Introduction to Molecular Dynamics with GROMACS

does, then a list of (optional) input/output files, and finally a list of available flags. The same informationis present with “man pdb2gmx” (if it was properly installed), and you can also use thewww.gromacs.org website: Choose “documentation” in the top menu, and then “manuals” in the leftsidemenu. There is both an extensive book describing everything and a shorter online reference to thecommands and file formats.For the first command we’ll help you. pdb2gmx should read its input from a coordinate file, and justfor fun we’ll pick a specific water model and force field. Typepdb2gmx -f 1LYD.pdb -water tip3pThe program will ask you for a force field - for this tutorial I suggest that you use “OPLS-AA/L”, andthen write out lots of information. The important thing to look for is if there were any errors or warnings(it will say at the end). Otherwise you are probably fine. Read the help text of pdb2gmx (-h) andtry to find out what the different files it just created do, and what options are available for pdb2gmx!The important output files are (find out yourself what they do)conf.gro topol.top posre.itpIf you issue the command several times, Gromacs will back up old files so their names start with ahash mark (#) - you can safely ignore (or erase) these for now.Adding solvent water around the proteinIf you started a simulation now it would be in vacuo. This was actually the norm 25-30 years ago,since it is much faster, but since we want this simulation to be accurate you are going to add solvent.However, before adding solvent we have to determine how big the simulation box should be, andwhat shape to use for it. There should be at least half a cutoff between the molecule and the box side,and preferably more. Due to the limited time you be cheap here and only use 0.5nm margin betweenthe protein and box. To further reduce the volume of the box you could explore the box type optionsand for instance test a “rhombic dodecahedron” box instead of a cubic one. The program you should useto alter the box is “editconf”. This time you will have to figure out the exact command(s) yourself,but the primary options to look at in the documentation are “–f”, “–d”, “–o”, and “–bt”.As an exercise, you could try a couple of different box sizes (the volume is written on the output!)and also other box shapes like “cubic” or “octahedron”. What effect does it have on the volume?The last step before the simulation is to add water in the box to solvate the protein. This is done byusing a small pre-equilibrated system of water coordinates that is repeated over the box, and overlappingwater molecules removed. The Lysozyme system will require roughly 6000 water molecules with 0.5nmto the box side, which increases the number of atoms significantly (from 2900 to over 20,000).GROMACS does not use a special pre-equilibrated system for TIP3P water since water coordinates canbe used with any model – the actual parameters are stored in the topology and force field. For this reasonthe program only provides pre-equilibrated coordinates for a small 216-water system using a watermodel called “SPC” (simple point charge). The best program to add solvent water is “genbox”, i.e.you will generate a box of solvent. You can select the SPC water coordinates to add with the flag “–csspc216.gro”, and for the rest of the choices you can use the documentation. Remember that yourtopology needs to be updated with the new water molecules! The easiest way to do that is to supply thetopology to the genbox command.Before you continue it is once more a good idea to look at this system in PyMOL. However, Py-MOL and most other programs cannot read the Gromacs gro file directly. There are two simple ways tocreate a PDB file for PyMOL: