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

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.

Synergy of Mutation-Induced Effects in Human Vitamin K Epoxide Reductase: Perspectives and Challenges for Allo-Network Modulator Design

The human Vitamin K Epoxide Reductase Complex (hVKORC1), a key enzyme transforming vitamin K into the form necessary for blood clotting, requires for its activation the reducing equivalents delivered by its redox partner through thiol-disulfide exchange reactions. The luminal loop (L-loop) is the principal mediator of hVKORC1 activation, and it is a region frequently harbouring numerous missense mutations. Four L-loop hVKORC1 mutants, suggested in vitro as either resistant (A41S, H68Y) or completely inactive (S52W, W59R), were studied in the oxidised state by numerical approaches (in silico). The DYNASOME and POCKETOME of each mutant were characterised and compared to the native protein, recently described as a modular protein composed of the structurally stable transmembrane domain (TMD) and the intrinsically disordered L-loop, exhibiting quasi-independent dynamics. The DYNASOME of mutants revealed that L-loop missense point mutations impact not only its folding and dynamics, but also those of the TMD, highlighting a strong mutation-specific interdependence between these domains. Another consequence of the mutation-induced effects manifests in the global changes (geometric, topological, and probabilistic) of the newly detected cryptic pockets and the alternation of the recognition properties of the L-loop with its redox protein. Based on our results, we postulate that (i) intra-protein allosteric regulation and (ii) the inherent allosteric regulation and cryptic pockets of each mutant depend on its DYNASOME; and (iii) the recognition of the redox protein by hVKORC1 (INTERACTOME) depend on their DYNASOME. This multifaceted description of proteins produces "omics" data sets, crucial for understanding the physiological processes of proteins and the pathologies caused by alteration of the protein properties at various "omics" levels. Additionally, such characterisation opens novel perspectives for the development of "allo-network drugs" essential for the treatment of blood disorders.

Effect of Soybean Protein Concentrate Preparation on Copy Numbers and Structural Characteristics of DNA from Genetically Modified Soybean

To regulate the degradation of transgenic DNA and lay theoretical foundations for the rational utilization of genetically modified (GM) products, variations in copy numbers and structural characteristics of DNA from GM soybean event GTS 40-3-2 during soybean protein concentrate (SPC) preparation were evaluated. Results showed that defatting and the first ethanol extraction were key procedures inducing DNA degradation. After these two procedures, copy numbers of the lectin and cp4 epsps targets decreased by more than 4 × 10⁸, occupying 36.88-49.30% of the total copy numbers from raw soybean. Atomic force microscopy images visually revealed the degradation of DNA that thinned and shortened during SPC preparation. Circular dichroism spectra suggested a lower helicity of DNA from defatted soybean kernel flour and a conformation transition of DNA from B-type to A-type after ethanol extraction. The fluorescence intensity of DNA decreased during SPC preparation, verifying the DNA damage along this preparation chain.

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune "bond" lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol's role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA's conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Overall, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

Folding and Intrinsic Disorder of the Receptor Tyrosine Kinase KIT Insert Domain Seen by Conventional Molecular Dynamics Simulations

The kinase insert domain (KID) of RTK KIT is the key recruitment region for downstream signalling proteins. KID, studied by molecular dynamics simulations as a cleaved polypeptide and as a native domain fused to KIT, showed intrinsic disorder represented by a set of heterogeneous conformations. The accurate atomistic models showed that the helical fold of KID is mainly sequence dependent. However, the reduced fold of the native KID suggests that its folding is allosterically controlled by the kinase domain. The tertiary structure of KID represents a compact array of highly variable α- and 3₁₀-helices linked by flexible loops playing a principal role in the conformational diversity. The helically folded KID retains a collapsed globule-like shape due to non-covalent interactions associated in a ternary hydrophobic core. The free energy landscapes constructed from first principles-the size, the measure of the average distance between the conformations, the amount of helices and the solvent-accessible surface area-describe the KID disorder through a collection of minima (wells), providing a direct evaluation of conformational ensembles. We found that the cleaved KID simulated with restricted N- and C-ends better reproduces the native KID than the isolated polypeptide. We suggest that a cyclic, generic KID would be best suited for future studies of KID f post-transduction effects.

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.

Molecular dynamics study of the recognition of ATP by nucleic acid aptamers

