Stabilizing salt-bridge enhances protein thermostability by reducing the heat capacity change of unfolding

Establishing knowledge on the sequence arrangement pattern of nucleated protein folding

The heat-tolerance mechanisms of (hyper)thermophilic proteins provide a unique opportunity to investigate the unsolved protein folding problem. In an attempt to determine whether the interval between residues in sequence might play a role in determining thermostability, we constructed a sequence interval-dependent value function to calculate the residue pair frequency. Additionally, we identified a new sequence arrangement pattern, where like-charged residues tend to be adjacently assembled, while unlike-charged residues are distributed over longer intervals, using statistical analysis of a large sequence database. This finding indicated that increasing the intervals between unlike-charged residues can increase protein thermostability, with the arrangement patterns of these charged residues serving as thermodynamically favorable nucleation points for protein folding. Additionally, we identified that the residue pairs K-E, R-E, L-V and V-V involving long sequence intervals play important roles involving increased protein thermostability. This work demonstrated a novel approach for considering sequence intervals as keys to understanding protein folding. Our findings of novel relationships between residue arrangement and protein thermostability can be used in industry and academia to aid the design of thermostable proteins.

In Silico Designing of an Industrially Sustainable Carbonic Anhydrase Using Molecular Dynamics Simulation

Carbonic anhydrase (CA) is a family of metalloenzymes that has the potential to sequestrate carbon dioxide (CO₂) from the environment and reduce pollution. The goal of this study is to apply protein engineering to develop a modified CA enzyme that has both higher stability and activity and hence could be used for industrial purposes. In the current study, we have developed an in silico method to understand the molecular basis behind the stability of CA. We have performed comparative molecular dynamics simulation of two homologous α-CA, one of thermophilic origin (Sulfurihydrogenibium sp.) and its mesophilic counterpart (Neisseria gonorrhoeae), for 100 ns each at 300, 350, 400, and 500 K. Comparing the trajectories of two proteins using different stability-determining factors, we have designed a highly thermostable version of mesophilic α-CA by introducing three mutations (S44R, S139E, and K168R). The designed mutant α-CA maintains conformational stability at high temperatures. This study shows the potential to develop industrially stable variants of enzymes while maintaining high activity.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Electrostatic effects on the folding stability of FKBP12

The roles of electrostatic interactions in protein folding stability have been a matter of debate, largely due to the complexity in the theoretical treatment of these interactions. We have developed computational methods for calculating electrostatic effects on protein folding stability. To rigorously test and further refine these methods, here we carried out experimental studies into electrostatic effects on the folding stability of the human 12-kD FK506 binding protein (FKBP12). This protein has a close homologue, FKBP12.6, with amino acid substitutions in only 18 of their 107 residues. Of the 18 substitutions, 8 involve charged residues. Upon mutating FKBP12 residues at these 8 positions individually into the counterparts in FKBP12.6, the unfolding free energy (ΔGu) of FKBP12 changed by -0.3 to 0.7 kcal/mol. Accumulating stabilizing substitutions resulted in a mutant with a 0.9 kcal/mol increase in stability. Additional charge mutations were grafted from a thermophilic homologue, MtFKBP17, which aligns to FKBP12 with 31% sequence identity over 89 positions. Eleven such charge mutations were studied, with ΔΔGu varying from -2.9 to 0.1 kcal/mol. The predicted electrostatic effects by our computational methods with refinements herein had a root-mean-square deviation of 0.9 kcal/mol from the experimental ΔΔGu values on 16 single mutations of FKBP12. The difference in ΔΔGu between mutations grafted from FKBP12.6 and those from MtFKBP17 suggests that more distant homologues are less able to provide guidance for enhancing folding stability.

Proposed ionic bond between Arg300 and Glu270 and Glu271 are not involved in inactivation of a mutant firefly luciferase (LRR)

