Identification of potential Abl kinase inhibitors using virtual screening and free energy calculations for the treatment of chronic myeloid leukemia

Abl kinase, particularly the Bcr-Abl fusion protein, is a critical driver of chronic myeloid leukemia (CML) and remain significant therapeutic target in hematologic malignancies. Although first-generation tyrosine kinase inhibitors (TKIs) such as Imatinib have revolutionized CML treatment, resistance due to kinase domain mutations (e.g., T315I gatekeeper mutation) and side effects highlight needs for another candidate inhibitors. In this study, we employed a combined virtual screening strategy to discover novel Abl kinase inhibitors from an extensive chemical database (∼670 million compounds). Initially, shape-based similarity (USR/USRCAT) and electrostatic potential filters were used to refine the candidate compounds, followed by structure-based molecular docking against the Abl kinase domain. Top-ranked candidates were evaluated using molecular dynamics (MD) simulations and binding free energy calculations, such as MM/GBSA and free energy perturbation (FEP), to confirm stability and binding affinity. Five candidate compounds emerged with binding energies comparable to or higher than known Abl kinase inhibitors, including Imatinib and Bafetinib. Finally, based on these calculations, we selected two compounds as candidates as Abl tyrosine kinase inhibitors. Overall, the results showed the effectiveness of combining ligand-based and structure-based drug design strategies to identify new Abl kinase inhibitor leads for improved the CML therapy.

Integrated Approaches for Discovery of Staphylococcus aureus Antimicrobial Agents: Virtual Screening, Molecular Docking, Molecular Dynamics Simulations, and Density Functional Theory

Due to the excessive use of antibiotics, Staphylococcus aureus has developed resistance to conventional antibiotics. This study primarily employs virtual screening methods to explore the binding mode, biological stability, electronic properties, and antimicrobial activity of the drugs. Lactate dehydrogenase (LDH) was chosen as the primary target, and virtual screening of approximately 3180 FDA-approved drugs was performed. On the basis of binding affinity scores, the top 12 molecules were shortlisted for further analysis through precise docking and MMGBSA calculations. Molecular docking simulations revealed that these compounds exhibit a strong affinity for the target protein 6BAZ, with Gliquidone demonstrating the highest binding affinity at -76.25 kcal/mol. The top three hit molecules were subjected to 100 ns molecular dynamics simulations, which confirmed the stability of the ligand-protein complexes through hydrophobic and hydrogen bonding interactions, corroborating the docking and MMGBSA findings. Density functional theory (B3LYP level, 6-31 + G (d, p) basis set) was applied to evaluate molecular geometry optimization and vibrational frequencies, offering valuable insights into the structure and stability of the drug molecules, which further supports their potential as lead compounds for LDH inhibition and establishes a strong basis for future drug development and optimization.

Identification of a TonB-Dependent Siderophore Receptor as a Novel Anti-Biofilm Target and Virtual Screening for Its Inhibitor in Pseudomonas fluorescens PF08

Pseudomonas fluorescens is a vital food spoilage bacterium that commonly spoils foods in the biofilm state. Uncovering the targets responsible for biofilm formation and disrupting their function is a promising way to control bacterial biofilms and food spoilage. In this work, using the combination of qRT-PCR and construction of the gene deletion strain, Δtdsr, TonB-dependent siderophore receptor D7M10_RS23410 was, for the first time, proven to play an essential part in the biofilm development of P. fluorescens. By utilizing structure-based virtual screening technology, a natural compound, adenosine monophosphate (AMP), with the highest binding activity to D7M10_RS23410, was obtained as an effective biofilm inhibitor. AMP significantly decreased the cell autoaggregation and biofilm biomass at sub-MIC concentrations (2.5, 1.25, and 0.625 mg/mL), mainly through inhibiting the generation of extracellular polymeric substances (EPS) in the biofilm matrix and promoting the cell motility. Furthermore, AMP was found to form hydrogen bonds with specific amino acid residues and stretched the protein structure of D7M10_RS23410, and this structural alteration undoubtedly interfered with the functionality of the D7M10_RS23410 protein.

Virtual screening: hope, hype, and the fine line in between

INTRODUCTION

