Thermodynamics of BPTI folding

Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Despite advances in artificial intelligence methods, protein folding remains in many ways an enigma to be solved. Accurate computation of protein folding energetics could help drive fields such as protein and drug design and genetic interpretation. However, the challenge of calculating the state functions governing protein folding from first-principles remains unaddressed. We present here a simple approach that allows us to accurately calculate the energetics of protein folding. It is based on computing the energy of the folded and unfolded states at different temperatures using molecular dynamics simulations. From this, two essential quantities (ΔH and ΔCp) are obtained and used to calculate the conformational stability of the protein (ΔG). With this approach, we have successfully calculated the energetics of two- and three-state proteins, representatives of the major structural classes, as well as small stability differences (ΔΔG) due to changes in solution conditions or variations in an amino acid residue.

Implications of a temperature-dependent heat capacity for temperature-gated ion channels

Temperature influences dynamics and state-equilibrium distributions in all molecular processes, and only a relatively narrow range of temperatures is compatible with life-organisms must avoid temperature extremes that can cause physical damage or metabolic disruption. Animals evolved a set of sensory ion channels, many of them in the family of transient receptor potential cation channels that detect biologically relevant changes in temperature with remarkable sensitivity. Depending on the specific ion channel, heating or cooling elicits conformational changes in the channel to enable the flow of cations into sensory neurons, giving rise to electrical signaling and sensory perception. The molecular mechanisms responsible for the heightened temperature-sensitivity in these ion channels, as well as the molecular adaptations that make each channel specifically heat- or cold-activated, are largely unknown. It has been hypothesized that a heat capacity difference (ΔCp) between two conformational states of these biological thermosensors can drive their temperature-sensitivity, but no experimental measurements of ΔCp have been achieved for these channel proteins. Contrary to the general assumption that the ΔCp is constant, measurements from soluble proteins indicate that the ΔCp is likely to be a function of temperature. By investigating the theoretical consequences for a linearly temperature-dependent ΔCp on the open-closed equilibrium of an ion channel, we uncover a range of possible channel behaviors that are consistent with experimental measurements of channel activity and that extend beyond what had been generally assumed to be possible for a simple two-state model, challenging long-held assumptions about ion channel gating models at equilibrium.

Fluorine-induced polarity increases inhibitory activity of BPTI towards chymotrypsin

Substituting the P1 position in bovine pancreatic trypsin inhibitor (BPTI) is known to heavily influence its inhibitory activity towards serine proteases. Side-chain fluorinated aliphatic amino acids have been shown to alter numerous properties of peptides and proteins and thus are of interest in the context of BPTI. In our study, we systematically investigated the site-specific incorporation of non-canonical amino acids into BPTI by microwave-assisted solid-phase peptide synthesis (SPPS). Inhibitor activity of the variants was tested towards the serine protease α-chymotrypsin. We observed enhanced inhibition of two fluorinated BPTIs compared to wild type and hydrocarbon variants. To further investigate the complexes, we performed X-ray structure analysis. Our findings underline the power fluorine offers as a tool in protein engineering to beneficially alter the effects on phenomena as protein-protein interactions.

Why small proteins tend to have high denaturation temperatures

Data indicate that small globular proteins (consisting of less than about 70 residues) tend to have high denaturation temperatures. This finding is analysed by comparing experimental denaturation enthalpy and entropy changes of a selected set of small proteins with values calculated on the basis of average and common properties of globular proteins. The conclusion is that the denaturation entropy change is smaller than expected, leading to an increase in denaturation temperature. The proposed molecular rationalization considers the existence of long-wavelength, low-frequency vibrational modes in the native state of small proteins due to their large surface-to-interior ratio. The effect of decreasing the conformational entropy gain associated with denaturation on thermal stability is directly verified by means of an already devised theoretical model [G. Graziano, Phys. Chem. Chem. Phys. 2010, 12, 14245-14252; 2014, 16, 21755-21767].

A Second Backbone: The Contribution of a Buried Asparagine Ladder to the Global and Local Stability of a Leucine-Rich Repeat Protein

Parallel β-sheet-containing repeat proteins often display a structural motif in which conserved asparagines form a continuous ladder buried within the hydrophobic core. In such "asparagine ladders", the asparagine side-chain amides form a repetitive pattern of hydrogen bonds with neighboring main-chain NH and CO groups. Although asparagine ladders have been thought to be important for stability, there is little experimental evidence to support such speculation. Here we test the contribution of a minimal asparagine ladder from the leucine-rich repeat protein pp32 to stability and investigate lattice rigidity and hydrogen bond character using solution nuclear magnetic resonance (NMR) spectroscopy. Point substitutions of the two ladder asparagines of pp32 are strongly destabilizing and decrease the cooperativity of unfolding. The chemical shifts of the ladder side-chain HZ protons are shifted significantly downfield in the NMR spectrum and have low temperature coefficients, indicative of strong hydrogen bonding. In contrast, the HE protons are shifted upfield and have temperature coefficients close to zero, suggesting an asymmetry in hydrogen bond strength along the ladder. Ladder NH2 groups have weak 1H-15N cross-peak intensities; 1H-15N nuclear Overhauser effect and 15N CPMG experiments show this to be the result of high rigidity. Hydrogen exchange measurements demonstrate that the ladder NH2 groups exchange very slowly, with rates approaching the global exchange limit. Overall, these results show that the asparagine side chains are held in a very rigid, nondynamic structure, making a significant contribution to the overall stability. In this regard, buried asparagine ladders can be considered "second backbones" within the cores of their elongated β-sheet host proteins.

