Thermodynamics of denaturation of mutants of barnase with disulfide crosslinks

Misfolding of a Single Disulfide Bonded Globular Protein into a Low-Solubility Species Conformationally and Biophysically Distinct from the Native One

In practice and despite Anfinsen’s dogma, the refolding of recombinant multiple SS-bonded proteins is famously difficult because misfolded species with non-native SS-bonds appear upon the oxidization of their cysteine residues. On the other hand, single SS-bond proteins are thought to be simple to refold because their cysteines have only one SS-bond partner. Here, we report that dengue 4 envelope protein domain 3 (DEN4 ED3), a single SS-bonded protein can be irreversibly trapped into a misfolded species through the formation of its sole intramolecular SS-bond. The misfolded species had a much lower solubility than the native one at pHs higher than about 7, and circular dichroism measurements clearly indicated that its secondary structure content was different from the native species. Furthermore, the peaks in the Heteronuclear Single Quantum Correlation spectroscopy (HSQC) spectrum of DEN4 ED3 from the supernatant fraction were sharp and well dispersed, reflecting the beta-sheeted native structure, whereas the spectrum of the precipitated fraction showed broad signals clustered near its center suggesting no or little structure and a strong tendency to aggregate. The two species had distinct biophysical properties and could interconvert into each other only by cleaving and reforming the SS-bond, strongly suggesting that they are topologically different. This phenomenon can potentially happen with any single SS-bonded protein, and our observation emphasizes the need for assessing the conformation and biophysical properties of bacterially produced therapeutic proteins in addition to their chemical purities.

"Cooperative collapse" of the denatured state revealed through Clausius-Clapeyron analysis of protein denaturation phase diagrams

Protein phase diagrams have a unique potential to identify the presence of additional thermodynamic states even when non-2-state character is not readily apparent from the experimental observables used to follow protein unfolding transitions. Two-state analysis of the von Willebrand factor A3 domain has previously revealed a discrepancy in the calorimetric enthalpy obtained from thermal unfolding transitions as compared with Gibbs-Helmholtz analysis of free energies obtained from the Linear Extrapolation Method (Tischer and Auton, Prot Sci 2013; 22(9):1147-60). We resolve this thermodynamic conundrum using a Clausius-Clapeyron analysis of the urea-temperature phase diagram that defines how Δ H and the urea m-value interconvert through the slope of c_m versus T, ( ∂ c m / ∂ T ) = Δ H / ( m T ) . This relationship permits the calculation of Δ H at low temperature from m-values obtained through iso-thermal urea denaturation and high temperature m-values from Δ H obtained through iso-urea thermal denaturation. Application of this equation uncovers sigmoid transitions in both cooperativity parameters as temperature is increased. Such residual thermal cooperativity of Δ H and the m-value confirms the presence of an additional state which is verified to result from a cooperative phase transition between urea-expanded and thermally-compact denatured states. Comparison of the equilibria between expanded and compact denatured ensembles of disulfide-intact and carboxyamidated A3 domains reveals that introducing a single disulfide crosslink does not affect the presence of the additional denatured state. It does, however, make a small thermodynamically favorable free energy (∼-13 ± 1 kJ/mol) contribution to the cooperative denatured state collapse transition as temperature is raised and urea concentration is lowered. The thermodynamics of this "cooperative collapse" of the denatured state retain significant compensations between the enthalpy and entropy contributions to the overall free energy.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Insights from molecular dynamics simulations for computational protein design

