Pressure-Dependent Properties of Elementary Hydrophobic Interactions: Ramifications for Activation Properties of Protein Folding

Behavior of enzymes under high pressure in food processing: mechanisms, applications, and developments

High pressure processing (HPP) offers the benefits of safety, uniformity, energy-efficient, and low waste, which is widely applied for microbial inactivation and shelf-life extension for foods. Over the past forty years, HPP has been extensively researched in the food industry, enabling the inactivation or activation of different enzymes in future food by altering their molecular structure and active site conformation. Such activation or inactivation of enzymes effectively hinders the spoilage of food and the production of beneficial substances, which is crucial for improving food quality. This paper reviews the mechanism in which high pressure affects the stability and activity of enzymes, concludes the roles of key enzymes in the future food processed using high pressure technologies. Moreover, we discuss the application of modified enzymes based on high pressure, providing insights into the future direction of enzyme evolution under complex food processing conditions (e.g. high temperature, high pressure, high shear, and multiple elements). Finally, we conclude with prospects of high pressure technology and research directions in the future. Although HPP has shown positive effects in improving the future food quality, there is still a pressing need to develop new and effective combined processing methods, upgrade processing modes, and promote sustainable lifestyles.

Isothermal Compressibility Perturbation as a Protein Design Principle for T1 Lipase Stability-Activity Trade-Off Counteracting

Given the widely existing stability-activity trade-off in enzyme evolution, it is still a goal to obtain enzymes embracing both high activity and stability. Herein, we employed an isothermal compressibility (β_T) perturbation engineering (ICPE) strategy to comprehensively understand the stability-activity seesaw-like mechanism. The stability and activity of mutants derived from ICPE uncovered a high Pearson correlation (r = 0.93) in a prototypical enzyme T1 lipase. The best variant A186L/L188M/A190Y exhibited a high T_m value up to 78.70 °C, catalytic activity of 474.04 U/mg, and a 73.33% increase in dimethyl sulfoxide resistance compared to the wild type, one of the highest comprehensive performances reported to date. The elastic activation mechanism mediated by conformational change with a Δβ_T range of -6.81 × 10^-6 to -1.90 × 10^-6 bar^-1 may account for the balancing of stability and activity to achieve better performing enzymes. The ICPE strategy deepens our understanding of stability-activity trade-off and boosts its applications in enzyme engineering.

Behavior of Proteins under Pressure from Experimental Pressure-Dependent Structures

Structure-based models are coarse-grained representations of the interactions responsible for the protein folding process. In their simplest form, they use only the native contact map of a given protein to predict the main features of its folding process by computer simulation. Given their limitations, these models are frequently complemented with sequence-dependent contributions or additional information. Specifically, to analyze the effect of pressure on the folding/unfolding transition, special forms of these interaction potentials are employed, which may a priori determine the outcome of the simulations. In this work, we have tried to keep the original simplicity of structure-based models. Therefore, we have used folded structures that have been experimentally determined at different pressures to define native contact maps and thus interactions dependent on pressure. Despite the apparently tiny structural differences induced by pressure, our simulation results provide different thermodynamic and kinetic behaviors, which roughly correspond to experimental observations (when there is a possible comparison) of two proteins used as benchmarks, hen egg-white lysozyme and dihydrofolate reductase. Therefore, this work shows the feasibility of using experimental native structures at different pressures to analyze the global effects of this physical property on the protein folding process.

Implicit water model within the Zimm-Bragg approach to analyze experimental data for heat and cold denaturation of proteins

Studies of biopolymer conformations essentially rely on theoretical models that are routinely used to process and analyze experimental data. While modern experiments allow study of single molecules in vivo, corresponding theories date back to the early 1950s and require an essential update to include the recent significant progress in the description of water. The Hamiltonian formulation of the Zimm-Bragg model we propose includes a simplified, yet explicit model of water-polypeptide interactions that transforms into the equivalent implicit description after performing the summation of solvent degrees of freedom in the partition function. Here we show that our model fits very well to the circular dichroism experimental data for both heat and cold denaturation and provides the energies of inter- and intra-molecular H-bonds, unavailable with other processing methods. The revealed delicate balance between these energies determines the conditions for the existence of cold denaturation and thus clarifies its absence in some proteins.