Implementation of an experimental and computational tool set to study protein-mAb interactions

This work focused on the development of a combined experimental and computational tool set to study protein-mAb interactions. A model protein library was first screened using cross interaction chromatography to identify proteins with the strongest retention. Fluorescence polarization was then employed to study the interactions and thermodynamics of the selected proteins-lactoferrin, pyruvate kinase, and ribonuclease B with the mAb. Binding affinities of lactoferrin and pyruvate kinase to the mAb were seen to be relatively salt insensitive in the range examined. Further, a strong entropic contribution was observed, suggesting the importance of hydrophobic interactions. On the other hand, ribonuclease B-mAb binding was seen to be enthalpically driven and salt sensitive, indicating the importance of electrostatic interactions. Protein-protein docking was then carried out and the results identified the CDR region on the mAb as an important binding site for all three proteins. The binding interfaces identified for the mAb-lactoferrin and mAb-pyruvate kinase systems were found to contain complementary hydrophobic and oppositely charged clusters on the interacting regions which were indicative of both hydrophobic and electrostatic interactions. On the other hand, the binding site on ribonuclease B was predominantly positively charged with minimal hydrophobicity. This resulted in an alignment with negatively charged clusters on the mAb, supporting the contention that these interactions were primarily electrostatic in nature. Importantly, these computational results were found to be consistent with the fluorescence polarization data and this combined approach may have utility in examining mAb-HCP interactions which can often complicate the downstream processing of biologics. © 2019 American Institute of Chemical Engineers.

BPTI folding revisited: switching a disulfide into methylene thioacetal reveals a previously hidden path

The folding mechanism of the model protein bovine pancreatic trypsin inhibitor was revisited. By switching the solvent exposed disulfide bond with methylene thioacetal we uncovered a hidden pathway in its folding mechanism. In addition, this moiety enhanced protein stability while fully maintaining the protein structure and biological function.

Bovine pancreatic trypsin inhibitor (BPTI) is a 58-residue protein that is stabilized by three disulfide bonds at positions 5–55, 14–38 and 30–51. Widely studied for about 50 years, BPTI represents a folding model for many disulfide-rich proteins. In the study described below, we replaced the solvent exposed 14–38 disulfide bond with a methylene thioacetal bridge in an attempt to arrest the folding pathway of the protein at its two well-known intermediates, N′ and N*. The modified protein was expected to be unable to undergo the rate-determining step in the widely accepted BPTI folding mechanism: the opening of the 14–38 disulfide bond followed by rearrangements that leads to the native state, N. Surprisingly, instead of halting BPTI folding at N′ and N*, we uncovered a hidden pathway involving a direct reaction between the N* intermediate and the oxidizing reagent glutathione (GSSG) to form the disulfide-mixed intermediate N*–SG, which spontaneously folds into N. On the other hand, N′ was unable to fold into N. In addition, we found that the methylene thioacetal bridge enhances BPTI stability while fully maintaining its structure and biological function. These findings suggest a general strategy for enhancing protein stability without compromising on function or structure, suggesting potential applications for future therapeutic protein production.

Can Local Probes Go Global? A Joint Experiment-Simulation Analysis of λ(6-85) Folding

The process of protein folding is known to involve global motions in a cooperative affair; the structure of most of the protein sequences is gained or lost over a narrow range of temperature, denaturant, or pressure perturbations. At the same time, recent simulations and experiments reveal a complex structural landscape with a rich set of local motions and conformational changes. We couple experimental kinetic and thermodynamic measurements with specifically tailored analysis of simulation data to isolate local versus global folding probes. We find that local probes exhibit lower melting temperatures, smaller surface area changes, and faster kinetics compared to global ones. We also see that certain local probes of folding match the global behavior more closely than others. Our work highlights the importance of using multiple probes to fully characterize protein folding dynamics by theory and experiment.

Applications of pressure perturbation calorimetry to study factors contributing to the volume changes upon protein unfolding

BACKGROUND

Pressure perturbation calorimetry (PPC) is a biophysical method that allows direct determination of the volume changes upon conformational transitions in macromolecules.

SCOPE OF THIS REVIEW

This review provides novel details of the use of PPC to analyze unfolding transitions in proteins. The emphasis is made on the data analysis as well as on the validation of different structural factors that define the volume changes upon unfolding. Four case studies are presented that show the application of these concepts to various protein systems.

MAJOR CONCLUSIONS

The major conclusions are: 1. Knowledge of the thermodynamic parameters for heat induced unfolding facilitates the analysis of the PPC profiles. 2. The changes in the thermal expansion coefficient upon unfolding appear to be temperature dependent.3.Substitutions on the protein surface have negligible effects on the volume changes upon protein unfolding. 4. Structural plasticity of proteins defines the position dependent effect of amino acid substitutions of the residues buried in the native state. 5. Small proteins have positive volume changes upon unfolding which suggests difference in balance between the cavity/void volume in the native state and the hydration volume changes upon unfolding as compared to the large proteins that have negative volume changes.

GENERAL SIGNIFICANCE

The information provided here gives a better understanding and deeper insight into the role played by various factors in defining the volume changes upon protein unfolding.

