The B- to A-DNA transition and the reorganization of solvent at the DNA surface

Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning

Conformational dynamics play a crucial role in determining the behavior of the biomolecules. Polarizable force fields, such as AMOEBA, can accurately capture electrostatic interactions underlying the conformational space. However, applying a polarizable force field in molecular dynamics (MD) simulations can be computationally expensive, especially in studying long-time-scale dynamics. To overcome this challenge, we incorporated the AMOEBA potential with Milestoning, an enhanced sampling method in this work. This integration allows us to efficiently sample the rare and important conformational states of a biomolecule by using many short and independent molecular dynamics trajectories with the AMOEBA force field. We applied this method to investigate the conformational dynamics of alanine dipeptide, DNA, and RNA A-B form conversion. Well-converged thermodynamic and kinetic properties were obtained, including the free energy difference, mean first passage time, and critical transitions between states. Our results demonstrate the power of integrating polarizable force fields with enhanced sampling methods in quantifying the thermodynamic and kinetic properties of biomolecules at the atomic level.

Base pair compositional variability influences DNA structural stability and tunes hydration thermodynamics and dynamics

DNA deformability and differential hydration are crucial determinants of biological processes ranging from genetic material packaging to gene expression; their associative details, however, remain inadequately understood. Herein, we report investigations of the dynamic and thermodynamic responses of the local hydration of a variety of base pair sequences. Leveraging in silico sampling and our in-house analyses, we first report the local conformational propensity of sequences that are either predisposed toward the canonical A- or B-conformations or are restrained to potential transitory pathways. It is observed that the transition from the unrestrained A-form to the B-form leads to lengthwise structural deformation. The insertion of intermittent -(CG)- base pairs in otherwise homogeneous -(AT)- sequences bears dynamical consequences for the vicinal hydration layer. Calculation of the excess (pair) entropy suggests substantially higher values of hydration water surrounding A conformations over the B- conformations. Applying the Rosenfeld approximation, we project that the diffusivity of water molecules proximal to canonical B conformation is least for the minor groove of the canonical B-conformation. We determine that structure, composition, and conformation specific groove dimension together influence the local hydration characteristics and, therefore, are expected to be important determinants of biological processes.

Structural properties and binding mechanism of DNA aptamers sensing saliva melatonin for diagnosis and monitoring of circadian clock and sleep disorders

Circadian desynchrony with the external light-dark cycle influences the rhythmic secretion of melatonin which is among the first signs of circadian rhythm sleep disorders. An accurate dim light melatonin onset (established indicator of circadian rhythm sleep disorders) measurement requires lengthy assays, and antibody affinities alterations, especially in patients with circadian rhythm disorders whose melatonin salivary levels vary significantly, making antibodies detection mostly inadequate. In contrast, aptamers with their numerous advantages (e.g., target selectivity, structural flexibility in tuning binding affinities, small size, etc.) can become preferable biorecognition molecules for salivary melatonin detection with high sensitivity and specificity. This study thoroughly characterizes the structural property and binding mechanism of a single-stranded DNA aptamer full sequence (MLT-C-1) and its truncated versions (MLT-A-2, MLT-A-4) to decipher its optimal characteristics for saliva melatonin detection. We use circular dichroism spectroscopy to determine aptamers' conformational changes under different ionic strengths and showed that aptamers display a hairpin loop structure where few base pairs in the stem play a significant role in melatonin binding and formation of aptamer stabilized structure. Through microscale thermophoresis, aptamers demonstrated a high binding affinity in saliva samples (MLT-C-1F K_d = 12.5 ± 1.7 nM; MLT-A-4F K_d = 11.2 ± 1.6 nM; MLT-A-2F K_d = 2.4 ± 2.8 nM; limit-of-detection achieved in pM, highest sensitivity attained for MLT-A-2F aptamer with the lowest detection limit of 1.35 pM). Our data suggest that aptamers are promising as biorecognition molecules and provide the baseline parameters for the development of an aptamer-based point-of-care diagnostic system for melatonin detection and accurate profiling of its fluctuations in saliva.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Structure-Activity Relationships in Nucleic-Acid-Templated Vectors Based on Peptidic Dynamic Covalent Polymers

The use of nucleic acids as templates, which can trigger the self-assembly of their own vectors represent an emerging, simple and versatile, approach toward the self-fabrication of tailored nucleic acids delivery vectors. However, the structure-activity relationships governing this complex templated self-assembly process that accompanies the complexation of nucleic acids remains poorly understood. Herein, the class of arginine-rich dynamic covalent polymers (DCPs) composed of different monomers varying the number and position of arginines were studied. The combinations that lead to nucleic acid complexation, in saline buffer, using different templates, from short siRNA to long DNA, are described. Finally, a successful peptidic DCP featuring six-arginine repeating unit that promote the safe and effective delivery of siRNA in live cancer cells was identified.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Thermodynamic and kinetic properties of a single base pair in A-DNA and B-DNA

