The optimization of protein-solvent interactions: thermostability and the role of hydrophobic and electrostatic interactions

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Modulation of Flexible Loops in Catalytic Cavities Reveals the Thermal Activation Mechanism of a Glycogen-Debranching Enzyme

Some thermophilic enzymes not only exhibit high thermostability at high temperatures but also have an activation effect by thermal incubation. However, the correlations between temperature-induced structural modulation and thermal activation are still unclear. In this study, we selected a thermophilic glycogen-debranching enzyme from Saccharolobus solfataricus STB09 (SsGDE), which was a promising starch-debranching enzyme with a thermal activation property at temperatures ranging from 50 to 70 °C, to explore the thermal activation mechanism. Molecular dynamics simulations were performed for SsGDE at 30, 50, or 70 °C to reveal the temperature dependence of structure modulation and catalytic function. The results revealed that four loops (loop1 313-337, loop2 399-418, loop3 481-513, and loop4 540-574) in SsGDE were reshaped, which made the catalytic cavity more open. The internal residues, including the catalytic triad Asp3631, Glu399, and Asp471, could be exposed, due to the structural modulation, to exert catalytic functions. We proposed that the thermal activation effect of SsGDE was closely associated with the temperature-induced modulation of the catalytic cavity, which paved the way for further engineering enzymes to achieve higher catalytic performance and stability.

Cryo-electron structures of the extreme thermostable enzymes Sulfur Oxygenase Reductase and Lumazine Synthase

Thermostable enzymes have the potential for use in a wide variety of biotechnological applications. Cryo-electron microscopy (cryo-EM) enables the imaging of biomolecules in their native aqueous environment. Here, we present high resolution cryo-EM structures of two thermostable enzymes that exhibit multimeric cage-like structures arranged into two different point-group symmetries. First, we determined the structure of the Sulfur Oxygenase Reductase (SOR) enzyme that catalyzes both the oxygenation and disproportionation of elemental sulfur in Archea and is composed of 24 homomeric units each of MW ≃ 35 kDa arranged in octahedral symmetry. The structure of SOR from Acidianus ambivalens (7X9W) was determined at 2.78 Å resolution. The active site of each subunit inside the central nanocompartment is composed of Fe3+ coordinated to two water molecules and the three amino acids (H86, H90 and E114). Second, we determined the structure of Lumazine Synthase (LS) from Aquifex aeolicus (7X7M) at 2.33 Å resolution. LS forms a cage-like structure consisting of 60 identical subunits each of MW ≃ 15 kDa arranged in a strict icosahedral symmetry. The LS subunits are interconnected by ion-pair network. Due to their thermostability and relatively easy purification scheme, both SOR and LS can serve as a model for the catalytic and structural characterization of biocatalysts as well as a benchmark for cryo-EM sample preparation, optimization of the acquisition parameters and 3D reconstruction.

Some Clues about Enzymes from Psychrophilic Microorganisms

Enzymes purified from psychrophilic microorganisms prove to be efficient catalysts at low temperatures and possess a great potential for biotechnological applications. The low-temperature catalytic activity has to come from specific structural fluctuations involving the active site region, however, the relationship between protein conformational stability and enzymatic activity is subtle. We provide a survey of the thermodynamic stability of globular proteins and their rationalization grounded in a theoretical approach devised by one of us. Furthermore, we provide a link between marginal conformational stability and protein flexibility grounded in the harmonic approximation of the vibrational degrees of freedom, emphasizing the occurrence of long-wavelength and excited vibrations in all globular proteins. Finally, we offer a close view of three enzymes: chloride-dependent α-amylase, citrate synthase, and β-galactosidase.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Why small proteins tend to have high denaturation temperatures

Data indicate that small globular proteins (consisting of less than about 70 residues) tend to have high denaturation temperatures. This finding is analysed by comparing experimental denaturation enthalpy and entropy changes of a selected set of small proteins with values calculated on the basis of average and common properties of globular proteins. The conclusion is that the denaturation entropy change is smaller than expected, leading to an increase in denaturation temperature. The proposed molecular rationalization considers the existence of long-wavelength, low-frequency vibrational modes in the native state of small proteins due to their large surface-to-interior ratio. The effect of decreasing the conformational entropy gain associated with denaturation on thermal stability is directly verified by means of an already devised theoretical model [G. Graziano, Phys. Chem. Chem. Phys. 2010, 12, 14245-14252; 2014, 16, 21755-21767].

Protein stability governed by its structural plasticity is inferred by physicochemical factors and salt bridges

