Molecular dynamics of an in vacuo model of duplex d(CGCGAATTCGCG) in the B-form based on the amber 3.0 force field.

The characteristics of 100 ps of molecular dynamics (MD) on the DNA dodecamer d(CGCGAATTCGCG) at 300 K are described and investigated. The simulation is based on an in vacuo model of the oligomer and the AMBER 3.0 force field configured in the manner of Singh, U. C., S. J. Weiner, and P. A. Kollman, (1985, Proc. Natl. Acad. Sci. USA. 82:755-759). The analysis of the results was carried out using the "curves, dials, and windows" procedure (Ravishanker, G., S. Swaminathan, D. L. Beveridge, R. Lavery, and H. Sklenar. 1989. J. Biomol. Struct. Dyn. 6:669-699). The results indicate this dynamical model to be a provisionally stable double helix which lies at approximately 3.2 A rms deviation from the canonical B-form. There is, however, a persistent nonplanarity in the base pair orientations which resemble that observed in canonical A-DNA. The major groove width is seen to narrow during the course of the simulation and the minor groove expands, contravariant to the alterations in groove width seen in the crystal structure of the native dodecamer (Drew, H. R., R. M. Wing, T. Takano, C. Broka, S. Tanaka, I. Itakura, and R. E. Dickerson, 1981. Proc. Natl. Acad. Sci. USA. 78:2179-2183). The propeller twist in the bases, the sequence dependence of the base pair roll and aspects of bending in the helix axis are in some degree of agreement with the crystal structure. The patterns in DNA bending are observed to follow Zhurkin theory (Zhurkin, V. B. 1985. J. Biomol. Struct. Dyn. 2:785-804.). The relationship between the dynamical model and structure in solution is discussed.     

Back to Units of Protein Folding

Abstract

In response to the criticism by A. Finkelstein (J. Biomol. Struct. Dyn. 20, 311–314, 2002) of our Communication (J. Biomol. Struct. Dyn. 20, 5–6, 2002) several issues are dealt with. Importance of the notion of elementary folding unit, its size and structure, and the necessity of further characterization of the units for the elucidation of the protein folding in vivo are discussed. The criticism (J. Biomol. Struct. Dyn. 20, 311–314, 2002) on the hierarchical protein folding is also briefly addressed.     

Sequence-Dependent Bending of DNA and Phasing of Nucleosomes

Abstract

Conformational analysis has revealed anisotropic flexibility of the B-DNA double helix: it bends most easily into the grooves, being the most rigid when bent in a perpendicular direction. This result implies that DNA in a nucleosome is curved by means of relatively sharp bends (“mini-kinks”) which are directed into the major and minor grooves alternatively and separated by 5–6 base pairs. The “mini-kink” model proved to be in keeping with the x-ray structure of the B-DNA dodecamer resolved later, which exhibits two “annealed kinks”, also directed into the grooves.

The anisotropy of B DNA is sequence-dependent: the pyrimidine-purine dimers (YR) favor bending into the minor groove, and the purine-pyrimidine dinucleotides (RY), into the minor one. The RR and YY dimers appear to be the most rigid dinucleotides. Thus, a DNA fragment consisting of the interchanging oligopurine and oligopyrmidine blocks 5–6 base pairs long should manifest a spectacular curvature in solution.

Similarly, a nucleotide sequence containing the RY and YR dimers separated by a half-pitch of the double helix is the most suitable for wrapping around the nucleosomal core. Analysis of the numerous examples demonstrating the specific alignment of nucleosomes on DNA confirms this concept. So, the sequence-dependent “mechanical” properties of the double helix influence the spatial arrangement of DNA in chromatin.     

A and Z canonical conformations in d(CnGCGn) crystals characterized by microFTIR and microRaman spectroscopies

