Protein Folding and Confinement: Inherent Structure Analysis of Chaperonin Action

Effects of Residual Composition and Distribution on the Structural Characteristics of the Protein

The effect of ratio and consecutive number of hydrophobic residues in the repeating unit of protein chains was investigated by MD simulation. The modified off-lattice HNP model was applied in this study. The protein chains constituted by different HNP ratios or different numbers of consecutively hydrophobic residues with the same chain length were simulated under a broad temperature range. We concluded that the proteins with higher ratio or larger number of sequentially hydrophobic residues present more orientated and compact structure under a certain low temperature. It is attributed to the lower non-bonded potential energy between H-H residual pairs, especially more hydrophobic residues in a procession among the protein chain. Considering the microscopic structure of the protein, more residue contacts are achieved with the proteins with higher ratios and sequential H residues under the low temperature. Meanwhile, with the ratio and consecutive number of H residues increasing, the distribution of stem length showed a transition from exponential decline to unimodal and even multiple peaks, indicating the specific ordered structure formed. These results provide an insight into 3D structural properties of proteins from their residue sequences, which has a primary structure at molecular level and, ultimately, a practical possibility of applying in biotechnological applications.

Inherent structure versus geometric metric for state space discretization

Inherent structure (IS) and geometry-based clustering methods are commonly used for analyzing molecular dynamics trajectories. ISs are obtained by minimizing the sampled conformations into local minima on potential/effective energy surface. The conformations that are minimized into the same energy basin belong to one cluster. We investigate the influence of the applications of these two methods of trajectory decomposition on our understanding of the thermodynamics and kinetics of alanine tetrapeptide. We find that at the microcluster level, the IS approach and root-mean-square deviation (RMSD)-based clustering method give totally different results. Depending on the local features of energy landscape, the conformations with close RMSDs can be minimized into different minima, while the conformations with large RMSDs could be minimized into the same basin. However, the relaxation timescales calculated based on the transition matrices built from the microclusters are similar. The discrepancy at the microcluster level leads to different macroclusters. Although the dynamic models established through both clustering methods are validated approximately Markovian, the IS approach seems to give a meaningful state space discretization at the macrocluster level in terms of conformational features and kinetics.

Conformational sampling of peptides in the presence of protein crowders from AA/CG-multiscale simulations

Macromolecular crowding is recognized as an important factor influencing folding and conformational dynamics of proteins and nucleic acids. Previous views of crowding have focused on the mostly entropic volume exclusion effect of crowding, but recent studies are indicating the importance of enthalpic effects, in particular, changes in electrostatic interactions due to crowding. Here, temperature replica exchange molecular dynamics simulations of trp-cage and melittin in the presence of explicit protein crowders are presented to further examine the effect of protein crowders on peptide dynamics. The simulations involve a three-component multiscale modeling scheme where the peptides are represented at an atomistic level, the crowder proteins at a coarse-grained level, and the surrounding aqueous solvent as implicit solvent. This scheme optimally balances a physically realistic description for the peptide with computational efficiency. The multiscale simulations were compared with simulations of the same peptides in different dielectric environments with dielectric constants ranging from 5 to 80. It is found that the sampling in the presence of the crowders resembles sampling with reduced dielectric constants between 10 and 40. Furthermore, diverse conformational ensembles are generated in the presence of crowders including partially unfolded states for trp-cage. These findings emphasize the importance of enthalpic interactions over volume exclusion effects in describing the effects of cellular crowding.