Small Neutral Crowding Solute Effects on Protein Folding Thermodynamic Stability and Kinetics.

Molecular crowding is a ubiquitous phenomenon in biological systems, with significant consequences on protein folding and stability. Small compounds, such as the osmolyte trimethylamine N-oxide (TMAO), can also present similar effects. To analyze the effects arising from these solute-like molecules, we performed a series of crowded coarse-grained folding simulations. Two well-known proteins were chosen, CI2 and SH3, modeled by the alpha-carbon-structure-based model. In the simulations, the crowding molecules were represented by low-sized neutral atom beads in different concentrations. The results show that a low level of the volume fraction occupied by neutral agents can change protein stability and folding kinetics for the two systems. However, the kinetics were shown to be unaffected in their respective folding temperatures. The faster kinetics correlates with changes in the folding route for systems at the same temperature, where non-native contacts were shown to be relevant. The transition states of the two systems with and without crowders are similar. In the case of SH3, there are differences in the structuring of two strands, which may be associated with the increase in its folding rate, in addition to the destabilization of the denatured ensemble. The present study also detected a crossover in the thermodynamic stability behavior, previously observed experimentally and theoretically. As the temperature increases, crowders change from destabilizing to stabilizing agents.

Modeling Chikungunya control strategies and Mayaro potential outbreak in the city of Rio de Janeiro.

Mosquito-borne diseases have become a significant health issue in many regions around the world. For tropical countries, diseases such as Dengue, Zika, and Chikungunya, became epidemic in the last decades. Health surveillance reports during this period were crucial in providing scientific-based information to guide decision making and resources allocation to control outbreaks. In this work, we perform data analysis of the last Chikungunya epidemics in the city of Rio de Janeiro by applying a compartmental mathematical model. Sensitivity analyses were performed in order to describe the contribution of each parameter to the outbreak incidence. We estimate the "basic reproduction number" for those outbreaks and predict the potential epidemic outbreak of the Mayaro virus. We also simulated several scenarios with different public interventions to decrease the number of infected people. Such scenarios should provide insights about possible strategies to control future outbreaks.

Rational Design of Chymotrypsin Inhibitor 2 by Optimizing Non-Native Interactions.

Rational design of proteins via mutagenesis is crucial for several biotechnological applications. A significant challenge of the computational strategies used to predict optimized mutations is to understand the influence of each amino acid during the folding process. In the present work, chymotrypsin inhibitor 2 (CI2) and several of its designed mutants have been simulated using a non-native hydrophobic and electrostatic potential as a structure-based Cα model. Through these simulations, we could identify the most critical folding stage to accelerate CI2 and also the charged residues responsible for providing its thermostability. The replacement of ionizable residues for hydrophobic ones tended to promote the formation of the CI2 secondary structure in the early transition state, which speeds up folding. However, this same replacement destabilized the native structure, and there was a decrease in the protein thermostability. Such a simple method proved to be capable of providing valuable information about thermodynamics and kinetics of CI2 and its mutations, thus being a fast alternative to the study of rational protein design.

Non-Native Cooperative Interactions Modulate Protein Folding Rates.

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.

TKSA-MC: A web server for rational mutation through the optimization of protein charge interactions.

The TKSAMC is a web server which calculates protein charge-charge interactions via the Tanford-Kirkwood Surface Accessibility model with the Monte Carlo method for sampling different protein protonation states. The optimization of charge-charge interactions via directed mutations has successfully enhanced the thermal stability of different proteins and could be a key to protein engineering improvement. The server presents the electrostatic free energy contribution of each polar-charged residue to the protein native state stability. Specific residues are suggested to be mutated for improving thermal stability. The choice of a residue is based on its fraction of side chain exposed to solvent and its positive free energy contribution, which tends to destabilize the protein native state. Any residue energy contribution can be shown as a function of pH condition. The web server is freely available at UNESP (São Paulo State University - DF/IBILCE): http://tksamc.df.ibilce.unesp.br and also on GitHub https://github.com/contessoto/tksamc.

NTL9 Folding at Constant pH: The Importance of Electrostatic Interaction and pH Dependence.

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.

Analyzing the effect of homogeneous frustration in protein folding.

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.