Two crystals d(C₂GCG₂) and d(C₅GCG₅) have been studied under microscope by Fourier transform ir spectroscopy and Raman spectroscopy. The x-ray diffraction study of the latter crystal had shown that the d(C₅GCG₅) sequence is the first DNA dodecamer known to adopt a canonical A conformation [N. Verdaguer, J. Aymami, D. Fernandez-Forner, I. Fita, M. Coll, T. Huynh-Dinh, J. Igolen, and J. A. Subirana (1991) Journal of Molecular Biology, Vol. 221, pp. 623–635]. Characteristic ir marker bands and Raman marker peaks of the A conformation have thus been obtained and are compared with previously proposed assignments correlated to fiber diffraction x-ray results obtained on polymers. The d(C₂GCG₂) sequence crystal had previously been studied in an intermediate form between B and Z [L. Urpi, J. P. Ridoux, J. Liquier, N. Verdagner, I. Fita, J. A. Subirana, F. Iglesias, T. Huynh-Dinh, J. Igolen, and E. Taillandier (1989) Nucleic Acids Research, Vol. 17, pp. 6669–6679]. In this paper we present results obtained from a crystal with this oligonucleotide in Z conformation. The effect of the crystallization conditions on the geometry of the obtained oligomer helix is discussed. The influence of the addition, to the central tetramer CGCG, of dC_n stretches (at the 5′ end) and dG_n stretches (at the 3′ end) of different lengths, on the conformational flexibility of the nucleic acid, is considered. © 1993 John Wiley & Sons, Inc.     

Dependencies of Substitution Steps Number on Hamming Distance Are Identical For One-parameter Discrete Models of Both Direct and Parallel Genetic Diversity

Abstract

With the help of previously introduced enumeration procedure (M.Yu. Shchelkanov, A.N. Yudin, A.V Antonov, N.S. Starikov, A.A. Vedenov, E.V. Karamov, J. Biomol. Struct. Dyn. 15, 217–229 (1997)) and probability distribution function for the enumeration after some substitution steps (M.Yu. Shchelkanov, L.A. Soinov, V.V. Zalunin, D.A. Gumennyi, A.N. Yudin, A.A. Natan, V.B. Kireev, E.V. Karamov, J. Biomol. Struct. Dyn. 15, N 4, (1998)) we have demonstrated that dependencies of replication acts number on Hamming distance are identical for one-parameter discrete models of both direct and parallel genetic diversity.     

Back to units of protein folding

In response to the criticism by A. Finkelstein (J Biomol Struct Dyn 20, 311-314, 2002) of our Communication (J Biomol Struct Dyn 20, 5-6, 2002) several issues are dealt with. Importance of the notion of elementary folding unit, its size and structure, and the necessity of further characterization of the units for the elucidation of the protein folding in vivo are discussed. The criticism (J Biomol Struct Dyn 20, 311-314, 2002) on the hierarchical protein folding is also briefly addressed.     

The effect of deuterium oxide on the stability of the collagen model peptides H‐(Pro‐Pro‐Gly)10‐OH,H‐(Gly‐Pro‐4(R)Hyp)9‐OH,and Type I collagen

The collagen triple helix has a larger accessible surface area per molecular mass than globular proteins, and therefore potentially more water interaction sites. The effect of deuterium oxide on the stability of collagen model peptides and Type I collagen molecules was analyzed by circular dichroism and differential scanning calorimetry. The transition temperatures (T_m) of the protonated peptide (Pro‐Pro‐Gly)₁₀ were 25.4 and 28.7°C in H₂O and D₂O, respectively. The increase of the T_m of (Pro‐Pro‐Gly)₁₀ measured calorimetrically at 1.0°C min^?1 in a low pH solution from the protonated to the deuterated solvent was 5.1°C. The increases of the T_m for (Gly‐Pro‐4(R)Hyp)₉ and pepsin‐extracted Type I collagen were measured as 4.2 and 2.2°C, respectively. These results indicated that the increase in the T_m in the presence of D₂O is comparable to that of globular proteins, and much less than reported previously for collagen model peptides [Gough and Bhatnagar, J Biomol Struct Dyn 1999, 17, 481–491]. These experimental results suggest that the interaction of water molecules with collagen is similar to the interaction of water with globular proteins, when the ratio of collagen to water is very small and collagen is monomerically dispersed in the solvent. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 93–101, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Defining the Structure of Irregular Nucleic Acids: Conventions and Principles