Double stranded DNA can adopt different forms, the so-called A-, B-, and Z-DNA, which play different biological roles. In this work, the thermodynamic and the kinetic parameters for the base-pair closing and opening in A-DNA and B-DNA were calculated by all-atom molecular dynamics simulations at different temperatures. The thermodynamic parameters of the base pair in B-DNA were in good agreement with the experimental results. The free energy barrier of breaking a single base stack results from the enthalpy increase ΔH caused by the disruption of hydrogen bonding and base-stacking interactions, as well as water and base interactions. The free energy barrier of base pair closing comes from the unfavorable entropy loss ΔS caused by the restriction of torsional angles and hydration. It was found that the enthalpy change ΔH and the entropy change ΔS for the base pair in A-DNA are much larger than those in B-DNA, and the transition rates between the opening and the closing state for the base pair in A-DNA are much slower than those in B-DNA. The large difference of the enthalpy and entropy change for forming the base pair in A-DNA and B-DNA results from different hydration in A-DNA and B-DNA. The hydration pattern observed around DNA is an accompanying process for forming the base pair, rather than a follow-up of the conformation.

Kinetics of the Conformational Transformation between B- and A-Forms in the Drew-Dickerson Dodecamer

Some DNA sequences in crystals and in complexes with proteins can exist in the forms intermediate between the B- and A-DNA. Based on this, it was implied that the B-to-A transition for any DNA molecule should go through these intermediate forms also in kinetics. More precisely, the helix parameter Slide has to change first, and the molecule should take the E-form. After that, the Roll parameter changes. In the present work, we simulated the kinetics of the B-A transition in the Drew-Dickerson dodecamer, a known B-philic DNA oligomer. We used the "sugar" coarse-grained model that reproduces ribose flexibility, preserves sequence specificity, employs implicit water and explicit ions, and offers the possibility to vary friction. As the control parameter of the transition, we chose the volume available for a counterion and considered the change from a large to a small volume. In the described system, the B-to-A conformational transformation proved to correspond to a first-order phase transition. The molecule behaves like a small cluster in the region of such a transition, jumping between the A- and B-forms in a wide range of available volumes. The viscosity of the solvent does not affect the midpoint of the transition but only the overall mobility of the system. All helix parameters change synchronously on average, we have not observed the sequence "Slide first, Roll later" in kinetics, and the E-DNA is not a necessary step for the transition between the B- and A-forms in the studied system. So, the existence of the intermediate DNA forms requires specific conditions, shifting the common balance of interactions: certain nucleotide sequence in specific solution or/and the interaction with some protein.

A New Look into the Mode of Action of Metal-Based Anticancer Drugs

The mode of action of Pt- and Pd-based anticancer agents (cisplatin and Pd₂Spm) was studied by characterising their impact on DNA. Changes in conformation and mobility at the molecular level in hydrated DNA were analysed by quasi-elastic and inelastic neutron scattering techniques (QENS and INS), coupled to Fourier transform infrared (FTIR) and microRaman spectroscopies. Although INS, FTIR and Raman revealed drug-triggered changes in the phosphate groups and the double helix base pairing, QENS allowed access to the nanosecond motions of the biomolecule’s backbone and confined hydration water within the minor groove. Distinct effects were observed for cisplatin and Pd₂Spm, the former having a predominant effect on DNA’s spine of hydration, whereas the latter had a higher influence on the backbone dynamics. This is an innovative way of tackling a drug’s mode of action, mediated by the hydration waters within its pharmacological target (DNA).

Hydrophobic catalysis and a potential biological role of DNA unstacking induced by environment effects

The main stabilizer of the DNA double helix is not the base-pair hydrogen bonds but coin-pile stacking of base pairs, whose hydrophobic cohesion, requiring abundant water, indirectly makes the DNA interior dry so that hydrogen bonds can exert full recognition power. We report that certain semihydrophobic agents depress the stacking energy (measurable in single-molecule experiments), leading to transiently occurring holes in the base-pair stack (monitorable via binding of threading intercalators). Similar structures observed in DNA complexes with RecA and Rad51, and previous observations of spontaneous strand exchange catalyzed in semihydrophobic model systems, make us propose that some hydrophobic protein residues may have roles in catalyzing homologous recombination. We speculate that hydrophobic catalysis is a general phenomenon in DNA enzymes.

