DOCK 4.0: search strategies for automated molecular docking of flexible molecule databases

Identification of novel antistaphylococcal hit compounds

Staphylococcus aureus is one of the most common nosocomial biofilm-forming pathogens worldwide that has developed resistance mechanisms against majority of the antibiotics. Therefore, the search of novel antistaphylococcal agents with unexploited mechanisms of action, especially with antibiofilm activity, is of great interest. Seryl-tRNA synthetase is recognized as a promising drug target for the development of antibacterials. We have carried out molecular docking of compounds with antistaphycoccal activity, which were earlier found by us using phenotypic screening, into synthetic site of S. aureus SerRS and found seven hit compounds with low inhibitory activity. Further, we have performed search of S. aureus SerRS inhibitors among compounds which were previously tested by us for inhibitory activity toward S. aureus ThrRS, that belong to the same class of aminoacyl-tRNA synthetases. Among them six hits were identified. We have selected four compounds for antibacterial study and found that the most active compound 1-methyl-3-(1H-imidazol-1-methyl-2-yl)-5-nitro-1H-indazole has MIC values toward S. aureus multidrug-resistant clinical isolates ranging from 78.12 to 156.2 µg/ml. However, this compound precipitated during anti-biofilm study. Therefore, we used 3-[N'-(2-hydroxy-3-methoxybenzylidene)hydrazino]-6-methyl-4H-[1,2,4]triazin-5-one with better solubility (ClogS value = 2.9) among investigated compounds toward SerRS for anti-biofilm study. It was found that this compound has a significant inhibitory effect on the growth of planktonic and biofilm culture of S. aureus 25923 with MIC value of 32 µg ml-1. At the same time, this compound does not reveal antibacterial activity toward Esherichia coli ATCC 47076. Therefore, this compound can be proposed as effective antiseptic toward multidrug-resistant biofilm-forming S. aureus isolates.

Implicit Solvent with Explicit Ions Generalized Born Model in Molecular Dynamics: Application to DNA

The ion atmosphere surrounding highly charged biomolecules, such as nucleic acids, is crucial for their dynamics, structure, and interactions. Here, we develop an approach for the explicit treatment of ions within an implicit solvent framework suitable for atomistic simulations of biomolecules. The proposed implicit solvent/explicit ions model, GBION, is based on a modified generalized Born (GB) model; it includes separate, modified GB terms for solute-ion and ion-ion interactions. The model is implemented in the AMBER package (version 24), and its performance is thoroughly investigated in atomistic molecular dynamics (MD) simulations of double-stranded DNA on a microsecond time scale. The aggregate characteristics of monovalent (Na+ and K+) and trivalent (Cobalt Hexammine, CoHex3+) counterion distributions around double-stranded DNA predicted by the model are in reasonable agreement with the experiment (where available), all-atom explicit water MD simulations, and the expectation from the Manning condensation theory. The radial distributions of monovalent cations around DNA are reasonably close to the ones obtained using the explicit water model: expressed in units of energy, the maximum deviations of local ion concentrations from the explicit solvent reference are within 1 kBT, comparable to the corresponding deviations expected between different established explicit water models. The proposed GBION model is able to simulate DNA fragments in a large volume of solvent with explicit ions with little additional computational overhead compared with the fully implicit GB treatment of ions. Ions simulated using the developed model explore conformational space at least 2 orders of magnitude faster than in the explicit solvent. These advantages allowed us to observe and explore an unexpected "stacking" mode of DNA condensation in the presence of trivalent counterions (CoHex3+) that was revealed by recent experiments.

Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms

In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.

Combining the Fragment Molecular Orbital and GRID Approaches for the Prediction of Ligand-Metalloenzyme Binding Affinity: The Case Study of hCA II Inhibitors

Polarization and charge-transfer interactions play an important role in ligand-receptor complexes containing metals, and only quantum mechanics methods can adequately describe their contribution to the binding energy. In this work, we selected a set of benzenesulfonamide ligands of human Carbonic Anhydrase II (hCA II)-an important druggable target containing a Zn2+ ion in the active site-as a case study to predict the binding free energy in metalloprotein-ligand complexes and designed specialized computational methods that combine the ab initio fragment molecular orbital (FMO) method and GRID approach. To reproduce the experimental binding free energy in these systems, we adopted a machine-learning approach, here named formula generator (FG), considering different FMO energy terms, the hydrophobic interaction energy (computed by GRID) and logP. The main advantage of the FG approach is that it can find nonlinear relations between the energy terms used to predict the binding free energy, explicitly showing their mathematical relation. This work showed the effectiveness of the FG approach, and therefore, it might represent an important tool for the development of new scoring functions. Indeed, our scoring function showed a high correlation with the experimental binding free energy (R2 = 0.76-0.95, RMSE = 0.34-0.18), revealing a nonlinear relation between energy terms and highlighting the relevant role played by hydrophobic contacts. These results, along with the FMO characterization of ligand-receptor interactions, represent important information to support the design of new and potent hCA II inhibitors.

Advancing Drug-Target Interaction prediction with BERT and subsequence embedding

Exploring the relationship between proteins and drugs plays a significant role in discovering new synthetic drugs. The Drug-Target Interaction (DTI) prediction is a fundamental task in the relationship between proteins and drugs. Unlike encoding proteins by amino acids, we use amino acid subsequence to encode proteins, which simulates the biological process of DTI better. For this research purpose, we proposed a novel deep learning framework based on Bidirectional Encoder Representation from Transformers (BERT), which integrates high-frequency subsequence embedding and transfer learning methods to complete the DTI prediction task. As the first key module, subsequence embedding allows to explore the functional interaction units from drug and protein sequences and then contribute to finding DTI modules. As the second key module, transfer learning promotes the model learn the common DTI features from protein and drug sequences in a large dataset. Overall, the BERT-based model can learn two kinds features through the multi-head self-attention mechanism: internal features of sequence and interaction features of both proteins and drugs, respectively. Compared with other methods, BERT-based methods enable more DTI-related features to be discovered by means of attention scores which associated with tokenized protein/drug subsequences. We conducted extensive experiments for the DTI prediction task on three different benchmark datasets. The experimental results show that the model achieves an average prediction metrics higher than most baseline methods. In order to verify the importance of transfer learning, we conducted an ablation study on datasets, and the results show the superiority of transfer learning. In addition, we test the scalability of the model on the dataset in unseen drugs and proteins, and the results of the experiments show that it is acceptable in scalability.

Drug repurposing for respiratory infections

Respiratory infections such as Coronavirus disease 2019 are a substantial worldwide health challenge, frequently resulting in severe sickness and death, especially in susceptible groups. Conventional drug development for respiratory infections faces obstacles such as extended timescales, substantial expenses, and the rise of resistance to current treatments. Drug repurposing is a potential method that has evolved to quickly find and reuse existing medications for treating respiratory infections. Drug repurposing utilizes medications previously approved for different purposes, providing a cost-effective and time-efficient method to tackle pressing medical needs. This chapter summarizes current progress and obstacles in repurposing medications for respiratory infections, focusing on notable examples of repurposed pharmaceuticals and their probable modes of action. The text also explores the significance of computational approaches, high-throughput screening, and preclinical investigations in identifying potential candidates for repurposing. The text delves into the significance of regulatory factors, clinical trial structure, and actual data in confirming the effectiveness and safety of repurposed medications for respiratory infections. Drug repurposing is a valuable technique for quickly increasing the range of treatments for respiratory infections, leading to better patient outcomes and decreasing the worldwide disease burden.

Prediction of drug-target binding affinity based on deep learning models

The prediction of drug-target binding affinity (DTA) plays an important role in drug discovery. Computerized virtual screening techniques have been used for DTA prediction, greatly reducing the time and economic costs of drug discovery. However, these techniques have not succeeded in reversing the low success rate of new drug development. In recent years, the continuous development of deep learning (DL) technology has brought new opportunities for drug discovery through the DTA prediction. This shift has moved the prediction of DTA from traditional machine learning methods to DL. The DL frameworks used for DTA prediction include convolutional neural networks (CNN), graph convolutional neural networks (GCN), and recurrent neural networks (RNN), and reinforcement learning (RL), among others. This review article summarizes the available literature on DTA prediction using DL models, including DTA quantification metrics and datasets, and DL algorithms used for DTA prediction (including input representation of models, neural network frameworks, valuation indicators, and model interpretability). In addition, the opportunities, challenges, and prospects of the application of DL frameworks for DTA prediction in the field of drug discovery are discussed.

Large-Scale Pretraining Improves Sample Efficiency of Active Learning-Based Virtual Screening

Virtual screening of large compound libraries to identify potential hit candidates is one of the earliest steps in drug discovery. As the size of commercially available compound collections grows exponentially to the scale of billions, active learning and Bayesian optimization have recently been proven as effective methods of narrowing down the search space. An essential component of those methods is a surrogate machine learning model that predicts the desired properties of compounds. An accurate model can achieve high sample efficiency by finding hits with only a fraction of the entire library being virtually screened. In this study, we examined the performance of a pretrained transformer-based language model and graph neural network in a Bayesian optimization active learning framework. The best pretrained model identifies 58.97% of the top-50,000 compounds after screening only 0.6% of an ultralarge library containing 99.5 million compounds, improving 8% over the previous state-of-the-art baseline. Through extensive benchmarks, we show that the superior performance of pretrained models persists in both structure-based and ligand-based drug discovery. Pretrained models can serve as a boost to the accuracy and sample efficiency of active learning-based virtual screening.

Molecular docking and metagenomics assisted mitigation of microplastic pollution

Microplastics, tiny, flimsy, and direct progenitors of principal and subsidiary plastics, cause environmental degradation in aquatic and terrestrial entities. Contamination concerns include irrevocable impacts, potential cytotoxicity, and negative health effects on mortals. The detection, recovery, and degradation strategies of these pollutants in various biota and ecosystems, as well as their impact on plants, animals, and humans, have been a topic of significant interest. But the natural environment is infested with several types of plastics, all having different chemical makeup, structure, shape, and origin. Plastic trash acts as a substrate for microbial growth, creating biofilms on the plastisphere surface. This colonizing microbial diversity can be glimpsed with meta-genomics, a culture-independent approach. Owing to its comprehensive description of microbial communities, genealogical evidence on unconventional biocatalysts or enzymes, genomic correlations, evolutionary profile, and function, it is being touted as one of the promising tools in identifying novel enzymes for the degradation of polymers. Additionally, computational tools such as molecular docking can predict the binding of these novel enzymes to the polymer substrate, which can be validated through in vitro conditions for its environmentally feasible applications. This review mainly deals with the exploration of metagenomics along with computational tools to provide a clearer perspective into the microbial potential in the biodegradation of microplastics. The computational tools due to their polymathic nature will be quintessential in identifying the enzyme structure, binding affinities of the prospective enzymes to the substrates, and foretelling of degradation pathways involved which can be quite instrumental in the furtherance of the plastic degradation studies.

In silico approaches to investigate enzyme immobilization: a comprehensive systematic review

Enzymes are popular catalysts with many applications, especially in industry. Biocatalyst usage on a large scale is facing some limitations, such as low operational stability, low recyclability, and high enzyme cost. Enzyme immobilization is a beneficial strategy to solve these problems. Bioinformatics tools can often correctly predict immobilization outcomes, resulting in a cost-effective experimental phase with the least time consumed. This study provides an overview of in silico methods predicting immobilization processes via a comprehensive systematic review of published articles till 11 December 2022. It also mentions the strengths and weaknesses of the processes and explains the computational analyses in each method that are required for immobilization assessment. In this regard, Web of Science and Scopus databases were screened to gain relevant publications. After screening the gathered documents (n = 3873), 60 articles were selected for the review. The selected papers have applied in silico procedures including only molecular dynamics (MD) simulations (n = 20), parallel tempering Monte Carlo (PTMC) and MD simulations (n = 3), MD and docking (n = 1), density functional theory (DFT) and MD (n = 1), only docking (n = 11), metal ion binding site prediction (MIB) server and docking (n = 2), docking and DFT (n = 1), docking and analysis of enzyme surfaces (n = 1), only DFT (n = 1), only MIB server (n = 2), analysis of an enzyme structure and surface (n = 12), rational design of immobilized derivatives (RDID) software (n = 3), and dissipative particle dynamics (DPD; n = 2). In most included studies (n = 51), enzyme immobilization was investigated experimentally in addition to in silico evaluation.

Computational insights into the stereo-selectivity of catechins for the inhibition of the cancer therapeutic target EGFR kinase

The epidermal growth factor receptor (EGFR) plays a crucial role in regulating cellular growth and survival, and its dysregulation is implicated in various cancers, making it a prime target for cancer therapy. Natural compounds known as catechins have garnered attention as promising anticancer agents. These compounds exert their anticancer effects through diverse mechanisms, primarily by inhibiting receptor tyrosine kinases (RTKs), a protein family that includes the notable member EGFR. Catechins, characterized by two chiral centers and stereoisomerism, demonstrate variations in chemical and physical properties due to differences in the spatial orientation of atoms. Although previous studies have explored the membrane fluidity effects and transport across cellular membranes, the stereo-selectivity of catechins concerning EGFR kinase inhibition remains unexplored. In this study, we investigated the stereo-selectivity of catechins in inhibiting EGFR kinase, both in its wild-type and in the prevalent L858R mutant. Computational analyses indicated that all stereoisomers, including the extensively studied catechin (-)-EGCG, effectively bound within the ATP-binding site, potentially inhibiting EGFR kinase activity. Notably, gallated catechins emerged as superior EGFR inhibitors to their non-gallated counterparts, revealing intriguing binding trends. The top four stereoisomers exhibiting high dock scores and binding energies with wild-type EGFR comprise (-)-CG (-)-GCG (+)-CG, and (-)-EGCG. To assess dynamic behavior and stability, molecular dynamics simulations over 100 ns were conducted for the top-ranked catechin (-)-CG and the widely investigated catechin (-)-EGCG with EGFR kinase. This study enhances our understanding of how the stereoisomeric nature of a drug influences inhibitory potential, providing insights that could guide the selection of specific stereoisomers for improved efficacy inexisting drugs.

Redocking the PDB

Molecular docking is a standard technique in structure-based drug design (SBDD). It aims to predict the 3D structure of a small molecule in the binding site of a receptor (often a protein). Despite being a common technique, it often necessitates multiple tools and involves manual steps. Here, we present the JAMDA preprocessing and docking workflow that is easy to use and allows fully automated docking. We evaluate the JAMDA docking workflow on binding sites extracted from the complete PDB and derive key factors determining JAMDA's docking performance. With that, we try to remove most of the bias due to manual intervention and provide a realistic estimate of the redocking performance of our JAMDA preprocessing and docking workflow for any PDB structure. On this large PDBScan22 data set, our JAMDA workflow finds a pose with an RMSD of at most 2 Å to the crystal ligand on the top rank for 30.1% of the structures. When applying objective structure quality filters to the PDBScan22 data set, the success rate increases to 61.8%. Given the prepared structures from the JAMDA preprocessing pipeline, both JAMDA and the widely used AutoDock Vina perform comparably on this filtered data set (the PDBScan22-HQ data set).

DOCK 6: Incorporating hierarchical traversal through precomputed ligand conformations to enable large-scale docking

To allow DOCK 6 access to unprecedented chemical space for screening billions of small molecules, we have implemented features from DOCK 3.7 into DOCK 6, including a search routine that traverses precomputed ligand conformations stored in a hierarchical database. We tested them on the DUDE-Z and SB2012 test sets. The hierarchical database search routine is 16 times faster than anchor-and-grow. However, the ability of hierarchical database search to reproduce the experimental pose is 16% worse than that of anchor-and-grow. The enrichment performance is on average similar, but DOCK 3.7 has better enrichment than DOCK 6, and DOCK 6 is on average 1.7 times slower. However, with post-docking torsion minimization, DOCK 6 surpasses DOCK 3.7. A large-scale virtual screen is performed with DOCK 6 on 23 million fragment molecules. We use current features in DOCK 6 to complement hierarchical database calculations, including torsion minimization, which is not available in DOCK 3.7.

Sulfonamide Derivatives: Recent Compounds with Potent Anti-alzheimer's Disease Activity

Facile synthetic procedures and broad spectrum of biological activities are special attributes of sulfonamides. Sulfonamide derivatives have demonstrated potential as a class of compounds for the treatment of Alzheimer's disease (AD). Recent sulfonamide derivatives have been reported as prospective anti-AD agents, with a focus on analogues that significantly inhibit the function of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) enzymes and exhibit remarkable antioxidant and anti-inflammatory properties, all of which are critical for the treatment of AD. Sulfonamide- mediated activation of nuclear factor erythroid 2-related factor 2 (NRF2), a key regulator of the endogenous antioxidant response, has also been suggested as a potential therapeutic approach in AD. Additionally, it has been discovered that a number of sulfonamide derivatives show selectivity for the β- and γ-secretase enzymes and a significant reduction of amyloid B (Aβ) aggregation, which have been implicated in AD. The comparative molecular docking of benzenesulfonamide and donepezil, an AD reference drug showed comparable anti-AD activities. These suggest that sulfonamide derivatives may represent a new class of drugs for the treatment of AD. Thus, the current review will focus on recent studies on the chemical synthesis and evaluation of the anti-AD properties, molecular docking, pharmacological profile, and structure-activity relationship (SAR) of sulfonamide derivatives, as well as their potential anti-AD mechanisms of action. This paper offers a thorough assessment of the state of the art in this field of study and emphasizes the potential of sulfonamide derivatives synthesized during the 2012-2023 period as a new class of compounds for the treatment of AD.

Exploring Scoring Function Space: Developing Computational Models for Drug Discovery

BACKGROUND

The idea of scoring function space established a systems-level approach to address the development of models to predict the affinity of drug molecules by those interested in drug discovery.

OBJECTIVE

Our goal here is to review the concept of scoring function space and how to explore it to develop machine learning models to address protein-ligand binding affinity.

METHODS

We searched the articles available in PubMed related to the scoring function space. We also utilized crystallographic structures found in the protein data bank (PDB) to represent the protein space.

RESULTS

The application of systems-level approaches to address receptor-drug interactions allows us to have a holistic view of the process of drug discovery. The scoring function space adds flexibility to the process since it makes it possible to see drug discovery as a relationship involving mathematical spaces.

CONCLUSION

The application of the concept of scoring function space has provided us with an integrated view of drug discovery methods. This concept is useful during drug discovery, where we see the process as a computational search of the scoring function space to find an adequate model to predict receptor-drug binding affinity.

Computer aided therapeutic tripeptide design, in alleviating the pathogenic proclivities of nocuous α-synuclein fibrils

Parkinson's disorder (PD) exacerbates neuronal degeneration of motor nerves, thereby effectuating uncoordinated movements and tremors. Aberrant alpha-synuclein (α-syn) is culpable of triggering PD, wherein cytotoxic amyloid aggregates of α-syn get deposited in motor neurons to instigate neuro-degeneration. Amyloid aggregates, typically rich in beta sheets are cardinal targets to mitigate their neurotoxic effects. In this analysis, owing to their interaction specificity, we formulated an efficacious tripeptide out of the aggregation-prone region of α-syn protein. With the help of a proficient computational pipeline, systematic peptide shortening and an adept molecular simulation platform, we formulated a tripeptide, VAV from α-syn structure based hexapeptide KISVRV. Indeed, the VAV tripeptide was able to effectively mitigate the α-syn amyloid fibrils' dynamic rate of beta-sheet formation. Additional trajectory analyses of the VAV- α-syn complex indicated that, upon its dynamic interaction, VAV efficiently altered the distinct pathogenic structural dynamics of α-syn, further advocating its potential in alleviating aberrant α-syn's amyloidogenic proclivities. Consistent findings from various computational analyses have led us to surmise that VAV could potentially re-alter the pathogenic conformational orientation of α-syn, essential to mitigate its cytotoxicity. Hence, VAV tripeptide could be an efficacious therapeutic candidate to efficiently ameliorate aberrant α-syn amyloid mediated neurotoxicity, eventually attenuating the nocuous effects of PD.Communicated by Ramaswamy H. Sarma.

Computer-Aided Drug Design and Drug Discovery: A Prospective Analysis

In the dynamic landscape of drug discovery, Computer-Aided Drug Design (CADD) emerges as a transformative force, bridging the realms of biology and technology. This paper overviews CADDs historical evolution, categorization into structure-based and ligand-based approaches, and its crucial role in rationalizing and expediting drug discovery. As CADD advances, incorporating diverse biological data and ensuring data privacy become paramount. Challenges persist, demanding the optimization of algorithms and robust ethical frameworks. Integrating Machine Learning and Artificial Intelligence amplifies CADDs predictive capabilities, yet ethical considerations and scalability challenges linger. Collaborative efforts and global initiatives, exemplified by platforms like Open-Source Malaria, underscore the democratization of drug discovery. The convergence of CADD with personalized medicine offers tailored therapeutic solutions, though ethical dilemmas and accessibility concerns must be navigated. Emerging technologies like quantum computing, immersive technologies, and green chemistry promise to redefine the future of CADD. The trajectory of CADD, marked by rapid advancements, anticipates challenges in ensuring accuracy, addressing biases in AI, and incorporating sustainability metrics. This paper concludes by highlighting the need for proactive measures in navigating the ethical, technological, and educational frontiers of CADD to shape a healthier, brighter future in drug discovery.

Analysis of some flavonoids for inhibitory mechanism against cancer target phosphatidylinositol 3-kinase (PI3K) using computational tool

Cancer has been one of the leading causes of mortality worldwide over the past few years. Some progress has been made in the development of more effective cancer therapeutics, resulting in improved survival rates. However, the desired outcome in the form of successful treatment is yet to be achieved. There is high demand for the development of innovative, inexpensive, and effective anticancer treatments using natural resources. Natural compounds have been increasingly discovered and used for cancer therapy owing to their high molecular diversity, novel biofunctionality, and minimal side effects. These compounds can be utilized as chemopreventive agents because they can efficiently inhibit cell growth, control cell cycle progression, and block several tumor-promoting signaling pathways. PI3K is an important upstream protein of the PI3K-Akt-mTOR pathway and a well-established cancer therapeutic target. This study aimed to explore the small molecules, natural flavonoids, viz. quercetin, luteolin, kaempferol, genistein, wogonin, daidzein, and flavopiridol for PI3Kγ kinase activity inhibition. In this study, the binding pose, interacting residues, molecular interactions, binding energies, and dissociation constants were investigated. Our results showed that these flavonoids bound well with PI3Kγ with adequate binding strength scores and binding energy ranging from (-8.19 to -8.97 Kcal/mol). Among the explored ligands, flavopiridol showed the highest binding energy of -8.97 Kcal/mol, dock score (-44.40), and dissociation constant term, p K d of 6.58 against PI3Kγ. Based on the above results, the stability of the most promising ligand, flavopiridol, against PI3Kγ was evaluated by molecular dynamics simulations for 200 ns, confirming the stable flavopiridol and PI3Kγ complex. Our study suggests that among the selected flavonoids specifically flavopiridol may act as potential inhibitors of PI3Kγ and could be a therapeutic alternative to inhibit the PI3Kγ pathway, providing new insights into rational drug discovery research for cancer therapy.

Novel multi-objective affinity approach allows to identify pH-specific μ-opioid receptor agonists

Opioids are essential pharmaceuticals due to their analgesic properties, however, lethal side effects, addiction, and opioid tolerance are extremely challenging. The development of novel molecules targeting the [Formula: see text]-opioid receptor (MOR) in inflamed, but not in healthy tissue, could significantly reduce these unwanted effects. Finding such novel molecules can be achieved by maximizing the binding affinity to the MOR at acidic pH while minimizing it at neutral pH, thus combining two conflicting objectives. Here, this multi-objective optimal affinity approach is presented, together with a virtual drug discovery pipeline for its practical implementation. When applied to finding pH-specific drug candidates, it combines protonation state-dependent structure and ligand preparation with high-throughput virtual screening. We employ this pipeline to characterize a set of MOR agonists identifying a morphine-like opioid derivative with higher predicted binding affinities to the MOR at low pH compared to neutral pH. Our results also confirm existing experimental evidence that NFEPP, a previously described fentanyl derivative with reduced side effects, and recently reported [Formula: see text]-fluorofentanyls and -morphines show an increased specificity for the MOR at acidic pH when compared to fentanyl and morphine. We further applied our approach to screen a >50K ligand library identifying novel molecules with pH-specific predicted binding affinities to the MOR. The presented differential docking pipeline can be applied to perform multi-objective affinity optimization to identify safer and more specific drug candidates at large scale.

Recent Developments in Ultralarge and Structure-Based Virtual Screening Approaches

Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.

Docking covalent targets for drug discovery: stimulating the computer-aided drug design community of possible pitfalls and erroneous practices

The continuous approval of covalent drugs in recent years for the treatment of diseases has led to an increased search for covalent agents by medicinal chemists and computational scientists worldwide. In the computational parlance, molecular docking which is a popular tool to investigate the interaction of a ligand and a protein target, does not account for the formation of covalent bond, and the increasing application of these conventional programs to covalent targets in early drug discovery practice is a matter of utmost concern. Thus, in this comprehensive review, we sought to educate the docking community about the realization of covalent docking and the existence of suitable programs to make their future virtual-screening events on covalent targets worthwhile and scientifically rational. More interestingly, we went beyond the classical description of the functionality of covalent-docking programs down to selecting the 'best' program to consult with during a virtual-screening campaign based on receptor class and covalent warhead chemistry. In addition, we made a highlight on how covalent docking could be achieved using random conventional docking software. And lastly, we raised an alert on the growing erroneous molecular docking practices with covalent targets. Our aim is to guide scientists in the rational docking pursuit when dealing with covalent targets, as this will reduce false-positive results and also increase the reliability of their work for translational research.

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

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

Computer-aided drug design in seeking viral capsid modulators

Approved or licensed antiviral drugs have limited applications because of their drug resistance and severe adverse effects. By contrast, by stabilizing or destroying the viral capsid, compounds known as capsid modulators prevent viral replication by acting on new targets and, therefore, overcoming the problem of clinical drug resistance. For example. computer-aided drug design (CADD) methods, using strategies based on structures of biological targets (structure-based drug design; SBDD), such as docking, molecular dynamics (MD) simulations, and virtual screening (VS), have provided opportunities for fast and effective development of viral capsid modulators. In this review, we summarize the application of CADD in the discovery, optimization, and mechanism prediction of capsid-targeting small molecules, providing new insights into antiviral drug discovery modalities. Teaser: Computer-aided drug design will accelerate the development of viral capsid regulators, which brings new hope for the treatment of refractory viral diseases.

In Silico Methodologies to Improve Antioxidants' Characterization from Marine Organisms

Marine organisms have been reported to be valuable sources of bioactive molecules that have found applications in different industrial fields. From organism sampling to the identification and bioactivity characterization of a specific compound, different steps are necessary, which are time- and cost-consuming. Thanks to the advent of the -omic era, numerous genome, metagenome, transcriptome, metatranscriptome, proteome and microbiome data have been reported and deposited in public databases. These advancements have been fundamental for the development of in silico strategies for basic and applied research. In silico studies represent a convenient and efficient approach to the bioactivity prediction of known and newly identified marine molecules, reducing the time and costs of "wet-lab" experiments. This review focuses on in silico approaches applied to bioactive molecule discoveries from marine organisms. When available, validation studies reporting a bioactivity assay to confirm the presence of an antioxidant molecule or enzyme are reported, as well. Overall, this review suggests that in silico approaches can offer a valuable alternative to most expensive approaches and proposes them as a little explored field in which to invest.

Design and Development of Novel Glitazones for Activation of PGC-1α Signaling Via PPAR-γ Agonism: A Promising Therapeutic Approach against Parkinson's Disease

Herein, we rationally designed and developed two novel glitazones (G1 and G2) to target peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PGC-1α) signaling through peroxisome proliferator-activated receptors (PPAR)-γ agonism as a therapeutic for Parkinson's disease (PD). The synthesized molecules were analyzed by mass spectrometry and NMR spectroscopy. The neuroprotective functionality of the synthesized molecules was assessed by a cell viability assay in lipopolysaccharide-intoxicated SHSY5Y neuroblastoma cell lines. The ability of these new glitazones to scavenge free radicals was further ascertained via a lipid peroxide assay, and pharmacokinetic properties were verified using in silico absorption, distribution, metabolism, excretion, and toxicity analyses. The molecular docking reports recognized the mode of interaction of the glitazones with PPAR-γ. The G1 and G2 exhibited a noticeable neuroprotective effect in lipopolysaccharide-intoxicated SHSY5Y neuroblastoma cells with the half-maximal inhibitory concentration value of 2.247 and 4.509 μM, respectively. Both test compounds prevented 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine-induced motor impairment in mice, as demonstrated by the beam walk test. Further, treating the diseased mice with G1 and G2 resulted in significant restoration of antioxidant enzymes glutathione and superoxide and reduced the intensity of lipid peroxidation inside the brain tissues. Histopathological analysis of the glitazones-treated mice brain revealed a reduced apoptotic region and a rise in the number of viable pyramidal neurons and oligodendrocytes. The study concluded that G1 and G2 showed promising results in treating PD by activating PGC-1α signaling in brain via PPAR-γ agonism. However, more extensive research is necessary for a better understanding of functional targets and signaling pathways.

