Example: pdb2adf transforms a PDB file to a QMMM input file

Download pdb2adf.run

#! /bin/sh

# This example shows how to use the utility pdb2adf,
# which creates an amsified ADF input file (ADF>=2020) from a PDB file,
# for a subsequent QM/MM calculation using ADF.

# -------------------------
# First create the PDB file
# -------------------------

cat << eor > chymotrypsin.pdb
HEADER    COMPLEX (SERINE PROTEASE/INHIBITOR)     12-MAR-97   1AFQ              
TITLE     CRYSTAL STRUCTURE OF BOVINE GAMMA-CHYMOTRYPSIN COMPLEXED              
TITLE    2 WITH A SYNTHETIC INHIBITOR                                           
REMARK
REMARK    Adaptation of original PDB file by M. Swart, March 2005
REMARK    only coordinates of GAMMA-CHYMOTRYPSIN are kept;
REMARK    rest has been deleted.
REMARK
ATOM      1  N   CYS A   1      13.717  20.021  22.754  1.00 13.46      PROA N  
ATOM      2  CA  CYS A   1      14.211  18.932  23.617  1.00 13.34      PROA C  
ATOM      3  C   CYS A   1      13.597  19.033  25.005  1.00 13.34      PROA C  
ATOM      4  O   CYS A   1      12.953  20.026  25.329  1.00 13.48      PROA O  
ATOM      5  CB  CYS A   1      15.734  19.018  23.753  1.00 13.44      PROA C  
ATOM      6  SG  CYS A   1      16.298  20.647  24.361  1.00 13.30      PROA S  
ATOM      7  N   GLY A   2      13.801  17.985  25.813  1.00 13.44      PROA N  
ATOM      8  CA  GLY A   2      13.369  17.952  27.214  1.00 13.65      PROA C  
ATOM      9  C   GLY A   2      11.904  18.088  27.631  1.00 13.87      PROA C  
ATOM     10  O   GLY A   2      11.669  18.375  28.799  1.00 13.63      PROA O  
ATOM     11  N   VAL A   3      10.947  17.887  26.732  1.00 14.18      PROA N  
ATOM     12  CA  VAL A   3       9.559  17.968  27.090  1.00 14.86      PROA C  
ATOM     13  C   VAL A   3       8.875  16.684  26.624  1.00 15.04      PROA C  
ATOM     14  O   VAL A   3       8.529  16.546  25.452  1.00 14.91      PROA O  
ATOM     15  CB  VAL A   3       8.861  19.211  26.437  1.00 15.00      PROA C  
ATOM     16  CG1 VAL A   3       7.403  19.299  26.880  1.00 15.08      PROA C  
ATOM     17  CG2 VAL A   3       9.585  20.486  26.805  1.00 15.27      PROA C  
ATOM     18  N   PRO A   4       8.754  15.691  27.519  1.00 15.39      PROA N  
ATOM     19  CA  PRO A   4       8.121  14.407  27.206  1.00 16.09      PROA C  
ATOM     20  C   PRO A   4       6.675  14.535  26.769  1.00 16.35      PROA C  
ATOM     21  O   PRO A   4       5.957  15.387  27.275  1.00 16.43      PROA O  
ATOM     22  CB  PRO A   4       8.219  13.635  28.527  1.00 15.91      PROA C  
ATOM     23  CG  PRO A   4       9.369  14.297  29.244  1.00 16.02      PROA C  
ATOM     24  CD  PRO A   4       9.166  15.742  28.928  1.00 15.56      PROA C  
ATOM     25  N   ALA A   5       6.262  13.690  25.827  1.00 16.87      PROA N  
ATOM     26  CA  ALA A   5       4.874  13.703  25.351  1.00 17.53      PROA C  
ATOM     27  C   ALA A   5       4.020  13.055  26.437  1.00 17.86      PROA C  
ATOM     28  O   ALA A   5       2.862  13.413  26.643  1.00 17.81      PROA O  
ATOM     29  CB  ALA A   5       4.740  12.936  24.027  1.00 17.34      PROA C  
ATOM     30  N   ILE A   6       4.615  12.104  27.143  1.00 18.37      PROA N  
ATOM     31  CA  ILE A   6       3.941  11.419  28.239  1.00 18.99      PROA C  
ATOM     32  C   ILE A   6       4.553  11.993  29.514  1.00 19.59      PROA C  
ATOM     33  O   ILE A   6       5.726  11.771  29.807  1.00 19.17      PROA O  
ATOM     34  CB  ILE A   6       4.190   9.909  28.190  1.00 19.29      PROA C  
ATOM     35  CG1 ILE A   6       3.631   9.335  26.886  1.00 19.66      PROA C  
ATOM     36  CG2 ILE A   6       3.552   9.232  29.399  1.00 19.38      PROA C  
ATOM     37  CD1 ILE A   6       3.977   7.887  26.674  1.00 20.59      PROA C  
ATOM     38  N   GLN A   7       3.760  12.742  30.265  1.00 20.52      PROA N  
ATOM     39  CA  GLN A   7       4.262  13.374  31.468  1.00 21.69      PROA C  
ATOM     40  C   GLN A   7       4.683  12.459  32.597  1.00 22.02      PROA C  
ATOM     41  O   GLN A   7       3.954  11.535  32.978  1.00 21.89      PROA O  
ATOM     42  CB  GLN A   7       3.259  14.392  31.997  1.00 22.77      PROA C  
ATOM     43  CG  GLN A   7       3.369  15.749  31.349  1.00 24.86      PROA C  
ATOM     44  CD  GLN A   7       2.467  16.774  32.004  1.00 25.86      PROA C  
ATOM     45  OE1 GLN A   7       1.660  17.417  31.337  1.00 27.24      PROA O  
ATOM     46  NE2 GLN A   7       2.601  16.934  33.325  1.00 26.82      PROA N  
ATOM     47  N   PRO A   8       5.898  12.675  33.125  1.00 22.41      PROA N  
ATOM     48  CA  PRO A   8       6.345  11.830  34.231  1.00 22.84      PROA C  
ATOM     49  C   PRO A   8       5.524  12.215  35.459  1.00 23.50      PROA C  
ATOM     50  O   PRO A   8       5.069  13.359  35.575  1.00 23.62      PROA O  
ATOM     51  CB  PRO A   8       7.821  12.219  34.389  1.00 22.58      PROA C  
ATOM     52  CG  PRO A   8       7.864  13.641  33.894  1.00 22.61      PROA C  
ATOM     53  CD  PRO A   8       6.972  13.583  32.678  1.00 22.34      PROA C  
ATOM     54  N   VAL A   9       5.267  11.244  36.323  1.00 24.00      PROA N  
ATOM     55  CA  VAL A   9       4.516  11.478  37.543  1.00 24.57      PROA C  
ATOM     56  C   VAL A   9       5.471  11.122  38.665  1.00 24.99      PROA C  
ATOM     57  O   VAL A   9       5.927   9.982  38.759  1.00 24.82      PROA O  
ATOM     58  CB  VAL A   9       3.273  10.580  37.613  1.00 24.64      PROA C  
ATOM     59  CG1 VAL A   9       2.596  10.725  38.969  1.00 24.88      PROA C  
ATOM     60  CG2 VAL A   9       2.308  10.935  36.488  1.00 24.78      PROA C  
ATOM     61  N   LEU A  10       5.827  12.119  39.464  1.00 25.78      PROA N  
ATOM     62  CA  LEU A  10       6.752  11.921  40.568  1.00 26.56      PROA C  
ATOM     63  C   LEU A  10       6.043  11.958  41.914  1.00 26.89      PROA C  
ATOM     64  O   LEU A  10       5.187  12.847  42.105  1.00 27.07      PROA O  
ATOM     65  CB  LEU A  10       7.857  12.973  40.501  1.00 26.93      PROA C  
ATOM     66  CG  LEU A  10       8.721  12.752  39.255  1.00 27.39      PROA C  
ATOM     67  CD1 LEU A  10       9.351  14.055  38.769  1.00 27.66      PROA C  
ATOM     68  CD2 LEU A  10       9.768  11.681  39.555  1.00 27.46      PROA C  
ATOM     69  OXT LEU A  10       6.329  11.066  42.743  1.00 27.55      PROA O  
TER      70      LEU A  10                                                      
END                                                                             
eor

