Kamlesh Sahu

Anesthetic Alterations of Collective Terahertz Oscillations in Tubulin Correlate with Clinical Potency: Implications for Anesthetic Action and Post-Operative Cognitive Dysfunction

Anesthesia blocks consciousness and memory while sparing non-conscious brain activities. While the exact mechanisms of anesthetic action are unknown, the Meyer-Overton correlation provides a link between anesthetic potency and solubility in a lipid-like, non-polar medium. Anesthetic action is also related to an anesthetic’s hydrophobicity, permanent dipole, and polarizability, and is accepted to occur in lipid-like, non-polar regions within brain proteins. Generally the protein target for anesthetics is assumed to be neuronal membrane receptors and ion channels, however new evidence points to critical effects on intra-neuronal microtubules, a target of interest due to their potential role in post-operative cognitive dysfunction (POCD). Here we use binding site predictions on tubulin, the protein subunit of microtubules, with molecular docking simulations, quantum chemistry calculations, and theoretical modeling of collective dipole interactions in tubulin to investigate the effect of a group of gases including anesthetics, non-anesthetics, and anesthetic/convulsants on tubulin dynamics. We found that these gases alter collective terahertz dipole oscillations in a manner that is correlated with their anesthetic potency. Understanding anesthetic action may help reveal brain mechanisms underlying consciousness, and minimize POCD in the choice and development of anesthetics used during surgeries for patients suffering from neurodegenerative conditions with compromised cytoskeletal microtubules.

Metabolic targeting of EGFRvIII/PDK1 axis in temozolomide resistant glioblastoma

Glioblastomas are characterized by amplification of EGFR. Approximately half of tumors with EGFR over-expression also express a constitutively active ligand independent EGFR variant III (EGFRvIII). While current treatments emphasize surgery followed by radiation and chemotherapy with Temozolomide (TMZ), acquired chemoresistance is a universal feature of recurrent GBMs. To mimic the GBM resistant state, we generated an in vitro TMZ resistant model and demonstrated that dichloroacetate (DCA), a metabolic inhibitor of pyruvate dehydrogenase kinase 1 (PDK1), reverses the Warburg effect. Microarray analysis conducted on the TMZ resistant cells with their subsequent treatment with DCA revealed PDK1 as its sole target. DCA treatment also induced mitochondrial membrane potential change and apoptosis as evidenced by JC-1 staining and electron microscopic studies. Computational homology modeling and docking studies confirmed DCA binding to EGFR, EGFRvIII and PDK1 with high affinity. In addition, expression of EGFRvIII was comparable to PDK1 when compared to EGFR in GBM surgical specimens supporting our in silico prediction data. Collectively our current study provides the first in vitro proof of concept that DCA reverses the Warburg effect in the setting of EGFRvIII positivity and TMZ resistance leading to GBM cytotoxicity, implicating cellular tyrosine kinase signaling in cancer cell metabolism.

Chicoric acid binds to two sites and decreases the activity of the YopH bacterial virulence factor

Chicoric acid (CA) is a phenolic compound present in dietary supplements with a large spectrum of biological properties reported ranging from antioxidant, to antiviral, to immunostimulatory properties. Due to the fact that chicoric acid promotes phagocytic activity and was reported as an allosteric inhibitor of the PTP1B phosphatase, we examined the effect of CA on YopH phosphatase from pathogenic bacteria, which block phagocytic processes of a host cell. We performed computational studies of chicoric acid binding to YopH as well as validation experiments with recombinant enzymes. In addition, we performed similar studies for caffeic and chlorogenic acids to compare the results. Docking experiments demonstrated that, from the tested compounds, only CA binds to both catalytic and secondary binding sites of YopH. Our experimental results showed that CA reduces activity of recombinant YopH phosphatase from Yersinia enterocolitica and human CD45 phosphatase. The inhibition caused by CA was irreversible and did not induce oxidation of catalytic cysteine. We proposed that inactivation of YopH induced by CA is involved with allosteric inhibition by interacting with essential regions responsible for ligand binding.

α-Synuclein dimer structures found from computational simulations.

Dimer formation is likely the first step in the oligomerization of α-synuclein in Lewy bodies. In order to prevent α-synuclein aggregation, knowledge of the atomistic structures of possible α-synuclein dimers and the interaction affinity between the dimer domains is a necessary prerequisite in the process of rational design of dimerization inhibitors. Using computational methodology, we have investigated several possible α-synuclein dimer structures, focusing on dimers formed from α-helical forms of the protein found when it is membrane-bound, and dimers formed from β-sheet conformations predicted by simulations. Structures and corresponding binding affinities for the interacting monomers in possible α-synuclein dimers, along with properties including the contributions from different interaction energies and the radii of gyration, were found through molecular docking followed by MD simulations and binding-energy calculations. We found that even though α-synuclein is highly charged, hydrophobic contributions play a significant role in stabilizing dimers.

Mathematical and computational modeling in biology at multiple scales

A variety of topics are reviewed in the area of mathematical and computational modeling in biology, covering the range of scales from populations of organisms to electrons in atoms. The use of maximum entropy as an inference tool in the fields of biology and drug discovery is discussed. Mathematical and computational methods and models in the areas of epidemiology, cell physiology and cancer are surveyed. The technique of molecular dynamics is covered, with special attention to force fields for protein simulations and methods for the calculation of solvation free energies. The utility of quantum mechanical methods in biophysical and biochemical modeling is explored. The field of computational enzymology is examined.

