Unveiling the Selectivity Mechanism of Type-I LRRK2 Inhibitors by Computational Methods: Insights from Binding Thermodynamics and Kinetics Simulation.

Understanding the selectivity mechanisms of inhibitors toward highly similar proteins is very important in new drug discovery. Developing highly selective targeting of leucine-rich repeat kinase 2 (LRRK2) kinases for the treatment of Parkinson's disease (PD) is challenging because of the similarity of the kinase ATP binding pocket. During the development of LRRK2 inhibitors, off-target effects on other kinases, especially TTK and JAK2 kinases, have been observed. As a result, significant time and resources have been devoted to improving the selectivity for the LRRK2 target. DNL201 is an LRRK2 kinase inhibitor entering phase I clinical studies. The experiments have shown that DNL201 significantly inhibits LRRK2 kinase activity, with >167-fold selectivity over JAK2 and TTK kinases. However, the potential mechanisms of inhibitor preferential binding to LRRK2 kinase are still not well elucidated. In this work, to reveal the underlying general selectivity mechanism, we carried out several computational methods and comprehensive analyses from both the binding thermodynamics and kinetics on two representative LRRK2 inhibitors (DNL201 and GNE7915) to LRRK2. Our results suggest that the structural and kinetic differences between the proteins may play a key role in determining the activity of the selective small-molecule inhibitor. The selectivity mechanisms proposed in this work could be helpful for the rational design of novel selective LRRK2 kinase inhibitors against PD.

Improving drug-target affinity prediction via feature fusion and knowledge distillation.

Rapid and accurate prediction of drug-target affinity can accelerate and improve the drug discovery process. Recent studies show that deep learning models may have the potential to provide fast and accurate drug-target affinity prediction. However, the existing deep learning models still have their own disadvantages that make it difficult to complete the task satisfactorily. Complex-based models rely heavily on the time-consuming docking process, and complex-free models lacks interpretability. In this study, we introduced a novel knowledge-distillation insights drug-target affinity prediction model with feature fusion inputs to make fast, accurate and explainable predictions. We benchmarked the model on public affinity prediction and virtual screening dataset. The results show that it outperformed previous state-of-the-art models and achieved comparable performance to previous complex-based models. Finally, we study the interpretability of this model through visualization and find it can provide meaningful explanations for pairwise interaction. We believe this model can further improve the drug-target affinity prediction for its higher accuracy and reliable interpretability.

Development, validation, and evaluation of a deep learning model to screen cyclin-dependent kinase 12 inhibitors in cancers.

Deep learning-based in silico alternatives have been demonstrated to be of significant importance in the acceleration of the drug discovery process and enhancement of success rates. Cyclin-dependent kinase 12 (CDK12) is a transcription-related cyclin-dependent kinase that may act as a biomarker and therapeutic target for cancers. However, currently, there is no high selective CDK12 inhibitor in clinical development and the identification of new specific CDK12 inhibitors has become increasingly challenging due to their similarity with CDK13. In this study, we developed a virtual screening workflow that combines deep learning with virtual screening tools and can be applied rapidly to millions of molecules. We designed a Transformer architecture Drug-Target Interaction (DTI) model with dual-branched self-supervised pre-trained molecular graph models and protein sequence models. Our predictive model produced satisfactory predictions for various targets, including CDK12, with several novel hits. We screened a large compound library consisting of 4.5 million drug-like molecules and recommended a list of potential CDK12 inhibitors for further experimental testing. In kinase assay, compared to the positive CDK12 inhibitor THZ531, the compounds CICAMPA-01, 02, 03 displayed more effective inhibition of CDK12, up to three times as much as THZ531. The compounds CICAMPA-03, 05, 04, 07 showed less inhibition of CDK13 compare to THZ531. In vitro, the IC50 of CICAMPA-01, 04, 05, 06, 09 was less than 3 µM in the HER2 positive CDK12 amplification breast cancer cell line BT-474. Overall, this study provides a highly efficient and end-to-end deep learning protocol, in conjunction with molecular docking, for discovering CDK12 inhibitors in cancers. Additionally, we disclose five novel CDK12 inhibitors. These results may accelerate the discovery of novel chemical-class drugs for cancer treatment.

