Electrostatic contributions to the stability of hyperthermophilic proteins

Protein thermal stability enhancement by designing salt bridges: a combined computational and experimental study

Protein thermal stability is an important factor considered in medical and industrial applications. Many structural characteristics related to protein thermal stability have been elucidated, and increasing salt bridges is considered as one of the most efficient strategies to increase protein thermal stability. However, the accurate simulation of salt bridges remains difficult. In this study, a novel method for salt-bridge design was proposed based on the statistical analysis of 10,556 surface salt bridges on 6,493 X-ray protein structures. These salt bridges were first categorized based on pairing residues, secondary structure locations, and Cα–Cα distances. Pairing preferences generalized from statistical analysis were used to construct a salt-bridge pair index and utilized in a weighted electrostatic attraction model to find the effective pairings for designing salt bridges. The model was also coupled with B-factor, weighted contact number, relative solvent accessibility, and conservation prescreening to determine the residues appropriate for the thermal adaptive design of salt bridges. According to our method, eight putative salt-bridges were designed on a mesophilic β-glucosidase and 24 variants were constructed to verify the predictions. Six putative salt-bridges leaded to the increase of the enzyme thermal stability. A significant increase in melting temperature of 8.8, 4.8, 3.7, 1.3, 1.2, and 0.7°C of the putative salt-bridges N437K–D49, E96R–D28, E96K–D28, S440K–E70, T231K–D388, and Q277E–D282 was detected, respectively. Reversing the polarity of T231K–D388 to T231D–D388K resulted in a further increase in melting temperatures by 3.6°C, which may be caused by the transformation of an intra-subunit electrostatic interaction into an inter-subunit one depending on the local environment. The combination of the thermostable variants (N437K, E96R, T231D and D388K) generated a melting temperature increase of 15.7°C. Thus, this study demonstrated a novel method for the thermal adaptive design of salt bridges through inference of suitable positions and substitutions.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Thermal adaptation of dihydrofolate reductase from the moderate thermophile Geobacillus stearothermophilus

The thermal melting temperature of dihydrofolate reductase from Geobacillus stearothermophilus (BsDHFR) is ~30 °C higher than that of its homologue from the psychrophile Moritella profunda. Additional proline residues in the loop regions of BsDHFR have been proposed to enhance the thermostability of BsDHFR, but site-directed mutagenesis studies reveal that these proline residues contribute only minimally. Instead, the high thermal stability of BsDHFR is partly due to removal of water-accessible thermolabile residues such as glutamine and methionine, which are prone to hydrolysis or oxidation at high temperatures. The extra thermostability of BsDHFR can be obtained by ligand binding, or in the presence of salts or cosolvents such as glycerol and sucrose. The sum of all these incremental factors allows BsDHFR to function efficiently in the natural habitat of G. stearothermophilus, which is characterized by temperatures that can reach 75 °C.

Comparative analysis of amino acid sequences from mesophiles and thermophiles in respective of carbon-nitrogen hydrolase family

A comparative study of amino acid sequence and physicochemical properties indicates the affiliation of protein from the nitrilase/cyanide hydratase family. This family contains nitrilases that break carbon-nitrogen bonds and appear to be involved in the reduction of organic nitrogen compounds and ammonia production. They all have distinct substrate specificity and include nitrilase, cyanide hydratases, aliphatic amidases, beta-alanine synthase, and a few other proteins with unknown molecular function. These sequences were analyzed for different physical and chemical properties and to relate these observed differences to the thermostability properties, phylogenetic tree construction and the evolutionary relationship among them. In this work, in silico analysis of amino acid sequences of mesophilic (15) and thermophilic (archaea, 15 and bacteria, 15) proteins has been done. The physiochemical properties of these three groups of nitrilase/cyanide hydratase family also differ in number of amino acids, molecular weight, pI values, positively charged ions, i.e. Arg + Lys, aliphatic index and grand average of hydropathacity (GRAVY). The amino acid Ala (1.37-fold) was found to be higher in mesophilic bacteria as compared to thermophilic bacteria but Lys and Phe were found to be significantly high (1.43 and 1.39-fold, respectively) in case of thermophilic bacteria. The amino acids Ala, Cys, Gln, His and Thr were found to be significantly higher (1.41, 1.6, 1.77, 1.44 and 1.29-fold, respectively) in mesophilic bacteria as compared to thermophilic archaea, where Glu, Leu and Val were found significantly high (1.22, 1.19 and 1.26-fold, respectively).

