Difference between revisions of "Gromacs/C3/Production-Run-for-Lysozyme/English"
Snehalathak (Talk | contribs) |
|||
| Line 14: | Line 14: | ||
'''Learning Objectives''' | '''Learning Objectives''' | ||
|| In this tutorial, we will, | || In this tutorial, we will, | ||
| − | * Complete | + | * Complete 1 nanosecond production run for '''lysozyme''' |
* About timescales of motions | * About timescales of motions | ||
* Time taken for simulation | * Time taken for simulation | ||
Revision as of 12:57, 26 May 2022
| Visual Cue | Narration |
| Slide Number 1
Title Slide |
Welcome to the spoken tutorial on Production Run for Lysozyme. |
| Slide Number 2
Learning Objectives |
In this tutorial, we will,
|
| Slide Number 3
Learning Objectives |
* Save data to various output files
|
| Slide Number 4
System and Software Requirement |
To record this tutorial, I am using
|
| Slide Number 5
Pre-requisites |
To follow this tutorial,
|
| Slide Number 6
Code Files |
* Files used in this tutorial are provided in the code files link.
|
| Show input files.
Show, copied, files in firsmd directory. |
Several files are provided with this tutorial.
If needed, copy the required files from the previous step. |
| Open Ubuntu 20.04 LTS app (Windows) or terminal (Linux).
Type cd Documents/firstmd and press Enter. |
Open a terminal and go to the working directory. |
| Type cd ls and press Enter. | I will also list the files. |
| Cursor on npt.gro and npt.tpr. | The structure and configuration files for MD simulation are npt.gro & npt.tpr.
At this point, the following has happened to the starting PDB file. |
| Show all the files in the folder. | The protein is solvated, charge neutralized with ions and energy minimized.
Then, in the initial MD step, temperature and pressure was equilibrated. This prepared the system for the Production run MD. I have all the files from all the steps in this folder. |
| Cursor on md.mdp and open in a text editor. | Open the mdp file provided for MD simulation in a text editor. |
| Cursor on the time. | Here, the timescale of simulation is 1 ns. |
| This allows us to follow changes in the conformation of the protein.
In energy minimization and equilibration steps, the time duration was shorter. | |
| Back to the text editor and close the mdp file. | Often, step size and length of simulation are adjusted.
This file specifies ½ million steps in 1 ns. Frequency with which structures are written to the trajectory is also adjustable. Let’s close the mdp file. |
| Slide Number 7
Timescale of Protein Motions |
This image shows the timescale of molecular motions we usually encounter.
The ns-ps time encompasses the lower end of protein motions for classical MD. |
| Slide Number 8
Timescale of Protein Motions |
* The microsecond or millisecond, ligand binding or folding is studied.
|
| Go to the terminal. | Go to the terminal. |
| Type, cd Documents/firstmd and press Enter. | Change to the firstmd directory, where the files are saved. |
| Type, gmx grompp -f md.mdp -c npt.gro -t npt.cpt -p topol.top -o md-1.tpr and press Enter. | Enter the command as seen here to assemble the configuration for the md run.
npt.gro is the starting structure for the production md run. |
| Cursor on md-1.tpr. | The tpr file is the output file from this step.
An -r flag is used if we need to restrain the protein conformation. This is often done when studying protein-ligand binding. Here we will not used it. |
| Cursor on the message. | While creating the file, few messages are seen on the terminal. |
| Type ls and press Enter. | Enter ls on the terminal to list the files.
Notice the output md-1.tpr, that is created. This prepared the system for the production md. |
| Press Ctrl+L to clear the screen. | As we did for the energy minimization step, we will enter the next command. |
| Type,
gmx mdrun -v -deffnm md-1 and press Enter. |
Enter the mdrun command as seen to start the production run.
Usually, molecular dynamics refers to equilibration followed by production run. |
| Show the output on the terminal. | This time, the verbose flag shows only a few details.
The number of steps and time needed to complete are seen on the screen. |
| Show the time taken for the process. | Process will stop when it reaches 500 thousand steps in 1 ns.
This process may take a day or so on many personal computers. |
| Press Ctrl+C. | You may run this process later on.
Use the provided data for further analysis. To abort the process, press Control and C keys together once. In a few minutes Gromacs stops the process. |
| Show files in the file manager. | Several files with md prefix are generated in the working directory. |
| Highlight and delete the files to be deleted. | If you had aborted the process, delete the incomplete files that are created. |
| If you have completed the process, do not delete the files.
A gro file will get created at the completion of the process. | |
| Show files to be copied from provided file.
Copy them to the working directory. |
Copy the provided files with the md-1 prefix to the working directory. |
| Type ls -ltr and press Enter. | This step creates additional file types with cpt and xtc extensions. |
| Cursor on md-1.xtc file. | The xtc file is a compressed trajectory file.
|
| Cursor on md.trr, md.log and md.edr files. | Trr, log and edr files are also created in this process. |
| Show the output files. | In case of errors, learner may refer to the log file for trouble shooting. |
| We have generated several files for further analysis.
Data analysis path varies depending on our aim or objective of simulation. | |
| Slide Number 9
Data Analysis |
Few examples of data analysis are,
|
| Open a web browser.
Go to, https://manual.gromacs.org/documentation/2018/user-guide/cmdline.html . |
Open a web browser.
Let’s go to the Gromacs manual site as seen here. |
| Cursor on the left panel. | On the left frame notice the various categories of topics for analysis.
Atomic interactions, kinetics and PCA are also options. These are useful to analyse flexible structures or ones such as molten globules. |
| Learner is encouraged to read further on this to know more. | |
| Go to the terminal.
Press Ctrl+L. |
Go to the terminal.
I will clear the terminal for clarity in the video. |
| Type, vmd md-1.gro md-1.xtc and press Enter. | Let’s open VMD to view the trajectory from the 1ns long simulation.
This time we will use the xtc file instead of the trr file to load the trajectory. Xtc is a compressed format of the simulation trajectory. Windows users must first open VMD. Then open the gro file, followed by the xtc trajectory file. |
| Click on Graphics, representation.
|
Hiding the water molecules from graphics display will be useful for clarity. |
| You may adjust the graphic representation of the protein to your desire.
I will change it to show as ribbons. | |
| In the VMD terminal enter the command pbc box. | Display the box, using the command as seen here in the vmd terminal. |
| Play the trajectory. | Pause the video and play the trajectory.
Notice the xtc file has 1002 frames. |
| Later on, you may load the trr file and notice, it has only 101 frames.
Xtc format has precision of 3 decimal places compared to 6 decimals for trr. Hence, even though xtc has 10 times more frames, the file sizes are comparable. | |
| Stop playing the trajectory.
Close the Window. |
Visualization alone is not sufficient to know the features of the data.
Hence let’s do few data analysis to see the features embedded. |
| Click on Extensions, Analysis and choose RMSD Trajectory Tool. | Click on Extensions, Analysis and choose RMSD Trajectory Tool.
A RMSD trajectory tool window opens. |
| Let’s calculate RMSD of the backbone atoms in the trajectory. | |
| In Selection modifiers, check the box for Backbone. | In Selection modifiers, check the box for Backbone.
Notice that, Trace or NH can also be selected for plotting. |
| Cursor over trajectory section. | From the Trajectory section, we can also select a range of frames for analysis.
I will retain the defaults and check the boxes for Plot and Save. |
| Cursor on trajrmsd.dat. | Notice that, the data will also get saved with the given filename.
You may change it if you desire. I will leave it as is. |
| Click on md-1.gro in mol section. | In the molecule section, click on md-1.gro file to select it. |
| Click on the RMSD.
Show the graph. |
Click on the RMSD button and the graph is seen.
X axis is the frame number and y axis is the RMSD in the graph. |
| Click on File in graphics window. | We can also export the data in many different formats as seen here. |
| Screenshot of the file manager and output file. | The data can also be written to a text file in two column format.
Learner may use a plotting program of your choice to plot it. |
| Close the RMSD Trajectory Tool window. | Let’s close this window. |
| Cursor on the graphics window. | Now, let’s display frame 50 and align it to md-1.gro. |
| Close the RMSD graphical window.
Go to representations and click on Create Rep. |
Let’s close this window.
Go to representations window and create a new representation. |
| Go to the Trajectory tab. | Go to the Trajectory tab. |
| Cursor on Draw multiple frames form, type, 50 and press Enter. | In the Draw multiple frames form, type 50 and press Enter.
This shows the 50th frame also. |
| Open the RMSD Trajectory Tool. | To align them, go back to the RMSD Trajectory Tool. |
| For Reference mol, choose Selected. | For Reference mol, choose Selected. |
| Type 50 for Reference frame. | Type 50 for Reference frame. |
| Uncheck the the box for Plot. | I will uncheck the box for Plot. |
| Click on md-1.gro in mol section. | In the molecule section, click on md-1.gro file to select it. |
| Click on ALIGN.
Show alignment in the graphical window. |
Now, click on ALIGN and notice the two aligned molecules. |
| Click on FIle, Save Visualization State.
In the filename form, type filename analysis.vmd and click on Ok. |
Click on FIle, Save Visualization State to save the state in VMD.
Enter file name analysis.vmd and retain the working directory to save the file. |
| Show the Options menu. | Pause the video and explore the Options menu.
Various file formats and options are available for user convenience. |
| Screenshot of energy command options to export to xvg file. | You may further, extract and plot more parameters. |
| Slide Number 9
Summary |
Now let’s summarize. In this tutorial, we,
|
| Slide Number 10
Summary |
* Loaded the xtc trajectory file
|
| Slide Number 11
Assignment |
For assignment activity, please do the following.
|
| Slide Number 12
Assignment |
Using the energy command,
|
| Slide Number 13
Spoken Tutorial Project |
This video summarises the Spoken Tutorial Project.
Please download and watch it. |
| Slide Number 14
Spoken Tutorial workshops |
We conduct workshops using spoken tutorials and give certificates.
Please write to us. |
| Slide Number 15
Forum for questions |
Post your timed queries in this forum. |
| Slide Number 16
Acknowledgment |
Spoken Tutorial Project is funded by MoE, Government of India. |
| This is Rani from IIT, Bombay. Thank you for joining. |
