Protein disulfide engineering

Improving the stability of proteins is an important goal in many biomedical and industrial applications. A logical approach is to emulate stabilizing molecular interactions found in nature. Disulfide bonds are covalent interactions that provide substantial stability to many proteins and conform to well-defined geometric conformations, thus making them appealing candidates in protein engineering efforts. Disulfide engineering is the directed design of novel disulfide bonds into target proteins. This important biotechnological tool has achieved considerable success in a wide range of applications, yet the rules that govern the stabilizing effects of disulfide bonds are not fully characterized. Contrary to expectations, many designed disulfide bonds have resulted in decreased stability of the modified protein. We review progress in disulfide engineering, with an emphasis on the issue of stability and computational methods that facilitate engineering efforts.     

Structural stability of human alpha-thrombin studied by disulfide reduction and scrambling

Human alpha-thrombin is a very important plasma serine protease, which is involved in physiologically vital processes like hemostasis, thrombosis, and activation of platelets. Knowledge regarding the structural stability of alpha-thrombin is essential for understanding its biological regulation. Here, we investigated the structural and conformational stability of alpha-thrombin using the techniques of disulfide reduction and disulfide scrambling. alpha-Thrombin is composed of a light A-chain (36 residues) and a heavy B-chain (259 residues) linked covalently by an inter-chain disulfide bond (Cys(1)-Cys(122)). The B-chain is stabilized by three intra-chain disulfide bonds (Cys(42)-Cys(58), Cys(168)-Cys(182), and Cys(191)-Cys(220)) (Chymotrypsinogen nomenclature). Upon reduction with dithiothreitol (DTT), alpha-thrombin unfolded in a 'sequential' manner with sequential reduction of Cys(168)-Cys(182) within the B-chain followed by the inter-chain disulfide, generating two distinct partially reduced intermediates, I-1 and I-2, respectively. Conformational stability of alpha-thrombin was investigated by the technique of disulfide scrambling. alpha-Thrombin denatures by scrambling its native disulfide bonds in the presence of denaturant [urea, guanidine hydrochloride (GdmCl) or guanidine thiocyanate (GdmSCN)] and a thiol initiator. During the process, cleavage of the inter-chain disulfide bond and release of the A-chain from B-chain was the foremost event. The three disulfides in the B-chain subsequently scrambled to form three major isomers (designated as X-Ba, X-Bb, and X-Bc). Complete denaturation of alpha-thrombin was observed at low concentrations of denaturants (0.5 M GdmSCN, 1.5 M GdmCl, or 3 M urea) indicating low conformational stability of the protease.     

Gain of function in an ERV/ALR sulfhydryl oxidase by molecular engineering of the shuttle disulfide

The ERV/ALR sulfhydryl oxidase domain is a versatile module adapted for catalysis of disulfide bond formation in various organelles and biological settings. Its four-helix bundle structure juxtaposes a Cys-X-X-Cys dithiol/disulfide motif with a bound flavin adenine dinucleotide (FAD) cofactor, enabling transfer of electrons from thiol substrates to non-thiol electron acceptors. ERV/ALR family members contain an additional di-cysteine motif outside the four-helix-bundle core. Although the location and context of this "shuttle" disulfide differs among family members, it is proposed to perform the same basic function of mediating electron transfer from substrate to the enzyme active site. We have determined by X-ray crystallography the structure of AtErv1, an ERV/ALR enzyme that contains a Cys-X4-Cys shuttle disulfide and oxidizes thioredoxin in vitro, and compared it to ScErv2, which has a Cys-X-Cys shuttle and does not oxidize thioredoxin at an appreciable rate. The AtErv1 shuttle disulfide is in a region of the structure that is disordered and thus apparently mobile and exposed. This feature may facilitate access of protein substrates to the shuttle disulfide. To test whether the shuttle disulfide region is modular and can confer on other enzymes oxidase activity toward new substrates, we generated chimeric enzyme variants combining shuttle disulfide and core elements from AtErv1 and ScErv2 and monitored oxidation of thioredoxin by the chimeras. We found that the AtErv1 shuttle disulfide region could indeed confer thioredoxin oxidase activity on the ScErv2 core. Remarkably, various chimeras containing the ScErv2 Cys-X-Cys shuttle disulfide were found to function efficiently as well. Since neither the ScErv2 core nor the Cys-X-Cys motif is therefore incapable of participating in oxidation of thioredoxin, we conclude that wild-type ScErv2 has evolved to repress activity on substrates of this type, perhaps in favor of a different, as yet unknown, substrate.     

