Stabilization of a Tetrameric Malate Dehydrogenase by Introduction of a Disulfide Bridge at the Dimer–Dimer Interface

Genome-Wide Characterization and Gene Expression Analyses of Malate Dehydrogenase (MDH) Genes in Low-Phosphorus Stress Tolerance of Chinese Fir (Cunninghamia lanceolata)

Malate dehydrogenase (MDH) genes play vital roles in developmental control and environmental stress tolerance in sessile plants by modulating the organic acid-malic acid level. However, MDH genes have not yet been characterized in gymnosperm, and their roles in nutrient deficiency are largely unexplored. In this study, 12 MDH genes were identified in Chinese fir (Cunninghamia lanceolata), namely, ClMDH-1, -2, -3, …, and -12. Chinese fir is one of the most abundant commercial timber trees in China, and low phosphorus has limited its growth and production due to the acidic soil of southern China. According to the phylogenetic analysis, MDH genes were classified into five groups, and Group 2 genes (ClMDH-7, -8, -9, and 10) were only found to be present in Chinese fir but not in Arabidopsis thaliana and Populus trichocarpa. In particular, the Group 2 MDHs also had specific functional domains-Ldh_1_N (malidase NAD-binding functional domain) and Ldh_1_C (malate enzyme C-terminal functional domain)-indicating a specific function of ClMDHs in the accumulation of malate. All ClMDH genes contained the conserved MDH gene characteristic functional domains Ldh_1_N and Ldh_1_C, and all ClMDH proteins exhibited similar structures. Twelve ClMDH genes were identified from eight chromosomes, involving fifteen ClMDH homologous gene pairs, each with a Ka/Ks ratio of <1. The analysis of cis-elements, protein interactions, and transcription factor interactions of MDHs showed that the ClMDH gene might play a role in plant growth and development, and in response to stress mechanisms. The results of transcriptome data and qRT-PCR validation based on low-phosphorus stress showed that ClMDH1, ClMDH6, ClMDH7, ClMDH2, ClMDH4, ClMDH5, ClMDH10 and ClMDH11 were upregulated under low-phosphorus stress and played a role in the response of fir to low-phosphorus stress. In conclusion, these findings lay a foundation for further improving the genetic mechanism of the ClMDH gene family in response to low-phosphorus stress, exploring the potential function of this gene, promoting the improvement of fir genetics and breeding, and improving production efficiency.

Enzymatic Characterization and Coenzyme Specificity Conversion of a Novel Dimeric Malate Dehydrogenase from Bacillus subtilis

Malate is an important material to various industrials and clinical applications. Bacillus subtilis is a widely used biocatalyst tool for chemical production. However, the specific enzymatic properties of malate dehydrogenase from Bacillus subtilis (BsMDH) remain largely unknown. In the present study, BsMDH was cloned, recombinantly expressed and purified to test its enzymatic properties. The molecular weight of single unit of BsMDH was 34,869.7 Da. Matrix-Assisted Laser-Desorption Ionization-Time-of-Flight Mass Spectrometry and gel filtration analysis indicated that the recombinant BsMDH could form dimers. The k_cat/K_m values of oxaloacetate and NADH were higher than those of malate and NAD⁺, respectively, indicating a better catalysis in the direction of malate synthesis than the reverse. Furthermore, six BsMDH mutants were constructed with the substitution of amino acids at the coenzyme binding site. Among them, BsMDH-T7 showed a greatly higher affinity and catalysis efficiency to NADPH than NADH with the degree of alteration of 2039, suggesting the shift of the coenzyme dependence from NADH to NADPH. In addition, BsMDH-T7 showed a relatively lower K_m value, but a higher k_cat and k_cat/K_m than NADPH-dependent MDHs from Thermus flavus and Corynebacterium glutamicum. Overall, these results indicated that BsMDH and BsMDH-T7 mutant might be promising enzymes for malate production.

Heterodimerizing helices as tools for nanoscale control of the organization of protein-protein and protein-quantum dots