Opportunities and challenges in application of artificial intelligence in pharmacology

Artificial intelligence (AI) is a machine science that can mimic human behaviour like intelligent analysis of data. AI functions with specialized algorithms and integrates with deep and machine learning. Living in the digital world can generate a huge amount of medical data every day. Therefore, we need an automated and reliable evaluation tool that can make decisions more accurately and faster. Machine learning has the potential to learn, understand and analyse the data used in healthcare systems. In the last few years, AI is known to be employed in various fields in pharmaceutical science especially in pharmacological research. It helps in the analysis of preclinical (laboratory animals) and clinical (in human) trial data. AI also plays important role in various processes such as drug discovery/manufacturing, diagnosis of big data for disease identification, personalized treatment, clinical trial research, radiotherapy, surgical robotics, smart electronic health records, and epidemic outbreak prediction. Moreover, AI has been used in the evaluation of biomarkers and diseases. In this review, we explain various models and general processes of machine learning and their role in pharmacological science. Therefore, AI with deep learning and machine learning could be relevant in pharmacological research.

Discovery of new PKN2 inhibitory chemotypes via QSAR-guided selection of docking-based pharmacophores

Serine/threonine-protein kinase N2 (PKN2) plays an important role in cell cycle progression, cell migration, cell adhesion and transcription activation signaling processes. In cancer, however, it plays important roles in tumor cell migration, invasion and apoptosis. PKN2 inhibitors have been shown to be promising in treating cancer. This prompted us to model this interesting target using our QSAR-guided selection of docking-based pharmacophores approach where numerous pharmacophores are extracted from docked ligand poses and allowed to compete within the context of QSAR. The optimal pharmacophore was sterically-refined, validated by receiver operating characteristic (ROC) curve analysis and used as virtual search query to screen the National Cancer Institute (NCI) database for new promising anti-PKN2 leads of novel chemotypes. Three low micromolar hits were identified with IC₅₀ values ranging between 9.9 and 18.6 µM. Pharmacological assays showed promising cytotoxic properties for active hits in MTT and wound healing assays against MCF-7 and PANC-1 cancer cells.