Several organisms, specifically microorganisms survive in a wide range of harsh environments including extreme temperature, pH, and salt concentration. We analyzed systematically a large number of protein sequences with their structures to understand their stability and to discriminate extremophilic proteins from their non-extremophilic orthologs. Our results highlighted that the strategy for the packing of the protein core was influenced by the environmental stresses through substitutive structural events through better ionic interaction. Statistical analysis showed that a significant difference in number and composition of amino acid exist among them. The negative correlation of pairwise sequence alignments and structural alignments indicated that most of the extremophile and non-extremophile proteins didn’t contain any association for maintaining their functional stability. A significant numbers of salt bridges were noticed on the surface of the extremostable proteins. The Ramachandran plot data represented more occurrences of amino acids being present in helix and sheet regions of extremostable proteins. We also found that a significant number of small nonpolar amino acids and moderate number of charged amino acids like Arginine and Aspartic acid represented more nonplanar Omega angles in their peptide bond. Thus, extreme conditions may predispose amino acid composition including geometric variability for molecular adaptation of extremostable proteins against atmospheric variations and associated changes under natural selection pressure. The variation of amino acid composition and structural diversifications in proteins play a major role in evolutionary adaptation to mitigate climate change.

Enhanced catalytic efficiency of CotA-laccase by DNA shuffling

Bacterial CotA-laccases exhibit higher activity in alkaline pH and salt concentration conditions compared to laccases from white-rot fungi. They are considered as green catalysts in decolorizing of industrial dyes. However, CotA-laccases are limited due to the low yield and catalytic efficiency as the spore-bound nature of CotA. A DNA shuffling strategy was applied to generate a random mutation library. To improve laccase activities, a mutant (T232P/Q367R 5E29) with two amino acid substitutions was identified. The catalytic efficiency of mutant 5E29 was 1.21 fold higher compared with that of the wild-type. The Km and k_cat values of 5E29 for SGZ were of 20.3 ± 1.3 µM and 7.6 ± 2.7 s^-1. The thermal stability was a slight enhancement. Indigo Carmine and Congo red were efficiently decolorized by using this mutant at pH 9.0. These results provide that 5E29 CotA-laccase is a good candidate for biotechnology applications under alkaline condition, with an effective decolorization capability.

The hydrogen-bond network around Glu160 contributes to the structural stability of chitosanase CsnA from Renibacterium sp. QD1

CsnA, a chitosanase from Renibacterium sp. QD1, has great potential for industrial applications due to its high yield and broad pH stability. In this study, a specific Glu160 in CsnA was identified by sequence alignment, and structural analysis and MD simulation predicted that Glu160 formed a hydrogen-bond network with Lys163 and Thr114. To evaluate the effect of the network, we constructed four mutants, including E160A, E160Q, K163A, and T114A, which partially or completely destroy this network. Characterization of these mutants demonstrated that the disruption of the network significantly decreased the enzyme thermostability. The underlying mechanisms responsible for the change of thermostability analyzed by circular dichroism spectroscopy revealed that the hydrogen-bond network conferred the structural stability of CsnA. Moreover, the length of the side chain of residue at 160 impacted conformational stability of the enzyme. Taken together, the hydrogen-bond network around Glu160 plays important roles in stabilization of CsnA.

Effect of heavy water on the conformational stability of globular proteins

It is well established from the experimental point of view that the native state of globular proteins is more stable in heavy water than in water. No robust explanation, however, has been provided up to now. The application of the theoretical approach, originally devised to rationalize the general occurrence of cold denaturation, indicates that the magnitude of the solvent-excluded volume effect is slightly smaller in heavy water than in water and cannot explain the observed protein stabilization. The latter has to be due to the strength of protein-water van der Waals attractions which are weaker in heavy water due to the smaller molecular polarizability of D₂ O compared with that of H₂ O molecules. Since protein-water van der Waals attractions occur more in the denatured than in the native state, this contribution leads to a stabilization of the latter through a destabilization of the former.

Density discriminates between thermophilic and mesophilic proteins

Despite an intense interest and a remarkable number of studies on the subject, the relationships between thermostability and (primary, secondary and tertiary) structure of proteins are still not fully understood. Here, comparing the protein density - defined by the ratio between the residue number and protein excluded volume - for a set of thermophilic/mesophilic pairs, we provide evidence that this property is connected to the optimal growth temperature. In particular, our results indicate that thermophilic proteins have - in general - a lower density with respect to the mesophilic counterparts, being such a correlation more pronounced for optimal growth temperature differences greater than 40°C. The effect of the protein thermostability changes on the molecular shape is also presented.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Rigidity versus flexibility: the dilemma of understanding protein thermal stability

