Electrostatic contributions to the stability of hyperthermophilic proteins

Rapid comparison of properties on protein surface

The mapping of physicochemical characteristics onto the surface of a protein provides crucial insights into its function and evolution. This information can be further used in the characterization and identification of similarities within protein surface regions. We propose a novel method which quantitatively compares global and local properties on the protein surface. We have tested the method on comparison of electrostatic potentials and hydrophobicity. The method is based on 3D Zernike descriptors, which provides a compact representation of a given property defined on a protein surface. Compactness and rotational invariance of this descriptor enable fast comparison suitable for database searches. The usefulness of this method is exemplified by studying several protein families including globins, thermophilic and mesophilic proteins, and active sites of TIM beta/alpha barrel proteins. In all the cases studied, the descriptor is able to cluster proteins into functionally relevant groups. The proposed approach can also be easily extended to other surface properties. This protein surface-based approach will add a new way of viewing and comparing proteins to conventional methods, which compare proteins in terms of their primary sequence or tertiary structure.

A computational approach for the rational design of stable proteins and enzymes: optimization of surface charge-charge interactions

The design of stable proteins and enzymes is not only of particular biotechnological importance, but also addresses some important fundamental questions. While there are a number of different options available for designing or engineering stable proteins, the field of computational design provides fast and universal methods for stabilizing proteins of interest. One of the successful computational design strategies focuses on stabilizing proteins through the optimization of charge-charge interactions on the protein surface. By optimizing surface interactions, it is possible to alleviate some of the challenges that accompany efforts to redesign the protein core. The rational design of surface charge-charge interactions also allows one to optimize only the interactions that are distant from binding sites or active sites, making it possible to increase stability without adversely affecting activity. The optimization of surface charge-charge interactions is discussed in detail along with the experimental evidence to demonstrate that this is a robust and universal approach to designing proteins with enhanced stability.

Characteristics of mutants designed to incorporate a new ion pair into the structure of a cold adapted subtilisin-like serine proteinase

Structural comparisons of VPR, a subtilisin-like serine proteinase from a psychrotrophic Vibrio species and a thermophilic homologue, aqualysin I, have led us to hypothesize about the roles of different residues in the temperature adaptation of the enzymes. Some of these hypotheses are now being examined by analysis of mutants of the enzymes. The selected substitutions are believed to increase the stability of the cold adapted enzyme based on structural analysis of the thermostable structure. We report here on mutants, which were designed to incorporate an ion pair into the structure of VPR. The residues Asp17 and Arg259 are assumed to form an ion pair in aqualysin I. The cold adapted VPR contains Asn (Asn15) and Lys (Lys257) at corresponding sites in its structure. In VPR, Asn 15 is located on a surface loop with its side group pointing towards the side chain of Lys257. By substituting Asn15 by Asp (N15D) it was considered feasible that a salt bridge would form between the oppositely charged groups. To mimic further the putative salt bridge from the thermophile enzyme the corresponding double mutant (N15D/K257R) was also produced. The N15D mutation increased the thermal stability of VPR by approximately 3 degrees C, both in T(50%) and T(m). Addition of the K257R mutation did not however, increase the stability of the double mutant any further. Despite this stabilization of the VPR mutants the catalytic activity (k(cat)) against the substrate Suc-AAPF-NH-Np was increased in the mutants. Molecular dynamics simulations on wild type and the two mutant proteins suggested that indeed a salt bridge was formed in both cases. Furthermore, a truncated form of the N15D mutant (N15DDeltaC) was produced, lacking a 15 residue long C-terminal extended sequence not present in the thermophilic enzyme. In wild type VPR this supposedly moveable, negatively charged arm on the protein molecule might interfere with the new salt bridge introduced as a result of the N15D mutation. Removal of the C-terminal arm improved the thermal stability (T(m) approximately +1.5 degrees C) of the truncated enzyme (VPRDeltaC) as compared to the wild type VPR. Introduction of the N15D substitution into VPRDeltaC improved the thermal stability further by about 3 degrees C, or to about the same extent as in the wild type. However, contrary to what was observed for the wild type, the introduction of the putative salt bridge did not affect the catalytic properties (k(cat)) of the C-terminal truncated enzyme.

Kinetic and thermodynamic characterization of dUTP hydrolysis by Plasmodium falciparum dUTPase

Deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase) catalyzes the hydrolysis of dUTP to dUMP and pyrophosphate and plays an important role in nucleotide metabolism and DNA replication controlling relative cellular levels of dTTP/dUTP, both of which can be incorporated into DNA. Isothermal titration calorimetry has been applied to the determination of the kinetic and thermodynamic parameters of the trimeric Plasmodium falciparum dUTPase, a potential drug target against malaria. The role of divalent ions in binding, and inhibition by different uridine derivatives has been assessed. When dUTP hydrolysis in the presence of EDTA was evaluated, a 105-fold decrease and a 12-fold increase of the k(cat) and Km values, respectively, were observed when compared with the dUTP.Mg2+ complex. Calculation of the activation energy, E(a), and the thermodynamic activation parameters showed that the energetic barrier was approximately 4-fold higher when Mg2+ was depleted. Other divalent ions such as Co2+ or Mn2+ can substitute the physiological cofactor, however the k(cat) was significantly reduced compared to dUTP.Mg2+. Binding and inhibition by dU, dUMP, dUDP, and alpha,beta-imido-dUTP were analysed by ITC and compared with data obtained by spectrophotometric methods and binding equilibrium studies. Product inhibition (Kip dUMP: 99.34 microM) was insignificant yet Ki values for dUDP and alpha,beta-imido-dUTP were in the low micromolar range. The effect of ionic strength on protein stability was also monitored. DSC analysis evidenced a slight increase in the unfolding temperature, Tm, with increasing salt concentrations. Moreover, the thermal unfolding pathway in the presence of salt fits adequately to an irreversible two-state model (N3-->3D).

Engineering hyperthermostability into a GH11 xylanase is mediated by subtle changes to protein structure

