Contribution of individual side-chains to the stability of BPTI examined by alanine-scanning mutagenesis

Free-Energy Landscape and Rate Estimation of the Aromatic Ring Flips in Basic Pancreatic Trypsin Inhibitors Using Metadynamics

Aromatic side chains (phenylalanine and tyrosine) of a protein flip by 180° around the Cβ-Cγ axis (χ2 dihedral of the side chain), producing two symmetry-equivalent states. The study of ring flip dynamics with nuclear magnetic resonance (NMR) experiments helps to understand local conformational fluctuations. Ring flips are categorized as slow (milliseconds and onward) or fast (nanoseconds to near milliseconds) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to estimate the flip rate and discriminate between slow and fast ring flips for eight individual aromatic side chains (F4, Y10, Y21, F22, Y23, F33, Y35, and F45) of the basic pancreatic trypsin inhibitor. Well-tempered metadynamics simulations were performed to estimate the ring-flipping free-energy surfaces for all eight aromatic residues. The results indicate that χ2 as a standalone collective variable (CV) is not sufficient to obtain computationally consistent results. Inclusion of a complementary CV, such as χ1(Cα-Cβ), solved the problem for most residues and enabled us to classify fast and slow ring flips. This indicates the importance of librational motions in ring flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events were observed for residues F22 and F33, indicating a possible role of friction effects in ring flipping. The results demonstrate the successful application of infrequent metadynamics to estimate ring flip rates and identify certain limitations of the approach.

Probing protein stability: towards a computational atomistic, reliable, affordable, and improvable model

We present an improved application of a recently proposed computational method designed to evaluate the change of free energy as a function of the average value of a suitably chosen collective variable in proteins. The method is based on a full atomistic description of the protein and its environment. The goal is to understand how the protein melting temperature changes upon single-point mutations, because the sign of the temperature variation will allow us to discriminate stabilizing vs. destabilizing mutations in protein sequences. In this refined application the method is based on altruistic well-tempered metadynamics, a variant of multiple-walkers metadynamics. The resulting metastatistics is then modulated by the maximal constrained entropy principle. The latter turns out to be especially helpful in free-energy calculations as it is able to alleviate the severe limitations of metadynamics in properly sampling folded and unfolded configurations. In this work we apply the computational strategy outlined above in the case of the bovine pancreatic trypsin inhibitor, a well-studied small protein, which is a reference for computer simulations since decades. We compute the variation of the melting temperature characterizing the folding-unfolding process between the wild-type protein and two of its single-point mutations that are seen to have opposite effect on the free energy changes. The same approach is used for free energy difference calculations between a truncated form of frataxin and a set of five of its variants. Simulation data are compared to in vitro experiments. In all cases the sign of the change of melting temperature is reproduced, under the further approximation of using an empirical effective mean-field to average out protein-solvent interactions.

Molecular Information Theory Meets Protein Folding

We propose an application of molecular information theory to analyze the folding of single domain proteins. We analyze results from various areas of protein science, such as sequence-based potentials, reduced amino acid alphabets, backbone configurational entropy, secondary structure content, residue burial layers, and mutational studies of protein stability changes. We found that the average information contained in the sequences of evolved proteins is very close to the average information needed to specify a fold ∼2.2 ± 0.3 bits/(site·operation). The effective alphabet size in evolved proteins equals the effective number of conformations of a residue in the compact unfolded state at around 5. We calculated an energy-to-information conversion efficiency upon folding of around 50%, lower than the theoretical limit of 70%, but much higher than human-built macroscopic machines. We propose a simple mapping between molecular information theory and energy landscape theory and explore the connections between sequence evolution, configurational entropy, and the energetics of protein folding.

Large-scale design and refinement of stable proteins using sequence-only models

Engineered proteins generally must possess a stable structure in order to achieve their designed function. Stable designs, however, are astronomically rare within the space of all possible amino acid sequences. As a consequence, many designs must be tested computationally and experimentally in order to find stable ones, which is expensive in terms of time and resources. Here we use a high-throughput, low-fidelity assay to experimentally evaluate the stability of approximately 200,000 novel proteins. These include a wide range of sequence perturbations, providing a baseline for future work in the field. We build a neural network model that predicts protein stability given only sequences of amino acids, and compare its performance to the assayed values. We also report another network model that is able to generate the amino acid sequences of novel stable proteins given requested secondary sequences. Finally, we show that the predictive model-despite weaknesses including a noisy data set-can be used to substantially increase the stability of both expert-designed and model-generated proteins.

Simulation analysis of selective alanine mutation effect on stability of human prion protein

Prion diseases are neurodegenerative disorders caused by spongiform degeneration of the brain. Understanding the fundamental mechanism of prion protein aggregation caused by mutations is very crucial to resolve the pathology of prion diseases. To help understand the roles of individual residues on the stability of the human prion protein, the computational method of free energy simulations based on atomistic molecular dynamics trajectories is applied to Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations of the human prion protein. The simulations show that all three alanine mutations destabilize the human prion protein. The calculated free energy change differences, ΔΔG, for the Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations are in good agreement with the experimental values. The significant destabilizing effects on the mutants relative to the wild-type protein arise from van der Waals terms. Furthermore, our free energy decomposition analysis shows that the major contribution to destabilizing the V180A and V209A mutants relative to the wild-type protein is originated from van der Waals interactions from residues near the mutation sites. In contrast, the contribution to destabilizing the F175A mutant is mainly caused by van der Waals interactions from residues near and far away from the mutation site. Our results show that the free energy simulation with a thermodynamic integration approach for selected alanine scanning mutations is beneficial for understanding the detailed mechanism of human prion protein destabilization, specific residues' role, and the hydrophobic effect on protein stability.Communicated by Ramaswamy H. Sarma.

Membrane Protein Stabilization Strategies for Structural and Functional Studies

Accounting for nearly two-thirds of known druggable targets, membrane proteins are highly relevant for cell physiology and pharmacology. In this regard, the structural determination of pharmacologically relevant targets would facilitate the intelligent design of new drugs. The structural biology of membrane proteins is a field experiencing significant growth as a result of the development of new strategies for structure determination. However, membrane protein preparation for structural studies continues to be a limiting step in many cases due to the inherent instability of these molecules in non-native membrane environments. This review describes the approaches that have been developed to improve membrane protein stability. Membrane protein mutagenesis, detergent selection, lipid membrane mimics, antibodies, and ligands are described in this review as approaches to facilitate the production of purified and stable membrane proteins of interest for structural and functional studies.