The role of fluctuations in protein thermostability has recently received considerable attention. In the current literature a dualistic picture can be found: thermostability seems to be associated with enhanced rigidity of the protein scaffold in parallel with the reduction of flexible parts of the structure. In contradiction to such arguments it has been shown by experimental studies and computer simulation that thermal tolerance of a protein is not necessarily correlated with the suppression of internal fluctuations and mobility. Both concepts, rigidity and flexibility, are derived from mechanical engineering and represent temporally insensitive features describing static properties, neglecting that relative motion at certain time scales is possible in structurally stable regions of a protein. This suggests that a strict separation of rigid and flexible parts of a protein molecule does not describe the reality correctly. In this work the concepts of mobility/flexibility versus rigidity will be critically reconsidered by taking into account molecular dynamics calculations of heat capacity and conformational entropy, salt bridge networks, electrostatic interactions in folded and unfolded states, and the emerging picture of protein thermostability in view of recently developed network theories. Last, but not least, the influence of high temperature on the active site and activity of enzymes will be considered.

On the mechanism of cold denaturation

A theoretical rationalization of the occurrence of cold denaturation for globular proteins was devised, assuming that the effective size of water molecules depends upon temperature [G. Graziano, Phys. Chem. Chem. Phys., 2010, 12, 14245-14252]. In the present work, it is shown that the latter assumption is not necessary. By performing the same type of calculations in water, 40% (by weight) methanol, methanol, and carbon tetrachloride, it emerges that cold denaturation occurs only in water due to the special temperature dependence of its density and the small size of its molecules. These two coupled factors determine the magnitude and the temperature dependence of the stabilizing term that measures the gain in configurational/translational entropy of water molecules upon folding of the protein. This term has to be contrasted with the destabilizing contribution measuring the loss in conformational entropy of the polypeptide chain upon folding.

Protein domain definition should allow for conditional disorder

Proteins are often classified in a binary fashion as either structured or disordered. However this approach has several deficits. Firstly, protein folding is always conditional on the physiochemical environment. A protein which is structured in some circumstances will be disordered in others. Secondly, it hides a fundamental asymmetry in behavior. While all structured proteins can be unfolded through a change in environment, not all disordered proteins have the capacity for folding. Failure to accommodate these complexities confuses the definition of both protein structural domains and intrinsically disordered regions. We illustrate these points with an experimental study of a family of small binding domains, drawn from the RNA polymerase of mumps virus and its closest relatives. Assessed at face value the domains fall on a structural continuum, with folded, partially folded, and near unstructured members. Yet the disorder present in the family is conditional, and these closely related polypeptides can access the same folded state under appropriate conditions. Any heuristic definition of the protein domain emphasizing conformational stability divides this domain family in two, in a way that makes no biological sense. Structural domains would be better defined by their ability to adopt a specific tertiary structure: a structure that may or may not be realized, dependent on the circumstances. This explicitly allows for the conditional nature of protein folding, and more clearly demarcates structural domains from intrinsically disordered regions that may function without folding.

Structural differences between thermophilic and mesophilic membrane proteins

The evolutionary adaptations of thermophilic water-soluble proteins required for maintaining stability at high temperature have been extensively investigated. Little is known about the adaptations in membrane proteins, however. Here, we compare many properties of mesophilic and thermophilic membrane protein structures, including side-chain burial, packing, hydrogen bonding, transmembrane kinks, loop lengths, hydrophobicity, and other sequence features. Most of these properties are quite similar between mesophiles and thermophiles although we observe a slight increase in side-chain burial and possibly a slight decrease in the frequency of transmembrane kinks in thermophilic membrane protein structures. The most striking difference is the increased hydrophobicity of thermophilic transmembrane helices, possibly reflecting more stringent hydrophobicity requirements for membrane partitioning at high temperature. In agreement with prior work examining transmembrane sequences, we find that thermophiles have an increase in small residues (Gly, Ala, Ser, and Val) and a strong suppression of Cys. We also find a relative dearth of most strongly polar residues (Asp, Asn, Glu, Gln, and Arg). These results suggest that in thermophiles, there is significant evolutionary pressure to offload destabilizing polar amino acids, to decrease the entropy cost of side chain burial, and to eliminate thermally sensitive amino acids.