Fluorine teams up with water to restore inhibitor activity to mutant BPTI

Fluorinated derivatives of aminobutyric acid engage in unique interactions with structural waters within the BPTI/trypsin interface and restore inhibitor activity.

Introducing fluorine into molecules has a wide range of effects on their physicochemical properties, often desirable but in most cases unpredictable. The fluorine atom imparts the C–F bond with low polarizability and high polarity, and significantly affects the behavior of neighboring functional groups, in a covalent or noncovalent manner. Here, we report that fluorine, present in the form of a single fluoroalkyl amino acid side chain in the P1 position of the well-characterized serine-protease inhibitor BPTI, can fully restore inhibitor activity to a mutant that contains the corresponding hydrocarbon side chain at the same site. High resolution crystal structures were obtained for four BPTI variants in complex with bovine β-trypsin, revealing changes in the stoichiometry and dynamics of water molecules in the S1 subsite. These results demonstrate that the introduction of fluorine into a protein environment can result in “chemical complementation” that has a significantly favorable impact on protein–protein interactions.

Infinitely dilute partial molar properties of proteins from computer simulation

A detailed understanding of temperature and pressure effects on an infinitely dilute protein's conformational equilibrium requires knowledge of the corresponding infinitely dilute partial molar properties. Established molecular dynamics methodologies generally have not provided a way to calculate these properties without either a loss of thermodynamic rigor, the introduction of nonunique parameters, or a loss of information about which solute conformations specifically contributed to the output values. Here we implement a simple method that is thermodynamically rigorous and possesses none of the above disadvantages, and we report on the method's feasibility and computational demands. We calculate infinitely dilute partial molar properties for two proteins and attempt to distinguish the thermodynamic differences between a native and a denatured conformation of a designed miniprotein. We conclude that simple ensemble average properties can be calculated with very reasonable amounts of computational power. In contrast, properties corresponding to fluctuating quantities are computationally demanding to calculate precisely, although they can be obtained more easily by following the temperature and/or pressure dependence of the corresponding ensemble averages.

[Conformational stability of serine proteinase inhibitor from the sea anemone Heteractis crispa]

The influence of different environmental values of the pH and temperature on the spatial organization of serine proteinase inhibitor from the sea anemone Heteractis crispa (=Radianthus macrodactylus) on the level of tertiary and secondary structure was studied by CD spectroscopy. The molecule InhVJ was shown to possess a high conformational thermo- and pH-stability. We determined the point of conformational thermotransition of polypeptide (70 degrees C) after which the molecule gets denaturational stable state with conservation of 80% proteinase inhibitory activity. The significant partial reversible changes of molecule spatial organization were established to occur at the level of tertiary structure in the process of acid-base titration in the range of pH 11.0-13.0. This can be explained by of ionization of tyrosine residues. The molecule InhVJ is conformationally stable at the low pH values (2.0). The quenching of tyrosine residues by acrylamide showed that two of these residues are accessible to the quencher in full, while the third part is available.

Brownian dynamics simulation of protein solutions: structural and dynamical properties

The study of solutions of biomacromolecules provides an important basis for understanding the behavior of many fundamental cellular processes, such as protein folding, self-assembly, biochemical reactions, and signal transduction. Here, we describe a Brownian dynamics simulation procedure and its validation for the study of the dynamic and structural properties of protein solutions. In the model used, the proteins are treated as atomically detailed rigid bodies moving in a continuum solvent. The protein-protein interaction forces are described by the sum of electrostatic interaction, electrostatic desolvation, nonpolar desolvation, and soft-core repulsion terms. The linearized Poisson-Boltzmann equation is solved to compute electrostatic terms. Simulations of homogeneous solutions of three different proteins with varying concentrations, pH, and ionic strength were performed. The results were compared to experimental data and theoretical values in terms of long-time self-diffusion coefficients, second virial coefficients, and structure factors. The results agree with the experimental trends and, in many cases, experimental values are reproduced quantitatively. There are no parameters specific to certain protein types in the interaction model, and hence the model should be applicable to the simulation of the behavior of mixtures of macromolecules in cell-like crowded environments.

Natural and engineered kallikrein inhibitors: an emerging pharmacopoeia

The kallikreins and kallikrein-related peptidases are serine proteases that control a plethora of developmental and homeostatic phenomena, ranging from semen liquefaction to skin desquamation and blood pressure. The diversity of roles played by kallikreins has stimulated considerable interest in these enzymes from the perspective of diagnostics and drug design. Kallikreins already have well-established credentials as targets for therapeutic intervention and there is increasing appreciation of their potential both as biomarkers and as targets for inhibitor design. Here, we explore the current status of naturally occurring kallikrein protease-inhibitor complexes and illustrate how this knowledge can interface with strategies for rational re-engineering of bioscaffolds and design of small-molecule inhibitors.

Temperature dependence on structure and dynamics of Bovine Pancreatic Trypsin Inhibitor (BPTI): a neutron scattering study