Deep2Full: Evaluating strategies for selecting the minimal mutational experiments for optimal computational predictions of deep mutational scan outcomes

Performing a complete deep mutational scan with all single point mutations may not be practical, and may not even be required, especially if predictive computational models can be developed. Computational models are however naive to cellular response in the myriads of assay-conditions. In a realistic paradigm of assay context-aware predictive hybrid models that combine minimal experimental data from deep mutational scans with structure, sequence information and computational models, we define and evaluate different strategies for choosing this minimal set. We evaluated the trivial strategy of a systematic reduction in the number of mutational studies from 85% to 15%, along with several others about the choice of the types of mutations such as random versus site-directed with the same 15% data completeness. Interestingly, the predictive capabilities by training on a random set of mutations and using a systematic substitution of all amino acids to alanine, asparagine and histidine (ANH) were comparable. Another strategy we explored, augmenting the training data with measurements of the same mutants at multiple assay conditions, did not improve the prediction quality. For the six proteins we analyzed, the bin-wise error in prediction is optimal when 50-100 mutations per bin are used in training the computational model, suggesting that good prediction quality may be achieved with a library of 500-1000 mutations.

Crystal structures of highly simplified BPTIs provide insights into hydration-driven increase of unfolding enthalpy

We report a thermodynamic and structural analysis of six extensively simplified bovine pancreatic trypsin inhibitor (BPTI) variants containing 19–24 alanines out of 58 residues. Differential scanning calorimetry indicated a two-state thermal unfolding, typical of a native protein with densely packed interior. Surprisingly, increasing the number of alanines induced enthalpy stabilization, which was however over-compensated by entropy destabilization. X-ray crystallography indicated that the alanine substitutions caused the recruitment of novel water molecules facilitating the formation of protein–water hydrogen bonds and improving the hydration shells around the alanine’s methyl groups, both of which presumably contributed to enthalpy stabilization. There was a strong correlation between the number of water molecules and the thermodynamic parameters. Overall, our results demonstrate that, in contrast to our initial expectation, a protein sequence in which over 40% of the residues are alanines can retain a densely packed structure and undergo thermal denaturation with a large enthalpy change, mainly contributed by hydration.

Evaluation of pharmacokinetics and the stability of daptomycin in serum at various temperatures

BACKGROUND

Daptomycin exhibits concentration-dependent antibacterial activity. By monitoring daptomycin serum concentrations, clinicians may be able to predict the effectiveness of treatments for infections more accurately. However, it has been reported that daptomycin concentrations in plasma samples stored at -20°C decrease approximately 25% after 4 weeks. The aim of this study was to evaluate the stability of daptomycin in serum at various temperatures.

METHODS

Daptomycin serum samples were prepared and stored at different temperatures. The stability of daptomycin under various conditions was evaluated by sequential measurements of concentration.

RESULTS

Although the loss of concentration of daptomycin in serum samples stored in freezers (-80°C and -20°C) was less than 10% after 168days (6 months), the concentrations in samples stored in a refrigerator (4°C) decreased by more than 70% over the same period. Furthermore, daptomycin concentrations in serum samples stored at close to body temperature (35°C, 37°C, and 39°C) decreased by more than 50% after only 24h.

CONCLUSIONS

The results of the present study demonstrate that the measurement of serum concentrations of daptomycin needs to be performed rapidly. Furthermore, the degradation of daptomycin in serum may be involved in its elimination from the living body.

The Structure-Forming Juncture in Oxidative Protein Folding: What Happens in the ER?

The folding of disulfide bond containing proteins proceeds in a biphasic manner. Initially, cysteines are oxidized to form disulfide bonds. Structure is largely absent during this phase. Next, when a minimally correct number of native linkages of disulfide bonds have been acquired, the biopolymer conformationally folds into the native, or a native-like, state. Thus, at the end of this "oxidative folding" process, a stable and biologically active protein is formed. This review focuses on dissecting the "structure-forming step" in oxidative protein folding. The ability to follow this pivotal step in protein maturation in somewhat detail is uniquely facilitated in "oxidative" folding scenarios. We review this step using bovine pancreatic Ribonuclease A as a model while recognizing the impact that this step has in subcellular trafficking and protein aggregation.

Reversible thermal unfolding of a yfdX protein with chaperone-like activity

yfdX proteins are ubiquitously present in a large number of virulent bacteria. A member of this family of protein in E. coli is known to be up-regulated by the multidrug response regulator. Their abundance in such bacteria suggests some important yet unidentified functional role of this protein. Here, we study the thermal response and stability of yfdX protein STY3178 from Salmonella Typhi using circular dichroism, steady state fluorescence, dynamic light scattering and nuclear magnetic resonance experiments. We observe the protein to be stable up to a temperature of 45 °C. It folds back to the native conformation from unfolded state at temperature as high as 80 °C. The kinetic measurements of unfolding and refolding show Arrhenius behavior where the refolding involves less activation energy barrier than that of unfolding. We propose a homology model to understand the stability of the protein. Our molecular dynamic simulation studies on this model structure at high temperature show that the structure of this protein is quite stable. Finally, we report a possible functional role of this protein as a chaperone, capable of preventing DTT induced aggregation of insulin. Our studies will have broader implication in understanding the role of yfdX proteins in bacterial function and virulence.

A many-body term improves the accuracy of effective potentials based on protein coevolutionary data

The study of correlated mutations in alignments of homologous proteins proved to be successful not only in the prediction of their native conformation but also in the development of a two-body effective potential between pairs of amino acids. In the present work, we extend the effective potential, introducing a many-body term based on the same theoretical framework, making use of a principle of maximum entropy. The extended potential performs better than the two-body one in predicting the energetic effect of 308 mutations in 14 proteins (including membrane proteins). The average value of the parameters of the many-body term correlates with the degree of hydrophobicity of the corresponding residues, suggesting that this term partly reflects the effect of the solvent.

Single, double and quadruple alanine substitutions at oligomeric interfaces identify hydrophobicity as the key determinant of human neutrophil alpha defensin HNP1 function

