Hydrogen bonding revisited: geometric selection as a principal determinant of DNA replication fidelity

A method to measure molecular hybridization

Fluorescence-based oligonucleotide probes have a great importance in research of molecular interactions. Molecular beacons (MBs) are special case of fluorescent probes that form a stem-loop shape, bringing together a fluorophore and quencher, thus emitting fluorescence only when hybridized to a complementary target. Here we describe a new method for the quantitation of MB hybridization based on the measurement of changes in free energy instead of the fluorescence intensity. The MB energy state can be measured by micro-fluorescence detection. The approach allowed to determine hybridization energy of the MB with target nucleotide directly from fluorescence spectra and distinguish the MB in unfolded and hybridized states. Moreover, the method enabled us to discriminate between DNA duplexes with perfect complementarity or a single-nucleotide mismatch, based on the first direct experimental prove of enthalpy-entropy compensation.

Base pair dynamics, electrostatics, and thermodynamics at the LTR-III quadruplex:duplex junction

G-quadruplexes (GQs) play key regulatory roles within the human genome and have also been identified to play similar roles in other eukaryotes, bacteria, archaea, and viruses. Human immunodeficiency virus 1, the etiological agent of acquired immunodeficiency syndrome, can form two GQs in its long terminal repeat (LTR) promoter region, each of which act to regulate viral gene expression in opposing manners. The major LTR GQ, called LTR-III, is a distinct hybrid GQ containing a 12-nucleotide duplex loop attached to the quadruplex motif. The resulting quadruplex:duplex junction (QDJ) has been hypothesized to serve as a selective drug targeting site. To better understand the dynamics of this QDJ, we performed conventional and enhanced-sampling molecular dynamics simulations using the Drude-2017 force field. We observed unbiased and reversible formation of additional base pairs in the QDJ, between Ade4:Thy14 and Gua3:Thy14. Both base pairs were electrostatically favored, but geometric constraints within the junction may drive the formation of, and preference for, the Ade4:Thy14 base pair. Finally, we demonstrated that the base pairs are separated only by small energy barriers that may enable transitions between both base-paired states. Together, these simulations provide new insights into the dynamics, electrostatics, and thermodynamics of the LTR-III QDJ.

Hydrogen bonding in duplex DNA probed by DNP enhanced solid-state NMR N-H bond length measurements

Numerous biological processes and mechanisms depend on details of base pairing and hydrogen bonding in DNA. Hydrogen bonds are challenging to quantify by X-ray crystallography and cryo-EM due to difficulty of visualizing hydrogen atom locations but can be probed with site specificity by NMR spectroscopy in solution and the solid state with the latter particularly suited to large, slowly tumbling DNA complexes. Recently, we showed that low-temperature dynamic nuclear polarization (DNP) enhanced solid-state NMR is a valuable tool for distinguishing Hoogsteen base pairs (bps) from canonical Watson-Crick bps in various DNA systems under native-like conditions. Here, using a model 12-mer DNA duplex containing two central adenine-thymine (A-T) bps in either Watson-Crick or Hoogsteen confirmation, we demonstrate DNP solid-state NMR measurements of thymine N3-H3 bond lengths, which are sensitive to details of N-H···N hydrogen bonding and permit hydrogen bonds for the two bp conformers to be systematically compared within the same DNA sequence context. For this DNA duplex, effectively identical TN3-H3 bond lengths of 1.055 ± 0.011 Å and 1.060 ± 0.011 Å were found for Watson-Crick A-T and Hoogsteen A (syn)-T base pairs, respectively, relative to a reference amide bond length of 1.015 ± 0.010 Å determined for N-acetyl-valine under comparable experimental conditions. Considering that prior quantum chemical calculations which account for zero-point motions predict a somewhat longer effective peptide N-H bond length of 1.041 Å, in agreement with solution and solid-state NMR studies of peptides and proteins at ambient temperature, to facilitate direct comparisons with these earlier studies TN3-H3 bond lengths for the DNA samples can be readily scaled appropriately to yield 1.083 Å and 1.087 Å for Watson-Crick A-T and Hoogsteen A (syn)-T bps, respectively, relative to the 1.041 Å reference peptide N-H bond length. Remarkably, in the context of the model DNA duplex, these results indicate that there are no significant differences in N-H···N A-T hydrogen bonds between Watson-Crick and Hoogsteen bp conformers. More generally, high precision measurements of N-H bond lengths by low-temperature DNP solid-state NMR based methods are expected to facilitate detailed comparative analysis of hydrogen bonding for a range of DNA complexes and base pairing environments.

Measuring the Enthalpy of an Individual Hydrogen Bond in a DNA Duplex with Nucleobase Isotope Editing and Variable-Temperature Infrared Spectroscopy

The level of interest in probing the strength of noncovalent interactions in DNA duplexes is high, as these weak forces dictate the range of suprastructures the double helix adopts under different conditions, in turn directly impacting the biological functions and industrial applications of duplexes that require making and breaking them to access the genetic code. However, few experimental tools can measure these weak forces embedded within large biological suprastructures in the native solution environment. Here, we develop experimental methods for detecting the presence of a single noncovalent interaction [a hydrogen bond (H-bond)] within a large DNA duplex in solution and measure its formation enthalpy (ΔH_f). We report that introduction of a H-bond into the TC2═O group from the noncanonical nucleobase 2-aminopurine produces an expected decrease ∼10 ± 0.76 cm^-1 (from ∼1720 cm^-1 in Watson-Crick to ∼1710 cm^-1 in 2-aminopurine), which correlates with an enthalpy of ∼0.93 ± 0.066 kcal/mol for this interaction.

The origin of genetic and metabolic systems: Evolutionary structuralinsights