Abstract

The algorithm “Curves”, that we have recently presented in this journal (J. Biolmol. Str. Dynam. 6, 63–91 (1988)), is updated to take into account the conventions developed at the Cambridge meeting on DNA curvature (September 1988) and extended to the calculation of local parameters. In addition, the principles which govern the choices made in establishing the Curves algorithm are compared with the approaches adopted by other authors.     

A molecular model for RecA-promoted strand exchange via parallel triple-stranded helices.

A number of studies have concluded that strand exchange between a RecA-complexed DNA single strand and a homologous DNA duplex occurs via a single-strand invasion of the minor groove of the duplex. Using molecular modeling, we have previously demonstrated the possibility of forming a parallel triple helix in which the single strand interacts with the intact duplex in the minor groove, via novel base interactions (Bertucat et al., J. Biomol. Struct. Dynam. 16:535-546). This triplex is stabilized by the stretching and unwinding imposed by RecA. In the present study, we show that the bases within this triplex are appropriately placed to undergo strand exchange. Strand exchange is found to be exothermic and to result in a triple helix in which the new single strand occupies the major groove. This structure, which can be equated to so-called R-form DNA, can be further stabilized by compression and rewinding. We are consequently able to propose a detailed, atomic-scale model of RecA-promoted strand exchange. This model, which is supported by a variety of experimental data, suggests that the role of RecA is principally to prepare the single strand for its future interactions, to guide a minor groove attack on duplex DNA, and to stabilize the resulting, stretched triplex, which intrinsically favors strand exchange. We also discuss how this mechanism can incorporate homologous recognition.     

Influence of the incorporation of a cyclohexenyl nucleic acid (CeNA) residue onto the sequence d(CGCGAATTCGCG)

Cyclohexene nucleic acids (CeNA), which are characterized by the presence of a cyclohexene moiety instead of a natural (deoxy)ribose sugar, are known to increase the thermal and enzymatic stability when incorporated in RNA oligonucleotides. As it has been demonstrated that even a single cyclohexenyl nucleoside, when incorporated in an oligonucleotide, can have a profound effect on the biological activity of the oligonucleotide, further research is warranted to study the complex of such oligonucleotides with target proteins. In order to analyse the influence of CeNA residues onto the helix conformation and hydration of natural nucleic acid structures, a cyclohexenyl-adenine building block (xAr) was incorporated into the Dickerson sequence CGCGA(xAr)TTCGCG. The crystal structure of this sequence determined to a resolution of 1.90 Å. The global helix belongs to the B-type family and shows a water spine, which is partially broken up by the apolar cyclohexene residue. The cyclohexene ring adopts the ²E-conformation allowing a better incorporation of the residue in the dodecamer sequence. The crystal packing is stabilized by cobalt hexamine residues and belongs to space group P222₁, never before reported for nucleic acids.     

Development of the multiple sequence approximation within the AGADIR model of α-helix formation: Comparison with Zimm-Bragg and Lifson-Roig formalisms