A grand challenge in the field of structural biology is to design and engineer proteins that exhibit targeted functions. Although much success on this front has been achieved, design success rates remain low, an ever-present reminder of our limited understanding of the relationship between amino acid sequences and the structures they adopt. In addition to experimental techniques and rational design strategies, computational methods have been employed to aid in the design and engineering of proteins. Molecular dynamics (MD) is one such method that simulates the motions of proteins according to classical dynamics. Here, we review how insights into protein dynamics derived from MD simulations have influenced the design of proteins. One of the greatest strengths of MD is its capacity to reveal information beyond what is available in the static structures deposited in the Protein Data Bank. In this regard simulations can be used to directly guide protein design by providing atomistic details of the dynamic molecular interactions contributing to protein stability and function. MD simulations can also be used as a virtual screening tool to rank, select, identify, and assess potential designs. MD is uniquely poised to inform protein design efforts where the application requires realistic models of protein dynamics and atomic level descriptions of the relationship between dynamics and function. Here, we review cases where MD simulations was used to modulate protein stability and protein function by providing information regarding the conformation(s), conformational transitions, interactions, and dynamics that govern stability and function. In addition, we discuss cases where conformations from protein folding/unfolding simulations have been exploited for protein design, yielding novel outcomes that could not be obtained from static structures.

Impact of physiochemical properties on pharmacokinetics of protein therapeutics

Physicochemical properties, such as molecular weight, size, partition coefficient, acid dissociation constant and solubility have a great impact on pharmacokinetics of traditional small molecule drugs and substantially used in development of small drugs. However, predicting pharmacokinetic fate (absorption, distribution, metabolism and elimination) of protein therapeutics from their physicochemical parameters is extremely difficult due to the macromolecular nature of therapeutic proteins and peptides. Their structural complexity and immunogenicity are other contributing factors that determine their biological fate. Therefore, to develop generalized strategies concerning development of therapeutic proteins and peptides are highly challenging. However, reviewing the literature, authors found that physiochemical properties, such as molecular weight, charge and structural modification are having great impact on pharmacokinetics of protein therapeutics and an attempt is made to provide the major findings in this manuscript. This manuscript will serve to provide some bases for developing protein therapeutics with desired pharmacokinetic profile.

Engineering a disulfide bond in the lid hinge region of Rhizopus chinensis lipase: increased thermostability and altered acyl chain length specificity

The key to enzyme function is the maintenance of an appropriate balance between molecular stability and structural flexibility. The lid domain which is very important for “interfacial activation” is the most flexible part in the lipase structure. In this work, rational design was applied to explore the relationship between lid rigidity and lipase activity by introducing a disulfide bond in the hinge region of the lid, in the hope of improving the thermostability of R. chinensis lipase through stabilization of the lid domain without interfering with its catalytic performance. A disulfide bridge between F95C and F214C was introduced into the lipase from R. chinensis in the hinge region of the lid according to the prediction of the “Disulfide by Design” algorithm. The disulfide variant showed substantially improved thermostability with an eleven-fold increase in the t_1/2 value at 60°C and a 7°C increase of T_m compared with the parent enzyme, probably contributed by the stabilization of the geometric structure of the lid region. The additional disulfide bond did not interfere with the catalytic rate (k_cat) and the catalytic efficiency towards the short-chain fatty acid substrate, however, the catalytic efficiency of the disulfide variant towards pNPP decreased by 1.5-fold probably due to the block of the hydrophobic substrate channel by the disulfide bond. Furthermore, in the synthesis of fatty acid methyl esters, the maximum conversion rate by RCLCYS reached 95% which was 9% higher than that by RCL. This is the first report on improving the thermostability of the lipase from R. chinensis by introduction of a disulfide bond in the lid hinge region without compromising the catalytic rate.

Contribution of disulfide bonds to stability, folding, and amyloid fibril formation: the PI3-SH3 domain case

AIMS

The failure of proteins to fold or to remain folded very often leads to their deposition into amyloid fibrils and is the origin of a variety of human diseases. Accordingly, mutations that destabilize the native conformation are associated with pathological phenotypes in several protein models. Protein backbone cyclization by disulfide bond crosslinking strongly reduces the entropy of the unfolded state and, usually, increases protein stability. The effect of protein cyclization on the thermodynamic and kinetics of folding has been extensively studied, but little is know on its effect on aggregation reactions.

RESULTS