Understanding the structural basis for protein thermostability is of considerable biological and biotechnological importance as exemplified by the industrial use of xylanases at elevated temperatures in the paper pulp and animal feed sectors. Here we have used directed protein evolution to generate hyperthermostable variants of a thermophilic GH11 xylanase, EvXyn11. The Gene Site Saturation Mutagenesis (GSSM) methodology employed assesses the influence on thermostability of all possible amino acid substitutions at each position in the primary structure of the target protein. The 15 most thermostable mutants, which generally clustered in the N-terminal region of the enzyme, had melting temperatures (Tm) 1-8 degrees C higher than the parent protein. Screening of a combinatorial library of the single mutants identified a hyperthermostable variant, EvXyn11TS, containing seven mutations. EvXyn11TS had a Tm approximately 25 degrees C higher than the parent enzyme while displaying catalytic properties that were similar to EvXyn11. The crystal structures of EvXyn11 and EvXyn11TS revealed an absence of substantial changes to identifiable intramolecular interactions. The only explicable mutations are T13F, which increases hydrophobic interactions, and S9P that apparently locks the conformation of a surface loop. This report shows that the molecular basis for the increased thermostability is extraordinarily subtle and points to the requirement for new tools to interrogate protein folding at non-ambient temperatures.

Evaluation of electrostatic interactions

Described here are several computational procedures for the analysis of electrostatic interactions in molecular complexes, all based on a continuum model of solvation. The first section describes how to compute the residual potential, a measure of how electrostatically complementary a ligand is for its receptor. The second procedure describes electrostatic component analysis, a method by which the electrostatic contribution to the binding free energy can be broken up into terms directly attributable to individual chemical groups. Finally, electrostatic affinity optimization is described. This procedure is particularly useful in determining what portions of a ligand are the most suboptimal, and thus provide the greatest opportunity for the design of improvements.

Cooperativity of complex salt bridges

The energetic contribution of complex salt bridges, in which one charged residue (anchor residue) forms salt bridges with two or more residues simultaneously, has been suggested to have importance for protein stability. Detailed analysis of the net energetics of complex salt bridge formation using double- and triple-mutant cycle analysis revealed conflicting results. In two cases, it was shown that complex salt bridge formation is cooperative, i.e., the net strength of the complex salt bridge is more than the sum of the energies of individual pairs. In one case, it was reported that complex salt bridge formation is anti-cooperative. To resolve these different findings, we performed analysis of the geometries of salt bridges in a representative set of structures from the PDB and found that over 87% of all complex salt bridges anchored by Arg/Lys have a geometry such that the angle formed by their Calpha atoms, Theta, is <90 degrees . This preferred geometry is observed in the two reported instances when the energetics of complex salt bridge formation is cooperative, while in the reported anti-cooperative complex salt bridge, Theta is close to 160 degrees . Based on these observations, we hypothesized that complex salt bridges are cooperative for Theta < 90 degrees and anti-cooperative for 90 degrees < Theta < 180 degrees . To provide a further experimental test for this hypothesis, we engineered a complex salt bridge with Theta = 150 degrees into a model protein, the activation domain of human procarboxypeptidase A2 (ADA2h). Experimentally derived stabilities of the ADA2h variants allowed us to show that the complex salt bridge in ADA2h is anti-cooperative.

Electrostatic effects in a network of polar and ionizable groups in staphylococcal nuclease

His121 and His124 are embedded in a network of polar and ionizable groups on the surface of staphylococcal nuclease. To examine how membership in a network affects the electrostatic properties of ionizable groups, the tautomeric state and the pK(a) values of these histidines were measured with NMR spectroscopy in the wild-type nuclease and in 13 variants designed to disrupt the network. In the background protein, His121 and His124 titrate with pK(a) values of 5.2 and 5.6, respectively. In the variants, where the network was disrupted, the pK(a) values range from 4.03 to 6.46 for His121, and 5.04 to 5.99 for His124. The largest decrease in a pK(a) was observed when the favorable Coulomb interaction between His121 and Glu75 was eliminated; the largest increase was observed when Tyr91 or Tyr93 was substituted with Ala or Phe. In all variants, the dominant tautomeric state at neutral pH was the N(epsilon2) state. At one level the network behaves as a rigid unit that does not readily reorganize when disrupted: crystal structures of the E75A or E75Q variants show that even when the pivotal Glu75 is removed, the overall configuration of the network was unaffected. On the other hand, a few key hydrogen bonds appear to govern the conformation of the network, and when these bonds are disrupted the network reorganizes. Coulomb interactions within the network report an effective dielectric constant of 20, whereas a dielectric constant of 80 is more consistent with the magnitude of medium to long-range Coulomb interactions in this protein. The data demonstrate that when structures are treated as static, rigid bodies, structure-based pK(a) calculations with continuum electrostatics method are not useful to treat ionizable groups in cases where pK(a) values are governed by short-range polar and Coulomb interactions.

Discrimination of mesophilic and thermophilic proteins using machine learning algorithms

Discriminating thermophilic proteins from their mesophilic counterparts is a challenging task and it would help to design stable proteins. In this work, we have systematically analyzed the amino acid compositions of 3075 mesophilic and 1609 thermophilic proteins belonging to 9 and 15 families, respectively. We found that the charged residues Lys, Arg, and Glu as well as the hydrophobic residues, Val and Ile have higher occurrence in thermophiles than mesophiles. Further, we have analyzed the performance of different methods, based on Bayes rules, logistic functions, neural networks, support vector machines, decision trees and so forth for discriminating mesophilic and thermophilic proteins. We found that most of the machine learning techniques discriminate these classes of proteins with similar accuracy. The neural network-based method could discriminate the thermophiles from mesophiles at the five-fold cross-validation accuracy of 89% in a dataset of 4684 proteins. Moreover, this method is tested with 325 mesophiles in Xylella fastidosa and 382 thermophiles in Aquifex aeolicus and it could successfully discriminate them with the accuracy of 91%. These accuracy levels are better than other methods in the literature and we suggest that this method could be effectively used to discriminate mesophilic and thermophilic proteins.

Contribution of charged groups to the enthalpic stabilization of the folded states of globular proteins

