Computational analysis of regulatory regions in human protein kinases

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Genome-wide characterization and evolutionary analysis of linker histones in castor bean (Ricinus communis)

H1s, or linker histones, are ubiquitous proteins in eukaryotic cells, consisting of a globular GH1 domain flanked by two unstructured tails. Whilst it is known that numerous non-allelic variants exist within the same species, the degree of interspecific and intraspecific variation and divergence of linker histones remain unknown. The conserved basic binding sites in GH1 and evenly distributed strong positive charges on the C-terminal domain (CTD) are key structural characters for linker histones to bind chromatin. Based on these features, we identified five linker histones from 13 GH1-containing proteins in castor bean (Ricinus communis), which were named as RcH1.1, RcH1.2a, RcH1.2b, RcH1.3, and RcH1.4 based on their phylogenetic relationships with the H1s from five other economically important Euphorbiaceae species (Hevea brasiliensis Jatropha curcas, Manihot esculenta Mercurialis annua, and Vernicia fordii) and Arabidopsis thaliana. The expression profiles of RcH1 genes in a variety of tissues and stresses were determined from RNA-seq data. We found three RcH1 genes (RcH1.1, RcH1.2a, and RcH1.3) were broadly expressed in all tissues, suggesting a conserved role in stabilizing and organizing the nuclear DNA. RcH1.2a and RcH1.4 was preferentially expressed in floral tissues, indicating potential involvement in floral development in castor bean. Lack of non-coding region and no expression detected in any tissue tested suggest that RcH1.2b is a pseudogene. RcH1.3 was salt stress inducible, but not induced by cold, heat and drought in our investigation. Structural comparison confirmed that GH1 domain was highly evolutionarily conserved and revealed that N- and C-terminal domains of linker histones are divergent between variants, but highly conserved between species for a given variant. Although the number of H1 genes varies between species, the number of H1 variants is relatively conserved in more closely related species (such as within the same family). Through comparison of nucleotide diversity of linker histone genes and oil-related genes, we found similar mutation rate of these two groups of genes. Using Tajima's D and ML-HKA tests, we found RcH1.1 and RcH1.3 may be under balancing selection.

New insights into structural determinants of prion protein folding and stability

Prions are the etiological agent of fatal neurodegenerative diseases called prion diseases or transmissible spongiform encephalopathies. These maladies can be sporadic, genetic or infectious disorders. Prions are due to post-translational modifications of the cellular prion protein leading to the formation of a β-sheet enriched conformer with altered biochemical properties. The molecular events causing prion formation in sporadic prion diseases are still elusive. Recently, we published a research elucidating the contribution of major structural determinants and environmental factors in prion protein folding and stability. Our study highlighted the crucial role of octarepeats in stabilizing prion protein; the presence of a highly enthalpically stable intermediate state in prion-susceptible species; and the role of disulfide bridge in preserving native fold thus avoiding the misfolding to a β-sheet enriched isoform. Taking advantage from these findings, in this work we present new insights into structural determinants of prion protein folding and stability.