A genetic screen to identify variants of bovine pancreatic trypsin inhibitor with altered folding energetics.

A genetic screening procedure has been developed to identify mutant forms of bovine pancreatic trypsin inhibitor (BPTI) that can fold to an active conformation but are inactivated more rapidly than the wild-type protein. Small cultures of Escherichia coli containing plasmids with mutagenized BPTI genes were grown in microtiter plates, lysed, and treated with dithiothreitol (DTT). Under these conditions, unfolding and inactivation of the wild-type protein has a half-time of about 10 hours. Variants of BPTI that are inactivated within 1 hour were identified by adding trypsin and a chromogenic substrate. Approximately 11,000 mutagenized clones were screened in this way and 75 clones that produce proteins that can fold but are inactivated by DTT were isolated. The genes coding for 68 "DTT-sensitive" mutant proteins were sequenced, and 25 different single amino acid substitutions at 15 of the 58 residues of the protein were identified. Most of the altered residues are largely buried in the core of the native wild-type structure and are highly conserved among proteins homologous to BPTI. These results indicate that a large fraction of the sequence of the protein contributes to the kinetic stability of the active conformation, but it also appears that substitutions can be tolerated at most sites without completely preventing folding. Because this genetic screen is based on changes in folding energetics, further studies of the isolated mutants are expected to provide information about the roles of the altered residues in folding and unfolding.     

Structural stability of disulfide mutants of basic pancreatic trypsin inhibitor: A molecular dynamics study

The structure and folding of basic pancreatic trypsin inhibitor (BPTI) has been studied extensively by experimental means. We report a computer simulation study of the structural stability of various disulfide mutants of BPTI, involving eight 250-psec molecular dynamics simulations of the proteins in water, with and without a phosphate counterion. The presence of the latter alters the relative stability of the single disulfide species [5–55] and [30–51]. This conclusion can explain results of mutational studies and the conservation of residues in homologues of BPTI, and suggests a possible role of ions in stabilizing one intermediate over another in unfolding or folding processes. © 1996 Wiley-Liss, Inc.     

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI. © 1995 Wiley-Liss, Inc.     

Formation of a native-like subdomain in a partially folded intermediate of bovine pancreatic trypsin inhibitor.

In the folding of bovine pancreatic trypsin inhibitor (BPTI), the single-disulfide intermediate [30-51] plays a key role. We have investigated a recombinant analog of [30-51] using a 2-dimensional nuclear magnetic resonance (2D-NMR). This recombinant analog, named [30-51]Ala, contains a disulfide bond between Cys-30 and Cys-51, but contains alanine in place of the other cysteines in BPTI to prevent the formation of other intermediates. By 2D-NMR, [30-51]Ala consists of 2 regions-one folded and one predominantly unfolded. The folded region resembles a previously characterized peptide model of [30-51], named P alpha P beta, that contains a native-like subdomain with tertiary packing. The unfolded region includes the first 14 N-terminal residues of [30-51] and is as unfolded as an isolated peptide containing these residues. Using protein dissection, we demonstrate that the folded and unfolded regions of [30-51]Ala are structurally independent. The partially folded structure of [30-51]Ala explains many of the properties of authentic [30-51] in the folding pathway of BPTI. Moreover, direct structural characterization of [30-51]Ala has revealed that a crucial step in the folding pathway of BPTI coincides with the formation of a native-like subdomain, supporting models for protein folding that emphasize the formation of cooperatively folded subdomains.     

Comparison of solution structures of mutant bovine pancreatic trypsin inhibitor proteins using two-dimensional nuclear magnetic resonance.