Hydrophobic base stacking is a major contributor to DNA double-helix stability. We report the discovery of specific unstacking effects in certain semihydrophobic environments. Water-miscible ethylene glycol ethers are found to modify structure, dynamics, and reactivity of DNA by mechanisms possibly related to a biologically relevant hydrophobic catalysis. Spectroscopic data and optical tweezers experiments show that base-stacking energies are reduced while base-pair hydrogen bonds are strengthened. We propose that a modulated chemical potential of water can promote “longitudinal breathing” and the formation of unstacked holes while base unpairing is suppressed. Flow linear dichroism in 20% diglyme indicates a 20 to 30% decrease in persistence length of DNA, supported by an increased flexibility in single-molecule nanochannel experiments in poly(ethylene glycol). A limited (3 to 6%) hyperchromicity but unaffected circular dichroism is consistent with transient unstacking events while maintaining an overall average B-DNA conformation. Further information about unstacking dynamics is obtained from the binding kinetics of large thread-intercalating ruthenium complexes, indicating that the hydrophobic effect provides a 10 to 100 times increased DNA unstacking frequency and an “open hole” population on the order of 10⁻² compared to 10⁻⁴ in normal aqueous solution. Spontaneous DNA strand exchange catalyzed by poly(ethylene glycol) makes us propose that hydrophobic residues in the L2 loop of recombination enzymes RecA and Rad51 may assist gene recombination via modulation of water activity near the DNA helix by hydrophobic interactions, in the manner described here. We speculate that such hydrophobic interactions may have catalytic roles also in other biological contexts, such as in polymerases.

Lessons learned in atomistic simulation of double-stranded DNA: Solvation and salt concerns [Article v1.0]

Nucleic acids are highly charged macromolecules sensitive to their surroundings of water, salt, and other biomolecules. Molecular dynamics simulations with accurate biomolecular force fields provide a detailed atomistic view into DNA and RNA that has been useful to study the structure and dynamics of these molecules and their biological relevance. In this work we study the Drew-Dickerson dodecamer duplex with the sequence d(GCGCAATTGCGC)₂ in three different salt concentrations and using different monvalent salt types to detect possible structural influence. Overall, the DNA shows no major structural changes regardless of amount or type of monovalent ions used. Our results show that only at very high salt conditions (5M) is a small structural effect observed in the DNA duplex, which mainly consist of narrowing of the grooves due to increased residence of ions. We also present the importance of sampling time to achieve a converged ensemble, which is of major relevance in any simulation to avoid biased or non-meaningful results.

Dehydrated DNA in B-form: ionic liquids in rescue

The functional B-conformation of DNA succumbs to the A-form at low water activity. Methods for room temperature DNA storage that rely upon ‘anhydrobiosis’, thus, often encounter the loss of DNA activity due to the B→A-DNA transition. Here, we show that ionic liquids, an emerging class of green solvents, can induce conformational transitions in DNA and even rescue the dehydrated DNA in the functional B-form. CD spectroscopic analyses not only reveal rapid transition of A-DNA in 78% ethanol medium to B-conformation in presence of ILs, but also the high resistance of IL-bound B-form to transit to A-DNA under dehydration. Molecular dynamics simulations show the unique ability of ILs to disrupt Na⁺ ion condensation and form ‘IL spine’ in DNA minor groove to drive the A→B transition. Implications of these findings range from the plausible use of ILs as novel anhydrobiotic DNA storage medium to a switch for modulating DNA conformational transitions.

Changes in Microenvironment Modulate the B- to A-DNA Transition

B- to A-DNA transition is known to be sensitive to the macroscopic properties of the solution, such as salt and ethanol concentrations. Microenvironmental effects on DNA conformational transition have been broadly studied. Providing an intuitive picture of how DNA responds to environmental changes is, however, still needed. Analyzing the chemical equilibrium of B-to-A DNA transition at critical concentrations, employing explicit-solvent simulations, is envisioned to help understand such microenvironmental effects. In the present study, free-energy calculations characterizing the B- to A-DNA transition and the distribution of cations were carried out in solvents with different ethanol concentrations. With the addition of ethanol, the most stable structure of DNA changes from the B- to A-form, in agreement with previous experimental observation. In 60% ethanol, a chemical equilibrium is found, showing reversible transition between B- and A-DNA. Analysis of the microenvironment around DNA suggests that with the increase of ethanol concentration, the cations exhibit a significant tendency to move toward the backbone, and mobility of water molecules around the major groove and backbone decreases gradually, leading eventually to a B-to-A transition. The present results provide a free-energy view of DNA microenvironment and of the role of cation motion in the conformational transition.

