[AMBER] FW: AMBER Digest, Vol 3348, Issue 1

From: Markus Metz <mmetz.omeros.com>
Date: Thu, 29 Apr 2021 19:09:23 +0000

Dear Dave:
Thank you very much for your message.
Getting the new amber version is not an option.
However, would it be possible to get ambertool18 or 19?
Best,
Markus


From: amber-request.ambermd.org <amber-request.ambermd.org>
Sent: Thursday, April 29, 2021 12:00 PM
To: amber.ambermd.org
Subject: AMBER Digest, Vol 3348, Issue 1

Send AMBER mailing list submissions to
amber.ambermd.org<mailto:amber.ambermd.org>

To subscribe or unsubscribe via the World Wide Web, visit
http://lists.ambermd.org/mailman/listinfo/amber<http://lists.ambermd.org/mailman/listinfo/amber>
or, via email, send a message with subject or body 'help' to
amber-request.ambermd.org<mailto:amber-request.ambermd.org>

You can reach the person managing the list at
amber-owner.ambermd.org<mailto:amber-owner.ambermd.org>

When replying, please edit your Subject line so it is more specific
than "Re: Contents of AMBER digest..."


AMBER Mailing List Digest

Today's Topics:

1. CPPTRAJ configuration is failing (Markus Metz)
2. Re: CPPTRAJ configuration is failing (David A Case)
3. Trouble Creating a pdb File (Khong, Megan)
4. ??: [ambermd.org??]Re: ??: [ambermd.org??]Re: Max relative
errors questions in parallel and cuda.serial test (???)
5. Questions about REMD extracting temperature trajectory (Ojet)
6. Re: Trouble Creating a pdb File (Gustaf Olsson)
7. Re: Questions about REMD extracting temperature trajectory
(Charo del Genio)
8. Re: Problems to install AmberTools20 (Miguel Garav?s)
9. Re: ??: [ambermd.org??]Re: ??: [ambermd.org??]Re: Max
relative errors questions in parallel and cuda.serial test
(David A Case)
10. Re: Problems to install AmberTools20 (David A Case)
11. Adventures installing AMBER 20 in Centos 7 (Azurmendi, Hugo)
12. cpptraj energy using 2 masks (German P. Barletta)
13. Re: Adventures installing AMBER 20 in Centos 7 (David A Case)


----------------------------------------------------------------------

Message: 1
Date: Wed, 28 Apr 2021 19:04:56 +0000
From: Markus Metz <mmetz.omeros.com<mailto:mmetz.omeros.com>>
Subject: [AMBER] CPPTRAJ configuration is failing
To: "amber.ambermd.org<mailto:amber.ambermd.org>" <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID: <72a0cb1919e6435dbc798e925db44225.omeros.com<mailto:72a0cb1919e6435dbc798e925db44225.omeros.com>>
Content-Type: text/plain; charset=WINDOWS-1252

Hello all:
I am trying to compile a parallel version of Amber.
I first downloaded amber18 and ambertools20 and performed the serial compilation.
Next, I wanted to compile the parallel version.

I went to $AMBERHOME/AmberTools/src and issued ./configure_mpich gnu.
That went fine only some warnings.

Next I went back to $AMBERHOME and issued ./configure -mpi gnu.
Here is where I get my error. Please see the output:

Checking for updates...
Checking for available patches online. This may take a few seconds...


Available AmberTools 20 patches:

No patches available

Available Amber 20 patches:

update.1, update.2, update.3, update.4, update.5, update.6, update.7, update.8, update.9, update.10,update.11
There are patches available. Do you want to apply them now? [y/N] (Recommended Y)
N
NOT updating your tree and continuing anyway.


AMBER_PREFIX=/usr/local/amber/amber18/
AMBER_SOURCE=/usr/local/amber/amber18/

Using the AmberTools miniconda installation in /usr/local/amber/amber18//miniconda
version 3.8.5

Obtaining the gnu compiler suite versions, e.g.:
gcc -v
The C version is 4.8.5
The Fortran version is 9.3.0

Testing the gcc compiler:
gcc -fPIC -D_FILE_OFFSET_BITS=64 -D_LARGEFILE_SOURCE -o testp testp.c
OK

Testing the g++ compiler:
g++ -fPIC -o testp testp.cpp
OK

Testing the gfortran compiler:
gfortran -fPIC -O0 -o testp testp.f
OK

