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.     

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.     

Response of dynamic structure to removal of a disulfide bond: normal mode refinement of C77A/C95A mutant of human lysozyme.

In order to investigate the response of dynamic structure to removal of a disulfide bond, the dynamic structure of human lysozyme has been compared to its C77A/C95A mutant. The dynamic structures of the wild type and mutant are determined by normal mode refinement of 1.5-A-resolution X-ray data. The C77A/C95A mutant shows an increase in apparent fluctuations at most residues. However, most of the change originates from an increase in the external fluctuations, reflecting the effect of the mutation on the quality of crystals. The effects of disulfide bond removal on the internal fluctuations are almost exclusively limited to the mutation site at residue 77. No significant change in the correlation of the internal fluctuations is found in either the overall or local dynamics. This indicates that the disulfide bond does not have any substantial role to play in the dynamic structure. A comparison of the wild-type and mutant coordinates suggests that the disulfide bond does not prevent the 2 domains from parting from each other. Instead, the structural changes are characteristic of a cavity-creating mutation, where atoms surrounding the mutation site move cooperatively toward the space created by the smaller alanine side chain. Although this produces tighter packing, more than half of the cavity volume remains unoccupied, thus destabilizing the native state.     

Changing the determinants of protein stability from covalent to non-covalent interactions by in vitro evolution: a structural and energetic analysis

The three disulfide bonds of the gene-3-protein of the phage fd are essential for the conformational stability of this protein, and it unfolds when they are removed by reduction or mutation. Previously, we used an iterative in vitro selection strategy to generate a stable and functional form of the gene-3-protein without these disulfides. It yielded optimal replacements for the disulfide bonds as well as several stabilizing second-site mutations. The best selected variant showed a higher thermal stability compared with the disulfide-bonded wild-type protein. Here, we investigated the molecular basis of this strong stabilization by solving the crystal structure of this variant and by analyzing the contributions to the conformational stability of the selected mutations individually. They could mostly be explained by improved side-chain packing. The R29W substitution alone increased the midpoint of the thermal unfolding transition by 14 deg and the conformational stability by about 25 kJ mol^− 1. This key mutation (i) removed a charged side chain that forms a buried salt bridge in the disulfide-containing wild-type protein, (ii) optimized the local packing with the residues that replace the C46-C53 disulfide and (iii) improved the domain interactions. Apparently, certain residues in proteins indeed play key roles for stability.     

The Humicola grisea Cel12A enzyme structure at 1.2 A resolution and the impact of its free cysteine residues on thermal stability

As part of a program to discover improved glycoside hydrolase family 12 (GH 12) endoglucanases, we have extended our previous work on the structural and biochemical diversity of GH 12 homologs to include the most stable fungal GH 12 found, Humicola grisea Cel12A. The H. grisea enzyme was much more stable to irreversible thermal denaturation than the Trichoderma reesei enzyme. It had an apparent denaturation midpoint (T(m)) of 68.7 degrees C, 14.3 degrees C higher than the T. reesei enzyme. There are an additional three cysteines found in the H. grisea Cel12A enzyme. To determine their importance for thermal stability, we constructed three H. grisea Cel12A single mutants in which these cysteines were exchanged with the corresponding residues in the T. reesei enzyme. We also introduced these cysteine residues into the T. reesei enzyme. The thermal stability of these variants was determined. Substitutions at any of the three positions affected stability, with the largest effect seen in H. grisea C206P, which has a T(m) 9.1 degrees C lower than that of the wild type. The T. reesei cysteine variant that gave the largest increase in stability, with a T(m) 3.9 degrees C higher than wild type, was the P201C mutation, the converse of the destabilizing C206P mutation in H. grisea. To help rationalize the results, we have determined the crystal structure of the H. grisea enzyme and of the most stable T. reesei cysteine variant, P201C. The three cysteines in H. grisea Cel12A play an important role in the thermal stability of this protein, although they are not involved in a disulfide bond.     

Crystal structures of two mutants that have implications for the folding of bovine pancreatic ribonuclease A.

