Identification of Selective Dual ROCK1 and ROCK2 Inhibitors Using Structure-Based Drug Design

Caught between a ROCK and a hard place: current challenges in structure-based drug design

The discipline of structure-based drug design (SBDD) is several decades old and it is tempting to think that the proliferation of experimental structures for many drug targets might make computer-aided drug design (CADD) straightforward. However, this is far from true. In this review, we illustrate some of the challenges that CADD scientists face every day in their work, even now. We use Rho-associated protein kinase (ROCK), and public domain structures and data, as an example to illustrate some of the challenges we have experienced during our project targeting this protein. We hope that this will help to prevent unrealistic expectations of what CADD can accomplish and to educate non-CADD scientists regarding the challenges still facing their CADD colleagues.

Discovery of a novel ROCK2 ATP competitive inhibitor by DNA-encoded library selection

Neurodegeneration disorders, such as Alzheimer's disease (AD), have garnered significant attention due to their impact on individuals and society as a whole. Understanding the mechanisms behind these disorders and developing effective therapy strategies is of utmost importance. One potential therapeutic target that has emerged is Rho-associated coiled-coil containing protein kinase 2 (ROCK2), as its accumulation and activity have been closely linked to memory loss. In this report, we present the findings of a recent discovery involving a new molecule that has the ability to competitively inhibit ROCK2 activity. This molecule was identified through the utilization of a DNA-encoded library (DEL) screening platform. Following selection against ROCK2, an off-DNA compound was synthesized and examined to ascertain its inhibitory properties, selectivity, mechanism of action, and binding mode analysis. From the screening, compound CH-2 has demonstrated an IC50 value of 28 nM against ROCK2, while exhibiting a 5-fold selectivity over ROCK1. Further analysis through molecular docking has provided insights into the specific binding modes of this compound. Our findings suggest that DEL selection offers a rapid method for identifying new inhibitors. Among these, the CH-2 compound shows promise as a potential ROCK2 inhibitor and warrants further investigation.

Target-based drug design strategies to overcome resistance to antiviral agents: opportunities and challenges

Viral infections have a major impact in human health. Ongoing viral transmission and escalating selective pressure have the potential to favor the emergence of vaccine- and antiviral drug-resistant viruses. Target-based approaches for the design of antiviral drugs can play a pivotal role in combating drug-resistant challenges. Drug design computational tools facilitate the discovery of novel drugs. This review provides a comprehensive overview of current drug design strategies employed in the field of antiviral drug resistance, illustrated through the description of a series of successful applications. These strategies include technologies that enhance compound-target affinity while minimizing interactions with mutated binding pockets. Furthermore, emerging approaches such as virtual screening, targeted protein/RNA degradation, and resistance analysis during drug design have been harnessed to curtail the emergence of drug resistance. Additionally, host targeting antiviral drugs offer a promising avenue for circumventing viral mutation. The widespread adoption of these refined drug design strategies will effectively address the prevailing challenge posed by antiviral drug resistance.

Discovery of dual rho-associated protein kinase 1 (ROCK1)/apoptosis signal-regulating kinase 1 (ASK1) inhibitors as a novel approach for non-alcoholic steatohepatitis (NASH) treatment

In the current study we suggest a novel approach to curb non-alcoholic steatohepatitis (NASH) progression, and we suggest privileged scaffolds for the design of novel compounds for this aim. NASH is an advanced form of non-alcoholic fatty liver disease that can further progress into fibrosis, cirrhosis, and hepatocellular carcinoma. It is a widely emerging disease affecting 25% of the global population and has no current approved treatments. Protein kinases are key regulators of cellular pathways, of which, Rho-associated protein kinase 1 (ROCK1) and apoptosis signal-regulating kinase 1 (ASK1) play an important role in the progression of NASH and they stand out as promising targets for NASH therapy. Interestingly, their kinase domains are found to be similar in sequence and topology; therefore, dual inhibition of ROCK1 and ASK1 is expected to be amenable and could achieve a more favourable outcome. To reach this goal, a training set of ROCK1 and ASK1 protein structures co-crystalized with type 1 (ATP-competitive) inhibitors was constructed to manually generate receptor-based pharmacophore models representing ROCK1 and ASK1 inhibitors' common pharmacophoric features. The models produced were assessed using a test set of both ROCK1 and ASK1 actives and decoys, and their performance was evaluated using different assessment metrics. The best pharmacophore model obtained, showing a Mathew's correlation coefficient (MCC) of 0.71, was then used to screen the ZINC purchasable database retrieving 6178 hits that were filtered accordingly using several medicinal chemistry and pharmacokinetics filters returning 407 promising compounds. To confirm that these compounds are capable of binding to the target kinases, they were subjected to molecular docking simulations at both protein structures. The results were then assessed individually and filtered, setting the spotlight on various privileged scaffolds that could be exploited as the nucleus for designing novel ROCK1/ASK1 dual inhibitors.

Recent advances in the development of Rho kinase inhibitors (2015-2021)

Rho-associated coiled-coil kinases (ROCKs) are key downstream effectors of small GTPases. ROCK plays a central role in diverse cellular events with accumulating evidence supporting the concept that ROCK is important in tumor development and progression. Numerous ROCK inhibitors have been investigated for their therapeutic potential in the treatment of cancers. In this article, we review recent research progress on ROCK inhibitors, especially those with potential for the treatment of cancers, reported in the literature from 2015 to 2021. Most ROCK inhibitors show potent in vitro and in vivo antitumor activities and have potential in the treatment of cancers.

In Vitro and In Silico Investigation of BCI Anticancer Properties and Its Potential for Chemotherapy-Combined Treatments

BACKGROUND

DUSP6 phosphatase serves as a negative regulator of MAPK kinases involved in numerous cellular processes. BCI has been identified as a potential allosteric inhibitor with anticancer activity. Our study was designed to test the anticancer properties of BCI in colon cancer cells, to characterize the effect of this compound on chemotherapeutics such as irinotecan and oxaliplatin activity, and to identify potential molecular targets for this inhibitor.

METHODS

BCI cytotoxicity, proapoptotic activity, and cell cycle distribution were investigated in vitro on three colon cancer cell lines (DLD1, HT-29, and Caco-2). In silico investigation was prepared to assess BCI drug-likeness and identify potential molecular targets.

RESULTS

The exposure of colorectal cancer cells with BCI resulted in antitumor effects associated with cell cycle arrest and induction of apoptosis. BCI exhibited strong cytotoxicity on DLD1, HT-29, and Caco-2 cells. BCI showed no significant interaction with irinotecan, but strongly attenuated the anticancer activity of oxaliplatin when administered together. Analysis of synergy potential further confirmed the antagonistic interaction between these two compounds. In silico investigation indicated CDK5 as a potential new target of BCI.

CONCLUSIONS

Our studies point to the anticancer potential of BCI but note the need for a precise mechanism of action.

In silico screening-based discovery of benzamide derivatives as inhibitors of Rho-associated kinase-1 (ROCK1)

As a pivotal node in modulating various cell behaviors, Rho-associated kinase-1 (ROCK1) has attracted significant attention as a promising therapeutic target in a variety of diseases. Benzamide has been widely reported as a ROCK1 inhibitors in recent years. To better understand its pharmacological properties and to explore its potential inhibitors, a series of ROCK1 inhibitors derived from N-methyl-4-(4-pyrazolidinyl) benzamides (MPBs) were investigated by using three-dimensional quantitative structure-activity relationship (3D-QSAR) models, pharmacophore models, molecular docking, and molecular dynamics (MD) simulation. The comparative Molecular Field Analysis (CoMFA) model (q2 = 0.616, R2 = 0.972, ONC = 4, and r2pred = 0.983) and the best Comparative Molecular Similarity Indices Analysis (CoMSIA) model (q2= 0.740, R2 = 0.982, ONC = 6, and r2pred = 0.824) exhibited reliable predictability with satisfactory validation parameters. In the subsequent virtual screening, VS03 and VS05 were identified to have superior predicted activities and higher docking scores, meanwhile they demonstrated to be reasonably stable in the binding pocket through MD simulations. These results provide a significant theoretical direction for the rational design and development of novel ROCK1 inhibitors.Communicated by Ramaswamy H. Sarma.

Structure-based virtual screening of ROCK1 inhibitors for the discovery of Enterovirus-A71 antivirals

Human enterovirus A71 (EV-A71) is the major causative agent of hand, foot, and mouth disease (HFMD), which may lead to neurological sequelae and even death. Although EV-A71 seriously threatens public health, there remains no efficient drug for the treatment of EV-A71 infection. We previously demonstrated that ROCK1 is a novel host dependency factor for EV-A71 replication and can serve as a target for the development of anti-EV-A71 therapeutics. In this study, we identified a subset of inhibitors with potential anti-EV-A71 activity by virtual screening using ROCK1 as a target. Among the hits, Dasabuvir, an HCV polymerase inhibitor, was found to have the best antiviral activity which is consistent with the ranking scores in Autodock Vina and iGEMDOCK. We found that Dasabuvir efficiently suppressed EV-A71 replication in a dose-dependent manner. Moreover, Dasabuvir not only efficiently suppressed the replication of EV-A71 in RD cells, but also in multiple cell lines, including HEK-293T, Caco-2, HT-29, HepG2, and Huh7. Besides, Dasabuvir alleviated the release of proinflammatory cytokines caused by EV-A71 infection. Notably, Dasabuvir also exhibited antiviral activity of CVA10, indicating it may have broad-spectrum antiviral activity against species Enteroviruses A. Hence, our results further confirm that ROCK1 can be a potential drug target and suggest Dasabuvir could be a clinical candidate for the treatment of EV-A71 infection.

Discovery of indoline-based derivatives as effective ROCK2 inhibitors for the potential new treatment of idiopathic pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic, progressive, and devastating lung disease with a median survival of only 3-5 years. Due to the lack of effective therapy, IPF threatens human health. Recently, increasing reports have indicated that Rho-associated coiled-coil protein kinases (ROCKs) play important roles in the development of IPF and might represent a novel target for the treatment of IPF. Herein, a new series of selective ROCK2 inhibitors based on indoline were designed and synthesized. Structural modification resulted in optimized compound 9b with an IC₅₀ value of 6 nM against ROCK2 and the inhibition of collagen gel contraction. Cellular assays demonstrated that 9b could significantly suppress the expression of collagen I and α-SMA, and inhibited ROCK signaling pathway. Oral administration of compound 9b (10 mg/kg) exerted more significant anti-pulmonary fibrosis effects than nintedanib (100 mg/kg) and KD025 (100 mg/kg) in a bleomycin-induced IPF rat model, suggesting that 9b could serve as a potential lead compound for the treatment of IPF.

Application and Study of ROCK Inhibitors in Pulmonary Fibrosis: Recent Developments and Future Perspectives

Rho-associated coiled-coil-containing kinases (ROCKs), serine/threonine protein kinases, were initially identified as downstream targets of the small GTP-binding protein Rho. Pulmonary fibrosis (PF) is a lethal disease with limited therapeutic options and a particularly poor prognosis. Interestingly, ROCK activation has been demonstrated in PF patients and in animal PF models, making it a promising target for PF treatment. Many ROCK inhibitors have been discovered, and four of these have been approved for clinical use; however, no ROCK inhibitors are approved for the treatment of PF patients. In this article, we describe ROCK signaling pathways and the structure-activity relationship, potency, selectivity, binding modes, pharmacokinetics (PKs), biological functions, and recently reported inhibitors of ROCKs in the context of PF. We will also focus our attention on the challenges to be addressed when targeting ROCKs and discuss the strategy of ROCK inhibitor use in the treatment of PF.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

