Example: pdb2adf transforms a PDB file to a QMMM input file¶
#! /bin/sh
# This example shows how to use the utiliy 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