Structural perturbations due to a series of mutations at the 30-51 disulfide bond of bovine pancreatic trypsin inhibitor have been explored using NMR. The mutants replaced cysteines at positions 30 and 51 by alanine at position 51 and alanine, threonine, or valine at position 30. Chemical shift changes occur in residues proximate to the site of mutation. NOE assignments were made using an automated procedure, NASIGN, which used information from the wild-type crystal structure. Intensity information was utilized by a distance geometry algorithm, VEMBED, to generate a series of structures for each protein. Statistical analyses of these structures indicated larger averaged structural perturbations than would be expected from crystallographic and other information. Constrained molecular dynamics refinement using AMBER at 900 K was useful in eliminating structural movements that were not a necessary consequence of the NMR data. In most cases, statistically significant movements are shown to be those greater than approximately 1 A. Such movements do not appear to occur between wild type and A30A51, a result confirmed by crystallography (Eigenbrot, C., Randal, M., & Kossiakoff, A.A., 1990, Protein Eng. 3, 591-598). Structural alterations in the T30A51 or V30A51 mutant proteins near the limits of detection occur in the beta-loop (residues 25-28) or C-terminal alpha-helix, respectively.     

Normal modes of vibration in bovine pancreatic trypsin inhibitor and its mechanical property

The normal mode analysis of conformational fluctuation is carried out for a small globular protein, bovine pancreatic trypsin inhibitor. Results are analyzed mainly to reveal the mechanical construction of the protein molecule. We take dihedral angles, including peptide omega angles, as independent variables for the normal mode analysis. There are 306 such angles in this molecule. Motions in modes with frequencies lower than 120 cm-1 are shown to involve atoms in the whole protein molecule, and spatial change of displacement vectors is continuous, i.e., those of atoms near in space are similar. To quantitate the observation of the continuity, a correlation function of direction vectors of atomic displacements is calculated. From this function we define a quantity that is interpreted as the wave length of an equivalent elastic plane wave. From this quantity we deduce effective Young's modulus for each mode. For the mode with the lowest frequency 4.4 cm-1, it turned out to be 0.8 x 10(9) dyn cm-2, the value two orders of magnitude softer than, for instance, alpha-helices. Prompted by this observation, the four lowest frequency modes and also the harmonic motions in the thermal equilibrium are analyzed further mainly to detect relatively rigid structural elements in the molecule. From this analysis emerges a mechanical picture of the protein molecule that is made up of relatively rigid elements held together by very soft parts.     

Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.     

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-Å resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.     

Binding of the bovine basic pancreatic trypsin inhibitor (Kunitz inhibitor) to human and bovine factor Xa. A thermodynamic study

The effect of pH and temperature on the apparent association equilibrium constant (Ka) for the binding of the bovine basic pancreatic trypsin inhibitor (BPTI, Kunitz inhibitor) to human and bovine factor Xa (Stuart-Prower factor; EC 3.4.21.6) has been investigated. Under all the experimental conditions, values of Ka for BPTI binding to human and bovine factor Xa are identical. On lowering the pH from 9.5 to 4.5, values of Ka (at 21.0 degrees C) for BPTI binding to human and bovine factor Xa decrease, thus reflecting the acidic pK shift of the His57 catalytic residue from 7.1, in the free enzyme, to 5.2, in the proteinase-inhibitor complex. At pH 8.0, values of the apparent thermodynamic parameters for BPTI binding to human and bovine factor Xa are: Ka = 2.1 x 10(5)M-1 (at 21.0 degrees C), delta G degree = -29.7 kJ/mol (at 21.0 degrees C), delta S degree = +161 entropy units (at 21.0 degrees C), and delta H degree = +17.6 kJ/mol (temperature-independent over the explored range, from 5.0 degrees C to 45.0 degrees C). Thermodynamics of BPTI binding to human and bovine factor Xa have been analysed in parallel with those of related serine (pro)enzyme/Kazal- and /Kunitz-type inhibitor systems. Considering the known molecular models, the observed binding behaviour of BPTI to human and bovine factor Xa was related to the inferred stereochemistry of the proteinase/inhibitor contact region.     

A protein contortionist: core mutations of GB1 that induce dimerization and domain swapping

Immunoglobulin-binding domain B1 of streptococcal protein G (GB1), a small (56 residues), stable, single-domain protein, is one of the most extensively used model systems in the area of protein folding and design. Recently, NMR and X-ray structures of a quintuple GB1 core mutant (L5V/A26F/F30V/Y33F/A34F) that showed an unexpected, intertwined tetrameric architecture were determined. Here, we report the NMR structure of another mutant, derived from the tetramer by reverting the single amino acid position F26 back to the wild-type sequence A26. The structure reveals a domain-swapped dimer that involves exchange of the second beta-hairpin. The resulting overall structure comprises an eight-stranded beta-sheet whose concave side is covered by two alpha helices. The dimer dissociates into a partially folded, monomeric species with a dissociation constant of 93(+/-10)microM.     

