Analogue and structure based approaches for modelling HIV-1 integrase inhibitors

A set of 220 inhibitors belonging to different structure classes and having HIV-1 integrase activity were collected along with their experimental pIC50 values. Geometries of all the inhibitors were fully optimized using B3LYP/6-31 + G(d) level of theory. These ligands were docked against 4 different HIV-1 integrase receptors (PDB IDs: 4LH5, 5KRS, 3ZSQ and 3ZSV). 30 docked poses were generated for all 220 inhibitors and ligand interaction of the first docked pose and the docked pose with the highest score were analysed. Residue GLU170 of 4LH5 receptor shows the highest number of interactions followed by ALA169, GLN168, HIS171 and ASP167 residues. Hydrogen bonding and stacking are mainly responsible for the interactions of these inhibitors with the receptor. We performed Molecular Dynamics (MD) simulation to observe the root-mean-square deviation (RMSD), for measure the average change of displacement between the atoms for a particular frame with respect to a reference and The Root Mean Square Fluctuation (RMSF) for characterization of local changes along the protein chain of the docked complexes. Analogue based models were generated to predict the pIC50 values for integrase inhibitors using various types of descriptors such as constitutional, geometrical, topological, quantum chemical and docking based descriptors. The best models were selected on the basis of statistical parameters and were validated by training and test set division. A few new inhibitors were designed on the basis of structure activity relationship and their pIC50 values were predicted using the generated models. All the designed new inhibitors a very high potential and may be used as potent inhibitors of HIV integrase. These models may be useful for further design and development of new and potent HIV integrase inhibitors.Communicated by Ramaswamy H. Sarma.

Artificially expanded genetic information systems (AEGISs) as potent inhibitors of the RNA-dependent RNA polymerase of the SARS-CoV-2

The recent outbreak of the SARS-CoV-2 infection has affected the lives and economy of more than 200 countries. The unavailability of virus-specific drugs has created an opportunity to identify potential therapeutic agents that can control the rapid transmission of this pandemic. Here, the mechanisms of the inhibition of the RNA-dependent RNA polymerase (RdRp), responsible for the replication of the virus in host cells, are examined by different ligands, such as Remdesivir (RDV), Remdesivir monophosphate (RMP), and several artificially expanded genetic information systems (AEGISs) including their different sequences by employing molecular docking, MD simulations, and MM/GBSA techniques. It is found that the binding of RDV to RdRp may block the RNA binding site. However, RMP would acquire a partially flipped conformation and may allow the viral RNA to enter into the binding site. The internal dynamics of RNA and RdRp may help RMP to regain its original position, where it may inhibit the RNA-chain elongation reaction. Remarkably, AEGISs are found to obstruct the binding site of RNA. It is shown that dPdZ, a two-nucleotide sequence containing P and Z would bind to RdRp very strongly and may occupy the positions of two nucleotides in the RNA strand, thereby denying access of the substrate-binding site to the viral RNA. Thus, it is proposed that the AEGISs may act as novel therapeutic candidates against the SARS-CoV-2. However, in vivo evaluations of their potencies and toxicities are needed before using them against COVID-19.Communicated by Ramaswamy H. Sarma.

Rational design of pincer-nickel complexes for catalytic cyanomethylation of benzaldehyde: A systematic DFT study

The current study dwells upon the efforts to computationally probe a phosphine-free pincer-nickel complex that would demonstrate an efficiency better than the reported phosphine-based pincer-nickel complex (^iPr2 POCN^Et2 )Ni(CH₂ CN) for cyanomethylation reaction. For this purpose, the mechanism of cyanomethylation of benzaldehyde was studied quantum mechanically for a series of 11 pincer-nickel complexes. The energetics of various intermediates and transition states involved in the catalytic cycle for each catalyst was compared with the corresponding energetics of the Miller's catalyst (^iPr2 POCN^Et2 )Ni(CH₂ CN) that is reported to accomplish the cyanomethylation at room temperature. While pincer complexes (^iPr4 NNN)Ni(CH₂ CN) and (^iPr4 NCN)Ni(CH₂ CN) containing strong σ-donating amines were found to fare poorly, pincer-nickel complexes (^iPr2 NCN)Ni(CH₂ CN) and (^dm PheboxNCN)Ni(CH₂ CN) based on weaker σ-donating imines had energetics more favorable than the reported efficient catalyst (^iPr2 POCN^Et2 )Ni(CH₂ CN). While strong trans-influencing C as the pincer central atom was found to be pivotal for lowering the cyanomethylation kinetics, presence of a poor trans-influencing N proved to be detrimental on the overall energetics.

