Combined strategies in structure-based virtual screening

A head-to-head comparison of MM/PBSA and MM/GBSA in predicting binding affinities for the CB1 cannabinoid ligands

CONTEXT

The substantial increase in the number of active and inactive-state CB1 receptor experimental structures has provided opportunities for CB1 drug discovery using various structure-based drug design methods, including the popular end-point methods for predicting binding free energies-Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA). In this study, we have therefore evaluated the performance of MM/PBSA and MM/GBSA in calculating binding free energies for CB1 receptor. Additionally, with both MM/PBSA and MM/GBSA being known for their highly individualized performance, we have evaluated the effects of various simulation parameters including the use of energy minimized structures, choice of solute dielectric constant, inclusion of entropy, and the effects of the five GB models. Generally, MM/GBSA provided higher correlations than MM/PBSA (rMM/GBSA = 0.433 - 0.652 vs. rMM/PBSA = 0.100 - 0.486) regardless of the simulation parameters, while also offering faster calculations. Improved correlations were observed with the use of molecular dynamics ensembles compared with energy minimized structures and larger solute dielectric constants. Incorporation of entropic terms led to unfavorable results for both MM/PBSA and MM/GBSA for a majority of the dataset, while the evaluation of the various GB models exerted a varying effect on both the datasets. The findings obtained in this study demonstrate the utility of MM/PBSA and MM/GBSA in predicting binding free energies for the CB1 receptor, hence providing a useful benchmark for their applicability in the endocannabinoid system as well as other G protein-coupled receptors.

METHODS

The study utilized the docked dataset (Induced Fit Docking with Glide XP scoring function) from Loo et al., consisting of 46 ligands-23 agonists and 23 antagonists. The equilibrated structures from Loo et al. were subjected to 30 ns production simulations using GROMACS 2018 at 300 K and 1 atm with the velocity rescaling thermostat and the Parinello-Rahman barostat. AMBER ff99SB*-ILDN was used for the proteins, General Amber Force Field (GAFF) was used for the ligands, and Slipids parameters were used for lipids. MM/PBSA and MM/GBSA binding free energies were then calculated using gmx_MMPBSA. The solute dielectric constant was varied between 1, 2, and 4 to study the effect of different solute dielectric constants on the performance of MM/PB(GB)SA. The effect of entropy on MM/PB(GB)SA binding free energies was evaluated using the interaction entropy module implemented in gmx_MMPBSA. Five GB models, GBHCT, GBOBC1, GBOBC2, GBNeck, and GBNeck2, were evaluated to study the effect of the choice of GB models in the performance of MM/GBSA. Pearson correlation coefficients were used to measure the correlation between experimental and predicted binding free energies.

Research on Artificial-Intelligence-Assisted Medicine: A Survey on Medical Artificial Intelligence

With the improvement of economic conditions and the increase in living standards, people's attention in regard to health is also continuously increasing. They are beginning to place their hopes on machines, expecting artificial intelligence (AI) to provide a more humanized medical environment and personalized services, thus greatly expanding the supply and bridging the gap between resource supply and demand. With the development of IoT technology, the arrival of the 5G and 6G communication era, and the enhancement of computing capabilities in particular, the development and application of AI-assisted healthcare have been further promoted. Currently, research on and the application of artificial intelligence in the field of medical assistance are continuously deepening and expanding. AI holds immense economic value and has many potential applications in regard to medical institutions, patients, and healthcare professionals. It has the ability to enhance medical efficiency, reduce healthcare costs, improve the quality of healthcare services, and provide a more intelligent and humanized service experience for healthcare professionals and patients. This study elaborates on AI development history and development timelines in the medical field, types of AI technologies in healthcare informatics, the application of AI in the medical field, and opportunities and challenges of AI in the field of medicine. The combination of healthcare and artificial intelligence has a profound impact on human life, improving human health levels and quality of life and changing human lifestyles.

TransGEM: a molecule generation model based on Transformer with gene expression data

MOTIVATION

It is difficult to generate new molecules with desirable bioactivity through ligand-based de novo drug design, and receptor-based de novo drug design is constrained by disease target information availability. The combination of artificial intelligence and phenotype-based de novo drug design can generate new bioactive molecules, independent from disease target information. Gene expression profiles can be used to characterize biological phenotypes. The Transformer model can be utilized to capture the associations between gene expression profiles and molecular structures due to its remarkable ability in processing contextual information.

RESULTS

We propose TransGEM (Transformer-based model from gene expression to molecules), which is a phenotype-based de novo drug design model. A specialized gene expression encoder is used to embed gene expression difference values between diseased cell lines and their corresponding normal tissue cells into TransGEM model. The results demonstrate that the TransGEM model can generate molecules with desirable evaluation metrics and property distributions. Case studies illustrate that TransGEM model can generate structurally novel molecules with good binding affinity to disease target proteins. The majority of genes with high attention scores obtained from TransGEM model are associated with the onset of the disease, indicating the potential of these genes as disease targets. Therefore, this study provides a new paradigm for de novo drug design, and it will promote phenotype-based drug discovery.

AVAILABILITY AND IMPLEMENTATION

The code is available at https://github.com/hzauzqy/TransGEM.

Feeling the Strain: Quantifying Ligand Deformation in Photosynthesis

Structural distortion of protein-bound ligands can play a critical role in enzyme function by tuning the electronic and chemical properties of the ligand molecule. However, quantifying these effects is difficult due to the limited resolution of protein structures and the difficulty of generating accurate structural restraints for nonprotein ligands. Here, we seek to quantify these effects through a statistical analysis of ligand distortion in chlorophyll proteins (CP), where ring deformation is thought to play a role in energy and electron transfer. To assess the accuracy of ring-deformation estimates from available structural data, we take advantage of the C2 symmetry of photosystem II (PSII), comparing ring-deformation estimates for equivalent sites both within and between 113 distinct X-ray and cryogenic electron microscopy PSII structures. Significantly, we find that several deformation modes exhibit considerable variability in predictions, even for equivalent monomers, down to a 2 Å resolution, to an extent that probably prevents their utilization in optical calculations. We further find that refinement restraints play a critical role in determining deformation values to resolution as low as 2 Å. However, for those modes that are well-resolved in the structural data, ring deformation in PSII is strongly conserved across all species tested from cyanobacteria to algae. These results highlight both the opportunities and limitations inherent in structure-based analyses of the bioenergetic and optical properties of CPs and other protein-ligand complexes.

Inhibitor design for TMPRSS2: insights from computational analysis of its backbone hydrogen bonds using a simple descriptor

Transmembrane protease serine 2 (TMPRSS2) is an important drug target due to its role in the infection mechanism of coronaviruses including SARS-CoV-2. Current understanding regarding the molecular mechanisms of known inhibitors and insights required for inhibitor design are limited. This study investigates the effect of inhibitor binding on the intramolecular backbone hydrogen bonds (BHBs) of TMPRSS2 using the concept of hydrogen bond wrapping, which is the phenomenon of stabilization of a hydrogen bond in a solvent environment as a result of being surrounded by non-polar groups. A molecular descriptor which quantifies the extent of wrapping around BHBs is introduced for this. First, virtual screening for TMPRSS2 inhibitors is performed by molecular docking using the program DOCK 6 with a Generalized Born surface area (GBSA) scoring function. The docking results are then analyzed using this descriptor and its relationship to the solvent-accessible surface area term ΔGsa of the GBSA score is demonstrated with machine learning regression and principal component analysis. The effect of binding of the inhibitors camostat, nafamostat, and 4-guanidinobenzoic acid (GBA) on the wrapping of important BHBs in TMPRSS2 is also studied using molecular dynamics. For BHBs with a large increase in wrapping groups due to these inhibitors, the radial distribution function of water revealed that certain residues involved in these BHBs, like Gln438, Asp440, and Ser441, undergo preferential desolvation. The findings offer valuable insights into the mechanisms of these inhibitors and may prove useful in the design of new inhibitors.