Ubiquitin and protein turnover in synapse function

Enduring modification of synapses is central to long-lasting neural circuit plasticity. Such adaptations include rapid posttranslational modification of existing synaptic proteins over periods of minutes and persisting changes in the abundance of synaptic proteins over hours to days. Recently, ubiquitination and protein degradation have emerged as additional mechanisms for modifying the function and molecular composition of synapses. These recent findings raise intriguing questions as to how enduring changes at synapses are accomplished in the face of robust, ongoing molecular turnover.     

Protein engineering of disulfide bonds in subtilisin BPN'

Five single-disulfide mutants were studied in subtilisin BPN', a cysteine-free, secreted serine protease from Bacillus amyloliquefaciens. The disulfides were engineered between residues 26-232, 29-119, 36-210, 41-80, and 148-243. These bonds connected a variety of secondary structural elements, located in buried or exposed positions at least 10 A from the catalytic Ser-221, and linked residues that were separated by 39 up to 206 amino acids. All disulfide bonds formed in the enzyme when the expressed protein was secreted from Bacillus subtilis, and the disulfides had only minor effects on the enzyme kinetics. Although these disulfide bonds varied by over 50-fold in their equilibrium constants for reduction with dithiothreitol, there was no correlation between the strength of the disulfide bond and the stability it imparted to the enzyme to irreversible inactivation. In some cases, the disulfide-bonded protein was stabilized greatly relative to its reduced counterpart. However, no disulfide mutant was substantially more stable than wild-type subtilisin BPN'. Some of these results can be rationalized by destabilizing effects of the cysteine mutations that disrupt interactions present in the folded enzyme structure. It is also possible that the rate of irreversible inactivation depends upon the kinetics and not the thermodynamics of unfolding and so the entropically stabilizing effect expected from a disulfide bond may not apply.     

Hsp90 structure and function studied by NMR spectroscopy

The molecular chaperone Hsp90 plays a crucial role in folding and maturation of regulatory proteins. Key aspects of Hsp90's molecular mechanism and its adenosine-5'-triphosphate (ATP)-controlled active cycle remain elusive. In particular the role of conformational changes during the ATPase cycle and the molecular basis of the interactions with substrate proteins are poorly understood. The dynamic nature of the Hsp90 machine designates nuclear magnetic resonance (NMR) spectroscopy as an attractive method to unravel both the chaperoning mechanism and interaction with partner proteins. NMR is particularly suitable to provide a dynamic picture of protein-protein interactions at atomic resolution. Hsp90 is rather a challenging protein for NMR studies, due to its high molecular weight and its structural flexibility. The recent technologic advances allowed overcoming many of the traditional obstacles. Here, we describe the different approaches that allowed the investigation of Hsp90 using state-of-the-art NMR methods and the results that were obtained. NMR spectroscopy contributed to understanding Hsp90's interaction with the co-chaperones p23, Aha1 and Cdc37. A particular exciting prospect of NMR, however, is the analysis of Hsp90 interaction with substrate proteins. Here, the ability of this method to contribute to the structural characterization of not fully folded proteins becomes crucial. Especially the interaction of Hsp90 with one of its natural clients, the tumour suppressor p53, has been intensively studied by NMR spectroscopy. This article is part of a Special Issue entitled: Heat Shock Protein 90 (HSP90).     

Thermodynamics of the folding of D-glyceraldehyde-3-phosphate dehydrogenase assisted by protein disulfide isomerase studied by microcalorimetry.