The weakness of firefly luciferase is its rapid inactivation. Many studies have been done to develop thermostable luciferases. One of these modifications was LRR mutant in which the Leu300 was substituted with Arg in the E(354)RR(356)Lampyris turkestanicus luciferase as template. LRR was more thermostable than the wild type but with only 0.02% activity. In this study, site-directed mutagenesis was used to change the proposed ionic bond between the Arg and two neighboring residues (Glu270 and Glu271), to understand if the induced interactions were responsible for inactivation in LRR. Our results showed that substitution of Glu270 and 271 with Ala removed the interactions but the activity of enzyme did not return. The E270A mutant was more active than LRR but the E271A and E270A/E271A mutants were inactive. Fluorescence and CD measurements showed that these mutations were accompanied by conformational changes. Extrinsic fluorescence measurement and obtained quenching data by KI and acrylamide also confirmed that the mutants were less compact than the LRR enzyme. In conclusion, in LRR, the interactions between Arg300 and Glu270 and Glu271 were not responsible for the enzyme inactivation and it is proposed that the enzyme inactivation is due to conformational changes of LRR mutant of firefly luciferase.

Structural features determining thermal adaptation of esterases

The adaptation of microorganisms to extreme living temperatures requires the evolution of enzymes with a high catalytic efficiency under these conditions. Such extremophilic enzymes represent valuable tools to study the relationship between protein stability, dynamics and function. Nevertheless, the multiple effects of temperature on the structure and function of enzymes are still poorly understood at the molecular level. Our analysis of four homologous esterases isolated from bacteria living at temperatures ranging from 10°C to 70°C suggested an adaptation route for the modulation of protein thermal properties through the optimization of local flexibility at the protein surface. While the biochemical properties of the recombinant esterases are conserved, their thermal properties have evolved to resemble those of the respective bacterial habitats. Molecular dynamics simulations at temperatures around the optimal temperatures for enzyme catalysis revealed temperature-dependent flexibility of four surface-exposed loops. While the flexibility of some loops increased with raising the temperature and decreased with lowering the temperature, as expected for those loops contributing to the protein stability, other loops showed an increment of flexibility upon lowering and raising the temperature. Preserved flexibility in these regions seems to be important for proper enzyme function. The structural differences of these four loops, distant from the active site, are substantially larger than for the overall protein structure, indicating that amino acid exchanges within these loops occurred more frequently thereby allowing the bacteria to tune atomic interactions for different temperature requirements without interfering with the overall enzyme function.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

Exploring Strong Interactions in Proteins with Quantum Chemistry and Examples of Their Applications in Drug Design

OBJECTIVES

Three strong interactions between amino acid side chains (salt bridge, cation-π, and amide bridge) are studied that are stronger than (or comparable to) the common hydrogen bond interactions, and play important roles in protein-protein interactions.

METHODS

Quantum chemical methods MP2 and CCSD(T) are used in calculations of interaction energies and structural optimizations.

RESULTS

The energies of three types of amino acid side chain interactions in gaseous phase and in aqueous solutions are calculated using high level quantum chemical methods and basis sets. Typical examples of amino acid salt bridge, cation-π, and amide bridge interactions are analyzed, including the inhibitor design targeting neuraminidase (NA) enzyme of influenza A virus, and the ligand binding interactions in the HCV p7 ion channel. The inhibition mechanism of the M2 proton channel in the influenza A virus is analyzed based on strong amino acid interactions.

CONCLUSION

(1) The salt bridge interactions between acidic amino acids (Glu- and Asp-) and alkaline amino acids (Arg+, Lys+ and His+) are the strongest residue-residue interactions. However, this type of interaction may be weakened by solvation effects and broken by lower pH conditions. (2) The cation- interactions between protonated amino acids (Arg+, Lys+ and His+) and aromatic amino acids (Phe, Tyr, Trp and His) are 2.5 to 5-fold stronger than common hydrogen bond interactions and are less affected by the solvation environment. (3) The amide bridge interactions between the two amide-containing amino acids (Asn and Gln) are three times stronger than hydrogen bond interactions, which are less influenced by the pH of the solution. (4) Ten of the twenty natural amino acids are involved in salt bridge, or cation-, or amide bridge interactions that often play important roles in protein-protein, protein-peptide, protein-ligand, and protein-DNA interactions.

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.

Thermostability of in vitro evolved Bacillus subtilis lipase A: a network and dynamics perspective

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure.

