Free Energy Landscape between Watson-Crick to Hoogsteen Base Pairing Transitions.

Genetic information is encoded in the DNA double helix, which, in its physiological milieu, is characterized by the iconical Watson-Crick nucleo-base pairing. Recent NMR relaxation experiments revealed the transient presence of an alternative, Hoogsteen (HG) base pairing pattern in naked DNA duplexes, and estimated its relative stability and lifetime. In contrast with DNA, such structures were not observed in RNA duplexes. Understanding HG base pairing is important because the underlying “breathing” motion between the two conformations can significantly modulate protein binding. However, a detailed mechanistic insight into the transition pathways and kinetics is still missing. We performed enhanced sampling simulation (with combined metadynamics and adaptive force-bias method) and Markov state modeling to obtain accurate free energy, kinetics, and the intermediates in the transition pathway between Watson-Crick and HG base pairs for both naked B-DNA and A-RNA duplexes. The Markov state model constructed from our unbiased MD simulation data revealed previously unknown complex extrahelical intermediates in the seemingly simple process of base flipping in B-DNA. Extending our calculation to A-RNA, for which HG base pairing is not observed experimentally, resulted in relatively unstable, single-hydrogen-bonded, distorted Hoogsteen-like bases. Unlike B-DNA, the transition pathway primarily involved base paired and intrahelical intermediates with transition timescales much longer than that of B-DNA. The seemingly obvious flip-over reaction coordinate (i.e., the glycosidic torsion angle) is unable to resolve the intermediates. Instead, a multidimensional picture involving backbone dihedral angles and distance between hydrogen bond donor and acceptor atoms is required to gain insight into the molecular mechanism.

Direct evidence for (G)O6···H2-N4(C)+ hydrogen bonding in transient G(syn)-C+ and G(syn)-m5C+ Hoogsteen base pairs in duplex DNA from cytosine amino nitrogen off-resonance R1ρ relaxation dispersion measurements.

The vibrational fingerprints of hydrogen-bonding associated with the adenine-thymine (A-T) Watson-Crick (WC) base pair have been identified in an infrared study of the A-T mimics 4-aminopyrimidine-1-methylthymine (4APM-1MT) and 4-aminopyrimidine-6-methyl-4-pyrimidinone (4APM-M4PMN) in the gas-phase. The IR vibrational spectra were measured via a double resonance scheme utilizing femtosecond multiphoton ionization. The changes in the molecular structure, anharmonic vibrational parameters, and the assignment of the observed vibrational spectra in the NH/CH stretch region were investigated by carrying out high-level theoretical calculations of the anharmonic spectra. The experimental observations and theoretical calculations indicate that the hydrogen bonds associated with WC base-pairing are relatively stronger than those associated with reverse WC (rWC) base pairing. This is manifested in a more pronounced red-shift of the H-bonded vibrational modes associated with the WC as compared with the rWC base-pairing. An analysis of the factors contributing to the anharmonicity of the vibrational modes associated with H-bonding reveals that the magnitude of the off-diagonal anharmonic coupling of the H-bonded -NH2 stretch and the -NH2 bend is much smaller in WC base-pairing than in the corresponding rWC base-pairing. The chemical and biological implications of these results, especially in the context of using vibrational spectroscopy as a tool for identifying the signatures of nucleotide base vibrations is addressed.

DNA predominantly contains Watson–Crick (WC) base pairs, but a non-negligible fraction of base pairs are in the Hoogsteen (HG) hydrogen bonding motif at any time. In HG, the purine is rotated ∼180° relative to the WC motif. The transitions between WC and HG may play a role in recognition and replication, but are difficult to investigate experimentally because they occur quickly, but only rarely. To gain insight into the mechanisms for this process, we performed transition path sampling simulations on a model nucleotide sequence in which an AT pair changes from WC to HG. This transition can occur in two ways, both starting with loss of hydrogen bonds in the base pair, followed by rotation around the glycosidic bond. In one route the adenine base converts from WC to HG geometry while remaining entirely within the double helix. The other route involves the adenine leaving the confines of the double helix and interacting with water. Our results indicate that this outside route is more probable. We used transition interface sampling to compute rate constants and relative free energies for the transitions between WC and HG. Our results agree with experiments, and provide highly detailed insights into the mechanisms of this important process.

