High-resolution crystal structure of the intramolecular d(TpA) thymine-adenine photoadduct and its mechanistic implications

A high-resolution crystal structure is reported for d(TpA)*, the intramolecular thymine–adenine photoadduct that is produced by direct ultraviolet excitation of the dinucleoside monophosphate d(TpA). It confirms the presence of a central 1,3-diazacyclooctatriene ring linking the remnants of the T and A bases, as previously deduced from heteronuclear NMR measurements by Zhao et al. (The structure of d(TpA)*, the major photoproduct of thymidylyl-(3′-5′)-deoxyadenosine. Nucleic Acids Res., 1996, 24, 1554–1560). Within the crystal, the d(TpA)* molecules exist as zwitterions with a protonated amidine fragment of the eight-membered ring neutralizing the charge of the internucleotide phosphate monoanion. The absolute configuration at the original thymine C5 and C6 atoms is determined as 5S,6R. This is consistent with d(TpA)* arising by valence isomerization of a precursor cyclobutane photoproduct with cis–syn stereochemistry that is generated by [2 + 2] photoaddition of the thymine 5,6-double bond across the C6 and C5 positions of adenine. This mode of photoaddition should be favoured by the stacked conformation of adjacent T and A bases in B-form DNA. It is probable that the primary photoreaction is mechanistically analogous to pyrimidine dimerization despite having a much lower quantum yield.     

The structure of d(TpA), the major photoproduct of thymidylyl-(3'5')-deoxyadenosine.

Irradiation of the dinucleotide TpdA and TA-containing oligonucleotides and DNA produces the TA* photoproduct which was proposed to be the [2+2] cyclo-addition adduct between the C5-C6 double bonds of the T and the A [Bose,S.N., Kumar,S., Davies,R.J.H., Sethi,S.K. and McCloskey,J.A. (1984) Nucleic Acids Res. 12, 7929-7947]. The proposed structure was based on a variety of spectroscopic and chemical degradation studies, and the assignment of a trans-syn-I stereochemistry was based on an extensive 1H-NMR and molecular modeling study of the dinucleotide adduct [Koning,T.M.G., Davies,R.J.H. and Kaptein,R. (1990) Nucleic Acids Res. 18, 277-284]. However, a number of properties of TA* are not in accord with the originally proposed structure, and prompted a re-evaluation of the structure. To assign the 13C spectrum and establish the bond connectivities of the TA* photoproduct of TpdA [d(TpA)*], 1H-13C heteronuclear multiple-quantum coherence (HMQC) and heteronuclear multiple bond correlation (HMBC) spectra were obtained. The 13C shifts and connectivities were found to be inconsistent with the originally proposed cyclobutane ring fusion between the thymine and adenine, but could be explained by a subsequent ring-expansion reaction to give an eight-membered ring valence isomer. The new structure for the d(TpA)* resolves the inconsistencies with the originally proposed structure, and could have a stereochemistry that arises from the anti, anti glycosyl conformation found in B form DNA.     

1H NMR assignment and melting temperature study of cis-syn and trans-syn thymine dimer containing duplexes of d(CGTATTATGC).d(GCATAATACG)