Role of protein interactions in stabilizing canonical DNA features in simulations of DNA in crowded environments

Background

Cellular environments are highly crowded with biological macromolecules resulting in frequent non-specific interactions. While the effect of such crowding on protein structure and dynamics has been studied extensively, very little is known how cellular crowding affects the conformational sampling of nucleic acids.

Results

The effect of protein crowding on the conformational preferences of DNA (deoxyribonucleic acid) is described from fully atomistic molecular dynamics simulations of systems containing a DNA dodecamer surrounded by protein crowders. From the simulations, it was found that DNA structures prefer to stay in B-like conformations in the presence of the crowders. The preference for B-like conformations results from non-specific interactions of crowder proteins with the DNA sugar-phosphate backbone. Moreover, the simulations suggest that the crowder interactions narrow the conformational sampling to canonical regions of the conformational space.

Conclusions

The overall conclusion is that crowding effects may stabilize the canonical features of DNA that are most important for biological function. The results are complementary to a previous study of DNA in reduced dielectric environments where reduced dielectric environments alone led to a conformational shift towards A-DNA. Such a shift was not observed here suggested that the reduced dielectric response of cellular environments is counteracted by non-specific interactions with protein crowders under in vivo conditions.

Electronic supplementary material

The online version of this article (10.1186/s13628-018-0048-y) contains supplementary material, which is available to authorized users.

Molecular Dynamics of Double Stranded Xylo-Nucleic Acid

Xylo-nucleic acid (XyloNA) is a synthetic analogue of ribo-nucleic acid (RNA), where the ribose sugar has been replaced by xylose. We present a molecular dynamics study of the conformational evolution of XyloNA double strand oligomers derived from A-RNA through the substitution of β-d-ribofuranose by β-d-xylofuranose and having lengths of 8, 16, and 29 base pairs, using a set of independent all-atom simulations performed at various time scales ranging from 55 to 100 ns, with one long 500 ns simulation of the 29-mer. In order to validate the robustness of XyloNA conformation, a set of simulations using various cutoff distances and solvation box dimensions has also been performed. These independent simulations reveal the uncoiling or elongation of the initial conformation to form an open ladder type transient state conformation and the subsequent formation of a highly flexible duplex with a tendency to coil in a left-handed fashion. The observed open ladder conformation is in line with recently obtained NMR data on the XyloNA 8-mer derived using 5'-d(GUGUACAC)-3'. The observed negative interbase pair twist leads to the observed highly flexible left-handed duplex, which is significantly less rigid than the stable left-handed dXyloNA duplex having a strong negative twist. A comparison between the xylo-analogues of DNA and RNA shows a clear distinction between the helical parameters, with implications for the pairing mechanism.

Cluster Plasmonics: Dielectric and Shape Effects on DNA-Stabilized Silver Clusters

This work investigates the effects of dielectric environment and cluster shape on electronic excitations of fluorescent DNA-stabilized silver clusters, AgN-DNA. We first establish that the longitudinal plasmon wavelengths predicted by classical Mie-Gans (MG) theory agree with previous quantum calculations for excitation wavelengths of linear silver atom chains, even for clusters of just a few atoms. Application of MG theory to AgN-DNA with 400-850 nm cluster excitation wavelengths indicates that these clusters are characterized by a collective excitation process and suggests effective cluster thicknesses of ∼2 silver atoms and aspect ratios of 1.5 to 5. To investigate sensitivity to the surrounding medium, we measure the wavelength shifts produced by addition of glycerol. These are smaller than reported for much larger gold nanoparticles but easily detectable due to narrower line widths, suggesting that AgN-DNA may have potential for fluorescence-reported changes in dielectric environment at length scales of ∼1 nm.

Computational Approach to Explore the B/A Junction Free Energy in DNA

Protein-DNA interactions induce conformational changes in DNA such as B- to A-form transitions at a local level. Such transitions are associated with a junction free energy cost at the boundary of two different conformations in a DNA molecule. In this study, we performed umbrella sampling simulations to find the free energy values of the B-A transition at the dinucleotide and trinucleotide level of DNA. Using a combination of dinucleotide and trinucleotide free energy costs obtained from simulations, we calculated the B/A junction free energy. Our study shows that the B/A junction free energy is 0.52 kcal mol(-1) for the A-philic GG step and 1.59 kcal mol(-1) for the B-philic AA step. This observation is in agreement with experimentally derived values. After excluding junction effects, we obtained an absolute free energy cost for the B- to A-form conversion for all the dinucleotide steps. These absolute free energies may be used for predicting the propensity of structural transitions in DNA.

