Relaxed complex scheme suggests novel inhibitors for the lyase activity of DNA polymerase beta

Accommodating receptor flexibility and free energy calculation to reduce false positive binders in the discovery of natural products blockers of SARS-COV-2 spike RBD-ACE2 interface

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-COV-2), which causes coronavirus disease-19 (COVID-19) has caused more than 2 million deaths around the globe. The high transmissibility rate of the disease is related to the strong interaction between the virus spike receptor-binding domain (RBD) and the human angiotensin-converting enzyme 2 (ACE2) as documented in several reports. In this study, using state-of-the-art computational methods, natural products were screened and their molecular mechanism to disrupt spike RBD-ACE2 recognition was evaluated. There is the sensitivity of results to receptor ensemble docking calculations. Binding free energy and MD simulation are important tools to evaluate the thermodynamics of binding stability and the capacity of top hits to disrupt RBD-ACE2 recognition. The free energy profiles provide a slight decrease in binding affinity of the virus-receptor interaction. Three flavonoids parvisoflavone B (3), alpinumisoflavone (5) and norisojamicin (2) were effective in blocking the viral entry by binding strongly at the spike RBD-ACE2 interface with the inhibition constant of 0.56, 0.78 and 0.93 μM, respectively. The same compounds demonstrated similar effect on free ACE2 protein. Compound (2), also demonstrated ability to bind strongly on free spike RBD. Well-tempered metadynamics established that parvisoflavone B (3) works by binding to three sites namely interface α, β and loop thereby inhibiting viral cell entry. Owing to their desirable pharmacokinetic properties, the presented top hit natural products are suggested for further SARS-COV-2 molecular targets and subsequent in vitro and in vivo evaluations.

Luteolin: a blocker of SARS-CoV-2 cell entry based on relaxed complex scheme, molecular dynamics simulation, and metadynamics

Natural products have served human life as medications for centuries. During the outbreak of COVID-19, a number of naturally derived compounds and extracts have been tested or used as potential remedies against COVID-19. Tetradenia riparia extract is one of the plant extracts that have been deployed and claimed to manage and control COVID-19 by some communities in Tanzania and other African countries. The active compounds isolated from T. riparia are known to possess various biological properties including antimalarial and antiviral. However, the underlying mechanism of the active compounds against SARS-CoV-2 remains unknown. Results in the present work have been interpreted from the view point of computational methods including molecular dynamics, free energy methods, and metadynamics to establish the related mechanism of action. Among the constituents of T. riparia studied, luteolin inhibited viral cell entry and was thermodynamically stable. The title compound exhibit residence time and unbinding kinetics of 68.86 ms and 0.014 /ms, respectively. The findings suggest that luteolin could be potent blocker of SARS-CoV-2 cell entry. The study shades lights towards identification of bioactive constituents from T. riparia against COVID-19, and thus bioassay can be carried out to further validate such observations.

Identification of Novel Cyclin A2 Binding Site and Nanomolar Inhibitors of Cyclin A2-CDK2 Complex

BACKGROUND

Given the diverse roles of cyclin A2 both in cell cycle regulation and in DNA damage response, identifying small molecule regulators of cyclin A2 activity carries significant potential to regulate diverse cellular processes in both ageing/neurodegeneration and in cancer.

OBJECTIVE

Based on cyclin A2's recently discovered role in DNA repair, we hypothesized that small molecule inhibitors that were predicted to bind to both cyclin A2 and CDK2 will be useful as a radiosensitizer of cancer cells. In this study, we used structure-based drug discovery to find inhibitors that target both cyclin A2 and CDK2.

METHODS

Molecular dynamics simulations were used to generate diverse binding pocket conformations for application of the relaxed complex scheme. We then used structure-based virtual screening to find potential dual cyclin A2 and CDK2 inhibitors. Based on a consensus score of docked poses from Glide and AutoDock Vina, we identified about 40 promising hit compounds, where all PAINS scaffolds were removed from consideration. A biochemical luminescence assay of cyclin A2-CDK2 function was used for experimental verification.

RESULTS

Four lead inhibitors of cyclin A2-CDK2 complex have been identified using a relaxed complex scheme virtual screen have been verified in a biochemical luminescence assay of cyclin A2- CDK2 function. Two of the four lead inhibitors had inhibitory concentrations in the nanomolar range.

CONCLUSION

The four cyclin A2-CDK2 complex inhibitors are the first reported inhibitors that were specifically designed not to target the cyclin A2-CDK2 protein-protein interface. Overall, our results highlight the potential of combined advanced computational tools and biochemical verification to discover novel binding scaffolds.

An assay for DNA polymerase β lyase inhibitors that engage the catalytic nucleophile for binding

DNA polymerase β (Pol β) repairs cellular DNA damage. When such damage is inflicted upon the DNA in tumor cells treated with DNA targeted antitumor agents, Pol β thus diminishes their efficacy. Accordingly, this enzyme has long been a target for antitumor therapy. Although numerous inhibitors of the lyase activity of the enzyme have been reported, none has yet proven adequate for development as a therapeutic agent. In the present study, we developed a new strategy to identify lyase inhibitors that critically engage the lyase active site primary nucleophile Lys72 as part of the binding interface. This involves a parallel evaluation of the effect of the inhibitors on the wild-type DNA polymerase β (Pol β) and Pol β modified with a lysine analogue at position 72. A model panel of five structurally diverse lyase inhibitors identified in our previous studies (only one of which has been published) with unknown modes of binding were used for testing, and one compound, cis-9,10-epoxyoctadecanoic acid, was found to have the desired characteristics. This finding was further corroborated by in silico docking, demonstrating that the predominant mode of binding of the inhibitor involves an important electrostatic interaction between the oxygen atom of the epoxy group and N^ε of the main catalytic nucleophile, Lys72. The strategy, which is designed to identify compounds that engage certain structural elements of the target enzyme, could find broader application for identification of ligands with predetermined sites of binding.

A computational study on the role of water and conformational fluctuations in Hsp90 in response to inhibitors

Molecular chaperone Heat Shock Protein 90 (Hsp90) represents an interesting chemotherapeutic target for cancer treatments as it plays a role in cancer proliferation. Thus, continued effort to identify novel inhibitors of this target is an important task. Drug design using computational approach has gained significant attention in recent years. This work aims to propose docking protocols to re-purpose FDA-approved drugs targeting Hsp90. Sensitivity of results to different docking protocols such apo, holo and receptor ensembles (relaxed complex) structures, the role of water and conformational changes of Hsp90, are described. We show that the protein conformation and water have effects on drug binding. Holo relaxed complex receptors ensembles improves the binding energy of ligands to the protein. We also compare and contrast structural stability of three drugs namely: ezetimibe, pitavastatin and vilazodon in the Hsp90 protein. The results obtained serves as a possible basis towards developing Hsp90 inhibitors.

Computational Characterization of Small Molecules Binding to the Human XPF Active Site and Virtual Screening to Identify Potential New DNA Repair Inhibitors Targeting the ERCC1-XPF Endonuclease

The DNA excision repair protein ERCC-1-DNA repair endonuclease XPF (ERCC1-XPF) is a heterodimeric endonuclease essential for the nucleotide excision repair (NER) DNA repair pathway. Although its activity is required to maintain genome integrity in healthy cells, ERCC1-XPF can counteract the effect of DNA-damaging therapies such as platinum-based chemotherapy in cancer cells. Therefore, a promising approach to enhance the effect of these therapies is to combine their use with small molecules, which can inhibit the repair mechanisms in cancer cells. Currently, there are no structures available for the catalytic site of the human ERCC1-XPF, which performs the metal-mediated cleavage of a DNA damaged strand at 5′. We adopted a homology modeling strategy to build a structural model of the human XPF nuclease domain which contained the active site and to extract dominant conformations of the domain using molecular dynamics simulations followed by clustering of the trajectory. We investigated the binding modes of known small molecule inhibitors targeting the active site to build a pharmacophore model. We then performed a virtual screening of the ZINC Is Not Commercial 15 (ZINC15) database to identify new ERCC1-XPF endonuclease inhibitors. Our work provides structural insights regarding the binding mode of small molecules targeting the ERCC1-XPF active site that can be used to rationally optimize such compounds. We also propose a set of new potential DNA repair inhibitors to be considered for combination cancer therapy strategies.

Successful Identification of Cardiac Troponin Calcium Sensitizers Using a Combination of Virtual Screening and ROC Analysis of Known Troponin C Binders

