Accelerated exchange of a buried water molecule in selectively disulfide-reduced bovine pancreatic trypsin inhibitor

Using magnetic relaxation dispersion (MRD), we have previously shown that the four internal water molecules in bovine pancreatic trypsin inhibitor (BPTI) exchange with bulk water on time scales between 10(-8) and 10(-4) s at room temperature. Because this exchange is controlled by the protein structure, internal water molecules can be used to probe rare conformational fluctuations. Here, we report (2)H and (17)O MRD data at three temperatures for wild-type BPTI and two BPTI variants where the 14-38 disulfide bond has been cleaved by a double Cys --> Ser mutation or by disulfide reduction and carboxamidomethylation. The MRD data show that the internal water molecules are conserved on disulfide cleavage. However, the exchange rate of the water molecule buried near the disulfide bond is enhanced by 2-4 orders of magnitude. The relation of water exchange to other dynamic processes in BPTI is discussed.     

Pressure-dependent changes in the solution structure of hen egg-white lysozyme

The "rules" governing protein structure and stability are still poorly understood. Important clues have come from proteins that operate under extreme conditions, because these clarify the physical constraints on proteins. One obvious extreme is pressure, but so far little is known of the behavior of proteins under pressure, largely for technical reasons. We have therefore developed new methodology for calculating structure change in solution with pressure, using NMR chemical shift changes, and we report the change in structure of lysozyme on going from 30 bar to 2000 bar, this being the first solution structure of a globular protein under pressure. The alpha-helical domain is compressed by approximately 1%, due to tighter packing between helices. The interdomain region is also compressed. By contrast, the beta-sheet domain displays very little overall compression, but undergoes more structural distortion than the alpha-domain. The largest volume changes tend to occur close to hydrated cavities. Because isothermal compressibility is related to volume fluctuation, this suggests that buried water molecules play an important role in conformational fluctuation at normal pressures, and are implicated as the nucleation sites for structural changes leading to pressure denaturation or channel opening.     

High-resolution structure of bovine pancreatic trypsin inhibitor with altered binding loop sequence

A mutant of bovine pancreatic trypsin inhibitor (BPTI) has been constructed and expressed in Escherichia coli in order to probe the kinetic and structural consequences of truncating the binding loop residues to alanine. In addition to two such mutations (Thr11Ala and Pro13Ala), it has a conservative Lys15Arg substitution at position P(1) and an unrelated Met52Leu change. In spite of the binding loop modification, the affinity for trypsin is only 30 times lower than that of the wild-type protein. At pH 7.5 the protein can be crystallized on the time-scale of hours, yielding very stable crystals of a new (tetragonal) form of BPTI. Conventional source X-ray data collected to 1.4 A at room temperature allowed anisotropic structure refinement characterized by R=0.1048. The structure reveals all 58 residues, including the complete C terminus, which is in a salt-bridge contact with the N terminus. The Cys14-Cys38 disulfide bridge is observed in two distinct chiralities. This bridge, together with an internal water molecule, contributes to the stabilization of the binding loop. The Ala mutations have only an insignificant and localized effect on the binding loop, which retains its wild-type conformation (maximum deviation of loop C(alpha) atoms of 0.7 A at Ala13). Four (instead of the typical three) additional water molecules are buried in an internal cleft and connected to the surface via a sulfate anion. Three more SO(4)(2-) anions are seen in the electron density, one of them located on a 2-fold axis. It participates in the formation of a dimeric structure between symmetry-related BPTI molecules, in which electrostatic and hydrogen bonding interactions resulting from the mutated Lys15Arg substitution are of central importance. This dimeric interaction involves direct recognition loop-recognition loop contacts, part of which are hydrophobic interactions of the patches created by the alanine mutations. Another 2-fold symmetric interaction between the BPTI molecules involves the formation of an antiparallel intermolecular beta-sheet that, together with the adjacent intramolecular beta-hairpin loops, creates a four-stranded structure.     

An Effective Solvation Term Based on Atomic Occupancies for Use in Protein Simulations

Abstract

