Protein engineering of a disulfide bond in a beta/alpha-barrel protein

Multi-template approach to modeling engineered disulfide bonds

The key issue for disulfide bond engineering is to select the most appropriate location in the protein. By surveying the structure of experimentally engineered disulfide bonds, we found about half of them that have geometry incompatible with any native disulfide bond geometry. To improve the current prediction methods that tend to apply either ideal geometrical or energetical criteria to single three-dimensional structures, we have combined a novel computational protocol with the usage of multiple protein structures to take into account protein backbone flexibility. The multiple structures can be selected from either independently determined crystal structures for identical proteins, models of nuclear magnetic resonance experiments, or crystal structures of homology-related proteins. We have validated our approach by comparing the predictions with known disulfide bonds. The accuracy of prediction for native disulfide bonds reaches 99.6%. In a more stringent test on the reported engineered disulfide bonds, we have obtained a success rate of 93%. Our protocol also determines the oxido-reduction state of a predicted disulfide bond and the corresponding mutational cost. From the energy ranking, the user can easily choose top predicted sites for mutagenesis experiments. Our method provides information about local stability of the engineered disulfide bond surroundings.

The effect of engineered disulfide bonds on the stability of Drosophila melanogaster acetylcholinesterase

Background

Acetylcholinesterase is irreversibly inhibited by organophosphate and carbamate insecticides allowing its use in biosensors for detection of these insecticides. Drosophila acetylcholinesterase is the most sensitive enzyme known and has been improved by in vitro mutagenesis. However, its stability has to be improved for extensive utilization.

Results

To create a disulfide bond that could increase the stability of the Drosophila melanogaster acetylcholinesterase, we selected seven positions taking into account first the distance between Cβ of two residues, in which newly introduced cysteines will form the new disulfide bond and second the conservation of the residues in the cholinesterase family. Most disulfide bonds tested did not increase and even decreased the stability of the protein. However, one engineered disulfide bridge, I327C/D375C showed significant stability increase toward denaturation by temperature (170 fold at 50°C), urea, organic solvent and provided resistance to protease degradation. The new disulfide bridge links the N-terminal domain (first 356 aa) to the C-terminal domain. The quantities produced by this mutant were the same as in wild-type flies.

Conclusion

Addition of a disulfide bridge may either stabilize or unstabilize proteins. One bond out of the 7 tested provided significant stabilisation.

Novel yeast cell-based assay to screen for inhibitors of human cytomegalovirus protease in a high-throughput format