Methane Clathrate Formation is Catalyzed and Kinetically Inhibited by the Same Molecule: Two Facets of Methanol

Here, we perform molecular dynamics simulations to provide atomic-level insights into the dual roles of methanol in enhancing and delaying the rate of methane clathrate hydrate nucleation. Consistent with experiments, we find that methanol slows clathrate hydrate nucleation above 250 K but promotes clathrate formation at temperatures below 250 K. We show that this behavior can be rationalized by the unusual temperature dependence of the methane-methanol interaction in an aqueous solution, which emerges due to the hydrophobic effect. In addition to its antifreeze properties at temperatures above 250 K, methanol competes with water to interact with methane prior to the formation of clathrate nuclei. Below 250 K, methanol encourages water to occupy the space between methane molecules favoring clathrate formation and it may additionally promote water mobility.

Biomolecular Condensates under Extreme Martian Salt Conditions

Biomolecular condensates formed by liquid-liquid phase separation (LLPS) are considered one of the early compartmentalization strategies of cells, which still prevail today forming nonmembranous compartments in biological cells. Studies of the effect of high pressures, such as those encountered in the subsurface salt lakes of Mars or in the depths of the subseafloor on Earth, on biomolecular LLPS will contribute to questions of protocell formation under prebiotic conditions. We investigated the effects of extreme environmental conditions, focusing on highly aggressive Martian salts (perchlorate and sulfate) and high pressure, on the formation of biomolecular condensates of proteins. Our data show that the driving force for phase separation of proteins is not only sensitively dictated by their amino acid sequence but also strongly influenced by the type of salt and its concentration. At high salinity, as encountered in Martian soil and similar harsh environments on Earth, attractive short-range interactions, ion correlation effects, hydrophobic, and π-driven interactions can sustain LLPS for suitable polypeptide sequences. Our results also show that salts across the Hofmeister series have a differential effect on shifting the boundary of immiscibility that determines phase separation. In addition, we show that confinement mimicking cracks in sediments and subsurface saline water pools in the Antarctica or on Mars can dramatically stabilize liquid phase droplets, leading to an increase in the temperature and pressure stability of the droplet phase.

Influence of Ionic Strength on Hydrophobic Interactions in Water: Dependence on Solute Size and Shape

Hydrophobicity is a phenomenon of great importance in biology, chemistry, and biochemistry. It is defined as the interaction between nonpolar molecules or groups in water and their low solubility. Hydrophobic interactions affect many processes in water, for example, complexation, surfactant aggregation, and coagulation. These interactions play a pivotal role in the formation and stability of proteins or biological membranes. In the present study, we assessed the effect of ionic strength, solute size, and shape on hydrophobic interactions between pairs of nonpolar particles. Pairs of methane, neopentane, adamantane, fullerene, ethane, propane, butane, hexane, octane, and decane were simulated by molecular dynamics in AMBER 16.0 force field. As a solvent, TIP3P and TIP4PEW water models were used. Potential of mean force (PMF) plots of these dimers were determined at four values of ionic strength, 0, 0.04, 0.08, and 0.40 mol/dm³, to observe its impact on hydrophobic interactions. The characteristic shape of PMFs with three extrema (contact minimum, solvent-separated minimum, and desolvation maximum) was observed for most of the compounds for hydrophobic interactions. Ionic strength affected hydrophobic interactions. We observed a tendency to deepen contact minima with an increase in ionic strength value in the case of spherical and spheroidal molecules. Additionally, two-dimensional distribution functions describing water density and average number of hydrogen bonds between water molecules were calculated in both water models for adamantane and hexane. It was observed that the density of water did not significantly change with the increase in ionic strength, but the average number of hydrogen bonds changed. The latter tendency strongly depends on the water model used for simulations.

