- Open Access
How far in-silico computing meets real experiments. A study on the structure and dynamics of spin labeled vinculin tail protein by molecular dynamics simulations and EPR spectroscopy
© Prasad Gajula et al.; licensee BioMed Central Ltd. 2013
- Published: 15 February 2013
Investigation of conformational changes in a protein is a prerequisite to understand its biological function. To explore these conformational changes in proteins we developed a strategy with the combination of molecular dynamics (MD) simulations and electron paramagnetic resonance (EPR) spectroscopy. The major goal of this work is to investigate how far computer simulations can meet the experiments.
Vinculin tail protein is chosen as a model system as conformational changes within the vinculin protein are believed to be important for its biological function at the sites of cell adhesion. MD simulations were performed on vinculin tail protein both in water and in vacuo environments. EPR experimental data is compared with those of the simulated data for corresponding spin label positions.
The calculated EPR spectra from MD simulations trajectories of selected spin labelled positions are comparable to experimental EPR spectra. The results show that the information contained in the spin label mobility provides a powerful means of mapping protein folds and their conformational changes.
The results suggest the localization of dynamic and flexible regions of the vinculin tail protein. This study shows MD simulations can be used as a complementary tool to interpret experimental EPR data.
- Electron Paramagnetic Resonance
- Molecular Dynamic Simulation
- Electron Paramagnetic Resonance Spectrum
- Spin Label
Rapid advances in computer technology have led to the development of successful molecular simulations of protein structural dynamics that are intrinsic to understand biological processes. These simulations have resulted in the development of novel models and methods that increasingly agree with experimental observations, and suggest new experiments providing insights into biological mechanisms. Used in combination with the information gained by sophisticated experimental techniques, molecular simulations can help us, to understand biological complexity at the atomic and molecular levels. Here, we emphasize such an approach that illustrates the potential of molecular dynamics simulations in analyzing experimental results determined by EPR spectroscopy on vinculin as an example.
Site directed spin labeling EPR spectroscopy has evolved as a powerful technique to investigate protein structure and conformational changes under physiological or near-physiological conditions [15–18]. The shape of continuous wave (cw) EPR spectra recorded at room temperature is sensitive to the re-orientational motion of the bound spin label side chain providing information on the motional restriction of the nitroxide due to sterical interaction with the secondary and tertiary structure [18–20]. In addition, solvent accessibility of the spin label side chain and polarity of the nitroxide microenvironment characterize the protein topology. Dipolar coupling between two spin labels incorporated into a protein report on intra-molecular distances [21, 23]. Details of this method are summarized in recent reviews by Bordgnon et al, and Klare et al. Here we report on MD simulations of spin labelled vinculin tail in order to analyze experimental spin label EPR spectra.
The initial coordinates for the vinculin tail MD simulations were obtained from the Protein Data Bank (PDB code: 1ST6). All MD simulations were performed with the GROMACS simulation suite for in water and in vacuo simulations. The force field ffG43a1 was used for the simulations in water and the force field ffG43b1 from GROMOS, which is integrated into GROMACS, was used for the simulations in vacuo. Several MD simulations were performed with different spin labeled sites for 30 ns at 300 K in water, whereas the simulations in vacuo were performed for 10 ns at 600 K. In order to assure that the spin label covers the maximum accessible conformational space within a relative short MD run of 10 ns, a high temperature of 600 K was applied with position restraints on all backbone atoms. Periodic boundary conditions were applied as mentioned elsewhere. For a detailed description of the MD simulations methodology and applications refer [25, 29, 30]. The vinculin tail protein holds two native cysteines at positions 950 and 972 on helix 3. Cysteines were genetically introduced at positions 901, 909 on helix 1 (H1); 922,927,934 on helix 2(H2); 957 on helix 3(H3); 984 on helix 4(H4); 1024, 1033 on helix 5(H5); and 1062 in the C-terminus. EPR experiments on these mutants were described in . Refer to [24–26] for calculating EPR spectra from MD simulations data.
Spin label dynamics in water and in vacuo
Comparison of simulated and experimental data
Inter helical distances
26.5 ± 1
8 ± 1
8 ± 1
7 ± 1
The reorientational angle distribution for the spin label at position 901 (on H1 helix) in water and in vacuo simulations is presented in Figure 3. Similar broad distributions are found for this position both in water and in vacuo simulations. This spin label is located near the flexible loop region and also the head group of the side chain is exposed to the surface of the protein. The spin label at position 909(H1) shows a narrow distribution of β in water simulations as well as in vacuo simulations. This indicates that the spin label mobility at this position is restricted. Inspection of the structure indicates that this spin label is buried between two helices. The reorientational dynamics of the spin label at position 922(H2) is characterized by a broad distribution both in water and in vacuo simulations. This indicates that the spin label at this position is mobile. Though part of the spin label linker is buried, the nitroxide shows a significant mobility because the more flexible dihedrals χ4 and χ5 are not restricted (Figure 2). The spin label at position 927(H2) shows a moderate distribution of β in water simulations when compared to the position 922(H2). However the distribution of β in vacuo simulations indicates that the spin label at this position is highly restricted due to contacts with neighbouring helix atoms. The spin label at position 934(H2) shows restricted mobility, both in water and in vacuo simulations with narrow β distributions. This spin label side chain is surrounded by bulky residues. Considering the overall mobility pattern, the H2 helix shows restricted mobility with low β distributions.
Helix H3 carries three spin labels located at positions 950, 957 and 972. The distribution plots indicate that the mobility is moderate in water simulations, whereas in vacuo simulations the mobility is significantly increased. A very interesting distribution pattern is found in vacuo simulations of 972R1 that indicates two rotameric states of the spin label that are equally distributed (data not shown). This spin label is located near the loop region and oriented between two helices. However a functional study has shown a perturbed conformation of protein due to the bound spin label at position 972. Helix H4 has only one spin label at position 984 showing very narrow distribution of β in water simulations and a larger distribution in vacuo simulations. This spin label location is near to the kink region of the helix. The narrow distribution of β in water simulations reveals that this spin label has strong interactions during the simulations, whereas in the vacuo simulations the spin label is mobile. This indicates helical movement in the in water simulations that brings the spin label into strong tertiary contacts. H5 helix holds 1024R1 and 1033R1. Both MD simulations in water and in vacuo indicate that 1024R1 is very mobile when compared to 1033R1. 1024R1 is a surface exposed site. Though the spin label at 1033 is partially surface exposed, its mobility is restricted by the presence of an arginine residue located at 987 in the adjacent helix. 1062R1 is located on a C-terminus end that naturally shows a broad distribution of β in both in water and in vacuo simulations.
From Figures 5 and 6, it is clear that high mobility is associated with residues 901, 922, 950, 972, 1024 and 1062. Among these spin labels, the one at position 1062 is in the C-terminal region. 901R1 that belongs to H1 helix is immediately next to the loop region. Therefore it shows a high degree of freedom in water simulations.
The structural dynamics within a subunit of the vinculin tail protein depicted by the relative motion between spin labels 901 and 957 in helices H1 and H3 respectively is shown in Figure 7. The distance distribution shows the mean distance of 27.2 Å with an overall width of about 2 Å. The distance values between spin labels at positions 922(H2) and 957(H3) reveal a smaller distribution. A study by Palmer et al  suggests that conformational change in the vinculin C-terminal may depend on a critical histidine residue at position 906 in H1. The conformational change triggered by the presence of this histidine may push/pull H1 with respect to H3. The distribution of distance between 901 and 1033 is similar to the distribution between 901 and 957. It is interesting that the distance distribution between 984R1 (H4) and 1033R1 (H5) shows two distinct maxima at 7Å and 9Å respectively. It indicates that one of the spin labels fluctuates around two conformations out of which the most dominant conformation is at 7Å. This is in agreement with the observation that the angle β of 1033 shows two distinct maxima. The existence of multiple conformations in dynamic equilibrium also raises questions regarding the internal mobility of the individual side chains as well as of larger structural domains. However, the MD simulations data show that both these spin labeled positions show significant motional restriction. The simulated inter-nitroxide mean distances are in-line with the results of EPR experiments  (cf. Table 1).
The MD simulations data together with EPR results indicate, while the majority of the spin labeled sites in the vinculin tail show considerable dynamics, perhaps the most significant ones are those belonging to helices H1, H3, H5 and naturally the C-terminal end residue 1062. (refer to [31–33]). This dynamics might be of importance for the understanding of the unfolding process of the vinculin tail bundle and its interaction with the membrane.
MD simulations were performed on the spin labeled vinculin tail domain both in vacuo and in water environment. The behavior of the spin label at various positions of the vinculin tail domain was analyzed by means of RMSF analysis and simulated EPR spectra. The results were comparable to EPR experiments. A correlation in the dynamics of the spin label mobility was found when mobility parameter values from MD simulations were compared with the EPR inverse line width data for the most of the spin-labeled sites of vinculin tail. To estimate the magnitude of helix displacements in the vinculin tail domain, distance distributions between pairs of spin label side chains were calculated and compared with experimental data. The MD simulations results in combination with EPR data show that the information contained in the spin label mobility provide a powerful means of mapping protein folds and their changes.
Department of Science and Technology is highly appreciated for the funding through DST Ramanujan Fellowship award No. SR/S2/RJN-22/2011. We are very much thankful to Wolfgang Ziegler and Christoph Abe for all the valuable suggestions and extended help. We sincerely acknowledge everyone who directly or indirectly helped us to carryout this work.
The publication costs for this article were funded by the corresponding author's sponsor institution, DST through Ramanujan fellowship award No. SR/S2/RJN-22/2011.
This article has been published as part of BMC Genomics Volume 14 Supplement 2, 2013: Selected articles from ISCB-Asia 2012. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/14/S2.
- Ziegler WH, Liddington RC, Critchley DR: The structure and regulation of vinculin. Trends Cell Biol. 2006, 16: 453-460. 10.1016/j.tcb.2006.07.004.View ArticlePubMedGoogle Scholar
- DeMali KA: Vinculin - a dynamic regulator of cell adhesion. TRENDS in Biochem Sci. 2004, 29: 11-10.1016/j.tibs.2003.11.004.View ArticleGoogle Scholar
- Ziegler WH: The cytoskeletal connection: understanding adaptor proteins. Cell Migration: Signalling and Mechanisms. Transl Res Biomed. Edited by: Entschladen F, Zänker KS. 2010, Basel, Karger, 2: 136-162.View ArticleGoogle Scholar
- Jockusch BM, Riidiger M: Crosstalk between cell adhesion molecules: vinculin as a paradigm for regulation by conformation. Trends Cell Biol. 1996, 6: 311-315. 10.1016/0962-8924(96)10022-2.View ArticlePubMedGoogle Scholar
- Izard T et al: Vinculin activation by talin through helical bundle conversion. Nature. 2004, 427: 171-175. 10.1038/nature02281.View ArticlePubMedGoogle Scholar
- Bakolitsa C et al: Structural basis for vinculin activation at sites of cell adhesion. Nature. 2004, 430: 583-6. 10.1038/nature02610.View ArticlePubMedGoogle Scholar
- Bakolitsa C, de Pereda JM, Bagshaw CR, Critchley DR, Liddington RC: Crystal structure of the vinculin tail suggests a pathway for activation. Cell. 1999, 99: 603-613. 10.1016/S0092-8674(00)81549-4.View ArticlePubMedGoogle Scholar
- Johnson RP, Niggli V, Durrer P, Craig SW: A conserved motif in the tail domain of vinculin mediates association with and insertion into acidic phospholipid bilayers. Biochem. 1998, 37: 10211-222. 10.1021/bi9727242.View ArticleGoogle Scholar
- Niggli V: Structural properties of lipid-binding sites in cytoskeletal proteins. Trends Biochem Science. 2001, 26: 604-11. 10.1016/S0968-0004(01)01927-2.View ArticleGoogle Scholar
- DeMali KA, Barlow CA, Burridge K: Recruitment of the Arp2/3 complex to vinculin: coupling membrane protrusion to matrix adhesion. J Cell Biol. 2002, 159: 881-91. 10.1083/jcb.200206043.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller GJ, Dunn SD, Ball EH: Interaction of the N- and C-terminal domains of vinculin: characterization and mapping studies. J Biol Chem. 2001, 276: 11729-11734. 10.1074/jbc.M008646200.View ArticlePubMedGoogle Scholar
- McGregor A, Blanchard AD, Rowe AJ, Critchley DR: Identification of the vinculin-binding site in the cytoskeletal protein a-actinin. Biochem J. 1994, 301: 225-233.PubMed CentralView ArticlePubMedGoogle Scholar
- Johnson RP, Craig SW: An intramolecular association between the head and tail domains of vinculin modulates talin binding. J Biol Chem. 1994, 269: 12611-12619.PubMedGoogle Scholar
- Johnson RP, Craig SW: F-actin binding site masked by the intramolecular association of vinculin head and tail domains. Nature. 1995, 373: 261-264. 10.1038/373261a0.View ArticlePubMedGoogle Scholar
- Altenbach C, Marti T, Khorana HG, Hubbell WL: Transmembrane Protein Structure: Spin Labeling of Bacteriorhodopsin Mutants. Science. 1990, 248: 1088-1092. 10.1126/science.2160734.View ArticlePubMedGoogle Scholar
- Altenbach C, Steinhoff H-J, Greenhalgh DA, Khorana HG, Hubbell WL: Factors That Determine the EPR-Spectra of Nitroxide Side-Chains in Spin-Labeled Proteins and Analysis by Molecular-Dynamics Simulation. Biophys. 1994, 66:Google Scholar
- Altenbach C, Yang K, Farrens DL, Farahbakhsh ZT, Khorana HG, Hubbell WL: Structural Features and Light-Dependent Changes in the Cytoplasmic Interhelical E-F Loop Region of Rhodopsin: A Site-Directed Spin-Labeling Study. Biochemistry. 1996, 35: 12470-12478. 10.1021/bi960849l.View ArticlePubMedGoogle Scholar
- Bordignon E, Steinhoff H-J: Membrane protein structure and dynamics studied by site-directed spin labeling ESR. Biological Magnetic Resonance 27 - ESR Spectroscopy in Membrane Biophysics. Edited by: Hemminga MA, Berliner LJ. 2007, 129-164.View ArticleGoogle Scholar
- Pfeiffer M, Rink T, Gerwert K, Steinhoff H-J: Site-directed Spin-labeling Reveals the orientation of the Amino Acid Side-chains in the E-F Loop of Bacteriorhodopsin. J Mol Biol. 1999, 287: 163-171. 10.1006/jmbi.1998.2593.View ArticlePubMedGoogle Scholar
- Abe Christoph, Franziska D, Prasad Gajula , Monique B, Vogel KP, Maurice Gl, Susanne I, Wolfgang HG, Steinhoff H-J: Monomeric and Dimeric Conformation of the Vinculin Tail Five-Helix Bundle in Solution Determined by EPR-Spectroscopy. Biophysical Journal. 2011, 101 (7): 1772-1780. 10.1016/j.bpj.2011.08.048.PubMed CentralView ArticlePubMedGoogle Scholar
- Gajula P, Borovykh IV, Beier C, Shkuropatova T, Gast P, Steinhoff H-J: Spin-labeled photosynthetic reaction centers from Rhodobacter sphaeroides studied by electron paramagnetic resonance spectroscopy and molecular dynamics simulations. Applied Magnetic Resonance. 2007, 31: 167-178. 10.1007/BF03166254.View ArticleGoogle Scholar
- Borovykh V, Gajula P, Huber M, Gast P, Stenhoff H-J: Distance between a native cofactor and a spin label in the reaction centre of Rhodobacter sphaeroides by a two-frequency pulsed electron paramagnetic resonance method and molecular dynamics simulations. Journal of magnetic resonance. 2006, 180: 178-185. 10.1016/j.jmr.2006.02.008.View ArticlePubMedGoogle Scholar
- Klare JP, Steinhoff H-J: Spin Labeling EPR. Photosynth Res. 2009, 102: 377-90. 10.1007/s11120-009-9490-7.View ArticlePubMedGoogle Scholar
- Steinhoff HJ, Hubbell WL: Calculation of electron paramagnetic resonance spectra from Brownian dynamics trajectories: application to nitroxide side chains in proteins. Biophys J. 1996, 71 (4): 2201-12. 10.1016/S0006-3495(96)79421-3.PubMed CentralView ArticlePubMedGoogle Scholar
- LaConte Leslie, Voelz Vincent, Nelson Wendy, Enz Michael, Thomas David: Molecular Dynamics Simulation of Site-Directed Spin Labeling: Experimental Validation in Muscle Fibers. Biophysical Journal. 2002, 83: 1854-1866. 10.1016/S0006-3495(02)73950-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Beier C, Steinhoff HJ: A structure-based simulation approach for electron paramagnetic resonance spectra using molecular and stochastic dynamics simulations. Biophys J. 2006, 91 (7): 2647-64. 10.1529/biophysj.105.080051.PubMed CentralView ArticlePubMedGoogle Scholar
- Prasad Gajula MNV: PhD Thesis. [http://repositorium.uni-osnabrueck.de/bitstream/urn:nbn:de:gbv:700-2008041631/2/E-Diss781_thesis.pdf]
- Palmer SM, Playford MP, Craig SW, Schaller MD, Campbell SL: Lipid binding to the tail domain of vinculin: specificity and the role of the N and C termini. J Biol Chem. 2009, 284: 7223-7231.PubMed CentralView ArticlePubMedGoogle Scholar
- Lindahl E, Hess B, van der Spoel D: GROMACS 3.0: A package for molecular simulation and trajectory analysis. J Mol Mod. 2001, 7: 306-317.Google Scholar
- van der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC: GROMACS: Fast, Flexible and Free. J Comp Chem. 2005, 26: 1701-1718. 10.1002/jcc.20291.View ArticleGoogle Scholar
- Shen K, Tolbert CE, Guilluy C, Swaminathan VS, Berginski ME, Burridge K, Superfine R, Campbell SL: The vinculin C-terminal hairpin mediates F-actin bundle formation, focal adhesion, and cell mechanical properties. J Biol Chem. 2011, 286: 45103-45115. 10.1074/jbc.M111.244293.PubMed CentralView ArticlePubMedGoogle Scholar
- Critchley DR et al: Cytoskeletal proteins talin and vinculin in integrin-mediated adhesion. Biochem Soc Trans. 2004, 32: 831-836.View ArticlePubMedGoogle Scholar
- Saunders RM, Holt MR, Jennings L, Sutton DH, Barsukov IL, Bobkov A, Liddington RC, Adamson EA, Dunn GA, *Critchley DR: *Role of vinculin in regulating focal adhesion turnover. Eur J Cell Biol. 2006, 85: 487-500. 10.1016/j.ejcb.2006.01.014.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.