HISNAPI: a bioinformatic tool for dynamic hot spot analysis in nucleic acid-protein interface with a case study.

Protein-nucleic acid interactions play essential roles in many biological processes, such as transcription, replication and translation. In protein-nucleic acid interfaces, hotspot residues contribute the majority of binding affinity toward molecular recognition. Hotspot residues are commonly regarded as potential binding sites for compound molecules in drug design projects. The dynamic property is a considerable factor that affects the binding of ligands. Computational approaches have been developed to expedite the prediction of hotspot residues on protein-nucleic acid interfaces. However, existing approaches overlook hotspot dynamics, despite their essential role in protein function. Here, we report a web server named Hotspots In silico Scanning on Nucleic Acid and Protein Interface (HISNAPI) to analyze hotspot residue dynamics by integrating molecular dynamics simulation and one-step free energy perturbation. HISNAPI is capable of not only predicting the hotspot residues in protein-nucleic acid interfaces but also providing insights into their intensity and correlation of dynamic motion. Protein dynamics have been recognized as a vital factor that has an effect on the interaction specificity and affinity of the binding partners. We applied HISNAPI to the case of SARS-CoV-2 RNA-dependent RNA polymerase, a vital target of the antiviral drug for the treatment of coronavirus disease 2019. We identified the hotspot residues and characterized their dynamic behaviors, which might provide insight into the target site for antiviral drug design. The web server is freely available via a user-friendly web interface at http://chemyang.ccnu.edu.cn/ccb/server/HISNAPI/ and http://agroda.gzu.edu.cn:9999/ccb/server/HISNAPI/.

Bioinformatics toolbox for exploring protein phosphorylation network.

A clear systematic delineation of the interactions between phosphorylation sites on substrates and their effector kinases plays a fundamental role in revealing cellular activities, understanding signaling modulation mechanisms and proposing novel hypotheses. The emergence of bioinformatics tools contributes to studying phosphorylation network. Some of them feature the visualization of network, enabling more effective trace of the underlying biological problems in a clear and succinct way. In this review, we aimed to provide a toolbox for exploring phosphorylation network. We first systematically surveyed 19 tools that are available for exploring phosphorylation networks, and subsequently comparatively analyzed and summarized these tools to guide tool selection in terms of functionality, data sources, performance, network visualization and implementation, and finally briefly discussed the application cases of these tools. In different scenarios, the conclusion on the suitability of a tool for a specific user may vary. Nevertheless, easily accessible bioinformatics tools are proved to facilitate biological findings. Hopefully, this work might also assist non-specialists, students, as well as computational scientists who aim at developing novel tools in the field of phosphorylation modification.

Cloud 3D-QSAR: a web tool for the development of quantitative structure-activity relationship models in drug discovery.

Effective drug discovery contributes to the treatment of numerous diseases but is limited by high costs and long cycles. The Quantitative Structure-Activity Relationship (QSAR) method was introduced to evaluate the activity of a large number of compounds virtually, reducing the time and labor costs required for chemical synthesis and experimental determination. Hence, this method increases the efficiency of drug discovery. To meet the needs of researchers to utilize this technology, numerous QSAR-related web servers, such as Web-4D-QSAR and DPubChem, have been developed in recent years. However, none of the servers mentioned above can perform a complete QSAR modeling and supply activity prediction functions. We introduce Cloud 3D-QSAR by integrating the functions of molecular structure generation, alignment, molecular interaction field (MIF) computing and results analysis to provide a one-stop solution. We rigidly validated this server, and the activity prediction correlation was R2 = 0.934 in 834 test molecules. The sensitivity, specificity and accuracy were 86.9%, 94.5% and 91.5%, respectively, with AUC = 0.981, AUCPR = 0.971. The Cloud 3D-QSAR server may facilitate the development of good QSAR models in drug discovery. Our server is free and now available at http://chemyang.ccnu.edu.cn/ccb/server/cloud3dQSAR/ and http://agroda.gzu.edu.cn:9999/ccb/server/cloud3dQSAR/.

