Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations.

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2'/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.     

The solution structure of a B-DNA undecamer comprising a portion of the specific target site for the cAMP receptor protein in the gal operon. Refinement on the basis of interproton distance data.

A restrained least squares refinement of the solution structure of the double-stranded DNA undecamer 5'd(AAGTGT-GACAT).5'd(ATGTCACACTT) comprising a portion of the specific target site of the cAMP receptor protein in the gal operon is presented. The structure is refined on the basis of both distance and planarity restraints, 2331 in all. The distance restraints comprise 150 interproton distances determined from pre-steady state nuclear Overhauser enhancement measurements and 2159 other interatomic distances derived from idealized geometry (i.e., distances between covalently bonded atoms, between atoms defining fixed bond angles, and between atoms defining hydrogen bonding in AT and GC base pairs). Two refinements were carried out and in both cases the final RMS difference between the experimental and calculated interproton distances was 0.2 A. The difference between the two refined structures is small (overall RMS difference of 0.23 A) and represents the error in the refined coordinates. Although the refined structures have an overall B-type conformation there are large variations in many of the local conformational parameters including backbone and glycosidic bond torsion angles, helical twist and propellor twist, base roll and base tilt angles.     

Relationship of DNA structure to internal dynamics: correlation of helical parameters from NOE-based NMR solution structures of d(GCGTACGC)(2) and d(CGCTAGCG)(2) with (13)C order parameters implies conformational coupling in dinucleotide units

The coupling between the conformational properties of double-stranded DNA and its internal dynamics has been examined. The solution structures of the isomeric DNA oligomers d(GCGTACGC)(2) (UM) and d(CGCTAGCG)(2) (CTSYM) were determined with (1)H NMR spectroscopy by utilizing distance restraints from total relaxation matrix analysis of NOESY cross-peak intensities in restrained molecular dynamics calculations. The root-mean-square deviation of the coordinates for the ensemble of structures was 0.13 A for UM and 0.49 A for CTSYM, with crystallographic equivalent R(c)=0.41 and 0.39 and sixth-root residual R(x)=0.11 and 0.10 for UM and CTSYM, respectively. Both UM and CTSYM are B-form with straight helical axes and show sequence-dependent variations in conformation. The internal dynamics of UM and CTSYM were previously determined by analysis of (13)C relaxation parameters in the context of the Lipari & Szabo model-free formalism. Helical parameters for the two DNA oligomers were examined for linear correlations with the order parameters (S(2)) of groups of (13)C spins in base-pairs and dinucleotide units of UM and CTSYM. Correlations were found for six interstrand base-pair parameters tip, y-displacement, inclination, buckle and stretch with various combinations of S(2) for atoms in Watson-Crick base-pairs and for two inter-base-pair parameters, rise and roll with various combinations of S(2) for atoms in dinucleotides. The correlations for the interstrand base-pair helical parameters indicate that the conformations of the deoxyribose residues of each strand are dynamically coupled. Also, the inter-base-pair separation has a profound effect on the local internal motions available to the DNA, supporting the idea that rise is a principal degree of freedom for DNA conformational variability. The correlations indicate collective atomic motions of spins that may represent specific motional modes in DNA, and that base sequence has a predictable effect on the relative order of groups of spins both in the bases and in the deoxyribose ring of the DNA backbone. These observations suggest that an important functional outcome of DNA base sequence is the modulation of both the conformation and dynamic behavior of the DNA backbone.     

Determining local conformational variations in DNA. Nuclear magnetic resonance structures of the DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC) generated using back-calculation of the nuclear Overhauser effect spectra, a distance geometry algorithm and constrained molecular dynamics

Two-dimensional nuclear magnetic resonance (n.m.r.) spectroscopy and a variety of computational techniques have been used to generate three-dimensional structures of the two DNA duplexes d(CGCCTAATCG) and d(CGTCACGCGC). The central six base-pairs in these two decamers contain all ten dinucleotide pairs in DNA and thus, represent a model system for investigating how the local structure of DNA varies with base sequence. Resonance assignments were made for the non-exchangeable base protons and most of the C-1'-C-4' sugar protons in both decamers. Three-dimensional structures were generated using a distance geometry algorithm and these initial structures were refined by optimizing the fit of back-calculated spectra against the experimental two-dimensional nuclear Overhauser effect (NOE) spectra. This back-calculation procedure consists of calculating NOE cross relaxation rates for a given structure by solution of the Bloch equations, and directly accounts for spin diffusion effects. Use of this refinement procedure eliminates some assumptions that have been invoked when generating structures of DNA oligomers from n.m.r. data. Constrained energy minimization and constrained quenched molecular dynamics calculation were also performed on both decamers to help generate energetically favorable structures consistent with the experimental data. Analysis of the local conformational parameters of helical twist, helical rise, propeller twist, displacement and the alpha, beta, gamma, epison and zeta backbone torsion angles in these structures shows that these parameters span a large range of values relative to the X-ray data of nucleic acids. However, the glycosidic and pseudorotation angles are quite well defined in these structures. The implications that these results have for determination of local structural variations of DNA in solution, such as those predicted by Callidine's rules, are discussed. Our results differ significantly from some previous studies on determining local conformations of nucleic acids and comparisons with these studies are made.     

DNA Modeller: an interactive program for modelling stacks of DNA base pairs on a microcomputer

DNA Modeller is a microcomputer program for interactively manipulatingup to 20 bp in a DNA double helical arrangement. It calculatesthe van der Waals and electrostatic energies of base-base interactionsusing the AMBER potential, minimizes the energy with respectto the pair (buckle, propeller, opening, shear, stretch, stagger)and step (tilt, roll, twist, shift, slide, rise) parameters,calculates lengths of the canonical hydrogen bonds between thecomplementary bases, and calculates interatomic distances betweenthe successive base pairs. Input/output files are simple listsof the step and pair parameters or lists of the atom specifications(N1, C2, etc.) and their Cartesian coordinates (compatible withthe Desktop Molecular Modeller ^*.mol files). The program issupplied with a readbrk utility which transforms PDB/NDB tothe^*.mol format readable by DNA Modeller. The DNA crystal structuresdeposited in the PDB or NDB databases can thus be analyzed,and their bases visualized and interactively manipulated. Inaddition, DNA Modeller can calculate the base pair and stepgeometrical parameters and interaction energies. A plotter utilitycreates wire mono or stereo pictures of the bases. This programis designed for IBM-compatible computers working under DOS orcan run as a DOS application under MS Windows 3.x or Merge (SCOUnix DOS emulator).     

Rigid spin-labeled nucleoside C: a nonperturbing EPR probe of nucleic acid conformation

Rigid spin-labeled nucleoside C, an analog of deoxycytidine that base-pairs with deoxyguanosine, was incorporated into DNA oligomers by chemical synthesis. Thermal denaturation experiments and circular dichroism (CD) measurements showed that C has a negligible effect on DNA duplex stability and conformation. Nucleoside C was incorporated into several positions within single-stranded DNA oligomers that can adopt two hairpin conformations of similar energy, each of which contains a four-base loop. The relative mobility of nucleotides in the alternating C/G hairpin loops, 5'-d(GCGC) and 5'-d(CGCG), was determined by electron paramagnetic resonance (EPR) spectroscopy. The most mobile nucleotide in the loop is the second one from the 5'-end, followed by the third, first and fourth nucleotides, consistent with previous NMR studies of DNA hairpin loops of different sequences. The EPR hairpin data were also corroborated by fluorescence spectroscopy using oligomers containing reduced C (C(f)), which is fluorescent. Furthermore, EPR spectra of duplex DNAs that contained C at the end of the helix showed features that indicated dipolar coupling between two spins. These data are consistent with end-to-end duplex stacking in solution, which was only observed when G was paired to C, but not when C was paired with A, C or T.     

Base sequence and helix structure variation in B and A DNA

The observed propeller twist in base-pairs of crystalline double-helical DNA oligomers improves the stacking overlap along each individual helix strand. But, as proposed by Calladine, it also leads to clash or steric hindrance between purines at adjacent base-pairs on opposite strands of the helix. This clash can be relieved by: (1) decreasing the local helix twist angle between base-pairs; (2) opening up the roll angle between base-pairs on the side on which the clash occurs; (3) separating purines by sliding base-pairs along their long axes so that the purines are partially pulled out of the stack (leading to equal but opposite alterations in main-chain torsion angle delta at the two ends of the base-pair); and (4) flattening the propeller twist of the offending base-pairs. Simple sum functions, sigma 1 through sigma 4, are defined, by which the expected local variation in helix twist, base roll angle, torsion angle delta and propeller twist may be calculated from base sequence. All four functions are quite successful in predicting the behavior of B DNA. Only the helix twist and base roll functions are applicable to A DNA, and the helix twist function begins to fail for an A helical RNA/DNA hybrid. Within these limits, the sequence-derived sum functions match the observed helix parameter variation quite closely, with correlation coefficients greater than 0.900 in nearly all cases. Implications of this sequence-derived helix parameter variation for repressor-operator interactions are considered.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Helical coherence of DNA in crystals and solution

