The Trp-cage miniprotein with single-site mutations: Studies of stability and dynamics using molecular dynamics

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.

Boosting the conformational sampling by combining replica exchange with solute tempering and well-sliced metadynamics

Methods that combine collective variable (CV) based enhanced sampling and global tempering approaches are used in speeding-up the conformational sampling and free energy calculation of large and soft systems with a plethora of energy minima. In this paper, a new method of this kind is proposed in which the well-sliced metadynamics approach (WSMTD) is united with replica exchange with solute tempering (REST2) method. WSMTD employs a divide-and-conquer strategy wherein high-dimensional slices of a free energy surface are independently sampled and combined. The method enables one to accomplish a controlled exploration of the CV-space with a restraining bias as in umbrella sampling, and enhance-sampling of one or more orthogonal CVs using a metadynamics like bias. The new hybrid method proposed here enables boosting the sampling of more slow degrees of freedom in WSMTD simulations, without the need to specify associated CVs, through a replica exchange scheme within the framework of REST2. The high-dimensional slices of the probability distributions of CVs computed from the united WSMTD and REST2 simulations are subsequently combined using the weighted histogram analysis method to obtain the free energy surface. We show that the new method proposed here is accurate, improves the conformational sampling, and achieves quick convergence in free energy estimates. We demonstrate this by computing the conformational free energy landscapes of solvated alanine tripeptide and Trp-cage mini protein in explicit water.

Protein aggregation due to nsSNP resulting in P56S VABP protein is associated with amyotrophic lateral sclerosis

Mutations in the gene encoding vesicle-associated membrane protein (VAPB) cause amyotrophic lateral sclerosis (ALS), a fatal neurodegenerative disorder. The VAPB gene is mapped to chromosome number 20 and can be found at cytogenetic location 20q13.33 of the chromosome. VAPB is seen to play a significant role in the unfolded protein response (UPR), which is a process that suppresses the accumulation of unfolded proteins in the endoplasmic reticulum. Earlier studies have reported two points; which we have analyzed in our study. Firstly, the mutation P56S in the VAPB is seen to increase the stability of the protein and secondly, the mutation P56S in VAPB is seen to interrupt the functioning of the gene and loses its ability to be involved in the activation of the IRE1/XBP1 pathway which leads to ALS. With correlation on the previous research studies on the stability of this protein, we carried out Molecular dynamics (MD) simulation. We analyzed the SNP results of 17 nsSNPs obtained from dbSNP using SIFT, polyphen, I-Mutant, SNP&GO, PhDSNP and Mutpred to predict the role of nsSNPs in VAPB. MD simulation is carried out and plots for RMSD, RMSF, Rg, SASA, H-bond and PCA are obtained to check and prove the stability of the wild type and the mutant protein structure. The protein is checked for its aggregation and the results obtained show changes in the protein structure that might result in the loss of function.

Folding dynamics of the Trp-cage miniprotein: evidence for a native-like intermediate from combined time-resolved vibrational spectroscopy and molecular dynamics simulations

Trp-cage is a synthetic 20-residue miniprotein which folds rapidly and spontaneously to a well-defined globular structure more typical of larger proteins. Due to its small size and fast folding, it is an ideal model system for experimental and theoretical investigations of protein folding mechanisms. However, Trp-cage's exact folding mechanism is still a matter of debate. Here we investigate Trp-cage's relaxation dynamics in the amide I' spectral region (1530-1700 cm(-1)) using time-resolved infrared spectroscopy. Residue-specific information was obtained by incorporating an isotopic label ((13)C═(18)O) into the amide carbonyl group of residue Gly11, thereby spectrally isolating an individual 310-helical residue. The folding-unfolding equilibrium is perturbed using a nanosecond temperature-jump (T-jump), and the subsequent re-equilibration is probed by observing the time-dependent vibrational response in the amide I' region. We observe bimodal relaxation kinetics with time constants of 100 ± 10 and 770 ± 40 ns at 322 K, suggesting that the folding involves an intermediate state, the character of which can be determined from the time- and frequency-resolved data. We find that the relaxation dynamics close to the melting temperature involve fast fluctuations in the polyproline II region, whereas the slower process can be attributed to conformational rearrangements due to the global (un)folding transition of the protein. Combined analysis of our T-jump data and molecular dynamics simulations indicates that the formation of a well-defined α-helix precedes the rapid formation of the hydrophobic cage structure, implying a native-like folding intermediate, that mainly differs from the folded conformation in the orientation of the C-terminal polyproline II helix relative to the N-terminal part of the backbone. We find that the main free-energy barrier is positioned between the folding intermediate and the unfolded state ensemble, and that it involves the formation of the α-helix, the 310-helix, and the Asp9-Arg16 salt bridge. Our results suggest that at low temperature (T ≪ Tm) a folding path via formation of α-helical contacts followed by hydrophobic clustering becomes more important.