Coarse-grained strategy for modeling protein stability in concentrated solutions. II: phase behavior

We use highly efficient transition-matrix Monte Carlo simulations to determine equilibrium unfolding curves and fluid phase boundaries for solutions of coarse-grained globular proteins. The model we analyze derives the intrinsic stability of the native state and protein-protein interactions from basic information about protein sequence using heteropolymer collapse theory. It predicts that solutions of low hydrophobicity proteins generally exhibit a single liquid phase near their midpoint temperatures for unfolding, while solutions of proteins with high sequence hydrophobicity display the type of temperature-inverted, liquid-liquid transition associated with aggregation processes of proteins and other amphiphilic molecules. The phase transition occurring in solutions of the most hydrophobic protein we study extends below the unfolding curve, creating an immiscibility gap between a dilute, mostly native phase and a concentrated, mostly denatured phase. The results are qualitatively consistent with the solution behavior of hemoglobin (HbA) and its sickle variant (HbS), and they suggest that a liquid-liquid transition resulting in significant protein denaturation should generally be expected on the phase diagram of high-hydrophobicity protein solutions. The concentration fluctuations associated with this transition could be a driving force for the nonnative aggregation that can occur below the midpoint temperature.     

Protein stability and dynamics in the pressure-temperature plane

The pressure-temperature stability diagram of proteins and the underlying assumptions of the elliptical shape of the diagram are discussed. Possible extensions, such as aggregation and fibril formation, are considered. An important experimental observation is the extreme pressure stability of the mature fibrils. Molecular origins of the diagram in terms of models of the partial molar volume of a protein focus on cavities and hydration. Changes in thermal expansivity, compressibility and heat capacity in terms of fluctuations of the enthalpy and volume change of the unfolding should also focus on these parameters. It is argued that the study of water-soluble polymers might further our understanding of the stability diagram. Whereas the role of water in protein behaviour is unquestioned, the role of cavities is less clear.     

Expanding the pressure technique: insights into protein folding from combined use of pressure and chemical denaturants

The fundamental principles derived from in vitro protein folding experiments have practical application in understanding the pathology of diseases of protein misfolding and for the development of industrial processes to produce proteins as pharmaceuticals and biotechnological reagents. High pressure as a tool to denature or disaggregate proteins offers a number of unique advantages. The emphasis of this review is on how low concentrations of chemical denaturants can be used in combination with high pressure to extend the range and scope of this useful technique. This approach has already been used in a number of studies, which are discussed here in the context of the questions they address. These include: the origin of the volume change observed on protein unfolding, pressure-induced formation of partially structured intermediates, pressure-induced dissociation of oligomeric and aggregated proteins, and the use of volume changes to probe the structure of the transition state. Wider use of hydrostatic pressure as a denaturation tool, facilitated by combination with chemical denaturants, is likely to bring significant advances to our understanding of protein structure, stability and folding, particularly in relation to proteins associated with the amyloid and prion diseases.     

Conformational fluctuations of proteins revealed by variable pressure NMR

With the high-resolution variable-pressure NMR spectroscopy, one can study conformational fluctuations of proteins in a much wider conformational space than hitherto explored by NMR and other spectroscopic techniques. This is because a protein in solution generally exists as a dynamic mixture of conformers mutually differing in partial molar volume, and pressure can select the population of a conformer according to its relative volume. In this review, we describe how variable-pressure NMR can be used to probe conformational fluctuations of proteins in a wide conformational space from the folded to the fully unfolded structures, with actual examples. Furthermore, the newly emerging technique "NMR snapshots" expresses amply fluctuating protein structures as changes in atomic coordinates. Finally, the concept of conformational fluctuation is extended to include intermolecular association leading to amyloidosis.     

How molecular modelling can better broaden the understanding of glycosylations

Glycosylations are among the most ubiquitous post-translational modifications (PTMs) in proteins, and the effects of their perturbations are seen in various diseases such as cancers, diabetes and arthritis to name a few. Yet they remain one of the most enigmatic aspects of protein structure and function. On the other hand, molecular modelling techniques have been rapidly bridging this knowledge gap since the last decade. In this review, we discuss how these techniques have proven to be indispensable for a better understanding of the role of glycosylations in glycoprotein structure and function.     

The assembly of ionic currents in a thalamic neuron. III. The seven-dimensional model

We replace our earlier three-dimensional ha-model of a thalamic neuron (Rose & Hindmarsh, Proc. R. Soc. Lond. B 237, 289-312 (1989b)) by a seven-dimensional ha-model. The stability and state diagrams for this seven-dimensional model are shown to be similar to those of the three-dimensional system with which it is compared. Two examples that illustrate how the seven-dimensional model can be related to experimental recordings are then discussed in detail. In each case we show the state and stability diagrams and the time courses of the different ionic currents during the burst response. We discuss how the various components of the state diagram could be determined experimentally. These experiments are illustrated by using our model to simulate the results of actual experiments. Finally the state and stability diagrams are transferred to a current-voltage diagram. This shows the advantage of the state diagram compared with the current-voltage diagram for representing the dynamics of the firing of a thalamic neuron.     

