QSAR AND DRUG DESIGN

Corwen Hansch

workable developed a methodology Quantitative Activity Relationship.

Structure In 1968,

Fraser Crum-Brown and

published an equation first general formulation of QSAR’S.

 The relationship b/w chemical structure and biological activity since before the turn of the century.

  On alkylation of the basic nitrogen atom different alkaloids produced the different biological effects

ammonium compouns quarternery

Biological activity is a function of the chemical structure

BA = f ( C) ….. (1)

 Richet discovered that the toxicity of organic compounds inversely proportional to water solubility.

So changing biological activity (ΔBA ) corresponds to the change in the chemical and physiochemical properties.

ΔBA= f (ΔC ) ------ (2)

  All the QSAR equation corresponds to eqn 2,because only difference in BA quantitatively correlates with changes in lipophilicity of the compound under investigation.

QSAR involves the derivation of mathematical formula which relates the biological activities of a group of compounds to their measurable physiochemical parameters.

 These parameters have major influence on drug’s activity.

Biological activity= function(parameters ) Activity is expressed as log(1/c),c is the minimum concentration required to cause a defined biological response.

QSAR based on Hammet’s constant uses eletronic properties as the descriptors of structures.

  The importance expressed as coefficient of lipophilicity octanol-water partition on biological activity.

This parameter provides a measure of site.

bioavailability to compounds which will determine the amt of compounds that reaches the target

  All these reveal that biological activity of a drug is a function of chemicalfeatures(lipophilicity,electro nic and steric) of the substituents and skeleton of the molecule.

Lipophilicity is the major factor governing transport,distribution and metabolism of the drug in biological system.

  Electronic and steric features influence the metabolism and pharmacodynamic process of the drug.

A major problem in QSAR studies arise because hydrophobic,electronic and steric effect overlaps and cannot be separated.

 

Parameters

Parameters assigned to various chemical grps that can be used to modify the structure of a drug.

The parameter is a measure of the potential contribution of chemical group to a particular property of the parent drug.

 Various parameters used are 1)

Lipophilic parameters Partitioncoefficient, Chromatographic parameters and π – Substitution constant.

 ii)Polarizability parameters : Molar refractivity(MR),Molarvolume, Parachor.

iii)Electronic parameters: Hammett Constant, Field and resonance parameters, parameters derived from spectroscopic data, Charge-transfer constant, Dipole chemical parameters.

moment, Quantum

 iv)Steric parameters: Taft’s steric constant,Vanderwaal’s radii.

v) Miscellaneous parameters : Molecular weight,Geometric parameters,Conformational entropies,Connectivity indices and other topological parameters .

   Lipophilic parameters Lipophilicity is defined as the partitioning of a compound between an aqueous and a nonqueous phase Partition coefficient (p) and lipophilic substituent(π).

 Partition Coefficient A drug passes through a no of biological membranes to reach its site of action.

The movement of drug through these membranes measured by P= C Organic / C aqueous……..(3)

 It is the ratio of concentration of the substance in organic and aqueous phase of a two compartment system under equillibrium conditions.

For easily ionisable drug P= C Organic / C aqueous(1-α ) Alpha is the degree of ionisation.

  The acuracy of correlation of drug activity wit coefficient will depend on the solvent system used as a model for the membrane.

Both pure and buffered soln used for the aqueous soln.

The n-octanol frequently polarity.

water system is chosen,because it appears to be a good mimic of lipid

 n-octanol more constistent results for drugs absorbed in GI tract less polar solvents such as olive oil give more consistent correlation of drugs crossing the blood brain barrier.

  n-octanol has pressure,allowing measurements.

a low vapour reproducible N-octanol is a UV transaparent over a large quantitative range detemination compound very easy.

making of the a

  More polar solvents such chloroform gives more consistent value for buccal absorption.

