The importance of charge-charge interactions in the thermal stability of proteins is widely known. pH and ionic strength play a crucial role in these electrostatic interactions, as well as in the arrangement of ionizable residues in each protein-folding stage. In this study, two coarse-grained models were used to evaluate the effect of pH and salt concentration on the thermal stability of a protein G variant (1PGB-QDD), which was chosen due to the quantity of experimental data exploring these effects on its stability. One of these coarse-grained models, the TKSA, calculates the electrostatic free energy of the protein in the native state via the Tanford-Kirkwood approach for each residue. The other one, CpHMD-SBM, uses a Coulomb screening potential in addition to the structure-based model C. Both models simulate the system in constant pH. The comparison between the experimental stability analysis and the computational results obtained by these simple models showed a good agreement. Through the TKSA method, the role of each charged residue in the protein's thermal stability was inferred. Using CpHMD-SBM, it was possible to evaluate salt and pH effects throughout the folding process. Finally, the computational pK values were calculated by both methods and presented a good level of agreement with the experiments. This study provides, to our knowledge, new information and a comprehensive description of the electrostatic contribution to protein G stability.
The folding process of the N-terminal domain of ribosomal protein L9 (NTL9) was investigated at constant-pH computer simulations. Evaluation of the role of electrostatic interaction during folding was carried out by including a Debye-Hückel potential into a Cα structure-based model (SBM). In this study, the charges of the ionizable residues and the electrostatic potential are susceptible to the solution conditions, such as pH and ionic strength, as well as to the presence of charged groups. Simulations were performed under different pHs, and the results were validated by comparing them with experimental values of pKa and with denaturation experiment data. Also, the free energy profiles, Φ-values, and folding routes were calculated for each condition. It was shown how charges vary along the folding under different pH, which is subject to different scenarios. This study reveals how simplified models can capture essential physical features, reproducing experimental results, and presenting the role of electrostatic interactions before, during, and after the transition state.
The energy landscape theory and the funnel description have had remarkable success in describing protein folding mechanisms and function. However, there are experimental results that are not understood using this approach. Among the puzzling examples are the α-spectrin results, in which the R15 domain folds 3 orders of magnitude more rapidly than the homologous R16 and R17, even though they are structurally very similar to each other. Such anomalous observations are usually attributed to the influence of internal friction on protein folding rates, but this is not a satisfactory explanation. In this study, this phenomenon is addressed by focusing on non-native interactions that could account for this effect. We carried out molecular dynamics simulations with structure-based C α models, in which the folding process of α-spectrin domains was investigated. The simulations take into account the hydrophobic and electrostatic contributions separately. The folding time results have shown qualitative agreement with the experimental data. We have also investigated mutations in R16 and R17, and the simulation folding time results correlate with the observed experimental ones. We suggest that the origin of the internal friction, at least in this case, might emerge from a cooperativity effect of these non-native interactions.
Constant pH molecular dynamics (MD) simulations sample protonation states on the fly according to the conformational environment and user specified pH conditions; however, the current accuracy is limited due to the use of implicitsolvent models or a hybrid solvent scheme. Here, we report the first GPU-accelerated implementation, parametrization, and validation of the all-atom continuous constant pH MD (CpHMD) method with particle-mesh Ewald (PME) electrostatics in the Amber22 pmemd.cuda engine. The titration parameters for Asp, Glu, His, Cys, and Lys were derived for the CHARMM c22 and Amber ff14sb and ff19sb force fields. We then evaluated the PME-CpHMD method using the asynchronous pH replicaexchange titration simulations with the c22 force field for six benchmark proteins, including BBL, hen egg white lysozyme (HEWL), staphylococcal nuclease (SNase), thioredoxin, ribonuclease A (RNaseA), and human muscle creatine kinase (HMCK). The root-mean-square deviation from the experimental pK a 's of Asp, Glu, His, and Cys is 0.76 pH units, and the Pearson's correlation coefficient for the pK a shifts with respect to model values is 0.80. We demonstrated that a finite-size correction or much enlarged simulation box size can remove a systematic error of the calculated pK a 's and improve agreement with experiment. Importantly, the simulations captured the relevant biology in several challenging cases, e.g., the titration order of the catalytic dyad Glu35/Asp52 in HEWL and the coupled residues Asp19/Asp21 in SNase, the large pK a upshift of the deeply buried catalytic Asp26 in thioredoxin, and the large pK a downshift of the deeply buried catalytic Cys283 in HMCK. We anticipate that PME-CpHMD will offer proper pH control to improve the accuracies of MD simulations and enable mechanistic studies of proton-coupled dynamical processes that are ubiquitous in biology but remain poorly understood due to the lack of experimental tools and limitation of current MD simulations.
The influence of soft, hydration-mediated ion-ion and ion-surface interactions on the differential capacitance of an electric double layer is investigated using Monte Carlo simulations and compared to various mean-field models. We focus on a planar electrode surface at physiological concentration of monovalent ions in a uniform dielectric background. Hydration-mediated interactions are modeled on the basis of Yukawa potentials that add to the Coulomb and excluded volume interactions between ions. We present a mean-field model that includes hydration-mediated anion-anion, anion-cation, and cation-cation interactions of arbitrary strengths. In addition, finite ion sizes are accounted for through excluded volume interactions, described either on the basis of the Carnahan-Starling equation of state or using a lattice gas model. Both our Monte Carlo simulations and mean-field approaches predict a characteristic double-peak (the so-called camel shape) of the differential capacitance; its decrease reflects the packing of the counterions near the electrode surface. The presence of hydration-mediated ion-surface repulsion causes a thin charge-depleted region close to the surface, which is reminiscent of a Stern layer. We analyze the interplay between excluded volume and hydration-mediated interactions on the differential capacitance and demonstrate that for small surface charge density our mean-field model based on the Carnahan-Starling equation is able to capture the Monte Carlo simulation results. In contrast, for large surface charge density the mean-field approach based on the lattice gas model is preferable.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.