In this study, we tested the possibility of creating complexes of two proteins by fusing them with heterodimerizing helices. We used the fluorescent proteins GFP and mCHERRY expressed with a His-tag as our model system. We added heterodimer-forming sequences at the C- or N- termini of the proteins, opposite to the His-tag position. Heterodimerization was tested for both helices at the C-terminus or at the N- terminus and C-terminus. We observed complex formation with a nanomolar dissociation constant in both cases that was higher by one order of magnitude than the Kds measured for helices alone. The binding of two C-terminal helices was accompanied by an increased enthalpy change. The binding between helices could be stabilized by introducing an additional turn of the helix with cysteine, which was capable of forming disulphide bridges. Covalently linked proteins were obtained using this strategy and observed using fluorescence cross-correlation spectroscopy. Finally, we demonstrated the formation of complexes of protein dimers and quantum dots.

Improving the thermostability of Serratia marcescens B4A chitinase via G191V site-directed mutagenesis

Chitinases with high thermostability are important for many industrial and biotechnological applications. This study was conducted to enhance the stability of Serratia marcescens B4A chitinase by site directed mutagenesis of G191 V. Further characterization showed that the thermal stability of the mutant showed marked increase of about 5 and 15 fold at 50 and 60 °C respectively, while the optimum temperature and pH was retained. Kinetic analysis showed decreased K_m and V_max of the mutant in comparison with the wild type chitinase of about 1.3 and 3 fold, respectively. Based on structural prediction, it was speculated that this replacement shortened an important loop concomitant with the extension of adjacent β sheets. Accordingly, a higher thermostability of G191 V up to 90 °C supporting the decreased flexibility of unfolded state was also indicated. Finally, a practical proof of kinetic and thermal stabilization of chitinase was provided through decreased flexibility and entropic stabilization of its surface loops.

Stabilization of Multimeric Proteins via Intersubunit Cyclization

Proteins with high catalytic efficiency and selectivity under mild conditions have long been appreciated by industrial and medicinal fields. These proteins, which are commonly multimeric, often possess low stability, impeding wider application. Currently, strategies to improve the stability of multimeric proteins concentrate on enhancing the interaction at internal interface of the subunits. In this report, we confirmed that the largely underestimated subunit terminal ends are as significant as the internal interface for protein stability. By connecting both the terminal ends and internal interface of subunits, the tetrameric Leifsonia alcohol dehydrogenase (LnADH) protein can been cyclized into a rigid form with significantly improved thermostability and resilience. The improvement in the temperature at which enzyme activity is reduced to 50% after a 15-min heat treatment (T₅₀¹⁵) and melting temperature (T_m ) of the modified protein was 18°C and 23.3°C, respectively, which is superior to the results achieved by normal protein engineering. Our study provided a novel strategy to effectively improve the stability of multimeric proteins, which is suitable not only for the short-chain dehydrogenase/reductase (SDR) family but also other classes of proteins with close terminal ends.IMPORTANCE Industrially interesting proteins are generally multimeric proteins; however, their applications are often restricted due to low stability caused by the natural tendency of subunit disassociation. Current approaches targeting this problem mainly focus on enhancing the internal interfaces of the subunits to avoid their disassociation. In this study, we identified and confirmed the external interface to be significant for improving the stability of multimeric proteins. By connecting the terminal ends and internal interface with disulfide bonds, we found that the multimeric protein LnADH cyclized into a robust monomeric-like form, resulting in superior thermostability compared to traditional protein engineering. This intersubunit cyclization approach is efficient and easy to perform, providing a novel method for engineering many important classes of multimeric proteins.

Rational mutagenesis by engineering disulphide bonds improves Kluyveromyces lactis beta-galactosidase for high-temperature industrial applications

Kluyveromyces lactis β-galactosidase (Kl-β-Gal) is one of the most important enzymes in the dairy industry. The poor stability of this enzyme limits its use in the synthesis of galactooligosaccharides (GOS) and other applications requiring high operational temperature. To obtain thermoresistant variants, a rational mutagenesis strategy by introducing disulphide bonds in the interface between the enzyme subunits was used. Two improved mutants, R116C/T270C and R116C/T270C/G818C, had increased half-lives at 45 °C compared to Kl-β-Gal (2.2 and 6.8 fold increases, respectively). Likewise, Tm values of R116C/T270C and R116C/T270C/G818C were 2.4 and 8.5 °C, respectively, higher than Kl-β-Gal Tm. Enrichment in enzymatically active oligomeric forms in these mutant variants also increased their catalytic efficiency, due to the reinforcement of the interface contacts. In this way, using an artificial substrate (p-nitrophenyl-β-D-galactopyranoside), the Vmax values of the mutants were ~1.4 (R116C/T270C) and 2 (R116C/T270C/G818C) fold higher than that of native Kl-β-Gal. Using the natural substrate (lactose) the Vmax for R116C/T270C/G818C almost doubled the Vmax for Kl-β-Gal. Validation of these mutant variants of the enzyme for their use in applications that depend on prolonged incubations at high temperatures was achieved at the laboratory scale by monitoring their catalytic activity in GOS synthesis.

Engineering disulfide bonds in Selenomonas ruminantium β-xylosidase by experimental and computational methods

Homotetrameric β-xylosidase from Selenomonas ruminantium (SXA) is one of the most efficient enzymes known for the hydrolysis of cell wall hemicellulose. SXA shows a rapid rate of activity loss at temperatures above 50°C. In this study, we have introduced two inter-subunit disulfide bridges with one point mutation. Lys237 was chosen to be replaced with cysteine since it interacts with the same residue in the opposite subunit. While pH optimum, temperature profile and catalytic efficiency of the mutated variant were similar to the native enzyme, the mutated enzyme showed about 40% increase in thermal stability at 55°C. Our results showed that introduction of a single residue mutation in structure of SXA results in appearance of two disulfide bonds at dimer-dimer interface of the enzyme. Coarse-grained molecular dynamics (CG-MD) simulations also proved lower amounts of root mean square fluctuation (RMSF) for position 237 and potential energy for mutated SXA. Based these results, we suggest that choosing a correct residue for mutation in multi subunit proteins results in multiple site conversions which equals to several simultaneous mutations. Furthermore, CG-MD simulation in agreement with experimental methods showed higher thermostability of mutated SXA which proved applicability of this simulation for thermostability analysis.

Hyperstabilization of Tetrameric Bacillus sp. TB-90 Urate Oxidase by Introducing Disulfide Bonds through Structural Plasticity

Bacillus sp. TB-90 urate oxidase (BTUO) is one of the most thermostable homotetrameric enzymes. We previously reported [Hibi, T., et al. (2014) Biochemistry 53, 3879-3888] that specific binding of a sulfate anion induced thermostabilization of the enzyme, because the bound sulfate formed a salt bridge with two Arg298 residues, which stabilized the packing between two β-barrel dimers. To extensively characterize the sulfate-binding site, Arg298 was substituted with cysteine by site-directed mutagenesis. This substitution markedly increased the protein melting temperature by ∼ 20 °C compared with that of the wild-type enzyme, which was canceled by reduction with dithiothreitol. Calorimetric analysis of the thermal denaturation suggested that the hyperstabilization resulted from suppression of the dissociation of the tetramer into the two homodimers. The crystal structure of R298C at 2.05 Å resolution revealed distinct disulfide bond formation between the symmetrically related subunits via Cys298, although the Cβ distance between Arg298 residues of the wild-type enzyme (5.4 Å apart) was too large to predict stable formation of an engineered disulfide cross-link. Disulfide bonding was associated with local disordering of interface loop II (residues 277-300), which suggested that the structural plasticity of the loop allowed hyperstabilization by disulfide formation. Another conformational change in the C-terminal region led to intersubunit hydrogen bonding between Arg7 and Asp312, which probably promoted mutant thermostability. Knowledge of the disulfide linkage of flexible loops at the subunit interface will help in the development of new strategies for enhancing the thermostabilization of multimeric proteins.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

Additional disulfide bonds in insulin: Prediction, recombinant expression, receptor binding affinity, and stability

The structure of insulin, a glucose homeostasis-controlling hormone, is highly conserved in all vertebrates and stabilized by three disulfide bonds. Recently, we designed a novel insulin analogue containing a fourth disulfide bond located between positions A10-B4. The N-terminus of insulin's B-chain is flexible and can adapt multiple conformations. We examined how well disulfide bond predictions algorithms could identify disulfide bonds in this region of insulin. In order to identify stable insulin analogues with additional disulfide bonds, which could be expressed, the Cβ cut-off distance had to be increased in many instances and single X-ray structures as well as structures from MD simulations had to be used. The analogues that were identified by the algorithm without extensive adjustments of the prediction parameters were more thermally stable as assessed by DSC and CD and expressed in higher yields in comparison to analogues with additional disulfide bonds that were more difficult to predict. In contrast, addition of the fourth disulfide bond rendered all analogues resistant to fibrillation under stress conditions and all stable analogues bound to the insulin receptor with picomolar affinities. Thus activity and fibrillation propensity did not correlate with the results from the prediction algorithm. Statement: A fourth disulfide bond has recently been introduced into insulin, a small two-chain protein containing three native disulfide bonds. Here we show that a prediction algorithm predicts four additional four disulfide insulin analogues which could be expressed. Although the location of the additional disulfide bonds is only slightly shifted, this shift impacts both stability and activity of the resulting insulin analogues.

Rational design, preparation, and characterization of a therapeutic enzyme mutant with improved stability and function for cocaine detoxification

Cocaine esterase (CocE) is known as the most efficient natural enzyme for cocaine hydrolysis. The major obstacle to the clinical application of wild-type CocE is the thermoinstability with a half-life of only ∼12 min at 37 °C. The previously designed T172R/G173Q mutant (denoted as enzyme E172-173) with an improved in vitro half-life of ∼6 h at 37 °C is currently in clinical trial Phase II for cocaine overdose treatment. Through molecular modeling and dynamics simulation, we designed and characterized a promising new mutant of E172-173 with extra L196C/I301C mutations (denoted as enzyme E196-301) to produce cross-subunit disulfide bonds that stabilize the dimer structure. The cross-subunit disulfide bonds were confirmed by X-ray diffraction. The designed L196C/I301C mutations have not only considerably extended the in vitro half-life at 37 °C to >100 days, but also significantly improved the catalytic efficiency against cocaine by ∼150%. In addition, the thermostable E196-301 can be PEGylated to significantly prolong the residence time in mice. The PEGylated E196-301 can fully protect mice from a lethal dose of cocaine (180 mg/kg, LD100) for at least 3 days, with an average protection time of ∼94h. This is the longest in vivo protection of mice from the lethal dose of cocaine demonstrated within all studies using an exogenous enzyme reported so far. Hence, E196-301 may be developed to become a more valuable therapeutic enzyme for cocaine abuse treatment, and it demonstrates that a general design strategy and protocol to simultaneously improve both the stability and function are feasible for rational protein drug design.

Novel covalently linked insulin dimer engineered to investigate the function of insulin dimerization

An ingenious system evolved to facilitate insulin binding to the insulin receptor as a monomer and at the same time ensure sufficient stability of insulin during storage. Insulin dimer is the cornerstone of this system. Insulin dimer is relatively weak, which ensures dissociation into monomers in the circulation, and it is stabilized by hexamer formation in the presence of zinc ions during storage in the pancreatic β-cell. Due to the transient nature of insulin dimer, direct investigation of this important form is inherently difficult. To address the relationship between insulin oligomerization and insulin stability and function, we engineered a covalently linked insulin dimer in which two monomers were linked by a disulfide bond. The structure of this covalent dimer was identical to the self-association dimer of human insulin. Importantly, this covalent dimer was capable of further oligomerization to form the structural equivalent of the classical hexamer. The covalently linked dimer neither bound to the insulin receptor, nor induced a metabolic response in vitro. However, it was extremely thermodynamically stable and did not form amyloid fibrils when subjected to mechanical stress, underlining the importance of oligomerization for insulin stability.

Subunit stabilization and polyethylene glycolation of cocaine esterase improves in vivo residence time