Structure of the complex formed by bovine trypsin and bovine pancreatic trypsin inhibitor. Crystal structure determination and stereochemistry of the contact region

The structure of the complex of bovine trypsin and bovine pancreatic trypsin inhibitor has been determined by crystal structure analysis at 2.8 Å resolution. The structure is closely similar to the model predicted from the structures of the components. The complex is a tetrahedral adduct with a covalent bond between the carbonyl carbon of Lys-15I of the inhibitor and the γ-oxygen of Ser-195 of the enzyme. The imidazole of His-57 is hydrogen-bonded to Asp-102 and the bound seryl γ-oxygen in accord with the histidine being charged. The negatively charged carbonyl oxygen of Lys-15I forms two hydrogen bonds with the amide nitrogens of Gly-193 and Ser-195. Protonation of the leaving group N-H of Ala-16I to form an acyl-complex requires a conformational change of the imidazole of His-57. The tetrahedral adduct is further stabilized by hydrogen bonds between groups at the leaving group side and inhibitor and enzyme, which would be weakened in the acyl-enzyme. The kinetic data of inhibitor-enzyme interaction are reconciled with the structural model, and relations between enzyme-inhibitor interaction and productive enzyme-substrate interaction are proposed.     

Insights into sucrose isomerization from crystal structures of the Pseudomonas mesoacidophila MX-45 sucrose isomerase,MutB

Three-dimensional structures of a sucrose isomerase from Pseudomonas mesoacidophila MX-45, forming mainly trehalulose have been solved to resolutions in the range 1.8–2.2 Å. Native and mutant complexes give, for the first time, a thorough insight into substrate binding and recognition, and product specificities of these enzymes. This study has pinpointed essential residues for binding the substrate sucrose, and hereby given detailed information on the interactions between the enzyme active site and glucosyl- and fructosyl moieties. Moreover, the 3-D structures revealed an aromatic clamp formed by two phenylalanines, which plays an essential role in recognition of the substrate and in controlling the reaction specificity.     

A highly destabilizing mutation, G37A, of the bovine pancreatic trypsin inhibitor retains the average native conformation but greatly increases local flexibility

A point mutation, G37A, on the surface of bovine pancreatic trypsin inhibitor (BPTI) destabilizes the protein by approximately 5 kcal/mol, which is very high for addition of one methyl group. In wild-type (WT) BPTI, Gly 37 HN is in an unusual NH-aromatic-NH network of interactions with the ring of Tyr 35 and the side chain HN of Asn 44. G37A was designed to disrupt this interaction, since the phi and psi backbone angles of G37 are not favorable for an amino acid containing a beta-carbon. Investigations of the structure and dynamics by NMR methods show that G37A retains the average WT structure. The NH-aromatic-NH interactions remain intact, as indicated by NOEs and the large upfield ring current shift (approximately 4 ppm) of A37 HN. The NMR structure, confirmed by molecular modeling calculations, requires phi and psi backbone angles that are highly destabilizing when alanine is in position 37. Although the average structure is essentially unchanged, the dynamics are altered dramatically. Many residues in the region of the mutation have increased flexibility, as probed by aromatic ring flip rates and native state hydrogen exchange. We conclude that a large fraction of the destabilization arises from maintaining A37 in a high-energy conformation. This suggests that disruption of the NH-aromatic-NH network is energetically very costly, and may involve other cooperatively linked interactions. The results illustrate the importance of the Gly-Gly sequence at positions 36 and 37 and the 37 HN-35 aromatic interaction to the stability, folding, and dynamics of the BPTI.     

Correlation between the amide proton exchange rates and the denaturation temperatures in globular proteins related to the basic pancreatic trypsin inhibitor.

