Posts

Showing posts from May, 2017

Fixing bugs in FF14SB port for Gromacs

Bugs? The newly released amber14sb made modifications for protein, DNA and RNA. This force field could be ported into gromacs as I already introduced previous. However, there are some small bugs in the improper dihedral angles of HID, HIE and HIP. So when generating tpr file, errors would be prompted. When run grompp -f npt.mdp -p topol.top -c protein.gro -o product.tpr -n index, errors occur: ERROR 1 [file topolff14sb_Protein_chain_B.itp, line 208041]:   No default Improper Dih. types The reason is that the parameters for some dihedral angles are missing. So we need to fix this. How to fix? Add the following two lines in the ffbonded.itp after the 730th lines: NA  CW  CC  CT       4      180.00     4.60240     2     ;;; HID force field added NA  CV  CC  CT       4      180.00     4.60240     2     ;;; HID force field added The atom type CV and CW here are similar and share the same parameter value. Save the file and run the command again, in general, errors won&

Install new amber force fields ports in Gromacs

Image
Amber force fields are popular for protein, DNA, RNA and lipid simulations. Recently updated amberff12, amberff14 and amberff15 were reported to perform better than amber99sb and amber99sb-ildn. Gromacs is a commonly used simulation software, which using OpenMM to enable GPU acceleration. Gromacs is fast and easy to use. A lot of people use amberff in Gromacs. So to use the newly release amberff12 and amber14sb in gromacs, we need to some work. Step 1. Install Gromacs 4.6.x or Gromacs 5.0.x or Gromacs 5.1.x Installation instructions could be find here: For 4.6.x http://www.gromacs.org/Documentation/Installation_Instructions_4.6 For 5.0.x and 5.1.x http://www.gromacs.org/Documentation/Installation_Instructions_5.0 Suppose you are user john, you installed a gromacs 4.6.7 in /home/john/applications/. After installation, you should source the GMXRC by: source  /home/john/applications/gromacs.4.6.7/GMXRC And add topology files to environments by: export GMXLIB=/home/john/ap