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




Thermodynamic 2D Description of H-bonding

That hydrogen bonding plays an important role in chemical processes, including enzymatic catalysis, is well known. This observation has inspired many approaches to take such interactions into account in QSAR and molecular modeling investigations. Up to now, the number of atoms that can participate in H-bonding has been used most often as the simplest descriptor of this interaction . Of course, such indicator variables are primitive, and do not reflect the influence of substituents near H-bond acceptors or donors.

On the basis of their own experimental measurements and published thermodynamic data on hydrogen bonding, Raevsky et al found that H-bond strength depends on the nature of substituents near the atoms that participate in bond formation [1,2]. The enthalpy values in measured H-bond complexes cover a wide range; they vary from 0.9 to 15.9 kcal/mol for "O-H…O"; from 0.1 to 19.1 kcal/mol for "O-H…N"; from 0.5 to 10.5 kcal/mol for "N-H…O", and from 0.5 to 11.5 kcal for "N-H…N".

Using a multiplicative approach to describe experimental thermodynamic parameters for about 14 thousand H-bond complexes, Raevsky et al. [3-5]: created uniform enthalpy (DH) and free energy (DG) H-bond scales:

D H = k1(kcal/mol)EaEd       (1)

D G = k2(kcal/mol)CaCd+ k3(kcal/mol)       (2)

where Ea is the enthalpy acceptor factor , Ed is the enthalpy donor factor , Ca is the free energy acceptor factor, and Cd is the free energy donor factor. Later Raevsky proposed also overall H-bond acceptor scale based on correlation of experimental property values with free energy acceptor factors to estimate an ability of different weak neutral and ionized acceptor atoms. H-bond factors can quantitatively characterize the acceptor or donor strength of any atom in any molecule.

Large data bases for these factor values as well as a program to estimate such values for atoms in organic molecules are contained in the program HYBOT (HYdrogen BOnd Thermodynamics) for Windows 95/98/2000/NT or for UNIX or for LINUX [6]. Besides HYBOT also calculates descriptors connected with steric and electrostatic interactions. Fig. below presents the results of calculation of molecular and atomic descriptors by means of HYBOT (version 2.3, 2003).

NN

Desriptor

Symbol

1

Molecular polarizability

α

2

Maximal positive atomic charge in a molecule

Max Q+

3

Maximal negative atomic charge in a molecule

Max Q-

4

Maximal enthalpy H-bond  acceptor factor

Max Ea

5

Maximal enthalpy H-bond  donor factor

Max Ed

6

Maximal free energy H-bond  acceptor factor

Max Ca

7

Maximal free energy H-bond  donor factor

Max Cd

8

Maximal overall H-bond  acceptor factor

Max Ca(o)

 

Sum of all positive atomic charges in a molecule

ΣQ+

10

Sum of all negative atomic charges in a molecule

ΣQ-

11

Sum of absolute values for all atomic charhes

Σ|Q|

12

Sum of all enthalpy H-bond acceptor factors

ΣEa

13

Sum of all free energy H-bond acceptor factors

ΣCa

14

Sum of all overall H-bond acceptor factors

ΣCa(o)

15

Sum of all enthalpy H-bond donor factors

ΣEd

16

Sum of all free energy H-bond factors

ΣCd

17

Sum of absolute values of all enthalpy H-bond acceptor and donor factors

Σ|E|

18

Sum of absolute values of all free energy H-bond  acceptor and donor factors

Σ|C|

19

Sum of absolute values of all overall H-bond acceptor factors and all free energy H-bond donor factors

ΣCa(o)+Σ|Cd|

References:

  1. Raevsky, O.A. and  Novikov, V.V., H-Bond Description by Multiplicative Apprtoach, Khim-Pharm. Z. (russ.),  16 (1982) 583; Chem.Abst. 97(1982) 33128d;
  2. Raevsky, O.A., Avidon V.V.  and  Novikov, V.P., H-bond descriptors for QSAR, Khim-Pharm. Z.(russ),  16 (1982) 968; Chem.Abst. 97 (1982) 181323d.
  3. Raevsky, O.A. in Hadzi, D. and Jerman Blazic, B. (eds.), QSAR Description of Molecular Structure, in "QSAR in Drug Design and Toxicology", Elsevier, Amsterdam, 1987, pp.31-36.
  4. Raevsky, O.A., Grigor’ev, V.Ju., Kireev, D.B. and Zefirov, N.S., Complete Thermodynamic Description of H-bonding in the framework of multiplicative approach,  Quant. Struct.-Act.Relat., 11 (1992) 49.
  5. Raevsky, O.A., Molecular Structure Descriptors in the Computer-Aided Design of Biologically Active Compounds, Russ. Chem. Rew , 68 (1999) 437.
  6. Raevsky, O.A., Grigor’ev, V.Ju. and Trepalin, S.V., HYBOT  program package, Registration by Russian State Patent Agency N 990090 of  26.02.99; Raevsky, O.A., Skvortsov, V.S., Grigor’ev, V.Ju. and Trepalin, HYBOT in UNIX/LINUX  , Registration by Russian State Patent Agency N 2002610496 of 05.02.02. ( raevsky@ipac.ac.ru)