Calcium-dependent cardiac muscle contraction is regulated by the protein complex troponin. Calcium binds to the N-terminal domain of troponin C (cNTnC) which initiates the process of contraction. Heart failure is a consequence of a disruption of this process. With the prevalence of this condition, a strong need exists to find novel compounds to increase the calcium sensitivity of cNTnC. Desirable are small chemical molecules that bind to the interface between cTnC and the cTnI switch peptide and exhibit calcium sensitizing properties by possibly stabilizing cTnC in an open conformation. To identify novel drug candidates, we employed a structure-based drug discovery protocol that incorporated the use of a relaxed complex scheme (RCS). In preparation for the virtual screening, cNTnC conformations were identified based on their ability to correctly predict known cNTnC binders using a receiver operating characteristics analysis. Following a virtual screen of the National Cancer Institute's Developmental Therapeutic Program database, a small number of molecules were experimentally tested using stopped-flow kinetics and steady-state fluorescence titrations. We identified two novel compounds, 3-(4-methoxyphenyl)-6,7-chromanediol (NSC600285) and 3-(4-methylphenyl)-7,8-chromanediol (NSC611817), that show increased calcium sensitivity of cTnC in the presence of the regulatory domain of cTnI. The effects of NSC600285 and NSC611817 on the calcium dissociation rate was stronger than that of the known calcium sensitizer bepridil. Thus, we identified a 3-phenylchromane group as a possible key pharmacophore in the sensitization of cardiac muscle contraction. Building on this finding is of interest to researchers working on development of drugs for calcium sensitization.

Regulation of USP37 Expression by REST-Associated G9a-Dependent Histone Methylation

The deubiquitylase (DUB) USP37 is a component of the ubiquitin system and controls cell proliferation by regulating the stability of the cyclin-dependent kinase inhibitor 1B, (CDKN1B/p27Kip1). The expression of USP37 is downregulated in human medulloblastoma tumor specimens. In the current study, we show that USP37 prevents medulloblastoma growth in mouse orthotopic models, suggesting that it has tumor-suppressive properties in this neural cancer. Here, we also report on the mechanism underlying USP37 loss in medulloblastoma. Previously, we observed that the expression of USP37 is transcriptionally repressed by the RE1 silencing transcription factor (REST), which requires chromatin remodeling factors for its activity. Genetic and pharmacologic approaches were employed to identify a specific role for G9a, a histone methyltransferase (HMT), in promoting methylation of histone H3 lysine-9 (H3K9) mono- and dimethylation, and surprisingly trimethylation, at the USP37 promoter to repress its gene expression. G9a inhibition also blocked the tumorigenic potential of medulloblastoma cells in vivo Using isogenic low- and high-REST medulloblastoma cells, we further showed a REST-dependent elevation in G9a activity, which further increased mono- and trimethylation of histone H3K9, accompanied by downregulation of USP37 expression. Together, these findings reveal a role for REST-associated G9a and histone H3K9 methylation in the repression of USP37 expression in medulloblastoma.Implications: Reactivation of USP37 by G9a inhibition has the potential for therapeutic applications in REST-expressing medulloblastomas. Mol Cancer Res; 15(8); 1073-84. ©2017 AACR.

Applications of computer-aided approaches in the development of hepatitis C antiviral agents

INTRODUCTION

Hepatitis C virus (HCV) is a global health problem that causes several chronic life-threatening liver diseases. The numbers of people affected by HCV are rising annually. Since 2011, the FDA has approved several anti-HCV drugs; while many other promising HCV drugs are currently in late clinical trials. Areas covered: This review discusses the applications of different computational approaches in HCV drug design. Expert opinion: Molecular docking and virtual screening approaches have emerged as a low-cost tool to screen large databases and identify potential small-molecule hits against HCV targets. Ligand-based approaches are useful for filtering-out compounds with rich physicochemical properties to inhibit HCV targets. Molecular dynamics (MD) remains a useful tool in optimizing the ligand-protein complexes and understand the ligand binding modes and drug resistance mechanisms in HCV. Despite their varied roles, the application of in-silico approaches in HCV drug design is still in its infancy. A more mature application should aim at modelling the whole HCV replicon in its active form and help to identify new effective druggable sites within the replicon system. With more technological advancements, the roles of computer-aided methods are only going to increase several folds in the development of next-generation HCV drugs.

New design of nucleotide excision repair (NER) inhibitors for combination cancer therapy

Many cancer chemotherapy agents act by targeting the DNA of cancer cells, causing substantial damage within their genome and causing them to undergo apoptosis. An effective DNA repair pathway in cancer cells can act in a reverse way by removing these drug-induced DNA lesions, allowing cancer cells to survive, grow and proliferate. In this context, DNA repair inhibitors opened a new avenue in cancer treatment, by blocking the DNA repair mechanisms from removing the chemotherapy-mediated DNA damage. In particular, the nucleotide excision repair (NER) involves more than thirty protein-protein interactions and removes DNA adducts caused by platinum-based chemotherapy. The excision repair cross-complementation group 1 (ERCC1)-xeroderma pigmentosum, complementation group A (XPA) protein (XPA-ERCC1) complex seems to be one of the most promising targets in this pathway. ERCC1 is over expressed in cancer cells and the only known cellular function so far for XPA is to recruit ERCC1 to the damaged point. Here, we build upon our recent advances in identifying inhibitors for this interaction and continue our efforts to rationally design more effective and potent regulators for the NER pathway. We employed in silico drug design techniques to: (1) identify compounds similar to the recently discovered inhibitors, but more effective at inhibiting the XPA-ERCC1 interactions, and (2) identify different scaffolds to develop novel lead compounds. Two known inhibitor structures have been used as starting points for two ligand/structure-hybrid virtual screening approaches. The findings described here form a milestone in discovering novel inhibitors for the NER pathway aiming at improving the efficacy of current platinum-based therapy, by modulating the XPA-ERCC1 interaction.

