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

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.

Variability of the DNA Backbone Geometry in DNA-Protein Complexes: Experimental Data Analysis

We have analyzed and compared the available experimental data (PDB) on the backbone geometry of the DNA in solution (NMR), in crystals (X-rays), and in complexes with proteins (X-rays and cryo-electron microscopy). The deoxyribose (pseudorotational angle τ₀) and ε/ζ (BI-BII transition in phosphates) flexibilities are practically the same in the four samples. The α/γ mobility is minimal in crystalline DNA: on the histograms, there is one canonical and one noncanonical t/t peak. The α/γ mobility increases in DNA solutions (three more noncanonical peaks) and is maximal in DNA-protein complexes (another additional peak). On a large amount of data, we have confirmed that the three main degrees of freedom of the sugar-phosphate backbone are "orthogonal": changes in any of the angles τ₀, (ζ-ε), and (γ-α) occur, as a rule, at a constant (usually canonical) value of any other. In the DNA-protein complexes, none of the geometrical parameters commonly used to distinguish the A and B forms of DNA, except for Zp and its simpler analog Zp', show an unambiguous correlation with τ₀. Proteins, binding to DNA, in 59% of cases change the local shape of the helix up to the characteristic of the A-form without switching the deoxyribose conformation from south to north. However, we have found simple local characteristics of one nucleotide that correlate with the angles τ₀ and (ζ-ε). These are the angles C3'C1'N* and C4'C3'P(2), respectively. They are orthogonal in DNA-protein complexes exactly as the pair τ₀ and (ζ-ε). Most characteristics of DNA in complexes with proteins are the same in X-ray and in cryo-EM data, except for the histogram for the angle τ₀. We offer a possible explanation for this difference. We also discuss the artifacts on the ε/ζ histogram for DNA in solutions caused by the currently used NMR refinement protocols.

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.

Multiple plasmid origin-of-transfer regions might aid the spread of antimicrobial resistance to human pathogens

Antimicrobial resistance poses a great danger to humanity, in part due to the widespread horizontal gene transfer of plasmids via conjugation. Modeling of plasmid transfer is essential to uncovering the fundamentals of resistance transfer and for the development of predictive measures to limit the spread of resistance. However, a major limitation in the current understanding of plasmids is the incomplete characterization of the conjugative DNA transfer mechanisms, which conceals the actual potential for plasmid transfer in nature. Here, we consider that the plasmid-borne origin-of-transfer substrates encode specific DNA structural properties that can facilitate finding these regions in large datasets and develop a DNA structure-based alignment procedure for typing the transfer substrates that outperforms sequence-based approaches. Thousands of putative DNA transfer substrates are identified, showing that plasmid mobility can be twofold higher and span almost twofold more host species than is currently known. Over half of all putative mobile plasmids contain the means for mobilization by conjugation systems belonging to different mobility groups, which can hypothetically link previously confined host ranges across ecological habitats into a robust plasmid transfer network. This hypothetical network is found to facilitate the transfer of antimicrobial resistance from environmental genetic reservoirs to human pathogens, which might be an important driver of the observed rapid resistance development in humans and thus an important point of focus for future prevention measures.

Concerted dynamics of metallo-base pairs in an A/B-form helical transition

Metal-mediated base pairs expand the repertoire of nucleic acid structures and dynamics. Here we report solution structures and dynamics of duplex DNA containing two all-natural C-HgII-T metallo base pairs separated by six canonical base pairs. NMR experiments reveal a 3:1 ratio of well-resolved structures in dynamic equilibrium. The major species contains two (N3)T-HgII-(N3)C base pairs in a predominantly B-form helix. The minor species contains (N3)T-HgII-(N4)C base pairs and greater A-form characteristics. Ten-fold different 1J coupling constants (15N,199Hg) are observed for (N3)C-HgII (114 Hz) versus (N4)C-HgII (1052 Hz) connectivities, reflecting differences in cytosine ionization and metal-bonding strengths. Dynamic interconversion between the two types of C-HgII-T base pairs are coupled to a global conformational exchange between the helices. These observations inspired the design of a repetitive DNA sequence capable of undergoing a global B-to-A-form helical transition upon adding HgII, demonstrating that C-HgII-T has unique switching potential in DNA-based materials and devices.

Free-Energy Profiles for A-/B-DNA Conformational Transitions in Isolated and Aggregated States from All-Atom Molecular Dynamics Simulation

