Drug-Receptor Interaction of Peptidic HIV-1 Protease: The Hydrophobic Effect-I
Abstract
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 ΔSASAdr 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, SASASd and SASASr well 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 ΔSASAdr 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, (ΣΔSASAdr)binding-site, (ΣΔSASAdr)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 ΣΔSASAdr. Among the binding site, Val has maximum interaction than Ile, as it has lower vale of Δπdr and ΔSASAdr.
1. Introduction
Interactions between the hydrophobic regions of a binding site and those of a complementary ligand are often observed to provide the driving force for binding [1]. In this work Δπdr and ΔSASAdr have been used to measure the extent of hydrophobic interaction between peptidic HIV-1 protease inhibitors (HIV-1-PRIs) and receptor proteins (binding and catalytic site) on the HIV-1 protease enzyme (HIV-1-PR). These approaches have been parametrized using πd (lipophilicity of the substituent of the drug), πr (lipophilicity of the substituent of the receptor), SASASd (solvent accessible surface area of the substituent of the drug) and SASASr (solvent accessible surface area of the substituent of the receptor). π and SASAS have been used to identify those groups/substituents that are more likely have a strong hydrophobic interaction with nonpolar regions [2]. These consider not only atom type but also the non-additive effects arising from the shape and extent of a nonpolar region. 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 as shown in Figure 1.
Drug-receptor interaction in term of Δπdr and ΔSASAdr are defined as:
where πd is the log P of the substituent of drug (peptidic HIV-1-PRI), πr the log P of the substituent of receptor (amino acids) and Δπdr is their change (higher - lower).
where SASASd is the SASAS of the substituent of drug (HIV-1-PRI), SASASr the SASA of the substituent of receptor and ΔSASAdr is their change (higher - lower).
Lower the values of Δπdr and ΔSASAdr higher will the interaction. The log P and SASA of the substituent have been evaluated by solving the following equations [3].
where log Psubs (π) is the log P of the substituent of the drug/ receptor, log Pderivt is the log P of the derivative and log PH is the log P of the lead compound, SASAsubs (or SASAS) is the SASA of the substituent of the drug/ receptor, SASAderivt is the SASA of the derivative and SASAH is the SASA of the lead compound. In this work, we have studied the hydrophobic interaction governing the drug-receptor interaction of fifty-one HIV-1-PRIs with their receptor on HIV-1-PR.
2. Materials and Methods:
The study materials of the present research work are fifty-one peptidic HIV-1-PRIs and receptor proteins. The inhibitors have been taken from literatures [4, 5] and are listed in Table 1 to 3 with their observed bioactivities in term of IC50 [6]. The logarithms of the inverse of IC50 have been used as biological end point (log1/C) in the study. Interaction sites on receptor proteins are amino acids. These are valine and isoleucine amino acids on binding site, while aspartic acid, threonine and glycine amino acids on catalytic site, and are shown in figure 2. 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 PM3 methods [7]. Log P is calculated using the atom-typing scheme of Ghose and Crippen [8], while SASA by conductor like screening model (COSMO) [9]. The majority of HIV research is done with cells and these studies tend to overestimate log P for animal system. A study of their database revealed that log P for cells is about 1 log unit higher than for whole organisms [10].
3. Results and Discussion:
Survey of literatures exposed that HIV-1-PR is a viral encoded homodimeric aspartyl protease with C2 symmetry [11, 12], and a catalytic triad of Asp-Thr-Gly contributed by each monomer that comprises the active site of the enzyme. Literatures also showed that all protease inhibitors bind to the protease binding site pocket that has a considerable number of hydrophobic residues and the X-rays structure indicates that the S2/S2’ pocket in HIV-1-PR are essentially hydrophobic [11, 12, 13]. 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 [13]. The importance of hydrophobic residues in the bonding pocket reconfirms the contribution of the hydrophobicity of inhibitors on anti-HIV activity. Since Val and Ile are hydrophobic in nature, they must be playing a decisive role in hydrophobic interaction. In this article, we have studied hydrophobic interaction governing the drug-receptor interaction of fifty-one peptidic HIV-1-PRIs with their receptor on HIV-1-PR. All the fifty-one inhibitors have been divided in three sets on the basis of their structural similarities (Figure 3). The first, second and third set comprises of eighteen, seventeen and sixteen inhibitors, respectively [4, 5].
3.1. Hydrophobicity/ Lipophilicity
Hydrophobicity or lipophilicity (ratio of octanol solubility to water solubility) measured through log P. The molecular lipophilic potential (MLP) was the first method designed to calculate the hydrophobic profile of a molecule in three dimensions [14, 15]. The development of the MLP was based on the finding that the partition coefficient (P) of a molecule, which represents its relative distribution over an octanol/water boundary, could be estimated from its chemical structure [16]. From the assumption that the log P is an additive property of the molecular fragments that make up a molecule, values for a wide variety of atom types and groups have been calculated. The receptor pocket may not be completely homogeneous (hydrophobic), so that log P does not perform a very good job for a large molecule with multiple position of substitution where as hydrophobicity of substituent (π) is more appropriate in such cases. For sake of simplicity each set has been studied separately as described below.
3.1.1. First Set
This set comprises of eighteen inhibitors (Table 1) [4]. The hydrophobicity of the each compound (log P), hydrophobicities of its substituents (πR, πX, πY, and πZ) and its interactions (Δπdr) with the substituents of receptor amino acids (πr) of this set have been evaluated and are presented in Table 4. A reference to the table shows that among the four substituents R, X, Y and Z, substituent-R has the highest value of πd and substituent-X has the lowest value. The decreasing order of substituent hydrophobicity has the sequence πR > πY > πZ > πX. Thus, substituent-R has maximum probability to interact hydrophobically with receptor amino acids: valine (Val), isoleucine (Ile), aspartic acid (Asp), threonine (Thr) and glycine (Gly). Val and Ile constituting the receptor amino acids of binding site, while Asp, Thr and Gly that of catalytic site of the enzyme, HIV-1-PR.The hydrophobic interactions (Δπdr) between substituent-R and receptor amino acid have been evaluated by solving Eq.1 and are also presented in Table 4. The decreasing order of Δπdr has the sequence Δπdr(Asp) > Δπdr(Thr) > Δπdr(Val) > Δπdr(Ile). A reference to these data shows that substituent-R has maximum interaction with the Ile amino acids as this interaction has lowest value of the Δπdr. The values of (ΣΔπdr)binding-site and (ΣΔπdr)catalytic-site (Table 10) shows that peptidic HIV-1-PRIs of this set interact with binding site rather than catalytic site as binding site have lower value of ΣΔπdr.Among the binding site, Ile has maximum interaction with the drug than Val, as it has lower value of Δπdr.
3.1.2. Second Set
This set comprises of seventeen inhibitors (Table 2) [5]. The log P, πR1, πR2, πR3, πr and Δπdr have been evaluated and are presented in Table 5. A reference to the table shows that among the three substituents R1, R2 and R3, substituent-R1 has the highest value of π. Thus, substituent-R1 has maximum probability to interact hydrophobically with receptor amino acids of the HIV-1-PR enzyme. The hydrophobic interactions (Δπdr) between substituent-R1 and receptor amino acid have been evaluated and are also presented in Table 5. The decreasing order of Δπdr has the order Δπdr(Asp) > Δπdr(Thr) > Δπdr(Val) > Δπdr(Ile). A reference to these data shows that substituent-R1 has maximum interaction with the Ile amino acids as this interaction has lowest value of the Δπdr. The values of (ΣΔπdr)binding-site and (ΣΔπdr)catalytic-site (Table 10) shows that inhibitors of this set interact with binding site rather than catalytic site as binding site have lower value of ΣΔπdr.Among the binding site, substituent-R1 has maximum interaction with the Ile than Val, as this interaction has lower value of Δπdr.
3.1.3. Third Set:
This set comprises of sixteen peptidic inhibitors (Table 3) [5]. The log P, πX and πr have been evaluated and are presented in Table 6. The hydrophobic interactions (Δπdr) between substituent-X and receptor amino acid have been evaluated and are also presented in Table 6. The decreasing order of Δπdr has the order Δπdr(Asp) > Δπdr(Thr) > Δπdr(Val) > Δπdr(Ile). A reference to these data shows that substituent-X has maximum interaction with the Ile amino acids as this interaction has lowest value of the Δπdr. The values of (ΣΔπdr)binding-site and (ΣΔπdr)catalytic-site (Table 10) shows that inhibitors of this set interact with binding site rather than catalytic site as binding site have lower value of ΣΔπdr.Among the binding site, substituent-X has maximum interaction with the Ile than Val, as this interaction has lower value of Δπdr.
The strength of the hydrophobic interaction is thus influenced not only by the polarity but also by the shape and extent of the exposed molecular surface [17]. It was also reported by Wang et al. that at least two additional factors are important in the binding of a compound to HIV-1-PR [16].The first one is the conformational inhibitor molecule, and the second one is the hydrophobic interaction between an inhibitor and the enzyme. As reported that HIV-1-PR has four hydrophilic pockets near its active sites, and it has also been shown that favorable hydrophobic interactions with these pockets are driblet for an inhibitor to achieve nanomolar potency [19]. The work of Wang et al. is also able to rationalize differences in binding affinity for enzyme-inhibitor complexes with largely hydrophobic binding sites.
3.2. Solvent accessible surface area
Solvent accessible surface area (SASA) is used to compute surface properties. SASA was first described by Lee and Richards in 1971 and is also known as Lee–Richard molecular surface [20]. Shrake and Ruplet in 1973 developed the rolling ball algorithm to calculate solvent accessible surface area [20]. SASA also provides a useful tool to gain insight into the overall extent of a hydrophobic region on a molecule or in the binding site of a protein, but unable for any real account of the particular atom types that make up the binding site or their positions relative to one another. Besides, it provides no means of assessing the shape of the binding site because it only calculates the relative accessibility of the contributing atoms [21]. Here, we studied hydrophobic interaction with respect to SASA and for the sake of simplicity each set has separately been described below.
3.2.1. First Set
The solvent accessible surface area (SASA) of each compound, solvent accessible surface area of its substituents (SASAR, SASAX, SASAY, and SASAZ) and its interactions (ΔSASAdr) with the substituents of receptor amino acids (SASAr) of the first set have been evaluated and are presented in Table 7. A reference to the table shows that among the four substituents R, X, Y and Z, substituent-R has the highest value of SASAd and substituent-X has the lowest value. The decreasing order of substituent hydrophobicity has the sequence: SASAR > SASAZ > SASAY > SASAX. Thus, substituent-R has maximum probability to interact with receptor amino acids. The hydrophobic interactions (ΔSASAdr) between substituent-R and receptor amino acid have been evaluated by solving Eq.2 and are also presented in Table 7. The decreasing order of ΔSASAdr has the sequence ΔSASAdr(Thr) > ΔSASAdr(Asp) > ΔSASAdr(Val) > ΔSASAdr(Ile). The values of (ΣΔSASAdr)binding-site and (ΣΔSASAdr)catalytic-site (Table 10) shows that peptidic HIV-1-PRIs interact with binding site rather than catalytic site as binding site have lower value of ΣΔSASA. A comparative study of ΔSASAdr reveals that Ile amino acid of binding sites has lowest values of ΔSASAdr than remaining amino acids. Thus, inhibitors of this set interact with Ile amino acid of binding site.
3.2.2. Second Set:
The SASA, SASAR1, SASAR2, SASAR3, SASAr and ΔSASAdr of each compound of second set have been evaluated and are presented in Table 8.
A reference to the table shows that among the three substituents R1, R2 and R3, substituent-R1 has the highest value of SASAd. Thus, substituent-R1 has maximum probability to interact hydrophobically with receptor amino acids of the HIV-1 protease enzyme. The hydrophobic interactions (ΔSASAdr) between substituent-R1 and receptor amino acid have been evaluated and are also presented in Table 8. The hydrophobic interactions in term of ΔSASAdr of substituent-R1 with respect to each amino acid have been evaluated and are also presented in Table-5. A reference to these data shows that among the five amino acids Ile has lowest value of the ΔSASAdr and Thr has the highest value of ΔSASAdr. The decreasing order of ΔSASAdr has the sequence ΔSASAdr(Thr) > ΔSASAdr(Val) > ΔSASAdr(Asp) > ΔSASAdr(Ile). The values of (ΣΔSASAdr)binding-site and (ΣΔSASAdr)catalytic-site (Table 10) shows that peptidic HIV-1-PRIs interact with binding site rather than catalytic site as binding site have lower value of ΣΔSASA. A comparative study of ΔSASAdr reveals that Ile amino acid of binding sites has lowest values of ΔSASAdr than remaining amino acids. Thus, peptidic HIV-1-PRIs of this set interact with Ile amino acid of binding site.
3.2.3. Third Set
The SASA, SASAX and SASAr of each compound of the third set have been evaluated and are presented in Table 9. The hydrophobic interactions interm of ΔSASAdr of substituent-X with respect to each amino acid have been evaluated and are also presented in Table 9. A reference to these data shows that among the five amino acids Ile has lowest value of the ΔSASAdr and Thr has the highest value of ΔSASAdr. The decreasing order of ΔSASAdr has the sequence: ΔSASAdr(Thr) > ΔSASAdr(Val) > ΔSASAdr(Asp) > ΔSASAdr(Ile). The values of (ΣΔSASAdr)binding-site and (ΣΔSASAdr)catalytic-site (Table 10) shows that peptidic HIV-1-PRIs interact with binding site rather than catalytic site as binding site have lower value of ΣΔSASA. A comparative study of ΔSASAdr reveals that Ile amino acid of binding sites has lowest values of ΔSASAdr than remaining amino acids. Thus, inhibitors of this set interact with Ile amino acid of binding site.
4. Conclusions
- Analysis of πd and SASAd shows that among the four substituents R, X, Y and Z of the first set, substituent-R has the highest value of both πd and SASAd. Among the three substituents R1, R2 and R3 of the second set, substituent-R1 also has the highest value of both πd and SASAd. Among the receptor amino acids (Val, Ile, Asp, Thr and Gly), Ile has the highest value of both πd and SASAd. Thus, πd, πr, SASASd and SASASr well 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 ΔSASAdr of all the three sets shows that they have the same order of hydrophobic interaction with respect to amino acids, Δπdr(Asp) > Δπdr(Thr) > Δπdr(Val) > Δπdr (Ile) and ΔSASAdr(Thr) > ΔSASAdr(Val) > ΔSASAdr(Asp) > ΔSASAdr(Ile). A comparative study of the values of (ΣΔπdr)binding-site, (ΣΔπdr)catalytic-site, (ΣΔSASAdr)binding-site, (ΣΔSASAdr)catalytic-site and shows that peptidic HIV-1-PRIs interact with binding site rather than catalytic site as binding site have lower value of ΣΔπdr and ΣΔSASAdr. Among the binding site, Val has maximum interaction than Ile, as it has lower vale of Δπdr and ΔSASAdr. Thus, Δπdr, ΔSASAdr, ΣΔπdr and ΣΔSASAdr well describe the extent of hydrophobic interaction of the drug with the receptor.
Acknowledgements
The authors gratefully acknowledge the financial support (Project No: C.S.T./D.3564 /11/2009) given by “Council of Science & Technology, U.P., INDIA”.
References
- Li, J.; Ma, X.; Guo, S.; Hou, C.; Shi, L.; Zhang, H.; Zheng, B.; Liao, C.; Yang, L.; Ye, L.; He, X. A Hydrophobic-Interaction-Based Mechanism Triggers Docking between the SARS-CoV-2 Spike and Angiotensin-Converting Enzyme 2. Global Challenges 2020, 4, 2000067 (1-8).[CrossRef] [PubMed]
- Cole, G. B.; Sharpe, S. Hydrophobic matching of HIV-1 Vpu transmembrane helix-helix interactions is optimized for subcellular location. Biochimica et Biophysica Acta (BBA) – Biomembranes, 2019, 1861, 183022 (1-8).[CrossRef] [PubMed]
- Sahu, V. K.; Khan, A. K. R.; Singh, R. K.; Singh, P. P. Drug-receptor interaction based quantitative structure-activity relationship of tetrahydroimidazodiazepinone.2009, 1243-1254.[CrossRef]
- Getman, D. P.;DeCrescenzo, G. A.;Heintz, R. M.;Reed, K. L.;Talley, J. J.;Bryant, M. L.;Clare, M.;Houseman, K. A.;Marr, J. J. Discovery of a novel class of potent HIV-1 protease inhibitors containing the (R)-(hydroxyethyl)urea isostere.J. Med. Chem. 1993, 36, 288–291.[CrossRef] [PubMed]
- Holloway,M. K.; Wai, J. M.;Halgren, T. A.;Fitzgerald, P. M. D.;Vacca, J. P.;Dorsey, B. D.;Levin, R. B.;Thompson, W. J.;Chen, L. J.; A priori prediction of activity for HIV-1 protease inhibitors employing energy minimization in the active site. J. Med. Chem. 1995, 38, 305–317.[CrossRef] [PubMed]
- Sahu, V. K.; Singh, R. K.; Singh, P. P. Extended Rule of Five and Prediction of biological activity of peptidic HIV-1-PR inhibitors. Universal Journal of Pharmacy and Pharmacology 2022, 1 (1), 20-42.[CrossRef]
- Stewart, J. J. P. Optimization of parameters for semiempirical methods I. Method. J. Comput. Chem. 1989, 10, 209-220.[CrossRef]
- Ghose, A. K.; Crippen, G. M. Atomic physicochemical parameters for 3-dimensional structure-directed quantitative structure-activity-relationships. 1. partition-coefficients as a measure of hydrophobicity. J. Comput. Chem. 1986, 7, 565-577.[CrossRef]
- Klamt, A.; Schuurmann, G. Cosmo-a new approach to dielectric screening in solvents with explicit expressions for the screening energy and its gradient. J. Chem. Soc.-Perkin Trans 1993, 2, 799-805.[CrossRef]
- Garg, R. Gupta, S. P.; Gao, H.; Babu, M. S.; Debnath, A. K.; Hansch, C. Comparative quantitative structure–activity relationship studies on anti-HIV drugs. Chem. Rev. 1999, 99, 3525-3601.[CrossRef] [PubMed]
- Pyring, D.; Lindberg, J.; Rosenquist, A.; Zuccarello, G.; Kvarnstrom, I.; Zhang, H.; Vrang, L.; Unge, T.; Classon, B.; Hallberg, A.; Samuelsson, B. Design and Synthesis of Potent C2-Symmetric Diol-Based HIV-1 Protease Inhibitors:Effects of Fluoro Substitution. J. Med. Chem. 2001, 44, 3083-3091.[CrossRef] [PubMed]
- Alterman, M.; Andersson, H. O.; Garg, N.; Ahlse´n, G.; Lo¨vgren, S.; Classon, B.; Danielson, U. H.; Kvarnstro¨m, I.; Vrang, L.; Unge, T.; Samuelsson, B.; Hallberg, A.Design and Fast Synthesis of C-Terminal Duplicated Potent C2-Symmetric P1/P1-Modified HIV-1 Protease Inhibitors. J. Med. Chem. 1999, 42, 3835-3844.[CrossRef] [PubMed]
- Ghosh, A. K.; Osswald, H. L.; Prato, G. Recent Progress in the Development of HIV-1 Protease Inhibitors for the Treatment of HIV/AIDS. J. Med. Chem. 2016, 59, 5172-5208.[CrossRef] [PubMed]
- Matthew, D. K.; Ricardo, L. M. A new method for estimating the importance of hydrophobic groups in the binding site of a protein. J. Med. Chem. 2005, 48, 1069-1078.[CrossRef] [PubMed]
- Audry, E.; Dubost, J. P.; Colleter, J. C.; Dallet, P. A new approach of structure-activity relationships-the potential of molecular lipophily. Eur. J. Med. Chem. 1986, 21, 71-72.
- Fujita, T.; Iwasa, J.; Hansch, C. New substituent constant pi derived from partition coefficients. J. Am. Chem. Soc. 1964, 86, 5175-5180.[CrossRef]
- Lee, S. H.; Rossky, P. J. A. Comparison of the structure and dynamics of liquid water at hydrophobic and hydrophilic surfacess A molecular dynamics simulation study. J. Chem. Phys. 1994, 100, 3334-3345.[CrossRef]
- Wang, S.; Milne, G. W. A.; Yan, X.; Posey, I. J.; Nicklaus, M. C.; Graham, L.; Rice, W. G. Discovery of novel, non-peptide HIV-1 protease inhibitors by pharmacophore searching. J. Med. Chem. 1996, 39, 2047–2054.[CrossRef] [PubMed]
- Ala, P.; Anton, E.D.; Bacheler, L.T.; Chang, C.H.; Garber, S.S.; Jadhav, P.K.; Woerner, F.J.Cyclic urea amides: HIV-1 protease inhibitors with low nanomolar potency against both wild type and protease inhibitor resistant mutants of HIV. J. Med. Chem. 1997, 40, 181-191.[CrossRef] [PubMed]
- Lee, B.; Richards, F. M. Interpretation of protein structures-estimation of static accessibility. J. Mol. Biol. 1971, 55, 379-400.[CrossRef]
- Shrake, A.; Rupley, J. A. Environment and exposure to solvent of protein atoms-lysozyme and insulin. J. Mol. Biol. 1973, 79, 351-371.[CrossRef]
Copyright
© 2025 by authors and Scientific Publications. This is an open access article and the related PDF distributed under the Creative Commons Attribution License which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Article Metrics
If you find this article cited by other articles, please click the button to add a citation.