Jmol/Visualizing membrane position: Difference between revisions
Eric Martz (talk | contribs) No edit summary |
No edit summary |
||
(12 intermediate revisions by one other user not shown) | |||
Line 17: | Line 17: | ||
You can download the pseudoatom-enhanced PDB file from the server, and then [[Help:Uploading_molecules|upload it to Proteopedia]] for use in a ''Molecular Scene''. | You can download the pseudoatom-enhanced PDB file from the server, and then [[Help:Uploading_molecules|upload it to Proteopedia]] for use in a ''Molecular Scene''. | ||
In the case of 5LiL, I had difficulty getting the pseudoatoms to show from a green link (even though they displayed in the [[SAT]]). In the PDB file, the pseudoatoms followed a MASTER record. I deleted the MASTER record and all CONECT records, and then I got the | In the case of 5LiL, I had difficulty getting the pseudoatoms to show from a green link (even though they displayed in the [[SAT]]). In the [[PDB file]], the pseudoatoms followed a MASTER record. I deleted the MASTER record and all CONECT records, and then I got the scene shown, using [[Image:5lil_opm2.pdb]]. | ||
* [http://FirstGlance.Jmol.Org/fg.htm?mol=http%3A//proteopedia.org/wiki/images/d/da/5lil_opm2.pdb View | * [http://FirstGlance.Jmol.Org/fg.htm?mol=http%3A//proteopedia.org/wiki/images/d/da/5lil_opm2.pdb View 5lil_opm2.pdb in FirstGlance in Jmol]: Use the ''Ligands+'' button to toggle visibility of the pseudoatoms. | ||
The pseudoatoms have residue name DUM (for dummy atoms). This can be used in the SAT to select them. The red pseudoatoms are oxygens, and the blue pseudoatoms are nitrogens. Thus, coloring them by ''chemical element'' ([[CPK]]) will make them red and blue. | |||
==Methods: Translucent Cylinder== | ==Methods: Translucent Cylinder== | ||
Line 26: | Line 28: | ||
Viewing the model in [[FirstGlance in Jmol]], you can use ''Hydrophobic/Polar'' (Views tab) to visualize the domain with a hydrophobic surface. Touching or clicking on atoms at the domain boundaries will identify those amino acids. If you wish, you can center a particular residue, then zoom in, and see atomic detail with ''Vines/Sticks'' (Views tab). | Viewing the model in [[FirstGlance in Jmol]], you can use ''Hydrophobic/Polar'' (Views tab) to visualize the domain with a hydrophobic surface. Touching or clicking on atoms at the domain boundaries will identify those amino acids. If you wish, you can center a particular residue, then zoom in, and see atomic detail with ''Vines/Sticks'' (Views tab). | ||
For 5LiL, I chose thr607 | For 5LiL, I chose thr607 and leu250. These will be used in the example commands below, but you must substitute appropriate amino acids for your model. | ||
5LiL is a homo-dimer of chains A and B. Since I did not specify a chain in the commands below, the position of, for example, thr607 will be the average position of all the atoms in thr607 in both chains. That average position will be at the axis of the rod-shaped molecule. Your molecule may require a different strategy. | |||
===Generate the Cylinder=== | ===Generate the Cylinder=== | ||
Line 32: | Line 36: | ||
To generate the cylinder in the [[SAT]], click on the button below JSmol ''Advanced: Open JSmol Console''. The command to generate the cylinder will have this form: | To generate the cylinder in the [[SAT]], click on the button below JSmol ''Advanced: Open JSmol Console''. The command to generate the cylinder will have this form: | ||
<pre>draw cyl1 cylinder diameter 100.0 color gray (thr607) (leu250)</pre> | <pre>draw cyl1 cylinder diameter 100.0 color gray (thr607) (leu250)</pre> | ||
Adapt this command to your situation, copy and paste it into the LOWER box of the Console, and click the ''Run'' button at the lower-left of the console box. | |||
*''cyl1'' is an arbitrary name for this cylinder. | |||
*You may adjust the diameter and color to your taste. (Gray is a [[DRuMS|standard color]] for hydrophobic. Yellow is sometimes used with the mnemonic ''butter'' in mind.) | |||
*If you wish to delete all cylinders before making a new one, the command is | |||
<pre>draw * delete</pre> | |||
*To make the cylinder translucent, the command is | |||
<pre>draw cyl1 color translucent -1</pre> | |||
Once you have the desired cylinder in the SAT, save the scene, insert a green link in the wikitext box, and save the page. Then test the green link. | |||
After displaying the cylinder with your new green link, if you then click a different green link, the cylinder may persist, unwanted. To fix this: | |||
# In the SAT, load the other scene. It doesn't matter if the cylinder is not visible in the SAT. | |||
# Open the command console. | |||
# Enter and Run the above ''delete'' command. | |||
# Save the scene. | |||
# Update the green link to use the newly saved version of the scene. | |||
# Test. | |||
===Color The Cylinder Ends=== | |||
Here are the commands to generate a translucent cylinder with a red top and a blue bottom: | |||
<pre>draw cyl1 cylinder diameter 100.0 color gray (thr607) (leu250) nofill mesh | |||
draw cyl1 color translucent -1 | |||
draw ctop circle diameter 100.0 color red (thr607) (leu250) | |||
draw ctop color translucent -1 | |||
draw cbot circle diameter 100.0 color blue (leu250) (thr607) | |||
draw cbot color translucent -1</pre> | |||
*Adding '''nofill mesh''' to the end of the first ''draw'' command removes the ends of the gray cylinder. | |||
*Next, the two ''circle'' commands create new cylinder ends (circles) colored red and blue. | |||
*Notice the reversed order of the amino acids in the two ''circle'' commands. The first position is the center of the circle. The second position makes the circle perpindicular to the axis between the two specified positions. | |||
If you have questions or need more details, don't hesitate to email {{Template:Contact}} | |||
</StructureSection> | </StructureSection> | ||
Eric Martz thanks Ángel Herráez for help with draw commands, and Gulalai Shah for inspiring discussion. | |||
== References == | == References == | ||
<references/> | <references/> | ||
[[Category: BioMolViz]] | |||
[[Category: Macromolecular Assemblies]] |
Latest revision as of 17:33, 8 November 2020
Here are two ways of visualizing the position of a lipid bilayer membrane around the integral membrane portion of a trans-membrane protein.
a hydrophobic vs. polar color scheme, the protein within the lipid bilayer is seen to have a largely hydrophobic surface. The example used here is 5lil. Methods: PseudoatomsAny protein structure that includes an integral membrane domain can have pseudoatoms added to the PDB file by the Orientations of Proteins in Membranes Server[1]. Structures in the Protein Data Bank are pre-processed; other structures may be uploaded to the server. The server removes hydrogen atoms from the model. You can download the pseudoatom-enhanced PDB file from the server, and then upload it to Proteopedia for use in a Molecular Scene. In the case of 5LiL, I had difficulty getting the pseudoatoms to show from a green link (even though they displayed in the SAT). In the PDB file, the pseudoatoms followed a MASTER record. I deleted the MASTER record and all CONECT records, and then I got the scene shown, using File:5lil opm2.pdb.
The pseudoatoms have residue name DUM (for dummy atoms). This can be used in the SAT to select them. The red pseudoatoms are oxygens, and the blue pseudoatoms are nitrogens. Thus, coloring them by chemical element (CPK) will make them red and blue. Methods: Translucent CylinderChoose boundary residuesFirst you must choose residues or atoms in the model that are closest to the ends of the desired cylinder. If you wish, you can use a pseudoatom-enhanced model to assist. Viewing the model in FirstGlance in Jmol, you can use Hydrophobic/Polar (Views tab) to visualize the domain with a hydrophobic surface. Touching or clicking on atoms at the domain boundaries will identify those amino acids. If you wish, you can center a particular residue, then zoom in, and see atomic detail with Vines/Sticks (Views tab). For 5LiL, I chose thr607 and leu250. These will be used in the example commands below, but you must substitute appropriate amino acids for your model. 5LiL is a homo-dimer of chains A and B. Since I did not specify a chain in the commands below, the position of, for example, thr607 will be the average position of all the atoms in thr607 in both chains. That average position will be at the axis of the rod-shaped molecule. Your molecule may require a different strategy. Generate the CylinderTo generate the cylinder in the SAT, click on the button below JSmol Advanced: Open JSmol Console. The command to generate the cylinder will have this form: draw cyl1 cylinder diameter 100.0 color gray (thr607) (leu250) Adapt this command to your situation, copy and paste it into the LOWER box of the Console, and click the Run button at the lower-left of the console box.
draw * delete
draw cyl1 color translucent -1 Once you have the desired cylinder in the SAT, save the scene, insert a green link in the wikitext box, and save the page. Then test the green link. After displaying the cylinder with your new green link, if you then click a different green link, the cylinder may persist, unwanted. To fix this:
Color The Cylinder EndsHere are the commands to generate a translucent cylinder with a red top and a blue bottom: draw cyl1 cylinder diameter 100.0 color gray (thr607) (leu250) nofill mesh draw cyl1 color translucent -1 draw ctop circle diameter 100.0 color red (thr607) (leu250) draw ctop color translucent -1 draw cbot circle diameter 100.0 color blue (leu250) (thr607) draw cbot color translucent -1
If you have questions or need more details, don't hesitate to email
|
|
Eric Martz thanks Ángel Herráez for help with draw commands, and Gulalai Shah for inspiring discussion.
ReferencesReferences
- ↑ 1.0 1.1 Lomize MA, Pogozheva ID, Joo H, Mosberg HI, Lomize AL. OPM database and PPM web server: resources for positioning of proteins in membranes. Nucleic Acids Res. 2012 Jan;40(Database issue):D370-6. doi: 10.1093/nar/gkr703., Epub 2011 Sep 2. PMID:21890895 doi:http://dx.doi.org/10.1093/nar/gkr703