Phosphorylation-Induced Conformational Dynamics and Inhibition of Janus Kinase 1 by Suppressors of Cytokine Signaling 1

Dynamic Interplay of Loop Motions Governs the Molecular Level Regulatory Dynamics in Spleen Tyrosine Kinase: Insights from Molecular Dynamics Simulations

The spleen tyrosine kinase (Syk) is a key regulator in immune cell signaling and is linked to various mechanisms in cancer and neurodegenerative diseases. Although most computational research on Syk focuses on novel drug design, the molecular-level regulatory dynamics remain unexplored. In this study, we utilized 5 × 1 μs all-atom molecular dynamics simulations of the Syk kinase domain, examining it in combinations of activation segment phosphorylated/unphosphorylated (at Tyr525, Tyr526) and the "DFG"-Asp protonated/deprotonated (at Asp512) states to investigate conformational variations and regulatory dynamics of various loops and motifs within the kinase domain. Our findings revealed that the formation and disruption of several electrostatic interactions among residues within and near the activation segment likely influenced its dynamics. The protein structure network analysis indicated that the N-terminal and C-terminal anchors were stabilized by connections with the nearby stable helical regions. The P-loop showed conformational variation characterized by movements toward and away from the conserved "HRD"-motif. Additionally, there was a significant correlation between the movement of the β3-αC loop and the P-loop, which controls the dimensions of the adenine-binding cavity of the C-spine region. Overall, understanding these significant motions of the Syk kinase domain enhances our knowledge of its functional regulatory mechanism and can guide future research.

Molecular dynamics simulations reveal phosphorylation-induced conformational dynamics of the fibroblast growth factor receptor 1 kinase

The Fibroblast Growth Factor Receptor1 (FGFR1) kinase wields exquisite control on cell fate, proliferation, differentiation, and homeostasis. An imbalance of FGFR1 signaling leads to several pathogeneses of diseases ranging from multiple cancers to allergic and neurodegenerative disorders. In this study, we investigated the phosphorylation-induced conformational dynamics of FGFR1 in apo and ATP-bound states via all-atom molecular dynamics simulations. All simulations were performed for 2 × 2 µs. We have also investigated the energetics of the binding of ATP to FGFR1 using the molecular mechanics Poisson-Boltzmann scheme. Our study reveals that the FGFR1 kinase can reach a fully active configuration through phosphorylation and ATP binding. A 3-10 helix formation in the activation loop signifies its rearrangement leading to stability upon ATP binding. The interaction of phosphorylated tyrosine (pTyr654) with positively charged residues forms strong salt-bridge interactions, driving the compactness of the structure. The dynamic cross-correlation map reveals phosphorylation enhances correlated motions and reduces anti-correlated motions between different domains. We believe that the mechanistic understanding of large-conformational changes upon the activation of the FGFR1 kinase will aid the development of novel targeted therapeutics.Communicated by Ramaswamy H. Sarma.

Microsecond dynamics of H10N7 influenza neuraminidase reveals the plasticity of loop regions and drug resistance due to the R292K mutation

At the beginning of the last century, multiple pandemics caused by influenza (flu) viruses severely impacted public health. Despite the development of vaccinations and antiviral medications to prevent and control impending flu outbreaks, unforeseen novel strains and continuously evolving old strains continue to represent a serious threat to human life. Therefore, the recently identified H10N7, for which not much data is available for rational structure-based drug design, needs to be further explored. Here, we investigated the structural dynamics of neuraminidase N7 upon binding of inhibitors, and the drug resistance mechanisms against the oseltamivir (OTV) and laninamivir (LNV) antivirals due to the crucial R292K mutation on the N7 using the computational microscope, molecular dynamics (MD) simulations. In this study, each system underwent long 2 × 1 μs MD simulations to answer the conformational changes and drug resistance mechanisms. These long time-scale dynamics simulations and free energy landscapes demonstrated that the mutant systems showed a high degree of conformational variation compared to their wildtype (WT) counterparts, and the LNV-bound mutant exhibited an extended 150-loop conformation. Further, the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculation and MM/GBSA free energy decomposition were used to characterize the binding of OTV and LNV with WT, and R292K mutated N7, revealing the R292K mutation as drug-resistant, facilitated by a decline in binding interaction and a reduction in the dehydration penalty. Due to the broader binding pocket cavity of the smaller K292 mutant residue relative to the wildtype, the drug carboxylate to K292 hydrogen bonding was lost, and the area surrounding the K292 residue was more accessible to water molecules. This implies that drug resistance could be reduced by strengthening the hydrogen bond contacts between N7 inhibitors and altered N7, creating inhibitors that can form a hydrogen bond to the mutant K292, or preserving the closed cavity conformations.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Development of chalcone-like derivatives and their biological and mechanistic investigations as novel influenza nuclear export inhibitors

Concerning the emergence of resistance to current anti-influenza drugs, our previous phenotypic-based screening study identified the compound A9 as a promising lead compound. This chalcone analog, containing a 2,6-dimethoxyphenyl moiety, exhibited significant inhibitory activity against oseltamivir-resistant strains (H1N1 pdm09), with an EC50 value of 1.34 μM. However, it also displayed notable cytotoxicity, with a CC50 value of 41.46 μM. Therefore, compound A9 was selected as a prototype structure for further structural optimization in this study. Initially, it was confirmed that the substituting the α,β-unsaturated ketone with pent-1,4-diene-3-one as a linker group significantly reduced the cytotoxicity of the final compounds. Subsequently, the penta-1,4-dien-3-one group was utilized as a privileged fragment for further structural optimization. Following two subsequent rounds of optimizations, we identified compound IIB-2, which contains a 2,6-dimethoxyphenyl- and 1,4-pentadiene-3-one moieties. This compound exhibited inhibitory effects on oseltamivir-resistant strains comparable to its precursor (compound A9), while demonstrating reduced toxicity (CC50 > 100 μM). Furthermore, we investigated its mechanism of action against anti-influenza virus through immunofluorescence, Western blot, and surface plasmon resonance (SPR) experiments. The results revealed that compound IIB-2 can impede virus proliferation by blocking the export of influenza virus nucleoprotein. Thusly, our findings further emphasize influenza nuclear export as a viable target for designing novel chalcone-like derivatives with potential inhibitory properties that could be explored in future lead optimization studies.