as It is a suitable model of the lipid constituents of biological membrane due to its long alkyl chain and the polar hydroxyl grp.

 The former parameter refers to the whole molecule while the latter related to substituted groups.

 Chromatographic Parameters When the solubility of a solute is considerably greater in 1 st phase than the other RP ,partition coefficient becomes difficult to determine experimentally.

Chromatographic parameters obtained from thin layer occassionaly used as substitutes for partition coefficient .

chromatography are

 Silica gel plate coated with hydrophobic phase is eluted with aqueous/organic solvent system of increasing water content.

 R f values are converted to R m which are the true measure of lipophilicity from the following equation values R m = Log(1/R f -1)

  Used as substitute for partition coefficient in QSAR investigations.

a)Compounds need not be pure.

b) Only traces of materials needed c) A wide range of hydrophilic and hydrophobic congeners can be investigated d) The measurement of practicaly insoluble analogs possess no problem e) No quantitative method for concentration determination needed f) The R m values.

values offer advantages compared to log P Several compounds simultaneously can be estimated

 The main disadvantages are a) Lack of precision and reproducibility b) Use of different organic solvent system renders the derivation of π and f related scales impossible.

Π Substituent constant/Lipophilic Substituent constant Lipophilic substituent constant are also known as hydrophobic constants (Π) The Π substituent constant defined by Hansch and Co-workers by the equation Π= log P RH - log P RX Where P RH and P RX are the partition coefficient of the std compound and its monosubstituted derivative respectively.

 Calculation of Log P value of m chlorotoluene is given as log P= log P = 2.13+0.71+0.56=3.40

log P= log P = 2.69+0.76 =3.45

log P= log P benzene toluene + π chlorobenzene + π Cl me- Cl + π + π me meta- Me = 2.84+0.51=3.35

logP exp = 3.28

 Polarizability parameters Molar Refractivity(MR) It is a measure of both the volume of a compound and how easily it is polarized It is expressed as MR = (n 2 -1)Mw (n 2 -1) d Where n is the refractive iondex M is the molecular weight d is the density

  The term Mw/d define a volume ,while the term (n 2 -1) / (n indicating how easily the 2 + 2) indicates a correction factor by substituent can be polarised.

This is particularly significant if the substituent has a π electron or lone pair of electrons.

Ligand –enzyme inter

 The significance of molar refractivity terms in QSAR equation of ligand enzyme interaction could be interpreted with the help of 3D structure bind in .The investigation shows that substituent modelled by MR polar areas while π,bind in hydrophobic space .

 The positive sign of MR in QSAR equation indicates that substituent binds to polar surface while a – ve sign or nonlinear relationship indicates steric hindrance at the binding sites.

 Electronic Parameters The distribution of electron in a drug molecule has considerable influence on the distribution and activity of the drug.

Polar and nonpolar drug unionized form are transported through the membrane than polar drugs in their ionized form.

  If the drug reaches the target site, the distributed electron will control the type of bond that it forms with the target site,which inturn affects its biological activity.

The Hammet Constant ( of the cpd.

σ ) He quantified the electronic effects of groups on the physiochemical properties The distribution of the elecrons within a molecule depends on the nature of the electron withdrawing and donating group found in the structure.

 Hammett used this concept to calculate Hammet constant(σ)for a variety of monosubstituted benzoic acids.

These constants were used to calculate the equillibrium constants for chemical reactions.

σ=log K Bx K B

 ie.,

σ

x

= log K

Bx

-log K

B pka= -log Ka σ x =pK B -pK Bx Where K B and K Bx constants for monosubstituted respectively.

are the equillibrium benzoic acids benzoic and acids

 Hammett substituent const is a measure of substituent and is the electron withdrawing or donating ability of a determined by comparing the dissociation of a series of substituted acid with that of the parent or unsubstituted acid

  A negative value of electron donor σ electron withdrawing x indicates that the substituent is acting as an and the +ve value indicates that it is acting as a group.

Hammett constant take into acount both resonance and inductive effect

   The values for σ for a particular substituent varies depending on their substitution meta or para .