DNA is derived from reverse transcription and its origin is related to reverse transcriptase, DNA polymerase and integrase. The gene structure originated from the evolution of the first RNA polymerase. Thus, an explanation of the origin of the genetic system must also explain the evolution of these enzymes. This paper proposes a polymer structure model, termed the stable complex evolution model, which explains the evolution of enzymes and functional molecules. Enzymes evolved their functions by forming locally tightly packed complexes with specific substrates. A metabolic reaction can therefore be considered to be the result of adaptive evolution in this way when a certain essential molecule is lacking in a cell. The evolution of the primitive genetic and metabolic systems was thus coordinated and synchronized. According to the stable complex model, almost all functional molecules establish binding affinity and specific recognition through complementary interactions, and functional molecules therefore have the nature of being auto-reactive. This is thermodynamically favorable and leads to functional duplication and self-organization. Therefore, it can be speculated that biological systems have a certain tendency to maintain functional stability or are influenced by an inherent selective power. The evolution of dormant bacteria may support this hypothesis, and inherent selectivity can be unified with natural selection at the molecular level.

Discovery, implications and initial use of semi-synthetic organisms with an expanded genetic alphabet/code

Much recent interest has focused on developing proteins for human use, such as in medicine. However, natural proteins are made up of only a limited number of canonical amino acids with limited functionalities, and this makes the discovery of variants with some functions difficult. The ability to recombinantly express proteins containing non-canonical amino acids (ncAAs) with properties selected to impart the protein with desired properties is expected to dramatically improve the discovery of proteins with different functions. Perhaps the most straightforward approach to such an expansion of the genetic code is through expansion of the genetic alphabet, so that new codon/anticodon pairs can be created to assign to ncAAs. In this review, I briefly summarize more than 20 years of effort leading ultimately to the discovery of synthetic nucleotides that pair to form an unnatural base pair, which when incorporated into DNA, is stably maintained, transcribed and used to translate proteins in Escherichia coli. In addition to discussing wide ranging conceptual implications, I also describe ongoing efforts at the pharmaceutical company Sanofi to employ the resulting 'semi-synthetic organisms' or SSOs, for the production of next-generation protein therapeutics. This article is part of the theme issue 'Reactivity and mechanism in chemical and synthetic biology'.

Kinetic Proofreading

Biochemistry and molecular biology rely on the recognition of structural complementarity between molecules. Molecular interactions must be both quickly reversible, i.e., tenuous, and specific. How the cell reconciles these conflicting demands is the subject of this article. The problem and its theoretical solution are discussed within the wider theoretical context of the thermodynamics of stochastic processes (stochastic thermodynamics). The solution-an irreversible reaction cycle that decreases internal error at the expense of entropy export into the environment-is shown to be widely employed by biological processes that transmit genetic and regulatory information. Expected final online publication date for the Annual Review of Biochemistry, Volume 91 is June 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Conformational dynamics during misincorporation and mismatch extension defined using a DNA polymerase with a fluorescent artificial amino acid

High-fidelity DNA polymerases select the correct nucleotide over the structurally similar incorrect nucleotides with extremely high specificity while maintaining fast rates of incorporation. Previous analysis revealed the conformational dynamics and complete kinetic pathway governing correct nucleotide incorporation using a high-fidelity DNA polymerase variant containing a fluorescent unnatural amino acid. Here we extend this analysis to investigate the kinetics of nucleotide misincorporation and mismatch extension. We report the specificity constants for all possible misincorporations and characterize the conformational dynamics of the enzyme during misincorporation and mismatch extension. We present free energy profiles based on the kinetic measurements and discuss the effect of different steps on specificity. During mismatch incorporation and subsequent extension with the correct nucleotide, the rates of the conformational change and chemistry are both greatly reduced. The nucleotide dissociation rate, however, increases to exceed the rate of chemistry. To investigate the structural basis for discrimination against mismatched nucleotides, we performed all atom molecular dynamics simulations on complexes with either the correct or mismatched nucleotide bound at the polymerase active site. The simulations suggest that the closed form of the enzyme with a mismatch bound is greatly destabilized due to weaker interactions with active site residues, nonideal base pairing, and a large increase in the distance from the 3'-OH group of the primer strand to the α-phosphate of the incoming nucleotide, explaining the reduced rates of misincorporation. The observed kinetic and structural mechanisms governing nucleotide misincorporation reveal the general principles likely applicable to other high-fidelity DNA polymerases.

Thermodynamic uncertainty relation to assess biological processes

We review the trade-offs between speed, fluctuations, and thermodynamic cost involved with biological processes in nonequilibrium states and discuss how optimal these processes are in light of the universal bound set by the thermodynamic uncertainty relation (TUR). The values of the uncertainty product Q of TUR, which can be used as a measure of the precision of enzymatic processes realized for a given thermodynamic cost, are suboptimal when the substrate concentration is at the Michaelis constant, and some of the key biological processes are found to work around this condition. We illustrate the utility of Q in assessing how close the molecular motors and biomass producing machineries are to the TUR bound, and for the cases of biomass production (or biological copying processes), we discuss how their optimality quantified in terms of Q is balanced with the error rate in the information transfer process. We also touch upon the trade-offs in other error-minimizing processes in biology, such as gene regulation and chaperone-assisted protein folding. A spectrum of Q recapitulating the biological processes surveyed here provides glimpses into how biological systems are evolved to optimize and balance the conflicting functional requirements.

Interplay between Position-Dependent Codon Usage Bias and Hydrogen Bonding at the 5' End of ORFeomes

