TarFisDock: a web server for identifying drug targets with docking approach

TarFisDock is a web-based tool for automating the procedure of searching for small molecule–protein interactions over a large repertoire of protein structures. It offers PDTD (potential drug target database), a target database containing 698 protein structures covering 15 therapeutic areas and a reverse ligand–protein docking program. In contrast to conventional ligand–protein docking, reverse ligand–protein docking aims to seek potential protein targets by screening an appropriate protein database. The input file of this web server is the small molecule to be tested, in standard mol2 format; TarFisDock then searches for possible binding proteins for the given small molecule by use of a docking approach. The ligand–protein interaction energy terms of the program DOCK are adopted for ranking the proteins. To test the reliability of the TarFisDock server, we searched the PDTD for putative binding proteins for vitamin E and 4H-tamoxifen. The top 2 and 10% candidates of vitamin E binding proteins identified by TarFisDock respectively cover 30 and 50% of reported targets verified or implicated by experiments; and 30 and 50% of experimentally confirmed targets for 4H-tamoxifen appear amongst the top 2 and 5% of the TarFisDock predicted candidates, respectively. Therefore, TarFisDock may be a useful tool for target identification, mechanism study of old drugs and probes discovered from natural products. TarFisDock and PDTD are available at

Glycan Epitopes and Potential Glycoside Antagonists of DC-SIGN Involved in COVID-19: In Silico Study

Glycosylation is an important post-translational modification that affects a wide variety of physiological functions. DC-SIGN (Dendritic Cell-Specific Intercellular adhesion molecule-3-Grabbing Non-integrin) is a protein expressed in antigen-presenting cells that recognizes a variety of glycan epitopes. Until now, the binding of DC-SIGN to SARS-CoV-2 Spike glycoprotein has been reported in various articles and is regarded to be a factor in systemic infection and cytokine storm. The mechanism of DC-SIGN recognition offers an alternative method for discovering new medication for COVID-19 treatment. Here, we discovered three potential pockets that hold different glycan epitopes by performing molecular dynamics simulations of previously reported oligosaccharides. The “EPN” motif, “NDD” motif, and Glu354 form the most critical pocket, which is known as the Core site. We proposed that the type of glycan epitopes, rather than the precise amino acid sequence, determines the recognition. Furthermore, we deduced that oligosaccharides could occupy an additional site, which adds to their higher affinity than monosaccharides. Based on our findings and previously described glycoforms on the SARS-CoV-2 Spike, we predicted the potential glycan epitopes for DC-SIGN. It suggested that glycan epitopes could be recognized at multiple sites, not just Asn234, Asn149 and Asn343. Subsequently, we found that Saikosaponin A and Liquiritin, two plant glycosides, were promising DC-SIGN antagonists in silico

Exploring the RING-catalyzed ubiquitin transfer mechanism by MD and QM/MM calculations.

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors

Targeting the Small- and Intermediate Conductance Ca2+- Activated Potassium Channels: The Drug Binding Pocket at the Channel/Calmodulin Interface

The small- and intermediate-conductance Ca 2+ -activated potassium (SK/IK) channels play important roles in the regulation of excitable cells in both the central nervous and cardiovascular systems. Evidence from animal models has implicated SK/IK channels in neurological conditions such as ataxia and alcohol use disorders. Further, genome-wide association studies have suggested that cardiovascular abnormalities such as arrhythmias and hypertension are associated with single nucleotide polymorphisms that occur within the genes encoding the SK/IK channels. The Ca 2+ sensitivity of the SK/IK channels stems from a constitutively bound Ca 2+ -binding protein: calmodulin. Small-molecule positive modulators of SK/IK channels have been developed over the past decade, and recent structural studies have revealed that the binding pocket of these positive modulators is located at the interface between the channel and calmodulin. SK/IK channel positive modulators can potentiate channel activity by enhancing the coupling between Ca 2+ sensing via calmodulin and mechanical opening of the channel. Here, we review binding pocket studies that have provided structural insight into the mechanism of action for SK/IK channel positive modulators. These studies lay the foundation for structure-based drug discovery efforts that can identify novel SK/IK channel positive modulators

Theoretical Study of the anti-NCP Molecular Mechanism of Traditional Chinese Medicine Lianhua-Qingwen Formula (LQF)