The meta and para sigma value are commonly used and indicated by a subscript m or p after the symbol , The o – substituent due to steric hindrance and other effects such as are unreliable intramolecular hydrogen bonding.

  Resonance or Mesomeric effect Inductive effect.

 3 Kinds of static or electical influence states predominates Electrical effect which is transmitted by the polarization of electrons from atom to another.

Hammett Constant

Inductive effect --- Intrinsic tendency of a substituent to withdraw or release electns.

For eg G releases e-ns ,disperses charge and stabilises cation.

G C+

 G withdraws e-ns intensifies the charge and destabilizes the cation.

G C +

Disadvantage: They only apply to substituent directly attached to benzene ring.

If there is more than one substituent the value are summarized as Σσ.

This constant biological distribution involved.

was is unsuccessful activity not sine the to the only relate electron factor

 Inductive substituent Constant : Hammett constant are a measure of both inductive and mesomeric effect.

The para substituent constant (σ equivalent meta constant (σ m inductive substituent constant (σ i p ) ) has a greater resonance component than the ) and the inductive const can be calculated from the σ i = ½ (3 σ p - σ m )

 It is used in aliphatic compound where influencing and influenced form group do not form a part of the conjugated system.

Taft’s substituent constant( σ * ) It is a measure of the polar effects of a substituent in aliphatic compound.

They are based on the hydrolysis of ester and calculated frm the following eqn

σ * = 1 log(k/k o ) B 2.48

ie, 0.403 log(k/k

o ) B

- log(k/k o ) - log(k/k o ) A A Where k is the rate constants for the hydrolysis of substituted compound.

K o those of methyl derivative .

B represents basic hydrolysis and A acid hydrolysis respectively.

 The factor 1/2.48 is a constant B is influenced by the polar effect so that by substracting the acid term from the basic term only the polar effect remain.

In taft’s substituent constant only the methyl group is the standard for which the constant is zero.

 That can be compared with other constant by writing the methyl group in the form CH 2 -H and identify it as the group for H.

Taft and Inductive constant are related as substituent σ * = 2.51 σ i

 Steric Substitution Constant The drug approaches the binding site to interact with the receptor.The bulk,size and shape of the drug may influence on this process.

A steric substitution constant is a measure of the bulkiness of the grp it represents and its effect on the closeness of the group between the drug and receptor site.

CHARTON’S STERIC CONSTANT

 Corrected Vanderwaal’s radius U in which the minimum Vanderwaal’s radii of the substituent grp radius for hydrogen(r vH ), (r v(min) is corrected for the corresponding

U= r v(min) -r vH= r v(min) -1.20

 Other parameters Molecular wt was used by Lien to improve the fit of parabolic Hansch equation .A more appropriate use of MW was demonstrated in the QSAR study of multidrug resistance of tumour cells,where MW stands for the dependence of biological activities on diffusion rate constant.

  The relationship bw MW and volume implies that 3√MW correponding to linear dimension of size should be better than log MW.

Indicator dummy variables variables or known as de-nova constant are used in linear multiple regression analysis .It is used to account for other structural features

1.Intramolecular

bonding,hydrogen donor & acceptor properties,ortho models etc.

hydrogen effects,cis/trans isomers,diff parent skeleton,diff test

 Quantitative Models This is the most popular mathematical approach to QSAR introduced by Hansch.

a) b) Transport of drug to its site of action.

The binding of drug to its target site .

Each of these stages depends on the chemical and physical properties of the drug and its target site.

In Hansch analysis these properties are described by the parameters which correlates biological activity.

 The most physiochemical steric parameters commonly parameters Hansch analysis are logp,Π, σ and are practically all the parameters used in Hansch analysis are linear free energy related ,so it is known as linear free energy approach or extra thermodynamical approach.

used for

 The relnship b/w P and drug activity depends on P values.

If range is low,straight line is obtained .

log1 /c = K 1 log p +K 2 Where K1 and K2 are consts.

 If the hydrophobicity values are limited to a small values then the eqn will be linear.

log1 /c = k 1 log p +k 2 σ+K 3 E 3 + K 4 Where K1,K2,K3 are constant obtained by least square procedure.

 C is the molar concentration that produces biological action.

The molecule which are too hydrophilic or lipophilic will not be able to cross the lipophilic or hydrophilic barriers respectively given as .

If the P value are higher range ,then the equation will be parabolic and



log1/c =-K1 (logp)

2

+K

2

logP+ K

3

σ+K

4

ε

S

+K

5 The constants K 1 -K 5 are obtained by least square method .

Here the graph of log 1/c against log P Have max values logP 0 ---- parabolic form

 Below the

P

0

– drug reluctant to enter the membrane.

 

Above

P 0-- membrane.

reluctant to leave the The analogues with partition coeff near this optimum value ---- most active and important to study.

HANCH RULES

 Selection of independent variables Different parameters 1.

log P,Π, σ,MR, steric 2.Parameters should be independent.

The intercorrelation coefficient should not be larger than 0.6-0.7

.

 All the reasonable parameters must be validated cross validation.

by appropriate statistical procedure ie either by a step wise regression analysis or

 The best equation is normally one with lower standard deviation and higher F value .

   If all the equations are equal simplest one is taken Number of variable should be at least 5 or 6 data point per variable to avoid chance correlations .

Model should be consistent with the known physical.organic and bio-medicinal chemistry of the process under consider consideration .

 APPLICATIONS OF HANCH ANALYSIS It is used to predict the activity of an yet unsynthesised analogue.

This chemist enables to the medicinal synthesise analogue which is worthy.

an

 Eg; The adrenergic activity of a series of analogues of β haloarylamine was observed.

It was found that only values not the steric factor.

pi and sigma only related to activity and Log 1/c =1.78pi+1.674

 Hansch analysis may also used

to give an indication of the importance of the influence of parameters on the mechanism by which a drug acts.

 The smaller the value of coefficient of σ relative to that of Π in tha above eqn shows that electronic effect do not play an important role in the action of the drug.

The accuracy of Hansch equation depends on 1) The no of analogues used

2) The greater the number, the higher the probability of obtaining an accurate Hansch equation.

3)The accuracy of biological data used in the derivation of the equation.

4)The choice of parameters.

 Free Wilson Analysis It is a true SAR model.This is an alternative approach to Hamett model.

Free and Wilson used a series of substituent constant which related biological activity to the presence of a specific functional group at a specific location on the parent molecule.

 The relationship between biological activity and the presence or absence of a substituent is expressed by the following equation Activity=A+G ij X ij Where A is average biological activity of the series,G ij is the contribution

of actiivity to the functional grp i in the jth position and X jth position.

ij the presence or absence of functional grp i in the Whole molecular descriptors Some of the natural

 The method is based upon an additive mathematical model in which a particular substituent in a particular molecular position is assumed to make an additive and constant contribution to biological activity of a molecule in a series of chemically related molecule.

Mathematical model or denovo approach are the synonyms model,Additive

for free wilsonmethod.

The substituent constant based on activities are used rather than physical properties

 Introduction substituent at a particular molecular position of always a lead particular to a quantitatively as similar effect on biological potency of the whole molecules and expressed by the eqn BA=µ+Σa i j

 Where Σa i occurs position µ= contribution of unsubstituted compounds j = contribution of substituted compounds i = number of position the substituent j = number of substituent at that

APPLICATION 1)Easy to apply 2)Simple a first method look on to derive the substituent contribution and to have their possible dependence on physiochemical properties different

  The substituent which cannot fulfill the principle of additivity can be recognized.

Substituent constants like were not considered and so this method is effective σ and Π when substituent available.

constants are not

 Disadvantages The structural variation is necessary in at least two different position of the substitution.

Large no of parameter is necessary to describe relatively few compounds.

Only a small no of new nalogues can be predicted.