In ordinary aqueous solution, B-DNA is the major structural form of DNA. After the addition of ethanol, DNA is thought to be aggregated/condensed in the A-form structure. However, there is uncertainty as to whether the B-to-A conformational change is connected to the aggregation/condensation steps. In this study, we performed all-atom molecular dynamics simulations and calculated the free-energy surface involved in the A/B conformational transition for isolated and aggregated Dickerson-Drew dodecamers (DDDs) in water and 85% ethanol environments. We found in the case of an isolated DDD, the overall free-energy profile is entirely downhill to give the B-DNA conformation in both water and 85% ethanol. However, in the aggregated state and 85% ethanol environment, there is a free-energy minimum associated with the A-DNA region in addition to the global B-DNA minimum, and there is a ∼3 kcal/mol free-energy barrier to the A-to-B conformational change. The molecular dynamics results suggest that aggregation of DNA is essential for forming A-DNA.

AMOEBA Polarizable Atomic Multipole Force Field for Nucleic Acids

The AMOEBA polarizable atomic multipole force field for nucleic acids is presented. Valence and electrostatic parameters were determined from high-level quantum mechanical data, including structures, conformational energy, and electrostatic potentials, of nucleotide model compounds. Previously derived parameters for the phosphate group and nucleobases were incorporated. A total of over 35 μs of condensed-phase molecular dynamics simulations of DNA and RNA molecules in aqueous solution and crystal lattice were performed to validate and refine the force field. The solution and/or crystal structures of DNA B-form duplexes, RNA duplexes, and hairpins were captured with an average root-mean-squared deviation from NMR structures below or around 2.0 Å. Structural details, such as base pairing and stacking, sugar puckering, backbone and χ-torsion angles, groove geometries, and crystal packing interfaces, agreed well with NMR and/or X-ray. The interconversion between A- and B-form DNAs was observed in ethanol-water mixtures at 328 K. Crystal lattices of B- and Z-form DNA and A-form RNA were examined with simulations. For the RNA tetraloop, single strand tetramers, and HIV TAR with 29 residues, the simulated conformational states, 3 J-coupling, nuclear Overhauser effect, and residual dipolar coupling data were compared with NMR results. Starting from a totally unstacked/unfolding state, the rCAAU tetranucleotide was folded into A-form-like structures during ∼1 μs molecular dynamics simulations.

Effect of water and ionic liquids on biomolecules

The remarkable progress in the field of ionic liquids (ILs) in the last two decades has involved investigations on different aspects of ILs in various conditions. The nontoxic and biocompatible nature of ILs makes them a suitable substance for the storage and application of biomolecules. In this regard, the aqueous IL solutions have attracted a large number of studies to comprehend the role of water in modulating various properties of biomolecules. Here, we review some of the recent studies on aqueous ILs that concern the role of water in altering the behavior of ILs in general and in case of biomolecules solvated in ILs. The different structural and dynamic effects caused by water have been highlighted. We discuss the different modes of IL interaction that are responsible for stabilization and destabilization of proteins and enzymes followed by examples of water effect on this. The role of water in the case of nucleic acid storage in ILs, an area which has mostly been underrated, also has been emphasized. Our discussions highlight the fact that the effects of water on IL behavior are not general and are highly dependent on the nature of the IL under consideration. Overall, we aim to draw attention to the significance of water dynamics in the aqueous IL solutions, a better understanding of which can help in developing superior storage materials for application purposes.

Comparison of the Binding Geometry of Free-Base and Hexacoordinated Cationic Porphyrins to A- and B-Form DNA

Although the transition from B-DNA to the A-form is essential for many biological concerns, the properties of this transition have not been resolved. The B to A equilibrium can be analyzed conveniently because of the significant changes in circular dichroism (CD) and absorption spectrum. CD and linear dichroism (LD) methods were used to examine the binding of water-soluble meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (TMPyP) and its derivatives, Co-TMPyP, with B- and A-calf thymus DNA. B- to A-transitions occurred when the physiological buffer was replaced with a water-ethanol mixture (∼80 v/v %), and the fluorescence emission spectra of TMPyP bound to DNA showed a different pattern under ethanol-water conditions and water alone. The featureless broad emission bands of TMPyP were split into two peaks near at 658 and 715 nm in the presence of DNA under an aqueous solution. In the case of an ethanol-water system, however, the emission bands are split in two peaks near at 648 and 708 nm and 656 and 715 nm with and without DNA, respectively. This may be due to a change in the solution polarity. On the basis of the CD and LD data, TMPyP interacts with B-DNA via intercalation at a low ratio under a low ionic strength, 1 mM sodium phosphate. On the other hand, the interaction with A-DNA (80 v/v % ethanol-water system) occurs in a nonintercalating manner. This difference might be because the structural conformations, such as the groove of A-DNA, are not as deep as in B-DNA and the bases are much more tilted. In the case of Co-TMPyP, porphyrin binds preferably via an outside self-stacking mode with B- and A-DNA.

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.