Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins

Conformational Dynamics, Energetics, and the Divergent Evolution of Allosteric Regulation: The Case of the Yeast MAPK Family

Allosteric mechanisms provide finely-tuned control over signalling proteins. Proteins of the same family may share high sequence identity and structural similarity but show distinct traits of allosteric control and evolutionary divergent regulation. Revealing the determinants of such properties may be important to understand the molecular bases of different regulatory pathways. Herein, we investigate whether and how evolutionarily-divergent traits of allosteric regulation in homologous proteins can be decoded in terms of internal dynamics and interaction networks that support functionally oriented conformations. In this framework, we start from the comparative analysis of the dynamics and energetics of the yeast MAP Kinases (MAPKs) Fus3 and Kss1 in their native basins. Importantly, distinctive dynamic and energetic stabilization features emerge, which can be related to the two proteins' differential ability to be phosphorylated and engage with the allosteric activator Ste5. We then expanded our study to other evolutionarily-related MAPKs. We show that the dynamical and energetical traits defining the distinct regulatory profiles of Fus3 and Kss1 can be traced along their evolutionary tree. Overall, our approach is able to reconnect (latent) allostery with the principal elements of protein structural stabilization and dynamics, showing how allosteric regulation was encrypted in MAPKs structure well before Ste5 appearance.

Effect of Main and Side Chains on the Folding Mechanism of the Trp-Cage Miniprotein

Proteins that do not fold into their functional native state have been linked to diseases. In this study, the influence of the main and side chains of individual amino acids on the folding of the tryptophan cage (Trp-cage), a designed 20-residue miniprotein, was analyzed. For this purpose, we calculated the solvation free energy (SFE) contributions of individual atoms by using the 3D-reference interaction site model with the atomic decomposition method. The mechanism by which the Trp-cage is stabilized during the folding process was examined by calculating the total energy, which is the sum of the conformational energy and SFE. The folding process of the Trp-cage resulted in a stable native state, with a total energy that was 62.4 kcal/mol lower than that of the unfolded state. The solvation entropy, which is considered to be responsible for the hydrophobic effect, contributed 31.3 kcal/mol to structural stabilization. In other words, the contribution of the solvation entropy accounted for approximately half of the total contribution to Trp-cage folding. The hydrophobic core centered on Trp6 contributed 15.6 kcal/mol to the total energy, whereas the solvation entropy contribution was 6.3 kcal/mol. The salt bridge formed by the hydrophilic side chains of Asp9 and Arg16 contributed 10.9 and 5.0 kcal/mol, respectively. This indicates that not only the hydrophobic core but also the salt bridge of the hydrophilic side chains gain solvation entropy and contribute to stabilizing the native structure of the Trp-cage.

Biological and thermodynamic stabilization of lipid-based delivery systems through natural biopolymers; controlled release and molecular dynamics simulations

Nowadays, the use of lipid-based nanocarriers for the targeted and controlled delivery of a variety of hydrophobic and hydrophilic bioactive-compounds and drugs has increased significantly. However, challenges such as thermodynamic instability, oxidation, and degradation of lipid membranes, as well as the unintended release of loaded compounds, have limited the use of these systems in the food and pharmaceutical industries. Therefore, the present study reviews the latest achievements in evaluating the characteristics, production methods, challenges, functional, and biological stabilization strategies of lipid-based carriers (including changes in formulation composition, structural modification, membrane-rigidity, and finally monolayer or multilayer coating with biopolymers) in different conditions, as well as molecular dynamics simulations. The scientists' findings indicate the effect of natural biopolymers (such as chitosan, calcium alginate, pectin, dextran, xanthan, caseins, gelatin, whey-proteins, zein, and etc.) in modifying the external structure of lipid-based carriers, improving thermodynamic stability and resistance of membranes to physicochemical and mechanical tensions. However, depending on the type of bioactive compound as well as the design and production goals of the delivery-system, selecting the appropriate biopolymer has a significant impact on the stability of vesicles and maintaining the bioaccessibility of the loaded-compounds due to the stresses caused by the storage-conditions, formulation, processing and gastrointestinal tract.

Computational Epitope Prediction and Design for Antibody Development and Detection

The design of optimized protein antigens is a fundamental step in the development of new vaccine candidates and in the detection of therapeutic antibodies. A fundamental prerequisite is the identification of antigenic regions that are most prone to interact with antibodies, namely, B-cell epitopes. Here, we describe an efficient structure-based computational method for epitope prediction, called MLCE. In this approach, all that is required is the 3D structure of the antigen of interest. MLCE can be applied to glycosylated proteins, facilitating the identification of immunoreactive versus immune-shielding carbohydrates.