Discovery of 4-aminophenylacetamide derivatives as intestine-specific farnesoid X receptor antagonists for the potential treatment of nonalcoholic steatohepatitis

Farnesoid X receptor (FXR) plays a key role in bile acid homeostasis, inflammation, fibrosis, lipid and glucose metabolism and is emerging as a promising therapeutic target for nonalcoholic steatohepatitis (NASH). Emerging evidence suggested that intestine-specific FXR antagonists exhibited remarkable metabolic improvements and slowed NASH progression. In this study, we discovered several potent FXR antagonists using a multistage ligand- and structure-based virtual screening approach. Notably, compound V023-9340, which possesses a 4-aminophenylacetamide scaffold, emerged as the most potent FXR antagonist with an IC50 value of 4.27 μM. In vivo, V023-9340 demonstrated selective accumulation in the intestine, substantially ameliorating high-fat diet (HFD)-induced NASH in mice by mitigating hepatic steatosis and inflammation. Mechanistic studies revealed that V023-9340 strongly inhibited intestinal FXR while concurrently feedback-activated hepatic FXR. Further structure-activity relationship optimization employing V023-9340 has resulted in the synthesis of a more efficacious compound V02-8 with an IC50 value of 0.89 μM, which exhibited a 4.8-fold increase in FXR antagonistic activity compared to V023-9340. In summary, 4-aminophenylacetamide derivative V023-9340 represented a novel intestine-specific FXR antagonist and showed improved effects against HFD-induced NASH in mice, which may serve as a promising lead in discovering potential therapeutic drugs for NASH treatment.

Proposition of In silico Pharmacophore Models for Malaria: A Review

In the field of medicinal chemistry, the concept of pharmacophore refers to the specific region of a molecule that possesses essential structural and chemical characteristics for binding to a receptor and eliciting biological activity. Understanding the pharmacophore is crucial for drug research and development, as it allows the design of new drugs. Malaria, a widespread disease, is commonly treated with chloroquine and artemisinin, but the emergence of parasite resistance limits their effectiveness. This study aims to explore computer simulations to discover a specific pharmacophore for Malaria, providing new alternatives for its treatment. A literature review was conducted, encompassing articles proposing a pharmacophore for Malaria, gathered from the "Web of Science" database, with a focus on recent publications to ensure up-to-date analysis. The selected articles employed diverse methods, including ligand-based and structurebased approaches, integrating molecular structure and biological activity data to yield comprehensive analyses. Affinity evaluation between the proposed pharmacophore and the target receptor involved calculating free energy to quantify their interaction. Multiple linear regression was commonly utilized, though it is sensitive to multicollinearity issues. Another recurrent methodology was the use of the Schrödinger package, employing tools such as the Phase module and the OPLS force field for interaction analysis. Pharmacophore model proposition allows threedimensional representations guiding the synthesis and design of new biologically active compounds, offering a promising avenue for discovering therapeutic agents to combat Malaria.

A novel in silico scaffold-hopping method for drug repositioning in rare and intractable diseases

In the field of rare and intractable diseases, new drug development is difficult and drug repositioning (DR) is a key method to improve this situation. In this study, we present a new method for finding DR candidates utilizing virtual screening, which integrates amino acid interaction mapping into scaffold-hopping (AI-AAM). At first, we used a spleen associated tyrosine kinase inhibitor as a reference to evaluate the technique, and succeeded in scaffold-hopping maintaining the pharmacological activity. Then we applied this method to five drugs and obtained 144 compounds with diverse structures. Among these, 31 compounds were known to target the same proteins as their reference compounds and 113 compounds were known to target different proteins. We found that AI-AAM dominantly selected functionally similar compounds; thus, these selected compounds may represent improved alternatives to their reference compounds. Moreover, the latter compounds were presumed to bind to the targets of their references as well. This new "compound-target" information provided DR candidates that could be utilized for future drug development.

High-throughput virtual screening approach of natural compounds as target inhibitors of plasmepsin-II

Plasmepsin II is a key enzyme in the life cycle of the Plasmodium falciparum parasite responsible for malaria, a disease that is causing deaths on a worldwide scale. Recently, plasmepsin II enzyme has gained much importance as an attractive drug target for the investigation of antimalarial drugs. In this sense, structure-based virtual screening have been utilized as tools in the process of discovering novel natural compounds based on quinoline as potential plasmepsin II inhibitors. Among the 58 quinoline derivatives isolated from different plants was screened by utilizing docking molecular, ADMET approaches, molecular dynamics simulation and MM-PBSA binding free energy. The first step in this work is building the 3 D structures of the plasmepsin II enzyme by using the SWISS-MODEL software. The optimized structures were subjected to virtual screening by Autodock Vina, an entity implicated in PyRx software. 21 were selected based on their binding affinity. The binding modes and interactions of the top-21 selected compounds were evaluated using AutoDock 4.2. Then, the pharmacokinetic proprieties and toxicity of these compounds were evaluated using ADMET analysis. Ten compounds were predicted to have ADMET characteristics with no side effects. Compounds M49 and M53 were found to be potential inhibitors. The stability of the selected two compounds was confirmed by MD simulation and MM/PBSA calculation during 200 ns. This study can be used to predict and to design new antimalarial drugs.Communicated by Ramaswamy H. Sarma.

Ligand binding affinity prediction with fusion of graph neural networks and 3D structure-based complex graph

Accurate prediction of protein-ligand binding affinity is pivotal for drug design and discovery. Here, we proposed a novel deep fusion graph neural networks framework named FGNN to learn the protein-ligand interactions from the 3D structures of protein-ligand complexes. Unlike 1D sequences for proteins or 2D graphs for ligands, the 3D graph of protein-ligand complex enables the more accurate representations of the protein-ligand interactions. Benchmark studies have shown that our fusion models FGNN can achieve more accurate prediction of binding affinity than any individual algorithm. The advantages of fusion strategies have been demonstrated in terms of expressive power of data, learning efficiency and model interpretability. Our fusion models show satisfactory performances on diverse data sets, demonstrating their generalization ability. Given the good performances in both binding affinity prediction and virtual screening, our fusion models are expected to be practically applied for drug screening and design. Our work highlights the potential of the fusion graph neural network algorithm in solving complex prediction problems in computational biology and chemistry. The fusion graph neural networks (FGNN) model is freely available in https://github.com/LinaDongXMU/FGNN.

A generalized protein-ligand scoring framework with balanced scoring, docking, ranking and screening powers

Applying machine learning algorithms to protein-ligand scoring functions has aroused widespread attention in recent years due to the high predictive accuracy and affordable computational cost. Nevertheless, most machine learning-based scoring functions are only applicable to a specific task, e.g., binding affinity prediction, binding pose prediction or virtual screening, suggesting that the development of a scoring function with balanced performance in all critical tasks remains a grand challenge. To this end, we propose a novel parameterization strategy by introducing an adjustable binding affinity term that represents the correlation between the predicted outcomes and experimental data into the training of mixture density network. The resulting residue-atom distance likelihood potential not only retains the superior docking and screening power over all the other state-of-the-art approaches, but also achieves a remarkable improvement in scoring and ranking performance. We emphatically explore the impacts of several key elements on prediction accuracy as well as the task preference, and demonstrate that the performance of scoring/ranking and docking/screening tasks of a certain model could be well balanced through an appropriate manner. Overall, our study highlights the potential utility of our innovative parameterization strategy as well as the resulting scoring framework in future structure-based drug design.

Discovery of a potent and selective allosteric inhibitor targeting the SHP2 tunnel site for RTK-driven cancer treatment