CADD, AI and ML in drug discovery: A comprehensive review

Computer-aided drug design (CADD) is an emerging field that has drawn a lot of interest because of its potential to expedite and lower the cost of the drug development process. Drug discovery research is expensive and time-consuming, and it frequently took 10-15 years for a drug to be commercially available. CADD has significantly impacted this area of research. Further, the combination of CADD with Artificial Intelligence (AI), Machine Learning (ML), and Deep Learning (DL) technologies to handle enormous amounts of biological data has reduced the time and cost associated with the drug development process. This review will discuss how CADD, AI, ML, and DL approaches help identify drug candidates and various other steps of the drug discovery process. It will also provide a detailed overview of the different in silico tools used and how these approaches interact.

HCovDock: an efficient docking method for modeling covalent protein-ligand interactions

Covalent inhibitors have received extensive attentions in the past few decades because of their long residence time, high binding efficiency and strong selectivity. Therefore, it is valuable to develop computational tools like molecular docking for modeling of covalent protein-ligand interactions or screening of potential covalent drugs. Meeting the needs, we have proposed HCovDock, an efficient docking algorithm for covalent protein-ligand interactions by integrating a ligand sampling method of incremental construction and a scoring function with covalent bond-based energy. Tested on a benchmark containing 207 diverse protein-ligand complexes, HCovDock exhibits a significantly better performance than seven other state-of-the-art covalent docking programs (AutoDock, Cov_DOX, CovDock, FITTED, GOLD, ICM-Pro and MOE). With the criterion of ligand root-mean-squared distance < 2.0 Å, HCovDock obtains a high success rate of 70.5% and 93.2% in reproducing experimentally observed structures for top 1 and top 10 predictions. In addition, HCovDock is also validated in virtual screening against 10 receptors of three proteins. HCovDock is computationally efficient and the average running time for docking a ligand is only 5 min with as fast as 1 sec for ligands with one rotatable bond and about 18 min for ligands with 23 rotational bonds. HCovDock can be freely assessed at http://huanglab.phys.hust.edu.cn/hcovdock/.

