Analysis of molecular dynamics simulations of 10-residue peptide, chignolin, using statistical mechanics: Relaxation mode analysis and three-dimensional reference interaction site model theory.
Ontology highlight
ABSTRACT:
SUBMITTER: Maruyama Y
PROVIDER: S-EPMC6975981 | biostudies-literature | 2019
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Statistical mechanics provides the link between microscopic properties of many-particle systems and macroscopic properties such as pressure and temperature. Observations of similar "microscopic" quantities exist for the motion of zooplankton, as well as many species of other social animals. Herein, we propose to take average squared velocities as the definition of the "ecological temperature" of a population under different conditions on nutrients, light, oxygen and others. We test the usefulness of this definition on observations of the crustacean zooplankton Daphnia pulicaria. In one set of experiments, D. pulicaria is infested with the pathogen Vibrio cholerae, the causative agent of cholera. We find that infested D. pulicaria under light exposure have a significantly greater ecological temperature, which puts them at a greater risk of detection by visual predators. In a second set of experiments, we observe D. pulicaria in cold and warm water, and in darkness and under light exposure. Overall, our ecological temperature is a good discriminator of the crustacean's swimming behavior.
| S-EPMC4535848 | biostudies-literature
Project description:Understanding the dominant factor in thermodynamic stability of proteins remains an open challenge. Kauzmann's hydrophobic interaction hypothesis, which considers hydrophobic interactions between nonpolar groups as the dominant factor, has been widely accepted for about sixty years and attracted many scientists. The hypothesis, however, has not been verified or disproved because it is difficult, both theoretically and experimentally, to quantify the solvent effects on the free energy change in protein folding. Here, we developed a computational method for extracting the dominant factor behind thermodynamic stability of proteins and applied it to a small, designed protein, chignolin. The resulting free energy profile quantitatively agreed with the molecular dynamics simulations. Decomposition of the free energy profile indicated that intramolecular interactions predominantly stabilized collapsed conformations, whereas solvent-induced interactions, including hydrophobic ones, destabilized them. These results obtained for chignolin were consistent with the site-directed mutagenesis and calorimetry experiments for globular proteins with hydrophobic interior cores.
| S-EPMC6435801 | biostudies-literature
Project description:We use molecular dynamics simulations of the Kremer-Grest (KG) bead-spring model of polymer chains of length between 10 and 500, and a closely related analogue that allows for chain crossing, to clearly delineate the effects of entanglements on the length-scale-dependent chain relaxation in polymer melts. We analyze the resulting trajectories using the Rouse modes of the chains and find that entanglements strongly affect these modes. The relaxation rates of the chains show two limiting effective monomeric frictions, with the local modes experiencing much lower effective friction than the longer modes. The monomeric relaxation rates of longer modes vary approximately inversely with chain length due to kinetic confinement effects. The time-dependent relaxation of Rouse modes has a stretched exponential character with a minimum of stretching exponent in the vicinity of the entanglement chain length. None of these trends are found in models that allow for chain crossing. These facts, in combination, argue for the confined motion of chains for time scales between the entanglement time and their ultimate free diffusion.
| S-EPMC4196748 | biostudies-other
Project description:The concept of allostery in which macromolecules switch between two different conformations is a central theme in biological processes ranging from gene regulation to cell signaling to enzymology. Allosteric enzymes pervade metabolic processes, yet a simple and unified treatment of the effects of allostery in enzymes has been lacking. In this work, we take a step toward this goal by modeling allosteric enzymes and their interaction with two key molecular players-allosteric regulators and competitive inhibitors. We then apply this model to characterize existing data on enzyme activity, comment on how enzyme parameters (such as substrate binding affinity) can be experimentally tuned, and make novel predictions on how to control phenomena such as substrate inhibition.
| S-EPMC5452729 | biostudies-literature
Project description:Music, while allowing nearly unlimited creative expression, almost always conforms to a set of rigid rules at a fundamental level. The description and study of these rules, and the ordered structures that arise from them, is the basis of the field of music theory. Here, I present a theoretical formalism that aims to explain why basic ordered patterns emerge in music, using the same statistical mechanics framework that describes emergent order across phase transitions in physical systems. I first apply the mean field approximation to demonstrate that phase transitions occur in this model from disordered sound to discrete sets of pitches, including the 12-fold octave division used in Western music. Beyond the mean field model, I use numerical simulation to uncover emergent structures of musical harmony. These results provide a new lens through which to view the fundamental structures of music and to discover new musical ideas to explore.
| S-EPMC6524979 | biostudies-literature
Project description:BackgroundAnalysis of known protein structures reveals that identical sequence fragments in proteins can adopt different secondary structure conformations. The extent of this conformational diversity is influenced by various factors like the intrinsic sequence propensity, sequence context and other environmental factors such as pH, site directed mutations or alteration of the binding ligands. Understanding the mechanism by which the environment affects the structural ambivalence of these peptides has potential implications for protein design and reliable local structure prediction algorithms. Identification of the structurally ambivalent sequence fragments and determining the rules which dictate their conformational preferences play an important role in understanding the conformational changes observed in misfolding diseases. However, a systematic classification of their intrinsic sequence patterns or a statistical analysis of their properties and sequence context in relation to the origin of their structural diversity have largely remained unexplored.ResultsIn this work, the conformational variability of α-helices is studied by mapping sequences from the non-redundant database to identical sequences across all classes of the SCOP (Structural Classification of Proteins) database. Some helices retain their conformations when mapped in the SCOP database while others exhibit a complete/partial switch to non-helical conformations. The results clearly depict the differences in the propensities of amino acids for the variable and conserved helices. Sequences flanking these ambivalent sequence fragments have anisotropic propensities at the N- and C-termini. This structural variability is depicted by molecular dynamics simulations in explicit solvent, which show that the short conserved helices retain their conformations while their longer counterparts fray into two or more shorter helices. Variable helices in the non-redundant database exhibit a trend of retaining helical conformations while their corresponding non-helical conformations in SCOP database show large deviations from their respective initial structures by adopting partial or full helical conformations. Partially ambivalent helices are also found to retain their respective conformations.ConclusionsAll sequence fragments which show structural diversity in different proteins of the non-redundant database are investigated. The final conformation of these ambivalent sequences are dictated by a fine tuning of their intrinsic sequence propensity and the anisotropic amino acid propensity of the flanking sequences. This analysis may unravel the connection between diverse secondary structures, which conserve the overall structural fold of the protein thus determining its function.
| S-EPMC2973962 | biostudies-literature
Project description:Molecular dynamics (MD) simulations and nuclear magnetic resonance spin-relaxation measurements provide detailed insights into ps-ns structural dynamics of proteins. An analysis of discrepancies between the two methods is presented for the B3 immunoglobulin-binding domain of streptococcal protein G. MD simulations using three MD force fields (OPLS-AA, AMBER ff99SB, and AMBER ff03) overestimate the flexibility of backbone N--H vectors at the borders of secondary structure and in loops when compared with experimentally determined backbone amide generalized order parameters (Hall and Fushman, J Am Chem Soc 2006; 12:7855-7870). Comparison with a previous study of residual dipolar coupling constants (Bouvignies et al., Proc Natl Acad Sci USA 2005;102:13885-13890) indicates that slower timescale motions do not account for the discrepancies. Structural analysis reveals that relative imbalance between the description of hydrogen bonding and other terms of modern force fields may be responsible for disagreement.
| S-EPMC2750878 | biostudies-literature
Project description:Since Jacob and Monod's characterization of the role of DNA elements in gene control, it has been recognized that the linear organization of genome structure is important for the regulation of gene transcription and hence the manifestation of phenotypes. Similarly, it has long been hypothesized that the spatial organization (in three dimensions evolving through time), as part of the epigenome, makes a significant contribution to the genotype-phenotype transition. Proximity ligation assays commonly known as chromosome conformation capture (3C) and 3C based methodologies (e.g., GCC, HiC and ChIA-Pet) are increasingly being incorporated into empirical studies to investigate the role that three-dimensional genome structure plays in the regulation of phenotype. The apparent simplicity of these methodologies-crosslink chromatin, digest, dilute, ligate, detect interactions-belies the complexity of the data and the considerations that should be taken into account to ensure the generation and accurate interpretation of reliable data. Here we discuss the probabilistic nature of these methodologies and how this contributes to their endogenous limitations.
| S-EPMC3899129 | biostudies-literature
Project description:Unicellular living organisms, such as bacteria and algae, propel themselves through a medium via cyclic strokes involving the motion of cilia and flagella. Dense populations of such "active particles" or "swimmers" exhibit a rich collective behavior at large scales. Starting with a minimal physical model of a stroke-averaged swimmer in a fluid, we derive a continuum description of a suspension of active organisms that incorporates fluid-mediated, long-range hydrodynamic interactions among the swimmers. Our work demonstrates that hydrodynamic interactions provide a simple, generic origin for several nonequilibrium phenomena predicted or observed in the literature. The continuum model derived here does not depend on the microscopic physical model of the individual swimmer. The details of the large-scale physics do, however, differ for "shakers" (particles that are active but not self-propelled, such as melanocytes) and "movers" (self-propelled particles), "pushers" (most bacteria) and "pullers" (algae like Chlamydomonas). Our work provides a classification of the large-scale behavior of all these systems.
| S-EPMC2747162 | biostudies-literature
Project description:Flocking is a typical example of emergent collective behavior, where interactions between individuals produce collective patterns on the large scale. Here we show how a quantitative microscopic theory for directional ordering in a flock can be derived directly from field data. We construct the minimally structured (maximum entropy) model consistent with experimental correlations in large flocks of starlings. The maximum entropy model shows that local, pairwise interactions between birds are sufficient to correctly predict the propagation of order throughout entire flocks of starlings, with no free parameters. We also find that the number of interacting neighbors is independent of flock density, confirming that interactions are ruled by topological rather than metric distance. Finally, by comparing flocks of different sizes, the model correctly accounts for the observed scale invariance of long-range correlations among the fluctuations in flight direction.
| S-EPMC3324025 | biostudies-literature