Despite their great success in recognizing small molecules in vitro, nucleic acid aptamers are rarely used in clinical settings. This is partially due to the lack of structure-based mechanistic information. In this work, atomistic molecular dynamics simulations are used to study the static and dynamic supramolecular structures relevant to the process of the wild-type (wt) nucleic acid aptamer recognition and binding of ATP. The effects brought about by mutation of key residues in the recognition site are also explored. The simulations reveal that the aptamer displays a high degree of rigidity and is structurally very little affected by the binding of ATP. Interaction energy decomposition shows that dispersion forces from π-stacking between ATP and the G6 and A23 nucleobases in the aptamer binding site plays a more important role in stabilizing the supramolecular complex, compared to hydrogen-bond interaction between ATP and G22. Moreover, metadynamics simulations show that during the association process, water molecules act as essential bridges connecting ATP with G22, which favors the dynamic stability of the complex. The calculations carried out on three mutated aptamer structures confirm the crucial role of the hydrogen bonds and π-stacking interactions for the binding affinity of the ATP nucleic acid aptamer.

Nanostructured Polyelectrolyte Complexes Based on Water-Soluble Thiacalix[4]Arene and Pillar[5]Arene: Self-Assembly in Micelleplexes and Polyplexes at Packaging DNA

Controlling the self-assembly of polyfunctional compounds in interpolyelectrolyte aggregates is an extremely challenging task. The use of macrocyclic compounds offers new opportunities in design of a new generation of mixed nanoparticles. This approach allows creating aggregates with multivalent molecular recognition, improved binding efficiency and selectivity. In this paper, we reported a straightforward approach to the synthesis of interpolyelectrolytes by co-assembling of the thiacalix[4]arene with four negatively charged functional groups on the one side of macrocycle, and pillar[5]arene with 10 ammonium groups located on both sides. Nanostructured polyelectrolyte complexes show effective packaging of high-molecular DNA from calf thymus. The interaction of co-interpolyelectrolytes with the DNA is completely different from the interaction of the pillar[5]arene with the DNA. Two different complexes with DNA, i.e., micelleplex- and polyplex-type, were formed. The DNA in both cases preserved its secondary structure in native B form without distorting helicity. The presented approach provides important advantage for the design of effective biomolecular gene delivery systems.

Enhanced sampling and free energy calculations for protein simulations

Molecular dynamics simulation is a powerful computational technique to study biomolecular systems, which complements experiments by providing insights into the structural dynamics relevant to biological functions at atomic scale. It can also be used to calculate the free energy landscapes of the conformational transitions to better understand the functions of the biomolecules. However, the sampling of biomolecular configurations is limited by the free energy barriers that need to be overcome, leading to considerable gaps between the timescales reached by MD simulation and those governing biological processes. To address this issue, many enhanced sampling methodologies have been developed to increase the sampling efficiency of molecular dynamics simulations and free energy calculations. Usually, enhanced sampling algorithms can be classified into methods based on collective variables (CV-based) and approaches which do not require predefined CVs (CV-free). In this chapter, the theoretical basis of free energy estimation is briefly reviewed first, followed by the reviews of the most common CV-based and CV-free methods including the presentation of some examples and recent developments. Finally, the combination of different enhanced sampling methods is discussed.

Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes

Association of peripheral proteins with lipid bilayers regulates membrane signaling and dynamics. Pleckstrin homology (PH) domains bind to phosphatidylinositol phosphate (PIP) molecules in membranes. The effects of local PIP enrichment on the interaction of PH domains with membranes is unclear. Molecular dynamics simulations allow estimation of the binding energy of GRP1 PH domain to PIP3-containing membranes. The free energy of interaction of the PH domain with more than two PIP3 molecules is comparable to experimental values, suggesting that PH domain binding involves local clustering of PIP molecules within membranes. We describe a mechanism of PH binding proceeding via an encounter state to two bound states which differ in the orientation of the protein relative to the membrane, these orientations depending on the local PIP concentration. These results suggest that nanoscale clustering of PIP molecules can control the strength and orientation of PH domain interaction in a concentration-dependent manner.

Computational Studies of Catalytic Loop Dynamics in Yersinia Protein Tyrosine Phosphatase Using Pathway Optimization Methods

Yersinia Protein Tyrosine Phosphatase (YopH) is the most efficient enzyme among all known PTPases and relies on its catalytic loop movements for substrate binding and catalysis. Fluorescence, NMR, and UV resonance Raman (UVRR) techniques have been used to study the thermodynamic and dynamic properties of the loop motions. In this study, a computational approach based on the pathway refinement methods nudged elastic band (NEB) and harmonic Fourier beads (HFB) has been developed to provide structural interpretations for the experimentally observed kinetic processes. In this approach, the minimum potential energy pathways for the loop open/closure conformational changes were determined by NEB using a one-dimensional global coordinate. Two dimensional data analyses of the NEB results were performed as an efficient method to qualitatively evaluate the energetics of transitions along several specific physical coordinates. The free energy barriers for these transitions were then determined more precisely using the HFB method. Kinetic parameters were estimated from the energy barriers using transition state theory and compared against experimentally determined kinetic parameters. When the calculated energy barriers are calibrated by a simple "scaling factor", as have been done in our previous vibrational frequency calculations to explain the ligand frequency shift upon its binding to protein, it is possible to make structural interpretations of several observed enzyme dynamic rates. For example, the nanosecond kinetics observed by fluorescence anisotropy may be assigned to the translational motion of the catalytic loop and microsecond kinetics observed in fluorescence T-jump can be assigned to the loop backbone dihedral angle flipping. Furthermore, we can predict that a Trp354 conformational conversion associated with the loop movements would occur on the tens of nanoseconds time scale, to be verified by future UVRR T-jump studies.

Determining selection free energetics from nucleotide pre-insertion to insertion in viral T7 RNA polymerase transcription fidelity control

An elongation cycle of a transcribing RNA polymerase (RNAP) usually consists of multiple kinetics steps, so there exist multiple kinetic checkpoints where non-cognate nucleotides can be selected against. We conducted comprehensive free energy calculations on various nucleotide insertions for viral T7 RNAP employing all-atom molecular dynamics simulations. By comparing insertion free energy profiles between the non-cognate nucleotide species (rGTP and dATP) and a cognate one (rATP), we obtained selection free energetics from the nucleotide pre-insertion to the insertion checkpoints, and further inferred the selection energetics down to the catalytic stage. We find that the insertion of base mismatch rGTP proceeds mainly through an off-path along which both pre-insertion screening and insertion inhibition play significant roles. In comparison, the selection against dATP is found to go through an off-path pre-insertion screening along with an on-path insertion inhibition. Interestingly, we notice that two magnesium ions switch roles of leave and stay during the dATP on-path insertion. Finally, we infer that substantial selection energetic is still required to catalytically inhibit the mismatched rGTP to achieve an elongation error rate ∼10-4 or lower; while no catalytic selection seems to be further needed against dATP to obtain an error rate ∼10-2.

Long Time-Scale Atomistic Simulations of the Structure and Dynamics of Transcription Factor-DNA Recognition

Recent years have witnessed an explosion of interest in computational studies of DNA binding proteins, including both coarse-grained and atomistic simulations of transcription factor-DNA recognition, to understand how these transcription factors recognize their binding sites on the DNA with such exquisite specificity. The present study performs microsecond time scale all-atom simulations of the dimeric form of the lactose repressor (LacI), both in the absence of any DNA and in the presence of both specific and nonspecific complexes, considering three different DNA sequences. We examine, specifically, the conformational differences between specific and nonspecific protein-DNA interactions, as well as the behavior of the helix-turn-helix motif of LacI when interacting with the DNA. Our simulations suggest that stable LacI binding occurs primarily to bent A-form DNA, with a loss of LacI conformational entropy and optimization of correlated conformational equilibria across the protein. In addition, binding to the specific operator sequence involves a slightly larger number of stabilizing DNA-protein hydrogen bonds (in comparison to nonspecific complexes), which may account for the experimentally observed specificity for this operator. In doing so, our simulations provide a detailed atomistic description of potential structural drivers for LacI selectivity.

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.

Conformational changes of DNA induced by a trans-azobenzene derivative via non-covalent interactions

In biological environments and in aqueous solution, DNA generally adopts the canonical B conformation. Recently, an azobenzene photoswitch containing a polyamine chain with three positive charges was shown to induce a reversible conformational transition between the A and B forms of DNA, the transition being triggered by trans-cis isomerization of the photoswitch upon non-covalent intercalation. It was proposed that, in its trans conformation, azobenzene stabilizes the A form of DNA. The structural details and the mechanism upon which trans-azobenzene induces the B-to-A DNA transition remain, however, unclear. In the present work, two possible intercalating modes of trans-azobenzene, from the minor groove and from the major groove, were investigated with all-atom molecular-dynamics simulations. Intercalation from the major groove was found to be the most probable binding mode due to favorable electrostatic and π-π stacking interactions. The free-energy profile associated with the B-to-A conformational transition reveals that intercalation from the major groove leads to a conformational change of DNA, showing a slight tendency to interconvert from B- to A-DNA, in agreement with the CD spectrum obtained from the experiment. However, the presence of only one interacting azobenzene is not sufficient to lead to a global conformational change to A-DNA. The present results are expected to serve in the design of DNA switches, which can induce reversible DNA conformational changes.

Probing helical transitions in a DNA duplex

