1 / 23

Behaviour of velocities in protein folding events

Behaviour of velocities in protein folding events. Aldo Rampioni, University of Groningen Leipzig, 17th May 2007. Plan of the talk. Questions that we want to address System studied: the ß-heptapeptide Definition of folding event Methodology used for the analysis Results of the analysis

pegeen
Download Presentation

Behaviour of velocities in protein folding events

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Behaviour of velocitiesinprotein folding events Aldo Rampioni, University of Groningen Leipzig, 17th May 2007

  2. Plan of the talk • Questions that we want to address • System studied: the ß-heptapeptide • Definition of folding event • Methodology used for the analysis • Results of the analysis • Final remarks

  3. Questions that we want to address: Do velocities show any correlation or cooperative behaviour during the protein folding event? Can this information be used to detect when the folding event occurs? Imagine being an amino-acid…

  4. ß-heptapeptide • Small peptide  fast simulations • 50 ns sufficient to generate an equilibrium distribution (multiple folding-unfolding events good statistics) Figures from Daura, X et al. PROTEINS: Struc. Func. Gen. 34 (1999)

  5. GROMACS 3.2.1 software package force field GROMOS96 43a1 The end groups were protonated -NH3+ and –COOH Solvent: methanol (962 molecules) [model B3 in J.Chem.Phys.112 (2000)] Temperature 340 K Time step of 2fs Twin-range cutoff of 0.8/1.4 nm for all non-bonded interactions Initial structure: helical fold (shown in figure) Ten 50-ns MD simulations were performed using: Simulation conditions Five 100-ns MD simulations: Same conditions as above, but starting from an unfolded conformation

  6. Definition of folding event(first trial) We used a criterion of similarity (RMSD) to group different structures (cluster algorithm) and build a dynamics on grapho. It is natural to define “folding event” each jump to the cluster representative of the folded structure.

  7. 2 For each pair of structures the RMSD was calculated after fitting the backbone atoms of residues 2 to 6. 1 Structures were extracted from the trajectories at regular time intervals for analysis 6 This process was iterated until all structures were assigned to a cluster. 3 Using the criterion of similarity of two structures RMSD<cutoff, the number of neighbours for each of the structures in the initial pool was determined. 5 All the structures belonging to this cluster were removed from the pool. 4 The structure with the highest number of neighbours was the centre of the first cluster. Cluster algorithm

  8. Choice of the cutoff

  9. Cluster analysis over 50 ns

  10. 3-2 1-1 2-3 4-4 5-5 Central structures of the five most populated clusters Blue time interval 10 ps Red time interval 50 ps

  11. Time series of cluster

  12. Transitions among the 5 most populated clusters over 50 ns The total number of transitions among all clusters is1224/322 * After switching 2 and 3 in the cluster numbering of the set got using 50 ps time interval

  13. Limits of this definition • The representative structure of cluster number 2 and 5 are very close to the folded structure, i.e. the jump from those clusters to the cluster number 1 is the last step of different folding paths • How to consider jumps to cluster number 1 followed by an immediate jump out?

  14. Definition of folding event(second trial) We simply used a criterion of similarity (RMSD) to the folded structure, introducing two thresholds: below the lower one we consider the peptide folded, above the higher we consider the peptide unfolded. We define “folding event” every time the RMSD pass from values higher than the upper threshold to values lower than the bottom threshold.

  15. Definiton of folding event F: 2<n<7 S: 7<n VF: n<3

  16. According to this definition we extracted from 1 s simulation: 49 VERYFAST folding events 42 FAST folding events 40 SLOW folding events These events have been aligned choosing as t0 the last time the RMSD is above the higher threshold

  17. Methods j = 1,…,N denotes the atom coordinate k = 1,…,T denotes the time i = 1,…,M denotes the trajectory is the ith trajectory is a slice of the matrix at time k the average is over the trajectories Covariance matrix at time k Time autocovariance

  18. Covariance matrices of the velocities of the backbone atoms between t0-500 and t0+500 ps

  19. RMSD RGYR PC1 PC2 CV1 CV2 RGYR PC1 PC2 CV2 RMSD CV1

  20. If the principal components of motions in cartesian space do not correlate with the order parameter (RMSD), there is no hope to see something looking at velocities in cartesian space Thus we chose to investigate some internal degree of freedom such as torsional angles

  21. Acknowledgments 28th of April, Zlotoryja, Poland Dr. Tsjerk Wassenaar, University of Utrecht, The Netherlands Prof. Alan Mark, University of Queensland, Australia A particular thank to Drs. Magdalena Siwko now…in Rampioni!!!

More Related