How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain

Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.

Sequence analysis and structure prediction of enoyl-CoA hydratase from Avicennia marina: implication of various amino acid residues on substrate-enzyme interactions

Enoyl-CoA hydratase catalyzes the hydration of 2-trans-enoyl-CoA into 3-hydroxyacyl-CoA. The present study focuses on the correlation between the functional and structural aspects of enoyl-CoA hydratase from Avicennia marina. We have used bioinformatics tools to construct and analyze 3D homology models of A. marina enoyl-CoA hydratase (AMECH) bound to different substrates and inhibitors and studied the residues involved in the ligand-enzyme interaction. Structural information obtained from the models was compared with those of the reported crystal structures. We observed that the overall folds were similar; however, AMECH showed few distinct structural changes which include structural variation in the mobile loop, formation and loss of certain interactions between the active site residues and substrates. Some changes were also observed within specific regions of the enzyme. Glu106 is almost completely conserved in sequences of the isomerases/hydratases including AMECH while Glu86 which is the other catalytic residue in most of the isomerases/hydratases is replaced by Gly and shows no interaction with the substrate. Asp114 is located within 4Å distance of the catalytic water which makes it a probable candidate for the second catalytic residue in AMECH. Another prominent feature of AMECH is the presence of structurally distinct mobile loop having a completely different coordination with the hydrophobic binding pocket of acyl portion of the substrate.

Molecular dynamics perspective on the protein thermal stability: a case study using SAICAR synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) Cα, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.

Importance of the ion-pair interactions in the OPEP coarse-grained force field: parametrization and validation

We have derived new effective interactions that improve the description of ion-pairs in the OPEP coarse-grained force field without introducing explicit electrostatic terms. The iterative Boltzmann inversion method was used to extract these potentials from all atom simulations by targeting the radial distribution function of the distance between the center of mass of the side-chains. The new potentials have been tested on several systems that differ in structural properties, thermodynamic stabilities and number of ion-pairs. Our modeling, by refining the packing of the charged amino-acids, impacts the stability of secondary structure motifs and the population of intermediate states during temperature folding/unfolding; it also improves the aggregation propensity of peptides. The new version of the OPEP force field has the potentiality to describe more realistically a large spectrum of situations where salt-bridges are key interactions.

The entropic nature of protein thermal stabilization

We performed thermodynamic analysis of temperature-induced unfolding of mesophilic and thermophilic proteins. It was shown that the variability in protein thermostability associated with pH-dependent unfolding or linked to the substitution of amino acid residues on the protein surface is evidence of the governing role of the entropy factor. Numerical values of conformational components in enthalpy, entropy and free energy which characterize protein unfolding in the "gas phase" were obtained. Based on the calculated absolute values of entropy and free energy, a model of protein unfolding is proposed in which the driving force is the conformational entropy of native protein, as an energy of the heat motion (T·S(NC)) increasing with temperature and acting as an factor devaluating the energy of intramolecular weak bonds in the transition state.

Isolation of a highly reactive β-sheet-rich intermediate of lysozyme in a solvent-free liquid phase

The thermal denaturation of solvent-free liquid lysozyme at temperatures in excess of 200 °C was studied by synchrotron radiation circular dichroism spectroscopy. Temperature-dependent changes in the secondary structure were used to map the equilibrium denaturation pathway and characterize a reactive β-sheet-rich unfolding intermediate that was stable in the solvent-free liquid phase under anhydrous conditions but which underwent irreversible aggregation in the presence of water. The unfolding intermediate had a transition temperature of 78 °C and was extremely stable to temperature, eventually reaching the fully denatured state at 178 °C. We propose that the three-stage denaturation pathway arises from the decreased stability of the native state due to the absence of any appreciable hydrophobic effect, along with an entropically derived stabilization of the reactive intermediate associated with molecular crowding in the solvent-free liquid.

Atomic structure of the 75 MDa extremophile Sulfolobus turreted icosahedral virus determined by CryoEM and X-ray crystallography

Sulfolobus turreted icosahedral virus (STIV) was isolated in acidic hot springs where it infects the archeon Sulfolobus solfataricus. We determined the STIV structure using near-atomic resolution electron microscopy and X-ray crystallography allowing tracing of structural polypeptide chains and visualization of transmembrane proteins embedded in the viral membrane. We propose that the vertex complexes orchestrate virion assembly by coordinating interactions of the membrane and various protein components involved. STIV shares the same coat subunit and penton base protein folds as some eukaryotic and bacterial viruses, suggesting that they derive from a common ancestor predating the divergence of the three kingdoms of life. One architectural motif (β-jelly roll fold) forms virtually the entire capsid (distributed in three different gene products), indicating that a single ancestral protein module may have been at the origin of its evolution.

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.

The effect of protein composition on hydration dynamics

Water dynamics at the surface of two homologous proteins with different thermal resistances is found to be unaffected by the different underlying amino-acid compositions, and when proteins are folded it responds similarly to temperature variations. Upon unfolding the water dynamics slowdown with respect to bulk decreases by a factor of two. Our findings are explained by the dominant topological perturbation induced by the protein on the water hydrogen bond dynamics.

Correlated electrostatic mutations provide a reservoir of stability in HIV protease

HIV protease, an aspartyl protease crucial to the life cycle of HIV, is the target of many drug development programs. Though many protease inhibitors are on the market, protease eventually evades these drugs by mutating at a rapid pace and building drug resistance. The drug resistance mutations, called primary mutations, are often destabilizing to the enzyme and this loss of stability has to be compensated for. Using a coarse-grained biophysical energy model together with statistical inference methods, we observe that accessory mutations of charged residues increase protein stability, playing a key role in compensating for destabilizing primary drug resistance mutations. Increased stability is intimately related to correlations between electrostatic mutations – uncorrelated mutations would strongly destabilize the enzyme. Additionally, statistical modeling indicates that the network of correlated electrostatic mutations has a simple topology and has evolved to minimize frustrated interactions. The model's statistical coupling parameters reflect this lack of frustration and strongly distinguish like-charge electrostatic interactions from unlike-charge interactions for of the most significantly correlated double mutants. Finally, we demonstrate that our model has considerable predictive power and can be used to predict complex mutation patterns, that have not yet been observed due to finite sample size effects, and which are likely to exist within the larger patient population whose virus has not yet been sequenced.

HIV is incurable because its enzymes evolve rapidly by developing resistance mutations to retroviral inhibitors. Most of these mutations work synergistically, but the biophysical basis behind their cooperation is not well understood. Our work addresses these important issues by bridging the gap between the statistical modeling of HIV protease subtype B sequences with the energetics of mutations involving charged amino acids by showing that electrostatic stability is intimately related to correlations. Moreover, we demonstrate that our statistical model has considerable predictive power and can be used to predict complex mutation patterns that have not yet been observed due to the finite sizes of the current sequence databases. In other words, as the database size increases, our model has the ability to predict the identities of the high probability mutations patterns, which are more likely to be observed. Knowing which currently unobserved mutations are more likely to be observed can be very advantageous in combating the disease.

Comparative analysis of zinc finger proteins involved in plant disease resistance

A meta-analysis was performed to understand the role of zinc finger domains in proteins of resistance (R) genes cloned from different crops. We analyzed protein sequences of seventy R genes of various crops in which twenty six proteins were found to have zinc finger domains along with nucleotide binding sites - leucine rice repeats (NBS-LRR) domains. We identified thirty four zinc finger domains in the R proteins of nine crops and were grouped into 19 types of zinc fingers. The size of individual zinc finger domain within the R genes varied from 11 to 84 amino acids, whereas the size of proteins containing these domains varied from 263 to 1305 amino acids. The biophysical analysis revealed that molecular weight of Pi54 zinc finger was lowest whereas the highest one was found in rice Pib zinc finger named as Transposes Transcription Factor (TTF). The instability (R(2) =0.95) and the aliphatic (R(2) =0.94) indices profile of zinc finger domains follows the polynomial distribution pattern. The pairwise identity analysis showed that the Lin11, Isl-1 & Mec-3 (LIM) zinc finger domain of rice blast resistance protein pi21 have 12.3% similarity with the nuclear transcription factor, X-box binding-like 1 (NFX) type zinc finger domain of Pi54 protein. For the first time, we reported that Pi54 (Pi-k(h)-Tetep), a rice blast resistance (R) protein have a small zinc finger domain of NFX type located on the C-terminal in between NBS and LRR domains of the R-protein. Compositional analysis depicted by the helical wheel diagram revealed the presence of a hydrophobic region within this domain which might help in exposing the LRR region for a possible R-Avr interaction. This domain is unique among all other cloned plant disease resistance genes and might play an important role in broad-spectrum nature of rice blast resistance gene Pi54.