The complex conformational change from B-DNA to Z-DNA requires inversion of helix-handedness. Multiple degrees of freedom are intricately coupled during this transition, and formulating an appropriate reaction coordinate that captures the underlying complexity would be problematic. In this contribution, we adopt an alternative approach, based on the potential energy landscape perspective, to construct a kinetic transition network. Microscopic insight into the B → Z transition is provided in terms of geometrically defined discrete paths consisting of local minima and the transition states that connect them. We find that the inversion of handedness can occur via two competing mechanisms, either involving stretched intermediates, or a B-Z junction, in agreement with previous predictions. The organisation of the free energy landscape further suggests that this process is likely to be slow under physiological conditions. Our results represent a key step towards decoding the more intriguing features of the B → Z transition, such as the role of ionic strength and negative supercoiling in reshaping the landscape.

Mechanism of mRNA-STAR domain interaction: Molecular dynamics simulations of Mammalian Quaking STAR protein

STAR proteins are evolutionary conserved mRNA-binding proteins that post-transcriptionally regulate gene expression at all stages of RNA metabolism. These proteins possess conserved STAR domain that recognizes identical RNA regulatory elements as YUAAY. Recently reported crystal structures show that STAR domain is composed of N-terminal QUA1, K-homology domain (KH) and C-terminal QUA2, and mRNA binding is mediated by KH-QUA2 domain. Here, we present simulation studies done to investigate binding of mRNA to STAR protein, mammalian Quaking protein (QKI). We carried out conventional MD simulations of STAR domain in presence and absence of mRNA, and studied the impact of mRNA on the stability, dynamics and underlying allosteric mechanism of STAR domain. Our unbiased simulations results show that presence of mRNA stabilizes the overall STAR domain by reducing the structural deviations, correlating the ‘within-domain’ motions, and maintaining the native contacts information. Absence of mRNA not only influenced the essential modes of motion of STAR domain, but also affected the connectivity of networks within STAR domain. We further explored the dissociation of mRNA from STAR domain using umbrella sampling simulations, and the results suggest that mRNA binding to STAR domain occurs in multi-step: first conformational selection of mRNA backbone conformations, followed by induced fit mechanism as nucleobases interact with STAR domain.

Electric-Field-Induced Protein Translocation via a Conformational Transition in SecDF: An MD Study

SecDF is an important component of the Sec protein translocation machinery embedded in the bacterial membrane, which is associated with many functions, such as stabilizing other Sec translocon components within the membrane, maintaining the transmembrane (TM) potential, and facilitating the ATP-independent stage of the translocation mechanism. Related studies suggest that SecDF undergoes functionally important conformational changes that involve mainly its P1-head domain and that these changes are coupled with the proton motive force (Δp). However, there still is not a clear understanding of how SecDF functions, its exact role in the translocation machinery, and how its function is related to Δp. Here, using all-atom molecular dynamics simulations combined with umbrella sampling, we study the P1-head conformational change and how it is coupled to the proton motive force. We report potentials of mean force along a root-mean-square-distance-based reaction coordinate obtained in the presence and absence of the TM electrical potential. Our results show that the interaction of the P1 domain dipole moment with the TM electrical field considerably lowers the free-energy barrier in the direction of F-form to I-form transition.

Hydration forces between aligned DNA helices undergoing B to A conformational change: In-situ X-ray fiber diffraction studies in a humidity and temperature controlled environment

Hydration forces between DNA molecules in the A- and B-Form were studied using a newly developed technique enabling simultaneous in situ control of temperature and relative humidity. X-ray diffraction data were collected from oriented calf-thymus DNA fibers in the relative humidity range of 98%-70%, during which DNA undergoes the B- to A-form transition. Coexistence of both forms was observed over a finite humidity range at the transition. The change in DNA separation in response to variation in humidity, i.e. change of chemical potential, led to the derivation of a force-distance curve with a characteristic exponential decay constant of∼2Å for both A- and B-DNA. While previous osmotic stress measurements had yielded similar force-decay constants, they were limited to B-DNA with a surface separation (wall-to-wall distance) typically>5Å. The current investigation confirms that the hydration force remains dominant even in the dry A-DNA state and at surface separation down to∼1.5Å, within the first hydration shell. It is shown that the observed chemical potential difference between the A and B states could be attributed to the water layer inside the major and minor grooves of the A-DNA double helices, which can partially interpenetrate each other in the tightly packed A phase. The humidity-controlled X-ray diffraction method described here can be employed to perform direct force measurements on a broad range of biological structures such as membranes and filamentous protein networks.

Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 μs molecular dynamics simulations

Z

Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 μs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.

