Re: DRMS Colvar

From: Giacomo Fiorin (
Date: Wed Sep 09 2020 - 14:55:47 CDT

Hi Ethan, the closest that comes to it in Colvars is the *distancePairs*
which returns the raw NxM matrix as a flattened vector. You could feed it
into a script function:
to do the squared difference from the reference distance matrix. Note that
it would be slow for large groups, because distancePairs was originally
meant to support pairs of few atoms (or copies thereof) such as those used

If you can afford multi-dimensional sampling you could also work with the
tilt and spinAngle functions:
which effectively measure the two relevant angles of relative orientation
(in addition to the COM-COM distance):


On Wed, Sep 9, 2020 at 3:28 PM Ethan Croitoru <> wrote:

> Hello,
> I would like to use the DMRS colvar defined in this paper, *Convergence
> and Sampling in Determining Free Energy Landscapes for Membrane Protein
> Association*, by Domanski et. al. (
> to do free energy
> calculations on a TMD helical dimer, like their GpA example. I know
> the paper used plumed for their colvars, but plumed only works with
> constant volume simulations on namd. So, I was wondering if there is a namd
> colvar that works similarly. So far I have thought of using distance
> colvars to pair residues from each helix, and use those distances to
> calculate the rmsd myself. Let me know what you think, any help would be
> appreciated.
> Thanks in advance,
> Ethan

Giacomo Fiorin
Associate Professor of Research, Temple University, Philadelphia, PA
Research collaborator, National Institutes of Health, Bethesda, MD

This archive was generated by hypermail 2.1.6 : Thu Dec 31 2020 - 23:17:14 CST