In this paper we use the results from all-atom molecular dynamics (MD) simulations of proteins and peptides to assess the individual contribution of charged atomic groups to the enthalpic stability of the native state of globular proteins and investigate how the distribution of charged atomic groups in terms of solvent accessibility relates to protein enthalpic stability. The contributions of charged groups is calculated using a comparison of nonbonded interaction energy terms from equilibrium simulations of charged amino acid dipeptides in water (the "unfolded state") and charged amino acids in globular proteins (the "folded state"). Contrary to expectation, the analysis shows that many buried, charged atomic groups contribute favorably to protein enthalpic stability. The strongest enthalpic contributions favoring the folded state come from the carboxylate (COO(-)) groups of either Glu or Asp. The contributions from Arg guanidinium groups are generally somewhat stabilizing, while N(+)(3) groups from Lys contribute little toward stabilizing the folded state. The average enthalpic gain due to the transfer of a methyl group in an apolar amino acid from solution to the protein interior is described for comparison. Notably, charged groups that are less exposed to solvent contribute more favorably to protein native-state enthalpic stability than charged groups that are solvent exposed. While solvent reorganization/release has favorable contributions to folding for all charged atomic groups, the variation in folded state stability among proteins comes mainly from the change in the nonbonded interaction energy of charged groups between the unfolded and folded states. A key outcome is that the calculated enthalpic stabilization is found to be inversely proportional to the excess charge density on the surface, in support of an hypothesis proposed previously.

Thermostability of salt bridges versus hydrophobic interactions in proteins probed by statistical potentials

The temperature dependence of the interactions that stabilize protein structures is a long-standing issue, the elucidation of which would enable the prediction and the rational modification of the thermostability of a target protein. It is tackled here by deriving distance-dependent amino acid pair potentials from four datasets of proteins with increasing melting temperatures (Tm). The temperature dependence of the interactions is determined from the differences in the shape of the potentials derived from the four datasets. Note that, here, we use an unusual dataset definition, which is based on the Tm values, rather than on the living temperature of the host organisms. Our results show that the stabilizing weight of hydrophobic interactions (between Ile, Leu, and Val) remains constant as the temperature increases, compared to the other interactions. In contrast, the two minima of the Arg--Glu and Arg--Asp salt bridge potentials show a significant Tm dependence. These two minima correspond to two geometries: the fork--fork geometry, where the side chains point toward each other, and the fork--stick geometry, which involves the N(epsilon) side chain atom of Arg. These two types of salt bridges were determined to be significantly more stabilizing at high temperature. Moreover, a preference for more-compact salt bridges is noticeable in heat-resistant proteins, especially for the fork--fork geometry. The Tm-dependent potentials that have been defined here should be useful for predicting thermal stability changes upon mutation.

Amino-acid interactions in psychrophiles, mesophiles, thermophiles, and hyperthermophiles: insights from the quasi-chemical approximation

We investigate the mechanisms used by proteins to maintain thermostability throughout a wide range of temperatures. We use the quasi-chemical approximation to estimate interaction strengths for psychrophiles, mesophiles, thermophiles, and hyperthermophiles. Our results highlight the importance of core packing in thermophilic stability. Although we observed an increase in the number of charged residues, the contribution of salt bridges appears to be relatively modest by comparison. We observed results consistent with a gradual loosening of structure in psychrophiles, including a weakening of almost all types of interactions.

Enhancing the thermal stability of mitochondrial cytochrome b5 by introducing a structural motif characteristic of the less stable microsomal isoform

Outer mitochondrial membrane cytochrome b5 (OM b5) is the most thermostable cytochrome b5 isoform presently known. Herein, we show that OM b5 thermal stability is substantially enhanced by swapping an apparently invariant motif in its heme-independent folding core with the corresponding motif characteristic of its less stable evolutionary relative, microsomal cytochrome b5 (Mc b5). The motif swap involved replacing two residues, Arg15 with His and Glu20 with Ser, thereby introducing a Glu11-His15-Ser20 H-bonding triad on the protein surface along with a His15/Trp22 pi-stacking interaction. The ferric and ferrous forms of the OM b5 R15H/E20S double mutant have thermal denaturation midpoints (Tm values) of approximately 93 degrees C and approximately 104 degrees C, respectively. A 15 degrees C increase in apoprotein Tm plays a key role in the holoprotein thermal stability enhancement, and is achieved by one of the most common natural mechanisms for stabilization of thermophilic versus mesophilic proteins: raising the unfolding free energy along the entire stability curve.

Carboxyl pK(a) values, ion pairs, hydrogen bonding, and the pH-dependence of folding the hyperthermophile proteins Sac7d and Sso7d

Sac7d and Sso7d are homologous, hyperthermophile proteins with a high density of charged surface residues and potential ion pairs. To determine the relative importance of specific amino acid side-chains in defining the stability and function of these Archaeal chromatin proteins, pK(a) values were measured for the acidic residues in both proteins using (13)C NMR chemical shifts. The stability of Sso7d enabled titrations to pH 1 under low-salt conditions. Two aspartate residues in Sso7d (D16 and D35) and a single glutamate residue (G54) showed significantly perturbed pK(a) values in low salt, indicating that the observed pH-dependence of stability was primarily due to these three residues. The pH-dependence of backbone amide NMR resonances demonstrated that perturbation of all three pK(a) values was primarily the result of side-chain to backbone amide hydrogen bonds. Few of the significantly perturbed acidic pK(a) values in Sac7d and Sso7d could be attributed to primarily ion pair or electrostatic interactions. A smaller perturbation of E48 (E47 in Sac7d) was ascribed to an ion pair interaction that may be important in defining the DNA binding surface. The small number (three) of significantly altered pK(a) values was in good agreement with a linkage analysis of the temperature, pH, and salt-dependence of folding. The linkage of the ionization of two or more side-chains to protein folding led to apparent cooperativity in the pH-dependence of folding, although each group titrated independently with a Hill coefficient near unity. These results demonstrate that the acid pH-dependence of protein stability in these hyperthermophile proteins is due to independent titration of acidic residues with pK(a) values perturbed primarily by hydrogen bonding of the side-chain to the backbone. This work demonstrates the need for caution in using structural data alone to argue the importance of ion pairs in stabilizing hyperthermophile proteins.

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.

Differences in amino acids composition and coupling patterns between mesophilic and thermophilic proteins

Thermophilic proteins show substantially higher intrinsic thermal stability than their mesophilic counterparts. Amino acid composition is believed to alter the intrinsic stability of proteins. Several investigations and mutagenesis experiment have been carried out to understand the amino acid composition for the thermostability of proteins. This review presents some generalized features of amino acid composition found in thermophilic proteins, including an increase in residue hydrophobicity, a decrease in uncharged polar residues, an increase in charged residues, an increase in aromatic residues, certain amino acid coupling patterns and amino acid preferences for thermophilic proteins. The differences of amino acids composition between thermophilic and mesophilic proteins are related to some properties of amino acids. These features provide guidelines for engineering mesophilic protein to thermophilic protein.

