Usage of pdb2adf

Short description

The program works interactively, and should be straightforwardly to use. However, for some of the stages in the output a short description is given below.

                                  P D B 2 A D F - program
                                     version 2008.01
                              Written by: Marcel Swart, 2008

                         This program uses AMBER parameter files
                            see: http://amber.scripps.edu

Do you want a logfile to be written (Y/n) ?

This option exists to create a logfile of what pdb2adf does. However, it should normally be used only for debugging purposes.

Ignoring atom on line:
ATOM    974   OH LYS A 128     -10.073 42.775 15.690 1.00 38.79     5AZU1065

This is a warning that the atom on that particular line is ignored, should normally occur only few times (less than ten). Depends also on how well the PDB file follows the PDB format rules.

Data Processed:
     Nat:       2519
    Nmol:        196
 NChains:          1

Information about what has been read on the PDB file: the total number of atoms (Nat), number of molecules/residues (Nmol) and number of protein chains (Nchains).

Please wait, making connection tables

At this point, the connections between the atoms are being made by looking at atom distances. It may take a while, depending on the size of the system.

Do you want to make separate files for each chain (Y/n) ?

You have the option to make different inputfiles for different protein chains, but you can also make one inputfile for all of them together.

Found the following terminal amino acid residues : (C-term) 128 (N-term) 1
Do you want to use these as terminal residues (Y/n) ?

Info is given about the C- and N-terminal residue of each chain. Reported for making sure they are chosen correctly. Note, if the C- and N-terminal residues are connected (rarely the case probably), enter N here.

Multiple AMBER options for HIS :
  0         Decide every time differently
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

Suggested option: 0

For a number of residues (His, Glu, Asp, Lys and Cys) there is more than one option available in the AMBER95 force field, depending on the protonation state (His, Glu, Asp and Lys) or the existence of a sulphur bridge/connection to a metal atom (Cys). The default is to choose a different option for the His and Cys residues, and use one option for Glu, Asp and Lys (fully charged). However, if wanted you can make a choice for all residues.