Melting transitions in biomembranes

We investigated melting transitions in native biological membranes containing their membrane proteins. The membranes originated from E. coli, B. subtilis, lung surfactant and nerve tissue from the spinal cord of several mammals. For some preparations, we studied the pressure, pH and ionic strength dependence of the transition. For porcine spine, we compared the transition of the native membrane to that of the extracted lipids. All preparations displayed melting transitions of 10–20° below physiological or growth temperature, independent of the organism of origin and the respective cell type. We found that the position of the transitions in E. coli membranes depends on the growth temperature.We discuss these findings in the context of the thermodynamic theory of membrane fluctuations close to transition that predicts largely altered elastic constants, an increase in fluctuation lifetime and in membrane permeability. We also discuss how to distinguish lipid melting from protein unfolding transitions. Since the feature of a transition slightly below physiological temperature is conserved even when growth conditions change, we conclude that the transitions are likely to be of major biological importance for the survival and the function of the cell.     

The emergence of protein complexes: quaternary structure, dynamics and allostery. Colworth Medal Lecture

All proteins require physical interactions with other proteins in order to perform their functions. Most of them oligomerize into homomers, and a vast majority of these homomers interact with other proteins, at least part of the time, forming transient or obligate heteromers. In the present paper, we review the structural, biophysical and evolutionary aspects of these protein interactions. We discuss how protein function and stability benefit from oligomerization, as well as evolutionary pathways by which oligomers emerge, mostly from the perspective of homomers. Finally, we emphasize the specificities of heteromeric complexes and their structure and evolution. We also discuss two analytical approaches increasingly being used to study protein structures as well as their interactions. First, we review the use of the biological networks and graph theory for analysis of protein interactions and structure. Secondly, we discuss recent advances in techniques for detecting correlated mutations, with the emphasis on their role in identifying pathways of allosteric communication.     

Characterization of the pressure-induced intermediate and unfolded state of red-shifted green fluorescent protein--a static and kinetic FTIR,UV/VIS and fluorescence spectroscopy study

The green fluorescence proteins (GFP) are widely used as reporters in molecular and cell biology. For their use it in high-pressure microbiology and biotechnology studies, their structural properties, thermodynamic parameters and stability diagrams have to be known. We investigated the pressure stability of the red-shifted green fluorescent protein (rsGFP) using Fourier-transform infrared spectroscopy, fluorescence and UV/Vis spectroscopy. We found that rsGFP does not unfold up to approximately 9kbar at room temperature. Its unique three-dimensional structure is held responsible for the high-pressure stability. At higher temperatures, its secondary structure collapses below 9kbar (e.g. the denaturation pressure at 58 degrees C is 7.8kbar). The analysis of the IR data shows that the pressure-denatured state contains more disordered structures at the expense of a decrease of intramolecular beta-sheets. As indicated by the large volume change of DeltaV degrees (u) approximately -250(+/-50)mlmol(-1) at 58 degrees C, this highly cooperative transition can be interpreted as a collapse of the beta-can structure of rsGFP. For comparison, the temperature-induced unfolding of rsGFP has also been studied. At high temperature (T(m)=78 degrees C), the unfolding resulted in the formation of an aggregated state. Contrary to the pressure-induced unfolding, the temperature-induced unfolding and aggregation of GFP is irreversible. From the FT-IR data, a tentative p,T-stability diagram for the secondary structure collapse of GFP has been obtained. Furthermore, changes in fluorescence and absorptivity were found which are not correlated to the secondary structural changes. The fluorescence and UV/Vis data indicate smaller conformational changes in the chromophore region at much lower pressures ( approximately 4kbar) which are probably accompanied by the penetration of water into the beta-can structure. In order to investigate also the kinetics of this initial step, pressure-jump relaxation experiments were carried out. The partial activation volumes observed indicate that the conformational changes in the chromophore region when passing the transition state are indeed rather small, thus leading to a comparably small volume change of -20 ml mol(-1) only. The use of the chromophore absorption and fluorescence band of rsGFP in using GFP as reporter for gene expression and other microbiological studies under high pressure conditions is thus limited to pressures of about 4kbar, which still exceeds the pressure range relevant for studies in vivo in micro-organisms, including piezophilic bacteria from deep-sea environments.     

Compressibility of protein transitions