Molecular dynamics studies on the NMR structures of rabbit prion protein wild type and mutants: surface electrostatic charge distributions

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species such as sheep and goats, cattle, deer and elk, and humans. But for rabbits, studies have shown that they have a low susceptibility to be infected by prion diseases. This paper does molecular dynamics (MD) studies of rabbit NMR structures (of the wild type and its two mutants of two surface residues), in order to understand the specific mechanism of rabbit prion proteins (RaPrP(C)). Protein surface electrostatic charge distributions are specially focused to analyze the MD trajectories. This paper can conclude that surface electrostatic charge distributions indeed contribute to the structural stability of wild-type RaPrP(C); this may be useful for the medicinal treatment of prion diseases.

Using first passage statistics to extract environmentally dependent amino acid correlations

In this work, we study the first passage statistics of amino acid primary sequences, that is the probability of observing an amino acid for the first time at a certain number of residues away from a fixed amino acid. By using this rich mathematical framework, we are able to capture the background distribution for an organism, and infer lengths at which the first passage has a probability that differs from what is expected. While many features of an organism's genome are due to natural selection, others are related to amino acid chemistry and the environment in which an organism lives, constraining the randomness of genomes upon which selection can further act. We therefore use this approach to infer amino acid correlations, and then study how these correlations vary across a wide range of organisms under a wide range of optimal growth temperatures. We find a nearly universal exponential background distribution, consistent with the idea that most amino acids are globally uncorrelated from other amino acids in genomes. When we are able to extract significant correlations, these correlations are reliably dependent on optimal growth temperature, across phylogenetic boundaries. Some of the correlations we extract, such as the enhanced probability of finding, for the first time, a cysteine three residues away from a cysteine or glutamic acid two residues away from an arginine, likely relate to thermal stability. However, other correlations, likely appearing on alpha helical surfaces, have a less clear physiochemical interpretation and may relate to thermal stability or unusual metabolic properties of organisms that live in a high temperature environment.

Structure of a His170Tyr mutant of thermostable pNPPase from Geobacillus stearothermophilus

Using directed evolution based on random mutagenesis and heat-treated selection, a thermostable His170Tyr mutant of Geobacillus stearothermophilus thermostable p-nitrophenylphosphatase (TpNPPase) was obtained. The temperature at which the His170Tyr mutant lost 50% of its activity (T1/2) was found to be 4.40 K higher than that of wild-type TpNPPase, and the melting temperature of the His170Tyr mutant increased by 2.39 K. The crystal structure of the His170Tyr mutant was then determined at 2.0 Å resolution in the presence of a sodium ion and a sulfate ion in the active site. The cap domain of chain B shows a half-closed conformation. The hydrophobic side chain of the mutated residue, the hydroxyphenyl group, forms a hydrophobic contact with the methyl group of Ala166. This hydrophobic interaction was found using the Protein Interactions Calculator (PIC) web server with an interaction distance of 4.6 Å, and might be a key factor in the thermostabilization of the His170Tyr mutant. This study potentially offers a molecular basis for both investigation of the catalytic mechanism and thermostable protein engineering.

Three-dimensional structure of an alkaline xylanase Xyn11A-LC from alkalophilic Bacillus sp. SN5 and improvement of its thermal performance by introducing arginines substitutions

The alkaline xylanase Xyn11A-LC from the alkalophilic Bacillus sp. SN5 was expressed in E. coli, purified and crystallized. The crystal structure was determined at a resolution of 1.49 Å. Xyn11A-LC has the β-jelly roll structure typical of family 11 xylanases. To improve its thermostability and thermophilicity, a mutant SB3 was constructed by introducing three arginines on the different sides of the protein surface. SB3 increased the optimum temperature by 5 °C. The wild type and SB3 had the half-lives of 22 and 68 min at 65 °C at pH 8.0 (Tris/HCl buffer), respectively. CD spectroscopy revealed that the melting temperature (T m) of the wild type and SB3 were 55.3 and 66.9 °C, respectively. These results showed that the introduction of arginines enhance the thermophilicity and thermostability of Xyn11A-LC.

SBION: A Program for Analyses of Salt-Bridges from Multiple Structure Files