Thermodynamics of the refolding of denatured D-glyceraldehyde 3-phosphate dehydrogenase (GAPDH) assisted by protein disulfide isomerase (PDI), a molecular chaperone, has been studied by isothermal microcalorimetry at different molar ratios of PDI/GAPDH and temperatures using two thermodynamic models proposed for chaperone-substrate binding and chaperone-assisted substrate folding, respectively. The binding of GAPDH folding intermediates to PDI is driven by a large favorable enthalpy decrease with a large unfavorable entropy reduction, and shows strong enthalpy-entropy compensation and weak temperature dependence of Gibbs free energy change. A large negative heat-capacity change of the binding, -156 kJ.mol(-1).K(-1), at all temperatures examined indicates that hydrophobic interaction is a major force for the binding. The binding stoichiometry shows one dimeric GAPDH intermediate per PDI monomer. The refolding of GAPDH assisted by PDI is a largely exothermic reaction at 15.0-25.0 degrees C. With increasing temperature from 15.0 to 37.0 degrees C, the PDI-assisted reactivation yield of denatured GAPDH upon dilution decreases. At 37.0 degrees C, the spontaneous reactivation, PDI-assisted reactivation and intrinsic molar enthalpy change during the PDI-assisted refolding of GAPDH are not detected.     

Conformation of the second disulfide loop in human transforming growth factor-alpha studied by two-dimensional NMR spectroscopy

The determination of the conformation of a cyclic heptadecapeptide derived from the second loop of human transforming growth factor-α, [Ala²¹]-hTGF-α-(16–32), is described. Two-dimensional homonuclear Hartmann–Hahn and rotating-frame cross-relaxation spectroscopy in H₂O were used to obtain the complete proton resonance assignments and the necessary distance constraints between nonbonded hydrogen atoms to derive a conformation without involving any energy minimization. The result is an ellipsoidal-shaped structure with a turn at Gly19 and a bend formed by residues 26–29, Gln-Glu-Asp-Lys. Comparison is made with the second loop of human epidermal growth factor and the results are discussed in terms of receptor binding and biological activity.     

Structure of the complement factor 5a receptor-ligand complex studied by disulfide trapping and molecular modeling

Complement factor 5a (C5a) is an anaphylatoxin that acts by binding to a G protein-coupled receptor, the C5aR. The relative orientation of this ligand-receptor pair is investigated here using the novel technique of disulfide trapping by random mutagenesis (DTRM) and molecular modeling. In the DTRM technique, an unpaired cysteine is introduced in the ligand, and a library of randomly mutagenized receptors is screened to identify mutants that introduce a cysteine at a position in the receptor that allows functional interactions with the ligand. By repeating this analysis at six positions of C5a, we identify six unique sets of intermolecular interactions for the C5a-C5aR complex, which are then compared with an independently developed computational three-dimensional model of the complex. This analysis reveals that the interface of the receptor N terminus with the cysteine-containing ligand molecules is selected from a variety of possible receptor conformations that exist in dynamic equilibrium. In contrast, DTRM identifies a single position in the second extracellular loop of the receptor that interacts specifically with a cysteine probe placed in the C-terminal tail of the C5a ligand.     

Mechanical design of proteins studied by single-molecule force spectroscopy and protein engineering

Mechanical unfolding and refolding may regulate the molecular elasticity of modular proteins with mechanical functions. The development of the atomic force microscopy (AFM) has recently enabled the dynamic measurement of these processes at the single-molecule level. Protein engineering techniques allow the construction of homomeric polyproteins for the precise analysis of the mechanical unfolding of single domains. alpha-Helical domains are mechanically compliant, whereas beta-sandwich domains, particularly those that resist unfolding with backbone hydrogen bonds between strands perpendicular to the applied force, are more stable and appear frequently in proteins subject to mechanical forces. The mechanical stability of a domain seems to be determined by its hydrogen bonding pattern and is correlated with its kinetic stability rather than its thermodynamic stability. Force spectroscopy using AFM promises to elucidate the dynamic mechanical properties of a wide variety of proteins at the single molecule level and provide an important complement to other structural and dynamic techniques (e.g., X-ray crystallography, NMR spectroscopy, patch-clamp).     

N-terminal insertion of alamethicin in channel formation studied using its covalent dimer N-terminally linked by disulfide bond