Technological advancements in virtual screening (VS) have rapidly accelerated its application in drug discovery, as reflected by the exponential growth in VS-related publications. However, a significant gap remains between the volume of computational predictions and their experimental validation. This discrepancy has led to a rise in the number of unverified 'claimed' hits which impedes the drug discovery efforts.

AREAS COVERED

This perspective examines the current VS landscape, highlighting essential practices and identifying critical challenges, limitations, and common pitfalls. Using case studies and practices, this perspective aims to highlight strategies that can effectively mitigate or overcome these challenges. Furthermore, the perspective explores common approaches for addressing pharmacodynamic and pharmacokinetic issues in optimizing VS hits.

EXPERT OPINION

VS has become a tried-and-true technique of drug discovery due to the rapid advances in computational methods and machine learning (ML) over the past two decades. Although each VS workflow varies depending on the chosen approach and methodology, integrated strategies that combine biological and in silico data have consistently yielded higher success rates. Moreover, the widespread adoption of ML has enhanced the integration of VS into the drug discovery pipeline. However, the absence of standardized evaluation criteria hinders the objective assessment of VS studies' success and the identification of optimal adoption methods.

In Silico Structural Insights and Potential Inhibitor Identification Based on the Benzothiazole Core for Targeting Leishmania major Pteridine Reductase 1

Leishmaniasis is reported as the second most common protozoonosis, with the highest prevalence and mortality rate. Among the Leishmania drug targets, Pteridine Reductase 1 of Leishmania major (LmPTR1) proved to be promising because Leishmania is auxotrophic for folates. Thus, this study employed a combination of ligand- and structure-based approaches to screen new benzothiazole compounds as LmPTR1 inhibitor candidates. Initially, a highly predictive quantitative structure-activity relationship (QSAR) model was able to identify the relevant hybrid descriptors, with an accuracy of over 93%. Insights into the mechanism of action indicated Phe113, His241, Leu188, Met183, and Leu226 as key residues. New commercially available compounds were screened using QSAR, docking, and pharmacokinetic properties as filters. Molecular dynamics, non-covalent interactions analysis, and quantum chemical calculation of binding enthalpy demonstrated that the lead compound (ZINC 72229720) forms a stable complex with LmPTR1, indicating it as a new promising LmPTR1 inhibitor.

A review on the screening methods for the discovery of natural antimicrobial peptides

Natural antimicrobial peptides (AMPs) are promising candidates for the development of a new generation of antimicrobials to combat antibiotic-resistant pathogens. They have found extensive applications in the fields of medicine, food, and agriculture. However, efficiently screening AMPs from natural sources poses several challenges, including low efficiency and high antibiotic resistance. This review focuses on the action mechanisms of AMPs, both through membrane and non-membrane routes. We thoroughly examine various highly efficient AMP screening methods, including whole-bacterial adsorption binding, cell membrane chromatography (CMC), phospholipid membrane chromatography binding, membrane-mediated capillary electrophoresis (CE), colorimetric assays, thin layer chromatography (TLC), fluorescence-based screening, genetic sequencing-based analysis, computational mining of AMP databases, and virtual screening methods. Additionally, we discuss potential developmental applications for enhancing the efficiency of AMP discovery. This review provides a comprehensive framework for identifying AMPs within complex natural product systems.

The potential of Mitragyna speciosa leaves as a natural source of antioxidants for disease prevention

Mitragyna speciosa is famous for its addictive effect. On the other hand, this plant has good potential as an antioxidant agent, and so far, it was not explicitly explained what the most contributing compound in the leaves to that activity is. This study has been conducted using several computational methods to determine which compounds are the most active in interacting with cytochrome P450, myeloperoxidase, and NADPH oxidase proteins. First, virtual screening was carried out based on molecular docking, followed by profiling the properties of adsorption, distribution, metabolism, excretion, and toxicity (ADMET); the second one is the molecular dynamics (MD) simulations for 100 ns. The virtual screening results showed that three compounds acted as inhibitors for each protein: (-)-epicatechin, sitogluside, and corynoxeine. The ADMET profiles of the three compounds exhibit good drug ability and toxicity. The trajectories study from MD simulations predicts that the complexes of these three compounds with their respective target proteins are stable. Furthermore, these compounds identified in this computational study can be a potential guide for future experiments aimed at assessing the antioxidant properties through in vitro testing.