The SRC homology 3 domain (SH3) of p85α subunit of bovine phosphatidyl-inositol-3'-kinase (PI3-SH3) domain is a small globular protein, whose folding and amyloid properties are well characterized. Here we describe the effect of polypeptide backbone cyclization on both processes.

INNOVATION

We show that a cyclized PI3-SH3 variant is more stable, folds faster, aggregates slower, and forms conformationally and functionally different amyloid fibrils than the wild-type domain.

CONCLUSION

Disulfide bridges may act as key molecular determinants of both productive protein folding and deleterious aggregation reactions.

All-atom model for stabilization of alpha-helical structure in peptides by hydrocarbon staples

Recent work has shown that the incorporation of an all-hydrocarbon "staple" into peptides can greatly increase their alpha-helix propensity, leading to an improvement in pharmaceutical properties such as proteolytic stability, receptor affinity, and cell permeability. Stapled peptides thus show promise as a new class of drugs capable of accessing intractable targets such as those that engage in intracellular protein-protein interactions. The extent of alpha-helix stabilization provided by stapling has proven to be substantially context dependent, requiring cumbersome screening to identify the optimal site for staple incorporation. In certain cases, a staple encompassing one turn of the helix (attached at residues i and i+4) furnishes greater helix stabilization than one encompassing two turns (i,i+7 staple), which runs counter to expectation based on polymer theory. These findings highlight the need for a more thorough understanding of the forces that underlie helix stabilization by hydrocarbon staples. Here we report all-atom Monte Carlo folding simulations comparing unmodified peptides derived from RNase A and BID BH3 with various i,i+4 and i,i+7 stapled versions thereof. The results of these simulations were found to be in quantitative agreement with experimentally determined helix propensities. We also discovered that staples can stabilize quasi-stable decoy conformations, and that the removal of these states plays a major role in determining the helix stability of stapled peptides. Finally, we critically investigate why our method works, exposing the underlying physical forces that stabilize stapled peptides.

Understanding the Folding Mechanism of an α-Helical Hairpin†

The alpha-helical hairpin is the fundamental building block of the widespread helix-turn-helix DNA binding motif. With two antiparallel helices connected by a reverse turn, the alpha-helical hairpin structure may be regarded as a "supersecondary structural element" and, therefore, could exhibit rather unique folding properties. So far, the folding mechanism of alpha-helical hairpins has not been studied in detail and remains elusive. Herein, we examine the effects of the turn, the hydrophobic cluster, and a disulfide cross-linker on the folding kinetics of a designed alpha-helical hairpin, Z34C, using an infrared temperature-jump (T-jump) method in conjunction with site-specific mutagenesis. Our results show that Z34C folds with an ultrafast rate ( approximately 4.0 x 10(5) s(-1)) and support a folding mechanism in which the rate-limiting step corresponds to the formation of the reverse turn. On the other hand, the hydrophobic cluster and the disulfide cross-linker appear to largely stabilize the native state but not the folding transition state.

Role of the disulphide bridge in folding, stability and function of porcine odorant binding protein: spectroscopic equilibrium studies on C63A/C155A double mutant

Porcine odorant binding protein (pOBP) contains a single disulphide bridge linking residues Cys63 and Cys155. In order to get information on the role played by this crosslink in determining the structural and functional properties of the protein, we substituted these two Cys residues with two Ala residues by site directed mutagenesis and investigated the changes in folding, stability and functional features, as detected by fluorescence and circular dichroism measurements. In particular, we studied both chemical and thermal unfolding/refolding processes under equilibrium conditions, the first induced by guanidinium hydrochloride and the second by raising the temperature from 15 to 90 degrees C. Chemical unfolding curves, as obtained from intrinsic fluorescence and far-UV circular dichroism data, can be fitted by a simple two-state cooperative sigmoidal function; however, their partial overlap (C(1/2)=0.57+/-0.05 from fluorescence and 0.66+/-0.03 from CD) suggests the formation of an intermediate, which lacks tertiary structural features. Thermal unfolding was found to be reversible if the protein was heated up to 65 degrees C, but irreversible above that temperature because of aggregation. The thermodynamic unfolding parameters of this double mutant protein, when compared to those of the wild type protein, clearly point out the important role played by the disulphide bridge on the stability and function of this protein family and probably of many other lipocalins.