The effects of cosolutes and crowding on the kinetics of protein condensate formation based on liquid-liquid phase separation: a pressure-jump relaxation study

Biomolecular assembly processes based on liquid-liquid phase separation (LLPS) are ubiquitous in the biological cell. To fully understand the role of LLPS in biological self-assembly, it is necessary to characterize also their kinetics of formation and dissolution. Here, we introduce the pressure-jump relaxation technique in concert with UV/Vis and FTIR spectroscopy as well as light microscopy to characterize the evolution of LLPS formation and dissolution in a time-dependent manner. As a model system undergoing LLPS we used the globular eye-lens protein γD-crystallin. As cosolutes and macromolecular crowding are known to affect the stability and dynamics of biomolecular condensates in cellulo, we extended our kinetic study by addressing also the impact of urea, the deep-sea osmolyte trimethylamine-N-oxide (TMAO) and a crowding agent on the transformation kinetics of the LLPS system. As a prerequisite for the kinetic studies, the phase diagram of γD-crystallin at the different solution conditions also had to be determined. The formation of the droplet phase was found to be a very rapid process and can be switched on and off on the 1-4 s timescale. Theoretical treatment using the Johnson-Mehl-Avrami-Kolmogorov model indicates that the LLPS proceeds via a diffusion-limited nucleation and growth mechanism at subcritical protein concentrations, a scenario which is also expected to prevail within biologically relevant crowded systems. Compared to the marked effect the cosolutes take on the stability of the LLPS region, their effect at biologically relevant concentrations on the phase transformation kinetics is very small, which might be a particular advantage in the cellular context, as a fast switching capability of the transition should not be compromised by the presence of cellular cosolutes.

The Cold-Unfolded State Is Expanded but Contains Long- and Medium-Range Contacts and Is Poorly Described by Homopolymer Models

Cold unfolding of proteins is predicted by the Gibbs-Helmholtz equation and is thought to be driven by a strongly temperature-dependent interaction of protein nonpolar groups with water. Studies of the cold-unfolded state provide insight into protein energetics, partially structured states, and folding cooperativity and are of practical interest in biotechnology. However, structural characterization of the cold-unfolded state is much less extensive than studies of thermally or chemically denatured unfolded states, in large part because the midpoint of the cold unfolding transition is usually below freezing. We exploit a rationally designed point mutation (I98A) in the hydrophobic core of the C-terminal domain of the ribosomal protein L9 that allows the cold denatured state ensemble to be observed above 0 °C at near neutral pH and ambient pressure in the absence of added denaturants. A combined approach consisting of paramagnetic relaxation enhancement measurements, analysis of small-angle X-ray scattering data, all-atom simulations, and polymer theory provides a detailed description of the cold-unfolded state. Despite a globally expanded ensemble, as determined by small-angle X-ray scattering, sequence-specific medium- and long-range interactions in the cold-unfolded state give rise to deviations from homopolymer-like behavior. Our results reveal that the cold-denatured state is heterogeneous with local and long-range intramolecular interactions that may prime the folded state and also demonstrate that significant long-range interactions are compatible with expanded unfolded ensembles. The work also highlights the limitations of homopolymer-based descriptions of unfolded states of proteins.

Pressure Sensitivity of SynGAP/PSD-95 Condensates as a Model for Postsynaptic Densities and Its Biophysical and Neurological Ramifications

Biomolecular condensates consisting of proteins and nucleic acids can serve critical biological functions, so that some condensates are referred as membraneless organelles. They can also be disease-causing, if their assembly is misregulated. A major physicochemical basis of the formation of biomolecular condensates is liquid-liquid phase separation (LLPS). In general, LLPS depends on environmental variables, such as temperature and hydrostatic pressure. The effects of pressure on the LLPS of a binary SynGAP/PSD-95 protein system mimicking postsynaptic densities, which are protein assemblies underneath the plasma membrane of excitatory synapses, were investigated. Quite unexpectedly, the model system LLPS is much more sensitive to pressure than the folded states of typical globular proteins. Phase-separated droplets of SynGAP/PSD-95 were found to dissolve into a homogeneous solution already at ten-to-hundred bar levels. The pressure sensitivity of SynGAP/PSD-95 is seen here as a consequence of both pressure-dependent multivalent interaction strength and void volume effects. Considering that organisms in the deep sea are under pressures up to about 1 kbar, this implies that deep-sea organisms have to devise means to counteract this high pressure sensitivity of biomolecular condensates to avoid harm. Intriguingly, these findings may shed light on the biophysical underpinning of pressure-related neurological disorders in terrestrial vertebrates.

