Change in backbone torsion angle distribution on protein folding

Robosample: A rigid-body molecular simulation program based on robot mechanics

BACKGROUND

Compared with all-atom molecular dynamics (MD), constrained MD methods allow for larger time steps, potentially reducing computational cost. For this reason, there has been continued interest in improving constrained MD algorithms to increase configuration space sampling in molecular simulations.

METHODS

Here, we introduce Robosample, a software package that implements high-performance constrained dynamics algorithms, originally developed for robotics, and applies them to simulations of biomolecular systems. As in the gMolmodel package developed by Spiridon and Minh in 2017, Robosample uses Constrained Dynamics Hamiltonian Monte Carlo (CDHMC) as a Gibbs sampling move - a type of Monte Carlo move where a subset of coordinates is allowed to change. In addition to the previously described Cartesian and torsional dynamics moves, Robosample implements spherical and cylindrical joints that can be distributed along the molecule by the user.

RESULTS

In alanine dipeptide simulations, the free energy surface is recovered by mixing fully flexible with torsional, cylindrical, or spherical dynamics moves. Ramachandran dynamics, where only the two key torsions are mobile, accelerate the slowest transition by an order of magnitude. We also show that simulations of a complex glycan cover significantly larger regions of the configuration space when mixed with constrained dynamics.

MAJOR CONCLUSIONS

Robosample is a tool of choice for efficient conformational sampling of large biomolecules.

GENERAL SIGNIFICANCE

Robosample is intended as a reliable and user-friendly simulation package for fast biomolecular sampling that does not require extensive expertise in mechanical engineering or in the statistical mechanics of reduced coordinates.

Protein phenotype diagnosis of autosomal dominant calmodulin mutations causing irregular heart rhythms

The life-threatening group of irregular cardiac rhythmic disorders also known as Cardiac Arrhythmias (CA) are caused by mutations in highly conserved Calmodulin (CALM/CaM) genes. Herein, we present a multidimensional approach to diagnose changes in phenotypic, stability, and Ca²⁺ ion binding properties of CA-causing mutations. Mutation pathogenicity was determined by diverse computational machine learning approaches. We further modeled the mutations in 3D protein structure and analyzed residue level phenotype plasticity. We have also examined the influence of torsion angles, number of H-bonds, and free energy dynamics on the stability, near-native simulation dynamic potential of residue fluctuations in protein structures, Ca²⁺ ion binding potentials, of CaM mutants. Our study recomends to use M-CAP method for measuring the pathogenicity of CA causing CaM variants. Interestingly, most CA-causing variants we analyzed, exists in either third (V/H-96, S/I-98, V-103) or fourth (G/V-130, V/E/H-132, H-134, P-136, G-141, and L-142) EF-hands located in carboxyl domains of the CaM molecule. We observed that the minor structural fluctuations caused by these variants are likely tolerable owing to the highly flexible nature of calmodulin's globular domains. However, our molecular docking results supports that these variants disturb the affinity of CaM toward Ca²⁺ ions and corroborate previous findings from functional studies. Taken together, these computational findings can explain the molecular reasons for subtle changes in structure, flexibility, and stability aspects of mutant CaM molecule. Our comprehensive molecular scanning approach demonstrates the utility of computational methods in quick preliminary screening of CA- CaM mutations before undertaking time consuming and complicated functional laboratory assays.

Using C' deviation to study structures of central amino acids in peptide fragments

In this investigation, we attempted to study the backbone geometry of amino acids in peptides using C' deviation. Diameters of distribution were used to describe the various atomic structures, and scatter graphs provided visual evaluation. The length of peptide fragments and the secondary structure of amino acids in the central position of the peptide fragments were also analyzed. The results showed that the atomic distribution of the central amino acids of five-residue peptide fragments was much more restricted than that of their corresponding three-residue peptide fragments. In identical three-residue fragments, atoms of central amino acids with different secondary structures, were distributed in distinct areas.

A simple model for polyproline II structure in unfolded states of alanine-based peptides

The striking similarity between observed circular dichroism spectra of nonprolyl homopolymers and that of regular left-handed polyproline II (P(II)) helices prompted Tiffany and Krimm to propose in 1968 that unordered peptides and unfolded proteins are built of P(II) segments linked by sharp bends. A large body of experimental evidence, accumulated over the past three decades, provides compelling evidence in support of the original hypothesis of Tiffany and Krimm. Of particular interest are the recent experiments of Shi et al. who find significant P(II) structure in a short unfolded alanine-based peptide. What is the physical basis for P(II) helices in peptide and protein unfolded states? The widely accepted view is that favorable chain-solvent hydrogen bonds lead to a preference for dynamical fluctuations about noncooperative P(II) helices in water. Is this preference simply a consequence of hydrogen bonding or is it a manifestation of a more general trend for unfolded states which are appropriately viewed as chains in a good solvent? The prevalence of closely packed interiors in folded proteins suggests that under conditions that favor folding, water-which is a better solvent for itself than for any polypeptide chain-expels the chain from its midst, thereby maximizing chain packing. Implicit in this view is a complementary idea: under conditions that favor unfolding, chain-solvent interactions are preferred and in a so-called good solvent, chain packing density is minimized. In this work we show that minimization of chain packing density leads to preferred fluctuations for short polyalanyl chains around canonical, noncooperative P(II)-like conformations. Minimization of chain packing is modeled using a purely repulsive soft-core potential between polypeptide atoms. Details of chain-solvent interactions are ignored. Remarkably, the simple model captures the essential physics behind the preference of short unfolded alanine-based peptides for P(II) helices. Our results are based on a detailed analysis of the potential energy landscape which determines the system's structural and thermodynamic preferences. We use the inherent structure formalism of Stillinger and Weber, according to which the energy landscape is partitioned into basins of attraction around local minima. We find that the landscape for the experimentally studied seven-residue alanine-based peptide is dominated by fluctuations about two noncooperative structures: the left-handed polyproline II helix and its symmetry mate.

Protein unfolding transitions in an intrinsically unstable annexin domain: molecular dynamics simulation and comparison with nuclear magnetic resonance data

Unfolding transitions of an intrinsically unstable annexin domain and the unfolded state structure have been examined using multiple approximately 10-ns molecular dynamics simulations. Three main basins are observed in the configurational space: native-like state, compact partially unfolded or intermediate compact state, and the unfolded state. In the native-like state fluctuations are observed that are nonproductive for unfolding. During these fluctuations, after an initial loss of approximately 20% of the core residue native contacts, the core of the protein transiently completely refolds to the native state. The transition from the native-like basin to the partially unfolded compact state involves approximately 75% loss of native contacts but little change in the radius of gyration or core hydration properties. The intermediate state adopts for part of the time in one of the trajectories a novel highly compact salt-bridge stabilized structure that can be identified as a conformational trap. The intermediate-to-unfolded state transition is characterized by a large increase in the radius of gyration. After an initial relaxation the unfolded state recovers a native-like topology of the domain. The simulated unfolded state ensemble reproduces in detail experimental nuclear magnetic resonance data and leads to a convincing complete picture of the unfolded domain.