Several approaches to the treatment of solvent effects based on continuum models are reviewed and a new method based on occupied atomic volumes (occupancies) is proposed and tested. The new method describes protein-water interactions in terms of atomic solvation parameters, which represent the solvation free energy per unit of volume. These parameters were determined for six different atoms types, using experimental free energies of solvation. The method was implemented in the GROMOS and PRESTO molecular simulation program suites. Simulations with the solvation term require 20-50% more CPU time than the corresponding vacuum simulations and are approximately 20 times faster than explicit water simulations. The method and parameters were tested by carrying out 200 ps simulations of BPTI in water, in vacuo, and with the solvation term. The performance of the solvation term was assessed by comparing the structures and energies from the solvation simulations with the equivalent quantities derived from several BPTI crystal structures and from the explicit water and vacuum simulations. The model structures were evaluated in terms of exposed total surface, buried and exposed polar surfaces, secondary structure preservation, number of hydrogen bonds, energy contributions, and positional deviations from BPTI crystal structures. Vacuum simulations produced unrealistic structures with respect to all criteria applied. The structures resulting from the simulations with explicit water were closer to the 5PTI crystal structure, although part of the secondary structure dissolved. The simulations with the effective solvation term produce structures that are normal according to all evaluations and in most respects are remarkably similar to the 5PTI crystal structure despite considerable positional fluctuations during the simulations. The segments where the model and crystal structures differ are known to be flexible and the observed difference may be physically realistic. The effective solvation term based on occupancies is not only very efficient in terms of computer time but also results in meaningful structural properties for BPTI. It may therefore be generally useful in molecular dynamics of macromolecules.     

Molecular dynamics simulation of solvated protein at high pressure.

We have completed a molecular dynamics simulation of protein (bovine pancreatic trypsin inhibitor, BPTI) in solution at high pressure (10 kbar). The structural and energetic effects of the application of high pressure to solvated protein are analyzed by comparing the results of the high-pressure simulation with a corresponding simulation at low pressure. The volume of the simulation cell containing one protein molecule plus 2943 water molecules decreases by 24.7% at high pressure. This corresponds to a compressibility for the protein solution of beta = 1.8 x 10(-2) kbar-1. The compressibility of the protein is estimated to be about one-tenth that of bulk water, while the protein hydration layer water is found to have a greater compressibility as compared to the bulk, especially for water associated with hydrophobic groups. The radius of gyration of BPTI decreases by 2% and there is a one third decrease in the protein backbone atomic fluctuations at high pressure. We have analyzed pressure effects on the hydration energy of the protein. The total hydration energy is slightly (4%) more favorable at high pressure even though the surface accessibility of the protein has decreased by a corresponding amount. Large pressure-induced changes in the structure of the hydration shell are observed. Overall, the solvation shell waters appear more ordered at high pressure; the pressure-induced ordering is greatest for nonpolar surface groups. We do not observe evidence of pressure-induced unfolding of the protein over the 100-ps duration of the high-pressure simulation. This is consistent with the results of high-pressure optical experiments on BPTI.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hydrogen exchange in BPTI variants that do not share a common disulfide bond.

Bovine pancreatic trypsin inhibitor (BPTI) is stabilized by 3 disulfide bonds, between cysteines 30-51, 5-55, and 14-38. To better understand the influence of disulfide bonds on local protein structure and dynamics, we have measured amide proton exchange rates in 2 folded variants of BPTI, [5-55]Ala and [30-51; 14-38]V5A55, which share no common disulfide bonds. These proteins resemble disulfide-bonded intermediates that accumulate in the BPTI folding pathway. Essentially the same amide hydrogens are protected from exchange in both of the BPTI variants studied here as in native BPTI, demonstrating that the variants adopt fully folded, native-like structures in solution. However, the most highly protected amide protons in each variant differ, and are contained within the sequences of previously studied peptide models of related BPTI folding intermediates containing either the 5-55 or the 30-51 disulfide bond.     

Cluster analysis of consensus water sites in thrombin and trypsin shows conservation between serine proteases and contributions to ligand specificity.