Biomolecular Phase Separation: From Molecular Driving Forces to Macroscopic Properties

Biological phase separation is known to be important for cellular organization, which has recently been extended to a new class of biomolecules that form liquid-like droplets coexisting with the surrounding cellular or extracellular environment. These droplets are termed membraneless organelles, as they lack a dividing lipid membrane, and are formed through liquid-liquid phase separation (LLPS). Elucidating the molecular determinants of phase separation is a critical challenge for the field, as we are still at the early stages of understanding how cells may promote and regulate functions that are driven by LLPS. In this review, we discuss the role that disorder, perturbations to molecular interactions resulting from sequence, posttranslational modifications, and various regulatory stimuli play on protein LLPS, with a particular focus on insights that may be obtained from simulation and theory. We finally discuss how these molecular driving forces alter multicomponent phase separation and selectivity.

Critical phenomena in the temperature-pressure-crowding phase diagram of a protein

In the cell, proteins fold and perform complex functions through global structural rearrangements. Function requires a protein to be at the brink of stability to be susceptible to small environmental fluctuations, yet stable enough to maintain structural integrity. These apparently conflicting behaviors are exhibited by systems near a critical point, where distinct phases merge-a concept beyond previous studies indicating proteins have a well-defined folded/unfolded phase boundary in the pressure-temperature plane. Here, by modeling the protein phosphoglycerate kinase (PGK) on the temperature (T), pressure (P), and crowding volume-fraction (ϕ) phase diagram, we demonstrate a critical transition where phases merge, and PGK exhibits large structural fluctuations. Above the critical point, the difference between the intermediate and unfolded phases disappears. When ϕ increases, the critical point moves to lower T _c. We verify the calculations with experiments mapping the T-P-ϕ space, which likewise reveal a critical point at 305 K and 170 MPa that moves to lower T _c as ϕ increases. Crowding places PGK near a critical line in its natural parameter space, where large conformational changes can occur without costly free energy barriers. Specific structures are proposed for each phase based on simulation.

Temperature, Hydrostatic Pressure, and Osmolyte Effects on Liquid-Liquid Phase Separation in Protein Condensates: Physical Chemistry and Biological Implications

Liquid-liquid phase separation (LLPS) of proteins and other biomolecules play a critical role in the organization of extracellular materials and membrane-less compartmentalization of intra-organismal spaces through the formation of condensates. Structural properties of such mesoscopic droplet-like states were studied by spectroscopy, microscopy, and other biophysical techniques. The temperature dependence of biomolecular LLPS has been studied extensively, indicating that phase-separated condensed states of proteins can be stabilized or destabilized by increasing temperature. In contrast, the physical and biological significance of hydrostatic pressure on LLPS is less appreciated. Summarized here are recent investigations of protein LLPS under pressures up to the kbar-regime. Strikingly, for the cases studied thus far, LLPSs of both globular proteins and intrinsically disordered proteins/regions are typically more sensitive to pressure than the folding of proteins, suggesting that organisms inhabiting the deep sea and sub-seafloor sediments, under pressures up to 1 kbar and beyond, have to mitigate this pressure-sensitivity to avoid unwanted destabilization of their functional biomolecular condensates. Interestingly, we found that trimethylamine-N-oxide (TMAO), an osmolyte upregulated in deep-sea fish, can significantly stabilize protein droplets under pressure, pointing to another adaptive advantage for increased TMAO concentrations in deep-sea organisms besides the osmolyte's stabilizing effect against protein unfolding. As life on Earth might have originated in the deep sea, pressure-dependent LLPS is pertinent to questions regarding prebiotic proto-cells. Herein, we offer a conceptual framework for rationalizing the recent experimental findings and present an outline of the basic thermodynamics of temperature-, pressure-, and osmolyte-dependent LLPS as well as a molecular-level statistical mechanics picture in terms of solvent-mediated interactions and void volumes.

