From: Giacomo Fiorin (giacomo.fiorin_at_gmail.com)
Date: Tue Jan 17 2017 - 20:47:39 CST
What you describe can't be a torsional angle, unless you define an
additional axis to represent the line where the two hypersurfaces (planes)
intersect.
You can, however, compute the angle of the quasi-rigid rotation (e.g.
orientationAngle for the small-molecule ligand) in a rotated frame that is
defined by the protein as fittingGroup. There is an example for the use of
fittingGroup here:
https://github.com/colvars/colvars/blob/master/examples/07_permeation.colvars.in
For proteins, this is also numerically more stable than computing the
inertia eigenvectors because the identity of the atoms is tracked by doing
the same optimal alignment done in computing minimal RMSD.
For more details, please read the doc:
http://colvars.github.io/colvars-refman-namd/colvars-refman-namd.html
On Tue, Jan 17, 2017 at 8:11 PM, Stefano Guglielmo <
stefano.guglielmo_at_unito.it> wrote:
> Dear NAMD users, is it possible to define a collective variable as a
> torsional angle between the principal axes of two molecules in a
> protein-ligand system?
>
> Thanks in advance
> Stefano.
>
> --
> Stefano GUGLIELMO PhD
> Assistant Professor of Medicinal Chemistry
> Department of Drug Science and Technology
> Via P. Giuria 9
> 10125 Turin, ITALY
> ph. +39 (0)11 6707678
>
>
-- Giacomo Fiorin Associate Professor of Research Institute for Computational Molecular Science (ICMS) College of Science and Technology, Temple University 1925 North 12th Street (035-07), Room 704D Philadelphia, PA 19122-1801 Phone: +1-215-204-4213 Scholar: http://goo.gl/Q3TBQU Personal: http://giacomofiorin.github.io/ Lab page: https://icms.cst.temple.edu/members.html *"As computer programmers we have a responsibility to make sure that we run the computers instead of the computers running us."* - Steve Oualline
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:02 CST