Effects of the difference in the unfolded-state ensemble on the folding of Escherichia coli dihydrofolate reductase

The unfolded state of a protein is an ensemble of a large number of conformations ranging from fully extended to compact structures. To investigate the effects of the difference in the unfolded-state ensemble on protein folding, we have studied the structure, stability, and folding of "circular" dihydrofolate reductase (DHFR) from Escherichia coli in which the N and C-terminal regions are cross-linked by a disulfide bond, and compared the results with those of disulfide-reduced "linear" DHFR. Equilibrium studies by circular dichroism, difference absorption spectra, solution X-ray scattering, and size-exclusion chromatography show that whereas the native structures of both proteins are essentially the same, the unfolded state of circular DHFR adopts more compact conformations than the unfolded state of the linear form, even with the absence of secondary structure. Circular DHFR is more stable than linear DHFR, which may be due to the decrease in the conformational entropy of the unfolded state as a result of circularization. Kinetic refolding measurements by stopped-flow circular dichroism and fluorescence show that under the native conditions both proteins accumulate a burst-phase intermediate having the same structures and both fold by the same complex folding mechanism with the same folding rates. Thus, the effects of the difference in the unfolded state of circular and linear DHFRs on the refolding reaction are not observed after the formation of the intermediate. This suggests that for the proteins with close termini in the native structure, early compaction of a protein molecule to form a specific folding intermediate with the N and C-terminal regions in close proximity is a crucial event in folding. If there is an enhancement in the folding reflecting the reduction in the breadth of the unfolded-state ensemble for circular DHFR, this acceleration must occur in the sub-millisecond time-range.

Dual effects of an extra disulfide bond on the activity and stability of a cold-adapted alpha-amylase

Chloride-dependent alpha-amylases constitute a well conserved family of enzymes thereby allowing investigation of the characteristics of each member to understand, for example, relevant properties required for environmental adaptation. In this context, we have constructed a double mutant (Q58C/A99C) of the cold-active and heat-labile alpha-amylase from the Antarctic bacterium Pseudoalteromonas haloplanktis, defined on the basis of its strong similarity with the mesophilic enzyme from pig pancreas. This mutant was characterized to understand the role of an extra disulfide bond specific to warm-blooded animals and located near the entrance of the catalytic cleft. We show that the catalytic parameters of the mutant are drastically modified and similar to those of the mesophilic enzyme. Calorimetric studies demonstrated that the mutant is globally stabilized (DeltaDeltaG = 1.87 kcal/mol at 20 degrees C) when compared with the wild-type enzyme, although the melting point (T(m)) was not increased. Moreover, fluorescence quenching experiments indicate a more compact structure for the mutated alpha-amylase. However, the strain imposed on the active site architecture induces a 2-fold higher thermal inactivation rate at 45 degrees C as well as the appearance of a less stable calorimetric domain. It is concluded that stabilization by the extra disulfide bond arises from an enthalpy-entropy compensation effect favoring the enthalpic contribution.

The cystine knot promotes folding and not thermodynamic stability in vascular endothelial growth factor

Cystine knots consist of three intertwined disulfide bridges and are considered major determinants of protein stability in proteins in which they occur. We questioned this function and observed that removal of individual disulfide bridges in human vascular endothelial growth factor (VEGF) does not reduce its thermodynamic stability but reduces its unexpected high thermal stability of 108 degrees C by up to 40 degrees C. In wild-type VEGF (deltaG(u,25)(0) = 5.1 kcal.mol(-1)), the knot is responsible for a large entropic stabilization of TdeltaS(u,25)(0) = -39.3 kcal mol(-1), which is compensated for by a deltaH(u,25)(0) of -34.2 kcal mol(-1). In the disulfide-deficient mutants, this entropic stabilization disappears, but instead of a decrease, we observe an increase in the thermodynamic stability by about 2 kcal.mol(-1). A detailed crystallographic analysis of the mutant structures suggests a role of the cystine knot motif in protein folding rather than in the stabilization of the folded state. When assuming that the sequential order of the disulfide bridge formation is conserved between VEGF and glycoprotein alpha-subunit, the crystal structure of the mutant C61A-C104A, which deviates by a root mean square deviation of more than 2.2 A from wild-type VEGF, identifies a true folding intermediate of VEGF.