Testing mixed C/Fortran compilation:
gcc -fPIC -D_FILE_OFFSET_BITS=64 -D_LARGEFILE_SOURCE -c -o testp.c.o testp.c
gfortran -fPIC -O0 -c -o testp.f.o testp.f
gcc -fPIC -D_FILE_OFFSET_BITS=64 -D_LARGEFILE_SOURCE -o testp testp.c.o testp.f.o -lgfortran -w
OK

Testing pointer size:
gcc -fPIC -D_FILE_OFFSET_BITS=64 -D_LARGEFILE_SOURCE -o test_pointer_size test_pointer_size.c
Detected 64 bit operating system.

Testing flex: OK

Testing bison: OK

Checking NetCDF...
Using bundled NetCDF library.
Using existing NetCDF in '/usr/local/amber/amber18/'

Checking for zlib: OK

Checking for libbz2: OK

Configuring fftw-3.3 (may be time-consuming)...OK
Compiling the FFTW3 interface (may be time-consuming)...OK

Configuring boost (may be time-consuming)...OK
Compiling boost (may be time-consuming)...OK

Configuring CPPTRAJ...
CPPTRAJ configure failed. Check '/usr/local/amber/amber18/AmberTools/src/cpptraj_config.log' for details.
Configure failed due to the errors above!
(amber) [amber18]$ more /usr/local/amber/amber18/AmberTools/src/cpptraj_config.log
fftw3 enabled.
fftw3 specified: /usr/local/amber/amber18/
sanderlib disabled.
netcdf specified: /usr/local/amber/amber18/
blas specified: /usr/local/amber/amber18/
lapack specified: /usr/local/amber/amber18/
arpack specified: /usr/local/amber/amber18/
Testing basic C++11 support: Not present
Testing fix for STDIO/MPI ordering:
Error: Test compile failed: mpicxx -Wall -O3 -fPIC -DMPICH_IGNORE_CXX_SEEK -o testp testp.cpp
Error: Check the output below for error messages:
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: warning: librt.so.1, needed by /usr/local/amber/amber18/lib/libmpicxx.so, not found (try using -rpath or -rpath-link)
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: warning: libpthread.so.0, needed by /usr/local/amber/amber18/lib/libmpicxx.so, not found (try using -rpath or -rpath-link)
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: warning: libgfortran.so.5, needed by /usr/local/amber/amber18/lib/libmpi.so, not found (try using -rpath or -rpath-link)
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `aio_error64.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `pthread_mutexattr_setpshared.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `aio_suspend64.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `clock_getres.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `pthread_mutexattr_init.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `pthread_mutexattr_destroy.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `aio_write64.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `clock_gettime.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `aio_return64.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `pthread_create.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `pthread_mutexattr_settype.GLIBC_2.2.5'
/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: /usr/local/amber/amber18/lib/libmpi.so: undefined reference to `aio_read64.GLIBC_2.2.5'

I confirmed that librt.so.1, libpthread.so.0 and libgfortran.so.5 are present on my system and are also accessible through the LD_LIBRARY_PATH variable.

I am not sure how to proceed. Any pointers from you are greatly appreciated?

Best wishes,
Markus




------------------------------

Message: 2
Date: Wed, 28 Apr 2021 16:24:51 -0400
From: David A Case <david.case.rutgers.edu<mailto:david.case.rutgers.edu>>
Subject: Re: [AMBER] CPPTRAJ configuration is failing
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<20210428202451.ynet2ydjx2ptzio5.gsappa209l01.rad.rutgers.edu<mailto:20210428202451.ynet2ydjx2ptzio5.gsappa209l01.rad.rutgers.edu>>
Content-Type: text/plain; charset=us-ascii; format=flowed

On Wed, Apr 28, 2021, Markus Metz wrote:

>I am trying to compile a parallel version of Amber.
>I first downloaded amber18 and ambertools20 and performed the serial
>compilation.

I think you are on the wrong pathway. Amber18 is only designed to be
compiled in conjunction with AmberTools18 or AmberTools19. The connections
between Amber and AmberTools change every two years.

(You can install AmberTools20 in a separate folder and use those
executables, setting AMBERHOME to the amber20 install folder. Then use the
pmemd from Amber18 when you need pmemd; pmemd doesn't make use of the
AMBERHOME variable. But this *use* is very different from building Amber18
in the first place.)

The straightforward path, if possible, is to upgrade to Amber20. But see
some more comments below.

>The C version is 4.8.5
>The Fortran version is 9.3.0

This combination is unlikely to work, even in the best circumstances. Try
to get C and Fortran versions at the same level (best to upgrade C to 9.3).

