An iterative compound screening contest method for identifying target protein inhibitors using the tyrosine-protein kinase Yes

Novel Glycomimetics Protect against Glycated Low-Density Lipoprotein-Induced Vascular Calcification In Vitro via Attenuation of the RAGE/ERK/CREB Pathway

Heparan sulphate (HS) can act as a co-receptor on the cell surface and alterations in this process underpin many pathological conditions. We have previously described the usefulness of mimics of HS (glycomimetics) in protection against β-glycerophosphate-induced vascular calcification and in the restoration of the functional capacity of diabetic endothelial colony-forming cells in vitro. This study aims to investigate whether our novel glycomimetic compounds can attenuate glycated low-density lipoprotein (g-LDL)-induced calcification by inhibiting RAGE signalling within the context of critical limb ischemia (CLI). We used an established osteogenic in vitro vascular smooth muscle cell (VSMC) model. Osteoprotegerin (OPG), sclerostin and glycation levels were all significantly increased in CLI serum compared to healthy controls, while the vascular calcification marker osteocalcin (OCN) was down-regulated in CLI patients vs. controls. Incubation with both CLI serum and g-LDL (10 µg/mL) significantly increased VSMC calcification vs. controls after 21 days, with CLI serum-induced calcification apparent after only 10 days. Glycomimetics (C2 and C3) significantly inhibited g-LDL and CLI serum-induced mineralisation, as shown by a reduction in alizarin red (AR) staining and alkaline phosphatase (ALP) activity. Furthermore, secretion of the osteogenic marker OCN was significantly reduced in VSMCs incubated with CLI serum in the presence of glycomimetics. Phosphorylation of cyclic AMP response element-binding protein (CREB) was significantly increased in g-LDL-treated cells vs. untreated controls, which was attenuated with glycomimetics. Blocking CREB activation with a pharmacological inhibitor 666-15 replicated the protective effects of glycomimetics, evidenced by elevated AR staining. In silico molecular docking simulations revealed the binding affinity of the glycomimetics C2 and C3 with the V domain of RAGE. In conclusion, these findings demonstrate that novel glycomimetics, C2 and C3 have potent anti-calcification properties in vitro, inhibiting both g-LDL and CLI serum-induced VSMC mineralisation via the inhibition of LDLR, RAGE, CREB and subsequent expression of the downstream osteogenic markers, ALP and OCN.

Virtual Screening Process: A Guide in Modern Drug Designing

Due to its capacity to drastically cut the cost and time necessary for experimental screening of compounds, virtual screening (VS) has grown to be a crucial component of drug discovery and development. VS is a computational method used in drug design to identify potential drugs from enormous libraries of chemicals. This approach makes use of molecular modeling and docking simulations to assess the small molecule's ability to bind to the desired protein. Virtual screening has a bright future, as high computational power and modern techniques are likely to further enhance the accuracy and speed of the process.

Discovery of a Hidden Trypanosoma cruzi Spermidine Synthase Binding Site and Inhibitors through In Silico, In Vitro, and X-ray Crystallography

In drug discovery research, the selection of promising binding sites and understanding the binding mode of compounds are crucial fundamental studies. The current understanding of the proteins-ligand binding model extends beyond the simple lock and key model to include the induced-fit model, which alters the conformation to match the shape of the ligand, and the pre-existing equilibrium model, selectively binding structures with high binding affinity from a diverse ensemble of proteins. Although methods for detecting target protein binding sites and virtual screening techniques using docking simulation are well-established, with numerous studies reported, they only consider a very limited number of structures in the diverse ensemble of proteins, as these methods are applied to a single structure. Molecular dynamics (MD) simulation is a method for predicting protein dynamics and can detect potential ensembles of protein binding sites and hidden sites unobservable in a single-point structure. In this study, to demonstrate the utility of virtual screening with protein dynamics, MD simulations were performed on Trypanosoma cruzi spermidine synthase to obtain an ensemble of dominant binding sites with a high probability of existence. The structure of the binding site obtained through MD simulation revealed pockets in addition to the active site that was present in the initial structure. Using the obtained binding site structures, virtual screening of 4.8 million compounds by docking simulation, in vitro assays, and X-ray analysis was conducted, successfully identifying two hit compounds.

Quantitative comparison of protein-protein interaction interface using physicochemical feature-based descriptors of surface patches