Purification and characterization of a thermostable glutamate dehydrogenase from a thermophilic bacterium isolated from a sterilization drying oven

Glutamate dehydrogenase from axenic bacterial cultures of a new microorganism, called GWE1, isolated from the interior of a sterilization drying oven, was purified by anion-exchange and molecular-exclusion liquid chromatography. The apparent molecular mass of the native enzyme was 250.5 kDa and was shown to be an hexamer with similar subunits of molecular mass 40.5 kDa. For glutamate oxidation, the enzyme showed an optimal pH and temperature of 8.0 and 70 degrees C, respectively. In contrast to other glutamate dehydrogenases isolated from bacteria, the enzyme isolated in this study can use both NAD(+) and NADP(+) as electron acceptors, displaying more affinity for NADP(+) than for NAD(+). No activity was detected with NADH or NADPH, 2-oxoglutarate and ammonia. The enzyme was exceptionally thermostable, maintaining more than 70% of activity after incubating at 100(o)C for more than five hours suggesting being one of the most thermoestable enzymes reported in the family of dehydrogenases.

Relationship between stability and flexibility in the most flexible region of Photinus pyralis luciferase

Firefly luciferase is a protein with a large N-terminal and a small C-terminal domain. B-factor analysis shows that its C-terminal is much more flexible than its N-terminal. Studies on hyperthermophile proteins have been shown that the increased thermal stability of hyperthermophile proteins is due to their enhanced conformational rigidity and the relationship between flexibility, stability and function in most of proteins is on debate. Two mutations (D474K and D476N) in the most flexible region of firefly luciferase were designed. Thermostability analysis shows that D476N mutation doesn't have any significant effect but D474K mutation destabilized protein. On the other hand, flexibility analysis using dynamic quenching and limited proteolysis demonstrates that D474K mutation became much more flexible than wild type although D476N doesn't have any significant difference. Intrinsic and ANS fluorescence studies demonstrate that D476N mutation is brought about by structural changes without significant effect on thermostability and flexibility. Molecular modeling reveals that disruption of a salt bridge between D(474) and K(445) accompanying with some H-bond deletion may be involved in destabilization of D474K mutant.

On the molecular origin of cold denaturation of globular proteins

A polypeptide chain can adopt very different conformations, a fundamental distinguishing feature of which is the water accessible surface area, WASA, that is a measure of the layer around the polypeptide chain where the center of water molecules cannot physically enter, generating a solvent-excluded volume effect. The large WASA decrease associated with the folding of a globular protein leads to a large decrease in the solvent-excluded volume, and so to a large increase in the configurational/translational freedom of water molecules. The latter is a quantity that depends upon temperature. Simple calculations over the -30 to 150 °C temperature range, where liquid water can exist at 1 atm, show that such a gain decreases significantly on lowering the temperature below 0 °C, paralleling the decrease in liquid water density. There will be a temperature where the destabilizing contribution of the polypeptide chain conformational entropy exactly matches the stabilizing contribution of the water configurational/translational entropy, leading to cold denaturation.

Discrimination of thermophilic and mesophilic proteins

Background

There is a considerable literature on the source of the thermostability of proteins from thermophilic organisms. Understanding the mechanisms for this thermostability would provide insights into proteins generally and permit the design of synthetic hyperstable biocatalysts.

Results

We have systematically tested a large number of sequence and structure derived quantities for their ability to discriminate thermostable proteins from their non-thermostable orthologs using sets of mesophile-thermophile ortholog pairs. Most of the quantities tested correspond to properties previously reported to be associated with thermostability. Many of the structure related properties were derived from the Delaunay tessellation of protein structures.

Conclusions

Carefully selected sequence based indices discriminate better than purely structure based indices. Combined sequence and structure based indices improve performance somewhat further. Based on our analysis, the strongest contributors to thermostability are an increase in ion pairs on the protein surface and a more strongly hydrophobic interior.

Structural characterization of the organic solvent-stable cholesterol oxidase from Chromobacterium sp. DS-1

Cholesterol oxidase is of significant commercial interest as it is widely used as a biosensor for the detection of cholesterol in clinical samples, blood serum and food. Increased stability of this enzyme with regards to temperature and different solvent conditions are of great importance to the reliability and versatility of its applications. We here report the crystal structure of the cholesterol oxidase of Chromobacterium sp. DS-1 (CHOLOX). In contrast to other previously characterized cholesterol oxidases, this enzyme retains high activity in organic solvents and detergents at temperatures above 85 degrees C despite its mesophilic origin. With the availability of one other homologous oxidase of known three-dimensional structure, a detailed comparison of its sequence and structure was performed to elucidate the mechanisms of stabilization. In contrast to factors that typically contribute to the stability of thermophilic proteins, the structure of CHOLOX exhibits a larger overall cavity volume, less charged residues and less salt bridge interactions. Moreover, the vast majority of residue substitutions were found on or near the protein's solvent exposed surface. We propose that the engineering of enhanced stability may also be accomplished through selective engineering of the protein periphery rather than by redesigning its entire core.

