This review outlines the role of electrostatics in computational molecular biophysics and its implication in altering wild-type characteristics of biological macromolecules, and thus the contribution of electrostatics to disease mechanisms. The work is not intended to review existing computational approaches or to propose further developments. Instead, it summarizes the outcomes of relevant studies and provides a generalized classification of major mechanisms that involve electrostatic effects in both wild-type and mutant biological macromolecules. It emphasizes the complex role of electrostatics in molecular biophysics, such that the long range of electrostatic interactions causes them to dominate all other forces at distances larger than several Angstroms, while at the same time, the alteration of short-range wild-type electrostatic pairwise interactions can have pronounced effects as well. Because of this dual nature of electrostatic interactions, being dominant at long-range and being very specific at short-range, their implications for wild-type structure and function are quite pronounced. Therefore, any disruption of the complex electrostatic network of interactions may abolish wild-type functionality and could be the dominant factor contributing to pathogenicity. However, we also outline that due to the plasticity of biological macromolecules, the effect of amino acid mutation may be reduced, and thus a charge deletion or insertion may not necessarily be deleterious.

Microtubules are key players in several stages of the cell cycle and are also involved in the transportation of cellular organelles. Microtubules are polymerized by α/β tubulin dimers with a highly dynamic feature, especially at the plus ends of the microtubules. Therefore, understanding the interactions among tubulins is crucial for characterizing microtubule dynamics. Studying microtubule dynamics can help researchers make advances in the treatment of neurodegenerative diseases and cancer. In this study, we utilize a series of computational approaches to study the electrostatic interactions at the binding interfaces of tubulin monomers. Our study revealed that among all the four types of tubulin-tubulin binding modes, the electrostatic attractive interactions in the α/β tubulin binding are the strongest while the interactions of α/α tubulin binding in the longitudinal direction are the weakest. Our calculations explained that due to the electrostatic interactions, the tubulins always preferred to form α/β tubulin dimers. The interactions between two protofilaments are the weakest. Thus, the protofilaments are easily separated from each other. Furthermore, the important residues involved in the salt bridges at the binding interfaces of the tubulins are identified, which illustrates the details of the interactions in the microtubule. This study elucidates some mechanistic details of microtubule dynamics and also identifies important residues at the binding interfaces as potential drug targets for the inhibition of cancer cells.

Janus tyrosine kinase 3 (JAK3) is primarily expressed in immune cells and is needed for signaling by the common gamma chain (γc) family of cytokines. Abnormal JAK3 signal transduction can manifest as hematological disorders, e.g., leukemia, severe combined immunodeficiency (SCID) and autoimmune disease states. While regulatory JAK3 phosphosites have been well studied, here a functional proteomics approach coupling a JAK3 autokinase assay to mass spectrometry revealed ten previously unreported autophosphorylation sites (Y105, Y190, Y238, Y399, Y633, Y637, Y738, Y762, Y824, and Y841). Of interest, Y841 was determined to be evolutionarily conserved across multiple species and JAK family members, suggesting a broader role for this residue. Phospho-substitution mutants confirmed that Y841 is also required for STAT5 tyrosine phosphorylation. The homologous JAK1 residue Y894 elicited a similar response to mutagenesis, indicating the shared importance for this site in JAK family members. Phospho-specific Y841-JAK3 antibodies recognized activated kinase from various T-cell lines and transforming JAK3 mutants. Computational biophysics analysis linked Y841 phosphorylation to enhanced JAK3 JH1 domain stability across pH environments, as well as to facilitated complementary electrostatic JH1 dimer formation. Interestingly, Y841 is not limited to tyrosine kinases, suggesting it represents a conserved ubiquitous enzymatic function that may hold therapeutic potential across multiple kinase families.

Overactive Janus kinases (JAKs) are known to drive leukemia, making them well-suited targets for treatment. We sought to identify new JAK-activating mutations and instead found a JAK1-inactivating pseudokinase mutation, V666G. In contrast to other pseudokinase mutations that canonically lead to an active kinase, the JAK1 V666G mutation led to under-activation seen by reduced phosphorylation. To understand the functional role of JAK1 V666G in modifying kinase activity we investigated its influence on other JAK kinases and within the Interleukin-2 pathway. JAK1 V666G not only inhibited its own activity, but its presence could inhibit other JAK kinases. These findings provide new insights into the potential of JAK1 pseudokinase to modulate its own activity, as well as of other JAK kinases. Thus, the features of the JAK1 V666 region in modifying JAK kinases can be exploited to allosterically inhibit overactive JAKs.

Janus tyrosine kinase (JAK) variants are known drivers for hematological disorders. With the full-length structure of mouse JAK1 being recently resolved, new observations on the localization of variants within closed, open, and dimerized JAK structures are possible. Full-length homology models of human wild-type JAK family members were developed using the Glassman et al. reported mouse JAK1 containing the V658F structure as a template. Many mutational sites related to proliferative hematological disorders reside in the JH2 pseudokinase domains facing the region important in dimerization of JAKs in both closed and open states. More than half of all JAK gain of function (GoF) variants are changes in polarity, while only 1.2% are associated with a change in charge. Within a JAK1-JAK3 homodimer model, IFNLR1 (PDB ID7T6F) and the IL-2 common gamma chain subunit (IL2Rγc) were aligned with the respective dimer implementing SWISS-MODEL coupled with ChimeraX. JAK3 variants were observed to encircle the catalytic site of the kinase domain, while mutations in the pseudokinase domain align along the JAK-JAK dimerization axis. FERM domains of JAK1 and JAK3 are identified as a hot spot for hematologic malignancies. Herein, we propose new allosteric surfaces for targeting hyperactive JAK dimers.

Janus Kinase 3 (JAK3) plays a key role in the development, proliferation, and differentiation of various immune cells. It regulates gene expression by phosphorylation of Signal Transducers and Activators of Transcriptions (STATs) via the JAK/STAT pathway. Recently, we found a new JAK3 phosphorylation site, tyrosine 841 (Y841). The results showed that pY841 helps the kinase domain flip around the pseudo kinase domain, which may cause JAK3 conformational changes. It also reduces the size of the cleft between the N-lobe and the C-lobe of the JAK3 kinase domain. However, pY841 was found to enlarge the cleft when ATP/ADP was bound to the kinase. The increase in the cleft size suggested that pY841 enhanced the elasticity of the kinase domain. For unphosphorylated JAK3 (JAK3-Y841), the binding forces between the kinase domain and ATP or ADP were similar. After phosphorylation of Y841, JAK3-pY841 exhibited more salt bridges and hydrogen bonds between ATP and the kinase than between ADP and the kinase. Consequently, the electrostatic binding force between ATP and the kinase was higher than that between ADP and the kinase. The result was that compared to ADP, ATP was more attractive to JAK3 when Y841 was phosphorylated. Therefore, JAK3-pY841 tended to bind ATP rather than ADP. This work provides new insights into the role of phosphorylation in kinase activation and ATP hydrolysis and sheds light on the importance of understanding the molecular mechanisms that regulate the kinase function.