Peptide-like and small-molecule inhibitors against Covid-19

Coronavirus disease strain (SARS-CoV-2) was discovered in 2019, and it is spreading very fast around the world causing the disease Covid-19. Currently, more than 1.6 million individuals are infected, and several thousand are dead across the globe because of Covid-19. Here, we utilized the in-silico approaches to identify possible protease inhibitors against SARS-CoV-2. Potential compounds were screened from the CHEMBL database, ZINC database, FDA approved drugs and molecules under clinical trials. Our study is based on 6Y2F and 6W63 co-crystallized structures available in the protein data bank (PDB). Seven hundred compounds from ZINC/CHEMBL databases and fourteen hundred compounds from drug-bank were selected based on positive interactions with the reported binding site. All the selected compounds were subjected to standard-precision (SP) and extra-precision (XP) mode of docking. Generated docked poses were carefully visualized for known interactions within the binding site. Molecular mechanics-generalized born surface area (MM-GBSA) calculations were performed to screen the best compounds based on docking scores and binding energy values. Molecular dynamics (MD) simulations were carried out on four selected compounds from the CHEMBL database to validate the stability and interactions. MD simulations were also performed on the PDB structure 6YF2F to understand the differences between screened molecules and co-crystallized ligand. We screened 300 potential compounds from various databases, and 66 potential compounds from FDA approved drugs. Cobicistat, ritonavir, lopinavir, and darunavir are in the top screened molecules from FDA approved drugs. The screened drugs and molecules may be helpful in fighting with SARS-CoV-2 after further studies.Communicated by Ramaswamy H. Sarma.

Synthesis of β- and γ-lactam fused dihydropyrazinones from Ugi adducts via a sequential ring construction strategy

A modular approach for the construction of β- and γ-lactam fused dihydropyrazinones from the readily available Ugi adducts has been described. The sequential construction of rings through base-mediated cycloisomerization followed by acid-mediated cyclization yielded β-lactam fused dihydropyrazinones. However, the Ugi-derived dihydropyrazinones afforded γ-lactam fused dihydropyrazinones under base-mediated cycloisomerization. The regioselectivity in the cycloisomerization reactions is explained on the basis of ring-strain. Substrate scope, limitations and mechanistic investigations through DFT-calculations have been explored.

pH-Regulated anion transport activities of bis(iminourea) derivatives across the cell and vesicle membrane

Recently, synthetic anion transporters have gained considerable attention because of their ability to disrupt cellular anion homeostasis and promote cell death. Herein, we report the development of bis(iminourea) derivatives as a new class of selective Cl- ion carrier. The bis(iminourea) derivatives were synthesized via a one-pot approach under mild reaction conditions. The presence of iminourea moieties suggests that the bis(iminourea) derivatives can be considered as unique guanidine mimics, indicating that the protonated framework could have much stronger anion recognition properties. The cooperative interactions of H+ and Cl- ions with these iminourea moieties results in the efficient transport of HCl across the lipid bilayer in an acidic environment. Under physiological conditions these compounds weakly transport Cl- ions via an antiport exchange mechanism. This pH-dependent gating/switching behavior (9-fold) within a narrow window could be due to the apparent pKa values (6.2-6.7) of the compounds within the lipid bilayer. The disruption of ionic homeostasis by the potent compounds was found to induce cell death.

Phosphine-free pincer-ruthenium catalyzed biofuel production: high rates, yields and turnovers of solventless alcohol alkylation

High TONs and TOFs are observed for the β-alkylation of alcohols using phosphine-free pincer-ruthenium catalysts at a very low base loading. Kinetic studies and DFT calculations were complementary and provide a clear understanding on the mechanism.

Selective and high yield transformation of glycerol to lactic acid using NNN pincer ruthenium catalysts

A sterically less hindered 2,6-bis(benzimidazol-2-yl)pyridine based pincer–ruthenium complex has been used here to accomplish the catalytic conversion of glycerol selectively to lactic acid in high yield.

Benchmarking of different molecular docking methods for protein-peptide docking

Background