>Error: Test compile failed: mpicxx -Wall -O3 -fPIC -DMPICH_IGNORE_CXX_SEEK -o testp testp.cpp
>Error: Check the output below for error messages:
>/anaconda3/envs/amber/bin/x86_64-conda_cos6-linux-gnu-ld: warning: librt.so.1, needed by /usr/local/amber/amber18/lib/libmpicxx.so, not found (try using -rpath or -rpath-link)

Oooh...looks like you also have conda environment with some sort of "amber"
folder in it! I don't know what this is. But you will most likely need to
deactivate any such conda evironment for building Amber.

...sorry to not be more encouraging, but I think you need to start afresh.

....dac




------------------------------

Message: 3
Date: Wed, 28 Apr 2021 21:14:36 +0000
From: "Khong, Megan" <MEGAN.KHONG.UCDENVER.EDU<mailto:MEGAN.KHONG.UCDENVER.EDU>>
Subject: [AMBER] Trouble Creating a pdb File
To: "amber.ambermd.org<mailto:amber.ambermd.org>" <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<CY4PR05MB354380E00726B8525E9D320B93409.CY4PR05MB3543.namprd05.prod.outlook.com<mailto:CY4PR05MB354380E00726B8525E9D320B93409.CY4PR05MB3543.namprd05.prod.outlook.com>>

Content-Type: text/plain; charset="iso-8859-1"

Hello, I am trying to create a pdb file using wat_md.rst as the input file. I am using AMBER 18. I have tried the command

ambpdb -p wat.prmtop < wat_md.rst

and I have been getting an error that says, "Error: could not read restart atoms/time". Has anyone encounter this problem before? Does anyone have any suggestions on how I can check the validity of my file?

Thank you.

-Megan Khong
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wat.prmtop
Type: application/octet-stream
Size: 214476 bytes
Desc: wat.prmtop
Url : http://lists.ambermd.org/mailman/private/amber/attachments/20210428/178ee20a/attachment-0002.obj<http://lists.ambermd.org/mailman/private/amber/attachments/20210428/178ee20a/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wat_md.rst
Type: application/octet-stream
Size: 63168 bytes
Desc: wat_md.rst
Url : http://lists.ambermd.org/mailman/private/amber/attachments/20210428/178ee20a/attachment-0003.obj<http://lists.ambermd.org/mailman/private/amber/attachments/20210428/178ee20a/attachment-0003.obj>

------------------------------

Message: 4
Date: Thu, 29 Apr 2021 00:57:46 +0000
From: ??? <guqin.shi.qilu-pharma.com<mailto:guqin.shi.qilu-pharma.com>>
Subject: [AMBER] ??: [ambermd.org??]Re: ??: [ambermd.org??]Re: Max
relative errors questions in parallel and cuda.serial test
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID: <1a10203f9a3f4624be2199f45243fad2.qilu-pharma.com<mailto:1a10203f9a3f4624be2199f45243fad2.qilu-pharma.com>>
Content-Type: text/plain; charset="utf-8"

Dear David,

Thanks for your follow-up. If the errors are with GB calculation, then that's quite serious because the most common post-processing tasks we will do is MMPB(GB)SA.

Unfortunately just a few days ago the server suffered a crash when we tried to re-allocate the disk space. I lost the amber20_installation directories. I will re-install everything and report the results.

David, what will be the potential way to solve this DPFP errors? Re-compiling with different ways? Install on a different workstation/cluster.....?

Thanks a lot,
Guqin

-----????-----
???: David A Case [mailto:david.case.rutgers.edu]
????: 2021?4?28? 19:58
???: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
??: [ambermd.org??]Re: [AMBER] ??: [ambermd.org??]Re: Max relative errors questions in parallel and cuda.serial test

On Wed, Apr 28, 2021, ??? wrote:
>
>Thought I need to bring this up again...I encountered large CUDA test
>error with DPFP precision mode. The maximum relative errors are usually
>at the level of e-01.

The errors you see are quite serious. But they are all in generalized Born calculations, so you are probably safe in running explicit solvent simulations.

I think you mentioned that you didn't see errors in SPFP. Could you go to $AMBERHOME/test and run "make test.cuda.serial.SPFP"? Do you see no errors there?

....thanks....dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org>
http://lists.ambermd.org/mailman/listinfo/amber<http://lists.ambermd.org/mailman/listinfo/amber>


***********????*************

??????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????

***********Business Email Disclaimer**************