Dynamic Arrangement of Ion Pairs and Individual Contributions to the Thermal Stability of the Cofactor-Binding Domain of Glutamate Dehydrogenase from Thermotoga maritima

The dynamics of a hyperthermophilic protein fragment in a water environment, as studied by performing molecular dynamics (MD) simulations at various temperatures, is compared to the dynamical behavior of a homologous mesophilic protein simulated under identical conditions. The effects on the stability of the spatial arrangement and mobility of the charged residues in solution were quantified by calculating free energy changes upon salt bridge formation in these proteins. Electrostatic free energy terms derived from a thermodynamic cycle were obtained by solving the linearized Poisson-Boltzmann equation for a series of protein conformations generated by MD simulations and placed subsequently in a continuum solvent medium. Our results show that the ion pairs are electrostatically stabilizing in most of the cases, but their individual contributions vary significantly. The greater contribution of the charged residues to the stability of the hyperthermophilic protein as compared with the mesophilic counterpart was evidenced only by the calculations that included conformations sampled at 343 and 373 K. The "dynamic" structure of the hyperthermophilic protein fragment simulated at elevated temperatures reveals an optimum placement of the ionizable residues within the protein structure as well as the role of their cooperative interactions in promoting thermal stability. The thermodynamic properties such as electrostatic free energy differences, configurational entropies, and specific heat capacities calculated in the dynamic context of the protein structure provided new insight into the mechanism of protein thermostabilization.

Comparative study of generalized born models: Born radii and peptide folding

In this study, we have implemented four analytical generalized Born (GB) models and investigated their performance in conjunction with the GROMOS96 force field. The four models include that of Still and co-workers, the HCT model of Cramer, Truhlar, and co-workers, a modified form of the AGB model of Levy and co-workers, and the GBMV2 model of Brooks and co-workers. The models were coded independently and implemented in the GROMOS software package and in TINKER. They were compared in terms of their ability to reproduce the results of Poisson-Boltzmann (PB) calculations and in their performance in the ab initio peptide folding of two peptides, one that forms a beta-hairpin in solution and one that forms an alpha-helix. In agreement with previous work, the GBMV2 model is most successful in reproducing PB results while the other models tend to underestimate the effective Born radii of buried atoms. In contrast, stochastic dynamics simulations on the folding of the two peptides, the C-terminus beta-hairpin of the B1 domain of protein G and the alanine-based alpha-helical peptide 3K(I), suggest that the simpler GB models are more effective in sampling conformational space. Indeed, the Still model used in conjunction with the GROMOS96 force field is able to fold the hairpin peptide to a native-like structure without the benefit of enhanced sampling techniques. This is due in part to the properties of the united-atom GROMOS96 force field which appears to be more flexible, and hence to sample more efficiently, than force fields such as OPLSAA. Our results suggest a general strategy which involves using different combinations of force fields and solvent models in different applications, for example, using GROMOS96 and a simple GB model in sampling and OPLSAA and a more accurate GB model in refinement. The fact that various methods have been implemented in a unified way should facilitate the testing and subsequent use of different methods to evaluate conformational free energies in different applications. Our results also bear on some general issues involved in peptide folding and structure prediction which are addressed in the Discussion.

Effect of Protein Charge on the Generation of Aggregation-Prone Conformers

This study describes how charge modification affects aggregation of ovalbumin, thereby distinguishing the role of conformational and electrostatic stability in the process. Ovalbumin variants were engineered using chemical methylation or succinylation to obtain a range of protein net charge from -1 to -26. Charge modification significantly affected the denaturation temperature. From urea-induced equilibrium denaturation studies, it followed that unfolding proceeded via an intermediate state. However, the heat-induced denaturation process could still be described as a two-state irreversible unfolding transition, suggesting that the occurrence of an intermediate has no influence on the kinetics of unfolding. By monitoring the aggregation kinetics, the net charge was found not to be discriminative in the process. It is concluded that the dominant factor determining ovalbumin aggregation propensity is the rate of denaturation and not electrostatic repulsive forces.

Electrostatic contributions to protein stability and folding energy

The ability to predict the thermal stability of proteins based on their corresponding sequence is a problem of great fundamental and practical importance. Here we report an approach for calculating the electrostatic contribution to protein stability based on the use of the semimacroscopic protein dipole Langevin dipole (PDLD/S) in its linear response approximation version for self-energy with a dielectric constant, (epsilon(p)) and an effective dielectric for charge-charge interactions (epsilon(eff)). The method is applied to the test cases of ubiquitin, lipase, dihydrofolate reductase and cold shock proteins with series of epsilon(p) and epsilon(eff). It is found that the optimal values of these dielectric constants lead to very promising results, both for the relative stability and the absolute folding energy. Consideration of the specific values of the optimal dielectric constants leads to an exciting conceptual description of the reorganization effect during the folding process. Although this description should be examined by further microscopic studies, the practical use of the current approach seems to offer a powerful tool for protein design and for studies of the energetics of protein folding.

Understanding antibody–antigen associations by molecular dynamics simulations: Detection of important intra- and inter-molecular salt bridges

1 NSec molecular dynamics (MD) simulation of anti-hen egg white antibody, HyHEL63 (HH63), complexed with HEL reveals important molecular interactions, not revealed in its X-ray crystal structure. These molecular interactions were predicted to be critical for the complex formation, based on structure-function studies of this complex and 3-other anti-HEL antibodies, HH8, HH10 and HH26, HEL complexes. All four antibodies belong to the same structural family, referred to here as HH10 family. Ala scanning results show that they recognize 'coincident epitopes'. 1 NSec explicit, with periodic boundary condition, MD simulation of HH63- HEL reveals the presence of functionally important saltbridges. Around 200 ps in vacuo and an additional 20 ps explicit simulation agree with the observations from 1 Nsec simulation. Intra-molecular salt-bridges predicted to play significant roles in the complex formation, were revealed during MD simulation. A very stabilizing saltbridge network, and another intra-molecular salt-bridge, at the binding site of HEL, revealed during the MD simulation, is proposed to predipose binding site geometry for specific binding. All the revealed saltbridges are present in one or more of the other three complexes and/or involve \"hot-spot\" epitope and paratope residues. Most of these charged epitope residues make large contribution to the binding free energy. The "hot spot" epitope residue Lys97Y, which significantly contributes to the free energy of binding in all the complexes, forms an intermolecular salt-bridge in several MD conformers. Our earlier computations have shown that this inter-molecular salt-bridge plays a significant role in determining specificity and flexibility of binding in the HH8-HEL and HH26-HEL complexes. Using a robust criterion of salt-bridge detection, this intermolecular salt-bridge was detected in the native structures of the HH8-HEL and HH26-HEL complexes, but was not revealed in the crystal structure of HH63-HEL complex. The electrostatic strength of this revealed saltbridge was very strong. During 1 Nsec MD simulation this salt-bridge networks with another inter-molecular salt-bridge to form an inter-molecular salt-bridge triad. Participation of Lys97Y in the formation of inter-molecular triad further validates the functional importance of Lys97Y in HH63-HEL associations. These results demonstrate that many important structural details of biomolecular interactions can be better understood when studied in a dynamic environment, and that MD simulations can complement and expand information obtained from static X-ray structure. This study also highlights "hot-spot" molecular interactions in HyHEL63-HEL complex.