Geometric Patterns for Neighboring Bases Near the Stacked State in Nucleic Acid Strands

Structural variation in base stacking has been analyzed frequently in isolated double helical contexts for nucleic acids, but not as often in nonhelical geometries or in complex biomolecular environments. In this study, conformations of two neighboring bases near their stacked state in any environment are comprehensively characterized for single-strand dinucleotide (SSD) nucleic acid crystal structure conformations. An ensemble clustering method is used to identify a reduced set of representative stacking geometries based on pairwise distances between select atoms in consecutive bases, with multiple separable conformational clusters obtained for categories divided by nucleic acid type (DNA/RNA), SSD sequence, stacking face orientation, and the presence or absence of a protein environment. For both DNA and RNA, SSD conformations are observed that are either close to the A-form, or close to the B-form, or intermediate between the two forms, or further away from either form, illustrating the local structural heterogeneity near the stacked state. Among this large variety of distinct conformations, several common stacking patterns are observed between DNA and RNA, and between nucleic acids in isolation or in complex with proteins, suggesting that these might be stable stacking orientations. Noncanonical face/face orientations of the two bases are also observed for neighboring bases in the same strand, but their frequency is much lower, with multiple SSD sequences across categories showing no occurrences of such unusual stacked conformations. The resulting reduced set of stacking geometries is directly useful for stacking-energy comparisons between empirical force fields, prediction of plausible localized variations in single-strand structures near their canonical states, and identification of analogous stacking patterns in newly solved nucleic acid containing structures.

Exploring Flexibility of Progesterone Receptor Ligand Binding Domain Using Molecular Dynamics

Progesterone receptor (PR), a member of nuclear receptor (NR) superfamily, plays a vital role for female reproductive tissue development, differentiation and maintenance. PR ligand, such as progesterone, induces conformation changes in PR ligand binding domain (LBD), thus mediates subsequent gene regulation cascades. PR LBD may adopt different conformations upon an agonist or an antagonist binding. These different conformations would trigger distinct transcription events. Therefore, the dynamics of PR LBD would be of general interest to biologists for a deep understanding of its structure-function relationship. However, no apo-form (non-ligand bound) of PR LBD model has been proposed either by experiments or computational methods so far. In this study, we explored the structural dynamics of PR LBD using molecular dynamics simulations and advanced sampling tools in both ligand-bound and the apo-forms. Resolved by the simulation study, helix 11, helix 12 and loop 895–908 (the loop between these two helices) are quite flexible in antagonistic conformation. Several residues, such as Arg899 and Glu723, could form salt-bridging interaction between helix 11 and helix 3, and are important for the PR LBD dynamics. And we also propose that helix 12 in apo-form PR LBD, not like other NR LBDs, such as human estrogen receptor α (ERα) LBD, may not adopt a totally extended conformation. With the aid of umbrella sampling and metadynamics simulations, several stable conformations of apo-form PR LBD have been sampled, which may work as critical structural models for further large scale virtual screening study to discover novel PR ligands for therapeutic application.

Study on the Application of the Combination of TMD Simulation and Umbrella Sampling in PMF Calculation for Molecular Conformational Transitions

Free energy calculations of the potential of mean force (PMF) based on the combination of targeted molecular dynamics (TMD) simulations and umbrella samplings as a function of physical coordinates have been applied to explore the detailed pathways and the corresponding free energy profiles for the conformational transition processes of the butane molecule and the 35-residue villin headpiece subdomain (HP35). The accurate PMF profiles for describing the dihedral rotation of butane under both coordinates of dihedral rotation and root mean square deviation (RMSD) variation were obtained based on the different umbrella samplings from the same TMD simulations. The initial structures for the umbrella samplings can be conveniently selected from the TMD trajectories. For the application of this computational method in the unfolding process of the HP35 protein, the PMF calculation along with the coordinate of the radius of gyration (R_g) presents the gradual increase of free energies by about 1 kcal/mol with the energy fluctuations. The feature of conformational transition for the unfolding process of the HP35 protein shows that the spherical structure extends and the middle α-helix unfolds firstly, followed by the unfolding of other α-helices. The computational method for the PMF calculations based on the combination of TMD simulations and umbrella samplings provided a valuable strategy in investigating detailed conformational transition pathways for other allosteric processes.

Structural Flexibility and Conformation Features of Cyclic Dinucleotides in Aqueous Solutions

