Morphs: Difference between revisions
Eric Martz (talk | contribs) |
Eric Martz (talk | contribs) |
||
(16 intermediate revisions by 3 users not shown) | |||
Line 24: | Line 24: | ||
* [[User:Luis E Ramirez-Tapia/T7 RNA polymerase|T7 RNA Polymerase]] | * [[User:Luis E Ramirez-Tapia/T7 RNA polymerase|T7 RNA Polymerase]] | ||
* [[Lipase lid morph]] | * [[Lipase lid morph]] | ||
* [[Hexokinase#Conformational_change_associated_with_substrate_binding]] | |||
* [[Calmodulin#Calmodulin_in_Motion]] | |||
* [[Human lactoferrin]] | |||
===Morphs Elsewhere=== | ===Morphs Elsewhere=== | ||
Line 29: | Line 32: | ||
* Database of Macromolecular Movements with Associated Tools for Flexibility and Geometric Analysis, [http://molmovdb.org molmovdb.org], developed by Mark Gerstein and coworkers at Yale University, USA.<ref name='flores'>The Database of Macromolecular Motions: | * Database of Macromolecular Movements with Associated Tools for Flexibility and Geometric Analysis, [http://molmovdb.org molmovdb.org], developed by Mark Gerstein and coworkers at Yale University, USA.<ref name='flores'>The Database of Macromolecular Motions: | ||
new features added at the decade mark. Flores, S. ''et al.'', Nucleic Acids Res., .34 (Database issue) D1–D6. 2006. [http://www.ncbi.nlm.nih.gov/pubmed/16381870 PubMed 16381870], [http://papers.gersteinlab.org/e-print/molmovdb-update-nar/preprint.pdf PDF].</ref><ref>Database of Macromolecular Movements with Associated Tools for Flexibility and Geometric Analysis, [http://molmovdb.org molmovdb.org], developed by Mark Gerstein and coworkers at Yale University, USA.</ref> | new features added at the decade mark. Flores, S. ''et al.'', Nucleic Acids Res., .34 (Database issue) D1–D6. 2006. [http://www.ncbi.nlm.nih.gov/pubmed/16381870 PubMed 16381870], [http://papers.gersteinlab.org/e-print/molmovdb-update-nar/preprint.pdf PDF].</ref><ref>Database of Macromolecular Movements with Associated Tools for Flexibility and Geometric Analysis, [http://molmovdb.org molmovdb.org], developed by Mark Gerstein and coworkers at Yale University, USA.</ref> | ||
* [[Protein Explorer]] developed by [[User:Eric Martz|Eric Martz]] at the University of Massachusetts, Amherst, USA.<ref name='pe'>Protein Explorer: easy yet powerful macromolecular visualization. Martz, E., Trends Biochem Sci. 27:107-9. 2002. [http://www.ncbi.nlm.nih.gov/pubmed/11852249 PubMed 11852249] developed by [[User:Eric Martz|Eric Martz]] at the University of Massachusetts, Amherst, USA.</ref> | * [[Protein Explorer]] developed by [[User:Eric Martz|Eric Martz]] at the University of Massachusetts, Amherst, USA.<ref name='pe'>Protein Explorer: easy yet powerful macromolecular visualization. Martz, E., Trends Biochem Sci. 27:107-9. 2002. [http://www.ncbi.nlm.nih.gov/pubmed/11852249 PubMed 11852249] developed by [[User:Eric Martz|Eric Martz]] at the University of Massachusetts, Amherst, USA.</ref> Note that Protein Explorer went out of service in the early 2000's because it was built around [[Chime]]. Chime was not open source, and its owners ceased to maintain it. | ||
* Database of the Morphit Pro server [http://morphit-pro.cmp.uea.ac.uk/MorphItPro/faces/faces/listMorphs.xhtml?faces-redirect=true]. and associated DynDom software [http://dyndom.cmp.uea.ac.uk/dyndom/dyndomDatabases.jsp#NRD] by Steven Hayward and coworkers at the University of East Anglia, UK. | |||
===Non-Morph Animations=== | ===Non-Morph Animations=== | ||
Line 36: | Line 40: | ||
==Why Morph?== | ==Why Morph?== | ||
<table align='right' border='0' width='285' cellpadding='10' bgcolor='#d0d0d0' hspace='8'><tr><td rowspan='2'> </td><td bgcolor='#e8e8e8'> | <table align='right' border='0' width='285' cellpadding='10' bgcolor='#d0d0d0' hspace='8'><tr><td rowspan='2'> </td><td bgcolor='#e8e8e8'> | ||
[[Image:Mage_hb.gif]]</td></tr><tr><td bgcolor='#e8e8e8'>Toggling between the carbonmonoxy and deoxy conformations of heme in hemoglobin. (This is NOT a morph.) Convergent stereo snapshots from a Kinemage. (Stops after 25 cycles; | [[Image:Mage_hb.gif]]</td></tr><tr><td bgcolor='#e8e8e8'>Toggling between the carbonmonoxy and deoxy conformations of heme in hemoglobin. (This is NOT a morph.) Convergent stereo snapshots from a Kinemage. (Stops after 25 cycles; Shift-Reload this page to restart the toggling.)</td></tr></table> | ||
The purpose of molecular morphing is to smooth the visual transition between two molecular conformations, making it easier to see and understand the structural differences between them. | The purpose of molecular morphing is to smooth the visual transition between two molecular conformations, making it easier to see and understand the structural differences between them. | ||
In contrast, predicting the actual trajectory through which a conformational change occurs is rarely, if ever, the goal of a morph. | In contrast, predicting the actual trajectory through which a conformational change occurs is rarely, if ever, the goal of a morph. | ||
Some proteins perform their functions without major conformational changes. On the other hand, some proteins must undergo major changes in secondary, tertiary, or quaternary structure in order to perform their functions. In quite a few cases, investigators have succeeded in obtaining empirically determined structures for a protein in two or more conformations. The challenge for visualization is then to be able to follow the changes in each region between the two conformations. | Some proteins perform their functions without major conformational changes. On the other hand, some proteins must undergo major changes in secondary, tertiary, or quaternary structure in order to perform their functions. In quite a few cases, investigators have succeeded in obtaining [[Empirical models|empirically determined structures]] for a protein in two or more conformations. The challenge for visualization is then to be able to follow the changes in each region between the two conformations. | ||
When the differences are small, simply toggling an image between the two states is adequate. [http://kinemage.biochem.duke.edu/ David Richardson's MAGE]<ref>[http://kinemage.biochem.duke.edu/ Kinemages] by David and Jane Richardson, Duke University, USA.</ref>, first available in 1992<ref name='history'>[http://history.molviz.org History of Macromolecular Visualization]</ref> supports visual toggling between macromolecular conformations. Hundreds of interactive molecular structure tutorials called kinemages take advantage of this capability. At the right are snapshots of hemoglobin toggled in MAGE. (MAGE is available in java applet form, and is an option for molecular displays in Proteopedia<ref>For sample pages that use the MAGE applet in Proteopedia, see [[Hemoglobin]] and [[Ribulose-1,5-bisphosphate carboxylase/oxygenase]]</ref>.) | When the differences are small, simply toggling an image between the two states is adequate. [http://kinemage.biochem.duke.edu/ David Richardson's MAGE]<ref>[http://kinemage.biochem.duke.edu/ Kinemages] by David and Jane Richardson, Duke University, USA.</ref>, first available in 1992<ref name='history'>[http://history.molviz.org History of Macromolecular Visualization]</ref> supports visual toggling between macromolecular conformations. Hundreds of interactive molecular structure tutorials called kinemages take advantage of this capability. At the right are snapshots of hemoglobin toggled in MAGE. (MAGE is available in java applet form, and is an option for molecular displays in Proteopedia<ref>For sample pages that use the MAGE applet in Proteopedia, see [[Hemoglobin]] and [[Ribulose-1,5-bisphosphate carboxylase/oxygenase]]</ref>.) | ||
Line 63: | Line 67: | ||
===True Movies=== | ===True Movies=== | ||
<table align='right' border='0' width='184' cellpadding='10' bgcolor='#d0d0d0' hspace='8'><tr><td rowspan='2'> </td><td bgcolor='#e8e8e8'> | <table align='right' border='0' width='184' cellpadding='10' bgcolor='#d0d0d0' hspace='8'><tr><td rowspan='2'> </td><td bgcolor='#e8e8e8'> | ||
[[Image:Lacrep_anim_small.gif]]</td></tr><tr><td bgcolor='#e8e8e8'>True movie of the morph shown in Jmol elsewhere near the top of this page. (Stops after 100 cycles; | [[Image:Lacrep_anim_small.gif]]</td></tr><tr><td bgcolor='#e8e8e8'>True movie of the morph shown in Jmol elsewhere near the top of this page. (Stops after 100 cycles; Shift-Reload this page to restart this movie.)</td></tr></table> | ||
A true movie is a series of static snaphots displayed sequentially in rapid succession. An example of a true movie of a morph is shown at right. Note that you '''cannot rotate the molecule with the mouse''', so you can view the morph from only the single perspective chosen by the author of the movie. Also, a separate movie must be provided for each change in rendering or coloring, while the rendering and coloring in Jmol can be easily be changed with script commands while displaying the same morph PDB file. | A true movie is a series of static snaphots displayed sequentially in rapid succession. An example of a true movie of a morph is shown at right. Note that you '''cannot rotate the molecule with the mouse''', so you can view the morph from only the single perspective chosen by the author of the movie. Also, a separate movie must be provided for each change in rendering or coloring, while the rendering and coloring in Jmol can be easily be changed with script commands while displaying the same morph PDB file. | ||
Line 102: | Line 106: | ||
==Morphing Methods== | ==Morphing Methods== | ||
All morphs begin with at least two, usually empirical, [[atomic coordinate files]] that represent different conformations of the same macromolecule. In some cases, more than two are available<ref name='vonrhein' />. Simple morphs then involve calculating a series of interpolated intermediate conformations using one of the methods below. | All morphs begin with at least two, usually [[Empirical models|empirical]], [[atomic coordinate files]] that represent different conformations of the same macromolecule. In some cases, more than two are available<ref name='vonrhein' />. Simple morphs then involve calculating a series of interpolated intermediate conformations using one of the methods below. | ||
Before making a new morph from published PDB files, search the [http://molmovdb.org/cgi-bin/movie.cgi Gallery of Morphs] where you will find thousands of morphs from jobs previously submitted to the [http://molmovdb.org Yale Morph Server] (see below). It is quite possible that your morph has already been made, and is already available there! | Before making a new morph from published PDB files, search the [http://molmovdb.org/cgi-bin/movie.cgi Gallery of Morphs] where you will find thousands of morphs from jobs previously submitted to the [http://molmovdb.org Yale Morph Server] (see below). It is quite possible that your morph has already been made, and is already available there! | ||
Line 120: | Line 124: | ||
====Proteopedia PyMOL Morpher==== | ====Proteopedia PyMOL Morpher==== | ||
'''<font color="magenta">Recommended</font>'''. Proteopedia provides a page http://proteopedia.org/cgi-bin/morph that sends a morph request to the [https://pymol.org/2/ PyMOL program] and automatically uploads the resulting multi-model file to proteopedia. An example of a morph made by the server is [http://proteopedia.org/wiki/index.php/Image:Morph_semet_apo-chaind.pdb_zn_complex-chaina.pdb_u4043.pdb here], and it is shown on this [http://proteopedia.org/wiki/index.php/User:Kristian_Koski/P4H#flexible_loops page]. | '''<font color="magenta">Recommended</font>'''. Proteopedia provides a page http://proteopedia.org/cgi-bin/morph that sends a morph request to the [https://pymol.org/2/ PyMOL program] and automatically uploads the resulting multi-model file to proteopedia. An example of a morph made by the server is [http://proteopedia.org/wiki/index.php/Image:Morph_semet_apo-chaind.pdb_zn_complex-chaina.pdb_u4043.pdb here], and it is shown when you click the green links on this [http://proteopedia.org/wiki/index.php/User:Kristian_Koski/P4H#flexible_loops page]. | ||
Advantages: | Advantages: | ||
Line 132: | Line 136: | ||
The [http://fatcat.godziklab.org FATCAT Structural Alignment Server] also produces "chemically possible" morphs. After an optionally flexible alignment (permitting twists at hinge points determined by FATCAT), a linear interpolation is done between the aligned models. Then "the intermediate structures are optimized by energy gradient minimization employing a reduced representation force field." | The [http://fatcat.godziklab.org FATCAT Structural Alignment Server] also produces "chemically possible" morphs. After an optionally flexible alignment (permitting twists at hinge points determined by FATCAT), a linear interpolation is done between the aligned models. Then "the intermediate structures are optimized by energy gradient minimization employing a reduced representation force field." | ||
Limitations: FATCAT does only one chain from each model. | Limitations: FATCAT does only one chain from each model. Chains are re-numbered so it is time consuming to relate the new numbers to the original sequence numbers. In the resulting morph, it appears that some changes may be made to at least one model even when a "rigid" alignment is specified. | ||
====Biomolecular Morphing by Kleywegt at Uppsala==== | ====Biomolecular Morphing by Kleywegt at Uppsala==== | ||
Line 139: | Line 143: | ||
====Yale Morph Server==== | ====Yale Morph Server==== | ||
<!--''SADLY, AS OF 2020, THE YALE MORPH SERVER AND GALLERY OF MOLECULAR MOVEMENTS HAS BEEN OFFLINE FOR A NUMBER OF YEARS.''--> | |||
''NOT WORKING AGAIN in January 2023. STILL NOT WORKING November 2024.'' | |||
'' | ''WORKING AGAIN in February, 2021. Their animation player doesn't work because it still requires Java, but you can download the morph PDB file and animate it in Proteopedia or Jmol. See [[#Animating Morph PDB Files|instructions below]].'' | ||
The [http://molmovdb.org Yale Morph Server] from Mark Gerstein's group. In addition to the chemically possible nature of the results, this server is very attractive because it totally automates both the creation of the multiple-model morph PDB data file, and also its visualization. The server does the structural alignment of the starting PDB files, and can handle modest differences in sequence between the initial and final PDB files. It then does a linear interpolation followed by some energy minimization to render each frame ''chemically possible''. | |||
Morph2 (version 2) of the Yale Morph Server allows you to pre-align (or have the server align) your structures, and gives you the option of minimizing, or not. | |||
===Linear Interpolation=== | ===Linear Interpolation=== | ||
In this method, a series of intermediate models is created in which each atom moves in a straight line from is initial position to its final position. Linear interpolation by itself (without any subsequent energy minimization) is relatively easy to calculate, and often suffices, but also has a number of limitations. When animated, the interpolated intermediate conformations often greatly help in visualizing the differences between the initial and final empirical models. However, bond lengths and angles become unrealistic, domains may artifactually shrink, expand, or distort, and chains may even pass through each other during the interpolated movements. When there are substantial changes in secondary structure or large movements of domains, these artifacts make the eye unable to follow the details of the conformational change, and linear interpolation becomes unsatisfactory. | In this method, a series of intermediate models is created in which each atom moves in a straight line from is initial position to its final position. Linear interpolation by itself (without any subsequent energy minimization) is relatively easy to calculate, and often suffices, but also has a number of limitations. When animated, the interpolated intermediate conformations often greatly help in visualizing the differences between the initial and final [[empirical models]]. However, bond lengths and angles become unrealistic, domains may artifactually shrink, expand, or distort, and chains may even pass through each other during the interpolated movements. When there are substantial changes in secondary structure or large movements of domains, these artifacts make the eye unable to follow the details of the conformational change, and linear interpolation becomes unsatisfactory. | ||
The presence of obvious artifacts of linear interpolation, such as objects passing through each other, can be useful because it reminds viewers that the morph's value is only to help visualize the details of the conformational change -- and is not intended to suggest the actual trajectory of change. An example is the morph of [[Recoverin, a calcium-activated myristoyl switch|Recoverin]], in which the N-linked myristoyl group passes through a protein chain. | The presence of obvious artifacts of linear interpolation, such as objects passing through each other, can be useful because it reminds viewers that the morph's value is only to help visualize the details of the conformational change -- and is not intended to suggest the actual trajectory of change. An example is the morph of [[Recoverin, a calcium-activated myristoyl switch|Recoverin]], in which the N-linked myristoyl group passes through a protein chain. | ||
Line 156: | Line 164: | ||
====MS-DOS morph2.exe==== | ====MS-DOS morph2.exe==== | ||
A very simple program, named '''morph2.exe''', performs linear interpolation and has been freely available since the late 1990's in the ''PDBTools'' package by Eric Martz.<ref>[http://www.umass.edu/microbio/rasmol/pdbtools.htm PDBTools] by Eric Martz include a simple linear interpolation morphing program. This program operates only in MS-DOS (under MS Windows), but the C source code is included. Available within ''Protein Explorer'' are [http://www.umass.edu/microbio/chime/morpher/morphmtd.htm step by step instructions] for making linear interpolation morphs.</ref>. This program, which was used for most of the morphs made by Martz (including the Lac repressor morph near the top of this page), requires that the two starting PDB files contain exactly the same atoms in exactly the same order. (Achieving this usually requires some hand editing of the PDB files.) Its advantages include that it is straightforward to include ligand in the morph, or at the end of the morph, and to include multiple protein and nucleic acid chains. [http://www.umass.edu/microbio/chime/morpher/morphmtd.htm Step by step instructions] are available within ''Protein Explorer''. While these instructions were written for Chime, rather than Jmol, everything there applies equally to Jmol except the last section on ''Playback Scripts''. Caution: check the results -- morph2.exe appears to fail with very large PDB files but its limits have not been defined. | A very simple program, named '''morph2.exe''', performs linear interpolation and has been freely available since the late 1990's in the ''PDBTools'' package by Eric Martz.<ref>[http://www.umass.edu/microbio/rasmol/pdbtools.htm PDBTools] by Eric Martz include a simple linear interpolation morphing program. This program operates only in MS-DOS (under MS Windows), but the C source code is included. Available within ''Protein Explorer'' are [http://www.umass.edu/microbio/chime/morpher/morphmtd.htm step by step instructions] for making linear interpolation morphs.</ref>. This program, which was used for most of the morphs made by Martz (including the Lac repressor morph near the top of this page), requires that the two starting PDB files contain exactly the same atoms in exactly the same order. (Achieving this usually requires some hand editing of the PDB files.) Its advantages include that it is straightforward to include ligand in the morph, or at the end of the morph, and to include multiple protein and nucleic acid chains. [http://www.umass.edu/microbio/chime/morpher/morphmtd.htm Step by step instructions] are available within ''Protein Explorer''. While these instructions were written for Chime, rather than Jmol, everything there applies equally to Jmol except the last section on ''Playback Scripts''. Caution: check the results -- morph2.exe appears to fail with very large PDB files but its limits have not been defined. | ||
===Combination of rigid body movement and linear interpolation=== | |||
To explore conformational changes on the fly within Jmol, the storymorph suite defines several Jmol functions to superimpose and morph related structures. The method is described in detail at [[Jmol/Storymorph]], together with some demonstrations and interactive examples. | |||
===Animating Morph PDB Files=== | ===Animating Morph PDB Files=== | ||
If you upload your multiple-model morph PDB file to Proteopedia, animating it is as simple as checking that option in Proteopedia's | If you upload your multiple-model morph PDB file to Proteopedia, animating it is as simple as checking that option in Proteopedia's [[Scene authoring tools]]. | ||
<blockquote> | <blockquote> | ||
For those who wish to make their own animation scripts in Jmol outside of Proteopedia (e.g. in the [http://bioinformatics.org/jmol-tutorials Jmol Tutorial-Authoring Template]), the script required to animate a multiple-model PDB file in Jmol is much simpler than what was needed in Chime. Here are the commands needed in Jmol: | For those who wish to make their own animation scripts in [[Jmol/Application|Jmol]] outside of Proteopedia (e.g. in the [http://bioinformatics.org/jmol-tutorials Jmol Tutorial-Authoring Template]), the script required to animate a multiple-model PDB file in Jmol is much simpler than what was needed in Chime. Here are the commands needed in Jmol: | ||
<pre> | <pre> | ||
anim mode palindrome; | anim mode palindrome; |