The twist, rise, slide, shift, tilt and roll between adjoining base pairs in DNA depend on the identity of the bases. The resulting dependence of the double helix conformation on the nucleotide sequence is important for DNA recognition by proteins, packaging and maintenance of genetic material, and other interactions involving DNA. This dependence, however, is obscured by poorly understood variations in the stacking geometry of the same adjoining base pairs within different sequence contexts. In this article, we approach the problem of sequence-dependent DNA conformation by statistical analysis of X-ray and NMR structures of DNA oligomers. We evaluate the corresponding helical coherence length—a cumulative parameter quantifying sequence-dependent deviations from the ideal double helix geometry. We find, e.g. that the solution structure of synthetic oligomers is characterized by 100–200 Å coherence length, which is similar to ~150 Å coherence length of natural, salmon-sperm DNA. Packing of oligomers in crystals dramatically alters their helical coherence. The coherence length increases to 800–1200 Å, consistent with its theoretically predicted role in interactions between DNA at close separations.     

Structural studies on ligand–DNA systems: A robust approach in drug design

Molecular docking, molecular mechanics, molecular dynamics and relaxation matrix simulation protocols have been extensively used to generate the structural details of ligand-receptor complexes in order to understand the binding interactions between the two entities. Experimental methods like NMR spectroscopy and X-ray crystallography are known to provide structural information about ligand-receptor complexes. In addition, fluorescence spectroscopy, circular dichroism (CD) spectroscopy and molecular docking have also been utilized to decode the phenomenon of the ligand-DNA interactions, with good correlation between experimental and computational results. The DNA binding affinity was demonstrated by analysing fluorescence spectral data. Structural rigidity of DNA upon ligand binding was identified by CD spectroscopy. Docking is carried out using the DNA-Dock program which results in the binding affinity data along with structural information like interatomic distances and H-bonding, etc. The complete structural analyses of various drug-DNA complexes have afforded results that indicate a specific DNA binding pattern of these ligands. It also exhibited that certain structural features of ligands can make a ligand to be AT- or GC-specific. It was also demonstrated that changing specificity from AT base pairs to GC base pairs further improved the DNA topoisomerase inhibiting activity in certain ligands. Thus, a specific molecular recognition signature encrypted in the structure of ligand can be decoded and can be effectively employed in designing more potent antiviral and antitumour agents.     

The use of BrCN for assembling modified DNA duplexes and DNA-RNA hybrids; comparison with water-soluble carbodiimide.

