Abstract: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.

Residue-specific free energy analysis in ligand bindings to JAK2

Computational Alanine Scanning with Interaction Entropy for Protein-ligand Binding Free Energies.

Computational Analysis for Residue-Specific CDK2-Inhibitor Bindings

Theoretical Exploring Selective-Binding Mechanisms of JAK3 by 3D-QSAR, Molecular Dynamics Simulation and Free Energy Calculation

Interaction Entropy for Computational Alanine Scanning in Protein–protein Binding

Interaction Entropy for Computational Alanine Scanning.

Amino Acid Residue E543 in Jak2 C618r is A Potential Therapeutic Target for Myeloproliferative Disorders Caused by Jak2 C618r Mutation

The Effects of R683S (G) Genetic Mutations on the JAK2 Activity, Structure and Stability.

Computational analyses of JAK1 kinase domain: subtle changes in the catalytic cleft influence inhibitor specificity.

A High-Throughput Percentage-of-binding Strategy to Measure Binding Energies in DNA-protein Interactions: Application to Genome-Scale Site Discovery.

Calculation of Relative Binding Free Energy for Mutations in Protein Complexes: The Alchemical Path

Binding Free Energy Estimation for Protein-Ligand Complex Based on MM-PBSA with Various Partial Charge Models.

Computational Analysis of Residue-Specific Binding Free Energies of Androgen Receptor to Ligands

Molecular Mechanism of Ligand Bindings to Zika Virus at SAM Site

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

Accurate and Efficient Calculation of Protein-Protein Binding Free Energy-Interaction Entropy with Residue Type-Specific Dielectric Constants.

A Computational Analysis of the Binding Affinities of FKBP12 Inhibitors Using the MM‐PB/SA Method

On achieving high accuracy and reliability in the calculation of relative protein–ligand binding affinities

Calculation of PARP1-ligand Binding Free Energy Using Alanine Scanning and Interaction Entropy Method

Predictions Of Binding Of A Diverse Set Of Ligands To Gelatinase-A By A Combination Of Molecular Dynamics And Continuum Solvent Models

Comparison of Methodologies for Absolute Binding Free Energy Calculations of Ligands to Intrinsically Disordered Proteins