![]() Next, load the coordinates file dna.rst7 and it should be loaded like an AMBER7 restart file. ![]() Load the file and in the VMD main windows it should display the number of atoms. It will recognize the file as an AMBER7 parm file. Launch VMD and load the topology file dna.prmtop. We can check our topology and coordinates file with VMD. We now have the topology file dna.prmtop that contains all the parameter definitions for each atom in our system and a coordinate file dna.rst7 that has all the xyz coordinates for each atom. This will generate a topology and coordinates file using the force field and force field modifications we loaded previously. Next we generate the topology file and the coordinates file of our system: > saveamberparm dna dna.prmtop dna.rst7 Around 5336 water molecules should be added. The tleap program will add at least 9.0 angstroms of water molecules between the solute (our DNA molecule) and the borders of the box. ![]() The command solaveoct will create a truncated octahedron periodic cell centered on the variable dna which contains our systems and add explicit water molecules using the OPC water model. Next, we solvate the system: > solvateoct dna OPCBOX 9.0 The ions are added close to the negative sites of the molecule. This command will add ions to the variabla dna, using Na+ type ions until it reaches a charge of ‘0’. At this step we can proceed and add counter ions to neutralize the charge. The structure has a -22 charge which is expected since we have 22 phosphate groups. WARNING: The unperturbed charge of the unit: -22.000000 is not zero. We now load our DNA structure with the command loadpdb and asign it to the variable dna: > dna = loadpdb A-dna.pdbĬheck our structure for possible errors: > check dna This will load the water model known as Optimum Point Charge ( OPC) and the ion parameters that have been shown to be compatible with that current water model. The current AMBER instruction are designed to load the water model and the ions model recommended: > source Since the molecule we are using has a negative charge, we load the ions force field so we can use them to net neutralize our system. In tleap, if we type the comand list, we can now see all the residues that we have loaded: > list Which means that we are loading the ff99bsc0 modifications for the alpha/gamma dihedrals plus the OL1 set of modifications for the epsilon/zeta dihedral and the OL4 modifications for the chi angel. For nucleic acids, it loads: OL15 force field for DNA (99bsc0-betaOL1-eps-zetaOL1-chiOL4) see Source the force field, which has the name 15 $ source 15Īs you can see, when the force field loads, there is a line which states what modifications are being loaded. Open tleap (we do not need to use the visual interface, so tleap is the terminal version): $ tleap A comprehensive review of this force field is available in these set of published articles ( link 1 and link 2). The recommended force field for nucleic acids (DNA) has the name OL15, which is named after the city of Olomouc where it was developed in the year 2015. We do that using the LEaP module of AmberTools. ![]() Next we need to generate the topology and the coordinate files for the system so we can run our simulation. Open the file using VMD so you can see the structure that we are going to use for our tutorial. The structure has the sequence d(CGCGAATTCGCG), which is known as the Drew-Dickerson dodecamer, or DDD, and is one of the most famous and studied short oligomers available. Generating the initial topology and coordinate filesĭownload the initial PDB structure file (A-dna.pdb) available here. The starting structure is a type A-DNA generated using the Nucleic Acids Builder tool available as part of the AmberTools collection of programs. After the simulation we will use CPPTRAJ to analyze the characteristics of the conformational switch during the simulation. In this tutorial, we will run a DNA molecule with 12 base pairs starting in the A-DNA configuration and study the conversion to the B-DNA form. This example consists of a full set of instructions that will guide the user on how to run a simulation using AMBER and performing the analysis with CPPTRAJ. ![]()
0 Comments
Leave a Reply. |
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |