In silico Studies on N- (Pyridin-2-yl) Thiobenzamides as NNRTIs against Wild and Mutant HIV-1 Strains
Anuradha Singh, Vishal Kumar Singh, Rajesh Verma, and Ramendra K. Singh*
Bioorganic Research Laboratory, Department of Chemistry,
University of Allahabad, Allahabad-211002, India
*Corresponding author: This email address is being protected from spambots. You need JavaScript enabled to view it.
ABSTRACT
In the present study, keeping the Lipinski’s Rule of Five in focus, a series of new 4-(4-benzenesulfonylamino)-N-(5-substituted-pyridin-2-yl)-thiobenzamides bearing different substituents at the C-4 position of benzenesulfonylamino ring have been designed as NNRTIs of wild-type (WT) and mutant HIV-1 strains. Molecules having drug-like character were further docked into the active domain of wild-type (WT) RT/1 with entry code (PDB: ID 3mec) and K103N/TYR181 mutant RT/2 complex (PDB: ID 3BGR) using Discovery Studio 2.5 software. Analysis of the docking results revealed that all molecules formed hydrogen bonds with active amino acids (Lys101, Lys103, Tyr181, and Tyr318) and exhibited π-stacking interactions with Tyr181, Tyr188, Phe227, and Trp229 present in the NNIBP with both WT and mutant HIV-1 RT. The designed ligands adopted ‘horseshoe/seahorse’ conformation inside the NNIBP like other second generation NNRTIs and formed more stable complexes (total interaction energy found in the range of (-) 54 – (-) 77 kcal/mol) with HIV-1 RT in comparison to etravirine and rilpivirine (-)61.43 and (-)50.23 Kcal/mol, respectively. Consequently, lower EC50 values were predicted for N- (Pyridin-2-yl) derivatives. Structure-activity relationships (SARs) are discussed in terms of a possible interaction with the RT binding site, depending on the nature of the substituent at ring A and ring C.
INTRODUCTION
Kappaphycus alvarezii (Doty ex P.C. Silva 1996) is an introduced-species, which has been mono-cultivated as the largest seaweed commodity not only in Indonesia, but also in other tropical Asian countries. Indonesia is currently leading the production of K. alvarezii with nearly 1.5 x 106 tons in 2009, driven by the increasing demand for κ-carrageenan (Meinita et al. 2012). Another potential seaweed is Padina australis (Hauck 1887), a native species of brown seaweed that is increasingly being cultivated and studied as source of alginate (Widyastuti 2012). In addition to carrageenan and alginate, K. alvarezii and P. australis are also rich sources of natural pigments that have other additional biological functions (i.e., antioxidant, anti-obesity, anti-inflammation, and anti-hyper cholesterol) as well as natural colorants, depending on the chemical properties of the pigment. However, while great interest has been expressed in recent decade to characterize the distribution, biosynthesis, and functions as well as application of the photosynthetic pigments of microalgae, little effort has been dedicated to characterize other pigments in those seaweeds......read more
REFERENCES
ABRAM ME, FERRIS AL, DAS K, QUINOÑES O, SHAO W, TUSKE S, ALVORD WG, ARNOLD E, HUGHES SH. 2014. Mutations in HIV-1 Reverse Transcriptase Affect the Errors Made in a Single Cycle of Viral Replication. J Virol 88(13):7589-7601.
BARRECA ML, RAO A, LUCA LD, IRACI N, MONFORTE AM, MAGA G, DE CLERCQ E, PANNECOUQUE C, BALZARINI J, CHIMIRRI A. 2007. Discovery of novel benzimidazoles as potent non-nucleoside reverse transcriptase inhibitors active against wild-type and mutant HIV-1 strains. Bioorg Med Chem Lett 17:1956–-1960.
BINDEWALD E, SKOLNICK J. 2005. A Scoring Function for Docking Ligands to Low-Resolution Protein Structures. J Comput Chem 26:374-383.
BÖHM HJ. 1994. The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure. J Comput Aided Mol Des 8:243-256.
BÖHM HJ. 1998. Prediction of binding constants of protein ligands: a fast method for the prioritization of hits obtained from de novo design or 3D database search programs. J Comput Aided Mol Des12:309-32.
BROOKS BR, BRUCCOLERI RE, OLAFSON BD, STATES DJ,
SWAMINATHAN S, KARPLUS M. 1983. CHARMM: a program for macromolecular energy, minimization, and dynamics calculations. J Comput Chem 4:187-217.
ChemDraw [Computer software]. Available at http://www.cambridgesoft.com/software/overview.aspx
DAS K, LEWI PJ, HUGHES SH, ARNOLD E. 2005. Crystallography
and the design of anti-AIDS drugs: conformational flexibility and positional adaptability are important in the design of non-nucleoside HIV-1 reverse transcriptase inhibitors. Prog Biophys Mol Biol 88:209-31.
DAS K, BAUMAN JD, CLARK ADJ, FRENKEL YV, LEWI PJ, SHATKIN AJ, HUGHES SH, ARNOLD E. 2008. High-resolution structures of HIV-1 reverse transcriptase/TMC-278 complexes: strategic flexibility explains potency against resistance mutations. Proc. Natl Acad Sci USA. 105:1466-1471.
DE CLERCQ E. 2004. Non-nucleoside reverse transcriptase inhibitors (NNRTIs): past, present, and future. Chem Biodiverse 1:44-64.
Discovery Studio (Version 2.5) [Computer software]. 25 Jul 2009. Available at http://accelrys.com
ESPOSITO F, CORONA A, TRAMONTANO E. 2012. HIV-1 Reverse Transcriptase Still Remains a New Drug Target: Structure, Function, Classical Inhibitors, and New Inhibitors with Innovative Mechanisms of Actions. Mol Biol Int p. 1-23.
GALLIVAN JP , DOUGHERTY DA. 1999. Cation-pi interactions in structural biology. Proc Natl Acad Sci USA 96:9459-9464.
GEHLHAARL DK, VERKHIVKER GM, REJTO PA, SHERMAN CJ, FOGEL DB, FOGEL LJ, FREER ST. 1995. Molecular recognition of the inhibitor AC1343 by HIV-l protease: conformationally flexible docking by evolutionary programming. Chemistry & Biology 2:317-324.
JAIN AN 2006. Scoring functions for protein-ligand docking. Current Protein & Peptide Science 7(5):407-420.
JONCKHEERE H, ANNE J, DE CLERCQ E. 2000. The HIV-1 reverse transcriptase (RT) process as target for RT inhibitors. Med Res Rev 20:129-154.
KUMARI G, SINGH RK. 2012. Highly Active Antiretroviral Therapy for treatment of HIV/AIDS patients: Current Status and Future Prospects and the Indian Scenario. HIV&AIDS Review 11:5-14.
KUMARI G, SINGH RK. 2013. Anti-HIV drug development: Structural features and limitations of Present day drugs and future challenges in the successful HIV/AIDS treatment. Curr. Pharm. Des. 19:1767-83.
LA REGINA G, COLUCCIA A, SILVESTRI R. 2010. Looking for an active conformation of the future HIV type-1 non-nucleoside reverse transcriptase inhibitors. Antivir Chem Chemother. 20:213-37.
LAGOS CF, CABALLERO J, GONZALEZ-NILO FD, DAVID PESSOA-MAHANA C, PEREZ-ACLE T. 2008. Docking and quantitative structure-activity relationship studies for the Bisphenylbenzimidazole Family of Non-Nucleoside Inhibitors of HIV-1 Reverse Transcriptase. Chem Biol Drug Des 72:360-369.
LANSDON EB, BRENDZA KM, HUNG M, WANG R, MUKUND S, JIN D, BIRKUS G, KUTTY N, LIU X. 2010. Crystal Structures of HIV-1 Reverse Transcriptase with Etravirine (TMC125) and Rilpivirine (TMC278): Implications for Drug Design. J Med Chem 53:4295-99.
LIPINSKI CA, LOMBARDO F, DOMINY BW, FEENEY PJ. 2001.Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 46:3-26.
MARTINS S, RAMOS MJ, FERNANDES PA. 2008. The current status of the NNRTI family of Antiretroviral used in the HAART regime against HIV infection. Curr Med Chem 15:1083-95.
MAYO SL, OLAFSON BD, GODDARD WA. 1990. Dreiding: a generic force forcefield for molecular simulation. J Phys Chem 94:8897-8909.
Molinspiration Cheminformatics. Available at http://www.molinspiration.com
MOMANY FA, RONE R. 1992. Validation of the general purpose QUANTA _3.2 ⁄ CHARMm_ force field. J Comput Chem 13:888-990.
MUEGGE I, MARTIN YC. 1999. A general and fast scoring function for protein-ligand interactions: a simplified potential approach. J Med Chem 42:791-804.
PENG Z, XUWANG C, DONGYUE L, ZENGJUN F, DE CLERCQ E, XINYONG L. 2013. HIV-1 NNRTIs: structural diversity, pharmacophore similarity, and implications for drug design. Med Res Rev 33:1-72.
RCSB Protein Data Bank. Retrieved from http://www.rcsb.org
SINGH RK, YADAV D, RAI D, KUMARI G, PANNECOUQUE C. 2010. Synthesis, structure- activity relationship and anti viral activity of 3'-N, N-dimethylamino-2', 3'-dideoxythymidine and its prodrugs. Eur J Med Chem 45:3787-93.
SINGH A, YADAV D, YADAV M, DHAMANAGE A, KULKARNI S, SINGH RK. 2015. Molecular modeling, synthesis and biological evaluation of N-heteroaryl compounds as reverse transcriptase inhibitors against HIV-1. Chem Biol Drug Des 85:336-347.
SINGH A, YADAV M, SRIVASTAVA R, SINGH N, KAUR R, GUPTA S K, SINGH RK. 2016. Design and anti-HIV activity of aryl sulphonamide as non-nucleoside reverse transcriptase inhibitors. Med Chem Res 25:2842-59.
[UNAIDS] Joint United Nations Programme on HIV and AIDS. Global AIDS Update 2016. Retrieved from http://www.unaids.org/en/resources/ on 20 Mar 2017.
VENKATACHALAM CM, JIANG X, OLDFIELD T, WALDMAN M. 2003. LigandFit: a novel method for the shape-directed rapid docking of ligands to protein active sites. J Mol Graphics Modell. 21:289-307.
YADAV M, SINGH A, SINGH RK. 2017. Docking Studies on Novel Bisphenylbenzimidazoles (BPBIs) as Non-nucleoside Inhibitors of HIV-1 Reverse Transcriptase. IJCB 56B(7):714-723
ZHAO X, LIU X, WANG Y, CHEN Z, KANG L, ZHANG H, LUO X, ZHU W, CHEN K, LI H, WANG X, JIANG H. 2008. An improved PMF scoring function for universally predicting the interactions of a ligand with protein, DNA, and RNA. J Chem Inf Model 48(7):1438-47.
ZHOU Z, LIN X, MADURA JD. 2006. HIV-1 RT nonnucleoside inhibitors and their interaction with RT for antiviral drug development. Infect. Disord Drug Targets 6:391-413.