Programming chemical kinetics: engineering dynamic reaction networks with DNA strand displacement

Over the last century, the silicon revolution has enabled us to build faster, smaller and more sophisticated computers. Today, these computers control phones, cars, satellites, assembly lines, and other electromechanical devices. Just as electrical wiring controls electromechanical devices, living organisms employ "chemical wiring" to make decisions about their environment and control physical processes. Currently, the big difference between these two substrates is that while we have the abstractions, design principles, verification and fabrication techniques in place for programming with silicon, we have no comparable understanding or expertise for programming chemistry.

In this thesis we take a small step towards the goal of learning how to systematically engineer prescribed non-equilibrium dynamical behaviors in chemical systems. We use the formalism of chemical reaction networks (CRNs), combined with mass-action kinetics, as our programming language for specifying dynamical behaviors. Leveraging the tools of nucleic acid nanotechnology (introduced in Chapter 1), we employ synthetic DNA molecules as our molecular architecture and toehold-mediated DNA strand displacement as our reaction primitive.

Abstraction, modular design and systematic fabrication can work only with well-understood and quantitatively characterized tools. Therefore, we embark on a detailed study of the "device physics" of DNA strand displacement (Chapter 2). We present a unified view of strand displacement biophysics and kinetics by studying the process at multiple levels of detail, using an intuitive model of a random walk on a 1-dimensional energy landscape, a secondary structure kinetics model with single base-pair steps, and a coarse-grained molecular model that incorporates three-dimensional geometric and steric effects. Further, we experimentally investigate the thermodynamics of three-way branch migration. Our findings are consistent with previously measured or inferred rates for hybridization, fraying, and branch migration, and provide a biophysical explanation of strand displacement kinetics. Our work paves the way for accurate modeling of strand displacement cascades, which would facilitate the simulation and construction of more complex molecular systems.

In Chapters 3 and 4, we identify and overcome the crucial experimental challenges involved in using our general DNA-based technology for engineering dynamical behaviors in the test tube. In this process, we identify important design rules that inform our choice of molecular motifs and our algorithms for designing and verifying DNA sequences for our molecular implementation. We also develop flexible molecular strategies for "tuning" our reaction rates and stoichiometries in order to compensate for unavoidable non-idealities in the molecular implementation, such as imperfectly synthesized molecules and spurious "leak" pathways that compete with desired pathways.

We successfully implement three distinct autocatalytic reactions, which we then combine into a de novo chemical oscillator. Unlike biological networks, which use sophisticated evolved molecules (like proteins) to realize such behavior, our test tube realization is the first to demonstrate that Watson-Crick base pairing interactions alone suffice for oscillatory dynamics. Since our design pipeline is general and applicable to any CRN, our experimental demonstration of a de novo chemical oscillator could enable the systematic construction of CRNs with other dynamic behaviors.

Probabilistic reasoning with a bayesian DNA device based on strand displacement

We present a computing model based on the DNA strand displacement technique which performs Bayesian inference. The model will take single stranded DNA as input data, representing the presence or absence of a specific molecular signal (evidence). The program logic encodes the prior probability of a disease and the conditional probability of a signal given the disease playing with a set of different DNA complexes and their ratios. When the input and program molecules interact, they release a different pair of single stranded DNA species whose relative proportion represents the application of Bayes? Law: the conditional probability of the disease given the signal. The models presented in this paper can empower the application of probabilistic reasoning in genetic diagnosis in vitro.

A coupled complex of T4 DNA replication helicase (gp41) and polymerase (gp43) can perform rapid and processive DNA strand-displacement synthesis

We have developed a coupled helicase–polymerase DNA unwinding assay and have used it to monitor the rate of double-stranded DNA unwinding catalyzed by the phage T4 DNA replication helicase (gp41). This procedure can be used to follow helicase activity in subpopulations in systems in which the unwinding-synthesis reaction is not synchronized on all the substrate-template molecules. We show that T4 replication helicase (gp41) and polymerase (gp43) can be assembled onto a loading site located near the end of a long double-stranded DNA template in the presence of a macromolecular crowding agent, and that this coupled “two-protein” system can carry out ATP-dependent strand displacement DNA synthesis at physiological rates (400 to 500 bp per sec) and with high processivity in the absence of other T4 DNA replication proteins. These results suggest that a direct helicase–polymerase interaction may be central to fast and processive double-stranded DNA replication, and lead us to reconsider the roles of the other replication proteins in processivity control.

