Receptor- and ligand-based study of fullerene analogues: comprehensive computational approach including quantum-chemical, QSAR and molecular docking simulations.

Fullerene and its derivatives have potential antiviral activity due to their specific binding interactions with biological molecules. In this study fullerene derivatives were investigated by the synergic combination of three approaches: quantum-mechanical calculations, protein-ligand docking and quantitative structure-activity relationship methods. The protein-ligand docking studies and improved structure-activity models have been able both to predict binding affinities for the set of fullerene-C60 derivatives and to help in finding mechanisms of fullerene derivative interactions with human immunodeficiency virus type 1 aspartic protease, HIV-1 PR. Protein-ligand docking revealed several important molecular fragments that are responsible for the interaction with HIV-1 PR. In addition, a density functional theory method has been utilized to identify the optimal geometries and predict physico-chemical parameters of the studied compounds. The 5-variable GA-MLRA based model showed the best predictive ability (r(2)training = 0.882 and r(2)test = 0.738), with high internal and external correlation coefficients.

[1]  Soumi Banerjee.  Molecular dynamics study of self-agglomeration of charged fullerenes in solvents. , 2013, The Journal of chemical physics.

[2]  A. Ferrarini,et al.  Monomeric fullerenes in lipid membranes: effects of molecular shape and polarity. , 2011, Langmuir : the ACS journal of surfaces and colloids.

[3]  Chi‐Huey Wong,et al.  HIV-1 protease: mechanism and drug discovery. , 2003, Organic & biomolecular chemistry.

[4]  Sanjay Mathur,et al.  Graphene oxide, but not fullerenes, targets immunoproteasomes and suppresses antigen presentation by dendritic cells. , 2013, Small.

[5]  L. Echegoyen,et al.  Electrochemistry of Fullerenes and Their Derivatives , 1998 .

[6]  Jerzy Leszczynski,et al.  Improved model for fullerene C60 solubility in organic solvents based on quantum-chemical and topological descriptors , 2011 .

[7]  Yuan-yuan Li,et al.  Protective effect of C60‐methionine derivate on lead‐exposed human SH‐SY5Y neuroblastoma cells , 2011, Journal of applied toxicology : JAT.

[8]  Alexandra Shulman-Peleg,et al.  Interaction of c(60)-fullerene and carboxyfullerene with proteins: docking and binding site alignment. , 2006, Bioconjugate chemistry.

[9]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[10]  N. Higashi,et al.  Immobilization and cleavage of DNA at cationic, self-assembledmonolayers containing C60 on gold , 1997 .

[11]  H. Okayama,et al.  Gene delivery by aminofullerenes: structural requirements for efficient transfection. , 2006, Chemistry, an Asian journal.

[12]  Xue Xu,et al.  Multidrug resistance protein P-glycoprotein does not recognize nanoparticle C60: experiment and modeling , 2012 .

[13]  Schmid,et al.  "Scaffold-Hopping" by Topological Pharmacophore Search: A Contribution to Virtual Screening. , 1999, Angewandte Chemie.

[14]  Serdar Durdagi,et al.  3D QSAR CoMFA/CoMSIA, molecular docking and molecular dynamics studies of fullerene-based HIV-1 PR inhibitors. , 2008, Bioorganic & medicinal chemistry letters.

[15]  R. Zhou,et al.  Amino acid analogues bind to carbon nanotube via π-π interactions: comparison of molecular mechanical and quantum mechanical calculations. , 2012, The Journal of chemical physics.

[16]  Yong Pan,et al.  Prediction of the net heat of combustion of organic compounds based on atom-type electrotopological state indices , 2009 .

[17]  A. Tropsha,et al.  Beware of q2! , 2002, Journal of molecular graphics & modelling.

[18]  P. Labute.  A widely applicable set of descriptors. , 2000, Journal of molecular graphics & modelling.

[19]  T. Luh,et al.  Carboxyfullerene Prevents Iron‐Induced Oxidative Stress in Rat Brain , 1999, Journal of neurochemistry.

[20]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[21]  K. Roy,et al.  Exploring quantitative structure–activity relationship studies of antioxidant phenolic compounds obtained from traditional Chinese medicinal plants , 2010 .

[22]  Ruhong Zhou,et al.  Binding preference of carbon nanotube over proline-rich motif ligand on SH3-domain: a comparison with different force fields. , 2013, The journal of physical chemistry. B.

[23]  Ya-Wen Chen,et al.  Fullerene derivatives protect against oxidative stress in RAW 264.7 cells and ischemia-reperfused lungs. , 2004, American journal of physiology. Regulatory, integrative and comparative physiology.

[24]  Zafar Iqbal,et al.  Single-walled Carbon Nanotubes Are a New Class of Ion Channel Blockers* , 2003, Journal of Biological Chemistry.

[25]  Johann Gasteiger,et al.  Deriving the 3D structure of organic molecules from their infrared spectra , 1999 .