Driving mechanisms of many biological functions in a cell include physical interactions of proteins. As protein-protein interactions (PPIs) are also important in disease development, protein-protein interactions are highlighted in the pharmaceutical industry as possible therapeutic targets in recent years. To understand the variety of protein-protein interactions in a proteome, it is essential to establish a method that can identify similarity and dissimilarity between protein-protein interactions for inferring the binding of similar molecules, including drugs and other proteins. In this study, we developed a novel method, protein-protein interaction-Surfer, which compares and quantifies similarity of local surface regions of protein-protein interactions. protein-protein interaction-Surfer represents a protein-protein interaction surface with overlapping surface patches, each of which is described with a three-dimensional Zernike descriptor (3DZD), a compact mathematical representation of 3D function. 3DZD captures both the 3D shape and physicochemical properties of the protein surface. The performance of protein-protein interaction-Surfer was benchmarked on datasets of protein-protein interactions, where we were able to show that protein-protein interaction-Surfer finds similar potential drug binding regions that do not share sequence and structure similarity. protein-protein interaction-Surfer is available at https://kiharalab.org/ppi-surfer.

Inverse Mixed-Solvent Molecular Dynamics for Visualization of the Residue Interaction Profile of Molecular Probes

To ensure efficiency in discovery and development, the application of computational technology is essential. Although virtual screening techniques are widely applied in the early stages of drug discovery research, the computational methods used in lead optimization to improve activity and reduce the toxicity of compounds are still evolving. In this study, we propose a method to construct the residue interaction profile of the chemical structure used in the lead optimization by performing “inverse” mixed-solvent molecular dynamics (MSMD) simulation. Contrary to constructing a protein-based, atom interaction profile, we constructed a probe-based, protein residue interaction profile using MSMD trajectories. It provides us the profile of the preferred protein environments of probes without co-crystallized structures. We assessed the method using three probes: benzamidine, catechol, and benzene. As a result, the residue interaction profile of each probe obtained by MSMD was a reasonable physicochemical description of the general non-covalent interaction. Moreover, comparison with the X-ray structure containing each probe as a ligand shows that the map of the interaction profile matches the arrangement of amino acid residues in the X-ray structure.

Screening for Inhibitors of Main Protease in SARS-CoV-2: In Silico and In Vitro Approach Avoiding Peptidyl Secondary Amides

In addition to vaccines, antiviral drugs are essential for suppressing COVID-19. Although several inhibitor candidates were reported for SARS-CoV-2 main protease, most are highly polar peptidomimetics with poor oral bioavailability and cell membrane permeability. Here, we conducted structure-based virtual screening and in vitro assays to obtain hit compounds belonging to a new chemical space, excluding peptidyl secondary amides. In total, 180 compounds were subjected to the primary assay at 20 μM, and nine compounds with inhibition rates of >5% were obtained. The IC50 of six compounds was determined in dose-response experiments, with the values on the order of 10-4 M. Although nitro groups were enriched in the substructure of the hit compounds, they did not significantly contribute to the binding interaction in the predicted docking poses. Physicochemical properties prediction showed good oral absorption. These new scaffolds are promising candidates for future optimization.

MERMAID: an open source automated hit-to-lead method based on deep reinforcement learning

The hit-to-lead process makes the physicochemical properties of the hit molecules that show the desired type of activity obtained in the screening assay more drug-like. Deep learning-based molecular generative models are expected to contribute to the hit-to-lead process. The simplified molecular input line entry system (SMILES), which is a string of alphanumeric characters representing the chemical structure of a molecule, is one of the most commonly used representations of molecules, and molecular generative models based on SMILES have achieved significant success. However, in contrast to molecular graphs, during the process of generation, SMILES are not considered as valid SMILES. Further, it is quite difficult to generate molecules starting from a certain molecule, thus making it difficult to apply SMILES to the hit-to-lead process. In this study, we have developed a SMILES-based generative model that can be generated starting from a certain molecule. This method generates partial SMILES and inserts it into the original SMILES using Monte Carlo Tree Search and a Recurrent Neural Network. We validated our method using a molecule dataset obtained from the ZINC database and successfully generated molecules that were both well optimized for the objectives of the quantitative estimate of drug-likeness (QED) and penalized octanol-water partition coefficient (PLogP) optimization. The source code is available at https://github.com/sekijima-lab/mermaid.

Vitexin isolated from Prosopis cineraria leaves induce apoptosis in K-562 leukemia cells via inhibition of the BCR-ABL-Ras-Raf pathway

OBJECTIVES

