Contribution of charged groups to the enthalpic stabilization of the folded states of globular proteins

Reverse Protonation of Buried Ion-Pairs in Staphylococcal Nuclease Mutants

Although buried titratable residues in protein cavities are often of major functional importance, it is generally challenging to understand their properties such as the ionization state and factors of stabilization based on experimental studies alone. A specific set of examples involve buried Glu-Lys pairs in a series of variants of Staphylococcal nuclease, for which recent structural and thermodynamic studies appeared to suggest that both the stability and the ionization state of the buried Glu-Lys pair are sensitive to its orientation (i.e., Glu23-Lys36 vs Lys23-Glu36). To further clarify the situation, especially ionization states of the buried Glu-Lys pairs, we have conducted extensive molecular dynamics simulations and free energy computations. Microsecond molecular dynamics simulations show that the hydration level of the cavity depends on the orientation of the buried ion-pair therein as well as its ionization state; free energy simulations recapitulate the relative stability of Glu23-Lys36 (EK) vs Lys23-Glu36 (KE) mutants measured experimentally, although the difference is similar in magnitude regardless of the ionization state of the Glu-Lys pair. A complementary set of free energy simulations strongly suggests that, in contrast to the original suggestion in the experimental analysis, the Glu and Lys residues prefer to adopt their charge-neutral rather than the ionized states. This result is consistent with the low dielectric constant computed for water in the cavity, which makes it difficult for the protein cavity to stabilize a pair of charged Glu-Lys residues, even with water penetration. The current study highlights the role of free energy simulations in understanding the ionization state of buried titratable residues and the relevant energetic contributions, forming the basis for the rational design of buried charge networks in proteins.

Evaluating the dynamics and electrostatic interactions of folded proteins in implicit solvents

Three implicit solvent models, namely GBMVII, FACTS, and SCPISM, were evaluated for their abilities to emulate an explicit solvent environment by comparing the simulated conformational ensembles, dynamics, and electrostatic interactions of the Src SH2 domain and the Lyn kinase domain. This assessment in terms of structural features in folded proteins expands upon the use of hydration energy as a metric for comparison. All-against-all rms coordinate deviation, average positional fluctuations, and ion-pair distance distribution were used to compare the implicit solvent models with the TIP3P explicit solvent model. Our study shows that the Src SH2 domains solvated with TIP3P, GBMVII, and FACTS sample similar global conformations. Additionally, the Src SH2 ion-pair distance distributions of solvent-exposed side chains corresponding to TIP3P, GBMVII, and FACTS do not differ substantially, indicating that GBMVII and FACTS are capable of modeling these electrostatic interactions. The ion-pair distance distributions of SCPISM are distinct from others, demonstrating that these electrostatic interactions are not adequately reproduced with the SCPISM model. On the other hand, for the Lyn kinase domain, a non-globular protein with bilobal structure and a large concavity on the surface, implicit solvent does not accurately model solvation to faithfully reproduce partially buried electrostatic interactions and lobe-lobe conformations. Our work reveals that local structure and dynamics of small, globular proteins are modeled well using FACTS and GBMVII. Nonetheless, global conformations and electrostatic interactions in concavities of multi-lobal proteins resulting from simulations with implicit solvent models do not match those obtained from explicit water simulations.

Role of key salt bridges in thermostability of G. thermodenitrificans EstGtA2: distinctive patterns within the new bacterial lipolytic enzyme subfamily XIII.2 [corrected]

Bacterial lipolytic enzymes were originally classified into eight different families defined by Arpigny and Jaeger (families I-VIII). Recently, the discovery of new lipolytic enzymes allowed for extending the original classification to fourteen families (I-XIV). We previously reported that G. thermodenitrificans EstGtA2 (access no. AEN92268) belonged to a novel group of bacterial lipolytic enzymes. Here we propose a 15(th) family (family XV) and suggest criteria for the assignation of protein sequences to the N' subfamily. Five selected salt bridges, hallmarks of the N' subfamily (E3/R54, E12/R37, E66/R140, D124/K178 and D205/R220) were disrupted in EstGtA2 using a combinatorial alanine-scanning approach. A set of 14 (R/K→A) mutants was produced, including five single, three double, three triple and three quadruple mutants. Despite a high tolerance to non-conservative mutations for folding, all the alanine substitutions were destabilizing (decreasing T m by 5 to 14°C). A particular combination of four substitutions exceeded this tolerance and prevents the correct folding of EstGtA2, leading to enzyme inactivation. Although other mutants remain active at low temperatures, the accumulation of more than two mutations had a dramatic impact on EstGtA2 activity at high temperatures suggesting an important role of these conserved salt bridge-forming residues in thermostability of lipolytic enzymes from the N' subfamily. We also identified a particular interloop salt bridge in EstGtA2 (D194/H222), located at position i -2 and i -4 residues from the catalytic Asp and His respectively which is conserved in other related bacterial lipolytic enzymes (families IV and XIII) with high tolerance to mutations and charge reversal. We investigated the role of residue identity at position 222 in controlling stability-pH dependence in EstGtA2. The introduction of a His to Arg mutation led to increase thermostability under alkaline pH. Our results suggest primary targets for optimization of EstGtA2 for specific biotechnological purposes.

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

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

Prediction of protein relative enthalpic stability from molecular dynamics simulations of the folded and unfolded states

For proteins of known structure, the relative enthalpic stability with respect to wild-type, ΔΔH(U), can be estimated by direct computation of the folded and unfolded state energies. We propose a model by which the change in stability upon mutation can be predicted from all-atom molecular dynamics simulations for the folded state and a peptide-based model for the unfolded state. The unfolding enthalpies are expressed in terms of environmental and hydration-solvent reorganization contributions that readily allow a residue-specific analysis of ΔΔH(U). The method is applied to estimate the relative enthalpic stability of variants with buried charged groups in T4 lysozyme. The predicted relative stabilities are in good agreement with experimental data. Environmental factors are observed to contribute more than hydration to the overall ΔΔH(U). The residue-specific analysis finds that the effects of burying charge are both localized and long-range. The enthalpy for hydration-solvent reorganization varies considerably among different amino-acid types, but because the variant folded state structures are similar to those of the wild-type, the hydration-solvent reorganization contribution to ΔΔH(U) is localized at the mutation site, in contrast to environmental contributions. Overall, mutation of apolar and polar amino acids to charged amino acids are destabilizing, but the reasons are complex and differ from site to site.