- Posts: 416
Philosophy of the martinilize
- Leliel
- Topic Author
- Visitor
Thanks a lot!
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
Please let us know if you get a reliable model of DNA based on Martini.
Leliel wrote: Hi,i am vary interested in martini FF and want to martinilise the DNA in my own philosophy.but i cant know how to get info in the atomic pdb files and change them into the cg .gro and .top files.So who can tell me how to do this? not a script but the core method.
Thanks a lot!
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
And is it right that once you settle down the position of every CG bead ,just use them to calculate the bond distance,angle,dihedral angle ...these parameters?
then how to parameterize the connection constant?like some k ?
If we get some progress use this model,i will surely inform you about this.
Again Thanks!
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
You would run a atomistic model simulation to extract the force constant of the bonds/angle/dihedrals ...
Leliel wrote: Thanks a lot!
And is it right that once you settle down the position of every CG bead ,just use them to calculate the bond distance,angle,dihedral angle ...these parameters?
then how to parameterize the connection constant?like some k ?
If we get some progress use this model,i will surely inform you about this.
Again Thanks!
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
I am an undergraduate student and these questions may be stupid ,but can you tell the procedure of constructing a martini model??
Thanks
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
You have to use an iterative procedure in which you will refine your topology at each step. You first need to collect all the experimental data you can on the molecule you want to parameterise. Partitioning data (between different environments) is fundamental. You then want (if possible) to run an atomistic simulation of the extract bonded terms. From this all you can build your first CG topology and run an MD simulation from which you can extract data to compare to the data you collected ... it is simple after you go through the process a few times :))
good luck.
Leliel wrote: Well,I have read most of your published paper,and on that it said that the nonbond interaction are parameterized with some thermodynamical datas and bond interaction are parameterized with the AT MD structure,but the papers shows little about the details of how to do this...
I am an undergraduate student and these questions may be stupid ,but can you tell the procedure of constructing a martini model??
Thanks
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
1.define the position of CG beads(use the collected beads' center of mass)
2.use the position of CG beads to calculate the distance of bonds,angles,dihedral angels.
3.compare these parameters to the experimental datas(or AT MD structure) using the iterative way to refine them.
that may be similar in the nonbond-interaction paramaters ?
And pls the name of book or chapter (or the website of the book)so i can get it and learn by myself.thanks in advance!
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
The book chapters I mention are:
- S.J. Marrink, M. Fuhrmans, H.J. Risselada, X. Periole. The MARTINI force field. In "Coarse graining of condensed phase and biomolecular systems", G. Voth ed., CRC press, Chapter 2, 2008.
- X. Periole, S.J. Marrink. The Martini coarse-grained force field. In "Methods in molecular biology", Vol 924, L. Monticelli & E. Salonen Eds., Springer, 2013, pp 533-565.
You can find a reprint on the martini webpage (cgmartini.nl) in the publication section.
Leliel wrote: Am i right to explain the procedure in this way?
1.define the position of CG beads(use the collected beads' center of mass)
2.use the position of CG beads to calculate the distance of bonds,angles,dihedral angels.
3.compare these parameters to the experimental datas(or AT MD structure) using the iterative way to refine them.
that may be similar in the nonbond-interaction paramaters ?
And pls the name of book or chapter (or the website of the book)so i can get it and learn by myself.thanks in advance!
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
thx!
Please Log in or Create an account to join the conversation.
- xavier
- Offline
- Admin
- Posts: 416
You have to get the tutorial and go through it. It will show you how to build atomistic and CG topologies to use in combination to go from an atomistic to a CG trajectory. This way you can extract the bonded terms (bonds, angles etc).
Leliel wrote: Thanks!i have already read these two chapters and get to know the whole 3step-recipe of building a new topology.But i still couldn't find out by what programs or tools you make these steps. Like which script or tool you will use to calculate the center of mass position or the distribution ?
thx!
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
i think you mention the "adaptive resolution scheme"section in Gromacs tutorial,right?
So i will install a "votca"package and use iterative boltzmann inversion method to get the parameters?
i will try to do it,thanks!
Please Log in or Create an account to join the conversation.
- Clement
- Offline
- Admin
- Posts: 211
To get the value of your bonded potentials (bond lengths, angles, dihedrals), you need to:
- run an atomistic (AA) simulation of the whole and/or some parts of your molecule, long enough to have converged values for the next step (for you to test; depending on the size of your molecule, can be nanoseconds to hundreds of nanoseconds).
- convert your AA trajectory into a CG trajectory; to be able to do that, you need a AA > CG mapping of your molecule, i.e. which AA atoms go in which CG bead (you need to write that by yourself; but it's the fun part of the parameterization). The reprocessing of the trajectory is done using this code: http://md.chem.rug.nl/cgmartini/index.php/downloads/tools/113-rt , following this tutorial: http://md.chem.rug.nl/cgmartini/index.php/tutorial/reverse-transformation . You'll find gazillions of posts concerning this reverse transformation on this Forum.
- extract the values you need: with the tools (g_bond, g_angle, etc.) provided by GROMACS (that's what we usually use, the rest being done with the same simulation package), you can measure distributions of bond length, angle, etc.; this will give you the first set of parameters for your trial topology.
In case you didn't see it already, these steps are explained in more details here: http://md.chem.rug.nl/cgmartini/index.php/user-platform/faq#create topos .
Good luck!
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
i have read through the tutorials of 1e. and 3. both mentions the resolution transformation.but it says that i the package is based upon gromacs3.3.1 ,but the version of gromacs on my computer is 4.6.2. so is it OK to install two version of gromacs? or the 4.6.2 can afford the g_fg2cg function??
Please Log in or Create an account to join the conversation.
- Clement
- Offline
- Admin
- Posts: 211
No problem to have multiple versions of GROMACS on your machine; you need to use the different executables (mdrun, g_fg2cg, etc.) giving the full path (i.e. /wherever/is/your/gromacs-3.3.1/bin/mdrun), or source it in your current bash session (i.e. source /wherever/is/your/gromacs-3.3.1/bin/GMXRC).
Please Log in or Create an account to join the conversation.
- Leliel
- Topic Author
- Visitor
And now i wanna extract some bond parameters using this tool.At hand i have the cg.gro file of an DNA,using a cpp file written by myself, and a .top and .itp files written artificially.
I summarize the steps of this recipe, could pls help me to validate this?
1) do a long-time scale-fine grained MD of this DNA,long enough to extract credible parameters(100ns or somewhat?)
2) use the fg2cg to extract a coarse grained trajectory of this MD.
3) make a index using the make_ndx ,mentioned which i would like to count.
4) use g_angle to extract the gold-standard angle and dihedral as a ref.
5) use my artificially .gro to run a MD ,and also extract some angle/dihedral ,and alter the parameters until finely accord to the gold-standard.
am i right??
as well i have some questions about this recipe's details. how can i put this 100ns MD? in the water(solution) or just do a MD only with this DNA?
thanks in advance,for helping me do a coarse-grained recipe step by step, enforce my knowledge.
Please Log in or Create an account to join the conversation.
- djurre
- Offline
- Admin
- Posts: 272
However, I would also like to make a warning: I won't be as easy as it sounds here. This protocol works very well for small molecules, but is much harder for larger biopolymers. For example the structures of proteins is not stable without an elastic network. A similar thing you might encounter for DNA. Second, the bases of DNA are very tightly packed. This might not be possible with regular of S-type Martini beads. These are only some of the problems you will face and have to solve.
In our lab, someone is working on DNA, and even thought he is one of the four smartest PhD students we have, he is still breaking his head on it. So be warned!
Please Log in or Create an account to join the conversation.