Allosteric mechanism of quinoline inhibitors for HIV RT‐associated RNase with MD simulation and dynamics fluctuation network

The human immunodeficiency virus (HIV) is a retrovirus which infects T lymphocyte of human body and causes immunodeficiency. Reverse transcriptase inhibitors (RTIs) can inhibit some functions of RT, preventing virus synthesis (double‐stranded DNA), so that HIV virus replication can be reduced. Experimental results indicate a series of benzimidazole‐based inhibitors which target HIV RT‐associated RNase to inhibit the reverse transcription of HIV virus. However, the allosteric mechanism is still unclear. Here, molecular dynamics simulations and dynamics fluctuation network analysis were used to reveal the binding mode between the inhibitors and RT‐associated RNase. The most active molecule has more hydrophobic and electrostatic interactions than the less active inhibitor. Dynamics correlation network analysis indicates that the most active inhibitor perturbs the network of RT‐associated RNase and decreases the correlation of nodes. 3D‐QSAR model suggests that two robust and reliable models were constructed and validated by independent test set. 3D‐QSAR model also shows that bulky negatively charged or hydrophilic substituent is favorable to bioactivity. These results reveal the allosteric mechanism of quinoline inhibitors and help to improve the bioactivity.

[1]  J. Chermann,et al.  Isolation of a T-lymphotropic retrovirus from a patient at risk for acquired immune deficiency syndrome (AIDS). , 1983, Science.

[2]  G. Klebe,et al.  Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. , 1994, Journal of medicinal chemistry.

[3]  R Brookmeyer,et al.  Identification of a reservoir for HIV-1 in patients on highly active antiretroviral therapy. , 1997, Science.

[4]  Hao Liu,et al.  Synergistic Allosteric Mechanism of Fructose-1, 6-bisphosphate and Serine for Pyruvate Kinase M2 via Dynamics Fluctuation Network Analysis , 2016, J. Chem. Inf. Model..

[5]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[6]  Dong Song,et al.  The IDP-Specific Force Field ff14IDPSFF Improves the Conformer Sampling of Intrinsically Disordered Proteins , 2017, J. Chem. Inf. Model..

[7]  Ray Luo,et al.  Synergistic Modification Induced Specific Recognition between Histone and TRIM24 via Fluctuation Correlation Network Analysis , 2016, Scientific Reports.

[8]  C. Tzeng,et al.  Synthesis and antibacterial evaluation of certain quinolone derivatives. , 2001, Journal of medicinal chemistry.

[9]  J. Moore,et al.  HIV-1 entry inhibitors: Evading the issue , 1999, Nature Medicine.

[10]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[11]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[12]  M. Karelson,et al.  Quantum-Chemical Descriptors in QSAR/QSPR Studies. , 1996, Chemical reviews.

[13]  Garland R. Marshall,et al.  VALIDATE: A New Method for the Receptor-Based Prediction of Binding Affinities of Novel Ligands , 1996 .

[14]  S. Free,et al.  A MATHEMATICAL CONTRIBUTION TO STRUCTURE-ACTIVITY STUDIES. , 1964, Journal of medicinal chemistry.

[15]  K A Johnson,et al.  Mechanism of inhibition of HIV-1 reverse transcriptase by nonnucleoside inhibitors , 1995, Science.

[16]  R. Philen,et al.  Synthsis and antitumor properties of bis(quinaldine) derivatives. , 1977, Journal of medicinal chemistry.

[17]  R. Jain,et al.  Substituted 4-methylquinolines as a new class of anti-tuberculosis agents. , 2003, Bioorganic & medicinal chemistry letters.

[18]  Hans-Joachim Böhm,et al.  The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[19]  M.G.B. Drew,et al.  The art of molecular dynamics simulation , 1996 .

[20]  Nouri Neamati,et al.  Application of CoMFA and CoMSIA 3D-QSAR and docking studies in optimization of mercaptobenzenesulfonamides as HIV-1 integrase inhibitors. , 2004, Journal of medicinal chemistry.

[21]  S. Nagamani,et al.  Pharmacophore filtering and 3D-QSAR in the discovery of new JAK2 inhibitors. , 2011, Journal of molecular graphics & modelling.

[22]  C. Lipinski Lead- and drug-like compounds: the rule-of-five revolution. , 2004, Drug discovery today. Technologies.

[23]  Asim Kumar Debnath,et al.  Application of 3D-QSAR techniques in anti-HIV-1 drug design--an overview. , 2005, Current pharmaceutical design.

[24]  V. Buchta,et al.  Quinaldine derivatives: preparation and biological activity. , 2005, Medicinal chemistry (Shariqah (United Arab Emirates)).

[25]  Y. Pommier,et al.  Integrase inhibitors to treat HIV/Aids , 2005, Nature Reviews Drug Discovery.

[26]  Tomozumi Imamichi,et al.  Action of anti-HIV drugs and resistance: reverse transcriptase inhibitors and protease inhibitors. , 2004, Current pharmaceutical design.

[27]  P A Kollman,et al.  Free energy calculations on dimer stability of the HIV protease using molecular dynamics and a continuum solvent model. , 2000, Journal of molecular biology.

[28]  Xin Xu,et al.  Structure-based rational quest for potential novel inhibitors of human HMG-CoA reductase by combining CoMFA 3D QSAR modeling and virtual screening. , 2007, Journal of combinatorial chemistry.

[29]  J. Dearden,et al.  QSAR modeling: where have you been? Where are you going to? , 2014, Journal of medicinal chemistry.

[30]  W F van Gunsteren,et al.  On the interpretation of biochemical data by molecular dynamics computer simulation. , 1992, European journal of biochemistry.

[31]  C. Hansch,et al.  p-σ-π Analysis. A Method for the Correlation of Biological Activity and Chemical Structure , 1964 .

[32]  H. Kubinyi,et al.  3D QSAR in drug design. , 2002 .

[33]  Hao Liu,et al.  Allosteric Autoinhibition Pathway in Transcription Factor ERG: Dynamics Network and Mutant Experimental Evaluations , 2017, J. Chem. Inf. Model..

[34]  J N Weinstein,et al.  Three-dimensional quantitative structure-activity relationship (QSAR) of HIV integrase inhibitors: a comparative molecular field analysis (CoMFA) study. , 1995, Journal of medicinal chemistry.

[35]  S. Sarafianos,et al.  Design, Synthesis, Biochemical, and Antiviral Evaluations of C6 Benzyl and C6 Biarylmethyl Substituted 2-Hydroxylisoquinoline-1,3-diones: Dual Inhibition against HIV Reverse Transcriptase-Associated RNase H and Polymerase with Antiviral Activities , 2014, Journal of medicinal chemistry.

[36]  Anil K. Saxena,et al.  Consensus Superiority of the Pharmacophore-Based Alignment, Over Maximum Common Substructure (MCS): 3D-QSAR Studies on Carbamates as Acetylcholinesterase Inhibitors , 2009, J. Chem. Inf. Model..

[37]  R. Cramer,et al.  Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. , 1988, Journal of the American Chemical Society.

[38]  R. Luo,et al.  Dynamics Correlation Network for Allosteric Switching of PreQ1 Riboswitch , 2016, Scientific Reports.

[39]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[40]  D. Jochmans Novel HIV-1 reverse transcriptase inhibitors. , 2008, Virus research.