Insights into the DNA stabilizing contributions of a bicyclic cytosine analogue: crystal structures of DNA duplexes containing 7,8-dihydropyrido [2,3-d]pyrimidin-2-one

The incorporation of the bicyclic cytosine analogue 7,8-dihydropyrido[2,3-d]pyrimidin-2-one (X) into DNA duplexes results in a significant enhancement of their stability (3–4 K per modification). To establish the effects of X on the local hydrogen-bonding and base stacking interactions and the overall DNA conformation, and to obtain insights into the correlation between the structure and stability of X-containing DNA duplexes, the crystal structures of [d(CGCGAATT-X-GCG)]₂ and [d(CGCGAAT-X-CGCG)]₂ have been determined at 1.9–2.9 Å resolutions. In all of the structures, the analogue X base pairs with the purine bases on the opposite strands through Watson–Crick and/or wobble type hydrogen bonds. The additional ring of the X base is stacked on the thymine bases at the 5′-side and overall exhibits greatly enhanced stacking interactions suggesting that this is a major contribution to duplex stabilization.     

Importance of the sphingosine base double-bond geometry for the structural and thermodynamic properties of sphingomyelin bilayers

The precise role of the sphingosine base trans double bond for the unique properties of sphingomyelins (SMs), one of the main lipid components in raftlike structures of biological membranes, has not been fully explored. Several reports comparing the hydration, lipid packing, and hydrogen-bonding behaviors of SM and glycerophospholipid bilayers found remarkable differences overall. However, the atomic interactions linking the double-bond geometry with these thermodynamic and structural changes remained elusive. A recent report on ceramides, which differ from SMs only by their hydroxyl headgroup, has shown that replacing the trans double bond of the sphingosine base by cis weakens the hydrogen-bonding potential of these lipids and thereby alters their biological activity. Based on data from extensive (a total 0.75 μs) atomistic molecular dynamics simulations of bilayers composed of all-trans, all-cis, and a trans/cis (4:1 ratio) racemic mixture of sphingomyelin lipids, here we show that the trans configuration allows for the formation of significantly more hydrogen bonds than the cis. The extra hydrogen bonds enabled tighter packing of lipids in the all-trans and trans/cis bilayers, thus reducing the average area per lipid while increasing the chain order and the bilayer thickness. Moreover, fewer water molecules access the lipid-water interface of the all-trans bilayer than of the all-cis bilayer. These results provide the atomic basis for the importance of the natural sphingomyelin trans double-bond conformation for the formation of ordered membrane domains.     

Replication of non-hydrogen bonded bases by DNA polymerases: A mechanism for steric matching

Recent experiments have presented evidence that Watson–Crick hydrogen bonds in a base pair are not absolute requirements for efficient synthesis of that pair by DNA polymerase enzymes. Here we examine quantitative steady-state kinetic data from several published studies involving poorly hydrogen-bonding DNA base analogues and adducts, and analyze the results in terms of solvation, hydrogen bonding, and steric effects. We propose a mechanism that can explain the surprising lack of hydrogen-bonding requirement accompanied by significant selectivity in pairing. This hypothesis makes use of steric matching, enforced both by the tightly confined polymerase active site and by the DNA backbone, as a chief factor determining nucleotide selection during DNA synthesis. The results also suggest that hydrogen bonds from bases to water (solvation) may be important in increasing the effective size of DNA bases, which may help prevent misinsertion of small bases opposite each other. © 1998 John Wiley & Sons, Inc. Biopoly 48: 3–17, 1998     

Structural diversity and isomorphism of hydrogen-bonded base interactions in nucleic acids

The wide structural diversity of RNA results in part from the diversity of non-Watson-Crick interactions between bases. To examine the repertoire of possible hydrogen bond interactions among bases, we computed databases of base-pairs and base-triples by systematically matching all possible hydrogen-bond donors and acceptors between bases and evaluating the geometries of each planar configuration. For base-pairs, we find 53 arrangements having at least two hydrogen bonds, including 23 pairs with protonated bases that have not previously been modeled. A comparison with experimentally observed base-pairs reveals an unexpected G:U pair recently observed in the ribosome. For base-triples, we find 840 arrangements in which the three bases are constrained by a total of at least three hydrogen bonds. Base-triples in particular exhibit a wide range of structural diversity, suggesting how compact or elongated nucleic acid structures may be constructed using different hydrogen-bonding patterns. Base-pair and base-triple conformations were systematically compared to identify structurally isomorphic combinations, and the experimentally observed arrangements within double and triple helices are among the most isomorphic. Unexpectedly, however, other combinations in the database are even more isomorphic, including several in which all-purine arrangements overlap with all-pyrimidine arrangements. These studies highlight some of the combinatoric and geometric versatility of base interactions and help provide a framework for analyzing and modeling isomorphic interactions and potentially for designing novel nucleic acid structures.     

