target based rational drug design l.
Skip this Video
Loading SlideShow in 5 Seconds..
Target Based Rational Drug Design PowerPoint Presentation
Download Presentation
Target Based Rational Drug Design

Loading in 2 Seconds...

play fullscreen
1 / 38

Target Based Rational Drug Design - PowerPoint PPT Presentation

  • Uploaded on

Target Based Rational Drug Design

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about 'Target Based Rational Drug Design' - Samuel

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
target based rational drug design

Target Based Rational Drug Design

Objectives: Current drug discovery efforts are based on rationally identifying chemical compounds that will bind to therapeutic target molecules such as proteins.  These efforts are based on the 3D structure of the target molecule and use a variety of computer-based molecular modeling techniques to exploit the 3D structural information.  Upon completion of these lectures, the student will be familiar with molecular modeling approaches that are used in target based drug design and how they facilitate bringing drugs to market.

goals of target based drug design
Goals of target-based drug design
  • Understand atomic details of drug binding strength and specificity

2) Identify or create novel molecules to bind to a selected target and elicit a biological response via de novo drug design or database searching techniques

3) Optimize the therapeutic index of an already available drug or lead compound

Computer-Aided Drug Design (CADD)

molecular modeling
Molecular Modeling

3D representation of molecules based on graphic or mathematical representations of chemical structures

Qualitative versus quantitative models

Relate 3D structure/conformation to energy via a mathematical equation

quantum mechanics versus molecular mechanics
Quantum mechanics versus molecular mechanics

Quantum mechanicsi) treat electrons explicitlyii) limited to 100 atoms

Molecular mechanicsi) atom is smallest particle ii) every atom interacts with every other atom iii) allows systems containing 10s of thousands to millions of atoms to be studied

molecular mechanics empirical energy function
Molecular mechanics(empirical energy function)

Determine total potential energy, Vtotal, as a function of conformation

Internal energybondsvalence anglesdihedrals

External energy

(nonbond, interaction)van der Waals (VDW, LJ)electrostatic

potential energy vs bond length
Potential energy vs. bond length

Similar equation for valence angles:

how can a mathematical model be related to the type of chemical being studied
How can a mathematical model be related to the type of chemical being studied?


single bond: bo of 1.54 Å,Kb of 100 kcal/mole

double bond: bo of 1.34 Å, Kb of 150 kcal/mole

triple bond: bo of 1.20 Å, Kb of 200 kcal/mole


Treatment of hydrogen bonding via partial atomic charges, where the sum of the electrostatic interactions between each atom (i.e. O with H, O with N, C with H and C with N) yields an interaction energy appropriate for a hydrogen bond (~-1 to -5 kcal/mol).

energy minimization
Energy Minimization

Force vs bond length

hypothetical energy surface for conformational changes in a protein
Hypothetical energy surface for conformational changes in a protein

Acetic acid: 1 degree of freedom leading to 2 conformations

Proteins: ~10 conformations/amino acid leading to 10n conformations where n is the number of amino acids. Note global minimum vs. local energy minima

Therefore: require methods to more rigorously sample conformational space that allow energy barriers to be overcome

information required for calculation of energy for minimization and for molecular dynamics
Information required for calculation of energy, for minimization and for molecular dynamics

Energy 1) atom positions

Energy minimization 1) atom positions 2) forces

Molecular dynamics simulation 1) atom positions 2) forces 3) atomic velocities: assign randomly to correspond to desired temperature

computer aided rational drug design lead identification
Computer-Aided Rational Drug Design:Lead Identification
  • Dock known chemicals from an in silico database into the receptor target site (database searching or screening).

2) Design of novel compounds to fit into receptor target site (de novo design).


Geometric arrangement of types of functional groups that are required for “activity” to the active site of HIV protease

final pharmacophore
Final Pharmacophore

Many molecular could fulfill the above pharmacophore: Identify such molecules via database screening or de novo design

  • a) Prepare a 3D lattice of grid points encompassing the binding site
  • b) Determine interaction energy of different types of functional groups with binding site at each grid point.
  • c) Select favorable interaction sites
  • d) Determine relative spatial orientations of the selected interaction sites
  • Types of functional groups and their positions yields a pharmacophore. Grid may also be used to calculate interaction energies of ligands with a receptor.
