Drug-Receptor Interaction of Peptidic HIV-1 Protease: Polar Effect-II
Abstract
Klopman described the chemical reaction of metal ions and base ions in term of softness, E‡n and E‡m, respectively. By simple modification of known methods, Singh et al. made it applicable for neutral Lewis acids (transition metal salts) and bases (organic molecules) and also extended its application to biological systems for site selectivity and to explain reaction mechanism (markovnikov and anti-markovnikov rule), ligand-receptor interaction of testosterones, estrogens and tetrahydroimidazobenzodiazepinone. In this study effective atomic softness E‡n(eff) and E‡m(eff), and their change ΔE‡nm have been used for site selectivity and polar interaction between 51 peptidic HIV-1 protease inhibitors and receptor amino acids. ΔE‡nm values derived from drug-receptor interaction show that when one moiety on receptor behaves as nucleophile (O of valine amino acid) at the same time maximum electrophilic site of the drug (C-atom of the maximum E‡n(eff) value) orient itself to come close the respective site and make maximum interaction, while when another moiety on receptor behaves as electrophilic site (C of isoleucine amino acid), at the same time maximum nucleophilic site of the drug (O-atom of the maximum E‡m(eff) value) also orient itself to come close the respective site and make maximum interaction.
1. Introduction
In our successive publications, we have studied pharmacokinetics followed by hydrophobic interaction of peptidic HIV-1 protease inhibitors [1, 2]. The present work describes the polar interaction based on effective atomic softness. Klopman provides a very convenient way to describe the chemical reactivity of a compound with the help of atomic softness values in terms of E‡n and E‡m [3]. This concept was based on the charge and frontier orbital controlled chemical reaction of perturbation theory.
where, E‡n = softness of an acid, E‡m = softness of a base, IP = ionization potential of atom, EA = electron affinity of atom, a2 = ¾, b2 = ¼, C = 1, ϵ = dielectric constant of the medium in which reaction is carried out, R = radius of atom whose softness is to be determined, q = charge on the atom, C = electron density, χ = q –(q – 1)√k, and k = 0.75, a and b are variational parameters. By a simple modification of known methods [4, 5, 6, 7], Singh et al. were calculated the values of IP, q and R and made it applicable for a neutral chemical system [8]. The softness values so derived by them are termed as “effective softness” and are designated by symbols E‡n(eff) for Lewis acid and E‡m(eff) for Lewis base [9, 10, 11]. In this study semiempirical method has been used to calculate effective atomic softness: E‡n(eff) and E‡m(eff), for 51 peptidic HIV-1 protease inhibitors (HIV-1-PRIs).
2. Materials and Methods
Fifty-one HIV-1-PRIs have been used as study material and are separately listed in Table 1, 2 and 3 with their observed biological activities in term of IC50 [12, 13, 14]. All the fifty-one inhibitors have been divided in three sets on the basis of their structural similarities (Figure 1, 2and 3). The first, second and third set comprises of eighteen, seventeen and sixteen inhibitors, respectively. The logarithms of the inverse of IC50 have been used as biological end point (log1/C) in the study. For solving the modified Klopman equations, the 3D modeling and geometry optimization of all the compounds have been performed with the help of CAChe Pro software of Fujitsu [15, 16]. The study is based on semiempirical PM3 method [17, 18]. The PM3 based calculations have been performed with MOPAC 2002 software associated with CAChe. The values of various parameters to solve modified Klopman softness have been calculated by softness calculator which was developed by Singh research group [4]. Singh et al. made Klopman equations applicable for neutral Lewis acids (transition metal salts) and bases (organic molecules) and also extended its application to biological systems for site selectivity and to explain reaction mechanism (markovnikov and anti-markovnikov rule), ligand-receptor interaction of testosterones, estrogens and tetrahydroimidazobenzodiazepinone [19, 20, 21, 22, 23].
3. Results and Discussion
Out of fifty-one compounds under study, the eighteen compounds (compound no. 1-18) have the parent skeleton of Figure 1, which has 25 sites. Out of remaining thirty-three, the seventeen compounds (compound no. 19-35) have the parent skeleton of Figure 2, which has 33 sites. While the remaining sixteen compounds (compound no. 36-51) have the parent skeleton of Figure 3, which has 29 sites [10, 11]. The effective softness values represented by E‡n(eff) describe the electrophilic character of compound, whereas E‡m(eff) describe the nucleophilic character of compound [19, 20, 21, 22, 23]. E‡n(eff), IP and EA of electron deficient carbon center (electrophilic site) have been evaluated at sixteen sites (C2, C3, C5, C6, C7, C9, C11, C12, C16, C17, C18, C19, C20, C21, C22 and C24) of compound no 1-18, twenty six sites (C1, C3, C5, C6, C7, C8, C9, C11, C12, C13, C14, C15, C16, C17, C18, C19, C20, C21, C22, C24, C25, C26, C27, C28, C29 and C30) of compound no 19-35 and twenty three sites (C1, C3, C5, C6, C7, C8, C9, C10, C11, C12, C14, C15, C16, C17, C18, C19, C20, C22, C24, C25, C26, C27 and C28) of compound no 36-51. While E‡m(eff), IP and EA of electron rich center (nucleophilic site) have been evaluated at nine sites (N1, N4, N8, N10, O13, N14, O15, O23 and O25) of compound no 1-18, seven sites (O2, N4, N10, O23, O31, O32 and O33) of compound no 19-35 and six sites (O2, N4, O13, O21, O23 and O29) of compound no 36-51. Softness parameter is a very dominating factor when correlated with the mechanism of action of a variety of known therapeutic agents and their pharmacoactivities. Because it included atomic radius (R), electron density (C), charge (q), ionization potential (IP), and electron affinity (EA) of atom. The incorporation of dielectric constant (ε) shows the effect of medium in which reaction is carried out. Correlation of this dominating factor, with mechanism of action of various therapeutic agents and their pharmacoactivities, would be very valuable in search of a new advance drug before its synthesis. This technique when applied saves time and resources with the limited facilities for a medicinal scientist. Effective softness values represented by E‡n(eff) describe the electrophilic character of an atom within a molecule. Higher the value of E‡n(eff) of a site greater will be the electrophilic character of that site within the molecule. Further, the site of highest E‡n(eff) value characterizes the susceptibility of the molecules toward the attack of nucleophile. The highest E‡n(eff) values of compounds are placed in Table 4. A reference to this table shows that, in case of compounds of group-A, the E‡n(eff) value of all the compounds is highest at position C11, except compounds 5, 6, 10 and 12, which has highest E‡n(eff) value at position C7. A close look at parent skeleton of compounds of group-A clearly indicates that the positions C7 and C11 are carbon to which the urea isostere and -CONH2 group is linked to it, respectively. In group-B, the E‡n(eff) value of compounds 19, 22, 23, 27, 31, 33, 34 and 35 is highest at position C21. While of compounds 20, 21, 25, 26, 28, 29 and 30, it is highest at position C22 and in compounds 24 and 32 it is highest at position C7. A close look at parent skeleton of compounds of group-B clearly indicates that the position C7, C21 and C22 is carbon of methylene, methyl and methyl moiety, respectively. In group-C, the E‡n(eff) value of compounds 37, 38, 40-45, 47-49 and 51 is highest at position C11. While of compounds 36 and 39, it is highest at position C12 and in compounds 46 and 50 it is highest at position C7. A close look at parent skeleton of compounds of group-C clearly indicates that the position C7, C11 and C12 is also carbon of methylene, methyl and methyl moiety, respectively. The examination of Table 1 indicates that there must be a relationship between effective softness, E‡n(eff), and observed biological activity. Of course, there exist relationship between E‡n and log 1/C but there no sequential rise or fall. In order to provide sequential relationship, we have divided the compounds into subgroups, group-A: subgroup-a, b, c and d; group-B: subgroup-e, f, g and h; and group-C: subgroup-i, j and k. A close look to these subgroups has shown that (i) In subgroup-a, compound 1, 2 and 16; in subgroup-b, compound 3, 5 and 17; in subgroup-c, compound 6, 7, 8, 9, 11 and 13; and in subgroup-d, compound 10, 14 and 18 show the direct relationship very clearly. Compound 4, 12 and 15 do not follow the sequential trend, (ii) In subgroup-e, compound 20, 24, 30, 31 and 32 and in subgroup-f, compound 22, 25, 26 and 27 show the direct relationship very clearly. While in subgroup-g, compound 21, 23 and 35 and in subgroup-h, compound 19, 33 and 34 shows the inverse relationship very clearly. Compound 28 and 29 do not follow the sequential trend, and (iii) In subgroup-i, compound 36, 38, 42, 43 and 50; and in subgroup-j, compound 40, 47 and 49 show the direct relationship very clearly. While in subgroup-l, compound 39, 41, 45, 46 and 48 shows the inverse relationship very clearly. Compound 28, 29, 37 and 44 do not follow the sequential trend.
Effective softness values represented by E‡m(eff) describe the nucleophilic character of an atom within a molecule. Higher the value of E‡m(eff) with negative sign of a site greater will be the nucleophilic character of that site within the molecule. Further, the site of highest E‡m(eff) value characterizes the susceptibility of the molecules toward the attack of electrophile. The highest E‡m(eff) values of compounds are also placed in Table 4. A reference to this table also shows that, in group-A, the highest E‡m(eff) value in all the eighteen compounds is associated with position O13. A close look at parent skeleton of compounds of group-A clearly indicates that the position O13 is oxygen of amidic moiety. In group-B, E‡m(eff) value of compounds 19-24, 28, 31, 33 and 34, it is highest at position O23. While of compounds 25-27, 29, 30, 32 and 35 is highest at position O32. A close look at parent skeleton of compounds of group-B clearly indicates that the positions O23 and O32 are oxygen of carboxylic and amidic moiety, respectively. In group-C, the E‡m(eff) value of compounds 36 and 37 it is highest at position O21. While of compounds 39-46, it is highest at position O13 and in compounds 38 and 47-51 it is highest at position C29. A close look at parent skeleton of compounds of group-C clearly indicates that the positions O13, O21 and O29 are oxygen of carboxylic, hydroxyl and carbonyl moiety, respectively. The examination of Table 1 also indicates that there must be a relationship between effective softness, E‡m(eff), and observed biological activity. Of course, there also exist relationship between E‡m(eff) and log 1/C, but there is also no sequential rise or fall. In order to provide sequential relationship, we have divided the compounds into same three groups, group-A: compound no. 1-18; group-B: compound no. 19-35 and group-C: compound no. 36-51, on the basis of their parent skeleton. But, each group has subdivided into different subgroups, group-A: subgroup-a, b, c and d; group-B: subgroup-e, f and g; and group-C: subgroup-h, i and j. A close look to these subgroups has shown that (i) In subgroup-a, compound 2, 8, 15 and 16; in subgroup-b, compound 5, 6, 13, 14, 17 and 18 show the direct relationship very clearly. In subgroup-c, compound 1, 3, 4, 7 and 11; and in subgroup-d, compound 9, 10 and 12 show the inverse relationship very clearly, (ii) In subgroup-e, compound 20, 28 and 30; in subgroup-f, compound 22, 27 and 34 and in subgroup-g, compound 19, 21, 23, 24, 29, 32, 33 and 35 show the inverse relationship very clearly. Compound 25, 26 and 31 do not follow the sequential trend, and (iii) In subgroup-h, compound 39, 40, 42 and 47; and in subgroup-j, compound 37, 38, 44, 48 and 50 show the direct relationship very clearly. While in subgroup-i, compound 41, 43, 45, 46 and 49 show the inverse relationship Compound 36 and 51 do not follow the sequential trend.
It is well established that the stability of the compound formed between nucleophile and electrophile depends upon the value of difference between softness values of E‡m(eff) of nucleophile, and softness values of E‡n(eff) of electrophile, ΔE‡nm represent the difference. The higher is the ΔE‡nm greater is the stability of the compound [24].
HIV-1 protease enzyme (HIV-1-PR) is a viral encoded homodimeric aspartyl protease with C2 symmetry [25, 26]. A catalytic triad of Asp-Thr-Gly contributed by each monomer comprises the active site of the enzyme. The amino acids constituting the binding site (S2/S2’ pocket) are Val-32, Ile-47, Ile-50, and Ile-84 in each monomeric polypeptidic unit of the protease enzyme [27, 28]. HIV inhibitors make initial contact with the receptor amino acids Val and Ile of binding site rather than Asp, Thr and Gly of catalytic site. The study has been made on interaction with Val and Ile of binding site. The nitrogen of amino and oxygen atom of carboxylic group of receptor protein have the nucleophilic character and can interact with the electrophilic center on the inhibitor, while the electron deficient carbon atoms can interact with the nucleophilic site on the inhibitor and vice versa. Thus, the highest E‡n(eff) value of compounds characterize the susceptibility of the molecules toward the attack of nucleophile, while the highest E‡m(eff) value of compounds characterize the susceptibility of the molecules toward the attack of electrophile. The softness values E‡n(eff) of at 5 sites of valine and 6 sites of isoleucine, while the softness values E‡m(eff) of 2 sites of both Val and isoleucine have also been evaluated and are presented in Table 5. When inhibitors are treated as acids and receptor proteins as base, the highest values of softness E‡n(eff) of inhibitors and highest values of E‡m(eff) of amino acids have been used for deriving ΔE‡nm values. While, when inhibitors are treated as bases and receptor proteins as acids, the highest values of softness E‡m(eff) of inhibitors and highest values of E‡n(eff) of amino acids have been used for deriving ΔE‡nm values. The highest values of softness and the ΔE‡nm derived from them are given in Table 6. A reference to this table 1 shows that in the former case the ΔE‡nm values in case of Val amino acid is higher than the Ile amino acid, while in later case the ΔE‡nm values in case of Ile amino acid is higher than the Val amino acid.
4. Conclusion
Softness parameter is a very dominating factor when correlated with the mechanism of action of a variety of known therapeutic agents and their pharmacoactivities. Because it included atomic radius (R), electron density (C), charge (q), ionization potential (IP), and electron affinity (EA) of atom. ΔE‡nm values derived from drug-receptor interaction show that when one moiety on receptor behaves as nucleophile (O of valine amino acid) at the same time maximum electrophilic site of the drug (C-atom of the maximum E‡n(eff) value) orient itself to come close the respective site and make maximum interaction, while when another moiety on receptor behaves as electrophilic site (C of isoleucine amino acid), at the same time maximum nucleophilic site of the drug (O-atom of the maximum E‡m(eff) value) also orient itself to come close the respective site and make maximum interaction.
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
- 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, 20-42.[CrossRef]
- Sahu, V. K.; Singh, R. K.; Singh, P. P. Drug-Receptor Interaction of Peptidic HIV-1 Protease: The Hydrophobic Effect-I. Online Journal of Microbiological Research. 2022, 1, 33-48.
- Klopman, G. Chemical Reactivity and The Concept of Charge and Frontier Controlled Reactions. J. Am. Chem. Soc. 1968, 90, 223-234.[CrossRef]
- Dewar, M. J. S.; Mortia, T. F. Ground States of Conjugated Molecules. XII. Improved Calculations for Coumpounds Containing Nitrogen and Oxygen. J. Am. Chem. Soc. 1968, 91, 796-802[CrossRef]
- Sanderson, R. T. Vhemical Bond and Bond Energy. Academic Press, New York, 1971, pp19-20 .
- Huheey, J. E. Inorganic Chemistry. Harper & Row, New York, 1972, pp71-72
- Pauling, L. The Nature of Chemical Bond, 3rd Edn. Cornell University Press, New York, 1960, pp513.
- Singh, P. P.; Srivastava, S. K.; Sivastava, A. K. A Matching Between Lewis Acids and Lewis Bases on The Basis of Quantitative Softness Values and Their Relation With The Stability of The Complexes-I. J. Inorg. Nucl. Chem. 1980, 42, 521-532.[CrossRef]
- Singh, P. P. Quantitative Softness Values and Their Novel Applications. J. Scient. Ind. Res. 1983, 42, 140-148.
- Singh, P. P; Naqvi, M. I.; Singh, N. B. New softness parameters and Hammett’s equation. Indian J. Chem. 1992, 33A, 586-589
- Singh, P. P.; Pasha, F. A.; Sivastava, H. K. Novel application of softness parameter for regioselectivity and reaction mechanism. Indian J. Chem. 2004, 43B, 983-991.
- 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]
- 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-3602.[CrossRef] [PubMed]
- Thiel, W. Fast Semiempirical Geometry Optimizations. J. Mol. Struct. 1988, 163, 415–429.[CrossRef]
- Komornicki A.; McIver, J. W. Rapid Geometry Optimization for Semiempirical Molecular Orbital Methods. Chem. Phys. Lett, 1971, 10, 303-306.[CrossRef]
- Stewart, J. J. P. Optimization of parameters for semiempirical methods I. Method. J. Comput. Chem. 1989, 10, 209-220.[CrossRef]
- Stewart, J. J. P. Optimization of parameters for semiempirical methods IV: extension of MNDO, AM1, and PM3 to more main group elements. J. Mol. Model. 2004, 10, 155-164.[CrossRef] [PubMed]
- Singh, P. P.; Pasha, F. A.; Sivastava, H. K. DFT Based Atomic Softness and Its Application in Site Selectivity. QSAR Comb. Sci. 2003, 22, 843-851.[CrossRef]
- Sivastava, H. K.; Pasha, F. A.; Singh, P. P. Atomic Softness-Based QSAR Study of Testosterone. Int. J. Quantum Chem. 2005, 103, 237-245 .[CrossRef]
- Pasha, F. A.; Sivastava, H. K.; Singh, P. P. Semiempirical QSAR study and ligand receptor interaction of estrogen. Mol. Diver. 2005, 9, 215-220 .[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. Int. J. Quantum Chem. 2009, 109, 1243-1254.[CrossRef]
- Sivastava, H. K.; Pasha, F. A.; Mishra, S. K.; Singh, P. P. Noval applications of atomic softness and QSAR study of testosterone derivatives. Med. Chem. Res. 2009, 103, 237-245.
- Sahu, V. K.; Khan, G.; Verma, R. N.; Singh, P. P. Complexes of Cobaltocene: A Effective Atomic Softness and Fukui Function Based Study. RJPBCS 2010, 1, 535-544.
- 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]
- Szeltner, Z.; Polga, R. L. Conformational Stability and Catalytic Activity of HIV-1 Protease Are Both Enhanced at High Salt Concentration. J. Biol. Chem. 1996, 271, 5458-5463.[CrossRef] [PubMed]
- Srivastava, H. K.; Choudhury, C.; Sastry, G. N. The Efficacy of Conceptual DFT Descriptors and Docking Scores on the QSAR Models of HIV Protease Inhibitors. Medicinal Chemistry 2012, 8, 811-825.[CrossRef] [PubMed]
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.