Increasing protein stability by altering long-range coulombic interactions

Revisiting structural organization of proteins at high temperature from a network perspective

Interactions between distantly placed amino acids in the primary chain (long-range) play a very crucial role in the formation and stabilization of the tertiary structure of a protein, while interactions between closely placed amino acids in the primary chain (short-range) mostly stabilize the secondary structures. Every protein needs to maintain marginal stability in order to perform its physiological functions in its native environment. The requirements for this stability in mesophilic and thermophilic proteins are different. Thermophilic proteins need to form more interactions as well as more stable interactions to survive in the extreme environment, they live in. Here, we aim to find out how the interacting amino acids in three-dimensional space are positioned in the primary chains in thermophilic and mesophilic. How does this arrangement help thermophiles to maintain their structural integrity at high temperatures? Working on a dataset of 1560 orthologous pairs we perceive that thermophiles are not only enriched with long-range interactions, they feature bigger connected clusters and higher network densities compared to their mesophilic orthologs, at higher interaction strengths between the amino acids. Moreover, we have observed the enrichment of different types of interactions at different secondary structural regions.

Charge reversal mutations in mesophilic-thermophilic orthologous protein pairs and their role in enhancing coulombic interaction energy

Proteins from thermophilic organisms are a matter of immense interest for decades because of its application in fields like de-novo protein design, thermostable variants of biocatalysts etc. Previous studies have found several sequence and structural adaptations related to thermal stability, while charge reversal study remains ignored. Here we address whether charge reversal mutations naturally occur in mesophilic-thermophilic/hyperthermophilic orthologous proteins. Do they contribute to thermal stability? Our systematic study on 1550 mesophilic-thermophilic/hyperthermophilic orthologous protein pairs with remarkable structural and topological similarity, shows gain in coulombic interaction energy in thermophilic/hyperthermophilic proteins at short range associated with partially exposed and buried charge reversal mutations, which may enhance thermostability. Our findings call forth its application in future protein engineering studies. Communicated by Ramaswamy H. Sarma.

Stability profile of the neuronal SNARE complex reflects its potency to drive fast membrane fusion

Soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs) form the SNARE complex to mediate most fusion events of the secretory pathway. The neuronal SNARE complex is featured by its high stability and half-zippered conformation required for driving robust and fast synaptic exocytosis. However, these two features seem to be thermodynamically mutually exclusive. In this study, we have employed temperature-dependent disassociation assays and single-molecule Förster resonance energy transfer (FRET) experiments to analyze the stability and conformation of the neuronal SNARE complex. We reclassified the amino acids of the SNARE motif into four sub-groups (core, core-side I and II, and non-contact). Our data showed that the core residues predominantly contribute to the complex stability to meet a basal requirement for SNARE-mediated membrane fusion, while the core-side residues exert an unbalanced effect on the N- and C-half bundle stability that determines the half-zippered conformation of the neuronal SNARE complex, which would accommodate essential regulations by complexins and synaptotagmins for fast Ca2+-triggered membrane fusion. Furthermore, our data confirmed a strong coupling of folding energy between the N- and C-half assembly of the neuronal SNARE complex, which rationalizes the strong potency of the half-zippered conformation to conduct robust and fast fusion. Overall, these results uncovered that the stability profile of the neuronal SNARE complex reflects its potency to drive fast and robust membrane fusion. Based on these results, we also developed a new parameter, the stability factor (Fs), to characterize the overall stability of the neuronal SNARE complex and resolved a linear correlation between the stability and inter-residue coulombic interactions of the neuronal SNARE complex, which would help rationally design artificial SNARE complexes and remold functional SNARE complexes with desirable stability.

pH and alcohol induced structural transition in Ntf2 a nuclear transport factor of Saccharomyces cerevisiae

