1
|
Gelenter MD, Chen KJ, Hong M. Off-resonance 13C- 2H REDOR NMR for site-resolved studies of molecular motion. JOURNAL OF BIOMOLECULAR NMR 2021; 75:335-345. [PMID: 34342847 PMCID: PMC8830769 DOI: 10.1007/s10858-021-00377-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Accepted: 07/19/2021] [Indexed: 06/06/2023]
Abstract
We introduce a 13C-2H Rotational Echo DOuble Resonance (REDOR) technique that uses the difference between on-resonance and off-resonance 2H irradiation to detect dynamic segments in deuterated molecules. By selectively inverting specific regions of the 2H magic-angle spinning (MAS) sideband manifold to recouple some of the deuterons to nearby carbons, we distinguish dynamic and rigid residues in 1D and 2D 13C spectra. We demonstrate this approach on deuterated GB1, H/D exchanged GB1, and perdeuterated bacterial cellulose. Numerical simulations reproduce the measured mixing-time and 2H carrier-frequency dependence of the REDOR dephasing of bacterial cellulose. Combining numerical simulations with experiments thus allow the extraction of motionally averaged quadrupolar couplings from REDOR dephasing values.
Collapse
Affiliation(s)
- Martin D Gelenter
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, MA, 02139, Cambridge, USA
| | - Kelly J Chen
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, MA, 02139, Cambridge, USA
| | - Mei Hong
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, MA, 02139, Cambridge, USA.
| |
Collapse
|
2
|
Meirovitch E, Liang Z, Freed JH. Protein dynamics in the solid-state from 2H NMR lineshape analysis. III. MOMD in the presence of Magic Angle Spinning. SOLID STATE NUCLEAR MAGNETIC RESONANCE 2018; 89:35-44. [PMID: 29208317 PMCID: PMC5772661 DOI: 10.1016/j.ssnmr.2017.11.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2017] [Revised: 11/11/2017] [Accepted: 11/13/2017] [Indexed: 06/07/2023]
Abstract
We report on a new approach to the analysis of dynamic NMR lineshapes from polycrystalline (i.e., macroscopically disordered) samples in the presence of Magic Angle Spinning (MAS). This is an application of the Stochastic Liouville Equation developed by Freed and co-workers for treating restricted (i.e., microscopically ordered) motions. The 2H nucleus in an internally-mobile C-CD3 moiety serves as a prototype probe. The acronym is 2H/MOMD/MAS, where MOMD stands for "microscopic-order-macroscopic-disorder." The key elements describing internal motions - their type, the local spatial restrictions, and related features of local geometry - are treated in MOMD generally, within their rigorous three-dimensional tensorial requirements. Based on this representation a single physically well-defined model of local motion has the capability of reproducing experimental spectra. There exist other methods for analyzing dynamic 2H/MAS spectra which advocate simple motional modes. Yet, to reproduce satisfactorily the experimental lineshapes, one has either to use unusual parameter values, or combine several simple motional modes. The multi-simple-mode reasoning assumes independence of the constituent modes, features ambiguity as different simple modes may be used, renders inter-system comparison difficult as the overall models differ, and makes possible model-improvement only by adding yet another simple mode, i.e., changing the overall model. 2H/MOMD/MAS is free of such limitations and inherently provides a clear physical interpretation. These features are illustrated. The advantage of 2H/MOMD/MAS in dealing with sensitive but hardly investigated slow-motional lineshapes is demonstrated by applying it to actual experimental data. The results differ from those obtained previously with a two-site exchange scheme that yielded unusual parameters.
Collapse
Affiliation(s)
- Eva Meirovitch
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat-Gan 52900, Israel.
| | - Zhichun Liang
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853-1301, USA
| | - Jack H Freed
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853-1301, USA
| |
Collapse
|
3
|
Kumar R, Ries A, Wengel J. Synthesis and Excellent Duplex Stability of Oligonucleotides Containing 2'-Amino-LNA Functionalized with Galactose Units. Molecules 2017; 22:molecules22050852. [PMID: 28531137 PMCID: PMC6153924 DOI: 10.3390/molecules22050852] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2017] [Revised: 05/17/2017] [Accepted: 05/17/2017] [Indexed: 02/07/2023] Open
Abstract
A convenient method for the preparation of oligonucleotides containing internally-attached galactose and triantennary galactose units has been developed based on click chemistry between 2′-N-alkyne 2′-amino-LNA nucleosides and azido-functionalized galactosyl building blocks. The synthesized oligonucleotides show excellent binding affinity and selectivity towards complementary DNA/RNA strands with an increase in the melting temperature of up to +23.5 °C for triply-modified variants.
Collapse
Affiliation(s)
- Rajesh Kumar
- Biomolecular Nanoscale Engineering Center, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark.
- Department of Chemistry, University of Delhi, Delhi 110007, India.
| | - Annika Ries
- Biomolecular Nanoscale Engineering Center, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark.
| | - Jesper Wengel
- Biomolecular Nanoscale Engineering Center, Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark.
| |
Collapse
|
4
|
Meirovitch E, Liang Z, Freed JH. Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective. J Phys Chem B 2015; 119:2857-68. [PMID: 25594631 DOI: 10.1021/jp511386b] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Deuterium line shape analysis of CD3 groups has emerged as a particularly useful tool for studying microsecond-millisecond protein motions in the solid state. The models devised so far consist of several independently conceived simple jump-type motions. They are comprised of physical quantities encoded in their simplest form; improvements are only possible by adding yet another simple motion, thereby changing the model. The various treatments developed are case-specific; hence comparison among the different systems is not possible. Here we develop a new methodology for (2)H NMR line shape analysis free of these limitations. It is based on the microscopic-order-macroscopic-disorder (MOMD) approach. In MOMD motions are described by diffusion tensors, spatial restrictions by potentials/ordering tensors, and geometric features by relative tensor orientations. Jump-type motions are recovered in the limit of large orientational potentials. Model improvement is accomplished by monitoring the magnitude, symmetry, and orientation of the various tensors. The generality of MOMD makes possible comparison among different scenarios. CD3 line shapes from the Chicken Villin Headpiece Subdomain and the Streptomyces Subtilisin Inhibitor are used as experimental examples. All of these spectra are reproduced by using rhombic local potentials constrained for simplicity to be given by the L = 2 spherical harmonics, and by axial diffusion tensors. Potential strength and rhombicity are found to be ca. 2-3 k(B)T. The diffusion tensor is tilted at 120° from the C-CD3 axis. The perpendicular (parallel) correlation times for local motion are 0.1-1.0 ms (3.3-30 μs). Activation energies in the 1.1-8.0 kcal/mol range are estimated. Future prospects include extension to the (2)H relaxation limit, application to the (15)N and (13)C NMR nuclei, and accounting for collective motions and anisotropic media.
Collapse
Affiliation(s)
- Eva Meirovitch
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University , Ramat-Gan 52900, Israel
| | | | | |
Collapse
|
5
|
Astakhova IK, Wengel J. Scaffolding along nucleic acid duplexes using 2'-amino-locked nucleic acids. Acc Chem Res 2014; 47:1768-77. [PMID: 24749544 DOI: 10.1021/ar500014g] [Citation(s) in RCA: 54] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
CONSPECTUS: Incorporation of chemically modified nucleotide scaffolds into nucleic acids to form assemblies rich in function is an innovative area with great promise for nanotechnology and biomedical and material science applications. The intrinsic biorecognition potential of nucleic acids combined with advanced properties of the locked nucleic acids (LNAs) provide opportunities to develop new nanomaterials and devices like sensors, aptamers, and machines. In this Account, we describe recent research on preparation and investigation of the properties of LNA/DNA hybrids containing functionalized 2'-amino-LNA nucleotides. By application of different chemical reactions, modification of 2'-amino-LNA scaffolds can be efficiently performed in high yields and with various tags, postsynthetically or during the automated oligonucleotide synthesis. The choice of a synthetic method for scaffolding along 2'-amino-LNA mainly depends on the chemical nature of the modification, its price, its availability, and applications of the product. One of the most useful applications of the product LNA/DNA scaffolds containing 2'-amino-LNA is to detect complementary DNA and RNA targets. Examples of these applications include sensing of clinically important single-nucleotide polymorphisms (SNPs) and imaging of nucleic acids in vitro, in cell culture, and in vivo. According to our studies, 2'-amino-LNA scaffolds are efficient within diagnostic probes for DNA and RNA targets and as therapeutics, whereas both 2'-amino- and isomeric 2'-α-l-amino-LNA scaffolds have promising properties for stabilization and detection of DNA nanostructures. Attachment of fluorescent groups to the 2'-amino group results in very high fluorescent quantum yields of the duplexes and remarkable sensitivity of the fluorescence signal to target binding. Notably, fluorescent LNA/DNA probes bind nucleic acid targets with advantages of high affinity and specificity. Thus, molecular motion of nanodevices and programmable self-assembly of chemically modified LNA/DNA nanomaterials can be followed by bright fluorescence signaling from the functionalized LNA units. Another appealing aspect of the amino-LNA scaffolds is specific targeting of nucleic acids and proteins for therapeutic applications. 2'-Amino-LNA/DNA conjugates containing peptide and polyaromatic hydrocarbon (PAH) groups are promising in this context as well as for advanced imaging and diagnostic purposes in vivo. For imaging applications, photostability of fluorescence dyes is of crucial importance. Chemically stable and photostable fluorescent PAH molecules attached to 2'-amino functionality of the 2'-amino-LNA are potent for in vitro and in vivo imaging of DNA and RNA targets. We believe that rational evolution of the biopolymers of Nature may solve the major challenges of the future material science and biomedicine. However, this requires strong scientific progress and efficient interdisciplinary research. Examples of this Account demonstrate that among other synthetic biopolymers, synthetic nucleic acids containing functionalized 2'-amino-LNA scaffolds offer great opportunities for material science, diagnostics, and medicine of the future.
Collapse
Affiliation(s)
- I. Kira Astakhova
- Nucleic Acid Center,
Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Jesper Wengel
- Nucleic Acid Center,
Department
of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| |
Collapse
|
6
|
Kinde-Carson MN, Ferguson C, Oyler NA, Harbison GS, Meints GA. Solid state 2H NMR analysis of furanose ring dynamics in DNA containing uracil. J Phys Chem B 2010; 114:3285-93. [PMID: 20151717 DOI: 10.1021/jp9091656] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions, such as uracil, in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Increased dynamic properties in lesion-containing DNAs have been suggested to assist recognition and specificity. Deuterium solid-state nuclear magnetic resonance (SSNMR) has been used to directly observe local dynamics of the furanose ring within a uracil:adenine (U:A) base pair and compared to a normal thymine:adenine (T:A) base pair. Quadrupole echo lineshapes, <T(1Z)>, and <T(2e)> relaxation data were collected, and computer modeling was performed. The results indicate that the relaxation times are identical within the experimental error, the solid lineshapes are essentially indistinguishable above the noise level, and our lineshapes are best fit with a model that does not have significant local motions. Therefore, U:A base pair furanose rings appear to have essentially identical dynamic properties as a normal T:A base pair, and the local dynamics of the furanose ring are unlikely to be the sole arbiter for uracil recognition and specificity in U:A base pairs.
Collapse
|
7
|
Meirovitch E, Shapiro YE, Polimeno A, Freed JH. Structural dynamics of bio-macromolecules by NMR: the slowly relaxing local structure approach. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2010; 56:360-405. [PMID: 20625480 PMCID: PMC2899824 DOI: 10.1016/j.pnmrs.2010.03.002] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Affiliation(s)
- Eva Meirovitch
- The Mina and Everard Goodman Faculty of Life Sciences, Bar–Ilan University, Ramat-Gan 52900 Israel
| | - Yury E. Shapiro
- The Mina and Everard Goodman Faculty of Life Sciences, Bar–Ilan University, Ramat-Gan 52900 Israel
| | - Antonino Polimeno
- Department of Physical Chemistry, University of Padua, 35131 Padua, Italy
| | - Jack H. Freed
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853-1301, U.S.A
| |
Collapse
|
8
|
Brut M, Estève A, Landa G, Renvez G, Djafari Rouhani M. The Static Modes: an alternative approach for the treatment of macro- and bio-molecular induced-fit flexibility. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2009; 28:17-25. [PMID: 19139941 DOI: 10.1140/epje/i2008-10397-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2008] [Accepted: 11/19/2008] [Indexed: 05/27/2023]
Abstract
We present a new competitive method for the atomic scale treatment of macromolecular flexibility called Static Mode method. This method is based on the "induced-fit" concept, i.e. it maps the intrinsic deformations of a macromolecule subject to diverse external excitations. The algorithm makes it possible to obtain a set of deformations, each one corresponding to a specific interaction on a specific molecular site, in terms of force constants contained in the energy model. In this frame, the docking problem can be expressed in terms of interaction sites between the two molecules, the molecular deformations being extracted from the pre-calculated Static Modes of each molecule. Some preliminary basic examples aimed at illustrating potential applications where macro- or bio-molecular flexibility is of key importance are given: flexibility inducing conformational changes in the case of furanose ring and flexibility for the characterization, including allostery, of poly(N-isopropylacrylamide)(P-NIPAM) active sites. We also discuss how this procedure allows "induced-fit" flexible molecular docking, beyond state-of-the-art semi-rigid methods.
Collapse
Affiliation(s)
- M Brut
- LAAS-CNRS, Université de Toulouse, 7 avenue du Colonel Roche, F-31077 Toulouse, France.
| | | | | | | | | |
Collapse
|
9
|
Echodu D, Goobes G, Shajani Z, Pederson K, Meints G, Varani G, Drobny G. Furanose dynamics in the HhaI methyltransferase target DNA studied by solution and solid-state NMR relaxation. J Phys Chem B 2008; 112:13934-44. [PMID: 18844399 PMCID: PMC2735271 DOI: 10.1021/jp801723x] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Both solid-state and solution NMR relaxation measurements are routinely used to quantify the internal dynamics of biomolecules, but in very few cases have these two techniques been applied to the same system, and even fewer attempts have been made so far to describe the results obtained through these two methods through a common theoretical framework. We have previously collected both solution 13C and solid-state 2H relaxation measurements for multiple nuclei within the furanose rings of several nucleotides of the DNA sequence recognized by HhaI methyltransferase. The data demonstrated that the furanose rings within the GCGC recognition sequence are very flexible, with the furanose rings of the cytidine, which is the methylation target, experiencing the most extensive motions. To interpret these experimental results quantitatively, we have developed a dynamic model of furanose rings based on the analysis of solid-state 2H line shapes. The motions are modeled by treating bond reorientations as Brownian excursions within a restoring potential. By applying this model, we are able to reproduce the rates of 2H spin-lattice relaxation in the solid and 13C spin-lattice relaxation in solution using comparable restoring force constants and internal diffusion coefficients. As expected, the 13C relaxation rates in solution are less sensitive to motions that are slower than overall molecular tumbling than to the details of global molecular reorientation, but are somewhat more sensitive to motions in the immediate region of the Larmor frequency. Thus, we conclude that the local internal motions of this DNA oligomer in solution and in the hydrated solid state are virtually the same, and we validate an approach to the conjoint analysis of solution and solid-state NMR relaxation and line shapes data, with wide applicability to many biophysical problems.
Collapse
Affiliation(s)
- Dorothy Echodu
- Department of Chemistry, University of Washington, Seattle, WA 98195-1700, USA
| | | | | | | | | | | | | |
Collapse
|
10
|
Duchardt E, Nilsson L, Schleucher J. Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study. Nucleic Acids Res 2008; 36:4211-9. [PMID: 18579564 PMCID: PMC2475628 DOI: 10.1093/nar/gkn375] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Using (13)C spin relaxation NMR in combination with molecular dynamic (MD) simulations, we characterized internal motions within double-stranded DNA on the pico- to nano-second time scale. We found that the C-H vectors in all cytosine ribose moieties within the Dickerson-Drew dodecamer (5'-CGCGAATTCGCG-3') are subject to high amplitude motions, while the other nucleotides are essentially rigid. MD simulations showed that repuckering is a likely motional model for the cytosine ribose moiety. Repuckering occurs with a time constant of around 100 ps. Knowledge of DNA dynamics will contribute to our understanding of the recognition specificity of DNA-binding proteins such as cytosine methyltransferase.
Collapse
Affiliation(s)
- Elke Duchardt
- Department of Medical Biochemistry and Biophysics, Umeå University, SE-901 87 Umeå, Sweden.
| | | | | |
Collapse
|
11
|
Pederson K, Meints GA, Shajani Z, Miller PA, Drobny GP. Backbone dynamics in the DNA HhaI protein binding site. J Am Chem Soc 2008; 130:9072-9. [PMID: 18570423 DOI: 10.1021/ja801243d] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.
Collapse
Affiliation(s)
- Kari Pederson
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| | | | | | | | | |
Collapse
|
12
|
Meints GA, Miller PA, Pederson K, Shajani Z, Drobny G. Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site. J Am Chem Soc 2008; 130:7305-14. [PMID: 18489097 DOI: 10.1021/ja075775n] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The dynamics of the furanose rings in the GCGC moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied by using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs selectively deuterated on the furanose rings of nucleotides within the 5'-GCGC-3' moiety indicated that all of these positions are structurally flexible. The furanose ring within the deoxycytidine that is the methylation target displays the largest-amplitude structural changes according to the observed deuterium NMR line shapes, whereas the furanose rings of nucleotides more remote from the methylation site have less-mobile furanose rings (i.e., with puckering amplitudes < 0.3 A). Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through the furanose ring. These NMR data indicate that the 5'-GCGC-3' is dynamic, with the largest-amplitude motions occurring nearest the methylation site. The inherent flexibility of this moiety in DNA makes the molecule more amenable to the large-amplitude structural rearrangements that must occur when the DNA binds to the HhaI methyltransferase.
Collapse
Affiliation(s)
- Gary A Meints
- Department of Chemistry, Missouri State University, Springfield, Missouri 65897, and Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| | | | | | | | | |
Collapse
|
13
|
Ivanova A, Rösch N. The structure of LNA:DNA hybrids from molecular dynamics simulations: the effect of locked nucleotides. J Phys Chem A 2007; 111:9307-19. [PMID: 17718546 DOI: 10.1021/jp073198j] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Locked nucleic acids (LNAs) exhibit a modified sugar fragment that is restrained to the C3'-endo conformation. LNA-containing duplexes are rather stable and have a more rigid structure than DNA duplexes, with a propensity for A-conformation of the double helix. To gain detailed insight into the local structure of LNA-modified DNA oligomers (as a foundation for subsequent exploration of the electron-transfer capabilities of such modified duplexes), we carried out molecular dynamics simulations on a set of LNA:DNA 9-mer duplexes and analyzed the resulting structures in terms of base step parameters and the conformations of the sugar residues. The perturbation introduced by a single locked nucleotide was found to be fairly localized, extending mostly to the first neighboring base pairs; such duplexes featured a B-type helix. With increasing degree of LNA modification the structure gradually changed; the duplex with one complete LNA strand assumed a typical A-DNA structure. The relative populations of the sugar conformations agreed very well with NMR data, lending credibility to the validity of the computational protocol.
Collapse
Affiliation(s)
- Anela Ivanova
- Department Chemie, Theoretische Chemie, Technische Universität München, 85747 Garching, Germany
| | | |
Collapse
|
14
|
Miller PA, Shajani Z, Meints GA, Caplow D, Goobes G, Varani G, Drobny GP. Contrasting views of the internal dynamics of the HhaI methyltransferase target DNA reported by solution and solid-state NMR spectroscopy. J Am Chem Soc 2007; 128:15970-1. [PMID: 17165714 DOI: 10.1021/ja066329n] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Solution and solid-state NMR have been used conjointly to probe the internal motions of a DNA dodecamer containing the recognition site for the HhaI methyltransferase. The results strongly suggest that ns-mus motions contribute to the functionally relevant dynamic properties of nucleic acids during DNA methylation.
Collapse
Affiliation(s)
- Paul A Miller
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | | | | | | | | | | | | |
Collapse
|
15
|
Banavali NK, Roux B. Free energy landscape of A-DNA to B-DNA conversion in aqueous solution. J Am Chem Soc 2005; 127:6866-76. [PMID: 15869310 DOI: 10.1021/ja050482k] [Citation(s) in RCA: 109] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.
Collapse
Affiliation(s)
- Nilesh K Banavali
- Department of Physiology, Biophysics, and Systems Biology, Weill Medical College of Cornell University, 1300 York Avenue, New York, NY 10021, USA.
| | | |
Collapse
|
16
|
Shishkin OV, Gorb L, Zhikol OA, Leszczynski J. Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides. J Biomol Struct Dyn 2005; 22:227-44. [PMID: 15317483 DOI: 10.1080/07391102.2004.10506998] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The molecular structure of different conformers of isolated canonical purine 2'-deoxyribonucleotides 2-deoxyadenosine-5'-phosphate (pdA) and 2'-deoxyguanosine-5'-phosphate (pdG) was optimized using the B3LYP/6-31G(d) method. The results of the calculations reveal that the geometrical parameters and relative stability of the conformers significantly depend on the nature of the nucleobase, its orientation, the conformation of the furanose ring, the charge of the phosphate group and the character of the intramolecular hydrogen bonds. Analysis of the electron density distribution in purine nucleotides reveals the existence of a number of intramolecular hydrogen bonds. In general, the south conformer has a lower energy at the anti orientation of the base, and both conformers occur as the most stable for the syn orientation of the nucleobases. The results of the calculations reveal that the geometry and relative energy of the conformers of purine DNTs may be easily tuned by the charge of the phosphate group.
Collapse
Affiliation(s)
- Oleg V Shishkin
- Institute for Scintillation Materials, National Academy of Science of Ukraine, 60 Lenina Ave., Kharkiv 61001, Ukraine.
| | | | | | | |
Collapse
|
17
|
Shishkin OV, Gorb L, Zhikol OA, Leszczynski J. Conformational analysis of canonical 2-deoxyribonucleotides. 1. Pyrimidine nucleotides. J Biomol Struct Dyn 2004; 21:537-54. [PMID: 14692798 DOI: 10.1080/07391102.2004.10506947] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
The molecular structure and relative stability of north and south conformers of 2'-deoxyribonucleotides containing pyrimidine nucleic acid bases ( 2'-deoxythymidilic (pdT), 2'-deoxycytidilic (pdC) acids and their mono- and dianions) have been obtained and analyzed at the DFT/B3LYP level using the standard 6-31G(d) basis set. We have revealed that, when the nucleobase moiety is incorporated into the nucleotides, it maintains a nonplanar and nonrigid conformation due to out-of-plane deformation of the amino group and pyrimidine ring. It has been demonstrated that an increase of negative charge of the phosphate group results in increase of amino group pyramidalization, discrimination between conformers with syn and anti orientation of base with respect to sugar, strengthening of intramolecular C-H.O hydrogen bonds leading to deformation and fixation of geometry of nucleotides, and weakening of phosphodiester bond. These results allow to make suggestions about sources of twist and buckle deformations of base pairs, mechanisms of repaire of DNA via change of base orientation, and conditions for breakage of the P-O bonds during hydrolysis.
Collapse
Affiliation(s)
- Oleg V Shishkin
- Institute for Scintillation Materials, National Academy of Science of Ukraine, 60 Lenina Ave, Kharkiv 61001, Ukraine.
| | | | | | | |
Collapse
|
18
|
Djuranovic D, Hartmann B. Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states. J Biomol Struct Dyn 2003; 20:771-88. [PMID: 12744707 DOI: 10.1080/07391102.2003.10506894] [Citation(s) in RCA: 67] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.
Collapse
Affiliation(s)
- D Djuranovic
- Laboratoire de Biochimie Theorique, CNRS UPR 9080, Institut de Biologie Physico-chimique, 13 rue P. et M. Curie, Paris 75005, France.
| | | |
Collapse
|
19
|
Freedberg DI. An alternative method for pucker determination in carbohydrates from residual dipolar couplings: a solution NMR study of the fructofuranosyl ring of sucrose. J Am Chem Soc 2002; 124:2358-62. [PMID: 11878992 DOI: 10.1021/ja010981v] [Citation(s) in RCA: 49] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A simple solution NMR method is presented for pucker determination of five-membered rings using only residual dipolar couplings obtained in a single liquid crystalline medium, DMPC/DHPC bicelles (DMPC = dimyristoylphosphatidylcholine; DHPC = dihexanoylphosphatidylcholine). The method was applied to determine the pucker of the fructofuranosyl ring of sucrose. The results indicate a fructofuranosyl pucker phase in the 20 degrees - 70 degrees range. The pucker phases are in agreement with those from previous NMR and optical spectroscopic studies and, importantly, do not rely on empirically parametrized Karplus curves. Furthermore, the analysis implies more than one stable pucker phase and rapid ring interconversion in this range. The present results suggest that using residual dipolar couplings alone can reveal multiple conformations present in solution. Hence, when a sufficient number of residual dipolar coupling constants is measured, the outcome is a robust, reliable, and independent route for determining carbohydrate and nucleic acid structure by NMR spectroscopy.
Collapse
Affiliation(s)
- Darón I Freedberg
- Laboratory of Biophysics, CBER/FDA, 29 Lincoln Drive, Building 29 HFM-419, Bethesda, Maryland 20892, USA.
| |
Collapse
|