Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study.

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 µs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the µs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the µs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-µs time scale, which could easily lead to inappropriate interpretation of the simulation data.

Understanding RNA Flexibility Using Explicit Solvent Simulations: The Ribosomal and Group I Intron Reverse Kink-Turn Motifs.

Reverse kink-turn is a recurrent elbow-like RNA building block occurring in the ribosome and in the group I intron. Its sequence signature almost matches that of the conventional kink-turn. However, the reverse and conventional kink-turns have opposite directions of bending. The reverse kink-turn lacks basically any tertiary interaction between its stems. We report unrestrained, explicit solvent molecular dynamics simulations of ribosomal and intron reverse kink-turns (54 simulations with 7.4 µs of data in total) with different variants (ff94, ff99, ff99bsc0, ff99χOL, and ff99bsc0χOL) of the Cornell et al. force field. We test several ion conditions and two water models. The simulations characterize the directional intrinsic flexibility of reverse kink-turns pertinent to their folded functional geometries. The reverse kink-turns are the most flexible RNA motifs studied so far by explicit solvent simulations which are capable at the present simulation time scale to spontaneously and reversibly sample a wide range of geometries from tightly kinked ones through flexible intermediates up to extended, unkinked structures. A possible biochemical role of the flexibility is discussed. Among the tested force fields, the latest χOL variant is essential to obtaining stable trajectories while all force field versions lacking the χ correction are prone to a swift degradation toward senseless ladder-like structures of stems, characterized by high-anti glycosidic torsions. The type of explicit water model affects the simulations considerably more than concentration and the type of ions.

Explicit Water Models Affect the Specific Solvation and Dynamics of Unfolded Peptides While the Conformational Behavior and Flexibility of Folded Peptides Remain Intact.

Conventional molecular dynamics simulations on 50 ns to 1 µs time scales were used to study the effects of explicit solvent models on the conformational behavior and solvation of two oligopeptide solutes: α-helical EK-peptide (14 amino acids) and a ß-hairpin chignolin (10 amino acids). The widely used AMBER force fields (ff99, ff99SB, and ff03) were combined with four of the most commonly used explicit solvent models (TIP3P, TIP4P, TIP5P, and SPC/E). Significant differences in the specific solvation of chignolin among the studied water models were identified. Chignolin was highly solvated in TIP5P, whereas reduced specific solvation was found in the TIP4P, SPC/E, and TIP3P models for kinetic, thermodynamic, and both kinetic and thermodynamic reasons, respectively. The differences in specific solvation did not influence the dynamics of structured parts of the folded peptide. However, substantial differences between TIP5P and the other models were observed in the dynamics of unfolded chignolin, stability of salt bridges, and specific solvation of the backbone carbonyls of EK-peptide. Thus, we conclude that the choice of water model may affect the dynamics of flexible parts of proteins that are solvent-exposed. On the other hand, all water models should perform similarly for well-structured folded protein regions. The merits of the TIP3P model include its high and overestimated mobility, which accelerates simulation processes and thus effectively increases sampling.

Changes in electrostatic surface potential of Na+/K+-ATPase cytoplasmic headpiece induced by cytoplasmic ligand(s) binding.

A set of single-tryptophan mutants of the Na(+)/K(+)-ATPase isolated, large cytoplasmic loop connecting transmembrane helices M4 and M5 (C45) was prepared to monitor effects of the natural cytoplasmic ligands (i.e., Mg(2+) and/or ATP) binding. We introduced a novel method for the monitoring of the changes in the electrostatic surface potential (ESP) induced by ligand binding, using the quenching of the intrinsic tryptophan fluorescence by acrylamide or iodide. This approach opens a new way to understanding the interactions within the proteins. Our experiments revealed that the C45 conformation in the presence of the ATP (without magnesium) substantially differed from the conformation in the presence of Mg(2+) or MgATP or in the absence of any ligand not only in the sense of geometry but also in the sense of the ESP. Notably, the set of ESP-sensitive residues was different from the set of geometry-sensitive residues. Moreover, our data indicate that the effect of the ligand binding is not restricted only to the close environment of the binding site and that the information is in fact transmitted also to the distal parts of the molecule. This property could be important for the communication between the cytoplasmic headpiece and the cation binding sites located within the transmembrane domain.