No small-molecule therapeutic is available to treat cocaine addiction, but enzyme-based therapy to accelerate cocaine hydrolysis in serum has gained momentum. Bacterial cocaine esterase (CocE) is the fastest known native enzyme that hydrolyzes cocaine. However, its lability at 37°C has limited its therapeutic potential. Cross-linking subunits through disulfide bridging is commonly used to stabilize multimeric enzymes. Herein we use structural methods to guide the introduction of two cysteine residues within dimer interface of CocE to facilitate intermolecular disulfide bond formation. The disulfide-crosslinked enzyme displays improved thermostability, particularly when combined with previously described mutations that enhance stability (T172R-G173Q). The newly modified enzyme yielded an extremely stable form of CocE (CCRQ-CocE) that retained greater than 90% of its activity after 41 days at 37°C, representing an improvement of more than 4700-fold over the wild-type enzyme. CCRQ-CocE could also be modified by polyethylene glycol (PEG) polymers, which improved its in vivo residence time from 24 to 72 h, as measured by a cocaine lethality assay, by self-administration in rodents, and by measurement of inhibition of cocaine-induced cardiovascular effects in rhesus monkeys. PEG-CCRQ elicited negligible immune response in rodents. Subunit stabilization and PEGylation has thus produced a potential protein therapeutic with markedly higher stability both in vitro and in vivo.

Engineered disulfide bonds restore chaperone-like function of DJ-1 mutants linked to familial Parkinson's disease

Loss-of-function mutations such as L166P, A104T, and M26I in the DJ-1 gene (PARK7) have been linked to autosomal-recessive early onset Parkinson's disease (PD). Cellular and structural studies of the familial mutants suggest that these mutations may destabilize the dimeric structure. To look for common dynamical signatures among the DJ-1 mutants, short MD simulations of up to 1000 ps were conducted to identify the weakest region of the protein (residues 38-70). In an attempt to stabilize the protein, we mutated residue Val 51 to cysteine (V51C) to make a symmetry-related disulfide bridge with the preexisting Cys 53 on the opposite subunit. We found that the introduction of this disulfide linkage stabilized the mutants A104T and M26I against thermal denaturation, improved their ability to scavenge reactive oxygen species (ROS), and restored a chaperone-like function of blocking alpha-synuclein aggregation. The L166P mutant was far too unstable to be rescued by introduction of the V51C mutation. The results presented here point to the possible development of pharmacological chaperones, which may eventually lead to PD therapeutics.

Novel vascular endothelial growth factor D variants with increased biological activity

Members of the vascular endothelial growth factor (VEGF) family play a pivotal role in angiogenesis and lymphangiogenesis. They are potential therapeutics to induce blood vessel formation in myocardium and skeletal muscle, when normal blood flow is compromised. Most members of the VEGF/platelet derived growth factor protein superfamily exist as covalently bound antiparallel dimers. However, the mature form of VEGF-D (VEGF-D(DeltaNDeltaC)) is predominantly a non-covalent dimer even though the cysteine residues (Cys-44 and Cys-53) forming the intersubunit disulfide bridges in the other members of the VEGF family are also conserved in VEGF-D. Moreover, VEGF-D bears an additional cysteine residue (Cys-25) at the subunit interface. Guided by our model of VEGF-D(DeltaNDeltaC), the cysteines at the subunit interface were mutated to study the effect of these residues on the structural and functional properties of VEGF-D(DeltaNDeltaC). The conserved cysteines Cys-44 and Cys-53 were found to be essential for the function of VEGF-D(DeltaNDeltaC). More importantly, the substitution of the Cys-25 at the dimer interface by various amino acids improved the activity of the recombinant VEGF-D(DeltaNDeltaC) and increased the dimer to monomer ratio. Specifically, substitutions to hydrophobic amino acids Ile, Leu, and Val, equivalent to those found in other VEGFs, most favorably affected the activity of the recombinant VEGF-D(DeltaNDeltaC). The increased activity of these mutants was mainly due to stabilization of the protein. This study enables us to better understand the structural determinants controlling the biological activity of VEGF-D. The novel variants of VEGF-D(DeltaNDeltaC) described here are potential agents for therapeutic applications, where induction of vascular formation is required.

Thermal stabilization of the protozoan Entamoeba histolytica alcohol dehydrogenase by a single proline substitution