Cyclic dinucleotides are able to trigger the innate immune system by activating STING. It was found that the binding affinity of asymmetric 2'3'-cGAMP to symmetric dimer of STING is 3 orders of magnitude higher than that of the symmetric 3'3'-cyclic dinucleotides. Such a phenomenon has not been understood yet. Here we show that the subtle changes in phosphodiester linkage of CDNs lead to their distinct structural properties which correspond to the varied binding affinities. 2'-5' and/or 3'-5' linked CDNs adopt specific while different types of ribose puckers and backbone conformations. That ribose conformations and base types have different propensities for anti or syn glycosidic conformations further affects the overall flexibility of CDNs. The counterbalance between backbone ring tension and electrostatic repulsion, both affected by the ring size, also contributes to the different flexibility of CDNs. Our calculations reveal that the free energy cost for 2'3'-cGAMP to adopt the STING-bound structure is smaller than that for 3'3'-cGAMP and cyclic-di-GMP. These findings may serve as a reference for design of CDN-analogues as vaccine adjuvants. Moreover, the cyclization pattern of CDNs closely related to their physiological roles suggests the importance of understanding structural properties in the study of protein-ligand interactions.

Elastic Energy Partitioning in DNA Deformation and Binding to Proteins

We study the elasticity of DNA based on local principal axes of bending identified from over 0.9-μs all-atom molecular dynamics simulations of DNA oligos. The calculated order parameters describe motion of DNA as an elastic rod. In 10 possible dinucleotide steps, bending about the two principal axes is anisotropic yet linearly elastic. Twist about the centroid axis is largely decoupled from bending, but DNA tends to overtwist for unbending beyond the typical range of thermal motion, which is consistent with experimentally observed twist-stretch coupling. The calculated elastic stiffness of dinucleotide steps yield sequence-dependent persistence lengths consistent with previous single-molecule experiments, which is further analyzed by performing coarse-grained simulations of DNA. Flexibility maps of oligos constructed from simulation also match with those from the precalculated stiffness of dinucleotide steps. These support the premise that base pair interaction at the dinucleotide-level is mainly responsible for the elasticity of DNA. Furthermore, we analyze 1381 crystal structures of protein-DNA complexes. In most structures, DNAs are mildly deformed and twist takes the highest portion of the total elastic energy. By contrast, in structures with the elastic energy per dinucleotide step greater than about 4.16 kBT (kBT: thermal energy), the major bending becomes dominant. The extensional energy of dinucleotide steps takes at most 35% of the total elastic energy except for structures containing highly deformed DNAs where linear elasticity breaks down. Such partitioning between different deformational modes provides quantitative insights into the conformational dynamics of DNA as well as its interaction with other molecules and surfaces.

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.

Computational investigation of curcumin, a natural polyphenol that inhibits the destabilization and the aggregation of human SOD1 mutant (Ala4Val)

Curcumin inhibits the aberrant aggregation in mutant SOD1 protein, thereby decreasing the propensity of β-sheets and the toxicity level.

Enhanced conformational sampling technique provides an energy landscape view of large-scale protein conformational transitions

A novel in silico approach (NMA–ITS) is introduced to rapidly and effectively sample the configuration space and give quantitative data for exploring the conformational changes of proteins.

RNA approaches the B-form in stacked single strand dinucleotide contexts

Duplex RNA adopts an A-form structure, while duplex DNA interconverts between the A- and B-forms depending on the environment. The C2'-endo sugar pucker seen in B-form DNA can occur infrequently in ribose sugars as well, but RNA is not understood to assume B-form conformations. Through analysis of over 45,000 stacked single strand dinucleotide (SSD) crystal structure conformations, this study demonstrates that RNA is capable of adopting a wide conformational range between the canonical A- and B-forms at the localized SSD level, including many B-form-like conformations. It does so through C2'-endo ribose conformations in one or both nucleotides, and B-form-like neighboring base stacking patterns. As chemical reactions on nucleic acids involve localized changes in chemical bonds, the understanding of how enzymes distinguish between DNA and RNA nucleotides is altered by the energetic accessibility of these rare B-form-like RNA SSD conformations. The existence of these conformations also has direct implications in parametrization of molecular mechanics energy functions used extensively to model nucleic acid behavior., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 65-82, 2016.

Mechanism of the αβ conformational change in F1-ATPase after ATP hydrolysis: free-energy simulations

