===== Molecular Visualization using VMD (Visual Molecular Dynamics) =====

VMD, a molecular graphics program, is developed in the Theoretical and Computational Biophysics group (Univ. of Illinois) together with other software (eg. NAMD) to simulate molecular systems. The group is lead by Klaus Schulten.

The tutorial will first have some insight into a raw PDB file entry. Then the tutorial will guide you through the topics given in the table of contents.

The structures used in this tutorial are Hen Egg Lysozyme (1HEW) for the basic parts, a Phosphotransferase from E. coli (1WCR), HSP90 (1YER and 2IOR) for alignment and electrostatic calculations/visualizations.

====  Exploring the RCSB Protein Data Bank (PDB) ====
**//Goal: getting familiar with the PDB Database //**

<code>
 * Open a browser and go to: http://www.rcsb.org
 * Search for the PDB ID: 1HEW
 * On the right part of the Structure Explorer page, click on Download_Files 
    -> PDB_File and save 1HEW.pdb as a PDB file (Text)
 * Open this file with a text editor and answer the following questions: 
</code>

Questions: 
  - What protein is this[[courses:hidden:vmd#Exploring the RCSB Protein Data Bank (PDB)|_?_]]
  - Does it contain any ligands? If yes, which[[courses:hidden:vmd#Exploring the RCSB Protein Data Bank (PDB)|_?_]]
  - What technique was used to determine its structure[[courses:hidden:vmd#Exploring the RCSB Protein Data Bank (PDB)|_?_]]
  - What is the resolution of the structure[[courses:hidden:vmd#Exploring the RCSB Protein Data Bank (PDB)|_?_]]
  - Where could we find more details about the structure determination procedure[[courses:hidden:vmd#Exploring the RCSB Protein Data Bank (PDB)|_?_]]




====  Protein and ligand visualization and selection of structure parts (1HEW) ====
**//Goal: identification of the nature of inter-molecular interactions//**

<code>
  * Open the VMD Program (vmd)
  * In VMD_Main window, go to: File 
     -> New_Molecule; click Browse button and browse for downloaded in Task 1 file 1HEW.pdb, 
     then click Load and close the Molecule_File_Browser window
  * In VMD_Main window, go to: Graphics -> Representations...
    * Type 'protein' (instead of 'all') in the Selected_Atoms field
    * Click on: Create_Rep button
    * For the new representation choose Secondary Structure under Coloring_Method and 
      New_Cartoon under Drawing_Method
  * Note that double click on representation list window (with Style, Color, Selection table) 
    switches on and off displaying of representations. Switch off the other represenations and 
    leave only the new one, and answer the following questions: 
</code>
Questions:
  - How many alpha-helices are present in the protein[[courses:hidden:vmd#Protein and ligand visualization (1HEW) and selection of structure parts|_?_]]
  - How many beta-sheets are present in the protein[[courses:hidden:vmd#Protein and ligand visualization (1HEW) and selection of structure parts|_?_]]
  - How many beta-strands form the beta-sheet? Are they parallel or anti-parallel[[courses:hidden:vmd#Protein and ligand visualization (1HEW) and selection of structure parts|_?_]]

====  Create Graphics representations ====
//**Explore the capabilities to visualize structural parts**//
<code>
  * In the Drawing_Method field, change from New_Cartoon to Surf
  * In the Coloring_Method field, change from Secondary Structure 
    to Backbone (protein backbone will be rendered in green and the 
    side chains in blue)
  * Click on Create_Rep button and replace 'protein' by 'resname NAG' in the Selected_Atoms field 
  * Change the Drawing_Method from Surf to Licorice and the Coloring_Method 
    from Backbone to Name, and answer the following question:
</code>

  - Does the inhibitor interact mostly with the protein backbone or side chains[[courses:hidden:vmd#Create Graphics representations|_?_]]

<code>
   * Click on representation number 2 (Surface-Backbone-protein) and change the Coloring Method 
     from Backbone to Res_Type (apolar residue are represented in white, neutral polar in green, 
     positively charged in blue, and negatively charged in red). It will give you an idea of 
     the chemical nature of the enzyme's active site.
</code>

====  Inspecting Protein Structure: The Ramachandran Plot ====
//**Goal: analysis of protein structure via the Ramachandran Plot **//

<code>
  *  In VMD_Main window, go to Extensions -> Analysis -> Ramachandran Plot
  *  In the Molecule menu select 1HEW.pdb
  *  Note that several points (on the right hand side of the graph) fall outside of the allowed conformational region
  *  Click on these points (yellow squares) to identify these residues and answer the following question:
</code>
  - Explain why these residues can adopt the so-called forbidden conformations[[courses:hidden:vmd#Inspecting Protein Structure: The Ramachandran Plot|_?_]]
<code>
  * Close the VMD program (File -> Quit)
</code>





====  Coloring Chains and Hydrogen patterns ====
//**Goal: identification of the nature of inter-molecular interactions **//

<code>
  * Download, save a molecule with PDB ID 1WCR (as in Task 1)
  * Open VMD and load this molecule - you should see the phosphotransferase enzyme IIA
    from Escherichia coli in the VMD graphical window
  * Color it by chain
    * Create a new represenation and 
      * change the Drawing_Method to HBonds 
      * and the Coloring_Method to ColorID . 
      * Set ColorID to 4 (yellow)
    * Relax the hydrogen-bond parameters a bit (distance cutoff: 3.5, and angle cutoff: 25)
    * Increase Line_Thickness to 6
    * Manipulate the model and answer the following questions: 
</code>
  - What type of hydrogen-bond pattern do you observe, i.e., intra- or interchain[[courses:hidden:vmd#Coloring Chains and Hydrogen patterns|_?_]]
  - Based on the hydrogen-bond pattern, what kind of interactions (hydrophilic or hydrophobic) do you think is mainly responsible for keeping the protein structure as a trimer (3 units)[[courses:hidden:vmd#Coloring Chains and Hydrogen patterns|_?_]]

==== Alignment of proteins ====
//**Goal: Compare structures from different organisms or different folds**//
<code>
    * Sometimes similar proteins exist in different species - they may perform the same function,
      but have different sequences. Therefore, it is useful to compare structures of these proteins. 
      To do that, one needs to align the proteins based on their sequence 
      and sometimes also taking the structure into account.
    * Download from RCSB the pdb entries **2IOR** and **1YER** and store them in a folder on the desktop.
</code>
   - How do these protein differ[[courses:hidden:vmd#Alignment of proteins|_?_]] Look for answers on the PDB webpages and/or in the downloaded files.
<code>
   * Load both structures into vmd
</code>
  - Can you see both proteins[[courses:hidden:vmd#Alignment of proteins|_?_]] Search for them by zooming and/or rotating (translation mode can be switched on by pressing 'T' and switched off by pressing 'R' on the keyboard). 
<code>
   * To align the proteins Go to the menu: Extensions->Analysis->Multiseq
   * You need to specify a folder for temporary files. Create an empty folder on the desktop and choose it
   * Next, in the window that appears, select Tools->Stamp Structural Alignment and press Ok.
   * If the center of your protein is shifted now, Select Display -> Reset view
   * Please note, while closing the panel of the Multiseq Analysis, the original representation is started again.
     Go to Display -> Representations and enable the Cartoon Representation again.
   * You can color the structures differently using the Graphical Representations to clearly see the differences.
   * Create an additional representation for each protein that selects the following atoms: 
     not (protein or water). Display this in Licorice represenation.
</code>
  - What region displays largest difference between the two structures[[courses:hidden:vmd#Alignment of proteins|_?_]] 
<code>
  * In the VMD Main window select the molecule of 1yer and right click on it to save its coordinates.
    Choose only the protein part (chain A) in 'Selected Atoms' and save it in a file called 1YER_aligned.pdb and repeat these steps also for the other molecule 2ior_aligned.pdb
  * Close the VMD program 
</code>

====  Visualization of electrostatic potentials ====
//**Goal: Visualize possible interaction sites on proteins**//

<code>
  * Start VMD again and use the aligned Proteins from the previous section (load the previously saved coordinates). 
    The proteins should show up as aligned entities.
  * Go to the PDB2PQR webserver: http://nbcr-222.ucsd.edu/pdb2pqr_2.1.1/
  * Upload one of the aligned coordinate files. 
    Choose the force field (AMBER) and naming scheme (AMBER), and click ’Submit’.
  * Save the output files (xxxx_aligned.in, xxxx_aligned.propka and xxxx_aligned.pqr) on disk.
  * Press "Run APBS" -> Launch
  * Again, save the output files, espcially the *.dx.gz grid files. 
    You might want to look at the program output and error logs if there are errors or warnings.
  * Repeat for both structures and unpack the *.dx.gz file.
</code>
In the respective pqr files, what are the net charges for the two proteins[[courses:hidden:vmd#Alignment of proteins|_?_]] 
<code>
  * To visualize the electrostatic potential, select in the vmd main window the molecule 1yer.
  * Go to File-> Load data in molecule. Select the corresponding dx grid file from your apbs output Press Load.
  * In Graphics Representation, change the standard representation to New Cartoon.
  * Create a new representation, Color by ID (blue), Select drawing method Isosurface. Select an Isovalue of 3 
    and draw as Solid surface.
  * Again Create a new representation. Color by ID (red), Select again drawing method Isosurface with an Isovalue of -3.
    also draw as Solid surface.
  * Repeat the steps with the other molecule. 
</code>
What are the differences between the electrostatic potentials of the two proteins? From which species are the respective proteins from[[courses:hidden:vmd#Alignment of proteins|_?_]] 
==== Tutor only ====

{{courses:hidden:visualizationstatesbasicvmdtutorial.zip|Visualization states from basic tutorial}}

{{:courses:hidden:hsp90_aligned_electrostatics.zip|Visualization state alignment and electrostatics}}