[26]  Davide Giust,et al.  [60]Fullerene-based monolayers as neuroprotective biocompatible hybrid materials. , 2011, Chemical Communications.

[27]  Andrew G Mercader,et al.  QSAR study of flavonoids and biflavonoids as influenza H1N1 virus neuraminidase inhibitors. , 2010, European journal of medicinal chemistry.

[28]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[29]  C. Wilkins,et al.  Synthesis of a fullerene derivative for the inhibition of HIV enzymes , 1993 .

[30]  Mark A. Ratner,et al.  6‐31G* basis set for third‐row atoms , 2001, J. Comput. Chem..

[31]  Anderson Coser Gaudio,et al.  BuildQSAR: A New Computer Program for QSAR Analysis , 2000 .

[32]  Maykel Pérez González,et al.  2D-autocorrelation descriptors for predicting cytotoxicity of naphthoquinone ester derivatives against oral human epidermoid carcinoma. , 2007, Bioorganic & medicinal chemistry.

[33]  Jarmo Huuskonen,et al.  Estimation of Aqueous Solubility for a Diverse Set of Organic Compounds Based on Molecular Topology , 2000, J. Chem. Inf. Comput. Sci..

[34]  Carsten Bingel,et al.  Cyclopropanierung von Fullerenen , 1993 .

[35]  T. Puzyn,et al.  Toward the development of "nano-QSARs": advances and challenges. , 2009, Small.

[36]  Stephen R. Wilson,et al.  [60]fullerene is a powerful antioxidant in vivo with no acute or subacute toxicity. , 2005, Nano letters.

[37]  Paola Gramatica,et al.  Structure/Response Correlations and Similarity/Diversity Analysis by GETAWAY Descriptors, 2. Application of the Novel 3D Molecular Descriptors to QSAR/QSPR Studies , 2002, J. Chem. Inf. Comput. Sci..

[38]  Eiichi Nakamura,et al.  Photoinduced biochemical activity of fullerene carboxylic acid , 1993 .

[39]  Zhen Cao,et al.  Molecular Dynamics Simulation of Fullerene C60 in Ethanol Solution , 2009 .

[40]  Igor V. Tetko,et al.  Neural Network Modeling for Estimation of Partition Coefficient Based on Atom-Type Electrotopological State Indices , 2000, J. Chem. Inf. Comput. Sci..

[41]  Magdalena Bacilieri,et al.  Prediction of the aqueous solvation free energy of organic compounds by using autocorrelation of molecular electrostatic potential surface properties combined with response surface analysis. , 2008, Bioorganic & medicinal chemistry.

[42]  Stephen R. Wilson,et al.  Trisamine C(60)-fullerene adducts inhibit neuronal nitric oxide synthase by acting as highly potent calmodulin antagonists. , 2002, Archives of biochemistry and biophysics.

[44]  Maurizio Prato,et al.  Excited-State Properties of C60 Fullerene Derivatives , 2000 .

[45]  V. Knight,et al.  A fullerene-paclitaxel chemotherapeutic: synthesis, characterization, and study of biological activity in tissue culture. , 2005, Journal of the American Chemical Society.

[46]  L. Ho,et al.  Local carboxyfullerene protects cortical infarction in rat brain , 2002, Neuroscience Research.

[47]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

[48]  Li Xiao,et al.  The water-soluble fullerene derivative "Radical Sponge" exerts cytoprotective action against UVA irradiation but not visible-light-catalyzed cytotoxicity in human skin keratinocytes. , 2006, Bioorganic & medicinal chemistry letters.

[49]  R. Luxenhofer,et al.  Neuronal uptake and intracellular superoxide scavenging of a fullerene (C60)-poly(2-oxazoline)s nanoformulation. , 2011, Biomaterials.

[50]  G. Miller,et al.  Nanobionics of Pharmacologically Active Derivatives of Fullerene C60 , 2003 .

[51]  Stephen R. Wilson,et al.  Anti-human immunodeficiency virus activity and cytotoxicity of derivatized buckminsterfullerenes , 1996 .

[52]  P. Roy,et al.  On Some Aspects of Variable Selection for Partial Least Squares Regression Models , 2008 .

[53]  Serdar Durdagi,et al.  Computational design of novel fullerene analogues as potential HIV-1 PR inhibitors: Analysis of the binding interactions between fullerene inhibitors and HIV-1 PR residues using 3D QSAR, molecular docking and molecular dynamics simulations. , 2008, Bioorganic & medicinal chemistry.

[54]  Jerzy Leszczynski,et al.  Multiplicative SMILES-based optimal descriptors: QSPR modeling of fullerene C60 solubility in organic solvents , 2007 .

[55]  Jerzy Leszczynski,et al.  Structure-toxicity relationships of nitroaromatic compounds , 2006, Molecular Diversity.

[56]  N. Monteiro-Riviere,et al.  Penetration of intact skin by quantum dots with diverse physicochemical properties. , 2006, Toxicological sciences : an official journal of the Society of Toxicology.