Conformational stability and order of Hoogsteen base pair induced by protein binding

NMR relaxation dispersion studies indicate that in canonical duplex DNA, Watson-Crick base pairs (bps) exist in dynamic equilibrium with short-lived low abundance excited state Hoogsteen bps. N1-methylated adenine (m1A) and guanine (m1G) are naturally occurring forms of damage that stabilize Hoogsteen bps in duplex DNA. NMR dynamic ensembles of DNA duplexes with m1A-T Hoogsteen bps reveal significant changes in sugar pucker and backbone angles in and around the Hoogsteen bp, as well as kinking of the duplex towards the major groove. Whether these structural changes also occur upon forming excited state Hoogsteen bps in unmodified duplexes remains to be established because prior relaxation dispersion probes provided limited information regarding the sugar-backbone conformation. Here, we demonstrate measurements of C3' and C4' spin relaxation in the rotating frame (R1ρ) in uniformly 13C/15N labeled DNA as sensitive probes of the sugar-backbone conformation in DNA excited states. The chemical shifts, combined with structure-based predictions using an automated fragmentation quantum mechanics/molecular mechanics method, show that the dynamic ensemble of DNA duplexes containing m1A-T Hoogsteen bps accurately model the excited state Hoogsteen conformation in two different sequence contexts. Formation of excited state A-T Hoogsteen bps is accompanied by changes in sugar-backbone conformation that allow the flipped syn adenine to form hydrogen-bonds with its partner thymine and this in turn results in overall kinking of the DNA toward the major groove. Results support the assignment of Hoogsteen bps as the excited state observed in canonical duplex DNA, provide an atomic view of DNA dynamics linked to formation of Hoogsteen bps, and lay the groundwork for a potentially general strategy for solving structures of nucleic acid excited states.

UV-light-induced cyclobutane pyrimidine dimers (CPDs) present a severe block to synthesis by replicative DNA polymerases (Pols), whereas Poleta promotes proficient and error-free replication through CPDs. Although the archael Dpo4, which, like Poleta, belongs to the Y family of DNA Pols, can also replicate through a CPD, it is much less efficient than Poleta. The x-ray crystal structure of Dpo4 complexed with either the 3'-thymine (T) or the 5' T of a cis-syn TT dimer has indicated that, whereas the 3' T of the dimer forms a Watson-Crick base pair with the incoming dideoxy ATP, the 5' T forms a Hoogsteen base pair with the dideoxy ATP in syn conformation. Based upon these observations, a similar mechanism involving Hoogsteen base pairing of the 5' T of the dimer with the incoming A has been proposed for Poleta. Here we examine the mechanisms of CPD bypass by Dpo4 and Poleta using nucleotide analogs that specifically disrupt the Hoogsteen or Watson-Crick base pairing. Our results show that both Dpo4 and Poleta incorporate dATP opposite the 5' T of the CPD via Watson-Crick base pairing and not by Hoogsteen base pairing. Furthermore, opposite the 3' T of the dimer, the two Pols differ strikingly in the mechanisms of dATP incorporation, with Dpo4 incorporating opposite an abasic-like intermediate and Poleta using the normal Watson-Crick base pairing. These observations have important implications for the mechanisms used for the inefficient vs. efficient bypass of CPDs by DNA Pols.

The insertion of a DNA base moiety at the end of a DNA duplex to form a Watson-Crick or wobble pair during DNA annealing or replication is a step of fundamental biological importance. Therefore, we investigated the energetics of a formation of the terminal G x C, G x T, and G x A base pairs in DNA containing a 5'-dangling G adjacent to the base insertion point using differential scanning calorimetry and computer simulations. The energies calculated along classical molecular dynamics trajectories in aqueous solution were analyzed in the framework of linear-response approximation (LRA) to obtain relative free energies for the base insertion and their electrostatic, van der Waals, and preorganization components. Using the generic set of LRA parameters, the calculated free energies disfavored the mispair formation by 2.5 (G x C --> G x T) and 1.7 (G x C --> G x A) kcal/mol, in reasonable agreement with the experimental free energy differences of 1.8 and 1.4 kcal/mol, respectively. The calculated preorganization components of these free energies of 0.6 (G x C --> G x T) and -0.1 (G x C --> G x A) kcal/mol show that electrostatic preorganization, which is an important source of DNA replication fidelity, plays a lesser role in the mispair destabilization in the absence of DNA polymerase.