Structure, stability, and folding of ribonuclease H1 from the moderately thermophilic Chlorobium tepidum: comparison with thermophilic and mesophilic homologues

Proteins from thermophilic organisms are able to function under conditions that render a typical mesophilic protein inactive. Pairwise comparisons of homologous mesophilic and thermophilic proteins can help to identify the energetic features of a protein's energy landscape that lead to such thermostability. Previous studies of bacterial ribonucleases H (RNases H) from the thermophile Thermus thermophilus and the mesophile Escherichia coli revealed that the thermostability arises in part from an unusually low change in heat capacity upon unfolding (DeltaC(p)) for the thermophilic protein [Hollien, J., and Marqusee, S. (1999) Biochemistry 38, 3831-3836]. Here, we have further examined how nearly identical proteins can adapt to different thermal constraints by adding a moderately thermophilic homologue to the previously characterized mesophilic and thermophilic pair. We identified a putative RNase H from Chlorobium. tepidum and demonstrated that it is an active RNase H and adopts the RNase H fold. The moderately thermophilic protein has a melting temperature (T(m)) similar to that of the mesophilic homologue yet also has a surprisingly low DeltaC(p), like the thermophilic homologue. This new RNase H folds through a pathway similar to that of the previously studied RNases H. These results suggest that lowering the DeltaC(p) may be a general strategy for achieving thermophilicity for some protein families and implicate the folding core as the major contributor to this effect. It should now be possible to design RNases H that display the desired thermophilic or mesophilic properties, as defined by their DeltaC(p) values, and therefore fine-tune the energy landscape in a predictable fashion.

Glutamic acid residues in the C-terminal extension of small heat shock protein 25 are critical for structural and functional integrity

Small heat shock proteins (sHsps) are intracellular molecular chaperones that prevent the aggregation and precipitation of partially folded and destabilized proteins. sHsps comprise an evolutionarily conserved region of 80-100 amino acids, denoted the alpha-crystallin domain, which is flanked by regions of variable sequence and length: the N-terminal domain and the C-terminal extension. Although the two domains are known to be involved in the organization of the quaternary structure of sHsps and interaction with their target proteins, the role of the C-terminal extension is enigmatic. Despite the lack of sequence similarity, the C-terminal extension of mammalian sHsps is typically a short, polar segment which is unstructured and highly flexible and protrudes from the oligomeric structure. Both the polarity and flexibility of the C-terminal extension are important for the maintenance of sHsp solubility and for complexation with its target protein. In this study, mutants of murine Hsp25 were prepared in which the glutamic acid residues in the C-terminal extension at positions 190, 199 and 204 were each replaced with alanine. The mutants were found to be structurally altered and functionally impaired. Although there were no significant differences in the environment of tryptophan residues in the N-terminal domain or in the overall secondary structure, an increase in exposed hydrophobicity was observed for the mutants compared with wild-type Hsp25. The average molecular masses of the E199A and E204A mutants were comparable with that of the wild-type protein, whereas the E190A mutant was marginally smaller. All mutants displayed markedly reduced thermostability and chaperone activity compared with the wild-type. It is concluded that each of the glutamic acid residues in the C-terminal extension is important for Hsp25 to act as an effective molecular chaperone.

Defining the role of salt bridges in protein stability