The preparation and spectroscopic characterization of duplex decamers containing site-specific cis-syn and trans-syn thymine dimers are described. Three duplex decamers, d(CGTATTATGC).d(GCATAATACG), d(CGTAT[c,s]TATGC).d(GCATAATACG), and d(CGTAT[t,s]TATGC).d(GCATAATACG), were prepared by solid-phase phosphoramidite synthesis utilizing cis-syn and trans-syn cyclobutane thymine dimer building blocks (Taylor et al., 1987; Taylor & Brockie, 1988). NMR spectra (500 MHz 2D 1H and 202 MHz 1D 31P) were obtained in "100%" D2O at 10 degrees C, and 1D exchangeable 1H spectra were obtained in 10% D2O at 10 degrees C. 1H NMR assignments for H5, H6, H8, CH3, H1', H2', and H2" were made on the basis of standard sequential NOE assignment strategies and verified in part by DQF COSY data. Comparison of the chemical shift data suggests that the helix structure is perturbed more to the 3'-side of the cis-syn dimer and more to the 5'-side of the trans-syn dimer. Thermodynamic parameters for the helix in equilibrium coil equilibrium were obtained by two-state, all or none, analysis of the melting behavior of the duplexes. Analysis of the temperature dependence of the T5CH3 1H NMR signal gave delta H = 44 +/- 4 kcal and delta S = 132 +/- 13 eu for the trans-syn duplex. Analysis of the concentration and temperature dependence of UV spectra gave delta H = 64 +/- 6 kcal and delta S = 178 +/- 18 eu for the parent duplex and delta H = 66 +/- 7 kcal and delta S = 189 +/- 19 eu for cis-syn duplex. It was concluded that photodimerization of the dTpdT unit to give the cis-syn product causes little perturbation of the DNA whereas dimerization to give the trans-syn product causes much greater perturbation, possibly in the form of a kink or dislocation at the 5'-side of the dimer.     

NMR studies of the exocyclic 1,N6-ethenodeoxyadenosine adduct (epsilon dA) opposite thymidine in a DNA duplex. Nonplanar alignment of epsilon dA(anti) and dT(anti) at the lesion site

Two-dimensional proton NMR studies are reported on the complementary d(C-A-T-G-T-G-T-A-C).d(G-T-A-C-epsilon A-C-A-T-G) nonanucleotide duplex (designated epsilon dA.dT 9-mer duplex) containing 1,N6-ethenodeoxyadenosine (epsilon dA), a carcinogen-DNA adduct, positioned opposite thymidine in the center of the helix. Our NMR studies have focused on the conformation of the epsilon dA.dT 9-mer duplex at neutral pH with emphasis on defining the alignment at the dT5.epsilon dA14 lesion site. The through-space NOE distance connectivities establish that both dT5 and epsilon dA14 adopt anti glycosidic torsion angles, are directed into the interior of the helix, and stack with flanking Watson-Crick dG4.dC15 and dG6.dC13 pairs. Furthermore, the d(G4-T5-G6).d(C13-epsilon A14-C15) trinucleotide segment centered about the dT5.epsilon dA14 lesion site adopts a right-handed helical conformation in solution. Energy minimization computations were undertaken starting from six different alignments of dT5(anti) and epsilon dA14(anti) at the lesion site and were guided by distance constraints defined by lower and upper bounds estimated from NOESY data sets on the epsilon dA.dT 9-mer duplex. Two families of energy-minimized structures were identified with the dT5 displaced toward either the flanking dG4.dC15 or the dG6.dC13 base pair. These structures can be differentiated on the basis of the observed NOEs from the imino proton of dT5 to the imino proton of dG4 but not dG6 and to the amino protons of dC15 but not dC13 that were not included in the constraints data set used in energy minimization. Our NMR data are consistent with a nonplanar alignment of epsilon dA14(anti) and dT5(anti) with dT5 displaced toward the flanking dG4.dC15 base pair within the d(G4-T5-G6).d(C13-epsilon A14-C15) segment of the epsilon dA.dT 9-mer duplex.     

NMR studies of abasic sites in DNA duplexes: deoxyadenosine stacks into the helix opposite acyclic lesions