Structural Studies of DNA Fragments: The G·T Wobble Base Pair in A,B and Z DNA; The G·A Base Pair in B-DNA

Abstract

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G·T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with 02 of T and 06 of G with N3 of T. The X-ray analyses establish that the G·T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G·A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Solution structure of DNA containing α-OH-PdG: the mutagenic adduct produced by acrolein

Acrolein is a cell metabolic product and a main component of cigarette smoke. Its reaction with DNA produces two guanine lesions γ-OH-PdG, a major adduct that is nonmutagenic in mammalian cells, and the positional isomer α-OH-PdG. We describe here the solution structure of a short DNA duplex containing a single α-OH-PdG lesion, as determined by solution NMR spectroscopy and restrained molecular dynamics simulations. The spectroscopic data show a mostly regular right-handed helix, locally perturbed at its center by the presence of the lesion. All undamaged residues of the duplex are in anti orientation, forming standard Watson–Crick base-pair alignments. Duplication of proton signals near the damaged site differentiates two enantiomeric duplexes, thus establishing the exocyclic nature of the lesion. At the lesion site, α-OH-PdG rotates to a syn conformation, pairing to its counter cytosine residue that is protonated at pH 5.9. Three-dimensional models produced by restrained molecular dynamics simulations show different hydrogen-bonding patterns between the lesion and its cytosine partner and identify further stabilization of α-OH-PdG in a syn conformation by intra-residue hydrogen bonds. We compare the α-OH-PdG•dC duplex structure with that of duplexes containing the analogous lesion propano-dG and discuss the implications of our findings for the mutagenic bypass of acrolein lesions.     

The hydrogen bonding of cytosine with guanine: calorimetric and 1H-nmr analysis of the molecular interactions of nucleic acid bases

The enthalpy of hydrogen-bond formation between guanine (G) and cytosine (C) in o-dichlorobenzene and in chloroform at 25°C has been determined by direct calorimetric measurement. We derivatized 2′-deoxyguanosine and 2′-deoxycytidine at the 5′- and 3′-hydroxyls with triisopropylsilyl groups; these groups increase the solubility of the nucleic acid bases in nonaqueous solvents. Such derivatization also prevents the ribose hydroxyls from forming hydrogen bonds. Consequently, hydrogen-bond formation in our system is primarily between the bases, and to a lesser extent, between base and solvent, and can be measured directly with calorimetry. To obtain the data on base-pair formation, we first took into account the contributions from self-association of each base, and where possible, have determined the ΔH of self-association. From isoperibolic titration calorimetry, our measured ΔH of C₂ formation in chloroform is ?1.7 kcal/mol of C. Our measured ΔH of C:G base-pair formation in o-dichlorobenzene is ?6.65 ± 0.32 kcal/mol. Since o-dichlorobenzene does not form hydrogen bonds, the ΔH of C:G base-pair formation in this solvent represents the ΔH of the hydrogen-bonding interaction of C with G in a nonassociating solvent. In contrast, our measured ΔH of C:G base-pair formation in chloroform is ?5.77 ± 0.20 kcal/mol; thus, the absolute value of the enthalpy of hydrogen bonding in the C:G base pair is greater in o-dichlorobenzene than in chloroform. Since chloroform is a solvent known to form hydrogen bonds, the decrease in enthalpic contribution to C:G base pairing in chloroform is due to the formation of hydrogen bonds between the bases and the solvent. The ΔH of hydrogen bonding of G with C reported here differs from previous indirect estimates: Our measurements indicate the ΔH is 50% less in magnitude than the ΔH based on spectroscopic measurements of the extent of interaction. We have also observed that the enthalpy of hydrogen bonding of C with G in chloroform is greater when G is in excess than when C is in excess. This increased heat is due to the formation of C:G_{n > 1} complexes that we have observed using ¹H-nmr. Although C:G₂ structures have previously been observed in triple-stranded polymeric nucleic acids, higher order structures have not been observed between C and G monomers in nonaqueous solvents until now. By using monomers as a model system to investigate hydrogen-bonding interactions in DNA and RNA, we have obtained the following results: A direct measurement of the ΔH of hydrogen bonding in the C:G complex in two nonaqueous solvents, and the first observation of C:G_{n > 1} complexes between monomers. These results reinforce the importance of hydrogen bonding in the stabilization of various nucleic acid secondary and tertiary structures.     

