Application of principal component analysis in protein unfolding: an all-atom molecular dynamics simulation study.
29 Oct 2007-Journal of Chemical Physics (American Institute of Physics)-Vol. 127, Iss: 16, pp 165103-165103
TL;DR: The cross-correlation matrix obtained from MD simulation trajectory provided important information regarding the anisotropy of backbone dynamics that leads to unfolding and was applied to give a new insight to protein dynamics.
Abstract: We have performed molecular dynamics (MD) simulation of the thermal denaturation of one protein and one peptide—ubiquitin and melittin. To identify the correlation in dynamics among various secondary structural fragments and also the individual contribution of different residues towards thermal unfolding, principal component analysis method was applied in order to give a new insight to protein dynamics by analyzing the contribution of coefficients of principal components. The cross-correlation matrix obtained from MD simulation trajectory provided important information regarding the anisotropy of backbone dynamics that leads to unfolding. Unfolding of ubiquitin was found to be a three-state process, while that of melittin, though smaller and mostly helical, is more complicated.
Citations
More filters
[...]
TL;DR: The findings related to the folding of ubiquitin are consistent, for the most part, with the folding principles derived from the simulation of fast-folding proteins, suggesting that these principles may be applicable to a wider range of proteins.
Abstract: Equilibrium molecular dynamics simulations, in which proteins spontaneously and repeatedly fold and unfold, have recently been used to help elucidate the mechanistic principles that underlie the folding of fast-folding proteins. The extent to which the conclusions drawn from the analysis of such proteins, which fold on the microsecond timescale, apply to the millisecond or slower folding of naturally occurring proteins is, however, unclear. As a first attempt to address this outstanding issue, we examine here the folding of ubiquitin, a 76-residue-long protein found in all eukaryotes that is known experimentally to fold on a millisecond timescale. Ubiquitin folding has been the subject of many experimental studies, but its slow folding rate has made it difficult to observe and characterize the folding process through all-atom molecular dynamics simulations. Here we determine the mechanism, thermodynamics, and kinetics of ubiquitin folding through equilibrium atomistic simulations. The picture emerging from the simulations is in agreement with a view of ubiquitin folding suggested from previous experiments. Our findings related to the folding of ubiquitin are also consistent, for the most part, with the folding principles derived from the simulation of fast-folding proteins, suggesting that these principles may be applicable to a wider range of proteins.
287 citations
[...]
TL;DR: The present data suggest that a family-centered point of view is necessary in the comparative analyses of cold- and warm-adapted enzymes, and enzymes belonging to the same family or superfamily have evolved similar structural and dynamic patterns to overcome the detrimental effects of low temperatures.
Abstract: The identification of molecular mechanisms underlying enzyme cold adaptation is a hot-topic both for fundamental research and industrial applications. In the present contribution, we review the last decades of structural computational investigations on cold-adapted enzymes in comparison to their warm-adapted counterparts. Comparative sequence and structural studies allow the definition of a multitude of adaptation strategies. Different enzymes carried out diverse mechanisms to adapt to low temperatures, so that a general theory for enzyme cold adaptation cannot be formulated. However, some common features can be traced in dynamic and flexibility properties of these enzymes, as well as in their intra- and inter-molecular interaction networks. Interestingly, the current data suggest that a family-centered point of view is necessary in the comparative analyses of cold- and warm-adapted enzymes. In fact, enzymes belonging to the same family or superfamily, thus sharing at least the three-dimensional fold and common features of the functional sites, have evolved similar structural and dynamic patterns to overcome the detrimental effects of low temperatures.
40 citations
Cites background or methods from "Application of principal component ..."
[...]
[...]
[...]
TL;DR: Conformers with similar dynamic and structural features, as revealed by LLS relaxation times, could be observed, in the absence of urea, in two ubiquitin mutants, L67S and L69S.
Abstract: The relaxation of long-lived states (LLS) corresponds to the slow return to statistical thermal equilibrium between symmetric and antisymmetric proton spin states. This process is remarkably sensitive to the presence of external spins and can be used to obtain information about partial unfolding of proteins. We detected the appearance of a destabilized conformer of ubiquitin when urea is added to the protein in its native state. This conformer shows increased mobility in the C-terminus, which significantly extends the lifetimes of proton LLS magnetisation in Ser-65. These changes could not be detected by conventional measurements of T(1) and T(2) relaxation times of protons, and would hardly be sensed by carbon-13 or nitrogen-15 relaxation measurements. Conformers with similar dynamic and structural features, as revealed by LLS relaxation times, could be observed, in the absence of urea, in two ubiquitin mutants, L67S and L69S.
39 citations
[...]
TL;DR: Trajectories obtained from molecular dynamics simulations of some large-ring cyclodextrins (LR-CDs) were analyzed by a statistical method, principal component analysis (PCA), thus enabling the monitoring of the dominant PCA modes for concerted motions of the macroring atoms in a lower-dimensions subspace.
Abstract: Trajectories obtained from molecular dynamics (MD) simulations of some large-ring cyclodextrins (LR-CDs) were analyzed by a statistical method, principal component analysis (PCA), thus enabling the monitoring of the dominant PCA modes for concerted motions of the macroring atoms in a lower-dimensions subspace. Earlier analyses of macrorings' conformational deformations based on examination of snapshots extracted from the MD simulation trajectories were further supported on more quantitative grounds. The first 10 lowest-indexed modes describe more than 90% of the total atomic motion in all cases, with more than 50% of the contribution coming from the two highest-eigenvalue principal components. Representative average geometries of the cyclodextrin macroring were also obtained.
33 citations
[...]
TL;DR: This study shows that molecular dynamics simulations combined with calculation of crowding effects provide an avenue for characterize the transition-state ensemble in atomic details.
Abstract: The four-helix bundle protein Rd-apocyt b562, a redesigned stable variant of apocytochrome b562, exhibits two-state folding kinetics. Its transition-state ensemble has been characterized by Φ-value analysis. To elucidate the molecular basis of the transition-state ensemble, we have carried out high-temperature molecular dynamics simulations of the unfolding process. In six parallel simulations, unfolding started with the melting of helix I and the C-terminal half of helix IV, and followed by helix III, the N-terminal half of helix IV and helix II. This ordered melting of the helices is consistent with the conclusion from native-state hydrogen exchange, and can be rationalized by differences in intrinsic helix propensity. Guided by experimental Φ-values, a putative transition-state ensemble was extracted from the simulations. The residue helical probabilities of this transition-state ensemble show good correlation with the Φ-values. To further validate the putative transition-state ensemble, the effect of macromolecular crowding on the relative stability between the unfolded ensemble and the transition-state ensemble was calculated. The resulting effect of crowding on the folding kinetics agrees well with experimental observations. This study shows that molecular dynamics simulations combined with calculation of crowding effects provide an avenue for characterize the transition-state ensemble in atomic details.
19 citations
References
More filters
[...]
TL;DR: The crystal structure of human erythrocytic ubiquitin has been refined at 1.8 A resolution using a restrained least-squares procedure and features a number of unusual secondary structural features, including a parallel G1 beta-bulge, two reverse Asx turns, and a symmetrical hydrogen-bonding region that involves the two helices and two of the reverse turns.
Abstract: The crystal structure of human erythrocytic ubiquitin has been refined at 1.8 A resolution using a restrained least-squares procedure. The crystallographic R-factor for the final model is 0.176. Bond lengths and bond angles in the molecule have root-mean-square deviations from ideal values of 0.016 A and 1.5 degrees, respectively. A total of 58 water molecules per molecule of ubiquitin are included in the final model. The last four residues in the molecule appear to have partial occupancy or large thermal motion. The overall structure of ubiquitin is extremely compact and tightly hydrogen-bonded; approximately 87% of the polypeptide chain is involved in hydrogen-bonded secondary structure. Prominent secondary structural features include three and one-half turns of alpha-helix, a short piece of 3(10)-helix, a mixed beta-sheet that contains five strands, and seven reverse turns. There is a marked hydrophobic core formed between the beta-sheet and alpha-helix. The molecule features a number of unusual secondary structural features, including a parallel G1 beta-bulge, two reverse Asx turns, and a symmetrical hydrogen-bonding region that involves the two helices and two of the reverse turns.
1,551 citations
[...]
TL;DR: An implementation of classical molecular dynamics on parallel computers of increased efficiency has enabled a simulation of protein folding with explicit representation of water for 1 microsecond, about two orders of magnitude longer than the longest simulation of a protein in water reported to date.
Abstract: An implementation of classical molecular dynamics on parallel computers of increased efficiency has enabled a simulation of protein folding with explicit representation of water for 1 microsecond, about two orders of magnitude longer than the longest simulation of a protein in water reported to date. Starting with an unfolded state of villin headpiece subdomain, hydrophobic collapse and helix formation occur in an initial phase, followed by conformational readjustments. A marginally stable state, which has a lifetime of about 150 nanoseconds, a favorable solvation free energy, and shows significant resemblance to the native structure, is observed; two pathways to this state have been found.
1,321 citations
[...]
TL;DR: Melittin from bee venom is water-soluble, yet integrates into membranes and lyses cells, and is referred to as a surface-active protein in the two crystal forms grown from concentrated salt solutions.
Abstract: Melittin from bee venom is water-soluble, yet integrates into membranes and lyses cells. Each melittin chain consists of 26 amino acid residues and in aqueous salt solutions it exists as a tetramer. We have determined the molecular structure of the tetramer in two crystal forms grown from concentrated salt solutions. In both crystal forms the melittin polypeptide is a bent alpha-helical rod, with the "inner" surface largely consisting of hydrophobic sidechains and the "outer" surface consisting of hydrophilic side chains. Thus, the helix is strongly amphiphilic. In the tetramer, four such helices contribute their hydrophobic side chains to the center of the molecule. The packing of melittin tetramers is also very similar in the two crystal forms: they are packed in planar layers with the outsides forming hydrophilic surfaces and the insides (the centers of melittin tetramers) forming a hydrophobic surface. We suggest that the surface activity of melittin can be rationalized in terms of these surfaces. The lytic activity of melittin can also be interpreted in terms of the molecular structure observed in the crystals: the hydrophobic inner surface of a melittin helix may integrate into the apolar region of a bilayer with the helix axis approximately parallel to the plane of the bilayer, and with the hydrophilic surface exposed to the aqueous phase. This integration would be expected to disrupt the bilayer because of melittin helix would penetrate only a short distance into it. Additionally, the integration of melittin from one side of a bilayer would produce a surface area difference across the bilayer, perhaps leading to lysis. In this view, melittin is distinct from membrane proteins that penetrate evenly into both leaflets of a bilayer or exactly halfway through a bilayer, and hence we refer to melittin as a surface-active protein.
365 citations