Novel Dual AChE and ROCK2 Inhibitor Induces Neurogenesis via PTEN/AKT Pathway in Alzheimer's Disease Model

Alzheimer's disease (AD) is a progressive and complex neurodegenerative disease. Acetylcholinesterase inhibitors (AChEIs) are a major class of drugs used in AD therapy. ROCK2, another promising target for AD, has been associated with the induction of neurogenesis via PTEN/AKT. This study aimed to characterize the therapeutic potential of a novel donepezil-tacrine hybrid compound (TA8Amino) to inhibit AChE and ROCK2 protein, leading to the induction of neurogenesis in SH-SY5Y cells. Experiments were carried out with undifferentiated and neuron-differentiated SH-SY5Y cells submitted to treatments with AChEIs (TA8Amino, donepezil, and tacrine) for 24 h or 7 days. TA8Amino was capable of inhibiting AChE at non-cytotoxic concentrations after 24 h. Following neuronal differentiation for 7 days, TA8Amino and donepezil increased the percentage of neurodifferentiated cells and the length of neurites, as confirmed by β-III-tubulin and MAP2 protein expression. TA8Amino was found to participate in the activation of PTEN/AKT signaling. In silico analysis showed that TA8Amino can stably bind to the active site of ROCK2, and in vitro experiments in SH-SY5Y cells demonstrate that TA8Amino significantly reduced the expression of ROCK2 protein, contrasting with donepezil and tacrine. Therefore, these results provide important information on the mechanism underlying the action of TA8Amino with regard to multi-target activities.

Identification of a Potent, Selective, and Brain-Penetrant Rho Kinase Inhibitor and its Activity in a Mouse Model of Huntington's Disease

The Rho kinase (ROCK) pathway is implicated in the pathogenesis of several conditions, including neurological diseases. In Huntington's disease (HD), ROCK is implicated in mutant huntingtin (HTT) aggregation and neurotoxicity, and members of the ROCK pathway are increased in HD mouse models and patients. To validate this mode of action as a potential treatment for HD, we sought a potent, selective, central nervous system (CNS)-penetrant ROCK inhibitor. Identifying a compound that could be dosed orally in mice with selectivity against other AGC kinases, including protein kinase G (PKG), whose inhibition could potentially activate the ROCK pathway, was paramount for the program. We describe the optimization of published ligands to identify a novel series of ROCK inhibitors based on a piperazine core. Morphing of the early series developed in-house by scaffold hopping enabled the identification of a compound exhibiting high potency and desired selectivity and demonstrating a robust pharmacodynamic (PD) effect by the inhibition of ROCK-mediated substrate (MYPT1) phosphorylation after oral dosing.

3D-printing magnesium-polycaprolactone loaded with melatonin inhibits the development of osteosarcoma by regulating cell-in-cell structures

Melatonin has been proposed as a potent anticarcinogen presents a short half-life for osteosarcoma (OS). Cell-in-cell (CIC) structures play a role in the development of malignant tumors by changing the tumor cell energy metabolism. This study developed a melatonin-loaded 3D printed magnesium-polycaprolactone (Mg-PCL) scaffold and investigated its effect and molecular mechanism on CIC in OS. Mg-PCL scaffold was prepared by 3D-printing and its characteristic was determined. The effect and molecular mechanism of Mg-PCL scaffold as well as melatonin-loaded Mg-PCL on OS growth and progression were investigated in vivo and in vitro. We found that melatonin receptor 1 (MT1) and CIC expressions were increased in OS tissues and cells. Melatonin treatment inhibit the key CIC pathway, Rho/ROCK, through the cAMP/PKA signaling pathway, interfering with the mitochondrial physiology of OS cells, and thus playing an anti-invasion and anti-metastasis role in OS. The Mg-PCL-MT could significantly inhibit distant organ metastasis of OS in the in vivo model. Our results showed that melatonin-loaded Mg-PCL scaffolds inhibited the proliferation, invasion and metastasis of OS cells through the CIC pathway. The Mg-PCL-MT could be a potential therapeutics for OS.

