System and method for improved computer drug design

• US 7,801,685 B2
• Filed: 08/19/2005
• Issued: 09/21/2010
• Est. Priority Date: 08/19/2004
• Status: Active Grant

• Current Assignee: Drug Design Methodologies LLC
• Original Assignee: Drug Design Methodologies LLC
• Inventors: Ho, Chris Meichung Wang
• Primary Examiner(s): Moran; Marjorie
• Assistant Examiner(s): Skibinsky; Anna
• Application Number:
  US11/208,189
• Publication Number:
  US 20060041414A1
• Time in Patent Office: 1,859 Days
• Field of Search: None
• US Class Current: 702/19
• CPC Class Codes: G16B 15/00 ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
• CPC Class Codes: G16B 15/30 Drug targeting using structural data; Docking or binding prediction
• CPC Class Codes: G16C 20/40 Searching chemical structures or physicochemical data
• CPC Class Codes: G16C 20/50 Molecular design, e.g. of drugs

Abstract

A system and method for computer-aided drug design employs adaptive sampling and iterative fitting of multi-atomic subunits. The iterative fitting is performed by successive perturbation of the location of the multi-atomic subunits in directions that reduces potential energy.

Claims

1. A method for computer-aided molecule design, the method comprising:

executing an adaptive sampling function and rapid minimization function, the rapid minimization function employing iterative fitting of multi-atomic subunits of a molecule, the multi-atomic subunits including subunits having a plurality of non-hydrogen atoms, the multi-atomic subunits being fitted by perturbing the location of the multi-atomic subunits, while keeping the coordinates of the constituent atoms of a given subunit fixed with respect to all of the other constituent atoms of that subunit, in a real-space direction that reduces potential energy;

wherein the rapid minimization function perturbs the real-space position of the multi-atomic subunits by calculating a scaled force vector for the subunits'"'"' constituent atoms, and then solving for new coordinates for the constituent atoms using root mean square fitting onto the new coordinates, while still keeping the coordinates of the constituent atoms of each subunit fixed with respect to one another, the scaled force vectors resulting from non-bonded interactions with a receptor and from bonded interactions with neighboring subunits, and the scaled force vectors being scaled according to molecular mechanics contributions to resolve steric violations first, geometries that stretch or compress bond lengths second, bond angle distortions third, and alterations in torsion angles last; and

wherein the adaptive sampling function and rapid minimization function are executed on a computer.

2. The method of claim 1 further comprising determining the optimal six-dimensional coordinates for units using RMS fitting.

3. The method of claim 1 wherein force vectors are also determined for electrostatic sources of potential energy.

4. The method of claim 1, wherein the multi-atomic subunits include only non-rotatable bonds.

5. The method of claim 2, wherein the multi-atomic subunits are bonded to each other only by rotatable bonds.

6. The method of claim 1, further comprising performing an iterative conformational searching function having minimization function as an element of conformational searching iterations.