ATP and magnesium drive conformational changes of the Na+/K+-ATPase cytoplasmic headpiece.

Conformational changes of the Na(+)/K(+)-ATPase isolated large cytoplasmic segment connecting transmembrane helices M4 and M5 (C45) induced by the interaction with enzyme ligands (i.e. Mg(2+) and/or ATP) were investigated by means of the intrinsic tryptophan fluorescence measurement and molecular dynamic simulations. Our data revealed that this model system consisting of only two domains retained the ability to adopt open or closed conformation, i.e. behavior, which is expected from the crystal structures of relative Ca(2+)-ATPase from sarco(endo)plasmic reticulum for the corresponding part of the entire enzyme. Our data revealed that the C45 is found in the closed conformation in the absence of any ligand, in the presence of Mg(2+) only, or in the simultaneous presence of Mg(2+) and ATP. Binding of the ATP alone (i.e. in the absence of Mg(2+)) induced open conformation of the C45. The fact that the transmembrane part of the enzyme was absent in our experiments suggested that the observed conformational changes are consequences only of the interaction with ATP or Mg(2+) and may not be related to the transported cations binding/release, as generally believed. Our data are consistent with the model, where ATP binding to the low-affinity site induces conformational change of the cytoplasmic part of the enzyme, traditionally attributed to E2-->E1 transition, and subsequent Mg(2+) binding to the enzyme-ATP complex induces in turn conformational change traditionally attributed to E1-->E2 transition.

Two C-terminal ankyrin repeats form the minimal stable unit of the ankyrin repeat protein p18INK4c.

Ankyrin repeat proteins (ARPs) appear to be abundant in organisms from all phyla, and play critical regulatory roles, mediating specific interactions with target biomolecules and thus ordering the sequence of events in diverse cellular processes. ARPs possess a non-globular scaffold consisting of repeating motifs named ankyrin (ANK) repeats, which stack on each other. The modular architecture of ARPs provides a new paradigm for understanding protein stability and folding mechanisms. In the present study, the stability of various C-terminal fragments of the ARP p18(INK4c) was investigated by all-atomic 450 ns molecular dynamics (MD) simulations in explicit water solvent. Only motifs with at least two ANK repeats made stable systems in the available timescale. All smaller fragments were unstable, readily losing their native fold and alpha-helical content. Since each non-terminal ANK repeat has two hydrophobic sides, we may hypothesize that at least one hydrophobic side must be fully covered and shielded from the water as a necessary, but not sufficient, condition to maintain ANK repeat stability. Consequently, at least two ANK repeats are required to make a stable ARP.

How the stabilization of INK4 tumor suppressor 3D structure evaluated by quantum chemical and molecular mechanics calculations corresponds well with experimental results: interplay of association enthalpy, entropy, and solvation effects.

The folding free energy of the INK4c tumor suppressor core, consisting of 10 helices, was determined as the sum of gas-phase interaction enthalpy, gas-phase interaction entropy, and dehydration and hydration free energy. The interaction energy and the hydration free energy were determined using the nonempirical density functional theory (DFT) method, augmented by a dispersion-energy correction term, the semiempirical density-functional tight-binding method covering the dispersion energy, and the density functional theory/conductor-like screening model (DFT/COSMO) procedure, whereas the interaction entropy was calculated with the empirical Cornell et al. force field. Alternatively, all contributions were evaluated consistently using empirical methods. All the values of the interaction energy of helix pairs are stabilizing, and the dominant stabilizing terms stem from the London dispersion energy and, in the case of charged systems, the electrostatic energy. The stabilization energy of the core, determined as the difference of the energy of the core and 10 separate helices, amounts to approximately 450 kcal/mol. Systematically, the difference in the hydration free energy of a helix pair and its separate components is smaller in magnitude than the interaction energy, and it is negative for some pairs while positive for others. The average total free energy of a core formation amounts to -29.6 kcal/mol (yielded by scaled quantum-chemical methods) and +13.9 kcal/mol (resulting from empirical methods). These values are considerably smaller than their single components, which are dominated by the interaction energy. The computationally predicted interval encloses the experimental value of the folding free energy (-2.8 kcal/mol).