HNP1 is a human alpha defensin that forms dimers and multimers governed by hydrophobic residues, including Tyr¹⁶, Ile²⁰, Leu²⁵, and Phe²⁸. Previously, alanine scanning mutagenesis identified each of these residues and other hydrophobic residues as important for function. Here we report further structural and functional studies of residues shown to interact with one another across oligomeric interfaces: I20A-HNP1 and L25A-HNP1, plus the double alanine mutants I20A/L25A-HNP1 and Y16A/F28A-HNP1, and the quadruple alanine mutant Y16A/I20A/L25A/F28A-HNP1. We tested binding to HIV-1 gp120 and HNP1 by surface plasmon resonance, binding to HIV-1 gp41 and HNP1 by fluorescence polarization, inhibition of anthrax lethal factor, and antibacterial activity using the virtual colony count assay. Similar to the previously described single mutant W26A-HNP1, the quadruple mutant displayed the least activity in all functional assays, followed by the double mutant Y16A/F28A-HNP1. The effects of the L25A and I20A single mutations were milder than the double mutant I20A/L25A-HNP1. Crystallographic studies confirmed the correct folding and disulfide pairing, and depicted an array of dimeric and tetrameric structures. These results indicate that side chain hydrophobicity is the critical factor that determines activity at these positions.

Analysis of amino acid contributions to protein solubility using short peptide tags fused to a simplified BPTI variant

Protein solubility is usually characterized in terms of a hydrophobicity scale, which refers to the free energy of transfer of a molecule from an aqueous to a nonpolar solution and is not a "solubility propensity scale" per se. Using a "host-guest" approach, we measured the effects of short poly-amino-acid tags (guests) on the solubility of a host protein, a simplified bovine pancreatic trypsin inhibitor (BPTI), to which they were fused at the C-terminus. We analyzed 10 amino acid types, representing the full range of biophysical properties (acidic, basic, polar, and hydrophobic). As anticipated, positively charged residues significantly increased the solubility of the model protein, at both pH 4.7 and 7.7, whereas very hydrophobic poly-Ile markedly reduced the solubility of BPTI. Poly-Asp and poly-Glu barely affected BPTI solubility at pH 4.7, but induced an eight to ten-fold increase at pH 7.7, attributable to the ionization of their side chains. Although Pro is the most soluble amino acid, poly-Pro did not affect the protein's solubility. The effects of the other tags on BPTI solubility ranged from none to an eight-fold increase. To ensure that the measured solubility values were context independent and could provide a "solubility propensity scale", we confirmed that the tags remained independent of the structure, thermal stability, and biochemical activity of the host protein. These observations suggest that this approach is valuable for measuring the solubility propensities of amino acids, which could eventually allow the calculation of a polypeptide's relative solubility from its amino acid sequence.

Oxidative folding and structural analyses of a Kunitz-related inhibitor and its disulfide intermediates: functional implications

Tick-derived protease inhibitor (TdPI) is a tight-binding Kunitz-related inhibitor of human tryptase β with a unique structure and disulfide-bond pattern. Here we analyzed its oxidative folding and reductive unfolding by chromatographic and disulfide analyses of acid-trapped intermediates. TdPI folds through a stepwise generation of heterogeneous populations of one-disulfide, two-disulfide, and three-disulfide intermediates, with a major accumulation of the nonnative three-disulfide species IIIa. The rate-limiting step of the process is disulfide reshuffling within the three-disulfide population towards a productive intermediate that oxidizes directly into the native four-disulfide protein. TdPI unfolds through a major accumulation of the native three-disulfide species IIIb and the subsequent formation of two-disulfide and one-disulfide intermediates. NMR characterization of the acid-trapped and further isolated IIIa intermediate revealed a highly disordered conformation that is maintained by the presence of the disulfide bonds. Conversely, the NMR structure of IIIb showed a native-like conformation, with three native disulfide bonds and increased flexibility only around the two free cysteines, thus providing a molecular basis for its role as a productive intermediate. Comparison of TdPI with a shortened variant lacking the flexible prehead and posthead segments revealed that these regions do not contribute to the protein conformational stability or the inhibition of trypsin but are important for both the initial steps of the folding reaction and the inhibition of tryptase β. Taken together, the results provide insights into the mechanism of oxidative folding of Kunitz inhibitors and pave the way for the design of TdPI variants with improved properties for biomedical applications.

Sequence determinants of thermodynamic stability in a WW domain--an all-beta-sheet protein

The stabilities of 66 sequence variants of the human Pin1 WW domain have been determined by equilibrium thermal denaturation experiments. All 34 residues composing the hPin1 WW three-stranded beta-sheet structure could be replaced one at a time with at least one different natural or non-natural amino acid residue without leading to an unfolded protein. Alanine substitutions at only four positions within the hPin1 WW domain lead to a partially or completely unfolded protein-in the absence of a physiological ligand. The side chains of these four residues form a conserved, partially solvent-inaccessible, continuous hydrophobic minicore comprising the N- and C-termini. Ala mutations at five other residues, three of which constitute the ligand binding patch on the concave side of the beta-sheet, significantly destabilize the hPin1 WW domain without leading to an unfolded protein. The remaining mutations affect protein stability only slightly, suggesting that only a small subset of side chain interactions within the hPin1 WW domain are mandatory for acquiring and maintaining a stable, cooperatively folded beta-sheet structure.

Brownian dynamics simulations of a wild type and mutants of bovine pancreatic trypsin inhibitors

Brownian dynamics simulations of a wild type and mutants of bovine pancreatic trypsin inhibitors were performed using the program brownian. The results of the simulations were consistent with experimentally determined stabilities of these proteins, and the computational times of the simulations were much less than those of molecular dynamics simulations. These results indicate that Brownian dynamics simulations can be useful for investigating structural features of proteins that have some mutants, such as drug-metabolizing enzymes.

Stable-isotope labeling using an inducible viral infection system in suspension-cultured plant cells