The Tyr92-Pro93 peptide group of bovine pancreatic ribonuclease A (RNase A) exists in the cis conformation in the native state. From unfolding/refolding kinetic studies of the disulfide-intact wild-type protein and of a variant in which Pro93 had been replaced by Ala, it had been suggested that the Tyr92-Ala93 peptide group also exists in the cis conformation in the native state. Here, we report the crystal structure of the P93A variant. Although there is disorder in the region of residues 92 and 93, the best structural model contains a cis peptide at this position, lending support to the results of the kinetics experiments. We also report the crystal structure of the C[40, 95]A variant, which is an analog of the major rate-determining three-disulfide intermediate in the oxidative folding of RNase A, missing the 40-95 disulfide bond. As had been detected by NMR spectroscopy, the crystal structure of this analog shows disorder in the region surrounding the missing disulfide. However, the global chain fold of the remainder of the protein, including the disulfide bond between Cys65 and Cys72, appears to be unaffected by the mutation.     

Variability of conformations at crystal contacts in BPTI represent true low-energy structures: correspondence among lattice packing and molecular dynamics structures.

The structures of five basic pancreatic trypsin inhibitor (BPTI) molecules are compared to establish the extent and nature of the conformational variability resulting from crystal packing effects. BPTI is an ideal system to evaluate such factors because of the availability of high resolution X-ray models of five different BPTI structures, each in a different crystal packing environment. Differences observed among the structures are found to be distributed throughout the molecule, although the regions that display most variability are associated with the loop structures (residues 14-17 and 24-29). The regions of structure that show the largest rms deviations from the mean of the five packing motifs correlate well with the presence of intermolecular contacts in the crystal lattice. For most of the molecules there is also a correspondence between a larger number of intermolecular contacts and systematically higher B-factors, although it is not apparent whether this is induced by the crystal contact or results from the fact that the contacts are made predominantly through surface loops. The conformational differences seen among the X-ray models constitute more than local shifts at the lattice contact surfaces, and in fact involve in some cases the making and breaking of intramolecular H-bonds. The magnitudes of the differences among packing models are significantly larger than those usually associated with changes induced by mutagenesis; for instance; the structural differences at the site of mutation observed on removing an internal disulfide from the molecule are significantly less than those associated with lattice contact effects. The crystal packing conformations are compared to representative structures of BPTI generated during a 96-psec molecular dynamics (MD) simulation. This comparison shows a high level of correspondence between the protein flexibility indicated by the X-ray and MD analyses, and specifically between those regions that are most variable. This suggests that the regions that show most variability among the crystal packing models are not artifacts of crystallization, but rather represent true low-energy conformers that have been preferentially selected by crystallization factors.     

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.     

Denaturant-dependent folding of bovine pancreatic trypsin inhibitor mutants with two intact disulfide bonds.

The equilibrium and kinetic behavior of the guanidine hydrochloride (Gdn-HCl) induced unfolding/refolding of four bovine pancreatic trypsin inhibitor (BPTI) mutants was examined by using ultraviolet difference spectroscopy. In three of the mutants, we replaced the buried 30-51 disulfide bond with alanine at position 51 and valine (Val30/Ala51), alanine (Ala30/Ala51), or threonine (Thr30/Ala51) at position 30. For the fourth mutant, the solvent-exposed 14-38 disulfide was substituted by a pair of alanines (Ala14/Ala38). All mutants retained the 5-55 disulfide. Experiments were performed under oxidizing conditions; thus, both the unfolded and folded forms retained two native disulfide bonds. Equilibrium experiments demonstrated that all four mutants were destabilized relative to wild-type BPTI. However, the stability of the 30-51 mutants increased with the hydrophobicity of the residue substituted at position 30. Kinetic experiments showed that all four mutants contained two minor slow refolding phases with characteristics of proline isomerization. The specific behavior of the phases depended on the location of the disulfide bonds. The major unfolding/refolding phase for each of the 30-51 mutants was more than an order of magnitude slower than for Ala14/Ala38 or for BPTI in which the 14-38 disulfide bond was specifically reduced and blocked with iodoacetamide [Jullien, M., & Baldwin, R. L. (1981) J. Mol. Biol. 145, 265-280]. Since this effect is independent of the stability of the protein, it is consistent with a model in which the proper docking of the interior residues of the protein is the rate-limiting step in the folding of these mutants.     

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 dynamic properties of the single disulfide-deficient α-amylase inhibitor [C45A/C73A]tendamistat: An NMR study