Src homology 2 domain-containing phosphatase 2 (SHP2) is a cytoplasmic protein tyrosine phosphatase (PTP) that regulates signal transduction of receptor tyrosine kinases (RTKs). Abnormal SHP2 activity is associated with tumorigenesis and metastasis. Because SHP2 contains multiple allosteric sites, identifying inhibitors at specific allosteric binding sites remains challenging. Here, we used structure-based virtual screening to directly search for the SHP2 "tunnel site" allosteric inhibitor. A novel hit (70) was identified as the SHP2 allosteric inhibitor with an IC₅₀ of 10.2 μM against full-length SHP2. Derivatization of hit compound 70 using molecular modeling-guided structure-based modification allowed the discovery of an effective and selective SHP2 inhibitor, compound 129, with 122-fold improved potency compared to the hit. Further studies revealed that 129 effectively inhibited signaling in multiple RTK-driven cancers and RTK inhibitor-resistant cancer cells. Remarkably, 129 was orally bioavailable (F = 55%) and significantly inhibited tumor growth in haematological malignancy. Taken together, compound 129 developed in this study may serve as a promising lead or candidate for cancers bearing RTK oncogenic drivers and SHP2-related diseases.

Discovery of the new alpha-glucosidase inhibitor with therapeutic potential in type 2 diabetes mellitus by a novel high-throughput virtual screening and free energy evaluation

Type 2 diabetes can cause a variety of complications, significantly affecting people's health. Given their ability to suppress carbohydrate digestion, alpha-glucosidase inhibitors are effective treatments for diabetes. However, the current approved glucosidase inhibitors' side effects of abdominal discomfort limit their use. We used the compound Pg3R from the natural fruit berry as a reference, screening against a large database of 22 million compounds to identify potential health-friendly alpha-glucosidase inhibitors. Ligand-based screening enables us to identify 3968 ligands that exhibit structural similarity compared to the natural compound. These lead hits were used for LeDock, and their binding free energies were evaluated by MM/GBSA. Among the top-scoring candidates, ZINC263584304 exhibited the strongest binding affinity to alpha-glucosidase, with a "low-fat" structural characteristic. Its recognition mechanism was further investigated by microsecond MD simulations and free energy landscapes, exhibiting novel conformational changes during the binding process. Our study provided a novel alpha-glucosidase inhibitor with the potential to treat type 2 diabetes.

Identification of a Family of Glycoside Derivatives Biologically Active against Acinetobacter baumannii and Other MDR Bacteria Using a QSPR Model

As the rate of discovery of new antibacterial compounds for multidrug-resistant bacteria is declining, there is an urge for the search for molecules that could revert this tendency. Acinetobacter baumannii has emerged as a highly virulent Gram-negative bacterium that has acquired multiple resistance mechanisms against antibiotics and is considered of critical priority. In this work, we developed a quantitative structure-property relationship (QSPR) model with 592 compounds for the identification of structural parameters related to their property as antibacterial agents against A. baumannii. QSPR mathematical validation (R2 = 70.27, RN = −0.008, a(R2) = 0.014, and δK = 0.021) and its prediction ability (Q2LMO = 67.89, Q2EXT = 67.75, a(Q2) = −0.068, δQ = 0.0, rm2¯ = 0.229, and Δrm2 = 0.522) were obtained with different statistical parameters; additional validation was done using three sets of external molecules (R2 = 72.89, 71.64 and 71.56). We used the QSPR model to perform a virtual screening on the BIOFACQUIM natural product database. From this screening, our model showed that molecules 32 to 35 and 54 to 68, isolated from different extracts of plants of the Ipomoea sp., are potential antibacterials against A. baumannii. Furthermore, biological assays showed that molecules 56 and 60 to 64 have a wide antibacterial activity against clinically isolated strains of A. baumannii, as well as other multidrug-resistant bacteria, including Staphylococcus aureus, Escherichia coli, Klebsiella pneumonia, and Pseudomonas aeruginosa. Finally, we propose 60 as a potential lead compound due to its broad-spectrum activity and its structural simplicity. Therefore, our QSPR model can be used as a tool for the investigation and search for new antibacterial compounds against A. baumannii.

Identification of LRRK2 Inhibitors through Computational Drug Repurposing

Parkinson's disease (PD) is the second most common neurodegenerative disorder that affects more than ten million people worldwide. However, the current PD treatments are still limited and alternative treatment strategies are urgently required. Leucine-rich repeat kinase 2 (LRRK2) has been recognized as a promising target for PD treatment. However, there are no approved LRRK2 inhibitors on the market. To rapidly identify potential drug repurposing candidates that inhibit LRRK2 kinase, we report a structure-based drug repurposing workflow that combines molecular docking, recursive partitioning model, molecular dynamics (MD) simulation, and molecular mechanics-generalized Born surface area (MM-GBSA) calculation. Thirteen compounds screened from our drug repurposing workflow were further evaluated through the experiment. The experimental results showed six drugs (Abivertinib, Aumolertinib, Encorafenib, Bosutinib, Rilzabrutinib, and Mobocertinib) with IC₅₀ less than 5 μM that were identified as potential LRRK2 kinase inhibitors. The most potent compound Abivertinib showed potent inhibitions with IC₅₀ toward G2019S mutation and wild-type LRRK2 of 410.3 nM and 177.0 nM, respectively. Our combination screening strategy had a 53% hit rate in this repurposing task. MD simulations and MM-GBSA free energy analysis further revealed the atomic binding mechanism between the identified drugs and G2019S LRRK2. In summary, the results showed that our drug repurposing workflow could be used to identify potent compounds for LRRK2. The potent inhibitors discovered in our work can be a starting point to develop more effective LRRK2 inhibitors.

Structural Analysis, Multi-Conformation Virtual Screening and Molecular Simulation to Identify Potential Inhibitors Targeting pS273R Proteases of African Swine Fever Virus

The African Swine Fever virus (ASFV) causes an infectious viral disease in pigs of all ages. The development of antiviral drugs primarily aimed at inhibition of proteases required for the proteolysis of viral polyproteins. In this study, the conformation of the pS273R protease in physiological states were investigated, virtually screened the multi-protein conformation of pS273R target proteins, combined various molecular docking scoring functions, and identified five potential drugs from the Food and Drug Administration drug library that may inhibit pS273R. Subsequent validation of the dynamic interactions of pS273R with the five putative inhibitors was achieved using molecular dynamics simulations and binding free energy calculations using the molecular mechanics/Poison-Boltzmann (Generalized Born) (MM/PB(GB)SA) surface area. These findings demonstrate that the arm domain and Thr159-Lys167 loop region of pS273R are significantly more flexible compared to the core structural domain, and the Thr159-Lys167 loop region can serve as a "gatekeeper" in the substrate channel. Leucovorin, Carboprost, Protirelin, Flavin Mononucleotide, and Lovastatin Acid all have Gibbs binding free energies with pS273R that were less than -20 Kcal/mol according to the MM/PBSA analyses. In contrast to pS273R in the free energy landscape, the inhibitor and drug complexes of pS273R showed distinct structural group distributions. These five drugs may be used as potential inhibitors of pS273R and may serve as future drug candidates for treating ASFV.

In silico prediction of boiling point, octanol-water partition coefficient, and retention time index of polycyclic aromatic hydrocarbons through machine learning

Polycyclic aromatic hydrocarbons (PAHs), a special class of persistent organic pollutants (POPs) with two or more aromatic rings, have received extensive attention owing to their carcinogenic, mutagenic, and teratogenic effects. Quantitative structure-property relationship (QSPR) is powerful chemometric method to correlate structural descriptors of PAHs with their physicochemical properties. In this manuscript, a QSPR study of PAHs was performed to predict their boiling point (bp), octanol-water partition coefficient (LogK_ow ), and retention time index (RI). In addition to traditional molecular descriptors, structural fingerprints play an important role in the correlation of the above properties. Three regression methods, partial least squares (PLS), multiple linear regression (MLR), and genetic function approximation (GFA), were used to establish QSPR models for each property of PAHs. The correlation coefficient (R² _test ) and root mean square error (RMSE) of best model were 0.980 and 24.39% (PLS), 0.979 and 35.80% (GFA), 0.926 and 22.90% (MLR) for bp, LogK_ow, and RI, respectively. The model proposed here can be used to estimate physicochemical properties and inform toxicity prediction of environmental chemicals.

