Molecular Dynamics Simulations of 2‐Amino‐6‐arylsulphonylbenzonitriles Analogues as HIV Inhibitors: Interaction Modes and Binding Free Energies

Molecular dynamics (MD) simulations in water environment were carried out on the HIV‐1 reverse transcriptase (RT), and its complexes with one representative of each of three series of inhibitors: 2‐amino‐6‐arylsulphonylbenzonitriles and their thio and sulphinyl congeners. Molecular Mechanics Generalized Born Surface Area (MM‐GBSA) was used to calculate the binding free energy based on the obtained MD trajectories. Calculated energies are correlated to activity. A comparison of interaction modes, binding free energy, contributions of the residues to the binding free energy and H‐bonds was carried out with the average structures. The results show that there exist different interaction modes between RT and ligands and different specific interactions with some residues. The higher binding affinity of the most potent inhibitor in the series of molecules under study is favoured by electrostatic interactions and solvation contribution.

[1]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[2]  A. D. Clark,et al.  Crystal structure of HIV‐1 reverse transcriptase in complex with a polypurine tract RNA:DNA , 2001, The EMBO journal.

[3]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[4]  K. Roy,et al.  Classical QSAR Modeling of HIV‐1 Reverse Transcriptase Inhibitor 2‐Amino‐6‐arylsulfonylbenzonitriles and Congeners , 2004 .

[5]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[6]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[7]  Yvonne Jones,et al.  Mechanism of inhibition of HIV-1 reverse transcriptase by non-nucleoside inhibitors , 1995, Nature Structural Biology.

[8]  T. Steitz,et al.  Crystal structure at 3.5 A resolution of HIV-1 reverse transcriptase complexed with an inhibitor. , 1992, Science.

[9]  J Andrew McCammon,et al.  Elucidating the inhibition mechanism of HIV-1 non-nucleoside reverse transcriptase inhibitors through multicopy molecular dynamics simulations. , 2009, Journal of molecular biology.

[10]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[11]  Bhyravabhotla Jayaram,et al.  Solvation Free Energy of Biomacromolecules: Parameters for a Modified Generalized Born Model Consistent with the AMBER Force Field , 1998 .

[12]  Y. Cheng,et al.  Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. , 1973, Biochemical pharmacology.

[13]  J. Doucet,et al.  QSAR models for 2-amino-6-arylsulfonylbenzonitriles and congeners HIV-1 reverse transcriptase inhibitors based on linear and nonlinear regression methods. , 2009, European journal of medicinal chemistry.

[14]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[15]  G L Verdine,et al.  Structure of a covalently trapped catalytic complex of HIV-1 reverse transcriptase: implications for drug resistance. , 1998, Science.

[16]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[17]  E. De Clercq,et al.  New anti-HIV agents and targets. , 2002, Medicinal research reviews.

[18]  A. D. Clark,et al.  Structure of unliganded HIV-1 reverse transcriptase at 2.7 A resolution: implications of conformational changes for polymerization and inhibition mechanisms. , 1996, Structure.

[19]  A. D. Clark,et al.  Crystal structure of human immunodeficiency virus type 1 reverse transcriptase complexed with double-stranded DNA at 3.0 A resolution shows bent DNA. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[20]  D I Stuart,et al.  2-Amino-6-arylsulfonylbenzonitriles as non-nucleoside reverse transcriptase inhibitors of HIV-1. , 2001, Journal of medicinal chemistry.

[21]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998, J. Comput. Chem..

[22]  Jörg Weiser,et al.  Approximate atomic surfaces from linear combinations of pairwise overlaps (LCPO) , 1999, J. Comput. Chem..

[23]  P A Kollman,et al.  Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution. , 1998, Journal of biomolecular structure & dynamics.

[24]  Yvonne Jones,et al.  High resolution structures of HIV-1 RT from four RT–inhibitor complexes , 1995, Nature Structural Biology.

[25]  Ruisheng Zhang,et al.  Receptor- and ligand-based 3D-QSAR study for a series of non-nucleoside HIV-1 reverse transcriptase inhibitors. , 2009, Bioorganic & medicinal chemistry.

[26]  D. Stuart,et al.  The structure of HIV-1 reverse transcriptase complexed with 9-chloro-TIBO: lessons for inhibitor design. , 1995, Structure.

[27]  M. P. Freitas.  MIA-QSAR modelling of anti-HIV-1 activities of some 2-amino-6-arylsulfonylbenzonitriles and their thio and sulfinyl congeners. , 2006, Organic & biomolecular chemistry.

[28]  Kunal Roy,et al.  QSAR modeling of HIV-1 reverse transcriptase inhibitor 2-amino-6-arylsulfonylbenzonitriles and congeners using molecular connectivity and E-state parameters. , 2004, Bioorganic & medicinal chemistry.

[29]  P A Kollman,et al.  Continuum solvent studies of the stability of RNA hairpin loops and helices. , 1998, Journal of biomolecular structure & dynamics.