Role of Computer-Aided Drug Design in Drug Development.

The introduction of computational techniques to pharmaceutical chemistry and molecular biology in the 20th century has changed the way people develop drugs [...].

AIMMS suite: a web server dedicated for prediction of drug resistance on protein mutation.

Drug resistance is one of the most intractable issues for successful treatment in current clinical practice. Although many mutations contributing to drug resistance have been identified, the relationship between the mutations and the related pharmacological profile of drug candidates has yet to be fully elucidated, which is valuable both for the molecular dissection of drug resistance mechanisms and for suggestion of promising treatment strategies to counter resistant. Hence, effective prediction approach for estimating the sensitivity of mutations to agents is a new opportunity that counters drug resistance and creates a high interest in pharmaceutical research. However, this task is always hampered by limited known resistance training samples and accurately estimation of binding affinity. Upon this challenge, we successfully developed Auto In Silico Macromolecular Mutation Scanning (AIMMS), a web server for computer-aided de novo drug resistance prediction for any ligand-protein systems. AIMMS can qualitatively estimate the free energy consequences of any mutations through a fast mutagenesis scanning calculation based on a single molecular dynamics trajectory, which is differentiated with other web services by a statistical learning system. AIMMS suite is available at http://chemyang.ccnu.edu.cn/ccb/server/AIMMS/.

In Silico Methods for Identification of Potential Active Sites of Therapeutic Targets.

Target identification is an important step in drug discovery, and computer-aided drug target identification methods are attracting more attention compared with traditional drug target identification methods, which are time-consuming and costly. Computer-aided drug target identification methods can greatly reduce the searching scope of experimental targets and associated costs by identifying the diseases-related targets and their binding sites and evaluating the druggability of the predicted active sites for clinical trials. In this review, we introduce the principles of computer-based active site identification methods, including the identification of binding sites and assessment of druggability. We provide some guidelines for selecting methods for the identification of binding sites and assessment of druggability. In addition, we list the databases and tools commonly used with these methods, present examples of individual and combined applications, and compare the methods and tools. Finally, we discuss the challenges and limitations of binding site identification and druggability assessment at the current stage and provide some recommendations and future perspectives.

PIIMS Server: A Web Server for Mutation Hotspot Scanning at the Protein-Protein Interface.

