Re: [gmx-users] Topology file for phosphothreonine (TPO) for OPLS force field
Thank you very much for your quick response, Justin. I guess I need a little bit more practise to make the parametrization myself. Do you have any clues on where I could find published parameters? Are there databases for this or something similar? Regards, Antoine Le 5/02/2010 18:23, Justin A. Lemkul a écrit : Antoine Delmotte wrote: Dear gmx-users, I would like to simulate the CDK9 protein (3BLH in PDB), but there is one residue (TPO, phosphothreonine), which Gromacs does not recognise. I found the PRODRG website as well as force field parameters for TPO for the Gromos96.1 force field. However, I would like to use the OPLS-AA force field. I therefore need a topology file for this force field, which cannot be provided by PRODRG (only useful for united atoms force fields). Does anyone have an idea about where I could get the phosphothreonine itp file for OPLS force field? Search the literature to see if anyone has published parameters for it. Otherwise, could someone give me clues or the link to a good tutorial about how to create this itp file myself? That's what Chapter 5 of the manual is for, provided you have reliable parameters (see below). Molecular dynamics is something completely new for me so, sorry if I ask stupid questions. It also seems I am not the first one to have this problem but I could not find any solution on the archive of the mailing list. If MD is completely new to you, and you are unable to find previously published and validated parameters, the parameterization of a new species is likely not a suitable task. Parameterization is an expert topic, and it can take many months or years to develop reliable parameters, even if you know what you are doing at the very start. Be forewarned: http://www.gromacs.org/Documentation/How-tos/Parameterization -Justin I thank you in advance. Best regards, Antoine Hotmail : une messagerie fiable avec la protection anti-spam performante de Microsoft Inscrivez-vous https://signup.live.com/signup.aspx?id=60969 -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] Topology file for phosphothreonine (TPO) for OPLS force field
Thank you very much for this, Bruce. It will help me, for sure. Could you just let me know where you found the parameters or how you generated them? Thanks in advance. Regards, Antoine Le 5/02/2010 22:39, Bruce D. Ray a écrit : On Fri, February 5, 2010 at 1:06:34 PM, Antoine Delmotte antoinedelmo...@hotmail.com wrote: I would like to simulate the CDK9 protein (3BLH in PDB), but there is one residue (TPO, phosphothreonine), which Gromacs does not recognise. I found the PRODRG website as well as force field parameters for TPO for the Gromos96.1 force field. However, I would like to use the OPLS-AA force field. I therefore need a topology file for this force field, which cannot be provided by PRODRG (only useful for united atoms force fields). Does anyone have an idea about where I could get the phosphothreonine itp file for OPLS force field? Otherwise, could someone give me clues or the link to a good tutorial about how to create this itp file myself? Molecular dynamics is something completely new for me so, sorry if I ask stupid questions. It also seems I am not the first one to have this problem but I could not find any solution on the archive of the mailing list. The charges probably need a little more work, but this might be a start towards an rtp file for phosphothreonine to add to ffoplsaa.rtp. [ TPO ] [ atoms ] Nopls_238 -0.500 1 Hopls_2410.300 1 CAopls_224B 0.140 1 HAopls_1400.060 1 CBopls_1580.205 2 HBopls_1400.060 2 OG1opls_447 -0.700 2 Popls_4451.795 3 O1Popls_446 -1.120 3 O2Popls_446 -1.120 3 O3Popls_446 -1.120 3 CG2opls_135 -0.180 4 HG21opls_1400.060 4 HG22opls_1400.060 4 HG23opls_1400.060 4 Copls_2350.500 5 Oopls_236 -0.500 5 [ bonds ] N H NCA CAHA CACB CA C CBHB CB OG1 CB CG2 OG1 P O1P P O2P P O3P P CG2 HG21 CG2 HG22 CG2 HG23 C O -C N [ dihedrals ] ; override some of the typebased dihedrals NCACB OG1dih_SER_THR_chi1_N_C_C_O CCACB OG1dih_SER_THR_chi1_CO_C_C_O CACB OG1 HG1dih_SER_THR_chi2_C_C_OH_HO [ impropers ] -CCA N Himproper_Z_N_X_Y CA+N C Oimproper_O_C_X_Y I hope that helps. -- Bruce D. Ray, Ph.D. Associate Scientist IUPUI Physics Dept. 402 N. Blackford St. Indianapolis, IN 46202-3273 -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php
Re: [gmx-users] Topology file for phosphothreonine (TPO) for OPLS force field
Thank you very much Fairuz. I will try this. Le 6/02/2010 1:20, fairuz zulkifli a écrit : to create itp file with OPLS force feild, you can try use GAMESS. -fairuz- -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/mailing_lists/users.php