Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

Interproton distance bounds from 2D NOE intensities: Effect of experimental noise and peak integration errors

Summary The effect of experimental and integration errors on the calculations in interproton distances from NOE intensities is examined. It is shown that NOE intensity errors can have a large impact on the distances determined. When multiple spin (spin diffusion) effects are significant, the calculated distances are often underestimated, even when using a complete relaxation matrix analysis. In this case, the bias of distances to smaller values is due to the random errors in the NOE intensities. We show here that accurate upper and lower bounds of the distances can be obtained if the intensity errors are properly accounted for in the complete relaxation matrix calculations, specifically the MARDIGRAS algorithm. The basic MARDIGRAS algorithm has been previously described [Borgias, B.A. and James, T.L. (1990) J. Magn. Reson., 87, 475–487]. It has been shown to provide reasonably good interproton distance bounds, but experimental errors can compromise the quality of the resulting restraints, especially for weak cross peaks. In a new approach introduced here, termed RANDMARDI (random error MARDIGRAS), errors due to random noise and integration errors are mimicked by the addition of random numbers from within a specified range to each input intensity. Interproton distances are then calculated for the modified intensity set using MARDIGRAS. The distribution of distances that define the upper and lower distance bounds is obtained by using N randomly modified intensity sets. RANDMARDI has been used in the solution structure determination of the interstrand cross-link (XL) formed between 4-hydroxymethyl-4,5,8-trimethylpsoralen (HMT) and the DNA oligomer d(5-GCGTACGC-3)₂ [Spielmann, H.P. et al. (1995) Biochemistry, 34, 12937–12953]. RANDMARDI generates accurate distance bounds from the experimental NOESY cross-peak intensities for the fixed (known) interproton distances in XL. This provides an independent internal check for the ability of RANDMARDI to accurately fit the experimental data. The XL structure determined using RANDMARDI-generated restrains is in good agreement with other biophysical data that indicate that there is no bend introduced into the DNA by the cross-link. In contrast, isolated spin-pair approximation calculations give distance restraints that, when applied in a restrained molecular dynamics protocol, produce a bent structure.Abbreviations NOE nuclear Overhauser effect - SD standard deviation - HMT 4-hydroxymethyl-4,5,8-trimethylpsoralen - XL psoralen-DNA interstrand cross-link     

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

High-resolution NMR structure of an AT-rich DNA sequence

We have determined, by proton NMR and complete relaxation matrix methods, the high-resolution structure of a DNA oligonucleotide in solution with nine contiguous AT base pairs. The stretch of AT pairs, TAATTATAATTATAATTA, is imbedded in a 27-nucleotide stem-and-loop construct, which is stabilized by terminal GC base pairs and an extraordinarily stable DNA loop GAA (Hirao et al., 1994, Nucleic Acids Res. 22, 576–582). The AT-rich sequence has three repeated TAATTA motifs, one in the reverse orientation. Comparison of the local conformations of the three motifs shows that the sequence context has a minor effect here: atomic RMSD between the three TAATTA fragments is 0.4–0.5 Å, while each fragment is defined within the RMSD of 0.3–0.4 Å. The AT-rich stem also contains a consensus sequence for the Pribnow box, TATAAT. The TpA, ApT, and TpTApA steps have characteristic local conformations, a combination of which determines a unique sequence-dependent pattern of minor groove width variation. All three TpA steps are locally bent in the direction compressing the major groove of DNA. These bends, however, compensate each other, because of their relative position in the sequence, so that the overall helical axis is essentially straight.     

NMR structures of loop B RNAs from the stem-loop IV domain of the enterovirus internal ribosome entry site: a single C to U substitution drastically changes the shape and flexibility of RNA

