System and method for structure-based drug design that includes accurate prediction of binding free energy

A method for building molecules for binding at a receptor site, comprising the steps of:

• (a) evaluating a receptor site for a molecular make up of at least a portion of the receptor site to which a molecule being grown will bind and generating at least a coordinate of at least a portion of the receptor site to which the molecule being grown will bind, and outputting, at least with respect to the molecular make up of the receptor site, the coordinate of the portion of the receptor site to which the molecule being grown will bind;
• (b) estimating the free energy of the molecule being grown using knowledge based potential data to estimate free energy and outputting the estimated free energy, comprising the substeps of:
• (1) selecting an interaction radius of a length that negates adverse affects from solvent entropy effects;
• (2) developing an accessible listing regarding known structures of molecule/receptor complexes that includes at least molecular fragment
• (3) developing an accessible listing regarding atom types based on classes of molecular interactions that include at least atom types predicted to interact in the molecule/receptor complex of the molecule being grown and the receptor site;
• (4) generating an accessible listing of molecular interactions for known structures of molecule/receptor complexes based on the selection at substep (1) and the listing developed at substeps (2) and (3);
• (5) selecting a reference state for the molecule being grown that negates the adverse effects from solvent and configurational entropy effects;
• (6) developing a quasichemical approximation for the molecule that is being grown;
• and
• (7) generating the estimated free energy of the molecule being grown based on the listing developed at substep (4), the selection at substep (5), and the quasichemical approximation developed at substep (6);
• and
• (c) building a molecule for binding to the receptor site using the outputs from steps (a) and (b), with the building step including building the molecule by selecting molecular fragments at orientations that will result in free energy estimates for the molecule that may be higher than a lowest free energy estimate possible for the molecule.

The method for estimating free energy in a molecule being grown for connection at a receptor site in a molecule/receptor complex, comprising the steps of:

• (a) selecting an interaction radius of a length that negates adverse affects from solvent entropy effects;
• (b) developing an accessible listing of known structures of molecule/receptor complexes that includes at least atom coordinates and corresponding chemical elements;
• (c) developing an accessible listing of atom types based on clases of molecular interactions that include at least atom types predicted to interact in the molecule/receptor complex of the molecule being grown and the receptor site;
• (d) generating an accessible listing of molecular interactions for known structures of molecule/receptor complexes based on the selection at step (a) and the listing developed at steps (b) and (c);
• (e) selecting a reference state for the molecule being grown that negates the adverse effects from solvent and configurational entropy effects;
• (f) developing a quasichemical approximation for the molecule that is being grown;
• and
• (g) generating the estimated free energy of the molecule being grown based on the listing developed at step (d), the selection at step (e), and the quasichemical approximation developed at step (f).

A method for building molecules for binding at a receptor site, comprising the steps of:

• (a) evaluating a receptor site for a molecular make up of at least a portion of the receptor site to which a molecule being grown will bind and generating at least a coordinate of at least a portion of the receptor site to which the molecule being grown will bind;
• (b) loading at least hydrogen ("H") atoms at least at the portion of the receptor site to which the molecule being grown will bind;
• (c) randomly selecting an H atom of the H atoms loading at the receptor site;
• (d) randomly selecting a molecular fragment from the predetermined group of molecular fragments and randomly selecting an H atom from the molecular fragment;
• (e) forming a bond at the selected H atoms selected at steps (c) and (d);
• (f) orienting the molecular fragment with respect to the bond formed at step (e) so that the molecule being grown has low free energy;
• (g) estimating free energy for the molecule being grown using knowledge-based potential;
• and
• (h) comparing the estimated free energy at step (g) with a free energy estimate of a molecule bound at the receptor site before the fragment was bonded to the receptor site at step (e) to determine if the estimated free energy at step (g) is more or less than a prior estimated free energy for the molecule before the fragment was bonded to he receptor site at site (e), and accepting the addition and orientation of the molecular fragment if the estimated free energy is less and accepting the addition and orientation of the first molecular fragment according to a probability determined in a predetermined manner if the estimated free energy is more.

A method for growing molecules for binding at a receptor site, comprising the steps of:

• (a) evaluating a receptor site for a molecular make up of at least a portion of the receptor site to which a molecule being grown will bind and generating at least a coordinate of at least a portion of the receptor site to which the molecule being grown will bind;
• (b) loading at least hydrogen ("H") atoms at least at the portion of the receptor site to which the molecule being grown will bind;
• (c) randomly selecting an H atom of the H atoms loading at the receptor site;
• (d) randomly selecting a molecular fragment from the predetermined group of molecular fragments and randomly selecting a H atom from the molecular fragment;
• (e) forming a bond at the selected H atoms selected at steps (c) and (d);
• (f) orienting the molecular fragment with respect to the bond formed at step (e) so that the molecule being grown has low free energy;
• (g) estimating free energy for the molecule being grown using knowledge-based potential;
• and
• (h) comparing the estimated free energy at step (g) with a free energy estimate of a molecule bound at the receptor site before the fragment was bonded to the receptor site at step (e) to determine if the estimated free energy at step (g) is more or less than a prior estimated free energy for the molecule before the fragment was bonded to the receptor site at step (e), and accepting the addition and orientation of the molecular fragment if the estimated free energy is less and accepting the addition and orientation of the first molecular fragment according to a probability determined in a predetermined manner if the estimated free energy is more.