Classification of Current Scoring Functions

Spectroscopic, electrochemical, and molecular docking studies of the interaction between the antihistamine drug desloratadine and dsDNA

Through the utilization of fluorescence spectroscopy, electrochemical, and molecular docking methods, this research investigates the interaction between the antihistamine drug desloratadine and calf thymus double-stranded DNA (ct-dsDNA). Deoxyguanosine (dGuo) and deoxyadenosine (dAdo) oxidation signals were diminished by incubation with varying concentrations of desloratadine, as determined by differential pulse voltammetry (DPV). This change was ascribed to desloratadine's binding mechanism to ct-dsDNA. The binding constant (Kb) between desloratadine and ct-dsDNA was determined to be 2.2 × 105 M-1 throughout electrochemical experiments. In order to further develop our comprehension of the interaction mechanism between desloratadine and ct-dsDNA, a series of spectroscopic experiments and molecular docking simulations were conducted. The Kb value was found to be 8.85 × 104 M-1 at a temperature of 25 °C by the use of fluorescence spectroscopic techniques. In summary, the utilization of electrochemical and spectroscopic techniques, alongside molecular docking investigations, has led to the prediction that desloratadine has the capability to interact with ct-dsDNA by groove binding.

SAnDReS 2.0: Development of machine-learning models to explore the scoring function space

Classical scoring functions may exhibit low accuracy in determining ligand binding affinity for proteins. The availability of both protein-ligand structures and affinity data make it possible to develop machine-learning models focused on specific protein systems with superior predictive performance. Here, we report a new methodology named SAnDReS that combines AutoDock Vina 1.2 with 54 regression methods available in Scikit-Learn to calculate binding affinity based on protein-ligand structures. This approach allows exploration of the scoring function space. SAnDReS generates machine-learning models based on crystal, docked, and AlphaFold-generated structures. As a proof of concept, we examine the performance of SAnDReS-generated models in three case studies. For all three cases, our models outperformed classical scoring functions. Also, SAnDReS-generated models showed predictive performance close to or better than other machine-learning models such as KDEEP, CSM-lig, and ΔVinaRF20. SAnDReS 2.0 is available to download at https://github.com/azevedolab/sandres.

HydraScreen: A Generalizable Structure-Based Deep Learning Approach to Drug Discovery

