2013 AMBER Tutorial with UMP and OMP
For additional Rizzo Lab tutorials see AMBER Tutorials.
In this tutorial, we will learn how to run a molecular dynamics simulation of a protein-ligand complex. We will then post-process that simulation by calculating structural fluctuations (with RMSD) and free energies of binding (MM-GBSA).
- 1 I. Introduction
- 2 II. Structural Preparation
- 3 III. Simulation using sander
- 4 IV. Simulation Analysis
- 5 V. Frequently Encountered Problems
Amber - Assisted Model Building with Energy Refinement - is a suite of about multiple programs for perform macromolecular simulations. Amber11, the current version of Amber, includes newly released functionality such as PMEMD, particle mesh Ewald MD and soft-core Thermodynamics Integration MD. For the tutorial, we are using the newest version AMBER12.
The Amber 12 Manual is the primary resource to get started with Amber12. (Tip: Using Adobe Acrobat to view the file, you can simply search the document for keywords such as the name of a simulation parameter, which saves much time.) In addition, Amber Tools User's Manual serves as another reference while using Amber tools.
Here are some programs in Amber
- LEaP: an preparing program for constructing new or modified systems in Amber. It consists of the functions of prep, link, edit, and parm for earlier version of Amber.
- ANTECHAMBER: in additional to LEap, this main Antechamber suite program is for preparing input files other than standard nucleic acids and proteins.
- SANDER: according to the Amber 12 manual, it is 'a basic energy minimizer and molecular dynamics program' that can be used to minimize, equilibrate and sample molecular conformations. And this is the program we mainly use in this tutorial to generate trajectory files of the molecular system.
- PMEMD: version of SANDER that has improved parallel scaling property and optimized speed.
- PTRAJ: an analysis program for processing trajectory files. One can use ptraj to rotate, translate the structures, evaluate geometrical features and so on.
There is a mailing list you could sign-up for, as an additional resource.
UMP and OMP
For information of the UMP-OMP system, see 2013 DOCK tutorial with Orotodine Monophosphate Decarboxylase.
While performing MD simulations, it is convenient to adopt a standard directory structure / naming scheme, so that files are easy to find / identify. For this tutorial, we will use something similar to the following:
~username/AMS536/AMBER-Tutorial/001.CHIMERA.MOL.PREP/ 002.ANTE.TLEAP/ 003.SANDER/ 004.PTRAJ/ 005.MMGBSA/
II. Structural Preparation
Preparation in Chimera
In this AMBER tutorial, we will use the same system with previous DOCK part. To begin with, we need three files under directory 001.CHIMERA.MOL.PREP.
Downloading the PDB file (1LOQ)
Since we need to edit the PDB before we use it in Chimera we should do this manually. Go to PDB homepage (http://www.rcsb.org/pdb/home/home.do ) enter the protein ID (1LOQ) in the search bar, click Download Files in the top-right of the webpage, then select PDB File (text). In the new window, save the file in Downloads.
Preparing the ligand and receptor in Chimera
In this section, we will create three new files and save them in the 001.CHIMERA.MOL.PREP/ folder:
1LOQ.dockprep.mol2 (complete system with hydrogens and charges) 1LOQ.receptor.noH.mol2 (the receptor alone, without hydrogens) 1LOQ.ligand.mol2 (the ligand alone)
To prepare these files, first copy the original PDB file into the 00.files/ folder and open it with VIM ($ vim 1LOQ.pdb). Because the residue name of the ligand (U) will give us some problems when assigning charges, change the residue name "U" to "LIG" starting at line 2082. Here is an example command that will change all instances of " U" into "LIG", while preserving the correct spacing:
For this command, g is short for global and c is short for check with the user before making the change.
Next, open up the PDB file (1LOQ.pdb) in Chimera. To delete water molecules and other ligands, click Tools -> Structure Editing -> Dock Prep. Check all boxes and click Okay to the end. Alternatively, the waters can be deleted manually by choosing Select -> Residue -> HOH, then go to Actions -> Atoms/Bonds -> Delete. Hydrogen atoms can be added manually by choosing Tools -> Structure Editing -> Add H.
Next, to add charges to the ligand and receptor, go to Select -> Residue -> LIG, then go to Tools -> Structure Editing -> Add Charge. Choose AMBER ff99SB as the charge model, click Okay, and when prompted chose AM1-BCC charges for the ligand, and make sure the Net Charge is set to -1. (You must consider the chemistry of the ligand when assigning a charge state).
Finally, save this file as 1LOQ.dockprep.mol2.
Generating the final files
To create the receptor file with no hydrogen atoms: Open 1LOQ.dockprep.mol2, click Select -> Chemistry -> Element -> H, then chose Actions -> Atoms/Bonds -> Delete. Save the file as 1LOQ.receptor.noH.mol2.
To create the ligand file: Open 1LOQ.dockprep.mol2, click Select -> Residue -> LIG, then click Select -> Invert, then chose Actions -> Atoms/Bonds -> Delete. Save the file as 1LOQ.ligand.mol2.
An antechamber input file for the next step requires all the atom names to be unique and it only uses the first 3 characters as the name. So if we use 1LOQ.lig.chimera.mol2 as the input file, it will cause errors ("H102" and "H103" will have the same name "H10") We can use Vim to rename the hydrogens, or any other atom with 4 character name. Make these edits and save as 1LOQ.ligand.mol2
An antechamber input file requires all the atom names to be unique. So if we use 1LOQ.ligand.mol2 as the input file, it will cause errors. The program can only recognize atom names of 3 characters ( In this case, H5' and H5 cannot be distinguished from each other. )
22 H5' 40.0697 36.0506 37.6716 H 1 LIG210 0.0761 23 H5 40.6060 37.0416 36.2883 H 1 LIG210 0.0349 24 H4' 38.2510 37.5673 38.1082 H 1 LIG210 0.0967 25 H3' 38.2723 38.5564 35.5613 H 1 LIG210 0.1052 26 H2' 40.3056 39.7322 35.6920 H 1 LIG210 0.1188 27 H1' 39.7587 40.5927 38.4795 H 1 LIG210 0.0753 28 H6 41.3847 41.7815 39.3543 H 1 LIG210 0.0322 29 H5 43.6843 42.3632 39.6188 H 1 LIG210 0.1693 30 H3 44.6515 39.3570 36.8437 H 1 LIG210 0.3481 31 HO3' 36.5607 38.9743 36.9843 H 1 LIG210 0.4373 32 HP3 39.3989 32.7036 35.9199 H 1 LIG210 0.4245 33 HO2' 39.5501 41.7087 35.4449 H 1 LIG210 0.4228
We need to manually rename the atoms.
22 H1 40.0697 36.0506 37.6716 H 1 LIG210 0.0761 23 H2 40.6060 37.0416 36.2883 H 1 LIG210 0.0349 24 H3 38.2510 37.5673 38.1082 H 1 LIG210 0.0967 25 H4 38.2723 38.5564 35.5613 H 1 LIG210 0.1052 26 H5 40.3056 39.7322 35.6920 H 1 LIG210 0.1188 27 H6 39.7587 40.5927 38.4795 H 1 LIG210 0.0753 28 H7 41.3847 41.7815 39.3543 H 1 LIG210 0.0322 29 H8 43.6843 42.3632 39.6188 H 1 LIG210 0.1693 30 H9 44.6515 39.3570 36.8437 H 1 LIG210 0.3481 31 HO1 36.5607 38.9743 36.9843 H 1 LIG210 0.4373 32 HP3 39.3989 32.7036 35.9199 H 1 LIG210 0.4245 33 HO2 39.5501 41.7087 35.4449 H 1 LIG210 0.4228
To begin with, go to 002.ANTE.TLEAP directory. To make sure we have access to the three programs that we want to run (antechamber, parmchk and tleap) and we are using the correct version of amber, we can use the which command, type:
which antechamber which parmchk which tleap
Your results should be similar to this:
/home/wjallen/AMS536/local/amber12/bin/antechamber /home/wjallen/AMS536/local/amber12/bin/parmchk /home/wjallen/AMS536/local/amber12/bin/tleap
Copy parameters of ions to your working directory from the following resource:
cp -r ~lingling/AMS536/AMBER_Tutorial/002.ANTE.TLEAP/rizzo_amber7.ionparms
Then we use antechamber to convert our input mol2 file into files ready for LEaP.Type command:
antechamber -i ../001.CHIMERA.MOL.PREP/1LOQ.lig.mol2 -fi mol2 -o 1LOQ.lig.ante.pdb -fo pdb
Here, -i input file name; -fi input file format; -o output file name; -fo output file format. You will have an output file:1LOQ.lig.ante.pdb
Similarly, we can use antechamber to change the fomat of 1LOQ.lig.mol2 file to prep file:
antechamber -i ../001.CHIMERA.MOL.PREP/1LOQ.lig.mol2 -fi mol2 -o 1LOQ.lig.ante.prep -fo prepi
You will get a set of output files:
ANTECHAMBER_AC.AC ANTECHAMBER_AC.AC0 ANTECHAMBER_BOND_TYPE.AC ANTECHAMBER_BOND_TYPE.AC0 ANTECHAMBER_PREP.AC ANTECHAMBER_PREP.AC0 ATOMTYPE.INF NEWPDB.PDB PREP.INF 1LOQ.lig.ante.prep