Medicinal Chemistry Unit 5
Medicinal Chemistry Unit 5
Medicinal Chemistry Unit 5
Learning Outcomes:
1.Student will learn about the different techniques of
drug designing.
2. Student will learn about the Relation of drug and its
Activity.
Quantitative Structure Activity
Relationships (QSAR)
CONTENTS
INTRODUCTION
HYDROPHOBICITY OF MOLECULE
HYDROPHOBICITY OF SUBSTITUENTS
ELECTRONIC EFFECT
STERIC EFFECT
HANSCH EQUATION
C R AIG PLOT
REFERENCES
INTRODUCTION
Drug designing....
Traditional
Lead compound
Analogue molecules designing new molecule
Eg; salicylic acid and aspirin
PHYSICOCHEMICAL PROPERTIES
Log (1/C)
. . .
.
.. . . .
0.78 3.82 Log P
HYDROPHOBICITY
Log (1/C)
. . .
log(1/c) = 0.075 log P + 230
(n= 42, r= 0.960 s= 0.159)
.
.. . . .
0.78 3.82 Log P
If the partition coefficient is the only
factor influencing biological activity, the
parabolic curve can expressed by the
equation
Log (1/C)
log(1/c) = -K1 (log P)2 + K2 log P + k3
The equation is
x= logPx-logPH
A possitive value shows that the substituent is more
hydrophobic than hydrogen
Example:
Cl
Log P(theory) = log P(benzene) +
2
O
NH2
meta c h l o r o b e n z a m i d e
Have an effect on how easily a drug can pass through the cell
membrane or how strongly it can interact with a binding site.
+
COO +
- H
COOH
[PhCO -2]
K H = Dissociationconstant = [PhCO 2H]
X= electron withdrawing group (e.g. NO2,)
X = electron
withdrawing X X
group CO2H CO2 + H
Charge is stabilised by X
Equilibrium shifts to right
KX > KH
X = log
K X = logK
X - logK H
KH
Positive value
X= electron donating group (e.g. CH3)
X X +
COOH COO
-
+ H
Charge destabilised
Equilibrium shifts to left
KX < KH
Negative value
EXAMPLES: p (NO2) m (NO2)
meta-Substitution
O
N
O
para-Substitution
O O O O O O O O
N N N N
e-withdrawing
(inductive +
resonance effects)
DRUG DRUG DRUG DRUG
o value depends on inductive and resonance effects
The constants σ,R and F can only be used for aromatic substituents
Aliphatic electronic substituents
The bulk, size and shape of a drug will influence how easily it can
approach and interact with binding site.
A bulky substituents may act like a shield and hinder the ideal
interaction between a drug and its binding site.
Disadvantages
(n 2 - 1) mol. wt.
MR = x
(n 2 - 2) density
B4 B3
O B
3
B2
C
B1
O H O C O
B
4
L
HANSCH EQUATION
Log 1 C - k1(logP)2 + k 2 logP + k 3 + k 4 Es + k 5
Craig Plot
Craig plot shows values for 2 different physicochemical
properties for various substituents
. . + 1.0
. .
CF 3 SO 2
. .. .. CN
.75
.50
NO 2
SF 5
.
SO 2 NH 2 CH 3S O 2
CF 3
..
CH 3 CO
.
CONH 2
. .25
OCF 3
. . -.8 -.4
CO 2 H
.4
Cl Br I
.
-2.0 -1.6 -1.2 .8 1.2 1.6 2.0
.
F
-
. .CH 3CONH +
. OCH 3
-.25 Me Et
t-B utyl
.
OH
. NH 2
-.50
NMe 2
-.75
-1.0
-
Allows an easy identification of suitable substituents for
a QSAR analysis which includes both relevant
properties
•This technique uses the same reaction conditions with the same
reaction vessels to produce a large range of analogues.
• natural products
•synthetic custom crafted organic molecules made in
small numbers
• analogues of known actives (analogue me-toos)
Advantages
• Specific reactants can be bound to specific beads
• Beads can be mixed and reacted in the same reaction vessel
• Products formed are distinctive for each bead and physically distinct
• Excess reagents can be used to drive reactions to completion
• Excess reagents and by products are easily removed
• Reaction intermediates are attached to bead and do not need to be isolated
and purified
• Individual beads can be separated to isolate individual products
• Polymeric support can be regenerated and re-used after cleaving the
product
• Automation is possible
1. SOLID PHASE TECHNIQUES
Requirements
• A resin bead or a functionalised surface to act as a solid support
• An anchor or linker
• A bond linking the substrate to the linker. The bond must be
stable to the reaction conditions used in the synthesis
• A means of cleaving the product from the linker at the end
• Protecting groups for functional groups not involved in the
synthesis
1. SOLID PHASE TECHNIQUES
For amines.
➢ Boc ( t-butoxycarbonyl )
➢ Fmoc (9-fluorenylmetoxy carbonyl)
Tmsec (2 [ trimethylsilyl ] ethoxycarbonyl)
Linker R H
HO2C NHBoc
O O
R2 H
O O O
NH2 coupling NH NHBoc
R H R H H
R2
O HF
O
Release from OH
aa1aa2aa3 aan NH2
solid support
HO2C aa1aa2aa3 aan NH2
Peptide
5. Identification of structures from mixed
combinatorial synthesis
5.1 Recursive Deconvolution:
• Method of identifying the active component in a mixture
• Quicker than separately synthesising all possible components
• Need to retain samples before each mix and split stage
Example:
Consider all 27 tripeptides synthesised by the mix and split strategy
from glycine, alanine and valine
Gly
Gly
Ala
Ala
Val
Val
Active
H 2N H
O NH2 HN
MeOS O
NH H
O Tryptophan
SOMe
HN
O
Lysine
NH2
NH2
5.2 Tagging
Example NH2
NH2 NH2 aa1
RCHBrCO2H amino acid(aa 1) R'NH2
NH2 NH R HN
Step 1 Tag 1 NH R
Step 2
O Br
O Br
NH2 NH2
NH2 aa2 aa2
aa3 NH2
aa2
O NR'COR"
6. Identification of struc tures LIGHT from
combinatorial synthesis LIGHT
Deprotection CO2H
NHX NH2 NH2 NHX NHX NHX NHX NHX NHX NHX
coupling
NHX NHX NHX NHX NHX NHX NHX NHX NHX NHX NHX NHX
6. Identification of structures from
combinatorial synthesis
6.2 Photolithography - example Y
Y Y Y
repeat
O O Me
Y fluorescent tag
O2N OMe
Target receptor
amino acids X= Nitroveratryloxycarbonyl
7. Planning a Combinatorial Synthesis
7.1 Aims
• To generate a large number of compounds
• To generate a diverse range of compounds
• Increase chances of finding a lead compound to fit a binding site
• Synthesis based on producing a molecular core or scaffold with
functionality attached
Centroid
or scaffold
Substituent
'arms' Binding groups
7. Planning a Combinatorial Syntheses
7.1 Aims
Target molecules should obey Lipinski’s ‘Rule of Five’ for oral
activity
Screen compound
library
RECEPTOR
BINDING
SITE
Privileged scaffolds
- scaffolds which are common in medicinal chemistry and
which are associated with a diverse range of activities
- benzodiazepines, hydantoins, benzenesulphonamide etc
7. Planning a Combinatorial Syntheses
7.2 Scaffolds - examples
R"
R O Me O
N R4
O R2 R3 R HO2C
N R1 4
R' R2
O
N
N N N
X Ar R1 R2 O R5 R3
R3 R
R3 R4
O
R5
• Good scaffolds
•
C N
R1 O C N
R6
Spider like
O
• Low molecular weight
R2 • Variety of synthetic routes available
Dipeptides
7. Planning a Combinatorial Syntheses
7.2 Scaffolds - poor examples
OR5
Spider like and small molecular weight - good
R4O O
OR1
points
R3O
OR2 But multiple OH groups
Glucose Difficult to vary R1-R5 independently
Me
Me
M.Wt. relatively high
Restricts no. of functional groups to keep MWt.<
R1CO 500
R2
Relatively few positions where substituents
Steroid easily added
O
R3
Tadpole like scaffold
H 2N
R2 Restricted region of
N
R
variability
Indole
VIRTUAL SCREENING IN
DRUG DISCOVERY
dsdht.wikispaces.com
Computational
Screening,
Experimental Screening
“Available” Validated
Compounds Targets
Choosing the right molecule
• Goal: to find a lead compound that can be optimized to give a drug candidate
− Estimates vary
• Reymond et al. suggest there are ~1 billion compounds with up to 13 heavy atoms
• High throughput screening allows large (up to 1 million) numbers of compounds to be tested
Screening Database
Virtual
Screening
Leads
HTS
Comb
Library Lead
Optimization
Virtual Screening
Depending upon structural and Bioactvity
data available :
3D
Structure
• One or more actives molecule known
perform similarity searching.
unknown Known
• Several active known try to identify a
common 3D pharmacophore and
Structure then do 3D database search.
Ligand
Based Based • Reasonable number of active and
inactive known train a machine
learning model.
Actives Known Actives and Inactive Known
• 3D structure of protein known use
protein ligand docking.
PhyChem property
Post Process descriptor calculation
and filtration
Cscore, MM/PBSA, Solvation Corrections
Apply protonation at
Potential Lead compounds pH 7.4
Drug Like Properties
Drug-like properties are an integral element of drug discovery projects.
Properties of interest to discovery scientists include the following:
• Structural properties
Hydrogen Bonding, Polar Surface area , Lipophilicity, Shape , Molecular
Weight, Reactivity, pka
• Physicochemical Properties
Solubility, Permeability, Chemical Stability
• Biochemical Properties
Metabolism(Phase 1 and 2) , Protein and tissue binding, transport
• Druglike
-Importantly, exclude chemically reactive functional groups ,‘promiscuous
inhibitors’, ‘frequent hitters’ and warheads
- MW < 500
- cloP < 5
- H-bond donors < 5
- Sum of N and O (H-bond acceptors) < 10
- Polar surface area < 140 A2
- Number of rotatable bonds <= 10
Filtering molecules using structuralproperties
Basic Washing –
• Removing Salts & Unwanted Elements
Filter out cationic atoms: Ca2+, Na+, etc.
Filter out metals:
Sc,Ti,V,Cr,Mn,Fe,Co,Ni,Cu,Zn,Y,Zr,Nb,Mo,Tc,Ru,Rh,Pd,Ag,Cd
Often the salt “filter” = keeping the largest molecule in the sdf entry.
• ALLOWED_ELEMENTS H, C, N, O, F, P, S, Cl, Br, I
• Check proper Atom Types by adding hydrogen and checks if O, N, C valences are
correct.
• Check formal charge
Filter out Reactives (false positives for proteins)
Filter out: Synthesis Intermediates, Chelators
‘warhead’ agents - functional groups which shows high reactivity to proteins due
which there is high attrition rate in drug development.
PAINS Filter
• PAINS = “Pan-Assay Interference Compounds”
• Problematic scaffolds – has cost their Institute time and $$
dsdht.wikispaces.com
Muchmore, SW et al. “Cheminformatic Tools for Medicinal Chemists” J. Med. Chem. (2010) 53, 4830 – 4841
Similarity Searching
What is it ??
Chemical, pharmacological or biological properties of two compounds
match.
The more the common features, the higher the similarity between two
molecules.
Chemical
The two structures on top are chemically similar to each other. This is reflected in their
common sub-graph, or scaffold: they share 14 atoms
Pharmacophore
The two structures above are less similar chemically (topologically) yet have the same
pharmacological activity, namely they both are Angiotensin-Converting Enzyme (ACE)
inhibitors
What is required for a similarity search ?
• A Database SQL or NoSQL ( Postgres, MySQL,
MongoDB) or flat file of descriptors eg: ChemFP
• Chemical Cartridge to generate fingerprints(descriptors)
for molecules ( RDKit, openbabel)
• Similarity function to calculate similarity( Jaccard, Dice,
Tversky) this can be written in c,c++ or python as a
function inside SQL databases.
2D fingerprints: molecules
represented as binary vectors
• Each bit in the bit string (binary vector) represents one molecular fragment. Typical length is ~1000 bits
• The bit string for a molecule records the presence (“1”) or absence (“0”) of each fragment in the molecule
− for a query substructure to be present in a database molecule each bit set to “1”
in the query must also be set to “1” in the database structure
• Similarity is based on determining the number of bits that are common to two structures
Calculate molecular similarity
Sequences/vectors of bits, or numeric values that can be compared by
distance functions, similarity metrics .
n B(x & y)
x yi
T (x, y)
E(x, y)
2
i B(x) B( y) B(x &y)
i1
3D based similarity
• Shape-based
ROCS (Rapid Overlay of Chemical Structures)
Silicos-it.com (Shape it)
− Flexible search
• In drug design, the term 'pharmacophore‘ refers to a set of features that is common to a
series of active molecules
• Hydrogen-bond donors and acceptors, positively and negatively charged groups, and
hydrophobic regions are typical features
It predicts ..
a) the pose of the molecule in the binding site
b) The binding affinity or score representing the strength of binding
Pose and Binding Site
• Binding Site (or “active site”)