Quantitative assessment of kinase selectivity based the water-ring network in protein binding sites using molecular dynamics simulations

Unraveling Extremely Damaging IRAK4 Variants and Their Potential Implications for IRAK4 Inhibitor Efficacy

Interleukin-1-receptor-associated kinase 4 (IRAK4) possesses a crucial function in the toll-like receptor (TLR) signaling pathway, and the dysfunction of this molecule could lead to various infectious and immune-related diseases in addition to cancers. IRAK4 genetic variants have been linked to various types of diseases. Therefore, we conducted a comprehensive analysis to recognize the missense variants with the most damaging impacts on IRAK4 with the employment of diverse bioinformatics tools to study single-nucleotide polymorphisms' effects on function, stability, secondary structures, and 3D structure. The residues' location on the protein domain and their conservation status were investigated as well. Moreover, docking tools along with structural biology were engaged in analyzing the SNPs' effects on one of the developed IRAK4 inhibitors. By analyzing IRAK4 gene SNPs, the analysis distinguished ten variants as the most detrimental missense variants. All variants were situated in highly conserved positions on an important protein domain. L318S and L318F mutations were linked to changes in IRAK4 secondary structures. Eight SNPs were revealed to have a decreasing effect on the stability of IRAK4 via both I-Mutant 2.0 and Mu-Pro tools, while Mu-Pro tool identified a decreasing effect for the G198E SNP. In addition, detrimental effects on the 3D structure of IRAK4 were also discovered for the selected variants. Molecular modeling studies highlighted the detrimental impact of these identified SNP mutant residues on the druggability of the IRAK4 ATP-binding site towards the known target inhibitor, HG-12-6, as compared to the native protein. The loss of important ligand residue-wise contacts, altered protein global flexibility, increased steric clashes, and even electronic penalties at the ligand-binding site interfaces were all suggested to be associated with SNP models for hampering the HG-12-6 affinity towards IRAK4 target protein. This given model lays the foundation for the better prediction of various disorders relevant to IRAK4 malfunction and sheds light on the impact of deleterious IRAK4 variants on IRAK4 inhibitor efficacy.

TWN-FS method: A novel fragment screening method for drug discovery

Fragment-based drug discovery (FBDD) is a well-established and effective method for generating diverse and novel hits in drug design. Kinases are suitable targets for FBDD due to their well-defined structure. Water molecules contribute to structure and function of proteins and also influence the environment within the binding pocket. Water molecules form a variety of hydrogen-bonded cyclic water-ring networks, collectively known as topological water networks (TWNs). Analyzing the TWNs in protein binding sites can provide valuable insights into potential locations and shapes for fragments within the binding site. Here, we introduce TWN-based fragment screening (TWN-FS) method, a novel screening method that suggests fragments through grouped TWN analysis within the protein binding site. We used this method to screen known CDK2, CHK1, IGF1R and ERBB4 inhibitors. Our findings suggest that TWN-FS method has the potential to effectively screen fragments. The TWN-FS method package is available on GitHub at https://github.com/pkj0421/TWN-FS.

TWN-RENCOD: A novel method for protein binding site comparison

Several diverse proteins possess similar binding sites. Protein binding site comparison provides valuable insights for the drug discovery and development. Binding site similarities are useful in understanding polypharmacology, identifying potential off-targets and repurposing of known drugs. Many binding site analysis and comparison methods are available today, however, these methods may not be adequate to explain variation in the activity of a drug or a small molecule against a number of similar proteins. Water molecules surrounding the protein surface contribute to structure and function of proteins. Water molecules form diverse types of hydrogen-bonded cyclic water-ring networks known as topological water networks (TWNs). Analysis of TWNs in binding site of proteins may improve understanding of the characteristics of binding sites. We propose TWN-based residue encoding (TWN-RENCOD), a novel binding site comparison method which compares the aqueous environment in binding sites of similar proteins. As compared to other existing methods, results obtained using our method correlated better with differences in wide range of activity of a known drug (Sunitinib) against nine different protein kinases (KIT, PDGFRA, VEGFR2, PHKG2, ITK, HPK1, MST3, PAK6 and CDK2).

Small Molecule Inhibitors of DYRK1A Identified by Computational and Experimental Approaches

Dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) is a protein kinase with diverse functions in cell regulation. Abnormal expression and activity of DYRK1A contribute to numerous human malignancies, Down syndrome, and Alzheimer's disease. Notably, DYRK1A has been proposed as a potential therapeutic target for the treatment of diabetes because of its key role in pancreatic β-cell proliferation. Consequently, DYRK1A is an attractive drug target for a variety of diseases. Here, we report the identification of several DYRK1A inhibitors using our in-house topological water network-based approach. All inhibitors were further verified by in vitro assay.

Pharmacophoric Site Identification and Inhibitor Design for Autotaxin

Autotaxin (ATX) is a potential drug target that is associated with inflammatory diseases and various cancers. In our previous studies, we have designed several inhibitors targeting ATX using computational and experimental approaches. Here, we have analyzed topological water networks (TWNs) in the binding pocket of ATX. TWN analysis revealed a pharmacophoric site inside the pocket. We designed and synthesized compounds considering the identified pharmacophoric site. Furthermore, we performed biological experiments to determine their ATX inhibitory activities. High potency of the designed compounds supports the predictions of the TWN analysis.

Topological Water Network Analysis Around Amino Acids

Water molecules play a key role in protein stability, folding, function and ligand binding. Protein hydration has been studied using free energy perturbation algorithms. However, the study of protein hydration without free energy calculation is also an active field of research. Accordingly, topological water network (TWN) analysis has been carried out instead of free energy calculation in the present work to investigate hydration of proteins. Water networks around 20 amino acids in the aqueous solution were explored through molecular dynamics (MD) simulations. These simulation results were compared with experimental observations. Water molecules from the protein data bank structures showed TWN patterns similar to MD simulations. This work revealed that TWNs are effected by the surrounding environment. TWNs could provide valuable clues about the environment around amino acid residues in the proteins. The findings from this study could be exploited for TWN-based drug discovery and development.

Computational insights into the role of α-strand/sheet in aggregation of α-synuclein

The α-synuclein is a major component of amyloid fibrils found in Lewy bodies, the characteristic intracellular proteinaceous deposits which are pathological hallmarks of neurodegenerative diseases such as Parkinson’s disease (PD) and dementia. It is an intrinsically disordered protein that may undergo dramatic structural changes to form amyloid fibrils. Aggregation process from α-synuclein monomers to amyloid fibrils through oligomeric intermediates is considered as the disease-causative toxic mechanism. However, mechanism underlying aggregation is not well-known despite several attempts. To characterize the mechanism, we have explored the effects of pH and temperature on the structural properties of wild-type and mutant α-synuclein using molecular dynamics (MD) simulation technique. MD studies suggested that amyloid fibrils can grow by monomer. Conformational transformation of the natively unfolded protein into partially folded intermediate could be accountable for aggregation and fibrillation. An intermediate α-strand was observed in the hydrophobic non-amyloid-β component (NAC) region of α-synuclein that could proceed to α-sheet and initiate early assembly events. Water network around the intermediate was analyzed to determine its influence on the α-strand structure. Findings of this study provide novel insights into possible mechanism of α-synuclein aggregation and promising neuroprotective strategy that could aid alleviate PD and its symptoms.

Design of a Novel and Selective IRAK4 Inhibitor Using Topological Water Network Analysis and Molecular Modeling Approaches

Protein kinases are deeply involved in immune-related diseases and various cancers. They are a potential target for structure-based drug discovery, since the general structure and characteristics of kinase domains are relatively well-known. However, the ATP binding sites in protein kinases, which serve as target sites, are highly conserved, and thus it is difficult to develop selective kinase inhibitors. To resolve this problem, we performed molecular dynamics simulations on 26 kinases in the aqueous solution, and analyzed topological water networks (TWNs) in their ATP binding sites. Repositioning of a known kinase inhibitor in the ATP binding sites of kinases that exhibited a TWN similar to interleukin-1 receptor-associated kinase 4 (IRAK4) allowed us to identify a hit molecule. Another hit molecule was obtained from a commercial chemical library using pharmacophore-based virtual screening and molecular docking approaches. Pharmacophoric features of the hit molecules were hybridized to design a novel compound that inhibited IRAK4 at low nanomolar levels in the in vitro assay.