Heat does not come in different colours: entropy-enthalpy compensation, free energy windows, quantum confinement, pressure perturbation calorimetry, solvation and the multiple causes of heat capacity effects in biomolecular interactions

Modern techniques in microcalorimetry allow us to measure directly the heat changes and associated thermodynamics for biomolecular processes in aqueous solution at reasonable concentrations. All these processes involve changes in solvation/hydration, and it is natural to assume that the heats for these processes should reflect, in some way, such changes in solvation. However, the interpretation of data is still somewhat ambiguous, since different non-covalent interactions may have similar thermodynamic signatures, and analysis is frustrated by large entropy-enthalpy compensation effects. Changes in heat capacity (Delta C(p)) have been related to changes in hydrophobic hydration and non-polar accessible surface areas, but more recent empirical and theoretical work has shown how this need not always be the case. Entropy-enthalpy compensation is a natural consequence of finite Delta C(p) values and, more generally, can arise as a result of quantum confinement effects, multiple weak interactions, and limited free energy windows, giving rise to thermodynamic homeostasis that may be of evolutionary and functional advantage. The new technique of pressure perturbation calorimetry (PPC) has enormous potential here as a means of probing solvation-related volumetric changes in biomolecules at modest pressures, as illustrated with preliminary data for a simple protein-inhibitor complex.

Disulfide bond effects on protein stability: designed variants of Cucurbita maxima trypsin inhibitor-V

Attempts to increase protein stability by insertion of novel disulfide bonds have not always been successful. According to the two current models, cross-links enhance stability mainly through denatured state effects. We have investigated the effects of removal and addition of disulfide cross-links, protein flexibility in the vicinity of a cross-link, and disulfide loop size on the stability of Cucurbita maxima trypsin inhibitor-V (CMTI-V; 7 kD) by differential scanning calorimetry. CMTI-V offers the advantage of a large, flexible, and solvent-exposed loop not involved in extensive intra-molecular interactions. We have uncovered a negative correlation between retention time in hydrophobic column chromatography, a measure of protein hydrophobicity, and melting temperature (T(m)), an indicator of native state stabilization, for CMTI-V and its variants. In conjunction with the complete set of thermodynamic parameters of denaturation, this has led to the following deductions: (1) In the less stable, disulfide-removed C3S/C48S (Delta Delta G(d)(50 degrees C) = -4 kcal/mole; Delta T(m) = -22 degrees C), the native state is destabilized more than the denatured state; this also applies to the less-stable CMTI-V* (Delta Delta G(d)(50 degrees C) = -3 kcal/mole; Delta T(m) = -11 degrees C), in which the disulfide-containing loop is opened by specific hydrolysis of the Lys(44)-Asp(45) peptide bond; (2) In the less stable, disulfide-inserted E38C/W54C (Delta Delta G(d)(50 degrees C) = -1 kcal/mole; Delta T(m) = +2 degrees C), the denatured state is more stabilized than the native state; and (3) In the more stable, disulfide-engineered V42C/R52C (Delta Delta G(d)(50 degrees C) = +1 kcal/mole; Delta T(m) = +17 degrees C), the native state is more stabilized than the denatured state. These results show that a cross-link stabilizes both native and denatured states, and differential stabilization of the two states causes either loss or gain in protein stability. Removal of hydrogen bonds in the same flexible region of CMTI-V resulted in less destabilization despite larger changes in the enthalpy and entropy of denaturation. The effect of a cross-link on the denatured state of CMTI-V was estimated directly by means of a four-state thermodynamic cycle consisting of native and denatured states of CMTI-V and CMTI-V*. Overall, the results show that an enthalpy-entropy compensation accompanies disulfide bond effects and protein stabilization is profoundly modulated by altered hydrophobicity of both native and denatured states, altered flexibility near the cross-link, and residual structure in the denatured state.

