Electrostatic stabilization of a thermophilic cold shock protein

Insights on Rigidity and Flexibility at the Global and Local Levels of Protein Structures and Their Roles in Homologous Psychrophilic, Mesophilic, and Thermophilic Proteins: A Computational Study

The rigidity and flexibility of homologous psychrophilic (P), mesophilic (M), and thermophilic (T) proteins have been investigated at the global and local levels in terms of "packing factors" and "atomic fluctuations" obtained from B-factors. For comparison of atomic fluctuations, correction of errors by considering errors in B-factors from all sources in a consolidated manner and conversion of the fluctuations to the same temperature have been suggested and validated. The results indicate no differences in the global values like the average packing factor among the three classes of protein homologues, but at local levels there are differences. A comparison of homologous protein triplets show that the average atomic fluctuations at a given temperature mainly obey the order P > M > T. Packing factors and the atomic fluctuations are anti-correlated, suggesting that altering the rigidity of the active site might be a potential strategy to make tailor-made psychrophilic or thermophilic proteins from their mesophilic homologues. The computer codes developed and used in this work are available at the link https://github.com/Munna-Sarkar/proteins-rigidity-flexibility.git.

Cold-Shock Domains-Abundance, Structure, Properties, and Nucleic-Acid Binding

Simple Summary

Proteins are composed of compact domains, often of known three-dimensional structure, and natively unstructured polypeptide regions. The abundant cold-shock domain is among the set of canonical nucleic acid-binding domains and conserved from bacteria to man. Proteins containing cold-shock domains serve a large variety of biological functions, which are mostly linked to DNA or RNA binding. These functions include the regulation of transcription, RNA splicing, translation, stability and sequestration. Cold-shock domains have a simple architecture with a conserved surface ideally suited to bind single-stranded nucleic acids. Because the binding is mostly by non-specific molecular interactions which do not involve the sugar-phosphate backbone, cold-shock domains are not strictly sequence-specific and do not discriminate reliably between DNA and RNA. Many, but not all functions of cold shock-domain proteins in health and disease can be understood based of the physical and structural properties of their cold-shock domains.

Abstract

The cold-shock domain has a deceptively simple architecture but supports a complex biology. It is conserved from bacteria to man and has representatives in all kingdoms of life. Bacterial cold-shock proteins consist of a single cold-shock domain and some, but not all are induced by cold shock. Cold-shock domains in human proteins are often associated with natively unfolded protein segments and more rarely with other folded domains. Cold-shock proteins and domains share a five-stranded all-antiparallel β-barrel structure and a conserved surface that binds single-stranded nucleic acids, predominantly by stacking interactions between nucleobases and aromatic protein sidechains. This conserved binding mode explains the cold-shock domains’ ability to associate with both DNA and RNA strands and their limited sequence selectivity. The promiscuous DNA and RNA binding provides a rationale for the ability of cold-shock domain-containing proteins to function in transcription regulation and DNA-damage repair as well as in regulating splicing, translation, mRNA stability and RNA sequestration.

The Role of Electrostatics and Folding Kinetics on the Thermostability of Homologous Cold Shock Proteins

Understanding which aspects contribute to the thermostability of proteins is a challenge that has persisted for decades, and it is of great relevance for protein engineering. Several types of interactions can influence the thermostability of a protein. Among them, the electrostatic interactions have been a target of particular attention. Aiming to explore how this type of interaction can affect protein thermostability, this paper investigated four homologous cold shock proteins from psychrophilic, mesophilic, thermophilic, and hyperthermophilic organisms using a set of theoretical methodologies. It is well-known that electrostatics as well as hydrophobicity are key-elements for the stabilization of these proteins. Therefore, both interactions were initially analyzed in the native structure of each protein. Electrostatic interactions present in the native structures were calculated with the Tanford-Kirkwood model with solvent accessibility, and the amount of hydrophobic surface area buried upon folding was estimated by measuring both folded and extended structures. On the basis of Energy Landscape Theory, the local frustration and the simplified alpha-carbon structure-based model were modeled with a Debye-Hückel potential to take into account the electrostatics and the effects of an implicit solvent. Thermodynamic data for the structure-based model simulations were collected and analyzed using the Weighted Histogram Analysis and Stochastic Diffusion methods. Kinetic quantities including folding times, transition path times, folding routes, and Φ values were also obtained. As a result, we found that the methods are able to qualitatively infer that electrostatic interactions play an important role on the stabilization of the most stable thermophilic cold shock proteins, showing agreement with the experimental data.