Salt-bridge and network salt-bridge are specific electrostatic interactions that contribute to the overall stability of proteins. In hierarchical protein folding model, these interactions play crucial role in nucleation process. The advent and growth of protein structure database and its availability in public domain made an urgent need for context dependent rapid analysis of salt-bridges. While these analyses on single protein is cumbersome and time-consuming, batch analyses need efficient software for rapid topological scan of a large number of protein for extracting details on (i) fraction of salt-bridge residues (acidic and basic). (ii) Chain specific intra-molecular salt-bridges, (iii) inter-molecular salt-bridges (protein-protein interactions) in all possible binary combinations (iv) network salt-bridges and (v) secondary structure distribution of salt-bridge residues. To the best of our knowledge, such efficient software is not available in public domain. At this juncture, we have developed a program i.e. SBION which can perform all the above mentioned computations for any number of protein with any number of chain at any given distance of ion-pair. It is highly efficient, fast, error-free and user friendly. Finally we would say that our SBION indeed possesses potential for applications in the field of structural and comparative bioinformatics studies.

Expression and properties of the highly alkalophilic phenylalanine ammonia-lyase of thermophilic Rubrobacter xylanophilus

The sequence of a phenylalanine ammonia-lyase (PAL; EC: 4.3.1.24) of the thermophilic and radiotolerant bacterium Rubrobacter xylanophilus (RxPAL) was identified by screening the genomes of bacteria for members of the phenylalanine ammonia-lyase family. A synthetic gene encoding the RxPAL protein was cloned and overexpressed in Escherichia coli TOP 10 in a soluble form with an N-terminal His6-tag and the recombinant RxPAL protein was purified by Ni-NTA affinity chromatography. The activity assay of RxPAL with l-phenylalanine at various pH values exhibited a local maximum at pH 8.5 and a global maximum at pH 11.5. Circular dichroism (CD) studies showed that RxPAL is associated with an extensive α-helical character (far UV CD) and two distinctive near-UV CD peaks. These structural characteristics were well preserved up to pH 11.0. The extremely high pH optimum of RxPAL can be rationalized by a three-dimensional homology model indicating possible disulfide bridges, extensive salt-bridge formation and an excess of negative electrostatic potential on the surface. Due to these properties, RxPAL may be a candidate as biocatalyst in synthetic biotransformations leading to unnatural l- or d-amino acids or as therapeutic enzyme in treatment of phenylketonuria or leukemia.

The conserved lid tryptophan, W211, potentiates thermostability and thermoactivity in bacterial thermoalkalophilic lipases

We hypothesize that aggregation of thermoalkalophilic lipases could be a thermostability mechanism. The conserved tryptophans (W211, W234) in the lid are of particular interest owing to their previous involvements in aggregation and thermostability mechanisms in many other proteins. The thermoalkalophilic lipase from Bacillus thermocatenulatus (BTL2) and its mutants (W211A, W234A) were expressed and purified to homogeneity. We found that, when aggregated, BTL2 is more thermostable than its non-aggregating form, showing that aggregation potentiates thermostability in the thermoalkalophilic lipase. Among the two lid mutants, the W211A lowered aggregation tendency drastically and resulted in a much less thermostable variant of BTL2, which indicated that W211 stabilizes the intermolecular interactions in BTL2 aggregates. Further thermoactivity and CD spectroscopy analyses showed that W211A also led to a strong decrease in the optimal and the melting temperature of BTL2, implying stabilization by W211 also to the intramolecular interactions. The other lid mutant W234A had no effects on these properties. Finally, we analyzed the molecular basis of these experimental findings in-silico using the dimer (PDB ID: 1KU0) and the monomer (PDB ID: 2W22) lipase structures. The computational analyses confirmed that W211 stabilized the intermolecular interactions in the dimer lipase and it is critical to the stability of the monomer lipase. Explicitly W211 confers stability to the dimer and the monomer lipase through distinct aromatic interactions with Y273-Y282 and H87-P232 respectively. The insights revealed by this work shed light not only on the mechanism of thermostability and its relation to aggregation but also on the particular role of the conserved lid tryptophan in the thermoalkalophilic lipases.

