A fully atomistic computer simulation study of cold denaturation of a β-hairpin

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations

Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.

Coherent X-ray Scattering Reveals Nanoscale Fluctuations in Hydrated Proteins

Hydrated proteins undergo a transition in the deeply supercooled regime, which is attributed to rapid changes in hydration water and protein structural dynamics. Here, we investigate the nanoscale stress-relaxation in hydrated lysozyme proteins stimulated and probed by X-ray Photon Correlation Spectroscopy (XPCS). This approach allows us to access the nanoscale dynamics in the deeply supercooled regime (T = 180 K), which is typically not accessible through equilibrium methods. The observed stimulated dynamic response is attributed to collective stress-relaxation as the system transitions from a jammed granular state to an elastically driven regime. The relaxation time constants exhibit Arrhenius temperature dependence upon cooling with a minimum in the Kohlrausch-Williams-Watts exponent at T = 227 K. The observed minimum is attributed to an increase in dynamical heterogeneity, which coincides with enhanced fluctuations observed in the two-time correlation functions and a maximum in the dynamic susceptibility quantified by the normalized variance χ_T. The amplification of fluctuations is consistent with previous studies of hydrated proteins, which indicate the key role of density and enthalpy fluctuations in hydration water. Our study provides new insights into X-ray stimulated stress-relaxation and the underlying mechanisms behind spatiotemporal fluctuations in biological granular materials.

Comparison of On-the-Fly Probability Enhanced Sampling and Parallel Tempering Combined with Metadynamics for Atomistic Simulations of RNA Tetraloop Folding

Atomistic simulations with reliable models are extremely useful in providing exquisitely detailed pictures of biomolecular phenomena that are not always accessible to experiments. One such biomolecular phenomenon is RNA folding, which often requires exhaustive simulations with combined advanced sampling techniques. In this work, we employed the multithermal-multiumbrella on-the-fly probability enhanced sampling (MM-OPES) technique and compared it against combined parallel tempering and metadynamics simulations. We found that MM-OPES simulations were successful in reproducing the free energy surfaces from combined parallel tempering and metadynamics simulations. Importantly, we also investigated a wide range of temperature sets (minimum and maximum temperatures) for MM-OPES simulations in order to identify some guidelines for deciding the temperature limits for an accurate and efficient exploration of the free energy landscapes. We found that most temperature sets yielded almost the same accuracy in reproducing the free energy surface at the ambient conditions as long as (i) the maximum temperature is reasonably high, (ii) the temperature at which we run the simulation is reasonably high (in our simulations, running temperature is defined as [minimum temperature + maximum temperature]/2), and (iii) the effective sample size at the temperature of interest is statistically reasonable. In terms of the computational cost, all MM-OPES simulations were nearly 4 times less costly than the combined parallel tempering and metadynamics simulations. We concluded that the demanding combined parallel tempering and metadynamics simulations can safely be replaced with approximately 4 times less costly MM-OPES simulations (with carefully selected temperature limits) to obtain the same information.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Stress-dependent flexibility of a full-length human monoclonal antibody: Insights from molecular dynamics to support biopharmaceutical development

After several decades of advancements in drug discovery, product development of biopharmaceuticals remains a time- and resource-consuming endeavor. One of the main reasons is associated to the lack of fundamental understanding of conformational dynamics of such biologic entities, and how they respond to various stresses encountered during manufacturing. In this work, we have studied the conformational dynamics of human IgG₁κ b12 monoclonal antibody (mAb) using molecular dynamics simulations. The hundreds of nanoseconds long trajectories reveal that b12 mAb is highly flexible. Its variable domains show greater conformational fluctuations than the constant domains. Additionally, it collapses towards a more globular shape in response to thermal stress, leading to decrease in the total solvent exposed surface area and radius of gyration. This behavior is more pronounced for the deglycosylated b12 mAb, and it appears to correlate with increase in inter-domain contacts between specific regions of the antibody. Conformational fluctuations also cause temporary formation and disruption of hydrophobic and charged patches on the antibody surface, which is particularly important for the prediction of CMC properties during development phases of antibody-based biotherapeutics. The insights gained through these simulations may help the development of biologic drugs, especially with regards to manufacturing processes where antibodies may undergo significant thermal stress.

Protein Stability Depends Critically on the Surface Hydrogen-Bonding Network: A Case Study of Bid Protein