The effects of disulfide bonds on the denatured state of barnase

The effects of engineered disulfide bonds on protein stability are poorly understood because they can influence the structure, dynamics, and energetics of both the native and denatured states. To explore the effects of two engineered disulfide bonds on the stability of barnase, we have conducted a combined molecular dynamics and NMR study of the denatured state of the two mutants. As expected, the disulfide bonds constrain the denatured state. However, specific extended beta-sheet structure can also be detected in one of the mutant proteins. This mutant is also more stable than would be predicted. Our study suggests a possible cause of the very high stability conferred by this disulfide bond: the wild-type denatured ensemble is stabilized by a nonnative hydrophobic cluster, which is constrained from occurring in the mutant due to the formation of secondary structure.

A study on the enthalpy-entropy compensation in protein unfolding

A large number of thermodynamic data including the free energy, enthalpy, entropy, and heat capacity changes were collected for the denaturation of various proteins. Regression indicated that remarkable enthalpy-entropy compensation occurred in protein unfolding, which meant that the change in enthalpy was almost compensated by a corresponding change in entropy resulting in a smaller net free energy change. This behavior was proposed to result from the water molecule reorganization, which contributed significantly to the enthalpy and entropy changes but little to the free energy change in protein unfolding. It turned out that the enthalpy-entropy compensation could provide novel insights into the problem of enthalpy and entropy convergence in protein unfolding.

Measuring denatured state energetics: deviations from random coil behavior and implications for the folding of iso-1-cytochrome c

The changes in the free energy of the denatured state of a set of yeast iso-1-cytochrome c variants with single surface histidine residues have been measured in 3 M guanidine hydrochloride. The thermodynamics of unfolding by guanidine hydrochloride is also reported. All variants have decreased stability relative to the wild-type protein. The free energy of the denatured state was determined in 3 M guanidine hydrochloride by evaluating the strength of heme-histidine ligation through determination of the pK(a) for loss of histidine binding to the heme. The data are corrected for the presence of the N-terminal amino group which also ligates to the heme under similar solution conditions. Significant deviations from random coil behavior are observed. Relative to a variant with a single histidine at position 26, residual structure of the order of -1.0 to -2.5 kcal/mol is seen for the other variants studied. The data explain the slower folding of yeast iso-1-cytochrome c relative to the horse protein. The greater number of histidines and the greater strength of ligation are expected to slow conversion of the histidine-misligated forms to the obligatory aquo-heme intermediate during the ligand exchange phase of folding. The particularly strong association of histidine residues at positions 54 and 89 may indicate regions of the protein with strong energetic propensities to collapse against the heme during early folding events, consistent with available data in the literature on early folding events for cytochrome c.

Instability, stabilization, and formulation of liquid protein pharmaceuticals

One of the most challenging tasks in the development of protein pharmaceuticals is to deal with physical and chemical instabilities of proteins. Protein instability is one of the major reasons why protein pharmaceuticals are administered traditionally through injection rather than taken orally like most small chemical drugs. Protein pharmaceuticals usually have to be stored under cold conditions or freeze-dried to achieve an acceptable shelf life. To understand and maximize the stability of protein pharmaceuticals or any other usable proteins such as catalytic enzymes, many studies have been conducted, especially in the past two decades. These studies have covered many areas such as protein folding and unfolding/denaturation, mechanisms of chemical and physical instabilities of proteins, and various means of stabilizing proteins in aqueous or solid state and under various processing conditions such as freeze-thawing and drying. This article reviews these investigations and achievements in recent years and discusses the basic behavior of proteins, their instabilities, and stabilization in aqueous state in relation to the development of liquid protein pharmaceuticals.