Virtual Screening Identifies Novel and Potent Inhibitors of Mycobacterium tuberculosis PknB with Antibacterial Activity

Mycobacterium tuberculosis protein kinase B (PknB) is essential to mycobacterial growth and has received considerable attention as an attractive target for novel anti-tuberculosis drug development. Here, virtual screening, validated by biological assays, was applied to select candidate inhibitors of M. tuberculosis PknB from the Specs compound library (www.specs.net). Fifteen compounds were identified as hits and selected for in vitro biological assays, of which three indoles (2, AE-848/42799159; 4, AH-262/34335013; 10, AP-124/40904362) inhibited growth of M. tuberculosis H37Rv with minimal inhibitory concentrations of 6.2, 12.5, and 6.2 μg/mL, respectively. Two compounds, 2 and 10, inhibited M. tuberculosis PknB activity in vitro, with IC₅₀ values of 14.4 and 12.1 μM, respectively, suggesting this to be the likely basis of their anti-tubercular activity. In contrast, compound 4 displayed anti-tuberculosis activity against M. tuberculosis H37Rv but showed no inhibition of PknB activity (IC₅₀ > 128 μM). We hypothesize that hydrolysis of its ethyl ester to a carboxylate moiety generates an active species that inhibits other M. tuberculosis enzymes. Molecular dynamics simulations of modeled complexes of compounds 2, 4, and 10 bound to M. tuberculosis PknB indicated that compound 4 has a lower affinity for M. tuberculosis PknB than compounds 2 and 10, as evidenced by higher calculated binding free energies, consistent with experiment. Compounds 2 and 10 therefore represent candidate inhibitors of M. tuberculosis PknB that provide attractive starting templates for optimization as anti-tubercular agents.

Functional dynamics of SARS-CoV-2 3C-like protease as a member of clan PA

SARS-CoV-2 3C-like protease (3CL^pro), a potential therapeutic target for COVID-19, consists of a chymotrypsin fold and a C-terminal α-helical domain (domain III), the latter of which mediates dimerization required for catalytic activation. To gain further understanding of the functional dynamics of SARS-CoV-2 3CL^pro, this review extends the scope to the comparative study of many crystal structures of proteases having the chymotrypsin fold (clan PA of the MEROPS database). First, the close correspondence between the zymogen-enzyme transformation in chymotrypsin and the allosteric dimerization activation in SARS-CoV-2 3CL^pro is illustrated. Then, it is shown that the 3C-like proteases of family Coronaviridae (the protease family C30), which are closely related to SARS-CoV-2 3CL^pro, have the same homodimeric structure and common activation mechanism via domain III mediated dimerization. The survey extended to order Nidovirales reveals that all 3C-like proteases belonging to Nidovirales have domain III, but with various chain lengths, and 3CL^pro of family Mesoniviridae (family C107) has the same homodimeric structure as that of C30, even though they have no sequence similarity. As a reference, monomeric 3C proteases belonging to the more distant family Picornaviridae (family C3) lacking domain III are compared with C30, and it is shown that the 3C proteases are rigid enough to maintain their structures in the active state.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01020-x.

MILCDock: Machine Learning Enhanced Consensus Docking for Virtual Screening in Drug Discovery

Molecular docking tools are regularly used to computationally identify new molecules in virtual screening for drug discovery. However, docking tools suffer from inaccurate scoring functions with widely varying performance on different proteins. To enable more accurate ranking of active over inactive ligands in virtual screening, we created a machine learning consensus docking tool, MILCDock, that uses predictions from five traditional molecular docking tools to predict the probability a ligand binds to a protein. MILCDock was trained and tested on data from both the DUD-E and LIT-PCBA docking datasets and shows improved performance over traditional molecular docking tools and other consensus docking methods on the DUD-E dataset. LIT-PCBA targets proved to be difficult for all methods tested. We also find that DUD-E data, although biased, can be effective in training machine learning tools if care is taken to avoid DUD-E's biases during training.

Identification of Active Compounds against Melanoma Growth by Virtual Screening for Non-Classical Human DHFR Inhibitors

Antifolates such as methotrexate (MTX) have been largely known as anticancer agents because of their role in blocking nucleic acid synthesis and cell proliferation. Their mechanism of action lies in their ability to inhibit enzymes involved in the folic acid cycle, especially human dihydrofolate reductase (hDHFR). However, most of them have a classical structure that has proven ineffective against melanoma, and, therefore, inhibitors with a non-classical lipophilic structure are increasingly becoming an attractive alternative to circumvent this clinical resistance. In this study, we conducted a protocol combining virtual screening (VS) and cell-based assays to identify new potential non-classical hDHFR inhibitors. Among 173 hit compounds identified (average logP = 3.68; average MW = 378.34 Da), two-herein, called C1 and C2-exhibited activity against melanoma cell lines B16 and A375 by MTT and Trypan-Blue assays. C1 showed cell growth arrest (39% and 56%) and C2 showed potent cytotoxic activity (77% and 51%) in a dose-dependent manner. The effects of C2 on A375 cell viability were greater than MTX (98% vs 60%) at equivalent concentrations and times. Our results indicate that the integrated in silico/in vitro approach provided a benchmark to identify novel promising non-classical DHFR inhibitors showing activity against melanoma cells.

The use of machine learning modeling, virtual screening, molecular docking, and molecular dynamics simulations to identify potential VEGFR2 kinase inhibitors

Targeting the signaling pathway of the Vascular endothelial growth factor receptor-2 is a promising approach that has drawn attention in the quest to develop novel anti-cancer drugs and cardiovascular disease treatments. We construct a screening pipeline using machine learning classification integrated with similarity checks of approved drugs to find new inhibitors. The statistical metrics reveal that the random forest approach has slightly better performance. By further similarity screening against several approved drugs, two candidates are selected. Analysis of absorption, distribution, metabolism, excretion, and toxicity, along with molecular docking and dynamics are performed for the two candidates with regorafenib as a reference. The binding energies of molecule1, molecule2, and regorafenib are - 89.1, - 95.3, and - 87.4 (kJ/mol), respectively which suggest candidate compounds have strong binding to the target. Meanwhile, the median lethal dose and maximum tolerated dose for regorafenib, molecule1, and molecule2 are predicted to be 800, 1600, and 393 mg/kg, and 0.257, 0.527, and 0.428 log mg/kg/day, respectively. Also, the inhibitory activity of these compounds is predicted to be 7.23 and 7.31, which is comparable with the activity of pazopanib and sorafenib drugs. In light of these findings, the two compounds could be further investigated as potential candidates for anti-angiogenesis therapy.

PyPLIF HIPPOS and Receptor Ensemble Docking Increase the Prediction Accuracy of the Structure-Based Virtual Screening Protocol Targeting Acetylcholinesterase

In this article, the upgrading process of the structure-based virtual screening (SBVS) protocol targeting acetylcholinesterase (AChE) previously published in 2017 is presented. The upgraded version of PyPLIF called PyPLIF HIPPOS and the receptor ensemble docking (RED) method using AutoDock Vina were employed to calculate the ensemble protein–ligand interaction fingerprints (ensPLIF) in a retrospective SBVS campaign targeting AChE. A machine learning technique called recursive partitioning and regression trees (RPART) was then used to optimize the prediction accuracy of the protocol by using the ensPLIF values as the descriptors. The best protocol resulting from this research outperformed the previously published SBVS protocol targeting AChE.

Discovery of Epipodophyllotoxin-Derived B2 as Promising XooFtsZ Inhibitor for Controlling Bacterial Cell Division: Structure-Based Virtual Screening, Synthesis, and SAR Study