We review the results of compressibility studies on proteins and low molecular weight compounds that model the hydration properties of these biopolymers. In particular, we present an analysis of compressibility changes accompanying conformational transitions of globular proteins. This analysis, in conjunction with experimental compressibility data on protein transitions, were used to define the changes in the hydration properties and intrinsic packing associated with native-to-molten globule, native-to-partially unfolded, and native-to-fully unfolded transitions of globular proteins. In addition, we discuss the molecular origins of predominantly positive changes in compressibility observed for pressure-induced denaturation transitions of globular proteins. Throughout this review, we emphasize the importance of compressibility data for characterizing protein transitions, while also describing how such data can be interpreted to gain insight into role that hydration and intrinsic packing play in modulating the stability of and recognition between proteins and other biologically important compounds.     

Experimental and theoretical studies of mechanical unfolding of different proteins

Mechanical properties of proteins are important for a wide range of biological processes including cell adhesion, muscle contraction, and protein translocation across biological membranes. It is necessary to reveal how proteins achieve their required mechanical stability under natural conditions in order to understand the biological processes and also to use the knowledge for constructing new biomaterials for medical and industrial purposes. In this connection, it is important to know how a protein will behave in response to various impacts. Theoretical and experimental works on mechanical unfolding of globular proteins will be considered in detail in this review.     

Characterization of the distribution of internal motions in the basic pancreatic trypsin inhibitor using a large number of internal NMR probes

The experimental observations described in this article indicated that a distribution of many different fluctuations is present in a globular protein. These fluctuations were characterized by observation of many natural internal probes such as the labile peptide protons and the aromatic side chains. The conditions which are necessary to get reactions of the internal probes have been discussed in detail. The structural interpretation of the data was facilitated by the development and the use of new NMR techniques which provided the identification of the resonances of all the labile peptide protons. With NOE measurements a distinction between correlated and uncorrelated exchange events was obtained. This enabled us to elucidate the exchange mechanism over a wide range of p2H and temperature and to classify different subsets of fluctuations with respect to their lifetimes. It was further demonstrated that a change of external conditions such as temperature, p2H or pressure can change the distribution of fluctuations in the protein. The mechanisms responsible for rotation of internal aromatic side chains were also found to change with temperature, and mechanistic aspects of these fluctuations were discussed. This demonstration of a manifold of spatial fluctuations in a small protein provides an impression on the kind of fluctuations which have to be expected for larger proteins. When studying protein reactions one should therefore consider the presence of a large number of different, transiently formed, spatial structures available for the partner in the reaction, which may pick out only that structure which will optimally perform a particular reaction with the highest efficiency.     

Understanding the role of hydrogen bonds in water dynamics and protein stability

The mechanisms of cold and pressure denaturation of proteins are a matter of debate, but it is commonly accepted that water plays a fundamental role in the process. It has been proposed that the denaturation process is related to an increase of hydrogen bonds among hydration water molecules. Other theories suggest that the causes of denaturation are the density fluctuations of surface water, or the destabilization of hydrophobic contacts as a consequence of water molecule inclusions inside the protein, especially at high pressures. We review some theories that have been proposed to give insight into this problem, and we describe a coarse-grained model of water that compares well with experiments for proteins’ hydration water. We introduce its extension for a homopolymer in contact with the water monolayer and study it by Monte Carlo simulations in an attempt to understand how the interplay of water cooperativity and interfacial hydrogen bonds affects protein stability.     

Exploring the temperature-pressure phase diagram of staphylococcal nuclease

The temperature dependence of the pressure-induced equilibrium unfolding of staphylococcal nuclease (Snase) was determined by fluorescence of the single tryptophan residue, FTIR absorption for the amide I' and tyrosine O-H bands, and small-angle X-ray scattering (SAXS). The results from these three techniques were similar, although the stability as measured by fluorescence was slightly lower than that measured by FTIR and SAXS. The resulting phase diagram exhibits the well-known curvature for heat and cold denaturation of proteins, due to the large decrease in heat capacity upon folding. The volume change for unfolding became less negative with increasing temperatures, consistent with a larger thermal expansivity for the unfolded state than for the folded state. Fluorescence-detected pressure-jump kinetics measurements revealed that the curvature in the phase diagram is due primarily to the rate constant for folding, indicating a loss in heat capacity for the transition state relative to the unfolded state. The similar temperature dependence of the equilibrium and activation volume changes for folding indicates that the thermal expansivities of the folded and transition states are similar. This, along with the fact that the activation volume for folding is positive over the temperature range examined, the nonlinear dependence of the folding rate constant upon temperature implicates significant dehydration in the rate-limiting step for folding of Snase.     

Interaction of the antimicrobial peptide gramicidin S with dimyristoyl--phosphatidylcholine bilayer membranes: a densitometry and sound velocimetry study