RecA tests homology at both pairing and strand exchange

RecA is a 38-kDa protein from Escherichia coli that polymerizes on single-stranded DNA, forming a nucleoprotein filament that pairs with homologous duplex DNA and carries out strand exchange in vitro. To observe the effects of mismatches on the kinetics of the RecA-catalyzed recombination reaction, we used assays based upon fluorescence energy transfer that can differentiate between the pairing and strand displacement phases. Oligonucleotide sequences that produced 2–14% mismatches in the heteroduplex product of strand exchange were tested, as well as completely homologous and heterologous sequences. The equilibrium constant for pairing decreased as the number of mismatches increased, which appeared to result from both a decrease in the rate of formation and an increase in the rate of dissociation of the intermediates. In addition, the rate of strand displacement decreased with increasing numbers of mismatches, roughly in proportion to the number of mismatches. The equilibrium constant for pairing and the rate constant for strand displacement both decreased 6-fold as the heterology increased to 14%. These results suggest that discrimination of homology from heterology occurs during both pairing and strand exchange.

Pilot study of the utility and acceptability of tampon sampling for the diagnosis of Neisseria gonorrhoeae and Chlamydia trachomatis infections by duplex realtime polymerase chain reaction in United Kingdom sex workers

We aimed to evaluate the acceptability of self-collected tampon samples for the screening of female sex workers for sexually transmitted infections. We recruited 65 sex workers, and 63 agreed to provide tampon samples. The tampon samples were processed by realtime polymerase chain reaction (PCR) targeting Neisseria gonorrhoeae and Chlamydia trachomatis. Urethral and endocervical swabs were also obtained from 61 of 63 participants and tested using culture (N. gonorrhoeae) and the BD ProbeTec strand displacement amplification (SDA) (C. trachomatis) assay. Tampon sampling was preferred by 95% of the women and all favoured being tested away from genitourinary medicine clinics; the most common reasons cited were avoidance of embarrassment (40%) and convenience (30%). Besides near-universal acceptability of tampon sampling, the tampon sampling-PCR approach described in this study appeared to have enhanced sensitivity compared with conventional testing, suggesting the possibility of a residual hidden burden of N. gonorrhoeae and/or C. trachomatis genital infections in UK female sex workers.

Repair of abasic sites in DNA

Repair of both normal and reduced AP sites is activated by AP endonuclease, which recognizes and cleaves a phosphodiester bond 5' to the AP site. For a short period of time an incised AP site is occupied by poly(ADP-ribose) polymerase and then DNA polymerase beta adds one nucleotide into the repair gap and simultaneously removes the 5'-sugar phosphate. Finally, the DNA ligase III/XRCC1 complex accomplishes repair by sealing disrupted DNA ends. However, long-patch BER pathway, which is involved in the removal of reduced abasic sites, requires further DNA synthesis resulting in strand displacement and the generation of a damage-containing flap that is later removed by the flap endonuclease. Strand-displacement DNA synthesis is accomplished by DNA polymerase delta/epsilon and DNA ligase I restores DNA integrity. DNA synthesis by DNA polymerase delta/epsilon is dependent on proliferating cell nuclear antigen, which also stimulates the DNA ligase I and flap endonuclease. These repair events are supported by multiple protein-protein interactions. (C) 2003 Elsevier B.V. All rights reserved.

Vulvovaginal-swab or first-catch urine specimen to detect Chlamydia trachomatis in women in a community setting?