Single-molecule portrait of DNA and RNA double helices

The composition and geometry of the genetic information carriers were described as double-stranded right helices sixty years ago. The flexibility of their sugar-phosphate backbones and the chemistry of their nucleotide subunits, which give rise to the RNA and DNA polymers, were soon reported to generate two main structural duplex states with biological relevance: the so-called A and B forms. Double-stranded (ds) RNA adopts the former whereas dsDNA is stable in the latter. The presence of flexural and torsional stresses in combination with environmental conditions in the cell or in the event of specific sequences in the genome can, however, stabilize other conformations. Single-molecule manipulation, besides affording the investigation of the elastic response of these polymers, can test the stability of their structural states and transition models. This approach is uniquely suited to understanding the basic features of protein binding molecules, the dynamics of molecular motors and to shedding more light on the biological relevance of the information blocks of life. Here, we provide a comprehensive single-molecule analysis of DNA and RNA double helices in the context of their structural polymorphism to set a rigorous interpretation of their material response both inside and outside the cell. From early knowledge of static structures to current dynamic investigations, we review their phase transitions and mechanochemical behaviour and harness this fundamental knowledge not only through biological sciences, but also for Nanotechnology and Nanomedicine.

Conformational preferences of DNA in reduced dielectric environments

The effect of reduced dielectric environments on the conformational sampling of DNA was examined through molecular dynamics simulations. Different dielectric environments were used to model one aspect of cellular environments. Implicit solvent based on the Generalized Born methodology was used to reflect different dielectric environments in the simulations. The simulation results show a tendency of DNA structures to favor noncanonical A-like conformations rather than canonical A- and B-forms as a result of the reduced dielectric environments. The results suggest that the reduced dielectric response in cellular environments may be sufficient to enhance the sampling of A-like DNA structures compared to dilute solvent conditions.

Balancing the interactions of ions, water, and DNA in the Drude polarizable force field

Recently we presented a first-generation all-atom Drude polarizable force field for DNA based on the classical Drude oscillator model, focusing on optimization of key dihedral angles followed by extensive validation of the force field parameters. Presently, we describe the procedure for balancing the electrostatic interactions between ions, water, and DNA as required for development of the Drude force field for DNA. The proper balance of these interactions is shown to impact DNA stability and subtler conformational properties, including the conformational equilibrium between the BI and BII states, and the A and B forms of DNA. The parametrization efforts were simultaneously guided by gas-phase quantum mechanics (QM) data on small model compounds and condensed-phase experimental data on the hydration and osmotic properties of biologically relevant ions and their solutions, as well as theoretical predictions for ionic distribution around DNA oligomer. In addition, fine-tuning of the internal base parameters was performed to obtain the final DNA model. Notably, the Drude model is shown to more accurately reproduce counterion condensation theory predictions of DNA charge neutralization by the condensed ions as compared to the CHARMM36 additive DNA force field, indicating an improved physical description of the forces dictating the ionic solvation of DNA due to the explicit treatment of electronic polarizability. In combination with the polarizable DNA force field, the availability of Drude polarizable parameters for proteins, lipids, and carbohydrates will allow for simulation studies of heterogeneous biological systems.

Sequence dependent free energy profiles of localized B- to A-form transition of DNA in water

DNA carries an inherent polymorphism, which surfaces under various external conditions. While B-form remains predominant under normal physiological conditions for most of the DNA sequences, low humidity and increased ion concentration cause B- to A-form transition. Certain proteins and molecules also sometimes cause local deformation of the DNA to the specific A-form. Previous experimental and computational studies focused on the overall B- to A-form transition. Here for the first time we investigated thermodynamics and mechanism of B- to A-form transition in water for various DNA sequences at a local dinucleotide base pair level. We introduced a new reaction coordinate Zp', based on the unique order parameter Zp, to drive B- to A-form transition locally and thereby calculate free energy profiles for the same for all the ten different dinucleotide steps embedded in a twelve base pair DNA. Results show that the trend of "A" and "B" philicity observed in experiment is preserved even at this local dinucleotide level, indicating its localized origin. Higher free energy cost obtained here is attributed to the cost of creating B∕A junctions along with formation of B->A transition at dimer level. We find that while water energetically stabilizes A-form for all the ten different dinucleotide steps to various extents, entropy acts against it. Therefore, we find that the stability of B-form DNA in water is entropic in origin. Mechanism of the conversion appears to be triggered by Slide; however, backbone parameters change concertedly.