# --------------------------------------------------------------------------------------
# then run program to create ADF inputfile
# The program works interactively. The input described here are answers to the questions
# that were asked interactively.
# In cases where the user agrees with the suggestion given by the program,
# the user can press the **Enter** key, which is shown here with an empty line.
# --------------------------------------------------------------------------------------

$AMSBIN/pdb2adf << eor
chymotrypsin.pdb




3 4 5
c
5
3 4 15 16
c


Y
1
1
17.5
eor

# --------------------------------------------------------------------------------------
# Questions asked were:
#  Q1: Please give name of PDB-file
#  A1: chymotrypsin.pdb
#  Q2: Do you want a logfile to be written (Y/n) ?
#  A2: Enter
#  Q3: .. 
#      Found the following terminal amino acid residues : (C-term)    10 (N-term)     1
#      Do you want to use these as terminal residues (Y/n) ?
#  A3: Enter
#  Q4: Multiple AMBER options for CYS :
#        0         Decide every time differently
#        1   CYS   Cysteine (SH)                 
#        2   CYM   Deprotonated Cysteine (S-)    
#        3   CYX   Cystine (S-S bridge)          
#      Suggested option: 0
#  A4: Enter
#  Q5: Multiple AMBER options for CYS    1 (    1) :
#      ...
#      Suggestion: 1
#  A5: Enter
#  Q6: ...
#      Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule
#           1: CYS    1         4: PRO    4         7: GLN    7        10: LEU   10
#           2: GLY    2         5: ALA    5         8: PRO    8
#           3: VAL    3         6: ILE    6         9: VAL    9
#      Give option number of molecules to be put in QM region (or 'c' to continue):
#      Note: by specifying a negative number a molecule is removed from the QM region
#  A6: 3 4 5
#  Q7:    ...
#      Give option number of molecules to be put in QM region (or 'c' to continue):
#      Note: by specifying a negative number a molecule is removed from the QM region
#  A7: c
#  Q8: Make a choice for the QM/MM treatment of VAL    3
#       0:  Put completely in QM region
#       1:  Cut off at C-alpha (put NH in QM region, CO in MM region)
#       2:  Cut off at C-alpha (put NH in MM region, CO in QM region)
#       3:  Cut off at C-alpha (put NH and CO in MM region)
#       4:  Cut off at C-alpha (put NH and CO in QM region,  sidechain in MM region)
#       5:  Put only part of sidechain in QM region
#      Suggestion: 2
#      Give choice:
#  A8: 5
#  Q9: Atoms belonging to molecule VAL     3
#         1: N     MM                6: HB    MM               11: CG2   MM               16: O     MM
#         2: H     MM                7: CG1   MM               12: HG21  MM
#         3: CA    MM                8: HG11  MM               13: HG22  MM
#         4: HA    MM                9: HG12  MM               14: HG23  MM
#         5: CB    MM               10: HG13  MM               15: C     MM
#      Give option number of atoms to be put in QM region (or 'c' to continue):
#      (Note: a range can be entered as 3-21, while a negative number removes an atom)
#  A9: 3 4 15 16
# Q10: ...
#      Give option number of atoms to be put in QM region (or 'c' to continue):
#      (Note: a range can be entered as 3-21, while a negative number removes an atom)
# A10: c
# Q11: Make a choice for the QM/MM treatment of PRO    4
#      ...
#      Suggestion: 2
# A11: Enter
# Q12: Make a choice for the QM/MM treatment of ALA    5
#      ...
#      Suggestion: 1
# A12: Enter
# Q13: Do you want to add solvent to your system (Y/n) ?
# A13: Y
# Q14: Solvent (box) available:
#         1:  HOH     HOH  Water molecule                     
#         2:  MOH     MOH  Methanol molecule                  
#         3:  CHL     CHL  Chloroform molecule
# A14: 1
# Q15: Box Shape options:
#       1  Spherical box
#       2  Cubic box
# A15: 1
# Q16: Give boxsize (def.:    16.71 Angs)
# A16: 17.5
# --------------------------------------------------------------------------------------

# ----------------------------
# now for checking correctness
# ----------------------------

echo ""
echo "ADF inputfile made by pdb2adf"
echo ""

cat chymotrypsin.pdb2adf

rm chymotrypsin.pdb chymotrypsin.pdb.log chymotrypsin.p2a.pdb chymotrypsin.pdb2adf