We established a novel strategy for preparing uniformly stable isotope-labeled proteins by using suspension-cultured plant cells and an inducible virus vector encoding the research target. By using this new method, we demonstrated the expression of three proteins, namely, Escherichia coli dihydrofolate reductase (DHFR), chicken calmodulin (CaM), and porcine protein kinase C-dependent protein phosphatase-1 inhibitor with a molecular mass of 17-kDa (CPI-17). In addition, we successfully expressed bovine pancreatic trypsin inhibitor (BPTI), which contains three pairs of disulfide bonds, as the soluble form. In the most efficient case, as little as 50 ml culture yielded 3-4 mg (15)N-labeled protein suitable for NMR experiments. The (1)H-(15)N HSQC spectra of all of these proteins clearly indicated that their structures were identical to those of their counterparts reported previously. Thus, the present results suggest that our novel protocol is a potential method for NMR sample preparation.

Crystal structure of an extensively simplified variant of bovine pancreatic trypsin inhibitor in which over one-third of the residues are alanines

We report the high-resolution crystal structures of an extensively simplified variant of bovine pancreatic trypsin inhibitor containing 20 alanines (BPTI-20st) and a reference single-disulfide-bonded variant (BPTI-[5,55]st) at, respectively, 1.39 and 1.09 A resolutions. The sequence was simplified based on the results of an alanine scanning experiment, as reported previously. The effects of the multiple alanine substitutions on the overall backbone structure were surprisingly small (C(alpha) atom RMSD of 0.53 A) being limited to small local structural perturbations. Both BPTI variants retained a wild-type level of trypsin inhibitory activity. The side-chain configurations of residues buried in the hydrophobic cores (<30% accessible surface area) were almost perfectly retained in both BPTI-20st and BPTI-[5,55]st, indicating that neither multiple alanine replacements nor the removal of the disulfide bonds affected their precise placements. However, the side chains of three partially buried residues (Q31, R20, and to some extent Y21) and several unburied residues rearranged into alternative dense-packing structures, suggesting some plasticity in their shape complementarity. These results indicate that a protein sequence simplified over its entire length can retain its densely packed, native side-chain structure, and suggest that both the design and fold recognition of natively folded proteins may be easier than previously thought.

Effect of phospholipids on conformational structure of bovine pancreatic trypsin inhibitor (BPTI) and its thermolabile mutants

The effects of negatively charged phosphatidylserine-prepared membranes (PS) and neutral phosphatidylcholine-prepared membranes (PC) on the structure of wild-type and mutant bovine pancreatic trypsin inhibitor (BPTI) at neutral pH were investigated. The presence of PC did not have any effect on the protein structure while PS induced a non-native structure in three mutant BPTI proteins. However, the negatively charged membrane did not have any effect on wild-type BPTI. The findings revealed that (i) elimination of some disulphide bonds results in dramatic change in protein structure, and, (ii) that this biochemical interaction is surface-driven and electrostatic interactions may play a very strong role in influencing the fore-stated changes in protein structure. Of further interest were the results obtained from investigating the possible role of PS fluidity and concentration in altering mutant. When the value of Gibbs free-energy change of unfolding (DeltaG(U)) was positive, various non-native structures were formed in a concentration-dependent manner. However, when the value of DeltaG(U) was negative, only two types of non-native structures were formed: one with high beta structure content at low PS fluidity state, and the other with a high alpha-helical content at high PS fluidity state. Our study reveals how particular combinations of phospholipid:protein interactions can induce a protein conformation transition from a native to a non-native one at neutral pH, especially when the native structure is predestabilized by amino acid substitutions. This revelation may open up opportunities to explore alternative ways in which phospholipids may play a role in protein mis-folding and the related pathologies.

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.

Relationship between local structural entropy and protein thermostability

We developed a technique to compute structural entropy directly from protein sequences. We explored the possibility of using structural entropy to identify residues involved in thermal stabilization of various protein families. Examples include methanococcal adenylate kinase, Ribonuclease HI and holocytochrome c(551). Our results show that the positions of the largest structural entropy differences between wild type and mutant usually coincide with the residues relevant to thermostability. We also observed a good linear relationship between the average structural entropy and the melting temperatures for adenylate kinase and its chimeric constructs. To validate this linear relationship, we compiled a large dataset comprised of 1153 sequences and found that most protein families still display similar linear relationships. Our results suggest that the multitude of interactions involved in thermal stabilization may be generalized into the tendency of proteins to maintain local structural conservation. The linear relationship between structural entropy and protein thermostability should be useful in the study of protein thermal stabilization.

Recent advances in our understanding of protein conformational stability from a pharmaceutical perspective

The marginal conformational stability of proteins has made them in some cases less than ideal candidates for pharmaceutical agents. Recent progress in our understanding of protein structure and stability has provided the opportunity to design the desired degree of stability into protein drug candidates. Modifications such as the optimisation of interior side-chain packing, the introduction of new ion-pairs, as well as the design of stabilising disulfide bridges and ligand binding sites, all offer the opportunity to produce proteins with enhanced stability properties.

Side-chain Contributions to Membrane Protein Structure and Stability

The molecular forces that stabilize membrane protein structure are poorly understood. To investigate these forces we introduced alanine substitutions at 24 positions in the B helix of bacteriorhodopsin and examined their effects on structure and stability. Although most of the results can be rationalized in terms of the folded structure, there are a number of surprises. (1) We find a remarkably high frequency of stabilizing mutations (17%), indicating that membrane proteins are not highly optimized for stability. (2) Helix B is kinked, with the kink centered around Pro50. The P50A mutation has no effect on stability, however, and a crystal structure reveals that the helix remains bent, indicating that tertiary contacts dominate in the distortion of this helix. (3) We find that the protein is stabilized by about 1kcal/mol for every 38A(2) of surface area buried, which is quite similar to soluble proteins in spite of their dramatically different environments. (4) We find little energetic difference, on average, in the burial of apolar surface or polar surface area, implying that van der Waals packing is the dominant force that drives membrane protein folding.

Electrostatic interactions in the reconstitution of an SH2 domain from constituent peptide fragments