One of the motive forces for F1-ATPase rotation is the conformational change of the catalytically active β subunit due to closing and opening motions caused by ATP binding and hydrolysis, respectively. The closing motion is accomplished in two steps: the hydrogen-bond network around ATP changes and then the entire structure changes via B-helix sliding, as shown in our previous study. Here, we investigated the opening motion induced by ATP hydrolysis using all-atom free-energy simulations, combining the nudged elastic band method and umbrella sampling molecular-dynamics simulations. Because hydrolysis requires residues in the α subunit, the simulations were performed with the αβ dimer. The results indicate that the large-scale opening motion is also achieved by the B-helix sliding (in the reverse direction). However, the sliding mechanism is different from that of ATP binding because sliding is triggered by separation of the hydrolysis products ADP and Pi. We also addressed several important issues: 1), the timing of the product Pi release; 2), the unresolved half-closed β structure; and 3), the ADP release mechanism. These issues are fundamental for motor function; thus, the rotational mechanism of the entire F1-ATPase is also elucidated through this αβ study. During the conformational change, conserved residues among the ATPase proteins play important roles, suggesting that the obtained mechanism may be shared with other ATPase proteins. When combined with our previous studies, these results provide a comprehensive view of the β-subunit conformational change that drives the ATPase.

Increasing the sampling efficiency of protein conformational transition using velocity-scaling optimized hybrid explicit/implicit solvent REMD simulation

The application of temperature replica exchange molecular dynamics (REMD) simulation on protein motion is limited by its huge requirement of computational resource, particularly when explicit solvent model is implemented. In the previous study, we developed a velocity-scaling optimized hybrid explicit/implicit solvent REMD method with the hope to reduce the temperature (replica) number on the premise of maintaining high sampling efficiency. In this study, we utilized this method to characterize and energetically identify the conformational transition pathway of a protein model, the N-terminal domain of calmodulin. In comparison to the standard explicit solvent REMD simulation, the hybrid REMD is much less computationally expensive but, meanwhile, gives accurate evaluation of the structural and thermodynamic properties of the conformational transition which are in well agreement with the standard REMD simulation. Therefore, the hybrid REMD could highly increase the computational efficiency and thus expand the application of REMD simulation to larger-size protein systems.

Structural investigation on the electrostatic loop of native and mutated SOD1 and their interaction with therapeutic compounds

The electrostatic loop of the native and mutated SOD1 protein with single point mutation in the loop is subjected to MD simulation. The structure and electrostatic properties of the native and mutated loops before/after interacting with small compounds are compared.

Single-Walled Carbon Nanotubes Modulate the B- to A-DNA Transition

We study the conformational equilibrium between B-to-A forms of ds-DNA adsorbed onto a single-walled carbon nanotube (SWNT) using free energy profile calculations based on all-atom molecular dynamics simulations. The potential of mean force (PMF) of the B-to-A transition of ds-DNA in the presence of an uncharged (10,0) carbon nanotube for two dodecamers with poly-AT or poly-GC sequences is calculated as a function of a root-mean-square-distance (ΔRMSD) difference metric for the B-to-A transition. The calculations reveal that in the presence of a SWNT DNA favors B-form DNA significantly in both poly-GC and poly-AT sequences. Furthermore, the poly-AT DNA:SWNT complex shows a higher energy penalty for adopting an A-like conformation than poly-GC DNA:SWNT by several kcal/mol. The presence of a SWNT on either poly-AT or poly-GC DNA affects the PMF of the transition such that the B form is favored by as much as 10 kcal/mol. In agreement with published data, we find a potential energy minimum between A and B-form DNA at ΔRMSD ≈ -1.5 Å and that the presence of the SWNT moves this minimum by as much as ΔRMSD = 3 Å.

Mapping central α-helix linker mediated conformational transition pathway of calmodulin via simple computational approach

The effects of intrinsic structural flexibility of calmodulin protein on the mechanism of its allosteric conformational transition are investigated in this article. Using a novel in silico approach, the conformational transition pathways of intact calmodulin as well as the isolated N- and C- terminal domains are identified and energetically characterized. It is observed that the central α-helix linker amplifies the structural flexibility of intact Ca(2+)-free calmodulin, which might facilitate the transition of the two domains. As a result, the global conformational transition of Ca(2+)-free calmodulin is initiated by the barrierless transition of two domains and proceeds through the barrier associated unwinding and bending of the central α-helix linker. The binding of Ca(2+) cations to calmodulin further increases the structural flexibility of the C-terminal domain and results in a downhill transition pathway of which all regions transit in a concerted manner. On the other hand, the separation of the N- and C-terminal domains from calmodulin protein loses the mediating function of central α-helix linker, leading to more difficult conformational transitions of both domains. The present study provides novel insights into the correlation of the integrity of protein, the structural flexibility, and the mechanism of conformational transition of proteinlike calmodulin.

All-atom polarizable force field for DNA based on the classical Drude oscillator model