Proton and phosphorus NMR studies are reported for two complementary nonanucleotide duplexes containing acyclic abasic sites. The first duplex, d(C-A-T-G-A-G-T-A-C).d(G-T-A-C-P-C-A-T-G), contains an acyclic propanyl moiety, P, located opposite a deoxyadenosine at the center of the helix (designated APP 9-mer duplex). The second duplex, d(C-A-T-G-A-G-T-A-C).d(G-T-A-C-E-C-A-T-G), contains a similarly located acyclic ethanyl moiety, E (designated APE 9-mer duplex). The ethanyl moiety is one carbon shorter than the natural carbon-phosphodiester backbone of a single nucleotide unit of DNA. The majority of the exchangeable and nonexchangeable base and sugar protons in both the APP 9-mer and APE 9-mer duplexes, including those at the abasic site, have been assigned by recording and analyzing two-dimensional phase-sensitive NOESY data sets in H2O and D2O solution between -5 and 5 degrees C. These spectroscopic observations establish that A5 inserts into the helix opposite the abasic site (P14 and E14) and stacks between the flanking G4.C15 and G6.C13 Watson-Crick base pairs in both the APP 9-mer and APE 9-mer duplexes. The helix is right-handed at and adjacent to the abasic site, and all glycosidic torsion angles are anti in both 9-mer duplexes. Proton NMR parameters for the APP 9-mer and APE 9-mer duplexes are similar to those reported previously for the APF 9-mer duplex (F = furan) in which a cyclic analogue of deoxyribose was embedded in an otherwise identical DNA sequence [Kalnik, M. W., Chang, C. N., Grollman, A. P., & Patel, D. J. (1988) Biochemistry 27, 924-931]. These proton NMR experiments demonstrate that the structures at abasic sites are very similar whether the five-membered ring is open or closed or whether the phosphodiester backbone is shortened by one carbon atom. Phosphorus spectra of the APP 9-mer and APE 9-mer duplexes (5 degrees C) indicate that the backbone conformation is similarly perturbed at three phosphodiester backbone torsion angles. These same torsion angles are also distorted in the APF 9-mer but assume a different conformation than those in the APP 9-mer and APE 9-mer duplexes.     

NMR studies of bipyrimidine cyclobutane photodimers

Cyclobutane-type photodimers of dinucleoside monophosphates dCpdT, dTpdC and dTpdT were prepared by ultraviolet irradiation in the presence of acetophenone as photosensitizer. The cytosine-containing derivatives were found to deaminate forming uracil products. Using one- and two-dimensional NMR, the photoproducts were characterized as cis-syn and trans-syn cyclobutane photodimers. On the basis of NOE data the structures of the cis-syn and trans-syn products of dUpdT were determined using distance-geometry and restrained-energy-minimization methods. The cis-syn structures showed (high-ANTI/SYN)/high-ANTI glycosidic linkages while the trans-syn structures were in the SYN-ANTI region. The backbone conformations of both structures were in fair agreement with the coupling-constant-data. The trans-syn structures were found to be very rigid and similar in all three products. For the three cis-syn structures more conformational freedom and more variation among the three structures was observed.     

Solution-state structure of the Dewar pyrimidinone photoproduct of thymidylyl-(3'----5')-thymidine

The preparation, spectroscopic investigation, structure determination, conformational analysis, and modeling of the Dewar pyrimidinone photoproduct of thymidylyl-(3'----5')-thymidine, previously referred to as TpT3 [Johns, H. E., Pearson, M. L., LeBlanc, J. C., & Heilleiner, C. W. (1964) J. Mol. Biol. 9, 503-524], is described. TpT3 was prepared in quantitative yield by photolysis of an aqueous solution of the (6-4) photoproduct of TpT with Pyrex-filtered medium-pressure mercury arc light. TpT3 was analyzed by FAB MS, IR, UV, and 1H, 13C, and 31P NMR spectroscopy. The spectroscopic data led to the conclusion that TpT3 results from the photoisomerization of the pyrimidinone ring of the (6-4) product of TpT to its Dewar valence isomer. Torsion angle and interproton distance information derived from coupling constants and NOE data was used to constrain ring conformation searches by utilizing the SYBYL molecular modeling program subroutine SEARCH. Sets of angles derived from the ring search procedure were then used to construct structures whose geometries were optimized by the energy-minimization subroutine MAXIMIN. A two-state model for the solution-state structure of the Dewar photoproduct was chosen which was energetically sound, fit the experimental coupling constants with an RMS deviation of 1.15 Hz, and was consistent with the NOE data. The model for the Dewar photoproduct was compared to a model for the (6-4) photoproduct and the TpT subunits of the Dickerson dodecamer structure by a least-squares fitting procedure. It was concluded that the Dewar photoproduct more closely resembles a B-form TpT unit than does the (6-4) photoproduct.     