Alamethicin is supposed to form helix-bundle-type channels by inserting the N terminus into bilayer lipid membranes under sufficient voltages. The N-terminal insertion has been studied with an alamethicin dimer (di-alm) N-terminally linked by a disulfide bond and by the asymmetric addition of dithiothreitol (DTT) and tetrathionate (TT) to the membrane. When di-alm was added to the cis-side membrane, it forms long-lasting channels with the lifetime τ of about 100 ms at cis-positive voltages. The lifetime was reduced to a few milliseconds by addition of DTT to the cis-side membrane, indicating that most of the channels were formed by the monomers (alm-SH) that resulted from the cleavage of the disulfide bond in di-alm. The succeeding addition of TT to the trans-side produced channels of τ=10-20 ms besides the channels of alm-SH. The results suggested that TT reacted with the N-terminal thiol group of alm-SH located at the trans-side of the membrane to alter the lifetime. The N-terminal insertion of alamethicin helices by voltage activation, therefore, was confirmed.     

N-terminal insertion of alamethicin in channel formation studied using its covalent dimer N-terminally linked by disulfide bond

Alamethicin is supposed to form helix-bundle-type channels by inserting the N terminus into bilayer lipid membranes under sufficient voltages. The N-terminal insertion has been studied with an alamethicin dimer (di-alm) N-terminally linked by a disulfide bond and by the asymmetric addition of dithiothreitol (DTT) and tetrathionate (TT) to the membrane. When di-alm was added to the cis-side membrane, it forms long-lasting channels with the lifetime tau of about 100 ms at cis-positive voltages. The lifetime was reduced to a few milliseconds by addition of DTT to the cis-side membrane, indicating that most of the channels were formed by the monomers (alm-SH) that resulted from the cleavage of the disulfide bond in di-alm. The succeeding addition of TT to the trans-side produced channels of tau=10-20 ms besides the channels of alm-SH. The results suggested that TT reacted with the N-terminal thiol group of alm-SH located at the trans-side of the membrane to alter the lifetime. The N-terminal insertion of alamethicin helices by voltage activation, therefore, was confirmed.     

Folding kinetics of staphylococcal nuclease studied by tryptophan engineering and rapid mixing methods

To monitor the development of tertiary structural contacts during folding, a unique tryptophan residue was introduced at seven partially buried locations (residues 15, 27, 61, 76, 91, 102 and 121) of a tryptophan-free variant of staphylococcal nuclease (P47G/P117G/H124L/W140H). Thermal unfolding measurements by circular dichroism indicate that the variants are destabilized, but maintain the ability to fold into a native-like structure. For the variants with Trp at positions 15, 27 and 61, the intrinsic fluorescence is significantly quenched in the native state due to close contact with polar side-chains that act as intramolecular quenchers. All other variants exhibit enhanced fluorescence under native conditions consistent with burial of the tryptophan residues in an apolar environment. The kinetics of folding was observed by continuous and stopped-flow fluorescence measurements over refolding times ranging from 100 micros to 10 s. The folding kinetics of all variants is quantitatively described by a mechanism involving a major pathway with a series of intermediate states and a minor parallel channel. The engineered tryptophan residues in the beta-barrel and the N-terminal part of the alpha-helical domain become partially shielded from the solvent at an early stage (<1 ms), indicating that this region undergoes a rapid collapse. For some variants, a major increase in fluorescence coincides with the rate-limiting step of folding on the 100 ms time scale, indicating that these tryptophan residues are buried only during the late stages of folding. Other variants exhibit a transient increase in fluorescence during the 10 ms phase followed by a decrease during the rate-limiting phase. These observations are consistent with burial of these probes in a collapsed, but loosely packed intermediate, followed by the rate-limiting formation of the densely packed native core, which brings the tryptophan residues into close contact with intramolecular quenchers.     

Linear free energy relationships in enzyme binding interactions studied by protein engineering.

