The Application of Some Computing Techniques in the Drug Design
Nazmul Islam1* and Savas Kaya2
1Department of Basic Science & Humanities/ Chemistry, Govt. Engineering College Ramgarh, India
2Department of Chemistry, Cumhuriyet University, Turkey
Submission: April 26, 2017; Published: June 02, 2017
*Corresponding author: Nazmul Islam, Department of Basic Science & Humanities/ Chemistry, Govt. Engineering College Ramgarh, Jharkhand, India, Tel: + 91.9432878737 Email: nazmul.islam786@gmail.com
How to cite this article: Nazmul I & Savas K. The Application of Some Computing Techniques in the Drug Design.Organic & Medicinal Chem IJ. 2017; 2(5): 555596. DOI: 10.19080/omcij.2016.02.555596
Abstract
Using global and local reactivity descriptors of DFT and MO theory, in this work, we obtained a nice correlation between theoretically computed and experimentally determined drug activities (MIC). This study once again reveals the fundamental statement of theoretical chemistry that "the structure of a molecule determines its electronic wave function which in turn determines many of its physicochemical properties”. It is reasonable to suppose that the biological activity of molecules its potency as a drug for examples will also, often enough, be dependent upon the molecular electronic wave function. It is possible in principle to determine any molecular property that is ultimately dependent upon the molecule's electronic ground state wave function. In practice this can sometimes be achieved using rigorous procedures.
Abbreviations: QSAR: Quantitative Structure Activity Relationship; QC: Quantum-Chemical; EHT: Extended Huckel Theory; CNDO: Complete Neglect of Differential Overlap; INDO: Intermediate Neglect of Differential Overlap; MINDO: Modified Intermediate Neglect of Differential Overlap; MNDO: Modified Neglect of Diatomic Overlap; AM1: Austin Model 1; PM3: Parametric Model 3; RHF: Restricted Hartree-Fock level.
Introduction
Years of cumulative research can result in the development of a clinically useful drug, providing either cures from a particular disease or symptomatic relief form a physiological disorder. A dug, with a diseased pharmacological activity may have associated with some undesirable side effects which limits its' use. Current research reveals that the function of a drug may associate with its structural features. Computational theoretical chemistry is a branch of chemistry that uses theoretical studies to assist in solving chemical problems. It uses the results of theoretical chemistry, incorporated into efficient computer programs, to calculate the structures and properties of molecules and solids. While its results normally complement the information obtained by chemical experiments, it can in some cases predict hitherto unobserved chemical phenomena. It is widely used in the design of new drugs and materials.
Now a day, the synthesis of novel pharmacologically active molecules with reduced toxicity is of prime interest. Recently, QSAR has gained importance in the field of pharmacological sciences [1]. Quantitative structure Activity relationships (QSAR) are predictive tools for a preliminary evaluation of the activity of chemical compounds by using computer-aided models. The computational chemist can assist the experimental chemist or it can challenge the experimental chemist to find entirely new chemical objects.The Hohenberg and Khon theorm based DFT [2-4] provide a major boost to the computational chemistry The performance of DFT method in description of structural, energetic and magnetic molecular properties has been reviewed quite substantially in recent time.
DFT methods are in general capable of generating a variety of isolated molecular properties [5-12]. The commercial exploitation of organic compounds as a medicinal drug for examples is likely to require, at some stage of its development, the determination of biological or chemical activities or properties related to the intended end use of the compound. It is desirable therefore to have at hand relatively straightforward and inexpensive procedures enabling the efficient and accurate prediction of a molecular activity or property especially when its direct measurement by experiment is, for one reason or another, to be avoided if at all possible. The procedures that are conventionally used for indirect determinations of activities make use of molecular “descriptors” which include suitable molecular properties and physical-organic constructs obtained from both experimental and computational sources. Molecular descriptors are ultimately related to molecular structure. Hence the relationships between activities and the descriptors on which they depend are generally known as Quantitative Structure- Activity Relationships (QSARs) or Quantitative Structure- Property Relationships (QSPRs) depending upon whether the property of interest is to be characterized as biological or nonbiological, respectively.
The structure of a molecule determines its electronic wave function which in turn determines many of its physicochemical properties. It is reasonable to suppose that the biological activity of molecules its potency as a drug for examples will also, often enough, be dependent upon the molecular electronic wave function. It is now possible to perform accurate ab initio calculations routinely on molecules of moderate size at reasonable cost. Thus, it would clearly be advantageous to have available a procedure for estimating a molecular activity or property from descriptors which are derived from, and in turn characterize, the electronic wave function of a molecule.
It is possible in principle to determine any molecular property that is ultimately dependent upon the molecule's electronic ground state wave function. In practice this can sometimes be achieved using rigorous procedures. It can, for instance, be achieved for certain properties using quantum-mechanical expectation values. But when rigorous procedures are difficult or virtually impossible then it is necessary to use either a less rigorous model calculation or to use the Quantitative structure- activity relationship or Quantitative structure-property relationship methods.
Quantitative structure-activity relationship (QSAR) techniques increase the probability of success and reduce time and cost involvement in drug discovery process [13,14]. In this work, a correlation between the DFT descriptors and the drug activity of the Benzo[f]quinazolin-1(2H)-ones has been discussed. The aim of the present work is to explore the efficacy of the DFT based descriptors to provide some relationship by which the biological activity of hypothetical molecule can be measure prior to their synthesis. This technique shall reduce the drug discovery coast, time and efforts. In recent year heterocyclic compounds analogues of benzothaizoles and derivatives have attracted strong interest due to their useful biological and pharmacological properties. The chemistry of quinazolinones is interested because of its biological significance. Many of them shows antifungal [15], antibacterial [16], anticancer[17] anti- inflammatory[18] anticonvulsant [19] and antiproliferative activities as well as inhibitory effects for thymidylate synthase and poly-(ADP-ribose) polymerase (PARP) [20]. Quinazolin- 1-one derivatives are also found to be tranquilizer, antiallergic agent, an antiulcer agent and antiasthmatic agent.
Quantum-chemical (QC) calculations are a key element in biological research [21]. When constantly tested for their range of validity QC methods provide a description of how molecules interact and form their three-dimensional shape, which in turn determines molecular function. They can aid the formulation of hypotheses that provide the connecting link between experimentally determined structures and biological function. QC calculations can be used to understand enzyme mechanisms, hydrogen bonding, polarization effects, spectra, ligand binding and other fundamental processes both in normal and aberrant biological contexts. The power of parallel computing and progress in computer algorithms are enlarging the domain of QC applications to ever more realistic models of biological macromolecules.
The key insight of chemistry is the relationship between molecular structure and molecular function. We use the details of molecular structure to predict the properties of molecules. Medicinal chemistry is a particularly glaring example of our use of structure-function relationships. There is a tremendous need to be able to quickly design new drugs for curing human disease. The rapid prediction of the activities of compounds for use as drugs and the discovery of new compounds is an important goal. Quantitative Structure Activity Relationships, or QSAR [22] enables us to predict the properties of compounds and are a quantitative expression of structure-function relationships. QSAR has been responsible for the rapid development of many new drugs.
An example of a QSAR study is the narcotic activity of esters, alcohols, ketones, and ethers with tadpole. In this study various organic compounds were added to water with swimming tadpoles. The swimming speed of the tadpoles was observed and the amount of the compound that was necessary to slow the tadpoles swimming was determined. A very effective compound has a very low concentration for the production of the desired effect. In QSAR studies we often like to have the more effective compounds have a higher "activity,” not a lower. Therefore, it is very common to transform the concentration for a desired effect, C, to an activity by:
A = log(1/C) (1)
The log(1/C) value increases with compound efficacy
The genesis of QSAR is from physical organic chemistry and linear free energy relationships. The first such studies were done by L. P. Hammett [23]. His goal was to uncover the effects of electronic structure on organic reactivity. A short discussion of his work will be instructive as we start to understand the foundations of QSAR. Hammett's first studies were to understand the effect of electron withdrawing and donating groups on the pKa's of substituted benzoic acids, Hammett first wanted to develop a descriptor that described inductive substituent effects. He compared the log Ka for a variety of substituted benzoic acids with the log KaH for unsubstitued benzoic acid to define the substituent constant σ.
σ = log Ka - log KaH (2)
σ the Hammett constants for meta and para substituents. He then postulated that other properties, other than acidity, would be likewise effected by the same substituent effects, and that these other properties would follow the relationship
log (property) = ρ σ + cst (3)
In other words, the expected function was a linear relationship with the ρ inductive constant.
Hammett obtained the equation from the relationship between Gibbs free energy and the equilibrium constant for a reaction,
Δ rG = - RT ln Keq (4)
Semi-empirical methods can be used to calculate quantum chemical descriptors [24]. A number of semi empirical methods have been developed over the last several decades. To name but some of the most popular extended Huckel theory (EHT) [25], complete neglect of differential overlap (CNDO) [26] intermediate neglect of differential overlap (lNDO),[27]modified INDO- (MINDO), [28]modified neglect of diatomic overlap (MNDO),[29] Austin model 1 (AM1), ),[30] and parametric model 3 (PM3) [31].
The density functional theoretical parameters
Given the electron density function p(r) in a chemical system (atom or molecule) and the energy functional E(p), the chemical potential, μ of that system having N number of electrons in equilibrium has been defined as the derivative of the energy with respect to the number of electrons at fixed molecular geometry
The chemical potential [32], μ, is given by
μ= [δE(p) / δp] (5)
The differential definition more appropriate to atomic system is
μ = [∂E / ∂N]v (6)
Then following Iczkowski and Margrave [33], Parr and Pearson [34] defined the electro negativity as the additive inverse of the chemical potential-
χ = - μ (7)
or, χ = - [∂E/∂N] v (8)
and the absolute hardness, η, as
η = 1/2[∂μ/∂N] v = 1/2[∂2E/∂N2)] v (9)
The softness is defined as inverse of hardness (S=1/η).
Parr and Pearson [34] invoking the calculus of finite difference approximation suggested the approximate and operational formulae of electro negativity and hardness as under
χ = (I + A)/2 (10)
η = (I - A)/2 (11)
where l and A are the first ionization potential and electron affinity of the chemical species. Pearson [35] proceeded further to evaluate I and A in terms of orbital energies of the highest occupied molecular orbital, HOMO and the lowest unoccupied molecular orbital, LUMO as
χ = (€HOMO -€LUMO)/2 (12)
η = (-€HOMO + -€LUMO)/2 (13)
The electrophilicity, w is a descriptor of reactivity that allows a quantitative classification of the global electrophilic nature of a molecule within a relative scale. Parr et al. [36] suggested that electro negativity squared divided by hardness measures the electrophilic power of a ligand its prosperity to "soak up” electrons.
Thus,
ω = μ2/2η (14)
It is further anticipated that electrophilicity,w should be related to electron affinity, because both w and electron affinity measures capacity of an agent to accept electrons. Electron affinity reflect capability of an agent to accept only one electron from the environment, whereas electrophilicity index measures the energy lowering of a ligand due to maximal electron flow between the donor and acceptor The electron flows may be either less or more than one. Thus the electrophilicity index provides the direct relationship between the rates of reaction and the electrophilic power of the inhibitors [37].
Fukui functions play a prominent role in the field known as conceptual Density Functional Theory (DFT). Parr and Yang [38] based on the original ideas of Fukui [39], introduced fukui function which reflect the response of a molecular system towards a change in the number of electrons (Ne) of the molecular system under consideration .The fukui functions is a measure of local reactivity and defined as:
f (r) = ( ∂ ρ (r) /∂N )v (15)
where ρ(r) is the electron density.
Parr et al. [38] gave us the statement "the preferred direction is the one with largest f(r) at the direction side" and their firm prediction was -
Governing electrophilic attack f- (r)= [ ∂ ρ (r) / ∂N]- v (16)
Governing nuclophilic attack f+ (r)= [ ∂ ρ (r) / ∂N] +v (17)
Governing nutral attack f 0(r)= [ ∂ ρ (r) / ∂N] 0v (18)
These three cases have μs> μr, μs< μr and μs~ μ'r.
A "frozen core" approximation reveals dρ = dρvalue in each case and therefore governing electrophilic attack, f- (r) ≈ ρ (HOMO) (r) (19 ),
governing nuclophilic attack, f+(r) ≈ ρ (LUMO) (r) (20 ),
and governing neutral attack, f 0(r) ≈ 1/2[ρ (HOMO) (r) + ρ (LUMO) (r)] (21).
Method of Computation
We have used semi-empirical computational procedure, AM1 [31] to compute the orbital energies (HOMO, LUMO). Side by side we have also adopted a sophisticated ab-initio quantum chemical calculation procedure to compute the orbital energies (HOMO, LUMO). To be more specific, the planar structure of each of the substituted biphenyls was drown on the programme ArgusLab4.0.1[40] and geometry optimization was carried out at the restricted Hartree-Fock level (RHF) using STO-3G minimal basis set to compute the orbital energies (HOMO, LUMO) using AM1 procedure[31]. The ab-initio values are obtained when geometry optimization of the planar structure of each of the substituted biphenyls were carried out at restricted HartreeFock level (RHF) using STO-3G minimal basis set on the Hyper Chem 8.06 software [41].
The Ionization energies and electron affinities are calculated using the Koopmans theorem. The hardness value was calculated using Parr and Pearson method [34]. The electrophilicity index was calculated using the Parr et al's formula [36]. HOMO and LUMO Eigen functions obtained in the AM1 calculations were used to compute the fukui functions using Parr et al formulae [38].
Results and Discussion
A look on the Tables 1a &1b reveal that, in both computational procedure-first principle and AM1, the ionization energy is maximum in the case of chloro substituted Benzo[f]quinazolin- 1(2H)-one and it is minimum in case of Benzo[f]quinazolin- 1(2H)-one. Further scrutiny of the table reveals that the electron affinity is maximum for NH2 substituted Benzo[f]quinazolin- 1(2H)-one.and minimum for Benzo[f]quinazolin-1(2H)-one The electro negativity is maximum for Cl substituted Benzo[f] quinazolin-1(2H)-one and minimum for OH Benzo[f]quinazolin- 1(2H)-one. The hardness is maximum for Cl substituted and minimum f Benzo[f]quinazolin-1(2H)-one or OH substituted. The chemical potential is maximum for OH substituted Benzo[f] quinazolin-1(2H)-one Benzo[f]quinazolin-1(2H)-one.andminimum for H Benzo[f]quinazolin-1(2H)-one.