Neural Networks in the Design of Molecules with Affinity to Selected Protein Domains

Drug design with machine learning support can speed up new drug discoveries. While current databases of known compounds are smaller in magnitude (approximately 108), the number of small drug-like molecules is estimated to be between 1023 and 1060. The use of molecular docking algorithms can help in new drug development by sieving out the worst drug-receptor complexes. New chemical spaces can be efficiently searched with the application of artificial intelligence. From that, new structures can be proposed. The research proposed aims to create new chemical structures supported by a deep neural network that will possess an affinity to the selected protein domains. Transferring chemical structures into SELFIES codes helped us pass chemical information to a neural network. On the basis of vectorized SELFIES, new chemical structures can be created. With the use of the created neural network, novel compounds that are chemically sensible can be generated. Newly created chemical structures are sieved by the quantitative estimation of the drug-likeness descriptor, Lipinski's rule of 5, and the synthetic Bayesian accessibility classifier score. The affinity to selected protein domains was verified with the use of the AutoDock tool. As per the results, we obtained the structures that possess an affinity to the selected protein domains, namely PDB IDs 7NPC, 7NP5, and 7KXD.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Innovative microbial disease biocontrol strategies mediated by quorum quenching and their multifaceted applications: A review

With the increasing resistance exhibited by undesirable bacteria to traditional antibiotics, the need to discover alternative (or, at least, supplementary) treatments to combat chemically resistant bacteria is becoming urgent. Quorum sensing (QS) refers to a novel bacterial communication system for monitoring cell density and regulation of a network of gene expression that is mediated by a group of signaling molecules called autoinducers (AIs). QS-regulated multicellular behaviors include biofilm formation, horizontal gene transfer, and antibiotic synthesis, which are demonstrating increasing pathogenicity to plants and aquacultural animals as well as contamination of wastewater treatment devices. To inhibit QS-regulated microbial behaviors, the strategy of quorum quenching (QQ) has been developed. Different quorum quenchers interfere with QS through different mechanisms, such as competitively inhibiting AI perception (e.g., by QS inhibitors) and AI degradation (e.g., by QQ enzymes). In this review, we first introduce different signaling molecules, including diffusible signal factor (DSF) and acyl homoserine lactones (AHLs) for Gram-negative bacteria, AIPs for Gram-positive bacteria, and AI-2 for interspecies communication, thus demonstrating the mode of action of the QS system. We next exemplify the QQ mechanisms of various quorum quenchers, such as chemical QS inhibitors, and the physical/enzymatic degradation of QS signals. We devote special attention to AHL-degrading enzymes, which are categorized in detail according to their diverse catalytic mechanisms and enzymatic properties. In the final part, the applications and advantages of quorum quenchers (especially QQ enzymes and bacteria) are summarized in the context of agricultural/aquacultural pathogen biocontrol, membrane bioreactors for wastewater treatment, and the attenuation of human pathogenic bacteria. Taken together, we present the state-of-the-art in research considering QS and QQ, providing theoretical evidence and support for wider application of this promising environmentally friendly biocontrol strategy.