Ntf2 is a nuclear envelope protein, which play a pivotal role in nucleocytoplasmic transport and mediates the nuclear import of RanGDP. It interacts with various nucleoporins along with Ran-GDP and part of a multicomponent system that assembles at the nuclear pore complex (NCP) during nuclear import. Here, we have described the biophysical characterization of Ntf2 from Saccharomyces cerevisiae. Recombinant Ntf2 showed increment in the β-sheet content as well as decrement in the α-helix content from pH-7.0 to pH-4.0. A subsequent decrease in the pH led to increment in the α-helical content along with decrement in β-sheet content. Intrinsic fluorescence studies demonstrated the unfolding of the protein below physiological pH. Ntf2 showed stabilization as well as phenomenal phase transition (β sheet to α helix) by increase in alcohol concentration from 10% to 70%. Further increase in alcohol concentration (90%) resulted in residual secondary structure in Ntf2 protein. Presence of ammonium sulfate also stabilizes the secondary structure of Ntf2 protein. The structural characterization reveals the flexibility and the stability of Ntf2 at various conditions. These structural alterations in Ntf2 protein probably occurs in the course of nucleocytoplasmic transport when it interacts with other proteins moving towards its final destination.

NADPH biosensor-based identification of an alcohol dehydrogenase variant with improved catalytic properties caused by a single charge reversal at the protein surface

Alcohol dehydrogenases (ADHs) are used in reductive biotransformations for the production of valuable chiral alcohols. In this study, we used a high-throughput screening approach based on the NADPH biosensor pSenSox and fluorescence-activated cell sorting (FACS) to search for variants of the NADPH-dependent ADH of Lactobacillus brevis (LbADH) with improved activity for the reduction of 2,5-hexanedione to (2R,5R)-hexanediol. In a library of approx. 1.4 × 10⁶ clones created by random mutagenesis we identified the variant LbADH^K71E. Kinetic analysis of the purified enzyme revealed that LbADH^K71E had a ~ 16% lowered K_M value and a 17% higher V_max for 2,5-hexanedione compared to the wild-type LbADH. Higher activities were also observed for the alternative substrates acetophenone, acetylpyridine, 2-hexanone, 4-hydroxy-2-butanone, and methyl acetoacetate. K71 is solvent-exposed on the surface of LbADH and not located within or close to the active site. Therefore, K71 is not an obvious target for rational protein engineering. The study demonstrates that high-throughput screening using the NADPH biosensor pSenSox represents a powerful method to find unexpected beneficial mutations in NADPH-dependent alcohol dehydrogenases that can be favorable in industrial biotransformations.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Networks of electrostatic and hydrophobic interactions modulate the complex folding free energy surface of a designed βα protein

The successful de novo design of proteins can provide insights into the physical chemical basis of stability, the role of evolution in constraining amino acid sequences, and the production of customizable platforms for engineering applications. Previous guanidine hydrochloride (GdnHCl; an ionic denaturant) experiments of a designed, naturally occurring βα fold, Di-III_14, revealed a cooperative, two-state unfolding transition and a modest stability. Continuous-flow mixing experiments in our laboratory revealed a simple two-state reaction in the microsecond to millisecond time range and consistent with the thermodynamic results. In striking contrast, the protein remains folded up to 9.25 M in urea, a neutral denaturant, and hydrogen exchange (HDX) NMR analysis in water revealed the presence of numerous high-energy states that interconvert on a time scale greater than seconds. The complex protection pattern for HDX corresponds closely with a pair of electrostatic networks on the surface and an extensive network of hydrophobic side chains in the interior of the protein. Mutational analysis showed that electrostatic and hydrophobic networks contribute to the resistance to urea denaturation for the WT protein; remarkably, single charge reversals on the protein surface restore the expected urea sensitivity. The roughness of the energy surface reflects the densely packed hydrophobic core; the removal of only two methyl groups eliminates the high-energy states and creates a smooth surface. The design of a very stable βα fold containing electrostatic and hydrophobic networks has created a complex energy surface rarely observed in natural proteins.

Increasing Enzyme Stability and Activity through Hydrogen Bond-Enhanced Halogen Bonds

The construction of more stable proteins is important in biomolecular engineering, particularly in the design of biologics-based therapeutics. We show here that replacing the tyrosine at position 18 (Y18) of T4 lysozyme with the unnatural amino acid m-chlorotyrosine (^mClY) increases both the thermal stability (increasing the melting temperature by ∼1 °C and the melting enthalpy by 3 kcal/mol) and the enzymatic activity at elevated temperatures (15% higher than that of the parent enzyme at 40 °C) of this classic enzyme. The chlorine of ^mClY forms a halogen bond (XB) to the carbonyl oxygen of the peptide bond at glycine 28 (G28) in a tight loop near the active site. In this case, the XB potential of the typically weak XB donor Cl is shown from quantum chemical calculations to be significantly enhanced by polarization via an intramolecular hydrogen bond (HB) from the adjacent hydroxyl substituent of the tyrosyl side chain, resulting in a distinctive synergistic HB-enhanced XB (or HeX-B for short) interaction. The larger halogens (bromine and iodine) are not well accommodated within this same loop and, consequently, do not exhibit the effects on protein stability or function associated with the HeX-B interaction. Thus, we have for the first time demonstrated that an XB can be engineered to stabilize and increase the activity of an enzyme, with the increased stabilizing potential of the HeX-B further extending the application of halogenated amino acids in the design of more stable protein therapeutics.

Probing Denatured State Conformational Bias in a Three-Helix Bundle, UBA(2), Using a Cytochrome c Fusion Protein

Previous work with the four-helix-bundle protein cytochrome c' from Rhodopseudomonas palustris using histidine-heme loop formation methods revealed fold-specific deviations from random coil behavior in its denatured state ensemble. To examine the generality of this finding, we extend this work to a three-helix-bundle polypeptide, the second ubiquitin-associated domain, UBA(2), of the human DNA excision repair protein. We use yeast iso-1-cytochrome c as a scaffold, fusing the UBA(2) domain at the N-terminus of iso-1-cytochrome c. We have engineered histidine into highly solvent accessible positions of UBA(2), creating six single histidine variants. Guanidine hydrochloride denaturation studies show that the UBA(2)-cytochrome c fusion protein unfolds in a three-state process with iso-1-cytochrome c unfolding first. Furthermore, engineered histidine residues in UBA(2) strongly destabilize the iso-1-cytochrome c domain. Equilibrium and kinetic histidine-heme loop formation measurements in the denatured state at 4 and 6 M guanidine hydrochloride show that loop stability decreases as the size of the histidine-heme loop increases, in accord with the Jacobson-Stockmayer equation. However, we observe that the His27-heme loop is both more stable than expected from the Jacobson-Stockmayer relationship and breaks more slowly than expected. These results show that the sequence near His27, which is in the reverse turn between helices 2 and 3 of UBA(2), is prone to persistent interactions in the denatured state. Therefore, consistent with our results for cytochrome c', this reverse turn sequence may help to establish the topology of this fold by biasing the conformational distribution of the denatured state.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Linking computation and experiments to study the role of charge-charge interactions in protein folding and stability

Over the past two decades there has been an increase in appreciation for the role of surface charge-charge interactions in protein folding and stability. The perception shifted from the belief that charge-charge interactions are not important for protein folding and stability to the near quantitative understanding of how these interactions shape the folding energy landscape. This led to the ability of computational approaches to rationally redesign surface charge-charge interactions to modulate thermodynamic properties of proteins. Here we summarize our progress in understanding the role of charge-charge interactions for protein stability using examples drawn from my own laboratory and touch upon unanswered questions.

Role of Proteome Physical Chemistry in Cell Behavior

We review how major cell behaviors, such as bacterial growth laws, are derived from the physical chemistry of the cell’s proteins. On one hand, cell actions depend on the individual biological functionalities of their many genes and proteins. On the other hand, the common physics among proteins can be as important as the unique biology that distinguishes them. For example, bacterial growth rates depend strongly on temperature. This dependence can be explained by the folding stabilities across a cell’s proteome. Such modeling explains how thermophilic and mesophilic organisms differ, and how oxidative damage of highly charged proteins can lead to unfolding and aggregation in aging cells. Cells have characteristic time scales. For example, E. coli can duplicate as fast as 2–3 times per hour. These time scales can be explained by protein dynamics (the rates of synthesis and degradation, folding, and diffusional transport). It rationalizes how bacterial growth is slowed down by added salt. In the same way that the behaviors of inanimate materials can be expressed in terms of the statistical distributions of atoms and molecules, some cell behaviors can be expressed in terms of distributions of protein properties, giving insights into the microscopic basis of growth laws in simple cells.