Influence of prior pH and thermal stresses on thermal tolerance of foodborne pathogens

Improper food processing is one of the major causes of foodborne illness. Accurate prediction of the thermal destruction rate of foodborne pathogens is therefore vital to ensure proper processing and food safety. When bacteria are subjected to pH and thermal stresses during growth, sublethal stresses can occur that may lead to differences in their subsequent tolerance to thermal treatment. As a preliminary study to test this concept, the current study evaluated the effect of prior pH and thermal stresses on thermal tolerance of Salmonella and Staphylococcus using a tryptic soy broth supplemented with yeast extract. Bacteria incubated at three pH values (6.0, 7.4, and 9.0) and four temperatures (15, 25, 35, and 45°C) for 24 hr were subjected to thermal treatments at 55, 60, and 65°C. At the end of each treatment time, bacterial suspensions were surface-plated on standard method agar for quantification of bacterial survival and further calculation of the thermal death decimal reduction time (D-value) and thermal destruction temperature (z-value). The effect of pH stress alone during the incubation on the thermal tolerance of both bacteria was generally insignificant. An increasing pattern of D-value was observed with the increment of thermal stress (incubation temperature). The bacteria incubated at 35°C required the highest z-value to reduce the 90% in D-values. Staphylococcus mostly displayed higher tolerance to thermal treatment than Salmonella. Although further research is needed to validate the current findings on food matrices, findings in this study clearly affirm that adaptation of bacteria to certain stresses may reduce the effectiveness of preservation procedures applied during later stage of food processing and storage.

Implementation of an experimental and computational tool set to study protein-mAb interactions

This work focused on the development of a combined experimental and computational tool set to study protein-mAb interactions. A model protein library was first screened using cross interaction chromatography to identify proteins with the strongest retention. Fluorescence polarization was then employed to study the interactions and thermodynamics of the selected proteins-lactoferrin, pyruvate kinase, and ribonuclease B with the mAb. Binding affinities of lactoferrin and pyruvate kinase to the mAb were seen to be relatively salt insensitive in the range examined. Further, a strong entropic contribution was observed, suggesting the importance of hydrophobic interactions. On the other hand, ribonuclease B-mAb binding was seen to be enthalpically driven and salt sensitive, indicating the importance of electrostatic interactions. Protein-protein docking was then carried out and the results identified the CDR region on the mAb as an important binding site for all three proteins. The binding interfaces identified for the mAb-lactoferrin and mAb-pyruvate kinase systems were found to contain complementary hydrophobic and oppositely charged clusters on the interacting regions which were indicative of both hydrophobic and electrostatic interactions. On the other hand, the binding site on ribonuclease B was predominantly positively charged with minimal hydrophobicity. This resulted in an alignment with negatively charged clusters on the mAb, supporting the contention that these interactions were primarily electrostatic in nature. Importantly, these computational results were found to be consistent with the fluorescence polarization data and this combined approach may have utility in examining mAb-HCP interactions which can often complicate the downstream processing of biologics. © 2019 American Institute of Chemical Engineers.

Biophysical and structural characterization of the thermostable WD40 domain of a prokaryotic protein, Thermomonospora curvata PkwA