Covalent linkages such as disulfide bonds are important for the stabilization of proteins. In the present NMR study we compare the structure and the dynamics of the single disulfide-deficient variant C45A/C73A of the α-amylase inhibitor tendamistat and the wild-type protein, which contains two disulfide bonds (C11-C27 and C45-C73). Complete proton assignment was achieved by standard homonuclear 2D techniques for the variant. Chemical shift differences, intra-strand NOE effects and protected amide proton were used to compare the connectivity of the secondary structure elements of variant and wild-type. Dynamic properties of the wild-type protein were studied by ¹³C_α heteronuclear NOE experiments with carbon in natural abundance. ¹⁵N isotope labeling was necessary to obtain the relaxation parameters of the variant, because of sample degradation. The ¹⁵N resonance assignment was achieved by a ¹⁵N 3D-NOESY-HMQC. Removal of the C45-C73 bond by the C45A/C73A mutation has no influence upon the β-barrel structure of tendamistat beside very local changes at the mutation site. The relaxation data revealed only subtle differences between variant and wild-type on a subnanosecond time scale. Only the N-terminus and G62 in the connecting loop between the anti-parallel β-sheets showed an increased mobility. The results are discussed in respect to thermodynamic stability and the secretion efficiency of tendamistat. Proteins 33:285–294, 1998. © 1998 Wiley-Liss, Inc.     

Investigating the effects of double mutation C30A/C75A on onconase structure: Studies at atomic resolution

The structure of onconase C30A/C75A double mutant has been determined at 1.12Å resolution. The structure has high structural homology to other onconase structures. The changes being results of mutation are relatively small, distributed asymmetrically around the two mutated positions, and they are observed not only in the mutation region but expanded to entire molecule. Different conformation of Lys31 side chain that influences the hydrogen bonding network around catalytic triad is probably responsible for lower catalytic efficiency of double mutant. The decrease in thermal stability observed for the onconase variant might be explained by a less dense packing as manifested by the increase of the molecular volume and the solvent accessible surface area. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 454–460, 2014.     

Electron spin resonance and fluorescence studies of the bound-state conformation of a model protein substrate to the chaperone SecB.

SecB is a homotetrameric, cytosolic chaperone that forms part of the protein translocation machinery in Escherichia coli. We have investigated the bound-state conformation of a model protein substrate of SecB, bovine pancreatic trypsin inhibitor (BPTI) as well as the conformation of SecB itself by using proximity relationships based on site-directed spin-labeling and pyrene fluorescence methods. BPTI is a 58-residue protein and contains three disulfide groups between residues 5 and 55, 14 and 38, as well as 30 and 51. Mutants of BPTI that contained only a single disulfide were reduced, and the free cysteines were labeled with either thiol-specific spin labels or pyrene maleimide. The relative proximity of the labeled residues was studied using either electron spin resonance spectroscopy or fluorescence spectroscopy. The data suggest that SecB binds a collapsed coil of reduced unfolded BPTI, which then undergoes a structural rearrangement to a more extended state upon binding to SecB. Binding occurs at multiple sites on the substrate, and the binding site on each SecB monomer accommodates less than 21 substrate residues. In addition, we have labeled four solvent-accessible cysteine residues in the SecB tetramer and have investigated their relative spatial arrangement in the presence and absence of the substrate protein. The electron spin resonance data suggest that these cysteine residues are in close proximity (15 A) when no substrate protein is bound but move away to a distance of greater than 20 A when SecB binds substrate. This is the first direct evidence of a conformational change in SecB upon binding of a substrate protein.     

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.     

In silico mutation of cysteine residues in the ligand-binding domain of an N-methyl-D-aspartate receptor

The precise nature of redox modulation of N-methyl-d-aspartate (NMDA) receptors is still unclear, although it is thought to be related to the formation and breaking of disulfide bonds. Recent structural data demonstrated the way in which disulfide bonds in the ligand-binding core of the NR1 subunit are arranged. However, the structures were not able to reconcile existing experimental data that examined the effects of mutating these cysteine residues. We have used molecular dynamics (MD) simulations of a series of in silico mutations to try and address this in terms of the current structure of the NR1 ligand-binding domain. A double mutation that removes the disulfide bridge between C744 and C798 gives rise to greater interlobe mobility which was predicted from the crystal structure information but, unexpectedly, also appears to predispose the receptor toward greater flexibility in the hinge region. Removal of the disulfide bond between C454 and C420 did not show any appreciable difference from the "wild-type" simulation, suggesting that removal of this would not change receptor properties, which is in agreement with experimental findings. Furthermore, the position of the C454 side chain could be characterized into discrete rotamers, which may reflect the observation of alternative density in the crystal structure for this residue. Simulations in which two of the disulfide bonds are removed via mutations to alanine (C420A and C436A) resulted in a tendency of the protein to adopt a partially closed conformation.     