Leukemia is one of the severe cancer types all around the globe. Even though some chemotherapeutic drugs are available for treating leukemia, they have various side effects. As an alternative approach, herbal drugs are focused on current research to overcome leukemia. The present work was conducted to investigate the antileukemic mechanism of active phytochemical vitexin, which was isolated from ethno-medicine (Prosopis cineraria leaf) used by traditional healers of West Bengal, India.

METHODS

Antiproliferative mechanisms of selected phyto-compound against K-562 cells were evaluated using cellular uptake, morphological changes, DNA fragmentation, mitochondrial membrane potential and signaling pathways analysis.

KEY FINDINGS

Vitexin exhibited cytotoxicity by reducing mitochondrial membrane potential (32.40%) and causing DNA fragmentation (84.15%). The western blotting study indicated inhibition of cell survival proteins (BCR, ABL, H-RAS, N-RAS, K-RAS and RAF) and expression of apoptotic proteins (p38, BAX and caspase-9) in leukemia cells upon treatment with vitexin.

CONCLUSIONS

Based on the results, presently investigated phyto-compound vitexin could be considered for developing safe and natural drugs to treat leukemia after conducting suitable preclinical and clinical trials.

Computer aided drug discovery review for infectious diseases with case study of anti-Chagas project

Neglected tropical diseases (NTDs) are parasitic and bacterial infections that are widespread, especially in the tropics, and cause health problems for about one billion people over 149 countries worldwide. However, in terms of therapeutic agents, for example, nifurtimox and benznidazole were developed in the 1960s to treat Chagas disease, but new drugs are desirable because of their side effects. Drug discovery takes 12 to 14 years and costs $2.6 billon dollars, and hence, computer aided drug discovery (CADD) technology is expected to reduce the time and cost. This paper describes our methods and results based on CADD, mainly for NTDs. An overview of databases, molecular simulation and pharmacophore modeling, contest-based drug discovery, and machine learning and their results are presented herein.

Identification of key interactions between SARS-CoV-2 main protease and inhibitor drug candidates

The number of cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection (COVID-19) has reached over 114,000. SARS-CoV-2 caused a pandemic in Wuhan, China, in December 2019 and is rapidly spreading globally. It has been reported that peptide-like anti-HIV-1 drugs are effective against SARS-CoV Main protease (Mpro). Due to the close phylogenetic relationship between SARS-CoV and SARS-CoV-2, their main proteases share many structural and functional features. Thus, these drugs are also regarded as potential drug candidates targeting SARS-CoV-2 Mpro. However, the mechanism of action of SARS-CoV-2 Mpro at the atomic-level is unknown. In the present study, we revealed key interactions between SARS-CoV-2 Mpro and three drug candidates by performing pharmacophore modeling and 1 μs molecular dynamics (MD) simulations. His41, Gly143, and Glu166 formed interactions with the functional groups that were common among peptide-like inhibitors in all MD simulations. These interactions are important targets for potential drugs against SARS-CoV-2 Mpro.

Integrating Ligand and Target-Driven Based Virtual Screening Approaches With in vitro Human Cell Line Models and Time-Resolved Fluorescence Resonance Energy Transfer Assay to Identify Novel Hit Compounds Against BCL-2

Antiapoptotic members of B-cell leukemia/lymphoma-2 (BCL-2) family proteins are one of the overexpressed proteins in cancer cells that are oncogenic targets. As such, targeting of BCL-2 family proteins raises hopes for new therapeutic discoveries. Thus, we used multistep screening and filtering approaches that combine structure and ligand-based drug design to identify new, effective BCL-2 inhibitors from a small molecule database (Specs SC), which includes more than 210,000 compounds. This database is first filtered based on binary “cancer-QSAR” model constructed with 886 training and 167 test set compounds and common 26 toxicity quantitative structure-activity relationships (QSAR) models. Predicted non-toxic compounds are considered for target-driven studies. Here, we applied two different approaches to filter and select hit compounds for further in vitro biological assays and human cell line experiments. In the first approach, a molecular docking and filtering approach is used to rank compounds based on their docking scores and only a few top-ranked molecules are selected for further long (100-ns) molecular dynamics (MD) simulations and in vitro tests. While docking algorithms are promising in predicting binding poses, they can be less prone to precisely predict ranking of compounds leading to decrease in the success rate of in silico studies. Hence, in the second approach, top-docking poses of each compound filtered through QSAR studies are subjected to initially short (1 ns) MD simulations and their binding energies are calculated via molecular mechanics generalized Born surface area (MM/GBSA) method. Then, the compounds are ranked based on their average MM/GBSA energy values to select hit molecules for further long MD simulations and in vitro studies. Additionally, we have applied text-mining approaches to identify molecules that contain “indol” phrase as many of the approved drugs contain indole and indol derivatives. Around 2700 compounds are filtered based on “cancer-QSAR” model and are then docked into BCL-2. Short MD simulations are performed for the top-docking poses for each compound in complex with BCL-2. The complexes are again ranked based on their MM/GBSA values to select hit molecules for further long MD simulations and in vitro studies. In total, seven molecules are subjected to biological activity tests in various human cancer cell lines as well as Time-Resolved Fluorescence Resonance Energy Transfer (TR-FRET) assay. Inhibitory concentrations are evaluated, and biological activities and apoptotic potentials are assessed by cell culture studies. Four molecules are found to be limiting the proliferation capacity of cancer cells while increasing the apoptotic cell fractions.