[57]  Jerzy Leszczynski,et al.  Choosing safe dispersing media for C60 fullerenes by using cytotoxicity tests on the bacterium Escherichia coli. , 2010, Journal of hazardous materials.

[58]  Roberto Todeschini,et al.  Structure/Response Correlations and Similarity/Diversity Analysis by GETAWAY Descriptors, 1. Theory of the Novel 3D Molecular Descriptors , 2002, J. Chem. Inf. Comput. Sci..

[59]  Fred Wudl,et al.  Inhibition of the HIV-1 protease by fullerene derivatives: model building studies and experimental verification , 1993 .

[60]  J. Leszczynski,et al.  Aconitum and Delphinium alkaloids of curare-like activity. QSAR analysis and molecular docking of alkaloids into AChBP. , 2010, European journal of medicinal chemistry.

[61]  V. Subramanian,et al.  Exploring the changes in the structure of α-helical peptides adsorbed onto a single walled carbon nanotube using classical molecular dynamics simulation. , 2010, The journal of physical chemistry. B.

[62]  Jerzy Leszczynski,et al.  Fullerene (C60) forms stable complex with nucleic acid base guanine , 2009 .

[63]  M. Shukla,et al.  Interaction of nucleic acid bases and Watson–Crick base pairs with fullerene: Computational study , 2010 .

[64]  Jan Balzarini,et al.  Synthesis and anti-HIV properties of new water-soluble bis-functionalized[60]fullerene derivatives. , 2003, Bioorganic & medicinal chemistry letters.

[65]  Jerzy Leszczynski,et al.  A Quantitative Structure‐Activity Relationship (QSAR) Study of the Antioxidant Activity of Flavonoids , 2005 .

[66]  Gerhard F. Ecker,et al.  Topological Distance Based 3D Descriptors for Use in QSAR and Diversity Analysis , 2004, J. Chem. Inf. Model..

[67]  Michael R Hamblin,et al.  Photodynamic therapy with fullerenes in vivo: reality or a dream? , 2011, Nanomedicine.

[68]  M. Prato,et al.  Anti-HIV properties of cationic fullerene derivatives. , 2005, Bioorganic & medicinal chemistry letters.

[69]  A. Mikos,et al.  Water-Soluble Fullerene (C60) Derivatives as Nonviral Gene-Delivery Vectors , 2008, Molecular pharmaceutics.

[70]  Petra Schneider,et al.  Distance phenomena in high‐dimensional chemical descriptor spaces: Consequences for similarity‐based approaches , 2009, J. Comput. Chem..

[71]  James Andrew McCammon,et al.  A virtual screening study of the acetylcholine binding protein using a relaxed-complex approach , 2009, Comput. Biol. Chem..

[72]  Yves Rubin,et al.  Sequential {open_quotes}double-Michael{close_quotes} additions of dienolates with C{sub 60}: Rapid access to sterically congested Buckminsterfullerene derivatives with defined stereochemistry , 1995 .

[73]  Kunal Roy,et al.  Some case studies on application of “rm2” metrics for judging quality of quantitative structure–activity relationship predictions: Emphasis on scaling of response data , 2013, J. Comput. Chem..

[74]  S. C. O'brien,et al.  C60: Buckminsterfullerene , 1985, Nature.

[75]  Jerzy Leszczynski,et al.  Using nano-QSAR to predict the cytotoxicity of metal oxide nanoparticles. , 2011, Nature nanotechnology.

[76]  Jerzy Leszczynski,et al.  SMILES‐based optimal descriptors: QSAR analysis of fullerene‐based HIV‐1 PR inhibitors by means of balance of correlations , 2009, J. Comput. Chem..

[77]  Jerzy Leszczynski,et al.  InChI-based optimal descriptors: QSAR analysis of fullerene[C60]-based HIV-1 PR inhibitors by correlation balance. , 2010, European journal of medicinal chemistry.

[78]  G L Kenyon,et al.  Optimizing the binding of fullerene inhibitors of the HIV-1 protease through predicted increases in hydrophobic desolvation. , 1998, Journal of medicinal chemistry.

[79]  A. Tropsha,et al.  Quantitative nanostructure-activity relationship modeling. , 2010, ACS nano.

[80]  Jerzy Leszczynski,et al.  QSAR modeling of measured binding affinity for fullerene-based HIV-1 PR inhibitors by CORAL , 2010 .

[81]  J. Huuskonen.  QSAR modeling with the electrotopological state indices: predicting the toxicity of organic chemicals. , 2003, Chemosphere.

[82]  J. Huff,et al.  HIV protease: a novel chemotherapeutic target for AIDS. , 1991, Journal of medicinal chemistry.

[83]  Ruhong Zhou,et al.  Plugging into proteins: poisoning protein function by a hydrophobic nanoparticle. , 2010, ACS nano.