Cluster analysis is presented as a technique for analyzing the conservation and chemistry of water sites from independent protein structures, and applied to thrombin, trypsin, and bovine pancreatic trypsin inhibitor (BPTI) to locate shared water sites, as well as those contributing to specificity. When several protein structures are superimposed, complete linkage cluster analysis provides an objective technique for resolving the continuum of overlaps between water sites into a set of maximally dense microclusters of overlapping water molecules, and also avoids reliance on any one structure as a reference. Water sites were clustered for ten superimposed thrombin structures, three trypsin structures, and four BPTI structures. For thrombin, 19% of the 708 microclusters, representing unique water sites, contained water molecules from at least half of the structures, and 4% contained waters from all 10. For trypsin, 77% of the 106 microclusters contained water sites from at least half of the structures, and 57% contained waters from all three. Water site conservation correlated with several environmental features: highly conserved microclusters generally had more protein atom neighbors, were in a more hydrophilic environment, made more hydrogen bonds to the protein, and were less mobile. There were significant overlaps between thrombin and trypsin conserved water sites, which did not localize to their similar active sites, but were concentrated in buried regions including the solvent channel surrounding the Na+ site in thrombin, which is associated with ligand selectivity. Cluster analysis also identified water sites conserved in thrombin but not trypsin, and vice versa, providing a list of water sites that may contribute to ligand discrimination. Thus, in addition to facilitating the analysis of water sites from multiple structures, cluster analysis provides a useful tool for distinguishing between conserved features within a protein family and those conferring specificity.     

Secretion incompetence of bovine pancreatic trypsin inhibitor expressed in Escherichia coli.

There is increasing evidence that protein folding and protein export are competing processes in prokaryotic cells. Virtually all secretion studies reported to date, however, have employed proteins that are relatively uncharacterized in terms of their folding behavior and three-dimensional structure. In contrast, the structural and biochemical parameters governing the folding of bovine pancreatic trypsin inhibitor (BPTI) and several of its mutants have been studied intensively. We therefore undertook a study of the secretion behavior in Escherichia coli of recombinant BPTI and its mutants. Wild-type BPTI and two well-characterized folding mutants (C14A, C38A)BPTI and (C30A, C51A)BPTI (missing the 14-38 and 30-51 disulfide bonds, respectively), were investigated by analyzing their expression fused to an E. coli signal sequence or to two synthetic IgG-binding domains of staphylococcal protein A. Both disulfide mutants are destabilized relative to wild-type BPTI and exhibit markedly altered folding kinetics: one (C14A, C38A) folds more slowly than wild-type BPTI and the other (C30A, C51A) unfolds more rapidly. Both mutants were observed to be exported 3-10 times more efficiently than the wild-type molecule. Moreover, the levels of unprocessed preprotein in the cytoplasm were severalfold higher for the wild-type fusion than for the fusion to the two folding mutants. Intracellular degradation of the BPTI moiety was also observed. These results are consistent with traffic of intracellular BPTI preproteins on at least three routes along the secretory pathway: (a) facile secretion of unfolded material, (b) intracellular folding leading to secretion blockage, and (c) degradation followed by export of truncated molecules. A novel feature of these findings is the implication that disulfide bonds can form in the bacterial cytoplasm and lead to secretion incompetence.     

Production of native, correctly folded bovine pancreatic trypsin inhibitor by Escherichia coli

A gene for bovine pancreatic trypsin inhibitor (BPTI) was fused to the coding sequence for the Escherichia coli alkaline phosphatase signal peptide and expressed in E. coli under the control of the alkaline phosphatase promoter. When induced in phosphate-depleted medium such cells produced a trypsin inhibitor that was indistinguishable from native, properly folded BPTI. In particular, the BPTI produced by E. coli had three disulfide bonds that appeared to be identical to those found in native BPTI, as assayed by sensitivity to iodoacetate, dithiothreitol, and urea. This expression/secretion system will make possible the production of variant BPTI molecules, thus allowing the perturbing effects of amino acid substitutions on BPTI folding, structure, and function to be assessed.     

Structure of tick anticoagulant peptide at 1.6 A resolution complexed with bovine pancreatic trypsin inhibitor