Entropy in bimolecular simulations: A comprehensive review of atomic fluctuations-based methods

Entropy of binding constitutes a major, and in many cases a detrimental, component of the binding affinity in biomolecular interactions. While the enthalpic part of the binding free energy is easier to calculate, estimating the entropy of binding is further more complicated. A precise evaluation of entropy requires a comprehensive exploration of the complete phase space of the interacting entities. As this task is extremely hard to accomplish in the context of conventional molecular simulations, calculating entropy has involved many approximations. Most of these golden standard methods focused on developing a reliable estimation of the conformational part of the entropy. Here, we review these methods with a particular emphasis on the different techniques that extract entropy from atomic fluctuations. The theoretical formalisms behind each method is explained highlighting its strengths as well as its limitations, followed by a description of a number of case studies for each method. We hope that this brief, yet comprehensive, review provides a useful tool to understand these methods and realize the practical issues that may arise in such calculations.

Clinical significance of a point mutation in DNA polymerase beta (POLB) gene in gastric cancer

Gastric cancer (GC) is a major cause of global cancer mortality. Genetic variations in DNA repair genes can modulate DNA repair capability and, consequently, have been associated with risk of developing cancer. We have previously identified a T to C point mutation at nucleotide 889 (T889C) in DNA polymerase beta (POLB) gene, a key enzyme involved in base excision repair in primary GCs. The purpose of this study was to evaluate the mutation and expression of POLB in a larger cohort and to identify possible prognostic roles of the POLB alterations in GC. Primary GC specimens and their matched normal adjacent tissues were collected at the time of surgery. DNA, RNA and protein samples were isolated from GC specimens and cell lines. Mutations were detected by PCR-RFLP/DHPLC and sequencing analysis. POLB gene expression was examined by RT-PCR, tissue microarray, Western blotting and immunofluorescence assays. The function of the mutation was evaluated by chemosensitivity, MTT, Transwell matrigel invasion and host cell reactivation assays. The T889C mutation was detected in 18 (10.17%) of 177 GC patients. And the T889C mutation was associated with POLB overexpression, lymph nodes metastases and poor tumor differentiation. In addition, patients with- the mutation had significantly shorter survival time than those without-, following postoperative chemotherapy. Furthermore, cell lines with T889C mutation in POLB gene were more resistant to the treatment of 5-fluorouracil, cisplatin and epirubicin than those with wild type POLB. Forced expression of POLB gene with T889C mutation resulted in enhanced cell proliferation, invasion and resistance to anticancer drugs, along with increased DNA repair capability. These results suggest that POLB gene with T889C mutation in surgically resected primary gastric tissues may be clinically useful for predicting responsiveness to chemotherapy in patients with GC. The POLB gene alteration may serve as a prognostic biomarker for GC.

Oxidatively induced DNA damage and its repair in cancer

Oxidatively induced DNA damage is caused in living organisms by endogenous and exogenous reactive species. DNA lesions resulting from this type of damage are mutagenic and cytotoxic and, if not repaired, can cause genetic instability that may lead to disease processes including carcinogenesis. Living organisms possess DNA repair mechanisms that include a variety of pathways to repair multiple DNA lesions. Mutations and polymorphisms also occur in DNA repair genes adversely affecting DNA repair systems. Cancer tissues overexpress DNA repair proteins and thus develop greater DNA repair capacity than normal tissues. Increased DNA repair in tumors that removes DNA lesions before they become toxic is a major mechanism for development of resistance to therapy, affecting patient survival. Accumulated evidence suggests that DNA repair capacity may be a predictive biomarker for patient response to therapy. Thus, knowledge of DNA protein expressions in normal and cancerous tissues may help predict and guide development of treatments and yield the best therapeutic response. DNA repair proteins constitute targets for inhibitors to overcome the resistance of tumors to therapy. Inhibitors of DNA repair for combination therapy or as single agents for monotherapy may help selectively kill tumors, potentially leading to personalized therapy. Numerous inhibitors have been developed and are being tested in clinical trials. The efficacy of some inhibitors in therapy has been demonstrated in patients. Further development of inhibitors of DNA repair proteins is globally underway to help eradicate cancer.