Higher order structural effects stabilizing the reverse Watson-Crick Guanine-Cytosine base pair in functional RNAs

The G:C reverse Watson-Crick (W:W trans) base pair, also known as Levitt base pair in the context of tRNAs, is a structurally and functionally important base pair that contributes to tertiary interactions joining distant domains in functional RNA molecules and also participates in metabolite binding in riboswitches. We previously indicated that the isolated G:C W:W trans base pair is a rather unstable geometry, and that dicationic metal binding to the Guanine base or posttranscriptional modification of the Guanine can increase its stability. Herein, we extend our survey and report on other H-bonding interactions that can increase the stability of this base pair. To this aim, we performed a bioinformatics search of the PDB to locate all the occurencies of G:C trans base pairs. Interestingly, 66% of the G:C trans base pairs in the PDB are engaged in additional H-bonding interactions with other bases, the RNA backbone or structured water molecules. High level quantum mechanical calculations on a data set of representative crystal structures were performed to shed light on the structural stability and energetics of the various crystallographic motifs. This analysis was extended to the binding of the preQ1 metabolite to a preQ1-II riboswitch.

Molecular dynamics simulations of double-stranded DNA in an explicit solvent model with the zero-dipole summation method

Molecular dynamics (MD) simulations of a double-stranded DNA with explicit water and small ions were performed with the zero-dipole summation (ZD) method, which was recently developed as one of the non-Ewald methods. Double-stranded DNA is highly charged and polar, with phosphate groups in its backbone and their counterions, and thus precise treatment for the long-range electrostatic interactions is always required to maintain the stable and native double-stranded form. A simple truncation method deforms it profoundly. On the contrary, the ZD method, which considers the neutralities of charges and dipoles in a truncated subset, well reproduced the electrostatic energies of the DNA system calculated by the Ewald method. The MD simulations using the ZD method provided a stable DNA system, with similar structures and dynamic properties to those produced by the conventional Particle mesh Ewald method.

Biological physics in México: Review and new challenges

Biological and physical sciences possess a long-standing tradition of cooperativity as separate but related subfields of science. For some time, this cooperativity has been limited by their obvious differences in methods and views. Biological physics has recently experienced a kind of revival (or better a rebirth) due to the growth of molecular research on animate matter. New avenues for research have been opened for both theoretical and experimental physicists. Nevertheless, in order to better travel for such paths, the contemporary biological physicist should be armed with a set of specialized tools and methods but also with a new attitude toward multidisciplinarity. In this review article, we intend to somehow summarize what has been done in the past (in particular, as an example we will take a closer look at the Mexican case), to show some examples of fruitful investigations in the biological physics area and also to set a proposal of new curricula for physics students and professionals interested in applying their science to get a better understanding of the physical basis of biological function.

Sole and stable RNA duplexes of G-rich sequences located in the 5'-untranslated region of protooncogenes

Guanine- (G-) rich nucleic acid sequences can form four-stranded structures called G-quadruplexes. It is widely held that the formation of a G-quadruplex in RNA is more feasible than in DNA because of the lack of a complementary strand in mRNA. Here, we analyzed sequences of 5'-untranslated regions of protooncogenes and surprisingly found that these regions showed an enrichment of not only guanine (G) but also cytosine (C) nucleotides. Since neighboring cytosine- (C-) rich regions can affect the formation and stability of a G-quadruplex structure, we further investigated the properties of DNA and RNA structures of G-rich and GC-rich regions. We selected typical GC-rich RNA sequences from protooncogenes and corresponding DNA sequences and investigated their structures. It was found that the GC-rich RNA sequences formed stable A-form duplexes as their major structure independent of the surrounding conditions, including the presence of different cations (Na(+), K(+), or Li(+)) or molecular crowding with 40 wt % poly(ethylene glycol) with an average molecular mass of 200 Da although there are a few exceptions in which only a combination of K(+) and molecular crowding induced a G-quadruplex structure of an extremely G-rich RNA sequence. In contrast, structural polymorphisms involving duplexes, G-quadruplexes, and i-motifs were observed for GC-rich DNA sequences depending on the surrounding factors. These results demonstrate the considerable structural and functional differences in GC-rich sequences of the genome (DNA) and transcriptosome (mRNA) with respect to the nucleic acid backbone. Moreover, it was suggested that structural study for a G-rich RNA sequence should be carried out under cell-mimicking condition where K(+) and crowding cosolutes exist.

Helical structure of xylose-DNA