Structural response to mutation at a protein-protein interface

We have crystallised three mutants of the barnase-barstar complex in which interactions across the interface have been deleted by simultaneous mutation of both residues involved in the interaction. Each mutant deletes a different type of interaction at the interface: the first complex bnHis102-->Ala-bsTyr29-->Phe (bn, barnase; bs, barstar), deletes a van der Waals packing interaction; the second complex, bnLys27-->Ala-bsThr42-->Ala, deletes a hydrogen bond; the third, bnLys27-->Ala-bsAsp35-->Ala, deletes a long-range charge-charge interaction. The contribution of each of these side-chains to the stability of the complex is known; the coupling energy between the deleted side-chains is also known. Despite each of the double mutants being significantly destabilised compared with the wild-type, the effects of mutation are local. Only small movements in the main-chain surrounding the sites of mutation and some larger movements of neighbouring side-chains are observed in the mutant complexes. The exact response to mutation is context-dependent and for the same mutant can vary depending upon the environment within the crystal. In some double mutant complexes, interfacial pockets, which are accessible to bulk solvent are formed, whereas interfacial cavities which are isolated from bulk solvent, are formed in others. In all double mutants, water molecules fill the created pockets and cavities. These water molecules mimic the deleted side-chains by occupying positions close to the non-carbon atoms of truncated side-chains and re-making many hydrogen bonds made by the truncated side-chains in the wild-type. It remains extremely difficult, however, to correlate energetic and structural responses to mutation because of unknown changes in entropy and entropy-enthalpy compensation.

Contribution of a conserved asparagine to the conformational stability of ribonucleases Sa, Ba, and T1

The contribution of hydrogen bonding by peptide groups to the conformational stability of globular proteins was studied. One of the conserved residues in the microbial ribonuclease (RNase) family is an asparagine at position 39 in RNase Sa, 44 in RNase T1, and 58 in RNase Ba (barnase). The amide group of this asparagine is buried and forms two similar intramolecular hydrogen bonds with a neighboring peptide group to anchor a loop on the surface of all three proteins. Thus, it is a good model for the hydrogen bonding of peptide groups. When the conserved asparagine is replaced with alanine, the decrease in the stability of the mutant proteins is 2.2 (Sa), 1.8 (T1), and 2.7 (Ba) kcal/mol. When the conserved asparagine is replaced by aspartate, the stability of the mutant proteins decreases by 1.5 and 1.8 kcal/mol for RNases Sa and T1, respectively, but increases by 0.5 kcal/mol for RNase Ba. When the conserved asparagine was replaced by serine, the stability of the mutant proteins was decreased by 2.3 and 1.7 kcal/mol for RNases Sa and T1, respectively. The structure of the Asn 39 --> Ser mutant of RNase Sa was determined at 1.7 A resolution. There is a significant conformational change near the site of the mutation: (1) the side chain of Ser 39 is oriented differently than that of Asn 39 and forms hydrogen bonds with two conserved water molecules; (2) the peptide bond of Ser 42 changes conformation in the mutant so that the side chain forms three new intramolecular hydrogen bonds with the backbone to replace three hydrogen bonds to water molecules present in the wild-type structure; and (3) the loss of the anchoring hydrogen bonds makes the surface loop more flexible in the mutant than it is in wild-type RNase Sa. The results show that burial and hydrogen bonding of the conserved asparagine make a large contribution to microbial RNase stability and emphasize the importance of structural information in interpreting stability studies of mutant proteins.

Conformational stability and thermodynamics of folding of ribonucleases Sa, Sa2 and Sa3