Molecular docking studies on protein-peptide interactions are a challenging and time-consuming task because peptides are generally more flexible than proteins and tend to adopt numerous conformations. There are several benchmarking studies on protein-protein, protein-ligand and nucleic acid-ligand docking interactions. However, a series of docking methods is not rigorously validated for protein-peptide complexes in the literature. Considering the importance and wide application of peptide docking, we describe benchmarking of 6 docking methods on 133 protein-peptide complexes having peptide length between 9 to 15 residues. The performance of docking methods was evaluated using CAPRI parameters like FNAT, I-RMSD, L-RMSD.

Result

Firstly, we performed blind docking and evaluate the performance of the top docking pose of each method. It was observed that FRODOCK performed better than other methods with average L-RMSD of 12.46 Å. The performance of all methods improved significantly for their best docking pose and achieved highest average L-RMSD of 3.72 Å in case of FRODOCK. Similarly, we performed re-docking and evaluated the performance of the top and best docking pose of each method. We achieved the best performance in case of ZDOCK with average L-RMSD 8.60 Å and 2.88 Å for the top and best docking pose respectively. Methods were also evaluated on 40 protein-peptide complexes used in the previous benchmarking study, where peptide have length up to 5 residues. In case of best docking pose, we achieved the highest average L-RMSD of 4.45 Å and 2.09 Å for the blind docking using FRODOCK and re-docking using AutoDock Vina respectively.

Conclusion

The study shows that FRODOCK performed best in case of blind docking and ZDOCK in case of re-docking. There is a need to improve the ranking of docking pose generated by different methods, as the present ranking scheme is not satisfactory. To facilitate the scientific community for calculating CAPRI parameters between native and docked complexes, we developed a web-based service named PPDbench (http://webs.iiitd.edu.in/raghava/ppdbench/).

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2449-y) contains supplementary material, which is available to authorized users.

Secondary amine salt catalyzed controlled activation of 2-deoxy sugar lactols towards alpha-selective dehydrative glycosylation

A new organocatalytic glycosylation method exploiting the lactol functionality has been disclosed. The catalytic generation of glycosyl oxacarbenium ions from lactols under forcible conditions via weakly Brønsted-acidic, readily available secondary amine salts affects the diastereoselective glycosylation of 2-deoxypyranoses and furanoses. This operationally simple iminium catalyzed activation of 2-deoxy hemi-acetals is a potential alternative to the existing cumbersome methods that need specialized handling. The mechanisms for this unique transformation and kinetic/thermodynamic effects have been discussed based on both experimental evidence and theoretical studies.

Trisulfides over disulfides: highly selective synthetic strategies, anti-proliferative activities and sustained H2S release profiles

Highly selective synthesis of trisulfides over disulfides is demonstrated along with their potential as anti-proliferative agents and sustained donors of H₂S.

One- and Two-Photon Uncaging: Carbazole Fused o-Hydroxycinnamate Platform for Dual Release of Alcohols (Same or Different) with Real-Time Monitoring

A one- and two-photon activated photoremovable protecting group (PRPG) was designed based on a carbazole fused o-hydroxycinnamate platform for the dual (same or different) release of alcohols. The mechanism for the dual release proceeds through a stepwise pathway and also monitors the first and second photorelease in real time by an increase in fluorescence intensity and color change, respectively. Further, its application in staining live neurons and ex vivo imaging with two-photon excitation is shown.

Metal-Free Sequential C(sp2)-H/OH and C(sp3)-H Aminations of Nitrosoarenes and N-Heterocycles to Ring-Fused Imidazoles

Hydrogen bond assisted ortho-selective C(sp²)-H amination of nitrosoarenes and subsequent α-C(sp³)-H functionalization of aliphatic amines is achieved under metal-free conditions. The annulation of nitrosoarenes and 2-hydroxy-C-nitroso compounds with N-heterocycles provides a facile excess to a wide range of biologically relevant ring-fused benzimidazoles and structurally novel polycyclic imidazoles, respectively. Nucleophilic aromatic hydrogen substitution (S_NArH) was found to be preferred over classical S_NAr reaction during the C(sp²)-H amination of halogenated nitrosoarenes.

A web server for analysis, comparison and prediction of protein ligand binding sites

BACKGROUND

One of the major challenges in the field of system biology is to understand the interaction between a wide range of proteins and ligands. In the past, methods have been developed for predicting binding sites in a protein for a limited number of ligands.