Fragment complementation has been used to delineate the essential recognition elements for stable folding in Src homology 2 (SH2) domains by using NMR spectroscopy, alanine scanning, and surface plasmon resonance. The unfolded 9-kD and 5-kD peptide fragments formed by limited proteolytic digestion of the N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3'-kinase fold into an active native-like structure on interaction with one another. The corresponding 5-kD fragment of the homologous Src protein, however, was not capable of structurally complementing the p85 9-kD fragment, indicating that fragment complementation among these SH2 domains is sensitive to the sequence differences between the Src and p85 domains. Partial complementation and folding activity could be recovered with hybrid sequences of these SH2 domains. Complete alanine scanning of the 5-kD p85 fragment was used to identify the sequence recognition elements required for complex formation. The alanine substitutions in the p85 5-kD fragment that abolished binding affinity with the cognate 9-kD fragment correlate well with highly conserved residues among SH2 domains that are either integrally involved in core packing or found at the interface between fragments. Surprisingly, however, mutation of a nonconserved surface-exposed aspartic acid to alanine was found to have a significant effect on complementation. A single additional mutation of arginine to aspartic acid allowed for recovery of native structure and increased the thermal stability of the designed Src-p85 chimera by 18 degrees C. This modification appears to relieve an unfavorable surface electrostatic interaction, demonstrating the importance of surface charge interactions in protein stability.

Screening for stable mutants with amino acid pairs substituted for the disulfide bond between residues 14 and 38 of bovine pancreatic trypsin inhibitor (BPTI)

We have developed a screening method to identify stable protein mutants from a large number of sequences using a cellular quality control system. This method was used to screen amino acid pairs substituted for the disulfide (S-S) bond between residues 14 and 38 of bovine pancreatic trypsin inhibitor. The mutants selected could be divided into two groups: one with mutation C14G and the other with mutation C38V. Although each mutation did not fully compensate for the destabilizing effect of removal of the S-S bond, these mutants have midpoint temperatures of thermal unfolding that are 12-17 degrees C higher than that of the C14A/C38A mutant. This fact indicates that these mutations are better substitutions for the S-S bond than C14A/C38A. The C14G mutants inhibited trypsin more strongly at 37 degrees C than did the C14A/C38A mutant, although bulky amino acids at position 14 largely diminished the inhibitory activity of the C38V mutants. Thermodynamic analysis indicated that the enthalpy of unfolding of the C14G and C38V mutant groups differed considerably, which suggests different stabilizing mechanisms in these two groups. Because renaturation of S-S bonds is often difficult in the large scale production of proteins, this method should provide a useful tool with which to increase the production of recombinant proteins by eliminating S-S bonds with minimum concomitant stability loss.

A helix initiation signal in T4 lysozyme identified by polyalanine mutagenesis

To better understand the relation between sequence and structure, and in an attempt to simplify the protein folding problem, a series of alanine substitutions was introduced into bacteriophage T4 lysozyme. In contrast to previous studies in this system, which were restricted to single alpha-helices, the present analysis included a helix-turn-helix region, a loop-helix region, and two alpha-helices that were well separated in the three-dimensional structure. It was shown previously that T4 lysozyme is very tolerant of alanine substitutions within alpha-helices, especially at solvent-exposed sites. The present study shows that the protein is also tolerant of such substitutions in turn and loop regions, although less than in helices. The results confirm that the structural information in the amino acid sequence is highly redundant. For example, the protein with the sequence 127AAAAAALAAAAWAAA141 folds normally, has melting temperature only 0.8 degrees C lower than wildtype, and has a crystal structure that is also very similar to wildtype. Polyalanine substitutions within turns or loops can, however, lead to differences in structure and in folding. In one example the triple substitution K35A/S36A/P37A caused this region of the molecule to change to a more helical conformation. In a second case the mutant with the sequence 34AAAAALAAAKAALAAA49, which spans a loop-helix region, had a dramatically altered thermal unfolding transition, suggesting that this region may tend to form a single, uninterrupted, helix. Substitution of Ala38 in the above construct with aspartic acid caused the unfolding to be more like wildtype, suggesting that residue 38, which is at a helix-capping position in the wildtype structure, provides an initiation signal that is essential in the polyalanine mutant for the correct formation of alpha-helix 39-50. In a typical protein, the information that codes for the 3D structure is presumably distributed over many amino acids. The present results suggest that in simplified sequences the key folding information may be restricted to a subset of critical residues, and so be more readily accessible to experimental analysis.

The Poisson-Boltzmann equation for biomolecular electrostatics: a tool for structural biology

Electrostatics plays a fundamental role in virtually all processes involving biomolecules in solution. The Poisson-Boltzmann equation constitutes one of the most fundamental approaches to treat electrostatic effects in solution. The theoretical basis of the Poisson-Boltzmann equation is reviewed and a wide range of applications is presented, including the computation of the electrostatic potential at the solvent-accessible molecular surface, the computation of encounter rates between molecules in solution, the computation of the free energy of association and its salt dependence, the study of pKa shifts and the combination with classical molecular mechanics and dynamics. Theoretical results may be used for rationalizing or predicting experimental results, or for suggesting working hypotheses. An ever-increasing body of successful applications proves that the Poisson-Boltzmann equation is a useful tool for structural biology and complementary to other established experimental and theoretical methodologies.

Combinatorial mutagenesis to restrict amino acid usage in an enzyme to a reduced set

We developed an effective strategy to restrict the amino acid usage in a relatively large protein to a reduced set with conservation of its in vivo function. The 213-residue Escherichia coli orotate phosphoribosyltransferase was subjected to 22 cycles of segment-wise combinatorial mutagenesis followed by 6 cycles of site-directed random mutagenesis, both coupled with a growth-related phenotype selection. The enzyme eventually tolerated 73 amino acid substitutions: In the final variant, 9 amino acid types (A, D, G, L, P, R, T, V, and Y) occupied 188 positions (88%), and none of 7 amino acid types (C, H, I, M, N, Q, and W) appeared. Therefore, the catalytic function associated with a relatively large protein may be achieved with a subset of the 20 amino acid. The converged sequence also implies simpler constituents for proteins in the early stage of evolution.

Robustness of hen lysozyme monitored by random mutations

