Atoms too fast/periodic cell too small with ABF protein-ligand

From: Francesco Pietra (chiendarret_at_gmail.com)
Date: Sun Jun 25 2017 - 05:56:01 CDT

Hello:

I am attempting a protein-ligand ABF, following the 2017 tutorial, by using
a 58.2ns problemless equilibrated system (ts=1.0 fs, bond restriction on
water only) in a TIP3P water box on a main pure-CPU cluster. Rather large
protein, organic ligand as accurately parameterized as I could, by fitting
torsions and water interaction.

I am experiencing immediate "atoms moving too fast" (two H atoms of the
ligand) when using a linux 4-core cpu desktop, or "periodic cell has become
too small" on a linux GPU workstation. At the moment I have no access to
the cluster.

I used ts=1.0 fs, i.e. no bond restriction, except for TIP3P water, as in
the equilibration.

As a possible cause, that I was unable to verify, is the setting of Euler
and polar angle in the colvars definition. That is , I used the following
values from "cv printframe"

>Main< (Conformation:5syf-SAA1) 28 % cv printframe
> 0 1.88926419197828e+01 5.22845307973583e+01
> -3.39612142392650e+01 9.26730407864645e+00 1.01135055762437e+02
> -5.05311208874754e+01 0.00000000000000e+00 0.00000000000000e+00
>

by discarding the first two, and using the directly measured intercenter
distance of 27.6A in place of 18.89 from printframe, i.e., as follows:

harmonic {
   colvars r
   forceConstant 0.0
   centers 27.6 # OK measured
}

harmonic {
   colvars Theta
   forceConstant 0.0
   centers 52.3 # from printframe
}

harmonic {
   colvars Phi
   forceConstant 0.0
   centers -40.0 # from printframe
}

harmonic {
   colvars Psi
   forceConstant 0.0
   centers 9.27 # from printframe
}

harmonic {
   colvars theta
   forceConstant 0.0
   centers 101.1 # from printframe
}

harmonic {
   colvars phi
   forceConstant 0.0
   centers -50.5 # from printframe
}

Should this assignment of colvars be correct, where to look for the causes
of the instability of the system?
I must confess that I am the first time at such an ABF beyond the tutorial

Thanks for advice

francesco pietra

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:23 CST