Apoglobin Stability Is the Major Factor Governing both Cell-free and in Vivo Expression of Holomyoglobin

Expression levels in animal muscle tissues and in Escherichia coli vary widely for naturally occurring mammalian myoglobins (Mb). To explore this variation, we developed an in vitro transcription and wheat germ extract-based translation assay to examine quantitatively the factors that govern expression of holoMb. We constructed a library of naturally occurring Mbs from two terrestrial and four deep-diving aquatic mammals and three distal histidine mutants designed to enhance apoglobin stability but decrease hemin affinity. A strong linear correlation is observed between cell-free expression levels of holo-metMb variants and their corresponding apoglobin stabilities, which were measured independently by guanidine HCl-induced unfolding titrations using purified proteins. In contrast, there is little dependence of expression on hemin affinity. Our results confirm quantitatively that deep diving mammals have highly stable Mbs that express to higher levels in animal myocytes, E. coli, and the wheat germ cell-free system than Mbs from terrestrial mammals. Our theoretical analyses show that the rate of aggregation of unfolded apoMb is very large, and as a result, the key factor for high level expression of holoMb, and presumably other heme proteins, is an ultra high fraction of folded, native apoglobin that is capable of rapidly binding hemin. This fraction is determined by the overall equilibrium folding constant and not hemin affinity. These results also demonstrate that the cell-free transcription/translation system can be used as a high throughput platform to screen for apoglobin stability without the need to generate large amounts of protein for in vitro unfolding measurements.

Positive selection in octopus haemocyanin indicates functional links to temperature adaptation

BACKGROUND

Octopods have successfully colonised the world's oceans from the tropics to the poles. Yet, successful persistence in these habitats has required adaptations of their advanced physiological apparatus to compensate impaired oxygen supply. Their oxygen transporter haemocyanin plays a major role in cold tolerance and accordingly has undergone functional modifications to sustain oxygen release at sub-zero temperatures. However, it remains unknown how molecular properties evolved to explain the observed functional adaptations. We thus aimed to assess whether natural selection affected molecular and structural properties of haemocyanin that explains temperature adaptation in octopods.

RESULTS

Analysis of 239 partial sequences of the haemocyanin functional units (FU) f and g of 28 octopod species of polar, temperate, subtropical and tropical origin revealed natural selection was acting primarily on charge properties of surface residues. Polar octopods contained haemocyanins with higher net surface charge due to decreased glutamic acid content and higher numbers of basic amino acids. Within the analysed partial sequences, positive selection was present at site 2545, positioned between the active copper binding centre and the FU g surface. At this site, methionine was the dominant amino acid in polar octopods and leucine was dominant in tropical octopods. Sites directly involved in oxygen binding or quaternary interactions were highly conserved within the analysed sequence.

CONCLUSIONS

This study has provided the first insight into molecular and structural mechanisms that have enabled octopods to sustain oxygen supply from polar to tropical conditions. Our findings imply modulation of oxygen binding via charge-charge interaction at the protein surface, which stabilize quaternary interactions among functional units to reduce detrimental effects of high pH on venous oxygen release. Of the observed partial haemocyanin sequence, residue 2545 formed a close link between the FU g surface and the active centre, suggesting a role as allosteric binding site. The prevalence of methionine at this site in polar octopods, implies regulation of oxygen affinity via increased sensitivity to allosteric metal binding. High sequence conservation of sites directly involved in oxygen binding indicates that functional modifications of octopod haemocyanin rather occur via more subtle mechanisms, as observed in this study.

Forces stabilizing proteins

The goal of this article is to summarize what has been learned about the major forces stabilizing proteins since the late 1980s when site-directed mutagenesis became possible. The following conclusions are derived from experimental studies of hydrophobic and hydrogen bonding variants. (1) Based on studies of 138 hydrophobic interaction variants in 11 proteins, burying a -CH2- group on folding contributes 1.1±0.5 kcal/mol to protein stability. (2) The burial of non-polar side chains contributes to protein stability in two ways: first, a term that depends on the removal of the side chains from water and, more importantly, the enhanced London dispersion forces that result from the tight packing in the protein interior. (3) Based on studies of 151 hydrogen bonding variants in 15 proteins, forming a hydrogen bond on folding contributes 1.1±0.8 kcal/mol to protein stability. (4) The contribution of hydrogen bonds to protein stability is strongly context dependent. (5) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (6) Polar group burial can make a favorable contribution to protein stability even if the polar group is not hydrogen bonded. (7) Hydrophobic interactions and hydrogen bonds both make large contributions to protein stability.

The role of surface electrostatics on the stability, function and regulation of human cystathionine β-synthase, a complex multidomain and oligomeric protein

Human cystathionine β-synthase (hCBS) is a key enzyme of sulfur amino acid metabolism, controlling the commitment of homocysteine to the transsulfuration pathway and antioxidant defense. Mutations in hCBS cause inherited homocystinuria (HCU), a rare inborn error of metabolism characterized by accumulation of toxic homocysteine in blood and urine. hCBS is a complex multidomain and oligomeric protein whose activity and stability are independently regulated by the binding of S-adenosyl-methionine (SAM) to two different types of sites at its C-terminal regulatory domain. Here we study the role of surface electrostatics on the complex regulation and stability of hCBS using biophysical and biochemical procedures. We show that the kinetic stability of the catalytic and regulatory domains is significantly affected by the modulation of surface electrostatics through noticeable structural and energetic changes along their denaturation pathways. We also show that surface electrostatics strongly affect SAM binding properties to those sites responsible for either enzyme activation or kinetic stabilization. Our results provide new insight into the regulation of hCBS activity and stability in vivo with implications for understanding HCU as a conformational disease. We also lend experimental support to the role of electrostatic interactions in the recently proposed binding modes of SAM leading to hCBS activation and kinetic stabilization.

An in silico method for designing thermostable variant of a dimeric mesophilic protein based on its 3D structure

Designing proteins with enhanced thermostability has been a major interest of protein engineering because of its potential industrial applications. Here, we have presented a computational method for designing dimeric thermostable protein based on rational mutations on a mesophilic protein. Experimental and structural data indicate that the surface stability of a protein is a major factor controlling denaturation of a protein and ion-pairs are most efficient in enhancing the stability of the surfaces of the monomers and the interface between them. Our mutation based strategy is to first identify several polar or charged residues on the protein surface, interacting weakly with the rest of the protein and then replacing the side-chains of suitable neighboring residues to increase the interaction between these two residues. In stabilizing the interface, mutation is done in the interface for forming an ion pairs between the monomers. Application of this design strategy to a homo-dimeric protein and a hetero-dimeric protein as examples has produced excellent results. In both the cases the designed mutated proteins including the individual monomers and the interfaces were found to be considerably more stable than the respective mesophilic proteins as judged by self-energies and residue-wise interaction patterns. This method is easily applicable to any multi-meric proteins.

Simplicity within the complexity: bilateral impact of DMSO on the functional and unfolding patterns of α-chymotrypsin

New understanding of the fundamental links between protein stability, conformational flexibility and function, can be gained through synergic studies on their catalytic and folding/unfolding properties under the influence of stabilizing/destabilizing additives. We explored an impact of dimethyl sulfoxide (DMSO), the moderate effector of multilateral action, on the kinetic (functional) and thermodynamic (thermal unfolding) patterns of a hydrolytic enzyme, α-chymotrypsin (α-CT), over a wide range of additive concentrations, 0-70% (v/v). Both the calorimetric and kinetic data exhibited rich behavior pointing to the complex interplay of global/local stability (and flexibility) patterns. The complex action of DMSO is explained through the negative and positive preferential solvation motifs that prevail for the extreme opposite, native-like and unfolded states, respectively, implying essential stabilization of compact domains by enhancement of interfacial water networks and destabilization of a flexible active site by direct binding of DMSO to the unoccupied specific positions intended for elongated polypeptide substrates.

