Re: [AMBER] MMPBSA error with receptor prmtop but not clear how to fix Amber20

From: Elvis Martis <elvis_bcp.elvismartis.in>
Date: Mon, 14 Dec 2020 17:05:02 +0530

OK.
Have you compiled AMBER20 using the legacy configure method or the new
cmake method?
If you have using the new cmake method of compilation, then I guess you
should find that file in amber20_src/AmberTools/src/sander/mdread2.F90.
if you have compiled using the legacy method, then you files should in the
path where I sent previously.
Best Regards
Elvis



On Mon, 14 Dec 2020 at 16:55, Vaibhav Dixit <vaibhavadixit.gmail.com> wrote:

> Dear All and Elivs,
> I'm afraid, I still can't find this file inside the AmberTools folder,
> there is no folder by the same sander.
> All tests passed while installing Amber20, thus I think there is some
> confusion regarding the directory structure and file locations in Amber18
> and Amber20 as well.
> Please let me know if you know the fix to this as well.
> thank you and best regards.
>
> (base) [exx.c107739 amber20]$ ls
> $AMBERHOME/AmberTools/src/sander/mdread2.F90
> ls: cannot access
> /home/exx/Downloads/amber20/AmberTools/src/sander/mdread2.F90: No such file
> or directory
> (base) [exx.c107739 amber20]$ cd AmberTools/src/
> (base) [exx.c107739 src]$ ls -arlth
> total 28K
> drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 cpptraj
> drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 mm_pbsa
> drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 FEW
> drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 pytraj
> drwxrwxr-x. 7 exx exx 4.0K Dec 4 15:25 ..
> -rw-r--r--. 1 exx exx 1.9K Dec 4 17:26 config.h
> drwxrwxr-x. 6 exx exx 4.0K Dec 4 17:45 .
>
> On Mon, Dec 14, 2020 at 7:46 PM Elvis Martis <elvis_bcp.elvismartis.in>
> wrote:
>
> > Hi
> > Sorry my bad. For AMBER20
> > you can look in
> > $AMBERHOME/AmberTools/src/sander/mdread2.F90
> > Start looking from line number 1508 (this should be same unless there
> were
> > changes made earlier)
> > this line starts with " else if ( gbsa == 2 ) then"
> >
> > Best Regards
> > Elvis
> >
> >
> >
> > On Mon, 14 Dec 2020 at 15:57, Vaibhav Dixit <vaibhavadixit.gmail.com>
> > wrote:
> >
> > > Dear All,
> > > I got this suggestion earlier last month also.
> > > But to my surprise, I can't find these specified folders and files in
> the
> > > expected locations.
> > > I ran a find command to search for the mdread.f file, but it is not
> found
> > > anywhere inside amber20 folder and subfolders.
> > > This command works and finds other files as shown below.
> > >
> > > Please do let me know what I'm missing to understand here.
> > > thanks
> > >
> > > (base) [exx.c107739 amber20]$ ls -arlth
> > > total 128K
> > > -rw-r--r--. 1 exx exx 7.5K Dec 6 2018 GNU_LGPL_v3
> > > -rwxr-xr-x. 1 exx exx 116 Dec 6 2018 amber-interactive.sh
> > > -rw-r--r--. 1 exx exx 1.1K Apr 28 2020 Makefile
> > > -rw-r--r--. 1 exx exx 37K Apr 28 2020 LICENSE
> > > drwxrwxr-x. 8 exx exx 4.0K Dec 4 15:17 licenses
> > > drwxrwxr-x. 2 exx exx 4.0K Dec 4 15:17 doc
> > > drwxrwxr-x. 19 exx exx 4.0K Dec 4 15:17 miniconda
> > > drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 share
> > > drwxrwxr-x. 9 exx exx 4.0K Dec 4 15:17 benchmarks
> > > drwxrwxr-x. 16 exx exx 4.0K Dec 4 15:17 dat
> > > drwxrwxr-x. 7 exx exx 4.0K Dec 4 15:25 AmberTools
> > > -rw-r--r--. 1 exx exx 1.9K Dec 4 17:26 config.h
> > > drwxrwxr-x. 8 exx exx 4.0K Dec 10 15:54 logs
> > > drwxrwxr-x. 157 exx exx 4.0K Dec 10 15:58 test
> > > -rwxr-xr-x. 1 exx exx 1.8K Dec 10 20:38 amber.sh
> > > -rwxr-xr-x. 1 exx exx 1.9K Dec 10 20:38 amber.csh
> > > drwxrwxr-x. 15 exx exx 4.0K Dec 10 21:50 .
> > > drwxr-xr-x. 2 exx exx 4.0K Dec 10 21:50 include
> > > drwxrwxr-x. 4 exx exx 4.0K Dec 10 21:50 lib64
> > > drwxr-xr-x. 4 exx exx 4.0K Dec 10 21:50 lib
> > > drwxrwxr-x. 3 exx exx 4.0K Dec 10 21:50 bin
> > > drwxr-xr-x. 24 exx exx 4.0K Dec 12 12:13 ..
> > > (base) [exx.c107739 amber20]$ cd AmberTools/
> > > (base) [exx.c107739 AmberTools]$ ls -arlth
> > > total 28K
> > > drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:17 benchmarks
> > > drwxrwxr-x. 6 exx exx 4.0K Dec 4 15:17 examples
> > > drwxrwxr-x. 3 exx exx 4.0K Dec 4 15:25 .pytest_cache
> > > drwxrwxr-x. 7 exx exx 4.0K Dec 4 15:25 .
> > > drwxrwxr-x. 6 exx exx 4.0K Dec 4 17:45 src
> > > drwxrwxr-x. 15 exx exx 4.0K Dec 10 21:50 ..
> > > drwxrwxr-x. 65 exx exx 4.0K Dec 10 21:59 test
> > > (base) [exx.c107739 AmberTools]$ cd ..
> > > (base) [exx.c107739 amber20]$ find . -name mdread.f
> > > (base) [exx.c107739 amber20]$ find . -name amber.sh
> > > ./amber.sh
> > > (base) [exx.c107739 amber20]$ find . -name am1bcc
> > > ./bin/am1bcc
> > > ./bin/wrapped_progs/am1bcc
> > >
> > > On Mon, Dec 14, 2020 at 7:14 PM Elvis Martis <elvis_bcp.elvismartis.in
> >
> > > wrote:
> > >
> > > > Hi
> > > > Here is a solution in a different thread, you will have to edit the
> > > > mdread.f file. Please have a look at link below to know how to edit
> the
> > > > file and where to find it.
> > > > Re: [AMBER] bad atom type: Br from Jason Swails on 2011-11-11 (Amber
> > > > Archive Nov 2011) (ambermd.org)
> > > > <http://archive.ambermd.org/201111/0370.html>
> > > > Best Regards
> > > > Elvis
> > > >
> > > >
> > > >
> > > > On Mon, 14 Dec 2020 at 15:34, Vaibhav Dixit <vaibhavadixit.gmail.com
> >
> > > > wrote:
> > > >
> > > > > Dear All,
> > > > > I'm getting the following error message while trying to run
> MMPBSA.py
> > > in
> > > > > Amber20.
> > > > > It says that the receptor mask overwritten with default, which I
> > don't
> > > > know
> > > > > why it is doing even after I gave a specific receptor mask (1st
> > residue
> > > > is
> > > > > my ligand in the trajectory). It also gives the bad atom type
> error,
> > > > that
> > > > > I got with Amber18 last month.
> > > > >
> > > > > >From the output files (attached), it is not clear what is the
> exact
> > > > nature
> > > > > of the error, thus I have no clue how to proceed to fix the same.
> > > > > Can you please help me understand why the dry prmtop is not
> > compatible?
> > > > > And how can I possibly attempt to get this analysis done?
> > > > > My mmpbsa input file is given below.
> > > > > Thanking you all in advance.
> > > > > Best regards
> > > > >
> > > > > mmpbsa input file:
> > > > > (base) [exx.c107739 mmpbsa-6701-tyr]$ more mmpbsa-decomp.in
> > > > > &general
> > > > > startframe=1, endframe=100, interval=5,
> > > > > verbose=2, keep_files=2, netcdf=1, receptor_mask=2-386,
> > > > > strip_mask=:387-16750,
> > > > > /
> > > > > &gb
> > > > > igb=5, saltcon=0.150,
> > > > > /
> > > > > &pb
> > > > > istrng=0.15, fillratio=4.0,
> > > > > /
> > > > > &decomp
> > > > > idecomp=1, print_res="56-57; 74-75; 167; 222; 225-227; 229-231;
> > > > 272-273;
> > > > > 275-277; 336; 368; 370; 372-373"
> > > > > dec_verbose=1,
> > > > > /
> > > > >
> > > > >
> > > > > Error message on the terminal
> > > > > (base) [exx.c107739 mmpbsa-6701-tyr]$ MMPBSA.py -O -i
> > mmpbsa-decomp.in
> > > > -o
> > > > > mmpbs-rouf-6701-tyr.dat -do decomp-rouf-6701-tyr.dat -sp
> > > > > ../rouf-6701-tyr-solv.prmtop -cp ../rouf-6701-tyr-dry.prmtop -lp
> > > > > ../tyr-dry.prmtop -y ../rouf-6701-tyr-solv-prod300.nc -rp
> > > > > ../rouf-6701-nolig-dry.prmtop
> > > > > Loading and checking parameter files for compatibility...
> > > > >
> > > > >
> > > >
> > >
> >
> */home/exx/Downloads/amber20/lib/python3.8/site-packages/MMPBSA_mods/main.py:603:
> > > > > UserWarning: receptor_mask overwritten with default*
> > > > >
> > > > > warnings.warn('receptor_mask overwritten with default\n')
> > > > > cpptraj found! Using /home/exx/Downloads/amber20/bin/cpptraj
> > > > > sander found! Using /home/exx/Downloads/amber20/bin/sander
> > > > > Preparing trajectories for simulation...
> > > > > 20 frames were processed by cpptraj for use in calculation.
> > > > >
> > > > > Running calculations on normal system...
> > > > >
> > > > > Beginning GB calculations with
> /home/exx/Downloads/amber20/bin/sander
> > > > > calculating complex contribution...
> > > > > * bad atom type: M1*
> > > > > File "/home/exx/Downloads/amber20/bin/MMPBSA.py", line 100, in
> > > <module>
> > > > > app.run_mmpbsa()
> > > > > File
> > > > >
> > > > >
> > > >
> > >
> >
> "/home/exx/Downloads/amber20/lib/python3.8/site-packages/MMPBSA_mods/main.py",
> > > > > line 218, in run_mmpbsa
> > > > > self.calc_list.run(rank, self.stdout)
> > > > > File
> > > > >
> > > > >
> > > >
> > >
> >
> "/home/exx/Downloads/amber20/lib/python3.8/site-packages/MMPBSA_mods/calculation.py",
> > > > > line 82, in run
> > > > > calc.run(rank, stdout=stdout, stderr=stderr)
> > > > > File
> > > > >
> > > > >
> > > >
> > >
> >
> "/home/exx/Downloads/amber20/lib/python3.8/site-packages/MMPBSA_mods/calculation.py",
> > > > > line 156, in run
> > > > > raise CalcError('%s failed with prmtop %s!' % (self.program,
> > > > > CalcError: /home/exx/Downloads/amber20/bin/sander failed with
> prmtop
> > > > > ../rouf-6701-tyr-dry.prmtop!
> > > > > Exiting. All files have been retained.
> > > > >
> > > > > --
> > > > >
> > > > > Regards,
> > > > >
> > > > > Dr. Vaibhav A. Dixit,
> > > > >
> > > > > Visiting Scientist at the Manchester Institute of Biotechnology
> > (MIB),
> > > > The
> > > > > University of Manchester, 131 Princess Street, Manchester M1 7DN,
> UK.
> > > > > AND
> > > > > Assistant Professor,
> > > > > Department of Pharmacy,
> > > > > ▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄
> > > > > Birla Institute of Technology and Sciences Pilani (BITS-Pilani),
> > > > > VidyaVihar Campus, street number 41, Pilani, Rajasthan 333031.
> > > > > India.
> > > > > Phone No. +91 1596 255652, Mob. No. +91-7709129400,
> > > > > Email: vaibhav.dixit.pilani.bits-pilani.ac.in,
> > vaibhavadixit.gmail.com
> > > > > http://www.bits-pilani.ac.in/pilani/vaibhavdixit/profile
> > > > > https://www.linkedin.com/in/vaibhav-dixit-b1a07a39/
> > > > >
> > > > > ORCID ID: https://orcid.org/0000-0003-4015-2941
> > > > >
> > > > >
> http://scholar.google.co.in/citations?user=X876BKcAAAAJ&hl=en&oi=sra
> > > > >
> > > > > P Please consider the environment before printing this e-mail
> > > > > _______________________________________________
> > > > > AMBER mailing list
> > > > > AMBER.ambermd.org
> > > > > http://lists.ambermd.org/mailman/listinfo/amber
> > > > >
> > > > _______________________________________________
> > > > AMBER mailing list
> > > > AMBER.ambermd.org
> > > > http://lists.ambermd.org/mailman/listinfo/amber
> > > >
> > >
> > >
> > > --
> > >
> > > Regards,
> > >
> > > Dr. Vaibhav A. Dixit,
> > >
> > > Visiting Scientist at the Manchester Institute of Biotechnology (MIB),
> > The
> > > University of Manchester, 131 Princess Street, Manchester M1 7DN, UK.
> > > AND
> > > Assistant Professor,
> > > Department of Pharmacy,
> > > ▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄
> > > Birla Institute of Technology and Sciences Pilani (BITS-Pilani),
> > > VidyaVihar Campus, street number 41, Pilani, Rajasthan 333031.
> > > India.
> > > Phone No. +91 1596 255652, Mob. No. +91-7709129400,
> > > Email: vaibhav.dixit.pilani.bits-pilani.ac.in, vaibhavadixit.gmail.com
> > > ​http://www.bits-pilani.ac.in/pilani/vaibhavdixit/profile
> > > https://www.linkedin.com/in/vaibhav-dixit-b1a07a39/
> > >
> > > ORCID ID: https://orcid.org/0000-0003-4015-2941
> > >
> > > http://scholar.google.co.in/citations?user=X876BKcAAAAJ&hl=en&oi=sra
> > >
> > > P Please consider the environment before printing this e-mail
> > > _______________________________________________
> > > AMBER mailing list
> > > AMBER.ambermd.org
> > > http://lists.ambermd.org/mailman/listinfo/amber
> > >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
>
>
> --
>
> Regards,
>
> Dr. Vaibhav A. Dixit,
>
> Visiting Scientist at the Manchester Institute of Biotechnology (MIB), The
> University of Manchester, 131 Princess Street, Manchester M1 7DN, UK.
> AND
> Assistant Professor,
> Department of Pharmacy,
> ▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄▄
> Birla Institute of Technology and Sciences Pilani (BITS-Pilani),
> VidyaVihar Campus, street number 41, Pilani, Rajasthan 333031.
> India.
> Phone No. +91 1596 255652, Mob. No. +91-7709129400,
> Email: vaibhav.dixit.pilani.bits-pilani.ac.in, vaibhavadixit.gmail.com
> ​http://www.bits-pilani.ac.in/pilani/vaibhavdixit/profile
> https://www.linkedin.com/in/vaibhav-dixit-b1a07a39/
>
> ORCID ID: https://orcid.org/0000-0003-4015-2941
>
> http://scholar.google.co.in/citations?user=X876BKcAAAAJ&hl=en&oi=sra
>
> P Please consider the environment before printing this e-mail
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Dec 14 2020 - 04:30:02 PST
Custom Search