We have studied the influence of temperature on the structure of BPTI in solution by small angle neutron scattering. We have investigated the variation of the radius of gyration and the modification of the shape of BPTI between ambient temperature and 368 K. Results have shown an increase of the radius of gyration from 10.9 A at ambient temperature up to 13.3 A at 368 K. Global and internal dynamics of BPTI in solution were studied by quasielastic neutron scattering. The analysis of neutron data in terms of intermediate scattering function reveals two relaxation times tau(1) and tau(2) related respectively to global translational diffusive motions and to internal motions of protein. Motions of protons belonging to lateral chains of residues located at the surface of the protein have been detected. The results are compared to the recently published results concerning the influence of pressure on structure and dynamics of BPTI in solution [Appavou MS et al. Biochimica et Biophysica Acta, 1764, 2006, pp 414-423].

A computational approach for the rational design of stable proteins and enzymes: optimization of surface charge-charge interactions

The design of stable proteins and enzymes is not only of particular biotechnological importance, but also addresses some important fundamental questions. While there are a number of different options available for designing or engineering stable proteins, the field of computational design provides fast and universal methods for stabilizing proteins of interest. One of the successful computational design strategies focuses on stabilizing proteins through the optimization of charge-charge interactions on the protein surface. By optimizing surface interactions, it is possible to alleviate some of the challenges that accompany efforts to redesign the protein core. The rational design of surface charge-charge interactions also allows one to optimize only the interactions that are distant from binding sites or active sites, making it possible to increase stability without adversely affecting activity. The optimization of surface charge-charge interactions is discussed in detail along with the experimental evidence to demonstrate that this is a robust and universal approach to designing proteins with enhanced stability.

Microcalorimetry of proteins and their complexes

Ultrasensitive microcalorimetric techniques for measuring the heat capacities of proteins in dilute solutions over a broad temperature range (DSC) and the heats of protein reactions at fixed temperatures (ITC) are described and the methods of working with these instruments are considered. Particular attention is paid to analyzing the thermal properties of individual proteins, their stability, the energetics of their folding, and their association with specific macromolecular partners. Use of these calorimetric methods is illustrated with examples of small compact globular proteins, small proteins having loose noncompact structure, multidomain proteins, and protein complexes, particularly with DNA.

Estimating the temperature dependence of peptide folding entropies and free enthalpies from total energies in molecular dynamics simulations

The temperature dependence of thermodynamic quantities, such as heat capacity, entropy and free enthalpy, may be obtained by using well-known equations that relate these quantities to the enthalpy of the molecular system of interest at a range of temperatures. In turn, the enthalpy of a molecular system can be estimated from molecular dynamics simulations of an appropriate model. To demonstrate this, we have investigated the temperature dependence of the enthalpy, heat capacity, entropy and free enthalpy of a system that consists of a beta-heptapeptide in methanol and have used the statistical mechanics relationships to describe the thermodynamics of the folding/unfolding equilibrium of the peptide. The results illustrate the power of current molecular simulation force fields and techniques in establishing the link between thermodynamic quantities and conformational distributions.

Mechanism of low density lipoprotein (LDL) release in the endosome: implications of the stability and Ca2+ affinity of the fifth binding module of the LDL receptor

Uptake of low density lipoproteins (LDL) by their receptor, LDLR, is the primary mechanism by which cells incorporate cholesterol from plasma. Mutations in LDLR lead to familial hypercholesterolemia, a common disease affecting 1 in 500 of the human population. LDLR is a modular protein that uses several small repeats to bind LDL. The repeats contain around 40 residues, including three disulfide bonds and a calcium ion. Repeat 5 (LR5) is critical for LDL and beta-migrating very low density lipoprotein binding. Based on the crystal structure of LDLR at endosomal pH (but close to extracellular calcium concentration), LR5 has been proposed to bind to the epidermal growth factor (EGF) precursor domain of LDLR in the endosome, thus releasing the LDL particles previously bound in extracellular conditions. We report here the conformational stability of LR5 as a function of temperature and calcium concentration under both extracellular and endosomal pH conditions. The repeat was very stable when it bore a bound calcium ion but was severely destabilized in the absence of calcium and even further destabilized at acidic versus neutral pH. The temperature and calcium concentration dependence of LR5 stability clearly indicate that under endosomal conditions the unfolded conformation of the repeat is largely dominant. We thus propose a new mechanism for LDL release in the endosome in which calcium depletion and decreased stability at acidic pH drives LR5 unfolding, which triggers LDL release from the receptor. Subsequent binding of LR5 to the EGF precursor domain, if it takes place at low calcium concentrations, would contribute to a further shifting of the equilibrium toward dissociation.

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.

Influence of pressure on structure and dynamics of bovine pancreatic trypsin inhibitor (BPTI): small angle and quasi-elastic neutron scattering studies