The emergence of phytopathogenic bacteria resistant to antibacterial agents has rendered previously manageable plant diseases intractable, highlighting the need for safe and environmentally responsible agrochemicals. Inhibition of bacterial cell division by targeting bacterial cell division protein FtsZ has been proposed as a promising strategy for developing novel antibacterial agents. We previously identified 4'-demethylepipodophyllotoxin (DMEP), a naturally occurring substance isolated from the barberry species Dysosma versipellis, as a novel chemical scaffold for the development of inhibitors of FtsZ from the rice blight pathogen Xanthomonas oryzae pv. oryzae (Xoo). Therefore, constructing structure-activity relationship (SAR) studies of DMEP is indispensable for new agrochemical discovery. In this study, we performed a structure-activity relationship (SAR) study of DMEP derivatives as potential XooFtsZ inhibitors through introducing the structure-based virtual screening (SBVS) approach and various biochemical methods. Notably, prepared compound B₂, a 4'-acyloxy DMEP analog, had a 50% inhibitory concentration of 159.4 µM for inhibition of recombinant XooFtsZ GTPase, which was lower than that of the parent DMEP (278.0 µM). Compound B₂ potently inhibited Xoo growth in vitro (minimum inhibitory concentration 153 mg L^-1) and had 54.9% and 48.4% curative and protective control efficiencies against rice blight in vivo. Moreover, compound B₂ also showed low toxicity for non-target organisms, including rice plant and mammalian cell. Given these interesting results, we provide a novel strategy to discover and optimize promising bactericidal compounds for the management of plant bacterial diseases.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Boosting Protein-Ligand Binding Pose Prediction and Virtual Screening Based on Residue-Atom Distance Likelihood Potential and Graph Transformer

The past few years have witnessed enormous progress toward applying machine learning approaches to the development of protein-ligand scoring functions. However, the robust performance and wide applicability of scoring functions remain a big challenge for increasing the success rate of docking-based virtual screening. Herein, a novel scoring function named RTMScore was developed by introducing a tailored residue-based graph representation strategy and several graph transformer layers for the learning of protein and ligand representations, followed by a mixture density network to obtain residue-atom distance likelihood potential. Our approach was resolutely validated on the CASF-2016 benchmark, and the results indicate that RTMScore can outperform almost all of the other state-of-the-art methods in terms of both the docking and screening powers. Further evaluation confirms the robustness of our approach that can not only retain its docking power on cross-docked poses but also achieve improved performance as a rescoring tool in larger-scale virtual screening.

Discovery of novel DprE1 inhibitors via computational bioactivity fingerprints and structure-based virtual screening

Decaprenylphosphoryl-β-D-ribose oxidase (DprE1) plays important roles in the biosynthesis of mycobacterium cell wall. DprE1 inhibitors have shown great potentials in the development of new regimens for tuberculosis (TB) treatment. In this study, an integrated molecular modeling strategy, which combined computational bioactivity fingerprints and structure-based virtual screening, was employed to identify potential DprE1 inhibitors. Two lead compounds (B2 and H3) that could inhibit DprE1 and thus kill Mycobacterium smegmatis in vitro were identified. Moreover, compound H3 showed potent inhibitory activity against Mycobacterium tuberculosis in vitro (MIC_Mtb = 1.25 μM) and low cytotoxicity against mouse embryo fibroblast NIH-3T3 cells. Our research provided an effective strategy to discover novel anti-TB lead compounds.

Synergy between machine learning and natural products cheminformatics: Application to the lead discovery of anthraquinone derivatives

Cheminformatics utilizing machine learning (ML) techniques have opened up a new horizon in drug discovery. This is owing to vast chemical space expansion with rocketing numbers of expected hits and lead compounds that match druggable macromolecular targets, in particular from natural compounds. Due to the natural products' (NP) structural complexity, uniqueness, and diversity, they could occupy a bigger space in pharmaceuticals, allowing the industry to pursue more selective leads in the nanomolar range of binding affinity. ML is an essential part of each step of the drug design pipeline, such as target prediction, compound library preparation, and lead optimization. Notably, molecular mechanic and dynamic simulations, induced docking, and free energy perturbations are essential in predicting best binding poses, binding free energy values, and molecular mechanics force fields. Those applications have leveraged from artificial intelligence (AI), which decreases the computational costs required for such costly simulations. This review aimed to describe chemical space and compound libraries related to NPs. High-throughput screening utilized for fractionating NPs and high-throughput virtual screening and their strategies, and significance, are reviewed. Particular emphasis was given to AI approaches, ML tools, algorithms, and techniques, especially in drug discovery of macrocyclic compounds and approaches in computer-aided and ML-based drug discovery. Anthraquinone derivatives were discussed as a source of new lead compounds that can be developed using ML tools for diverse medicinal uses such as cancer, infectious diseases, and metabolic disorders. Furthermore, the power of principal component analysis in understanding relevant protein conformations, and molecular modeling of protein-ligand interaction were also presented. Apart from being a concise reference for cheminformatics, this review is a useful text to understand the application of ML-based algorithms to molecular dynamics simulation and in silico absorption, distribution, metabolism, excretion, and toxicity prediction.

Drugsniffer: An Open Source Workflow for Virtually Screening Billions of Molecules for Binding Affinity to Protein Targets

The SARS-CoV2 pandemic has highlighted the importance of efficient and effective methods for identification of therapeutic drugs, and in particular has laid bare the need for methods that allow exploration of the full diversity of synthesizable small molecules. While classical high-throughput screening methods may consider up to millions of molecules, virtual screening methods hold the promise of enabling appraisal of billions of candidate molecules, thus expanding the search space while concurrently reducing costs and speeding discovery. Here, we describe a new screening pipeline, called drugsniffer, that is capable of rapidly exploring drug candidates from a library of billions of molecules, and is designed to support distributed computation on cluster and cloud resources. As an example of performance, our pipeline required ∼40,000 total compute hours to screen for potential drugs targeting three SARS-CoV2 proteins among a library of ∼3.7 billion candidate molecules.

Identification of allosteric fingerprints of alpha-synuclein aggregates in matrix metalloprotease-1 and substrate-specific virtual screening with single molecule insights

Alpha-synuclein (aSyn) has implications in pathological protein aggregations in neurodegeneration. Matrix metalloproteases (MMPs) are broad-spectrum proteases and cleave aSyn, leading to aggregation. Previous reports showed that allosteric communications between the two domains of MMP1 on collagen fibril and fibrin depend on substrates, activity, and ligands. This paper reports quantification of allostery using single molecule measurements of MMP1 dynamics on aSyn-induced aggregates by calculating Forster Resonance Energy Transfer (FRET) between two dyes attached to the catalytic and hemopexin domains of MMP1. The two domains of MMP1 prefer open conformations that are inhibited by a single point mutation E219Q of MMP1 and tetracycline, an MMP inhibitor. A two-state Poisson process describes the interdomain dynamics, where the two states and kinetic rates of interconversion between them are obtained from histograms and autocorrelations of FRET values. Since a crystal structure of aSyn-bound MMP1 is unavailable, binding poses were predicted by molecular docking of MMP1 with aSyn using ClusPro. MMP1 dynamics were simulated using predicted binding poses and compared with the experimental interdomain dynamics to identify an appropriate pose. The selected aSyn-MMP1 binding pose near aSyn residue K45 was simulated and analyzed to define conformational changes at the catalytic site. Allosteric residues in aSyn-bound MMP1 exhibiting strong correlations with the catalytic motif residues were compared with allosteric residues in free MMP1, and aSyn-specific residues were identified. The allosteric residues in aSyn-bound MMP1 are K281, T283, G292, G327, L328, E329, R337, F343, G345, N346, Y348, G353, Q354, D363, Y365, S366, S367, F368, P371, R372, V374, K375, A379, F391, A394, R399, M414, F419, V426, and C466. Shannon entropy was defined to quantify MMP1 dynamics. Virtual screening was performed against a site on selected aSyn-MMP1 binding poses, which showed that lead molecules differ between free MMP1 and substrate-bound MMP1. Also, identifying aSyn-specific allosteric residues in MMP1 enabled further selection of lead molecules. In other words, virtual screening needs to take substrates into account for potential substrate-specific control of MMP1 activity in the future. Molecular understanding of interactions between MMP1 and aSyn-induced aggregates may open up the possibility of degrading aggregates by targeting MMPs.