A human ether-á-go-go-related (hERG) ion channel atomistic model generated by long supercomputer molecular dynamics simulations and its use in predicting drug cardiotoxicity

Acquired cardiac long QT syndrome (LQTS) is a frequent drug-induced toxic event that is often caused through blocking of the human ether-á-go-go-related (hERG) K(+) ion channel. This has led to the removal of several major drugs post-approval and is a frequent cause of termination of clinical trials. We report here a computational atomistic model derived using long molecular dynamics that allows sensitive prediction of hERG blockage. It identified drug-mediated hERG blocking activity of a test panel of 18 compounds with high sensitivity and specificity and was experimentally validated using hERG binding assays and patch clamp electrophysiological assays. The model discriminates between potent, weak, and non-hERG blockers and is superior to previous computational methods. This computational model serves as a powerful new tool to predict hERG blocking thus rendering drug development safer and more efficient. As an example, we show that a drug that was halted recently in clinical development because of severe cardiotoxicity is a potent inhibitor of hERG in two different biological assays which could have been predicted using our new computational model.

A Refined Model of the HCV NS5A protein bound to daclatasvir explains drug-resistant mutations and activity against divergent genotypes

Many direct-acting antiviral agents (DAAs) that selectively block hepatitis C virus (HCV) replication are currently under development. Among these agents is Daclatasvir, a first-in-class inhibitor targeting the NS5A viral protein. Although Daclatasvir is the most potent HCV antiviral molecule yet developed, its binding location and mode of binding remain unknown. The drug exhibits a low barrier to resistance mutations, particularly in genotype 1 viruses, but its efficacy against other genotypes is unclear. Using state-of-the-art modeling techniques combined with the massive computational power of Blue Gene/Q, we identified the atomic interactions of Daclatasvir within NS5A for different HCV genotypes and for several reported resistant mutations. The proposed model is the first to reveal the detailed binding mode of Daclatasvir. It also provides a tool to facilitate design of second generation drugs, which may confer less resistance and/or broader activity against HCV.

Detailed computational study of the active site of the hepatitis C viral RNA polymerase to aid novel drug design

The hepatitis C virus (HCV) RNA polymerase, NS5B, is a leading target for novel and selective HCV drug design. The enzyme has been the subject of intensive drug discovery aimed at developing direct acting antiviral (DAA) agents that inhibit its activity and hence prevent the virus from replicating its genome. In this study, we focus on one class of NS5B inhibitors, namely nucleos(t)ide mimetics. Forty-one distinct nucleotide structures have been modeled within the active site of NS5B for the six major HCV genotypes. Our comprehensive modeling protocol employed 287 different molecular dynamics simulations combined with the molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) methodology to rank and analyze these structures for all genotypes. The binding interactions of the individual compounds have been investigated and reduced to the atomic level. The present study significantly refines our understanding of the mode of action of NS5B-nucleotide-inhibitors, identifies the key structural elements necessary for their activity, and implements the tools for ranking the potential of additional much needed novel inhibitors of NS5B.

Protein flexibility in docking and surface mapping

Structure-based drug design has become an essential tool for rapid lead discovery and optimization. As available structural information has increased, researchers have become increasingly aware of the importance of protein flexibility for accurate description of the native state. Typical protein-ligand docking efforts still rely on a single rigid receptor, which is an incomplete representation of potential binding conformations of the protein. These rigid docking efforts typically show the best performance rates between 50 and 75%, while fully flexible docking methods can enhance pose prediction up to 80-95%. This review examines the current toolbox for flexible protein-ligand docking and receptor surface mapping. Present limitations and possibilities for future development are discussed.

Discovery of small molecule inhibitors that interact with γ-tubulin

Recent studies have shown an overexpression of γ-tubulin in human glioblastomas and glioblastoma cell lines. As the 2-year survival rate for glioblastoma is very poor, potential benefit exists for discovering novel chemotherapeutic agents that can inhibit γ-tubulin, which is known to form a ring complex that acts as a microtubule nucleation center. We present experimental evidence that colchicine and combretastatin A-4 bind to γ-tubulin, which are to our knowledge the first drug-like compounds known to interact with γ-tubulin. Molecular dynamics simulations and docking studies were used to analyze the hypothesized γ-tubulin binding domain of these compounds. The suitability of the potential binding modes was evaluated and suggests the subsequent rational design of novel targeted inhibitors of γ-tubulin.