Ribonucleases Sa, Sa2, and Sa3 are three small, extracellular enzymes produced by different strains of Streptomyces aureofaciens with amino acid sequences that are 50% identical. We have studied the unfolding of these enzymes by heat and urea to determine the conformational stability and its dependence on temperature, pH, NaCl, and the disulfide bond. All three of the Sa ribonucleases unfold reversibly by a two-state mechanism with melting temperatures, Tm, at pH 7 of 48.4 degrees C (Sa), 41.1 degrees C (Sa2), and 47.2 degrees C (Sa3). The Tm values are increased in the presence of 0.5 M NaCl by 4.0 deg. C (Sa), 0.1 deg. C (Sa2), and 7.2 deg. C (Sa3). The Tm values are decreased by 20.0 deg. C (Sa), 31.5 deg. C (Sa2), and 27.0 deg. C (Sa3) when the single disulfide bond in the molecules is reduced. We compare these results with similar studies on two other members of the microbial ribonuclease family, RNase T1 and RNase Ba (barnase), and with a member of the mammalian ribonuclease family, RNase A. At pH 7 and 25 degrees C, the conformational stabilities of the ribonucleases are (kcal/mol): 2.9 (Sa2), 5.6 (Sa3), 6.1 (Sa), 6.6 (T1), 8.7 (Ba), and 9.2 (A). Our analysis of the stabilizing forces suggests that the hydrophobic effect contributes from 90 to 110 kcal/mol and that hydrogen bonding contributes from 70 to 105 kcal/mol to the stability of these ribonucleases. Thus, we think that the hydrophobic effect and hydrogen bonding make large but comparable contributions to the conformational stability of these proteins.

Nine hydrophobic side chains are key determinants of the thermodynamic stability and oligomerization status of tumour suppressor p53 tetramerization domain

The contribution of almost each amino acid side chain to the thermodynamic stability of the tetramerization domain (residues 326-353) of human p53 has been quantitated using 25 mutants with single-residue truncations to alanine (or glycine). Truncation of either Leu344 or Leu348 buried at the tetramer interface, but not of any other residue, led to the formation of dimers of moderate stability (8-9 kcal/mol of dimer) instead of tetramers. One-third of the substitutions were moderately destabilizing (<3.9 kcal/mol of tetramer). Truncations of Arg333, Asn345 or Glu349 involved in intermonomer hydrogen bonds, Ala347 at the tetramer interface or Thr329 were more destabilizing (4.1-5.7 kcal/mol). Strongly destabilizing (8.8- 11.7 kcal/mol) substitutions included those of Met340 at the tetramer interface and Phe328, Arg337 and Phe338 involved peripherally in the hydrophobic core. Truncation of any of the three residues involved centrally in the hydrophobic core of each primary dimer either prevented folding (Ile332) or allowed folding only at high protein concentration or low temperature (Leu330 and Phe341). Nine hydrophobic residues per monomer constitute critical determinants for the stability and oligomerization status of this p53 domain.

Glutamine, alanine or glycine repeats inserted into the loop of a protein have minimal effects on stability and folding rates

Natural proteins can contain flexible regions in their polypeptide chain. We have investigated the effects of glycine, alanine and glutamine repeats on the stability and folding of a protein by inserting stretches of 7 to 13 residues into a suitable position in a model system, the chymotrypsin inhibitor-2 (CI2). This folds by residues (1-40) docking with residues (41-64) to form a folding nucleus. The peptides GQ4GM, GQ6GM, GQ8GM, GQ10GM, GA2SA4SA2GM and G3SG4SG3M were inserted after residue 40. The stability of the mutant proteins changes only weakly with chain length and nature of insertion, suggesting that the presence of unstructured polypeptide chains in a protein does not have a great energetic penalty. This has implications in catalysis, for example, where floppy regions have been noted in active sites, and in DNA transcription where activators, transcription factors and intermediary proteins all show long repeats of glycine/serine and/or glutamine, which are thought to be important for function. We find that the rate of folding is very insensitive to the length of the linker. The changes in rate are close to those predicted from polymer theory for the loss of configuration entropy on closing a loop. This implies that all the diffusion steps are relatively rapid.