Molecular Mechanics and Dynamics of an Abasic Frameshift in DNA and Comparison to NMR Data

Abstract

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide:

5′d(C1 P1 G2 P2 G3 P3 dr4 P4 G5 P5 G6 P6 C7)3′

3′d(G13P12C12PllCll P10 C10 P9 C9 P8 G8)5′

where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2′ endo or C3′ endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4.

In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium.

Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.     

Crystal and molecular structure of a DNA duplex containing the carcinogenic lesion O6-methylguanine

The crystal and molecular structure of the first DNA duplex containing the carcinogenic lesion O6MeG has been determined to a resolution of 1.9 A and refined to an R factor of 19%. (d[CGC-(O6Me)GCG])2 crystallizes in the left-handed Z DNA form and has crystal parameters and conformational features similar to those of the parent sequence [d(CG)3]2. The methyl groups on O6 of G4 and G10 have C5-C6-O6-O6Me torsion angles of 73 degrees and 56 degrees, respectively, and protrude onto the major groove surface. The base-pairing conformation for the methylated G.C base pairs is of the Watson-Crick type as opposed to a wobble-type conformation that had been proposed in a B DNA fragment. As in other Z DNA structures, a spine of hydration is seen in the minor groove.     

Adenine photodimerization in deoxyadenylate sequences: elucidation of the mechanism through structural studies of a major d(ApA) photoproduct.

The mechanism of the photodimerization of adjacent adenine bases on the same strand of DNA has been elucidated by determining the structure of one of the two major photoproducts that are formed by UV irradiation of the deoxydinucleoside monophosphate d(ApA). The photoproduct, denoted d(ApA)*, corresponds to a species of adenine photodimer first described by P?rschke (P?rschke, D. (1973) J.Am.Chem.Soc. 95, 8440-8446). From a detailed examination of its chemical and spectroscopic properties, including comparisons with the model compound N-cyano-N1-(1-methylimidazol-5-yl)formamidine, it is deduced that d(ApA)* contains a deoxyadenosine unit covalently linked through its C(8) position to C(4) of an imidazole N(1) deoxyribonucleoside moiety bearing an N-cyanoformamidino substituent at C(5). On treatment with acid, d(ApA)* is degraded with high specificity to 8-(5-amino-imidazol-4-yl)adenine whose identity has been confirmed by independent chemical synthesis. It is concluded that the primary event in adenine photodimerization entails photoaddition of the N(7)-C(8) double bond of the 5'-adenine across the C(6) and C(5) positions of the 3'-adenine. The azetidine species thus generated acts as a common precursor to both types of d(ApA) photoproduct which are formed from it by competing modes of azetidine ring fission.     

Antigen structural requirements for recognition by a cyclobutane thymine dimer-specific monoclonal antibody.

A monoclonal antibody (TDM-2) specific to a UV-induced cyclobutane pyrimidine dimer (T[cis-syn]T) has previously been established; however,the immunization had used UV-irradiated calf-thymus DNA containing a heterogeneous mixture of photoproduct sites. We investigated here the structural requirements of antigen recognition by the antibody using chemically synthesized antigen analogs. TDM-2 bound with cis-syn,but not trans-syn thymine dimer,and could bind strongly with four nucleotide analogs in which the cis-syn pyrimidine dimer was located in the center. Antigen analogs containing abasic linkers at the 5'- or 3'-side of the cis-syn cyclobutane pyrimidine dimer were synthesized and tested for binding to TDM-2. The results indicated that TDM-2 recognizes not only the cyclobutane ring but also both the 5'- and 3'-side nucleosides of the cyclobutane dimer. Furthermore,it was proved that either the 5'- or 3'-side phosphate group at a cyclobutane dimer site was absolutely required for the affinity to TDM-2. The antibody showed a strong binding to single stranded DNA but indicated little binding to double stranded DNA.     

