# Using parameters of a drug from supplementary information of paper **URL:** \ **Category:** User discussions\ **Tags:** pdb2gmx, forcefield\ **Created:** [January 27, 2022, 8:08pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558 "2022-01-27T20:08:17Z")\ **Posts on this page:** 13\ **Page:** 1
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [January 27, 2022, 8:08pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/1 "2022-01-27T20:08:17Z")
GROMACS version:2016.3 GROMACS modification: No Here post your question Hello, I wanted to model the drug-bound c-Src kinase with the dasatinib bound drug (RCSB: 3G5D) in GROMACS. For dasatinib parameters, I used the parameter values provided in Supplementary Information of “Role of Desolvation in Thermodynamics and Kinetics of Ligand Binding to a Kinase” (Figure attached). I have made an rtp file, and a modified pdb file ([github gist](https://gist.github.com/shaunak-badani/34012f2512fbf514e28d1fe80362b924)) to name the atoms in dasatinib as per the parameters (residue name 1N1), which I made from the extent of my understanding in order to successfully build my system in gromacs. I wanted to use the OPLS-AA force field for the rest of the system. ``` Is this the correct methodology? How do I pass this rtp file to pdb2gmx, in order for it to pick up drug parameters? ``` Kind regards, Shaunak ![Parameters_SI](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/5/5cbbf3e87b7cda137b00e8a5c711757dbd74220e.png) ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [January 28, 2022, 1:37pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/2 "2022-01-28T13:37:02Z")
Update : I tried doing the following : 1. Copied the rtp file to “/usr/local/share/gromacs/top/oplsaa.ff”, made an atp file containing atom types and their masses. ([gist](https://gist.github.com/shaunak-badani/34012f2512fbf514e28d1fe80362b924)) 2. Ran pdb2gmx using OPLS-AA force field and TIP4P water. It gave me the following error : Residue 1 named 1N1 of a molecule in the input file was mapped to an entry in the topology database, but the atom CA used in that entry is not found in the input file. Perhaps your atom and/or residue naming needs to be fixed. There is no CA entry in 1N1, residue database, or in the input file. How do I interpret this error? Regards, Shaunak. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [January 28, 2022, 1:53pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/3 "2022-01-28T13:53:23Z")
You will need to provide the content of the 1N1 `.rtp` entry as well as the complete screen output of `pdb2gmx`. ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [January 28, 2022, 2:10pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/4 "2022-01-28T14:10:50Z")
Respected sir, Thank you for your reply. I was not able to upload files earlier, hence I provided a link to a github gist in the previous post. Attached is the `dasatinib.rtp` file I made. I realize I did not provide the `[bonded types]` directive, but it assumed some type for now. I wish to understand how gromacs does the modelling, so I can provide the parameters in the right way. [dasatinib\_atp.dat](https://gromacs.bioexcel.eu/uploads/short-url/hsn2xmmMjN1b5gjE2D5RgFuYh3H.dat) (486 Bytes) [dasatinib\_rtp.dat](https://gromacs.bioexcel.eu/uploads/short-url/vxMlQgzxbTLBxX2wUYlN9ibARH7.dat) (13.4 KB) pdb file (final\_drug\_bound\_1.pdb): [final\_drug\_bound\_1\_pdb.dat](https://gromacs.bioexcel.eu/uploads/short-url/Rys2pKfSe34PIyIfyfNGJnci8w.dat) (385.5 KB) All the files are in dat format because of uploading constraints Screenshot : ![Screenshot from 2022-01-28 19-38-33](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/f/f0fb02847e8ca9be0dddfc52f0ea7df75adef598.png) Command executed : gmx pdb2gmx -f final\_drug\_bound\_1.pdb -o 3g5d\_processed.gro Any help would be greatly appreciated. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [January 28, 2022, 6:46pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/5 "2022-01-28T18:46:00Z")
Please provide the **full** screen output from `pdb2gmx`. The error is not necessarily that informative. Everything that came before it is. There is a huge amount of diagnostic information that `pdb2gmx` prints that may be useful in pinpointing the problem. ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [January 29, 2022, 5:32am UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/6 "2022-01-29T05:32:01Z")
Sir, Here is the full screen output from the command executed : ![final](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/3/33eff94e63807a63fd4d07bd87d4035b6453ed36.jpeg) ![final_2](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/b/b7770a7e4ce67ce89babeb9c8090a99d6a5d3775.png) Kind regards, Shaunak. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [January 29, 2022, 1:38pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/7 "2022-01-29T13:38:53Z")
Are you specifying 1N1 as Protein in `residuetypes.dat`? If so, don’t. `pdb2gmx` is trying to apply NH3+ and COO- patches to it as if it were a protein residue. That’s why you’re getting a fatal error. In the future, please simply copy and paste the text rather than providing a series of screenshots with tiny text that is hard to read. ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [January 29, 2022, 5:31pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/8 "2022-01-29T17:31:40Z")
Respected Mr. Lemkul Apologies for the format of the data. I fixed the residuetype of 1N1 in residuetypes.dat to ‘Other’, and made some changes to the pdb. Attached is the pdb I now pass to pdb2gmx (GOL removed, TIP renamed to SOL) : [final\_drug\_bound\_1\_pdb.dat](https://gromacs.bioexcel.eu/uploads/short-url/9Ul1PUUfxIBvGFPBJVluZLXubbo.dat) (384.4 KB) When I execute the same command, it doesn’t give an error on 1N1, but rather related to terminal residues, even after using the -ter option (GLY: 0 - NH3+, LEU: 0 - COO-) : [pdb2gmx\_out\_dump.dat](https://gromacs.bioexcel.eu/uploads/short-url/8JBYhgb4GkNvaFtIeQxuk1IKjpR.dat) (10.7 KB) I also tried modelling the terminal residues using CHARMM-GUI, but It provided an extra two atoms at the end of the protein, at the amino acid LEU, OT1 and OT2. I believe this is the COO- end for the protein. When I tried passing it through pdb2gmx, it said “OT1 was not found in residue LEU”, which is fair, considering the termini were added by CHARMM-GUI. Kindly guide me as to how I should proceed. And thank you for your support so far. Kind regards, Shaunak. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [January 31, 2022, 11:08am UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/9 "2022-01-31T11:08:13Z")
> [@shaunak](#): > > When I execute the same command, it doesn’t give an error on 1N1, but rather related to terminal residues, even after using the -ter option (GLY: 0 - NH3+, LEU: 0 - COO-) The error has nothing to do with the termini, it’s incorrect nomenclature: ``` ------------------------------------------------------- Program: gmx pdb2gmx, version 2021.5 Source file: src/gromacs/gmxpreprocess/pdb2gmx.cpp (line 790) Fatal error: Atom OH2 in residue HO4 3 was not found in rtp entry HO4 with 4 atoms while sorting atoms. . For more information and tips for troubleshooting, please check the GROMACS website at http://www.gromacs.org/Documentation/Errors ------------------------------------------------------- ``` You’ve chosen TIP4P as your water model, which requires atoms OW, HW1, HW2, and MW. You’ve supplied atoms named OH2, H1, and H2 (e.g. CHARMM TIP3P nomenclature). Therefore, `pdb2gmx` fails to find the atoms it needs. I will also note that you may have issues with the topology you’re trying to generate anyway, because `pdb2gmx` will try to follow the `.rtp` entry for TIP4P (HO4), which does not specify a virtual site construction. You get this from tip4p.itp, but that means you should not include the water molecules in what you pass to `pdb2gmx` (though in this case, you have to generate the MW sites yourself). ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [February 1, 2022, 5:36am UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/10 "2022-02-01T05:36:16Z")
Respected Mr. Lemkul, Thank you for your input. I ended up removing the initial water molecules and solvating the system later. The `pdb2gmx` command ran successfully. However, after specifying the box size, solvating the system (`editconf, solvate`), I ran `grompp` to generate a `.tpr` file, which gave me an error `Atomtype C4 not found`. Here is the output dump for `grompp` : [grompp\_dump.dat](https://gromacs.bioexcel.eu/uploads/short-url/bQecNfcIpX99yerJpjTy4aGnepB.dat) (3.1 KB) I understand this error as the paper has introduced new atomtypes for the drug. When I searched online, it said I would have to make an `itp` file and specify (mass, sigma\_LJ, epsilon\_LJ and the charge) for each atom. 1. I have already provided LJ parameters in the rtp file, so is an itp file needed in this case? To phrase the question differently, how does one introduce a new atom type in a force field? 2. How can I estimate the charge value for each atom? That is not provided in the parameters of the paper. I could copy the charges from a different force field, but would that be the correct approach? 3. Also, how do I assign the `chargegroup` field for all atoms in the rtp file? I had allotted each atom its own chargegroup initially to see if gromacs accepted the file, which it did. But its not the right parameterization for the molecule. Kind regards, Shaunak. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [February 1, 2022, 11:53am UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/11 "2022-02-01T11:53:19Z")
> [@shaunak](#): > > I have already provided LJ parameters in the rtp file, so is an itp file needed in this case? To phrase the question differently, how does one introduce a new atom type in a force field? There are no LJ parameters in `.rtp` files. You need to add them to `ffnonbonded.itp` under `[atomtypes]` or in an auxiliary `.itp` file in the system topology that defines the new parameters > [@shaunak](#): > > How can I estimate the charge value for each atom? That is not provided in the parameters of the paper. I could copy the charges from a different force field, but would that be the correct approach? Third time this week I’ve said this :) You cannot mix and match different force fields. A force field is a self-consistent entity and each has different parametrization strategies. If you need to develop parameters for any new species, you need to do some thorough reading about how to do this for your chosen force field. There are tools available for many force fields to aid in this process, but it is extremely odd to me that a published study provided new LJ parameters for something but not the charges. The two go hand in hand. > [@shaunak](#): > > Also, how do I assign the `chargegroup` field for all atoms in the rtp file? I had allotted each atom its own chargegroup initially to see if gromacs accepted the file, which it did. But its not the right parameterization for the molecule. Charge groups are irrelevant in modern versions of GROMACS. ---
**Author:** ![shaunak](https://avatars.discourse-cdn.com/v4/letter/s/34f0e0/32.png) [@shaunak](https://gromacs.bioexcel.eu/u/shaunak)\ **Post date:** [February 3, 2022, 11:33am UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/12 "2022-02-03T11:33:08Z")
Respected Mr. Lemkul, Thank you for your input. I generated an `itp` file and grompp worked. Apologies for stating that the charges weren’t present in the SI, because they were. In an itp file, each atom type is given a charge. But in the parameters provided, there were multiple atoms of 1 type with different charges, hence I took an arbitrary value of charge for the atom type and placed it in the itp file. I am banking on the fact that the charges specified in the `rtp` file override the ones in the `itp` file, as stated in the [gromacs documentation](https://manual.gromacs.org/documentation/2019/reference-manual/topologies/pdb2gmx-input-files.html#residue-database). I ran into another error during grompp : `Proper Dih. multiplicity can not be perturbed grompp.` This error occurred because gromacs assumed I was providing parameters for [proper dihedrals of the periodic type](https://manual.gromacs.org/documentation/2019/reference-manual/functions/bonded-interactions.html#proper-dihedrals-periodic-type). I wanted to specify the use of Ryckaert-Bellemans dihedrals, which is function type 3. I did provide that in `dasatinib.rtp` and ran all of the steps again. I had to manually change the function from `1` to `3` in the `itp` file which gromacs generates after executing the pdb2gmx, editconf and solvate commands. How do I specify which functions to use without manually fiddling with the itp file generated by gromacs? Nonetheless, I have successfully generated the gro and tpr files, and just need to verify if the system is modelled correctly. Thank you so much for your help and support. Kind regards, Shaunak. ---
**Author:** ![jalemkul](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/jalemkul/32/18_2.png) [@jalemkul](https://gromacs.bioexcel.eu/u/jalemkul)\ **Post date:** [February 3, 2022, 6:48pm UTC](https://gromacs.bioexcel.eu/t/using-parameters-of-a-drug-from-supplementary-information-of-paper/3558/13 "2022-02-03T18:48:20Z")
> [@shaunak](#): > > How do I specify which functions to use without manually fiddling with the itp file generated by gromacs? The third column of the `[bondedtypes]` directive of the `.rtp` file sets the dihedral function type.