Multiple AMBER options for CYS    3 (    3) :
  1   CYS   Cysteine (SH)
  2   CYM   Deprotonated Cysteine (S-)
  3   CYX   Cystine (S-S bridge)

  Connections and Nearest Atoms for SG CYS    3 SG  ( P2A #   41 PDB#   20 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.82      38      19  CB  CYS    3  CB       1       3.79    2382     980  O   HOH  151  O
 2    2.02     461     193  SG  CYS   26  SG       2       3.80      22       0  HC  GLN    2
                                                   3       4.04    2391     983  O   HOH  154  O
                                                   4       4.15     509     206  O   GLN   28  O
                                                   5       4.18     522       0  HA  PHE   29
Suggestion: 3

The options for Cys3 are given, with information about the atoms bonded to the SG sulphur atom (on the left), as well as the closest five non-bonded atoms (on the right). This information may help you decide which choice to make for this particular residue. Also given (on the bottom) is the suggested choice, which is based, in this case, on the presence of a sulphur bridge.

Multiple AMBER options for HIS   46 (   46) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS   46 ND1 ( P2A #  844 PDB#  347 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.37     843     346  CG  HIS   46  CG       1       2.62    2166       0  H1  MET  121
 2    1.33     846     349  CE  HIS   46  CE1      2       3.23    2080     863  ND  HIS  117  ND1
 3    2.04    2318     959  CU  CU   130  CU       3 HB    3.33    2163     900  S   MET  121  SD
                                                   4       3.40    2164     901  CT  MET  121  CE
                                                   5       3.57    2082     865  CE  HIS  117  CE1

  Connections and Nearest Atoms for NE HIS   46 NE2 ( P2A #  848 PDB#  350 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.32     846     349  CE  HIS   46  CE1      1 HB    2.70     162      67  O   ASN   10  O
 2    1.37     850     348  CD  HIS   46  CD2      2       2.83     814       0  H1  MET   44
                                                   3       3.23    2166       0  H1  MET  121
                                                   4       3.52     822     332  O   MET   44  O
                                                   5       3.74     813     334  CT  MET   44  CG
Suggestion: 2

For His residues, the information is given for both the delta- and the epsilon nitrogen atoms. Also indicated (by HB) is the presence of a hydrogen bond with another atom. The definition used here is that two atoms are hydrogen bonded if they are both non-carbon/non-hydrogen atoms, and the distance between them is less than the sum of the van der Waals radii of the atoms. It is a simple definition, but seems to be effective. In this case, as the N(delta) is bonded to copper, the proton should be attached to the N(epsilon).

Making choice for which molecules should be QM, which MM

Now we come to the part where the division in the QM and MM systems is made.

Residues belonging to chain  0
Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule
     1: ALA    1        28: GLN   28        55: ASP   55        82: ALA   82       109: ALA  109
     2: GLN    2        29: PHE   29        56: LYS   56        83: HIS   83       110: TYR  110
etc

All molecules/residues belonging to chain 0 are given, with an option number.

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

Here you are asked to enter the option numbers of the residues you want to put in the QM system.

Putting GLY   45 in QM region
Putting HIS   46 in QM region

In this case, Gly45 and His46 have been put in the QM system.

Make a choice for the QM/MM treatment of GLY   45
 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:

A choice should be made for where to cut-off the QM system. Normally this is done at the C(alpha) position, and you should simply choose the Suggestion.

Solvent molecules (SOL/HOH) belonging to this chain:
    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15   16   17   18   19   20
   21   22   23   24   25   26   27   28   29   30   31   32   33   34   35   36   37   38   39   40
   41   42   43   44   45   46   47   48   49   50   51   52   53   54   55   56   57   58   59   60
   61   62   63   64   65   66

Give the number of the molecule to be put in QM region (or 'c' to continue):

Also water molecules can be put in the QM system.

Box Shape options:
 1  Spherical box
 2  Cubic box
Make a choice:

Type of box to be used.

Maximum atomic distance (Angs) from center       25.62
Give boxsize (def.:    28.62 Angs)

Size of box to be used to put a layer of solvent molecules around the system. Max. dist. is the maximal distance of any protein atom from the center of mass of the protein. Usually you should choose a boxsize at least 6 Angstrom larger (so at least two solvent molecules are surrounding the system).

...

Using BOXSIZE value of  30.0000
Adding atoms for box     1  Added (Box):     0 (Total):     0  Excl. (1):   648  Excl. (2):     0
Adding atoms for box     2  Added (Box):     9 (Total):     9  Excl. (1):   639  Excl. (2):     0
Adding atoms for box    63  Added (Box):     3 (Total):  7635  Excl. (1):   645  Excl. (2):     0
Adding atoms for box    64  Added (Box):     0 (Total):  7635  Excl. (1):   648  Excl. (2):     0
Writing inputfile for chain   1

A total amount of 7635 atoms (2545 water molecules) has been added.

Inputfile(s) written, everything processed, work has been done.
Thank you for using the PDB2ADF program.

================================
Normal ending of PDB2ADF program
================================

ADF inputfile(s) have been written, the PDB-file has been processes. Everything is done.

An example on protein structure

The idea of this example is to make an adf-input file using a PDB of an azurin (1DYZ.pdb.txt). The result of this example should be that in the adf-input file the active site of azurin (Figure 1) is in the QM part, and the rest of the protein is in the MM part, and that the solvent water is added (in a box), which is also in the MM part.

/scm-uploads/doc.2019/QMMM/_images/azurin.png

Figure 1: the active site of azurin

Usage of pdb2adf

The program works interactively. Given below in bold are the parts that the user has to type. In cases where the user agrees with the suggestion given by the program, the user can press the Enter key indicated with Enter.

                                  P D B 2 A D F - program
                                     version 2008.01
                              Written by: Marcel Swart, 2008

                         This program uses AMBER parameter files
                            see: http://amber.scripps.edu

Please give name of PDB-file

1DYZ.pdb.txt

Do you want a logfile to be written (Y/n) ?

Enter

 read fragments

Data Processed:
     Nat:       2519
    Nmol:        196
 NChains:          1

Please wait, making connection tables
Now finding nearby atoms
Assigning chain ID to all residues
Completing residues for which only option is available

Found the following terminal amino acid residues : (C-term)   129 (N-term)     1
Do you want to use these as terminal residues (Y/n) ?

Enter

Refinding nearby atoms (including atoms added in residue completion)

Multiple AMBER options for HIS :
  0         Decide every time differently
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

Suggested option: 0

Enter

Using 0: Decide every time differently

Multiple AMBER options for GLU :
  0         Decide every time differently
  1   GLU   Glutamic acid (COO-)
  2   GLH   Neutral Glutamic acid (COOH)

Suggested option: 1

Enter

Using 17   GLU   Glutamic acid (COO-)

Multiple AMBER options for ASP :
  0         Decide every time differently
  1   ASP   Aspartic acid (COO-)
  2   ASH   Neutral Aspartatic acid (COOH)

Suggested option: 1

Enter

Using 18   ASP   Aspartic acid (COO-)

Multiple AMBER options for LYS :
  0         Decide every time differently
  1   LYS   Charged Lysine (NH3+)
  2   LYN   Neutral Lysine (NH2)

Suggested option: 1

Enter

Using 19   LYS   Charged Lysine (NH3+)

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

Enter

Using 0: Decide every time differently
 - - - - - - - - - - - - - - - - - - - - - - - - - - -
           Making Choices for Chain      0
 - - - - - - - - - - - - - - - - - - - - - - - - - - -

Multiple AMBER options for CYS    3 (    3) :
  1   CYS   Cysteine (SH)
  2   CYM   Deprotonated Cysteine (S-)
  3   CYX   Cystine (S-S bridge)

  Connections and Nearest Atoms for SG CYS    3 SG  ( P2A #   41 PDB#   20 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.82      38      19  CB  CYS    3  CB       1       3.79    2382     980  O   HOH  151  O
 2    2.02     461     193  SG  CYS   26  SG       2       3.80      22       0  HC  GLN    2
                                                   3       4.04    2391     983  O   HOH  154  O
                                                   4       4.15     509     206  O   GLN   28  O
                                                   5       4.18     522       0  HA  PHE   29
Suggestion: 3

Enter

Multiple AMBER options for CYS   26 (   26) :
  1   CYS   Cysteine (SH)
  2   CYM   Deprotonated Cysteine (S-)
  3   CYX   Cystine (S-S bridge)

  Connections and Nearest Atoms for SG CYS   26 SG  ( P2A #  461 PDB#  193 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.82     458     192  CB  CYS   26  CB       1       3.41     522       0  HA  PHE   29
 2    2.02      41      20  SG  CYS    3  SG       2       3.43     411     168  O   ASP   23  O
                                                   3       3.60    2322     960  O   HOH  131  O
                                                   4       3.91     403     169  CB  ASP   23  CB
                                                   5       4.15     387       0  HC  VAL   22
Suggestion: 3

Enter

Multiple AMBER options for HIS   32 (   32) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS   32 ND1 ( P2A #  581 PDB#  244 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.39     580     243  CG  HIS   32  CG       1       3.41     545       0  HC  THR   30
 2    1.33     583     246  CE  HIS   32  CE1      2       3.43      76      33  O   ALA    5  O
                                                   3       3.58      90      40  OH  THR    6  OG1
                                                   4       3.99      91       0  HO  THR    6
                                                   5       4.17      68       0  H   ALA    5

  Connections and Nearest Atoms for NE HIS   32 NE2 ( P2A #  585 PDB#  247 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.31     583     246  CE  HIS   32  CE1      1       2.86     544       0  HC  THR   30
 2    1.37     587     245  CD  HIS   32  CD2      2       3.00     545       0  HC  THR   30
                                                   3       3.14    1677       0  HO  SER   94
                                                   4       3.42     542     229  CT  THR   30  CG2
                                                   5       3.65    1676     688  OH  SER   94  OG
Suggestion: 1

3

Multiple AMBER options for HIS   35 (   35) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS   35 ND1 ( P2A #  649 PDB#  271 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.38     648     270  CG  HIS   35  CG       1       2.46     682       0  H   GLY   37
 2    1.32     651     273  CE  HIS   35  CE1      2       2.69    1604       0  H1  GLY   89
                                                   3       3.31     681     282  N   GLY   37  N
                                                   4       3.56    1602     653  CT  GLY   89  CA
                                                   5       3.67     152       0  H1  ASN   10

  Connections and Nearest Atoms for NE HIS   35 NE2 ( P2A #  653 PDB#  274 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.33     651     273  CE  HIS   35  CE1      1 HB    2.91     822     332  O   MET   44  O
 2    1.37     655     272  CD  HIS   35  CD2      2       3.24     814       0  H1  MET   44
                                                   3       3.24     850     348  CD  HIS   46  CD2
                                                   4       3.34    1593       0  H1  GLY   88
                                                   5       3.75     848     350  NE  HIS   46  NE2
Suggestion: 2

3

Multiple AMBER options for HIS   46 (   46) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS   46 ND1 ( P2A #  844 PDB#  347 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.37     843     346  CG  HIS   46  CG       1       2.62    2166       0  H1  MET  121
 2    1.33     846     349  CE  HIS   46  CE1      2       3.23    2080     863  ND  HIS  117  ND1
 3    2.04    2318     959  CU  CU   130  CU       3 HB    3.33    2163     900  S   MET  121  SD
                                                   4       3.40    2164     901  CT  MET  121  CE
                                                   5       3.57    2082     865  CE  HIS  117  CE1

  Connections and Nearest Atoms for NE HIS   46 NE2 ( P2A #  848 PDB#  350 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.32     846     349  CE  HIS   46  CE1      1 HB    2.70     162      67  O   ASN   10  O
 2    1.37     850     348  CD  HIS   46  CD2      2       2.83     814       0  H1  MET   44
                                                   3       3.23    2166       0  H1  MET  121
                                                   4       3.52     822     332  O   MET   44  O
                                                   5       3.74     813     334  CT  MET   44  CG
Suggestion: 2

Enter

Multiple AMBER options for HIS   83 (   83) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS   83 ND1 ( P2A # 1494 PDB#  613 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.39    1493     612  CG  HIS   83  CG       1       2.67    1317       0  HC  VAL   73
 2    1.33    1496     615  CE  HIS   83  CE1      2       3.63    1315     542  CT  VAL   73  CG2
                                                   3       3.74    1310       0  HC  VAL   73
                                                   4       3.82    1316       0  HC  VAL   73
                                                   5       3.86    1313       0  HC  VAL   73

  Connections and Nearest Atoms for NE HIS   83 NE2 ( P2A # 1498 PDB#  616 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.32    1496     615  CE  HIS   83  CE1      1       3.09    1313       0  HC  VAL   73
 2    1.38    1500     614  CD  HIS   83  CD2      2       3.44    1317       0  HC  VAL   73
                                                   3       3.88    2385     981  O   HOH  152  O
                                                   4       3.93    1311     541  CT  VAL   73  CG1
                                                   5       4.03    1309     540  CT  VAL   73  CB
Suggestion: 2

3

Multiple AMBER options for CYS  112 (  112) :
  1   CYS   Cysteine (SH)
  2   CYM   Deprotonated Cysteine (S-)
  3   CYX   Cystine (S-S bridge)

  Connections and Nearest Atoms for SG CYS  112 SG  ( P2A # 2001 PDB#  828 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.82    1998     827  CB  CYS  112  CB       1       2.53     858       0  H   ASN   47
 2    2.14    2318     959  CU  CU   130  CU       2       2.65    2023       0  H   PHE  114
                                                   3       3.00    2028       0  HC  PHE  114
                                                   4       3.29     868       0  H   ASN   47
                                                   5       3.39    2027       0  HC  PHE  114
Suggestion: 2

Enter

Multiple AMBER options for HIS  117 (  117) :
  1   HID   Histidine Delta Hydrogen
  2   HIE   Histidine Epsilon Hydrogen
  3   HIP   Histidine E & D Hydrogens

  Connections and Nearest Atoms for ND HIS  117 ND1 ( P2A # 2080 PDB#  863 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.37    2079     862  CG  HIS  117  CG       1       2.82    2028       0  HC  PHE  114
 2    1.34    2082     865  CE  HIS  117  CE1      2       3.23     844     347  ND  HIS   46  ND1
 3    1.99    2318     959  CU  CU   130  CU       3       3.26    2031       0  HA  PHE  114
                                                   4       3.27     832     340  O   GLY   45  O
                                                   5       3.43     846     349  CE  HIS   46  CE1

  Connections and Nearest Atoms for NE HIS  117 NE2 ( P2A # 2084 PDB#  866 )
      Dist  P2A Nr  PDB Nr  Label                 Near     Dist  P2A Nr  PDB Nr  Label
 1    1.31    2082     865  CE  HIS  117  CE1      1       2.57     209       0  H1  MET   13
 2    1.37    2086     864  CD  HIS  117  CD2      2       2.65    2031       0  HA  PHE  114
                                                   3 HB    2.74    2406     988  O   HOH  159  O
                                                   4       3.34    2030     841  CA  PHE  114  CD1
                                                   5       3.41     204       0  H1  MET   13
Suggestion: 2

Enter

 - - - - - - - - - - - - - - - - - - - - - - - - - - -
           Making Choices for Chain      1
 - - - - - - - - - - - - - - - - - - - - - - - - - - -

Completing residues with multiple options available, and solvent molecules

Checking positions of newly added atoms

Making choice for which molecules should be QM, which MM

Residues belonging to chain  0
Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule    Option  Molecule
     1: ALA    1        28: GLN   28        55: ASP   55        82: ALA   82       109: ALA  109
     2: GLN    2        29: PHE   29        56: LYS   56        83: HIS   83       110: TYR  110
     3: CYS    3        30: THR   30        57: GLN   57        84: THR   84       111: PHE  111
     4: GLU    4        31: MET   31        58: ALA   58        85: LYS   85       112: CYS  112
     5: ALA    5        32: HIS   32        59: VAL   59        86: VAL   86       113: SER  113
     6: THR    6        33: LEU   33        60: ALA   60        87: ILE   87       114: PHE  114
     7: VAL    7        34: LYS   34        61: THR   61        88: GLY   88       115: PRO  115
     8: GLU    8        35: HIS   35        62: ASP   62        89: GLY   89       116: GLY  116
     9: SER    9        36: VAL   36        63: GLY   63        90: GLY   90       117: HIS  117
    10: ASN   10        37: GLY   37        64: MET   64        91: GLU   91       118: TRP  118
    11: ASP   11        38: LYS   38        65: GLY   65        92: SER   92       119: ALA  119
    12: ALA   12        39: MET   39        66: ALA   66        93: ASP   93       120: MET  120
    13: MET   13        40: ALA   40        67: GLY   67        94: SER   94       121: MET  121
    14: GLN   14        41: LYS   41        68: LEU   68        95: VAL   95       122: LYS  122
    15: TYR   15        42: VAL   42        69: ALA   69        96: THR   96       123: GLY  123
    16: ASN   16        43: ALA   43        70: GLN   70        97: PHE   97       124: THR  124
    17: VAL   17        44: MET   44        71: ASP   71        98: ASP   98       125: LEU  125
    18: LYS   18        45: GLY   45        72: TYR   72        99: VAL   99       126: LYS  126
    19: GLU   19        46: HIS   46        73: VAL   73       100: SER  100       127: LEU  127
    20: ILE   20        47: ASN   47        74: LYS   74       101: LYS  101       128: GLY  128
    21: VAL   21        48: LEU   48        75: ALA   75       102: ILE  102       129: SER  129
    22: VAL   22        49: VAL   49        76: GLY   76       103: ALA  103       130: CU   130
    23: ASP   23        50: LEU   50        77: ASP   77       104: ALA  104
    24: LYS   24        51: THR   51        78: THR   78       105: GLY  105
    25: SER   25        52: LYS   52        79: ARG   79       106: GLU  106
    26: CYS   26        53: ASP   53        80: VAL   80       107: ASN  107
    27: LYS   27        54: ALA   54        81: ILE   81       108: TYR  108

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

45 46 112 117 121 130

Putting GLY   45 in QM region
Putting HIS   46 in QM region
Putting CYS  112 in QM region
Putting HIS  117 in QM region
Putting MET  121 in QM region
Putting CU   130 in QM region

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

c

Make a choice for the QM/MM treatment of GLY   45
 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:

Enter

Make a choice for the QM/MM treatment of HIS   46
 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: 1
Give choice:

Enter

Make a choice for the QM/MM treatment of CYS  112
 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: 3
Give choice:

Enter

Make a choice for the QM/MM treatment of HIS  117
 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: 3
Give choice:

Enter

Make a choice for the QM/MM treatment of MET  121
 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: 3
Give choice:

Enter

Make a choice for the QM/MM treatment of CU   130
 0:  Put completely in QM region
 1:  Put only part of molecule in QM region

Suggestion: 0
Give choice:

Enter

Total formal charge on molecule CU    130      2.0000

Solvent molecules (SOL/HOH) belonging to this chain:
    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15   16   17   18   19   20
   21   22   23   24   25   26   27   28   29   30   31   32   33   34   35   36   37   38   39   40
   41   42   43   44   45   46   47   48   49   50   51   52   53   54   55   56   57   58   59   60
   61   62   63   64   65   66

Give the number of the molecule to be put in QM region (or 'c' to continue):

c

Residues belonging to chain  1

Do you want to add solvent to your system (Y/n) ?

Enter

Solvent (box) available:
   1:  HOH     HOH  Water molecule
   2:  MOH     MOH  Methanol molecule
   3:  CHL     CHL  Chloroform molecule

1

Reading contents of solvent box p2abox.HOH

Box Shape options:
 1  Spherical box
 2  Cubic box
Make a choice:

1

Writing inputfile for chain   0

Using total charge   1.0 and total spin   1.0

Maximum atomic distance (Angs) from center       25.62
Give boxsize (def.:    28.62 Angs)

30.0

Using BOXSIZE value of  30.0000
Adding atoms for box     1  Added (Box):     0 (Total):     0  Excl. (1):   648  Excl. (2):     0
Adding atoms for box     2  Added (Box):     9 (Total):     9  Excl. (1):   639  Excl. (2):     0
Adding atoms for box     3  Added (Box):     3 (Total):    12  Excl. (1):   645  Excl. (2):     0
Adding atoms for box     4  Added (Box):     0 (Total):    12  Excl. (1):   648  Excl. (2):     0
Adding atoms for box     5  Added (Box):     6 (Total):    18  Excl. (1):   642  Excl. (2):     0
Adding atoms for box     6  Added (Box):   228 (Total):   246  Excl. (1):   420  Excl. (2):     0
Adding atoms for box     7  Added (Box):   219 (Total):   465  Excl. (1):   429  Excl. (2):     0
Adding atoms for box     8  Added (Box):     9 (Total):   474  Excl. (1):   639  Excl. (2):     0
Adding atoms for box     9  Added (Box):     0 (Total):   474  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    10  Added (Box):   225 (Total):   699  Excl. (1):   423  Excl. (2):     0
Adding atoms for box    11  Added (Box):   216 (Total):   915  Excl. (1):   432  Excl. (2):     0
Adding atoms for box    12  Added (Box):     6 (Total):   921  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    13  Added (Box):     0 (Total):   921  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    14  Added (Box):     6 (Total):   927  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    15  Added (Box):    12 (Total):   939  Excl. (1):   636  Excl. (2):     0
Adding atoms for box    16  Added (Box):     0 (Total):   939  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    17  Added (Box):    12 (Total):   951  Excl. (1):   636  Excl. (2):     0
Adding atoms for box    18  Added (Box):   210 (Total):  1161  Excl. (1):   438  Excl. (2):     0
Adding atoms for box    19  Added (Box):   219 (Total):  1380  Excl. (1):   429  Excl. (2):     0
Adding atoms for box    20  Added (Box):     3 (Total):  1383  Excl. (1):   645  Excl. (2):     0
Adding atoms for box    21  Added (Box):   216 (Total):  1599  Excl. (1):   417  Excl. (2):    15
Adding atoms for box    22  Added (Box):   381 (Total):  1980  Excl. (1):     3  Excl. (2):   264
Adding atoms for box    23  Added (Box):   261 (Total):  2241  Excl. (1):     3  Excl. (2):   384
Adding atoms for box    24  Added (Box):   183 (Total):  2424  Excl. (1):   423  Excl. (2):    42
Adding atoms for box    25  Added (Box):   189 (Total):  2613  Excl. (1):   426  Excl. (2):    33
Adding atoms for box    26  Added (Box):   186 (Total):  2799  Excl. (1):     3  Excl. (2):   459
Adding atoms for box    27  Added (Box):   351 (Total):  3150  Excl. (1):     3  Excl. (2):   294
Adding atoms for box    28  Added (Box):   222 (Total):  3372  Excl. (1):   420  Excl. (2):     6
Adding atoms for box    29  Added (Box):     9 (Total):  3381  Excl. (1):   639  Excl. (2):     0
Adding atoms for box    30  Added (Box):   162 (Total):  3543  Excl. (1):   429  Excl. (2):    57
Adding atoms for box    31  Added (Box):   219 (Total):  3762  Excl. (1):   426  Excl. (2):     3
Adding atoms for box    32  Added (Box):     6 (Total):  3768  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    33  Added (Box):     6 (Total):  3774  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    34  Added (Box):   219 (Total):  3993  Excl. (1):   426  Excl. (2):     3
Adding atoms for box    35  Added (Box):   216 (Total):  4209  Excl. (1):   432  Excl. (2):     0
Adding atoms for box    36  Added (Box):     6 (Total):  4215  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    37  Added (Box):   219 (Total):  4434  Excl. (1):   426  Excl. (2):     3
Adding atoms for box    38  Added (Box):   279 (Total):  4713  Excl. (1):     6  Excl. (2):   363
Adding atoms for box    39  Added (Box):   231 (Total):  4944  Excl. (1):     0  Excl. (2):   417
Adding atoms for box    40  Added (Box):   195 (Total):  5139  Excl. (1):   432  Excl. (2):    21
Adding atoms for box    41  Added (Box):   231 (Total):  5370  Excl. (1):   414  Excl. (2):     3
Adding atoms for box    42  Added (Box):   324 (Total):  5694  Excl. (1):     0  Excl. (2):   324
Adding atoms for box    43  Added (Box):   408 (Total):  6102  Excl. (1):     6  Excl. (2):   234
Adding atoms for box    44  Added (Box):   204 (Total):  6306  Excl. (1):   435  Excl. (2):     9
Adding atoms for box    45  Added (Box):     6 (Total):  6312  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    46  Added (Box):   177 (Total):  6489  Excl. (1):   435  Excl. (2):    36
Adding atoms for box    47  Added (Box):   219 (Total):  6708  Excl. (1):   429  Excl. (2):     0
Adding atoms for box    48  Added (Box):     6 (Total):  6714  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    49  Added (Box):     0 (Total):  6714  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    50  Added (Box):     3 (Total):  6717  Excl. (1):   645  Excl. (2):     0
Adding atoms for box    51  Added (Box):     6 (Total):  6723  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    52  Added (Box):     0 (Total):  6723  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    53  Added (Box):     9 (Total):  6732  Excl. (1):   639  Excl. (2):     0
Adding atoms for box    54  Added (Box):   222 (Total):  6954  Excl. (1):   426  Excl. (2):     0
Adding atoms for box    55  Added (Box):   213 (Total):  7167  Excl. (1):   426  Excl. (2):     9
Adding atoms for box    56  Added (Box):     6 (Total):  7173  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    57  Added (Box):     3 (Total):  7176  Excl. (1):   645  Excl. (2):     0
Adding atoms for box    58  Added (Box):   219 (Total):  7395  Excl. (1):   423  Excl. (2):     6
Adding atoms for box    59  Added (Box):   219 (Total):  7614  Excl. (1):   429  Excl. (2):     0
Adding atoms for box    60  Added (Box):     6 (Total):  7620  Excl. (1):   642  Excl. (2):     0
Adding atoms for box    61  Added (Box):     0 (Total):  7620  Excl. (1):   648  Excl. (2):     0
Adding atoms for box    62  Added (Box):    12 (Total):  7632  Excl. (1):   636  Excl. (2):     0
Adding atoms for box    63  Added (Box):     3 (Total):  7635  Excl. (1):   645  Excl. (2):     0
Adding atoms for box    64  Added (Box):     0 (Total):  7635  Excl. (1):   648  Excl. (2):     0
Writing inputfile for chain   1
There are no atoms in this chain, ignoring it

Inputfile(s) written, everything processed, work has been done.
Thank you for using the PDB2ADF program.

================================
Normal ending of PDB2ADF program
================================

Contents of the 1DYZ.pdb2adf file generated by pdb2adf

The file is not given completely, since it contains more than 9000 atoms. Note that the old style QM/MM input format is used if the environment variable SCM_PDB2ADF is set to OLD. This example uses the default NEWQMMM format.

#! /bin/sh

$ADFBIN/adf << eor
TITLE QM/MM calculation setup by pdb2adf: M.Swart, 2008

Symmetry NOSYM

EPRINT
  SFO NOEIG NOOVL
END

XC
 GGA Becke-Perdew
END

GEOMETRY
 iterations 100
 CONVERGE grad=1.0e-3 rad=0.1
END

BASIS
 type TZP
 core small
END

SCF
 diis ok=0.01
 Converge 1.0e-5 1.0e-5
 Iterations 99
END

CHARGE   1.0  1.0
UNRESTRICTED

ATOMS
     1 N     -1.1930    25.6890    17.1840
     2 H     -0.3133    25.1929    17.1970
     3 H     -1.3738    25.1438    18.0148
     4 H     -1.5170    24.8559    16.7138
     5 C     -1.4820    27.1340    16.8960
     6 H     -2.1350    27.2082    16.0264
     7 C     -2.1950    27.7860    18.0880
     8 H     -1.5602    27.7210    18.9717
     9 H     -2.3971    28.8331    17.8627
    10 H     -3.1350    27.2677    18.2776
    11 C     -0.1820    27.8790    16.5880
    12 O      0.8890    27.4920    17.0690
    13 N     -0.2890    28.9420    15.7940
    14 H     -1.1936    29.2105    15.4339
    15 C      0.8750    29.7430    15.4220
    16 H      0.5616    30.5606    14.7728
    17 C      1.5270    30.3290    16.6860
    18 H      0.7563    30.7529    17.3298
    19 H      2.0544    29.5402    17.2224
    20 C      2.5320    31.4280    16.4240
    21 H      2.9095    31.3405    15.4052
    22 H      2.0516    32.3982    16.5508
    23 C      3.7140    31.3600    17.3720
    24 O      4.5190    30.4170    17.3220
    25 N      3.8270    32.3530    18.2440
    26 H      4.5950    32.3600    18.9000
    27 H      3.1446    33.0976    18.2480
    28 C      1.8460    28.8170    14.6770
    29 O      1.4330    28.0490    13.7940
    30 N      3.1200    28.8760    15.0450
    31 H      3.4036    29.5187    15.7707
    32 C      4.1380    28.0400    14.4360
    33 H      3.6966    27.4628    13.6235
    34 C      5.2490    28.9240    13.8830
    35 H      5.5773    29.6209    14.6541
    36 H      6.0893    28.3021    13.5744
    37 S      4.6970    29.8870    12.4450
    38 C      4.6800    27.1010    15.5090
    39 O      5.8890    27.0060    15.7280
    40 N      3.7690    26.4380    16.2040
    41 H      2.7886    26.5676    15.9987
    42 C      4.1640    25.5320    17.2520
    43 H      5.0753    25.0110    16.9583
    44 C      4.4480    26.3010    18.5380
    45 H      4.7889    25.6090    19.3080
    46 H      5.2208    27.0470    18.3527
    47 C      3.2450    27.0350    19.0990
    48 H      2.6884    27.4975    18.2839
    49 H      2.6007    26.3297    19.6240
    50 C      3.6030    28.1430    20.0880
    51 O      4.7940    28.3110    20.4290
    52 O      2.6740    28.8670    20.5030
    53 C      3.0820    24.5140    17.4740
    54 O      1.9190    24.7350    17.1470
    55 N      3.4810    23.3770    18.0170
    56 H      4.4563    23.2437    18.2430
    57 C      2.5550    22.3010    18.3030
    58 H      1.5462    22.7037    18.3941
    59 C      2.5590    21.2690    17.1740
    60 H      3.5602    20.8520    17.0658
    61 H      1.8565    20.4695    17.4093
    62 H      2.2632    21.7490    16.2412
    63 C      2.9700    21.6650    19.6140
    64 O      4.1160    21.8010    20.0550
    65 N      2.0180    21.0240    20.2660
    66 H      1.0839    20.9887    19.8835
    67 C      2.2640    20.3560    21.5330
    68 H      3.3189    20.4450    21.7924
    69 C      1.4300    21.0100    22.6680
    70 H      0.3783    21.0286    22.3823
    71 C      1.5670    20.2160    23.9790
    72 H      2.6136    20.1931    24.2828
    73 H      0.9735    20.6946    24.7580
    74 H      1.2111    19.1971    23.8262
    75 O      1.8900    22.3550    22.8810
    76 H      1.3759    22.7617    23.5824
    77 C      1.9150    18.8760    21.4070
    78 O      0.8730    18.5150    20.8640
...
  9746 O     31.1328    34.4612    22.6903
  9747 H     31.8908    34.5740    22.1167
  9748 H     30.6706    35.2981    22.6446
END

QMMM
  FORCE_FIELD_FILE $ADFRESOURCES/ForceFields/amber95.ff
  ELSTAT_COUPLING_MODEL=1

  LINK_BONDS
      691 -   689   1.3320 H H1
      698 -   711   1.3861 H H1
     1648 -  1646   1.3362 H H1
     1648 -  1654   1.3927 H H1
     1710 -  1708   1.3471 H H1
     1710 -  1723   1.3951 H H1
     1778 -  1776   1.3310 H H1
     1778 -  1791   1.3799 H H1
  SUBEND
  ! ---------------------------------------------------
  ! use NEWQMMM developer implementation
  !
  ! NOTE:
  ! the rest of QMMM block is only for use with NEWQMMM
  ! ---------------------------------------------------
  NEWQMMM
  QMMM_PREOPTIMIZE

  QM_INFO
      689    H1  LI   -0.055275  GLY    45  N       691    687    690
      691  CT    QM   -0.025200  GLY    45  CA      689    694    692    693
      692  H1    QM    0.069800  GLY    45  HA2     691
      693  H1    QM    0.069800  GLY    45  HA3     691
      694  C     QM    0.597300  GLY    45  C       691    695    696
      695  O     QM   -0.567900  GLY    45  O       694
      696  N     QM   -0.415700  HIS    46  N       698    694    697
      697  H     QM    0.271900  HIS    46  H       696
      698  CT    QM   -0.058100  HIS    46  CA      696    700    711    699
      699  H1    QM    0.136000  HIS    46  HA      698
      700  CT    QM   -0.007400  HIS    46  CB      698    703    701    702
      701  HC    QM    0.036700  HIS    46  HB2     700
      702  HC    QM    0.036700  HIS    46  HB3     700
      703  CC    QM    0.186800  HIS    46  CG      700    704    709
      704  NB    QM   -0.543200  HIS    46  ND1     703    705   1915
      705  CR    QM    0.163500  HIS    46  CE1     704    707    706
      706  H5    QM    0.143500  HIS    46  HE1     705
      707  NA    QM   -0.279500  HIS    46  NE2     705    709    708
      708  H     QM    0.333900  HIS    46  HE2     707
      709  CW    QM   -0.220700  HIS    46  CD2     703    707    710
      710  H4    QM    0.186200  HIS    46  HD2     709
      711    H1  LI   -0.055275  HIS    46  C       698    712    713
     1646    H1  LI   -0.055275  CYS   112  N      1648   1644   1647
     1648  CT    QM    0.035000  CYS   112  CA     1646   1650   1654   1649
     1649  H1    QM    0.048000  CYS   112  HA     1648
     1650  CT    QM   -0.736000  CYS   112  CB     1648   1653   1651   1652
     1651  H1    QM    0.244000  CYS   112  HB3    1650
     1652  H1    QM    0.244000  CYS   112  HB2    1650
     1653  SH    QM   -0.736000  CYS   112  SG     1650   1915
     1654    H1  LI   -0.055275  CYS   112  C      1648   1655   1656
     1708    H1  LI   -0.055275  HIS   117  N      1710   1706   1709
     1710  CT    QM   -0.058100  HIS   117  CA     1708   1712   1723   1711
     1711  H1    QM    0.136000  HIS   117  HA     1710
     1712  CT    QM   -0.007400  HIS   117  CB     1710   1715   1713   1714
     1713  HC    QM    0.036700  HIS   117  HB2    1712
     1714  HC    QM    0.036700  HIS   117  HB3    1712
     1715  CC    QM    0.186800  HIS   117  CG     1712   1716   1721
     1716  NB    QM   -0.543200  HIS   117  ND1    1715   1717   1915
     1717  CR    QM    0.163500  HIS   117  CE1    1716   1719   1718
     1718  H5    QM    0.143500  HIS   117  HE1    1717
     1719  NA    QM   -0.279500  HIS   117  NE2    1717   1721   1720
     1720  H     QM    0.333900  HIS   117  HE2    1719
     1721  CW    QM   -0.220700  HIS   117  CD2    1715   1719   1722
     1722  H4    QM    0.186200  HIS   117  HD2    1721
     1723    H1  LI   -0.055275  HIS   117  C      1710   1724   1725
     1776    H1  LI   -0.055275  MET   121  N      1778   1774   1777
     1778  CT    QM   -0.023700  MET   121  CA     1776   1780   1791   1779
     1779  H1    QM    0.088000  MET   121  HA     1778
     1780  CT    QM    0.034200  MET   121  CB     1778   1783   1781   1782
     1781  HC    QM    0.024100  MET   121  HB2    1780
     1782  HC    QM    0.024100  MET   121  HB3    1780
     1783  CT    QM    0.001800  MET   121  CG     1780   1786   1784   1785
     1784  H1    QM    0.044000  MET   121  HG2    1783
     1785  H1    QM    0.044000  MET   121  HG3    1783
     1786  S     QM   -0.273700  MET   121  SD     1783   1787
     1787  CT    QM   -0.053600  MET   121  CE     1786   1788   1789   1790
     1788  H1    QM    0.068400  MET   121  HE1    1787
     1789  H1    QM    0.068400  MET   121  HE2    1787
     1790  H1    QM    0.068400  MET   121  HE3    1787
     1791    H1  LI   -0.055275  MET   121  C      1778   1792   1793
     1915  CU    QM    2.000000  CU    130  CU      704   1653   1716
  SUBEND

  QMMM_INFO
        1  N3    MM    0.141400  ALA     1  N         5      2      3      4
        2  H     MM    0.199700  ALA     1  H1        1
        3  H     MM    0.199700  ALA     1  H2        1
        4  H     MM    0.199700  ALA     1  H3        1
        5  CT    MM    0.096200  ALA     1  CA        1      7     11      6
        6  HP    MM    0.088900  ALA     1  HA        5
        7  CT    MM   -0.059700  ALA     1  CB        5      8      9     10
        8  HC    MM    0.030000  ALA     1  HB1       7
        9  HC    MM    0.030000  ALA     1  HB2       7
       10  HC    MM    0.030000  ALA     1  HB3       7
       11  C     MM    0.616300  ALA     1  C         5     12     13
       12  O     MM   -0.572200  ALA     1  O        11
       13  N     MM   -0.415700  GLN     2  N        15     11     14
       14  H     MM    0.271900  GLN     2  H        13
       15  CT    MM   -0.003100  GLN     2  CA       13     17     28     16
       16  H1    MM    0.085000  GLN     2  HA       15
       17  CT    MM   -0.003600  GLN     2  CB       15     20     18     19
       ...
     9746  OW    MM   -0.834000  HOH  2545  O      9747   9748
     9747  HW    MM    0.417000  HOH  2545  H1     9746
     9748  HW    MM    0.417000  HOH  2545  H2     9746
  SUBEND
END

eor

An example on solvent shell run

The idea of this example is to make an adf-input file using a PDB file of water (hoh.pdb.txt), in the solvent methanol. The water molecule in the adf-input file should be in the QM part, and the solvent methanol (in a box) is in MM part.

Contents of the hoh.pdb file

TITLE     PDB-FILE CORRESPONDING TO pdb2adf-GENERATED ADF-INPUTFILE
REMARK    Written by M. Swart, March 2005
HETATM    1 H1   HOH     1       1.716  26.282  11.239  1.00  0.00      1DYZ H
HETATM    2 O    HOH     1       2.439  25.795  11.634  1.00  0.00      1DYZ O
HETATM    3 H2   HOH     1       3.140  26.440  11.729  1.00  0.00      1DYZ H
END   

Usage of pdb2adf

The program works interactively. Given below in bold are the parts that the user has to type. In cases where the user agrees with the suggestion given by the program, the user can press the Enter key indicated with Enter.

                                  P D B 2 A D F - program
                                     version 2008.01
                              Written by: Marcel Swart, 2008

                         This program uses AMBER parameter files
                            see: http://amber.scripps.edu
Please give name of PDB-file

hoh.pdb.txt

Do you want a logfile to be written (Y/n) ?

Enter

 read fragments

Data Processed:
     Nat:          3
    Nmol:          1
 NChains:          0

Please wait, making connection tables
Now finding nearby atoms
Assigning chain ID to all residues
Completing residues for which only option is available

Refinding nearby atoms (including atoms added in residue completion)

 - - - - - - - - - - - - - - - - - - - - - - - - - - -
           Making Choices for Chain      0
 - - - - - - - - - - - - - - - - - - - - - - - - - - -

Completing residues with multiple options available, and solvent molecules

Checking positions of newly added atoms

Making choice for which molecules should be QM, which MM

Residues belonging to chain  0

Solvent molecules (SOL/HOH) belonging to this chain:
    1

Give the number of the molecule to be put in QM region (or 'c' to continue):

1

Putting HOH    1 in QM region

Give the number of the molecule to be put in QM region (or 'c' to continue):

c

Do you want to add solvent to your system (Y/n) ?

Enter

Solvent (box) available:
   1:  HOH     HOH  Water molecule
   2:  MOH     MOH  Methanol molecule
   3:  CHL     CHL  Chloroform molecule

2

Reading contents of solvent box p2abox.MOH

Box Shape options:
 1  Spherical box
 2  Cubic box
Make a choice:

1

Writing inputfile for chain   0

Using total charge   0.0 and total spin   0.0

Maximum atomic distance (Angs) from center        0.92
Give boxsize (def.:    15.00 Angs)

14.0

Using BOXSIZE value of  14.0000
Adding atoms for box     1  Added (Box):    84 (Total):    84  Excl. (1):   660  Excl. (2):     6
Adding atoms for box     2  Added (Box):   102 (Total):   186  Excl. (1):   642  Excl. (2):     6
Adding atoms for box     3  Added (Box):   102 (Total):   288  Excl. (1):   642  Excl. (2):     6
Adding atoms for box     4  Added (Box):   108 (Total):   396  Excl. (1):   642  Excl. (2):     0
Adding atoms for box     5  Added (Box):   120 (Total):   516  Excl. (1):   630  Excl. (2):     0
Adding atoms for box     6  Added (Box):    96 (Total):   612  Excl. (1):   654  Excl. (2):     0
Adding atoms for box     7  Added (Box):   108 (Total):   720  Excl. (1):   642  Excl. (2):     0
Adding atoms for box     8  Added (Box):   102 (Total):   822  Excl. (1):   642  Excl. (2):     6

Inputfile(s) written, everything processed, work has been done.
Thank you for using the PDB2ADF program.

================================
Normal ending of PDB2ADF program
================================

Contents of the hoh.pdb2adf file generated by pdb2adf

The file is not given completely, since it contains more than 800 atoms. Note that the old style QM/MM input format is used if the environment variable SCM_PDB2ADF is set to OLD. This example uses the default NEWQMMM format.

#! /bin/sh

$ADFBIN/adf << eor
TITLE QM/MM calculation setup by pdb2adf: M.Swart, 2008

Symmetry NOSYM

EPRINT
  SFO NOEIG NOOVL
END

XC
 GGA Becke-Perdew
END

GEOMETRY
 iterations 100
 CONVERGE grad=1.0e-3 rad=0.1
END

BASIS
 type TZP
 core small
END

SCF
 diis ok=0.01
 Converge 1.0e-5 1.0e-5
 Iterations 99
END

CHARGE   0.0

ATOMS
     1 O      2.4390    25.7950    11.6340
     2 H      1.7160    26.2820    11.2390
     3 H      3.1400    26.4400    11.7290
     4 C    -10.0667    22.2493    11.7437
     5 H    -10.2077    21.5053    10.9597
     6 H    -10.5047    21.8683    12.6667
     7 H    -10.5167    23.2103    11.4977
     8 O     -8.7387    22.3983    12.0617
     9 H     -8.3007    22.6943    11.2607
    10 C     -0.2827    19.0253     2.2847
    11 H     -0.5357    18.2063     2.9567
    12 H      0.7633    19.2913     2.4407
    13 H     -0.9267    19.8753     2.5107
    14 O     -0.4997    18.6373     0.9467
    15 H      0.1123    17.9313     0.7287
....
   823 H      5.4711    27.9401    19.5645
   824 O      5.5611    28.7181    17.7095
   825 H      5.2631    27.8621    17.3935
END

QMMM
  FORCE_FIELD_FILE $ADFRESOURCES/ForceFields/amber95.ff
  ELSTAT_COUPLING_MODEL=1
  ! ---------------------------------------------------
  ! use NEWQMMM developer implementation
  !
  ! NOTE:
  ! the rest of QMMM block is only for use with NEWQMMM
  ! ---------------------------------------------------
  NEWQMMM
  QMMM_PREOPTIMIZE

  QM_INFO
        1  OW    QM   -0.834000  HOH     1  O         2      3
        2  HW    QM    0.417000  HOH     1  H1        1
        3  HW    QM    0.417000  HOH     1  H2        1
  SUBEND

  QMMM_INFO
        1  OW    QM   -0.834000  HOH     1  O         2      3
        2  HW    QM    0.417000  HOH     1  H1        1
        3  HW    QM    0.417000  HOH     1  H2        1
        4  CT    MM    0.116600  MOH     1  C1        5      6      7      8
        5  H1    MM    0.037200  MOH     1  HC1       4
        6  H1    MM    0.037200  MOH     1  HC2       4
        7  H1    MM    0.037200  MOH     1  HC3       4
        8  OH    MM   -0.649700  MOH     1  O1        4      9
        9  HO    MM    0.421500  MOH     1  HO1       8
       10  CT    MM    0.116600  MOH     2  C1       11     12     13     14
       11  H1    MM    0.037200  MOH     2  HC1      10
       12  H1    MM    0.037200  MOH     2  HC2      10
       13  H1    MM    0.037200  MOH     2  HC3      10
       14  OH    MM   -0.649700  MOH     2  O1       10     15
       15  HO    MM    0.421500  MOH     2  HO1      14
....
      823  H1    MM    0.037200  MOH   137  HC3     820
      824  OH    MM   -0.649700  MOH   137  O1      820    825
      825  HO    MM    0.421500  MOH   137  HO1     824
  SUBEND
END
eor