We have studied the influence of pressure on structure and dynamics of a small protein belonging to the enzymatic catalysis: the bovine pancreatic trypsin inhibitor (BPTI). Using a copper-beryllium high-pressure cell, we have performed small angle neutron scattering (SANS) experiment on NEAT spectrometer at HMI (Berlin, Germany). In the SANS configuration, the evolution of the radius of gyration and of the shape of the protein under pressures up to 6,000 bar has been studied. When increasing pressure from atmospheric pressure up to 6,000 bar, the pressure effects on the global structure of BPTI result on a reduction of the radius of gyration from 13.4 A down to 12.0 A. Between 5,000 and 6,000 bar, some transition already detected by FTIR [N. Takeda, K. Nakano, M. Kato, Y. Taniguchi, Biospectroscopy, 4, 1998, pp. 209-216] is observed. The pressure effect is not reversible because the initial value of the radius of gyration is not recovered after pressure release. By extending the range of wave-vectors to high q, we have observed a change of the form factor (shape) of the BPTI under pressure. At atmospheric pressure BPTI exhibits an ellipsoidal form factor that is characteristic of the native state. When the pressure is increased from atmospheric pressure up to 6,000 bar, the protein keeps its ellipsoidal shape. The parameters of the ellipsoid vary and the transition detected between 5,000 and 6,000 bar in the form factor of BPTI is in agreement with the FTIR results. After pressure release, the form factor of BPTI is characteristic of an ellipsoid of revolution with a semi-axis a, slightly elongated with respect to that of the native one, indicating that the pressure-induced structural changes on the protein are not reversible. The global motions and the internal dynamics of BPTI protein have been investigated in the same pressure range by quasi-elastic neutron scattering experiments on IN5 time-of-flight spectrometer at ILL (Grenoble, France). The diffusion coefficients D and the internal relaxation times <tau(2)> of BPTI deduced from the analysis of the intermediate scattering functions show a slowing down of protein dynamics when increasing pressure.

pH dependence thermal stability of a chymotrypsin inhibitor from Schizolobium parahyba seeds

The thermal stability of a Schizolobium parahyba chymotrypsin inhibitor (SPCI) as a function of pH has been investigated using fluorescence, circular dichroism, and differential scanning calorimetry (DSC). The thermodynamic parameters derived from all methods are remarkably similar and strongly suggest the high stability of SPCI under a wide range of pH. The transition temperature (T(m)) values ranging from 57 to 85.3 degrees C at acidic, neutral, and alkaline pH are in good agreement with proteins from mesophilic and thermophilic organisms and corroborate previous data regarding the thermal stability of SPCI. All methods gave transitions curves adequately fitted to a two-state model of the unfolding process as judged by the cooperative ratio between the van't Hoff and the calorimetric enthalpy energies close to unity in all of the pH conditions analyzed, except at pH 3.0. Thermodynamic analysis using all these methods reveals that SPCI is thermally a highly stable protein, over the wide range of pH from 3.0 to 8.8, exhibiting high stability in the pH region of 5.0-7.0. The corresponding maximum stabilities, DeltaG(25), were obtained at pH 7.0 with values of 15.4 kcal mol(-1) (combined fluorescence and circular dichroism data), and 15.1 kcal mol(-1) (DSC), considering a DeltaC(p) of 1.72 +/- 0.24 kcal mol(-1) K(-1). The low histidine content ( approximately 1.7%) and the high acidic residue content ( approximately 22.5%) suggests a flat pH dependence of thermal stability in the region 2.0-8.8 and that the decrease in thermal stability at low pH can be due to the differences in pK values of the acidic groups.

Dynamics of protein and peptide hydration

Biological processes often involve the surfaces of proteins, where the structural and dynamic properties of the aqueous solvent are modified. Information about the dynamics of protein hydration can be obtained by measuring the magnetic relaxation dispersion (MRD) of the water (2)H and (17)O nuclei or by recording the nuclear Overhauser effect (NOE) between water and protein protons. Here, we use the MRD method to study the hydration of the cyclic peptide oxytocin and the globular protein BPTI in deeply supercooled solutions. The results provide a detailed characterization of water dynamics in the hydration layer at the surface of these biomolecules. More than 95% of the water molecules in contact with the biomolecular surface are found to be no more than two-fold motionally retarded as compared to bulk water. In contrast to small nonpolar molecules, the retardation factor for BPTI showed little or no temperature dependence, suggesting that the exposed nonpolar residues do not induce clathrate-like hydrophobic hydration structures. New NOE data for oxytocin and published NOE data for BPTI were analyzed, and a mutually consistent interpretation of MRD and NOE results was achieved with the aid of a new theory of intermolecular dipolar relaxation that accounts explicitly for the dynamic perturbation at the biomolecular surface. The analysis indicates that water-protein NOEs are dominated by long-range dipolar couplings to bulk water, unless the monitored protein proton is near a partly or fully buried hydration site where the water molecule has a long residence time.

Identification of a Residue Critical for Maintaining the Functional Conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.

DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution

The thermal properties of two forms of the Drosophila melanogaster HMG-D protein, with and without its highly basic 26 residue C-terminal tail (D100 and D74) and the thermodynamics of their non-sequence-specific interaction with linear DNA duplexes were studied using scanning and titration microcalorimetry, spectropolarimetry, fluorescence anisotropy and FRET techniques at different temperatures and salt concentrations. It was shown that the C-terminal tail of D100 is unfolded at all temperatures, whilst the state of the globular part depends on temperature in a rather complex way, being completely folded only at temperatures close to 0 degrees C and unfolding with significant heat absorption at temperatures below those of the gross denaturational changes. The association constant and thus Gibbs energy of binding for D100 is much greater than for D74 but the enthalpies of their association are similar and are large and positive, i.e. DNA binding is a completely entropy-driven process. The positive entropy of association is due to release of counterions and dehydration upon forming the protein/DNA complex. Ionic strength variation showed that electrostatic interactions play an important but not exclusive role in the DNA binding of the globular part of this non-sequence-specific protein, whilst binding of the positively charged C-terminal tail of D100 is almost completely electrostatic in origin. This interaction with the negative charges of the DNA phosphate groups significantly enhances the DNA bending. An important feature of the non-sequence-specific association of these HMG boxes with DNA is that the binding enthalpy is significantly more positive than for the sequence-specific association of the HMG box from Sox-5, despite the fact that these proteins bend the DNA duplex to a similar extent. This difference shows that the enthalpy of dehydration of apolar groups at the HMG-D/DNA interface is not fully compensated by the energy of van der Waals interactions between these groups, i.e. the packing density at the interface must be lower than for the sequence-specific Sox-5 HMG box.