But we have surprisingly noted that the electro negativity value of all compounds studied by first principle on the minimal basis set STO-3G is negative. These do not agree with our normal experiences. We have thus adopted semi-empirical AM1 calculation on the same compounds. The semi-empirical results are discussed on the foregoing sections.

A look on Table 2 reveals that the centre 2C has the maximum f+ value and the centre12N has the minimum f+ value. Also it is found that the centrellC has the maximum f- value and the centre 15O has the minimum f- value in case of Benzo[f]quinazolin-1(2H)-one. A look on Table 3 reveals that the centre 4C has the maximum f+ value and the centre...17C has the minimum f+ value. Also it is found that the centre 3C has the maximum f- value and the centre 6C has the minimum f-value in case of Cl substitution Benzo[f]quinazolin-1(2H)-one.


A look on Table 4 reveals that the centre 17N has the maximum f+ value and the centre 11C has the minimum f+ value. Also it is found that the centre 3C has the maximum f- value and the centre 1C has the minimum f-value in case of NH2 substitution Benzo[f]quinazolin-1(2H)-one. A look on Table 5 reveals that the centre 2C has the maximum f+ value and the centre 8C has the minimum f+ value. Also it is found that the centre 3C has the maximum f- value and the centre 1C has the minimum f-value in case of CH3 substitution Benzo[f]quinazolin-1(2H)-one .


A deeper look on the Tables 3-6 reveals that the fukui functions on the different atoms of the molecules studied depends upon the substituents. The nucleophilic fukui functions of the 12 N atom varies in the following order- Benzo[f]quinazolin- 1(2H)-one OH->Cl>CH3>H>NH2 Benzo[f]quinazolin-1(2H)-one .The electrophilic fukui function of the 12 N atom varies in the following order- H>NH2l>CH3>Cl>OH Benzo[f]quinazolin-1(2H)- one .The nucleophilic fukui function of the 14 N atom varies in the following order- Cl>CH3>OH>H>NH2 Benzo[f]quinazolin- 1(2H)-one. The electrophilic fukui function of the 14 N atom varies in the following order- H>Cl>CH3>OH>NH2 Benzo[f] quinazolin-1(2H)-one. The nucleophilic fukui function on the 15 O atom varies in the following order-OH-Benzo[f]quinazolin- 1(2H)-one>Cl-Benzo[f]quinazolin-1(2H)-one>CH3- Benzo[f] quinazolin-1(2H)-one >H- Benzo[f]quinazolin-1(2H)-one>NH2- Benzo[f]quinazolin-1(2H)-one.
The electrophilic fukui function on the 15 O atom varies in the following order-CH3 Benzo[f]quinazolin-1(2H)-one >Cl Benzo[f]quinazolin-1(2H)-one >OH Benzo[f]quinazolin-1(2H)- one >NH2>H Benzo[f]quinazolin-1(2H)-one. The nucleophilic fukui function on the 16N atom varies in the following order-
Cl- Benzo[f]quinazolin-1(2H)-one> Benzo[f]quinazolin- 1(2H)-one>OH- Benzo[f]quinazolin-1(2H)-one >CH3- Benzo[f] quinazolin-1(2H)-one> NH2- Benzo[f]quinazolin-1(2H)-one. The electrophilic fukui function on the 16N atom varies in the following order-NH2- Benzo[f]quinazolin-1(2H)-one>OH- Benzo[f]quinazolin-1(2H)-one>CH3 Benzo[f]quinazolin-1(2H)- one >Cl- Benzo[f]quinazolin-1(2H)-one> Benzo[f]quinazolin- 1(2H)-one.
We have found a linear correlation between the global reactivity parameters of conceptual density functional theory and their drug activities.
We have correlated the global reactivity parameters of conceptual density functional theory with the drug activities of the Benzo[f]quinazolin-1(2H)-ones. We have found excellent correlation between the global hardness of the molecules and their drug activity.
Our suggested ansatz for the drug activity are‑
Drug Activity = ‐0.1289rηab‐initio + 1.1378 (22]
Drug Activity =‐0.1919 χab‐initio +0.3446 (23)
Drug Activity = 0.1919 μab-initio +.3446 (24)
Drug Activity = 0.081ωab‐initio + 1.3171 (25)
Drug Activity = 0.0086Sab‐initio + 0.3268 (26)
And using AM1 reactivity parameters
Drug Activity = ‐6.00386ηAM1 + 1.0789 (2 7)
Drug Activity =2.2924 χAM1 +0.1481 (28)
Drug Activity = ‐2.2924 μAM1 ‐0.1481 (29)
Drug Activity = 0.0894ωAM1 + 0.2604 (30)
Drug Activity = 0.0688SAM1 + 1.0789 (31)




