Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Synergizing drug repurposing and target identification for neurodegenerative diseases

Despite dedicated research efforts, the absence of disease-curing remedies for neurodegenerative diseases (NDDs) continues to jeopardize human society and stands as a challenge. Drug repurposing is an attempt to find new functionality of existing drugs and take it as an opportunity to discourse the clinically unmet need to treat neurodegeneration. However, despite applying this approach to rediscover a drug, it can also be used to identify the target on which a drug could work. The primary objective of target identification is to unravel all the possibilities of detecting a new drug or repurposing an existing drug. Lately, scientists and researchers have been focusing on specific genes, a particular site in DNA, a protein, or a molecule that might be involved in the pathogenesis of the disease. However, the new era discusses directing the signaling mechanism involved in the disease progression, where receptors, ion channels, enzymes, and other carrier molecules play a huge role. This review aims to highlight how target identification can expedite the whole process of drug repurposing. Here, we first spot various target-identification methods and drug-repositioning studies, including drug-target and structure-based identification studies. Moreover, we emphasize various drug repurposing approaches in NDDs, namely, experimental-based, mechanism-based, and in silico approaches. Later, we draw attention to validation techniques and stress on drugs that are currently undergoing clinical trials in NDDs. Lastly, we underscore the future perspective of synergizing drug repurposing and target identification in NDDs and present an unresolved question to address the issue.

Perspectives on Ligand/Protein Binding Kinetics Simulations: Force Fields, Machine Learning, Sampling, and User-Friendliness

Computational techniques applied to drug discovery have gained considerable popularity for their ability to filter potentially active drugs from inactive ones, reducing the time scale and costs of preclinical investigations. The main focus of these studies has historically been the search for compounds endowed with high affinity for a specific molecular target to ensure the formation of stable and long-lasting complexes. Recent evidence has also correlated the in vivo drug efficacy with its binding kinetics, thus opening new fascinating scenarios for ligand/protein binding kinetic simulations in drug discovery. The present article examines the state of the art in the field, providing a brief summary of the most popular and advanced ligand/protein binding kinetics techniques and evaluating their current limitations and the potential solutions to reach more accurate kinetic models. Particular emphasis is put on the need for a paradigm change in the present methodologies toward ligand and protein parametrization, the force field problem, characterization of the transition states, the sampling issue, and algorithms' performance, user-friendliness, and data openness.

Predicting Protein-Ligand Binding and Unbinding Kinetics with Biased MD Simulations and Coarse-Graining of Dynamics: Current State and Challenges

The prediction of drug-target binding and unbinding kinetics that occur on time scales between milliseconds and several hours is a prime challenge for biased molecular dynamics simulation approaches. This Perspective gives a concise summary of the theory and the current state-of-the-art of such predictions via biased simulations, of insights into the molecular mechanisms defining binding and unbinding kinetics as well as of the extraordinary challenges predictions of ligand kinetics pose in comparison to binding free energy predictions.

In silico prediction of siRNA ionizable-lipid nanoparticles In vivo efficacy: Machine learning modeling based on formulation and molecular descriptors

In silico prediction of the in vivo efficacy of siRNA ionizable-lipid nanoparticles is desirable as it can save time and resources dedicated to wet-lab experimentation. This study aims to computationally predict siRNA nanoparticles in vivo efficacy. A data set containing 120 entries was prepared by combining molecular descriptors of the ionizable lipids together with two nanoparticles formulation characteristics. Input descriptor combinations were selected by an evolutionary algorithm. Artificial neural networks, support vector machines and partial least squares regression were used for QSAR modeling. Depending on how the data set is split, two training sets and two external validation sets were prepared. Training and validation sets contained 90 and 30 entries respectively. The results showed the successful predictions of validation set log (siRNA dose) with R_val ²= 0.86-0.89 and 0.75-80 for validation sets one and two, respectively. Artificial neural networks resulted in the best R_val ² for both validation sets. For predictions that have high bias, improvement of R_val ² from 0.47 to 0.96 was achieved by selecting the training set lipids lying within the applicability domain. In conclusion, in vivo performance of siRNA nanoparticles was successfully predicted by combining cheminformatics with machine learning techniques.

Therapeutic drug repositioning with special emphasis on neurodegenerative diseases: Threats and issues

Drug repositioning or repurposing is the process of discovering leading-edge indications for authorized or declined/abandoned molecules for use in different diseases. This approach revitalizes the traditional drug discovery method by revealing new therapeutic applications for existing drugs. There are numerous studies available that highlight the triumph of several drugs as repurposed therapeutics. For example, sildenafil to aspirin, thalidomide to adalimumab, and so on. Millions of people worldwide are affected by neurodegenerative diseases. According to a 2021 report, the Alzheimer's disease Association estimates that 6.2 million Americans are detected with Alzheimer's disease. By 2030, approximately 1.2 million people in the United States possibly acquire Parkinson's disease. Drugs that act on a single molecular target benefit people suffering from neurodegenerative diseases. Current pharmacological approaches, on the other hand, are constrained in their capacity to unquestionably alter the course of the disease and provide patients with inadequate and momentary benefits. Drug repositioning-based approaches appear to be very pertinent, expense- and time-reducing strategies for the enhancement of medicinal opportunities for such diseases in the current era. Kinase inhibitors, for example, which were developed for various oncology indications, demonstrated significant neuroprotective effects in neurodegenerative diseases. This review expounds on the classical and recent examples of drug repositioning at various stages of drug development, with a special focus on neurodegenerative disorders and the aspects of threats and issues viz. the regulatory, scientific, and economic aspects.

Ligand Unbinding Pathway and Mechanism Analysis Assisted by Machine Learning and Graph Methods

We present two methods to reveal protein-ligand unbinding mechanisms in biased unbinding simulations by clustering trajectories into ensembles representing unbinding paths. The first approach is based on a contact principal component analysis for reducing the dimensionality of the input data, followed by identification of unbinding paths and training a machine learning model for trajectory clustering. The second approach clusters trajectories according to their pairwise mean Euclidean distance employing the neighbor-net algorithm, which takes into account input data bias in the distances set and is superior to dendrogram construction. Finally, we describe a more complex case where the reaction coordinate relevant for path identification is a single intraligand hydrogen bond, highlighting the challenges involved in unbinding path reaction coordinate detection.

Public Data Set of Protein-Ligand Dissociation Kinetic Constants for Quantitative Structure-Kinetics Relationship Studies

Protein-ligand binding affinity reflects the equilibrium thermodynamics of the protein-ligand binding process. Binding/unbinding kinetics is the other side of the coin. Computational models for interpreting the quantitative structure-kinetics relationship (QSKR) aim at predicting protein-ligand binding/unbinding kinetics based on protein structure, ligand structure, or their complex structure, which in principle can provide a more rational basis for structure-based drug design. Thus far, most of the public data sets used for deriving such QSKR models are rather limited in sample size and structural diversity. To tackle this problem, we have compiled a set of 680 protein-ligand complexes with experimental dissociation rate constants (k off), which were mainly curated from the references accumulated for updating our PDBbind database. Three-dimensional structure of each protein-ligand complex in this data set was either retrieved from the Protein Data Bank or carefully modeled based on a proper template. The entire data set covers 155 types of protein, with their dissociation kinetic constants (k off) spanning nearly 10 orders of magnitude. To the best of our knowledge, this data set is the largest of its kind reported publicly. Utilizing this data set, we derived a random forest (RF) model based on protein-ligand atom pair descriptors for predicting k off values. We also demonstrated that utilizing modeled structures as additional training samples will benefit the model performance. The RF model with mixed structures can serve as a baseline for testifying other more sophisticated QSKR models. The whole data set, namely, PDBbind-koff-2020, is available for free download at our PDBbind-CN web site (http://www.pdbbind.org.cn/download.php).

Determination of Multidirectional Pathways for Ligand Release from the Receptor: A New Approach Based on Differential Evolution

Steered molecular dynamics (SMD) simulation is a powerful method in computer-aided drug design as it can be used to access the relative binding affinity with high precision but with low computational cost. The success of SMD depends on the choice of the direction along which the ligand is pulled from the receptor-binding site. In most simulations, the unidirectional pathway was used, but in some cases, this choice resulted in the ligand colliding with the complex surface of the exit tunnel. To overcome this difficulty, several variants of SMD with multidirectional pulling have been proposed, but they are not completely devoid of disadvantages. Here, we have proposed to determine the direction of pulling with a simple scoring function that minimizes the receptor–ligand interaction, and an optimization algorithm called differential evolution is used for energy minimization. The effectiveness of our protocol was demonstrated by finding expulsion pathways of Huperzine A and camphor from the binding site of Torpedo California acetylcholinesterase and P450cam proteins, respectively, and comparing them with the previous results obtained using memetic sampling and random acceleration molecular dynamics. In addition, by applying this protocol to a set of ligands bound with LSD1 (lysine specific demethylase 1), we obtained a much higher correlation between the work of pulling force and experimental data on the inhibition constant IC50 compared to that obtained using the unidirectional approach based on minimal steric hindrance.

Applications of machine learning in computer-aided drug discovery

Machine learning (ML) has revolutionised the field of structure-based drug design (SBDD) in recent years. During the training stage, ML techniques typically analyse large amounts of experimentally determined data to create predictive models in order to inform the drug discovery process. Deep learning (DL) is a subfield of ML, that relies on multiple layers of a neural network to extract significantly more complex patterns from experimental data, and has recently become a popular choice in SBDD. This review provides a thorough summary of the recent DL trends in SBDD with a particular focus on de novo drug design, binding site prediction, and binding affinity prediction of small molecules.

Exploring the thermodynamic, kinetic and inhibitory mechanisms of 5-iTU targeting mitotic kinase haspin by integrated molecular dynamics

As a human mitotic kinase, haspin is considered as a promising target for various diseases including cancers. However, no inhibitors targeting haspin have entered clinical trials presently. 5-iTU (5-iodotubercidin) is a useful and classical chemical probe for the investigation of haspin activity, but its inhibitory mechanism remains unclear. In this study, integrated molecular dynamics (MD) of conventional MD, extended adaptive biasing force (eABF), random acceleration MD and well-tempered metadynamics were applied to investigate the thermodynamic and kinetic features of 5-iTU and three derivatives targeting haspin. To emphasize the importance of gatekeeper Phe605, two haspin mutants (F605Y and F605T) were also built. The results showed that the binding affinity of 5-iTU and haspin was highest in all wild type (WT) systems, relying on the strong halogen aromatic π interaction between 5-iTU and gatekeeper Phe605. Gatekeeper mutations, because of damage to this interaction, led to the rearrangement of water distributions at the binding site and the decrease of 5-iTU residence times. Additionally, compared with the smaller 5-fTU, 5-iTU dissociated from WT haspin with more difficulty through distinct unbinding pathways. These findings will provide crucial guidance for the design and development of novel haspin inhibitors and the rational modification of existing inhibitors.

Rapid Quantification of Protein-Ligand Binding via 19F NMR Lineshape Analysis

Fluorine incorporation is ideally suited to many NMR techniques, and incorporation of fluorine into proteins and fragment libraries for drug discovery has become increasingly common. Here, we use one-dimensional 19F NMR lineshape analysis to quantify the kinetics and equilibrium thermodynamics for the binding of a fluorine-labeled Src homology 3 (SH3) protein domain to four proline-rich peptides. SH3 domains are one of the largest and most well-characterized families of protein recognition domains and have a multitude of functions in eukaryotic cell signaling. First, we showe that fluorine incorporation into SH3 causes only minor structural changes to both the free and bound states using amide proton temperature coefficients. We then compare the results from lineshape analysis of one-dimensional 19F spectra to those from two-dimensional 1H-15N heteronuclear single quantum coherence spectra. Their agreement demonstrates that one-dimensional 19F lineshape analysis is a robust, low-cost, and fast alternative to traditional heteronuclear single quantum coherence-based experiments. The data show that binding is diffusion limited and indicate that the transition state is highly similar to the free state. We also measured binding as a function of temperature. At equilibrium, binding is enthalpically driven and arises from a highly positive activation enthalpy for association with small entropic contributions. Our results agree with those from studies using different techniques, providing additional evidence for the utility of 19F NMR lineshape analysis, and we anticipate that this analysis will be an effective tool for rapidly characterizing the energetics of protein interactions.

In Silico Prediction of the Dissociation Rate Constants of Small Chemical Ligands by 3D-Grid-Based VolSurf Method

Accumulated evidence suggests that binding kinetic properties—especially dissociation rate constant or drug-target residence time—are crucial factors affecting drug potency. However, quantitative prediction of kinetic properties has always been a challenging task in drug discovery. In this study, the VolSurf method was successfully applied to quantitatively predict the k_off values of the small ligands of heat shock protein 90α (HSP90α), adenosine receptor (AR) and p38 mitogen-activated protein kinase (p38 MAPK). The results showed that few VolSurf descriptors can efficiently capture the key ligand surface properties related to dissociation rate; the resulting models demonstrated to be extremely simple, robust and predictive in comparison with available prediction methods. Therefore, it can be concluded that the VolSurf-based prediction method can be widely applied in the ligand-receptor binding kinetics and de novo drug design researches.