The scripts and input files that accompany this demo can be found in the
demos/protocol_captures directory of the Rosetta weekly releases.
This is a protocol capture, and represents the protocol at a fixed point in time. It may not work with the current version of Rosetta.
KEYWORDS: MEMBRANES UTILITIES
Visualizing the position and orientation of a biomolecule in the membrane bilayer (embedding) is an important part of analyzing membrane models from Rosetta. This visualization allows the user to distingish conformations with poor embeddings from conformations with native-like/reasonable embeddings.
This application displays a set of two parallel planes in PyMOL, representing the position and geometry of the membrane bilayer. During a simulation, Rosetta extracts the membrane center, normal and thickness from the MEM residue, calculates the position of the planes, and sends this information to PyMOL in real time. PyMOL then uses this information to draw two CGO plane objects representing the membrane.
This tool is part of a standalone application and can also be used in combination with any JD2-supported Rosetta application.
Publication describing the method:
-show_simulation_in_pymol 0 flag with any Rosetta Membrane Framework application
Two inputs are required for using the standalone visualization app:
Steps for generating these inputs are found below. A set of example inputs can also be found in example_inputs/. Here, 1c3w is used as an example:
PDB File: If an output model is not already available from Rosetta, generate a PDB file where the membrane protein structure is transformed into PDB coordinates (z-axis is membrane normal). This can be done either by downloading the transformed PDB directly from the PDBTM website (http://pdbtm.enzim.hu/) or by downloading a PDB file from the PDB and running it through the PPM server (http://opm.phar.umich.edu/server.php).
Span File: Generate a spanfile from the PDB structure using the spanfile_from_pdb application described in the MP_spanfile-from-pdb protocol capture in Rosetta/demos/protocol_captures/2014. An example commandline using 1c3w is also provided here:
Rosetta/main/source/bin/spanfile_from_pdb.linuxgccrelease -database /path/to/db -in:file:s example_inputs/1c3w_tr.pdb
For this example, this command will produce 1 output files: = 1c3w_tr.span: Spanfile containing predicted trans-membrane spans
Here, we describe the steps required to run the MP_PyMOLViewer protocol. As an example, all steps use the PDB 1c3w:
Required Options: Options (flags) needed to run this application. A file with these flags, pymol_flags, is also provided for 1c3w in this demo:
flags descriptions -------------------------------------------------------------------------------------------------- -in:file:s <pdbfile> Input PDB Structure: PDB file for protein structure -membrane_new:setup:spanfiles Spanfile describing spanning topology of starting structure for full symmetric structure -show_simulation_in_pymol 0 Use the PyMOL viewer to visualize membrane planes for structures -keep_pymol_simulation_history 1 Keep pymol frames for making movies/replaying simulations (optional)
Startup the PyMOL PyRosetta Session:
Run the PyMOLPyrosettaServer.py script using the following command line in the pymol window:
Once run, a message should appear in the PyMOL terminal window indicating the server was initialized successfully.
Run Rosetta application:
From the regular terminal, run the standalone application or other membrane framework apps from the command line:
Rosetta/main/source/bin/view_membrane_protein.linuxgccrelease -database /path/to/db @pymol_flags
Within ~10-20 seconds, 2 parallel planes (PyMOL object entitled membrane_planes) will appear in the PyMOL session.
Note: In earlier versions of PyMOL, small lines will be drawn to the MEM residue as this is a HETATM in the PDB. These can be turned off manually by hiding the line view for the MEM residue.
The following outputs will be generated from the standalone pymol viewer application. A version of these outputs are also provided in the example_outputs/ directory:
1c3w.pse: Example pymol session file including membrane planes objects
1c3w.png: Example image of PDB 1c3w (bacteriorhodopsin) embedded in the membrane (from session file)
1c3w_tr_0001.pdb: Output PDB file containing membrane position information in MEM residue
Note: Rosetta will also output a score file (score.sc) which is not needed for this analysis.
Baugh EH, Lyskov S, Weitzner BD, Gray JJ (2011) Real-Time PyMOL Visualization for Rosetta and PyRosetta. PLoS ONE 6: e21931.
DeLano W (n.d.) The PyMOL Manual: Compiled Graphics Objects (CGOs) and Molscript Ribbons. Available: http://pymol.sourceforge.net/newman/user/toc.html.