We have computed the drug activity of the Benzo[f] quinazolin-1(2H)-ones using the above ansatz 1,2 and 4 and 6,7 and 9 and presented them in Tables 6-12. As the chemical potential is the negative counterpart of electro negativity and the softness is the inverse of hardness only, we have not considered them for computing the activity value. A deep scrutiny of the Tables 6-12 reveal that the theoretical activity value correlate well with the experimental values. In the case of hardness, this correlation is the best among the reactivity parameters in both ab-initio and AM1 computed results. We have drown Figures 1 & 2 to explore the correlation between theoretical computed and experimentally determined drug activity value of the Benzo[f] quinazolin-1(2H)-one using ab-initio and semi empirical hardness values respectively. From Figures 1 & 2 we can see that the activity data of the four compound correlate excellently with the experimental values. In the case of AM1 hardness data excellent correlation found in two values and other values fairly correlate with the experimental data.








Conclusion
Using global and local reactivity descriptors of DFT and MO theory, we have found a nice correlation between theoretically computed and experimentally determined drug activities (MIC). The fukui functions and local softness's at different sites changes with the change in substituents. This study once again reveals the fundamental statement of theoretical chemistry that "the structure of a molecule determines its electronic wave function which in turn determines many of its physicochemical properties”. It is reasonable to suppose that the biological activity of molecules its potency as a drug for examples will also, often enough, be dependent upon the molecular electronic wave function. It is possible in principle to determine any molecular property that is ultimately dependent upon the molecule's electronic ground state wave function. In practice this can sometimes be achieved using rigorous procedures.
References
- Smeyers Y G, L Bouniam, N J Smeyers, A Ezzamarty, A Hernandez Laguna et al. (19980 Quantum mechanical and QSAR study of some α-arylpropionic acids as anti-inflammatory agents. Eur J Med Chem 33: 103.
- Parr RG, W Yang (1989) Density Functional Theory of Atoms and Molecules. Oxford University Press, New York, USA.
- Kohn W, AD Becke, RG Parr (1996) Density Functional Theory of Electronic Structure. J Phys Chem 100(31): 12974.
- Hohenberg P, W Kohn (1964) Phys Rev B136: 864.
- Chattaraj PK, A Cedillo, R G Parr (1991) J Phys Chem 103: 7645.
- Ayers PW, RG Parr (2000) Variational Principles for Describing Chemical Reactions: The Fukui Function and Chemical Hardness Revisited. J Am Chem Soc 122(9): 2010-2018.
- De Proft F, J M L Martin, P Geerlings (1996) On the performance of density functional methods for describing atomic populations, dipole moments and infrared intensities. Chem Phys Let 250: 393-401.
- Geerlings P, De Proft F, Martin JML, (1996) In Theoretical and Computational Chemistry; Seminario. J Ed Elsveir Amsterdam 4: 773.
- De Proft F, Martin J M L, Geerlings P (1996) Calculation of molecular electrostatic potentials and Fukui functions using density functional methods. Chem Phys Let 256: 400-408.
- De Proft F, Geerlings P (1997) Calculation of ionization energies, electron affinities, electronegativities, and hardnesses using density functional methods. J Chem Phys 106: 3270.
- Geerlings P, De Proft F, Langenaeker W (1996) Density Functional Theory : A Source of Chemical Concepts and a Cost-Effective Methodology for Their Calculation. Adv Quantum Chem 33: 303-328.
- Parr R G, R A Donnelly, Levy M, Palke W E (1978) Electronegativity: The density functional viewpoint. J Chem Phys 68(8): 3801.
- Hansch C, PG Sammes, Taylor J B (1990) Computers and the medicinal chemist; in:Comprehensive Medicinal Chemistry, vol. 4, Eds. Pergamon Press, Oxford, pp. 33
- Franke R (1984) Theoretical Drug Design Methods. Elsevier Amsterdam.
- Sawhney S N, Tower R K, Singh S P (1980) Indian J Chem 19: 415.
- Ana B, Boteanu, Farmacia S (1971) 19: 683.
- Dempcy R O, Skibo E B (1993) The synthesis and structural characterization of way-120,491; a novel potassium channel activator Bioorg Med Chem Lett 1: 39.
- Bhargava P N, Singh G C J (1961) Indian Chem Soc 38: 77.
- Gujral M L, Saxena PN, Kohli R P (1957) Indian J Med Res 45: 201.
- Griffin J, Srinivasan S, Bowman K, Kalvert H A, Curtin N J, et al. (1998) B T J Med Chem 41: 5247.
- Carloni P, Alber F (Eds.) (2002) Quantum Medicinal Chemistry, WILEY- VCH Verlag GmbH & Co. KGaA, Wenham.
- Karelson M (2000) Molecular descriptors in QSAR/QSPR, Wiley Inter Science. New York, USA.
- Hammett LP (1935) Some Relations between Reaction Rates and Equilibrium Constants. chem Rev 17:125-136.
- Hasanein A A (1996) Computational method in quantum chemistry. world scientific.
- Hoffmann R J (1964) Chem Phys 40: 2480.
- Pople JA, Santry DP, Segal GA (1965) J Chem Phys 43: S29.
- Pople JA, Beveridge DL, Dobosh PA (1967) Approximate Self-Consistent Molecular-Orbital Theory. V. Intermediate Neglect of Differential Overlap. J Chem Phys 47: 2026.
- Dewar MJS, Kiopman G (1967) Ground states of .sigma.-bonded molecules. I. Semiempirical S.C.F. molecular orbital. J Am Chem Soc 89: 3089-3098.
- Dewar MJS,Thiel W (1977) Ground states of molecules. 38. The MNDO method. Approximations and parameters. J Am Chem Soc 99: 48994907.
- Dewar MJS, Zoebisch EG (1988) J Mol Struct THEOCHEM 49: 1.
- StewartJ JP (1990) Reply to "Comments on a comparison of AM1 with the recently developed PM3 method”. J ComputChem11: 543.
- Gyftpoulous EP, Hatsopoulos GN (1968) Quantum-Thermodynamic Definition Of Electronegativity . Proc Natl Acad Sci 60(3): 786-793.
- Iczkowski RP, Margrave JL (1961) Electronegativity. J Am Chem Soc 83: 3547-3551.
- Parr RG, Pearson RG (1983) Absolute hardness: companion parameter to absolute electronegativity. J Am Chem Soc105: 7512-7516.
- Pearson RG (1986) beta-Glucuronidase from Escherichia coli as a genefusion marker. Proc Natl Acad Sci 83: 8447-8451.
- Parr RG, Szentpaly LV, Liu S (1999) Electrophilicity Index. J Am Chem Soc 121:1922.
- Maynard AT, Covell DG (2001) Reactivity of Zinc Finger Cores: Analysis of Protein Packing and Electrostatic Screening. J Am ChemSoc123: 1047-1058.
- Parr RG, YangW (1984) Density functional approach to the frontier- electron theory of chemical reactivity. J Am Chem Soc 106: 4049.
- Fukui K, Yonezawa T, Shingu H (1952) A Molecular Orbital Theory of Reactivity in Aromatic Hydrocarbons. J Chem Phys 20: 722.
- Stewart J J P (1989) Optimization of parameters for semiempirical methods I. Method. J Comput Chem 10: 209-220.
- Thompson M A (ArgusLab Version 4.0.1, Planaria Software LLC Seattle, WA. http://www.planaria-software.com.