Modeling furanose ring dynamics in DNA

Off-resonance 13C-2H REDOR NMR for site-resolved studies of molecular motion

We introduce a ¹³C-²H Rotational Echo DOuble Resonance (REDOR) technique that uses the difference between on-resonance and off-resonance ²H irradiation to detect dynamic segments in deuterated molecules. By selectively inverting specific regions of the ²H 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 ¹³C spectra. We demonstrate this approach on deuterated GB1, H/D exchanged GB1, and perdeuterated bacterial cellulose. Numerical simulations reproduce the measured mixing-time and ²H 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.

Protein dynamics in the solid-state from 2H NMR lineshape analysis. III. MOMD in the presence of Magic Angle Spinning

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 ²H nucleus in an internally-mobile C-CD₃ moiety serves as a prototype probe. The acronym is ²H/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 ²H/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. ²H/MOMD/MAS is free of such limitations and inherently provides a clear physical interpretation. These features are illustrated. The advantage of ²H/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.

Synthesis and Excellent Duplex Stability of Oligonucleotides Containing 2'-Amino-LNA Functionalized with Galactose Units

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.

Protein dynamics in the solid state from 2H NMR line shape analysis: a consistent perspective

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.

Scaffolding along nucleic acid duplexes using 2'-amino-locked nucleic acids

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.

Solid state 2H NMR analysis of furanose ring dynamics in DNA containing uracil

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.

The Static Modes: an alternative approach for the treatment of macro- and bio-molecular induced-fit flexibility

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.

Furanose dynamics in the HhaI methyltransferase target DNA studied by solution and solid-state NMR relaxation

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.

Cytosine ribose flexibility in DNA: a combined NMR 13C spin relaxation and molecular dynamics simulation study

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.

Backbone dynamics in the DNA HhaI protein binding site

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.

Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site

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.

The structure of LNA:DNA hybrids from molecular dynamics simulations: the effect of locked nucleotides

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.

Contrasting views of the internal dynamics of the HhaI methyltransferase target DNA reported by solution and solid-state NMR spectroscopy

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.

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

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.

Conformational analysis of canonical 2-deoxyribonucleotides. 2. Purine nucleotides

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.

Conformational analysis of canonical 2-deoxyribonucleotides. 1. Pyrimidine nucleotides

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.

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

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.

An alternative method for pucker determination in carbohydrates from residual dipolar couplings: a solution NMR study of the fructofuranosyl ring of sucrose

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.