We propose HydraScreen, a deep-learning framework for safe and robust accelerated drug discovery. HydraScreen utilizes a state-of-the-art 3D convolutional neural network designed for the effective representation of molecular structures and interactions in protein-ligand binding. We designed an end-to-end pipeline for high-throughput screening and lead optimization, targeting applications in structure-based drug design. We assessed our approach using established public benchmarks based on the CASF-2016 core set, achieving top-tier results in affinity and pose prediction (Pearson's r = 0.86, RMSE = 1.15, Top-1 = 0.95). We introduced a novel approach for interaction profiling, aimed at detecting potential biases within both the model and data sets. This approach not only enhanced interpretability but also reinforced the impartiality of our methodology. Finally, we demonstrated HydraScreen's ability to generalize effectively across novel proteins and ligands through a temporal split. We also provide insights into potential avenues for future development aimed at enhancing the robustness of machine learning scoring functions. HydraScreen (accessible at http://hydrascreen.ro5.ai/paper) provides a user-friendly GUI and a public API, facilitating the easy-access assessment of protein-ligand complexes.

Do electrostatic interactions make a difference in physics-based AutoDock4 scoring function?

Physics-based scoring function AutoDock4 is one of the most successfully applied tools in the area of structure-based drug design. However, current scoring functions are still far from being perfect. In a recent work highlighting the strengths and deficiencies of current scoring functions, we discovered that the residual error of ΔGbind predictions made by AutoDock4 is highly correlated to the presence of formally charged fragments in a ligand. In this work, we study how the use of the high-quality atomic charges, applied for contemporary force fields calculation, affects the quality of the experimental ΔGbind prediction by means of AutoDock4. We initially expected that the previously found discrepancy could be attributed to the Gasteiger charges used within AutoDock4. We show that AutoDock4 is, surprisingly, not sensitive to the charges used, and the use of QC-derived atomic charges does not lead to any statistical improvements. We also briefly discuss the role of the explicit empirical hydrogen bond term along with the electrostatic term.

Geometry Optimization Algorithms in Conjunction with the Machine Learning Potential ANI-2x Facilitate the Structure-Based Virtual Screening and Binding Mode Prediction

Structure-based virtual screening utilizes molecular docking to explore and analyze ligand-macromolecule interactions, crucial for identifying and developing potential drug candidates. Although there is availability of several widely used docking programs, the accurate prediction of binding affinity and binding mode still presents challenges. In this study, we introduced a novel protocol that combines our in-house geometry optimization algorithm, the conjugate gradient with backtracking line search (CG-BS), which is capable of restraining and constraining rotatable torsional angles and other geometric parameters with a highly accurate machine learning potential, ANI-2x, renowned for its precise molecular energy predictions reassembling the wB97X/6-31G(d) model. By integrating this protocol with binding pose prediction using the Glide, we conducted additional structural optimization and potential energy prediction on 11 small molecule-macromolecule and 12 peptide-macromolecule systems. We observed that ANI-2x/CG-BS greatly improved the docking power, not only optimizing binding poses more effectively, particularly when the RMSD of the predicted binding pose by Glide exceeded around 5 Å, but also achieving a 26% higher success rate in identifying those native-like binding poses at the top rank compared to Glide docking. As for the scoring and ranking powers, ANI-2x/CG-BS demonstrated an enhanced performance in predicting and ranking hundreds or thousands of ligands over Glide docking. For example, Pearson's and Spearman's correlation coefficients remarkedly increased from 0.24 and 0.14 with Glide docking to 0.85 and 0.69, respectively, with the addition of ANI-2x/CG-BS for optimizing and ranking small molecules binding to the bacterial ribosomal aminoacyl-tRNA receptor. These results suggest that ANI-2x/CG-BS holds considerable potential for being integrated into virtual screening pipelines due to its enhanced docking performance.

Rational Design of Drugs Targeting G-Protein-Coupled Receptors: Ligand Search and Screening

G protein-coupled receptors (GPCRs) are transmembrane proteins that participate in many physiological processes and represent major pharmacological targets. Recent advances in structural biology of GPCRs have enabled the development of drugs based on the receptor structure (structure-based drug design, SBDD). SBDD utilizes information about the receptor-ligand complex to search for suitable compounds, thus expanding the chemical space of possible receptor ligands without the need for experimental screening. The review describes the use of structure-based virtual screening (SBVS) for GPCR ligands and approaches for the functional testing of potential drug compounds, as well as discusses recent advances and successful examples in the application of SBDD for the identification of GPCR ligands.

Water Network-Augmented Two-State Model for Protein-Ligand Binding Affinity Prediction

Water network rearrangement from the ligand-unbound state to the ligand-bound state is known to have significant effects on the protein-ligand binding interactions, but most of the current machine learning-based scoring functions overlook these effects. In this study, we endeavor to construct a comprehensive and realistic deep learning model by incorporating water network information into both ligand-unbound and -bound states. In particular, extended connectivity interaction features were integrated into graph representation, and graph transformer operator was employed to extract features of the ligand-unbound and -bound states. Through these efforts, we developed a water network-augmented two-state model called ECIFGraph::HM-Holo-Apo. Our new model exhibits satisfactory performance in terms of scoring, ranking, docking, screening, and reverse screening power tests on the CASF-2016 benchmark. In addition, it can achieve superior performance in large-scale docking-based virtual screening tests on the DEKOIS2.0 data set. Our study highlights that the use of a water network-augmented two-state model can be an effective strategy to bolster the robustness and applicability of machine learning-based scoring functions, particularly for targets with hydrophilic or solvent-exposed binding pockets.

An ensemble-based approach to estimate confidence of predicted protein-ligand binding affinity values

When designing a machine learning-based scoring function, we access a limited number of protein-ligand complexes with experimentally determined binding affinity values, representing only a fraction of all possible protein-ligand complexes. Consequently, it is crucial to report a measure of confidence and quantify the uncertainty in the model's predictions during test time. Here, we adopt the conformal prediction technique to evaluate the confidence of a prediction for each member of the core set of the CASF 2016 benchmark. The conformal prediction technique requires a diverse ensemble of predictors for uncertainty estimation. To this end, we introduce ENS-Score as an ensemble predictor, which includes 30 models with different protein-ligand representation approaches and achieves Pearson's correlation of 0.842 on the core set of the CASF 2016 benchmark. Also, we comprehensively investigate the residual error of each data point to assess the normality behavior of the distribution of the residual errors and their correlation to the structural features of the ligands, such as hydrophobic interactions and halogen bonding. In the end, we provide a local host web application to facilitate the usage of ENS-Score. All codes to repeat results are provided at https://github.com/miladrayka/ENS_Score.

Enhancing Generalizability in Protein-Ligand Binding Affinity Prediction with Multimodal Contrastive Learning

Improving the generalization ability of scoring functions remains a major challenge in protein-ligand binding affinity prediction. Many machine learning methods are limited by their reliance on single-modal representations, hindering a comprehensive understanding of protein-ligand interactions. We introduce a graph-neural-network-based scoring function that utilizes a triplet contrastive learning loss to improve protein-ligand representations. In this model, three-dimensional complex representations and the fusion of two-dimensional ligand and coarse-grained pocket representations converge while distancing from decoy representations in latent space. After rigorous validation on multiple external data sets, our model exhibits commendable generalization capabilities compared to those of other deep learning-based scoring functions, marking it as a promising tool in the realm of drug discovery. In the future, our training framework can be extended to other biophysical- and biochemical-related problems such as protein-protein interaction and protein mutation prediction.

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.

Prediction of protein-ligand binding affinity via deep learning models

Accurately predicting the binding affinity between proteins and ligands is crucial in drug screening and optimization, but it is still a challenge in computer-aided drug design. The recent success of AlphaFold2 in predicting protein structures has brought new hope for deep learning (DL) models to accurately predict protein-ligand binding affinity. However, the current DL models still face limitations due to the low-quality database, inaccurate input representation and inappropriate model architecture. In this work, we review the computational methods, specifically DL-based models, used to predict protein-ligand binding affinity. We start with a brief introduction to protein-ligand binding affinity and the traditional computational methods used to calculate them. We then introduce the basic principles of DL models for predicting protein-ligand binding affinity. Next, we review the commonly used databases, input representations and DL models in this field. Finally, we discuss the potential challenges and future work in accurately predicting protein-ligand binding affinity via DL models.

Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Software and Databases for Protein-Protein Docking

Protein-protein interactions (PPIs) provide valuable insights for understanding the principles of biological systems and for elucidating causes of incurable diseases. One of the techniques used for computational prediction of PPIs is protein-protein docking calculations, and a variety of software has been developed. This chapter is a summary of software and databases used for protein-protein docking.

Tetraenone A: A New β-Ionone Derivative from Tetraena aegyptia

In this study, the chemical investigation of Tetraena aegyptia (Zygophyllaceae) led to the identification of a new megastigmene derivative, tetraenone A ((2S, 5R, 6R, 7E)-2-hydroxy-5,6-dihydro-β-ionone) (1), along with (3S, 5R, 6S, 7E)-3-hydroxy-5,6-epoxy-5,6-dihydro-β-ionone- (2), 3,4-dihydroxy-cinnamyl alcohol-4-glucoside (3), 3β,19α-dihydroxy-ursan-28-oic acid (4), quinovic acid (5), p-coumaric acid (6), and ferulic acid (7), for the first time. The chemical structures of 1-7 were confirmed by analysis of their 1D and 2D NMR and HRESIMS spectra and by their comparison with the relevant literature. The absolute configurations of 1 and 2 were assigned based on NOESY interactions and ECD spectra. Conformational analysis showed that 1 existed exclusively in one of the two theoretically possible chair conformers with a predominant s-trans configuration for the 3-oxobut-1-en-1-yl group with the ring, while the half-chair conformer had a pseudo-axial hydroxy group that was predominant over the other half-chair conformation. Boat conformations were not among the most stable conformations, and the s-trans isomerism was in favor of s-cis configuration. In silico investigation revealed that 1 and 2 had more favorable binding interactions with Mpro rather than with TMPRSS2. Accordingly, molecular dynamic simulations were performed on the complexes of compounds 1 and 2 with Mpro to explore the stability of their interaction with the target protein structure. Compounds 1 and 2 might offer a possible starting point for developing covalent inhibitors of Mpro of SARS-CoV-2.

Ensemble of local and global information for Protein-Ligand Binding Affinity Prediction

Accurately predicting protein-ligand binding affinities is crucial for determining molecular properties and understanding their physical effects. Neural networks and transformers are the predominant methods for sequence modeling, and both have been successfully applied independently for protein-ligand binding affinity prediction. As local and global information of molecules are vital for protein-ligand binding affinity prediction, we aim to combine bi-directional gated recurrent unit (BiGRU) and convolutional neural network (CNN) to effectively capture both local and global molecular information. Additionally, attention mechanisms can be incorporated to automatically learn and adjust the level of attention given to local and global information, thereby enhancing the performance of the model. To achieve this, we propose the PLAsformer approach, which encodes local and global information of molecules using 3DCNN and BiGRU with attention mechanism, respectively. This approach enhances the model's ability to encode comprehensive local and global molecular information. PLAsformer achieved a Pearson's correlation coefficient of 0.812 and a Root Mean Square Error (RMSE) of 1.284 when comparing experimental and predicted affinity on the PDBBind-2016 dataset. These results surpass the current state-of-the-art methods for binding affinity prediction. The high accuracy of PLAsformer's predictive performance, along with its excellent generalization ability, is clearly demonstrated by these findings.

An overview of recent advances and challenges in predicting compound-protein interaction (CPI)

Compound-protein interactions (CPIs) are critical in drug discovery for identifying therapeutic targets, drug side effects, and repurposing existing drugs. Machine learning (ML) algorithms have emerged as powerful tools for CPI prediction, offering notable advantages in cost-effectiveness and efficiency. This review provides an overview of recent advances in both structure-based and non-structure-based CPI prediction ML models, highlighting their performance and achievements. It also offers insights into CPI prediction-related datasets and evaluation benchmarks. Lastly, the article presents a comprehensive assessment of the current landscape of CPI prediction, elucidating the challenges faced and outlining emerging trends to advance the field.

Anatomy of Potency Predictions Focusing on Structural Analogues with Increasing Potency Differences Including Activity Cliffs

Potency predictions are popular in compound design and optimization but are complicated by intrinsic limitations. Moreover, even for nonlinear methods, activity cliffs (ACs, formed by structural analogues with large potency differences) represent challenging test cases for compound potency predictions. We have devised a new test system for potency predictions, including AC compounds, that is based on partitioned matched molecular pairs (MMP) and makes it possible to monitor prediction accuracy at the level of analogue pairs with increasing potency differences. The results of systematic predictions using different machine learning and control methods on MMP-based data sets revealed increasing prediction errors when potency differences between corresponding training and test compounds increased, including large prediction errors for AC compounds. At the global level, these prediction errors were not apparent due to the statistical dominance of analogue pairs with small potency differences. Test compounds from such pairs were accurately predicted and determined the observed global prediction accuracy. Shapley value analysis, an explainable artificial intelligence approach, was applied to identify structural features determining potency predictions using different methods. The analysis revealed that numerical predictions of different regression models were determined by features that were shared by MMP partner compounds or absent in these compounds, with opposing effects. These findings provided another rationale for accurate predictions of similar potency values for structural analogues and failures in predicting the potency of AC compounds.

Applying molecular docking to pesticides

Pesticide creation is related to the development of sustainable agricultural and ecological safety, and molecular docking technology can effectively help in pesticide innovation. This paper introduces the basic theory behind molecular docking, pesticide databases, and docking software. It also summarizes the application of molecular docking in the pesticide field, including the virtual screening of lead compounds, detection of pesticides and their metabolites in the environment, reverse screening of pesticide targets, and the study of resistance mechanisms. Finally, problems with the use of molecular docking technology in pesticide creation are discussed, and prospects for the future use of molecular docking technology in new pesticide development are discussed. © 2023 Society of Chemical Industry.

Multi-shelled ECIF: improved extended connectivity interaction features for accurate binding affinity prediction

Motivation

Extended connectivity interaction features (ECIF) is a method developed to predict protein-ligand binding affinity, allowing for detailed atomic representation. It performed very well in terms of Comparative Assessment of Scoring Functions 2016 (CASF-2016) scoring power. However, ECIF has the limitation of not being able to adequately account for interatomic distances.

Results

To investigate what kind of distance representation is effective for P-L binding affinity prediction, we have developed two algorithms that improved ECIF's feature extraction method to take distance into account. One is multi-shelled ECIF, which takes into account the distance between atoms by dividing the distance between atoms into multiple layers. The other is weighted ECIF, which weights the importance of interactions according to the distance between atoms. A comparison of these two methods shows that multi-shelled ECIF outperforms weighted ECIF and the original ECIF, achieving a CASF-2016 scoring power Pearson correlation coefficient of 0.877.

Availability and implementation

All the codes and data are available on GitHub (https://github.com/koji11235/MSECIFv2).

Molecular Docking and Molecular Dynamics Studies Reveal the Anticancer Potential of Medicinal-Plant-Derived Lignans as MDM2-P53 Interaction Inhibitors

The interaction between the tumor suppressor protein p53 and its negative regulator, the MDM2 oncogenic protein, has gained significant attention in cancer drug discovery. In this study, 120 lignans reported from Ferula sinkiangensis and Justicia procumbens were assessed for docking simulations on the active pocket of the MDM2 crystal structure bound to Nutlin-3a. The docking analysis identified nine compounds with higher docking scores than the co-crystallized reference. Subsequent AMDET profiling revealed satisfactory pharmacokinetic and safety parameters for these natural products. Three compounds, namely, justin A, 6-hydroxy justicidin A, and 6'-hydroxy justicidin B, were selected for further investigation due to their strong binding affinities of -7.526 kcal/mol, -7.438 kcal/mol, and -7.240 kcal/mol, respectively, which surpassed the binding affinity of the reference inhibitor Nutlin-3a (-6.830 kcal/mol). To assess the stability and reliability of the binding of the candidate hits, a molecular dynamics simulation was performed over a duration of 100 ns. Remarkably, the thorough analysis demonstrated that all the hits exhibited stable molecular dynamics profiles. Based on their effective binding to MDM2, favorable pharmacokinetic properties, and molecular dynamics behavior, these compounds represent a promising starting point for further refinement. Nevertheless, it is essential to synthesize the suggested compounds and evaluate their activity through in vitro and in vivo experiments.

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.

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

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

Protein-ligand binding affinity prediction exploiting sequence constituent homology

MOTIVATION

Molecular docking is a commonly used approach for estimating binding conformations and their resultant binding affinities. Machine learning has been successfully deployed to enhance such affinity estimations. Many methods of varying complexity have been developed making use of some or all the spatial and categorical information available in these structures. The evaluation of such methods has mainly been carried out using datasets from PDBbind. Particularly the Comparative Assessment of Scoring Functions (CASF) 2007, 2013, and 2016 datasets with dedicated test sets. This work demonstrates that only a small number of simple descriptors is necessary to efficiently estimate binding affinity for these complexes without the need to know the exact binding conformation of a ligand.

RESULTS

The developed approach of using a small number of ligand and protein descriptors in conjunction with gradient boosting trees demonstrates high performance on the CASF datasets. This includes the commonly used benchmark CASF2016 where it appears to perform better than any other approach. This methodology is also useful for datasets where the spatial relationship between the ligand and protein is unknown as demonstrated using a large ChEMBL-derived dataset.

AVAILABILITY AND IMPLEMENTATION

Code and data uploaded to https://github.com/abbiAR/PLBAffinity.

Solvated interaction energy: from small-molecule to antibody drug design

Scoring functions are ubiquitous in structure-based drug design as an aid to predicting binding modes and estimating binding affinities. Ideally, a scoring function should be broadly applicable, obviating the need to recalibrate and refit its parameters for every new target and class of ligands. Traditionally, drugs have been small molecules, but in recent years biologics, particularly antibodies, have become an increasingly important if not dominant class of therapeutics. This makes the goal of having a transferable scoring function, i.e., one that spans the range of small-molecule to protein ligands, even more challenging. One such broadly applicable scoring function is the Solvated Interaction Energy (SIE), which has been developed and applied in our lab for the last 15 years, leading to several important applications. This physics-based method arose from efforts to understand the physics governing binding events, with particular care given to the role played by solvation. SIE has been used by us and many independent labs worldwide for virtual screening and discovery of novel small-molecule binders or optimization of known drugs. Moreover, without any retraining, it is found to be transferrable to predictions of antibody-antigen relative binding affinities and as accurate as functions trained on protein-protein binding affinities. SIE has been incorporated in conjunction with other scoring functions into ADAPT (Assisted Design of Antibody and Protein Therapeutics), our platform for affinity modulation of antibodies. Application of ADAPT resulted in the optimization of several antibodies with 10-to-100-fold improvements in binding affinity. Further applications included broadening the specificity of a single-domain antibody to be cross-reactive with virus variants of both SARS-CoV-1 and SARS-CoV-2, and the design of safer antibodies by engineering of a pH switch to make them more selective towards acidic tumors while sparing normal tissues at physiological pH.

Molecular docking in organic, inorganic, and hybrid systems: a tutorial review

Molecular docking simulation is a very popular and well-established computational approach and has been extensively used to understand molecular interactions between a natural organic molecule (ideally taken as a receptor) such as an enzyme, protein, DNA, RNA and a natural or synthetic organic/inorganic molecule (considered as a ligand). But the implementation of docking ideas to synthetic organic, inorganic, or hybrid systems is very limited with respect to their use as a receptor despite their huge popularity in different experimental systems. In this context, molecular docking can be an efficient computational tool for understanding the role of intermolecular interactions in hybrid systems that can help in designing materials on mesoscale for different applications. The current review focuses on the implementation of the docking method in organic, inorganic, and hybrid systems along with examples from different case studies. We describe different resources, including databases and tools required in the docking study and applications. The concept of docking techniques, types of docking models, and the role of different intermolecular interactions involved in the docking process to understand the binding mechanisms are explained. Finally, the challenges and limitations of dockings are also discussed in this review.

Graphical abstract

Recent Advances in Computer-Aided Structure-Based Drug Design on Ion Channels

Ion channels play important roles in fundamental biological processes, such as electric signaling in cells, muscle contraction, hormone secretion, and regulation of the immune response. Targeting ion channels with drugs represents a treatment option for neurological and cardiovascular diseases, muscular degradation disorders, and pathologies related to disturbed pain sensation. While there are more than 300 different ion channels in the human organism, drugs have been developed only for some of them and currently available drugs lack selectivity. Computational approaches are an indispensable tool for drug discovery and can speed up, especially, the early development stages of lead identification and optimization. The number of molecular structures of ion channels has considerably increased over the last ten years, providing new opportunities for structure-based drug development. This review summarizes important knowledge about ion channel classification, structure, mechanisms, and pathology with the main focus on recent developments in the field of computer-aided, structure-based drug design on ion channels. We highlight studies that link structural data with modeling and chemoinformatic approaches for the identification and characterization of new molecules targeting ion channels. These approaches hold great potential to advance research on ion channel drugs in the future.

Designing highly potent compounds using a chemical language model

Compound potency prediction is a major task in medicinal chemistry and drug design. Inspired by the concept of activity cliffs (which encode large differences in potency between similar active compounds), we have devised a new methodology for predicting potent compounds from weakly potent input molecules. Therefore, a chemical language model was implemented consisting of a conditional transformer architecture for compound design guided by observed potency differences. The model was evaluated using a newly generated compound test system enabling a rigorous assessment of its performance. It was shown to predict known potent compounds from different activity classes not encountered during training. Moreover, the model was capable of creating highly potent compounds that were structurally distinct from input molecules. It also produced many novel candidate compounds not included in test sets. Taken together, the findings confirmed the ability of the new methodology to generate structurally diverse highly potent compounds.

Bifenthrin disrupts cytochrome c oxidase activity and reduces mitochondrial DNA copy number through oxidative damage in pool barb (Puntius sophore)

Bifenthrin (BF), a synthetic pyrethroid is used worldwide for both agricultural and non-agricultural purposes due to its high insecticidal activity and low toxicity in mammals. However, its improper usage implies a possible risk to aquatic life. The Study was aimed to correlate the association of BF toxicity with mitochondrial DNA copy number variation in edible fish Punitus sophore. The 96-h LC 50 of BF in P. sophore was 3.4 μg/L, fish was treated with sub-lethal doses (0.34 μg/L,0.68 μg/L) of BF for 15 days. The activity and expression level of cytochrome c oxidase (Mt-COI) were measured to assess mitochondrial dysfunction caused by BF. Results showed BF reduced the level of Mt-COI mRNA in treated groups, hindered complex IV activity and increased ROS generation leading to oxidative damage. mtDNAcn was decreased in the muscle, brain and liver after BF treatment. Furthermore, BF induced neurotoxicity in brain and muscle cells through the inhibition of AchE activity. The treated groups showed elevated level of malondialdehyde (MDA) and an imbalance of antioxidant enzymes activity. Molecular docking and simulation analysis also predicted that BF binds to the active sites of the enzyme and restricts the fluctuation of active sites' residues. Hence, outcome of the study suggests reduction of mtDNAcn could be a potential biomarker to assess Bifenthrin induced toxicity in aquatic ecosystem.

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Identification of Catechins' Binding Sites in Monomeric Aβ42 through Ensemble Docking and MD Simulations

The assembly of the amyloid-β peptide (Aβ) into toxic oligomers and fibrils is associated with Alzheimer's disease and dementia. Therefore, disrupting amyloid assembly by direct targeting of the Aβ monomeric form with small molecules or antibodies is a promising therapeutic strategy. However, given the dynamic nature of Aβ, standard computational tools cannot be easily applied for high-throughput structure-based virtual screening in drug discovery projects. In the current study, we propose a computational pipeline-in the framework of the ensemble docking strategy-to identify catechins' binding sites in monomeric Aβ₄₂. It is shown that both hydrophobic aromatic interactions and hydrogen bonding are crucial for the binding of catechins to Aβ₄₂. Additionally, it has been found that all the studied ligands, especially EGCG, can act as potent inhibitors against amyloid aggregation by blocking the central hydrophobic region of Aβ. Our findings are evaluated and confirmed with multi-microsecond MD simulations. Finally, it is suggested that our proposed pipeline, with low computational cost in comparison with MD simulations, is a suitable approach for the virtual screening of ligand libraries against Aβ.

Exploring the therapeutic mechanisms of Gleditsiae Spina acting on pancreatic cancer via network pharmacology, molecular docking and molecular dynamics simulation

Pancreatic cancer is one of the most aggressive tumors and also has a low survival rate. The dried spines of Gleditsia sinensis Lam are known as "Gleditsiae Spina" and they mostly contain flavonoids, phenolic acids, terpenoids, steroids, and other chemical components. In this study, the potential active components and molecular mechanisms of Gleditsiae Spina for treating pancreatic cancer were systematically revealed by network pharmacology, molecular docking and molecular dynamics simulations (MDs). RAC-alpha serine/threonine-protein kinase (AKT1), cellular tumor antigen p53 (TP53), tumor necrosis factor α (TNFα), interleukin-6 (IL6) and vascular endothelial growth factor A (VEGFA) were common targets of Gleditsiae Spina, human cytomegalovirus infection signaling pathway, AGE-RAGE signaling pathway in diabetic complications, and MAPK signaling pathway were critical pathways of fisetin, eriodyctiol, kaempferol and quercetin in the treatment of pancreatic cancer. Molecular dynamics simulations (MDs) results showed that eriodyctiol and kaempferol have long-term stable hydrogen bonds and high binding free energy for TP53 (-23.64 ± 0.03 kcal mol-1 and -30.54 ± 0.02 kcal mol-1, respectively). Collectively, our findings identify active components and potential targets in Gleditsiae Spina for the treatment of pancreatic cancer, which may help to explore leading compounds and potential drugs for pancreatic cancer.

Design of Novel Letrozole Analogues Targeting Aromatase for Breast Cancer: Molecular Docking, Molecular Dynamics, and Theoretical Studies on Gold Nanoparticles

The use of aromatase inhibitors is an established therapy for estrogen-dependent breast cancer in postmenopausal women. However, the only commercially available aromatase inhibitor, letrozole, is not highly selective; in addition to aromatase, it has an affinity for binding to desmolase, an enzyme involved in steroidogenesis, which explains the main side effects. Therefore, we designed new compounds based on the structure of letrozole. More than five thousand compounds were constructed based on the letrozole structure. Then, these compounds were screened for their binding ability toward the target protein, aromatase. Quantum docking, Glide docking, and ADME studies showed 14 new molecules with docking scores of ≤-7 kcal/mol, compared to the docking score of -4.109 kcal/mol of the reference, letrozole. Moreover, molecular dynamics (MD) and post-MD MM-GBSA calculations were calculated for the top three compounds, and the results supported in their interaction's stability. Finally, the density-functional theory (DFT) study applied to the top compound to study the interaction with gold nanoparticles revealed the most stable position for the interaction with the gold nanoparticles. The results of this study confirmed that these newly designed compounds could be useful starting points for lead optimization. Further in vitro and in vivo studies are recommended for these compounds to verify these promising results experimentally.

An artificial neural network model to predict structure-based protein-protein free energy of binding from Rosetta-calculated properties

The prediction of the free energy (ΔG) of binding for protein-protein complexes is of general scientific interest as it has a variety of applications in the fields of molecular and chemical biology, materials science, and biotechnology. Despite its centrality in understanding protein association phenomena and protein engineering, the ΔG of binding is a daunting quantity to obtain theoretically. In this work, we devise a novel Artificial Neural Network (ANN) model to predict the ΔG of binding for a given three-dimensional structure of a protein-protein complex with Rosetta-calculated properties. Our model was tested using two data sets, and it presented a root-mean-square error ranging from 1.67 kcal mol^-1 to 2.45 kcal mol^-1, showing a better performance compared to the available state-of-the-art tools. Validation of the model for a variety of protein-protein complexes is showcased.

GB-score: Minimally designed machine learning scoring function based on distance-weighted interatomic contact features

In recent years, thanks to advances in computer hardware and dataset availability, data-driven approaches (like machine learning) have become one of the essential parts of the drug design framework to accelerate drug discovery procedures. Constructing a new scoring function, a function that can predict the binding score for a generated protein-ligand pose during docking procedure or a crystal complex, based on machine and deep learning has become an active research area in computer-aided drug design. GB-Score is a state-of-the-art machine learning-based scoring function that utilizes distance-weighted interatomic contact features, PDBbind-v2019 general set, and Gradient Boosting Trees algorithm to the binding affinity prediction. The distance-weighted interatomic contact featurization method used the distance between different ligand and protein atom types for numerical representation of the protein-ligand complex. GB-Score attains Pearson's correlation 0.862 and RMSE 1.190 on the CASF-2016 benchmark test in the scoring power metric. GB-Score's codes are freely available on the web at https://github.com/miladrayka/GB_Score.

Application of computational methods for class A GPCR Ligand discovery

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

Predicting Potent Compounds Using a Conditional Variational Autoencoder Based upon a New Structure-Potency Fingerprint

Prediction of the potency of bioactive compounds generally relies on linear or nonlinear quantitative structure-activity relationship (QSAR) models. Nonlinear models are generated using machine learning methods. We introduce a novel approach for potency prediction that depends on a newly designed molecular fingerprint (FP) representation. This structure-potency fingerprint (SPFP) combines different modules accounting for the structural features of active compounds and their potency values in a single bit string, hence unifying structure and potency representation. This encoding enables the derivation of a conditional variational autoencoder (CVAE) using SPFPs of training compounds and apply the model to predict the SPFP potency module of test compounds using only their structure module as input. The SPFP-CVAE approach correctly predicts the potency values of compounds belonging to different activity classes with an accuracy comparable to support vector regression (SVR), representing the state-of-the-art in the field. In addition, highly potent compounds are predicted with very similar accuracy as SVR and deep neural networks.

Persistent Path-Spectral (PPS) Based Machine Learning for Protein-Ligand Binding Affinity Prediction

Molecular descriptors are essential to quantitative structure activity/property relationship (QSAR/QSPR) models and machine learning models. Here we propose persistent path-spectral (PPS), PPS-based molecular descriptors, and PPS-based machine learning model for the prediction of the protein-ligand binding affinity, for the first time. For the graph, simplicial complex, and hypergraph representation of molecular structures and interactions, the path-Laplacian can be constructed and the derived path-spectral naturally gives a quantitative description of molecules. Further, by introducing the filtration process of the representation, the persistent path-spectral can be derived, which gives a multiscale characterization of molecules. Molecular descriptors from the persistent path-spectral attributes then are combined with the machine learning model, in particular, the gradient boosting tree, to form our PPS-ML model. We test our model on three most commonly used data sets, i.e., PDBbind-v2007, PDBbind-v2013, and PDBbind-v2016, and our model can achieve competitive results.

Machine Learning Scoring Functions for Drug Discovery from Experimental and Computer-Generated Protein-Ligand Structures: Towards Per-Target Scoring Functions

In recent years, machine learning has been proposed as a promising strategy to build accurate scoring functions for computational docking finalized to numerically empowered drug discovery. However, the latest studies have suggested that over-optimistic results had been reported due to the correlations present in the experimental databases used for training and testing. Here, we investigate the performance of an artificial neural network in binding affinity predictions, comparing results obtained using both experimental protein-ligand structures as well as larger sets of computer-generated structures created using commercial software. Interestingly, similar performances are obtained on both databases. We find a noticeable performance suppression when moving from random horizontal tests to vertical tests performed on target proteins not included in the training data. The possibility to train the network on relatively easily created computer-generated databases leads us to explore per-target scoring functions, trained and tested ad-hoc on complexes including only one target protein. Encouraging results are obtained, depending on the type of protein being addressed.

Docking is not enough: 17-trifluoromethylphenyl trinor PGF2α is only a very weak ligand of neurokinin-1 receptor

17-trifluoromethylphenyl trinor prostaglandin F_2α (17-CF₃PTPGF_2α) was reported recently to exhibit in vitro and in vivo anticancer activity. Based solely on the results of in silico molecular docking, it was claimed that this compound is NK1 receptor (NK1R) antagonist and that its activity is through this receptor. In this contribution we show that 17-CF₃PTPGF_2α is only a very weak NK1R ligand (IC₅₀ > 200 μM). In connection with that we discuss the issue of this compound's molecular target. Finally, we briefly narrate on the proper use of molecular docking in biomedical research.

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.

Molecular Recognition Patterns between Vitamin B12 and Proteins Explored through STD-NMR and In Silico Studies

Ligand-receptor molecular recognition is the basis of biological processes. The Saturation Transfer Difference-NMR (STD-NMR) technique has been recently used to gain qualitative and quantitative information about physiological interactions at an atomic resolution. The molecular recognition patterns between the cyanocobalamin (CNBL)/aqua cobalamin (OHBL) and different plant and animal proteins were investigated via STD-NMR supplemented by molecular docking. This study demonstrates that myoglobin has the highest binding affinity and that gluten has the lowest affinity. Casein also shows a higher binding affinity for cyanocobalamin when compared with that of plant-based proteins. STD-NMR results showed the moderate binding capability of casein with both CNBL and OHBL. Computer simulation confirmed the recognition mode in theory and was compared with the experiments. This work is beneficial for understanding the binding affinity and biological action of cyanocobalamin and will attract researchers to use NMR technology to link the chemical and physiological properties of nutrients.

Hierarchical analysis of the target-based scoring function modification for the example of selected class A GPCRs

Computational methods, especially molecular docking-based calculations, have become indispensable in the modern drug discovery workflow. The constantly increasing chemical space requires fast, robust but most of all highly predictive methods to search for new bioactive agents. Thus, the scoring function (SF) is a useful and broadly applied energy-based element of docking software, allowing quick and effective evaluation of a ligand's propensity to bind to selected protein targets. Despite many spectacular successes of molecular docking applications in virtual screening (VS), the obtained results are often far from ideal, leading to incorrect selection of hit molecules and poor pose prediction. In our study we focused on docking calculation for the selected class A G-protein coupled receptors (GPCRs), with experimentally determined 3D structures and a sufficient set of known ligands with affinity values reported in the ChEMBL database. Our goal is to investigate how much the energy-based scoring function for this particular target class changes when changing from the default to the re-estimated weighting scheme on the specified energy terms in the SF definition. Additionally, we want to verify if indeed more accurate results are obtained when considering different levels of the biological hierarchy, namely: the whole class A GPCRs, sub-subfamilies, or just the individual proteins while applying default or specifically designed weighting coefficients. The performed calculation and evaluation factor values suggest a significant improvement of docking results for the designed SF definition. This individual approach improves the accuracy of binding affinity prediction and active compound recognition. The designed scoring function for classes, sub-subfamilies, or proteins leads to a significant improvement of molecular docking performance, especially at the level of individual proteins. Our results show that to increase the efficiency and predictive power of molecular docking calculations applied in classical VS, the strategy based on the individual approach for scoring function definition for selected proteins should be considered.

Leveraging scaffold information to predict protein-ligand binding affinity with an empirical graph neural network

Protein-ligand binding affinity prediction is an important task in structural bioinformatics for drug discovery and design. Although various scoring functions (SFs) have been proposed, it remains challenging to accurately evaluate the binding affinity of a protein-ligand complex with the known bound structure because of the potential preference of scoring system. In recent years, deep learning (DL) techniques have been applied to SFs without sophisticated feature engineering. Nevertheless, existing methods cannot model the differential contribution of atoms in various regions of proteins, and the relationship between atom properties and intermolecular distance is also not fully explored. We propose a novel empirical graph neural network for accurate protein-ligand binding affinity prediction (EGNA). Graphs of protein, ligand and their interactions are constructed based on different regions of each bound complex. Proteins and ligands are effectively represented by graph convolutional layers, enabling the EGNA to capture interaction patterns precisely by simulating empirical SFs. The contributions of different factors on binding affinity can thus be transparently investigated. EGNA is compared with the state-of-the-art machine learning-based SFs on two widely used benchmark data sets. The results demonstrate the superiority of EGNA and its good generalization capability.

Comprehensive Survey of Consensus Docking for High-Throughput Virtual Screening

The rapid advances of 3D techniques for the structural determination of proteins and the development of numerous computational methods and strategies have led to identifying highly active compounds in computer drug design. Molecular docking is a method widely used in high-throughput virtual screening campaigns to filter potential ligands targeted to proteins. A great variety of docking programs are currently available, which differ in the algorithms and approaches used to predict the binding mode and the affinity of the ligand. All programs heavily rely on scoring functions to accurately predict ligand binding affinity, and despite differences in performance, none of these docking programs is preferable to the others. To overcome this problem, consensus scoring methods improve the outcome of virtual screening by averaging the rank or score of individual molecules obtained from different docking programs. The successful application of consensus docking in high-throughput virtual screening highlights the need to optimize the predictive power of molecular docking methods.

In Vitro and In Silico Studies for the Identification of Potent Metabolites of Some High-Altitude Medicinal Plants from Nepal Inhibiting SARS-CoV-2 Spike Protein

Despite ongoing vaccination programs against COVID-19 around the world, cases of infection are still rising with new variants. This infers that an effective antiviral drug against COVID-19 is crucial along with vaccinations to decrease cases. A potential target of such antivirals could be the membrane components of the causative pathogen, SARS-CoV-2, for instance spike (S) protein. In our research, we have deployed in vitro screening of crude extracts of seven ethnomedicinal plants against the spike receptor-binding domain (S1-RBD) of SARS-CoV-2 using an enzyme-linked immunosorbent assay (ELISA). Following encouraging in vitro results for Tinospora cordifolia, in silico studies were conducted for the 14 reported antiviral secondary metabolites isolated from T. cordifolia-a species widely cultivated and used as an antiviral drug in the Himalayan country of Nepal-using Genetic Optimization for Ligand Docking (GOLD), Molecular Operating Environment (MOE), and BIOVIA Discovery Studio. The molecular docking and binding energy study revealed that cordifolioside-A had a higher binding affinity and was the most effective in binding to the competitive site of the spike protein. Molecular dynamics (MD) simulation studies using GROMACS 5.4.1 further assayed the interaction between the potent compound and binding sites of the spike protein. It revealed that cordifolioside-A demonstrated better binding affinity and stability, and resulted in a conformational change in S1-RBD, hence hindering the activities of the protein. In addition, ADMET analysis of the secondary metabolites from T. cordifolia revealed promising pharmacokinetic properties. Our study thus recommends that certain secondary metabolites of T. cordifolia are possible medicinal candidates against SARS-CoV-2.

Systematic Investigation of Docking Failures in Large-Scale Structure-Based Virtual Screening

In recent years, large-scale structure-based virtual screening has attracted increasing levels of interest for identification of novel compounds corresponding to potential drug targets. It is critical to understand the strengths and weaknesses of docking algorithms to increase the success rate in practical applications. Here, we systematically investigated the docking successes and failures of two representative docking programs: UCSF DOCK 3.7 and AutoDock Vina. DOCK 3.7 performed better in early enrichment on the Directory of Useful Decoys: Enhanced (DUD-E) data set, although both docking methods were roughly comparable in overall enrichment performance. DOCK 3.7 also showed superior computational efficiency. Intriguingly, the Vina scoring function showed a bias toward compounds with higher molecular weights. Both the tested docking approaches yielded incorrectly predicted ligand binding poses caused by the limitations of torsion sampling. Based on a careful analysis of docking results from six representative cases, we propose the reasons underlying docking failures; furthermore, we provide a few solutions, representing practical guidance for large-scale virtual screening campaigns and future docking algorithm development.

Biomolecular Topology: Modelling and Analysis

With the great advancement of experimental tools, a tremendous amount of biomolecular data has been generated and accumulated in various databases. The high dimensionality, structural complexity, the nonlinearity, and entanglements of biomolecular data, ranging from DNA knots, RNA secondary structures, protein folding configurations, chromosomes, DNA origami, molecular assembly, to others at the macromolecular level, pose a severe challenge in their analysis and characterization. In the past few decades, mathematical concepts, models, algorithms, and tools from algebraic topology, combinatorial topology, computational topology, and topological data analysis, have demonstrated great power and begun to play an essential role in tackling the biomolecular data challenge. In this work, we introduce biomolecular topology, which concerns the topological problems and models originated from the biomolecular systems. More specifically, the biomolecular topology encompasses topological structures, properties and relations that are emerged from biomolecular structures, dynamics, interactions, and functions. We discuss the various types of biomolecular topology from structures (of proteins, DNAs, and RNAs), protein folding, and protein assembly. A brief discussion of databanks (and databases), theoretical models, and computational algorithms, is presented. Further, we systematically review related topological models, including graphs, simplicial complexes, persistent homology, persistent Laplacians, de Rham-Hodge theory, Yau-Hausdorff distance, and the topology-based machine learning models.

Harnessing systematic protein-ligand interaction fingerprints for drug discovery

Determining protein-ligand interaction characteristics and mechanisms is crucial to the drug discovery process. Here, we review recent progress and successful applications of a systematic protein-ligand interaction fingerprint (IFP) approach for investigating proteome-wide protein-ligand interactions for drug development. Specifically, we review the use of this IFP approach for revealing polypharmacology across the kinome, predicting promising targets from which to design allosteric inhibitors and covalent kinase inhibitors, uncovering the binding mechanisms of drugs of interest, and demonstrating resistant mechanisms of specific drugs. Together, we demonstrate that the IFP strategy is efficient and practical for drug design research for protein kinases as targets and is extensible to other protein families.

Prediction of protein-protein interaction sites in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) participate in many biological processes by interacting with other proteins, including the regulation of transcription, translation, and the cell cycle. With the increasing amount of disorder sequence data available, it is thus crucial to identify the IDP binding sites for functional annotation of these proteins. Over the decades, many computational approaches have been developed to predict protein-protein binding sites of IDP (IDP-PPIS) based on protein sequence information. Moreover, there are new IDP-PPIS predictors developed every year with the rapid development of artificial intelligence. It is thus necessary to provide an up-to-date overview of these methods in this field. In this paper, we collected 30 representative predictors published recently and summarized the databases, features and algorithms. We described the procedure how the features were generated based on public data and used for the prediction of IDP-PPIS, along with the methods to generate the feature representations. All the predictors were divided into three categories: scoring functions, machine learning-based prediction, and consensus approaches. For each category, we described the details of algorithms and their performances. Hopefully, our manuscript will not only provide a full picture of the status quo of IDP binding prediction, but also a guide for selecting different methods. More importantly, it will shed light on the inspirations for future development trends and principles.

Computational methods for calculation of protein-ligand binding affinities in structure-based drug design

Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.