Raman study of the thermal behaviour and conformational stability of basic pancreatic trypsin inhibitor

We have studied the thermal denaturation of native basic pancreatic trypsin inhibitor (BPTI) by monitoring the Raman bands in the 4000-400 cm(-1) range. In agreement with results obtained by calorimetry, a cooperative melting transition is observed starting at 75 degrees C. This transition is found to involve predominantly the unfolding of helical structures accompanied by beta-aggregation, loss of hydrophobic interactions between side chains and changes in CSSC dihedral angles. However, salt bridge breaking starts near 40 degrees C, as deduced from the nu(s)(COO(-)) band and from the bands close to 1320 and 1345 cm(-1) which for the first time have been shown to be due largely to vibrations of the arginine guanidyl group in BPTI. The thermal stability is, hence, attributable to cooperative contributions from hydrophobic and backbone hydrogen bond interactions as well as from disulfide bonds.

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.

Unfolding of a leucine zipper is not a simple two-state transition

Temperature-induced unfolding of the leucine zipper, an alpha-helical, double-stranded, coiled-coil, was studied by circular dichroism spectroscopy, spectrofluorimetry and heat capacity scanning calorimetry. It is shown that this process does not represent a simple two-state transition, as previously believed, but consists of several stages. The first transition starts at the very beginning of heating from 0 degrees C and proceeds with significant heat absorption and decrease of ellipticity. This transition does not depend on the concentration of protein and is sensitive to modification of the N terminus; it is therefore associated with unfolding or fraying of this part of the leucine zipper. The second transition takes place at a considerably higher temperature; it is more pronounced than the first one and does not depend on the concentration of protein, i.e. it is unimolecular. This transition is sensitive to modification of both termini of the leucine zipper and affects the optical properties of a tryptophan residue placed in the central part of the zipper. It therefore involves the whole dimer but does not result in its dissociation, presumably being associated with some repacking of the coiled-coil. This second transition is followed at higher temperatures by the concentration-dependent cooperative unfolding/dissociation of the two strands. The enthalpy and entropy of the temperature-induced structural changes of the leucine zipper that take place before its cooperative unfolding/dissociation comprises almost 40% of the total enthalpy and entropy of unfolding of the completely folded coiled-coil, the state in which it appears to be below 0 degrees C. Comparison of the total enthalpy of leucine zipper unfolding with that of a single-stranded alpha-helix shows that their temperature-dependence correlates with the extent of intramolecular non-polar contacts and allows an assessment of the enthalpy of hydrogen bonding in alpha-helices, which appears to be about 3.3kJmol(-1) at 20 degrees C.

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.

Selection of potent chymotrypsin and elastase inhibitors from M13 phage library of basic pancreatic trypsin inhibitor (BPTI)

The combinatorial approach offered by phage display has proved to be powerful in obtaining novel variants of canonical inhibitors of serine proteinases that show new binding patterns. We applied this strategy to search for variants of basic pancreatic trypsin inhibitor (BPTI) that would be strong inhibitors of two serine proteinases: bovine alpha-chymotrypsin and porcine pancreatic elastase. BPTI only moderately inhibits the first and does not inhibit the second enzyme. A representative library of 3.2 x 10(4) BPTI variants, randomized at P(1), P(1)', P(2)' and P(3)' positions of the proteinase binding loop, was displayed on the surface of phage M13. After four to five rounds of selection on the target proteinase consensus sequences of the inhibitor binding loop were obtained. In both cases, the variants selected differed from BPTI at two to four positions, with a strong preference for selection of hydrophobic residues. Nevertheless, five of these variants expressed in a free form appeared to be correctly folded, stable proteins, and did not aggregate during thermal denaturation. The midpoints of the thermal unfolding curves of these variants were lowered by 5-20 degrees C as compared to BPTI. The expressed variants proved to be new potent inhibitors of the target enzymes with association constants up to 6.9 x 10(9) M(-1) and 3.7 x 10(10) M(-1) for elastase and chymotrypsin, respectively. Thus, the inhibitory properties of BPTI were improved by as much as 7 x 10(6)-fold towards elastase and 420-fold towards chymotrypsin.

Are proteins well-packed?

The average packing density inside proteins is as high as in crystalline solids. Does this mean proteins are well-packed? We go beyond average densities, and look at the full distribution functions of free volumes inside proteins. Using a new and rigorous Delaunay triangulation method for parsing space into empty and filled regions, we introduce formal definitions of interior and surface packing densities. Although proteins look like organic crystals by the criterion of average density, they look more like liquids and glasses by the criterion of their free volume distributions. The distributions are broad, and the scalings of volume-to-surface, volume-to-cluster-radius, and numbers of void versus volume show that the interiors of proteins are more like randomly packed spheres near their percolation threshold than like jigsaw puzzles. We find that larger proteins are packed more loosely than smaller proteins. And we find that the enthalpies of folding (per amino acid) are independent of the packing density of a protein, indicating that van der Waals interactions are not a dominant component of the folding forces.