We investigated the robustness of hen lysozyme by using random mutant libraries. Six random mutant libraries containing 1, 1.5, 2, 3, 5 and 14 amino acid mutations per hen lysozyme were systematically constructed by varying the concentrations of Mg(2+) and Mn(2+) on polymerase chain reaction. The mutated genes from the six libraries were cloned to a yeast expression vector and a total of 4000 clones were screened on the basis of lysis activity and ELISA employing monoclonal antibody that recognized only lysozyme with native conformation. About 80% of the clones with an average of two amino acid mutations retained active structure. Almost all clones with an average of five mutations lost active structure. On the other hand, 80% of the clones with an average of two amino acid mutations retained both gross conformation and active structure and 24% of the clones with an average of 14 amino acid mutations retained gross conformation. These results show that gross conformation is robust against mutations and so is active structure to a lesser extent.

NMR structures of two variants of bovine pancreatic trypsin inhibitor (BPTI) reveal unexpected influence of mutations on protein structure and stability

Here we determined NMR solution structures of two mutants of bovine pancreatic trypsin inhibitor (BPTI) to reveal structural reasons of their decreased thermodynamic stability. A point mutation, A16V, in the solvent-exposed loop destabilizes the protein by 20 degrees C, in contrast to marginal destabilization observed for G, S, R, L or W mutants. In the second mutant introduction of eight alanine residues at proteinase-contacting sites (residues 11, 13, 17, 18, 19, 34, 37 and 39) provides a protein that denatures at a temperature about 30 degrees C higher than expected from additive behavior of individual mutations. In order to efficiently determine structures of these variants, we applied a procedure that allows us to share data between regions unaffected by mutation(s). NOAH/DYANA and CNS programs were used for a rapid assignment of NOESY cross-peaks, structure calculations and refinement. The solution structure of the A16V mutant reveals no conformational change within the molecule, but shows close contacts between V16, I18 and G36/G37. Thus, the observed 4.3kcal/mol decrease of stability results from a strained local conformation of these residues caused by introduction of a beta-branched Val side-chain. Contrary to the A16V mutation, introduction of eight alanine residues produces significant conformational changes, manifested in over a 9A shift of the Y35 side-chain. This structural rearrangement provides about 6kcal/mol non-additive stabilization energy, compared to the mutant in which G37 and R39 are not mutated to alanine residues.

A partially buried site in homologous HPr proteins is not optimized for stability

The energetic consequences of site-specific replacement of a residue at a partially buried site in the two homologous HPr proteins from Escherichia coli and Bacillus subtilis is described. We determined previously that the replacement of a partially buried Lys residue with Glu at position 49 in E.coli HPr increased the conformational stability of the protein substantially because the side-chain of the latter residue could act as a hydrogen-bond acceptor. Here, we extend this analysis to other side-chains with different chemical properties and abilities to form hydrogen bonds to compare the properties of this position in the backgrounds of two different homologous HPr proteins. We find that the variants with polar residues that can form a tertiary hydrogen bond with a nearby site in the protein are more stable than either hydrophobic residues or polar residues that become buried yet are incapable of forming a new hydrogen bond. Furthermore, the protein with the wild-type residue in each HPr variant is not among the most stable of the proteins studied. These results suggest a general strategy for designing variants in which the overall stability of a protein can be modulated in a defined fashion.

Positive contribution of hydration structure on the surface of human lysozyme to the conformational stability