WD40 proteins belong to a big protein family with members identified in every eukaryotic proteome. However, WD40 proteins were only reported in a few prokaryotic proteomes. Using WDSP (http://wu.scbb.pkusz.edu.cn/wdsp/), a prediction tool, we identified thousands of prokaryotic WD40 proteins, among which few proteins have been biochemically characterized. As shown in our previous bioinformatics study, a large proportion of prokaryotic WD40 proteins have higher intramolecular sequence identity among repeats and more hydrogen networks, which may indicate better stability than eukaryotic WD40s. Here we report our biophysical and structural study on the WD40 domain of PkwA from Thermomonospora curvata (referred as tPkwA-C). We demonstrated that the stability of thermophilic tPkwA-C correlated to ionic strength and tPkwA-C exhibited fully reversible unfolding under different denaturing conditions. Therefore, the folding kinetics was also studied through stopped-flow circular dichroism spectra. The crystal structure of tPkwA-C was further resolved and shed light on the key factors that stabilize its beta-propeller structure. Like other WD40 proteins, DHSW tetrad has a significant impact on the stability of tPkwA-C. Considering its unique features, we proposed that tPkwA-C should be a great structural template for protein engineering to study key residues involved in protein-protein interaction of a WD40 protein.

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.

Mechanics of Protein Adaptation to High Temperatures

Inspired by Somero's corresponding state principle that relates protein enhanced thermal stability with mechanical rigidity, we deployed state of the art computational techniques (based on atomistic steered molecular dynamics and Hamiltonian-replica exchange simulations) to study the in silico realization of mechanical and thermal unfolding of two homologous Csp proteins that have evolved to thrive in different thermal environments. By complementing recent single-molecule experiments, we unambiguously show that, for these homologues whose structures are very similar, the increased thermal resistance of the thermophilic variant is not associated with an increased mechanical stability. Our approach provides microscopic insights that are otherwise inaccessible to experimental techniques, and explains why the protein weak spots for thermal and mechanical denaturation are distinct.

Unusual dimerization of a BcCsp mutant leads to reduced conformational dynamics

Cold shock proteins (Csp) constitute a family of ubiquitous small proteins that act as RNA-chaperones to avoid cold-induced termination of translation. All members contain two subdomains composed of 2 and 3 β-strands, respectively, which are connected by a hinge loop and fold into a β-barrel. Bacillus caldolyticus Csp (BcCsp) is one of the most studied members of the family in terms of its folding, function, and structure. This protein has been described as a monomer in solution, although a recent crystal structure showed dimerization via domain swapping (DS). In contrast, other cold shock proteins of the same fold are known to dimerize in a nonswapped arrangement. Hypothesizing that reducing the size of the hinge loop may promote swapping as in several other DS proteins with different folds we deleted two residues from these region (BcCsp∆36-37), leading to a protein in monomer-dimer equilibrium with similar folding stability to that of the wild-type. Strikingly, the crystal structure of BcCsp∆36-37 revealed a nonswapped dimer with its interface located at the nucleic acid-binding surface, showing that the deletion led to structural consequences far from the perturbation site. Concomitantly, circular dichroism experiments on BcCsp∆36-37 demonstrated that binding of the oligonucleotide hexathymidine disrupts the dimer. Additionally, HDXMS shows a protective effect on the protein structure upon dimerization, where the resulting interactions between ligand-binding surfaces in the dimer reduced the extent of exchange throughout the whole protein. Our work provides evidence of the complex interplay between conformational dynamics, deletions, and oligomerization within the Csp protein family.

DATABASES

Structural data are available in the Protein Data Bank under accession number 5JX4.

Thermodynamics of globular proteins

The analysis of temperature-induced unfolding of proteins in aqueous solutions was performed. Based on the data of thermodynamic parameters of protein unfolding and using the method of semi-empirical calculations of hydration parameters at reference temperature 298 K, we obtained numerical values of enthalpy, free energy, and entropy which characterize the unfolding of proteins in the 'gas phase'. It was shown that specific values of the energy of weak intramolecular bonds (∆H^int), conformational free energy (∆G^conf) and entropy (∆S^conf) are the same for proteins with molecular weight 7-25 kDa. Using the energy value (∆H^int) and the proposed approach for estimation of the conformational entropy of native protein (S^NC), numerical values of the absolute free energy (G^NC) were obtained.

Simple approach for ranking structure determining residues

Mutating residues has been a common task in order to study structural properties of the protein of interest. Here, we propose and validate a simple method that allows the identification of structural determinants; i.e., residues essential for preservation of the stability of global structure, regardless of the protein topology. This method evaluates all of the residues in a 3D structure of a given globular protein by ranking them according to their connectivity and movement restrictions without topology constraints. Our results matched up with sequence-based predictors that look up for intrinsically disordered segments, suggesting that protein disorder can also be described with the proposed methodology.

A theoretical method to compute sequence dependent configurational properties in charged polymers and proteins

A general formalism to compute configurational properties of proteins and other heteropolymers with an arbitrary sequence of charges and non-uniform excluded volume interaction is presented. A variational approach is utilized to predict average distance between any two monomers in the chain. The presented analytical model, for the first time, explicitly incorporates the role of sequence charge distribution to determine relative sizes between two sequences that vary not only in total charge composition but also in charge decoration (even when charge composition is fixed). Furthermore, the formalism is general enough to allow variation in excluded volume interactions between two monomers. Model predictions are benchmarked against the all-atom Monte Carlo studies of Das and Pappu [Proc. Natl. Acad. Sci. U. S. A. 110, 13392 (2013)] for 30 different synthetic sequences of polyampholytes. These sequences possess an equal number of glutamic acid (E) and lysine (K) residues but differ in the patterning within the sequence. Without any fit parameter, the model captures the strong sequence dependence of the simulated values of the radius of gyration with a correlation coefficient of R(2) = 0.9. The model is then applied to real proteins to compare the unfolded state dimensions of 540 orthologous pairs of thermophilic and mesophilic proteins. The excluded volume parameters are assumed similar under denatured conditions, and only electrostatic effects encoded in the sequence are accounted for. With these assumptions, thermophilic proteins are found-with high statistical significance-to have more compact disordered ensemble compared to their mesophilic counterparts. The method presented here, due to its analytical nature, is capable of making such high throughput analysis of multiple proteins and will have broad applications in proteomic studies as well as in other heteropolymeric systems.

Mutagenesis of Trichoderma reesei endoglucanase I: impact of expression host on activity and stability at elevated temperatures

BACKGROUND

Trichoderma reesei is a key cellulase source for economically saccharifying cellulosic biomass for the production of biofuels. Lignocellulose hydrolysis at temperatures above the optimum temperature of T. reesei cellulases (~50°C) could provide many significant advantages, including reduced viscosity at high-solids loadings, lower risk of microbial contamination during saccharification, greater compatibility with high-temperature biomass pretreatment, and faster rates of hydrolysis. These potential advantages motivate efforts to engineer T. reesei cellulases that can hydrolyze lignocellulose at temperatures ranging from 60-70°C.

RESULTS

A B-factor guided approach for improving thermostability was used to engineer variants of endoglucanase I (Cel7B) from T. reesei (TrEGI) that are able to hydrolyze cellulosic substrates more rapidly than the recombinant wild-type TrEGI at temperatures ranging from 50-70°C. When expressed in T. reesei, TrEGI variant G230A/D113S/D115T (G230A/D113S/D115T Tr_TrEGI) had a higher apparent melting temperature (3°C increase in Tm) and improved half-life at 60°C (t1/2 = 161 hr) than the recombinant (T. reesei host) wild-type TrEGI (t1/2 = 74 hr at 60°C, Tr_TrEGI). Furthermore, G230A/D113S/D115T Tr_TrEGI showed 2-fold improved activity compared to Tr_TrEGI at 65°C on solid cellulosic substrates, and was as efficient in hydrolyzing cellulose at 60°C as Tr_TrEGI was at 50°C. The activities and stabilities of the recombinant TrEGI enzymes followed similar trends but differed significantly in magnitude depending on the expression host (Escherichia coli cell-free, Saccharomyces cerevisiae, Neurospora crassa, or T. reesei). Compared to N.crassa-expressed TrEGI, S. cerevisiae-expressed TrEGI showed inferior activity and stability, which was attributed to the lack of cyclization of the N-terminal glutamine in Sc_TrEGI and not to differences in glycosylation. N-terminal pyroglutamate formation in TrEGI expressed in S. cerevisiae was found to be essential in elevating its activity and stability to levels similar to the T. reesei or N. crassa-expressed enzyme, highlighting the importance of this ubiquitous modification in GH7 enzymes.

CONCLUSION

Structure-guided evolution of T. reesei EGI was used to engineer enzymes with increased thermal stability and activity on solid cellulosic substrates. Production of TrEGI enzymes in four hosts highlighted the impact of the expression host and the role of N-terminal pyroglutamate formation on the activity and stability of TrEGI enzymes.

Evolutionary aspects of enzyme dynamics

The role of evolutionary pressure on the chemical step catalyzed by enzymes is somewhat enigmatic, in part because chemistry is not rate-limiting for many optimized systems. Herein, we present studies that examine various aspects of the evolutionary relationship between protein dynamics and the chemical step in two paradigmatic enzyme families, dihydrofolate reductases and alcohol dehydrogenases. Molecular details of both convergent and divergent evolution are beginning to emerge. The findings suggest that protein dynamics across an entire enzyme can play a role in adaptation to differing physiological conditions. The growing tool kit of kinetics, kinetic isotope effects, molecular biology, biophysics, and bioinformatics provides means to link evolutionary changes in structure-dynamics function to the vibrational and conformational states of each protein.

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.

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.

Computational tools for rational protein engineering of aldolases

In this mini-review we describe the different strategies for rational protein engineering and summarize the computational tools available. Computational tools can either be used to design focused libraries, to predict sequence-function relationships or for structure-based molecular modelling. This also includes de novo design of enzymes. Examples for protein engineering of aldolases and transaldolases are given in the second part of the mini-review.

Interaction between soluble and membrane-embedded potassium channel peptides monitored by Fourier transform infrared spectroscopy

Recent studies have explored the utility of Fourier transform infrared spectroscopy (FTIR) in dynamic monitoring of soluble protein-protein interactions. Here, we investigated the applicability of FTIR to detect interaction between synthetic soluble and phospholipid-embedded peptides corresponding to, respectively, a voltage-gated potassium (Kv) channel inactivation domain (ID) and S4–S6 of the Shaker Kv channel (KV1; including the S4–S5 linker “pre-inactivation” ID binding site). KV1 was predominantly α-helical at 30°C when incorporated into dimyristoyl-l-α-phosphatidylcholine (DMPC) bilayers. Cooling to induce a shift in DMPC from liquid crystalline to gel phase reversibly decreased KV1 helicity, and was previously shown to partially extrude a synthetic S4 peptide. While no interaction was detected in liquid crystalline DMPC, upon cooling to induce the DMPC gel phase a reversible amide I peak (1633 cm⁻¹) consistent with novel hydrogen bond formation was detected. This spectral shift was not observed for KV1 in the absence of ID (or vice versa), nor when the non-inactivating mutant V7E ID was applied to KV1 under similar conditions. Alteration of salt or redox conditions affected KV1-ID hydrogen bonding in a manner suggesting electrostatic KV1-ID interaction favored by a hairpin conformation for the ID and requiring extrusion of one or more KV1 domains from DMPC, consistent with ID binding to S4–S5. These findings support the utility of FTIR in detecting reversible interactions between soluble and membrane-embedded proteins, with lipid state-sensitivity of the conformation of the latter facilitating control of the interaction.

Bluues server: electrostatic properties of wild-type and mutated protein structures

MOTIVATION

Electrostatic calculations are an important tool for deciphering many functional mechanisms in proteins. Generalized Born (GB) models offer a fast and convenient computational approximation over other implicit solvent-based electrostatic models. Here we present a novel GB-based web server, using the program Bluues, to calculate numerous electrostatic features including pKa-values and surface potentials. The output is organized allowing both experts and beginners to rapidly sift the data. A novel feature of the Bluues server is that it explicitly allows to find electrostatic differences between wild-type and mutant structures.

AVAILABILITY

The Bluues server, examples and extensive help files are available for non-commercial use at URL: http://protein.bio.unipd.it/bluues/.

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.

Relative stabilities of conserved and non-conserved structures in the OB-fold superfamily

The OB-fold is a diverse structure superfamily based on a beta-barrel motif that is often supplemented with additional non-conserved secondary structures. Previous deletion mutagenesis and NMR hydrogen exchange studies of three OB-fold proteins showed that the structural stabilities of sites within the conserved beta-barrels were larger than sites in non-conserved segments. In this work we examined a database of 80 representative domain structures currently classified as OB-folds, to establish the basis of this effect. Residue-specific values were obtained for the number of Calpha-Calpha distance contacts, sequence hydrophobicities, crystallographic B-factors, and theoretical B-factors calculated from a Gaussian Network Model. All four parameters point to a larger average flexibility for the non-conserved structures compared to the conserved beta-barrels. The theoretical B-factors and contact densities show the highest sensitivity. Our results suggest a model of protein structure evolution in which novel structural features develop at the periphery of conserved motifs. Core residues are more resistant to structural changes during evolution since their substitution would disrupt a larger number of interactions. Similar factors are likely to account for the differences in stability to unfolding between conserved and non-conserved structures.

Protein stabilization by the rational design of surface charge-charge interactions

The design of proteins with increased stability has many important applications in biotechnology. In recent years, strategies involving directed evolution, sequence-based design, or computational design have proven successful for generating stabilized proteins. A brief overview of the various methods that have been used to increase protein stability is presented, followed by a detailed example of how the rational design of surface charge-charge interactions has provided a robust method for protein stabilization.

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

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

Destabilization of psychrotrophic RNase HI in a localized fashion as revealed by mutational and X-ray crystallographic analyses

The Arg97 --> Gly and Asp136 --> His mutations stabilized So-RNase HI from the psychrotrophic bacterium Shewanella oneidensis MR-1 by 5.4 and 9.7 degrees C, respectively, in T(m), and 3.5 and 6.1 kJ x mol(-1), respectively, in DeltaG(H2O). These mutations also stabilized the So-RNase HI derivative (4x-RNase HI) with quadruple thermostabilizing mutations in an additive manner. As a result, the resultant sextuple mutant protein (6x-RNase HI) was more stable than the wild-type protein by 28.8 degrees C in T(m) and 27.0 kJ x mol(-1) in DeltaG(H2O). To analyse the effects of the mutations on the protein structure, the crystal structure of the 6x-RNase HI protein was determined at 2.5 A resolution. The main chain fold and interactions of the side-chains of the 6x-RNase HI protein were basically identical to those of the wild-type protein, except for the mutation sites. These results indicate that all six mutations independently affect the protein structure, and are consistent with the fact that the thermostabilizing effects of the mutations are roughly additive. The introduction of favourable interactions and the elimination of unfavourable interactions by the mutations contribute to the stabilization of the 6x-RNase HI protein. We propose that So-RNase HI is destabilized when compared with its mesophilic and thermophilic counterparts in a localized fashion by increasing the number of amino acid residues unfavourable for protein stability.

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.

Cooperativity of complex salt bridges

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

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

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

Computational methods for biomolecular electrostatics

An understanding of intermolecular interactions is essential for insight into how cells develop, operate, communicate, and control their activities. Such interactions include several components: contributions from linear, angular, and torsional forces in covalent bonds, van der waals forces, as well as electrostatics. Among the various components of molecular interactions, electrostatics are of special importance because of their long range and their influence on polar or charged molecules, including water, aqueous ions, and amino or nucleic acids, which are some of the primary components of living systems. Electrostatics, therefore, play important roles in determining the structure, motion, and function of a wide range of biological molecules. This chapter presents a brief overview of electrostatic interactions in cellular systems, with a particular focus on how computational tools can be used to investigate these types of interactions.

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

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

Folding and stability of the isolated Greek key domains of the long-lived human lens proteins gammaD-crystallin and gammaS-crystallin

The transparency of the eye lens depends on the high solubility and stability of the lens crystallin proteins. The monomeric gamma-crystallins and oligomeric beta-crystallins have paired homologous double Greek key domains, presumably evolved through gene duplication and fusion. Prior investigation of the refolding of human gammaD-crystallin revealed that the C-terminal domain folds first and nucleates the folding of the N-terminal domain. This result suggested that the human N-terminal domain might not be able to fold on its own. We constructed and expressed polypeptide chains corresponding to the isolated N- and C-terminal domains of human gammaD-crystallin, as well as the isolated domains of human gammaS-crystallin. Both circular dichroism and fluorescence spectroscopy indicated that the isolated domains purified from Escherichia coli were folded into native-like monomers. After denaturation, the isolated domains refolded efficiently at pH 7 and 37 degrees C into native-like structures. The in vitro refolding of all four domains revealed two kinetic phases, identifying partially folded intermediates for the Greek key motifs. When subjected to thermal denaturation, the isolated N-terminal domains were less stable than the full-length proteins and less stable than the C-terminal domains, and this was confirmed in equilibrium unfolding/refolding experiments. The decrease in stability of the N-terminal domain of human gammaD-crystallin with respect to the complete protein indicated that the interdomain interface contributes of 4.2 kcal/mol to the overall stability of this very long-lived protein.

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

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

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.

A proteomic approach toward the selection of proteins with enhanced intrinsic conformational stability

A detailed understanding of the molecular basis of protein folding and stability determinants partly relies on the study of proteins with enhanced conformational stability properties, such as those from thermophilic organisms. In this study, we set up a methodology aiming at identifying the subset of cytosolic hyperstable proteins using Sulfurispharea sp., a hyperthermophilic archaeon, able to grow between 70 and 97 degrees C, as a model organism. We have thermally and chemically perturbed the cytosolic proteome as a function of time (up to 96 h incubation at 90 degrees C), and proceeded with analysis of the remaining proteins by combining one- and two-dimensional gel electrophoresis, liquid chromatography fractionation, and protein identification by N-terminal sequencing and mass spectrometry methods. In total, 14 proteins with enhanced stabilities which are involved in key cellular processes such as detoxification, nucleic acid processing, and energy metabolism were identified including a superoxide dismutase, a peroxiredoxin, and a ferredoxin. We demonstrate that these proteins are biologically active after extensive thermal treatment of the proteome. The relevance of these and other targets is discussed in terms of the organism's ecology. This work thus illustrates an experimental approach aimed at mining a proteome for hyperstable proteins, a valuable tool for target selection in protein stability and structural studies.

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

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

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.

The Mechanism of TC230′s Thermostability: A Molecular Dynamics Simulation Study

The quasielastic neutron scattering index beta and the modulus of a protein's quasi-electric dipole moment were utilized to quantitate the thermostability of wildtype TC23O and its mutants. Charged residues Arg314, Glu246, Glu291, and some prolines near the C-terminus of the sequence (Pro228, Pro296, and Pro308) were identified to be critical for the thermostability of wildtype TC23O according to these two criteria. By analyzing the molecular conformation changes during the simulation, it was demonstrated how the mutant P228S was destabilized by disrupting two salt-bridges Asp116OD1-Lys215N and Glu210OE1-Lys217N at an adjacent beta-turn. The destabilization of P296S also shown to be intimate correlated with the break down of ion pair Lys188N-Glu291OE1. The sensitivity of its electrostatic network to the local structure is an important feature. It reveals that the 'proline effect' and electrostatic interactions together influences the thermostability of TC23O a lot.