In Silico Drug Repositioning to Target the SARS-CoV-2 Main Protease as Covalent Inhibitors Employing a Combined Structure-Based Virtual Screening Strategy of Pharmacophore Models and Covalent Docking

The epidemic caused by the SARS-CoV-2 coronavirus, which has spread rapidly throughout the world, requires urgent and effective treatments considering that the appearance of viral variants limits the efficacy of vaccines. The main protease of SARS-CoV-2 (Mpro) is a highly conserved cysteine proteinase, fundamental for the replication of the coronavirus and with a specific cleavage mechanism that positions it as an attractive therapeutic target for the proposal of irreversible inhibitors. A structure-based strategy combining 3D pharmacophoric modeling, virtual screening, and covalent docking was employed to identify the interactions required for molecular recognition, as well as the spatial orientation of the electrophilic warhead, of various drugs, to achieve a covalent interaction with Cys145 of Mpro. The virtual screening on the structure-based pharmacophoric map of the SARS-CoV-2 Mpro in complex with an inhibitor N3 (reference compound) provided high efficiency by identifying 53 drugs (FDA and DrugBank databases) with probabilities of covalent binding, including N3 (Michael acceptor) and others with a variety of electrophilic warheads. Adding the energy contributions of affinity for non-covalent and covalent docking, 16 promising drugs were obtained. Our findings suggest that the FDA-approved drugs Vaborbactam, Cimetidine, Ixazomib, Scopolamine, and Bicalutamide, as well as the other investigational peptide-like drugs (DB04234, DB03456, DB07224, DB7252, and CMX-2043) are potential covalent inhibitors of SARS-CoV-2 Mpro.

Approaches for discovery of small-molecular antivirals targeting to influenza A virus PB2 subunit

Influenza is an acute respiratory infectious disease caused by influenza virus, leading to huge morbidity and mortality in humans worldwide. Despite the availability of antivirals in the clinic, the emergence of resistant strains calls for antivirals with novel mechanisms of action. The PB2 subunit of the influenza A virus polymerase is a promising target because of its vital role in the 'cap-snatching' mechanism. In this review, we summarize the technologies and medicinal chemistry strategies for hit identification, hit-to-lead and lead-to-candidate optimization, and current challenges in PB2 inhibitor development, as well as offering insights for the fight against drug resistance.

Investigation of Unprecedented Sites and Proposition of New Ligands for Programmed Cell Death Protein I through Molecular Dynamics with Probes and Virtual Screening

Cancer immunotherapy has attracted increasing attention over the last few years. Programmed cell death protein 1 (PD-1) promotes self-tolerance and inhibits immune responses by modulating the T-cell function. The interaction between PD-1 and programmed cell death ligand-1 (PD-L1) leads to immune exhaustion, protecting cancer cells from destruction. Here, we computationally designed a novel ligand named 1508 that binds to an unprecedented PD-1 cavity identified by MixMD and defined by amino acid residues Lys78 to Val97. We showed through a set of MD simulations totaling 12.5 μs that ligand 1508 establishes frequent cation-π and hydrogen bonding interactions with amino acid residues Lys78 and Arg86, respectively, and stabilizes the PD-1 C'D loop in a conformation that does not favor PD-1-PD-L1 complex formation. This study highlights the power of MixMD in exposing new cavities prone to protein-protein complex inhibition and establishes the basis for the design of new molecules that target the PD-1 C'D cavity as an alternative for exploring the modulation of the PD-1-PD-L1 complex in cancer therapy.

Computer-Aided Drug Design (CADD) to De-Orphanize Marine Molecules: Finding Potential Therapeutic Agents for Neurodegenerative and Cardiovascular Diseases

Computer-aided drug design (CADD) techniques allow the identification of compounds capable of modulating protein functions in pathogenesis-related pathways, which is a promising line on drug discovery. Marine natural products (MNPs) are considered a rich source of bioactive compounds, as the oceans are home to much of the planet's biodiversity. Biodiversity is directly related to chemodiversity, which can inspire new drug discoveries. Therefore, natural products (NPs) in general, and MNPs in particular, have been used for decades as a source of inspiration for the design of new drugs. However, NPs present both opportunities and challenges. These difficulties can be technical, such as the need to dive or trawl to collect the organisms possessing the compounds, or biological, due to their particular marine habitats and the fact that they can be uncultivable in the laboratory. For all these difficulties, the contributions of CADD can play a very relevant role in simplifying their study, since, for example, no biological sample is needed to carry out an in-silico analysis. Therefore, the amount of natural product that needs to be used in the entire preclinical and clinical study is significantly reduced. Here, we exemplify how this combination between CADD and MNPs can help unlock their therapeutic potential. In this study, using a set of marine invertebrate molecules, we elucidate their possible molecular targets and associated therapeutic potential, establishing a pipeline that can be replicated in future studies.

Screening S protein - ACE2 blockers from natural products: Strategies and advances in the discovery of potential inhibitors of COVID-19

The Coronavirus disease, 2019 (COVID-19) is caused by severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), which poses a major threat to human life and health. Given its continued development, limiting the spread of COVID-19 in the population remains a challenging task. Currently, multiple therapies are being tried around the world to deal with SARS-CoV-2 infection, and a variety of studies have shown that natural products have a significant effect on COVID-19 patients. The combination of SARS-CoV-2 S protein with Angiotensin converting enzyme II(ACE2) of host cell to promote membrane fusion is an initial critical step for SARS-CoV-2 infection. Therefore, screening natural products that inhibit the binding of SARS-CoV-2 S protein and ACE2 also provides a feasible strategy for the treatment of COVID-19. Establishment of high throughput screening model is an important basis and key technology for screening S protein-ACE2 blockers. Based on this, the molecular structures of SARS-CoV-2 and ACE2 and their processes in the life cycle of SARS-CoV-2 and host cell infection were firstly reviewed in this paper, with emphasis on the methods and techniques of screening S protein-ACE2 blockers, including Virtual Screening (VS), Surface Plasmon Resonance (SPR), Biochromatography, Biotin-avidin with Enzyme-linked Immunosorbent assay and Gene Chip Technology. Furthermore, the technical principle, advantages and disadvantages and application scope were further elaborated. Combined with the application of the above screening technologies in S protein-ACE2 blockers, a variety of natural products, such as flavonoids, terpenoids, phenols, alkaloids, were summarized, which could be used as S protein-ACE2 blockers, in order to provide ideas for the efficient discovery of S protein-ACE2 blockers from natural sources and contribute to the development of broad-spectrum anti coronavirus drugs.

Accelerating Molecular Dynamics Enrichments of High-Affinity Ligands for Proteins

Molecular docking algorithms are used to seek the most active compounds from a pool of ligands. In principle, molecular dynamics (MD) simulations with accurate physical potentials and sampling could yield better enrichments, but they are computationally expensive. Here, we describe a method called MELD-Bracket that utilizes biased replica exchange ladders in MD in order to compete different ligands against each other within a fast bracket style "binding tournament". MELD-Bracket finds best-binders rapidly when ligands are well separated in their binding affinities.

In Silico End-to-End Protein-Ligand Interaction Characterization Pipeline: The Case of SARS-CoV-2