NMR characterization of three-disulfide variants of lysozyme, C64A/C80A, C76A/C94A, and C30A/C115A--a marginally stable state in folded proteins

Our earlier NMR study showed that a two-disulfide variant of hen lysozyme containing intra-alpha-domain disulfide bridges, C6-C127 and C30-C115, is partially folded, with the alpha domain tightly folded to the nativelike conformation and the beta domain flexible or unfolded. With a view that the formation of a third disulfide bridge is a key for the accomplishment of the overall chain fold, three-dimensional structures of three-disulfide variants of hen lysozyme lacking one disulfide bridge (C64A/C80A, C76A/C94A, and C30A/C115A) were studied in detail using NMR spectroscopy. Amide hydrogen exchange rates were measured to estimate the degree of conformational fluctuation in a residue-specific manner. The structure of C76A/C94A was found to be quite similar to that of the wild type, except for the peptide segment of residues 74-78. The structure of C64A/C80A was considerably disordered in the entire region of the loop (residues 62-79). Further, it was found that a network of hydrogen bonds within the beta sheet and the 3(10) helix in the beta domain were disrupted and fluctuating. In C30A/C115A, the D helix was unstructured and the interface of the B helix with the D helix was significantly perturbed. However, the structural disorder generated in the hydrophobic core of the alpha domain was prevented by the C helix from propagating toward the beta domain. A marginally stable state in folded proteins is discussed based on the structures remaining in each three-disulfide variant.     

Crystal structure of heat-labile enterotoxin from Escherichia coli with increased thermostability introduced by an engineered disulfide bond in the A subunit.

Cholera toxin (CT) produced by Vibrio cholerae and heat-labile enterotoxin (LT-I), produced by enterotoxigenic Escherichia coli, are AB5 heterohexamers with an ADP-ribosylating A subunit and a GM1 receptor binding B pentamer. These toxins are among the most potent mucosal adjuvants known and, hence, are of interest both for the development of anti-diarrheal vaccines against cholera or enterotoxigenic Escherichia coli diarrhea and also for vaccines in general. However, the A subunits of CT and LT-I are known to be relatively temperature sensitive. To improve the thermostability of LT-I an additional disulfide bond was introduced in the A1 subunit by means of the double mutation N40C and G166C. The crystal structure of this double mutant of LT-I has been determined to 2.0 A resolution. The protein structure of the N40C/G166C double mutant is very similar to the native structure except for a few local shifts near the new disulfide bond. The introduction of this additional disulfide bond increases the thermal stability of the A subunit of LT-I by 6 degrees C. The enhancement in thermostability could make this disulfide bond variant of LT-I of considerable interest for the design of enterotoxin-based vaccines.     

Intriguing conformation changes associated with the trans/cis isomerization of a prolyl residue in the active site of the DsbA C33A mutant

Escherichia coli DsbA belongs to the thioredoxin family and catalyzes the formation of disulfide bonds during the folding of proteins in the bacterial periplasm. It active site (C30-P31-H32-C33) consists of a disulfide bridge that is transferred to newly translocated proteins. The work reported here refers to the DsbA mutant termed C33A that retains, towards reduced unfolded thrombin inhibitor, an activity comparable with the wild-type enzyme. Besides, C33A is also able to form a stable covalent complex with DsbB, the membrane protein responsible for maintaining DsbA in its active form. We have determined the crystal structure of C33A at 2.0 angstroms resolution. Although the general architecture of wt DsbA is conserved, we observe the trans/cis isomerization of P31 in the active site and further conformational changes in the so-called "peptide binding groove" region. Interestingly, these modifications involve residues that are specific to DsbA but not to the thioredoxin family fold. The C33A crystal structure exhibits as well a hydrophobic ligand bound close to the active site of the enzyme. The structural analysis of C33A may actually explain the peculiar behavior of this mutant in regards with its interaction with DsbB and thus provides new insights for understanding the catalytic cycle of DsbA.