NMR and molecular modeling studies on two glycopeptides from the carbohydrate-protein linkage region of connective tissue proteoglycans.

Complete 1H and 13C NMR assignments are reported for two glycopeptides representing the carbohydrate-protein linkage region of connective tissue proteoglycans. These glycopeptides are the octasaccharide hexapeptide, Ser(GlcpAbeta(1-->3) Galpbeta(1-->3)Galpbeta(1-->4)Xylpbeta)-Gly-Ser-Gly-Se r (GlcpAbeta(1-->3)Galpbeta(1-->3)Galpbeta(1-->4)Xylp beta)-Gly (1), and the tetrasaccharide dipeptide, Ser(GlcpAbeta(1-->3)Galpbeta(1-->3)Galpbeta(1-->4)X ylpbeta)-Gly (2). The vicinal coupling constant data show that the monosaccharide residues adopt4 C 1 chair conformations. Distance geometry/simulated annealing calculations using 2D NOESY derived distance constraints yielded a single family of structures for the tetrasaccharide moiety, with well defined interglycosidic linkage conformations. The straight phi torsion angles of the glycosidic C1'-O1 bonds showed a strict preference for the -sc range whereas the psi torsion angles (O1-Cn) exhibited dependence upon the interglycosidic linkage position (-ac for beta(1-->3) linkage, +ac for beta(1-->4) linkage). The predominant conformation about the glycopeptide bond is straight phi = -sc and psi = +ac. The presence of strong daN (i, i+1) NOE contacts, and the general absence of dNN (i, i+1) contacts (except for a weak Ser-5/Gly-6 dNN contact) and the dbN (i, i+1) contacts (except for Ser-1/Gly-2) in the ROESY spectrum, suggest that the backbone for 1 is predominantly in an extended conformation. A comparison of the ROESY data for 1 with those obtained from the unglycosylated hexapeptide (3) of the same sequence suggests that glycosylation has only a marginal influence on the backbone conformation of the hexapeptide.     

NMR studies of the exocyclic 1,N6-ethenodeoxyadenosine adduct (epsilon dA) opposite deoxyguanosine in a DNA duplex. Epsilon dA(syn).dG(anti) pairing at the lesion site

Proton NMR studies are reported on the complementary d(C-A-T-G-G-G-T-A-C).d(G-T-A-C-epsilon A-C-A-T-G) nonanucleotide duplex (designated epsilon dA.dG 9-mer duplex), which contains exocyclic adduct 1,N6-ethenodeoxyadenosine positioned opposite deoxyguanosine in the center of the helix. The present study focuses on the alignment of dG5 and epsilon dA14 at the lesion site in the epsilon dA.dG 9-mer duplex at neutral pH. This alignment has been characterized by monitoring the NOEs originating from the NH1 proton of dG5 and the H2, H5, and H7/H8 protons of epsilon dA14 in the central d(G4-G5-G6).d(C13-epsilon A14-C15) trinucleotide segment of the epsilon dA.dG 9-mer duplex. These NOE patterns establish that epsilon dA14 adopts a syn glycosidic torsion angle that positions the exocyclic ring toward the major groove edge while all the other bases including dG5 adopt anti glycosidic torsion angles. We detect a set of intra- and interstrand NOEs between protons (exchangeable and nonexchangeable) on adjacent residues in the d(G4-G5-G6).d(C13-epsilon A14-C15) trinucleotide segment which establish formation of right-handed helical conformations on both strands and stacking of the dG5(anti).epsilon dA14(syn) pair between stable dG4.dC15 and dG6.dC13 pairs. The energy-minimized conformation of the central d(G4-G5-G6).d(C13-epsilon A14-C15) segment establishes that the dG5(anti).epsilon dA14(syn) alignment is stabilized by two hydrogen bonds from the NH1 and NH2-2 of dG5(anti) to N9 and N1 of epsilon dA14(syn), respectively.(ABSTRACT TRUNCATED AT 250 WORDS)     