Computational approaches for identifying potential inhibitors on targeting protein interactions in drug discovery

In the era of big data, the interplay of artificial and human intelligence is the demanding job to address the concerns involving exchange of decisions between both sides. Drug discovery is one of the key sources of the big data, which involves synergy among various computational methods to achieve a clinical success. Rightful acquisition, mining and analysis of the data related to ligand and targets are crucial to accomplish reliable outcomes in the entire process. Novel designing and screening tactics are necessary to substantiate a potent and efficient lead compounds. Such methods are emphasized and portrayed in the current review targeting protein-ligand and protein-protein interactions involved in various diseases with potential applications.

A prospective compound screening contest identified broader inhibitors for Sirtuin 1

Potential inhibitors of a target biomolecule, NAD-dependent deacetylase Sirtuin 1, were identified by a contest-based approach, in which participants were asked to propose a prioritized list of 400 compounds from a designated compound library containing 2.5 million compounds using in silico methods and scoring. Our aim was to identify target enzyme inhibitors and to benchmark computer-aided drug discovery methods under the same experimental conditions. Collecting compound lists derived from various methods is advantageous for aggregating compounds with structurally diversified properties compared with the use of a single method. The inhibitory action on Sirtuin 1 of approximately half of the proposed compounds was experimentally accessed. Ultimately, seven structurally diverse compounds were identified.

Molecular Dynamics Simulation reveals the mechanism by which the Influenza Cap-dependent Endonuclease acquires resistance against Baloxavir marboxil

Baloxavir marboxil (BXM), an antiviral drug for influenza virus, inhibits RNA replication by binding to RNA replication cap-dependent endonuclease (CEN) of influenza A and B viruses. Although this drug was only approved by the FDA in October 2018, drug resistant viruses have already been detected from clinical trials owing to an I38 mutation of CEN. To investigate the reduction of drug sensitivity by the I38 mutant variants, we performed a molecular dynamics (MD) simulation on the CEN-BXM complex structure to analyze variations in the mode of interaction. Our simulation results suggest that the side chain methyl group of I38 in CEN engages in a CH-pi interaction with the aromatic ring of BXM. This interaction is abolished in various I38 mutant variants. Moreover, MD simulation on various mutation models and binding free energy prediction by MM/GBSA method suggest that the I38 mutation precludes any interaction with the aromatic ring of BXA and thereby reduces BXA sensitivity.

Improved Method of Structure-Based Virtual Screening via Interaction-Energy-Based Learning

Virtual screening is a promising method for obtaining novel hit compounds in drug discovery. It aims to enrich potentially active compounds from a large chemical library for further biological experiments. However, the accuracy of current virtual screening methods is insufficient. In this study, we develop a new virtual screening method named Similarity of Interaction Energy VEctor Score (SIEVE-Score), in which protein-ligand interaction energies are extracted to represent docking poses for machine learning. SIEVE-Score offers substantial improvements compared to other state-of-the-art virtual screening methods, namely, other machine-learning-based scoring functions, interaction fingerprints, and docking software, for the enrichment factor 1% results on the Directory of Useful Decoys, Enhanced (DUD-E). The screening results are also human-interpretable in the form of important interactions for distinguishing between active and inactive compounds. The source code is available at https://github.com/sekijima-lab/SIEVE-Score .

Homology Model and Docking-Based Virtual Screening for Ligands of Human Dyskerin as New Inhibitors of Telomerase for Cancer Treatment

