|
|
Line 1: |
Line 1: |
− | =Introduction=
| |
− | AMBER is a molecular dynamics program that can be run on your protein/ligand complex to ensure that the interactions between the two structures are stable. DOCK shows us how the two interact with each other at one point in time. AMBER looks at those interactions over time to ensure that forces will not occur which will push the ligand out of the binding site as the complex naturally moves. This tutorial will again be working with PDB #4s0v
| |
| | | |
− | ==Setting Up Your Environment==
| |
− | Just as with DOCK you should set up for directory structure at this point to keep everything organized and easy to find. We will be creating a new structure which looks like:
| |
− |
| |
− | =Structure=
| |
− | Before starting the analysis it's best to download a new protein/ligand complex from the PDB and isolate both the protein and ligand structures. Follow the steps in the *[http://ringo.ams.stonybrook.edu/index.php/2023_DOCK_tutorial_1_with_PDBID_4S0V] tutorial to do this. The inputs we need are the isolated protein with NO hydrogens and NO charges; and the ligand with hydrogens and charges. In other words, once you isolate the protein structure in Chimera, save it with a filename such as, 4s0v_protein_for_AMBER.pdb. Then isolate the ligand structure, add hydrogens and re-do whatever protonation changes you made in the *[http://ringo.ams.stonybrook.edu/index.php/2023_DOCK_tutorial_1_with_PDBID_4S0V] tutorial. Once the protonation state is correct, add charges and save the file as 4s0v_ligand_for_AMBER.mol2.
| |
− |
| |
− | Once these two files have been generated, scp them over to the 001.structure directory on Seawulf.
| |
− |
| |
− | =Force Field Parameters=
| |
− | AMBER needs force field parameters to run correctly. This only needs to be done for the ligand since information about the protein has been done over previous years. To generate the necessary files for AMBER we will be working on the command line on Seawulf, please cd into your 002.forceFieldParameters directory.
| |
− |
| |
− | To generate the parameters for the ligand we need to run the following slurm script:
| |
− | #!/bin/bash
| |
− | #
| |
− | #SBATCH --job-name=4s0v_AMBER_parameters
| |
− | #SBATCH --output=parameters_output.txt
| |
− | #SBATCH --ntasks-per-node=24
| |
− | #SBATCH --nodes=6
| |
− | #SBATCH --time=48:00:00
| |
− | #SBATCH -p long-24core
| |
− |
| |
− | module load amber/16
| |
− |
| |
− | antechamber -i ../001.structure/4s0v_ligand_for_AMBER.pdb -fi pdb -o 4s0v_ligand_antechamber.mol2 -fo mol2 -at gaff2 -c bcc -rn LIG -nc 0
| |
− |
| |
− | The important parameter in the above command is the -nc option at the end. This is telling antechamber what the total charge is on your ligand. This number needs to be the same as the number used in the previous step when you added charges to the structure in Chimera.
| |
− |
| |
− | Once this has completed running you will see multiple new files in your directory including, 4s0v_ligand_antechamber.mol2. This is the file with the parameters we just generated and what we need to use to generate the next file we need, a frcmod file, which will contain modified force field parameters. To generate the frcmod file, run the command:
| |
− | parmchk2 -i 4s0v_ligand_antechamber.mol2 -f mol2 -o 4s0v_ligand.am1bcc.frcmod
| |
− |
| |
− | If you get an error when running this line, try typing:
| |
− | module load amber/16
| |
− |
| |
− | and then running the command again. Once it's done running you will see the 4s0v_ligand.am1bcc.frcmod file in your directory.
| |
− |
| |
− | =TLeap Implemenation=
| |
− | Remembering that the point of running AMBER on our complex is to investigate the dynamics of the entire complex over time. In this next step we bring together the separate ligand and protein files so AMBER can be run on the whole comlex. This is done with a program called tleap. Again we will be working on the command line in Seawulf, please cd into your 003.tleap directory.
| |
− |
| |
− | To create the input file:
| |
− | vi 4s0v_tleap.in
| |
− |
| |
− | and type the following lines into the input file:
| |
− | #!/usr/bin/sh
| |
− | ###load protein force field
| |
− | source leaprc.protein.ff14SB
| |
− |
| |
− | ###load GAFF force field (for our ligand)
| |
− | source leaprc.gaff
| |
− |
| |
− | ###load TIP3P (water) force field
| |
− | source leaprc.water.tip3p
| |
− |
| |
− | ###load ions frcmod for the tip3p model
| |
− | loadamberparams frcmod.ionsjc_tip3p
| |
− |
| |
− | ###needed so we can use igb=8 model
| |
− | set default PBradii mbondi3
| |
− |
| |
− | ###load protein pdb file
| |
− | rec=loadpdb ../001.structure/4s0v_protein_for_AMBER.pdb
| |
− |
| |
− | #bond rec.649.SG rec.762.SG
| |
− | #bond rec.454.SG rec.472.SG
| |
− | #bond rec.444.SG rec.447.SG
| |
− | #bond rec.328.SG rec.339.SG
| |
− | #bond rec.385.SG rec.394.SG
| |
− | ###load ligand frcmod/mol2
| |
− | loadamberparams ../002.forceFieldParameters/4s0v_ligand.am1bcc.frcmod
| |
− | lig=loadmol2 ../002.forceFieldParameters/4s0v_ligand_antechamber.mol2
| |
− |
| |
− | ###create gase-phase complex
| |
− | gascomplex= combine {rec lig}
| |
− |
| |
− | ###write gas-phase pdb
| |
− | savepdb gascomplex 4s0v.gas.complex.pdb
| |
− |
| |
− | ###write gase-phase toplogy and coord files for MMGBSA calc
| |
− | saveamberparm gascomplex 4s0v.complex.parm7 4s0v.gas.complex.rst7
| |
− | saveamberparm rec 4s0v.gas.receptor.parm7 4s0v.gas.receptor.rst7
| |
− | saveamberparm lig 4s0v.gas.ligand.parm7 4s0v.gas.ligand.rst7
| |
− |
| |
− | ###create solvated complex (albeit redundant)
| |
− | solvcomplex= combine {rec lig}
| |
− |
| |
− | ###solvate the system
| |
− | solvateoct solvcomplex TIP3PBOX 12.0
| |
− |
| |
− | ###Neutralize system
| |
− | addions solvcomplex Cl- 0
| |
− | addions solvcomplex Na+ 0
| |
− |
| |
− | #write solvated pdb file
| |
− | savepdb solvcomplex 4s0v.wet.complex.pdb
| |
− |
| |
− | ###check the system
| |
− | charge solvcomplex
| |
− | check solvcomplex
| |
− |
| |
− | ###write solvated toplogy and coordinate file
| |
− | saveamberparm solvcomplex 4s0v.wet.complex.parm7 4s0v.wet.complex.rst7
| |
− | quit
| |
− |
| |
− | We should again run this input file using a slurm script:
| |
− | #!/bin/bash
| |
− | #
| |
− | #SBATCH --job-name=4s0v_tleap
| |
− | #SBATCH --output=tleap_output.txt
| |
− | #SBATCH --ntasks-per-node=24
| |
− | #SBATCH --nodes=6
| |
− | #SBATCH --time=48:00:00
| |
− | #SBATCH -p long-24core
| |
− |
| |
− | module load amber/16
| |
− |
| |
− | tleap -f 4s0v_tleap.in
| |
− |
| |
− | =Complex Minimization and Equilibration=
| |