Structure of a thermophilic cyanobacterial b6f-type Rieske protein

The `Rieske protein' PetC is one of the key subunits of the cytochrome b(6)f complex. Its Rieske-type [2Fe-2S] cluster participates in the photosynthetic electron-transport chain. Overexpression and careful structure analysis at 2.0 Å resolution of the extrinsic soluble domain of PetC from the thermophilic cyanobacterium Thermosynechococcus elongatus BP-1 enabled in-depth spectroscopic and structural characterization and suggested novel structural features. In particular, both the protein structure and the positions of the internal water molecules unexpectedly showed a higher similarity to eukaryotic PetCs than to other prokaryotic PetCs. The structure also revealed a deep pocket on the PetC surface which is oriented towards the membrane surface in the whole complex. Its surface properties suggest a binding site for a hydrophobic compound and the complete conservation of the pocket-forming residues in all known PetC sequences indicates the functional importance of this pocket in the cytochrome b(6)f complex.

Formyl-coenzyme A (CoA):oxalate CoA-transferase from the acidophile Acetobacter aceti has a distinctive electrostatic surface and inherent acid stability

Bacterial formyl-CoA:oxalate CoA-transferase (FCOCT) and oxalyl-CoA decarboxylase work in tandem to perform a proton-consuming decarboxylation that has been suggested to have a role in generalized acid resistance. FCOCT is the product of uctB in the acidophilic acetic acid bacterium Acetobacter aceti. As expected for an acid-resistance factor, UctB remains folded at the low pH values encountered in the A. aceti cytoplasm. A comparison of crystal structures of FCOCTs and related proteins revealed few features in UctB that would distinguish it from nonacidophilic proteins and thereby account for its acid stability properties, other than a strikingly featureless electrostatic surface. The apparently neutral surface is a result of a "speckled" charge decoration, in which charged surface residues are surrounded by compensating charges but do not form salt bridges. A quantitative comparison among orthologs identified a pattern of residue substitution in UctB that may be a consequence of selection for protein stability by constant exposure to acetic acid. We suggest that this surface charge pattern, which is a distinctive feature of A. aceti proteins, creates a stabilizing electrostatic network without stiffening the protein or compromising protein-solvent interactions.

Electrostatic contribution of surface charge residues to the stability of a thermophilic protein: benchmarking experimental and predicted pKa values

Optimization of the surface charges is a promising strategy for increasing thermostability of proteins. Electrostatic contribution of ionizable groups to the protein stability can be estimated from the differences between the pKa values in the folded and unfolded states of a protein. Using this pKa-shift approach, we experimentally measured the electrostatic contribution of all aspartate and glutamate residues to the stability of a thermophilic ribosomal protein L30e from Thermococcus celer. The pKa values in the unfolded state were found to be similar to model compound pKas. The pKa values in both the folded and unfolded states obtained at 298 and 333 K were similar, suggesting that electrostatic contribution of ionizable groups to the protein stability were insensitive to temperature changes. The experimental pKa values for the L30e protein in the folded state were used as a benchmark to test the robustness of pKa prediction by various computational methods such as H++, MCCE, MEAD, pKD, PropKa, and UHBD. Although the predicted pKa values were affected by crystal contacts that may alter the side-chain conformation of surface charged residues, most computational methods performed well, with correlation coefficients between experimental and calculated pKa values ranging from 0.49 to 0.91 (p<0.01). The changes in protein stability derived from the experimental pKa-shift approach correlate well (r = 0.81) with those obtained from stability measurements of charge-to-alanine substituted variants of the L30e protein. Our results demonstrate that the knowledge of the pKa values in the folded state provides sufficient rationale for the redesign of protein surface charges leading to improved protein stability.

Design of thermostable luciferases through arginine saturation in solvent-exposed loops

In most bioluminescence systems the oxidation of luciferin and production of light is catalyzed by luciferases. Protein engineering studies have shown that thermostable proteins from thermophilic organisms have a higher frequency of Arg, especially in exposed states. To further clarify the arginine saturation effect on thermostability of firefly luciferase, some of hydrophobic solvent-exposed residues in Lampyris turkestanicus luciferase are changed to arginine. All of these residues are located at the surface loops of L.turkestanicus luciferase. Starting with a luciferase mutant (E³⁵⁴Q/Arg³⁵⁶), single mutation (-Q³⁵R, -I¹⁸²R, -I²³²R and -L(300)R), double mutation (-Q³⁵R/I²³²R) and triple mutation (-Q³⁵R/I²³²R/I¹⁸²R) are added. Bioluminescence emission spectra indicate that substitution of Arg by these residues, do not effect on the maximum wavelength of emission spectrum. It should be noted, introduction of double mutation (-Q³⁵R/I²³²R) and triple mutation (-Q³⁵R/I²³²R/I¹⁸²R) were kept specific activity of firefly luciferase. By addition of positively charged residue, some specific mutations (-I²³²R, -Q³⁵R/I²³²R and -Q³⁵R/I²³²R/I¹⁸²R) showed that optimum temperature of activity was increased to 40°C which are 12 and 15°C higher than E³⁵⁴Q/Arg³⁵⁶ and wild-type luciferases, respectively. Also, after 40 min incubation of enzymes at 40°C, the relative remaining activity of wild type was only 5%, whereas for -I²³²R, -Q³⁵R/I²³²R and -Q³⁵R/I²³²R/I¹⁸²R was 60, 80 and 80% of original activity, respectively.

Analysis of electrostatic interactions in the denatured state ensemble of the N-terminal domain of L9 under native conditions

The pH dependence of protein stability is defined by the difference in the number of protons bound to the folded state and to the denatured state ensemble (DSE) as a function of pH. In many cases, the protonation behavior can be described as the sum of a set of independently titrating residues; in this case, the pH dependence of stability reflects differences in folded and DSE pK(a)'s. pH dependent stability studies have shown that there are energetically important interactions involving charged residues in the DSE of the N-terminal domain of L9 (NTL9), which affect significantly the stability of the protein. The DSE of wild type NTL9 cannot be directly characterized under native conditions because of its high stability. A destabilized double mutant of NTL9, V3AI4A, significantly populates the folded state and the DSE in the absence of denaturant. The two states are in slow exchange on the nuclear magnetic resonance time scale, and diffusion measurements indicate that the DSE is compact. The DSE pK(a)'s of all of the acidic residues were directly determined. The DSE pK(a) of Asp8 and Asp23 are depressed relative to model compounds values. Use of the mutant DSE pK(a)'s together with known native state pK(a)'s leads to a significantly improved agreement between the measured pH dependent stability and that predicted by the Tanford-Wyman linkage relationship. An analysis of the literature suggests that DSE interactions involving charged residues are relatively common and should be considered in discussions of protein stability.

Simulating electrostatic energies in proteins: perspectives and some recent studies of pKas, redox, and other crucial functional properties

Electrostatic energies provide what is arguably the most effective tool for structure-function correlation of biological molecules. Here, we provide an overview of the current state-of-the-art simulations of electrostatic energies in macromolecules, emphasizing the microscopic perspective but also relating it to macroscopic approaches. We comment on the convergence issue and other problems of the microscopic models and the ways of keeping the microscopic physics while moving to semi-macroscopic directions. We discuss the nature of the protein dielectric "constants" reiterating our long-standing point that the dielectric "constants" in semi-macroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering different functional properties including pK(a)'s, redox potentials, ion and proton channels, enzyme catalysis, ligand binding, and protein stability. We emphasize the microscopic overcharging approach for studying pK(a) 's of internal groups in proteins and give a demonstration of power of this approach. We also emphasize recent advances in coarse grained models with a physically based electrostatic treatment and provide some examples including further directions in treating voltage activated ion channels.

Contribution of hydrophobic interactions to protein stability

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

A rigidifying salt-bridge favors the activity of thermophilic enzyme at high temperatures at the expense of low-temperature activity

BACKGROUND

Thermophilic enzymes are often less active than their mesophilic homologues at low temperatures. One hypothesis to explain this observation is that the extra stabilizing interactions increase the rigidity of thermophilic enzymes and hence reduce their activity. Here we employed a thermophilic acylphosphatase from Pyrococcus horikoshii and its homologous mesophilic acylphosphatase from human as a model to study how local rigidity of an active-site residue affects the enzymatic activity.

METHODS AND FINDINGS

Acylphosphatases have a unique structural feature that its conserved active-site arginine residue forms a salt-bridge with the C-terminal carboxyl group only in thermophilic acylphosphatases, but not in mesophilic acylphosphatases. We perturbed the local rigidity of this active-site residue by removing the salt-bridge in the thermophilic acylphosphatase and by introducing the salt-bridge in the mesophilic homologue. The mutagenesis design was confirmed by x-ray crystallography. Removing the salt-bridge in the thermophilic enzyme lowered the activation energy that decreased the activation enthalpy and entropy. Conversely, the introduction of the salt-bridge to the mesophilic homologue increased the activation energy and resulted in increases in both activation enthalpy and entropy. Revealed by molecular dynamics simulations, the unrestrained arginine residue can populate more rotamer conformations, and the loss of this conformational freedom upon the formation of transition state justified the observed reduction in activation entropy.

CONCLUSIONS

Our results support the conclusion that restricting the active-site flexibility entropically favors the enzymatic activity at high temperatures. However, the accompanying enthalpy-entropy compensation leads to a stronger temperature-dependency of the enzymatic activity, which explains the less active nature of the thermophilic enzymes at low temperatures.

In silico characterization of thermostable lipases

Thermostable lipases are of high priority for industrial applications as they are endowed with the capability of carrying out diversified reactions at elevated temperatures. Extremophiles are their potential source. Sequence and structure annotation of thermostable lipases can elucidate evolution of lipases from their mesophilic counterparts with enhanced thermostability hence better industrial potential. Sequence analysis highlighted the conserved residues in bacterial and fungal thermostable lipases. Higher frequency of AXXXA motif and poly Ala residues in lid domain of thermostable Bacillus lipases were distinguishing characteristics. Comparison of amino acid composition among thermostable and mesostable lipases brought into light the role of neutral, charged and aromatic amino acid residues in enhancement of thermostability. Structural annotation of thermostable lipases with that of mesostable lipases revealed some striking features which are increment of gamma turns in thermostable lipases; being first time reported in our paper, longer beta strands, lesser beta-branched residues in helices, increase in charged-neutral hydrogen bonding pair, hydrophobic-hydrophobic contact and differences in the N-cap and C-cap residues of the α helices. Conclusively, it can be stated that subtle changes in the arrangement of amino acid residues in the tertiary structure of lipases contributes to enhanced thermostability.

Propensities of aromatic amino acids versus leucine and proline to induce residual structure in the denatured-state ensemble of iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X=Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.

Increasing protein stability: importance of DeltaC(p) and the denatured state

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site-directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the beta-turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28 degrees C higher than the wild-type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, DeltaC(p), by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.

Amino acid transport in thermophiles: characterization of an arginine-binding protein in Thermotoga maritima. 2. Molecular organization and structural stability

ABC transport systems provide selective passage of metabolites across cell membranes and typically require the presence of a soluble binding protein with high specificity to a specific ligand. In addition to their primary role in nutrient gathering, the binding proteins associated with bacterial transport systems have been studied for their potential to serve as design scaffolds for the development of fluorescent protein biosensors. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the physicochemical properties of a hyperthermophilic binding protein from Thermotoga maritima. We demonstrated preferential binding for the polar amino acid arginine and experimentally monitored the significant stabilization achieved upon binding of ligand to protein. The effect of temperature, pH, and detergent was also studied to provide a more complete picture of the protein dynamics. A protein structure model was obtained and molecular dynamic experiments were performed to investigate and couple the spectroscopic observations with specific secondary structural elements. The data determined the presence of a buried beta-sheet providing significant stability to the protein under all conditions investigated. The specific amino acid residues responsible for arginine binding were also identified. Our data on dynamics and stability will contribute to our understanding of bacterial binding protein family members and their potential biotechnological applications.

Increasing protein stability by improving beta-turns

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-turns in nine proteins with 66-370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 +/- 0.3; Range = 0.3-1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.