Experiments on mutants of tyrosyl-tRNA synthetase have shown that there can be linear free energy relationships (LFERs) between changes in activation free energies and changes in binding energies when groups are deleted that bind to non-reacting parts of the substrate (Fersht et al., 1986, 1987). It has now been proposed (Straub and Karplus, 1990) that such LFERs can occur for the mutation of hydrogen bonding groups only for the limiting examples of Br?nsted beta of 0, 1 or infinity, and that fractional values of beta are not permissible. The reasoning behind this is that the energy of a hydrogen bond is not linear with distance and the (false) premise that an LFER requires that there is a linear relationship between bond energy and distance. We show from a simple model how LFERs can arise for binding interactions and how they can give fractional values of beta, in accord with experimental evidence. An LFER occurs between binding and catalysis when a set of interactions exists in which each member contributes to the binding energy of the transition state the same fraction of the binding energy it contributes to the products (both relative to the ground state).     

Cell Regulation by Phosphotyrosine-Targeted Ubiquitin Ligases

Three classes of E3 ubiquitin ligases, members of the Cbl, Hakai, and SOCS-Cul5-RING ligase families, stimulate the ubiquitination of phosphotyrosine-containing proteins, including receptor and nonreceptor tyrosine kinases and their phosphorylated substrates. Because ubiquitination frequently routes proteins for degradation by the lysosome or proteasome, these E3 ligases are able to potently inhibit tyrosine kinase signaling. Their loss or mutational inactivation can contribute to cancer, autoimmunity, or endocrine disorders, such as diabetes. However, these ligases also have biological functions that are independent of their ubiquitination activity. Here we review relevant literature and then focus on more-recent developments in understanding the structures, substrates, and pathways through which the phosphotyrosine-specific ubiquitin ligases regulate diverse aspects of cell biology.     

Exploring the folding pathway of green fluorescent protein through disulfide engineering

We have introduced two disulfide crosslinks into the loop regions on opposite ends of the beta barrel in superfolder green fluorescent protein (GFP) in order to better understand the nature of its folding pathway. When the disulfide on the side opposite the N/C‐termini is formed, folding is 2× faster, unfolding is 2000× slower, and the protein is stabilized by 16 kJ/mol. But when the disulfide bond on the side of the termini is formed we see little change in the kinetics and stability. The stabilization upon combining the two crosslinks is approximately additive. When the kinetic effects are broken down into multiple phases, we observe Hammond behavior in the upward shift of the kinetic m‐value of unfolding. We use these results in conjunction with structural analysis to assign folding intermediates to two parallel folding pathways. The data are consistent with a view that the two fastest transition states of folding are "barrel closing" steps. The slower of the two phases passes through an intermediate with the barrel opening occurring between strands 7 and 8, while the faster phase opens between 9 and 4. We conclude that disulfide crosslink‐induced perturbations in kinetics are useful for mapping the protein folding pathway.     

The role of disulfide bond in the amyloidogenic state of beta(2)-microglobulin studied by heteronuclear NMR

beta(2)-Microglobulin (beta2-m) is a major component of dialysis-related amyloid fibrils. Although recombinant beta2-m forms needle-like fibrils by in vitro extension reaction at pH 2.5, reduced beta2-m, in which the intrachain disulfide bond is reduced, cannot form typical fibrils. Instead, thinner and flexible filaments are formed, as shown by atomic force microscopy images. To clarify the role of the disulfide bond in amyloid fibril formation, we characterized the conformations of the oxidized (intact) and reduced forms of beta2-m in the acid-denatured state at pH 2.5, as well as the native state at pH 6.5, by heteronuclear NMR. [(1)H]-(15)N NOE at the regions between the two cysteine residues (Cys25-Cys80) revealed a marked difference in the pico- and nanosecond time scale dynamics between that the acid-denatured oxidized and reduced states, with the former showing reduced mobility. Intriguingly, the secondary chemical shifts, DeltaCalpha, DeltaCO, and DeltaHalpha, and (3)J(HNHalpha) coupling constants indicated that both the oxidized and reduced beta2-m at pH 2.5 have marginal alpha-helical propensity at regions close to the C-terminal cysteine, although it is a beta-sheet protein in the native state. The results suggest that the reduced mobility of the denatured state is an important factor for the amylodogenic potential of beta2-m, and that the marginal helical propensity at the C-terminal regions might play a role in modifying this potential.