Molecular dynamics simulations to investigate the domain swapping mechanism of human cystatin C

Human cystatin C (HCC), one of the amyloidgenic proteins, has been proved to form a dimeric structure via a domain swapping process and then cause amyloid deposits in the brains of patients suffering from Alzheimer's disease. HCC monomer consists of a core with a five-stranded antiparallel beta-sheet (beta region) wrapped around a central helix. The connectivity of these secondary structures is: (N)-beta1-alpha-beta2-L1-beta3-AS-beta4-L2-beta5-(C). In this study, various molecular dynamics simulations were conducted to investigate the conformational changes of the monomeric HCC at different temperatures (300 and 500 K) and pH levels (2, 4, and 7) to gain insight into the domain swapping mechanism. The results show that high temperature (500 K) and low pH (pH 2) will trigger the domain swapping process of HCC. We further proposed that the domain swapping mechanism of HCC follows four steps: (1) the alpha-helix moves away from the beta region; (2) the contacts between beta2 and beta3-AS disappear; (3) the beta2-L1-beta3 hairpin unfolds following the so-called "zip-up" mechanism; and finally (4) the HCC dimer is formed. Our study shows that high temperature can accelerate the unfolding of HCC and the departure of the alpha-helix from the beta-region, especially at low pH value. This is attributed to the fact that that low pH results in the protonation of the side chains of Asp, Glu, and His residues, which further disrupts the following four salt-bridge interactions stabilizing the alpha-beta interface of the native structure: Asp15-Arg53 (beta1-beta2), Glu21/20-Lys54 (helix-beta2), Asp40-Arg70 (helix-AS), and His43-Asp81 (beta2-AS).

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.

Thermal stability of isocitrate dehydrogenase from Archaeoglobus fulgidus studied by crystal structure analysis and engineering of chimers

Isocitrate dehydrogenase from Archaeoglobus fulgidus (AfIDH) has an apparent melting temperature (T(m)) of 98.5 degrees C. To identify the structural features involved in thermal stabilization of AfIDH, the structure was solved to 2.5 A resolution. AfIDH was strikingly similar to mesophilic IDH from Escherichia coli (EcIDH) and displayed almost the same number of ion pairs and ionic networks. However, two unique inter-domain networks were present in AfIDH; one three-membered ionic network between the large and the small domain and one four-membered ionic network between the clasp and the small domain. The latter ionic network was presumably reduced in size when the clasp domain of AfIDH was swapped with that of EcIDH and the T (m) decreased by 18 degrees C. Contrarily, EcIDH was only stabilized by 4 degrees C by the clasp domain of AfIDH, a result probably due to the introduction of a unique inter-subunit aromatic cluster in AfIDH that may strengthen the dimeric interface in this enzyme. A unique aromatic cluster was identified close to the N-terminus of AfIDH that could provide additional stabilization of this region. Common and unique heat adaptive traits of AfIDH with those recently observed for hyperthermophilic IDH from Aeropyrum pernix (ApIDH) and Thermotoga maritima (TmIDH) are discussed herein.

Mechanisms for stabilisation and the maintenance of solubility in proteins from thermophiles

Background

The database of protein structures contains representatives from organisms with a range of growth temperatures. Various properties have been studied in a search for the molecular basis of protein adaptation to higher growth temperature. Charged groups have emerged as key distinguishing factors for proteins from thermophiles and mesophiles.

Results

A dataset of 291 thermophile-derived protein structures is compared with mesophile proteins. Calculations of electrostatic interactions support the importance of charges, but indicate that increases in charge contribution to folded state stabilisation do not generally correlate with the numbers of charged groups. Relative propensities of charged groups vary, such as the substitution of glutamic for aspartic acid sidechains. Calculations suggest an energetic basis, with less dehydration for longer sidechains. Most other properties studied show weak or insignificant separation of proteins from moderate thermophiles or hyperthermophiles and mesophiles, including an estimate of the difference in sidechain rotameric entropy upon protein folding. An exception is increased burial of alanine and proline residues and decreased burial of phenylalanine, methionine, tyrosine and tryptophan in hyperthermophile proteins compared to those from mesophiles.

Conclusion

Since an increase in the number of charged groups for hyperthermophile proteins is separable from charged group contribution to folded state stability, we hypothesise that charged group propensity is important in the context of protein solubility and the prevention of aggregation. Accordingly we find some separation between mesophile and hyperthermophile proteins when looking at the largest surface patch that does not contain a charged sidechain. With regard to our observation that aromatic sidechains are less buried in hyperthermophile proteins, further analysis indicates that the placement of some of these groups may facilitate the reduction of folding fluctuations in proteins of the higher growth temperature organisms.

Role of the Charge–Charge Interactions in Defining Stability and Halophilicity of the CspB Proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.

Electrostatic interactions modulate the conformation of collagen I

The pH- and electrolyte-dependent charging of collagen I fibrils was analyzed by streaming potential/streaming current experiments using the Microslit Electrokinetic Setup. Differential scanning calorimetry and circular dichroism spectroscopy were applied in similar electrolyte solutions to characterize the influence of electrostatic interactions on the conformational stability of the protein. The acid base behavior of collagen I was found to be strongly influenced by the ionic strength in KCl as well as in CaCl(2) solutions. An increase of the ionic strength with KCl from 10(-4) M to 10(-2) M shifts the isoelectric point (IEP) of the protein from pH 7.5 to 5.3. However, a similar increase of the ionic strength in CaCl(2) solutions shifts the IEP from 7.5 to above pH 9. Enhanced thermal stability with increasing ionic strength was observed by differential scanning calorimetry in both electrolyte systems. In line with this, circular dichroism spectroscopy results show an increase of the helicity with increasing ionic strength. Better screening of charged residues and the formation of salt bridges are assumed to cause the stabilization of collagen I with increasing ionic strength in both electrolyte systems. Preferential adsorption of hydroxide ions onto intrinsically uncharged sites in KCl solutions and calcium binding to negatively charged carboxylic acid moieties in CaCl(2) solutions are concluded to shift the IEP and influence the conformational stability of the protein.