This e-mail and any attachments are meant for the intended recipient only and may contain information belonging to Qilu Pharma that is privileged, confidential, proprietary, and/or otherwise protected or prohibited from disclosure. If you are not the correct recipient or received this e-mail erroneously, please inform the sender immediately and delete this mail from your system. Qilu Pharma state no liability for any damage to your system and data caused by this email. Unless this email is related to the business with the company, otherwise any views or opinions presented in this email are solely from the sender. Qilu Pharma has the right to monitor the sending and receiving of the e-mail.

------------------------------

Message: 5
Date: Thu, 29 Apr 2021 14:26:12 +0800 (CST)
From: Ojet <ojet0501.163.com<mailto:ojet0501.163.com>>
Subject: [AMBER] Questions about REMD extracting temperature
trajectory
To: amber.ambermd.org<mailto:amber.ambermd.org>
Message-ID: <71de993.442a.1791c4e6cde.Coremail.ojet0501.163.com<mailto:71de993.442a.1791c4e6cde.Coremail.ojet0501.163.com>>
Content-Type: text/plain; charset=GBK

Dear ambers expertsI did REMD according to the tutorial. After the simulation is over, an error is reported when using ? extract_300k_traj.x? script to extract the temperature trajectory:
remdtraj requested in trajin but this NETCDF trajectory does not contain temperatures! Trajectory will be processed as replica traj.



The ?extract_300k_traj.x? script is as follows:

#!/bin/csh

ptraj 5EK4_1.prmtop << EOF

trajin remd.mdcrd.001 remdtraj remdtrajtemp 300

trajout remd.300K.mdcrd nobox

go

EOF



The following are the key input files when I perform REMD:

&cntrl

irest=0, ntx=1,

nstlim=500, dt=0.002,

irest=0, ntt=3, gamma_ln=1.0,

temp0=XXXXX, ig=RANDOM_NUMBER,

ntc=2, ntf=2, nscm=1000,

ntb=0, igb=5,

cut=999.0, rgbmax=999.0,

ntpr=100, ntwx=1000, ntwr=10000,

nmropt=1,

numexchg=50000,

ntr=1,

restraint_wt=50.0,

restraintmask=':1-347,374-390.CA',

/

&wt TYPE='END'

/

DISANG=5EK4_chir.dat



REMD simulations are run for a temperature range between 290 - 472K. The starting structure of the 30 replicas of REMD is the 5EKC protein with 30 different structures. 50 ns for each replica.



remd.groupfile:

-O -rem 1 -remlog rem.log -i remd.mdin.001 -o remd.mdout.001 -c 5EK4_1.inpcrd -r remd.rst.001 -x remd.mdcrd.001 -inf remd.mdinfo.001 -p 5EK4_1.prmtop -ref 5EK4_1.inpcrd

-O -rem 1 -remlog rem.log -i remd.mdin.002 -o remd.mdout.002 -c 5EK4_2.inpcrd -r remd.rst.002 -x remd.mdcrd.002 -inf remd.mdinfo.002 -p 5EK4_1.prmtop -ref 5EK4_2.inpcrd

??

-O -rem 1 -remlog rem.log -i remd.mdin.028 -o remd.mdout.028 -c 5EK4_28.inpcrd -r remd.rst.028 -x remd.mdcrd.028 -inf remd.mdinfo.028 -p 5EK4_1.prmtop -ref 5EK4_28.inpcrd

-O -rem 1 -remlog rem.log -i remd.mdin.029 -o remd.mdout.029 -c 5EK4_29.inpcrd -r remd.rst.029 -x remd.mdcrd.029 -inf remd.mdinfo.029 -p 5EK4_1.prmtop -ref 5EK4_29.inpcrd

-O -rem 1 -remlog rem.log -i remd.mdin.030 -o remd.mdout.030 -c 5EK4_30.inpcrd -r remd.rst.030 -x remd.mdcrd.030 -inf remd.mdinfo.030 -p 5EK4_1.prmtop -ref 5EK4_30.inpcrd



The command to run REMD is:

srun -n 30 pmemd.cuda.MPI -ng 30 -groupfile remd.groupfile



The output file is as follows:

rem.log:

......

# exchange 50000

1 -1.00 483.93 -4919.86 472.00 472.00 0.00 -1

2 1.01 425.03 -5615.09 421.90 428.80 0.28 -1

3 -1.00 357.92 -6665.98 363.30 363.30 0.26 -1

4 -1.00 431.49 -5213.03 443.00 443.00 0.29 -1

