An Importance and Advancement of QSAR Parameters in Modern Drug Design: A Review
An Importance and Advancement of QSAR Parameters in Modern Drug Design: A Review
ABSTRACT
From the last 30 years QSAR study become important to help in designing of drugs. Many software and hardware are used in
discovery of drugs. QSAR techniques employ powerful computers to generate large data bases result. Quantities structure
activity relationship and quantative structure property relationship is essentially computerised statistical methods which help to
explain the observed variance in the biological effect of certain classes of compound. It assumes that a biological activity of
certain compound is a function of various physico chemical properties are favourable to the concern activity. QSAR involve
themathematical and statistical analysis of SAR data which help to minimize the molecular modification. The ultimate
objective of the present study is to understand the force governing the activity of a particular compound. QSAR is thus a
scientific achievement and an economic necessity to reduce empiricism in drug designing. The purpose of this study is to
describe the techniques like ANN, multiple linearanalysis, partial least square method. The study further extent the types of
QSAR like 1D , 2D, 3D,4D,5D etc used in designing of molecule with the help of different software. QSAR is thus a scientific
achievement and an economic necessity to reduce an empiricism in drug design to ensure that every drug synthesized and
pharmacologically tested should be as meaningful as possible. The uses of different parameter are important to calculate
biological activity of molecule.
Key words: QSAR, Drug design, QSAR PARAMETER, MLR, ANN, Free energy model, Quantum mechanical
method.
INTRODUCTION
The introduction of Hansch method in 1964 enables chemist to describe SAR study in quantative terms.1during past
decade, QSAR started to develop from a more advance stage. Various methods used in QSAR analysis can be explain
that how the molecule is biological active or not. The activity of the molecules depends upon the different parameters
like hydration energy , molecular volume of molecule, parachore, Log P values etc.2 Physical organic chemistry deals
with characterisation of the structure and prediction of the properties, the descriptor for which are usually found
experimentally. If some property depends on the set of selective descriptor,3 the ordering of the structure will parallel
the ordering of properties. In other term the structure information is coded in these properties. The good correlation of
physico chemical properties with a particular set of indices may help in understanding the contribution of these
invariants in determining the property.4
QSAR PARAMETERS
Mainly three types of physico chemical parameters are used in qsar study. Hydrophobic parameters, electronic
paramerter, steric paramerter.5Hydrophobic parameter like partition coefficients log P , Pi substituent’s constant , Rm
chromatographic parameters log Rm, solubility parameter δ, elution time in HPLC log K and parachore [P] are used
in study.6-7
Electronicparameters can be divided in two parts 1. Experimental parameters 2.Theoretical parameters. Experimental
parameters are generally Ionization constant Pka and sigma substituent’s constants σ and σ2 and spectrochemical
shifts are described in term of ppm. Resonance effect R and field effect F and Ionization potential calculated in term of
I.8-9
The theoretical quantum mechanical indices10 MO indices are generally studied as atomic charge densities ϵ, atomic
net charge q, QT, Q, σ etc. and superdelocalisibilitySrn and energy of molecular orbital’sby ELEMO andEHOMO and other
parameters like N, NH, F[A] are also used to describe and calculate the data for result . The third type of parameters are
steric parameters like Taft's steric substituent’s constants ES andvandersWaals radii ϒ and inter atomic distances B, L ,
molar refractivity MR and molar volume MV are used for calculation of value to obtained biological
activity.11Biological activity reflects the fundamental physico chemical properties of the bioactive compound. The
lipophilicity,polarity,chargesdistributions are main parameter to calculate biological activity.
Major problems in QSAR study arises because hydrophobic, electronic and steric effect overlap and cannot be separated
easily. The parameters which are used to obtain such relation can be divided in to two parts. Those which describe
mainly the physical properties of Skelton such as water solubility partition coefficient, chromatographic Rf value,
molecular weight, surface tension etc. And another parameters which describe the chemical properties like charge
densities, dipole moments, electron donor acceptor properties, Hammetts constant, Taft's steric constant.12-13
QSAR METHODS
There are various methods to calculate the qsar parameters. The introduction the hansch methods.In 1964 enables
chemist to describe the SAR studied in quantative term. There are various methods like free energy model and other
statistical methods, Pattern recognition, topological methods14,quantum mechanical methods and molecular
modelling15 are used for qsar analysis. The free energy model is used to describe the Hansch method for linear free
energy relationship calculation and other method is free wilson mathematical model.
The statistical methods are discriminate analysis, principal component analysis, Factor analysis, cluster analysis and
combined multivariate analysis. There are two way to get a quantative information about SAR, one may use QSAR
methods based on linear free energy relationship which relate the biological activity of a molecule with contributions
from various free energy related physicochemical parameter of the substituent’s. In other approach mathematical model
are used other than linear free energy relationship.
SUBSTITUENT CONSTANT
When any dose of drug is given to patients than the molecules change its confirmation from Site of administration to
site of action. Each confirmation has its differentenergy. The pharmacokinetics and pharmacodynamics property are
based on Lipophilic (Log P)16, electronic sigma and steric featureEs of the drug molecule. The biological activity BA
can be calculated by the equation
there were no optimum lipophilicity in each series, compound with infinite biological activity would result if only their
lipophilicity were high enough. In series of compounds where biological activity is dependent mainly upon
lipophilicity, one can not go on increasing the biological activity indefinitely by increasing lipophilicity of the
compound. Activity rises to a maximum (log Po) and then decline.
CHROMATOGRAPHIC PARAMETERS
Other parameter Rm is used to describe the lipophilicity of drug molecule. In 1965 Boyce and milborrow suggested the
use of Rm value19 from reverse phase thin layer chromatography (TLC)
as alternate lipophilicity parameter in QSAR.
Rm = log (1/Rf -1)
However Rm values cannot be regarded as true equilibrium parameters. Usually silica gel plates are impregnated with
liquid paraffin’s , silicone oil, ethyl oleate or n -octanol as stationary phase. While mobile phase may consist of mixture
of polar solvents like methonol, ethynolaceton with water or aqueous buffer solution. There are many advantages of
using Rm value instead of partition coefficient.
Molarrefractivity is a volume term, but is also proportional to electron polarisibility25. Hence its calculation in qsar is
difficult. If the confirmation dependent steric effect are eliminated molar refractivity may be used as a measure of
dispersion and dipole induced dipole forces. Molar refractivity and parachore may be regarded as corrected molar
volumes. The lipophilicity however not be correlated with molar volume, molar refractivity and parachore if polar
compound are induced.
The capacity factor (log k) determined by HPLC is an indicator of lipophilicity.
K =tr - to
to
wheretr and to are retention time of solute and unretened compound for a series eluted by
methanol–water mobile phase. Centrifugal counter cartography is used to measure partition coefficient.
ELECTRONICS PARAMETER
The polar charecters of drug are calculated by electronic paramete26. σm is called hammett constant for meta
substituents derived from ionisation of benzoic acid.27σpis known is hammett constant for para substituents derived
from ionisation of benzoic acid. σp- Hammet constant used when there is direct conjugation between substituents and
reaction centre derived from aniline and phenols. σ P+is calledH.C. brown constant derived from solvolysis of di
methyl phenyl carbonyl chlorides. σ1 constant describing solely polar effect. σR constant describing solely mesmeric
effect28.σ* is Taft's polar substituents constants29 derived from hydrolysis of aliphatic esters. Sigma parameter is
haemolytic constant for substituents interacting with a free radical reaction. F and R are field and resonance
components derived for linear combination of σm and σp values. The main commonly used electronic parameter is
Hammett substituents constantswhich can be obtained from the dissociation constant KX and KH of the benzoic acid.
The substituents constant σis depends on Delta G; free energy arises due to dissociation of benzoic acid. Hammetts told
that and electron withdrawing group attached to aromatic ring of benzoic acid would enhance the acid strength of
carboxylic group.
MR= (n2-1) MW
(n2+2) d
Where n= index of refration at the sodium D line
mw= molecular wt of compound.
d= density.
Larger the value of MR greeter the steric contribution.
If the substance is in liquid form the than molar refractivity value can be calculated by using the Lorentz-Lorentz
equation
MR= MW(n2-1)
d (n2+2) ( cm3/mol )
Molecular connectivity indexes are the third type of steric parameter. The degree of branching in given compound is
described by molerconncetivity indexes. Connectivity indexes can be calculated by writing Skelton formula without the
hydrogen atom.
Insert Formula x X=E (deltai) etc.
Effect of steric and electronic parameters on lipophilicity: - Due to inductive effect overall lipophilicity32is affected and
electron withdrawing groups increases Pie value when a hydrogen bonding group is involve. The nitro group and
hydroxyl group present at aromatic Skelton , the electron withdrawing inductive effect of the phenyl ring and the nitro
group make the non-bonded electrons on the hydroxyl group less available for hydrogen bonding. Theσ and Es is a
function of QSAR activity. Hence compound having substituents with almost the same value of these parameters may
be considered as isomeric bioisosters.
Hansch model and free wilson model are included in liener free energy method.
According to hansch model there are two types of dependent liener and nonliner .
Hansch suggested two processfirstly the journey of drug from the point of entry in the body and secondly drug receptor
interaction in the body. Biological activity will be as follows
MOLECULAR MODELLING
It is a technique which describe the actual shape, confirmation, size of the molecules which intract with receptors. It is
the branch of science that help in designing of drugs with the help of computers and software’s. The three
dimensionalstructural confirmation analysis to correlate physiochemical paramerts with the biological activity. The
softwares used for molecular modeling35. aregenrally molecular modeling pro, amber, Hydra, SYBL, FRODO etc.
Thesesoftware’s are available for calculation of molecular properties for search of new lead compounds.
Toplissdecision tree method help to select a limited no. of substituents which will give good relation between π, σ,and
Es. If the biological activity will increases it can be attributed to positive value of π, σ.
REGRATION ANALYSIS
The multiple linearregreation method (MLR) 36 is used to determine the biological active molecule from the large data
pool. MLR is the method used to predict the relationship between two of more explanatory variables by creating a
linear equation. Inqsar correlation coefficient r is calculated with the help of graph. Ploting a graph is prefer method to
obtained biological activity. The term correlation coefficientr , number of compound utilised n, standard deviation S,
and statical validity F is calculated in this method. A regression equation is derived to obtain high value of
regressioncoefficient r which indicate the significant of regreation analysis while low value of regreation analysis r
indicate the substituent constant is not important. For a better correlation large no of compound must be used. The
value of r must be related to n number of compound for example if r = 0.87 for n=10 is better correlation with r = 0.98
and n=3.
The value of statndard deviation S is large means larger is the accuracy with which the expected activity of new
compound can be judge.
The term regreation parameter r 2value help us to understand whether other parameters should be remain counted or
uncounted. Greater the value of r2 lesser the value of variance data that remains unaccounted by the equation.
The term F evaluate in qsar for the statical validity of a particular equation.
PURPOSE OF QSAR
QSAR techniques gives hints about the activity of molecules so it save the time and coast of synthesis of chemical
molecule. It saves time and effort in clinical trials. It provides the information about the biological activity of the
molecules from the large data of molecules so it save time and coast of synthetic chemist.
APPLICATION OF QSAR
QSAR help to know the biological active compounds in a series of molecules. It also gives idea about the toxic
substances. It also for cast the biological activity of molecules. It also provides idea about pharmacological response of
the molecule in better way. It can be also used to know the surface active agents, purfume, dye and fine chemicals.
QSAR provides idea about selection of proper substituents. It is also provide idea about drug receptor interaction. It also
provide idea about correlation between various type of parameters and pharmacokinetics feature of drug.
QSAR STEPS
Structural Entry and molecular modelling
Descriptor Generation
Feature selection
Construct model
MLRA or CNN
Model validation
CONCLUSION
QSAR is basically used to determine the biological activity of drug but now a day’s 3D qsar is very important to
correlate the structure feature with the biological activity and its toxic effect if any.
Various newly parameters are used and values calculated with the help of different qsarsoftware’are important to
calculate biological activity. Now a day’s 4D qsar, 5D qsar, 6D qsar are used but 3D qsar is the base of all model for
researcher. Various parameters like electronic, steric etc are used to search the most potent drugs for disease. The role
of use of physicochemical parameters is not limited because hydrophobic parameters, electronic parameters, and steric
parameters are very important and base of to calculate newly parameters.
REFERENCES
[1.] Hansch Crown ,fujita Toshio application of Hanschs model to Guaeanolide derivatives, Journal of american
chemical society 1964, 86,8 1616-1626.
[2.] doi.org/10.1021/jacs.0c07761
[3.] Robert C. Putnam ,Parachore and other properties of organotins, Canadian journal of chemistry, 1986,44,12 1343-
1350.
[4.] https://doi.org/10.1139/v66-202
[5.] Artemcherkasov, Eugene N. Muratov, denis F. Qsarmodeling where have you been where are you going to? J. of
med . Chem. 2014, 57(12) 4977-5010.
[6.] doi: 10.1021/jm4004285
[7.] sylviacorolina, alcazaralay, mariaangela, physicochemical properties modification and application of starch. Food
science technology 2015 , 35 (2) 215-236.
[8.] https://doi.org/10.1590/1678-457X.6749
[9.] kapoor Y , kumar K. Qsar in drug design science forcast 2019 2 (2). 1017.
[10.] ToshirFujita ,JunkichiIwasa and crown hansch . J. American chemical society 1964 86 (23) 5175-5180
[11.] doi: 10.1584/jpestics.M17-01
[12.] Xinliang YU, Xveyewang, Hangluwang , prediction of solubility parameters QSAR and combinatorial science
2005 (2) 156-161.
[13.] https://doi.org/10.1007/978-94-007-6064-6_4
[14.] O.G. mekenyan and G.D. veith . The electronic factor in qsar SAR and QSAR in environmental research 2006 2
(1-2) 129-143.
[15.] https://doi.org/10.1080/10629369408028844
[16.] CrownHansch, A leo and R.W. Taft A survey of Hammetts constantchem Rev.. 1991- 91(2) 165-195.
[17.] DOI: 10.1021/cr00002a004
[18.] R. CarboDorca L Amat. E Besalu. X Girones. D robert. Quantum mechanical origin of qsarJ.of molecular
structure ( Theochem) 504 (2000)181-228.
[19.] https://doi.org/10.1016/j.jmgm
[20.] ZinehAlmi, Salah Belaidi, TouhamiLanez. International letter of chemistry , physics and astronomy. 2014,
37,113-124.
[21.] https://doi.org/10.18052/www.scipress.com/ILCPA.37.113
[22.] Mehroldshapar, sharminEsmaipoor . Asian journal of green chemistry, 2017(2) 116-129.
[23.] DOI.10.22631/AJGC.2017.94413.1010
[24.] Depankardutta , Taftssubstitutents constants, journal of physical organic chemistry 1991 (4) 96-100.
[25.] https://doi.org/10.1021/jo01345a067
[26.] VladmirA ,polyulin, Eugene V, Radchenkonikolai, S. Zefirov, journal of chemical information and computer
science. 2000, 40 (3) 659-667.
[27.] https://doi.org/10.1021/ci9604605
[28.] Maral aminpour ,carlomontemagno, jack A. Molecule 2019,24 1693 (1-30).
[29.] https://doi.org/10.3390/molecules24091693
[30.] Grossy G. Chavanieu A. Molecular lipophilicity : A predominentdiscripter for qsar. Chemogenomics and chemical
genetics. 2011, 153-170.
[31.] https://doi.org/10.1021/ci025516b
[32.] CatlenC ,bannan , Calabro G. Dausy Y. Kyer. David L. Modeley. Calculating partition coffecient of small
molecules in octanol water and cyclohexane water. J. Chem. Theory and computation. 2016. (12) 4015-4024.
[33.] https://doi.org/10.1021/acs.jctc.6b00449
[34.] MeladiyahI ,Jumina J . mubarika S. Mustfa M. Biological activity quantativestaructure activity relationship and
molecular docking of xanthone derivatives as anticancer drugsJ. Med. Chem. 2001, 44, 5, 672–680.
[35.] https://doi.org/10.1021/jm000955s
[36.] Apostolov S Vastag G .Nakomcic J. StojiljkovicI ,charamatographic descriptor in QSAR study of barbeturates .
journal of liquid chromatography and related technology. 2019 Vol. (42) issue 7 194-203.
[37.] https://doi.org/10.1002/jssc.201900292
[38.] GozalbesR .Pinada A QSAR based solubility model for drug like compound. Bioorganic and medicinal chemistry.
2010. 18 (19) , 7078-7084.
[39.] DOI: 10.1016/j.bmc.2010.08.003
[40.] Anderson C. Yuji T calculation of molecular surface area with numerical factors. Computers and chemistry. 1992,
16 (4) 227-284.
[41.] https://doi.org/10.1021/jm980313t
[42.] shukla R.K. , kumaratul , shrivastavaurvashi. Density refractive index and molar refrectivity of binary liquid.
Arabian journal of chemistry. 2016. 9 (2) 1357-1367.
[43.] DOI:10.1016/j.arabjc.2012.02.013
[44.] V.S. kuzmin, S.B. katser calculation of vanderwaals volume of organic molecules. Bullatin of russian academy of
science. 1992, 41, 720-727.
[45.] DOI: 10.4018/978-1-4666-3604-0.ch059
[46.] Controversies over research embryonic stem cell S Rajpoot, D Paul, J Kaur - International Journal of Chemical
and Pharmaceutical 2019, 6 (2) 1-8.
[47.] Exner, conception and significance of the parachore. Nature 1962. 196., 890-891.
[48.] DOI: 10.1038/196890a0
[49.] Tian T, Declan S. Dale H, Lu Hua Li, electrocnicpolarisibility as the fundamental variable. Nano letter 2020. 20
(2) 841-851DOI: 10.1021/acs.nanolett.9b02982
[50.] Sahdev A.K., B. Sethi, sumaL.Rawat, Amrita singh, neetianand, A review article role of qsar significance and
uses in drug design. 2018 , 5 (1) 426-435
[51.] Gerrtschurmann , do hammett constant model electronic properties in qsar , scince of total environment1991. 110,
221-235.
[52.] Yagyesh K and Kapil K Drug Design: An Overview, Medicinal & Analytical Chemistry International Journal
2019, 3 (2) 000136.
[53.] DOI: 10.23880/macij-16000136
[54.] Future trends in fiber optics communicationS Rajpoot, P Singh, S Solanki, SJ YasinInternational Journal on
Cybernetics & Informatics 2017, 6 (1), 23-28
[55.] DOI: 10.5121/ijci.2017.6203
[56.] Review Article Review on Stem Cells: Basics Classification and Applications*S Rajpoot, G Tewari - Int. J. Pharm.
Sci. Rev, 2018. (12), 48-52
[57.] Kari N, Helmi N, Koch A, Erich K Tafts equations in the light of NBD, computational and theoretical chemistry,
2012, 981, 52-58.
[58.] https://doi.org/10.1021/ja00839a047
[59.] shankar P. Mitra , drug receptor interaction pharmacological binding and therodynamics A review. Journal of
surface science and technology 2009, 25 (3) 103-152.
[60.] Lemont B. Kier, the meaning of molecular connectivity: a biomolecular accessibility model CroaticachemicaActa,
2002, 75 (2) . 371-382.
[61.] https://hrcak.srce.hr/127526
[62.] TizianaGinex, javier V. Enrich G. Lipophilicity in drug design an overview. Future medicinal chemistry. 2019 11
(10), 111-116.DOI: 10.4155/fmc-2018-0435
[63.] Catlin C, Gaetano C, Daisy Y, calculating partition cofficent of small molecule in octanol water.
[64.] J. Chem. Theory Comput. 2016, 12, 8, 4015–4024 doi.org/10.1021/acs.jctc.6b00449
[65.] J.chem. Theory. 2016. 12 (8) 4015-4024.doi: 10.1021/acs.jctc.6b00449
[66.] Hansch Crown, A method for the correlation of biological activity and chemical structure. JACS.1964, 86, 1616.
[67.] https://doi.org/10.1248/cpb.40.127