Analysis of the three-dimensional structures of two closely related thermophilic and hyperthermophilic alcohol dehydrogenases (ADHs) from the respective microorganisms Entamoeba histolytica (EhADH1) and Thermoanaerobacter brockii (TbADH) suggested that a unique, strategically located proline residue (Pro275) at the center of the dimerization interface might be crucial for maintaining the thermal stability of TbADH. To assess the contribution of Pro275 to the thermal stability of the ADHs, we applied site-directed mutagenesis to replace Asp275 of EhADH1 with Pro (D275P-EhADH1) and conversely Pro275 of TbADH with Asp (P275D-TbADH). The results indicate that replacing Asp275 with Pro significantly enhances the thermal stability of EhADH1 (DeltaT(1/2) <or= +10 degrees C), whereas the reverse mutation in the thermophilic TbADH (P275D-TbADH) reduces the thermostability of the enzyme (DeltaT(1/2) <or= -18.8 degrees C). Analysis of the crystal structures of the thermostabilized mutant D275P-EhADH1 and the thermocompromised mutant P275D-TbADH suggest that a proline residue at position 275 thermostabilized the enzymes by reducing flexibility and by reinforcing hydrophobic interactions at the dimer-dimer interface of the tetrameric ADHs.

Design of Disulfide-linked Thioredoxin Dimers and Multimers Through Analysis of Crystal Contacts

Disulfide bonds play an important role in protein stability and function. Here, we describe a general procedure for generating disulfide-linked dimers and multimers of proteins of known crystal structures. An algorithm was developed to predict sites in a protein compatible with intermolecular disulfide formation with neighboring molecules in the crystal lattice. A database analysis was carried out on 46 PDB coordinates to verify the general applicability of this algorithm to predict intermolecular disulfide linkages. On the basis of the predictions from this algorithm, mutants were constructed and characterized for a model protein, thioredoxin. Of the five mutants, as predicted, in solution four formed disulfide-linked dimers while one formed polymers. Thermal and chemical denaturation studies on these mutant thioredoxins showed that three of the four dimeric mutants had similar stability to wild-type thioredoxin while one had lower stability. Three of the mutant dimers crystallized readily (in four to seven days) in contrast to the wild-type protein, which is particularly difficult to crystallize and takes more than a month to form diffraction-quality crystals. In two of the three cases, the structure of the dimer was exactly as predicted by the algorithm, while in the third case the relative orientation of the monomers in the dimer was different from the predicted one. This methodology can be used to enhance protein crystallizability, modulate the oligomerization state and to produce linear chains or ordered three-dimensional protein arrays.

Structural Stability of Covalently Linked GroES Heptamer: Advantages in the Formation of Oligomeric Structure

In order to understand how inter-subunit association stabilizes oligomeric proteins, a single polypeptide chain variant of heptameric co-chaperonin GroES (tandem GroES) was constructed from Escherichia coli heptameric GroES by linking consecutively the C-terminal of one subunit to the N-terminal of the adjacent subunit with a small linker peptide. The tandem GroES (ESC7) showed properties similar to wild-type GroES in structural aspects and co-chaperonin activity. In unfolding and refolding equilibrium experiments using guanidine hydrochloride (Gdn-HCl) as a denaturant at a low protein concentration (50 microg ml(-1)), ESC7 showed a two-state transition with a greater resistance toward Gdn-HCl denaturation (Cm=1.95 M) compared to wild-type GroES (Cm=1.1 M). ESC7 was found to be about 10 kcal mol(-1) more stable than the wild-type GroES heptamer at 50 microg ml(-1). Kinetic unfolding and refolding experiments of ESC7 revealed that the increased stability was mainly attributed to a slower unfolding rate. Also a transient intermediate was detected in the refolding reaction. Interestingly, at the physiological GroES concentration (>1 mg ml(-1)), the free energy of unfolding for GroES heptamer exceeded that for ESC7. These results showed that at low protein concentrations (<1 mg ml(-1)), the covalent linking of subunits contributes to the stability but also complicates the refolding kinetics. At physiological concentrations of GroES, however, the oligomeric state is energetically preferred and the advantages of covalent linkage are lost. This finding highlights a possible advantage in transitioning from multi-domain proteins to oligomeric proteins with small subunits in order to improve structural and kinetic stabilities.

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.

Protein disulfides and protein disulfide oxidoreductases in hyperthermophiles