The structure of tick anticoagulant peptide (TAP) has been determined by X-ray crystallography at 1.6 A resolution complexed with bovine pancreatic trypsin inhibitor (BPTI). The TAP-BPTI crystals are tetragonal, a = b = 46.87, c = 50.35 A, space group P41, four complexes per unit cell. The TAP molecules are highly dipolar and form an intermolecular helical array along the c-axis with a diameter of about 45 A. Individual TAP units interact in a head-to-tail fashion, the positive end of one molecule associating with the distal negative end of another, and vice versa. The BPTI molecules have a uniformly distributed positively charged surface that interacts extensively through 14 hydrogen bonds and two hydrogen bonded salt bridges with the helical groove around the helical TAP chains. Comparing the structure of TAP in TAP-BPTI with TAP bound to factor Xa(Xa) suggests a massive reorganization in the N-terminal tetrapeptide and the first disulfide loop of TAP (Cys5T-Cys15T) upon binding to Xa. The Tyr1(T)OH atom of TAP moves 14.2 A to interact with Asp189 of the S1 specificity site, Arg3(T)CZ moves 5.0 A with the guanidinium group forming a cation-pi-electron complex in the S4 subsite of Xa, while Lys7(T)NZ differs in position by 10.6 A in TAP-BPTI and TAP-Xa, all of which indicates a different pre-Xa-bound conformation for the N-terminal of TAP in its native state. In contrast to TAP, the BPTI structure of TAP-BPTI is practically the same as all those of previously determined structures of BPTI, only arginine and lysine side-chain conformations showing significant differences.     

The heme redox center of chloroplast cytochrome f is linked to a buried five-water chain.

The crystal structure of the 252-residue lumen-side domain of reduced cytochrome f, a subunit of the proton-pumping integral cytochrome b6f complex of oxygenic photosynthetic membranes, was determined to a resolution of 1.96 A from crystals cooled to -35 degrees. The model was refined to an R-factor of 15.8% with a 0.013-A RMS deviation of bond lengths from ideality. Compared to the structure of cytochrome f at 20 degrees, the structure at -35 degrees has a small change in relative orientation of the two folding domains and significantly lower isotropic temperature factors for protein atoms. The structure revealed an L-shaped array of five buried water molecules that extend in two directions from the N delta 1 of the heme ligand His 25. The longer branch extends 11 A within the large domain, toward Lys 66 in the prominent basic patch at the top of the large domain, which has been implicated in the interaction with the electron acceptor, plastocyanin. The water sites are highly occupied, and their temperature factors are comparable to those of protein atoms. Virtually all residues that form hydrogen bonds with the water chain are invariant among 13 known cytochrome f sequences. The water chain has many features that optimize it as a proton wire, including insulation from the protein medium. It is suggested that this chain may function as the lumen-side exit port for proton translocation by the cytochrome b6f complex.     

Thermodynamics of substrate binding to the chaperone SecB

The thermodynamics of binding of unfolded polypeptides to the chaperone SecB was investigated in vitro by isothermal titration calorimetry and fluorescence spectroscopy. The substrates were reduced and carboxamidomethylated forms of RNase A, BPTI, and alpha-lactalbumin. SecB binds both fully unfolded RNase A and BPTI as well as compact, partially folded disulfide intermediates of alpha-lactalbumin, which have 40-60% of native secondary structure. The heat capacity changes observed on binding the reduced and carboxamidomethylated forms of alpha-lactalbumin, BPTI, and RNase A were found to be -0.10, -0.29, and -0.41 kcal mol(-1) K(-1), respectively, and suggest that between 7 and 29 residues are buried upon substrate binding to SecB. In all cases, binding occurs with a stoichiometry of one polypeptide chain per monomer of SecB. There is no evidence for two separate types of binding sites for positively charged and hydrophobic ligands. Spectroscopic and proteolysis protection studies of the binding of SecB to poly-L-Lys show that binding of highly positively charged peptide ligands to negatively charged SecB leads to charge neutralization and subsequent aggregation of SecB. The data are consistent with a model where SecB binds substrate molecules at an exposed hydrophobic cleft. SecB aggregation in the absence of substrate is prevented by electrostatic repulsion between negatively charged SecB tetramers.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

Pressure-Dependent Structure Changes in Barnase on Ligand Binding Reveal Intermediate Rate Fluctuations