Molecular mechanics and dynamics of an abasic frameshift in DNA and comparison to NMR data

In a previous publication (Ph. Cuniasse, L.C. Sowers, R. Eritja, B. Kaplan, M.F. Goodman, J.A.H. Cognet, M. Le Bret, W. Guschlbauer and G.V. Fazakerley, Biochemistry 28, 2018 (1989), we determined by two dimensional NMR studies and molecular mechanics calculations the three-dimensional structure of a non-selfcomplementary oligonucleotide: [sequence; see text] where dr, at the center of the first strand, is a model abasic site. In order to explain all the results arising from NMR measurements, we found that an equilibrium between two conformations was necessary. These conformations differ mainly by the sugar pucker of G5 which is C2' endo or C3' endo. The latter is stabilized by addition of counterions between phosphate residues P3 and P4. In this paper, we have constructed systematically, all possible structures as a function of torsion angles delta of dr4 and of G5 by molecular mechanics in the presence or absence of counterions. Since these conformations were not forced with NMR distance measurements, this method allows detailed comparisons between all possible conformations and NMR data. Maps of contour lines of the potential energy, of fits to NMR distance measurements, and of helical twist as a function of torsion angles delta of dr4 and of G5 unravel the difficulties associated with the study of the G5 sugar pucker conformation equilibrium. Sugar puckers and proton distances are very sensitive criteria to monitor molecular dynamics. Relying on these experimental criteria, we have tested many molecular dynamics preparation phases and we propose a new warm-up and equilibration procedure for molecular dynamics. Thus we show with a 290 ps molecular dynamic run that G5 is in conformational equilibrium and that all NMR data are well reproduced.     

Experimental and theoretical studies of the conformational perturbations induced by an abasic site

The three-dimensional structure of the natural undecamer duplex d(CGCACACACGC). d(GCGTGTGTGCG) has been determined by the combined use of NMR spectroscopy and restrained molecular dynamics (rMD) and also by molecular mechanics calculations using the JUMNA program without experimental distance constraints. Both procedures have also been used to model the abasic structure d(CGCACOCACGC).d(GCGTGTGTGCG), where 'O' indicates a modified abasic site: 3-hydroxy-2-(hydroxymethyl) tetrahydrofuran. For the natural duplex, 134 interproton distances have been obtained by complete relaxation matrix analysis of the NOESY cross-peaks intensities, using MARDIGRAS software. These distances along with 100 torsion angles for sugar ring and additional data derived from canonical A and B-DNA, have been used for structures refinement by restrained molecular dynamics. Comparison of the natural oligomer with the abasic structure obtained earlier by NMR/rMD (Y. Coppel, N. Berthet, C. Coulombeau, Ce. Coulombeau, J. Garcia and J. Lhomme, Biochemistry 36, 4817-4830, 1997) confirms that the creation of an abasic site, in this sequence context, leads to marked helix kinking. It is also shown that the JUMNA procedure is capable of reproducing the overall structural features of the natural and damaged DNA conformations without the use of experimental constraints.     

QM/MM study of energy storage and molecular rearrangements due to the primary event in vision

The energy storage and the molecular rearrangements due to the primary photochemical event in rhodopsin are investigated by using quantum mechanics/molecular mechanics hybrid methods in conjunction with high-resolution structural data of bovine visual rhodopsin. The analysis of the reactant and product molecular structures reveals the energy storage mechanism as determined by the detailed molecular rearrangements of the retinyl chromophore, including rotation of the (C11-C12) dihedral angle from -11 degrees in the 11-cis isomer to -161 degrees in the all-trans product, where the preferential sense of rotation is determined by the steric interactions between Ala-117 and the polyene chain at the C13 position, torsion of the polyene chain due to steric constraints in the binding pocket, and stretching of the salt bridge between the protonated Schiff base and the Glu-113 counterion by reorientation of the polarized bonds that localize the net positive charge at the Schiff-base linkage. The energy storage, computed at the ONIOM electronic-embedding approach (B3LYP/6-31G*:AMBER) level of theory and the S0-->S1 electronic-excitation energies for the dark and product states, obtained at the ONIOM electronic-embedding approach (TD-B3LYP/6-31G*//B3LYP/6-31G*:AMBER) level of theory, are in very good agreement with experimental data. These results are particularly relevant to the development of a first-principles understanding of the structure-function relations in prototypical G-protein-coupled receptors.     

Determination of nucleic acid backbone conformation by 1H NMR.

In DNA or RNA duplexes, the six-bond C3'-O3'-P-O5'-C5'-C4'-C3' backbone linkage connecting adjacent residues contains six torsion angles (epsilon, zeta, alpha, beta, gamma, delta) but only four protons. This seriously limits the ability to define the backbone conformation by NMR using purely 1H-1H distance geometry (DG) methods. The problem is further compounded by the inability to assign two of the four backbone protons, namely the poorly resolved H5' and H5' protons, and invariably leads to DG structures with poorly defined backbone conformations. We have developed and tested a reliable method to constrain the beta, gamma, and epsilon (and indirectly alpha and zeta) backbone torsion angles by lower-bound NOE distances to unassigned H5'/H5' resonances combined with either 1H line widths or the conservative use of sigma J measurements; the method relies only on 1H 2-D NMR data, does not involve any structural assumptions, and leads to much improved backbone convergence among DG structures. The C4'-C5' torsion angle gamma is constrained by lower-bound NOE distances from H2' and from H6/H8 to any H5'/H5', as well as by sigma JH4, coupling measurements in the 3.9-4.4 ppm region; delta is constrained by H1'-H4' NOE distances and by H3'-H4' and H3'-H2' J couplings in COSY data; epsilon is partially constrained by H3' line width and/or further constrained by subtracting the minimum possible sigma JH3'-H from the observed sigma JH3' (COSY) to arrive at the maximum possible JH3'-P, which is then converted to H3'-P distance bounds. The angle beta is partially constrained via H5'-P and H5'-P distance bounds consistent with the maximum H5'-P and H5'-P J couplings derived from the observed H5' and H5' line widths, while alpha and zeta are indirectly constrained by lower distance bounds on the observed (n)H1' to (n + 1)H5'/H5' NOEs combined with the prior partial constraints on beta, gamma, delta, and epsilon. The combined effects of these additional constraints in determining distance geometry structures have been demonstrated using a 12-base duplex, [d(GCCGTTAACGGC)]2. Coordinate RMSDs per atom between structures refined with these constraints from random-embedded DG structures, from ideal A-DNA, and from B-DNA starting structures were less than 0.4 A for the central 8 base pairs indicating good convergence. All backbone angles for the central 8 base pairs are very well constrained with less than 10 degrees variation in any of the 48 torsion angles.     

The unusual conformation of cross-strand disulfide bonds is critical to the stability of β-hairpin peptides

The cross-strand disulfides (CSDs) found in β-hairpin antimicrobial peptides (β-AMPs) show a unique disulfide geometry that is characterized by unusual torsion angles and a short Cα-Cα distance. While the sequence and disulfide bond connectivity of disulfide-rich peptides is well studied, much less is known about the disulfide geometry found in CSDs and their role in the stability of β-AMPs. To address this, we solved the nuclear magnetic resonance (NMR) structure of the β-AMP gomesin (Gm) at 278, 298, and 310 K, examined the disulfide bond geometry of over 800 disulfide-rich peptides, and carried out extensive molecular dynamics (MD) simulation of the peptides Gm and protegrin. The NMR data suggests Cα-Cα distances characteristic for CSDs are independent of temperature. Analysis of disulfide-rich peptides from the Protein Data Bank revealed that right-handed and left-handed rotamers are equally likely in CSDs. The previously reported preference for right-handed rotamers was likely biased by restricting the analysis to peptides and proteins solved using X-ray crystallography. Furthermore, data from MD simulations showed that the short Cα-Cα distance is critical for the stability of these peptides. The unique disulfide geometry of CSDs poses a challenge to biomolecular force fields and to retain the stability of β-hairpin fold over long simulation times, restraints on the torsion angles might be required.     

DNA sequence modulates the conformation of the food mutagen 2-amino-3-methylimidazo[4,5-f]quinoline in the recognition sequence of the NarI restriction enzyme

The conformations of C8-dG adducts of 2-amino-3-methylimidazo[4,5-f]quinoline (IQ) positioned in the C-X1-G, G-X2-C, and C-X3-C contexts in the C-G1-G2-C-G3-C-C recognition sequence of the NarI restriction enzyme were compared, using the oligodeoxynucleotides 5'-d(CTCXGCGCCATC)-3'.5'-d(GATGGCGCCGAG)-3', 5'-d(CTCGXCGCCATC)-3'.5'-d(GATGGCGCCGAG)-3', and 5'-d(CTCGGCXCCATC)-3'.5'-d(GATGGCGCCGAG)-3' (X is the C8-dG adduct of IQ). These were the NarIIQ1, NarIIQ2, and NarIIQ3 duplexes, respectively. In each instance, the glycosyl torsion angle chi for the IQ-modified dG was in the syn conformation. The orientations of the IQ moieties were dependent upon the conformations of torsion angles alpha' [N9-C8-N(IQ)-C2(IQ)] and beta' [C8-N(IQ)-C2(IQ)-N3(IQ)], which were monitored by the patterns of 1H NOEs between the IQ moieties and the DNA in the three sequence contexts. The conformational states of IQ torsion angles alpha' and beta' were predicted from the refined structures of the three adducts obtained from restrained molecular dynamics calculations, utilizing simulated annealing protocols. For the NarIIQ1 and NarIIQ2 duplexes, the alpha' torsion angles were predicted to be -176 +/- 8 degrees and -160 +/- 8 degrees , respectively, whereas for the NarIIQ3 duplex, torsion angle alpha' was predicted to be 159 +/- 7 degrees . Likewise, for the NarIIQ1 and NarIIQ2 duplexes, the beta' torsion angles were predicted to be -152 +/- 8 degrees and -164 +/- 7 degrees , respectively, whereas for the NarIIQ3 duplex, torsion angle beta' was predicted to be -23 +/- 8 degrees . Consequently, the conformations of the IQ adduct in the NarIIQ1 and NarIIQ2 duplexes were similar, with the IQ methyl protons and IQ H4 and H5 protons facing outward in the minor groove, whereas in the NarIIQ3 duplex, the IQ methyl protons and the IQ H4 and H5 protons faced into the DNA duplex, facilitating the base-displaced intercalated orientation of the IQ moiety [Wang, F., Elmquist, C. E., Stover, J. S., Rizzo, C. J., and Stone, M. P. (2006) J. Am. Chem. Soc. 128, 10085-10095]. In contrast, for the NarIIQ1 and NarIIQ2 duplexes, the IQ moiety remained in the minor groove. These sequence-dependent differences suggest that base-displaced intercalation of the IQ adduct is favored when both the 5'- and 3'-flanking nucleotides in the complementary strand are guanines. These conformational differences may correlate with sequence-dependent differences in translesion replication.     

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 3. The magnitude of torsional angle epsilon in d(TpA) from CCOP and HCOP NMR coupling constants.

Carbon-13 NMR spectra of the deoxyribonucleotide d(TpA), 3',5'-cyclic AMP and 3',5'-cyclic dAMP were measured. It is shown that the different substitution of C2' in deoxyribonucleotides versus ribonucleotides does not affect the vicinal C2'-C3'-O3'-P coupling to a measurable extent. Therefore, the same set of Karplus parameters may be used for the C2'-C3-O3'-P couplings in ribonucleotides and in deoxyribonucleotides. Vicinal carbon-phosphorus and proton-phosphorus coupling constants are used to calculate the magnitude of the torsion angle epsilon (C4'-C3'-O3'-P), which amounts to 195(0) in the trans conformer and to 261(0) in the gauche(-) conformer.