Nuclear magnetic resonance was used to measure the hydrogen-deuterium exchange rates for individual interior amide protons in a group of small globular proteins related to the basic pancreatic trypsin inhibitor (BPTI). These proteins include two homologous proteins and seven chemical modifications of BPTI. It was previously shown that the spatial structure of BPTI is preserved in all these related proteins. The exchange rates for corresponding amide protons in the different proteins were found to vary by a factor of as much as 5 X 104. The proton exchange is correlated with the thermal stability of the proteins, i.e. the lower the denaturation temperature, the faster the NH exchange. Further evidence that the exchange of interior amide protons is promoted by global fluctuations of the protein structures comes from the observation that the order of the relative exchange rates for the individual protons is the same in all the different species. This is the third in a series of three papers on nuclear magnetic resonance studies of labile protons in BPTI-related proteins. A detailed interpretation of the data will be given in a forthcoming paper.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Exploring the function of Tyr83 in bacteriorhodopsin: features of the Y83F and Y83N mutants.

Tyrosine-83, a residue which is conserved in all halobacterial retinal proteins, is located at the extracellular side in helix C of bacteriorhodopsin. Structural studies indicate that its hydroxyl group is hydrogen bonded to Trp189 and possibly to Glu194, a residue which is part of the proton release complex (PRC) in bacteriorhodopsin. To elucidate the role of Tyr83 in proton transport, we studied the Y83F and Y83N mutants. The Y83F mutation causes an 11 nm blue shift of the absorption spectrum and decreases the size of the absorption changes seen upon dark adaptation. The light-induced fast proton release, which accompanies formation of the M intermediate, is observed only at pH above 7 in Y83F. The pK(a) of the PRC in M is elevated in Y83F to about 7.3 (compared to 5.8 in WT). The rate of the recovery of the initial state (the rate of the O --> BR transition) and light-induced proton release at pH below 7 is very slow in Y83F (ca. 30 ms at pH 6). The amount of the O intermediate is decreased in Y83F despite the longer lifetime of O. The Y83N mutant shows a similar phenotype in respect to proton release. As in Y83F, the recovery of the initial state is slowed several fold in Y83N. The O intermediate is not seen in this mutant. The data indicate that the PRC is functional in Y83F and Y83N but its pK(a) in M is increased by about 1.5 pK units compared to the WT. This suggests that Tyr83 is not the main source for the proton released upon M formation in the WT; however, Tyr83 is involved in the proton release affecting the pK(a) of the PRC in M and the rate of proton transport from Asp85 to PRC during the O --> bR transition. Both the Y83F and the Y83N mutations lead to a greatly decreased functionality of the pigment at high pH because most of the pigment is converted into the inactive P480 species, with a pK(a) 8-9.     

Folding in vitro of bovine pancreatic trypsin inhibitor in the presence of proteins of the endoplasmic reticulum.

The rates of folding and disulfide bond formation in reduced BPTI were measured in vitro in the presence and absence of total protein from the endoplasmic reticulum. The rates were increased substantially by the endoplasmic reticulum proteins, but only to the extent expected from the known content and activity of protein-disulfide-isomerase. No effects of added ATP or Ca2+ were observed, even though protein-disulfide-isomerase binds Ca2+ tightly.     

Kinetic roles and conformational properties of the non-native two-disulphide intermediates in the refolding of bovine pancreatic trypsin inhibitor.

The most productive folding pathway of reduced bovine pancreatic trypsin inhibitor (BPTI) proceeds through the disulphide intermediates (30-51), (30-51, 5-14), and (30-51, 5-38); these are important kinetic intermediates in folding, even though the latter pair contain non-native disulphide bonds. Analogues of these intermediates have been prepared by protein engineering methods and their conformational properties examined by circular dichroism and 1H-nuclear magnetic resonance. The (30-51), (30-51, 5-14) and (30-51, 5-38) analogues exhibit comparable degrees of stable structure, which cannot include those portions of the polypeptide chain involving Cys5, Cys14 and Cys38. These properties are consistent with the roles of (30-51, 5-14) and (30-51, 5-38) in the folding pathway of BPTI, which demand that they exhibit a considerable degree of conformational flexibility in part of the molecule.