Presented is a first generation atomistic force field (FF) for DNA in which electronic polarization is modeled based on the classical Drude oscillator formalism. The DNA model is based on parameters for small molecules representative of nucleic acids, including alkanes, ethers, dimethylphosphate, and the nucleic acid bases and empirical adjustment of key dihedral parameters associated with the phosphodiester backbone, glycosidic linkages, and sugar moiety of DNA. Our optimization strategy is based on achieving a compromise between satisfying the properties of the underlying model compounds in the gas phase targeting quantum mechanical (QM) data and reproducing a number of experimental properties of DNA duplexes in the condensed phase. The resulting Drude FF yields stable DNA duplexes on the 100-ns time scale and satisfactorily reproduce (1) the equilibrium between A and B forms of DNA and (2) transitions between the BI and BII substates of B form DNA. Consistency with the gas phase QM data for the model compounds is significantly better for the Drude model as compared to the CHARMM36 additive FF, which is suggested to be due to the improved response of the model to changes in the environment associated with the explicit inclusion of polarizability. Analysis of dipole moments associated with the nucleic acid bases shows the Drude model to have significantly larger values than those present in CHARMM36, with the dipoles of individual bases undergoing significant variations during the MD simulations. Additionally, the dipole moment of water was observed to be perturbed in the grooves of DNA.

Cytosine unstacking and strand slippage at an insertion-deletion mutation sequence in an overhang-containing DNA duplex

Base unstacking in template strands, when accompanied by strand slippage, can result in deletion mutations during strand extension by nucleic acid polymerases. In a GCCC mutation hot-spot sequence, which was previously identified to have a 50% probability of causing such mutations during DNA replication by a Y-family polymerase, a single-base deletion mutation could result from such unstacking of any one of its three template cytosines. In this study, the intrinsic energetic differences in unstacking among these three cytosines in a solvated DNA duplex overhang model were examined using umbrella sampling molecular dynamics simulations. The free energy profiles obtained show that cytosine unstacking grows progressively more unfavorable as one moves inside the duplex from the 5′-end of the overhang template strand. Spontaneous strand slippage occurs in response to such base unstacking in the direction of both the major and minor grooves for all three cytosines. Unrestrained simulations run from three distinct strand-slipped states and one non-strand-slipped state suggest that a more duplexlike environment can help stabilize strand slippage. The possible underlying reasons and biological implications of these observations are discussed in the context of nucleic acid replication active site dynamics.

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.

Energetic and structural details of the trigger-loop closing transition in RNA polymerase II

An evolutionarily conserved element in RNA polymerase II, the trigger loop (TL), has been suggested to play an important role in the elongation rate, fidelity of selection of the matched nucleoside triphosphate (NTP), catalysis of transcription elongation, and translocation in both eukaryotes and prokaryotes. In response to NTP binding, the TL undergoes large conformational changes to switch between distinct open and closed states to tighten the active site and avail catalysis. A computational strategy for characterizing the conformational transition pathway is presented to bridge the open and closed states of the TL. Information from a large number of independent all-atom molecular dynamics trajectories from Hamiltonian replica exchange and targeted molecular dynamics simulations is gathered together to assemble a connectivity map of the conformational transition. The results show that with a cognate NTP, TL closing should be a spontaneous process. One major intermediate state is identified along the conformational transition pathway, and the key structural features are characterized. The complete pathway from the open TL to the closed TL provides a clear picture of the TL closing.

The free energy landscape of pseudorotation in 3'-5' and 2'-5' linked nucleic acids

The five-membered furanose ring is a central component of the chemical structure of biological nucleic acids. The conformations of the furanose ring can be analytically described using the concept of pseudorotation, and for RNA and DNA they are dominated by the C2'-endo and C3'-endo conformers. While the free energy difference between these two conformers can be inferred from NMR measurements, a free energy landscape of the complete pseudorotation cycle of nucleic acids in solution has remained elusive. Here, we describe a new free energy calculation method for molecular dynamics (MD) simulations using the two pseudorotation parameters directly as the collective variables. To validate our approach, we calculated the free energy surface of ribose pseudorotation in guanosine and 2'-deoxyguanosine. The calculated free energy landscape reveals not only the relative stability of the different pseudorotation conformers, but also the main transition path between the stable conformations. Applying this method to a standard A-form RNA duplex uncovered the expected minimum at the C3'-endo state. However, at a 2'-5' linkage, the minimum shifts to the C2'-endo conformation. The free energy of the C3'-endo conformation is 3 kcal/mol higher due to a weaker hydrogen bond and a reduced base stacking interaction. Unrestrained MD simulations suggest that the conversion from C3'-endo to C2'-endo and vice versa is on the nanosecond and microsecond time scale, respectively. These calculations suggest that 2'-5' linkages may enable folded RNAs to sample a wider spectrum of their pseudorotation conformations.