Temperature-Controlled Liquid-Liquid Phase Separation of Disordered Proteins

The liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) is a commonly observed phenomenon within the cell, and such condensates are also highly attractive for applications in biomaterials and drug delivery. A better understanding of the sequence-dependent thermoresponsive behavior is of immense interest as it will aid in the design of protein sequences with desirable properties and in the understanding of cellular response to heat stress. In this work, we use a transferable coarse-grained model to directly probe the sequence-dependent thermoresponsive phase behavior of IDPs. To achieve this goal, we develop a unique knowledge-based amino acid potential that accounts for the temperature-dependent effects on solvent-mediated interactions for different types of amino acids. Remarkably, we are able to distinguish between more than 35 IDPs with upper or lower critical solution temperatures at experimental conditions, thus providing direct evidence that incorporating the temperature-dependent solvent-mediated interactions to IDP assemblies can capture the difference in the shape of the resulting phase diagrams. Given the success of the model in predicting experimental behavior, we use it as a high-throughput screening framework to scan through millions of disordered sequences to characterize the composition dependence of protein phase separation.

Narrow equilibrium window for complex coacervation of tau and RNA under cellular conditions

The mechanism that leads to liquid-liquid phase separation (LLPS) of the tau protein, whose pathological aggregation is implicated in neurodegenerative disorders, is not well understood. Establishing a phase diagram that delineates the boundaries of phase co-existence is key to understanding whether LLPS is an equilibrium or intermediate state. We demonstrate that tau and RNA reversibly form complex coacervates. While the equilibrium phase diagram can be fit to an analytical theory, a more advanced model is investigated through field theoretic simulations (FTS) that provided direct insight into the thermodynamic driving forces of tau LLPS. Together, experiment and simulation reveal that tau-RNA LLPS is stable within a narrow equilibrium window near physiological conditions over experimentally tunable parameters including temperature, salt and tau concentrations, and is entropy-driven. Guided by our phase diagram, we show that tau can be driven toward LLPS under live cell coculturing conditions with rationally chosen experimental parameters.

Factors Promoting the Formation of Clathrate-Like Ordering of Water in Biomolecular Structure at Ambient Temperature and Pressure

Clathrate hydrate forms when a hydrophobic molecule is entrapped inside a water cage or cavity. Although biomolecular structures also have hydrophobic patches, clathrate-like water is found in only a limited number of biomolecules. Also, while clathrate hydrates form at low temperature and moderately higher pressure, clathrate-like water is observed in biomolecular structure at ambient temperature and pressure. These indicate presence of other factors along with hydrophobic environment behind the formation of clathrate-like water in biomolecules. In the current study, we presented a systematic approach to explore the factors behind the formation of clathrate-like water in biomolecules by means of molecular dynamics simulation of a model protein, maxi, which is a naturally occurring nanopore and has clathrate-like water inside the pore. Removal of either confinement or hydrophobic environment results in the disappearance of clathrate-like water ordering, indicating a coupled role of these two factors. Apart from these two factors, clathrate-like water ordering also requires anchoring groups that can stabilize the clathrate-like water through hydrogen bonding. Our results uncover crucial factors for the stabilization of clathrate-like ordering in biomolecular structure which can be used for the development of new biomolecular structure promoting clathrate formation.

Coarse-grained residue-based models of disordered protein condensates: utility and limitations of simple charge pattern parameters