Although the energetic balance of forces stabilizing proteins has been established qualitatively over the last decades, quantification of the energetic contribution of particular interactions still poses serious problems. The reasons are the strong cooperativity and the interdependence ofnoncovalent interactions. Salt bridges are a typical example. One expects that ionizable side chains frequently form ion pairs in innumerable crystal structures. Since electrostatic attraction between opposite charges is strong per se, salt bridges can intuitively be regarded as an important factor stabilizing the native structure. Is that really so? In this chapter we critically reassess the available methods to delineate the role ofelectrostatic interactions and salt bridges to protein stability, and discuss the progress and the obstacles in this endeavor. The basic problem is that formation of salt bridges depends on the ionization properties of the participating groups, which is significantly influenced by the protein environment. Furthermore, salt bridges experience thermal fluctuations, continuously break and re-form, and their lifespan in solution is governed by the flexibility of the protein. Finally, electrostatic interactions are long-range and might be significant in the unfolded state, thus seriously influencing the energetic profile. Elimination of salt bridges by protonation/deprotonation at extreme pH or by mutation provides only rough energetic estimates, since there is no way to account for the nonadditive response of the protein moiety. From what we know so far, the strength of electrostatic interactions is strongly context-dependent, yet it is unlikely that salt bridges are dominant factors governing protein stability. Nevertheless, proteins from thermophiles and hyperthermophiles exhibit more, and frequently networked, salt bridges than proteins from the mesophilic counterparts. Increasing the thermal (not the thermodynamic) stability of proteins by optimization of charge-charge interactions is a good example for an evolutionary solution utilizing physical factors.

UPF201 archaeal specific family members reveal structural similarity to RNA-binding proteins but low likelihood for RNA-binding function

We have determined X-ray crystal structures of four members of an archaeal specific family of proteins of unknown function (UPF0201; Pfam classification: DUF54) to advance our understanding of the genetic repertoire of archaea. Despite low pairwise amino acid sequence identities (10-40%) and the absence of conserved sequence motifs, the three-dimensional structures of these proteins are remarkably similar to one another. Their common polypeptide chain fold, encompassing a five-stranded antiparallel beta-sheet and five alpha-helices, proved to be quite unexpectedly similar to that of the RRM-type RNA-binding domain of the ribosomal L5 protein, which is responsible for binding the 5S- rRNA. Structure-based sequence alignments enabled construction of a phylogenetic tree relating UPF0201 family members to L5 ribosomal proteins and other structurally similar RNA binding proteins, thereby expanding our understanding of the evolutionary purview of the RRM superfamily. Analyses of the surfaces of these newly determined UPF0201 structures suggest that they probably do not function as RNA binding proteins, and that this domain specific family of proteins has acquired a novel function in archaebacteria, which awaits experimental elucidation.

Salt bridges in the hyperthermophilic protein Ssh10b are resilient to temperature increases

A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 +/- 0.4 kJ/mol at 298 K and 2.2 +/- 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 +/- 0.2 kJ/mol at 298 K and 5.9 +/- 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability.

Teaching noncovalent interactions using protein molecular evolution

Noncovalent interactions and physicochemical properties of amino acids are important topics in biochemistry courses. Here, we present a computational laboratory where the capacity of each of the 20 amino acids to maintain different noncovalent interactions are used to investigate the stabilizing forces in a set of proteins coming from organisms adapted to different environments. Using protein sequence and structure information it is possible to evaluate the noncovalent contributions to the stabilization of a given protein fold. As a case study, we use the protein lumazine synthase from three different organisms adapted to live in extreme temperatures: one psychrophilic (optimal growth temperature, 0-20 °C), one mesophilic (optimal growth temperature, 20-50 °C), and one thermophilic (optimal growth temperature, 80-110 °C). We found that this computational laboratory for biochemistry and molecular biology courses enhances student amino acid noncovalent interaction understanding and how these interactions are involved in protein stability.

Characterization and crystallography of recombinant 7S globulins of Adzuki bean and structure-function relationships with 7S globulins of various crops

The recombinant proteins Adzuki 7S1, Adzuki 7S2, and Adzuki 7S3 were prepared through the Escherichia coli expression systems of three kinds of adzuki bean cDNAs. The recombinant proteins exhibited intrinsic thermal stabilities, surface hydrophobicities, and solubilities, although the homology of their amino acid sequences ranged from 95-98%. To understand why these individual proteins exhibited different properties, their three-dimensional structures were elucidated. The three proteins were successfully crystallized, and the three-dimensional structures of Adzuki 7S1 and Adzuki 7S3 were determined. The properties and structures of these two proteins were comprehensively compared with those of recombinant 7S globulins (soybean beta-conglycinins beta and alpha'c and mungbean 8Salpha) reported previously. It was likely that cavity sizes, hydrogen bonds, salt bridges, hydrophobic interactions, and lengths of loops determine the thermal stabilities of 7S globulins, and results indicated that cavity sizes strongly contribute to such stability. Surface hydrophobicity was also found to be determined not only by distributions of hydrophobic residues on the molecular surface. Furthermore, solubility at neutral and weak alkaline pH values at mu = 0.08 was found to be dominantly influenced by the electrostatic surface potentials.

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.