RESULTS

In order to address this problem, we developed a web server named 'LPIcom' to facilitate users in understanding protein-ligand interaction. Analysis, comparison and prediction modules are available in the "LPIcom' server to predict protein-ligand interacting residues for 824 ligands. Each ligand must have at least 30 protein binding sites in PDB. Analysis module of the server can identify residues preferred in interaction and binding motif for a given ligand; for example residues glycine, lysine and arginine are preferred in ATP binding sites. Comparison module of the server allows comparing protein-binding sites of multiple ligands to understand the similarity between ligands based on their binding site. This module indicates that ATP, ADP and GTP ligands are in the same cluster and thus their binding sites or interacting residues exhibit a high level of similarity. Propensity-based prediction module has been developed for predicting ligand-interacting residues in a protein for more than 800 ligands. In addition, a number of web-based tools have been integrated to facilitate users in creating web logo and two-sample between ligand interacting and non-interacting residues.

CONCLUSIONS

In summary, this manuscript presents a web-server for analysis of ligand interacting residue. This server is available for public use from URL http://crdd.osdd.net/raghava/lpicom .

Modeling anti-HIV compounds: the role of analogue-based approaches

There has been a tremendous progress in the development of anti-HIV therapies since the discovery of the HIV virus. Computer aided drug design in general and analogue-based approaches in particular have played an important role in the process of HIV drug discovery. Structure-based approaches also have played a vital role in this process. There are a large number of studies reported in the literature where QSAR methodology was employed to study the structural requirements for inhibition against various HIV targets like reverse transcriptase, protease, entry and integrase. The current review focuses on those studies and provides a detailed description on the QSAR methodology, descriptors, statistical significance and important findings. This review categorizes the reported QSAR studies on the basis of chemical scaffolds against a particular target. In reverse transcriptase category, QSAR studies on HEPT, TIBO, DABO, DAPY, DATA, AASBN, pyridone and DATZD derivatives have been reviewed. Cyclic urea, fullerene, AHPBA and dihydropyrone derivatives were considered in protease inhibitors category. In addition, QSAR studies on styrylquinoline, carboxylic acid, MBSA and chalcone derivatives were reviewed in integrase inhibitors category. QSAR studies on entry inhibitors like piperidine, benzyl piperidine, benzyl pyrazole, pyrrole and diazepane urea have also been reviewed.

Molecular dynamics investigation on a series of HIV protease inhibitors: assessing the performance of MM-PBSA and MM-GBSA approaches

The binding free energies (ΔG(Bind)) obtained from molecular mechanics with Poisson-Boltzmann surface area (MM-PBSA) or molecular mechanics with Generalized Born surface area (MM-GBSA) calculations using molecular dynamics (MD) trajectories are the most popular procedures to measure the strength of interactions between a ligand and its receptor. Several attempts have been made to correlate the ΔG(Bind) and experimental IC(50) values in order to observe the relationship between binding strength of a ligand (with its receptor) and its inhibitory activity. The duration of MD simulations seems very important for getting acceptable correlation. Here, we are presenting a systematic study to estimate the reasonable MD simulation time for acceptable correlation between ΔG(Bind) and experimental IC(50) values. A comparison between MM-PBSA and MM-GBSA approaches is also presented at various time scales. MD simulations (10 ns) for 14 HIV protease inhibitors have been carried out by using the Amber program. MM-PBSA/GBSA based ΔG(Bind) have been calculated and correlated with experimental IC(50) values at different time scales (0-1 to 0-10 ns). This study clearly demonstrates that the MM-PBSA based ΔG(Bind) (ΔG(Bind)-PB) values provide very good correlation with experimental IC(50) values (quantitative and qualitative) when MD simulation is carried out for a longer time; however, MM-GBSA based ΔG(Bind) (ΔG(Bind)-GB) values show acceptable correlation for shorter time of simulation also. The accuracy of ΔG(Bind)-PB increases and ΔG(Bind)-GB remains almost constant with the increasing time of simulation.

Calculation of transition dipole moment in fluorescent proteins--towards efficient energy transfer

Förster Resonance Energy Transfer (FRET) between fluorescent proteins (FPs) is widely used to construct fluorescent sensor proteins, to study intracellular protein-protein interactions and to monitor conformational changes in multidomain proteins. Although FRET depends strongly on the orientation of the transition dipole moments (TDMs) of the donor and acceptor fluorophores, this orientation dependence is currently not taken into account in FRET sensor design. Similarly, studies that use FRET to derive structural constrains typically assume a κ(2) of 2/3 or use the TDM of green fluorescent protein, as this is the only FP for which the TDM has been determined experimentally. Here we used time-dependent density functional theory (TD-DFT) methods to calculate the TDM for a comprehensive list of commonly used fluorescent proteins. The method was validated against higher levels of calculation. Validation with model compounds and the experimentally determined TDM of GFP shows that the TDM is mostly determined by the structure of the π-conjugated fluorophore and is insensitive to non-conjugated side chains or the protein surrounding. Our calculations not only provide TDM for most of the currently used FPs, but also suggest an empirical rule that can be used to obtain the TDMs for newly developed fluorescent proteins in the future.

Analogue-based approaches in anti-cancer compound modelling: the relevance of QSAR models

BACKGROUND

QSAR is among the most extensively used computational methodology for analogue-based design. The application of various descriptor classes like quantum chemical, molecular mechanics, conceptual density functional theory (DFT)- and docking-based descriptors for predicting anti-cancer activity is well known. Although in vitro assay for anti-cancer activity is available against many different cell lines, most of the computational studies are carried out targeting insufficient number of cell lines. Hence, statistically robust and extensive QSAR studies against 29 different cancer cell lines and its comparative account, has been carried out.

RESULTS

The predictive models were built for 266 compounds with experimental data against 29 different cancer cell lines, employing independent and least number of descriptors. Robust statistical analysis shows a high correlation, cross-validation coefficient values, and provides a range of QSAR equations. Comparative performance of each class of descriptors was carried out and the effect of number of descriptors (1-10) on statistical parameters was tested. Charge-based descriptors were found in 20 out of 39 models (approx. 50%), valency-based descriptor in 14 (approx. 36%) and bond order-based descriptor in 11 (approx. 28%) in comparison to other descriptors. The use of conceptual DFT descriptors does not improve the statistical quality of the models in most cases.

CONCLUSION

Analysis is done with various models where the number of descriptors is increased from 1 to 10; it is interesting to note that in most cases 3 descriptor-based models are adequate. The study reveals that quantum chemical descriptors are the most important class of descriptors in modelling these series of compounds followed by electrostatic, constitutional, geometrical, topological and conceptual DFT descriptors. Cell lines in nasopharyngeal (2) cancer average R2 = 0.90 followed by cell lines in melanoma cancer (4) with average R2 = 0.81 gave the best statistical values.

High efficiency, heterosis, and homeostasis in mitochondria of wheat

Mitochondria of a wheat hybrid 31MS x 28 when assayed polarographically exceeded mitochondria of parents of the hybrid in ADP:O ratios and in respiratory control when they utilized alpha-ketoglutarate, malate, and succinate. The ADP:O ratios of the hybrid mitochondria for the three substrates were 5.8 +/- 0.2, 3.8 +/- 0.15, and 3.4 +/- 0.2, respectively. High P:O ratios were also observed in manometric assays. Mixtures of parental mitochondria exceeded averages of the parents, i.e., they exhibited complementation and approached activities of hybrid mitochondria. Complementation was observed only with NAD-linked substrates. Dinitrophenol uncoupled all mitochondria, but uncoupling was smallest in hybrid mitochondria. Hybrid mitochondria had the greatest ATPase activity. Mixtures exhibited complementation in ATPase activity. ATPase activity was stimulated least in hybrid mitochondria and in parental mixtures. The results were discussed from the point of view that the high P:) ratios gave us direct evidence that heterosis is indeed characterized by high efficiency. It is also pointed out that because of consistent repeatability and their specificity to hybrids, P:O ratios studied in rapidly isolated mitochondria of hybrids may give us revealing information on supramaximal efficiency of oxidative phosphorylation.

Can DFT methods correctly and efficiently predict the coordination number of copper(I) complexes? A case study