Both cyanogen bromide (BrCN) and 1-ethyl-3-(3'-dimethylaminopropyl) carbodiimide may be used as coupling reagents for the template-directed assembly of DNA duplexes containing the sugar-phosphate backbone modification. Both reagents show similar ligation site structure-specific trend. Practical recommendations are given for selection of the condensing reagent depending on the properties of the duplex. Based on 31P NMR spectroscopy data, a scheme is suggested for BrCN activation of the nucleotide phosphomonoester group. Using both condensing reagents, we studied the condensation of oligonucleotides containing ribo-segments (from mononucleotide residue to full sequence) on the DNA template. Efficiency of the chemical ligation of RNA oligomers was shown to be much lower than that of DNA analogues. The coupling yield depends on the position of the RNA segment in the hybrid duplexes and on the position of the phosphate group in the nick.     

Solution conformation on bovine growth hormone releasing factor by 1H NMR and molecular modeling.

The structure of bovine growth hormone releasing factor (bGHRF) consisting of 44 amino acids has been studied in CD and 1H nuclear magnetic resonance (NMR) spectroscopy in conjunction with molecular modeling. Since bGHRF does not have an ordered structure in water alone, a 30% 2,2,2-trifluoroethanol (TFE) aqueous solvent was used to induce considerable alpha-helical structures, which corresponds to a helical content of approximately 62% as determined by circular dichroism (CD). The secondary structure was obtained from nuclear Overhauser enhancement and 3J(HN alpha) coupling constant in 30% TFE solution. Three-dimensional structures consistent with NMR data were generated by using distance geometry calculation. A set of 267 interproton distances derived from nuclear Overhauser effect correlation spectroscopy (NOESY) experiments and coupling constants were used. From the initial random conformations, 50 distance geometry structures with minimal violations were selected for further refinement. The 14 best structures were obtained after simulated annealing calculation with energy minimization. The structure of bGHRF in 30% TFE solution was characterized by one alpha-helix (residues 8-19), two poorly constrained helices (residues 23-27 and residues 31-34) and a beta I(III)-turn fragment (residues 20-23; phi(i+1) = -53.1 degrees, psi(i+1) = -19.6 degrees, phi(i+2) = -59.9 degrees, psi(i+2) = -20.6 degrees) connected by the segments of less defined structures in N-terminal and omega-shaped flexible C-terminal determined from NOESY cross peaks between helical segment (residues 14-18) and tail fragment (residues 42-44). The obtained structure will play an important role toward the understanding of the structural and functional role of the GHRF.     

Upper limit for the variances of some helical parameters in DNA double helix

Assuming that the realistic DNA chains are random sequences of purines and pyrimidines and, by using the Dickerson sum functions, it is shown that there exists an upper limit for the variance of helical twist angles, the base-plane roll angles and the other helical parameters respectively in the DNA sequences. The estimates of the variances of all the above helical parameters for the DNA sequences in the Los Alamos data base have been performed and found to be in good agreement with the theoretical results obtained in this paper.     

Molecular dynamics simulations of the d(CCAACGTTGG)(2) decamer: influence of the crystal environment

Molecular dynamics (MD) simulations of the DNA duplex d(CCAACGTTGG)(2) were used to study the relationship between DNA sequence and structure. Two crystal simulations were carried out; one consisted of one unit cell containing two duplexes, and the other of two unit cells containing four duplexes. Two solution simulations were also carried out, one starting from canonical B-DNA and the other starting from the crystal structure. For many helicoidal parameters, the results from the crystal and solution simulations were essentially identical. However, for other parameters, in particular, alpha, gamma, delta, (epsilon - zeta), phase, and helical twist, differences between crystal and solution simulations were apparent. Notably, during crystal simulations, values of helical twist remained comparable to those in the crystal structure, to include the sequence-dependent differences among base steps, in which values ranged from 20 degrees to 50 degrees per base step. However, in the solution simulations, not only did the average values of helical twist decrease to approximately 30 degrees per base step, but every base step was approximately 30 degrees, suggesting that the sequence-dependent information may be lost. This study reveals that MD simulations of the crystal environment complement solution simulations in validating the applicability of MD to the analysis of DNA structure.     

Conformational variation of the central CG site in d(ATGACGTCAT)2 and d(GAAAACGTTTTC)2. An NMR, molecular modelling and 3D-homology investigation.

The determination of the solution structure of two self-complementary oligomers d(ATGACGTCAT)2 (CG10) and d(GAAAACGTTTTC)2 (CG12), both containing the 5'-pur-ACGT-pyr-3' sequence, is reported. The impact of the base context on the conformation of the central CpG site has been examined by a combined approach of: (a) 2D 1H-NMR and 31P-NMR; (b) molecular mechanics under experimental constraints; (c) back-calculations of NOESY spectra and iterative refinements of distances; and (d) 3D-homology search of the central tetrad ACGT within the complete oligonucleotides. A full NMR study of each fragment is achieved by means of standard 2D experiments: NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Sugar phase angle, epsilon-zeta difference angle and NOE-derived distances are input as experimental constraints to generate molecular models by energy minimization with the help of jumna. The morass program is used to iteratively refine the structures obtained. The similarity of the two ACGTs within the whole oligonucleotides is investigated. Both the decamer and the dodecamer adopt a B-like DNA conformation. However, the helical parameters within this conformational type are significantly different in CG12 and CG10. The central CpG step conformation is not locked by its nearest environment (5'A and 3'T) as seen from the structural analysis of ACGT in the two molecules. In CG12, despite the presence of runs of A-T pairs, CpG presents a high twist of 43 degrees and a sugar phase at the guanine of about 180 degrees, previously observed in other ACGT-containing-oligomers. Conversely, ACGT in CG10 exhibits strong inclinations, positive rolls, a flat profile of sugar phase, twist and glycosidic angles, as a result of the nucleotide sequence extending beyond the tetrad. The structural specificity of CG10 and its flexibility (as reflected by its energy) are tentatively related to the process of recognition of the cyclic AMP response element by its cognate protein.     

A simple spectral-driven procedure for the refinement of DNA structures by NMR spectroscopy.

We have developed a simple and quantitative procedure (SPEDREF) for the refinement of DNA structures using experimental two-dimensional nuclear Overhauser effect (2D NOE) data. The procedure calculates the simulated 2D NOE spectrum using the full matrix relaxation method on the basis of a molecular model. The volume of all NOE peaks is measured and compared between the experimental and the calculated spectra. The difference of the experimental and simulated volumes is minimized by a conjugated gradient procedure to adjust the interproton distances in the model. An agreement factor (analogous to the crystallographic R-factor) is used to monitor the progress of the refinement. The procedure is an The agreement is considered to be complete when several parameters, including the R-factor, the energy associated with the molecule, the local conformation (as judged by the sugar pseudorotation), and the global conformation (as judged by the helical x-displacement), are refined to their respective convergence. With the B-DNA structure of d(CGATCG) as an example, we show that DNA structure may be refined to produce calculated NOE spectra that are in excellent agreement with the experimental 2D NOE spectra. This is judged to be effective by the low R-factor of approximately 15%. Moreover, we demonstrate that not only are NOE data very powerful in providing details of the local structure but, with appropriate weighting of the NOE constraints, the global structure of the DNA double helix can also be determined, even when starting with a grossly different model. The reliability and limitations of a DNA structure as determined by NMR spectroscopy are discussed.     

On the structure of poly d(A-S4T).

The helical twist of poly d(A-s4T) was determined from the periodicity of the cleavage patterns of the double stranded polydeoxynucleotide adsorbed on calcium phosphate and found to be 14 bp per turn. Both cleavage patterns and 31P NMR spectra indicate a mononucleotide structure rather than an alternating B DNA like poly d(A-T). The failure of nucleosome formation excludes a B type structure. The discrepancy of the mononucleotide structure found in 31P NMR spectra and the dinucleotide structure given by X ray fiber diffraction is explained by an alternating tilt of the planes of the base pairs (base roll) as a consequence of a strong propeller twist. The importance of interstrand stacking interactions of adjacent 4-thiothymidines for the helical stability is discussed.     

Solid-state NMR and molecular dynamics simulations reveal the oligomeric ion-channels of TM2-GABAA stabilized by intermolecular hydrogen bonding

The second transmembrane (TM2) domain of GABA_A receptor forms the inner-lining surface of chloride ion-channel and plays important roles in the function of the receptor protein. In this study, we report the first structure of TM2 in lipid bilayers determined using solid-state NMR and MD simulations. The interatomic ¹³C-¹⁵N distances measured from REDOR magic angle spinning experiments on multilamellar vesicles, containing a TM2 peptide site specifically labeled with ¹³C′ and ¹⁵N isotopes, were used to determine the secondary structure of the peptide. The ¹⁵N chemical shift and ¹H-¹⁵N dipolar coupling parameters measured from PISEMA experiments on mechanically aligned phospholipid bilayers, containing a TM2 peptide site specifically labeled with ¹⁵N isotopes, under static conditions were used to determine the membrane orientation of the peptide. Our results reveal that the TM2 peptide forms an alpha helical conformation with a tilted transmembrane orientation, which is unstable as a monomer but stable as pentameric oligomers as indicated by MD simulations. Even though the peptide consists of a number of hydrophilic residues, the transmembrane folding of the peptide is stabilized by intermolecular hydrogen bondings between the side chains of Ser and Thr residues as revealed by MD simulations. The results also suggest that peptide-peptide interactions in the tilted transmembrane orientation overcome the hydrophobic mismatch between the peptide and bilayer thickness.     

DNA structures from phosphate chemical shifts

For B-DNA, the strong linear correlation observed by nuclear magnetic resonance (NMR) between the ³¹P chemical shifts (δP) and three recurrent internucleotide distances demonstrates the tight coupling between phosphate motions and helicoidal parameters. It allows to translate δP into distance restraints directly exploitable in structural refinement. It even provides a new method for refining DNA oligomers with restraints exclusively inferred from δP. Combined with molecular dynamics in explicit solvent, these restraints lead to a structural and dynamical view of the DNA as detailed as that obtained with conventional and more extensive restraints. Tests with the Jun-Fos oligomer show that this δP-based strategy can provide a simple and straightforward method to capture DNA properties in solution, from routine NMR experiments on unlabeled samples.