We determined changes in the volume and adiabatic compressibility of large multi- and unilamellar vesicles composed of dimyristoylphosphatidylcholine containing various concentrations of the antimicrobial peptide gramicidin S (GS) by applying densitometry and sound velocimetry. Gramicidin S incorporation was found to progressively decrease the phase transition temperature of DMPC vesicles as well as to decrease the degree of cooperativity of the main phase transition and to increase the volume compressibility of the vesicles. GS probably enhanced thermal fluctuations at the region of main phase transition and provide more freedom of rotational movement for the phospholipid hydrocarbon chains. The ability of GS to increase the membrane compressibility and to decrease the phase transition temperature is evidence for regions of distorted membrane structure around incorporated gramicidin S molecules. At relatively high GS concentration (10 mol%), more significant changes of specific volume and compressibility appear. This might suggest changes in the integrity of the lipid bilayer upon interaction with high concentrations of GS.     

Phase diagram for cell cytoplasm from the calf lens

We have determined the phase diagram for cytoplasmic homogenate isolated from the cells of calf lenses. In the intact lens, this material is associated with the reversible opacification known as cold cataract. The coexistence curve we find resembles that found in concentrated high polymer solutions. Two coexisting phases are found below the transition temperature. Both these phases have similar protein compositions but different net protein concentrations. The constituent proteins have the following molecular weights: 17,000, 20,000, 22,000, 26,000, 28,000, 31,000, 36,000, 52,000 and 55,000 daltons. Scanning electron microscopy shows that the phase separation in the cytoplasmic homogenate is accompanied by a subtle structural reorganization which consists of increased spatial density fluctuations of dimensions comparable to the wavelength of light.     

Water at Interface with Proteins

Water is essential for the activity of proteins. However, the effect of the properties of water on the behavior of proteins is only partially understood. Recently, several experiments have investigated the relation between the dynamics of the hydration water and the dynamics of protein. These works have generated a large amount of data whose interpretation is debated. New experiments measure the dynamics of water at low temperature on the surface of proteins, finding a qualitative change (crossover) that might be related to the slowing down and stop of the protein’s activity (protein glass transition), possibly relevant for the safe preservation of organic material at low temperature. To better understand the experimental data several scenarios have been discussed. Here, we review these experiments and discuss their interpretations in relation with the anomalous properties of water. We summarize the results for the thermodynamics and dynamics of supercooled water at an interface. We consider also the effect of water on protein stability, making a step in the direction of understanding, by means of Monte Carlo simulations and theoretical calculations, how the interplay of water cooperativity and hydrogen bonds interfacial strengthening affects the protein cold denaturation.     

Effects of glycerol on the compaction and stability of the wild type and mutated rabbit muscle creatine kinase

All organisms have developed detect, repair, regulation, and stabilization mechanisms to survive from cellular and molecular damage induced by diverse stresses. Among them, the accumulation of osmolytes is a common mechanism evolved by cells to maintain cell volume and stabilize macromolecules against various environmental stresses. The molecular mechanisms by which osmolytes stabilize proteins and prevent aggregation have been well-established. However, little is known about the effects of osmolytes on mutated or damaged proteins. In this research, we investigated the effects of glycerol on the activity, structure, and stability of the wild type (WT) and D54G CK under normal and extreme (high temperature) conditions. It was found that glycerol had similar effects on the suppression of the aggregation during the refolding of both proteins. Under native conditions, the effect of glycerol on the mutated protein was more obvious than on the WT protein. Glycerol could efficiently force the mutated protein to fold to a state close to the WT protein, and thus stabilize the native state of the mutated protein. Glycerol could also protect both the WT and mutated proteins against heat-induced denaturation. However, the change in the transition free energy of heat-induced inactivation of the WT protein was larger than that of the mutated protein. These results suggested that glycerol might have differential effects on the changes of the chemical potential and the transition free energy of the WT and mutated proteins.     

Pressure-temperature phase diagrams of biomolecules

The pressure-temperature phase diagram of various biomolecules is reviewed. Special attention is focused on the elliptic phase diagram of proteins. The phenomenological thermodynamic theory describing this diagram explains the heat, cold and pressure denaturations in a unified picture. The limitations and possible developments of this theory are discussed as well. It is pointed out that a more complex diagram can be obtained when the intermolecular interactions are also taken into account. In this case metastable states appear on the pressure-temperature (p-T) diagram due to intermolecular interactions. Pressure-temperature phase diagrams of other biopolymers are also discussed. While the p-T diagrams of helix-coil transition of nucleic acids and of gel-liquid crystal transition of lipid bilayers are non-elliptical, those of gelatinization of starch and of phase separation of some synthetic polymers show an elliptic profile, similar to that of proteins. Finally, the p-T diagram of bacterial inactivation is shown to be elliptic. From the point of view of basic science, this fact shows that the key factor of inactivation should be the protein type, and from the viewpoint of practical applications, it serves as the theoretical basis of pressure treatment of biosystems.