In this work we present the development of the multiple sequence approximation (AGADIRms) and the standard one-sequence approximation (AGADIRIs) within the framework of AGADIR's α-helix formation model. The extensive comparison between these new formulations and the original one [AGADIR; v. Muñoz and L. Serrano (1994),Nat. Struct. Biol., Vol. 1, pp. 399–409] indicates that the standard one-sequence approximation is virtually identical to the multiple sequence approximation, while the previously used residue partition function approximation [Muñoz and Serrano (1994); (1995), J. Mol. Biol., Vol. 245, pp. 275–296] is less precise. The calculations of the average helical content performed with AGADIR are precise for peptides of less than 30 residues and progressively diverge from the multiple sequence formulation for longer peptides. The helicity distribution of heteropolypeptides with less than 50% average helical content is also well described, while those of quasi-homopolymers with high helical content tend to be flattened. These inaccuracies lead to an underestimation of 0.017 kcal/mol for the mean-residue enthalpic contribution in AGADIR, as compared to AGADIRms and AGADIRIs. The other energy contributions to α-helix stability are not affected by the original statistical approximation. We also discuss the particularities of the model for α-helix formation utilized in AGADIR and compare it with the classical Zimm-Bragg and Lifson-Roig theories. Moreover, we develop the mathematical relationships between the basic AGADIR energy contributions and helix nucleation and elongation, which permit the quantitative comparison between formalisms. Remarkably, the comparison between AGADIRms and the Lifson-Roig formalism shows that, despite the differences on treating helix/coil cooperativity, both theories give virtually identical results when an equivalent set of parameters is used. This indicates that the helix/coil transition is a solid theory independent of the particularities of the model for α-helix formation. © 1997 John Wiley & Sons, Inc. Biopoly 41: 495–509, 1997.     

Multistep Brownian dynamics: application to short wormlike chains

Brownian dynamics simulations of short wormlike chains are carried out using the method of Ermak and McCammon [(1978) J. Chem. Phys. 69 , 1352–1360]. Following Hagerman and Zimm [(1981) Biopolymers 20 , 1481–1502], the wormlike chain is modeled as a string of beads. In each simulation, the dynamic evolution of an ensemble of 100 randomly generated chains is calculated for a period of from 3 to 200 ns. Two different “experiments,” fluorescence depolarization and dynamic light scattering, were performed in these simulations. Since we are primarily interested in the bending motions and not the torsional motions in this work, we have placed the transition moments along the local symmetry axis of the wormlike chain in the fluorescence depolarization “experiment.” As predicted by the Barkley and Zimm theory [(1979) J. Chem. Phys. 70 , 2991–3008], a considerable amount of rapid bending motion was detected by fluorescence depolarization, though not as much as predicted by theory. We conclude that these differences are primarily due to differences between the model used in the theory and the simulations. The light-scattering experiment was found to be insensitive to internal motion in the low scattering angle limit.     

Symmetry and structure of RNA and DNA triple helices

Despite wide interest in nucleic acid triple helices, there has beenno stereochemically satisfactory structure of an RNA triple helixin atomic detail. An RNA triplex structure has previously been proposed based on fiber diffraction and molecular modeling [S. Arnott and P. J. Bond (1973) Nature New Biology, Vol. 244. pp. 99–101; S. Arnott. P. J. Bond. E. Seising, and P. J. C. Smith (1976) Nucleic Acids Research, Vol. 3. pp.2459–2470], but it has nonallowed close contacts at every triplet and is therefore not stereochemically acceptable. We propose here a new modelfor an RNA triple helix in which the three chains have identical backbone conformations and are symmetry related. There are no short contacts. The modeling employs a novel geometrical approach using the linked atom least squares [P. J. C. Smith and S. Arnott (1978) Acta Crystallographica, Vol. A34, pp. 3–11] program and is not based on energy minimization. In general, the method leads to a range of possible structures rather than a unique structure. In the present case, however, the constraints resulting from theintroduction of a third strand limit the possible structures to a very small range of conformation space. This method was used previously to obtain a model for DNA triple helices [G. Raghunathan, H. T. Miles, and V. Sasisekharan (1993) Biochemistry, Vol. 32, pp. 455–462], subsequently confirmed by fiber-type x-ray diffraction of oligomeric crystals [K. Liu. H. T. Miles. K. D. Parris, and V. Sasisekharan (1994) Nature Structural Biology, Vol. 1. pp. 11–12]. The above triple helices have Watson–Crick–Hoogsteen [K. Hoogsteen (1963) Acta Crystallographica, Vol. 16. pp. 907–916] pairing of the three bases. The same modeling method was used to investigate the feasibility of three-dimensional structures based on the three possible alternative hydrogen-bonding schemes: Watson–Crick–reverse Hoogsteen, Donogue [J. Donohue (1953) Proceeding of the national Academy of Science USA, Vol. 39, pp. 470–475] (reverse Watson–Crick)–Hoogsteen, and Donohue–reverse Hoogsteen. We found that none of these can occur in either RNA or DNA helices because they give rise only to structures with prohibitively short contacts between backbone and base atoms in the same chain. © 1995 John Wiley & Sons, Inc.     