Host-Guest Binding Free Energies à la Carte: An Automated OneOPES Protocol

Estimating absolute binding free energies from molecular simulations is a key step in computer-aided drug design pipelines, but the agreement between computational results and experiments is still very inconsistent. Both the accuracy of the computational model and the quality of the statistical sampling contribute to this discrepancy, yet disentangling the two remains a challenge. In this study, we present an automated protocol based on OneOPES, an enhanced sampling method that exploits replica exchange and can accelerate several collective variables to address the sampling problem. We apply this protocol to 37 host-guest systems. The simplicity of setting up the simulations and producing well-converged binding free energy estimates without the need to optimize simulation parameters provides a reliable solution to the sampling problem. This, in turn, allows for a systematic force field comparison and ranking according to the correlation between simulations and experiments, which can inform the selection of an appropriate model. The protocol can be readily adapted to test more force field combinations and study more complex protein-ligand systems, where the choice of an appropriate physical model is often based on heuristic considerations rather than systematic optimization.

DockThor-VS: A Free Platform for Receptor-Ligand Virtual Screening

The DockThor-VS platform (https://dockthor.lncc.br/v2/) is a free protein-ligand docking server conceptualized to facilitate and assist drug discovery projects to perform docking-based virtual screening experiments accurately and using high-performance computing. The DockThor docking engine is a grid-based method designed for flexible-ligand and rigid-receptor docking. It employs a multiple-solution genetic algorithm and the MMFF94S molecular force field scoring function for pose prediction. This engine was engineered to handle highly flexible ligands, such as peptides. Affinity prediction and ranking of protein-ligand complexes are performed with the linear empirical scoring function DockTScore. The main steps of the ligand and protein preparation are available on the DockThor Portal, making it possible to change the protonation states of the amino acid residues, and include cofactors as rigid entities. The user can also customize and visualize the main parameters of the grid box. The results of docking experiments are automatically clustered and ordered, providing users with a diverse array of meaningful binding modes. The platform DockThor-VS offers a user-friendly interface and powerful algorithms, enabling researchers to conduct virtual screening experiments efficiently and accurately. The DockThor Portal utilizes the computational strength of the Brazilian high-performance platform SDumont, further amplifying the efficiency and speed of docking experiments. Additionally, the web server facilitates and enhances virtual screening experiments by offering curated structures of potential targets and compound datasets, such as proteins related to COVID-19 and FDA-approved drugs for repurposing studies. In summary, DockThor-VS is a dynamic and evolving solution for docking-based virtual screening to be applied in drug discovery projects.

Increased throughput in methods for simulating protein ligand binding and unbinding

By incorporating full flexibility and enabling the quantification of crucial parameters such as binding free energies and residence times, methods for investigating protein-ligand binding and unbinding via molecular dynamics provide details on the involved mechanisms at the molecular level. While these advancements hold promise for impacting drug discovery, a notable drawback persists: their relatively time-consuming nature limits throughput. Herein, we survey recent implementations which, employing a blend of enhanced sampling techniques, a clever choice of collective variables, and often machine learning, strive to enhance the efficiency of new and previously reported methods without compromising accuracy. Particularly noteworthy is the validation of these methods that was often performed on systems mirroring real-world drug discovery scenarios.

Consensus holistic virtual screening for drug discovery: a novel machine learning model approach

In drug discovery, virtual screening is crucial for identifying potential hit compounds. This study aims to present a novel pipeline that employs machine learning models that amalgamates various conventional screening methods. A diverse array of protein targets was selected, and their corresponding datasets were subjected to active/decoy distribution analysis prior to scoring using four distinct methods: QSAR, Pharmacophore, docking, and 2D shape similarity, which were ultimately integrated into a single consensus score. The fine-tuned machine learning models were ranked using the novel formula "w_new", consensus scores were calculated, and an enrichment study was performed for each target. Distinctively, consensus scoring outperformed other methods in specific protein targets such as PPARG and DPP4, achieving AUC values of 0.90 and 0.84, respectively. Remarkably, this approach consistently prioritized compounds with higher experimental PIC50 values compared to all other screening methodologies. Moreover, the models demonstrated a range of moderate to high performance in terms of R2 values during external validation. In conclusion, this novel workflow consistently delivered superior results, emphasizing the significance of a holistic approach in drug discovery, where both quantitative metrics and active enrichment play pivotal roles in identifying the best virtual screening methodology.Scientific contributionWe presented a novel consensus scoring workflow in virtual screening, merging diverse methods for enhanced compound selection. We also introduced 'w_new', a groundbreaking metric that intricately refines machine learning model rankings by weighing various model-specific parameters, revolutionizing their efficacy in drug discovery in addition to other domains.

A combination of virtual screening, molecular dynamics simulation, MM/PBSA, ADMET, and DFT calculations to identify a potential DPP4 inhibitor

DPP4 inhibitors can control glucose homeostasis by increasing the level of GLP-1 incretins hormone due to dipeptidase mimicking. Despite the potent effects of DPP4 inhibitors, these compounds cause unwanted toxicity attributable to their effect on other enzymes. As a result, it seems essential to find novel and DPP4 selective compounds. In this study, we introduce a potent and selective DPP4 inhibitor via structure-based virtual screening, molecular docking, molecular dynamics simulation, MM/PBSA calculations, DFT analysis, and ADMET profile. The screened compounds based on similarity with FDA-approved DPP4 inhibitors were docked towards the DPP4 enzyme. The compound with the highest docking score, ZINC000003015356, was selected. For further considerations, molecular docking studies were performed on selected ligands and FDA-approved drugs for DPP8 and DPP9 enzymes. Molecular dynamics simulation was run during 200 ns and the analysis of RMSD, RMSF, Rg, PCA, and hydrogen bonding were performed. The MD outputs showed stability of the ligand-protein complex compared to available drugs in the market. The total free binding energy obtained for the proposed DPP4 inhibitor was more negative than its co-crystal ligand (N7F). ZINC000003015356 confirmed the role of the five Lipinski rule and also, have low toxicity parameter according to properties. Finally, DFT calculations indicated that this compound is sufficiently soft.

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.

The role and potential of computer-aided drug discovery strategies in the discovery of novel antimicrobials

Antimicrobial resistance (AMR) has become more of a concern in recent decades, particularly in infections associated with global public health threats. The development of new antibiotics is crucial to ensuring infection control and eradicating AMR. Although drug discovery and development are essential processes in the transformation of a drug candidate from the laboratory to the bedside, they are often very complicated, expensive, and time-consuming. The pharmaceutical sector is continuously innovating strategies to reduce research costs and accelerate the development of new drug candidates. Computer-aided drug discovery (CADD) has emerged as a powerful and promising technology that renews the hope of researchers for the faster identification, design, and development of cheaper, less resource-intensive, and more efficient drug candidates. In this review, we discuss an overview of AMR, the potential, and limitations of CADD in AMR drug discovery, and case studies of the successful application of this technique in the rapid identification of various drug candidates. This review will aid in achieving a better understanding of available CADD techniques in the discovery of novel drug candidates against resistant pathogens and other infectious agents.

Virtual Screening Process: A Guide in Modern Drug Designing

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

Scipion-Chem: An Open Platform for Virtual Drug Screening

Virtual drug screening (VDS) tackles the problem of drug discovery by computationally reducing the number of potential pharmacological molecules that need to be tested experimentally to find a new drug. To do so, several approaches have been developed through the years, typically focusing on either the physicochemical characteristics of the receptor structure (structure-based virtual screening) or those of the potential ligands (ligand-based virtual screening). Scipion is a workflow engine well suited for structural studies of biological macromolecules. Here, we present Scipion-chem, a new branch oriented to VDS. A total of 11 plugins have already been integrated from the most common programs used in the field. They can be used through the Scipion graphical user interface to execute and analyze typical VDS tasks. In addition, we have developed several consensus protocols that combine results from the different integrated programs to generate more robust predictions. Backstage, Scipion also facilitates the interoperability of those different software packages while tracking all of the intermediate files, parameters, and user decisions. In summary, in this article, we present Scipion-chem. This accessible, interoperable, and traceable platform provides the user with all of the tools to carry out a successful VDS workflow. Scipion-chem is openly available at https://github.com/scipion-chem.

Deep Learning with Geometry-Enhanced Molecular Representation for Augmentation of Large-Scale Docking-Based Virtual Screening

Structure-based virtual screening has been a crucial tool in drug discovery for decades. However, as the chemical space expands, the existing structure-based virtual screening techniques based on molecular docking and scoring struggle to handle billion-entry ultralarge libraries due to the high computational cost. To address this challenge, people have resorted to machine learning techniques to enhance structure-based virtual screening for efficiently exploring the vast chemical space. In those cases, compounds are usually treated as sequential strings or two-dimensional topology graphs, limiting their ability to incorporate three-dimensional structural information for downstream tasks. We herein propose a novel deep learning protocol, GEM-Screen, which utilizes the geometry-enhanced molecular representation of the compounds docking to a specific target and is trained on docking scores of a small fraction of a library through an active learning strategy to approximate the docking outcome for yet nontraining entries. This protocol is applied to virtual screening campaigns against the AmpC and D4 targets, demonstrating that GEM-Screen enriches more than 90% of the hit scaffolds for AmpC in the top 4% of model predictions and more than 80% of the hit scaffolds for D4 in the same top-ranking size of library. GEM-Screen can be used in conjunction with traditional docking programs for docking of only the top-ranked compounds to avoid the exhaustive docking of the whole library, thus allowing for discovering top-scoring compounds from billion-entry libraries in a rapid yet accurate fashion.

Improving drug discovery with a hybrid deep generative model using reinforcement learning trained on a Bayesian docking approximation

Generative approaches to molecular design are an area of intense study in recent years as a method to generate new pharmaceuticals with desired properties. Often though, these types of efforts are constrained by limited experimental activity data, resulting in either models that generate molecules with poor performance or models that are overfit and produce close analogs of known molecules. In this paper, we reduce this data dependency for the generation of new chemotypes by incorporating docking scores of known and de novo molecules to expand the applicability domain of the reward function and diversify the compounds generated during reinforcement learning. Our approach employs a deep generative model initially trained using a combination of limited known drug activity and an approximate docking score provided by a second machine learned Bayes regression model, with final evaluation of high scoring compounds by a full docking simulation. This strategy results in molecules with docking scores improved by 10-20% compared to molecules of similar size, while being 130 × faster than a docking only approach on a typical GPU workstation. We also show that the increased docking scores correlate with (1) docking poses with interactions similar to known inhibitors and (2) result in higher MM-GBSA binding energies comparable to the energies of known DDR1 inhibitors, demonstrating that the Bayesian model contains sufficient information for the network to learn to efficiently interact with the binding pocket during reinforcement learning. This outcome shows that the combination of the learned latent molecular representation along with the feature-based docking regression is sufficient for reinforcement learning to infer the relationship between the molecules and the receptor binding site, which suggest that our method can be a powerful tool for the discovery of new chemotypes with potential therapeutic applications.

Enhancement by pyrazolones of colistin efficacy against mcr-1-expressing E. coli: an in silico and in vitro investigation

Owing to the emergence of antibiotic resistance, the polymyxin colistin has been recently revived to treat acute, multidrug-resistant Gram-negative bacterial infections. Positively charged colistin binds to negatively charged lipids and damages the outer membrane of Gram-negative bacteria. However, the MCR-1 protein, encoded by the mobile colistin resistance (mcr) gene, is involved in bacterial colistin resistance by catalysing phosphoethanolamine (PEA) transfer onto lipid A, neutralising its negative charge, and thereby reducing its interaction with colistin. Our preliminary results showed that treatment with a reference pyrazolone compound significantly reduced colistin minimal inhibitory concentrations in Escherichia coli expressing mcr-1 mediated colistin resistance (Hanpaibool et al. in ACS Omega, 2023). A docking-MD combination was used in an ensemble-based docking approach to identify further pyrazolone compounds as candidate MCR-1 inhibitors. Docking simulations revealed that 13/28 of the pyrazolone compounds tested are predicted to have lower binding free energies than the reference compound. Four of these were chosen for in vitro testing, with the results demonstrating that all the compounds tested could lower colistin MICs in an E. coli strain carrying the mcr-1 gene. Docking of pyrazolones into the MCR-1 active site reveals residues that are implicated in ligand-protein interactions, particularly E246, T285, H395, H466, and H478, which are located in the MCR-1 active site and which participate in interactions with MCR-1 in ≥ 8/10 of the lowest energy complexes. This study establishes pyrazolone-induced colistin susceptibility in E. coli carrying the mcr-1 gene, providing a method for the development of novel treatments against colistin-resistant bacteria.

Advances in the discovery of new chemotypes through ultra-large library docking

INTRODUCTION

The size and complexity of virtual screening libraries in drug discovery have skyrocketed in recent years, reaching up to multiple billions of accessible compounds. However, virtual screening of such ultra-large libraries poses several challenges associated with preparing the libraries, sampling, and pre-selection of suitable compounds. The utilization of artificial intelligence (AI)-assisted screening approaches, such as deep learning, poses a promising countermeasure to deal with this rapidly expanding chemical space. For example, various AI-driven methods were recently successfully used to identify novel small molecule inhibitors of the SARS-CoV-2 main protease (M^pro).

AREAS COVERED

This review focuses on presenting various kinds of virtual screening methods suitable for dealing with ultra-large libraries. Challenges associated with these computational methodologies are discussed, and recent advances are highlighted in the example of the discovery of novel M^pro inhibitors targeting the SARS-CoV-2 virus.

EXPERT OPINION

With the rapid expansion of the virtual chemical space, the methodologies for docking and screening such quantities of molecules need to keep pace. Employment of AI-driven screening compounds has already been shown to be effective in a range from a few thousand to multiple billion compounds, furthered by de novo generation of drug-like molecules without human interference.

Identification of Potential Cytochrome P450 3A5 Inhibitors: An Extensive Virtual Screening through Molecular Docking, Negative Image-Based Screening, Machine Learning and Molecular Dynamics Simulation Studies

Cytochrome P450 3A5 (CYP3A5) is one of the crucial CYP family members and has already proven to be an important drug target for cardiovascular diseases. In the current study, the PubChem database was screened through molecular docking and high-affinity molecules were adopted for further assessment. A negative image-based (NIB) model was used for a similarity search by considering the complementary shape and electrostatics of the target and small molecules. Further, the molecules were segregated into active and inactive groups through six machine learning (ML) matrices. The active molecules found in each ML model were used for in silico pharmacokinetics and toxicity assessments. A total of five molecules followed the acceptable pharmacokinetics and toxicity profiles. Several potential binding interactions between the proposed molecules and CYP3A5 were observed. The dynamic behavior of the selected molecules in the CYP3A5 was explored through a molecular dynamics (MD) simulation study. Several parameters obtained from the MD simulation trajectory explained the stability of the protein–ligand complexes in dynamic states. The high binding affinity of each molecule was revealed by the binding free energy calculation through the MM-GBSA methods. Therefore, it can be concluded that the proposed molecules might be potential CYP3A5 molecules for therapeutic application in cardiovascular diseases subjected to in vitro/in vivo validations.

DOCKSTRING: Easy Molecular Docking Yields Better Benchmarks for Ligand Design

The field of machine learning for drug discovery is witnessing an explosion of novel methods. These methods are often benchmarked on simple physicochemical properties such as solubility or general druglikeness, which can be readily computed. However, these properties are poor representatives of objective functions in drug design, mainly because they do not depend on the candidate compound's interaction with the target. By contrast, molecular docking is a widely applied method in drug discovery to estimate binding affinities. However, docking studies require a significant amount of domain knowledge to set up correctly, which hampers adoption. Here, we present dockstring, a bundle for meaningful and robust comparison of ML models using docking scores. dockstring consists of three components: (1) an open-source Python package for straightforward computation of docking scores, (2) an extensive dataset of docking scores and poses of more than 260,000 molecules for 58 medically relevant targets, and (3) a set of pharmaceutically relevant benchmark tasks such as virtual screening or de novo design of selective kinase inhibitors. The Python package implements a robust ligand and target preparation protocol that allows nonexperts to obtain meaningful docking scores. Our dataset is the first to include docking poses, as well as the first of its size that is a full matrix, thus facilitating experiments in multiobjective optimization and transfer learning. Overall, our results indicate that docking scores are a more realistic evaluation objective than simple physicochemical properties, yielding benchmark tasks that are more challenging and more closely related to real problems in drug discovery.

Development of an Automatic Pipeline for Participation in the CELPP Challenge

The prediction of how a ligand binds to its target is an essential step for Structure-Based Drug Design (SBDD) methods. Molecular docking is a standard tool to predict the binding mode of a ligand to its macromolecular receptor and to quantify their mutual complementarity, with multiple applications in drug design. However, docking programs do not always find correct solutions, either because they are not sampled or due to inaccuracies in the scoring functions. Quantifying the docking performance in real scenarios is essential to understanding their limitations, managing expectations and guiding future developments. Here, we present a fully automated pipeline for pose prediction validated by participating in the Continuous Evaluation of Ligand Pose Prediction (CELPP) Challenge. Acknowledging the intrinsic limitations of the docking method, we devised a strategy to automatically mine and exploit pre-existing data, defining—whenever possible—empirical restraints to guide the docking process. We prove that the pipeline is able to generate predictions for most of the proposed targets as well as obtain poses with low RMSD values when compared to the crystal structure. All things considered, our pipeline highlights some major challenges in the automatic prediction of protein–ligand complexes, which will be addressed in future versions of the pipeline.

Characterizing the stabilization effects of stabilizers in protein-protein systems with end-point binding free energy calculations

Drug design targeting protein-protein interactions (PPIs) associated with the development of diseases has been one of the most important therapeutic strategies. Besides interrupting the PPIs with PPI inhibitors/blockers, increasing evidence shows that stabilizing the interaction between two interacting proteins may also benefit the therapy, such as the development of various types of molecular glues/stabilizers that mostly work by stabilizing the two interacting proteins to regulate the downstream biological effects. However, characterizing the stabilization effect of a stabilizer is usually hard or too complicated for traditional experiments since it involves ternary interactions [protein-protein-stabilizer (PPS) interaction]. Thus, developing reliable computational strategies will facilitate the discovery/design of molecular glues or PPI stabilizers. Here, by fully analyzing the energetic features of the binary interactions in the PPS ternary complex, we systematically investigated the performance of molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) methods on characterizing the stabilization effects of stabilizers in 14-3-3 systems. The results show that both MM/PBSA and MM/GBSA are powerful tools in distinguishing the stabilizers from the decoys (with area under the curves of 0.90-0.93 for all tested cases) and are reasonable for ranking protein-peptide interactions in the presence or absence of stabilizers as well (with the average Pearson correlation coefficient of ~0.6 at a relatively high dielectric constant for both methods). Moreover, to give a detailed picture of the stabilization effects, the stabilization mechanism is also analyzed from the structural and energetic points of view for individual systems containing strong or weak stabilizers. This study demonstrates a potential strategy to accelerate the discovery of PPI stabilizers.

Discovery of Novel BRD4 Ligand Scaffolds by Automated Navigation of the Fragment Chemical Space

Fragment-based drug discovery (FBDD) is a very effective hit identification method. However, the evolution of fragment hits into suitable leads remains challenging and largely artisanal. Fragment evolution is often scaffold-centric, meaning that its outcome depends crucially on the chemical structure of the starting fragment. Considering that fragment screening libraries cover only a small proportion of the corresponding chemical space, hits should be seen as probes highlighting privileged areas of the chemical space rather than actual starting points. We have developed an automated computational pipeline to mine the chemical space around any specific fragment hit, rapidly finding analogues that share a common interaction motif but are structurally novel and diverse. On a prospective application on the bromodomain-containing protein 4 (BRD4), starting from a known fragment, the platform yields active molecules with nonobvious scaffold changes. The procedure is fast and inexpensive and has the potential to uncover many hidden opportunities in FBDD.