To date, many experiments have revealed that the functional balance between hemagglutinin (HA) and neuraminidase (NA) plays a crucial role in viral mobility, production, and transmission. However, whether and how HA and NA maintain balance at the sequence level needs further investigation. Here, we applied principal component analysis and hierarchical clustering analysis on thousands of HA and NA sequences of A/H1N1 and A/H3N2. We discovered significant coevolution between HA and NA at the sequence level, which is closely related to the type of host species and virus epidemic years. Furthermore, we propose a sequence-to-sequence transformer model (S2STM), which mainly consists of an encoder and a decoder that adopts a multi-head attention mechanism for establishing the mapping relationship between HA and NA sequences. The training results reveal that the S2STM can effectively realize the “translation” from HA to NA or vice versa, thereby building a relationship network between them. Our work combines unsupervised and supervised machine learning methods to identify the sequence matching between HA and NA, which will advance our understanding of IAVs’ evolution and also provide a novel idea for sequence analysis methods.

p38α is a key serine/threonine kinase that can enable atypical auto-activation through Zap70 phosphorylation and initiate T cell receptor signaling. The auto-activation plays an important role in autoimmune diseases. Although the classical activation mechanism of p38α has been studied in-depth, the atypical activation mechanism of Y323 phosphorylation-induced p38α auto-activation remains largely unexplained, especially the regulatory effects of phosphorylation on different sites (Y323 vs T180). From the X-ray experimental data, we identified the inactive and active states of p38α using principal component analysis. To understand the auto-activation process and the internal driving mechanism, a computational paradigm that couples the targeted molecular dynamics simulations, the String Method, and the umbrella sampling strategy were employed to generate the conformational landscape of p38α, including p38α T180–Y323, p38α T180–pY323, and p38α pT180–pY323 systems (pT180/pY323: phosphorylated T180/Y323). We explored that pY323 could change the conformational distribution and promote the conformational transition of p38α from the inactive state to the active state. Auto-activation of p38α is regulated by pY323 through destabilization of the hydrophobic core structure and aided by R173. This study will further explain the conformational transition of p38α induced by Y323 phosphorylation and provide insights into the universal molecular auto-activation mechanism of the p38 subfamily at the atomic level.

Lymphokine-activated killer T-cell-originated protein kinase (TOPK) is a potential target for cancer therapy. To explore the micromechanism, we proposed the N-terminal premodel (NTPM) of the TOPK monomer via homology modeling and molecular dynamic simulations and analyzed the conformational dynamics by Markov state model analysis. The electronegative insert (ENI) motif of the NTPM can be opened with a small probability under wild type, regulated by the so-called “N–C” interaction zone consisting of the N-terminal head, the coil between β3-strand and αC-helix, and the ENI motif. Glutamate substitution at threonine residue 9 or tyrosine residue 74 promotes the closed–open transition, revealing the details of phosphorylation. Allosteric effects induce functionally relevant structural changes, such as increased structural flexibility and active sites, which are thought to be necessary for further activation or binding. These findings provide rational structural templates for designing state-dependent inhibitors and give insight into the molecular regulatory mechanisms of TOPK monomers.