SARS-CoV-2 triggered a worldwide pandemic disease, COVID-19, for which an effective treatment has not yet been settled. Among the most promising targets to fight this disease is SARS-CoV-2 main protease (Mpro), which has been extensively studied in the last few months. There is an urgency for developing effective computational protocols that can help us tackle these key viral proteins. Hence, we have put together a robust and thorough pipeline of in silico protein-ligand characterization methods to address one of the biggest biological problems currently plaguing our world. These methodologies were used to characterize the interaction of SARS-CoV-2 Mpro with an α-ketoamide inhibitor and include details on how to upload, visualize, and manage the three-dimensional structure of the complex and acquire high-quality figures for scientific publications using PyMOL (Protocol 1); perform homology modeling with MODELLER (Protocol 2); perform protein-ligand docking calculations using HADDOCK (Protocol 3); run a virtual screening protocol of a small compound database of SARS-CoV-2 candidate inhibitors with AutoDock 4 and AutoDock Vina (Protocol 4); and, finally, sample the conformational space at the atomic level between SARS-CoV-2 Mpro and the α-ketoamide inhibitor with Molecular Dynamics simulations using GROMACS (Protocol 5). Guidelines for careful data analysis and interpretation are also provided for each Protocol.

Rational Design of Anti-Angiogenic Peptides to Inhibit VEGF/VEGFR2 Interactions for Cancer Therapeutics

BACKGROUND

Angiogenesis is a critical physiological process that plays a key role in tumor progression, metastatic dissemination, and invasion. In the last two decades, the vascular endothelial growth factor (VEGF) signaling pathway has been the area of extensive researches. VEGF executes its special effects by binding to vascular endothelial growth factor receptors (VEGFRs), particularly VEGFR-2.

OBJECTIVE

The inhibition of VEGF/VEGFR2 interaction is known as an effective cancer therapy strategy. The current study pointed to design and model an anti-VEGF peptide based on VEGFR2 binding regions.

METHOD

The large-scale peptide mutation screening was used to achieve a potent peptide with high binding affinity to VEGF for possible application in inhibition of VEGF/VEGFR2 interaction. The AntiCP and Peptide Ranker servers were used to generate the possible peptides library with anticancer activities and prediction of peptides bioactivity. Then, the interaction of VEGF and all library peptides were analyzed using Hex 8.0.0 and ClusPro tools. A number of six peptides with favorable docking scores were achieved. All of the best docking scores of peptides in complexes with VEGF were evaluated to confirm their stability, using molecular dynamics simulation (MD) with the help of the GROMACS software package.

RESULTS

As a result, two antiangiogenic peptides with 13 residues of PepA (NGIDFNRDFFLGL) and PepC (NGIDFNRDKFLFL) were achieved and introduced to inhibit VEGF/VEGFR2 interactions.

CONCLUSIONS

In summary, this study provided new insights into peptide-based therapeutics development for targeting VEGF signaling pathway in tumor cells. PepA and PepC are recommended as potentially promising anticancer agents for further experimental evaluations.

Screening of β1- and β2-Adrenergic Receptor Modulators through Advanced Pharmacoinformatics and Machine Learning Approaches

Cardiovascular diseases (CDs) are a major concern in the human race and one of the leading causes of death worldwide. β-Adrenergic receptors (β1-AR and β2-AR) play a crucial role in the overall regulation of cardiac function. In the present study, structure-based virtual screening, machine learning (ML), and a ligand-based similarity search were conducted for the PubChem database against both β1- and β2-AR. Initially, all docked molecules were screened using the threshold binding energy value. Molecules with a better binding affinity were further used for segregation as active and inactive through ML. The pharmacokinetic assessment was carried out on molecules retained in the above step. Further, similarity searching of the ChEMBL and DrugBank databases was performed. From detailed analysis of the above data, four compounds for each of β1- and β2-AR were found to be promising in nature. A number of critical ligand-binding amino acids formed potential hydrogen bonds and hydrophobic interactions. Finally, a molecular dynamics (MD) simulation study of each molecule bound with the respective target was performed. A number of parameters obtained from the MD simulation trajectories were calculated and substantiated the stability between the protein-ligand complex. Hence, it can be postulated that the final molecules might be crucial for CDs subjected to experimental validation.

The impact of cross-docked poses on performance of machine learning classifier for protein-ligand binding pose prediction

Structure-based drug design depends on the detailed knowledge of the three-dimensional (3D) structures of protein-ligand binding complexes, but accurate prediction of ligand-binding poses is still a major challenge for molecular docking due to deficiency of scoring functions (SFs) and ignorance of protein flexibility upon ligand binding. In this study, based on a cross-docking dataset dedicatedly constructed from the PDBbind database, we developed several XGBoost-trained classifiers to discriminate the near-native binding poses from decoys, and systematically assessed their performance with/without the involvement of the cross-docked poses in the training/test sets. The calculation results illustrate that using Extended Connectivity Interaction Features (ECIF), Vina energy terms and docking pose ranks as the features can achieve the best performance, according to the validation through the random splitting or refined-core splitting and the testing on the re-docked or cross-docked poses. Besides, it is found that, despite the significant decrease of the performance for the threefold clustered cross-validation, the inclusion of the Vina energy terms can effectively ensure the lower limit of the performance of the models and thus improve their generalization capability. Furthermore, our calculation results also highlight the importance of the incorporation of the cross-docked poses into the training of the SFs with wide application domain and high robustness for binding pose prediction. The source code and the newly-developed cross-docking datasets can be freely available at https://github.com/sc8668/ml_pose_prediction and https://zenodo.org/record/5525936 , respectively, under an open-source license. We believe that our study may provide valuable guidance for the development and assessment of new machine learning-based SFs (MLSFs) for the predictions of protein-ligand binding poses.

Native or Non-Native Protein-Protein Docking Models? Molecular Dynamics to the Rescue

Molecular docking excels at creating a plethora of potential models of protein-protein complexes. To correctly distinguish the favorable, native-like models from the remaining ones remains, however, a challenge. We assessed here if a protocol based on molecular dynamics (MD) simulations would allow distinguishing native from non-native models to complement scoring functions used in docking. To this end, the first models for 25 protein-protein complexes were generated using HADDOCK. Next, MD simulations complemented with machine learning were used to discriminate between native and non-native complexes based on a combination of metrics reporting on the stability of the initial models. Native models showed higher stability in almost all measured properties, including the key ones used for scoring in the Critical Assessment of PRedicted Interaction (CAPRI) competition, namely the positional root mean square deviations and fraction of native contacts from the initial docked model. A random forest classifier was trained, reaching a 0.85 accuracy in correctly distinguishing native from non-native complexes. Reasonably modest simulation lengths of the order of 50-100 ns are sufficient to reach this accuracy, which makes this approach applicable in practice.

Discovery of potential inhibitors targeting the kinase domain of polynucleotide kinase/phosphatase (PNKP): Homology modeling, virtual screening based on multiple conformations, and molecular dynamics simulation

In recent years, the level of interest has been increased in developing the DNA-repair inhibitors, to enhance the cytotoxic effects in the treatment of cancers. Polynucleotide kinase/phosphatase (PNKP) is a critical human DNA repair enzyme that repairs DNA strand breaks by catalyzing the restoration of 5'-phosphate and 3'-hydroxyl termini that are required for subsequent processing by DNA ligases and polymerases. PNKP is the only protein that repairs the 3'-hydroxyl group and 5'-phosphate group, which depicts PNKP as a potential therapeutic target. Besides, PNKP is the only DNA-repair enzyme that contains the 5'-kinase activity, therefore, targeting this kinase domain would motivate the development of novel PNKP-specific inhibitors. However, there are neither crystal structures of human PNKP nor the kinase inhibitors reported so far. Thus, in this present study, a sequential molecular docking-based virtual screening with multiple PNKP conformations integrating homology modeling, molecular dynamics simulation, and binding free energy calculation was developed to discover novel PNKP kinase inhibitors, and the top-scored molecule was finally submitted to molecular dynamics simulation to reveal the binding mechanism between the inhibitor and PNKP. Taken together, the current study could provide some guidance for the molecular docking based-virtual screening of novel PNKP kinase inhibitors.

Methodological approaches for the analysis of transmembrane domain interactions: A systematic review

