$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Voor een functionele test, een DNA-uitvoering van de bimoleculaire katalytische reactie (dwz, A + B-> B + C) is gemaakt. De prestaties van plasmide afgeleide poorten vergeleken met gates samengesteld uit synthetische DNA. Katalytische reacties zijn een goede test voor de poort zuiverheid omdat een defecte poort kan irreversibel vangen een katalysator, waardoor een onevenredig effect op de hoeveelheid geproduceerd product 18,19. Tegelijkertijd wordt een klein lek reactie resulteert in de untriggered afgifte van het katalytische signaal lineair versterkt, wat leidt tot een onevenredige foutsignaal. Experimentele gegevens voor plasmide afgeleid, gesynthetiseerd poorten zijn getoond in figuur 8B en 8C respectievelijk. In de experimenten wordt het gehalte van signaal onderdeel A vast terwijl de hoeveelheid van het katalytische signaal B wordt gevarieerd. Signaal C wordt gebruikt voor het uitlezen van de voortgang van de reactie zonder onderbreking van de katalytischecyclus. Katalyse kan worden waargenomen in de data sinds reacties naderen voltooiing zelfs hoeveelheden katalysator B veel kleiner dan de hoeveelheid A. Omdat SDS niet toegevoegd aan experimenten gedaan met het systeem gesynthetiseerd wordt reactiesnelheid (die kunnen worden beïnvloed door de toevoeging van SDS) niet vergeleken en de analytische nadruk in plaats daarvan op katalytische omzet (als hieronder beschreven).
Verdere analyse van de katalytische omzet van deze reactie werd uitgevoerd. De omzet is gedefinieerd als de hoeveelheid signaal C geproduceerd voor elke katalysator B op een bepaald moment. Specifiek, werd de omzet berekend op basis van onze experimentele gegevens die door het lek afgetrokken signaal C te delen door het aanvankelijke bedrag van katalysator B toegevoegd. Voor een ideale katalytisch systeem, moet deze omzet getal lineair toenemen met de tijd en onafhankelijk van de hoeveelheid katalysator zolang het substraat niet beperkend. In een echt systeem, kan defecte poorten kat uitschakelen alysts, en de omzet zal een maximale waarde, zelfs indien niet alle beschikbare substraat wordt omgezet naar het product te bereiken. De maximale omzet waarde geeft aan hoeveel substraten (signaal A) een katalysator (signaal B) kan omzetten voordat hij geïnactiveerd. Hier wordt opgemerkt dat de gesynthetiseerde systeem afwijkt van het ideale lineaire toename van de omzet veel eerder dan het plasmide afgeleide systeem doet, aangeeft vastlegging van de katalysator door een ongewenste nevenreactie (Figuur 8D). De omzet vergelijking is alleen zichtbaar voor lage concentraties, omdat bij hoge concentraties van katalysatoren, zullen alle poorten worden geactiveerd en laat het signaal C. Het circuit lekkage wordt ook vergeleken en wordt opgemerkt dat de verhouding leksignaal gebruik plasmide afgeleide gates ongeveer 8% lager dan die gebruikt gesynthetiseerde poorten na 10 h reactie (figuur 8E).
iles / ftp_upload / 53.087 / 53087fig1.jpg "/>
Figuur 1. (A) CRNs dienen als een normatief programmeertaal. DNA reactienetwerken kan worden gemanipuleerd om de dynamiek van een formele CRN benaderen (B) DNA uitvoering van een voorbeeld chemische instructie. A + B-> B + C. DNA-strengen worden getekend als lijnen met pijlen aan het 3 'uiteinde en * geeft complementariteit. Alle signaal strengen A (, groen), B (, oranje) en C (, red) worden bestond uit een steunpunt domein (gelabeld als ta, tb en tc) en een identiteit domein (aangeduid als a, b en c). De bimoleculaire reactie A + B-> B + C vereist twee multi-stranded complexen Doe AB en Vork voor Christus, en vier extra strengen , , en . De reactie verloopt door zeven stappen strengverplaatsing, waarbij elke stap starts met steunpunt bindend. (C) Reporter strategie. De reactie wordt gevolgd met behulp van een verslaggever die de onderste streng wordt gelabeld met een fluorofoor (rode stip) en de bovenste streng is een quencher (zwarte stip) bevestigd. Door de co-lokalisatie van de fluorofoor en de quencher, wordt de fluorescentie gedoofd reporter in het intacte reporter. Het signaal C kan de bovenste streng van de reporter te vervangen, wat leidt tot een toename van fluorescentie. (Dit cijfer is aangepast van 29 Ref.) Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2. (A) NdsDNA gates uit bacterieel plasmide DNA. Meerdere exemplaren van de dubbelstrengs ndsDNA poort template worden gekloond in een plasmide. De gekloneerde plasmiden zijnn getransformeerd in E. coli-cellen en kolonies op de plaat zijn sequentie geverifieerd. Nadat de sequentie is bevestigd, wordt plasmide-DNA geamplificeerd en gewonnen. Tenslotte wordt het dubbelstrengs plasmide verwerkt tot de gewenste ndsDNA poorten door middel van enzymatische behandeling. (B) Enzymatische behandeling van ndsDNA poorten. De PvuII restrictie-enzym wordt gebruikt om het hek te openen uit het plasmide. De vrijgekomen poorten zijn verder verwerkt met behulp van enzymen nicking: Nb.BsrDI wordt gebruikt om namen te genereren voor Word lid van AB (Panel i); Nt.BstNBI wordt gebruikt om nicks voor Fork BC (Panel ii) te genereren. Beperking en nicking sites worden aangeduid als kleurgecodeerde dozen. (C) Sequence uitzicht op de poort template van Join AB (Panel i) en Vork BC (Panel ii). Het PvuII restrictieplaats (gemarkeerd paarse doos) is aan beide uiteinden van de ndsDNA poorten. De Nb.BsrDI en NT.BstNBI Nicking sites zijn gemarkeerd in rood en zwarte dozen, respectievelijk. De locaties van de snede zijn gemarkeerd met pijlpunten. Sequence N elk nucleotide. (Dit cijfer is gewijzigd met toestemming van Ref. 29) Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3. (A) PCR van een DNA matrijs poort. Een DNA-poort template bevat de ndsDNA poort sequenties in het centrum (een blauwe regio) en spacer sequenties aan beide uiteinden (zwarte regio's; deze twee eindsequenties orthogonaal). Primers kunnen binden aan de spacersequenties van de template poort, en het genereren van vier overlappende DNA-fragmenten door middel van PCR (overlappende sequenties kleurcode in de figuur). (B) Gibson samenstel. De vier versterkte DNA Fragments worden vervolgens samengevoegd tot gelineariseerde plasmide ruggengraat door middel Gibson montagemethode 43. (Dit cijfer is gewijzigd met toestemming van Ref. 29) Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4. Circuit prestaties met verschillende hoeveelheden enzym. (A) een vereenvoudigde voorstelling van de poort reporter, hulp- draden en strengen signaal voor het corresponderende experimenten. (B) Kinetiek experimenten met plasmide afgeleide Word AB verwerkt met verschillende hoeveelheden enzym . ik. 10 eenheden PvuII-HF en 45 eenheden van Nb.BsrDI per 1 ug plasmide; ik ik. 10 eenheden PvuII-HF en 4 eenheden van Nb.BsrDI per 1 ug plasmide; iii. 4 eenheden PvuII-HF en 4 eenheden van Nb.BsrDI per 1 ug plasmide. Alle extra strengen waren 2x (1x = 10 nM). De poort complex was 1,5x, en de experimenten werden uitgevoerd bij 35 ° C in 1 x TAE / Mg2 +. (Dit cijfer is gewijzigd met toestemming van Ref. 29) Klik hier om een grotere versie van deze figuur te bekijken.