Immortality is one of the main features of cancer cells. Tumor cells have an unlimited replicative potential, principally due to the holoenzyme telomerase. Telomerase is composed mainly by dyskerin (DKC1), a catalytic retrotranscriptase (hTERT) and an RNA template (hTR). The aim of this work is to develop new inhibitors of telomerase, selecting the interaction between hTR⁻DKC1 as a target. We designed two models of the human protein DKC1: homology and ab initio. These models were evaluated by different procedures, revealing that the homology model parameters were the most accurate. We selected two hydrophobic pockets contained in the PUA (pseudouridine synthase and archaeosine transglycosylase) domain, using structural and stability analysis. We carried out a docking-based virtual screen on these pockets, using the reported mutation K314 as the center of the docking. The hDKC1 model was tested against a library of 450,000 drug-like molecules. We selected the first 10 molecules that showed the highest affinity values to test their inhibitory activity on the cell line MDA MB 231 (Monroe Dunaway Anderson Metastasis Breast cancer 231), obtaining three compounds that showed inhibitory effect. These results allowed us to validate our design and set the basis to continue with the study of telomerase inhibitors for cancer treatment.

QEX: target-specific druglikeness filter enhances ligand-based virtual screening

Druglikeness is a useful concept for screening drug candidate compounds. We developed QEX, which is a new druglikeness index specific to individual targets. QEX is an improvement of the quantitative estimate of druglikeness (QED) method, which is a popular quantitative evaluation method of druglikeness proposed by Bickerton et al. QEX models the physicochemical properties of compounds that act on each target protein based on the concept of QED modeling physicochemical properties from information on US Food and Drug Administration-approved drugs. The result of the evaluation of PubChem assay data revealed that QEX showed better performance than the original QED did (the area under the curve value of the receiver operating characteristic curve improved by 0.069-0.236). We also present the c-Src inhibitor filtering results of the QEX constructed using Src family kinase inhibitors as a case study. QEX distinguished the inhibitors and non-inhibitors better than QED did. QEX works efficiently even when datasets of inactive compounds are unavailable. If both active and inactive compounds are present, QEX can be used as an initial filter to enhance the screening ability of conventional ligand-based virtual screenings.

Compound property enhancement by virtual compound synthesis

During drug discovery, drug candidates are narrowed down over several steps to develop pharmaceutical products. The theoretical chemical space in such steps is estimated to be [Formula: see text]. To cover that space, extensive virtual compound libraries have been developed; however, the compilation of extensive libraries comes at large computational cost. Thus, to reduce the computational cost, researchers have constructed custom-made virtual compound libraries that focus on target diseases. In this study, we develop a system that generates virtual compound libraries from input compounds. When a user inputs a compound, the system recursively applies virtual synthetic reaction rules to the compound to improve its properties. The synthetic pathway can also be traced by the user because the reaction rules in this system are based on real organic synthesis reactions. This system has useful functions for effective drug design, such as structural preservation, allowing the substructures necessary for potency to be maintained. In this paper, to confirm the effect of directional reaction sets, we applied the reaction sets to 100 compounds. Moreover, to confirm that the system can reproduce real synthetic pathways, the synthetic pathways of Ibuprofen and Ofloxacin were explored by inputting isobutyl benzene and 7,8-difluoro-2,3-dihydro-3-methyl-4H-benzoxazine. This application is available at the following URL: http://enh.sekijima-lab.org .

Exploring the selectivity of inhibitor complexes with Bcl-2 and Bcl-XL: A molecular dynamics simulation approach

B-cell lymphoma 2 (Bcl-2) family proteins are potential drug targets in cancer and have a relatively flat and flexible binding site. ABT-199 is one of the most promising selective Bcl-2 inhibitors, and A-1155463 selectively inhibits Bcl-XL. Although the amino acid sequences of the binding sites of these two inhibitors are similar, the inhibitors selectively bind the target protein. In order to determine the origin of the selectivity of these inhibitors, we conducted molecular dynamics simulations using protein-inhibitor modeling. We confirmed that ASP103 of Bcl-2 is a key residue and that hydrogen bonding between ASP103 and ABT-199 confers the Bcl-2 selectivity of this inhibitor. For Bcl-XL selectivity, the secondary structure of α-helix 3 is a key factor. PHE105, SER106, and LEU108 in the loose α-helix 3 interact with A-1155463 to confer Bcl-XL selectivity. These findings provide important insights into the molecular mechanisms of selective inhibitors of Bcl-2 family proteins.