Disulfide bonds are required for the stability and function of a large number of proteins. Recently, the results from genome analysis have suggested an important role for disulfide bonds concerning the structural stabilization of intracellular proteins from hyperthermophilic Archaea and Bacteria, contrary to the conventional view that structural disulfide bonds are rare in proteins from Archaea. A specific protein, known as protein disulfide oxidoreductase (PDO) is recognized as a potential key player in intracellular disulfide-shuffling in hyperthermophiles. The structure of this protein shows a combination of two thioredoxin-related units with low sequence identity which together, in tandem-like manner, form a closed protein domain. Each of these units contains a distinct CXXC active site motif. Due to their estimated conformational energies, both sites are likely to have different redox properties. The observed structural and functional characteristics suggest a relation to eukaryotic protein disulfide isomerase. Functional studies have revealed that both the archaeal and bacterial forms of this protein show oxidative and reductive activity and are able to isomerize protein disulfides. The physiological substrates and reduction systems, however, are to date unknown. The variety of active site disulfides found in PDOs from hyperthermophiles is puzzling. Nevertheless, the catalytic function of any PDO is expected to be correlated with the redox properties of its active site disulfides CXXC and with the distinct nature of its redox environment. The residues around the two active sites form two grooves on the protein surface. In analogy to a similar groove in thioredoxin, both grooves are suggested to constitute the substrate binding sites of PDO. The direct neighbourhood of the grooves and the different redox properties of both sites may favour sequential reactions in protein disulfide shuffling, like reduction followed by oxidation. A model for peptide binding by PDO is proposed to be derived from the analysis of crystal packing contacts mimicking substrate binding interactions. It is assumed, that PDO enzymes in hyperthermophilic Archaea and Bacteria may be part of a complex system involved in the maintenance of protein disulfide bonds. The regulation of disulfide bond formation may be dependent on a distinct interplay of thermodynamic and kinetic effects, including functional asymmetry and substrate-mediated protection of the active sites, in analogy to the situation in protein disulfide isomerase. Numerous questions related to the function of PDO enzymes in hyperthermophiles remain unanswered to date, but can probably successfully be studied by a number of approaches, such as first-line genetic and in vivo studies.

Structure-stability-activity relationship in covalently cross-linked N-carbamoyl D-amino acid amidohydrolase and N-acylamino acid racemase

N-Acylamino acid racemase (NAAAR) and N-carbamoyl-D-amino-acid amidohydrolase (D-NCAase) are important biocatalysts for producing enantiopure alpha-amino acids. NAAAR forms an octameric assembly and displays induced fit movements upon substrate binding, while D-NCAase is a tetramer that does not change conformation in the presence of a ligand. To investigate the effects of introducing potentially stabilizing S-S bridges in these different multimeric enzymes, cysteine residues predicted to form inter or intra-subunit disulfide bonds were introduced by site-directed mutagenesis. Inter-subunit S-S bonds were formed in two NAAAR variants (A68C-D72C and P60C-Y100C) and two d-NCAase variants (A302C and P295C-F304C). Intra-subunit S-S bonds were formed in two additional NAAAR variants (E149C-A182C and V265C). Crystal structures of NAAARs variants show limited deviations from the wild-type overall tertiary structure. An apo A68C-D72C subunit differs from the wild-type enzyme, in which it has an ordered lid loop, resembling ligand-bound NAAAR. The structures of A222C and A302C D-NCAases are nearly identical to the wild-type enzyme. All mutants with inter-subunit bridges had increases in thermostability. Compared with the wild-type enzyme, A68C-D72C NAAAR showed similar kcat/Km ratios, whereas mutant D-NCAases demonstrated increased kcat/Km ratios at high temperatures (A302C: 4.2-fold at 65 degrees C). Furthermore, molecular dynamic simulations reveal that A302C substantially sustains the fine-tuned catalytic site as temperature increases, achieving enhanced activity.

Directed evolution of enzyme stability

Modern enzyme development relies to an increasing extent on strategies based on diversity generation followed by screening for variants with optimised properties. In principle, these directed evolution strategies might be used for optimising any enzyme property, which can be screened for in an economically feasible way, even if the molecular basis of that property is not known. Stability is an interesting property of enzymes because (1) it is of great industrial importance, (2) it is relatively easy to screen for, and (3) the molecular basis of stability relates closely to contemporary issues in protein science such as the protein folding problem and protein folding diseases. Thus, engineering enzyme stability is of both commercial and scientific interest. Here, we review how directed evolution has contributed to the development of stable enzymes and to new insight into the principles of protein stability. Several recent examples are described. These examples show that directed evolution is an effective strategy to obtain stable enzymes, especially when used in combination with rational or semi-rational engineering strategies. With respect to the principles of protein stability, some important lessons to learn from recent efforts in directed evolution are (1) that there are many structural ways to stabilize a protein, which are not always easy to rationalize, (2) that proteins may very well be stabilized by optimizing their surfaces, and (3) that high thermal stability may be obtained without forfeiture of catalytic performance at low temperatures.

