If it is indeed imaging, you'll want to add 'autoimage' *before* your
'rms' command. Also note that if your intent is to speed up commands
like 'autoimage' etc, you will want to put the 'strip' command before
everything else, e.g.:
trajin md_simulation_tac-comb.nc
strip :WAT,Cl-
autoimage
rms first :1-526.CA,C,N out rmsd-tac.dat
atomicfluct out rmsf-tac.dat :
1-526.CA,C,N byres
Hope this helps,
-Dan
On Fri, Jul 8, 2016 at 6:28 AM, Saman Yousuf ali
<saman.yousufali64.yahoo.com> wrote:
> Dear vlad,
> Thank you for your help. My rmsd y-axis is in angstrom. I have run backbone atoms rmsd. Script is pasted below.
> trajin md_simulation_tac-comb.nc
> rms first :1-526.CA,C,N out rmsd-tac.dat
> atomicfluct out rmsf-tac.dat :1-526.CA,C,N byres
> strip :WAT
> strip :Cl-
> Best Regards, Saman Yousuf AliJunior Research Fellow,
> | Lab No. P-133, Computational Chemistry Laboratory
> Dr. Panjwani Center for Molecular Medicine & Drug Research,
> International Center for Chemical & Biological Sciences,
> University of Karachi – 75270.Karachi-Pakistan.
>
> Contact No:
> Office (92-21) 111222292 (Ext 309)
> Email ID: saman.yousufali64.yahoo.com
> saman.ali.iccs.edu
>
> |
>
>
>
> On Friday, July 8, 2016 5:17 AM, Vlad Cojocaru <vlad.cojocaru.mpi-muenster.mpg.de> wrote:
>
>
> This looks to me like an imaging issue with that frame (or those few frames) .... check the frame(s), you should be able to see in a visualization software what happens. If it is an imaging issue, simply remove it (them) from the processed trajectory.... If for some reason you really need exactly those frames affected, try some manual imaging procedures ... or re-run the simulation (5 ns should be doable)
>
> Just a remark ... if your y axis is nm, the rmsd looks to quite large during the entire simulation.... How did you calculate the RMSD ? What parts of the system did you fit ? Is this RMSD for the fitted parts or for other parts ?
>
> Best
> Vlad
>
>
>
> On 07/08/2016 01:42 PM, Saman Yousuf ali wrote:
>
> Dear zoran,
> This rmsd plot is generated after 5ns production run. I have performed protein ligand complex md simulation. When I reimaged my mdcrd after production it gives me following error for 1ns trajectory; Warning: autoimage: Frame 20 imaging failed, box lengths are zero. While rest of proceed normally without this warning. you see in attached file that the rmsd jumps to near20 angstroms. This is a sign of something seriouslywrong with the trajectory.
> Best Regards, Sarah
>
>
> On Friday, July 8, 2016 4:17 AM, zoran matovic <zmatovic.kg.ac.rs> wrote:
>
>
> Hi,
> What kind of MD (heating, equil., production)? What kind of molecule is the
> "complex"? What was the reference structure you have used for RMSD? What's
> up to me you have nice equilibrium after 5 ns.
>
> cheers
> zoran
>
> Dear all,
> I have performed 5ns md simulation using amber14. Rmsd plot of one my
> simulated complex shows unusua 20 angstrom fluctuation at starting ns. I am
> unable to understand the reason. Rmsd plot is attached below. Kindly help
> me. Best Regards, Sarah
>
>
>
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
>
>
>
>
> _______________________________________________AMBER mailing listAMBER@ambermd.orghttp://lists.ambermd.org/mailman/listinfo/amber
>
> --
> Dr. Vlad Cojocaru
> Computational Structural Biology Laboratory
> Department of Cell and Developmental Biology
> Max Planck Institute for Molecular Biomedicine
> Röntgenstrasse 20, 48149 Münster, Germany
> Tel: +49-251-70365-324; Fax: +49-251-70365-399
> Email: vlad.cojocaru[at]mpi-muenster.mpg.de
> http://www.mpi-muenster.mpg.de/43241/cojocaru
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
--
-------------------------
Daniel R. Roe, PhD
Department of Medicinal Chemistry
University of Utah
30 South 2000 East, Room 307
Salt Lake City, UT 84112-5820
http://home.chpc.utah.edu/~cheatham/
(801) 587-9652
(801) 585-6208 (Fax)
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Jul 08 2016 - 08:00:03 PDT