Synthetic biology and systems chemistry demonstrate a growing interest in modified nucleotides to achieve an enzymatically stable artificial nucleic acid. A potential candidate system is xylose-DNA, in which the 2'-deoxy-beta-D-ribo-furanose is substituted by 2'-deoxy-beta-D-xylo-furanose. We present here the helical structure and conformational analysis of xylose-DNA on the basis of 35 ns MD simulations of a 29-base-pair DNA duplex. Starting from a right-handed xylose-DNA helix, we observe a remarkable conformational transition from right- to left-handed helix. The left-handed xylose-DNA is highly dynamic, involving screwing and unscrewing motion of the helix. The sugar pucker induced helical changes influence the backbone to adopt the backbone angles for xylose-DNA while retaining the Watson-Crick base pairing and stacking interactions. The results demonstrate the chiral orthogonality of the ribose and xylose based episomes. As far as stability and compactness of information storage is concerned, the ribose based natural DNA is unsurpassed.

Spectroscopic and molecular dynamics evidence for a sequential mechanism for the A-to-B transition in DNA

The A-to-B form transition has been examined in three DNA duplexes, d(CGCGAATTCGCG)(2), d(CGCGAATTGCGC), and d(CGCAAATTTCGC), using circular dichroism spectroscopy, ultraviolet resonance Raman (UVRR) spectroscopy, and molecular dynamics (MD) simulation. Circular dichroism spectra confirm that these molecules adopt the A form under conditions of reduced water activity. UVRR results, obtained under similar conditions, suggest that the transition involves a series of intermediate forms between A and B. Cooperative and distinct transitions were observed for the bases and the sugars. Independent MD simulations on d(CGCGAATTCGCG)(2) show a spontaneous change from the A to B form in aqueous solution and describe a kinetic model that agrees well with UVRR results. Based on these observations, we predict that the mechanism of the transition involves a series of A/B hybrid forms and is sequential in nature, similar to previous crystallographic studies of derivatized duplexes. A simulation in which waters were restrained in the major groove of B DNA shows a rapid, spontaneous change from B to A at reduced water activity. These results indicate that a quasiergodic sampling of the solvent distribution may be a problem in going from B to A at reduced water activity in the course of an MD simulation.

Insights into structure, dynamics and hydration of locked nucleic acid (LNA) strand-based duplexes from molecular dynamics simulations

Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3′-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA–DNA, LNA–LNA, LNA–RNA, RNA–DNA and RNA–RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA–DNA and LNA–RNA duplexes are found to be similar to regular RNA–DNA and RNA–RNA duplexes, whereas the LNA–LNA duplex is found to have its helix partly unwound and does not resemble RNA–RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA–DNA and RNA–RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3′-endo conformations more efficiently than RNA. The LNA–LNA duplex has lesser backbone flexibility compared to the RNA–RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.

Comparison of protein force fields for molecular dynamics simulations

In the context of molecular dynamics simulations of proteins, the term "force field" refers to the combination of a mathematical formula and associated parameters that are used to describe the energy of the protein as a function of its atomic coordinates. In this review, we describe the functional forms and parameterization protocols of the widely used biomolecular force fields Amber, CHARMM, GROMOS, and OPLS-AA. We also summarize the ability of various readily available noncommercial molecular dynamics packages to perform simulations using these force fields, as well as to use modern methods for the generation of constant-temperature, constant-pressure ensembles and to treat long-range interactions. Finally, we finish with a discussion of the ability of these force fields to support the modeling of proteins in conjunction with nucleic acids, lipids, carbohydrates, and/or small molecules.

The role of electrostatics in the B to A transition of DNA: from solution to assembly

On the basis of a wealth of published experimental data and computer simulations, we build a simple physical model that allows us to rationalize the A to B transition of DNA in solution and in aggregates. In both cases we find that the electrostatic interactions are strong enough, alone, to induce the transition independently of other energetic contributions, e.g. those related to hydration. On the basis of this analysis we conclude that in ethanol/water mixtures, the effect responsible for the transition is the reduction of dielectric constant in the mixture. This is manifested in electrostatic self-energy terms that include the interaction of phosphate charges with condensed counterions. But in dense aggregates, electrostatics plays a dual role, giving rise to two competing effects. In the absence of groove localized counterions the electrostatic self-energy favours the B form, and the electrostatic interaction energy between neighbouring DNA favours the A form. However, the addition of enough counterions localized in the narrow groove reverses this. In dry aggregates of DNA both terms, in most cases, conspire to keep DNA in the A form. The analysis gives a broad picture of the B to A transition and sets a number of new research goals, particularly concerning simulations that may test our simple model for aggregates.

Nucleic acid solvation: from outside to insight