Role of disulfide bonds in conformational stability and folding of 5'-deoxy-5'-methylthioadenosine phosphorylase II from the hyperthermophilic archaeon Sulfolobus solfataricus

Sulfolobus solfataricus 5'-deoxy-5'-melthylthioadenosine phosphorylase II (SsMTAPII), is a hyperthermophilic hexameric protein with two intrasubunit disulfide bonds (C138-C205 and C200-C262) and a CXC motif (C259-C261). To get information on the role played by these covalent links in stability and folding, the conformational stability of SsMTAPII and C262S and C259S/C261S mutants was studied by thermal and guanidinium chloride (GdmCl)-induced unfolding and analyzed by fluorescence spectroscopy, circular dichroism, and SDS-PAGE. No thermal unfolding transition of SsMTAPII can be obtained under nonreducing conditions, while in the presence of the reducing agent Tris-(2-carboxyethyl) phosphine (TCEP), a Tm of 100°C can be measured demonstrating the involvement of disulfide bridges in enzyme thermostability. Different from the wild-type, C262S and C259S/C261S show complete thermal denaturation curves with sigmoidal transitions centered at 102°C and 99°C respectively. Under reducing conditions these values decrease by 4°C and 8°C respectively, highlighting the important role exerted by the CXC disulfide on enzyme thermostability. The contribution of disulfide bonds to the conformational stability of SsMTAPII was further assessed by GdmCl-induced unfolding experiments carried out under reducing and nonreducing conditions. Thermal unfolding was found to be reversible if the protein was heated in the presence of TCEP up to 90°C but irreversible above this temperature because of aggregation. In analogy, only chemical unfolding carried out in the presence of reducing agents resulted in a reversible process suggesting that disulfide bonds play a role in enzyme denaturation. Thermal and chemical unfolding of SsMTAPII occur with dissociation of the native hexameric state into denatured monomers, as indicated by SDS-PAGE.

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.

Thermophilic proteins: insight and perspective from in silico experiments

Proteins from thermophilic and hyperthermophilic organisms are stable and function at high temperatures (50-100 °C). The importance of understanding the microscopic mechanisms underlying this thermal resistance is twofold: it is key for acquiring general clues on how proteins maintain their fold stable and for targeting those medical and industrial applications that aim at designing enzymes that can work under harsh conditions. In this tutorial review we first provide the general background of protein thermostability by specifically focusing on the structural and thermodynamic peculiarities; next, we discuss how computational studies based on Molecular Dynamics simulations can broaden and refine our knowledge on such special class of proteins.

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 inter-subunit interactions to the thermostability of Pyrococcus furiosus citrate synthase

Using citrate synthase from the hyperthermophile Pyrococcus furiosus (PfCS) as our test molecule, we show through guanidine hydrochloride-induced unfolding that the dimer separates into folded, but inactive, monomers before individual subunit unfolding takes place. Given that forces across the dimer interface are vital for thermostability, a robust computational method was derived that uses the University of Houston Brownian Dynamics (UHBD) program to calculate both the hydrophobic and electrostatic contribution to the dimerisation energy at 100°C. The results from computational and experimental determination of the lowered stability of interface mutants were correlated, being both of the same order of magnitude and placing the mutant proteins in the same order of stability. This computational method, optimised for hyperthermophilic molecules and tested in the laboratory, after further testing on other examples, could be of widespread use in the prediction of thermostabilising mutations in other oligomeric proteins for which dissociation is the first step in unfolding.

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.

Engineering the stability and the activity of a glycoside hydrolase