Screening for chlamydia in women is widely recommended. We evaluated the performance of two nucleic acid amplification tests for detecting Chlamydia trachomatis in self-collected vulvovaginal-swab and first-catch urine specimens from women in a community setting and a strategy for optimizing the sensitivity of an amplified enzyme immunoassay on vulvovaginal-swab specimens. We tested 2,745 paired vulvovaginal-swab and urine specimens by PCR (Roche Cobas) or strand displacement amplification (SDA; Becton Dickinson). There were 146 women infected with chlamydia. The assays detected 97.3% (95% confidence interval [CI], 93.1 to 99.2%) of infected patients with vulvovaginal-swab specimens and 91.8% (86.1 to 95.7%) with urine specimens. We tested 2,749 vulvovaginal-swab specimens with both a nucleic acid amplification test and a polymer conjugate-enhanced enzyme immunoassay with negative-gray-zone testing. The relative sensitivities obtained after retesting specimens in the negative gray zone were 74.3% (95% CI, 62.8 to 83.8%) with PCR and 58.3% (95% CI, 46.1 to 69.8%) with SDA. In community settings, both vulvovaginal-swab and first-catch urine specimens from women are suitable substrates for nucleic acid amplification tests, but enzyme immunoassays, even after negative-gray-zone testing, should not be used in screening programs.

Epidemiological, social, diagnostic and economic evaluation of population screening for genital chlamydial infection

OBJECTIVES: To investigate epidemiological, social, diagnostic and economic aspects of chlamydia screening in non-genitourinary medicine settings. METHODS: Linked studies around a cross-sectional population-based survey of adult men and women invited to collect urine and (for women) vulvovaginal swab specimens at home and mail these to a laboratory for testing for Chlamydia trachomatis. Specimens were used in laboratory evaluations of an amplified enzyme immunoassay (PCE EIA) and two nucleic acid amplification tests [Cobas polymerase chain reaction (PCR), Becton Dickinson strand displacement amplification (SDA)]. Chlamydia-positive cases and two negative controls completed a risk factor questionnaire. Chlamydia-positive cases were invited into a randomised controlled trial of partner notification strategies. Samples of individuals testing negative completed psychological questionnaires before and after screening. In-depth interviews were conducted at all stages of screening. Chlamydia transmission and cost-effectiveness of screening were investigated in a transmission dynamic model. SETTING AND PARTICIPANTS: General population in the Bristol and Birmingham areas of England. In total, 19,773 women and men aged 16-39 years were randomly selected from 27 general practice lists. RESULTS: Screening invitations reached 73% (14,382/19,773). Uptake (4731 participants), weighted for sampling, was 39.5% (95% CI 37.7, 40.8%) in women and 29.5% (95% CI 28.0, 31.0%) in men aged 16-39 years. Chlamydia prevalence (219 positive results) in 16-24 year olds was 6.2% (95% CI 4.9, 7.8%) in women and 5.3% (95% CI 4.4, 6.3%) in men. The case-control study did not identify any additional factors that would help target screening. Screening did not adversely affect anxiety, depression or self-esteem. Participants welcomed the convenience and privacy of home-sampling. The relative sensitivity of PCR on male urine specimens was 100% (95% CI 89.1, 100%). The combined relative sensitivities of PCR and SDA using female urine and vulvovaginal swabs were 91.8% (86.1, 95.7, 134/146) and 97.3% (93.1, 99.2%, 142/146). A total of 140 people (74% of eligible) participated in the randomised trial. Compared with referral to a genitourinary medicine clinic, partner notification by practice nurses resulted in 12.4% (95% CI -3.7, 28.6%) more patients with at least one partner treated and 22.0% (95% CI 6.1, 37.8%) more patients with all partners treated. The health service and patients costs (2005 prices) of home-based postal chlamydia screening were 21.47 pounds (95% CI 19.91 pounds, 25.99) per screening invitation and 28.56 pounds (95% CI 22.10 pounds, 30.43) per accepted offer. Preliminary modelling found an incremental cost-effectiveness ratio (2003 prices) comparing screening men and women annually to no screening in the base case of 27,000 pounds/major outcome averted at 8 years. If estimated screening uptake and pelvic inflammatory disease incidence were increased, the cost-effectiveness ratio fell to 3700 pounds/major outcome averted. CONCLUSIONS: Proactive screening for chlamydia in women and men using home-collected specimens was feasible and acceptable. Chlamydia prevalence rates in men and women in the general population are similar. Nucleic acid amplification tests can be used on first-catch urine specimens and vulvovaginal swabs. The administrative costs of proactive screening were similar to those for opportunistic screening. Using empirical estimates of screening uptake and incidence of complications, screening was not cost-effective.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.