Structure of a Hyperthermophilic Archaeal Homing Endonuclease, I-Tsp061I: Contribution of Cross-domain Polar Networks to Thermostability

A novel LAGLIDADG-type homing endonuclease (HEase), I-Tsp061I, from the hyperthermophilic archaeon Thermoproteus sp. IC-061 16 S rRNA gene (rDNA) intron was characterized with respect to its structure, catalytic properties and thermostability. It was found that I-Tsp061I is a HEase isoschizomer of the previously described I-PogI and exhibits the highest thermostability among the known LAGLIDADG-type HEases. Determination of the crystal structure of I-Tsp061I at 2.1 A resolution using the multiple isomorphous replacement and anomalous scattering method revealed that the overall fold is similar to that of other known LAGLIDADG-type HEases, despite little sequence similarity between I-Tsp061I and those HEases. However, I-Tsp061I contains important cross-domain polar networks, unlike its mesophilic counterparts. Notably, the polar network Tyr6-Asp104-His180-107O-HOH12-104O-Asn177 exists across the two packed alpha-helices containing both the LAGLIDADG catalytic motif and the GxxxG hydrophobic helix bundle motif. Another important structural feature is the salt-bridge network Asp29-Arg31-Glu182 across N and C-terminal domain interface, which appears to contribute to the stability of the domain/domain packing. On the basis of these structural analyses and extensive mutational studies, we conclude that such cross-domain polar networks play key roles in stabilizing the catalytic center and domain packing, and underlie the hyperthermostability of I-Tsp061I.

Thermodynamic stability of a cold-adapted protein, type III antifreeze protein, and energetic contribution of salt bridges

A thermodynamic analysis of a cold-adapted protein, type III anti-freeze protein (AFP), was carried out. The results indicate that the folding equilibrium of type III AFP is a reversible, unimolecular, two-state process with no populated intermediates. Compared to most mesophilic proteins whose folding is two-state, the psychrophilic type III AFP has a much lower thermodynamic stability at 25 degrees C, approximately 3 kcal/mol, and presents a remarkably downshifted stability-temperature curve, reaching a maximum of 5 kcal/mol around 0 degrees C. Type III AFPs contain few and non-optimally distributed surface charges relative to their mesophilic homologs, the C-terminal domains of sialic acid synthases. We used thermodynamic double mutant cycles to evaluate the energetic role of every surface salt bridge in type III AFP. Two isolated salt bridges provided no contribution to stability, while the Asp36-Arg39 salt bridge, involved in a salt bridge network with the C-terminal carboxylate, had a substantial contribution (approximately 1 kcal/mol). However, this contribution was more than counteracted by the destabilizing effect of the Asp36 carboxylate itself, whose removal led to a net 30% increase in stability at 25 degrees C. This study suggests that type III AFPs may have evolved for a minimally acceptable stability at the restricted, low temperature range (around 0 degrees C) at which AFPs must function. In addition, it indicates that salt bridge networks are used in nature also for the stability of psychrophilic proteins, and has led to a type III AFP variant of increased stability that could be used for biotechnological purposes.

Side Chain Dynamics and Alternative Hydrogen Bonding in the Mechanism of Protein Thermostabilization

To elucidate the mechanism of protein thermostabilization, the thermodynamic properties of small monomeric proteins from mesophilic and thermophilic organisms have been analyzed. Molecular dynamics simulations were employed in the study of dynamic features of charged and polar side chains of amino acid residues. The basic conclusion has been made: surface charged and polar side chains with high conformational mobility can form alternative hydrogen bonded (H-bonded) donor-acceptor pairs. The correlation between the quantitative content of alternative H-bonds per residue and the temperature of maximal thermostability of proteins has been found. The proposed mechanism of protein thermostabilization suggests continuous disruption of the primary H-bonds and formation of alternative ones, which maintain constant the enthalpy value in the native state and prevent a rapid increase of the conformational entropy with the rising temperature. The analysis of the results show that the more residues located in the N- and C-terminal regions and in the extended loops that are capable of forming alternative longer-range H-bonded pairs, the higher the protein thermostability.

Structural analysis and prediction of protein mutant stability using distance and torsion potentials: Role of secondary structure and solvent accessibility

Analyzing the factors behind protein stability is a key research topic in molecular biology, and has direct implications on protein structure prediction and protein-protein interactions. We have analyzed protein stability upon point mutations using a distance-dependant pair potential representing mainly through-space interactions, and torsion angle potential representing mainly neighboring effects as a basic statistical mechanical setup for the analysis. The synergetic effect of accessible surface area and secondary structure preferences was used as a classifier for the potentials. In addition, short-, medium-, and long-range interactions of the protein environment were also analyzed. Two datasets of point mutations were taken for the comparison of theoretically predicted stabilizing energy values with experimental DeltaDeltaG and DeltaDeltaGH(2)O from thermal and chemical denaturation experiments. These include 1538 and 1603 mutations, respectively, and contain 101 proteins that share a wide range of sequence identity. The resulting force fields were carefully evaluated with different statistical tests. Results show a maximum correlation of 0.87 with a standard error of 0.71 kcal/mol between predicted and measured DeltaDeltaG values and a prediction accuracy of 85.3% (stabilizing or destabilizing) for all mutations together. A correlation of 0.77 (more than 80% prediction accuracy with a standard error of 0.95 kcal/mol) each for the test dataset of split-sample validation and fivefold crossvalidation was obtained and a correlation of 0.70 (77.4% prediction accuracy with a standard error of 1.17 kcal/mol) was shown by the jackknife test. The same model was implemented, and the results were analyzed for mutations with DeltaDeltaGH(2)O. A correlation of 0.78 (standard error 0.96 kcal/mol) was observed with a prediction efficiency of 84.65%. This model can be used for the future prediction of protein structural stability together with various experimental techniques.

Optimization of electrostatics as a strategy for cold-adaptation: a case study of cold- and warm-active elastases

