Browse
Submit Data
Databases
API
Help

Dataset Information

11 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Molecular Dynamics Simulations of the Temperature Induced Unfolding of Crambin Follow the Arrhenius Equation.

ABSTRACT: Molecular dynamics simulations have been used extensively to model the folding and unfolding of proteins. The rates of folding and unfolding should follow the Arrhenius equation over a limited range of temperatures. This study shows that molecular dynamic simulations of the unfolding of crambin between 500K and 560K do follow the Arrhenius equation. They also show that while there is a large amount of variation between the simulations the average values for the rate show a very high degree of correlation.

SUBMITTER: Dalby A

PROVIDER: S-EPMC4629273 | biostudies-literature | 2015

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Molecular Dynamics Simulations of the Temperature Induced Unfolding of Crambin Follow the Arrhenius Equation.

Dalby Andrew A Shamsir Mohd Shahir MS

F1000Research 20150820

Molecular dynamics simulations have been used extensively to model the folding and unfolding of proteins. The rates of folding and unfolding should follow the Arrhenius equation over a limited range of temperatures. This study shows that molecular dynamic simulations of the unfolding of crambin between 500K and 560K do follow the Arrhenius equation. They also show that while there is a large amount of variation between the simulations the average values for the rate show a very high degree of co ...[more]

PMID: 26539292

Similar Datasets

High temperature unfolding of Bacillus anthracis amidase-03 by molecular dynamics simulations.

Project description:The stability of amidase-03 structure (a cell wall hydrolase protein) from Bacillus anthracis was studied using classical molecular dynamics (MD) simulation. This protein (GenBank accession number: NP_844822) contains an amidase-03 domain which is known to exhibit the catalytic activity of N-acetylmuramoyl-L-alanine amidase (digesting MurNAc-Lalanine linkage of bacterial cell wall). The amidase-03 enzyme has stability at high temperature due to the core formed by the combination of several secondary structure elements made of beta-sheets. We used root-mean-square-displacement (RMSD) of the simulated structure from its initial state to demonstrate the unfolding of the enzyme using its secondary structural elements. Results show that amidase-03 unfolds in transition state ensemble (TSE). The data suggests that alpha-helices unfold before beta-sheets from the core during simulation.

| S-EPMC2737499 | biostudies-literature

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding.

Project description:We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.

| S-EPMC1301720 | biostudies-other

Temperature-induced unfolding of epidermal growth factor (EGF): insight from molecular dynamics simulation.

Project description:Thermal disruption of protein structure and function is a potentially powerful therapeutic vehicle. With the emerging nanoparticle-targeting and femtosecond laser technology, it is possible to deliver heating locally to specific molecules. It is therefore important to understand how fast a protein can unfold or lose its function at high temperatures, such as near the water boiling point. In this study, the thermal damage of EGF was investigated by combining the replica exchange (136 replicas) and conventional molecular dynamics simulations. The REMD simulation was employed to rigorously explore the free-energy landscape of EGF unfolding. Interestingly, besides the native and unfolded states, we also observed a distinct molten globule (MG) state that retained substantial amount of native contacts. Based on the understanding that which the unfolding of EGF is a three-state process, we have examined the unfolding kinetics of EGF (N-->MG and MG-->D) with multiple 20-ns conventional MD simulations. The Arrhenius prefactors and activation energy barriers determined from the simulation are within the range of previously studied proteins. In contrast to the thermal damage of cells and tissues which take place on the time scale of seconds to hours at relatively low temperatures, the denaturation of proteins occur in nanoseconds when the temperature of heat bath approaches the boiling point.

| S-EPMC2923262 | biostudies-literature

Molecular dynamics simulations using temperature-enhanced essential dynamics replica exchange.

Project description:Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.

| S-EPMC1877756 | biostudies-literature

Evaluating the Arrhenius equation for developmental processes.

Project description:The famous Arrhenius equation is well suited to describing the temperature dependence of chemical reactions but has also been used for complicated biological processes. Here, we evaluate how well the simple Arrhenius equation predicts complex multi-step biological processes, using frog and fruit fly embryogenesis as two canonical models. We find that the Arrhenius equation provides a good approximation for the temperature dependence of embryogenesis, even though individual developmental intervals scale differently with temperature. At low and high temperatures, however, we observed significant departures from idealized Arrhenius Law behavior. When we model multi-step reactions of idealized chemical networks, we are unable to generate comparable deviations from linearity. In contrast, we find the two enzymes GAPDH and β-galactosidase show non-linearity in the Arrhenius plot similar to our observations of embryonic development. Thus, we find that complex embryonic development can be well approximated by the simple Arrhenius equation regardless of non-uniform developmental scaling and propose that the observed departure from this law likely results more from non-idealized individual steps rather than from the complexity of the system.

| S-EPMC8377445 | biostudies-literature

High temperature unfolding simulations of the TRPZ1 peptide.

Project description:We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 mus with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.

| S-EPMC2480669 | biostudies-literature

Unfolding of the amyloid ?-peptide central helix: mechanistic insights from molecular dynamics simulations.