Understanding how proteins retain structural stability is not only of fundamental importance in biophysics but also critical to industrial production of antibodies and vaccines. Protein stability is known to depend mainly on two effects: internal hydrophobicity and H-bonding between the protein surface and solvent. A challenging task is to identify their individual contributions to a protein. Here, we investigate the structural stability of the apoptotic Bid protein in solutions containing various concentrations of guanidinium hydrochloride and urea using a combination of recently developed methods including the QTY (glutamine, threonine, and tyrosine) code and electron spin resonance-based peak-height analysis. We show that when the internal hydrophobicity of Bid is broken down using the QTY code, the surface H-bonding alone is sufficient to retain the structural stability intact. When the surface H-bonding is disrupted, Bid becomes sensitive to the temperature-dependent internal hydrophobicity such that it exhibits a reversible cold unfolding above water's freezing point. Using the combined approach, we show that the free-energy contributions of the two effects can be more reliably obtained. The surface H bonds are more important than the other effect in determining the structural stability of Bid protein.

Temperature effect on the SARS-CoV-2: A molecular dynamics study of the spike homotrimeric glycoprotein

Rapid spread of SARS-CoV-2 virus have boosted the need of knowledge about inactivation mechanisms to minimize the impact of COVID-19 pandemic. Recent studies have shown that SARS-CoV-2 virus can be disabled by heating, the exposure time for total inactivation depending on the reached temperature (e.g. more than 45 min at 329 K or less than 5 min at 373 K. In spite of recent crystallographic structures, little is known about the molecular changes induced by the temperature. Here, we unravel the molecular basis of the effect of the temperature over the SARS-CoV-2 spike glycoprotein, which is a homotrimer with three identical monomers, by executing atomistic molecular dynamics (MD) simulations at 298, 310, 324, 338, 358 and 373 K. Furthermore, both the closed down and open up conformational states, which affect the accessibility of receptor binding domain, have been considered. Our results suggest that the spike homotrimer undergoes drastic changes in the topology of the hydrogen bonding interactions and important changes on the secondary structure of the receptor binding domain (RBD), while electrostatic interactions (i.e. salt bridges) are mainly preserved. The proposed inactivation mechanism has important implications for engineering new approaches to fight the SARS-CoV-2 coronavirus, as for example, cleaving or reorganizing the hydrogen bonds through chaotropic agents or nanoparticles with local surface resonant plasmon effect.

Towards molecular movies with X-ray photon correlation spectroscopy

In this perspective article we highlight research opportunities and challenges in probing structural dynamics of molecular systems using X-ray Photon Correlation Spectroscopy (XPCS). The development of new X-ray sources, such as 4th generation storage rings and X-ray free-electron lasers (XFELs), provides promising new insights into molecular motion. Employing XPCS at these sources allows to capture a very broad range of timescales and lengthscales, spanning from femtoseconds to minutes and atomic scales to the mesoscale. Here, we discuss the scientific questions that can be addressed with these novel tools for two prominent examples: the dynamics of proteins in biomolecular condensates and the dynamics of supercooled water. Finally, we provide practical tips for designing and estimating feasibility of XPCS experiments as well as on detecting and mitigating radiation damage.

Hydrophilic Interactions Dominate the Inverse Temperature Dependence of Polypeptide Hydration Free Energies Attributed to Hydrophobicity

We address the association of the hydrophobic driving forces in protein folding with the inverse temperature dependence of protein hydration, wherein stabilizing hydration effects strengthen with increasing temperature in a physiological range. All-atom calculations of the free energy of hydration of aqueous deca-alanine conformers, holistically including backbone and side-chain interactions together, show that attractive peptide-solvent interactions and the thermal expansion of the solvent dominate the inverse temperature signatures that have been interpreted traditionally as the hydrophobic stabilization of proteins in aqueous solution. Equivalent calculations on a methane solute are also presented as a benchmark for comparison. The present study calls for a reassessment of the forces that stabilize folded protein conformations in aqueous solutions and of the additivity of hydrophobic/hydrophilic contributions.

Statistical and molecular dynamics (MD) simulation approach to investigate the role of intrinsically disordered regions of shikimate dehydrogenase in microorganisms surviving at different temperatures

Hyperthermophiles, a subset of prokaryotes that thrive in adverse temperatures, potentially utilize the protein molecular biosystem for maintaining thermostability in a wide range of temperatures. Recent studies revealed that these organisms have smaller proportions of intrinsically disordered proteins. In this study, we performed sequence and structural analysis to investigate the maintenance of protein conformation and their stability at different temperatures. The sequence analysis reveals the higher proportion of charged amino acids are responsible for preventing the helix formation and, hence, become disordered regions. For structural analysis, we chose shikimate dehydrogenase from four species, namely Listeria monocytogenes, Escherichia coli, Thermus thermophilus, and Methanopyrus kandleri, and evaluated the protein adaptation at 283 K, 300 K, and 395 K temperatures. From this investigation, we found more residues of shikimate dehydrogenase prefer an order-to-disorder transition at 395 K only for hyperthermophilic species. The solvent-accessible surface area (SASA) and hydrogen-bond analysis revealed that the tertiary conformation and the number of hydrogen bonds for hyperthermophilic shikimate dehydrogenase are highly preserved at 395 K, compared to 300 K. Our simulation results conjointly provide shikimate dehydrogenase of hyperthermophile which resists high temperatures through stronger protein tertiary conformations.

Pressure dependence of vibrational optical activity of model biomolecules. A computational study

Change of molecular properties with pressure is an attracting means to regulate molecular reactivity or biological activity. However, the effect is usually small and so far explored rather scarcely. To obtain a deeper insight and estimate the sensitivity of vibrational optical activity spectra to pressure-induced conformational changes, we investigate small model molecules. The Ala-Ala dipeptide, isomaltose disaccharide and adenine-uracil dinucleotide were chosen to represent three different biomolecular classes. The pressure effects were modeled by molecular dynamics and density functional theory simulations. The dinucleotide was found to be the most sensitive to the pressure, whereas for the disaccharide the smallest changes are predicted. Pressure-induced relative intensity changes in vibrational circular dichroism and Raman optical activity spectra are predicted to be 2-3-times larger than for non-polarized IR and Raman techniques.

Protein thermal stability

Proteins, in general, fold to a well-organized three-dimensional structure in order to function. The stability of this functional shape can be perturbed by external environmental conditions, such as temperature. Understanding the molecular factors underlying the resistance of proteins to the thermal stress has important consequences. First of all, it can aid the design of thermostable enzymes able to perform efficient catalysis in the high-temperature regime. Second, it is an essential brick of knowledge required to decipher the evolutionary pathways of life adaptation on Earth. Thanks to the development of atomistic simulations and ad hoc enhanced sampling techniques, it is now possible to investigate this problem in silico, and therefore provide support to experiments. After having described the methodological aspects, the chapter proposes an extended discussion on two problems. First, we focus on thermophilic proteins, a perfect model to address the issue of thermal stability and molecular evolution. Second, we discuss the issue of how protein thermal stability is affected by crowded in vivo-like conditions.

Differences in thermal structural changes and melting between mesophilic and thermophilic dihydrofolate reductase enzymes

The thermal resistance of two homolog enzymes is investigated, with an emphasis on their local stability and flexibility, and on the possible implications regarding their reactivity.

Molecular simulations by generalized-ensemble algorithms in isothermal-isobaric ensemble

Generalized-ensemble algorithms are powerful techniques for investigating biomolecules such as protein, DNA, lipid membrane, and glycan. The generalized-ensemble algorithms were originally developed in the canonical ensemble. On the other hand, not only temperature but also pressure is controlled in experiments. Additionally, pressure is used as perturbation to study relationship between function and structure of biomolecules. For this reason, it is important to perform efficient conformation sampling based on the isothermal-isobaric ensemble. In this article, we review a series of the generalized-ensemble algorithms in the isothermal-isobaric ensemble: multibaric-multithermal, pressure- and temperature-simulated tempering, replica-exchange, and replica-permutation methods. These methods achieve more efficient simulation than the conventional isothermal-isobaric simulation. Furthermore, the isothermal-isobaric generalized-ensemble simulation samples conformations of biomolecules from wider range of temperature and pressure. Thus, we can estimate physical quantities more accurately at any temperature and pressure values. The applications to the biomolecular system are also presented.

The Cold Thermal Response of an Amyloid Oligomer Differs from Typical Globular Protein Cold Denaturation

In contrast with the general behavior of folded proteins, the cold thermal response of amyloid assemblies is difficult to elicit with simple models. We exploit exhaustive simulations to evaluate the thermal response of a barrel-shaped model amyloid oligomer, with a distinct hydrophobic core akin to that of folded proteins. Cumulative thermal data over the range of 210-483 K indicate a sharp inflection and rise in structural stability as the temperature is decreased below the melting temperature of the water model. This is not commensurate with the equilibrium free energy profile obtained with core packing as the order parameter. However, energetic analyses and the size of their fluctuations indicate the crucial role of hydration in mediating structural transitions, beyond the expected temperature-dependent hydrophobic effect. Structural ordering of the hydration layer over bulk water is maximized at the transition and vanishes at high temperatures. This is a first direct demonstration of the microscopic influence of hydration water on the low-temperature response of an amyloid assembly close to the cryo-regime.

Computational Investigation of the Effect of Pressure on Protein Stability

Previous studies show parabolic or elliptical regions of protein stability in the pressure-temperature ( P, T) plane. The construction of stability diagrams requires accessing a sufficiently broad ( P, T) range, which is often frustrated by ice formation in experiments and sampling challenges in simulations. We perform a fully atomistic computational study of the miniprotein Trp-cage over the range of temperatures 210 ≤ T ≤ 420 K and pressures P ≤ 5 kbar and construct the corresponding stability diagram. At ambient temperature, pressure shifts the conformational states toward unfolding. Below 250 K, the native fold's stability depends nonmonotonically on pressure. While cold unfolding and thermal denaturation differ significantly at ambient pressure, they exhibit progressive similarity at elevated pressures. At ambient pressure, cold denaturation is an enthalpically driven process that preserves significant elements of Trp-cage's secondary structure. In contrast, cold unfolding at elevated pressures involves a more substantial loss of secondary and tertiary structure, similar to thermal denaturation.

Hydration shell differentiates folded and disordered states of a Trp-cage miniprotein, allowing characterization of structural heterogeneity by wide-line NMR measurements

Hydration properties of folded and unfolded/disordered miniproteins were monitored in frozen solutions by wide-line 1H-NMR. The amount of mobile water as function of T (-80 °C < T < 0 °C) was found characteristically different for folded (TC5b), semi-folded (pH < 3, TCb5(H+)) and disordered (TC5b_N1R) variants. Comparing results of wide-line 1H-NMR and molecular dynamics simulations we found that both the amount of mobile water surrounding proteins in ice, as well as their thaw profiles differs significantly as function of the compactness and conformational heterogeneity of their structure. We found that (i) at around -50 °C ~50 H2Os/protein melt (ii) if the protein is well-folded then this amount of mobile water remains quasi-constant up to -20 °C, (iii) if disordered then the quantity of the lubricating mobile water increases with T in a constant manner up to ~200 H2Os/protein by reaching -20 °C. Especially in the -55 °C ↔ -15 °C temperature range, wide-line 1H-NMR detects the heterogeneity of protein fold, providing the size of the hydration shell surrounding the accessible conformers at a given temperature. Results indicate that freezing of protein solutions proceeds by the gradual selection of the enthalpically most favored states that also minimize the number of bridging waters.

Influence of glycerol on the cooling effect of pair hydrophobicity in water: relevance to proteins’ stabilization at low temperature

Glycerol reduces the cooling effect of pair hydrophobicity (reduction of hydrophobicity with decreasing temperature) in water.

Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core

The primary driving force for protein folding is the formation of a well-packed, anhydrous core. However, recently, the crystal structure of an antifreeze protein, maxi, has been resolved where the core of the protein is filled with water, which apparently contradicts the existing notion of protein folding. Here, we have performed standard molecular dynamics (MD) simulation, replica exchange MD (REMD) simulation, and umbrella sampling using TIP4P water at various temperatures (300, 260, and 240 K) to explore the origin of this unusual structural feature. It is evident from standard MD and REMD simulations that the protein is found to be stable at 240 K in its unusual state. The core of protein has two layers of semi-clathrate water separating the methyl groups of alanine residues from different helical strands. However, with increasing temperature (260 and 300 K), the stability decreases as the core becomes dehydrated, and methyl groups of alanine are tightly packed driven by hydrophobic interactions. Calculation of the potential of mean force by an umbrella sampling technique between a pair of model hydrophobes resembling maxi protein at 240 K shows the stabilization of second solvent-separated minima (SSM), which provides a thermodynamic rationale of the unusual structural feature in terms of weakening of the hydrophobic interaction. Because the stabilization of SSMs is implicated for cold denaturation, it suggests that the maxi protein is so designed by nature where the cold denatured-like state becomes the biologically active form as it works near or below the freezing point of water.

How do glycerol and dimethyl sulphoxide affect local tetrahedral structure of water around a nonpolar solute at low temperature? Importance of preferential interaction

Glycerol and dimethyl sulphoxide (DMSO) have vital roles in cryoprotection of living cells, tissues, etc. The above action has been directly linked with disruption of hydrogen (H-) bond structure and dynamics of water by these cosolvents at bulk region and around various complex units, such as peptide, amino acid, protein, and lipid membrane. However, the disruption of the local structure of the water solvent around a purely hydrophobic solute is still not studied extensively. The latter is also important in the context of stabilization of protein from cold denaturation. Through all-atom molecular dynamics simulation, we have investigated the comparative effect of glycerol and DMSO on the orientational order of water around a nonpolar solute at -5 °C. A steady reduction of the tetrahedral order of water is observed at bulk (>10 Å distance from the solute) and solute interface (<5.5 Å distance from the solute) with increasing the cosolvent concentration. Contrasting roles of glycerol and DMSO have been evidenced. While DMSO affects the H-bond structure of the interfacial water more than that of the bulk water, glycerol affects the water structure almost uniformly at all regions around the solute. Furthermore, while glycerol helps to retain water molecules at the interface, DMSO significantly reduces the water content in that region. We have put forward a plausible mechanism for these contrasting roles of these cosolvents. The solute-cosolvent hydrophobic-interaction-induced orientational alignment of an interfacial cosolvent molecule determines whether the involvement of the cosolvent molecules in H-bonding with solvent water in the interface is akin to the bulk region or not.

Development of isothermal-isobaric replica-permutation method for molecular dynamics and Monte Carlo simulations and its application to reveal temperature and pressure dependence of folded, misfolded, and unfolded states of chignolin

We developed a two-dimensional replica-permutation molecular dynamics method in the isothermal-isobaric ensemble. The replica-permutation method is a better alternative to the replica-exchange method. It was originally developed in the canonical ensemble. This method employs the Suwa-Todo algorithm, instead of the Metropolis algorithm, to perform permutations of temperatures and pressures among more than two replicas so that the rejection ratio can be minimized. We showed that the isothermal-isobaric replica-permutation method performs better sampling efficiency than the isothermal-isobaric replica-exchange method and infinite swapping method. We applied this method to a β-hairpin mini protein, chignolin. In this simulation, we observed not only the folded state but also the misfolded state. We calculated the temperature and pressure dependence of the fractions on the folded, misfolded, and unfolded states. Differences in partial molar enthalpy, internal energy, entropy, partial molar volume, and heat capacity were also determined and agreed well with experimental data. We observed a new phenomenon that misfolded chignolin becomes more stable under high-pressure conditions. We also revealed this mechanism of the stability as follows: TYR2 and TRP9 side chains cover the hydrogen bonds that form a β-hairpin structure. The hydrogen bonds are protected from the water molecules that approach the protein as the pressure increases.

The fold preference and thermodynamic stability of α-synuclein fibrils is encoded in the non-amyloid-β component region

The strain-dependent synucleinopathies may be partially imprinted in the fold-dependent thermodynamic properties of non-amyloid-β component (NAC) fibrils.

Insilico direct folding of thrombin-binding aptamer G-quadruplex at all-atom level

The reversible folding of the thrombin-binding DNA aptamer G-quadruplexes (GQs) (TBA-15) starting from fully unfolded states was demonstrated using a prolonged time scale (10-12 μs) parallel tempering metadynamics (PTMetaD) simulation method in conjunction with a modified version of the AMBER bsc1 force field. For unbiased descriptions of the folding free energy landscape of TBA-15, this force field was minimally modified. From this direct folding simulation using the modified bsc1 force field, reasonably converged free energy landscapes were obtained in K+-rich aqueous solution (150 mM), providing detailed atomistic pictures of GQ folding mechanisms for TBA-15. This study found that the TBA folding occurred via multiple folding pathways with two major free energy barriers of 13 and 15 kcal/mol in the presence of several intermediate states of G-triplex variants. The early formation of these intermediates was associated with a single K+ ion capturing. Interestingly, these intermediate states appear to undergo facile transitions among themselves through relatively small energy barriers.

Water structure around hydrophobic amino acid side chain analogs using different water models

The water structure around hydrophobic groups governs various biochemical processes. There is an ongoing debate on whether water molecules near hydrophobic groups are more ordered with greater participation in water-water hydrogen bonding with respect to water in the pure bulk state. The water structure around six different hydrophobic amino acid side chain analog molecules has been studied in pure water using molecular dynamics simulations. The analysis of water tetrahedral order parameter and the number of hydrogen bonds formed by the individual water molecules in the first hydration shell of the hydrophobic analogs provide evidence that both ordering and hydrogen bonds involving water molecules are to some extent reduced in the hydrophobic hydration shell. It is revealed that the water tetrahedrality in the outer part of the first hydrophobic hydration shell is equivalent to bulk water for all the water models except for the TIP4P-2005 model which shows marginally higher tetrahedrality. However, irrespective of the model employed, water tetrahedrality has always been found to be reduced in the inner part of the first hydration shell, which eventually makes the overall water tetrahedrality in the first hydrophobic hydration shell marginally lower than that observed for pure bulk water. Importantly, it is noticed that the decrease in water structuring exhibits solute size dependencies. Around a small solute like methane, the water tetrahedral ordering or hydrogen bonding propensity is quite similar to that of the bulk state. The effect, reduction in water structuring, is however more pronounced for relatively larger solutes.

Effects of high pressure freezing (HPF) on denaturation of natural actomyosin extracted from prawn (Metapenaeus ensis)

Effects of protein denaturation caused by high pressure freezing, involving Pressure-Factors (pressure, time) and Freezing-Factors (temperature, phase transition, recrystallization, ice crystal types), are complicated. In the current study, the conformation and functional changes of natural actomyosin (NAM) under pressure assisted freezing (PAF, ^{100,150,300,400,500MPa}P_-20°C/25min), pressure shift freezing (PSF, ^200MPaP_-20°C/25min), and immersion freezing (^0.1MPaP_-20°C/5min) after pressure was released to 0.1MPa, as compared to normal immersion freezing process (IF, ^0.1MPaP_-20°C/30min). Results indicated that PSF (^200MPaP_-20°C/30min) could reduce the denaturation of frozen NAM and a pressure of 300MPa was the critical point to induce such a denaturation. During the periods of B→D in PSF or B→C→D in PAF, the generation and growth of ice crystals played an important role on changing the secondary and tertiary structure of the treated NAM.

Comparison of force fields for Alzheimer's A β42: A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three-dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN-NMR and CHARMM22*, in their ability to model Aβ42 , an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN-NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.

Computational investigation of cold denaturation in the Trp-cage miniprotein

The functional native states of globular proteins become unstable at low temperatures, resulting in cold unfolding and impairment of normal biological function. Fundamental understanding of this phenomenon is essential to rationalizing the evolution of freeze-tolerant organisms and developing improved strategies for long-term preservation of biological materials. We present fully atomistic simulations of cold denaturation of an α-helical protein, the widely studied Trp-cage miniprotein. In contrast to the significant destabilization of the folded structure at high temperatures, Trp-cage cold denatures at 210 K into a compact, partially folded state; major elements of the secondary structure, including the α-helix, are conserved, but the salt bridge between aspartic acid and arginine is lost. The stability of Trp-cage's α-helix at low temperatures suggests a possible evolutionary explanation for the prevalence of such structures in antifreeze peptides produced by cold-weather species, such as Arctic char. Although the 310-helix is observed at cold conditions, its position is shifted toward Trp-cage's C-terminus. This shift is accompanied by intrusion of water into Trp-cage's interior and the hydration of buried hydrophobic residues. However, our calculations also show that the dominant contribution to the favorable energetics of low-temperature unfolding of Trp-cage comes from the hydration of hydrophilic residues.

Contribution of Water to Pressure and Cold Denaturation of Proteins

The mechanisms of cold and pressure denaturation of proteins are matter of debate and are commonly understood as due to water-mediated interactions. Here, we study several cases of proteins, with or without a unique native state, with or without hydrophilic residues, by means of a coarse-grain protein model in explicit solvent. We show, using Monte Carlo simulations, that taking into account how water at the protein interface changes its hydrogen bond properties and its density fluctuations is enough to predict protein stability regions with elliptic shapes in the temperature-pressure plane, consistent with previous theories. Our results clearly identify the different mechanisms with which water participates to denaturation and open the perspective to develop advanced computational design tools for protein engineering.

On the Effect of Sodium Chloride and Sodium Sulfate on Cold Denaturation

Both sodium chloride and sodium sulfate are able to stabilize yeast frataxin, causing an overall increase of its thermodynamic stability curve, with a decrease in the cold denaturation temperature and an increase in the hot denaturation one. The influence of low concentrations of these two salts on yeast frataxin stability can be assessed by the application of a theoretical model based on scaled particle theory. First developed to figure out the mechanism underlying cold denaturation in water, this model is able to predict the stabilization of globular proteins provided by these two salts. The densities of the salt solutions and their temperature dependence play a fundamental role.