Adaptation to both high and low temperatures requires proteins with special properties. While organisms living at or close to the boiling point of water need to have proteins with increased stability, other properties are required at temperatures close to the freezing point of water. Indeed, it has been shown that enzymes adapted to cold environments are less resistant to heat with a concomitant increased activity as compared to their warm-active counter-parts. Several recent studies have pointed in the direction that electrostatic interactions play a central role in temperature adaptation, and in this study we investigate the role such interactions have in adaptation of elastase from Atlantic salmon and pig. Molecular dynamics (MD) simulations have been used to generate structural ensembles at 283 and 310 K of the psychrophilic and mesophilic elastase, and a total of eight 12 ns simulations have been carried out. Even though the two homologues have a highly similar three-dimensional structure, the location and number of charged amino acids are very different. Based on the simulated structures we find that very few salt-bridges are stable throughout the simulations, and provide little stabilization/destabilization of the proteins as judged by continuum electrostatic calculations. However, the mesophilic elastase is characterized by a greater number of salt-bridges as well as a putative salt-bridge network close to the catalytic site, indicating a higher rigidity of the components involved in the catalytic cycle. In addition, subtle differences are also found in the electrostatic potentials in the vicinity of the catalytic residues, which may explain the increased catalytic efficiency of the cold-adapted elastase.

Modeling electrostatic effects in proteins

Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biological molecules. This review considers the current state of simulations of electrostatic energies in macromolecules as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielectric 'constants' in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quantitative descriptions of electrostatic effects in proteins and thus to quantitative descriptions of the function of proteins.

The crystal structure of a hyperthermostable subfamily II isocitrate dehydrogenase from Thermotoga maritima

Isocitrate dehydrogenase (IDH) from the hyperthermophile Thermotoga maritima (TmIDH) catalyses NADP+- and metal-dependent oxidative decarboxylation of isocitrate to alpha-ketoglutarate. It belongs to the beta-decarboxylating dehydrogenase family and is the only hyperthermostable IDH identified within subfamily II. Furthermore, it is the only IDH that has been characterized as both dimeric and tetrameric in solution. We solved the crystal structure of the dimeric apo form of TmIDH at 2.2 A. The R-factor of the refined model was 18.5% (R(free) 22.4%). The conformation of the TmIDH structure was open and showed a domain rotation of 25-30 degrees compared with closed IDHs. The separate domains were found to be homologous to those of the mesophilic mammalian IDHs of subfamily II and were subjected to a comparative analysis in order to find differences that could explain the large difference in thermostability. Mutational studies revealed that stabilization of the N- and C-termini via long-range electrostatic interactions were important for the higher thermostability of TmIDH. Moreover, the number of intra- and intersubunit ion pairs was higher and the ionic networks were larger compared with the mesophilic IDHs. Other factors likely to confer higher stability in TmIDH were a less hydrophobic and more charged accessible surface, a more hydrophobic subunit interface, more hydrogen bonds per residue and a few loop deletions. The residues responsible for the binding of isocitrate and NADP+ were found to be highly conserved between TmIDH and the mammalian IDHs and it is likely that the reaction mechanism is the same.

A thermodynamic comparison of HPr proteins from extremophilic organisms

A thermodynamic stability study of five histidine-containing phosphocarrier protein (HPr) homologues derived from organisms inhabiting diverse environments is described. These HPr homologues are from Bacillus subtilis (Bs), Streptococcus thermophilus (St), Bacillus staerothermophilus (Bst), Bacillus halodurans (Bh), and Oceanobacillus iheyensis (Oi). Analyses of solvent and thermal denaturation experiments provide the cardinal thermodynamic parameters, like deltaG, deltaH, deltaS, T(m), and deltaC(p), that characterize the conformational stability for each homologue. The homologue from Bacillus staerothermophilus (BstHPr) was established as the most thermostable homologue and also the homologue with highest deltaG at all temperatures. A good correlation between habitat temperature of the organism and thermal stability of the protein is also seen. Stability curves (deltaG vs T) for every homologue are also reported; these reveal very similar deltaC(p) and temperature of maximum stability (T(S)) values for all HPr homologues. Stability curves show that the higher thermal stability of some homologues is not a result of change in curvature of the curve or a shift to higher temperature, but rather a displacement of the stability curves to higher deltaG values. Stability curves also allowed estimation of deltaG at habitat temperature of the organisms, and we find good agreement between homologues. Electrostatic contributions to stability of each homologue were investigated by measuring stability as a function of varying pH and NaCl concentration, and our results suggest that most HPr homologues share similar electrostatic contributions to stability.

Probing the structural basis for the difference in thermostability displayed by family 10 xylanases

Thermostability is an important property of industrially significant hydrolytic enzymes: understanding the structural basis for this attribute will underpin the future biotechnological exploitation of these biocatalysts. The Cellvibrio family 10 (GH10) xylanases display considerable sequence identity but exhibit significant differences in thermostability; thus, these enzymes represent excellent models to examine the structural basis for the variation in stability displayed by these glycoside hydrolases. Here, we have subjected the intracellular Cellvibrio mixtus xylanase CmXyn10B to forced protein evolution. Error-prone PCR and selection identified a double mutant, A334V/G348D, which confers an increase in thermostability. The mutant has a Tm 8 degrees C higher than the wild-type enzyme and, at 55 degrees C, the first-order rate constant for thermal inactivation of A334V/G348D is 4.1 x 10(-4) min(-1), compared to a value of 1.6 x 10(-1) min(-1) for the wild-type enzyme. The introduction of the N to C-terminal disulphide bridge into A334V/G348D, which increases the thermostability of wild-type CmXyn10B, conferred a further approximately 2 degrees C increase in the Tm of the double mutant. The crystal structure of A334V/G348D showed that the introduction of Val334 fills a cavity within the hydrophobic core of the xylanase, increasing the number of van der Waals interactions with the surrounding aromatic residues, while O(delta1) of Asp348 makes an additional hydrogen bond with the amide of Gly344 and O(delta2) interacts with the arabinofuranose side-chain of the xylose moiety at the -2 subsite. To investigate the importance of xylan decorations in productive substrate binding, the activity of wild-type CmXyn10B, the mutant A334V/G348D, and several other GH10 xylanases against xylotriose and xylotriose containing an arabinofuranose side-chain (AX3) was assessed. The enzymes were more active against AX3 than xylotriose, providing evidence that the arabinose side-chain makes a generic contribution to substrate recognition by GH10 xylanases.