Structural studies of DNA fragments: the G.T wobble base pair in A, B and Z DNA; the G.A base pair in B-DNA

The crystal structures of five double helical DNA fragments containing non-Watson-Crick complementary base pairs are reviewed. They comprise four fragments containing G.T base pairs: two deoxyoctamers d(GGGGCTCC) and d(GGGGTCCC) which crystallise as A type helices; a deoxydodecamer d(CGCGAATTTGCG) which crystallises in the B-DNA conformation; and the deoxyhexamer d(TGCGCG), which crystallises as a Z-DNA helix. In all four duplexes the G and T bases form wobble base pairs, with bases in the major tautomer forms and hydrogen bonds linking N1 of G with O2 of T and O6 of G with N3 of T. The X-ray analyses establish that the G.T wobble base pair can be accommodated in the A, B or Z double helix with minimal distortion of the global conformation. There are, however, changes in base stacking in the neighbourhood of the mismatched bases. The fifth structure, d(CGCGAATTAGCG), contains the purine purine mismatch G.A where G is in the anti and A in the syn conformation. The results represent the first direct structure determinations of base pair mismatches in DNA fragments and are discussed in relation to the fidelity of replication and mismatch recognition.     

Infrared spectra and conformation of monodisperse Oligo-γ-Benzyl-L-Glutamates in solution

The ir spectra of σ–nitrophenylthio-tetra- and hexa-γ-benzyl-L -glutamate ethylamides were measured at different concentrations in chloroform and ethylene dichloride. The molar extinction coefficients of two bands, each for the amide I and A modes, were observed as indicating the content of β structure and the fraction of hydrogen bonds; they were analyzed for elucidating the hydrogen-bonded state of peptide residues in the σ and β conformations of oligopeptides. While the content of β structure of the tetra-peptide increases with increasing concentration, the hexapeptide is in the β conformation above the critical concentration only. The fraction of hydrogen bonds remains finite even at infinite dilution or below the critical concentration, indicating the intramolecular hydrogen-bonding in the σ—form. As the fundamental structure of folded forms having only intramolecular hydrogen bonds, the 2₇ ribbon is most likely. With increasing concentration or above the critical concentration, the extended forms are stabilized by the intermolecular hydrogen bonds between residues of the β-form. The β-form is present only when intermolecular hydrogen bonds link two residues in an antiparallel way. Possible structures of the oligopeptides in the σ– and β conformations in the two solvents are described briefly.     

NMR structure of the DNA decamer duplex containing double T*G mismatches of cis-syn cyclobutane pyrimidine dimer: implications for DNA damage recognition by the XPC-hHR23B complex

The cis–syn cyclobutane pyrimidine dimer (CPD) is a cytotoxic, mutagenic and carcinogenic DNA photoproduct and is repaired by the nucleotide excision repair (NER) pathway in mammalian cells. The XPC–hHR23B complex as the initiator of global genomic NER binds to sites of certain kinds of DNA damage. Although CPDs are rarely recognized by the XPC–hHR23B complex, the presence of mismatched bases opposite a CPD significantly increased the binding affinity of the XPC–hHR23B complex to the CPD. In order to decipher the properties of the DNA structures that determine the binding affinity for XPC–hHR23B to DNA, we carried out structural analyses of the various types of CPDs by NMR spectroscopy. The DNA duplex which contains a single 3′ T·G wobble pair in a CPD (CPD/GA duplex) induces little conformational distortion. However, severe distortion of the helical conformation occurs when a CPD contains double T·G wobble pairs (CPD/GG duplex) even though the T residues of the CPD form stable hydrogen bonds with the opposite G residues. The helical bending angle of the CPD/GG duplex was larger than those of the CPD/GA duplex and properly matched CPD/AA duplex. The fluctuation of the backbone conformation and significant changes in the widths of the major and minor grooves at the double T·G wobble paired site were also observed in the CPD/GG duplex. These structural features were also found in a duplex that contains the (6–4) adduct, which is efficiently recognized by the XPC–hHR23B complex. Thus, we suggest that the unique structural features of the DNA double helix (that is, helical bending, flexible backbone conformation, and significant changes of the major and/or minor grooves) might be important factors in determining the binding affinity of the XPC–hHR23B complex to DNA.     