Project description:Alzheimer's disease (AD) pathogenesis is associated with formation of amyloid fibrils caused by polymerization of the amyloid ?-peptide (A?), which is a process that requires unfolding of the native helical structure of A?. According to recent experimental studies, stabilization of the A? central helix is effective in preventing A? polymerization into toxic assemblies. To uncover the fundamental mechanism of unfolding of the A? central helix, we performed molecular dynamics simulations for wild-type (WT), V18A/F19A/F20A mutant (MA), and V18L/F19L/F20L mutant (ML) models of the A? central helix. It was quantitatively demonstrated that the stability of the ?-helical conformation of both MA and ML is higher than that of WT, indicating that the ?-helical propensity of the three nonpolar residues (18, 19, and 20) is the main factor for the stability of the whole A? central helix and that their hydrophobicity plays a secondary role. WT was found to completely unfold by a three-step mechanism: 1) loss of ?-helical backbone hydrogen bonds, 2) strong interactions between nonpolar sidechains, and 3) strong interactions between polar sidechains. WT did not completely unfold in cases when any of the three steps was omitted. MA and ML did not completely unfold mainly due to the lack of the first step. This suggests that disturbances in any of the three steps would be effective in inhibiting the unfolding of the A? central helix. Our findings would pave the way for design of new drugs to prevent or retard AD.

| S-EPMC3049775 | biostudies-literature

Decoupling the Arrhenius equation <i>via</i> mechanochemistry.

Project description:Mechanochemistry continues to reveal new possibilities in chemistry including the opportunity for "greening" reactions. Nevertheless, a clear understanding of the energetic transformations within mechanochemical systems remains elusive. We employed a uniquely modified ball mill and strategically chosen Diels-Alder reactions to evaluate the role of several ball-milling variables. This revealed three different energetic regions that we believe are defining characteristics of most, if not all, mechanochemical reactors. Relative to the locations of a given ball mill's regions, activation energy determines whether a reaction is energetically easy (Region I), challenging (Region II), or unreasonable (Region III) in a given timeframe. It is in Region II, that great sensitivity to mechanochemical conditions such as vial material and oscillation frequency emerge. Our unique modifications granted control of reaction vessel temperature, which in turn allowed control of the locations of Regions I, II, and III for our mill. Taken together, these results suggest envisioning vibratory mills (and likely other mechanochemical methodologies) as molecular-collision facilitating devices that act upon molecules occupying a thermally-derived energy distribution. This unifies ball-milling energetics with solution-reaction energetics via a common tie to the Arrhenius equation, but gives mechanochemistry the unique opportunity to influence either half of the equation. In light of this, we discuss a strategy for translating solvent-based reaction conditions to ball milling conditions. Lastly, we posit that the extra control via frequency factor grants mechanochemistry the potential for greater selectivity than conventional solution reactions.

| S-EPMC5609516 | biostudies-literature

Modeling temperature entrainment of circadian clocks using the Arrhenius equation and a reconstructed model from Chlamydomonas reinhardtii.

Project description:Endogenous circadian rhythms allow living organisms to anticipate daily variations in their natural environment. Temperature regulation and entrainment mechanisms of circadian clocks are still poorly understood. To better understand the molecular basis of these processes, we built a mathematical model based on experimental data examining temperature regulation of the circadian RNA-binding protein CHLAMY1 from the unicellular green alga Chlamydomonas reinhardtii, simulating the effect of temperature on the rates by applying the Arrhenius equation. Using numerical simulations, we demonstrate that our model is temperature-compensated and can be entrained to temperature cycles of various length and amplitude. The range of periods that allow entrainment of the model depends on the shape of the temperature cycles and is larger for sinusoidal compared to rectangular temperature curves. We show that the response to temperature of protein (de)phosphorylation rates play a key role in facilitating temperature entrainment of the oscillator in Chlamydomonas reinhardtii. We systematically investigated the response of our model to single temperature pulses to explain experimentally observed phase response curves.

| S-EPMC3388193 | biostudies-literature

Polarizable molecular dynamics simulations of ionic liquids: Influence of temperature control.

Project description:Ionic liquids are an interesting class of soft matter with viscosities of one or two orders of magnitude higher than that of water. Unfortunately, classical, non-polarizable molecular dynamics (MD) simulations of ionic liquids result in too slow dynamics and demonstrate the need for explicit inclusion of polarizability. The inclusion of polarizability, here via the Drude oscillator model, requires amendments to the employed thermostat, where we consider a dual Nosé-Hoover thermostat, as well as a dual Langevin thermostat. We investigate the effects of the choice of a thermostat and the underlying parameters such as the masses and force constants of the Drude particles on static and dynamic properties of ionic liquids. Here, we show that Langevin thermostats are not suitable for investigating the dynamics of ionic liquids. Since polarizable MD simulations are associated with high computational costs, we employed a self-developed graphics processing unit enhanced code within the MD program CHARMM to keep the overall computational effort reasonable.

| S-EPMC7610910 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data