Glycosidases, the enzymes responsible in nature for the catabolism of carbohydrates, are well-studied catalysts widely used in industrial biotransformations and oligosaccharide synthesis, which are also attractive targets for drug development. Glycosidases from hyperthermophilic organisms (thriving at temperatures > 85 °C) are also interesting models to understand the molecular basis of protein stability and to produce robust tools for industrial applications. Here, we review the results obtained in the last two decades by our group on a β-glycosidase from the hyperthermophilic Archaeon Sulfolobus solfataricus. Our findings will be presented in the general context of the stability of proteins from hyperthermophiles and of the chemo-enzymatic synthesis of oligosaccharides.

Impact of distant charge reversals within a robust beta-barrel protein pore

Among all beta-barrel pores, staphylococcal alpha-hemolysin (alphaHL), a heptameric transmembrane protein of known high-resolution crystal structure, features a high stability in planar lipid bilayers under a wide range of harsh experimental conditions. Here, we employed single-channel electrical recordings and standard protein engineering to explore the impact of two distant charge reversals within the interior of the beta-barrel part of the pore. The charge reversals were replacements of lysines with aspartic acids. A charge reversal within the structurally stiff region of the beta barrel near the pore constriction reduced the open-state current of the pore, but produced a quiet pore, showing current fluctuation-free channel behavior. In contrast, a charge reversal on the trans entrance, within the structurally flexible glycine-rich turn of the beta barrel, increased the open-state current and produced gating activity of the pore in the form of large-amplitude and frequent current fluctuations. Remarkably, cumulative insertion of the two distant charge reversals resulted in a large-amplitude permanent blockade of the beta barrel, as judged by both single-channel and macroscopic current measurements. The results from this work suggest that these distant charge reversals are energetically coupled, producing different impacts on the ionic transport, the unitary conductance and the open-state probability of the pore.

Changes in non-core regions stabilise plastocyanin from the thermophilic cyanobacterium Phormidium laminosum

We report a theoretical investigation on the different stabilities of two plastocyanins. The first one belongs to the thermophilic cyanobacterium Phormidium laminosum and the second one belongs to its mesophilic relative Synechocystis sp. These proteins share the same topology and secondary-structure elements; however, the melting temperatures of their oxidised species differ by approximately 15 K. Long-time-scale molecular dynamics simulations, performed at different temperatures, show that the thermophilic protein optimises a set of intramolecular interactions (interstrand hydrogen bonding, salt bridging and hydrophobic clustering) within the region that comprises the strands beta 5 and beta 6, loop L5 and the helix. This region exhibits most of the differences in the primary sequence between the two proteins and, in addition, it is involved in the interaction with known physiological partners. Further work is in progress to unveil the specific structural features responsible for the different thermal stability of the two proteins.

Shape and evolution of thermostable protein structure

Organisms evolved at high temperatures must maintain their proteins' structures in the face of increased thermal disorder. This challenge results in differences in residue utilization and overall structure. Focusing on thermostable/mesostable pairs of homologous structures, we have examined these differences using novel geometric measures: specifically burial depth (distance from the molecular surface to each atom) and travel depth (distance from the convex hull to the molecular surface that avoids the protein interior). These along with common metrics like packing and Wadell Sphericity are used to gain insight into the constraints experienced by thermophiles. Mean travel depth of hyperthermostable proteins is significantly less than that of their mesostable counterparts, indicating smaller, less numerous and less deep pockets. The mean burial depth of hyperthermostable proteins is significantly higher than that of mesostable proteins indicating that they bury more atoms further from the surface. The burial depth can also be tracked on the individual residue level, adding a finer level of detail to the standard exposed surface area analysis. Hyperthermostable proteins for the first time are shown to be more spherical than their mesostable homologues, regardless of when and how they adapted to extreme temperature. Additionally, residue specific burial depth examinations reveal that charged residues stay unburied, most other residues are slightly more buried and Alanine is more significantly buried in hyperthermostable proteins.

Effective approach for calculations of absolute stability of proteins using focused dielectric constants