Prediction of adsorption capacity and biodegradability of polybrominated diphenyl ethers in soil

Polybrominated diphenyl ethers (PBDEs) are widely used brominated flame retardants with strong toxicity concerns. Understanding the behaviors of PBDEs in soil is essential to evaluate their environmental impact. However, the limited, incoherent, and inaccurate data has challenged predicting the adsorption capacity and biodegradability of all 209 PBDE congeners in the soil. Moreover, there are minimal studies regarding the interactions between adsorption and biodegradation behaviors of PBDEs in the soil. Herein, in this study, we adopted quantitative structure-property relationship (QSAR) modeling to predict the adsorption behavior of 209 PBDE congeners by estimating their organic carbon-water partition coefficient (K_OC) values. In addition, the biodegradability of commonly occurring PBDE congeners was evaluated by analyzing their affinity to extracellular enzymes responsible for biodegradation using molecular docking. The results highlight that the degree of bromination plays a significant role in both the absorption and biodegradation of PBDEs in the soil due to compound stability and molecular geometry. Our findings help to advance the knowledge on PBDE behaviors in the soil and facilitate PBDE remediation associated with a soil environment.

Computational Approaches to Designing Antiviral Drugs against COVID-19: A Comprehensive Review

The global impact of the COVID-19 pandemic caused by SARS-CoV-2 necessitates innovative strategies for the rapid development of effective treatments. Computational methodologies, such as molecular modelling, molecular dynamics simulations, and artificial intelligence, have emerged as indispensable tools in the drug discovery process. This review aimed to provide a comprehensive overview of these computational approaches and their application in the design of antiviral agents for COVID-19. Starting with an examination of ligand-based and structure-based drug discovery, the review has delved into the intricate ways through which molecular modelling can accelerate the identification of potential therapies. Additionally, the investigation extends to phytochemicals sourced from nature, which have shown promise as potential antiviral agents. Noteworthy compounds, including gallic acid, naringin, hesperidin, Tinospora cordifolia, curcumin, nimbin, azadironic acid, nimbionone, nimbionol, and nimocinol, have exhibited high affinity for COVID-19 Mpro and favourable binding energy profiles compared to current drugs. Although these compounds hold potential, their further validation through in vitro and in vivo experimentation is imperative. Throughout this exploration, the review has emphasized the pivotal role of computational biologists, bioinformaticians, and biotechnologists in driving rapid advancements in clinical research and therapeutic development. By combining state-of-the-art computational techniques with insights from structural and molecular biology, the search for potent antiviral agents has been accelerated. The collaboration between these disciplines holds immense promise in addressing the transmissibility and virulence of SARS-CoV-2.