Identification of LRRK2 Inhibitors through Computational Drug Repurposing.

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

Binding Thermodynamics and Dissociation Kinetics Analysis Uncover the Key Structural Motifs of Phenoxyphenol Derivatives as the Direct InhA Inhibitors and the Hotspot Residues of InhA.

Given the current epidemic of multidrug-resistant tuberculosis, there is an urgent need to develop new drugs to combat drug-resistant tuberculosis. Direct inhibitors of the InhA target do not require activation and thus can overcome drug resistance caused by mutations in drug-activating enzymes. In this work, the binding thermodynamic and kinetic information of InhA to its direct inhibitors, phenoxyphenol derivatives, were explored through multiple computer-aided drug design (CADD) strategies. The results show that the van der Waals interactions were the main driving force for protein-ligand binding, among which hydrophobic residues such as Tyr158, Phe149, Met199 and Ile202 have high energy contribution. The AHRR pharmacophore model generated by multiple ligands demonstrated that phenoxyphenol derivatives inhibitors can form pi-pi stacking and hydrophobic interactions with InhA target. In addition, the order of residence time predicted by random acceleration molecular dynamics was consistent with the experimental values. The intermediate states of these inhibitors could form hydrogen bonds and van der Waals interactions with surrounding residues during dissociation. Overall, the binding and dissociation mechanisms at the atomic level obtained in this work can provide important theoretical guidance for the development of InhA direct inhibitors with higher activity and proper residence time.

Molecular Modeling Study on the Interaction Mechanism between the LRRK2 G2019S Mutant and Type I Inhibitors by Integrating Molecular Dynamics Simulation, Binding Free Energy Calculations, and Pharmacophore Modeling.

Leucine-rich repeat kinase 2 (LRRK2) has been reported in the pathogenesis of Parkinson's disease (PD). G2019S mutant is the most common pathogenic mutation in LRRK2-related PD patients. Inhibition of LRRK2 kinase activity is proposed to be a new therapeutic approach for PD treatment. Therefore, understanding the molecular basis of the interaction between LRRK2 and its inhibitors will be valuable for the discovery and design of LRRK2 inhibitors. However, the structure of human LRRK2 in complex with the inhibitor has not been determined, and the inhibitory mechanism underlying LRRK2 still needs to be further investigated. In this study, molecular dynamics (MD) simulation combined with the molecular mechanics generalized born surface area (MM-GBSA) binding free energy calculation and pharmacophore modeling methods was employed to explore the critical residues in LRRK2 for binding of inhibitors and to investigate the general structural features of the inhibitors with diverse scaffolds. The results from MD simulations suggest that the hinge region residues Glu1948 and Ala1950 play a significant role in maintaining the intermolecular hydrogen bond interaction with the G2019S LRRK2 protein and inhibitor. The strong hinge hydrogen bond with an occupancy rate of more than 95% represents the high activity of LRRK2 inhibitors, and the hydrogen bond interaction with the kinase catalytic loop region could compromise selectivity. Further pharmacophore modeling reveals that the high activity LRRK2 inhibitor should have one aromatic ring, one hydrogen bond acceptor, and one hydrogen bond donor. Hence, the obtained results can provide valuable information to understand the interactions of LRRK2 inhibitors at the atomic level that will be helpful in designing potent inhibitors of LRRK2.

Virtual Screening and Biological Activity Evaluation of New Potent Inhibitors Targeting LRRK2 Kinase Domain.

