Image_logo.bmp (1885 bytes) MOL ECULAR PR OPERTIES O PTIMIZATION PROJECT




Thermodynamic 3D  Description of  H-bonding.

Most modern procedures for three-dimensional (3D) computer-assisted drug design use force-field calculations to take H-bonding into account. These calculations are used in many research areas including conformational analysis, pharmacophore identification, ligand docking, de novo ligand design, comparative molecular field analysis (CoMFA), and the identification of favorable binding sites from molecular interaction fields .

On the basis of a sophisticated and  parameterized method called GRID, Goodford et al. [1-3] made efforts to improve the calculation of the H-bond contribution to the interaction energy of ligand-macromolecular complexes. The GRID method is parameterized by fitting experimental X-ray data from protein and small-molecule crystals, and is designed to calculate the interactions of a probe (a small molecule such as water or ammonia) and a macromolecular system. An 8-6 function was adopted, and was found to give satisfactory result [2]:

Er = C/r8 - D/r6     (1),

where C = -3Emr8A8/mol; D = -4Emrm6A6/mol ; r is the distance between the acceptor atom and the donor heavy atom; Em is optimum H-bond enthalpy in kcal/mole; rm is the optimum H-bond length.

However there are only three fixed optimum H-bond lengths in the GRID framework: 2.8 A for “O…O”, 3.0 A for “N…O” and 3.2 A for “N…N” . And the wide intervals of enthalpy values for each those type of H-bond complexes (see previous Chapter)  permit to suppose any dependence between energies and distances even in the case of optimum arrangement of atoms participating in H-bonding. Hence, to create a realistic platform for the quantitative 3D description of H-bonding, one must recognize that optimum H-bond energies (Em) depend on optimum H-bond distances (rm) for the different types of H-bonds in complex formation. This our study is presented in [4].

First of all  we selected from Cambridge Structural Database 58 "ideal" H-bond complexes where the angles of the donor heavy atom, hydrogen and acceptor atoms were in the range of 173° to 187° and angles of the acceptor nucleus, the lone pair of electrons and the hydrogen atom were in the range of 170° to 190°. There were 13 such “O-H…O” complexes, 19 ”O-H…N” or “N‑H…O” complexes, and 26 “N-H…N” complexes among the 58 "ideal" representatives.

Because there are significant differences between the covalent and van der Waals radii of oxygen and nitrogen atoms, we decided to look for relationships between optimum energy (Em) and optimum distances (rm) in those three subsets. Such relationships are not expected to be linear because energy values approach zero as distances increase. In this study we used a sigmoid function:

Em = k4/(1 + 10k5+k6rm)    (2)

Values of k4 in eq. (2) may be fixed, and are limited by the maximum values for the enthalpies. Rearranging eq. (2) and simplifying the constants gives eq. (3):

rm = k7log[(k4-Em)/(Em)] + k8      (3)

Using eq. (3), we can now calculate rm from Em values. Thus, for each specific pair of atoms participating in an H-bond, its H-bonding potential can be calculated. The potentials fot 58 “ideal” complexes are presented on the fig.:

Hence, a scheme to calculate the H-bonding energies of atoms in a three-dimension arrangement can be carried out by following steps:

  1. Use HYBOT’s acceptor and donor enthalpy factors to estimate the optimum H-bonding energy.

  2. Estimate the optimum H-bonding distance corresponding to the calculated energy optimum by means of eq. (3) .

  3. Use eq. (1) or other functions to take into account the distance and angle deviations of the atoms involved.

An example of such calculation for a model system is presented in [4]. This represents the interactions of bases in double helix of RNA in A form.

These calculations were carried out by means of 3D HYBOT in UNIX, a new program, by considering only H-bonding interactions and Sybyl . We believe it is possible that this approach to H-bond calculations will not be too difficult to integrate with existing and future force-field schemes to produce a superior method to estimate intermolecular interactions in any three dimensional framework.

3-D hydrogen bond descriptors.

 

NN

Descriptor

Symbol

1

Van der Waals’acceptor surface area in A2  which is proportional to H-bond enthalpy factors  of acceptor atoms. n is number of acceptors in the molecule of interest. . So is a surface sphere with a radius of 1.36 A (Osp3)

2

Van der Waals’acceptor surface area in A2 which is proportional to H-bond overall free energy factors Ca(o) of acceptor atoms.

3

Van der Waals donor surface area in A2 which is proportional to H-bond enthalpy factors of donor atoms.n is number of donors in a molecule, , SH is a surface sphere with a radius of 1.08 A  (H atom).

 

4

Van der Waals donor surface area in A2 which is proportional to H-bond free energy factors of donor atoms.

 

5

Surface area around a molecule in A2 where interactions of acceptor atoms of a molecule with a H-bond donor probe have been optimumly placed and which is proportional to product of H-bond enthalpy factor absolute values of those atoms. Ed(probe) is the enthalpy factor of the probe H-bond donor, , Srm is the surface area of sphere with a radius of rm = 2.45 A for the strongest H-bonding

6

Surface area around a molecule in A2 where interactions of acceptor atoms of a molecule with H-bond donor probe have been optimumly placed and which is proportional to product of H-bond (overall) free energy factor absolute values of those atoms.

7

Surface area in A2 around a molecule where interactions of donor atoms of a molecule with H-bond acceptor probe have been optimumly placed and which is proportional to product of its H-bond enthalpy factor absolute values.

8

Surface area in A2 around a molecule where interactions of donor atoms of a molecule with H-bond acceptor probe have been optimumly placed and which is proportional to product of its H-bond free energy factor absolute values.

9

Surface integral for enthalpy values (kcal/M*A2) of interactions between acceptor atoms of a molecule and a donor probe on the surface OEASA .

10

Surface integral for enthalpy values (kcal/M*A2) for interactions between donor atoms of a molecule and an acceptor probe on the surface OEDSA

 

References:

1.     Goodford, P.J., J.Med.Chem., 28 (1985) 849.

2.     Wade, R.C., Clark, K. and Goodford, P.J., J.Med.Chem., 36 (1993) 140,148.

3.      Goodford, P.J., J.Chemometrics, 10 (1996) 107.

4.     Raevsky O.A. and Skvortsov, V.S., 3D Hydrogen Bond Thermodynamics (HYBOT) Potentials in Molecular Modelling, Comp.-Aided Mol.Des., 16, 1 (2002).