5 0.99 310.04 -7672.94 316.80 311.30 0.25 -1

6 -1.00 381.14 -6215.04 395.00 395.00 0.27 -1

7 0.99 457.53 -5069.34 457.60 450.20 0.29 -1

8 0.99 308.95 -7850.23 305.80 300.50 0.24 -1

9 1.01 310.87 -7623.60 311.30 316.80 0.25 -1

10 -1.00 435.48 -5500.29 435.90 435.90 0.29 -1

11 0.99 422.83 -5597.49 428.80 421.90 0.28 -1

12 -1.00 320.61 -7401.61 322.30 322.30 0.25 -1

13 1.01 448.75 -5117.91 450.20 457.60 0.28 -1

14 1.01 300.06 -7902.66 300.50 305.80 0.25 -1

15 -1.00 407.83 -6045.13 401.60 401.60 0.27 -1

16 -1.00 411.97 -5993.03 408.30 408.30 0.28 -1

17 -1.00 342.54 -7093.49 345.30 345.30 0.26 -1

18 -1.00 374.51 -6618.57 369.50 369.50 0.27 -1

19 -1.00 406.04 -5767.74 415.00 415.00 0.28 -1

20 -1.00 347.37 -7103.56 339.40 339.40 0.26 -1

21 -1.00 328.74 -7117.47 333.60 333.60 0.25 -1

22 -1.00 354.91 -6798.13 357.20 357.20 0.26 -1

23 -1.00 377.97 -6340.68 382.10 382.10 0.28 -1

24 0.99 294.19 -7967.01 295.20 290.00 0.25 -1

25 1.01 291.39 -7976.23 290.00 295.20 0.25 -1

26 -1.00 326.46 -7318.99 327.90 327.90 0.25 -1

27 -1.00 380.62 -6394.33 375.70 375.70 0.26 -1

28 -1.00 384.89 -6096.66 388.50 388.50 0.28 -1

29 -1.00 349.34 -6909.41 351.20 351.20 0.26 -1

30 -1.00 468.11 -5049.65 465.00 465.00 0.32 -1



remd.mdout.010 :

......

NSTEP = 25000000 TIME(PS) = 50000.000 TEMP(K) = 433.81 PRESS = 0.0

Etot = 1429.7036 EKtot = 6800.4443 EPtot = -5370.7407

BOND = 1641.3302 ANGLE = 4534.2074 DIHED = 5377.9737

1-4 NB = 1573.3763 1-4 EEL = 17498.4438 VDWAALS = -3220.7094

EELEC = -28666.4125 EGB = -4444.1999 RESTRAINT = 335.2497

EAMBER (non-restraint) = -5705.9904

TEMP0 = 435.9000 REPNUM = 10 EXCHANGE# = 50000

------------------------------------------------------------------------------

NMR restraints: Bond = 0.000 Angle = 0.000 Torsion = 18.340












------------------------------

Message: 6
Date: Thu, 29 Apr 2021 06:40:35 +0000
From: Gustaf Olsson <gustaf.olsson.lnu.se<mailto:gustaf.olsson.lnu.se>>
Subject: Re: [AMBER] Trouble Creating a pdb File
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID: <FAEE80C5-629A-4D72-A76E-56D4F22DC2BA.lnu.se<mailto:FAEE80C5-629A-4D72-A76E-56D4F22DC2BA.lnu.se>>
Content-Type: text/plain; charset="us-ascii"

As the RST file is a binary file (NetCDF, ncrst) file, the suggested command for ambpdb may not work as intended.

A quick search revealed this issue being adress in an earlier discussion:

The file redirect '<' only works with ASCII restarts - if this is a
NetCDF restart use the '-c' flag instead:

$AMBERHOME/bin/ambpdb -p cram.prmtop -c min_qmmm.rst > min_qmmm.rst.pdb

So for you files this would mean:

ambpdb -p wat.prmtop -c wat_md.rst > pdbfile.pdb

Producing a PDB file.


Best regards
// Gustaf

On 28 Apr 2021, at 23:14, Khong, Megan <MEGAN.KHONG.UCDENVER.EDU<mailto:MEGAN.KHONG.UCDENVER.EDU<mailto:MEGAN.KHONG.UCDENVER.EDU%3cmailto:MEGAN.KHONG.UCDENVER.EDU>>> wrote:

Hello, I am trying to create a pdb file using wat_md.rst as the input file. I am using AMBER 18. I have tried the command

ambpdb -p wat.prmtop < wat_md.rst