Leucine rich repeat kinase 2 (LRRK2) has been reported in the pathogenesis of Parkinson's disease (PD). Inhibition of LRRK2 kinase activity is a therapeutic approach that may provide new treatments for PD. In this study, novel LRRK2 inhibitors were identified by performing a docking-based virtual screening (VS). Due to the absence of a crystal structure of LRRK2, homology modeling was adopted to model human LRRK2 kinase domain that binds the inhibitor. Next, a docking-based virtual screening protocol was applied to identify LRRK2 small molecule inhibitors targeting the ATP binding pocket. A total of 28 compounds were selected and subjected to LRRK2 kinase inhibition assay. As a result, two small molecules with novel skeleton, compounds LY2019-005 and LY2019-006, were identified as potential LRRK2 kinase inhibitors with the IC50 of these two compounds against the wild-type and G2019S mutant LRRK2 kinase being 424.40 ± 1.31 nM, 378.80 ± 1.20 nM and 1526.00 ± 0.87 nM, 1165.00 ± 1.18 nM, respectively. Molecular dynamics (MD) simulation was carried out to reveal the binding mode of the newly identified compound LY2019-005 to the LRRK2 kinase domain. The binding modes indicate that the important hydrogen bond between hinge region (such as Ala1950) and inhibitor is crucial for the inhibition activity. In summary, our study provides a highly efficient way to discover LRRK2 inhibitors, and we find two highly efficient novel LRRK2 inhibitors, which could be helpful for the development of potential drugs targeting LRRK2 in PD therapy.

MolAICal: a soft tool for 3D drug design of protein targets by artificial intelligence and classical algorithm.

Deep learning is an important branch of artificial intelligence that has been successfully applied into medicine and two-dimensional ligand design. The three-dimensional (3D) ligand generation in the 3D pocket of protein target is an interesting and challenging issue for drug design by deep learning. Here, the MolAICal software is introduced to supply a way for generating 3D drugs in the 3D pocket of protein targets by combining with merits of deep learning model and classical algorithm. The MolAICal software mainly contains two modules for 3D drug design. In the first module of MolAICal, it employs the genetic algorithm, deep learning model trained by FDA-approved drug fragments and Vinardo score fitting on the basis of PDBbind database for drug design. In the second module, it uses deep learning generative model trained by drug-like molecules of ZINC database and molecular docking invoked by Autodock Vina automatically. Besides, the Lipinski's rule of five, Pan-assay interference compounds (PAINS), synthetic accessibility (SA) and other user-defined rules are introduced for filtering out unwanted ligands in MolAICal. To show the drug design modules of MolAICal, the membrane protein glucagon receptor and non-membrane protein SARS-CoV-2 main protease are chosen as the investigative drug targets. The results show MolAICal can generate the various and novel ligands with good binding scores and appropriate XLOGP values. We believe that MolAICal can use the advantages of deep learning model and classical programming for designing 3D drugs in protein pocket. MolAICal is freely for any nonprofit purpose and accessible at https://molaical.github.io.

Probing the Molecular Mechanism of Rifampin Resistance Caused by the Point Mutations S456L and D441V on Mycobacterium tuberculosis RNA Polymerase through Gaussian Accelerated Molecular Dynamics Simulation.

Rifampin is the first-line antituberculosis drug, with Mycobacterium tuberculosis RNA polymerase as the molecular target. Unfortunately, M. tuberculosis strains that are resistant to rifampin have been identified in clinical settings, which limits its therapeutic effects. In clinical isolates, S531L and D516V (in Escherichia coli) are two common mutated codons in the gene rpoB, corresponding to S456L and D441V in M. tuberculosis However, the resistance mechanism at the molecular level is still elusive. In this work, Gaussian accelerated molecular dynamics simulations were performed to uncover the resistance mechanism of rifampin due to S456L and D441V mutations at the atomic level. The binding free energy analysis revealed that the reduction in the ability of two mutants to bind rifampin is mainly due to a decrease in electrostatic interaction, specifically, a decrease in the energy contribution of the R454 residue. R454 acts as an anchor and forms stable hydrogen bond interaction with rifampin, allowing rifampin to be stably incorporated in the center of the binding pocket. However, the disappearance of the hydrogen bond between R454 and the mutated residues increases the flexibility of the side chain of R454. The conformation of R454 changes, and the hydrogen bond interaction between it and rifampin is disrupted. As result, the rifampin molecule moves to the outside of the pocket, and the binding affinity decreases. Overall, these findings can provide useful information for understanding the drug resistance mechanism of rifampin and also can give theoretical guidance for further design of novel inhibitors to overcome the drug resistance.