The study of protein-protein interactions (PPI) has proven fundamental for the understanding of the most relevant cell processes. Any protein domain can participate in PPI, including transmembrane (TM) segments that can establish interactions with other TM domains (TMDs). However, the hydrophobic nature of TMDs and the environment they occupy complicates the study of intramembrane PPI, which demands the use of specific approaches and techniques. In this review, we will explore some of the strategies available to study intramembrane PPI in vitro, in vivo, and, in silico, focusing on those techniques that could be carried out in a standard molecular biology laboratory regarding its previous experience with membrane proteins.

Current trends in computer aided drug design and a highlight of drugs discovered via computational techniques: A review

Computer-aided drug design (CADD) is one of the pivotal approaches to contemporary pre-clinical drug discovery, and various computational techniques and software programs are typically used in combination, in a bid to achieve the desired outcome. Several approved drugs have been developed with the aid of CADD. On SciFinder®, we evaluated more than 600 publications through systematic searching and refining, using the terms, virtual screening; software methods; computational studies and publication year, in order to obtain data concerning particular aspects of CADD. The primary focus of this review was on the databases screened, virtual screening and/or molecular docking software program used. Furthermore, we evaluated the studies that subsequently performed molecular dynamics (MD) simulations and we reviewed the software programs applied, the application of density functional theory (DFT) calculations and experimental assays. To represent the latest trends, the most recent data obtained was between 2015 and 2020, consequently the most frequently employed techniques and software programs were recorded. Among these, the ZINC database was the most widely preferred with an average use of 31.2%. Structure-based virtual screening (SBVS) was the most prominently used type of virtual screening and it accounted for an average of 57.6%, with AutoDock being the preferred virtual screening/molecular docking program with 41.8% usage. Following the screening process, 38.5% of the studies performed MD simulations to complement the virtual screening and GROMACS with 39.3% usage, was the popular MD software program. Among the computational techniques, DFT was the least applied whereby it only accounts for 0.02% average use. An average of 36.5% of the studies included reports on experimental evaluations following virtual screening. Ultimately, since the inception and application of CADD in pre-clinical drug discovery, more than 70 approved drugs have been discovered, and this number is steadily increasing over time.

A multi-conformational virtual screening approach based on machine learning targeting PI3Kγ

Nowadays, more and more attention has been attracted to develop selective PI3Kγ inhibitors, but the unique structural features of PI3Kγ protein make it a very big challenge. In the present study, a virtual screening strategy based on machine learning with multiple PI3Kγ protein structures was developed to screen novel PI3Kγ inhibitors. First, six mainstream docking programs were chosen to evaluate their scoring power and screening power; CDOCKER and Glide show satisfactory reliability and accuracy against the PI3Kγ system. Next, virtual screening integrating multiple PI3Kγ protein structures was demonstrated to significantly improve the screening enrichment rate comparing to that with an individual protein structure. Last, a multi-conformational Naïve Bayesian Classification model with the optimal docking programs was constructed, and it performed a true capability in the screening of PI3Kγ inhibitors. Taken together, the current study could provide some guidance for the docking-based virtual screening to discover novel PI3Kγ inhibitors.

Structural and Functional Diversity of Resistance-Nodulation-Cell Division Transporters

Multidrug resistant (MDR) bacteria are a global threat with many common infections becoming increasingly difficult to eliminate. While significant effort has gone into the development of potent biocides, the effectiveness of many first-line antibiotics has been diminished due to adaptive resistance mechanisms. Bacterial membrane proteins belonging to the resistance-nodulation-cell division (RND) superfamily play significant roles in mediating bacterial resistance to antimicrobials. They participate in multidrug efflux and cell wall biogenesis to transform bacterial pathogens into "superbugs" that are resistant even to last resort antibiotics. In this review, we summarize the RND superfamily of efflux transporters with a primary focus on the assembly and function of the inner membrane pumps. These pumps are critical for extrusion of antibiotics from the cell as well as the transport of lipid moieties to the outer membrane to establish membrane rigidity and stability. We analyze recently solved structures of bacterial inner membrane efflux pumps as to how they bind and transport their substrates. Our cumulative data indicate that these RND membrane proteins are able to utilize different oligomerization states to achieve particular activities, including forming MDR pumps and cell wall remodeling machineries, to ensure bacterial survival. This mechanistic insight, combined with simulated docking techniques, allows for the design and optimization of new efflux pump inhibitors to more effectively treat infections that today are difficult or impossible to cure.

Applications of Virtual Screening in Bioprospecting: Facts, Shifts, and Perspectives to Explore the Chemo-Structural Diversity of Natural Products

Natural products are continually explored in the development of new bioactive compounds with industrial applications, attracting the attention of scientific research efforts due to their pharmacophore-like structures, pharmacokinetic properties, and unique chemical space. The systematic search for natural sources to obtain valuable molecules to develop products with commercial value and industrial purposes remains the most challenging task in bioprospecting. Virtual screening strategies have innovated the discovery of novel bioactive molecules assessing in silico large compound libraries, favoring the analysis of their chemical space, pharmacodynamics, and their pharmacokinetic properties, thus leading to the reduction of financial efforts, infrastructure, and time involved in the process of discovering new chemical entities. Herein, we discuss the computational approaches and methods developed to explore the chemo-structural diversity of natural products, focusing on the main paradigms involved in the discovery and screening of bioactive compounds from natural sources, placing particular emphasis on artificial intelligence, cheminformatics methods, and big data analyses.

In silico Approaches for the Design and Optimization of Interfering Peptides Against Protein-Protein Interactions

Large contact surfaces of protein-protein interactions (PPIs) remain to be an ongoing issue in the discovery and design of small molecule modulators. Peptides are intrinsically capable of exploring larger surfaces, stable, and bioavailable, and therefore bear a high therapeutic value in the treatment of various diseases, including cancer, infectious diseases, and neurodegenerative diseases. Given these promising properties, a long way has been covered in the field of targeting PPIs via peptide design strategies. In silico tools have recently become an inevitable approach for the design and optimization of these interfering peptides. Various algorithms have been developed to scrutinize the PPI interfaces. Moreover, different databases and software tools have been created to predict the peptide structures and their interactions with target protein complexes. High-throughput screening of large peptide libraries against PPIs; "hotspot" identification; structure-based and off-structure approaches of peptide design; 3D peptide modeling; peptide optimization strategies like cyclization; and peptide binding energy evaluation are among the capabilities of in silico tools. In the present study, the most recent advances in the field of in silico approaches for the design of interfering peptides against PPIs will be reviewed. The future perspective of the field and its advantages and limitations will also be pinpointed.

Discovery of novel IDO1 inhibitors via structure-based virtual screening and biological assays

Indoleamine 2,3-dioxygenase 1 (IDO1) is a heme-containing enzyme that catalyzes the first and rate-limiting step in catabolism of tryptophan via the kynurenine pathway, which plays a pivotal role in the proliferation and differentiation of T cells. IDO1 has been proven to be an attractive target for many diseases, such as breast cancer, lung cancer, colon cancer, prostate cancer, etc. In this study, docking-based virtual screening and bioassays were conducted to identify novel inhibitors of IDO1. The cellular assay demonstrated that 24 compounds exhibited potent inhibitory activity against IDO1 at micromolar level, including 8 compounds with IC₅₀ values below 10 μM and the most potent one (compound 1) with IC₅₀ of 1.18 ± 0.04 μM. Further lead optimization based on similarity searching strategy led to the discovery of compound 28 as an excellent inhibitor with IC₅₀ of 0.27 ± 0.02 μM. Then, the structure-activity relationship of compounds 1, 2, 8 and 14 analogues is discussed. The interaction modes of two compounds against IDO1 were further explored through a Python Based Metal Center Parameter Builder (MCPB.py) molecular dynamics simulation, binding free energy calculation and electrostatic potential analysis. The novel IDO1 inhibitors of compound 1 and its analogues could be considered as promising scaffold for further development of IDO1 inhibitors.