The effect of net charge on the solubility, activity, and stability of ribonuclease Sa

The net charge and isoelectric pH (pI) of a protein depend on the content of ionizable groups and their pK values. Ribonuclease Sa (RNase Sa) is an acidic protein with a pI = 3.5 that contains no Lys residues. By replacing Asp and Glu residues on the surface of RNase Sa with Lys residues, we have created a 3K variant (D1K, D17K, E41K) with a pI = 6.4 and a 5K variant (3K + D25K, E74K) with a pI = 10.2. We show that pI values estimated using pK values based on model compound data can be in error by >1 pH unit, and suggest how the estimation can be improved. For RNase Sa and the 3K and 5K variants, the solubility, activity, and stability have been measured as a function of pH. We find that the pH of minimum solubility varies with the pI of the protein, but that the pH of maximum activity and the pH of maximum stability do not.

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.

Inhibition of six serine proteinases of the human coagulation system by mutants of bovine pancreatic trypsin inhibitor

A series of 12 bovine pancreatic trypsin inhibitor variants mutated in the P(4) and P(3) positions of the canonical binding loop containing additional K15R and M52L mutations were used to probe the role of single amino acid substitutions on binding to bovine trypsin and to the following human proteinases involved in blood clotting: plasmin, plasma kallikrein, factors X(a) and XII(a), thrombin, and protein C. The mutants were expressed in Escherichia coli as fusion proteins with the LE1413 hydrophobic polypeptide and purified from inclusion bodies; these steps were followed by CNBr cleavage and oxidative refolding. The mutants inhibited the blood-clotting proteinases with association constants in the range of 10(3)-10(10) m(-)(1). Inhibition of plasma kallikrein, factors X(a) and XII(a), thrombin, and protein C could be improved by up to 2 orders of magnitude by the K15R substitution. The highest increase in the association constant for P(3) mutant was measured for factor XII(a); P13S substitution increased the K(a) value 58-fold. Several other substitutions at P(3) resulted in about 10-fold increase for factor X(a), thrombin, and protein C. The cumulative P(3) and P(1) effects on K(a) values for the strongest mutant compared with the wild type bovine pancreatic trypsin inhibitor were in the range of 2.2- (plasmin) to 4,000-fold (factors XII(a) and X(a)). The substitutions at the P(4) site always caused negative effects (a decrease in the range from over 1,000- to 1.3-fold) on binding to all studied enzymes, including trypsin. Thermal stability studies showed a very large decrease of the denaturation temperature (about 22 degrees C) for all P(4) mutants, suggesting that substitution of the wild type Gly-12 residue leads to a change in the binding loop conformation manifesting itself in non-optimal binding to the proteinase active site.

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.

The structure of denatured bovine pancreatic trypsin inhibitor (BPTI)

In the presence of denaturant and thiol initiator, the native bovine pancreatic trypsin inhibitor (BPTI) denatures by shuffling its native disulfide bonds and converts to a mixture of scrambled isomers. The extent of denaturation is evaluated by the relative yields of the scrambled and native species of BPTI. BPTI is an exceedingly stable molecule and can be effectively denatured only by guanidine thiocyanate (GdmSCN) at concentrations higher than 3-4 M. The denatured BPTI consists of at least eight fractions of scrambled isomers. Their composition varies under increasing concentrations of GdmSCN. In the presence of 6 M GdmSCN, the most predominant fraction of scrambled BPTI accounts for 56% of the total structure of denatured BPTI. Structural analysis reveals that this predominant fraction contains the bead-form isomer of scrambled BPTI, bridged by three pairs of neighboring cysteines, Cys5-Cys14, Cys30-Cys38 and Cys51-Cys55. The extreme conformational stability of BPTI has important implications in its distinctive folding pathway.

Pressure response of protein backbone structure. Pressure-induced amide 15N chemical shifts in BPTI

The effect of pressure on amide 15N chemical shifts was studied in uniformly 15N-labeled basic pancreatic trypsin inhibitor (BPTI) in 90%1H2O/10%2H2O, pH 4.6, by 1H-15N heteronuclear correlation spectroscopy between 1 and 2,000 bar. Most 15N signals were low field shifted linearly and reversibly with pressure (0.468 +/- 0.285 ppm/2 kbar), indicating that the entire polypeptide backbone structure is sensitive to pressure. A significant variation of shifts among different amide groups (0-1.5 ppm/2 kbar) indicates a heterogeneous response throughout within the three-dimensional structure of the protein. A tendency toward low field shifts is correlated with a decrease in hydrogen bond distance on the order of 0.03 A/2 kbar for the bond between the amide nitrogen atom and the oxygen atom of either carbonyl or water. The variation of 15N shifts is considered to reflect site-specific changes in phi, psi angles. For beta-sheet residues, a decrease in psi angles by 1-2 degrees/2 kbar is estimated. On average, shifts are larger for helical and loop regions (0.553 +/- 0.343 and 0.519 +/- 0.261 ppm/2 kbar, respectively) than for beta-sheet (0.295 +/- 0.195 ppm/2 kbar), suggesting that the pressure-induced structural changes (local compressibilities) are larger in helical and loop regions than in beta-sheet. Because compressibility is correlated with volume fluctuation, the result is taken to indicate that the volume fluctuation is larger in helical and loop regions than in beta-sheet. An important aspect of the volume fluctuation inferred from pressure shifts is that they include motions in slower time ranges (less than milliseconds) in which many biological processes may take place.