Evolutionary conservation of residues in vertebrate DNA polymerase N conferring low fidelity and bypass activity.

POLN is a nuclear A-family DNA polymerase encoded in vertebrate genomes. POLN has unusual fidelity and DNA lesion bypass properties, including strong strand displacement activity, low fidelity favoring incorporation of T for template G and accurate translesion synthesis past a 5S-thymine glycol (5S-Tg). We searched for conserved features of the polymerase domain that distinguish it from prokaryotic pol I-type DNA polymerases. A Lys residue (679 in human POLN) of particular interest was identified in the conserved 'O-helix' of motif 4 in the fingers sub-domain. The corresponding residue is one of the most important for controlling fidelity of prokaryotic pol I and is a nonpolar Ala or Thr in those enzymes. Kinetic measurements show that K679A or K679T POLN mutant DNA polymerases have full activity on nondamaged templates, but poorly incorporate T opposite template G and do not bypass 5S-Tg efficiently. We also found that a conserved Tyr residue in the same motif not only affects sensitivity to dideoxynucleotides, but also greatly influences enzyme activity, fidelity and bypass. Protein sequence alignment reveals that POLN has three specific insertions in the DNA polymerase domain. The results demonstrate that residues have been strictly retained during evolution that confer unique bypass and fidelity properties on POLN.

3CAPS – a structural AP-site analogue as a tool to investigate DNA base excision repair

Abasic sites (AP-sites) are frequent DNA lesions, arising by spontaneous base hydrolysis or as intermediates of base excision repair (BER). The hemiacetal at the anomeric centre renders them chemically reactive, which presents a challenge to biochemical and structural investigation. Chemically more stable AP-site analogues have been used to avoid spontaneous decay, but these do not fully recapitulate the features of natural AP-sites. With its 3′-phosphate replaced by methylene, the abasic site analogue 3CAPS was suggested to circumvent some of these limitations. Here, we evaluated the properties of 3CAPS in biochemical BER assays with mammalian proteins. 3CAPS-containing DNA substrates were processed by APE1, albeit with comparably poor efficiency. APE1-cleaved 3CAPS can be extended by DNA polymerase β but repaired only by strand displacement as the 5′-deoxyribophosphate (dRP) cannot be removed. DNA glycosylases physically and functionally interact with 3CAPS substrates, underlining its structural integrity and biochemical reactivity. The AP lyase activity of bifunctional DNA glycosylases (NTH1, NEIL1, FPG), however, was fully inhibited. Notably, 3CAPS-containing DNA also effectively inhibited the activity of bifunctional glycosylases on authentic substrates. Hence, the chemically stable 3CAPS with its preserved hemiacetal functionality is a potent tool for BER research and a potential inhibitor of bifunctional DNA glycosylases.

Opening of a monomer-monomer interface of the trimeric bacteriophage T4-coded GP45 sliding clamp is required for clamp loading onto DNA

The replication system of bacteriophage T4 uses a trimeric ring-shaped processivity clamp (gp45) to tether the replication polymerase (gp43) to the template-primer DNA. This ring is placed onto the DNA by an ATPase-driven clamp-loading complex (gp44/62) where it then transfers, in closed form, to the polymerase. It generally has been assumed that one of the functions of the loading machinery is to open the clamp to place it around the DNA. However, the mechanism by which this occurs has not been fully defined. In this study we design and characterize a double-mutant gp45 protein that contains pairs of cysteine residues located at each monomer-monomer interface of the trimeric clamp. This mutant protein is functionally equivalent to wild-type gp45. However, when all three monomer-monomer interfaces are tethered by covalent crosslinks formed (reversibly or irreversibly) between the cysteine pairs these closed clamps can no longer be loaded onto the DNA nor onto the polymerase, effectively eliminating processive strand-displacement DNA synthesis. Analysis of the individual steps of the clamp-loading process shows that the ATPase-dependent interactions between the clamp and the clamp loader that precede DNA binding are hyperstimulated by the covalently crosslinked ring, suggesting that binding of the closed ring induces a futile, ATP-driven, ring-opening cycle. These findings and others permit further characterization and ordering of the steps involved in the T4 clamp-loading process.