Biomolecular condensates undergirded by phase separations of proteins and nucleic acids serve crucial biological functions. To gain physical insights into their genetic basis, we study how liquid-liquid phase separation (LLPS) of intrinsically disordered proteins (IDPs) depends on their sequence charge patterns using a continuum Langevin chain model wherein each amino acid residue is represented by a single bead. Charge patterns are characterized by the "blockiness" measure κ and the "sequence charge decoration" (SCD) parameter. Consistent with random phase approximation (RPA) theory and lattice simulations, LLPS propensity as characterized by critical temperature Tcr* increases with increasingly negative SCD for a set of sequences showing a positive correlation between κ and -SCD. Relative to RPA, the simulated sequence-dependent variation in Tcr* is often-though not always-smaller, whereas the simulated critical volume fractions are higher. However, for a set of sequences exhibiting an anti-correlation between κ and -SCD, the simulated Tcr*'s are quite insensitive to either parameter. Additionally, we find that blocky sequences that allow for strong electrostatic repulsion can lead to coexistence curves with upward concavity as stipulated by RPA, but the LLPS propensity of a strictly alternating charge sequence was likely overestimated by RPA and lattice models because interchain stabilization of this sequence requires spatial alignments that are difficult to achieve in real space. These results help delineate the utility and limitations of the charge pattern parameters and of RPA, pointing to further efforts necessary for rationalizing the newly observed subtleties.

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.

Pressure-Induced Dissolution and Reentrant Formation of Condensed, Liquid-Liquid Phase-Separated Elastomeric α-Elastin

We investigated the combined effects of temperature and pressure on liquid-liquid phase separation (LLPS) phenomena of α-elastin up to the multi-kbar regime. FT-IR spectroscopy, CD, UV/Vis absorption, phase-contrast light and fluorescence microscopy techniques were employed to reveal structural changes and mesoscopic phase states of the system. A novel pressure-induced reentrant LLPS was observed in the intermediate temperature range. A molecular-level picture, in particular on the role of hydrophobic interactions, hydration, and void volume in controlling LLPS phenomena is presented. The potential role of the LLPS phenomena in the development of early cellular compartmentalization is discussed, which might have started in the deep sea, where pressures up to the kbar level are encountered.

Molecular recognition and packing frustration in a helical protein

Biomolecular recognition entails attractive forces for the functional native states and discrimination against potential nonnative interactions that favor alternate stable configurations. The challenge posed by the competition of nonnative stabilization against native-centric forces is conceptualized as frustration. Experiment indicates that frustration is often minimal in evolved biological systems although nonnative possibilities are intuitively abundant. Much of the physical basis of minimal frustration in protein folding thus remains to be elucidated. Here we make progress by studying the colicin immunity protein Im9. To assess the energetic favorability of nonnative versus native interactions, we compute free energies of association of various combinations of the four helices in Im9 (referred to as H1, H2, H3, and H4) by extensive explicit-water molecular dynamics simulations (total simulated time > 300 μs), focusing primarily on the pairs with the largest native contact surfaces, H1-H2 and H1-H4. Frustration is detected in H1-H2 packing in that a nonnative packing orientation is significantly stabilized relative to native, whereas such a prominent nonnative effect is not observed for H1-H4 packing. However, in contrast to the favored nonnative H1-H2 packing in isolation, the native H1-H2 packing orientation is stabilized by H3 and loop residues surrounding H4. Taken together, these results showcase the contextual nature of molecular recognition, and suggest further that nonnative effects in H1-H2 packing may be largely avoided by the experimentally inferred Im9 folding transition state with native packing most developed at the H1-H4 rather than the H1-H2 interface.

Biomolecules need to recognize one another with high specificity: promoting “native” functional intermolecular binding events while avoiding detrimental “nonnative” bound configurations; i.e., “frustration”—the tendency for nonnative interactions—has to be minimized. Folding of globular proteins entails a similar discrimination. To gain physical insight, we computed the binding affinities of helical structures of the protein Im9 in various native or nonnative configurations by atomic simulations, discovering that partial packing of the Im9 core is frustrated. This frustration is overcome when the entire core of the protein is assembled, consistent with experiment indicating no significant kinetic trapping in Im9 folding. Our systematic analysis thus reveals a subtle, contextual aspect of biomolecular recognition and provides a general approach to characterize folding frustration.