Organotin Antifouling Compounds and Sex-Steroid Nuclear Receptor Perturbation: Some Structural Insights

Organotin compounds (OTCs) are a commercially important group of organometallic compounds of tin used globally as polyvinyl chloride stabilizers and marine antifouling biocides. Worldwide use of OTCs has resulted in their ubiquitous presence in ecosystems across all the continents. OTCs have metabolic and endocrine disrupting effects in marine and terrestrial organisms. Thus, harmful OTCs (tributyltin) have been banned by the International Convention on the Control of Harmful Antifouling Systems since 2008. However, continued manufacturing by non-member countries poses a substantial risk for animal and human health. In this study, structural binding of common commercial OTCs, tributyltin (TBT), dibutyltin (DBT), monobutyltin (MBT), triphenyltin (TPT), diphenyltin (DPT), monophenyltin (MPT), and azocyclotin (ACT) against sex-steroid nuclear receptors, androgen receptor (AR), and estrogen receptors (ERα, ERβ) was performed using molecular docking and MD simulation. TBT, DBT, DPT, and MPT bound deep within the binding sites of AR, ERα, and Erβ, showing good dock score, binding energy and dissociation constants that were comparable to bound native ligands, testosterone and estradiol. The stability of docking complex was shown by MD simulation of organotin/receptor complex with RMSD, RMSF, Rg, and SASA plots showing stable interaction, low deviation, and compactness of the complex. A high commonality (50-100%) of interacting residues of ERα and ERβ for the docked ligands and bound native ligand (estradiol) indicated that the organotin compounds bound in the same binding site of the receptor as the native ligand. The results suggested that organotins may interfere with the natural steroid/receptor binding and perturb steroid signaling.