The coordination number of various experimentally known Cu(I) compounds is studied using density functional theory. Various basis sets are tested, aiming to establish a reliable level for prediction of the coordination number of these and other Cu(I) complexes. It is found that most levels exhibit correct trends, namely, the bulkier ligands demonstrate larger preference for coordination of two ligands. Proper absolute values are obtained when dispersion corrections are also included in the calculations. It is concluded that the fairly small modified 6-31+G* basis set due to Pulay represents a good compromise between accuracy and efficiency, followed by Balabanov and Peterson's all-electron aug-cc-pVDZ basis set. The overall energy is decomposed into various components whose relative contribution to the overall tendency of forming a complex with a particular coordination is examined. It is shown that two opposing contributions play a major role: the interaction energy of the ligand being added and the deformation energy of the copper's coordination sphere prior to the ligand addition. The former being a stabilizing contribution, leads to higher coordination numbers while the later, a destabilizing contribution, is shown to favor lower coordination numbers.

Formation of Diethyl 2-Amino-1-cyclopentenylphosphonates: A Simple Synthesis with a Unique Mechanism

Reaction of diethyl 5-chloro-1-pentynylphosphonate with primary and secondary amines produces novel 2-amino-1-cyclopentenylphosphonates, 1, in excellent isolated yields (79-88%). Calculations supported by experimental facts point to a two-step mechanism: an initial amine addition to give a zwitterionic intermediate followed by cyclization and proton transfer. Calculations rule out the formation of an enamine as an intermediate.

Honeybees foraging response in genetically diversified opium poppy

Studies were carried out on honeybees foraging on plant flowers. Results showed significantly higher foraging response of honeybees (Apis mellifera) in genetically divergent narcotic plant opium poppy (Papaver somniferum). Of the 18 mutants and two locally adapted cultivars of diverse genotypes screened, eight revealed significantly greater foraging response manifesting honeybee's preference towards specific plant morphotypes. The number of flower bloom did not correspond to number of foraging bees in both mutant and cultivar plant types of opium poppy. The genotype specific foraging response of honeybees could be attributed to physico-chemical properties of opium poppy flowers. This could have implications for the development of opium alkaloid fortified honeys for novel pharmaceuticals and isolation of natural spray compounds to attract honeybee pollinators for promoting crossing and sustainable hybridity in crops.

Comparative QSAR study of phenol derivatives with the help of density functional theory

Quantum chemical reactivity descriptors based QSAR study of 50 phenol derivatives is presented in this paper. Four different methods have been employed to certify the reliability of QSAR study. The molecular weight, hardness, chemical potential, total energy, and electrophilicity index provide valuable information and have a significant role in the assessment of the toxicity of phenols. The first model has been drawn up with the help of AM1 calculations and in this model the correlation coefficient r2 is 0.88 and the cross-validation coefficient r(cv)2 is 0.78. Second and third models have been designed with the PM3 and PM5 calculations, respectively. The values of correlation coefficient r2 and cross-validation coefficient r(cv)2 in the second case are 0.85 and .070, while in the third case they are 0.85 and 0.71. Finally, the DFT calculations have been made for the same series of compounds by using a B88-PW91 GGA energy functional with the DZVP basis set. The DFT models have a higher predictive power than AM1, PM3, and PM5 methods, and the reliability of this model is clear from its correlation coefficient r2 0.91 and cross-validation coefficient r(cv)2 0.88. This study is also helpful in determining the effect of any particular phenol derivative of this series over Tetrahymena pyriformis.

Semiempirical QSAR study and ligand receptor interaction of estrogens

Softness values E(n)+/+ of estrogen derivatives and softness values E(m)+/+ of receptor lysine, histidine, tyrosine and cysteine have been evaluated by Klopman equation. The required parameters for the solution of Klopman equation have been calculated with the help of PM3 method. The difference deltaE(nm)+/+ between E(n)+/+ and E(m)+/+ has been derived for QSAR study. The estrogen derivatives have been divided into four different sets on the basis of their structural similarities, and their biological activity taken from literature in terms of relative binding affinity (RBA). The QSAR study shows that, deltaE(nm)+/+ values provide good relationship with biological activity.

DFT-based QSAR study of testosterone and its derivatives

QSAR study of derivatives of testosterone has been made with the help of quantum mechanical parameters such as Absolute Hardness (eta) and Electronegativity (chi). These two parameters have been derived with the help of density functional theory. The 3-D modeling and geometry optimization of all the compounds have been done with the help of PCMODEL software and semiempirical PM3 calculations performed with the help of WinMOPAC-7.21 software. The absolute hardness provides valuable information due to maximum hardness principle and used in development of QSAR. The information provided by electronegativity is not as clear as in case of absolute hardness.