Thank you for the help.
Converting to xyz format also does not work:
$ babel -i pdb EPE.pdb -o xyz EPE.xyz
33
EPE.pdb
N -0.82000 -0.33000 0.25800
C -1.66900 0.55700 1.06100
C -1.69300 0.07200 2.51100
N -0.34200 0.14700 3.07900
C 0.50700 -0.74000 2.27600
C 0.53100 -0.25500 0.82600
C -0.41600 -0.44100 4.42300
C 0.96600 -0.39500 5.07600
O 1.40200 0.96100 5.17400
C -0.74600 0.25800 -1.08500
C 0.12800 -0.62100 -1.98000
S 0.21900 0.10400 -3.64000
O -1.01900 -0.00100 -4.32700
O 0.98600 1.30100 -3.64600
O 1.11000 -0.87000 -4.39700
H -2.68300 0.55100 0.65900
H -1.27100 1.57100 1.02400
H -2.04400 -0.95900 2.54300
H -2.36700 0.70200 3.09300
H 1.52000 -0.73400 2.67700
H 0.10900 -1.75400 2.31300
H 1.20400 -0.88500 0.24400
H 0.88100 0.77500 0.79400
H -0.74800 -1.47600 4.34800
H -1.12300 0.12400 5.02900
H 1.67300 -0.96100 4.47000
H 0.91100 -0.83200 6.07300
H 2.27500 0.94500 5.59000
H -1.74800 0.32300 -1.50900
H -0.31300 1.25600 -1.02200
H 1.13100 -0.68700 -1.55700
H -0.30400 -1.62000 -2.04400
H 1.19400 -0.52900 -5.29800
==============================
*** Open Babel Error in OpenAndSetFormat
Cannot open EPE.xyz
1 molecule converted
Then I tried making a new folder, copy the pdb, and run the babel again, still the same, not work. So I think I already have the permission for the folders, but the babel does not work.
Then I used the antechamber. Even though the “no version information available” appears in the last two lines, an EPE.mol2 file was generated!
$ antechamber -fi pdb -fo mol2 -i EPE.pdb -o EPE.mol2
Welcome to antechamber 21.0: molecular input file processor.
acdoctor mode is on: check and diagnose problems in the input file.
The atom type is set to gaff; the options available to the -at flag are
gaff, gaff2, amber, bcc, and sybyl.
-- Check Format for pdb File --
Status: pass
-- Check Unusual Elements --
Status: pass
-- Check Open Valences --
Status: pass
-- Check Geometry --
for those bonded
for those not bonded
Status: pass
-- Check Weird Bonds --
Status: pass
-- Check Number of Units --
Status: pass
acdoctor mode has completed checking the input file.
/bin/bash: /home/lanselibai/miniconda3/envs/AmberTools21/lib/libtinfo.so.6: no version information available (required by /bin/bash)
/bin/bash: /home/lanselibai/miniconda3/envs/AmberTools21/lib/libtinfo.so.6: no version information available (required by /bin/bash)
So I can now try to use the generated mol2 file for acpype, but got the “mopac.sh: not found” error, which I mentioned in the 1st thread.:
$ acpype -i EPE.mol2
============================================================================
| ACPYPE: AnteChamber PYthon Parser interfacE v. 0 0 Rev: 0 (c) 2021 AWSdS |
============================================================================
==> ... charge set to 0
WARNING: Old version of antechamber. Strongly consider upgrading to AmberTools
WARNING: Setting mopac for antechamber
==> Executing Antechamber...
++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Welcome to antechamber 21.0: molecular input file processor.
acdoctor mode is on: check and diagnose problems in the input file.
The atom type is set to gaff; the options available to the -at flag are
gaff, gaff2, amber, bcc, and sybyl.
-- Check Format for mol2 File --
Status: pass
Info: Ignoring Mol2 record type (@<TRIPOS>SUBSTRUCTURE).
Info: Finished reading file (EPE.mol2); atoms read (33), bonds read (33).
Info: Determining atomic numbers from atomic symbols which are case sensitive.
-- Check Unusual Elements --
Status: pass
-- Check Open Valences --
Status: pass
-- Check Geometry --
for those bonded
for those not bonded
Status: pass
-- Check Weird Bonds --
Status: pass
-- Check Number of Units --
Status: pass
acdoctor mode has completed checking the input file.
Running: /home/lanselibai/miniconda3/envs/AmberTools21/bin/bondtype -j full -i ANTECHAMBER_BOND_TYPE.AC0 -o ANTECHAMBER_BOND_TYPE.AC -f ac
/bin/bash: /home/lanselibai/miniconda3/envs/AmberTools21/lib/libtinfo.so.6: no version information available (required by /bin/bash)
Running: /home/lanselibai/miniconda3/envs/AmberTools21/bin/atomtype -i ANTECHAMBER_AC.AC0 -o ANTECHAMBER_AC.AC -p gaff
/bin/bash: /home/lanselibai/miniconda3/envs/AmberTools21/lib/libtinfo.so.6: no version information available (required by /bin/bash)
Info: Total number of electrons: 128; net charge: 0
Running: /home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh
sh: 1: /home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh: not found
/home/lanselibai/miniconda3/envs/AmberTools21/bin/wrapped_progs/antechamber: Fatal Error!
Cannot properly run "/home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh".
++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
ERROR: Antechamber failed
++++++++++start_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
/bin/sh: 1: -i: not found
++++++++++end_quote+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
ERROR: Parmchk failed
ERROR: Tleap failed
==> Removing temporary files...
ACPYPE FAILED: [Errno 2] No such file or directory: 'EPE_AC.prmtop'
Total time of execution: less than a second
After looking at the 4th thread, I copy /home/lanselibai/miniconda3/envs/AmberTools21/amber.sh
to the last line of /home/lanselibai/.bashrc
.
Then I restart the WSL and run acpype -i EPE.mol2
, still the mopac
problem
Running: /home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh
sh: 1: /home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh: not found
/home/lanselibai/miniconda3/envs/AmberTools21/bin/wrapped_progs/antechamber: Fatal Error!
Cannot properly run "/home/lanselibai/miniconda3/envs/AmberTools21/bin/mopac.sh".