and I have been getting an error that says, "Error: could not read restart atoms/time". Has anyone encounter this problem before? Does anyone have any suggestions on how I can check the validity of my file?

Thank you.

-Megan Khong
<wat.prmtop><wat_md.rst>_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org<mailto:AMBER.ambermd.org%3cmailto:AMBER.ambermd.org>>
http://lists.ambermd.org/mailman/listinfo/amber<http://lists.ambermd.org/mailman/listinfo/amber>

-------------- next part --------------
A non-text attachment was scrubbed...
Name: pdbfile.pdb
Type: chemical/x-pdb
Size: 106352 bytes
Desc: pdbfile.pdb
Url : http://lists.ambermd.org/mailman/private/amber/attachments/20210429/28f2f652/attachment-0001.pdb<http://lists.ambermd.org/mailman/private/amber/attachments/20210429/28f2f652/attachment-0001.pdb>

------------------------------

Message: 7
Date: Thu, 29 Apr 2021 09:45:46 +0100
From: Charo del Genio <the.paraw.gmail.com<mailto:the.paraw.gmail.com>>
Subject: Re: [AMBER] Questions about REMD extracting temperature
trajectory
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID: <97c60f7b-edd8-0eae-dc06-bed0f3103493.gmail.com<mailto:97c60f7b-edd8-0eae-dc06-bed0f3103493.gmail.com>>
Content-Type: text/plain; charset=utf-8; format=flowed

On 29/04/2021 07:26, Ojet wrote:
> Dear ambers expertsI did REMD according to the tutorial. After the simulation is over, an error is reported when using ? extract_300k_traj.x? script to extract the temperature trajectory:
> remdtraj requested in trajin but this NETCDF trajectory does not contain temperatures! Trajectory will be processed as replica traj.
> The ?extract_300k_traj.x? script is as follows:
> #!/bin/csh
> ptraj 5EK4_1.prmtop << EOF
> trajin remd.mdcrd.001 remdtraj remdtrajtemp 300
> trajout remd.300K.mdcrd nobox
> go
> EOF

Hi, I already replied before. Please try using cpptraj rather than ptraj. In fact, I thought ptraj wasn't even any longer within the distribution...


Cheers,

Charo