In this work we measured ¹H NMR chemical shifts for the ribonuclease barnase at pressures from 3 MPa to 200 MPa, both free and bound to d(CGAC). Shift changes with pressure were used as restraints to determine the change in structure with pressure. Free barnase is compressed by ∼0.7%. The largest changes are on the ligand-binding face close to Lys-27, which is the recognition site for the cleaved phosphate bond. This part of the protein also contains the buried water molecules. In the presence of d(CGAC), the compressibility is reduced by ∼70% and the region of structural change is altered: the ligand-binding face is now almost incompressible, whereas changes occur at the opposite face. Because compressibility is proportional to mean square volume fluctuation, we conclude that in free barnase, volume fluctuation is largest close to the active site, but when the inhibitor is bound, the fluctuations become much smaller and are located mainly on the opposite face. The timescale of the fluctuations is nanoseconds to microseconds, consistent with the degree of ordering required for the fluctuations, which are intermediate between rapid uncorrelated side-chain dynamics and slow conformational transitions. The high-pressure technique is therefore useful for characterizing motions on this relatively inaccessible timescale.     

Hydrogen exchange of primary amide protons in basic pancreatic trypsin inhibitor: evidence for NH2 group rotation in buried asparagine side chains

Hydrogen-deuterium exchange is measured for the buried primary amide groups of Asn-43 and Asn-44 in bovine pancreatic trypsin inhibitor. Amide protons trans and cis to the amide carbonyl oxygen (HE and HZ, respectively) exchange at indistinguishable rates. Uncorrelated exchange of HE and HZ is established for both residues by following the nuclear Overhauser enhancement from HE to HZ during the deuterium exchange. The exchange of Asn-43 and Asn-44 side-chain protons differs qualitatively from exchange of primary amide groups in fully solvated model compounds, for which HE generally exchanges faster than HZ. The equal rates for the buried primary amide HE and HZ in BPTI are not a consequence of coupled exchange. The data indicate rapid rotation around the CO-NH2 bond for both Asn-43 and Asn-44 and suggest considerable lability of intramolecular hydrogen bonds. The side chain of Asn-43 has all of its polar atoms integrated into the very stable hydrogen-bonded structure of the protein. Asn-44 is hydrogen-bonded to side chains and to a buried water molecule. Solvent isotope exchange is several orders of magnitude more restricted by protein secondary and tertiary structure than the CO-NH2 rotation, indicating that N delta H2 groups flip many times before hydrogen isotope exchange occurs.     

Correlation between calculated local stability and hydrogen exchange rates in proteins

The attempt is made to find new correlations between local structural characteristics of proteins and the hydrogen exchange rates of their individual main-chain amides, and to relate such correlations to possible mechanisms of hydrogen exchange. It is found that in bovine pancreatic trypsin inhibitor (BPTI) the surface area buried by a particular residue and its neighbors correlates with the exchange rate of the main-chain amide of that residue. As the area buried by a particular fragment can be associated with the stabilization of the protein structure by this fragment, the correlation suggests a role for the energetics of the local unfolding in the mechanism of hydrogen exchange. Calculations based on the assumption that the exchange mechanism involves local unfolding lead to quantitative agreement between the calculated and experimentally measured exchange rates for 80% of the amides of BPTI that are buried or hydrogen bonded to the main-chain or to internal water molecules. The same degree of correlation is found between the calculated exchange rates and partial exchange data for ribonuclease S, hen lysozyme and cytochrome c. A similarly strong correlation is found between calculated exchange rates and the exchange rates of ribonuclease A determined by neutron diffraction in the crystal. The criteria of correlation are, however, less stringent in this case because of the experimental errors, which are larger than for solution data. It is suggested that the observed correlation be used for predictions of hydrogen exchange rates in proteins.     

Prediction of a 12-residue loop in bovine pancreatic trypsin inhibitor: effects of buried water