Nucleic acids are polyanionic molecules that were historically considered to be solely surrounded by a shell of water molecules and a neutralizing cloud of monovalent and divalent cations. In this respect, recent experimental and theoretical reports demonstrate that water molecules within complex nucleic acid structures can display very long residency times, and assist drug binding and catalytic reactions. Finally, anions can also bind to these polyanionic systems. Many of these recent insights are provided by state-of-the-art molecular dynamics simulations of nucleic acid systems, which will be described together with relevant methodological issues.

Theoretical study of large conformational transitions in DNA: the B<-->A conformational change in water and ethanol/water

We explore here the possibility of determining theoretically the free energy change associated with large conformational transitions in DNA, like the solvent-induced B<-->A conformational change. We find that a combination of targeted molecular dynamics (tMD) and the weighted histogram analysis method (WHAM) can be used to trace this transition in both water and ethanol/water mixture. The pathway of the transition in the A-->B direction mirrors the B-->A pathway, and is dominated by two processes that occur somewhat independently: local changes in sugar puckering and global rearrangements (particularly twist and roll) in the structure. The B-->A transition is found to be a quasi-harmonic process, which follows closely the first spontaneous deformation mode of B-DNA, showing that a physiologically-relevant deformation is in coded in the flexibility pattern of DNA.

The dynamics of the B-A transition of natural DNA double helices

The dynamics of the B-A transition of DNA double helices with different GC contents and various chain lengths has been characterized by an electric field pulse technique. The field-induced B-A reaction is separated from orientation effects using the magic angle technique. Amplitudes reflecting the B-A reaction are observed selectively in the limited range of ethanol contents, where CD spectra demonstrate the B-A transition. The maximum amplitude appears at 1-2% higher ethanol content than the center of the B-A transition observed by CD because electric field pulses induce a relatively large perturbation from the A- toward the B-form. The relaxation curves measured after pulse termination reflect a spectrum of up to three relaxation processes. For DNA's with approximately 50% GC, the main part of the amplitude ( approximately 75%) is associated with time constants of approximately 2 micros, and another major component appears with time constants of 50-100 micros. These relaxation effects have been observed for DNA samples with 859, 2629, 7160, and 48501 bp. The time constant associated with the main amplitude increases with decreasing GC content from approximately 2 micros at 50% GC to approximately 3 mus at 41% GC and approximately 10 micros at 0% GC at the center of the B-A transition. Model calculations on the kinetics of cooperative linear Ising lattices predict the appearance of a distinct maximum of the mean relaxation time at the center of the transition. The absence of such maximum in our experimental data indicates a low cooperativity of the B-A transition with a nucleation parameter of approximately 0.1. The rate of the B-A transition is lower by approximately 3 orders of magnitude than that predicted by molecular dynamics simulations.

Hydration and conformational transitions in DNA, RNA, and mixed DNA-RNA triplexes studied by gravimetry and FTIR spectroscopy

We have studied by gravimetric measurements and FTIR spectroscopy the hydration of duplexes and triplexes formed by combinations of dA(n), dT(n), rA(n), and rU(n) strands. Results obtained on hydrated films show important differences in their hydration and in the structural transitions which can be induced by varying the water content of the samples. The number of water molecules per nucleotide (w/n) measured at high relative humidity (98% R.H.) is found to be 21 for dA(n).dT(n) and 15 for rA(n).rU(n). Addition of a third rU(n) strand does not change the number of water molecules per nucleotide: w/n=21 for rU(n)*dA(n).dT(n) and w/n=15 for rU(n)*rA(n).rU(n). On the contrary, the addition of a third dT(n) strand changes the water content but in a different way, depending whether the duplex is DNA or RNA. Thus, a loss of four water molecules per nucleotide is measured for dT(n)*dA(n).dT(n) while an increase of two water molecules per nucleotide is observed for dT(n)*rA(n).rU(n). The final hydration is the same for both triplexes (w/n=17). The desorption profiles obtained by gravimetry and FTIR spectroscopy are similar for the rA(n).rU(n) duplex and the rU(n)*rA(n).rU(n) triplex. On the contrary, the desorption profiles of the dA(n).dT(n) duplex and the triplexes formed with it (rU(n)*dA(n).dT(n) and dT(n)*dA(n).dT(n)) are different from each other. This is correlated with conformational transitions induced by varying the hydration content of the different structures, as shown by FTIR spectroscopy. Modifications of the phosphate group hydration and of the sugar conformation (S to N type repuckering) induced by decrease of the water content are observed in the case of triplexes formed on the dA(n).dT(n) duplex.