Primary structure and properties of the cathepsin G/chymotrypsin inhibitor from the larval hemolymph of Apis mellifera

A member of the Ascaris inhibitor family exhibiting anti-cathepsin G and anti-chymotrypsin activity was purified from the larval hemolymph of the honey bee (Apis mellifera). Three forms of the inhibitor, designated AMCI 1-3, were isolated using gel filtration and anion-exchange chromatographies followed by reverse-phase HPLC. The amino-acid analyses indicated that AMCI-1 and AMCI-2 have an identical composition whereas AMCI-3 is shorter by two residues (Thr, Arg). All three forms contain as many as 10 cysteine residues and lack tryptophan, tyrosine, and histidine. The sequence of the isoinhibitors showed that the major form (AMCI-1) consisting of 56 amino-acid residues was a single-chain protein of molecular mass 5972 Da, whereas the other two forms were two-chain proteins with a very high residue identity. The AMCI-2 appeared to be derived from AMCI-1, as a result of the Lys24-Thr25 peptide bond splitting, while AMCI-3 was truncated at its N-terminus by the dipeptide Thr25-Arg26. The association constants for the binding of bovine alpha-chymotrypsin to all purified forms of the inhibitor were high and nearly identical, ranging from 4.8 x 10(10) M-1 for AMCI-1 to 2.7 x 10(9) M-1 for AMCI-3. The sensitivity of cathepsin G to inhibition by each inhibitor was different. Only the association constant for the interaction of this enzyme with AMCI-1 was high (2 x 10(8) M-1) whereas those for AMCI-2 and AMCI-3 were significantly lower, and appeared to be 3.7 x 10(7) M-1 and 4.5 x 10(6) M-1, respectively. The reactive site of the inhibitor, as identified by cathepsin G degradation and chemical modification, was found to be at Met30-Gln31. A search in the Protein Sequence Swiss-Prot databank revealed a significant degree of identity (44%) between the primary structure of AMCI and the trypsin isoinhibitor from Ascaris sp (ATI). On the basis of the cysteine residues alignment, the position of the reactive site as well as some sequence homology, the cathepsin G/chymotrypsin inhibitor from larval hemolymph of the honey bee may be considered to be a member of the Ascaris inhibitor family.

Effect of glycerol on the interactions and solubility of bovine pancreatic trypsin inhibitor

The effects of additives used to stabilize protein structure during crystallization on protein solution phase behavior are poorly understood. Here we investigate the effect of glycerol and ionic strength on the solubility and strength of interactions of the bovine pancreatic trypsin inhibitor. These two variables are found to have opposite effects on the intermolecular forces; attractions increase with [NaCl], whereas repulsions increase with glycerol concentration. These changes are mirrored in bovine pancreatic trypsin inhibitor solubility where the typical salting out behavior for NaCl is observed with higher solubility found in buffers containing glycerol. The increased repulsions induced by glycerol can be explained by a number of possible mechanisms, all of which require small changes in the protein or the solvent in its immediate vicinity. Bovine pancreatic trypsin inhibitor follows the same general phase behavior as other globular macromolecules where a robust correlation between protein solution second virial coefficient and solubility has been developed. This study extends previous reports of this correlation to solution conditions involving nonelectrolyte additives.

Kinetic folding pathway of a three-disulfide mutant of bovine pancreatic ribonuclease A missing the [40-95] disulfide bond

The oxidative refolding pathway of a three-disulfide mutant of bovine pancreatic ribonuclease A (RNase A) from the fully reduced unfolded form to the native state has been studied by using oxidized and reduced dithiothreitol as the redox reagents at pH 8.0 and 25 degrees C. This mutant was prepared by replacing Cys40 and Cys95 in RNase A with alanines while maintaining the other three native disulfide bonds to mimic one of the two major three-disulfide intermediates (des-[40-95]) observed in the regeneration of wild-type RNase A. The kinetics of refolding of this mutant were measured by quenching the regeneration reaction at various times with a rapid blocking reagent, 2-aminoethyl methanethiosulfonate (AEMTS), fractionating the disulfide intermediates by using cation-exchange HPLC, and analyzing the time course of each group of disulfide species. It was found that the disulfide intermediates formed during regeneration reach a steady-state distribution after a short period of preequilibration similar to that in the regeneration of wild-type RNase A. The experimental data acquired under different redox conditions were fit to a kinetic model with a steady-state treatment. The fitted results indicate that this mutant refolds through a rate-determining step which involves the oxidation of certain two-disulfide species to form a putative three-disulfide species which proceeds rapidly to the native protein. A rough estimation suggests that this pathway could constitute no more than 5% of the major pathway leading to the formation of des-[40-95] (the major three-disulfide intermediate formed) in the regeneration of wild-type RNase A. Several kinetic constants pertaining to the oxidation and reduction of various disulfide intermediates were compared with those obtained in the regeneration studies of wild-type RNase A to gain further understanding about the folding pathways of RNase A. Comparisons are also given for the oxidative refolding studies of several other three disulfide bond proteins, suggesting that the formation of a large number of disulfide-bonded intermediates during oxidative refolding is probably a common feature for most proteins.