Using Protein Design To Dissect the Effect of Charged Residues on Metal Binding and Protein Stability

Ca2+ controls biological processes by interacting with proteins with different affinities, which are largely influenced by the electrostatic interaction from the local negatively charged ligand residues in the coordination sphere. We have developed a general strategy for rationally designing stable Ca2+- and Ln3+-binding proteins that retain the native folding of the host protein. Domain 1 of cluster differentiation 2 (CD2) is the host for the two designed proteins in this study. We investigate the effect of local charge on Ca2+-binding affinity based on the folding properties and metal-binding affinities of the two proteins that have similarly located Ca2+-binding sites with two shared ligand positions. While mutation and Ca2+ binding do not alter the native structure of the protein, Ca2+ binding specifically induced changes around the designed Ca2+-binding site. The designed protein with a -5 charge at the binding sphere displays a 14-, 20-, and 12-fold increase in the binding affinity for Ca2+, Tb3+, and La3+, respectively, compared to the designed protein with a -3 charge, which suggests that higher local charges are preferred for both Ca2+ and Ln3+ binding. The localized charged residues significantly decrease the thermal stability of the designed protein with a -5 charge, which has a T(m) of 41 degrees C. Wild-type CD2 has a T(m) of 61 degrees C, which is similar to the designed protein with a -3 charge. This decrease is partially restored by Ca2+ binding. The effect on the protein stability is modulated by the environment and the secondary structure locations of the charged mutations. Our study demonstrates the capability and power of protein design in unveiling key determinants to Ca2+-binding affinity without the complexities of the global conformational changes, cooperativity, and multibinding process found in most natural Ca2+-binding proteins.

Contribution of Electrostatic Interactions, Compactness and Quaternary Structure to Protein Thermostability: Lessons from Structural Genomics of Thermotoga maritima

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.

Salting the charged surface: pH and salt dependence of protein G B1 stability

This study shows significant effects of protein surface charges on stability and these effects are not eliminated by salt screening. The stability for a variant of protein G B1 domain was studied in the pH-range of 1.5-11 at low, 0.15 M, and 2 M salt. The variant has three mutations, T2Q, N8D, and N37D, to guarantee an intact covalent chain at all pH values. The stability of the protein shows distinct pH dependence with the highest stability close to the isoelectric point. The stability is pH-dependent at all three NaCl concentrations, indicating that interactions involving charged residues are important at all three conditions. We find that 2 M salt stabilizes the protein at low pH (protein net charge is +6 and total number of charges is 6) but not at high pH (net charge is <or=-6 and total number of charges is >or=18). Furthermore, 0.15 M salt slightly decreases the stability of the protein over the pH range. The results show that a net charge of the protein is destabilizing and indicate that proteins contain charges for reasons other than improved stability. Salt seems to reduce the electrostatic contributions to stability under conditions with few total charges, but cannot eliminate electrostatic effects in highly charged systems.

Checking the conformational stability of cystatin C and its L68Q variant by molecular dynamics studies: why is the L68Q variant amyloidogenic?

Human L68Q cystatin C is one of the known human amyloidogenic proteins. In its native state it is a monomer with alpha/beta structure. Experimental evidence suggests that L68Q variant associates into dimeric intermediates and that the dimers subsequently self-assemble to form amyloid deposits and insoluble fibrils. Details of the pathway of L68Q mutant amyloid formation are unclear; however, different experimental approaches with resolutions at molecular level have provided some clues. Probably, the stability and flexibility of monomeric L68Q variant play essential roles in the early steps of amyloid formation; thus, it is necessary to characterize early conformational changes of L68Q cystatin C monomers. In this paper, we demonstrate the possibility that the differences between the monomeric forms of wild-type (wt) cystatin C and its L68Q variant are responsible for higher tendency of the L68Q cystatin C amyloidogenesis. We started our studies with the simulations of wt and L68Q cystatin C monomers. Nanosecond time scale molecular dynamics simulations at 308K were performed using AMBER7.0 program. The results show that the structure of the L68Q monomer was changed, relative to the wt cystatin C structure. The results support earlier speculation that the L68Q point mutation would easily lead to dimer formation.

All-atom contact potential approach to protein thermostability analysis

In this paper we use all-atom potential energy to define and analyze the inter-residue contacts in mesophilic and thermophilic proteins. Fifteen families of proteins are selected and each family has two representative proteins with greatly different preferred environmental temperatures. We find that both the number and energy of the contacts defined in this way show stronger correlations with the preferred temperatures of proteins than other factors used before. We also find that the charged-polar and charged-nonpolar residue contacts not only have larger contact numbers but also have lower single contact energies. Furthermore, the most important is that most of the thermophilic proteins have more charged-polar and charged-nonpolar residue contacts than their mesophilic counterparts. This suggests that they may play an important role in the thermostability of proteins, except usual charged-charged and nonpolar-nonpolar residue contacts. Charged residues may exert their profound influence by forming contacts not only with other charged residues but also with polar or nonpolar residues, thus further increasing the strength of contact network and then the thermostability of proteins.

Intersubunit Disulfide Interactions Play a Critical Role in Maintaining the Thermostability of Glucose-6-phosphate Dehydrogenase from the Hyperthermophilic Bacterium Aquifex aeolicus

Proteins from thermophilic microorganisms are stabilized by various mechanisms to preserve their native folded states at higher temperatures. A thermostable glucose-6-phosphate dehydrogenase (tG6PDH) from the hyperthermophilic bacterium Aquifex aeolicus was expressed as a recombinant protein in Escherichia coli. The A. aeolicus G6PDH is a homodimer exhibiting remarkable thermostability (t1/2 = 24 hr at 90 degrees C). Based on homology modeling and upon comparison of its structure with human G6PDH, it was predicted that cysteine 184 of one subunit could form a disulfide bond with cysteine 352 of the other subunit resulting in reinforced intersubunit interactions that hold the dimer together. Site-directed mutagenesis was performed on tG6PDH to convert C184 and C352 to serines. The tG6PDH double mutant exhibited a dramatic decrease in the half-life from 24 hr to 3 hr at 90 degrees C. The same decrease in half-life was also found when either C184 or C352 was mutated to serine. The result indicates that C184 and C352 may play a crucial role in strengthening the dimer interface through disulfide bond formation, thereby contributing to the thermal stability of the enzyme.