Mechanism of DNA Recognition by the Restriction Enzyme EcoRV

EcoRV, a restriction enzyme in Escherichia coli, destroys invading foreign DNA by cleaving it at the center step of a GATATC sequence. In the EcoRV-cognate DNA crystallographic complex, a sharp kink of 50° has been found at the center base-pair step (TA). Here, we examine the interplay between the intrinsic propensity of the cognate sequence to kink and the induction by the enzyme by performing all-atom molecular dynamics simulations of EcoRV unbound and interacting with three DNA sequences: the cognate sequence, GATATC (TA); the non-cognate sequence, GAATTC (AT); and with the cognate sequence methylated on the first adenine GA_CH3TATC (TA-CH₃). In the unbound EcoRV, the cleft between the two C-terminal subdomains is found to be open. Binding to AT narrows the cleft and forms a partially bound state. However, the intrinsic bending propensity of AT is insufficient to allow tight binding. In contrast, the cognate TA sequence is easier to bend, allowing specific, high-occupancy hydrogen bonds to form in the complex. The absence of cleavage for this methylated sequence is found to arise from the loss of specific hydrogen bonds between the first adenine of the recognition sequence and Asn185. On the basis of the results, we suggest a three-step recognition mechanism. In the first step, EcoRV, in an open conformation, binds to the DNA at a random sequence and slides along it. In the second step, when the two outer base pairs, GAxxTC, are recognized, the R loops of the protein become more ordered, forming strong hydrogen-bonding interactions, resulting in a partially bound EcoRV-DNA complex. In the third step, the flexibility of the center base pair is probed, and in the case of the full cognate sequence the DNA bends, the complex strengthens and the protein and DNA interact more closely, allowing cleavage.     

Ned Seeman and the prediction of amino acid-basepair motifs mediating protein-nucleic acid recognition

Fifty years ago, the first atomic-resolution structure of a nucleic acid double helix, the mini-duplex (ApU)₂, revealed details of basepair geometry, stacking, sugar conformation, and backbone torsion angles, thereby superseding earlier models based on x-ray fiber diffraction, including the original DNA double helix proposed by Watson and Crick. Just 3 years later, in 1976, Ned Seeman, John Rosenberg, and Alex Rich leapt from their structures of mini-duplexes and H-bonding motifs between bases in small-molecule structures and transfer RNA to predicting how proteins could sequence specifically recognize double helix nucleic acids. They proposed interactions between amino acid side chains and nucleobases mediated by two hydrogen bonds in the major or minor grooves. One of these, the arginine-guanine pair, emerged as the most favored amino acid-base interaction in experimental structures of protein-nucleic acid complexes determined since 1986. In this brief review we revisit the pioneering work by Seeman et al. and discuss the importance of the arginine-guanine pairing motif.     

Conformational changes of the phenyl and naphthyl isocyanate-DNA adducts during DNA replication and by minor groove binding molecules

DNA lesions produced by aromatic isocyanates have an extra bulky group on the nucleotide bases, with the capability of forming stacking interaction within a DNA helix. In this work, we investigated the conformation of the 2′-deoxyadenosine and 2′-deoxycytidine derivatives tethering a phenyl or naphthyl group, introduced in a DNA duplex. The chemical modification experiments using KMnO₄ and 1-cyclohexyl-3 -(2-morpholinoethyl) carbodiimide metho-p-toluenesulfonate have shown that the 2′-deoxycytidine lesions form the base pair with guanine while the 2′-deoxyadenosine lesions have less ability of forming the base pair with thymine in solution. Nevertheless, the kinetic analysis shows that these DNA lesions are compatible with DNA ligase and DNA polymerase reactions, as much as natural DNA bases. We suggest that the adduct lesions have a capability of adopting dual conformations, depending on the difference in their interaction energies between stacking of the attached aromatic group and base pairing through hydrogen bonds. It is also presented that the attached aromatic groups change their orientation by interacting with the minor groove binding netropsin, distamycin and synthetic polyamide. The nucleotide derivatives would be useful for enhancing the phenotypic diversity of DNA molecules and for exploring new non-natural nucleotides.     

