Jmol/Visualizing membrane position
From Proteopedia
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 Image: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.
References
- ↑ 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