Protein unfolding mechanisms and their effects on folding experiments

In this review, I discuss the various methods researchers use to unfold proteins in the lab in order to understand protein folding both in vitro and in vivo. The four main techniques, chemical-, heat-, pressure- and force-denaturation, produce distinctly different unfolded conformational ensembles. Recent measurements have revealed different folding kinetics from different unfolding mechanisms. Thus, comparing these distinct unfolded ensembles sheds light on the underlying free energy landscape of folding.

Theoretical Insights into the Biophysics of Protein Bi-stability and Evolutionary Switches

Deciphering the effects of nonsynonymous mutations on protein structure is central to many areas of biomedical research and is of fundamental importance to the study of molecular evolution. Much of the investigation of protein evolution has focused on mutations that leave a protein’s folded structure essentially unchanged. However, to evolve novel folds of proteins, mutations that lead to large conformational modifications have to be involved. Unraveling the basic biophysics of such mutations is a challenge to theory, especially when only one or two amino acid substitutions cause a large-scale conformational switch. Among the few such mutational switches identified experimentally, the one between the G_A all-α and G_B α+β folds is extensively characterized; but all-atom simulations using fully transferrable potentials have not been able to account for this striking switching behavior. Here we introduce an explicit-chain model that combines structure-based native biases for multiple alternative structures with a general physical atomic force field, and apply this construct to twelve mutants spanning the sequence variation between G_A and G_B. In agreement with experiment, we observe conformational switching from G_A to G_B upon a single L45Y substitution in the GA98 mutant. In line with the latent evolutionary potential concept, our model shows a gradual sequence-dependent change in fold preference in the mutants before this switch. Our analysis also indicates that a sharp G_A/G_B switch may arise from the orientation dependence of aromatic π-interactions. These findings provide physical insights toward rationalizing, predicting and designing evolutionary conformational switches.

The biological functions of globular proteins are intimately related to their folded structures and their associated conformational fluctuations. Evolution of new structures is an important avenue to new functions. Although many mutations do not change the folded state, experiments indicate that a single amino acid substitution can lead to a drastic change in the folded structure. The physics of this switch-like behavior remains to be elucidated. Here we develop a computational model for the relevant physical forces, showing that mutations can lead to new folds by passing through intermediate sequences where the old and new folds occur with varying probabilities. Our approach helps provide a general physical account of conformational switching in evolution and mutational effects on conformational dynamics.

Shedding light on the hydrophobicity puzzle

A general theory of hydrophobic hydration and pairwise hydrophobic interaction has been developed in the last years. The main ingredient is the recognition that: (a) cavity creation (necessary to insert a solute molecule into water) causes a solvent-excluded volume effect that leads to a loss in the translational entropy of water molecules; (b) the merging of two cavities (necessary to form the contact minimum configuration of two nonpolar molecules) causes a decrease in the solvent-excluded volume effect and so an increase in the translational entropy of water molecules. The performance of the theoretical approach is tested by reproducing both the hydration thermodynamics of xenon and the thermodynamics associated with the formation of the contact minimum configuration of two xenon atoms, over a large temperature range.

Steric interactions lead to collective tilting motion in the ribosome during mRNA-tRNA translocation

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNA–tRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/P–P/E and P/P–E/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/P–pe/E intermediate, where the 30S head is rotated ∼18°. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a ∼10° tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specific steric features are perturbed. These calculations show that head tilting may be attributed to specific steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

During protein elongation, the translocation of mRNA and tRNA molecules across the 30S ribosomal subunit is associated with large-scale motions of the 30S head domain. Here the authors carry out MD simulations to probe the associated steric interactions and identify novel tilting motions during the late stages of translocation.

Thermodynamics of Aβ16-21 dissociation from a fibril: Enthalpy, entropy, and volumetric properties