Designing of the N-ethyl-4-(pyridin-4-yl)benzamide based potent ROCK1 inhibitors using docking, molecular dynamics, and 3D-QSAR

Rho-associated kinase-1 (ROCK1) has been recognized for its pivotal role in heart diseases, different types of malignancy, and many neurological disorders. Hyperactivity of ROCK phosphorylates the protein kinase-C (PKC), which ultimately induces smooth muscle cell contraction in the vascular system. Inhibition of ROCK1 has been shown to be a promising therapy for patients with cardiovascular disease. In this study, we have conducted molecular modeling techniques such as docking, molecular dynamics (MD), and 3-Dimensional structure-activity relationship (3D-QSAR) on a series of N-ethyl-4-(pyridin-4-yl)benzamide-based compounds. Docking and MD showed critical interactions and binding affinities between ROCK1 and its inhibitors. To establish the structure-activity relationship (SAR) of the compounds, 3D-QSAR techniques such as Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA) were used. The CoMFA (q ² = 0.774, r ² = 0.965, ONC = 6, and r p r e d 2 = 0.703) and CoMSIA (q ² = 0.676, r ² = 0.949, ONC = 6, and r p r e d 2 = 0.548) both models have shown reasonable external predictive activity, and contour maps revealed favorable and unfavorable substitutions for chemical group modifications. Based on the contour maps, we have designed forty new compounds, among which, seven compounds exhibited higher predictive activity (pIC₅₀). Further, we conducted the MD study, ADME/Tox, and SA score prediction using the seven newly designed compounds. The combination of docking, MD, and 3D-QSAR studies helps to understand the coherence modification of existing molecules. Our study may provide valuable insight into the development of more potent ROCK1 inhibitors.

Discovery of a phenylpyrazole amide ROCK inhibitor as a tool molecule for in vivo studies

Structure-activity relationship optimization on a series of phenylpyrazole amides led to the identification of a dual ROCK1 and ROCK2 inhibitor (25) which demonstrated good potency, kinome selectivity and favorable pharmacokinetic profiles. Compound 25 was selected as a tool molecule for in vivo studies including evaluating hemodynamic effects in telemeterized mice, from which moderate decreases in blood pressure were observed.

Computational systematic selectivity of the Fasalog inhibitors between ROCK-I and ROCK-II kinase isoforms in Alzheimer's disease

Human Rho-associated coiled-coil forming kinase (ROCK) is a class of essential neurokinases that consists of two structurally conserved isoforms ROCK-I and ROCK-II; they have been revealed to play distinct roles in the pathogenesis of Alzheimer's disease (AD) and other neurological disorders. Selective targeting of the two kinase isoforms with small-molecule inhibitors is a great challenge due to the surprisingly high homology in kinase domain (92 %) and the full identity in kinase active site (100 %). Here, we describe a computational protocol to systematically profile the selectivity of Fasudil and its 25 analogs (termed as Fasalogs) between the two kinase isoforms. It is suggested that the substitution of Fasudil's 1,4-diazepane moiety with rigid ring such as Ripasudil and Dimehtylfasudil would render the resulting inhibitors of ROCK-II over ROCK-I (II-o-I) selectivity, while the substitution with long, flexible group such as H-89 and BDBM92607 tends to have I-o-II selectivity. Structural analysis reveals that the inhibitor affinity is not only determined by the identical active site, but also contributed from the non-identical first and second shells of the site as well as other non-conserved kinase regions, which can indirectly influence the active site and inhibitor binding through allosteric effect. A further kinase assay basically confirms the computational findings, which also exhibits a good consistence with theoretical selectivity over 10 tested samples (Rp = 0.89). In particular, the Fasalog compounds Dimehtylfasudil and H-89 are identified as II-o-I and I-o-II selective inhibitors. They can be considered as promising lead molecular entities to develop new specific ROCK isoform-selective Fasalog inhibitors.