Amylomaltase of Pyrobaculum aerophilum IM2 produces thermoreversible starch gels

Amylomaltases are 4-alpha-glucanotransferases (EC 2.4.1.25) of glycoside hydrolase family 77 that transfer alpha-1,4-linked glucans to another acceptor, which can be the 4-OH group of an alpha-1,4-linked glucan or glucose. The amylomaltase-encoding gene (PAE1209) from the hyperthermophilic archaeon Pyrobaculum aerophilum IM2 was cloned and expressed in Escherichia coli, and the gene product (PyAMase) was characterized. PyAMase displays optimal activity at pH 6.7 and 95 degrees C and is the most thermostable amylomaltase described to date. The thermostability of PyAMase was reduced in the presence of 2 mM dithiothreitol, which agreed with the identification of two possible cysteine disulfide bridges in a three-dimensional model of PyAMase. The kinetics for the disproportionation of malto-oligosaccharides, inhibition by acarbose, and binding mode of the substrates in the active site were determined. Acting on gelatinized food-grade potato starch, PyAMase produced a thermoreversible starch product with gelatin-like properties. This thermoreversible gel has potential applications in the food industry. This is the first report on an archaeal amylomaltase.

Rational engineering of enzyme stability

During the past 15 years there has been a continuous flow of reports describing proteins stabilized by the introduction of mutations. These reports span a period from pioneering rational design work on small enzymes such as T4 lysozyme and barnase to protein design, and directed evolution. Concomitantly, the purification and characterization of naturally occurring hyperstable proteins has added to our understanding of protein stability. Along the way, many strategies for rational protein stabilization have been proposed, some of which (e.g. entropic stabilization by introduction of prolines or disulfide bridges) have reasonable success rates. On the other hand, comparative studies and efforts in directed evolution have revealed that there are many mutational strategies that lead to high stability, some of which are not easy to define and rationalize. Recent developments in the field include increasing awareness of the importance of the protein surface for stability, as well as the notion that normally a very limited number of mutations can yield a large increase in stability. Another development concerns the notion that there is a fundamental difference between the "laboratory stability" of small pure proteins that unfold reversibly and completely at high temperatures and "industrial stability", which is usually governed by partial unfolding processes followed by some kind of irreversible inactivation process (e.g. aggregation). Provided that one has sufficient knowledge of the mechanism of thermal inactivation, successful and efficient rational stabilization of enzymes can be achieved.

Large improvement in the thermal stability of a tetrameric malate dehydrogenase by single point mutations at the dimer-dimer interface

The stability of tetrameric malate dehydrogenase from the green phototrophic bacterium Chloroflexus aurantiacus (CaMDH) is at least in part determined by electrostatic interactions at the dimer-dimer interface. Since previous studies had indicated that the thermal stability of CaMDH becomes lower with increasing pH, attempts were made to increase the stability by removal of (excess) negative charge at the dimer-dimer interface. Mutation of Glu165 to Gln or Lys yielded a dramatic increase in thermal stability at pH 7.5 (+23.6 -- + 23.9 degrees C increase in apparent t(m)) and a more moderate increase at pH 4.4 (+4.6 -- + 5.4 degrees C). The drastically increased stability at neutral pH was achieved without forfeiture of catalytic performance at low temperatures. The crystal structures of the two mutants showed only minor structural changes close to the mutated residues, and indicated that the observed stability effects are solely due to subtle changes in the complex network of electrostatic interactions in the dimer-dimer interface. Both mutations reduced the concentration dependency of thermal stability, suggesting that the oligomeric structure had been reinforced. Interestingly, the two mutations had similar effects on stability, despite the charge difference between the introduced side-chains. Together with the loss of concentration dependency, this may indicate that both E165Q and E165K stabilize CaMDH to such an extent that disruption of the inter-dimer electrostatic network around residue 165 no longer limits kinetic thermal stability.