The ability to predict the absolute stability of proteins based on their corresponding sequence and structure is a problem of great fundamental and practical importance. In this work, we report an extensive, refinement and validation of our recent approach (Roca et al., FEBS Lett 2007;581:2065-2071) for predicting absolute values of protein stability DeltaG(fold). This approach employs the semimacroscopic protein dipole Langevin dipole method in its linear response approximation version (PDLD/S-LRA) while using the best fitted values of the dielectric constants epsilon'(p) and epsilon'(eff) for the self energy and charge-charge interactions, respectively. The method is validated on a diverse set of 45 proteins. It is found that the best fitted values of both dielectric constants are around 40. However, the self energy of internal residues and the charge-charge interactions of Lys have to be treated with care, using a somewhat lower values of epsilon'(p) and epsilon'(eff). The predictions of DeltaG(fold) reported here, have an average error of only 1.8 kcal/mole compared to the observed values, making our method very promising for estimating protein stability. It also provides valuable insight into the complex electrostatic phenomena taking place in folded proteins.

A novel scoring function for discriminating hyperthermophilic and mesophilic proteins with application to predicting relative thermostability of protein mutants

Background

The ability to design thermostable proteins is theoretically important and practically useful. Robust and accurate algorithms, however, remain elusive. One critical problem is the lack of reliable methods to estimate the relative thermostability of possible mutants.

Results

We report a novel scoring function for discriminating hyperthermophilic and mesophilic proteins with application to predicting the relative thermostability of protein mutants. The scoring function was developed based on an elaborate analysis of a set of features calculated or predicted from 540 pairs of hyperthermophilic and mesophilic protein ortholog sequences. It was constructed by a linear combination of ten important features identified by a feature ranking procedure based on the random forest classification algorithm. The weights of these features in the scoring function were fitted by a hill-climbing algorithm. This scoring function has shown an excellent ability to discriminate hyperthermophilic from mesophilic sequences. The prediction accuracies reached 98.9% and 97.3% in discriminating orthologous pairs in training and the holdout testing datasets, respectively. Moreover, the scoring function can distinguish non-homologous sequences with an accuracy of 88.4%. Additional blind tests using two datasets of experimentally investigated mutations demonstrated that the scoring function can be used to predict the relative thermostability of proteins and their mutants at very high accuracies (92.9% and 94.4%). We also developed an amino acid substitution preference matrix between mesophilic and hyperthermophilic proteins, which may be useful in designing more thermostable proteins.

Conclusions

We have presented a novel scoring function which can distinguish not only HP/MP ortholog pairs, but also non-homologous pairs at high accuracies. Most importantly, it can be used to accurately predict the relative stability of proteins and their mutants, as demonstrated in two blind tests. In addition, the residue substitution preference matrix assembled in this study may reflect the thermal adaptation induced substitution biases. A web server implementing the scoring function and the dataset used in this study are freely available at http://www.abl.ku.edu/thermorank/.

Inverse tuning of metal binding affinity and protein stability by altering charged coordination residues in designed calcium binding proteins

Ca(2+ )binding proteins are essential for regulating the role of Ca(2+ )in cell signaling and maintaining Ca(2+ )homeostasis. Negatively charged residues such as Asp and Glu are often found in Ca(2+ )binding proteins and are known to influence Ca(2+ )binding affinity and protein stability. In this paper, we report a systematic investigation of the role of local charge number and type of coordination residues in Ca(2+ )binding and protein stability using de novo designed Ca(2+ )binding proteins. The approach of de novo design was chosen to avoid the complications of cooperative binding and Ca(2+)-induced conformational change associated with natural proteins. We show that when the number of negatively charged coordination residues increased from 2 to 5 in a relatively restricted Ca(2+)-binding site, Ca(2+ )binding affinities increased by more than 3 orders of magnitude and metal selectivity for trivalent Ln(3+ )over divalent Ca(2+ )increased by more than 100-fold. Additionally, the thermal transition temperatures of the apo forms of the designed proteins decreased due to charge repulsion at the Ca(2+ )binding pocket. The thermal stability of the proteins was regained upon Ca(2+ )and Ln(3+ )binding to the designed Ca(2+ )binding pocket. We therefore observe a striking tradeoff between Ca(2+)/Ln(3+ )affinity and protein stability when the net charge of the coordination residues is varied. Our study has strong implications for understanding and predicting Ca(2+)-conferred thermal stabilization of natural Ca(2+ )binding proteins as well as for designing novel metalloproteins with tunable Ca(2+ )and Ln(3+ )binding affinity and selectivity.PACS codes: 05.10.-a.