Protein adaptations in archaeal extremophiles

Extremophiles, especially those in Archaea, have a myriad of adaptations that keep their cellular proteins stable and active under the extreme conditions in which they live. Rather than having one basic set of adaptations that works for all environments, Archaea have evolved separate protein features that are customized for each environment. We categorized the Archaea into three general groups to describe what is known about their protein adaptations: thermophilic, psychrophilic, and halophilic. Thermophilic proteins tend to have a prominent hydrophobic core and increased electrostatic interactions to maintain activity at high temperatures. Psychrophilic proteins have a reduced hydrophobic core and a less charged protein surface to maintain flexibility and activity under cold temperatures. Halophilic proteins are characterized by increased negative surface charge due to increased acidic amino acid content and peptide insertions, which compensates for the extreme ionic conditions. While acidophiles, alkaliphiles, and piezophiles are their own class of Archaea, their protein adaptations toward pH and pressure are less discernible. By understanding the protein adaptations used by archaeal extremophiles, we hope to be able to engineer and utilize proteins for industrial, environmental, and biotechnological applications where function in extreme conditions is required for activity.

Engineering and kinetic stabilization of the therapeutic enzyme Anabeana variabilis phenylalanine ammonia lyase

Anabeana variabilis phenylalanine ammonia lyase has just recently been discovered and introduced in clinical trials of phenylketonuria enzyme replacement therapy for its outstanding kinetic properties. In the present study, kinetic stabilization of this therapeutically important enzyme has been explored by introduction of a disulfide bond into the structure. Site-directed mutagenesis was performed with quick-change PCR method. Recombinant wild-type and mutated enzymes were expressed in Escherichia coli, and his-tagged proteins were affinity purified. Formation of disulfide bond was confirmed by Ellman's method, and then chemical unfolding, kinetic behavior, and thermal inactivation of mutated enzyme were compared with the wild type. Based on our results, the Q292C mutation resulted in a significant improvement in kinetic stability and resistance against chemical unfolding of the enzyme while kinetic parameters and pH profile of enzyme activity were remained unaffected. The results of the present study provided an insight towards designing phenylalanine ammonia lyases with higher stability.

A multi-factors rational design strategy for enhancing the thermostability of Escherichia coli AppA phytase

Despite recent advances in our understanding of the importance of protein surface properties for protein thermostability,there are seldom studies on multi-factors rational design strategy, so a more scientific, simple and effective rational strategy is urgent for protein engineering. Here, we first attempted to use a three-factors rational design strategy combining three common structural features, protein flexibility, protein surface, and salt bridges. Escherichia coli AppA phytase was used as a model enzyme to improve its thermostability. Moreover, the structure and enzyme features of the thermostable mutants designed by our strategy were analyzed roundly. For the single mutants, two (Q206E and Y311K), in five exhibited thermostable property with a higher success rate of prediction (40 %). For the multiple mutants, the themostable sites were combined with another site, I427L, we obtained by directed evolution, Q206E/I427L, Y311K/I427L, and Q206E/Y311K/I427L, all exhibited thermostable property. The Y311K/I427L doubled thermostability (61.7 %, and was compared to 30.97 % after being heated at 80 °C for 10 min) and catalytic efficiency (4.46 was compared to 2.37) improved more than the wild-type AppA phytase almost without hampering catalytic activity. These multi-factors of rational design strategy can be applied practically as a thermostabilization strategy instead of the conventional single-factor approach.

Stability mechanisms of a thermophilic laccase probed by molecular dynamics

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T_opt (∼350 K) from 300 K, this change correlated with a beginning “unzipping” of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F⁻ (but not Cl⁻) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F⁻ intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes.

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.

Residue-specific annotation of disorder-to-order transition and cathepsin inhibition of a propeptide-like crammer from D. melanogaster