Disclosing the Mechanism of Spontaneous Aggregation and Template-Induced Misfolding of the Key Hexapeptide (PHF6) of Tau Protein Based on Molecular Dynamics Simulation.

The microtubule-associated protein tau is critical for the development and maintenance of the nervous system. Tau dysfunction is associated with a variety of neurodegenerative diseases called tauopathies, which are characterized by neurofibrillary tangles formed by abnormally aggregated tau protein. Studying the aggregation mechanism of tau protein is of great significance for elucidating the etiology of tauopathies. The hexapeptide 306VQIVYK311 (PHF6) of R3 has been shown to play a vital role in promoting tau aggregation. In this study, long-term all-atom molecular dynamics simulations in explicit solvent were performed to investigate the mechanisms of spontaneous aggregation and template-induced misfolding of PHF6, and the dimerization at the early stage of nucleation was further specifically analyzed by the Markov state model (MSM). Our results show that PHF6 can spontaneously aggregate to form multimers enriched with ß-sheet structure and the ß-sheets in multimers prefer to exist in a parallel way. It is observed that PHF6 monomer can be induced to form a ß-sheet structure on either side of the template but in a different way. In detail, the ß-sheet structure is easier to form on the left side but does not extend well, but on the right side, the monomer can form the extended ß-sheet structure. Furthermore, MSM analysis shows that the formation of dimer mainly occurs in three steps. First, the separated monomers collide with each other at random orientations, and then a dimer with short ß-sheet structure at the N-terminal forms; finally, ß-sheets elongate to form an extended parallel ß-sheet dimer. During these processes, multiple intermediate states are identified and multiple paths can form a parallel ß-sheet dimer from the disordered coil structure. Moreover, the residues I308, V309, and Y310 play an essential role in the dimerization. In a word, our results uncover the aggregation and misfolding mechanism of PHF6 from the atomic level, which can provide useful theoretical guidance for rational design of effective therapeutic drugs against tauopathies.

Conformation Transition of Intracellular Part of Glucagon Receptor in Complex With Agonist Glucagon by Conventional and Accelerated Molecular Dynamics Simulations.

The inactive conformations of glucagon receptor (GCGR) are widely reported by crystal structures that support the precision structure for drug discovery of type 2 diabetes. The previous study shows that the intracellular part is open in the glucagon-bound GCGR (glu-GCGR) and closed in the apo-GCGR by accelerated molecular dynamics (aMD) simulations. However, the crystal structure of GCGR in complex with partial agonist shows that the intracellular part is closed in the inactive conformation. To understand the differences between the studies of aMD simulations and crystal structure, the 2,500 ns conventional molecular dynamics (cMD) simulations are performed on the simulated model of glu-GCGR. The result shows that the transmembrane helices (TMH) 6 of glu-GCGR is outward ~4 Å to drive the intracellular part of glu-GCGR open until ~390 ns cMD simulations. The (TMH) 6 of glu-GCGR becomes closed after ~490 ns cMD simulations, which are consistent with the crystal structure of GCGR in complex with the partial agonist. To further elucidate the activation mechanism of GCGR deeply, the simulated models of glu-GCGR, apo-GCGR, and antagonist-bound GCGR (ant-GCGR) are constructed to perform 10 of parallel 300 ns aMD simulations, respectively. The results show that both of glu-GCGR and apo-GCGR can generate the open conformations of the intracellular part. But the glu-GCGR has the higher percentage of open conformations than apo-GCGR. The ant-GCGR is restricted to generate the open conformations of the intracellular part by antagonist MK-0893. It indicates that the glu-GCGR, apo-GCGR, and ant-GCGR can be distinguished by the aMD simulated method. Free energy landscape shows that the open conformations of the intracellular part of GCGR are in intermediate state. Our results show that aMD simulations enhance the space samplings of open conformations of GCGR via adding extra boost potential. It indicates that the aMD simulations are an effective way for drug discovery of GCGR.