Potential Disruption of Systemic Hormone Transport by Tobacco Alkaloids Using Computational Approaches

Tobacco/nicotine is one of the most toxic and addictive substances and continues to pose a significant threat to global public health. The harmful effects of smoking/nicotine affect every system in the human body. Nicotine has been associated with effects on endocrine homeostasis in humans such as the imbalance of gonadal steroid hormones, adrenal corticosteroid hormones, and thyroid hormones. The present study was conducted to characterize the structural binding interactions of nicotine and its three important metabolites, cotinine, trans-3'-hydroxycotinine, and 5'-hydroxycotinine, against circulatory hormone carrier proteins, i.e., sex-hormone-binding globulin (SHBG), corticosteroid-binding globulin (CBG), and thyroxine-binding globulin (TBG). Nicotine and its metabolites formed nonbonded contacts and/or hydrogen bonds with amino acid residues of the carrier proteins. For SHBG, Phe-67 and Met-139 were the most important amino acid residues for nicotine ligand binding showing the maximum number of interactions and maximum loss in ASA. For CBG, Trp-371 and Asn-264 were the most important amino acid residues, and for TBG, Ser-23, Leu-269, Lys-270, Asn-273, and Arg-381 were the most important amino acid residues. Most of the amino acid residues of carrier proteins interacting with nicotine ligands showed a commonality with the interacting residues for the native ligands of the proteins. Taken together, the results suggested that nicotine and its three metabolites competed with native ligands for binding to their carrier proteins. Thus, nicotine and its three metabolites may potentially interfere with the binding of testosterone, estradiol, cortisol, progesterone, thyroxine, and triiodothyronine to their carrier proteins and result in the disbalance of their transport and homeostasis in the blood circulation.

A molecular evolution algorithm for ligand design in DOCK

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Implications of the Essential Role of Small Molecule Ligand Binding Pockets in Protein-Protein Interactions

Protein-protein interactions (PPIs) and protein-metabolite interactions play a key role in many biochemical processes, yet they are often viewed as being independent. However, the fact that small molecule drugs have been successful in inhibiting PPIs suggests a deeper relationship between protein pockets that bind small molecules and PPIs. We demonstrate that 2/3 of PPI interfaces, including antibody-epitope interfaces, contain at least one significant small molecule ligand binding pocket. In a representative library of 50 distinct protein-protein interactions involving hundreds of mutations, >75% of hot spot residues overlap with small molecule ligand binding pockets. Hence, ligand binding pockets play an essential role in PPIs. In representative cases, evolutionary unrelated monomers that are involved in different multimeric interactions yet share the same pocket are predicted to bind the same metabolites/drugs; these results are confirmed by examples in the PDB. Thus, the binding of a metabolite can shift the equilibrium between monomers and multimers. This implicit coupling of PPI equilibria, termed "metabolic entanglement", was successfully employed to suggest novel functional relationships among protein multimers that do not directly interact. Thus, the current work provides an approach to unify metabolomics and protein interactomics.

Docking cyclic peptides formed by a disulfide bond through a hierarchical strategy

MOTIVATION

Cyclization is a common strategy to enhance the therapeutic potential of peptides. Many cyclic peptide drugs have been approved for clinical use, in which the disulfide-driven cyclic peptide is one of the most prevalent categories. Molecular docking is a powerful computational method to predict the binding modes of molecules. For protein-cyclic peptide docking, a big challenge is considering the flexibility of peptides with conformers constrained by cyclization.

RESULTS

Integrating our efficient peptide 3D conformation sampling algorithm MODPEP2.0 and knowledge-based scoring function ITScorePP, we have proposed an extended version of our hierarchical peptide docking algorithm, named HPEPDOCK2.0, to predict the binding modes of the peptide cyclized through a disulfide against a protein. Our HPEPDOCK2.0 approach was extensively evaluated on diverse test sets and compared with the state-of-the-art cyclic peptide docking program AutoDock CrankPep (ADCP). On a benchmark dataset of 18 cyclic peptide-protein complexes, HPEPDOCK2.0 obtained a native contact fraction of above 0.5 for 61% of the cases when the top prediction was considered, compared with 39% for ADCP. On a larger test set of 25 cyclic peptide-protein complexes, HPEPDOCK2.0 yielded a success rate of 44% for the top prediction, compared with 20% for ADCP. In addition, HPEPDOCK2.0 was also validated on two other test sets of 10 and 11 complexes with apo and predicted receptor structures, respectively. HPEPDOCK2.0 is computationally efficient and the average running time for docking a cyclic peptide is about 34 min on a single CPU core, compared with 496 min for ADCP. HPEPDOCK2.0 will facilitate the study of the interaction between cyclic peptides and proteins and the development of therapeutic cyclic peptide drugs.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/hpepdock/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of novel inhibitors for mycobacterial polyketide synthase 13 via in silico drug screening assisted by the parallel compound screening with genetic algorithm-based programs

Identifying small compounds capable of inhibiting Mycobacterium tuberculosis polyketide synthase 13 (Pks13), in charge of final step of mycolic acid biosynthesis, could lead to the development of a novel antituberculosis drug. This study screened for lead compounds capable of targeting M. tuberculosis Pks13 from a chemical library comprising 154,118 compounds through multiple in silico docking simulations. The parallel compound screening (PCS), conducted via two genetic algorithm-based programs was applied in the screening strategy. Out of seven experimentally validated compounds, four compounds showed inhibitory effects on the growth of the model mycobacteria (Mycobacterium smegmatis). Subsequent docking simulation of analogs of the promising leads with the assistance of PCS resulted in the identification of three additional compounds with potent antimycobacterial effects (compounds A1, A2, and A5). Further, molecular dynamics simulation predicted stable interaction between M. tuberculosis Pks13 active site and compound A2, which showed potent antimycobacterial activity comparable to that of isoniazid. The present study demonstrated the efficacy of in silico structure-based drug screening through PCS in antituberculosis drug discovery.

Identification of dual-targeted Mycobacterium tuberculosis aminoacyl-tRNA synthetase inhibitors using machine learning

