I am trying to pull apart two dimers of tubulin protein joined together, to form two separate dimers by steered molecular dynamics using Pull code. The objective is to separate the dimers by pulling along one axis only (Y axis). I ran the simulation for 25 ps and it accomplished well but while visualizing the trajectory in VMD, in the very beginning (first frame), dimers are at distance which is much greater than expected. I expected that in first frame they should not have separated too much (i.e. in just 1 ps, they are at greater distance) I guess this is unusual. Please let me know where I am doing wrong? I have used the following pull code:<br /><br />pull            = umbrella<br />pull_geometry   = distance<br />pull_dim        = N N Y<br />pull_start      = yes          <br />pull_init1      = 0.0<br />pull_ngroups    = 1<br />pull_group0     = dimer_ref               ; reference group (based on index file)<br />pull_group1     = dimer_pull           ; pull group<br />pull_rate1      = 0.1           <br />pull_k1         = 1000          <br />pull_nstfout    = 10           <br />pull_nstxout    = 10<br /><br />best regards,<br />Shilpi Chaurasia