A double-hexamer archaeal minichromosome maintenance protein is an ATP-dependent DNA helicase

The minichromosome maintenance (MCM) proteins are essential for DNA replication in eukaryotes. Thus far, all eukaryotes have been shown to contain six highly related MCMs that apparently function together in DNA replication. Sequencing of the entire genome of the thermophilic archaeon Methanobacterium thermoautotrophicum has allowed us to identify only a single MCM-like gene (ORF Mt1770). This gene is most similar to MCM4 in eukaryotic cells. Here we have expressed and purified the M. thermoautotrophicum MCM protein. The purified protein forms a complex that has a molecular mass of ≈850 kDa, consistent with formation of a double hexamer. The protein has an ATP-independent DNA-binding activity, a DNA-stimulated ATPase activity that discriminates between single- and double-stranded DNA, and a strand-displacement (helicase) activity that can unwind up to 500 base pairs. The 3′ to 5′ helicase activity requires both ATP hydrolysis and a functional nucleotide-binding site. Moreover, the double hexamer form is the active helicase. It is therefore likely that an MCM complex acts as the replicative DNA helicase in eukaryotes and archaea. The simplified replication machinery in archaea may provide a simplified model for assembly of the machinery required for initiation of eukaryotic DNA replication.

A quantitative method for evaluating the stabilities of nucleic acids

We report a general method for screening, in solution, the impact of deviations from canonical Watson-Crick composition on the thermodynamic stability of nucleic acid duplexes. We demonstrate how fluorescence resonance energy transfer (FRET) can be used to detect directly free energy differences between an initially formed “reference” duplex (usually a Watson-Crick duplex) and a related “test” duplex containing a lesion/alteration of interest (e.g., a mismatch, a modified, a deleted, or a bulged base, etc.). In one application, one titrates into a solution containing a fluorescently labeled, FRET-active, reference duplex, an unlabeled, single-stranded nucleic acid (test strand), which may or may not compete successfully to form a new duplex. When a new duplex forms by strand displacement, it will not exhibit FRET. The resultant titration curve (normalized fluorescence intensity vs. logarithm of test strand concentration) yields a value for the difference in stability (free energy) between the newly formed, test strand-containing duplex and the initial reference duplex. The use of competitive equilibria in this assay allows the measurement of equilibrium association constants that far exceed the magnitudes accessible by conventional titrimetric techniques. Additionally, because of the sensitivity of fluorescence, the method requires several orders of magnitude less material than most other solution methods. We discuss the advantages of this method for detecting and characterizing any modification that alters duplex stability, including, but not limited to, mutagenic lesions. We underscore the wide range of accessible free energy values that can be defined by this method, the applicability of the method in probing for a myriad of nucleic acid variations, such as single nucleotide polymorphisms, and the potential of the method for high throughput screening.

The 3′→5′ exonuclease of DNA polymerase δ can substitute for the 5′ flap endonuclease Rad27/Fen1 in processing Okazaki fragments and preventing genome instability

Many DNA polymerases (Pol) have an intrinsic 3′→5′ exonuclease (Exo) activity which corrects polymerase errors and prevents mutations. We describe a role of the 3′→5′ Exo of Pol δ as a supplement or backup for the Rad27/Fen1 5′ flap endonuclease. A yeast rad27 null allele was lethal in combination with Pol δ mutations in Exo I, Exo II, and Exo III motifs that inactivate its exonuclease, but it was viable with mutations in other parts of Pol δ. The rad27-p allele, which has little phenotypic effect by itself, was also lethal in combination with mutations in the Pol δ Exo I and Exo II motifs. However, rad27-p Pol δ Exo III double mutants were viable. They exhibited strong synergistic increases in CAN1 duplication mutations, intrachromosomal and interchromosomal recombination, and required the wild-type double-strand break repair genes RAD50, RAD51, and RAD52 for viability. Observed effects were similar to those of the rad27-null mutant deficient in the removal of 5′ flaps in the lagging strand. These results suggest that the 3′→5′ Exo activity of Pol δ is redundant with Rad27/Fen1 for creating ligatable nicks between adjacent Okazaki fragments, possibly by reducing the amount of strand-displacement in the lagging strand.