Background: The most serious challenge in the treatment of tuberculosis is the multidrug resistance of Mycobacterium tuberculosis to existing antibiotics. As a strategy to overcome resistance we used a multitarget drug design approach. The purpose of the work was to discover dual-targeted inhibitors of mycobacterial LeuRS and MetRS with machine learning. Methods: The artificial neural networks were built using module nnet from R 3.6.1. The inhibitory activity of compounds toward LeuRS and MetRS was investigated in aminoacylation assays. Results: Using a machine-learning approach, we identified dual-targeted inhibitors of LeuRS and MetRS among 2-(quinolin-2-ylsulfanyl)-acetamide derivatives. The most active compound inhibits MetRS and LeuRS with IC₅₀ values of 33 μm and 23.9 μm, respectively. Conclusion: 2-(Quinolin-2-ylsulfanyl)-acetamide scaffold can be useful for further research.

Drug Discovery by Automated Adaptation of Chemical Structure and Identity

Computer-aided drug design offers the potential to dramatically reduce the cost and effort required for drug discovery. While screening-based methods are valuable in the early stages of hit identification, they are frequently succeeded by iterative, hypothesis-driven computations that require recurrent investment of human time and intuition. To increase automation, we introduce a computational method for lead refinement that combines concerted dynamics of the ligand/protein complex via molecular dynamics simulations with integrated Monte Carlo-based changes in the chemical formula of the ligand. This approach, which we refer to as ligand-exchange Monte Carlo molecular dynamics, accounts for solvent- and entropy-based contributions to competitive binding free energies by coupling the energetics of bound and unbound states during the ligand-exchange attempt. Quantitative comparison of relative binding free energies to reference values from free energy perturbation, conducted in vacuum, indicates that ligand-exchange Monte Carlo molecular dynamics simulations sample relevant conformational ensembles and are capable of identifying strongly binding compounds. Additional simulations demonstrate the use of an implicit solvent model. We speculate that the use of chemical graphs in which exchanges are only permitted between ligands with sufficient similarity may enable an automated search to capture some of the benefits provided by human intuition during hypothesis-guided lead refinement.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Scoring Functions for Protein-Ligand Binding Affinity Prediction using Structure-Based Deep Learning: A Review

The rapid and accurate in silico prediction of protein-ligand binding free energies or binding affinities has the potential to transform drug discovery. In recent years, there has been a rapid growth of interest in deep learning methods for the prediction of protein-ligand binding affinities based on the structural information of protein-ligand complexes. These structure-based scoring functions often obtain better results than classical scoring functions when applied within their applicability domain. Here we review structure-based scoring functions for binding affinity prediction based on deep learning, focussing on different types of architectures, featurization strategies, data sets, methods for training and evaluation, and the role of explainable artificial intelligence in building useful models for real drug-discovery applications.

Maximizing the integration of virtual and experimental screening in hit discovery

INTRODUCTION

Experimental and virtual screening contributes to the discovery of more than 50% of clinical candidates. Considering the similar concept and goals, early-phase drug discovery would benefit from the effective integration of these approaches.

AREAS COVERED

After reviewing the recent trends in both experimental and virtual screening, the authors discuss different integration strategies from parallel, focused, sequential, and iterative screening. Strategic considerations are demonstrated in a number of real-life case studies.

EXPERT OPINION

Experimental and virtual screening are complementary approaches that should be integrated in lead discovery settings. Virtual screening can access extremely large synthetically feasible chemical space that can be effectively searched on GPU clusters or cloud architectures. Experimental screening provides reliable datasets by quantitative HTS applications, and DNA-encoded libraries (DEL) have enlarged the chemical space covered by these technologies. These developments, together with the use of artificial intelligence methods, represent new options for their efficient integration. The case studies discussed here demonstrate the benefits of complementary strategies, such as focused and iterative screening.

Machine-Learning- and Knowledge-Based Scoring Functions Incorporating Ligand and Protein Fingerprints

We propose a novel machine-learning-based scoring function for drug discovery that incorporates ligand and protein structural information into a knowledge-based PMF score. Molecular docking, a simulation method for structure-based drug design (SBDD), is expected to reduce the enormous costs associated with conventional experimental methods in terms of rational drug discovery. Molecular docking has two main purposes: to predict ligand-binding structures for target proteins and to predict protein-ligand binding affinity. Currently available programs of molecular docking offer an accurate prediction of ligand binding structures for many systems. However, the accurate prediction of binding affinity remains challenging. In this study, we developed a new scoring function that incorporates fingerprints representing ligand and protein structures as descriptors in the PMF score. Here, regression analysis of the scoring function was performed using the following machine learning techniques: least absolute shrinkage and selection operator (LASSO) and light gradient boosting machine (LightGBM). The results on a test data set showed that the binding affinity delivered by the newly developed scoring function has a Pearson correlation coefficient of 0.79 with the experimental value, which surpasses that of the conventional scoring functions. Further analysis provided a chemical understanding of the descriptors that contributed significantly to the improvement in prediction accuracy. Our approach and findings are useful for rational drug discovery.

Consensus scoring evaluated using the GPCR-Bench dataset: Reconsidering the role of MM/GBSA

The recent availability of large numbers of GPCR crystal structures has provided an unprecedented opportunity to evaluate their performance in virtual screening protocols using established benchmarking datasets. In this study, we evaluated the ability of MM/GBSA in consensus scoring-based virtual screening enrichment together with nine classical scoring functions, using the GPCR-Bench dataset consisting of 24 GPCR crystal structures and 254,646 actives and decoys. While the performance of consensus scoring was modest overall, combinations which included MM/GBSA performed relatively well compared to combinations of classical scoring functions. Combinations of MM/GBSA and good-performing scoring functions provided the highest proportion of improvements, with improvements observed in 32% and 19% of all combinations across all targets at the EF_1% and EF_5% levels respectively. Combinations of MM/GBSA and poor-performing scoring functions still outperformed classical scoring functions, with improvements observed in 26% and 17% of all combinations at the EF_1% and EF_5% levels. In comparison, only 14-22% and 6-11% of combinations of classical scoring functions produced improvements at EF_1% and EF_5% respectively. Efforts to improve performance by increasing the number of scoring functions in consensus scoring to three were mostly ineffective. We also observed that consensus scoring performed better for individual scoring functions possessing initially low enrichment factors, potentially implying their benefits are more relevant in such scenarios. Overall, this study demonstrated the first implementation of MM/GBSA in consensus scoring using the GPCR-Bench dataset and could provide a valuable benchmark of the performance of MM/GBSA in comparison to classical scoring functions in consensus scoring for GPCRs.

Validation of an allosteric binding site of Src kinase identified by unbiased ligand binding simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often "cryptic" without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a v-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein-ligand systems.