Influence of local sequence context on damaged base conformation in human DNA polymerase ι: molecular dynamics studies of nucleotide incorporation opposite a benzo[a]pyrene-derived adenine lesion

Human DNA polymerase ι is a lesion bypass polymerase of the Y family, capable of incorporating nucleotides opposite a variety of lesions in both near error-free and error-prone bypass. With undamaged templating purines polymerase ι normally favors Hoogsteen base pairing. Polymerase ι can incorporate nucleotides opposite a benzo[a]pyrene-derived adenine lesion (dA*); while mainly error-free, the identity of misincorporated bases is influenced by local sequence context. We performed molecular modeling and molecular dynamics simulations to elucidate the structural basis for lesion bypass. Our results suggest that hydrogen bonds between the benzo[a]pyrenyl moiety and nearby bases limit the movement of the templating base to maintain the anti glycosidic bond conformation in the binary complex in a 5′-CAGA*TT-3′ sequence. This facilitates correct incorporation of dT via a Watson−Crick pair. In a 5′-TTTA*GA-3′ sequence the lesion does not form these hydrogen bonds, permitting dA* to rotate around the glycosidic bond to syn and incorporate dT via a Hoogsteen pair. With syn dA*, there is also an opportunity for increased misincorporation of dGTP. These results expand our understanding of the versatility and flexibility of polymerase ι and its lesion bypass functions in humans.     

Structural basis of DNA flexibility

It has been recently shown by us, on the basis of crystal structure database that the flexibility of B-DNA double helices depends significantly on their base sequence. Our model building studies further indicated that the existence of bifurcated cross-strand hydrogen bonds between successive base pairs is possibly the main factor behind the sequence directed DNA flexibility. These cross-strand hydrogen bonds are, of course, weaker than the usual Watson-Crick hydrogen bonds and their bond geometry is characterized by relatively larger bond lengths and smaller bond angles. We have tried to improve our model structures by incorporating non-planarity of the amino groups in DNA bases due to the presence of lone pair electrons at the nitrogen atoms. Energy minimization studies have been carried out by using different quantum chemical methods, whereby it is found that in all cases of N-H....O type cross-strand hydrogen bonds, the bond geometry improves significantly. In the cases of N-H....N type hydrogen bonds, however, no such consistent improvements can be noticed. Perhaps the true picture would emerge only if all the other interactions present in the DNA macromolecule could be appropriately taken into account.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

How alcohol chain-length and concentration modulate hydrogen bond formation in a lipid bilayer

Molecular dynamics simulations are used to measure the change in properties of a hydrated dipalmitoylphosphatidylcholine bilayer when solvated with ethanol, propanol, and butanol solutions. There are eight oxygen atoms in dipalmitoylphosphatidylcholine that serve as hydrogen bond acceptors, and two of the oxygen atoms participate in hydrogen bonds that exist for significantly longer time spans than the hydrogen bonds at the other six oxygen atoms for the ethanol and propanol simulations. We conclude that this is caused by the lipid head group conformation, where the two favored hydrogen-bonding sites are partially protected between the head group choline and the sn-2 carbonyl oxygen. We find that the concentration of the alcohol in the ethanol and propanol simulations does not have a significant influence on the locations of the alcohol/lipid hydrogen bonds, whereas the concentration does impact the locations of the butanol/lipid hydrogen bonds. The concentration is important for all three alcohol types when the lipid chain order is examined, where, with the exception of the high-concentration butanol simulation, the alcohol molecules having the longest hydrogen-bonding relaxation times at the favored carbonyl oxygen acceptor sites also have the largest order in the upper chain region. The lipid behavior in the high-concentration butanol simulation differs significantly from that of the other alcohol concentrations in the order parameter, head group rotational relaxation time, and alcohol/lipid hydrogen-bonding location and relaxation time. This appears to be the result of the system being very near to a phase transition, and one occurrence of lipid flip-flop is seen at this concentration.