Drosophila melanogaster crammer is a novel cathepsin inhibitor involved in long-term memory formation. A molten globule-to-ordered structure transition is required for cathepsin inhibition. This study reports the use of alanine scanning to probe the critical residues in the two hydrophobic cores and the salt bridges of crammer in the context of disorder-to-order transition and cathepsin inhibition. Alanine substitution of the aromatic residues W9, Y12, F16, Y20, Y32, and W53 within the hydrophobic cores, and charged residues E8, R28, R29, and E67 in the salt bridges considerably decrease the ability of crammer to inhibit Drosophila cathepsin B (CTSB). Far-UV circular dichroism (CD), intrinsic fluorescence, and nuclear magnetic resonance (NMR) spectroscopies show that removing most of the aromatic and charged side-chains substantially reduces thermostability, alters pH-dependent helix formation, and disrupts the molten globule-to-ordered structure transition. Molecular modeling indicates that W53 in the hydrophobic Core 2 is essential for the interaction between crammer and the prosegment binding loop (PBL) of CTSB; the salt bridge between R28 and E67 is critical for the appropriate alignment of the α-helix 4 toward the CTSB active cleft. The results of this study show detailed residue-specific dissection of folding transition and functional contributions of the hydrophobic cores and salt bridges in crammer, which have hitherto not been characterized for cathepsin inhibition by propeptide-like cysteine protease inhibitors. Because of the involvements of cathepsin inhibitors in neurodegenerative diseases, these structural insights can serve as a template for further development of therapeutic inhibitors against human cathepsins.

Relationship between Escherichia coli AppA phytase's thermostability and salt bridges

In order to study on the relationship between Escherichia coli AppA phytase's thermostability and salt bridges, and indicate an effective technical route of which factor to think about and where to modify at AppA for enhancing its thermostability, a salt bridge subtraction mutant E31Q and a salt bridge addition mutant Q307D were constructed by site-directed mutagenesis. The residual activities of the wild-type AppA phytase, E31Q and Q307D were 31.42%, 17.46%, and 40.57%, respectively, after being heated at 80°C for 10 min. The salt bridge subtraction mutant E31Q showed 13.96% thermostability decreasement, and the salt bridge addition mutant Q307D showed 9.15% thermostability enhancement than the wild-type both without the pH and temperature optimum changed. It proved salt bridges play a key role in E. coli AppA phytase's thermostability and the α/β-domain of AppA may be sensitive to heat. Salt bridges and the α/β-domain of AppA should have high priority to think about to enhance AppA's thermostability for commercial application. Besides, molecular dynamics simulation was used for salt bridges analysis.

Promiscuous contacts and heightened dynamics increase thermostability in an engineered variant of the engrailed homeodomain

A thermostabilized variant (UVF) of the engrailed homeodomain (EnHD) was previously engineered by Mayo and co-workers. The melting temperature of the non-natural, designed protein is 50°C higher than the natural wild-type protein (>99 vs. 52°C), and the two proteins share 22% sequence identity. We have performed extensive (1 μs) all-atom, explicit solvent molecular dynamics simulations of the wild-type and engineered proteins to investigate their structural and dynamic properties at room temperature and at 100°C. Our simulations are in good agreement with nuclear magnetic resonance data available for the two proteins [nuclear Overhauser effect crosspeaks (NOEs), J-coupling constants and order parameters for EnHD; and NOEs for UVF], showing that we reproduce the backbone dynamics and side chain packing in the native state of both proteins. UVF was more dynamic at room temperature than EnHD, with respect to both its backbone and side chain motion. When the temperature was raised, the thermostable protein maintained this mobility while retaining its native conformation. EnHD, on the other hand, was unable to maintain its more rigid native structure at higher temperature and began to unfold. Heightened protein dynamics leading to promiscuous and dynamically interchangeable amino acid contacts makes UVF more tolerant to increasing temperature, providing a molecular explanation for heightened thermostability of this protein.

A study on the effect of surface lysine to arginine mutagenesis on protein stability and structure using green fluorescent protein