The x-ray conformations of 5-, 7-, 9-, and 12-residue loops in bovine pancreatic trypsin inhibitor (BPTI) were predicted by the use of multiple independent Monte Carlo simulating annealing (MCSA) runs starting from random conformations. Four buried water molecules interacted with a 12-residue loop that started at residue 8 and ended at residue 19, and that included the binding region. The final conformation at the end of an MCSA run was characterized. Solvation free energy based on the solvent accessible surface area was included in the energy function at low simulated annealing temperatures. Conformational states were interactively separated by a recently developed algorithm. Computed loops were characterized in terms of total energy, and backbone and side chain root mean square deviations (RMSDs) between computed native loop conformations and the x-ray conformation. The 12-residue loop was computed with and without buried water [called WL12(8-19) and L12(8-19), respectively]. The backbone was reliably and reproducibly computed to within 1.1 A in L12(8-19) and 0.9 A in WL12(8-19). L12(8-19) required significantly more MCSA runs to achieve the same level of reproducibility as WL12(8-19). Based on the size of the cluster of low energy native loop conformations, and the computational effort, WL12(8-19) had greater entropy. In calculations of 7-, 9-, and 12-residue loops without buried water, the effects of buried water became obvious in the 12-residue loop calculation, which interacted with all four buried water molecules. Nearly all conformations of the native loop conformer had a hydrogen bond between the Lys 15 side chain and the backbone of Gly 12, Pro 13, and Cys 14, which may have implications in the rate of exchange of buried water with bulk solvent and in protein folding. The present version of MCSA program was more efficient than earlier versions.     

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 of BPTI deduced from the analysis of the intermediate scattering functions show a slowing down of protein dynamics when increasing pressure.     

Binding of bovine pancreatic trypsin inhibitor to trypsinogen: spectroscopic and volumetric studies

We have investigated the binding of bovine pancreatic trypsin inhibitor (BPTI) to bovine trypsinogen by combining ultrasonic velocimetry, high precision densimetry, and fluorescence spectroscopy. We report the changes in volume, adiabatic compressibility, van't Hoff enthalpy, entropy, and free energy that accompany the association of the two proteins at 25 degrees C and pH 8.0. We have used the measured changes in volume and compressibility in conjunction with available structural data to characterize the binding-induced changes in the hydration properties and intrinsic packing of the two proteins. Our estimate reveals that 110 +/- 40 water molecules become released to the bulk from the hydration shells of BPTI and trypsinogen. Furthermore, we find that the intrinsic coefficient of adiabatic compressibility of the two proteins decreases by 14 +/- 2%, which is suggestive of the binding-induced rigidification of the proteins' interior. BPTI-trypsinogen association is an entropy-driven event which proceeds with an unfavorable change in enthalpy. The favorable change in entropy results from partial compensation between two predominant terms. Namely, a large favorable change in hydrational entropy slightly prevails over a close in magnitude but opposite in sign change in configurational entropy. The reduction in configurational entropy and, consequently, protein dynamics is consistent with the observed decrease in intrinsic compressibility. In general, results of this work emphasize the vital role that water plays in modulating protein recognition events.     

Free‐energy function based on an all‐atom model for proteins

We have developed a free‐energy function based on an all‐atom model for proteins. It comprises two components, the hydration entropy (HE) and the total dehydration penalty (TDP). Upon a transition to a more compact structure, the number of accessible configurations arising from the translational displacement of water molecules in the system increases, leading to a water‐entropy gain. To fully account for this effect, the HE is calculated using a statistical‐mechanical theory applied to a molecular model for water. The TDP corresponds to the sum of the hydration energy and the protein intramolecular energy when a fully extended structure, which possesses the maximum number of hydrogen bonds with water molecules and no intramolecular hydrogen bonds, is chosen as the standard one. When a donor and an acceptor (e.g., N and O, respectively) are buried in the interior after the break of hydrogen bonds with water molecules, if they form an intramolecular hydrogen bond, no penalty is imposed. When a donor or an acceptor is buried with no intramolecular hydrogen bond formed, an energetic penalty is imposed. We examine all the donors and acceptors for backbone‐backbone, backbone‐side chain, and side chain‐side chain intramolecular hydrogen bonds and calculate the TDP. Our free‐energy function has been tested for three different decoy sets. It is better than any other physics‐based or knowledge‐based potential function in terms of the accuracy in discriminating the native fold from misfolded decoys and the achievement of high Z‐scores. Proteins 2009. © 2009 Wiley‐Liss, Inc.