Codon usage bias exerts control over a wide variety of molecular processes. The positioning of synonymous codons within coding sequences (CDSs) dictates protein expression by mechanisms such as local translation efficiency, mRNA Gibbs free energy, and protein cotranslational folding. In this work, we explore how codon usage affects the position-dependent content of hydrogen bonding, which in turn influences energy requirements for unwinding double-stranded DNA (dsDNA). We categorized codons according to their hydrogen bond content and found differential effects on hydrogen bonding encoded by codon variants. The specific positional disposition of codon variants within CDSs creates a ramp of hydrogen bonding at the 5' end of the ORFeome in Escherichia coli CDSs occupying the first position of operons are subjected to selective pressure that reduces their hydrogen bonding compared to internal CDSs, and highly transcribed CDSs demand a lower maximum capacity of hydrogen bonds per codon, suggesting that the energetic requirement for unwinding the dsDNA in highly transcribed CDSs has evolved to be minimized in E. coli Subsequent analysis of over 14,000 ORFeomes showed a pervasive ramp of hydrogen bonding at the 5' end in Bacteria and Archaea that positively correlates with the probability of mRNA secondary structure formation. Both the ramp and the correlation were not found in Fungi The position-dependent hydrogen bonding might be part of the mechanism that contributes to the coordination between transcription and translation in Bacteria and Archaea A Web-based application to analyze the position-dependent hydrogen bonding of ORFeomes has been developed and is publicly available (https://juanvillada.shinyapps.io/hbonds/).IMPORTANCE Redundancy of the genetic code creates a vast space of alternatives to encode a protein. Synonymous codons exert control over a variety of molecular and physiological processes of cells mainly through influencing protein biosynthesis. Recent findings have shown that synonymous codon choice affects transcription by controlling mRNA abundance, mRNA stability, transcription termination, and transcript biosynthesis cost. In this work, by analyzing thousands of Bacteria, Archaea, and Fungi genomes, we extend recent findings by showing that synonymous codon choice, corresponding to the number of hydrogen bonds in a codon, can also have an effect on the energetic requirements for unwinding double-stranded DNA in a position-dependent fashion. This report offers new perspectives on the mechanism behind the transcription-translation coordination and complements previous hypotheses on the resource allocation strategies used by Bacteria and Archaea to manage energy efficiency in gene expression.

A Deep Dive into DNA Base Pairing Interactions Under Water

Base pairing plays a pivotal role in DNA functions and replication fidelity. But while the complementarity between Watson-Crick matched bases is generally believed to arise from the different number of hydrogen bonds in G|C pairs versus A|T, the energetics of these interactions are heavily renormalized by the aqueous solvent. Employing large-scale Monte Carlo simulations, we have extracted the solvent contribution to the free energy for canonical and some noncanonical and stacked base pairs. For all of them, the solvent's contribution to the base pairing free energy is exclusively destabilizing. While the direct hydrogen bonding interactions in the G|C pair is much stronger than A|T, the thermodynamic resistance produced by the solvent also pushes back much stronger against G|C compared to A|T, generating an only ∼1 kcal/mol free energy difference between them. We have profiled the density of water molecules in the solvent adjacent to the bases and observed a "freezing" behavior where waters are recruited into the gap between the bases to compensate for the unsatisfied hydrogen bonds between them. A very small number of water molecules that are associated with the Watson-Crick donor/acceptor atoms turn out to be responsible for the majority of the solvent's thermodynamic resistance to base pairing. The absence or presence of these near-field waters can be used to enhance fidelity during DNA replication.

Thermodynamic Cost, Speed, Fluctuations, and Error Reduction of Biological Copy Machines

Due to large fluctuations in cellular environments, transfer of information in biological processes without regulation is error-prone. The mechanistic details of error-reducing mechanisms in biological copying processes have been a subject of active research; however, how error reduction of a process is balanced with its thermodynamic cost and dynamical properties remain largely unexplored. Here, we study the error reducing strategies in light of the recently discovered thermodynamic uncertainty relation (TUR) that sets a physical bound to the cost-precision trade-off for dissipative processes. We found that the two representative copying processes, DNA replication by the exonuclease-deficient T7 DNA polymerase and mRNA translation by the E. coli ribosome, reduce the error rates to biologically acceptable levels while also optimizing the processes close to the physical limit dictated by TUR.

On the Cooperativity Effect in Watson and Crick and Wobble Pairs for a Halouracil Series and Its Potential Quantitative Application Studied through Surface-Enhanced Raman Spectroscopy

The nature of the cooperativity effect of hydrogen bonds in Watson and Crick and wobble base pairs formed with thymine, uracil, and its 5-halogenated derivatives (5-fluoro, -chloro, and -bromouracil) has been studied through SERS and by using chemometric tools to process data and extract relevant information. Remarkable differences between the two kinds of pairs were clearly observed, and the behavior correlated to the withdrawing character of different substituents at the 5-position of uracil was verified. Multivariate analyses have also unveiled information about the pair's stability, and a stronger cooperativity effect seems to rule the Watson and Crick pairs when compared to wobble pairs. Defined patterns in the behavior of Watson and Crick pairs allowed the design of an indirect methodology for quantifying 5-bromouracil using a partial least squares (PLS) method with variable selection. Limit of detection (LOD) values of 0.037 and 0.112 mmol L^-1 in the absence and presence of structurally similar interferences were reached, while its direct surface-enhanced Raman spectroscopy (SERS) quantification is only possible at ∼45 mmol L^-1.

Local energy decomposition analysis of hydrogen-bonded dimers within a domain-based pair natural orbital coupled cluster study

The local energy decomposition (LED) analysis allows for a decomposition of the accurate domain-based local pair natural orbital CCSD(T) [DLPNO-CCSD(T)] energy into physically meaningful contributions including geometric and electronic preparation, electrostatic interaction, interfragment exchange, dynamic charge polarization, and London dispersion terms. Herein, this technique is employed in the study of hydrogen-bonding interactions in a series of conformers of water and hydrogen fluoride dimers. Initially, DLPNO-CCSD(T) dissociation energies for the most stable conformers are computed and compared with available experimental data. Afterwards, the decay of the LED terms with the intermolecular distance (r) is discussed and results are compared with the ones obtained from the popular symmetry adapted perturbation theory (SAPT). It is found that, as expected, electrostatic contributions slowly decay for increasing r and dominate the interaction energies in the long range. London dispersion contributions decay as expected, as r-6. They significantly affect the depths of the potential wells. The interfragment exchange provides a further stabilizing contribution that decays exponentially with the intermolecular distance. This information is used to rationalize the trend of stability of various conformers of the water and hydrogen fluoride dimers.

Variants of sequence family B Thermococcus kodakaraensis DNA polymerase with increased mismatch extension selectivity

Fidelity and selectivity of DNA polymerases are critical determinants for the biology of life, as well as important tools for biotechnological applications. DNA polymerases catalyze the formation of DNA strands by adding deoxynucleotides to a primer, which is complementarily bound to a template. To ensure the integrity of the genome, DNA polymerases select the correct nucleotide and further extend the nascent DNA strand. Thus, DNA polymerase fidelity is pivotal for ensuring that cells can replicate their genome with minimal error. DNA polymerases are, however, further optimized for more specific biotechnological or diagnostic applications. Here we report on the semi-rational design of mutant libraries derived by saturation mutagenesis at single sites of a 3’-5’-exonuclease deficient variant of Thermococcus kodakaraensis DNA polymerase (KOD pol) and the discovery for variants with enhanced mismatch extension selectivity by screening. Sites of potential interest for saturation mutagenesis were selected by their proximity to primer or template strands. The resulting libraries were screened via quantitative real-time PCR. We identified three variants with single amino acid exchanges—R501C, R606Q, and R606W—which exhibited increased mismatch extension selectivity. These variants were further characterized towards their potential in mismatch discrimination. Additionally, the identified enzymes were also able to differentiate between cytosine and 5-methylcytosine. Our results demonstrate the potential in characterizing and developing DNA polymerases for specific PCR based applications in DNA biotechnology and diagnostics.

Proofreading of DNA polymerase: a new kinetic model with higher-order terminal effects

The fidelity of DNA replication by DNA polymerase (DNAP) has long been an important issue in biology. While numerous experiments have revealed details of the molecular structure and working mechanism of DNAP which consists of both a polymerase site and an exonuclease (proofreading) site, there were quite a few theoretical studies on the fidelity issue. The first model which explicitly considered both sites was proposed in the 1970s and the basic idea was widely accepted by later models. However, all these models did not systematically investigate the dominant factor on DNAP fidelity, i.e. the higher-order terminal effects through which the polymerization pathway and the proofreading pathway coordinate to achieve high fidelity. In this paper, we propose a new and comprehensive kinetic model of DNAP based on some recent experimental observations, which includes previous models as special cases. We present a rigorous and unified treatment of the corresponding steady-state kinetic equations of any-order terminal effects, and derive analytical expressions for fidelity in terms of kinetic parameters under bio-relevant conditions. These expressions offer new insights on how the higher-order terminal effects contribute substantially to the fidelity in an order-by-order way, and also show that the polymerization-and-proofreading mechanism is dominated only by very few key parameters. We then apply these results to calculate the fidelity of some real DNAPs, which are in good agreements with previous intuitive estimates given by experimentalists.

Controlled keto–enol tautomerism of coumarin containing β-ketodithioester by its encapsulation in cucurbit[7]uril

CB7 shifts the tautomeric equilibrium of CAM2, from the enol- to keto-form, whereas β-CD maintains the enol form.

DNA and RNA folds in transcription complex as evidenced by iodine-125 radioprobing

Abstract Folding of the DNA and RNA strands in an arrested T7 RNA polymerase (RNAP) transcription complex was studied by radioprobing, a novel method based on an analysis of the strand breaks produced by decay of the iodine-125 incorporated at the C5 position of cytosine. (125)I-labeled cytosines were incorporated into transcripts at different positions relative to the site of the arrest. The intensities of the DNA breaks inversely correlate with the distances from the (125)I decay site, so the radioprobing data provide information about the spatial RNA/DNA folding during transcription. We found that the yield of DNA strand breaks is significantly higher in the template than the non-template strand. This is consistent with local opening of the DNA duplex and formation of a hybrid between RNA and the template DNA strand. Our data demonstrate that the RNA-DNA hybrid has a nonuni form A-like structure. When the (125)I is incorporatedseven nucleotides from the active center of RNAP, the yield of strand breaks is substantially lower than if (125)I is positioned at the ends of the hybrid. Consequently, the DNA and RNA strands are located closer to each other at the ends of the hybrid and somehow separated in the middle. Surprisingly, the (125)I-induced breaks were detected in both DNA strands upstream from the transcription "bubble" indicating that DNA and RNA are closely associated outside the RNAP cleft. Thus, radioprobing data imply that the RNA/DNA fold in the complex with T7 RNAP is more complicated than had been anticipated by the existing models. Based on the present data, we suggest a sterically feasible model explaining how formation of the long RNA-DNA hybrid can result in the initiation-to-elongation switch in the T7 transcription complex. According to this model, the topological linkage between the DNA and RNA strands provides the necessary stability for the elongation complex, while permitting movement of the polymerase along the DNA duplex.

Kinetic selection vs. free energy of DNA base pairing in control of polymerase fidelity

What is the free energy source enabling high-fidelity DNA polymerases (pols) to favor incorporation of correct over incorrect base pairs by 10(3)- to 10(4)-fold, corresponding to free energy differences of ΔΔGinc∼ 5.5-7 kcal/mol? Standard ΔΔG° values (∼0.3 kcal/mol) calculated from melting temperature measurements comparing matched vs. mismatched base pairs at duplex DNA termini are far too low to explain pol accuracy. Earlier analyses suggested that pol active-site steric constraints can amplify DNA free energy differences at the transition state (kinetic selection). A recent paper [Olson et al. (2013)J Am Chem Soc135:1205-1208] used Vent pol to catalyze incorporations in the presence of inorganic pyrophosphate intended to equilibrate forward (polymerization) and backward (pyrophosphorolysis) reactions. A steady-state leveling off of incorporation profiles at long reaction times was interpreted as reaching equilibrium between polymerization and pyrophosphorolysis, yielding apparent ΔG° = -RTlnKeq, indicating ΔΔG° of 3.5-7 kcal/mol, sufficient to account for pol accuracy without need of kinetic selection. Here we perform experiments to measure and account for pyrophosphorolysis explicitly. We show that forward and reverse reactions attain steady states far from equilibrium for wrong incorporations such as G opposite T. Therefore,[Formula: see text]values obtained from such steady-state evaluations ofKeqare not dependent on DNA properties alone, but depend largely on constraints imposed on right and wrong substrates in the polymerase active site.

The expanded genetic alphabet

All biological information, since the last common ancestor of all life on Earth, has been encoded by a genetic alphabet consisting of only four nucleotides that form two base pairs. Long-standing efforts to develop two synthetic nucleotides that form a third, unnatural base pair (UBP) have recently yielded three promising candidates, one based on alternative hydrogen bonding, and two based on hydrophobic and packing forces. All three of these UBPs are replicated and transcribed with remarkable efficiency and fidelity, and the latter two thus demonstrate that hydrogen bonding is not unique in its ability to underlie the storage and retrieval of genetic information. This Review highlights these recent developments as well as the applications enabled by the UBPs, including the expansion of the evolution process to include new functionality and the creation of semi-synthetic life that stores increased information.

Altered minor-groove hydrogen bonds in DNA block transcription elongation by T7 RNA polymerase

DNA transcription depends upon the highly efficient and selective function of RNA polymerases (RNAPs). Modifications in the template DNA can impact the progression of RNA synthesis, and a number of DNA adducts, as well as abasic sites, arrest or stall transcription. Nonetheless, data are needed to understand why certain modifications to the structure of DNA bases stall RNA polymerases while others are efficiently bypassed. In this study, we evaluate the impact that alterations in dNTP/rNTP base-pair geometry have on transcription. T7 RNA polymerase was used to study transcription over modified purines and pyrimidines with altered H-bonding capacities. The results suggest that introducing wobble base-pairs into the DNA:RNA heteroduplex interferes with transcriptional elongation and stalls RNA polymerase. However, transcriptional stalling is not observed if mismatched base-pairs do not H-bond. Together, these studies show that RNAP is able to discriminate mismatches resulting in wobble base-pairs, and suggest that, in cases of modifications with minor steric impact, DNA:RNA heteroduplex geometry could serve as a controlling factor for initiating transcription-coupled DNA repair.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

The nature of the transition mismatches with Watson-Crick architecture: the G*·T or G·T* DNA base mispair or both? A QM/QTAIM perspective for the biological problem

This study provides the first accurate investigation of the tautomerization of the biologically important guanine*·thymine (G*·T) DNA base mispair with Watson-Crick geometry, involving the enol mutagenic tautomer of the G and the keto tautomer of the T, into the G·T* mispair (∆G = .99 kcal mol(-1), population = 15.8% obtained at the MP2 level of quantum-mechanical theory in the continuum with ε = 4), formed by the keto tautomer of the G and the enol mutagenic tautomer of the T base, using DFT and MP2 methods in vacuum and in the weakly polar medium (ε = 4), characteristic for the hydrophobic interfaces of specific protein-nucleic acid interactions. We were first able to show that the G*·T↔G·T* tautomerization occurs through the asynchronous concerted double proton transfer along two antiparallel O6H···O4 and N1···HN3 H-bonds and is assisted by the third N2H···O2 H-bond, that exists along the entire reaction pathway. The obtained results indicate that the G·T* base mispair is stable from the thermodynamic point of view complex, while it is dynamically unstable structure in vacuum and dynamically stable structure in the continuum with ε = 4 with lifetime of 6.4·10(-12) s, that, on the one side, makes it possible to develop all six low-frequency intermolecular vibrations, but, on the other side, it is by three orders less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. One of the more significant findings to emerge from this study is that the short-lived G·T* base mispair, which electronic interaction energy between the bases (-23.76 kcal mol(-1)) exceeds the analogical value for the G·C Watson-Crick nucleobase pair (-20.38 kcal mol(-1)), "escapes from the hands" of the DNA replication machinery by fast transforming into the G*·T mismatch playing an indirect role of its supplier during the DNA replication. So, exactly the G*·T mismatch was established to play the crucial role in the spontaneous point mutagenesis.

Recognition of Watson-Crick base pairs: constraints and limits due to geometric selection and tautomerism

The natural bases of nucleic acids have a strong preference for one tautomer form, guaranteeing fidelity in their hydrogen bonding potential. However, base pairs observed in recent crystal structures of polymerases and ribosomes are best explained by an alternative base tautomer, leading to the formation of base pairs with Watson-Crick-like geometries. These observations set limits to geometric selection in molecular recognition of complementary Watson-Crick pairs for fidelity in replication and translation processes.

Stochastic adhesion of hydroxylated atomic force microscopy tips to supported lipid bilayers

This work reports results of an atomic force microscopy (AFM) study of adhesion force between hydroxylated AFM tips and supported lipid bilayers (SLBs) of phosphatidylcholine in phosphate buffer saline solution at neutral pH. Silicon nitride AFM probes were hydroxylated by treatment in water vapor plasma and used in force spectroscopy measurements of adhesion force on SLBs with control of contact loading force and residence time. The measurements showed a stochastic behavior of adhesion force that was attributed to stochastic formation of hydrogen bonds between the hydroxyl groups on the AFM tip and oxygen atoms from the phosphate groups of the phosphatidylcholine molecules. Analysis of a large number of force curves revealed a very low probability of hydrogen bond formation, a probability that increased with the increase of contact loading force and residence time. The variance and mean values of adhesion force showed a linear dependence on each other, which indicated that hydrogen bond formation obeyed the Poisson distribution of probability. This allowed for the quantitative determination of the rupture force per hydrogen bond of about 40 pN and showed the absence of other nonspecific interaction forces.

Structural insights into DNA replication without hydrogen bonds

The genetic alphabet is composed of two base pairs, and the development of a third, unnatural base pair would increase the genetic and chemical potential of DNA. d5SICS-dNaM is one of the most efficiently replicated unnatural base pairs identified to date, but its pairing is mediated by only hydrophobic and packing forces, and in free duplex DNA it forms a cross-strand intercalated structure that makes its efficient replication difficult to understand. Recent studies of the KlenTaq DNA polymerase revealed that the insertion of d5SICSTP opposite dNaM proceeds via a mutually induced-fit mechanism, where the presence of the triphosphate induces the polymerase to form the catalytically competent closed structure, which in turn induces the pairing nucleotides of the developing unnatural base pair to adopt a planar Watson-Crick-like structure. To understand the remaining steps of replication, we now report the characterization of the prechemistry complexes corresponding to the insertion of dNaMTP opposite d5SICS, as well as multiple postchemistry complexes in which the already formed unnatural base pair is positioned at the postinsertion site. Unlike with the insertion of d5SICSTP opposite dNaM, addition of dNaMTP does not fully induce the formation of the catalytically competent closed state. The data also reveal that once synthesized and translocated to the postinsertion position, the unnatural nucleobases again intercalate. Two modes of intercalation are observed, depending on the nature of the flanking nucleotides, and are each stabilized by different interactions with the polymerase, and each appear to reduce the affinity with which the next correct triphosphate binds. Thus, continued primer extension is limited by deintercalation and rearrangements with the polymerase active site that are required to populate the catalytically active, triphosphate bound conformation.

The excision of 3' penultimate errors by DNA polymerase I and its role in endonuclease V-mediated DNA repair

Deamination of adenine can occur spontaneously under physiological conditions, and is enhanced by exposure of DNA to ionizing radiation, UV light, nitrous acid, or heat, generating the highly mutagenic lesion of deoxyinosine in DNA. Such DNA lesions tends to generate A:T to G:C transition mutations if unrepaired. In Escherichia coli, deoxyinosine is primarily removed through a repair pathway initiated by endonuclease V (endo V). In this study, we compared the repair of three mutagenic deoxyinosine lesions of A-I, G-I, and T-I using E. coli cell-free extracts as well as reconstituted protein system. We found that 3'-5' exonuclease activity of DNA polymerase I (pol I) was very important for processing all deoxyinosine lesions. To understand the nature of pol I in removing damaged nucleotides, we systemically analyzed its proofreading to 12 possible mismatches 3'-penultimate of a nick, a configuration that represents a repair intermediate generated by endo V. The results showed all mismatches as well as deoxyinosine at the 3' penultimate site were corrected with similar efficiency. This study strongly supports for the idea that the 3'-5' exonuclease activity of E. coli pol I is the primary exonuclease activity for removing 3'-penultimate deoxyinosines derived from endo V nicking reaction.

What Stabilizes the LinPn Inorganic Double Helices?

Recently, inorganic double-helical structures based on simple lithium and phosphorus salts have been demonstrated. We have analyzed the nature of bonding in these double helices of LinPn (n = 7-9). Similar to DNA, noncovalent interactions and cooperativity play an important role in stabilizing the inorganic double helices. Cooperativity imparts an additional stabilization of 4.5 to 10.1 kcal/mol per Li-P pair. In addition, the ionicity of Li-P units further augments the stability of these inorganic double-helix structures in contrast with the canonical base pairs in DNA, where noncovalent interactions determine the duplex stability. Unwinding is shown to be unfavorable, and cleavage of a few edge Li-P bonds leads to spontaneous self-healing into the intact double-helix tract.

The energetic difference between synthesis of correct and incorrect base pairs accounts for highly accurate DNA replication

To better understand the energetics of accurate DNA replication, we directly measured ΔG(o) for the incorporation of a nucleotide into elongating dsDNA in solution (ΔG(o)(incorporation)). Direct measurements of the energetic difference between synthesis of correct and incorrect base pairs found it to be much larger than previously believed (average ΔΔG(o)(incorporation) = 5.2 ± 1.34 kcal mol(-1)). Importantly, these direct measurements indicate that ΔΔG(o)(incorporation) alone can account for the energy required for highly accurate DNA replication. Evolutionarily, these results indicate that the earliest polymerases did not have to evolve sophisticated mechanisms to replicate nucleic acids; they may only have had to take advantage of the inherently more favorable ΔG(o) for polymerization of correct nucleotides. These results also provide a basis for understanding how polymerases replicate DNA (or RNA) with high fidelity.

DNA polymerase delta in DNA replication and genome maintenance

The eukaryotic genome is in a constant state of modification and repair. Faithful transmission of the genomic information from parent to daughter cells depends upon an extensive system of surveillance, signaling, and DNA repair, as well as accurate synthesis of DNA during replication. Often, replicative synthesis occurs over regions of DNA that have not yet been repaired, presenting further challenges to genomic stability. DNA polymerase δ (pol δ) occupies a central role in all of these processes: catalyzing the accurate replication of a majority of the genome, participating in several DNA repair synthetic pathways, and contributing structurally to the accurate bypass of problematic lesions during translesion synthesis. The concerted actions of pol δ on the lagging strand, pol ϵ on the leading strand, associated replicative factors, and the mismatch repair (MMR) proteins results in a mutation rate of less than one misincorporation per genome per replication cycle. This low mutation rate provides a high level of protection against genetic defects during development and may prevent the initiation of malignancies in somatic cells. This review explores the role of pol δ in replication fidelity and genome maintenance.

Thermodynamic and mechanistic insights into translesion DNA synthesis catalyzed by Y-family DNA polymerase across a bulky double-base lesion of an antitumor platinum drug

To determine how the Y-family translesion DNA polymerase η (Polη) processes lesions remains fundamental to understanding the molecular origins of the mutagenic translesion bypass. We utilized model systems employing a DNA double-base lesion derived from 1,2-GG intrastrand cross-links of a new antitumor Pt(II) complex containing a bulky carrier ligand, namely [PtCl(2)(cis-1,4-dach)] (DACH=diaminocyclohexane). The catalytic efficiency of Polη for the insertion of correct dCTP, with respect to the other incorrect nucleotides, opposite the 1,2-GG cross-link was markedly reduced by the DACH carrier ligand. This reduced efficiency of Polη to incorporate the correct dCTP could be due to a more extensive DNA unstacking and deformation of the minor groove induced in the DNA by the cross-link of bulky [PtCl(2)(cis-1,4-dach)]. The major products of the bypass of this double-base lesion produced by [PtCl(2)(cis-1,4-dach)] by Polη resulted from misincorporation of dATP opposite the platinated G residues. The results of the present work support the thesis that this misincorporation could be due to sterical effects of the bulkier 1,4-DACH ligand hindering the formation of the Polη-DNA-incoming nucleotide complex. Calorimetric analysis suggested that thermodynamic factors may contribute to the forces that governed enhanced incorporation of the incorrect dATP by Polη as well.

2'-Deoxyriboguanylurea, the primary breakdown product of 5-aza-2'-deoxyribocytidine, is a mutagen, an epimutagen, an inhibitor of DNA methyltransferases and an inducer of 5-azacytidine-type fragile sites

5-Aza-2′-deoxycytidine (5azaC-dR) has been employed as an inhibitor of DNA methylation, a chemotherapeutic agent, a clastogen, a mutagen, an inducer of fragile sites and a carcinogen. However, its effects are difficult to quantify because it rapidly breaks down in aqueous solution to the stable compound 2′-deoxyriboguanylurea (GuaUre-dR). Here, we used a phosphoramidite that permits the introduction of GuaUre-dR at defined positions in synthetic oligodeoxynucleotides to demonstrate that it is a potent inhibitor of human DNA methyltransferase 1 (hDNMT1) and the bacterial DNA methyltransferase (M.EcoRII) and that it is a mutagen that can form productive base pairs with either Guanine or Cytosine. Pure GuaUre-dR was found to be an effective demethylating agent and was able to induce 5azaC-dR type fragile sites FRA1J and FRA9E in human cells. Moreover, we report that demethylation associated with C:G → G:C transversion and C:G → T:A transition mutations was observed in human cells exposed to pure GuaUre-dR. The data suggest that most of the effects attributed to 5azaC-dR are exhibited by its stable primary breakdown product.

Expanding the horizon of the thymine isostere biochemistry: unique cyclobutane dimers formed by photoreaction between a thymine and a toluene residue in the dinucleotide framework

Substituted toluenyl groups are considered as close isosteres of the thymine residue. They can be recognized by DNA polymerases as if they were thymine. These toluene derivatives are generally inert toward radical additions, including the [2+2] photo-cycloadditions, due to the stable structure of the aromatic ring and are usually used as solvents for radical reactions. Surprisingly, after incorporating toluene into the dinucleotide framework, we found that the UV excited thymine residue readily dimerizes with the toluenyl moiety through a [2+2] photo-addition reaction. Furthermore, the reaction site on the toluenyl moiety is not the C5=C6 bond, as commonly observed in cyclobutane pyrimidine dimers, but the C4=C5 or C3=C4 instead. Such a reaction pattern suggests that in the stacked structure, it is one of these bonds, not the C5=C6, that is close to the thymine C5=C6 bond. A similar structural feature is found in DNA duplex with a thymine replaced by a 2,4-difluorotoluene. Our results argue that although the substituted toluenyl moieties closely mimic the size and shape of the thymine residue, their more hydrophobic nature determines that they stack on DNA bases differently from the natural thymine residue and likely cause local conformational changes in duplex DNA.

KlenTaq polymerase replicates unnatural base pairs by inducing a Watson-Crick geometry

Many candidate unnatural DNA base pairs have been developed, but some of the best-replicated pairs adopt intercalated structures in free DNA that are difficult to reconcile with known mechanisms of polymerase recognition. Here we present crystal structures of KlenTaq DNA polymerase at different stages of replication for one such pair, dNaM-d5SICS, and show that efficient replication results from the polymerase itself, inducing the required natural-like structure.

Amino acid templating mechanisms in selection of nucleotides opposite abasic sites by a family a DNA polymerase

Cleavage of the N-glycosidic bond that connects the nucleobase to the backbone in DNA leads to abasic sites, the most frequent lesion under physiological conditions. Several DNA polymerases preferentially incorporate an A opposite this lesion, a phenomenon termed "A-rule." Accordingly, KlenTaq, the large fragment of Thermus aquaticus DNA polymerase I, incorporates a nucleotide opposite an abasic site with efficiencies of A > G > T > C. Here we provide structural insights into constraints of the active site during nucleotide selection opposite an abasic site. It appears that these confines govern the nucleotide selection mainly by interaction of the incoming nucleotide with Tyr-671. Depending on the nucleobase, the nucleotides are differently positioned opposite Tyr-671 resulting in different alignments of the functional groups that are required for bond formation. The distances between the α-phosphate and the 3'-primer terminus increases in the order A < G < T, which follows the order of incorporation efficiency. Additionally, a binary KlenTaq structure bound to DNA containing an abasic site indicates that binding of the nucleotide triggers a remarkable rearrangement of enzyme and DNA template. The ability to resolve the stacking arrangement might be dependent on the intrinsic properties of the respective nucleotide contributing to nucleotide selection. Furthermore, we studied the incorporation of a non-natural nucleotide opposite an abasic site. The nucleotide was often used in studying stacking effects in DNA polymerization. Here, no interaction with Tyr-761 as found for the natural nucleotides is observed, indicating a different reaction path for this non-natural nucleotide.

Sequence context effect on strand slippage in natural DNA primer-templates

Strand slippage has been found to occur in primer-templates containing a templating thymine, cytosine, and guanine, leading to the formation of misaligned structures with a single-nucleotide bulge. If remained in the active site of low-fidelity polymerases during DNA replication, these misaligned structures can ultimately bring about deletion mutations. In this study, we performed NMR investigations on primer-template models containing a templating adenine. Similar to our previous results on guanine, adenine templates are also less prone to strand slippage than pyrimidine templates. Misalignment occurs only in primer-templates that form a terminal C·G or G·C base pair. Together with our previous findings on thymine, cytosine, and guanine templates, the present study reveals strand slippage can occur in any kind of natural templating bases during DNA replication, providing insights into the origin of mutation hotspots in natural DNA sequences. In addition to the type of incoming base upon misincorporation, the propensity of strand slippage in primer-templates depends also on the type of templating base, its upstream and downstream bases.

Importance of steric effects on the efficiency and fidelity of transcription by T7 RNA polymerase

DNA-dependent RNA polymerases such as T7 RNA polymerase (T7 RNAP) perform the transcription of DNA into mRNA with high efficiency and high fidelity. Although structural studies have provided a detailed account of the molecular basis of transcription, the relative importance of factors like hydrogen bonds and steric effects remains poorly understood. We report herein the first study aimed at systematically probing the importance of steric and electrostatic effects on the efficiency and fidelity of DNA transcription by T7 RNAP. We used synthetic nonpolar analogues of thymine with sizes varying in subangstrom increments to probe the steric requirements of T7 RNAP during the elongation mode of transcription. Enzymatic assays with internal radiolabeling were performed to compare the efficiency of transcription of modified DNA templates with a natural template containing thymine as a reference. Furthermore, we analyzed effects on the fidelity by measuring the composition of RNA transcripts by enzymatic digestion followed by two-dimensional thin layer chromatography separation. Our results demonstrate that hydrogen bonds play an important role in the efficiency of transcription but, interestingly, do not appear to be required for faithful transcription. Steric effects (size and shape variations) are found to be significant both in insertion of a new RNA base and in extension beyond it.

Solution structure and stability of the DNA undecamer duplexes containing oxanine mismatch

Solution structures of DNA duplexes containing oxanine (Oxa, O) opposite a cytosine (O:C duplex) and opposite a thymine (O:T duplex) have been solved by the combined use of ¹H NMR and restrained molecular dynamics calculation. One mismatch pair was introduced into the center of the 11-mer duplex of [d(GTGACO₆CACTG)/d(CAGTGX₁₇GTCAC), X = C or T]. ¹H NMR chemical shifts and nuclear Overhauser enhancement (NOE) intensities indicate that both the duplexes adopt an overall right-handed B-type conformation. Exchangeable resonances of C₁₇ 4-amino proton of the O:C duplex and of T₁₇ imino proton of O:T duplex showed unusual chemical shifts, and disappeared with temperature increasing up to 30°C, although the melting temperatures were >50°C. The O:C mismatch takes a wobble geometry with positive shear parameter where the Oxa ring shifted toward the major groove and the paired C₁₇ toward the minor groove, while, in the O:T mismatch pair with the negative shear, the Oxa ring slightly shifted toward the minor groove and the paired T₁₇ toward the major groove. The Oxa mismatch pairs can be wobbled largely because of no hydrogen bond to the O1 position of the Oxa base, and may occupy positions in the strands that optimize the stacking with adjacent bases.

DNA analysis based on the local structural disruption to the duplexes carrying a luminous lanthanide complex

The two DNA conjugates (split probes) carrying a metal chelator form an integrated luminous lanthanide (Ln(3+): Tb(3+) or Eu(3+)) complex on the complementary template DNA (target). The luminous property of this Ln(3+) complex has been used for DNA assay. The intensity of the luminescence was affected by the local structural disruption caused by one-base mispairing around the complex. Among the mispairings systematically introduced around the Ln(3+) center, vicinal mispairings to the center decreased the emission intensity more. This would be a novel nucleobase-discriminating principle, in which the split probes bind the target tightly, yet still retain sequence selectivity.