Water molecules make a hydration structure with the network of hydrogen bonds, covering on the surface of proteins. To quantitatively estimate the contribution of the hydration structure to protein stability, a series of hydrophilic mutant human lysozymes (Val to Ser, Tyr, Asp, Asn, and Arg) modified at three different positions on the surface, which are located in the alpha-helix (Val-110), the beta-sheet (Val-2), and the loop (Val-74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by x-ray crystallography at 100 K, respectively. The introduced polar residues made hydrogen bonds with protein atoms and/or water molecules, sometimes changing the hydration structure around the mutation site. Changes in the stability of the mutant proteins can be evaluated by a unique equation that considers the conformational changes resulting from the substitutions. Using this analysis, the relationship between the changes in the stabilities and the hydration structures for mutant human lysozymes substituted on the surface could be quantitatively estimated. The analysis indicated that the hydration structure on protein surface plays an important role in determining the conformational stability of the protein.

Toward development of a screen to identify randomly encoded, foldable sequences

The ability to identify sequences in a randomly encoded polypeptide library that are capable of acquiring unique and stably folded structures would be valuable in the examination of protein-folding issues. The quality control system of the yeast secretory pathway prevents the release of incompletely folded polypeptides. Earlier work has shown that this feature can be used in a screen to identify mutations that increase the stability of a protein. We sought to extend this strategy for use with random sequence libraries by combining a quality-control system-based screen with generic tag-based immunodetection that can be applied to any sequence. To test this method, we screened a library encoding random mutations in a bovine pancreatic trypsin inhibitor variant containing a small generic tag. Initial on-plate screening resulted in a large number of false positives: sequences that were secreted but not foldable. These false positives were excluded successfully in additional screening steps that used a liquid-culture secretion screen and a gel electrophoresis assay. Three positive clones were obtained that showed midpoint thermal denaturation temperatures 10-16 degrees C higher than the original bovine pancreatic trypsin inhibitor variant. Thus, this multistep screening method may be useful for finding novel, foldable sequences.

Thermodynamics of single peptide bond cleavage in bovine pancreatic trypsin inhibitor (BPTI)

A major goal of this paper was to estimate a dynamic range of equilibrium constant for the opening of a single peptide bond in a model protein, bovine pancreatic trypsin inhibitor (BPTI). Ten mutants of BPTI containing a single Xaa-->Met substitution introduced in different parts of the molecule were expressed in Escherichia coli. The mutants were folded, purified to homogeneity, and cleaved with cyanogen bromide to respective cleaved forms. Conformation of the intact mutants was similar to the wildtype, as judged from their circular dichroism spectra. Substantial conformational changes were observed on the chemical cleavage of three single peptide bonds--Met46-Ser, Met49-Cys, and Met53-Thr--located within the C-terminal helix. Cleavage of those peptide bonds caused a significant destabilization of the molecule, with a drop of the denaturation temperature by 56.4 degrees C to 68 degrees C at pH 4.3. Opening of the remaining seven peptide bonds was related to a 10.8 degrees C to 39.4 degrees C decrease in T(den). Free energies of the opening of 10 single peptide bonds in native mutants (Delta G(op,N)) were estimated from the thermodynamic cycle that links denaturation and cleavage free energies. To calculate those values, we assumed that the free energy of opening of a single peptide bond in the denatured state (Delta G(op,D)) was equal to -2.7 kcal/mole, as reported previously. Calculated Delta G(op,N) values in BPTI were in the range from 0.2 to 10 kcal/mole, which was equivalent to a >1 million-fold difference in equilibrium constants. The values of Delta G(op,N) were the largest for peptide bonds located in the C-terminal helix and significantly lower for peptide bonds in the beta-structure or loop regions. It appears that opening constants for single peptide bonds in various proteins span across 33 orders of magnitude. Typical equilibrium values for a single peptide bond opening in a protein containing secondary structure elements fall into negligibly low values, from 10(-3) to 10(-8), and are efficient to ensure stability against proteolysis.

Analysis of serine proteinase-inhibitor interaction by alanine shaving

We analyzed the energetic importance of residues surrounding the hot spot (the P(1) position) of bovine pancreatic trypsin inhibitor (BPTI) in interaction with two proteinases, trypsin and chymotrypsin, by a procedure called molecular shaving. One to eight residues of the structural epitope, composed of two extended and exposed loops, were mutated to alanine(s). Although truncation of the side chains of residues surrounding the P(1) position to methyl groups caused a decrease in Delta G(den) values up to 6.4 kcal mole(-1), it did not influence the overall conformation of the inhibitor. We found that the replacement of up to six residues with alanines was fully additive at the level of protein stability. To analyze the influence of the structural epitope on the association energy, we determined association constants for BPTI variants and both enzymes and applied the additivity analysis. Shaving of two binding loops led to a progressive drop in the association energy, more pronounced for trypsin (decrease up to 9.6 kcal mole(-1)) than chymotrypsin (decrease up to 3.5 kcal mole(-1)). In the case of extensively mutated variants interacting with chymotrypsin, the association energies agreed very well with the values calculated from single mutational effects. However, when P(1)-neighboring residues were shaved to alanine(s), their contribution to the association energy was not fully removed because of the presence of methyl groups and main chain-main chain intermolecular hydrogen bonds. Moreover, the hot spot had a different contribution to the complex stability in the fully shaved BPTI variant compared with the wild type, which was caused by perturbations of the P(1)-S(1) electrostatic interaction.

Correlative motions and memory effects in molecular dynamics simulations of molecules: principal components and rescaled range analysis suggest that the motions of native BPTI are more correlated than those of its mutants

In this work MD simulations of the native bovine pancreatic trypsin inhibitor (BPTI) and 16 mutants were done in vacuum in order to study memory effects in the mutants using principal component analysis (PCA) and the rescaled range analysis (Hurst exponents). Both PCA and the rescaled range analysis support our previous proposition, based on PCA of lysozyme, that the motions of a native protein are more correlated than those of mutants. The methods are compared, the nature and applications of the rule and the role of the long-range correlations in MD time series (i.e. memory) are discussed in the context of collective motions.

Mutational analysis of hydrogen bonding residues in the BPTI folding pathway

Nine BPTI variants with replacements that remove one or more hydrogen bonds from the native protein were constructed, and the folding pathways of these proteins were determined by isolating and identifying the disulfide-bonded intermediates that accumulated during unfolding and refolding. The forward and reverse rate constants for the individual steps in the folding pathways for each protein were measured, providing a detailed description of the energetic effects of the substitutions. The native forms of eight of the nine variants were measurably destabilized, by 1-7 kcal/mol (1 cal=4.184 J), with an average effect of 1.6 kcal/mol per hydrogen bond removed. The folding pathways for the variants were found to be similar to that previously described for the wild-type protein, with the kinetically preferred mechanism involving intramolecular rearrangements of intermediates with two disulfide bonds. Some of the substitutions, however, significantly destabilized the major intermediates and broadened the distribution of species with one or two disulfide bonds, thus identifying residues that play important roles in stabilizing the normal intermediates and defining specificity in the folding process. The kinetic data also suggest that one residue, Asn43, may play a distinctive role in defining the BPTI folding mechanism. Replacement of this residue with either Gly or Ala appeared to stabilize the major transition states for folding and unfolding. In the native protein, the side-chain of Asn43 participates directly in the hydrogen bonding pattern of the central beta-sheet, and the kinetic behavior of the Asn43 variants suggests that the major energy barriers in folding and unfolding may be due in part to the steric constraints imposed by this structural element, together with those imposed by the chemical transition states for thiol-disulfide exchange.

Are the parameters of various stabilization factors estimated from mutant human lysozymes compatible with other proteins?

The various factors which contribute to protein stability have been extensively examined using mutant proteins, but the same kinds of substitutions have given different results depending on the substitution sites. Recently, the contributions of some stabilization factors have been quantitatively derived as parameters by a unique equation, considering the conformational changes due to the mutations using mutant human lysozymes [Funahashi et al. (1999) Protein ENG: 12, 841-850]. To evaluate these parameters estimated from the mutant human lysozymes, stability-structure datasets for the mutant T4 lysozymes were selected. The stabilities for the mutant T4 lysozymes could be roughly estimated using these parameters. Notable differences between the estimated and experimental stabilities were caused by the uncertainty in part of the structures due to some Arg and Lys residues fluctuating on the surface of the T4 lysozyme. Excluding these atoms from the estimation gave a good correlation between the estimated and experimental stabilities. These results suggest that the parameters of the various stabilization factors derived from the mutant human lysozymes are compatible with the mutant T4 lysozymes, although they should be improved with respect to some points using more information.

Substitutions at the P(1) position in BPTI strongly affect the association energy with serine proteinases

The role of the S(1) subsite in trypsin, chymotrypsin and plasmin has been examined by measuring the association with seven different mutants of bovine pancreatic trypsin inhibitor (BPTI); the mutants contain Gly, Ala, Ser, Val, Leu, Arg, and Trp at the P(1) position of the reactive site. The effects of substitutions at the P(1) position on the association constants are very large, comprising seven orders of magnitude for trypsin and plasmin, and over five orders for chymotrypsin. All mutants showed a decrease of the association constant to the three proteinases in the same order: Ala>Gly>Ser>Arg>Val>Leu>Trp. Calorimetric and circular dichroism methods showed that none of the P1 substitutions, except the P1-Val mutant, lead to destabilisation of the binding loop conformation. The X-ray structure of the complex formed between bovine beta-trypsin and P(1)-Leu BPTI showed that the P(1)-Leu sterically conflicts with the side-chain of P(3)-Ile, which thereby is forced to rotate approximately 90 degrees. Ile18 (P(3)) in its new orientation, in turn interacts with the Tyr39 side-chain of trypsin. Introduction of a large side-chain at the P1' position apparently leads to a cascade of small alterations of the trypsin-BPTI interface that seem to destabilise the complex by it adopting a less optimized packing and by tilting the BPTI molecule up to 15 degrees compared to the native trypsin-BPTI complex.

Folding of bovine pancreatic trypsin inhibitor (BPTI) variants in which almost half the residues are alanine

Recent studies indicate that a fraction of the information contained in an amino acid sequence may be sufficient for specifying a native protein structure. An earlier alanine-scanning experiment conducted on bovine pancreatic trypsin inhibitor (BPTI; 58 residues) suggested that if cumulative mutations have additive effects on protein stability, a native protein structure could be built from BPTI sequences that contained many alanine residues distributed throughout the protein. To test this hypothesis, we designed and produced six BPTI mutants containing from 21 to 29 alanine residues. We found that the melting temperature of mutants containing up to 27 alanine residues (48 % of the total number of residues) could be predicted quite well by the sum of the change in melting temperature for the single mutations. Additionally, these same mutants folded into a native-like structure, as judged by their cooperative thermal denaturation curves and heteronuclear multiple quantum correlation (HMQC) NMR spectra. A BPTI mutant containing 22 alanine residues was further shown by 2D and 3D-NMR to fold into a structure very similar to that of native BPTI, and to be a functional trypsin inhibitor. These results provide insight into the extent to which native protein structure and function can be achieved with a highly simplified amino acid sequence.

Disulfide bonds and protein folding

The applications of disulfide-bond chemistry to studies of protein folding, structure, and stability are reviewed and illustrated with bovine pancreatic ribonuclease A (RNase A). After surveying the general properties and advantages of disulfide-bond studies, we illustrate the mechanism of reductive unfolding with RNase A, and discuss its application to probing structural fluctuations in folded proteins. The oxidative folding of RNase A is then described, focusing on the role of structure formation in the regeneration of the native disulfide bonds. The development of structure and conformational order in the disulfide intermediates during oxidative folding is characterized. Partially folded disulfide species are not observed, indicating that disulfide-coupled folding is highly cooperative. Contrary to the predictions of "rugged funnel" models of protein folding, misfolded disulfide species are also not observed despite the potentially stabilizing effect of many nonnative disulfide bonds. The mechanism of regenerating the native disulfide bonds suggests an analogous scenario for conformational folding. Finally, engineered covalent cross-links may be used to assay for the association of protein segments in the folding transition state, as illustrated with RNase A.

Remarkable thermal stability of doubly intramolecularly cross-linked hen lysozyme

In order to examine how a protein can be effectively stabilized, two intramolecular cross-links, Glu35-Trp108 and Lys1-His15, which have few unfavorable interactions in the folded state, were simultaneously introduced into hen lysozyme. Both of the intramolecularly cross-linked lysozymes, 35-108 CL and 1-15 CL, containing cross-links Glu35-Trp108 and Lys1-His15, respectively, showed increases in thermal stability of 13.9 and 5.2 degrees C, respectively, over that of wild type, at pH 2.7. On the other hand, a doubly cross-linked lysozyme showed an increase in thermal stability of 20.8 degrees C over that of wild type, under identical conditions. Since the sum of the differences in denaturation temperature between wild type and each of the cross-linked lysozymes was nearly equal to that between wild type and the doubly cross-linked lysozyme, we suggest that the efficient stabilization of the lysozyme molecule was the direct result of the double intramolecular cross-links.

Do proteins learn to evolve? The Hopfield network as a basis for the understanding of protein evolution

Correlations between amino-acid residues can be observed in sets of aligned protein sequences, and the analysis of their statistical and evolutionary significance and distribution has been thoroughly investigated. In this paper, we present a model based on such covariations in protein sequences in which the pairs of residues that have mutual influence combine to produce a system analogous to a Hopfield neural network. The emergent properties of such a network, such as soft failure and the connection between network architecture and stored memory, have close parallels in known proteins. This model suggests that an explanation for observed characters of proteins such as the diminution of function by substitutions distant from the active site, the existence of protein folds (superfolds) that can perform several functions based on one architecture, and structural and functional resilience to destabilizing substitutions might derive from their inherent network-like structure. This model may also provide a basis for mapping the relationship between structure, function and evolutionary history of a protein family, and thus be a powerful tool for rational engineering.

Increasing protein stability by altering long-range coulombic interactions

It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.

Early events in the disulfide-coupled folding of BPTI

Recent studies of the refolding of reduced bovine pancreatic trypsin inhibitor (BPTI) have shown that a previously unidentified intermediate with a single disulfide is formed much more rapidly than any other one-disulfide species. This intermediate contains a disulfide that is present in the native protein (between Cys14 and 38), but it is thermodynamically less stable than the other two intermediates with single native disulfides. To characterize the role of the [14-38] intermediate and the factors that favor its formation, detailed kinetic and mutational analyses of the early disulfide-formation steps were carried out. The results of these studies indicate that the formation of [14-38] from the fully reduced protein is favored by both local electrostatic effects, which enhance the reactivities of the Cys14 and 38 thiols, and conformational tendencies that are diminished by the addition of urea and are enhanced at lower temperatures. At 25 degrees C and pH 7.3, approximately 35% of the reduced molecules were found to initially form the 14-38 disulfide, but the majority of these molecules then undergo intramolecular rearrangements to generate non-native disulfides, and subsequently the more stable intermediates with native disulfides. Amino acid replacements, other than those involving Cys residues, were generally found to have only small effects on either the rate of forming [14-38] or its thermodynamic stability, even though many of the same substitutions greatly destabilized the native protein and other disulfide-bonded intermediates. In addition, those replacements that did decrease the steady-state concentration of [14-38] did not adversely affect further folding and disulfide formation. These results suggest that the weak and transient interactions that are often detected in unfolded proteins and early folding intermediates may, in some cases, not persist or promote subsequent folding steps.