Due to the good clinical efficacy in treating Novel Coronavirus Pneumonia (NCP) resulted from SARS-CoV-2, as the traditional Chinese medicine(TCM) prescription, Lianhua Qingwen Formula (LQF) was composed into the Diagnosis and Treatment Programs of 2019 New Coronavirus Pneumonia (from fourth to seventh editions) formulated by the National Health Commission of China. Aiming to prevent and treat viral influenza, LQF was patented from 2003 in China, and passed the Phase II clinical trial by FDA in the United States in 2015. However, the molecular mechanism of LQF anti SARS-CoV-2 pneumonia is still not clear. It is shown that the docking scores of three components in LQF including Rutin, Forsythoside E, and Hyperoside to main protease of SARS-CoV-2 are very large as -9.1, -9.0 and -8.7 kcal/mol, respectively, which are even better than those of Lopinavir at -7.3 kcal/mol. Importantly, the binding modes between active compounds and protein were verified via molecular dynamics (MD) simulation and calculation all the binding free energies at MM-PBSA level. Note that these donor-acceptor systems were stabilized by non-polar interactions including hydrogen bonds and hydrophobic interactions. At last, from the constructed component-target-pathway network, it is shown that the components in LQF are related important pathways to improve the human immunity such as T cell, B cell receptor signaling, natural killer cell mediated cytotoxicity, as well as anti inflammatory pathways including Fc epsilon RI, ErbB, MAPK signaling and so on. The present investigation represents the first report on the molecular mechanism of LQF as NCP inhibitor</div

Article Discovery of a Novel CCR5 Antagonist Lead Compound Through Fragment Assembly

molecule

Targeting the Small- and Intermediate-Conductance Ca2+-Activated Potassium Channels: The Drug-Binding Pocket at the Channel/Calmodulin Interface

The small- and intermediate-conductance Ca2+-activated potassium (SK/IK) channels play important roles in the regulation of excitable cells in both the central nervous and cardiovascular systems. Evidence from animal models has implicated SK/IK channels in neurological conditions such as ataxia and alcohol use disorders. Further, genome-wide association studies have suggested that cardiovascular abnormalities such as arrhythmias and hypertension are associated with single nucleotide polymorphisms that occur within the genes encoding the SK/IK channels. The Ca2+ sensitivity of the SK/IK channels stems from a constitutively bound Ca2+-binding protein: calmodulin. Small-molecule positive modulators of SK/IK channels have been developed over the past decade, and recent structural studies have revealed that the binding pocket of these positive modulators is located at the interface between the channel and calmodulin. SK/IK channel positive modulators can potentiate channel activity by enhancing the coupling between Ca2+ sensing via calmodulin and mechanical opening of the channel. Here, we review binding pocket studies that have provided structural insight into the mechanism of action for SK/IK channel positive modulators. These studies lay the foundation for structure-based drug discovery efforts that can identify novel SK/IK channel positive modulators. © 2014 S. Karger AG, Base

Structural insights into the potency of SK channel positive modulators.

Small-conductance Ca2+-activated K+ (SK) channels play essential roles in the regulation of cellular excitability and have been implicated in neurological and cardiovascular diseases through both animal model studies and human genetic association studies. Over the past two decades, positive modulators of SK channels such as NS309 and 1-EBIO have been developed. Our previous structural studies have identified the binding pocket of 1-EBIO and NS309 that is located at the interface between the channel and calmodulin. In this study, we took advantage of four compounds with potencies varying over three orders of magnitude, including 1-EBIO, NS309, SKS-11 (6-bromo-5-methyl-1H-indole-2,3-dione-3-oxime) and SKS-14 (7-fluoro-3-(hydroxyimino)indolin-2-one). A combination of x-ray crystallographic, computational and electrophysiological approaches was utilized to investigate the interactions between the positive modulators and their binding pocket. A strong trend exists between the interaction energy of the compounds within their binding site calculated from the crystal structures, and the potency of these compounds in potentiating the SK2 channel current determined by electrophysiological recordings. Our results further reveal that the difference in potency of the positive modulators in potentiating SK2 channel activity may be attributed primarily to specific electrostatic interactions between the modulators and their binding pocket