Computational study of paroxetine-like inhibitors reveals new molecular insight to inhibit GRK2 with selectivity over ROCK1

The G-protein coupled receptor kinase 2 (GRK2) regulates the desensitization of beta-adrenergic receptors (β-AR), and its overexpression has been implicated in heart failure. Hence, the inhibition of GRK2 is considered to be an important drug target for the treatment of heart failure. Due to the high sequence similarity of GRK2 with the A, G, and C family (AGC family) of kinases, the inhibition of GRK2 also leads to the inhibition of AGC kinases such as Rho-associated coiled-coil kinase 1 (ROCK1). Therefore, unraveling the mechanisms to selectively inhibit GRK2 poses an important challenge. We have performed molecular docking, three dimensional quantitative structure activity relationship (3D-QSAR), molecular dynamics (MD) simulation, and free energy calculations techniques on a series of 53 paroxetine-like compounds to understand the structural properties desirable for enhancing the inhibitory activity for GRK2 with selectivity over ROCK1. The formation of stable hydrogen bond interactions with the residues Phe202 and Lys220 of GRK2 seems to be important for selective inhibition of GRK2. Electropositive substituents at the piperidine ring and electronegative substituents near the amide linker between the benzene ring and pyrazole ring showed a higher inhibitory preference for GRK2 over ROCK1. This study may be used in designing more potent and selective GRK2 inhibitors for therapeutic intervention of heart failure.

Potently inhibiting cancer cell migration with novel 3H-pyrazolo[4,3-f]quinoline boronic acid ROCK inhibitors

Rho-associated protein kinases (ROCKs) are ubiquitously expressed in most adult tissues, and are involved in modulating the cytoskeleton, protein synthesis and degradation pathways, synaptic function, and autophagy to list a few. A few ROCK inhibitors, such as fasudil and netarsudil, are approved for clinical use. Here we present a new ROCK inhibitor, boronic acid containing HSD1590, which is more potent than netarsudil at binding to or inhibiting ROCK enzymatic activities. This compound exhibits single digit nanomolar binding to ROCK (Kds < 2 nM) and subnanomolar enzymatic inhibition profile (ROCK2 IC50 is 0.5 nM for HSD1590. Netarsudil, an FDA-approved drug, inhibited ROCK2 with IC50 = 11 nM under similar conditions). Whereas netarsudil was cytotoxic to breast cancer cell line, MDA-MB-231 (greater than 80% growth inhibition at concentrations greater than 5 μM), HSD1590 displayed low cytotoxicity to MDA-MB-231. Interestingly, at 1 μM HSD1590 inhibited the migration of MDA-MB-231 whereas netarsudil did not.

Accelerating Drug Discovery by Early Protein Drug Target Prediction Based on a Multi-Fingerprint Similarity Search

In this continuing work, we have updated our recently proposed Multi-fingerprint Similarity Search algorithm (MuSSel) by enabling the generation of dominant ionized species at a physiological pH and the exploration of a larger data domain, which included more than half a million high-quality small molecules extracted from the latest release of ChEMBL (version 24.1, at the time of writing). Provided with a high biological assay confidence score, these selected compounds explored up to 2822 protein drug targets. To improve the data accuracy, samples marked as prodrugs or with equivocal biological annotations were not considered. Notably, MuSSel performances were overall improved by using an object-relational database management system based on PostgreSQL. In order to challenge the real effectiveness of MuSSel in predicting relevant therapeutic drug targets, we analyzed a pool of 36 external bioactive compounds published in the Journal of Medicinal Chemistry from October to December 2018. This study demonstrates that the use of highly curated chemical and biological experimental data on one side, and a powerful multi-fingerprint search algorithm on the other, can be of the utmost importance in addressing the fate of newly conceived small molecules, by strongly reducing the attrition of early phases of drug discovery programs.