Re: constant pulling velocity in SMD

From: JC Gumbart (
Date: Thu Apr 26 2012 - 18:15:42 CDT

The velocity is in A/timestep. Assuming you're using a 1 fs timestep, that means you're pulling at the rate of 100 A/ns, 10x what you expect.

How well SMD controls the speed depends on the size of the force constant. See the tutorial Forcing Substrates Through Channels here:

On Apr 26, 2012, at 11:33 AM, Jiao, Dian NMN (-EXP) wrote:

> Hi, I was trying to pull this metal ion through an ion channel along Z-axis in SMD simulation. SMDvel was set at 0.0001, which means it will travel around 10 A in 1ns. However, my ion travels a whole lot faster. It is first placed near the top of the channel (extracellular surface), it goes from 15.233 to 7.753 (z direction) in merely 50 fs. It doesn't slow down until it hits this spot where the binding between the ion and the lumen facing residues (5 Glu) is very strong. This position is fairly close to the bottom of the channel (intracellular surface). The ion seems to be stablized there and never gets past this position throughout 5ps simulation. Isn't SMDvel supposed to control the speed of the ion moving? Why is it still going so fast? What do I do to make it easy for the ion to escape the binding of that glutamic ring and keeps going down? Thanks in advance.
> Dian

This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:21:55 CST