database searching
Database searching

Identify novel lead compounds from databases of known compounds

Chemical databasesCambridge Structural DatabaseChemical AbstractsNational Cancer InstitutePharmaceutical CompaniesSoP Computer-Aided Drug Design Center


1) Determine solvent accessible (Connelly) surface of the binding pocket by rolling a sphere the size of a water molecule along the surface


2) Generate "negative" image of receptor based on spheres that is complementary to the receptor surface generated in step 1

3) Determine sphere-sphere distances of the negative receptor image

4) Convert sphere-sphere distances to possible atom-atom distances

5) Compare atom-atom distances with actual atom distances of chemical compounds in a database

6) Select ligands with greatest overlap for further studies (perform conformational search and energy minimization to obtain low energy structure of ligand in binding site).

7) Calculate ligand-receptor interaction energies (via GRID, see below) of selected ligands

de novo drug design groupbuild
De Novo Drug Design: Groupbuild

Build novel compounds that are complementary to a target binding site on a protein via “random” combination of small molecular fragments.

Example fragments: formic acid, formaldehyde, formamide, amine, benzene, cyclohexane, cyclopentane, ethane, ethylene, water, methanol, methane, sulfone, thiophene

groupbuild structure generation
Groupbuild Structure Generation

1) Define scoring grids for the binding siteElectrostaticVDW (steric)

2) Structure generationa) Predock a "core" (starting) fragment into the binding siteb) Structure build up (extension of compound)c) Randomly select one of top 25% of the new structuresd) Iterate over steps B and C (i.e. keep building) until termination criteria are fulfilled

3) Selection of built structures for synthesis and analysis

analyses of built structures
Analyses of built structures
  • Visual examination of structures for chemical feasibility

B) Identify specific positions of certain functional groups etc. that my be related to a known pharmacophore

C) Database screening for similar compounds in chemical databases (avoid synthesis!).

limitation of database screening and de novo design
Limitation of database screening and de novo design

Rigid geometry of receptor and ligand

Ligands often treated as flexible

Multiple conformations of receptor can partially overcome rigid representation

Inherent assumptions and simplifications in molecular mechanics

steps involved in drug receptor interactions
Steps involved in drug-receptor interactions

1) diffusion controlled encounter rate

2) initial Michaelis complex

3) desolvation of both inhibitor and binding site

4) conformational changes of both inhibitor and binding site upon binding

5) correct orientation between drug and receptor binding site

All the above cannot be calculated, therefore, need a trick to make the problem computational feasible

how is an alchemical perturbation performed
How is an alchemical perturbation performed?

Single bond: bo of 1.54 Å, K of 100 kcal/mole Double bond: bo of 1.34 Å, K of 150 kcal/mole

Gradually change the equilibrium and force constants. For example the equilibrium bond length of 1.54 goes to 1.34 Å in 10 steps of 0.02 Å and the force constant of 100 goes to 200 kcal/mole in 10 steps of 10 kcal/mole.

Sum over the “change” in energy in each step to get the total free energy difference between C-C and C=C.

data for binding to a hypothetical receptor
Data for binding to a hypothetical receptor

G = Geq,2 – Geq,1 = Greceptor - Gsolution

free energy component analysis
Free Energy Component Analysis

Since energy calculations involve contributions from every atom in the system, the individual contributions of functional groups on the drug to drug-receptor interactions or solvation energies can be estimated. This allows for a detailed understanding of the relationship of structural changes to changes in binding to be determined. From this information atomic details of drug solubility and drug-receptor interactions are obtained that can be used to better interpret experimental data and make predictions on how the change a drug's structure to improve it's biological activity.

lead compound optimization via free energy perturbation calculation
Lead compound optimization via free energy perturbation calculation

Improved interpretation of experimental data allowing for rational design of compounds Solvation vs. receptor binding Contributions of amino acids and/or drug functional groups to binding

Prediction of modifications to improve activity