Aurintricarboxylic acid structure modifications lead to reduction of inhibitory properties against virulence factor YopH and higher cytotoxicity.

Yersinia sp. bacteria owe their viability and pathogenic virulence to the YopH factor, which is a highly active bacterial protein tyrosine phosphatase. Inhibition of YopH phosphatase results in the lack of Yersinia sp. pathogenicity. We have previously described that aurintricarboxylic acid inhibits the activity of YopH at nanomolar concentrations and represents a unique mechanism of YopH inactivation due to a redox process. This work is a continuation of our previous studies. Here we show that modifications of the structure of aurintricarboxylic acid reduce the ability to inactivate YopH and lead to higher cytotoxicity. In the present paper we examine the inhibitory properties of aurintricarboxylic acid analogues, such as eriochrome cyanine R (ECR) and pararosaniline. Computational docking studies we report here indicate that ATA analogues are not precluded to bind in the YopH active site and in all obtained binding conformations ECR and pararosaniline bind to YopH active site. The free binding energy calculations show that ECR has a stronger binding affinity to YopH than pararosaniline, which was confirmed by experimental YopH enzymatic activity studies. We found that ATA analogues can reversibly reduce the enzymatic activity of YopH, but possess weaker inhibitory properties than ATA. The ATA analogues induced inactivation of YopH is probably due to oxidative mechanism, as pretreatment with catalase prevents from inhibition. We also found that ATA analogues significantly decrease the viability of macrophage cells, especially pararosaniline, while ATA reveals only slight effect on cell viability.

Abstract 41: DCA bind to EGFR/EGFRvIII/PDK1 and affect the proliferation and growth in TMZ resistant glioblastoma model

Glioblastomas are characterized by amplification of epidermal growth factor receptor (EGFR). Approximately half of the GBM tumors with EGFR over-expression also express a constitutively active ligand independent EGFR variant III (EGFRvIII). This phenotype represents a tumor specific target correlating with a high growth potential. Despite advances in treatment, whether surgery, chemotherapeutics or molecular targeted therapies, very few have shown promise. Temozolomide (TMZ), as standard of care alkylation therapy, demonstrates improved survival when administered with concomitant radiotherapy. Unfortunately, acquired TMZ chemo-resistance is observed in more than 90% of recurrent GBMs. To understand the mechanisms in the context of TMZ resistance, we generated an in vitro TMZ resistant model by continuous exposure of U373 cells constitutively expressing EGFRvIII (U373vIII) to 150uM TMZ for 6 months (U373vIIIR). Dicholoroacetate (DCA) is a known metabolic inhibitor of pyruvate dehydrogenase kinase 1 (PDK1). Its application to cancer has recently been revived in the realm of metabolic oncology, where reversal of the Warburg effect has resulted in decreased tumor growth. Treatment of U373, U373vIII, and U373vIIIR GBM cell lines with DCA was found to induce cytotoxicity and reduced cell survival across all cell lines. Further, micro array studies conducted on U373vIII or U373vIIIR and their subsequent treatment with DCA revealed that PDK1 is the sole target in TMZ resistant tumors expressing EGFRvIII. Additionally, we demonstrated that 1mM DCA induced mitochondrial membrane potential change as evidenced in JC-1 staining and electron microscopic studies confirming mitochondrial apoptosis. Consistent with our previous findings that EGFR interacts with PDK1, our computational analysis revealed that DCA aligns itself to the binding sites of both EGFR and EGFRvIII with strong binding affinity apart from its binding to PDK1. Clinically, expression of EGFRvIII correlated with PDK1 when compared to EGFR in GBM surgical specimens. Collectively, our studies demonstrate that principles of metabolic oncology and DCA are active in GBM treatment despite resistance providing new insight into the development of alternative treatment in TMZ failure. Citation Format: Kiran Kumar Velpula, Kamlesh Sahu, Jack Tuszynski, Maheedhara R. Guda, Swapna Asuthkar, Sarah E. Martin, Justin D. Lathia, Andrew J. Tsung. DCA bind to EGFR/EGFRvIII/PDK1 and affect the proliferation and growth in TMZ resistant glioblastoma model. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 41.

Theoretical investigation of the interactions in binding pocket of Reverse Transcriptase.

Interactions in proteins have been studied using several chemical information techniques including quantum chemical methods that are applied to truncated systems composed of the ligand molecule and the surrounding amino acids of the receptor. In this work we adopt an approach to study these interactions accounting for as many as possible explicit solvent molecules and without the need of a fragmented calculation. Furthermore, we embed our quantum chemical calculations within a molecular dynamics framework that enables a fundamentally fast system for quantum molecular dynamic simulations (QCMD). Central to this new system for QCMD is the tight binding QC system, newly developed in our laboratories, and which combined with the MD paradigm results in an ultra accelerated QCMD method for protein–ligand interaction evaluations. We have applied our newly developed method to the Nevirapine (NVP)–Reverse Transcriptase (RT) system. We show how the proposed method leads us to new findings. The advanced QCMD was applied to a system of RT with NVP and it has led to the knowledge of specific groups and atoms that interact with surrounding amino acids of RT and help in drug binding. The information derived from this calculation may be used in designing drugs for NVP resistant virus strains that have binding capability like NVP.