The structures and stabilities of cyclic DNA octamers of different sequences have been studied by NMR and CD spectroscopy and by restrained molecular dynamics. At low oligonucleotide concentrations, some of these molecules form stable monomeric structures consisting of a short stem of two base pairs connected by two mini-loops of two residues. To our knowledge, these dumbbell-like structures are the smallest observed to date. The relative stabilities of these cyclic dumbbells have been established by studying their melting transitions. Dumbbells made up purely of GC stems are more stable than those consisting purely of AT base pairs. The order of the base pairs closing the loops also has an important effect on the stabilities of these structures. The NMR data indicate that there are significant differences between the solution structures of dumbbells with G-C base pairs in the stem compared to those with A-T base pairs. In the case of dumbbells with G-C base pairs, the residues in the stem form a short segment of a BDNA helix stabilized by two Watson-Crick base pairs. In contrast, in the case of d<pCATTCATT>, the stem is formed by two A-T base pairs with the glycosidic angles of the adenine bases in a syn conformation, most probably forming Hoogsteen base pairs. Although the conformations of the loop residues are not very well defined, the thymine residues at the first position of the loop are observed to fold back into the minor groove of the stem.

Measuring the strength of the hydrogen bonds between DNA base pairs is of vital importance for understanding how our genetic code is physically accessed and recognized in cells, particularly during replication and transcription. Therefore, it is important to develop probes for these key hydrogen bonds (H-bonds) that dictate events critical to cellular function, such as the localized melting of DNA. The vibrations of carbonyl bonds are well-known probes of their H-bonding environment, and their signals can be observed with infrared (IR) spectroscopy. Yet, pinpointing a single bond of interest in the complex IR spectrum of DNA is challenging due to the large number of carbonyl signals that overlap with each other. Here, we develop a method using isotope editing and infrared (IR) spectroscopy to isolate IR signals from the thymine (T) C2=O carbonyl. We use solvatochromatic studies to show that the TC2=O signal’s position in the IR spectrum is sensitive to the H-bonding capacity of the solvent. Our results indicate that C2=O of a single T base within DNA duplexes experiences weak H-bonding interactions. This finding is consistent with the existence of a third, noncanonical CH···O H-bond between adenine and thymine in both Watson–Crick and Hoogsteen base pairs in DNA.

The synthetic deoxydodecamer d(C-G-C-G-A-A-T-T-A-G-C-G) was analyzed by x-ray diffraction methods, and the structure was refined to a residual error of R = 0.17 at 2.5-A resolution (2 sigma data) with 83 water molecules located. The sequence crystallizes as a full turn of a B-DNA helix and contains 2 purine X purine (G.A) base pairs and 10 Watson-Crick base pairs. The analysis shows conclusively that adenine is in the syn orientation with respect to the sugar moiety whereas guanine adopts the usual trans orientation. Nitrogen atoms of both bases are involved in hydrogen bonding with the N-1 of guanine 2.84 A from the N-7 of adenine and the N-6 of adenine within 2.74 A of the O-6 of guanine. The C-1'...C-1' separation is 10.7 A close to that for standard Watson-Crick base pairs. The incorporation of the purine.purine base pairs at two steps in the dodecamer causes little perturbation of either the local or the global conformation of the double helix. Comparison of the structural features with those of the G.T wobble pair and the standard G.C pair suggests a rationale for the differential enzymatic repair of the two types of base-pair mismatches.

We found that Hoogsteen base pairs were stabilized by molecular crowding and a histone H3-mimicking peptide, which was not observed for Watson-Crick base pairs. Our findings demonstrate that the type of DNA base pair is critical for the interaction between DNA and histones.

