Next: Turning Energy Terms On Up: Energy Function Previous: Non-crystallographic Symmetry Interactions

The Explicit Hydrogen-Bond Term

The explicit hydrogen-bond energy term $E_{HBON}$ is used in certain parameter sets, in particular the ones that deal with nucleic acids. In the more recent protein parameter files, this term is zero, and the effect of hydrogen bonds is taken implicitly into account by appropriate parameterization of the partial charges and van der Waals parameters. The explicit hydrogen-bond term has the form

$\displaystyle E_{HBON}$	$\textstyle =$	$\displaystyle \sum_{i<j} (\frac{A}{r^{a}_{AD}} - \frac{B}{r^{b}_{AD}})$	(4.25)
		$\displaystyle cos^m(\theta_{A-H-D}) cos^{n} (\theta_{AA-A-H}) SW(r_{AD}^{2},r_{h_{on}}^{2}, r_{h_{off}}^{2})$
		$\displaystyle SW[cos^{2} (\theta_{A-D-H}), cos^{2}(\theta_{h_{on}}), cos^{2} (\theta_{h_{off}})]$

involving atoms AA, A, H, D (acceptor antecedent, acceptor, hydrogen, and donor heavy atom), where

and

are positive integers and

is the switching function defined in Eq. 4.9.

The switching function SW maintains a continuous and smooth function to the cutoff point so that derivatives are well defined at the cutoff point. The sum in Eq. 4.25 extends over all pairs of donors and acceptors (DONOr and ACCEptor statements; see Section 3.1.1) for which the donor-acceptor distance and the donor-hydrogen-acceptor angle satisfy the criteria specified by DCUT and ACUT (see Section 3.2.1). AEXP assigns the exponent, REXP assigns the exponent, HAEX assigns the exponent, and AAEX assigns the exponent; and are defined by the HBON statement (see Section 3.2.1). The remaining switching parameters are provided by the hbonds statement. The computational time for searching for hydrogen bonds that satisfy the cutoff criteria is reduced by introducing an approximation, storing the atomic pair indices () that satisfy the cutoff criteria in a list that is updated only when any atom has moved more than the value specified by TOLErance.

Xplor-NIH 2024-09-13