Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.     

The crystal structure at 1.5 angstroms resolution of an RNA octamer duplex containing tandem G.U basepairs

The crystal structure of the RNA octamer, 5'-GGCGUGCC-3' has been determined from x-ray diffraction data to 1.5 angstroms resolution. In the crystal, this oligonucleotide forms five self-complementary double-helices in the asymmetric unit. Tandem 5'GU/3'UG basepairs comprise an internal loop in the middle of each duplex. The NMR structure of this octameric RNA sequence is also known, allowing comparison of the variation among the five crystallographic duplexes and the solution structure. The G.U pairs in the five duplexes of the crystal form two direct hydrogen bonds and are stabilized by water molecules that bridge between the base of guanine (N2) and the sugar (O2') of uracil. This contrasts with the NMR structure in which only one direct hydrogen bond is observed for the G.U pairs. The reduced stability of the r(CGUG)2 motif relative to the r(GGUC)2 motif may be explained by the lack of stacking of the uracil bases between the Watson-Crick and G.U pairs as observed in the crystal structure.     

Conformational analysis of a covalently cross-linked Watson-Crick base pair model

Low-temperature NMR experiments and molecular modeling have been used to characterize the conformational behavior of a covalently cross-linked DNA base pair model. The data suggest that Watson-Crick or reverse Watson-Crick hydrogen bonding geometries have similar energies and can interconvert at low temperatures. This low-temperature process involves rotation about the crosslink CH(2)C(5') (psi) carbon-carbon bond, which is energetically preferred over the alternate CH(2)N(3) (phi) carbon-nitrogen bond rotation.     

Sheared Aanti.Aanti base pairs in a destabilizing 2 x 2 internal loop: the NMR structure of 5'(rGGCAAGCCU)2

The 5'(rGGCAAGCCU)(2) duplex contains tandem A.A pairs. The three-dimensional structure of the 5'(rGGCAAGCCU)(2) duplex was modeled by molecular dynamics and energy minimization with NMR-derived distance and dihedral angle restraints. Although the 5'(rCAAG)(2) loop is thermodynamically destabilizing by 1.1 kcal/mol, the tandem A.A pairs adopt a predominant conformation: a sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment similar to that observed in the crystal structure of the P4-P6 domain of the Tetrahymena thermophila intron [Cate, J. H., Gooding, A. R., Podell, E., Zhou, K., Golden, B. L., Kundrot, C. E., Cech, T. R., and Doudna, J. A. (1996) Science 273, 1678-1685]. The NMR-derived structure of the 5'(rGGCAAGCCU)(2) duplex exhibits cross-strand hydrogen bonds from N3 of A4 to an amino hydrogen of A5 and from the 2' oxygen of the A4 sugar to the other amino hydrogen of A5. An intrastrand hydrogen bond is formed from the 2' OH hydrogen of A4 to O5' of A5. The cross-strand A5 bases are stacked. The Watson-Crick G-C regions are essentially A-form. The sheared anti-anti (A.A trans Hoogsteen/Sugar-edge) alignment provides potential contact sites for tertiary interactions and, therefore, is a possible target site for therapeutics. Thus, thermodynamically destabilizing internal loops can be preorganized for tertiary interactions or ligand binding.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

NMR solution structure of a DNA dodecamer containing a transplatin interstrand GN7-CN3 cross-link.

The DNA duplex d(CTCTCG*AGTCTC).d(GAGAC-TC*GAGAG) containing a single trans- diammine-dichloroplatinum(II) interstrand cross-link (where G* and C* represent the platinated bases) has been studied by two-dimensional NMR. All the exchangeable and non-exchangeable proton resonance lines were assigned (except H5'/H5") and the NOE intensities were transformed into distances via the RELAZ program. By combining the NOESY and COSY data (330 constraints) and NMR-constrained molecular mechanics using JUMNA, a solution structure of the cross-linked duplex has been determined. The duplex is distorted over two base pairs on each side of the interstrand cross-link and exhibits a slight bending of its axis ( approximately 20 degrees ) towards the minor groove. The platinated guanine G* adopts a syn conformation. The rotation results in a Hoogsteen-type pairing between the complementary G(6)* and C(19)* residues which is mediated by the platinum moiety and is stabilized by a hydrogen bond between O6(G(6)*) and N4H(C(19)*). The rise between the cross-linked residues and the adjacent residues is increased owing to the interaction between these adjacent residues and the ammine groups of the platinum moiety. These results are discussed in relation to the slow rate of closure of the monofunctional adducts into interstrand cross-links.     

The crystal structure of the octamer [r(guauaca)dC]2 with six Watson-Crick base-pairs and two 3' overhang residues

The crystal structure of an alternating RNA octamer, r(guauaca)dC (RNA bases are in lower case while the only DNA base is in upper case), with two 3' overhang residues one of them a terminal deoxycytosine and the other a ribose adenine, has been determined at 2.2 A resolution. The refined structure has an Rwork 18.6% and Rfree 26.8%. There are two independent duplexes (molecules I and II) in the asymmetric unit cell, a = 24.95, b = 45.25 and c = 73.67 A, with space group P2(1)2(1)2(1). Instead of forming a blunt end duplex with two a+.c mispairs and six Watson-Crick base-pairs, the strands in the duplex slide towards the 3' direction forming a two-base overhang (radC) and a six Watson-Crick base-paired duplex. The duplexes are bent (molecule I, 20 degrees; molecule II, 25 degrees) and stack head-to-head to form a right-handed superhelix. The overhang residues are looped out and the penultimate adenines of the two residues at the top end (A15) are anti and at the bottom (A7) end are syn. The syn adenine bases form minor groove A*(G.C) base triples with C8-H...N2 hydrogen bonds. The anti adenine in molecule II also forms a triple and a different C2-H...N3 hydrogen bond, while the other anti adenine in molecule I does not, it stacks on the looped out overhang base dC. The 3' terminal deoxycytosines form two stacked hemiprotonated trans d(C.C)+ base-pairs and the pseudo dyad related molecules form four consecutive deoxyribose and ribose zipper hydrogen bonds in the minor groove.     

Theoretical studies of cis-Pt(II)-diammine binding to duplex DNA

The binding of cis-Pt(II) diammine (cis-DP) to double-stranded DNA was studied with several kinked conformations that can accommodate the formation of a square planar complex. Molecular mechanics (MM) calculations were performed to optimize the molecular fit. These results were combined with quantum mechanical (QM) calculations to ascertain the relative energetics of ligand binding through water vs direct binding of the phosphate to the ammine and platinum, and to guide the selection of DNA conformations to model complex formation. Based on QM and MM calculations, models are proposed that may be characterized by several general features. A structure involving hydrogen bonding between each ammine and distinct adjacent phosphate groups, referred to as closed conformation (CC), has already been reported. This is also found in the crystal structure of small dimers. We report alternative conformations that may be important in platination of duplex DNA. They are characterized by an intermediate conformation (IC), involving hydrogen bonding between one ammine and phosphate group, and an open conformation (OC), without ammine phosphate hydrogen bonding. The IC and OC can be stabilized by water bridges in the space between the ammine and the phosphate groups. Sugar puckers alternate from the type C(2')-endo or C(1')-exo (S), to the type C(3')-endo or C(2')-exo (N), with intermediate types near O(1')-endo (O). In general, the sugar puckers alternate from S to N to S through the platinated region (3'-TpG*pG*p-5'), with the complexed strand exhibiting, (3')-S*-N*-S-(5') alternation, while the complementary strand shows either (3')-S*-N*-S-(5') or (3')-S*-N*-O-(5') alternation. In both the OC and IC, a hydrogen bond is found between the ammine and O4(T) on thymine (T) at the (3') end, adjacent to the complex site. There is a continuous range of backbone conformations through the platinated region which relate the OC to the IC. The models presented suggest that the dynamics of the binding of the cis-Pt(II)-diammines to adjacent N7(G) in double-stranded DNA may encompass several conformational possibilities, and that water bridges may play a roll in supporting open and intermediate conformations. Proton-proton distances are reported to assist in the experimental determination of conformations.     

Effect of flanking residues on CA and AA dinucleotides: some rationale

Effect of flanking base pairs on CA and AA dinluceotide step-geometry has been studied using molecular dynamics method. Sixteen dodecameric sequences are constructed for each doublet with all possible bases at their 5' and 3' positions along with their complementary sequences. Structural parameter, formation of Extra Watson Crick bifurcated hydrogen bonds along or across the strands and effect of sodium ions are studied for these sequences. It is found that geometry of CA doublet step is perturbed by the neighboring base pairs, which might be due to inherent flexibility of the step. Flexible character of CA step is reflected in its low bifurcated hydrogen bond formation capability and lower preference of sodium ions to enter in minor or major grooves. AA step on the other hand is quite rigid according to different structural parameters and respond much less to environmental changes due to formation of strong Extra Watson-Crick hydrogen bonds.     

Molecular dynamic simulations of cisplatin- and oxaliplatin-d(GG) intrastand cross-links reveal differences in their conformational dynamics

Mismatch repair proteins, DNA damage-recognition proteins and translesion DNA polymerases discriminate between Pt-GG adducts containing cis-diammine ligands (formed by cisplatin (CP) and carboplatin) and trans-RR-diaminocyclohexane ligands (formed by oxaliplatin (OX)) and this discrimination is thought to be important in determining differences in the efficacy, toxicity and mutagenicity of these platinum anticancer agents. We have postulated that these proteins recognize differences in conformation and/or conformational dynamics of the DNA containing the adducts. We have previously determined the NMR solution structure of OX-DNA, CP-DNA and undamaged duplex DNA in the 5'-d(CCTCAGGCCTCC)-3' sequence context and have shown the existence of several conformational differences in the vicinity of the Pt-GG adduct. Here we have used molecular dynamics simulations to explore differences in the conformational dynamics between OX-DNA, CP-DNA and undamaged DNA in the same sequence context. Twenty-five 10 ns unrestrained fully solvated molecular dynamics simulations were performed starting from two different DNA conformations using AMBER v8.0. All 25 simulations reached equilibrium within 4 ns, were independent of the starting structure and were in close agreement with previous crystal and NMR structures. Our data show that the cis-diammine (CP) ligand preferentially forms hydrogen bonds on the 5' side of the Pt-GG adduct, while the trans-RR-diaminocyclohexane (OX) ligand preferentially forms hydrogen bonds on the 3' side of the adduct. In addition, our data show that these differences in hydrogen bond formation are strongly correlated with differences in conformational dynamics, specifically the fraction of time spent in different DNA conformations in the vicinity of the adduct, for CP- and OX-DNA adducts. We postulate that differential recognition of CP- and OX-GG adducts by mismatch repair proteins, DNA damage-recognition proteins and DNA polymerases may be due, in part, to differences in the fraction of time that the adducts spend in a conformation favorable for protein binding.     

2D NMR study of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) cross-linked by the antitumor-active dirhodium(II,II) unit at the cytosine-adenine step

The 2D NMR analysis in solution of the DNA duplex d(CTCTC*A*ACTTCC).d(GGAAGTTGAGAG) binding to the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+ showed that an unprecedented intrastrand adduct, dsII, is formed with the dirhodium unit cross-linking in the major groove residues C5 and A6 (indicated with asterisks), also corroborated by enzyme digestion studies. Formation of the dirhodium complex dsII destabilizes significantly the duplex as indicated by the substantial decrease in its melting temperature (DeltaTm = -22.9 degrees C). The reduced thermal stability of dsII is attributed to the decreased stacking of the bases and the complete disruption and/or weakening of the hydrogen bonds within the base pairs in the immediate vicinity of the metalation site (C5.G20 and A6.T19), but the effects due to the metal binding are more severe for the base pairs in the 5' direction to the lesion site. The NMR spectroscopic data indicate that Watson-Crick hydrogen bonding is completely disrupted for the C5.G20 site and considerably weakened for A6.T19. In dsII, the bases C5 and A6 bind to eq positions of the dirhodium unit cis-[Rh2(mu-O2CCH3)2(eta1-O2CCH3)]+, which retains one monodentate and two bridging acetate groups, presumably due to steric reasons. Binding of A6 takes place via N7, whereas binding of the C5 base takes place via the exocyclic N4 site, resulting in the anti-cytosine rotamer with respect to site N3 in its metal-stabilized rare iminooxo form.     

Explicit solvent molecular dynamics simulation of duplex formed by the modified oligonucleotide with alternating phosphate/phosphonate internucleoside linkages and its natural counterpart

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT(11).dA(11) duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5'-O-PO(2)-CH(2)-O-3' (5'PC3') or 5'-O-CH(2)-PO(2)-O-3' (5'CP3') type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5'PC3' and ggg, tgg for the 5'CP3', respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.     

Nuclear magnetic resonance spectroscopy and molecular modeling reveal that different hydrogen bonding patterns are possible for G.U pairs: one hydrogen bond for each G.U pair in r(GGCGUGCC)(2) and two for each G.U pair in r(GAGUGCUC)(2)

G.U pairs occur frequently and have many important biological functions. The stability of symmetric tandem G.U motifs depends both on the adjacent Watson-Crick base pairs, e.g., 5'G > 5'C, and the sequence of the G.U pairs, i.e., 5'-UG-3' > 5'-GU-3', where an underline represents a nucleotide in a G.U pair [Wu, M., McDowell, J. A., and Turner, D. H. (1995) Biochemistry 34, 3204-3211]. In particular, at 37 degrees C, the motif 5'-CGUG-3' is less stable by approximately 3 kcal/mol compared with other symmetric tandem G.U motifs with G-C as adjacent pairs: 5'-GGUC-3', 5'-GUGC-3', and 5'-CUGG-3'. The solution structures of r(GAGUGCUC)(2) and r(GGCGUGCC)(2) duplexes have been determined by NMR and restrained simulated annealing. The global geometry of both duplexes is close to A-form, with some distortions localized in the tandem G.U pair region. The striking discovery is that in r(GGCGUGCC)(2) each G.U pair apparently has only one hydrogen bond instead of the two expected for a canonical wobble pair. In the one-hydrogen-bond model, the distance between GO6 and UH3 is too far to form a hydrogen bond. In addition, the temperature dependence of the imino proton resonances is also consistent with the different number of hydrogen bonds in the G.U pair. To test the NMR models, U or G in various G.U pairs were individually replaced by N3-methyluridine or isoguanosine, respectively, thus eliminating the possibility of hydrogen bonding between GO6 and UH3. The results of thermal melting studies on duplexes with these substitutions support the NMR models.     

Wobble dG X dT pairing in right-handed DNA: solution conformation of the d(C-G-T-G-A-A-T-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

We report below on features of the three-dimensional structure of the d(C-G-T-G-A-A-T-T-C-G-C-G) self-complementary duplex (designated 12-mer GT) containing symmetrical G X T mismatches in the interior of the helix. The majority of the base and sugar protons in the 12-mer GT duplex were assigned by two-dimensional nuclear Overhauser effect (NOESY) spectra in H2O and D2O solution. A set of 92 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds for one symmetrical half of the 12-mer GT duplex were estimated from NOESY data sets recorded as a function of mixing time. These experimental distances combined with nucleotide bond length parameters were embedded into Cartesian space; several trial structures were refined to minimize bond geometry and van der Waals and chirality error. Confidence in this approach is based on the similarity of the refined structures for the solution conformation of the 12-mer GT duplex. The G and T bases pair through two imino-carbonyl hydrogen bonds, and stacking is maintained between the G X T wobble pair and adjacent Watson-Crick G X C pairs. The experimental distance information is restricted to base and sugar protons, and hence structural features such as base pair overlap, glycosidic torsion angles, and sugar pucker are well-defined by this combination of NMR and distance geometry methods. By contrast, we are unable to define the torsion angles about the bonds C3'-O3'-P-O5'-C5'-C4' in the backbone of the nucleic acid.     

Molecular mechanics calculations on a triple stranded DNA involving C+.G-T and T.A+-C mismatched base triples

We have carried out molecular modeling of a triple stranded pyrimidine(Y). purine(R): pyrimidine(Y) (where ':' refers to Watson-Crick and '.' to Hoogsteen bonding) DNA, formed by a homopurine (d-TGAGGAAAGAAGGT) and homo-pyrimidine (d-CTCCTTTCTTCC). Molecular mechanics calculations using NMR constraints have provided a detailed three dimensional structure of the triplex. The entire stretches of purine and the pyrimidine nucleotides have a conformation close to B-DNA. The three strands are held by the canonical C+.G:C and T.A:T hydrogen bonds. The structure also contains two mismatch C+.G-T and T.A+-C base triples which have been characterized for the first time. In the A+-C base-pair of the T.A+-C triple, both hydrogen donors are situated on the purine (A+(1N) and A+(6N)). We observe a unique hydrogen bonding interaction scheme in case of C+.G-T where one acceptor, G(60), is bonded to three donors (C+(3NH), C+(4NH2) and T(3NH)). Though the C+.G-T base triple is less stable than C+.G:C, it is significantly more stable than T.A:T. On the other hand, T.A+-C is as stable as the T.A:T base triad.     

dC-dG alternating oligonucleotides: thermodynamic and kinetic aspects of the B-Z transformation

The alternating cytosine-guanine oligodeoxyribonucleotides (dCdG)n, (dGdC)n, (dCdG)ndC (n = 3,4), (dGdC)7 and dG(dCdG)3 have been studied by UV and CD spectroscopy at different temperatures and NaCl concentrations. The analysis of the melting data, assuming an all-or-none model, reveals that in the B-conformation the 5'G/C3' stacking interactions are enthalpically favoured with respect to the 5'C/G3' one. The CD investigation of the B-Z equilibrium shows that the Z-conformation is enthalpically stabilized, while the B-conformation is entropically favoured, in the range of NaCl concentration considered (1 to 5 M). The kinetic data for the B-Z transformation, obtained with a salt-jump technique for the hexamer (dCdG)3, support a mechanism by which the Watson-Crick hydrogen bonds are broken before the bases flip over separately and eventually stack, reforming the H-bonds, in the new helix.     

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen bonded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel triplex, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41--O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees , helical twist of the average structure from this simulation had a value of 36 degrees , while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel triplex, it was energetically comparable to the parallel triplex. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel triplex by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically. Together these results indicate that the parallel C.G*G triplex with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.     

Herpes simplex virus 1 primase employs watson-crick hydrogen bonding to identify cognate nucleoside triphosphates

We utilized NTP analogues containing modified bases to probe the mechanism of NTP selection by the primase activity of the herpes simplex virus 1 helicase-primase complex. Primase readily bound NTP analogues of varying base shape, hydrophobicity, and hydrogen-bonding capacity. Remarkably, primase strongly discriminated against incorporating virtually all of the analogues, even though this enzyme misincorporates natural NTPs at frequencies as high as 1 in 7. This included analogues with bases much more hydrophobic than a natural base (e.g., 4- and 7-trifluoromethylbenzimidazole), a base of similar hydrophobicity as a natural base but with the Watson-Crick hydrogen-bonding groups in unusual positions (7-beta-d-guanine), bases shaped almost identically to the natural bases (4-aminobenzimidazole and 4,6-difluorobenzimidazole), bases shaped very differently than a natural base (e.g., 5- and 6-trifluoromethylbenzimidazole), and bases capable of forming just one Watson-Crick hydrogen bond with the template base (purine and 4-aminobenzimidazole). The only analogues that primase readily polymerized into primers (ITP and 3-deaza-ATP) were those capable of forming Watson-Crick hydrogen bonds with the template base. Thus, herpes primase appears to require the formation of Watson-Crick hydrogen bonds in order to efficiently polymerize a NTP. In contrast to primase's narrow specificity for NTP analogues, the DNA-dependent NTPase activity associated with the herpes primase-helicase complex exhibited very little specificity with respect to NTPs containing unnatural bases. The implications of these results with respect to the mechanism of the helicase-primase and current fidelity models are discussed.     

Intramolecular CH···O hydrogen bonds in the AI and BI DNA-like conformers of canonical nucleosides and their Watson-Crick pairs. Quantum chemical and AIM analysis

The aim of this work is to cast some light on the H-bonds in double-stranded DNA in its AI and BI forms. For this purpose, we have performed the MP2 and DFT quantum chemical calculations of the canonical nucleoside conformers, relative to the AI and BI DNA forms, and their Watson-Crick pairs, which were regarded as the simplest models of the double-stranded DNA. Based on the atoms-in-molecules analysis (AIM), five types of the CH···O hydrogen bonds, involving bases and sugar, were detected numerically from 1 to 3 per a conformer: C2'H···O5', C1'H···O2, C6H···O5', C8H···O5', and C6H···O4'. The energy values of H-bonds occupy the range of 2.3-5.6 kcal/mol, surely exceeding the kT value (0.62 kcal/mol). The nucleoside CH···O hydrogen bonds appeared to "survive" turns of bases against the sugar, sometimes in rather large ranges of the angle values, pertinent to certain conformations, which points out to the source of the DNA lability, necessary for the conformational adaptation in processes of its functioning. The calculation of the interactions in the dA·T nucleoside pair gives evidence, that additionally to the N6H···O4 and N1···N3H canonical H-bonds, between the bases adenine and thymine the third one (C2H···O2) is formed, which, though being rather weak (about 1 kcal/mol), satisfies the AIM criteria of H-bonding and may be classified as a true H-bond. The total energy of all the CH···O nontraditional intramolecular H-bonds in DNA nucleoside pairs appeared to be commensurable with the energy of H-bonds between the bases in Watson-Crick pairs, which implies their possible important role in the DNA shaping.     

Human DNA primase uses Watson-Crick hydrogen bonds to distinguish between correct and incorrect nucleoside triphosphates

Human DNA primase synthesizes short RNA primers that DNA polymerase alpha further elongates. Primase readily misincorporates the natural NTPs and will generate a wide variety of mismatches. In contrast, primase exhibited a remarkable resistance to polymerizing NTPs containing unnatural bases. This included bases whose shape was almost identical to the natural bases (4-aminobenzimidazole and 4,6-difluorobenzimidazole), bases shaped very differently than a natural base [e.g., 5- and 6-(trifluoromethyl)benzimidazole], bases much more hydrophobic than a natural base [e.g., 4- and 7-(trifluoromethyl)benzimidazole], bases of similar hydrophobicity as a natural base but with the Watson-Crick hydrogen-bonding groups in unusual positions (7-beta-D-guanine), and bases capable of forming only one Watson-Crick hydrogen bond with the template base (purine and 4-aminobenzimidazole). Primase only polymerized NTP analogues containing bases capable of forming hydrogen bonds between the equivalent of both N-1 and the exocyclic group at C-6 of a purine NTP (2-fluoroadenine, 2-chloroadenine, 3-deazaadenine, and hypoxanthine) and N-3 and the exocyclic group at C-4 of a pyrimidine. These data indicate that human primase requires the formation of Watson-Crick hydrogen bonds in order to polymerize a NTP, a situation very different than what is observed with some DNA polymerases. The implications of these results with respect to current theories of how polymerases discriminate between right and wrong (d)NTPs are discussed.