The protease encoded by the human cytomegalovirus (HCMV) is an attractive target for antiviral drug development because of its essential function in viral replication. We describe here a cellular assay in the yeast Saccharomyces cerevisiae for the identification of small molecule inhibitors of HCMV protease by conditional growth in selective medium. In this system, the protease cleavage sequence is inserted into the N-(5'-phosphoribosyl)anthranilate isomerase (Trp1p), a yeast protein essential for cell proliferation in the absence of tryptophan. Coexpression of HCMV protease with the engineered Trp1p substrate in yeast cells results in site-specific cleavage and functional inactivation of the Trp1p enzyme, thereby leading to an arrest of cell proliferation. This growth arrest can be suppressed by the addition of validated HCMV protease inhibitors. The growth selection system presented here provides the basis for a high-throughput screen to identify HCMV protease inhibitors that are active in eukaryotic cells.

Cooperativity of the oxidization of cysteines in globular proteins

Based on the 639 non-homologous proteins with 2910 cysteine-containing segments of well-resolved three-dimensional structures, a novel approach has been proposed to predict the disulfide-bonding state of cysteines in proteins by constructing a two-stage classifier combining a first global linear discriminator based on their amino acid composition and a second local support vector machine classifier. The overall prediction accuracy of this hybrid classifier for the disulfide-bonding state of cysteines in proteins has scored 84.1% and 80.1%, when measured on cysteine and protein basis using the rigorous jack-knife procedure, respectively. It shows that whether cysteines should form disulfide bonds depends not only on the global structural features of proteins but also on the local sequence environment of proteins. The result demonstrates the applicability of this novel method and provides comparable prediction performance compared with existing methods for the prediction of the oxidation states of cysteines in proteins.

Prediction of the disulfide-bonding state of cysteines in proteins based on dipeptide composition

In this paper, a novel approach has been introduced to predict the disulfide-bonding state of cysteines in proteins by means of a linear discriminator based on their dipeptide composition. The prediction is performed with a newly enlarged dataset with 8114 cysteine-containing segments extracted from 1856 non-homologous proteins of well-resolved three-dimensional structures. The oxidation of cysteines exhibits obvious cooperativity: almost all cysteines in disulfide-bond-containing proteins are in the oxidized form. This cooperativity can be well described by protein's dipeptide composition, based on which the prediction accuracy of the oxidation form of cysteines scores as high as 89.1% and 85.2%, when measured on cysteine and protein basis using the rigorous jack-knife procedure, respectively. The result demonstrates the applicability of this new relatively simple method and provides superior prediction performance compared with existing methods for the prediction of the oxidation states of cysteines in proteins.

Adaptation of Class-13 α-Amylases to Diverse Living Conditions

There are currently 35 available nonredundant molecular structures of class-13 alpha-amylases (EC 3.2.1.1), mostly from microbial organisms living under a wide range of environmental conditions. One of the most recent additions has been the first alpha-amylase structure of a hyperthermophilic archaeon [Linden et al., J. Biol. Chem. 2003, 278, 9875-9884]. The structure has been used for comparative analyses with a representative set of three alpha-amylases from thermophilic, mesophilic and psychrophilic sources to identify molecular parameters for environmental adaptation. Our analysis supports generally observed trends such as an increase in structural compactness as well as an increase in salt bridges in order to cope with high-temperature conditions. The two representative thermophilic structures used in this comparative study have independently evolved di-metal centres--not present in the mesophilic and psychrophilic structures--in the vicinity of the active site. These observations may provide impetus for the design of alpha-amylases with improved molecular properties to enhance their utility in biotechnological processes.

Stabilization of a (betaalpha)8-barrel protein by an engineered disulfide bridge

The aim of this study was to increase the stability of the thermolabile (betaalpha)8-barrel enzyme indoleglycerol phosphate synthase from Escherichia coli by the introduction of disulfide bridges. For the design of such variants, we selected two out of 12 candidates, in which newly introduced cysteines potentially form optimal disulfide bonds. These variants avoid short-range connections, substitutions near catalytic residues, and crosslinks between the new and the three parental cysteines. The variant linking residues 3 and 189 fastens the N-terminus to the (betaalpha)8-barrel. The rate of thermal inactivation at 50 degrees C of this variant with a closed disulfide bridge is 65-fold slower than that of the reference dithiol form, but only 13-fold slower than that of the parental protein. The near-ultraviolet CD spectrum, the reactivity of parental buried cysteines with Ellman's reagent as well as the decreased turnover number indicate that the protein structure is rigidified. To confirm these data, we have solved the X-ray structure to 2.1-A resolution. The second variant was designed to crosslink the terminal modules betaalpha1 and betaalpha8. However, not even the dithiol form acquired the native fold, possibly because one of the targeted residues is solvent-inaccessible in the parental protein.

Amino acid neighbours and detailed conformational analysis of cysteines in proteins

Here we present an investigation of the contacts that cysteines make with residues in their three-dimensional environment and a comprehensive analysis of the conformational features of 351 disulphide bridges in 131 non-homologous single-chain protein structures. Upstream half-cystines preferentially have downstream neighbours, whereas downstream half-cystines have mainly upstream neighbours. Non-disulphide bridged cysteines (free cysteines) have no preference for upstream or downstream neighbours. Free cysteines have more contacts to non-polar residues and fewer contacts to polar/charged residues than half-cystines, which correlates with our observation that free cysteines are more buried than half-cystines. Free cysteines prefer to be located in alpha-helices while no clear preference is observed for half-cystines. Histidine and methionine are preferentially seen nearby free cysteines. Tryptophan is found preferentially nearby half-cystines. We have merged sequential and spatial information, and highly interesting novel patterns have been discovered. The number of cysteines per protein is typically an even number, peaking at four. The number of residues separating two half-cystines is preferentially 11 and 16. Left-handed and right-handed disulphide bridges display different conformational parameters. Here we present side chain torsion angle information based on a 5-12 times larger number of disulphide bridges than has previously been published. Considering the importance of cysteines for maintaining the 3D-structural scaffold of proteins, it is essential to have as accurate information as possible concerning the packing and conformational preferences. The present work may provide key information for engineering the protein environment around cysteines.

The hyperthermostable indoleglycerol phosphate synthase from Thermotoga maritima is destabilized by mutational disruption of two solvent-exposed salt bridges

The recombinantly expressed protein indoleglycerol phosphate synthase from the hyperthermophilic bacterium Thermotoga maritima (tIGPS) was purified and characterized with respect to oligomerization state, catalytic properties and thermostability. This enzyme from the biosynthetic pathway of tryptophan is a monomer in solution. In contrast to IGPS from the hyperthermophilic archaeon Sulfolobus solfataricus, tIGPS shows high catalytic activity at room temperature and only weak product inhibition. In order to test the hypothesis that salt bridges in a critical context contribute to the high thermostability of tIGPS, two solvent-exposed salt bridges were selected, based on its three-dimensional structure, for individual disruption by site-directed mutagenesis. The first salt bridge fixes the N terminus to the core of the protein, and the second serves as a clamp between helices alpha1 and alpha8, which are widely separated in sequence but adjacent in the (betaalpha)8-barrel. Kinetics of irreversible heat inactivation reveal that the salt bridge crosslinking helices alpha1 and alpha8 stabilizes tIGPS more strongly than that tethering the N terminus.

Beta-turn propensities as paradigms for the analysis of structural motifs to engineer protein stability

The thermodynamic stability of a protein provides an experimental metric for the relationship of protein sequence and native structure. We have investigated an approach based on an analysis of the structural database for stability engineering of an immunoglobulin variable domain. The most frequently occurring residues in specific positions of beta-turn motifs were predicted to increase the folding stability of mutants that were constructed by site-directed mutagenesis. Even in positions in which different residues are conserved in immunoglobulin sequences, the predictions were confirmed. Frequently, mutants with increased beta-turn propensities display increased folding cooperativities, suggesting pronounced effects on the unfolded state independent of the expected effect on conformational entropy. We conclude that structural motifs with predominantly local interactions can serve as templates with which patterns of sequence preferences can be extracted from the database of protein structures. Such preferences can predict the stability effects of mutations for protein engineering and design.

Computer-aided modeling of structure stabilizing disulfide bonds in recombinant human interferon-gamma

We present a general search algorithm for possible insertion sites of disulfide bonds in proteins based on the coordinates of the solved X-ray or NMR structure, allowing the insertion of disulfide bonds with a minimum of conformational tension and backbone rearrangements. The FORTRAN 77 program "Ssuitable' was written for this purpose. This methodological approach was applied to recombinant human interferon-gamma (rhu-IFN-gamma), a cytokine of great pharmaceutical interest with a wide variety of biological activities including antiviral, antiproliferative and immunomodulatory effects. A model based on the C alpha-coordinates obtained from the Brookhaven data base was built. Four different insertion sites were selected in the model, connecting the two subunits of the homodimer. The thermodynamic stability of rhu-IFN-gamma is low, limiting its clinical application. We expect that the insertion of additional new disulfide bonds will enhance the thermodynamic stability as well as protect the protein against proteolytic degradation.

Site-directed mutants designed to test back-door hypotheses of acetylcholinesterase function

The location of the active site of the rapid enzyme, acetylcholinesterase, near the bottom of a deep and narrow gorge indicates that alternative routes may exist for traffic of substrate, products or solute into and out of the gorge. Molecular dynamics suggest the existence of a shutter-like back door near Trp84, a key- residue in the binding site for acetylcholine, in the Torpedo californica enzyme. The homology of the omega loop, bearing Trp84, with the lid which sequesters the substrate in neutral lipases displaying structural homology with acetylcholinesterase, suggests a flap-like back door. Both possibilities were examined by site-directed mutagenesis. The shutter-like back door was tested by generating a salt bridge which might impede opening of the shutter. The flap-like back door was tested by de novo insertion of a disulfide bridge which tethered the omega loop to the body of the enzyme. Neither type of mutation produced significant changes in catalytic activity, thus failing to provide experimental support for either back door model. Molecular dynamics revealed, however, substantial mobility of the omega loop in the immediate vicinity of Trp84, even when the loop was tethered, supporting the possibility that access to the active site, involving limited movement of a segment of the loop, is indeed possible.

Entropic effects of disulphide bonds on protein stability

To measure the thermodynamic consequences of the reduction in the number of polypeptide-chain conformations that accompanies protein folding, we developed a method called loop permutation analysis. In this approach, the stabilizing contributions of three engineered disulphide bonds were compared in extended and circularly permutated mutants of phage T4 lysozyme. The observed differences in disulphide contributions, although qualitatively consistent with theoretical estimates, were not solely proportional to the differences in loop length. These findings suggest that in addition to the length of the chain, the polypeptide sequence may influence the energetic consequences of conformational restrictions.

Detection of an intermediate in the folding of the (beta alpha)8-barrel N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli

We have used thermodynamic and kinetic techniques to monitor the guanidinium chloride induced (GdmCl-induced) denaturation of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli (ePRAI). Although CD-monitored equilibrium denaturation curves are consistent with cooperative unfolding of the protein centered at 1.45 M GdmCl, fluorescence readings drop by over 25% in the region preceding the CD-monitored transition, suggesting non-two-state behavior. Kinetics experiments measure a slow relaxation rate with negative fluorescence amplitude when protein is diluted from 0 to 0.5 M GdmCl, corroborating results from equilibrium conditions. Detection of several unfolding and refolding rates in final GdmCl concentrations from 0 to 5.0 M indicates the presence of at least one intermediate along unfolding and refolding pathways. GdmCl dependence of the relaxation rates can be explained most easily by a nonsequential mechanism for ePRAI unfolding, though a sequential mechanism cannot be ruled out. The data corroborate the fragment complementation studies of Eder and Kirschner [Eder, J., & Kischner, K. (1992) Biochemistry 31, 3617-3625], which are consistent with unfolding of the C-terminal portion of a yeast-derived PRAI in its folding intermediate. In ePRAI, such partial unfolding would expose W391 to quenching by solvent molecules; W356, ePRAI's other tryptophan, is buried in the hydrophobic core and is unlikely to be affected by local changes in structure. A C-terminally unfolded folding intermediate has been demonstrated in the folding of tryptophan synthase (alpha-subunit), a related beta alpha-barrel enzyme.(ABSTRACT TRUNCATED AT 250 WORDS)

Design of interchain disulfide bonds in the framework region of the Fv fragment of the monoclonal antibody B3

The Fv fragments are the smallest units of antibodies that retain the specific antigen binding characteristics of the whole molecule and are being used for the diagnosis and therapy of human diseases. These are noncovalently associated heterodimers of the heavy (VH) and the light (VL) chain variable domains, which, without modification, tend to dissociate, unfold, and/or nonspecifically aggregate. The fragment is usually stabilized by producing it as a single chain recombinant molecule in which the two chains are linked by means of a short polypeptide linker. An alternative strategy is to connect the two chains by means of an interchain disulfide bond. We used molecular graphics and other modeling tools to identify two possible interchain disulfide bond sites in the framework region of the Fv fragment of the monoclonal mouse antibody (mAb) B3. The mAb B3 binds to many human cancer cells and is being used in the development of a new anticancer agent. The two sites identified are VH44-VL105 and VH111-VL48. (VH44-VL100 and VH105-VL43 in the numbering scheme of Kabat et al., "Sequence of Proteins of Immunological Interest," U.S. DHHS, NIH publication No. 91-3242, 1991). This design was recently tested using the chimeric protein composed of a truncated form of Pseudomonas exotoxin and the Fv fragment of mAb B3 with the engineered disulfide bond at VH44-VL105 (Brinkmann et al., Proc. Natl. Acad. Sci. U.S.A. 90:7538, 1993). The chimeric toxin was found to be just as active as the corresponding single chain counterpart and considerably more stable. Because these disulfide bond sites are in the framework region, they can be located from sequence alignment alone. We expect that the disulfide bond at these sites will stabilize the Fv fragment of most antibodies and the antigen-specific portion of the T-cell receptors, which are homologous.

Disulfide bonds and the stability of globular proteins

An understanding of the forces that contribute to stability is pivotal in solving the protein-folding problem. Classical theory suggests that disulfide bonds stabilize proteins by reducing the entropy of the denatured state. More recent theories have attempted to expand this idea, suggesting that in addition to configurational entropic effects, enthalpic and native-state effects occur and cannot be neglected. Experimental thermodynamic evidence is examined from two sources: (1) the disruption of naturally occurring disulfides, and (2) the insertion of novel disulfides. The data confirm that enthalpic and native-state effects are often significant. The experimental changes in free energy are compared to those predicted by different theories. The differences between theory and experiment are large near 300 K and do not lend support to any of the current theories regarding the stabilization of proteins by disulfide bonds. This observation is a result of not only deficiencies in the theoretical models but also from difficulties in determining the effects of disulfide bonds on protein stability against the backdrop of numerous subtle stabilizing factors (in both the native and denatured states), which they may also affect.

Protein engineering for unusual environments

The ability to use proteins in unusual or non-natural environments greatly expands their potential applications in biotechnology. Because natural selection has neither maximized the stability of proteins nor optimized them to function under unusual conditions, there is considerable room for their improvement by protein engineering. Significant advances reported within the past year include a dramatic demonstration of a protein's ability to tolerate changes in its amino acid sequence, large increases in protein stability, and the use of random mutagenesis to obtain novel enzymatic properties. Approaches using random or site-directed mutagenesis have been successful in generating proteins able to function in an extended range of environments.

Disulfide bond contribution to protein stability: positional effects of substitution in the hydrophobic core of the two-stranded alpha-helical coiled-coil

To investigate the positional effect of the disulfide bond on the structure and stability of a two-stranded alpha-helical coiled-coil, an interchain disulfide bond was systematically introduced into the hydrophobic core of a de novo designed model coiled-coil at the N-terminus (position 2), C-terminus (position 33), and nonterminal positions a (positions 9, 16, 23, and 30) and d (positions 5, 12, 19, and 26). The rate of formation of a disulfide bond is faster at position d compared to at the corresponding position a under nondenaturing conditions, suggesting that position d is more suitable for engineering a disulfide bond. The structure and stability of the reduced and oxidized coiled-coils were determined by circular dichroism studies in the absence and presence of guanidine hydrochloride. Our results demonstrate that the improvement of protein stability by introduction of a disulfide bond is very relevant to its location and the most effective disulfide bonds are those that can be introduced in the hydrophobic core without any disruption of the protein structure. The disulfide bond at position d with near-optimal geometry does not perturb the coiled-coil structure and makes the largest contribution to coiled-coil stability. In contrast, the inappropriate geometry of the disulfide bond at nonterminal position a introduces a high strain energy on the disulfide bond which disrupts the coiled-coil structure. At positions a, the closer the disulfide bridge is to the center of the coiled-coil, the larger the disruption on the coiled-coil structure and the smaller the contribution the disulfide bond makes to coiled-coil stability. The computer modeling results also suggest that an insertion of an interchain disulfide bond at position a in the GCN4 leucine zipper X-ray structure has a higher potential energy than insertion at position d. The energy-minimized coiled-coil structure with an interchain disulfide bond at position a has a larger root mean square difference from the X-ray structure of GCN4 than the coiled-coil with a disulfide bond at position d. Because interhelical interactions are common in globular proteins as well as coiled-coils, the results obtained in this study will have general utility for selecting the sites for engineering disulfide bonds between alpha-helices.