Articles

When a drug interacts with its receptor, the nonpolar substituent of drug and receptor proteins attract each other because they have opposite magnitude with respect to each other. X-rays structure studies reflected that the S2/S2’ pocket in HIV-1 protease enzyme are essentially hydrophobic. The residues that make up these pockets are Val-32, Ile-47, Ile-50, and Ile-84 in each monomeric polypeptidic unit of the protease enzyme. Δπ_dr and ΔSASA_dr have been used to measure the extent of hydrophobic interaction between peptidic protease inhibitors and receptor proteins (binding site: valine‒isoleucine; and catalytic site: glycine‒aspartic acid‒threonine) on the HIV-1 protease enzyme. For measurement of hydrophobic interaction, the molecular modeling and geometry optimization of all the inhibitors and the receptor amino acids have been carried out with CAChe Pro software by opting semiempirical PM3 methods. Log P was calculated using the atom-typing scheme of Ghose and Crippen, while solvent accessible surface area by conductor likes screening model. π_d, π_r, SASA_Sd and SASA_Srwell describe the hydrophobicities of the substituents and play the effective role for site selectivity for interaction of the drug with the receptor. Comparative study of values of Δπ_dr and ΔSASA_dr show the order of hydrophobic interaction with respect to amino acids: Asp > Thr > Val > Ile and Thr > Val > Asp > Ile, respectively. Further, comparative study of the values of (ΣΔπ_dr)_binding-site, (ΣΔπ_dr)_{catalytic-site}, (ΣΔSASA_dr)_binding-site, (ΣΔSASA_dr)_{catalytic-site} shows that peptidic HIV-1-PRIs interact with binding site rather than catalytic site as binding site have lower value of ΣΔπ_dr and ΣΔSASA_dr. Among the binding site, Val has maximum interaction than Ile, as it has lower vale of Δπ_dr and ΔSASA_dr.

In this research work, we have applied “Lipinski’s RO5” for pharmacokinetics (PK) study and to predict the activity of peptidic HIV-1 protease inhibitors. Peptidic HIV-1-PRIs have been taken from literature with their observed biological activities (OBAs) in term of IC50. The logarithms of the inverse of IC50 have been used as biological end point o(log1/C) in the study. For calculation of physicochemical parameters, the molecular modeling and geometry optimization of all the derivatives have been carried out with CAChe Pro software using semiempirical PM3 method. Prediction of the biological activity of the inhibitors has shown that the best QSAR model is constructed from pharmacokinetic properties, molecular weight and hydrogen bond acceptor. This also proved that these properties play important role to describe the PKs of the drugs. On the basis of the derived models one can build up a theoretical basis to access the biological activity of the compounds of the same series.