--
Dr. Charo I. del Genio
Senior Lecturer in Statistical Physics
Applied Mathematics Research Centre (AMRC)
Design Hub
Coventry University Technology Park
Coventry CV1 5FB
UK
https://charodelgenio.weebly.com<https://charodelgenio.weebly.com>
------------------------------
Message: 8
Date: Thu, 29 Apr 2021 12:43:12 +0200
From: Miguel Garav?s <mgaravis.gmail.com<mailto:mgaravis.gmail.com>>
Subject: Re: [AMBER] Problems to install AmberTools20
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<CAJVSO568D8BFDxV0pfLnDy1jJuZ2-wE=V=MyrY8UtOvpV9yTFQ.mail.gmail.com<mailto:CAJVSO568D8BFDxV0pfLnDy1jJuZ2-wE=V=MyrY8UtOvpV9yTFQ.mail.gmail.com>>
Content-Type: text/plain; charset="UTF-8"
Thanks again.
I did a deep cleaning and tried to install it again. The result is still
the same.
In case it helps, I paste here some of the lines concerning existing errors:
-- Failed to find Perl modules: Chemistry::Mol (missing:
PERLMODULES_CHEMISTRY_MOL_MODULE)
-- Could NOT find Boost (missing: Boost_INCLUDE_DIR thread system
program_options iostreams regex timer chrono filesystem graph)
-- Cannot find PLUMED. You will still be able to load it at runtime. If
you want to link it at build time, set PLUMED_ROOT to where you installed
it.
CMake Error at cmake/LibraryUtils.cmake:75 (message):
Could not determine whether
"/Applications/Xcode.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX11.1.sdk/usr/lib/libbz2.tbd"
is a static or shared library, it does not have a known suffix.
The OSX version I have is MacOS Catalina 10.15.7 (19H15)
The version of XCode is: 12.4 (12D4e)
And the clang version is clang-1200.0.32.29
MGC
El mi?, 28 abr 2021 a las 18:41, David A Case (<david.case.rutgers.edu<mailto:david.case.rutgers.edu>>)
escribi?:
> On Wed, Apr 28, 2021, Miguel Garav?s wrote:
> >
> >I can confirm that update.11 is applied. It appears in the list after
> >typing the command you indicated.
>
> OK, we have to search for something else.
>
> First, I think it's worth trying a "clean" build: re-download the
> AmberTools20.tar.bz2 file (unless you still have it on your computer),
> un-tar and try a "run_cmake" job again. It's possible that some cached
> file
> is interfering with what you did before.
>
> If that doesn't help, please let us know what version of OSX you have, and
> which version of Xcode (open /Applications/Xcode.app and look at the "about
> Xcode menu item.) Also, type "clang --version" at a terminal prompt, and
> report the result.
>
> ...thanks...dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org<mailto:AMBER.ambermd.org>
> http://lists.ambermd.org/mailman/listinfo/amber<http://lists.ambermd.org/mailman/listinfo/amber>
>
------------------------------
Message: 9
Date: Thu, 29 Apr 2021 08:09:35 -0400
From: David A Case <david.case.rutgers.edu<mailto:david.case.rutgers.edu>>
Subject: Re: [AMBER] ??: [ambermd.org??]Re: ??: [ambermd.org??]Re:
Max relative errors questions in parallel and cuda.serial test
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<20210429120935.guzxnqmkym4ff4hn.rulx-presler-l.rad.rutgers.edu<mailto:20210429120935.guzxnqmkym4ff4hn.rulx-presler-l.rad.rutgers.edu>>
Content-Type: text/plain; charset=utf-8; format=flowed
On Thu, Apr 29, 2021, ??? wrote:
>
>Thanks for your follow-up. If the errors are with GB calculation, then
>that's quite serious because the most common post-processing tasks we will
>do is MMPB(GB)SA.
Do you use pmemd.cuda for such tasks? Is it really significantly faster
than pmemd or sander on a CPU? I ask because getting single frame energies
involves a lot of initialization tasks that reduce the speed advantage of
the GB calculation itself.
>
>David, what will be the potential way to solve this DPFP
>errors? Re-compiling with different ways? Install on a different
>workstation/cluster.....?
Both of the above are good ideas. It's certainly possible that you are
seeing a hardware problem, so trying a different GPU might help.
Changing compilers or OS might also help identify the problem. Right now,
we don't have any way to reproduce the problem at our end.
...regards....dac
------------------------------
Message: 10
Date: Thu, 29 Apr 2021 08:21:53 -0400
From: David A Case <david.case.rutgers.edu<mailto:david.case.rutgers.edu>>
Subject: Re: [AMBER] Problems to install AmberTools20
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<20210429122153.4nhifvuntn5vozjk.rulx-presler-l.rad.rutgers.edu<mailto:20210429122153.4nhifvuntn5vozjk.rulx-presler-l.rad.rutgers.edu>>
Content-Type: text/plain; charset=iso-8859-1; format=flowed
On Thu, Apr 29, 2021, Miguel Garav?s wrote:
>
>I did a deep cleaning and tried to install it again. The result is still
>the same.
>
>CMake Error at cmake/LibraryUtils.cmake:75 (message):
> Could not determine whether
>"/Applications/Xcode.app/Contents/Developer/Platforms/MacOSX.platform/Developer/SDKs/MacOSX11.1.sdk/usr/lib/libbz2.tbd"
> is a static or shared library, it does not have a known suffix.
>
>The OSX version I have is MacOS Catalina 10.15.7 (19H15)
>The version of XCode is: 12.4 (12D4e)
>And the clang version is clang-1200.0.32.29
I don't have any more suggestions for now. AmberTools21 will be released
today or tomorrow, so you might try that when it is available.
If you have access (via a friend?) to any other OSX machine, you might try
an installation there. If it works, one could try to see what difference
there is between the two computers.
If time is critical, you could try using the configure option in place of
cmake. See Section 2.5 of the Amber 2020 Reference Manual. Or, you could
install the conda package if that meets your needs.
It seems that very few Amber developers (besides me!) use OSX now, but I
think a lot of users do. We haven't previously seen any such reports since
update.11 was released last October. But maybe someone on the list will
chime in with ideas here.
....dac
------------------------------
Message: 11
Date: Thu, 29 Apr 2021 15:59:40 +0000
From: "Azurmendi, Hugo" <Hugo.Azurmendi.fda.hhs.gov<mailto:Hugo.Azurmendi.fda.hhs.gov>>
Subject: [AMBER] Adventures installing AMBER 20 in Centos 7
To: "amber.ambermd.org<mailto:amber.ambermd.org>" <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID: <D63291C4-E657-4E5F-A13A-59603C5398BE.fda.gov<mailto:D63291C4-E657-4E5F-A13A-59603C5398BE.fda.gov>>
Content-Type: text/plain; charset="utf-8"
Hello there.
I just finished installing AMBER 20 on a new Linux and found a few problems that I think developers should be aware of. I did a fresh install of Centos 7 installing the recommended packages for centos with yum. The problems were:
1. For parallel I followed the instructions for Fedora to install openmpi as follows (not sure all these are required, but it worked):
? yum install dnf
? dnf install openmpi openmpi-devel environment-modules
? source /etc/profile.d/modules.sh (to load modules)
? Add line 'module load mpi/openmpi-x86_64' to .bashrc/.cshrc.
At this point I could start ?make test.parallel?, however it would stall (without messages) after a few tests. I finally could get it to go all the way through by taking out (renaming the corresponding $AMBERHOME/test/testname) the tests for dhfr and rism3d. I haven?t investigated further this issue.
1. The initial installation and test I did it in bash, deciding to switch to csh later. I sourced to amber.csh (located in the amber20 folder) in my .cshrc but didn?t work:
* It warned that the shell was not csh, which it was. I fixed replacing ( "$0" != '-tcsh' && "$0" != '-csh' ) for ( "$0" !~ '*csh' ) in the ?if? line of amber.csh.
* AMBERHOME was set to the path of the .cshrc instead of the path to amber.csh. I commented the setenv line in amber.csh and set AMBERHOME in my .cshrc in the old fashion.
Hope that helps to improve future releases. Looking forward to comments if I did something wrong.
Hugo
------------------------------
Message: 12
Date: Thu, 29 Apr 2021 12:01:07 -0600
From: "German P. Barletta" <pbarletta.gmail.com<mailto:pbarletta.gmail.com>>
Subject: [AMBER] cpptraj energy using 2 masks
To: ?amber.ambermd.org<mailto:?amber.ambermd.org>? <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<CAJo28majJB5FPaOQ6Uikved8zSphPoGCBCJvosNw0oAbYMhygg.mail.gmail.com<mailto:CAJo28majJB5FPaOQ6Uikved8zSphPoGCBCJvosNw0oAbYMhygg.mail.gmail.com>>
Content-Type: text/plain; charset="UTF-8"
Dear amber users and developers,
I'd like to calculate the nonbonded (LJ + Coulombic) energy terms for the
interaction between a ligand (residue 164) and a protein (residues 1-163)
but if I use the mask ":1-164", cpptraj will also include the interactions
of the protein and the ligand with themselves.
gromacs tutorials
<http://www.mdtutorials.com/gmx/complex/09_analysis.html<http://www.mdtutorials.com/gmx/complex/09_analysis.html>> show
a convenient way to do this, but how can I do the same with cpptraj?
Thanks!
------------------------------
Message: 13
Date: Thu, 29 Apr 2021 14:29:49 -0400
From: David A Case <david.case.rutgers.edu<mailto:david.case.rutgers.edu>>
Subject: Re: [AMBER] Adventures installing AMBER 20 in Centos 7
To: AMBER Mailing List <amber.ambermd.org<mailto:amber.ambermd.org>>
Message-ID:
<20210429182949.fea2p6kkmhu7opjr.ufa33kbb210dl.rad.rutgers.edu<mailto:20210429182949.fea2p6kkmhu7opjr.ufa33kbb210dl.rad.rutgers.edu>>
Content-Type: text/plain; charset=us-ascii; format=flowed
On Thu, Apr 29, 2021, Azurmendi, Hugo wrote:
>I just finished installing AMBER 20 on a new Linux and found a few problems
>that I think developers should be aware of. I did a fresh install of Centos
>7 installing the recommended packages for centos with yum. The problems
>were:
>
> 1. For parallel I followed the instructions for Fedora
Can you say a little more about exactly what you did? Did you make any
changes to run_cmake other than setting -DMPI_TRUE? What value did you set
for DO_PARALLEL, and what was the exact command you used to run the tests?
* It warned that the shell was not csh
I can believe that the amber.csh may not work in all cases, but it would be
more helpful if you could quote the exact error message, since that would
tell us what the value of the "$0" variable is in your case.
...thx...dac
------------------------------
_______________________________________________
AMBER mailing list
AMBER.ambermd.org<mailto:AMBER.ambermd.org>
http://lists.ambermd.org/mailman/listinfo/amber<http://lists.ambermd.org/mailman/listinfo/amber>
End of AMBER Digest, Vol 3348, Issue 1
**************************************
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Apr 29 2021 - 12:30:02 PDT
Custom Search