The 5'-untranslated region of positive-strand RNA viruses harbors many cis-acting RNA structural elements that are important for various viral processes such as replication, translation, and packaging of new virions. Among these is loop B RNA of the stem-loop IV domain within the internal ribosomal entry site (IRES) of enteroviruses, including Poliovirus type 1 (PV1). Studies on PV1 have shown that specific recognition of loop B by the first KH (hnRNP K homology) domain of cellular poly(rC)-binding protein 2 (PCBP2) is essential for efficient translation of the viral mRNA. Here we report the NMR solution structures of two representative sequence variants of enteroviral loop B RNA. The two RNA variants differ at only one position (C vs U) within a six-nucleotide asymmetric internal loop sequence that is the binding site for the PCBP2 KH1 domain. Surprisingly, the two RNAs are drastically different in the overall shape and local dynamics of the bulge region. The RNA with the 5'-AUCCCU bulge sequence adopts an overall L shape. Its bulge nucleotides, especially the last four, are highly flexible and not very well defined by NMR. The RNA with the 5'-AUUCCU bulge sequence adopts an overall U shape, and its bulge sequence exhibits only limited flexibility. A detailed analysis of the two RNA structures and their dynamic properties, as well as available sequence data and known KH domain-RNA complex structures, not only provides insights into how loop B RNA might be recognized by the PCBP2 KH1 domain but also suggests a possible correlation between structural flexibility and pre-existing structural features for protein recognition.     

Visualization of the microwave beam generated by a plasma relativistic microwave amplifier

A method based on the detection of emission of a dielectric screen with metal microinclusions in open air is applied to visualize the transverse structure of a high-power microwave beam. In contrast to other visualization techniques, the results obtained in this work provide qualitative information not only on the electric field strength, but also on the structure of electric field lines in the microwave beam cross section. The interpretation of the results obtained with this method is confirmed by numerical simulations of the structure of electric field lines in the microwave beam cross section by means of the CARAT code.     

Sequence dependent modulating effect of camptothecin on the DNA-cleaving activity of the calf thymus type I topoisomerase.

High-resolution mapping of topol cleavages in the regions of human DNA including the oncogene c-Ha-ras and p53, has revealed three kinds of topol cleavage sites: cleavage sites not affected by camptothecin; cleavage sites reinforced only in the presence of camptothecin, and cleavage sites which weaken in the presence of camptothecin. Statistical analysis of sequences revealed certain nucleotide or dinucleotide preferences for three groups studied. The preferences in camptothecin-reduced sites predominate upstream from the cleavage point, whereas in camptothecin-induced sites the situation is reversed. The influence of camptothecin on cleavage sites induced by two molecular forms of topol has been also studied.     

Solution structure of Syrian hamster prion protein rPrP(90-231)

NMR has been used to refine the structure of Syrian hamster (SHa) prion protein rPrP(90-231), which is commensurate with the infectious protease-resistant core of the scrapie prion protein PrPSc. The structure of rPrP(90-231), refolded to resemble the normal cellular isoform PrPC spectroscopically and immunologically, has been studied using multidimensional NMR; initial results were published [James et al. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 10086-10091]. We now report refinement with better definition revealing important structural and dynamic features which can be related to biological observations pertinent to prion diseases. Structure refinement was based on 2778 unambiguously assigned nuclear Overhauser effect (NOE) connectivities, 297 ambiguous NOE restraints, and 63 scalar coupling constants (3JHNHa). The structure is represented by an ensemble of 25 best-scoring structures from 100 structures calculated using ARIA/X-PLOR and further refined with restrained molecular dynamics using the AMBER 4.1 force field with an explicit shell of water molecules. The rPrP(90-231) structure features a core domain (residues 125-228), with a backbone atomic root-mean-square deviation (RMSD) of 0.67 A, consisting of three alpha-helices (residues 144-154, 172-193, and 200-227) and two short antiparallel beta-strands (residues 129-131 and 161-163). The N-terminus (residues 90-119) is largely unstructured despite some sparse and weak medium-range NOEs implying the existence of bends or turns. The transition region between the core domain and flexible N-terminus, i.e., residues 113-128, consists of hydrophobic residues or glycines and does not adopt any regular secondary structure in aqueous solution. There are about 30 medium- and long-range NOEs within this hydrophobic cluster, so it clearly manifests structure. Multiple discrete conformations are evident, implying the possible existence of one or more metastable states, which may feature in conversion of PrPC to PrPSc. To obtain a more comprehensive picture of rPrP(90-231), dynamics have been studied using amide hydrogen-deuterium exchange and 15N NMR relaxation times (T1 and T2) and 15N{1H} NOE measurements. Comparison of the structure with previous reports suggests sequence-dependent features that may be reflected in a species barrier to prion disease transmission.