The relation between DNA polymerase fidelity and base pairing stability is investigated by using DNA primer-template duplexes that contain a common 9-base template sequence but have either correct (A.T) or incorrect (G.T, C.T, T.T) base pairs at the primer 3' terminus. Thermal melting and enzyme kinetic measurements are compared for each kind of terminus. Analysis of melting temperatures finds that differences between the free energy changes upon dissociation (delta delta Go) are only 0.2, 0.3, and 0.4 kcal.mol-1 (1 cal = 4.18 J) for terminal A.T compared to G.T, C.T, and T.T mispairs, respectively, at 37 degrees C. We show that enthalpy changes are directly correlated with entropy changes for normal and abnormal base pairs in DNA in aqueous solution and that delta delta Go values are small because of near cancellation of corresponding enthalpy and entropy components. The kinetics of elongating primer termini are measured with purified Drosophila DNA polymerase alpha. The matched A.T terminus is found to be extended approximately 200 times faster than a G.T mismatch and 1400 and 2500 times faster than C.T and T.T mismatches, respectively. Enzymatic discrimination against elongating mismatched termini is based mainly on Km rather than Vmax differences. From Km at 37 degrees C, we find delta delta Go values of 2.6-3.7 kcal.mol-1, about an order of magnitude greater than indicated by melting data. A similar measurement of nucleotide insertion kinetics has previously found rates of forming A.T base pairs to be 500 times greater than G.T mispairs and 20,000 times greater than C.T and T.T mispairs. Here also, Km differences are mainly responsible for discrimination and indicate even larger delta delta Go values (4.3-4.9 kcal.mol-1). Thus, free energy differences between correct and incorrect base pairs in the active site cleft of polymerase appear to be greater than 10 times as large as in aqueous medium. We explore the idea that a binding cleft that snugly fits correct base pairs and excludes water at the active site may amplify base-pair free energy differences by reducing entropy differences and increasing enthalpy differences sufficiently to account for nucleotide insertion and extension fidelity.

There is rapid development of drug resistant microorganisms in the environment due to drug abuse and poor environmental sanitation among other factors. Strains of Pseudomonas putida are inhabitants of soil and freshwater that exhibit amazing ability to metabolize a wide range of biogenic and xenobiotic compounds. Some strains are however known to possess plasmid bearing genes that encodes for antibiotic resistant factor which can be transferred to other microorganisms. To inhibit the growth of prokaryotes, antibiotics are prepared to target the ribosomal RNA. In order to make possible discovery of improved antibiotics that can target the ribosomal RNA of P. putida, the secondary and tertiary structures of P. putida strain TB3 16s rRNA is being described. P. putida strain TB3 was isolated following standard procedures from kolanut waste (kolanut husk and testa) obtained from Owena, Ondo State, one of the regions in Nigeria where kolanut is largely cultivated for commercial purposes. The isolates were identified using molecular techniques and Sanger sequencing to identify labeled chain of terminating nucleotides that are incorporated by DNA polymerase during replication of the templates. The secondary structure of the rRNA sequence of the isolate was determined by the interaction between its bases, including hydrogen bonding, base stacking and loops defined as R = r1, r2 . . . ; rN with the indices 1< i < N numbering the nucleotides in the 5l to 3l direction and base pairs were denoted by ri:rj. Nearest- neighbor model was used for the secondary structure prediction and the minimum free energy of the structure was estimated by summing individual energy contributions from base pair stacking, hairpins, bulges, internal loops and multi-branch loops using the function rnafold on Matlab software. An assortment of idiosyncratic chemical and structural properties was observed on the base pairs: they have high affinity for metal ions, they present conformational flexibility to different environments and they are almost thermodynamically unchanging as Watson-Crick base pairs. With the development of new drug discovery technologies, targeting RNA for better antibiotic development is emerging as a new frontier in drug discovery

We present an efficient method based on an extension of metadynamics for exploring complex free energy landscapes (FELs). The method employs two-step metadynamics simulations. In the first step, rapid metadynamics simulations using broad and tall Gaussians are performed to identify a free energy pathway (FEP) connecting the two states of interest. The FEP is then divided into a series of independent subphase spaces that comprise selected discrete images of the system. Using appropriate collective variables (CVs) chosen according to the FEP, the accurate FEL of each subphase space is separately calculated in subsequent divide-and-conquer metadynamics simulations with narrow and low Gaussians. Finally, all FELs calculated in each subphase space are merged to obtain the full FEL. We show that the method greatly improves the performance of the metadynamics approach. In particular, we are able to efficiently model chemical systems with complex FELs, such as chemical reactions at the air/water interface. We demonstrate the performance of this method on two model reactions: the hydrolysis of formaldehyde in the gas phase and at the air/water interface.