. Figuur 5 Flow grafiek van de kinetiek experimenten Blue. Materialen te voegen aan cuvette (0,875 ml synthetische kwarts cel). Referentietabel 14 specifieke volumes te voegen voor kinetische experiment A + B -> B + C. Groen: Aanwijzingen van een spectrofluorimeter (aangeduid als SPEX). Rood: Het mengen van instructies.53087fig5large.jpg "target =" _ blank "> Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6. Enzym dissociatie en circuit gedrag. (A) Een vereenvoudigde weergave van de poort reporter, hulp- draden en strengen signaal voor het corresponderende experimenten. (B) Kinetiek van de experimenten plasmide afgeleide uzelf AB gebruik 80 ° C verwarmen inactivatie (groene sporen), 0,15% natriumdodecylsulfaat (SDS) (rood), en een controle zonder warmte-inactivatie of toevoeging van SDS (blauw). De standaard concentratie 1x = 10 nM, en alle ondersteunende strengen en ingang B waren 2x. De poort complex was 1,5x, en de experimenten werden uitgevoerd bij 35 ° C in 1 x Tris-acetaat-EDTA-buffer bevattende 12,5 mM Mg2 + (1 x TAE / Mg2 ). (Dit cijfer is aangepast van 29 Ref.) Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7. Reporter kalibratie. (A) Reporter C kinetiek. De reporter concentratie aan 3x (1x = 50 nM), en de initiële concentratie van het signaal C wordt aangegeven in de figuur. (B) De fluorescentie niveaus van signaal C aan de meting eindpunt (40 min) toont een lineair verband met de beginconcentratie signaal C. In een kwantificering voorbeeld van Fork BC poort (groene stippellijn), de fluorescentie waarde van Fork BC was maatregeld als 3 x 10 6 (au), wat overeenkomt met 25 nm (0.5x) op basis van de kalibratiecurve. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 8. Bimoleculaire katalytische reactie kinetiek (A + B-> B + C). (A) Een vereenvoudigde weergave van de poort reporter, hulp- draden en strengen signaal voor het corresponderende experimenten. Experimenten werden uitgevoerd in 1 x Tris-acetaat-EDTA-buffer bevattende 12,5 mM Mg2 + (1 x TAE / Mg2 +). Alle gate complexen waren 75 nM concentratie (1,5 x) en ondersteunende strengen waren 100 nM concentratie (2x). Kinetics gegevens plasmide afgeleide poorten en gegevens gesynthetiseerde poorten zijn getoond in (B) en (C) , respectievelijk. Signal was 50 nM (1x). Verschillende hoeveelheden signaal (katalysator) werd ingevoerd in het systeem, en de reactie werd beproefd bij 35 ° C. (D) Plasmide-afgeleide poorten hadden hogere omzet dan gesynthetiseerd DNA poorten wanneer lage hoeveelheden input toegevoegd. (E) Omvang lekkage. Het diagram geeft de verhouding van de uiteindelijke lekkage naar het eindsignaal (C B0 = 0 / C B0 = 1) aan de eindpunten (10 uur). (Dit cijfer is aangepast van 29 Ref.) Klik hier om een grotere versie van deze figuur te bekijken.
| Poort Templates | Opeenvolgingen | Lengte (nt) |
| JoinAB | TCTAGTTCGATCAGAGCGTTATTACCAGTAGTCGATTGCTCAGCTGCTACATTGCTTCTACGAGTCATCCTTCCACCATTGCACCTTAGAGTCCGAATCCTACCATTGCTTAACCGAGTCTCACAACCAGCTGTCATTATGGACTTGACACACAGATTACACGGGAAAGTTGC | 173 |
| FORKBC | TCTAGTTCGATCAGAGCGTTATTACCAGTAGTCGATTGCTCAGCTGCCATCATAAGAGTCACCATACCCACATTGCCACATCGAGTCCCTTTTCCACCATTGCACCTTAGAGTCCGAATCCTACCATTGCTTAACCGAGTCTCACAACCAGCTGTCATTATGGACTTGACACACAGATTACACGGGAAAGTTGC | 194 |
Tabel 1. Sequences van ndsDNA poort templates.
| Poort | Streng | Lengte van de onderste streng (nt) |
| JoinAB | JoinAB-Bottom, , , | 87 |
| ForkBC | ForkBC-Bottom, , , , |
| 108 |
Tabel 2. Strands bestaande Doe AB en Vork voor Christus. (Deze tabel is gewijzigd van Ref. 29)
| Domein | Sequentie | Lengte (nt) |
| ta | CTGCTA | 6 |
| tb | TTCCAC | 6 |
| tc | TACCCA | 6 |
| tr | TCCTAC | 6 |
| tq | AACCAG | 6 |
| een | CATTGCTTCTACGAGTCATCC | 21 |
| b | CATTGCACCTTAGAGTCCGAA | 21 |
| c | CATTGCCACATCGAGTCCCTT | 21 |
| r | CATTGCTTAACCGAGTCTCAC | 21 |
| ik | CTGCCATCATAAGAGTCACCA | 21 |
jove_content "fo: keep-together.within-page =" always "fo: keep-met-previous.within-page =." always "> Tabel 3 Domein niveau sequenties voor de uitvoering van de chemische reactie A + B -> B + C . (Deze tabel is gewijzigd van Ref. 29) | Primer streng | Opeenvolgingen | Lengte (nt) |
| Forward primer-1 | AAGAGAGACCACATGGTCCTTCTTGAGTTTGTAACAG CGTTATTACCAGTAGTCGATTGC | 60 |
| Reverse primer-1 | ACTACTATTTACTAATCCCATTGCGTGTTCTTATT TAATCTGTGTGTCAAGTCCATAATG | 60 |
| Forward primer-2 | AATAAGAACACGCAATGGGATTAGTAAATAGTAGT CGTTATTACCAGTAGTCGATTGC | 58 |
| Reverse primer-2 | GCGAAACTAGCTTGTGGTGATATTGTCTCGTGTGT TAATCTGTGTGTCAAGTCCATAATG | 60 |
| Forward primer-3 | ACACACGAGACAATATCACCACAAGCTAGTTTCGC CGTTATTACCAGTAGTCGATTGC | 58 |
| Reverse primer-3 | ACATTGTACGCCTAAATCATCAAGAATAATTGTTG TAATCTGTGTGTCAAGTCCATAATG | 60 |
| Forward primer-4 | CAACAATTATTCTTGATGATTTAGGCGTACAATGT CGTTATTACCAGTAGTCGATTGC | 58 |
| Reverse primer-4 | GAGCGCAGCGAGTCAGTGAGCGAGGAAGCCTGCAG TAATCTGTGTGTCAAGTCCATAATG | 60 |
Tabel 4. Primer sequenties voor de PCR van ndsDNA gate templates.
| Reagens | Volume voor 1x reactie (pl) |
| High-copy plasmide backbone (~ 300 ng / ul) | 10 |
| PvuII-HF (20.000 eenheden / ml) | 2 |
| PstI-HF (20.000 eenheden / ml) | 2 |
| 10x Cut slimme buffer | 2 |
| H 2 O | 4 |
| Totaal volume | 20 d> |
Tabel 5. Protocol voor plasmide ruggengraat digest.
| Reagens | Volume voor 1x reactie (pl) |
| DNA vector (~ 50 ng / ul) | 1 |
| PCR geamplificeerde fragment-1 (~ 50 ng / ul) | 1 |
| PCR geamplificeerde fragment-2 (~ 50 ng / ul) | 1 |
| PCR geamplificeerde fragment-3 (~ 50 ng / ul) | 1 |
| PCR geamplificeerde fragment-4 (~ 50 ng / ul) | 1 |
| 2x Gibson Vergadering meester Mix | 5 |
| Totaal volume | 10 |
s "fo: keep-met-previous.within-page =" always ">
Tabel 6 Protocol voor Gibson montage.. | Reagens | Volume voor 1x reactie (pl) |
| Plasmide DNA (~ 1 ug / ul concentratie) | 1000 |
| PvuII-HF (20.000 eenheden / ml) | 200 |
| 10x Cut slimme buffer | 133,3 |
| Totaal volume | 1333,3 |
Tabel 7. Protocol voor ndsDNA hekken geplaatst plasmide verteren met restrictie-enzym PvuII-HF.
| Reagens | Volume (pl) |
| Join poorten (~ 5 ug / ul concentratie) | 150 |
| Nb.BsrDI (10.000 eenheden / ml) | 300 |
| 10x Cut slimme buffer | 50 |
| Totaal volume | 500 |
Tabel 8. Protocol mee poorten verteren met jatten enzym Nb.BsrDI.
| Reagens | Volume (pl) |
| Vork poorten (~ 5 ug / ul concentratie) | 150 |
| Nt.BstNBI (10.000 eenheden / ml) | 600 |
| 10x NEB buffer 3.1 | 83.3 |
| Totaal volume | 833,3 |
Tabel 9 Protocol voor de vork poorten verteren met jatten enzym Nt.BstNBI.
| Streng | Domein | Sequentie | Lengte (nt) |
| JoinAB-Bottom | tq * r * tr * b * tb * a * ta * | CTGGTT GTGAGACTCGGTTAAGCAATG GTAGGA TTCGGACTCTAAGGTGCAATG GTGGAA GGATGACTCGTAGAAGCAATG TAGCAG | 87 |
| FORKBC-Bottom | tq * r * tr * b * tb * c * tc * i * | CTGGTT GTGAGACTCGGTTAAGCAATG GTAGGA TTCGGACTCTAAGGTGCAATG GTGGAA AAGGGACTCGATGTGGCAATG TGGGTA TGGTGACTCTTATGATGGCAG | 108 |
| | | CTGCTA CATTGCTTCTACGAGTCATCC | 27 |
| | ta een | TTCCAC CATTGCACCTTAGAGTCCGAA | 27 |
| | tb b | TACCCA CATTGCCACATCGAGTCCCTT | 27 |
| | tc c | CATTGCTTCTACGAGTCATCC TTCCAC | 27 |
| | een tb | CATTGCACCTTAGAGTCCGAA TCCTAC | 27 |
| | b tr | CATTGCTTAACCGAGTCTCAC AACCAG | 27 |
| | r tq | CTGCCATCATAAGAGTCACCA | 21 |
| |
| ik | TCCTAC CATTGCTTAACCGAGTCTCAC | 27 |
| | tr r | CTGCCATCATAAGAGTCACCA TACCCA | 27 |
| | i tc | CATTGCCACATCGAGTCCCTT TCCTAC | 27 |
| | c tr | CATTGCCACATCGAGTCCCTT TTCCAC | 27 |
| | c tb | CATTGCACCTTAGAGTCCGAA TCCTAC | 27 |
| | b tr | CTGCCATCATAAGAGTCACCA TTCCAC | 27 |
| | i tb | CATTGCACCTTAGAGTCCGAA TTCCAC | 27 |
| | b tb | CATTGCACCTTAGAGTCCGAA TACCCA | 27 |
| | b tc | CATTGCCACATCGAGTCCCTT TCCTAC | 27 |
| | c tr | CATTGCACCTTAGAGTCCGAA TCCTAC | 27 |
| ROX- | b tr | / 56-ROXN / AAGGGACTCGATGTGGCAATG TGGGTA | 27 |
| -RQ | c * tc * | CATTGCCACATCGAGTCCCTT / 3IAbRQSp / | 21 |
| - TAMRA | c | CTGGTT GTGAGACTCGGTTAAGCAATG / 36-TAMTSp / | 27 |
| RQ- | tq * r * | / 5IAbRQ / CATTGCTTAACCGAGTCTCAC | 21 |
| r | | |
. Tabel 10 Strand sequenties voor de uitvoering van de chemische reactie A + B -.> B + C (. Deze tabel is gewijzigd van Ref 29)
| Reagens | Volume (pl) | Eindconcentratie |
| ROX- 100 uM | 10 | 10 uM (1x) |
| -RQ bij 100 μ; M | 13 | 13 uM (1.3x) |
| 10x TAE met 125 mm Mg 2+ | 10 | 1x TAE met 12,5 mM Mg 2+ |
| H 2 O | 67 | - |
| Totaal volume | 100 | 10 uM (1x) |
Tabel 11. Protocol voor het assembleren Reporter C.
| Reagens | Volume (pl) | Eindconcentratie |
| H 2 O | 514 | - |
| 10x TAE met 125 mm Mg 2+ | 60 | 1x TAE met 12,5 mM Mg 2+ |
| PolyT bij 300 μ; M | 2 | 1 uM |
| Reporter C bij 10 uM | 9 | 150 nM (3x) |
| 10% SDS | 9 | 0,15% |
| 5 uM | 6 | 50 Nm (1x) |
| Totaal volume | 600 | - |
Tabel 12. Protocol voor het kalibreren van Reporter C. De volumes die hier voor een totaal reactievolume van 600 pi (overeenkomend met het gebruik van een 0,875 ml synthetisch kwarts cel), maar kan worden aangepast om te werken met verschillende grootte cellen.
| Reagens | Volume (pl) | Final concentratie |
| H 2 O | 493 | - |
| 10x TAE met 125 mm Mg 2+ | 60 | 1x TAE met 12,5 mM Mg 2+ |
| polyT bij 300 uM | 2 | 1 uM |
| Reporter C bij 10 uM | 9 | 150 nM (3x) |
| bij 100 uM | 3 | 10x |
| bij 100 uM | 3 | 10x |
| 100 uM | 3 | 10x |
| 10% SDS | 9 | 0,15% |
| Vork BC bij ~ 1 uM (concentratie onbekend) | 15 | ~ 0,5x |
| bij 100 uM | 3 | 10x |
| Totaal volume | 600 | - |
Tabel 13. Protocol voor het kalibreren van Fork BC. De volumes die hier voor een totaal reactievolume van 600 pi, maar kan worden aangepast om te werken met verschillende grootte cellen.
| Reagens | Volume (pl) | Eindconcentratie |
| H 2 O | 407,2 | - | |
| 10x TAE met 125 mm Mg 2+ | 52.8 | 12,5 mM Mg2 | |
| polyT bij 300 uM | 2 | 1 uM | |
| Reporter C bij 10 uM | 9 | 150 nM (3x) | |
| bij 10 uM | 6 | 100 nM (2x) | |
| bij 10 uM | 6 | 100 nM (2x) | |
| 10 uM | 6 | 100 nM (2x) | |
| |
| 6 | 100 nM (2x) | |
| 10% SDS | 9 | 0,15% | |
| Join AB 1 uM | 45 | 75 nm (1,5x) | |
| Vork BC 1 uM | 45 | 75 nm (1,5x) | |
| bij 10 uM | 3 | 50 Nm (1x) | |
| 10 uM | 3 | 50 Nm (1x) | |
| Totaal volume | 600 | - | |
bij 10 uM
Tabel 14. Protocol voor een chemische reactie A + B-> B + C. De volumes die hier voor een totaal reactievolume van 600 pi, maar kan worden aangepast om te werken met verschillende grootte cellen.
| Gesynthetiseerd poorten | | Plasmide afgeleid poorten | | | |
| Beschrijving | Kosten | Join poorten | | Vork poorten | |
| PAGE gezuiverde lange streng (100 nt, diende als de bodem strengen van een poort) | ~ $ 75 | Beschrijving </ strong> | Kosten | Beschrijving | Kosten |
| PAGE gezuiverde korte streng (~ 30 nt, diende als top strengen van een poort) | ~ $ 185 | Poort template | ~ $ 100 | Poort template | ~ $ 100 |
| Totaal | ~ $ 260 | Plasmide extractiekit | ~ $ 26 | Plasmide extractiekit | ~ $ 26 |
| | Restrictie-enzym (PvuII-HF) | ~ $ 11 | Restrictie-enzym (PvuII-HF) | ~ $ 11 |
| | Inkeping enzym (Nt.BsrDI, Join poorten) | ~ $ 29 | Inkeping enzym (Nt.BstNBI, Fork poorten) | ~ $ 62 |
| | Totaal | ~ $ 166 | Totaal | ~ $ 199 |
fo:.. keep-met-previous.within-page = "always"> Tabel 15 Kosten vergelijking tussen plasmide afgeleid poorten en gesynthetiseerd poorten (. Deze tabel is gewijzigd van Ref 29)
| Gesynthetiseerd poorten | Plasmide afgeleid poorten | |
| Verwerking | Verwerkingstijd | Verwerking | Verwerkingstijd |
| Gloeien | 1 hr | Klonen | 5 uur |
| PAGE zuivering | 2 hr | Plasmide-extractie | 2 hr |
| Totaal | 3 hr | Twee stappen enzymdigestie | 0,5 hr |
| | Ethanolprecipitatie | 1 hr |
| | Totaal | 8,5 hr |
Tabel 16. Verwerking tijd vergelijking tussen plasmide afgeleid poorten en synthetische poorten. (Deze tabel is gewijzigd van Ref. 29)