Two positively charged basic amino acids, arginine and lysine, are mostly exposed to protein surface, and play important roles in protein stability by forming electrostatic interactions. In particular, the guanidinium group of arginine allows interactions in three possible directions, which enables arginine to form a larger number of electrostatic interactions compared to lysine. The higher pKa of the basic residue in arginine may also generate more stable ionic interactions than lysine. This paper reports an investigation whether the advantageous properties of arginine over lysine can be utilized to enhance protein stability. A variant of green fluorescent protein (GFP) was created by mutating the maximum possible number of lysine residues on the surface to arginines while retaining the activity. When the stability of the variant was examined under a range of denaturing conditions, the variant was relatively more stable compared to control GFP in the presence of chemical denaturants such as urea, alkaline pH and ionic detergents, but the thermal stability of the protein was not changed. The modeled structure of the variant indicated putative new salt bridges and hydrogen bond interactions that help improve the rigidity of the protein against different chemical denaturants. Structural analyses of the electrostatic interactions also confirmed that the geometric properties of the guanidinium group in arginine had such effects. On the other hand, the altered electrostatic interactions induced by the mutagenesis of surface lysines to arginines adversely affected protein folding, which decreased the productivity of the functional form of the variant. These results suggest that the surface lysine mutagenesis to arginines can be considered one of the parameters in protein stability engineering.

Highly conserved salt bridge stabilizes rigid signal patch at extracellular loop critical for surface expression of acid-sensing ion channels

Acid-sensing ion channels (ASICs) are non-selective cation channels activated by extracellular acidosis associated with many physiological and pathological conditions. A detailed understanding of the mechanisms that govern cell surface expression of ASICs, therefore, is critical for better understanding of the cell signaling under acidosis conditions. In this study, we examined the role of a highly conserved salt bridge residing at the extracellular loop of rat ASIC3 (Asp(107)-Arg(153)) and human ASIC1a (Asp(107)-Arg(160)) channels. Comprehensive mutagenesis and electrophysiological recordings revealed that the salt bridge is essential for functional expression of ASICs in a pH sensing-independent manner. Surface biotinylation and immunolabeling of an extracellular epitope indicated that mutations, including even minor alterations, at the salt bridge impaired cell surface expression of ASICs. Molecular dynamics simulations, normal mode analysis, and further mutagenesis studies suggested a high stability and structural constrain of the salt bridge, which serves to separate an adjacent structurally rigid signal patch, important for surface expression, from a flexible gating domain. Thus, we provide the first evidence of structural requirement that involves a stabilizing salt bridge and an exposed rigid signal patch at the destined extracellular loop for normal surface expression of ASICs. These findings will allow evaluation of new strategies aimed at preventing excessive excitability and neuronal injury associated with tissue acidosis and ASIC activation.

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

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

Solution structure of the dimerization domain of the eukaryotic stalk P1/P2 complex reveals the structural organization of eukaryotic stalk complex

The lateral ribosomal stalk is responsible for the kingdom-specific binding of translation factors and activation of GTP hydrolysis during protein synthesis. The eukaryotic stalk is composed of three acidic ribosomal proteins P0, P1 and P2. P0 binds two copies of P1/P2 hetero-dimers to form a pentameric P-complex. The structure of the eukaryotic stalk is currently not known. To provide a better understanding on the structural organization of eukaryotic stalk, we have determined the solution structure of the N-terminal dimerization domain (NTD) of P1/P2 hetero-dimer. Helix-1, -2 and -4 from each of the NTD-P1 and NTD-P2 form the dimeric interface that buries 2200 A² of solvent accessible surface area. In contrast to the symmetric P2 homo-dimer, P1/P2 hetero-dimer is asymmetric. Three conserved hydrophobic residues on the surface of NTD-P1 are replaced by charged residues in NTD-P2. Moreover, NTD-P1 has an extra turn in helix-1, which forms extensive intermolecular interactions with helix-1 and -4 of NTD-P2. Truncation of this extra turn of P1 abolished the formation of P1/P2 hetero-dimer. Systematic truncation studies suggest that P0 contains two spine-helices that each binds one copy of P1/P2 hetero-dimer. Modeling studies suggest that a large hydrophobic cavity, which can accommodate the loop between the spine-helices of P0, can be found on NTD-P1 but not on NTD-P2 when the helix-4 adopts an ‘open’ conformation. Based on the asymmetric properties of NTD-P1/NTD-P2, a structural model of the eukaryotic P-complex with P2/P1:P1/P2 topology is proposed.