Protein-protein interactions (PPIs) play vital roles in regulating biological processes, such as cellular and signaling pathways. Hotspots are certain residues located at protein-protein interfaces that contribute more in protein-protein binding than other residues. Research on the mutational effects of hotspots is important for understanding basic aspects of protein association. Hence, various computational tools have been developed to explore the impact of mutation hotspots, which will allow a better understanding of the forces that drive PPIs. However, tools that may provide comprehensive substitutions at hotspots are still rare. Hence, there is a strong need for a new free web server to explore mutational effects of hotspots. Herein we introduce a web server named PIIMS that integrates molecular dynamics simulation and one-step free energy perturbation. It contains two main computational functions: (1) computational alanine scanning analysis to identify hotspots and (2) full mutation scanning analysis to evaluate the effects of hotspot mutations. We rigidly validated its ability to predict binding free energy changes by using large and diverse datasets including 1,341 mutations from 50 PPIs with the correlation coefficient R = 0.75. The difference from the existing tools is that PIIMS can perform further evaluation of hotspot residues with regard to their different mutations. The PIIMS web server (accessible at http://chemyang.ccnu.edu.cn/ccb/server/PIIMS/index.php) is free and open to all users without login requirements.

Discovery of N-substituted-3-phenyl-1,6-naphthyridinone derivatives bearing quinoline moiety as selective type II c-Met kinase inhibitors against VEGFR-2.

New N-substituted-3-phenyl-1,6-naphthyridinone derivatives are designed and synthesized, based on structural modification of our previously reported compound 3. Extensive enzyme-based SAR studies and PK evaluation led to the discovery of compound 4r, with comparable c-Met potency to that of Cabozantinib and high VEGFR-2 selectivity, while Cabozantinib displayed no VEGFR-2 selectivity. More importantly, at oral doses of 45 mg/kg (Q.D.), compound 4r exhibits significant tumor growth inhibition (93%) in a U-87MG human gliobastoma xenograft model. The promising selectivity against VEGFR-2 and excellent tumor growth inhibition of compound 4r suggest that it could be used as a new lead molecule for further discovery of selective type II c-Met inhibitors.

Separation and detection of cyproconazole enantiomers and its stereospecific recognition with chiral stationary phase by high-performance liquid chromatography.

Cyproconazole, a chiral triazole fungicide, has been diffusely used and analyzed. The development of an effective analytical method for cyproconazole enantiomers can support their residual monitoring and risk assessment. In the present study, the absolute configuration of the cyproconazole enantiomers was confirmed by electronic circular dichroism and time-dependent density functional theory. The enantioseparation parameters were optimized by the response surface methodology using the Box-Behnken design on Lux Cellulose-2. The elution order of (2S,3R)-(+)-, (2S,3S)-(+)-, (2R,3S)-(-)-, and (2R,3R)-(-)-cyproconazole was simulated with molecular docking. The enantiomers were completely separated primarily via halogen bond and hydrogen bond interactions with the chiral stationary phases. The mean recoveries of the cyproconazole enantiomers in the four matrices were 71.8-102.0% with intraday relative standard deviations (RSDs) of 0.3-11.9% and interday RSDs of 0.9-10.6%. The results showed the chiral recognition mechanism clearly and confirmed that the method was accurate and convenient for the simultaneous detection of cyproconazole enantiomers in environmental and food matrices.

ACFIS: a web server for fragment-based drug discovery.

In order to foster innovation and improve the effectiveness of drug discovery, there is a considerable interest in exploring unknown 'chemical space' to identify new bioactive compounds with novel and diverse scaffolds. Hence, fragment-based drug discovery (FBDD) was developed rapidly due to its advanced expansive search for 'chemical space', which can lead to a higher hit rate and ligand efficiency (LE). However, computational screening of fragments is always hampered by the promiscuous binding model. In this study, we developed a new web server Auto Core Fragment in silico Screening (ACFIS). It includes three computational modules, PARA_GEN, CORE_GEN and CAND_GEN. ACFIS can generate core fragment structure from the active molecule using fragment deconstruction analysis and perform in silico screening by growing fragments to the junction of core fragment structure. An integrated energy calculation rapidly identifies which fragments fit the binding site of a protein. We constructed a simple interface to enable users to view top-ranking molecules in 2D and the binding mode in 3D for further experimental exploration. This makes the ACFIS a highly valuable tool for drug discovery. The ACFIS web server is free and open to all users at http://chemyang.ccnu.edu.cn/ccb/server/ACFIS/.

Fluorescence sensing of ADP over ATP and PPi in 100% aqueous solution.

An anthracene-bridged dinuclear zinc(ii)-dipicolylamine complex was found to show high selectivity for ADP with a significant fluorescence enhancement over ATP, PPi and other common analytes in 100% aqueous solution. This complex can be used for fluorescence detection of ADP in living cells and for monitoring the activity of kinases.

Coumarin-furo[2,3-d]pyrimidone hybrid molecules targeting human liver cancer cells: synthesis, anticancer effect, EGFR inhibition and molecular docking studies.

The design, synthesis and investigation of antitumor activities of some coumarin-furo[2,3-d]pyrimidone hybrid molecules are reported. In vitro, HepG2 cells were used to investigate the cytotoxicity of 6a-n and 10a-n. The results demonstrated that coupling a furopyrimidone scaffold with coumarin through a hydrazide linker can effectively improve their synergistic anticancer activity. The coumarin-furo[2,3-d]pyrimidone combination 10a exhibited significant inhibitory activity against HepG2 cells with IC50 = 7.72 ± 1.56 µM, which is better than those of gefitinib and sorafenib. It is worth mentioning that the coumarin-furo[2,3-d]pyrimidone combination 10a showed excellent inhibition of the EGFR enzymatic activity with IC50 = 1.53 µM and 90% inhibition at 10 µM concentration. In silico investigation predicts the possibility of direct binding between the new coumarin-furo[2,3-d]pyrimidone hybrid molecules and the EGFR. The results suggest that coumarin-furo[2,3-d]pyrimidone hybrid molecules are potential antitumor agents targeting human liver cancer cells.

Novel Coumarin-furo[2,3-d]pyrimidinone hybrid derivatives as anticancer agents: Synthesis, biological evaluation and molecular docking.

A series of coumarin-furo[2,3-d]pyrimidinone hybrid derivatives were synthesized, characterized by HR-MS, 1H NMR and 13C NMR. All synthesized compounds were evaluated for antiproliferative activities against hepatic carcinoma (HepG2) and cervical carcinoma (Hela) cell lines in vitro, and results shown that most of the compounds exhibited potent antitumor activity. Moreover, compound 3i, 8d and 8i were selected to induce apoptosis in HepG2 cells, and it displayed a significant concentration-dependent. Further, transwell migration assay was used to detect the most potent compound 8i, and the results revealed that 8i can significantly inhibit HepG2 cells migration and invasion. In addition, kinase activity assay showed compound 8i may be a multi-target inhibitor, which 8i has an inhibition rate of 40-20% on RON, ABL, GSK3α and so on ten different kinases at the concentration 1 µmol/L. At the same time, molecular docking studies revealed the possible binding modes of compounds 3i, 8d and 8i with kinase recepteur d'origine nantais (RON). A comparative molecular field analysis (CoMFA) model was established from 3D-QSAR study that guide us to a more bulkly and electro-positive Y group at the C-2 position of furo[2,3-d]pyrimidinone ring was preferable for the bioactivity improvement of our compounds. Our preliminary research indicated that the coumarin skeleton introducing to the furo[2,3-d]pyrimidine system had a significantly influence on the biological activities.

Comparison of the molecular mechanisms of Fuzi Lizhong Pill and Huangqin decoction in the treatment of the cold and heat syndromes of ulcerative colitis based on network pharmacology.

OBJECTIVE: The aim of this study was to illuminate the similarities and differences of two prescriptions as "cold" and "heat" drugs for treating ulcerative colitis (UC) with the simultaneous occurrence of heat and cold syndrome via network pharmacology. METHODS: (1) Active compounds of Fuzi-Lizhong Pill (FLP) and Huangqin Decoction (HQT) were retrieved from the TCMSP database, and their common active compounds were compared using the Venn diagram. (2) Potential proteins targeted to three sets of compounds either (i) shared by FLP and HQT, (ii) unique to FLP or (iii) unique to HQT were screened from the STP, STITCH and TCMSP databases, and three corresponding core compound sets were identified in Herb-Compound-Target (H-C-T) networks. (3) Targets related to UC were identified from the DisGeNET and GeneCards databases and compared with the FLP-HQT common targets to identify potential targets of FLP-HQT compounds related to UC. (4) Three potential target sets were imported into the STRING database for protein‒protein interaction (PPI) analysis, and three core target sets were defined. (5) The binding capabilities and interacting modes between core compounds and key targets were verified by molecular docking via Discovery Studio 2019 and molecular dynamics (MD) simulations via Amber 2018. (6) The target sets were enriched for KEGG pathways using the DAVID database. RESULTS: (1) FLP and HQT included 95 and 113 active compounds, respectively, with 46 common compounds, 49 FLP-specific compounds and 67 HQT-specific compounds. (2) 174 targets of FLP-HQT common compounds, 168 targets of FLP-specific compounds, and 369 targets of HQT-specific compounds were predicted from the STP, STITCH and TCMSP databases; six core compounds specific to FLP and HQT were screened in the FLP-specific and HQT-specific H-C-T networks, respectively. (3) 103 targets overlapped from the 174 predicted targets and the 4749 UC-related targets; two core compounds for FLP-HQT were identified from the FLP-HQT H-C-T network. (4) 103 FLP-HQT-UC common targets, 168 of FLP-specific targets and 369 of HQT-specific targets had shared core targets (AKT1, MAPK3, TNF, JUN and CASP3) based on the PPI network analysis. (5) Molecular docking demonstrated that naringenin, formononetin, luteolin, glycitein, quercetin, kaempferol and baicalein of FLP and HQT play a critical role in treating UC; meanwhile, MD simulations revealed the stability of protein‒ligand interactions. (6) The enriched pathways indicated that most targets were related to anti-inflammatory, immunomodulatory and other pathways. Compared with the pathways identified using traditional methods, FLP-specific pathways included the PPAR signaling pathway and the bile secretion pathway, and HQT-specific pathways included the vascular smooth muscle contraction pathway and the natural killer cell-mediated cytotoxicity pathway etc. CONCLUSION: In this study, we clarified the common mechanisms of FLP and HQT in treating UC and their specific mechanisms in treating cold and heat syndrome in UC through compound, target and pathway distinction and a literature comparison based on network pharmacology; these results provide a new perspective on the detailed mechanism of "multidrugs and single-disease" thought in traditional Chinese medicine.

In silico Methods for Identification of Potential Therapeutic Targets.

At the initial stage of drug discovery, identifying novel targets with maximal efficacy and minimal side effects can improve the success rate and portfolio value of drug discovery projects while simultaneously reducing cycle time and cost. However, harnessing the full potential of big data to narrow the range of plausible targets through existing computational methods remains a key issue in this field. This paper reviews two categories of in silico methods-comparative genomics and network-based methods-for finding potential therapeutic targets among cellular functions based on understanding their related biological processes. In addition to describing the principles, databases, software, and applications, we discuss some recent studies and prospects of the methods. While comparative genomics is mostly applied to infectious diseases, network-based methods can be applied to infectious and non-infectious diseases. Nonetheless, the methods often complement each other in their advantages and disadvantages. The information reported here guides toward improving the application of big data-driven computational methods for therapeutic target discovery.

Discovery of 3-pyrazolyl-substituted pyrazolo[1,5-a]pyrimidine derivatives as potent TRK inhibitors to overcome clinically acquired resistance.

Several secondary tropomyosin receptor kinase (TRK) mutations located in the solvent front, xDFG, and gatekeeper regions, are a common cause of clinical resistance. Mutations in the xDFG motif in particular limit sensitivity to second-generation TRK inhibitors, which represent an unmet clinical need. We designed a series of 3-pyrazolyl-substituted pyrazolo[1,5-a]pyrimidine derivatives toward these secondary mutations using ring-opening and scaffold-hopping strategies. Compound 5n was the most potent, with IC50 values of 2.3 nM, 0.4 nM, and 0.5 nM against TRKAG667C, TRKAF589L, and TRKAG595R, compared to selitrectinib with IC50 values of 12.6 nM, 5.8 nM, and 7.6 nM, respectively (approximately 5.4, 14.5, and 15.2-fold increases). Furthermore, 5n displayed favorable pharmacokinetic properties and satisfactory antitumor efficacy (tumor growth inhibition of 97% at 30 mg/kg and 73% at 100 mg/kg) in TRKAWT and TRKAG667C xenograft mouse models. Collectively, 5n is a promising TRK inhibitor lead compound for overcoming clinically acquired resistance to second-generation inhibitors, particularly for resistant tumors harboring the TRKAG667C mutation in the xDFG motif.

Green synthesis, structure optimization and biological evalution of Rhopaladins' analog 2-styryl-5-oxopyrrolidine-2- carboxamide RPDPRH on CaSki cells.

We have synthesized Rhopaladins' analog (2E,4E)-4-chlorobenzylidene-2-(4-chlorostyryl)-N-cyclohexyl-1-(4-fluorophenyl)-5-oxopyrrolidine-2-carboxamide (RPDPRH) via a highly facile, inexpensive and green approach and verified the structural superiority of compound RPDPRH through molecular docking. Moreover, we further detected the anti-proliferation, apoptosis and HPV E6/E7 effects of RPDPRH on CaSki cells. Finally, we confirmed that compared with the previous compound (E)-N-(tert-butyl)-2-(4-chlorobenzoyl)-4-(4-fluorobenzylidene)-1-isopropyl-5-oxopyrrolidine-2-carboxamide (RPDPB), RPDPRH could better inhibit proliferation, induce apoptosis, and down-regulate HPV E6/E7 mRNA expression on Caski cells. And preliminary RT-PCR experiments have demonstrated that RPDPRH also could affect the expression of Bcl-2, Bax and Caspase-3 mRNA in Caski cells. In summary, RPDPRH has potential as an effective agent against cervical cancer and will play an important role in our subsequent research.

Protocol for hit-to-lead optimization of compounds by auto in silico ligand directing evolution (AILDE) approach.

Hit-to-lead (H2L) optimization is crucial for drug design, which has become an increasing concern in medicinal chemistry. A virtual screening strategy of auto in silico ligand directing evolution (AILDE) has been developed to yield promising lead compounds rapidly and efficiently. The protocol includes instructions for fragment compound library construction, conformational sampling by molecular dynamics simulation, ligand modification by fragment growing, as well as the binding free energy prediction. For complete details on the use and execution of this protocol, please refer to Wu et al. (2020).

Discovery of Next-Generation Tropomyosin Receptor Kinase Inhibitors for Combating Multiple Resistance Associated with Protein Mutation.

Tropomyosin receptor kinase (TRK) inhibition is an effective therapeutic approach for treatment of a variety of cancers. Despite the use of first-generation TRK inhibitor (TRKI) larotrectinib (1) resulting in significant therapeutic response in patients, acquired resistance develops invariably. The emergence of secondary mutations occurring at the solvent-front, xDFG, and gatekeeper regions of TRK represents a common mechanism for acquired resistance. However, xDFG mutations remain insensitive to second-generation macrocyclic TRKIs selitrectinib (3) and repotrectinib (4) designed to overcome the resistance mediated by solvent-front and gatekeeper mutations. Here, we report the structure-based drug design and discovery of a next-generation TRKI. The structure-activity relationship studies culminated in the identification of a promising drug candidate 8 that showed excellent in vitro potency on a panel of TRK mutants, especially TRKAG667C in the xDFG motif, and improved in vivo efficacy than 1 and 3 in TRK wild-type and mutant fusion-driven tumor xenograft models, respectively.

Auto In Silico Ligand Directing Evolution to Facilitate the Rapid and Efficient Discovery of Drug Lead.

Motivated by the growing demand for reducing the chemical optimization burden of H2L, we developed auto in silico ligand directing evolution (AILDE, http://chemyang.ccnu.edu.cn/ccb/server/AILDE), an efficient and general approach for the rapid identification of drug leads in accessible chemical space. This computational strategy relies on minor chemical modifications on the scaffold of a hit compound, and it is primarily intended for identifying new lead compounds with minimal losses or, in some cases, even increases in ligand efficiency. We also described how AILDE greatly reduces the chemical optimization burden in the design of mesenchymal-epithelial transition factor (c-Met) kinase inhibitors. We only synthesized eight compounds and found highly efficient compound 5g, which showed an ∼1,000-fold improvement in in vitro activity compared with the hit compound. 5g also displayed excellent in vivo antitumor efficacy as a drug lead. We believe that AILDE may be applied to a large number of studies for rapid design and identification of drug leads.