Nuclear Overhauser Data and Stereochemical Considerations Suggest that Netropsin Binds Symmetrically Within the Minor Groove of Poly(dA)·Poly(dT), Forming Hydrogen Bonds with Both Strands of the Double Helix

Abstract

Sarma et al. (J. Biomol. Str. and Dynam. 2, 1085 (1985) have proposed, on the basis of nuclear magnetic resonance experiments on the complex of netropsin with poly(dA)·poly(dT), that the drug molecule lies asymmetrically along the dA side of the minor groove and makes hydrogen bonds only with the dA strand. If the crystal structure analyses of B-DNA (Fratini et al., J. Biol. Chem. 257, 14686 (1982)) and of its complex with netropsin (Kopka et al., J. Mol. Biol. 183, 553 (1985)) are any guide, this off-center, wide-groove model is stereochemically unlikely. More to the point, the off-center model is unnecessary to explain the observed nmr data. All of the nuclear Overhauser and other observations are fully explained by the structure seen in the x-ray crystal analysis, in which netropsin sits squarely centered within the minor groove, making bifurcated hydrogen bonds with both strands.     

Structure of (dG)n · (dC)n Under Superhelical Stress and Acid pH

Abstract

We have recently shown that a (GA)_n · (TC)_n tract undergoes a sharp structural transition under superhelical stress (V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 2, 327 (1985)). Unlike the well studied transitions to the cruciform and to the Z form, this novel transition was strongly pH-dependent. We have found the (dG)_n · (dC)_n insert to undergo a pH-dependent structural transition similar to that of the (GA)_n · (TC)_n tract. These new data meet our earlier expectations and disagree with the data of D.E. Pulleyblank, D.B. Haniford and A.R. Morgan, Cell 42, 271 (1985). We conclude that a novel DNA structure (the H-form) is typical of homopurine-homopyrimidine mirror repeats (the H palindromes) under superhelical stress and/or acid pH. In the H-form the homopyrimidine strand forms a hairpin while half of the homopurine strand interacts with the hairpin forming a triplex, the other half of the homopurine strand being unstructured (V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 2, 3,667 (1986)).     

ENGINEERING DNA TOPOLOGY WITH LOCKED NUCLEOSIDES: A STRUCTURAL STUDY

DNA dodecamers modified with nucleotide building blocks based on a bicyclo[3.1.0]hexane system that effectively “locks” the ribose template into an RNA-like or “North” (N) conformation were analyzed by various biophysical techniques including high field nuclear magnetic resonance (NMR). Replacement of either one or both of the center thymidines in the Dickerson Drew dodecamer (CGCGAAT*T*CGCG) caused a progressive shift in the bending propensity of the double helix as shown by a newly developed rapid technique that compares the residual dipolar coupling (RDC) values of the modified duplexes with those previously determined for the native DNA.     

Improved Procedure for the Reduction of N - 1 Content in Synthetic Oligonucleotides

Abstract

By incorporating a “capping step” at the start of an oligonucleotide synthesis (“pre-cap”) and following a “SUP” work-up protocol with ammonium hydroxide, an overall improvement is observed in the quality of oligonucleotides synthesized on a large scale on controlled pore glass support (CPG). Rationalization of these results is provided.