- Open Access
TMM@: a web application for the analysis of transmembrane helix mobility
BMC Bioinformatics volume 8, Article number: 232 (2007)
To understand the mechanism by which a protein transmits a signal through the cell membrane, an understanding of the flexibility of its transmembrane (TM) region is essential. Normal Mode Analysis (NMA) has become the method of choice to investigate the slowest motions in macromolecular systems. It has been widely used to study transmembrane channels and pumps. It relies on the hypothesis that the vibrational normal modes having the lowest frequencies (also named soft modes) describe the largest movements in a protein and are the ones that are functionally relevant. In particular NMA can be used to study dynamics of TM regions, but no tool making this approach available for non-experts, has been available so far.
We developed the web-application TMM@ (TransMembrane α-helical Mobility analyzer). It uses NMA to characterize the propensity of transmembrane α-helices to be displaced. Starting from a structure file at the PDB format, the server computes the normal modes of the protein and identifies which helices in the bundle are the most mobile. Each analysis is performed independently from the others and results can be visualized using only a web browser. No additional plug-in or software is required. For users who would like to further analyze the output data with their favourite software, raw results can also be downloaded.
We built a novel and unique tool, TMM@, to study the mobility of transmembrane α-helices. The tool can be applied to for example membrane transporters and provides biologists studying transmembrane proteins with an approach to investigate which α-helices are likely to undergo the largest displacements, and hence which helices are most likely to be involved in the transportation of molecules in and out of the cell.
α-helical transmembrane (TM) proteins represent approximately 20–30% of all open reading-frames in the genome of complex organisms. They are involved in many biological processes such as sight, smell, muscle contraction, photosynthesis, etc . Their signalling function is most often achieved by movements of the helices constituting the transmembrane bundle; the movements can be of different nature, involving the whole bundle like in the case of the mechanosensitive channel  or individual helices displacements such as those accomplished by the Ca-ATPase to transport calcium ions through the sarcoplasmic reticulum membrane . Even though the structural changes involved are now pretty well understood for these two proteins, there are still many for which the available structural information is not sufficient to understand the mechanism(s) by which signals are transmitted along the transmembrane region.
Molecular modeling is the approach of choice to study protein dynamics. Molecular Dynamics (MD) is perhaps the most widely used technique but it is computer demanding and simulations of slow protein dynamics still take months even on high-performance computers. Approximate Normal Mode Analysis (NMA) is the method of choice to investigate the slowest motions in macromolecular systems . Because of its very modest requirements in terms of computer power, it is especially useful for large biomolecular assemblies. It has been shown for many proteins that the vibrational normal modes having the lowest frequencies describe the largest movements in a protein and are the ones that are functionally relevant [4–9]. In recent years, it has become available to non-specialists thanks to the development of several web applications enabling traditional NMA approaches [10–15] for the study of large amplitude movements of protein domains.
NMA has proven useful in identifying mobile helices in TM bundles and describing, for example, the structural modifications accompanying the transport of calcium by the Ca-ATPase or the opening/closing mechanism of the mechanosensitive channel (MscL) [14, 16, 17]. Here we present a novel unique web application, the TransMembrane α-helical Movement Analyzer (TMM@) that analyzes the mobility of α-helices in transmembrane bundles. TMM@ performs a calculation of the normal modes and analyzes the displacement of the TM α-helices by calculating the overlap between the modes and displacement vectors describing movements with relevance to the transport function. The outcome is a ranking of the TM helices according to their mobility. For proteins where it is believed that the transport function involves helices movements, the characterization of their mobility will help elucidating the way the signal is transmitted through the cell membrane .
The architecture of TMM@ is similar to that of WEBnm@ . The web-interface is written using the DTML language of the Zope  webserver. The analysis core, written in Python, runs under the BIAZ application server . The calculation of the normal modes and the analysis of the displacements of the helices are implemented in Python and make extensive use of MMTK . The normal modes are calculated on our server (AMD Opteron, 2 core, 2.4 Ghz, 8GB ram). The calculation of the whole set of normal modes for a protein of 1000 residues (i.e. 3000 modes) is performed in about 10 minutes. Overlap plots are generated using the R  package and the Rpy library . We have embedded Jmol , an interactive molecular viewer (Java applet), into the web page to visualize the protein. VMD  (1.8.1) is used to produce VMD state files which can be downloaded if a more detailed visualization is preferred.
1. Normal modes calculations
A normal mode analysis (NMA) consists of the diagonalization of the matrix of the second derivatives of the energy with respect to the displacements of the atoms, in mass-weighted coordinates (Hessian matrix). The eigenvectors of the Hessian matrix are the normal modes, and the associated eigenvalues are the squares of the associated frequencies. We use the approximate normal modes calculation method developed by Hinsen  and implemented in the MMTK package . This method represents the low-frequency domain motions very well at negligible computational cost. The force field used has been described in reference . It uses only the Cα atoms of the protein; each atom is assigned the weight of the whole residue it represents.
Briefly, the functional form of the force field is
V (r) is the harmonic pair potential describing the interaction between the Cα atoms:
where is the pair distance vector (Ri - Rj) in the input configuration and k is the pair force constant:
2. Identification of trans-membrane α-helices bundle
TMM@ uses DSSP  and its own filter algorithm to produce a list of all α-helices present in the submitted protein structure and to identify the TM bundle. The filter algorithm makes use of the following structural properties: helix length, distance between helices, hydrophobicity, and the angle between helices. As the filter algorithm is based on empirical parameters, we recommend that each user review and correct if necessary the list suggested by TMM@.
3. Defining α-helical mobility
The projection of a normal mode vector onto a displacement vector defines the contribution of the normal mode to the given displacement. In TMM@ we define four different movements of relevance for the transport function: (i) rotation and (ii) translation of individual helices around and along their axis, respectively, (iii) slide of the α-helices perpendicular to the helix axis towards/away from the centre of the bundle, and (iv) tilt of helices perpendicular to the helix axis away from the centre of bundle, and (v) rotation of the helices around the bundle axis. The axis of a α-helix is defined as the principal axis of inertia of the Cα-atoms of the amino acids forming the helix, the axis of the bundle is defined as the principal axis of inertia of the Cα atoms of all helices in the bundle. The rotation vector on each Cα-atom of the α-helix is calculated as the cross-product between a unit vector collinear to the helix axis and the distance vector between the Cα-atom and the centre of mass of the α-helix. The translation vectors of the α-helices have a component for each Cα-atom, collinear to the axis. The rotation of the bundle is defined by the cross product between the axis of the bundle and the distance vector between the bundle centre and the helix centre. The slide vector is the cross product between the bundle rotation vector and the helix axis. The tilt vector is calculated in the same way as the slide vector, but with decreasing magnitude for residues closer to the centre of the helix, and opposite direction on the other side of the centre. Hence, we tilt the helix around the centre of the helix, directly away and towards the bundle axis.
The projections are defined by
p i = d·e i (4)
where e i is the normal mode vector of mode i, and d is the displacement vector (i.e. rotation, translation, slide or tilt of individual helices, bundle rotation). This satisfies the relation
because the normal mode vectors form a basis of configuration space (N is the number of atoms). Thus pi2 is interpreted as the contribution of mode i to the motion described by d. For each helix, the calculation of the cumulative overlap of one given displacement vector and all modes thus yields a curve that increases from 0 to 1 (y axis) when it is computed over all modes (on the x axis). The modes are ranked following increasing frequencies. Frequency is inversely proportional to amplitude. If α-helix H1 shows a plot (e.g. grey dotted line on Fig. 1) that reaches an asymptotic behaviour for fewer modes than the plot of another helix H2 (e.g. black plain line on Fig. 1), it means that the movement of H1 following d is of higher amplitude than H2 following d. As a consequence H1 is considered to be more mobile than H2. This is illustrated on Figure 1 where d is the rotation of each helix around its own axis.
Song et al recently described a new method to evaluate the overlap between a set of normal modes and a given transconformation. They point out that the global conformation displacement is a finite motion, while the normal mode motions are infinitesimal motions. Therefore, for large conformational difference, the global direction may have little to do with the initial direction. They thus propose to use an infinitesimal version of the transconformation instead of the global transconformation, and show that it is more reliable than the usual difference vector, especially when the displacement is large. However, our program is meant to evaluate the contribution of the modes to helices displacements that we define ourselves, using only one structure. We thus avoid the problems mentioned by Song et al.
III. Results and Discussion
TMM@ can be used by anybody with access to an internet browser; all results are presented within the webpages of TMM@ and no additional programs are needed on the computer of the user. The user starts by defining which structure is to be analyzed; from the main page, he/she is offered the possibility to upload his/her own protein structure file (in the Protein Data Bank format) or specify a PDB code (a local copy of the PDB database is maintained on our server). A third field allows the user to enter a job ID to continue working on a previously started calculation. Pressing the submit button will launch the calculation of the normal modes. Then TMM@ identifies all α-helices present in the protein and lists them on the next page where it also suggests which ones are TM helices. The user is offered the possibility to correct the definition of the TM α-helices with the aid of a Jmol applet. Alternatively the user can download a vmd state file  for a more detailed visualization of the predicted TM helices. For each TM helix in the bundle approved by the user, TMM@ defines the displacement vectors and calculates the overlap between each of them and the set of normal modes. The overlaps are then plotted using the R package and the Rpy library. The cumulative squared overlaps are plotted against mode numbers, four plots are drawn (available in the PDF and PNG formats); one for each type of displacement. In addition, the user can retrieve the raw results in a text file. Each curve on the plot corresponds to a TM α-helix, which is identified with a unique identifier, colour and line type. Displaying all helices on the same plot helps comparing their mobility. It is reasonable to believe that the most mobile helices will be involved in the structural modifications accompanying, for example, the transport of a molecule or ion through a tight TM bundle. In many cases, the transport function is indeed a dynamical process during which the protein undergoes structural rearrangements.
We have tested TMM@ on more than 20 transmembrane proteins taken from the PDB, and representing different families. In what follows, we describe two examples and use them to illustrate the different steps performed by TMM@: the SERCA1 calcium ATPase and the mechanosensitive channel (MscL). Snapshots of the application are given in Figure 2 and Figure 3, respectively. The calcium pump (Ca-ATPase) transports calcium from the cytoplasmic side (outside cell) to the lumenal side (inside cell). There is however no obvious channel in the protein leading to the lumenal side and it has been shown that the ion transport implies movements of the α-helices. After uploading the x-ray structure (1su4) of the E1Ca form of Ca-ATPase (Figure 2a) the normal modes will be computed. The filter algorithm outputs a list of 45 α-helices, in which 13 of these are suggested to be in the TM region (Figure 2b). However, since we know that only 10 of them are TM helices, we wish to correct the definition given by TMM@. The Jmol applet provided on the page aids in this work. The predicted TM helices of the Ca-ATPase are listed in Table 1. The next step is to submit the list of TM helices so that the overlap calculation will be performed, resulting in 5 plots (Figure 2c and 2d) described above. The overlap plot for the rotation of each helix around its individual axis (Figure 2d) shows that less than 300 modes are enough to describe 60% of the rotation of helices 3, 4 and 9 (M1-M3). Similarly, 60% of the rotation of helices 41 and 43 (M9 and M10) can be described with only 300 modes. Conversely, over 600 modes are needed to describe 60% of the rotation for helices 12, 31, 32, 34, 38 (M4 to M8). This means that helices number M1 to M3 can undergo larger amplitude displacements than helices M4 to M8. This leads to the hypothesis that they might play a role in the uptake/release of calcium ions since it is known that the ion transport requires displacements of the TM α-helices. This result is in agreement with the available x-ray structures of the calcium pump and in particular ref.  where Toyoshima et al. describe movements of M1 to M3, and in particular that M2 and M1 are pulled towards the cytoplasm by one and two turns of a α-helix, respectively. These complex movements are a combination of the rotation and translation around and along α-helical axes that TMM@ investigates opportunely.
Figure 3 (a–d) shows the resulting overlap plots for MscL (1msl). As expected, our calculations reproduce the symmetry of the molecule, i.e. all M2 helices have the same mobility in rotation and translation (Fig. 3a and 3b) around/along their own axis, the same applies to M1 helices. Therefore the plots show superimposed lines for symmetrically equivalent helices. For a bundle without that type of symmetry, like in the Ca-ATPase which contains 10 TM helices, the corresponding plots should contain 10 distinct lines. TMM@ also identifies the iris-like movement described by others [2, 16]. The plots on Figure 3c and 3d indeed reveal that a few low-frequency modes (modes #1, #7 and #20) describe a rotation of the helices around the bundle axis. The same type of movement was identified for the calcium pump (Figure 2c and ref. ).
We have successfully implemented and developed a unique tool for analysing the mobility of α-helical TM segments in proteins. We have tested it on a number of TM α-helical proteins and have compared the results, whenever possible, with existing structural data. In these cases, the helices that TMM@ identifies as being the most mobile are known experimentally to be involved in the protein function. We thus believe that this approach has a strong predictive power.
Using normal mode calculations as the basis of the tool makes it reliable, robust and fast. Providing the service in a user friendly web interface will make it easy to use, even for non-specialists. It thus provides biologists studying transmembrane proteins a unique tool for determining which helices undergo the largest displacements, and hence which might be involved in the transportation of molecules in and out of the cell. Considering that trans-membrane proteins are of vital importance for cell life, TMM@ might be a tool with great value. TMM@ is available from the website of the Norwegian Bioinformatics Platform http://www.bioinfo.no/tools/tmma.
Availability and requirements
Project name: TMM@
Project home page: http://www.bioinfo.no/tools/tmma
Operating system: Platform independent; tested on Windows XP (Firefox 2.0, MS Explorer 6.0 and 7.0), Mac OSX (Firefox 2.0, Safari 2.0), CentOS and Ubuntu Linux (Firefox 2.0, Opera 9.0)
Programming language: Python
Other requirements: Java
Curran AR, Engelman DM: Sequence motifs, polar interactions and conformational changes in helical membrane proteins. Curr Opin Struct Biol. 2003, 13 (4): 412-417. 10.1016/S0959-440X(03)00102-7.
Sukharev S, Anishkin A: Mechanosensitive channels: what can we learn from 'simple' model systems?. Trends Neurosci. 2004, 27 (6): 345-351. 10.1016/j.tins.2004.04.006.
Toyoshima C, Nomura H, Tsuda T: Lumenal gating mechanism revealed in calcium pump crystal structures with phosphate analogues. Nature. 2004, 432 (7015): 361-368. 10.1038/nature02981.
Hayward S: Normal mode analysis of biological molecules. Computational biochemistry and biophysics. Edited by: Becker OM, MacKerell AD, Roux B, Watanabe M. 2001, New-York , Marcel Dekker, Inc., 153-168.
Levitt M, Sander C, Stern PS: Protein normal-mode dynamics: trypsin inhibitor, crambin, ribonuclease and lysozyme. J Mol Biol. 1985, 181 (3): 423-447. 10.1016/0022-2836(85)90230-X.
Marques O, Sanejouand YH: Hinge-bending motion in citrate synthase arising from normal mode calculations. Proteins. 1995, 23 (4): 557-560. 10.1002/prot.340230410.
Mouawad L, Perahia D: Motions in hemoglobin studied by normal mode analysis and energy minimization: evidence for the existence of tertiary T-like, quaternary R-like intermediate structures. J Mol Biol. 1996, 258 (2): 393-410. 10.1006/jmbi.1996.0257.
Hayward S, Kitao A, Berendsen HJ: Model-free methods of analyzing domain motions in proteins from simulation: a comparison of normal mode analysis and molecular dynamics simulation of lysozyme. Proteins. 1997, 27 (3): 425-437. 10.1002/(SICI)1097-0134(199703)27:3<425::AID-PROT10>3.0.CO;2-N.
Ma J: Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. Structure (Camb). 2005, 13 (3): 373-380. 10.1016/j.str.2005.02.002.
Suhre K, Sanejouand YH: ElNemo: a normal mode web server for protein movement analysis and the generation of templates for molecular replacement. Nucleic Acids Res. 2004, 32 (Web Server issue): W610-4. 10.1093/nar/gkh368.
Echols N, Milburn D, Gerstein M: MolMovDB: analysis and visualization of conformational change and structural flexibility. Nucleic Acids Res. 2003, 31 (1): 478-482. 10.1093/nar/gkg104.
Wako H, Kato M, Endo S: ProMode: a database of normal mode analyses on protein molecules with a full-atom model. Bioinformatics. 2004, 20 (13): 2035-2043. 10.1093/bioinformatics/bth197.
Hollup SM, Salensminde G, Reuter N: WEBnm@: a web application for normal mode analyses of proteins. BMC Bioinformatics. 2005, 6 (1): 52-10.1186/1471-2105-6-52.
Yang LW, Liu X, Jursa CJ, Holliman M, Rader AJ, Karimi HA, Bahar I: iGNM: a database of protein functional motions based on Gaussian Network Model. Bioinformatics. 2005, 21 (13): 2978-2987. 10.1093/bioinformatics/bti469.
Delarue M, Lindahl E: Normal mode calculation and visualisation using Pymol. 2004, [http://lorentz.immstr.pasteur.fr/nma/]
Valadie H, Lacapere JJ, Sanejouand YH, Etchebest C: Dynamical properties of the MscL of Escherichia coli: a normal mode analysis. J Mol Biol. 2003, 332 (3): 657-674. 10.1016/S0022-2836(03)00851-9.
Shrivastava IH, Bahar I: Common mechanism of pore opening shared by five different potassium channels. Biophys J. 2006, 90 (11): 3929-3940. 10.1529/biophysj.105.080093.
Reuter N, Hinsen K, Lacapere JJ: Transconformations of the SERCA1 Ca-ATPase: a normal mode study. Biophys J. 2003, 85 (4): 2186-2197.
Zope Open Source web application server. [http://www.zope.org/]
Saelensminde G: The Biaz application server. 2003, [http://www.bioinfo.no/tools/Biaz]
Hinsen K: The Molecular Modeling Toolkit : a new approach to molecular simulations. J Comput Chem. 2000, 21 (2): 79-85. 10.1002/(SICI)1096-987X(20000130)21:2<79::AID-JCC1>3.0.CO;2-B.
R. [http://www.r-project.org/]Copyright (C) 1989, 1991 Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
Moreira W, Warnes GR: RPy. 2003, [http://rpy.sourceforge.net/]
Gezelter JD: Jmol - http://jmol.sourceforge.net/. 2006, [http://jmol.sourceforge.net/]
Humphrey W, Dalke A, Schulten K: VMD - Visual Molecular Dynamics. J Molec Graphics. 1996, 14 (1): 33-38. 10.1016/0263-7855(96)00018-5.
Hinsen K: Analysis of domain motions by approximate normal mode calculations. Proteins. 1998, 33: 417-429. 10.1002/(SICI)1097-0134(19981115)33:3<417::AID-PROT10>3.0.CO;2-8.
Hinsen K, Petrescu AJ, Dellerue S, Bellissent-Funel MC, Kneller GR: Harmonicity in slow protein dynamics. Chem Phys. 2000, 261: 25-37. 10.1016/S0301-0104(00)00222-6.
Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22 (12): 2577-2637. 10.1002/bip.360221211.
Song G, Jernigan RL: An enhanced elastic network model to represent the motions of domain-swapped proteins. Proteins. 2006, 63 (1): 197-209. 10.1002/prot.20836.
Toyoshima C, Mizutani T: Crystal structure of the calcium pump with a bound ATP analogue. Nature. 2004, 430 (6999): 529-535. 10.1038/nature02680.
Funding for this work was provided by FUGE (Norwegian functional genomics program of the Norwegian Research Council) through the technology platform for bioinformatics. We wish to thank Gisle Saelensminde and Svenn Helge. Grindhaug for technical assistance.
LS designed TMM@ and wrote the code. IJ was involved in the planning of TMM@ and contributed to the writing of the manuscript. NR supervised the project, and edited the manuscript. This work is a truly collaborative effort of all three authors. All authors read and approved the final manuscript.
About this article
Cite this article
Skjaerven, L., Jonassen, I. & Reuter, N. TMM@: a web application for the analysis of transmembrane helix mobility. BMC Bioinformatics 8, 232 (2007). https://doi.org/10.1186/1471-2105-8-232
- Normal Mode
- Normal Mode Analysis
- Helix Axis
- Mechanosensitive Channel