Bicyclo-DNA mimics with enhanced protein binding affinities: insights from molecular dynamics simulations

DNA-protein interactions occur at all levels of DNA expression and replication and are crucial determinants for the survival of a cell. Several modified nucleotides have been utilized to manipulate these interactions and have implications in drug discovery. In the present article, we evaluated the binding of bicyclo-nucleotides (generated by forming a methylene bridge between C1' and C5' in sugar, leading to a bicyclo system with C2' axis of symmetry at the nucleotide level) to proteins. We utilized four ssDNA-protein complexes with experimentally known binding free energies and investigated the binding of modified nucleotides to proteins via all-atom explicit solvent molecular dynamics (MD) simulations (200 ns), and compared the binding with control ssDNA-protein systems. The modified ssDNA displayed enhanced binding to proteins as compared to the control ssDNA, as seen by means of MD simulations followed by MM-PBSA calculations. Further, the Delphi-based electrostatic estimation revealed that the high binding of modified ssDNA to protein might be related to the enhanced electrostatic complementarity displayed by the modified ssDNA molecules in all the four systems considered for the study. The improved binding achieved with modified nucleotides can be utilized to design and develop anticancer/antisense molecules capable of targeting proteins or ssRNAs.Communicated by Ramaswamy H. Sarma.

Molecular dynamics and protein frustration analysis of human fused in Sarcoma protein variants in Amyotrophic Lateral Sclerosis type 6: An In Silico approach

Amyotrophic lateral sclerosis (ALS) is the most frequent adult-onset motor neuron disorder. The disease is characterized by degeneration of upper and lower motor neurons, leading to death usually within five years after the onset of symptoms. While most cases are sporadic, 5%-10% of cases can be associated with familial inheritance, including ALS type 6, which is associated with mutations in the Fused in Sarcoma (FUS) gene. This work aimed to evaluate how the most frequent ALS-related mutations in FUS, R521C, R521H, and P525L affect the protein structure and function. We used prediction algorithms to analyze the effects of the non-synonymous single nucleotide polymorphisms and performed evolutionary conservation analysis, protein frustration analysis, and molecular dynamics simulations. Most of the prediction algorithms classified the three mutations as deleterious. All three mutations were predicted to reduce protein stability, especially the mutation R521C, which was also predicted to increase chaperone binding tendency. The protein frustration analysis showed an increase in frustration in the interactions involving the mutated residue 521C. Evolutionary conservation analysis showed that residues 521 and 525 of human FUS are highly conserved sites. The molecular dynamics results indicate that protein stability could be compromised in all three mutations. They also affected the exposed surface area and protein compactness. The analyzed mutations also displayed high flexibility in most residues in all variants, most notably in the interaction site with the nuclear import protein of FUS.

Mechanism of collagen folding propagation studied by Molecular Dynamics simulations

Collagen forms a characteristic triple helical structure and plays a central role for stabilizing the extra-cellular matrix. After a C-terminal nucleus formation folding proceeds to form long triple-helical fibers. The molecular details of triple helix folding process is of central importance for an understanding of several human diseases associated with misfolded or unstable collagen fibrils. However, the folding propagation is too rapid to be studied by experimental high resolution techniques. We employed multiple Molecular Dynamics simulations starting from unfolded peptides with an already formed nucleus to successfully follow the folding propagation in atomic detail. The triple helix folding was found to propagate involving first two chains forming a short transient template. Secondly, three residues of the third chain fold on this template with an overall mean propagation of ~75 ns per unit. The formation of loops with multiples of the repeating unit was found as a characteristic misfolding event especially when starting from an unstable nucleus. Central Gly→Ala or Gly→Thr substitutions resulted in reduced stability and folding rates due to structural deformations interfering with folding propagation.

The extracellular matrix is stabilized by collagen, a fibrillar protein structure, which represents the most abundant protein of the human body. Collagen consists of three peptide chains that form an elongated triple helix with a repeating and largely conserved sequence pattern of two proline (or hydroxyproline) residues followed by a glycine. Several human diseases are associated with mutations in collagen. The folding propagation is the most critical step in the collagen structure formation and not well understood. We have used multiple Molecular Dynamics simulations to specifically investigate the mechanism of triple helix propagation and how it is affected by mutations. The folding propagation was found to involve first two chains forming a short transient template followed by three residues of the third chain to fold on this template. Additional simulations were used to characterize misfolding events such as loop formation and the effect of glycine substitutions on collagen folding.