Here, we provide insights into the thermodynamic properties of A β16-21 dissociation from an amyloid fibril using all-atom molecular dynamics simulations in explicit water. An umbrella sampling protocol is used to compute potentials of mean force (PMF) as a function of the distance ξ between centers-of-mass of the A β16-21 peptide and the preformed fibril at nine temperatures. Changes in the enthalpy and the entropic energy are determined from the temperature dependence of these PMF(s) and the average volume of the simulation box is computed as a function of ξ. We find that the PMF at 310 K is dominated by enthalpy while the entropic energy does not change significantly during dissociation. The volume of the system decreases during dissociation. Moreover, the magnitude of this volume change also decreases with increasing temperature. By defining dock and lock states using the solvent accessible surface area (SASA), we find that the behavior of the electrostatic energy is different in these two states. It increases (unfavorable) and decreases (favorable) during dissociation in lock and dock states, respectively, while the energy due to Lennard-Jones interactions increases continuously in these states. Our simulations also highlight the importance of hydrophobic interactions in accounting for the stability of A β16-21.

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.

Native contact density and nonnative hydrophobic effects in the folding of bacterial immunity proteins

The bacterial colicin-immunity proteins Im7 and Im9 fold by different mechanisms. Experimentally, at pH 7.0 and 10°C, Im7 folds in a three-state manner via an intermediate but Im9 folding is two-state-like. Accordingly, Im7 exhibits a chevron rollover, whereas the chevron arm for Im9 folding is linear. Here we address the biophysical basis of their different behaviors by using native-centric models with and without additional transferrable, sequence-dependent energies. The Im7 chevron rollover is not captured by either a pure native-centric model or a model augmented by nonnative hydrophobic interactions with a uniform strength irrespective of residue type. By contrast, a more realistic nonnative interaction scheme that accounts for the difference in hydrophobicity among residues leads simultaneously to a chevron rollover for Im7 and an essentially linear folding chevron arm for Im9. Hydrophobic residues identified by published experiments to be involved in nonnative interactions during Im7 folding are found to participate in the strongest nonnative contacts in this model. Thus our observations support the experimental perspective that the Im7 folding intermediate is largely underpinned by nonnative interactions involving large hydrophobics. Our simulation suggests further that nonnative effects in Im7 are facilitated by a lower local native contact density relative to that of Im9. In a one-dimensional diffusion picture of Im7 folding with a coordinate- and stability-dependent diffusion coefficient, a significant chevron rollover is consistent with a diffusion coefficient that depends strongly on native stability at the conformational position of the folding intermediate.

In order to fold correctly, a globular protein must avoid being trapped in wrong, i.e., nonnative conformations. Thus a biophysical account of how attractive nonnative interactions are bypassed by some amino acid sequences but not others is key to deciphering protein structure and function. We examine two closely related bacterial immunity proteins, Im7 and Im9, that are experimentally known to fold very differently: Whereas Im9 folds directly, Im7 folds through a mispacked conformational intermediate. A simple model we developed accounts for their intriguingly different folding kinetics in terms of a balance between the density of native-promoting contacts and the hydrophobicity of local amino acid sequences. This emergent principle is extensible to other biomolecular recognition processes.

Theoretical perspectives on nonnative interactions and intrinsic disorder in protein folding and binding

The diverse biological functions of intrinsically disordered proteins (IDPs) have markedly raised our appreciation of protein conformational versatility, whereas the existence of energetically favorable yet functional detrimental nonnative interactions underscores the physical limitations of evolutionary optimization. Here we survey recent advances in using biophysical modeling to gain insight into experimentally observed nonnative behaviors and IDP properties. Simulations of IDP interactions to date focus mostly on coupled folding-binding, which follows essentially the same organizing principle as the local-nonlocal coupling mechanism in cooperative folding of monomeric globular proteins. By contrast, more innovative theories of electrostatic and aromatic interactions are needed for the conceptually novel but less-explored 'fuzzy' complexes in which the functionally bound IDPs remain largely disordered.