PyTMs: a useful PyMOL plugin for modeling common post-translational modifications

Hypoxia increases persulfide and polysulfide formation by AMP kinase dependent cystathionine gamma lyase phosphorylation

Hydropersulfide and hydropolysulfide metabolites are increasingly important reactive sulfur species (RSS) regulating numerous cellular redox dependent functions. Intracellular production of these species is known to occur through RSS interactions or through translational mechanisms involving cysteinyl t-RNA synthetases. However, regulation of these species under cell stress conditions, such as hypoxia, that are known to modulate RSS remain poorly understood. Here we define an important mechanism of increased persulfide and polysulfide production involving cystathionine gamma lyase (CSE) phosphorylation at serine 346 and threonine 355 in a substrate specific manner, under acute hypoxic conditions. Hypoxic phosphorylation of CSE occurs in an AMP kinase dependent manner increasing enzyme activity involving unique inter- and intramolecular interactions within the tetramer. Importantly, both cellular hypoxia and tissue ischemia result in AMP Kinase dependent CSE phosphorylation that regulates blood flow in ischemic tissues. Our findings reveal hypoxia molecular signaling pathways regulating CSE dependent persulfide and polysulfide production impacting tissue and cellular response to stress.

An allosteric switch between the activation loop and a c-terminal palindromic phospho-motif controls c-Src function

Autophosphorylation controls the transition between discrete functional and conformational states in protein kinases, yet the structural and molecular determinants underlying this fundamental process remain unclear. Here we show that c-terminal Tyr 530 is a de facto c-Src autophosphorylation site with slow time-resolution kinetics and a strong intermolecular component. On the contrary, activation-loop Tyr 419 undergoes faster kinetics and a cis-to-trans phosphorylation switch that controls c-terminal Tyr 530 autophosphorylation, enzyme specificity, and strikingly, c-Src non-catalytic function as a substrate. In line with this, we visualize by X-ray crystallography a snapshot of Tyr 530 intermolecular autophosphorylation. In an asymmetric arrangement of both catalytic domains, a c-terminal palindromic phospho-motif flanking Tyr 530 on the substrate molecule engages the G-loop of the active kinase adopting a position ready for entry into the catalytic cleft. Perturbation of the phospho-motif accounts for c-Src dysfunction as indicated by viral and colorectal cancer (CRC)-associated c-terminal deleted variants. We show that c-terminal residues 531 to 536 are required for c-Src Tyr 530 autophosphorylation, and such a detrimental effect is caused by the substrate molecule inhibiting allosterically the active kinase. Our work reveals a crosstalk between the activation and c-terminal segments that control the allosteric interplay between substrate- and enzyme-acting kinases during autophosphorylation.

Histone H3 E50K mutation confers oncogenic activity and supports an EMT phenotype

Sequencing of human patient tumors has identified recurrent missense mutations in genes encoding core histones. We report that mutations that convert histone H3 amino acid 50 from a glutamate to a lysine (H3E50K) support an oncogenic phenotype in human cells. Expression of H3E50K is sufficient to transform human cells as evidenced by a dramatic increase in cell migration and invasion, and a statistically significant increase in proliferation and clonogenicity. H3E50K also increases the invasive phenotype in the context of co-occurring BRAF mutations, which are present in patient tumors characterized by H3E50K. H3E50 lies on the globular domain surface in a region that contacts H4 within the nucleosome. We find that H3E50K perturbs proximal H3 post-translational modifications globally and dysregulates gene expression, activating the epithelial to mesenchymal transition. Functional studies using S. cerevisiae reveal that, while yeast cells that express H3E50K as the sole copy of histone H3 show sensitivity to cellular stressors, including caffeine, H3E50K cells display some genetic interactions that are distinct from the characterized H3K36M oncohistone yeast model. Taken together, these data suggest that additional histone H3 mutations have the potential to be oncogenic drivers and function through distinct mechanisms that dysregulate gene expression.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

PyProtif: a PyMol plugin to retrieve and visualize protein motifs for structural studies

Proteins often possess several motifs and the ones with similar motifs were found to have similar biochemical properties and thus related biological functions. Thereby, multiple databases were developed to store information on such motifs in proteins. For instance, PDBsum stores the results of Promotif's generated structural motifs and Pfam stores pre-computed patterns of functional domains. In addition to the fact that all this stored information is extremely useful, we can further augment its importance if we ought to integrate these motifs into visualization software. In this work, we have developed PyProtif, a plugin for the PyMOL molecular visualization program, which automatically retrieves protein structural and functional motifs from different databases and integrates them in PyMOL for visualization and analyses. Through an expendable menu and a user-friendly interface, the plugin grants the users the ability to study simultaneously multiple proteins and to select and manipulate each motif separately. Thus, this plugin will be of great interest for structural, evolutionary and classification studies of proteins.

Disease-associated patterns of acetylation stabilize tau fibril formation

Assembly of tau into beta-sheet-rich amyloids dictates the pathology of a diversity of diseases. Lysine acetylation has been proposed to drive tau amyloid assembly, but no direct mechanism has emerged. Using tau fragments, we identify patterns of acetylation that flank amyloidogenic motifs on the tau fragments that promote rapid fibril assembly. We determined a 3.9 Å cryo-EM amyloid fibril structure assembled from an acetylated tau fragment uncovering how lysine acetylation can mediate gain-of-function interactions. Comparison of the structure to an ex vivo tauopathy fibril reveals regions of structural similarity. Finally, we show that fibrils encoding disease-associated patterns of acetylation are active in cell-based tau aggregation assays. Our data uncover the dual role of lysine residues in limiting tau aggregation while their acetylation leads to stabilizing pro-aggregation interactions. Design of tau sequence with specific acetylation patterns may lead to controllable tau aggregation to direct folding of tau into distinct amyloid folds.

14-3-3 phosphorylation inhibits 14-3-3θ's ability to regulate LRRK2 kinase activity

LRRK2 mutations are among the most common genetic causes for Parkinson's disease (PD), and toxicity is associated with increased kinase activity. 14-3-3 proteins are key interactors that regulate LRRK2 kinase activity. Phosphorylation of the 14-3-3θ isoform at S232 is dramatically increased in human PD brains. Here we investigate the impact of 14-3-3θ phosphorylation on its ability to regulate LRRK2 kinase activity. Both wildtype and the non-phosphorylatable S232A 14-3-3θ mutant reduced the kinase activity of wildtype and G2019S LRRK2, whereas the phosphomimetic S232D 14-3-3θ mutant had minimal effects on LRRK2 kinase activity, as determined by measuring autophosphorylation at S1292 and T1503 and Rab10 phosphorylation. However, wildtype and both 14-3-3θ mutants similarly reduced the kinase activity of the R1441G LRRK2 mutant. 14-3-3θ phosphorylation did not promote global dissociation with LRRK2, as determined by co-immunoprecipitation and proximal ligation assays. 14-3-3s interact with LRRK2 at several phosphorylated serine/threonine sites, including T2524 in the C-terminal helix, which can fold back to regulate the kinase domain. Interaction between 14-3-3θ and phosphorylated T2524 LRRK2 was important for 14-3-3θ's ability to regulate kinase activity, as wildtype and S232A 14-3-3θ failed to reduce the kinase activity of G2019S/T2524A LRRK2. Molecular modeling showed that 14-3-3θ phosphorylation causes a partial rearrangement of its canonical binding pocket, thus affecting the interaction between 14-3-3θ and the C-terminus of LRRK2. We conclude that 14-3-3θ phosphorylation destabilizes the interaction of 14-3-3θ with LRRK2 at T2524, which consequently promotes LRRK2 kinase activity.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as ₂₁₆PTPP₂₁₉, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

nNOS-induced tyrosine nitration of TRKB impairs BDNF signaling and restrains neuronal plasticity

Nitric oxide (NO) has been long recognized as an important modulator of neural plasticity, but characterization of the molecular mechanisms involved - specially the guanylyl cyclase-independent ones - has been challenging. There is evidence that NO could modify BDNF-TRKB signaling, a key mediator of neuronal plasticity. However, the mechanism underlying the interplay of NO and TRKB remains unclear. Here we show that NO induces nitration of the tyrosine 816 in the TRKB receptor in vivo and in vitro, and that post-translational modification inhibits TRKB phosphorylation and binding of phospholipase Cγ1 (PLCγ1) to this same tyrosine residue. Additionally, nitration triggers clathrin-dependent endocytosis of TRKB through the adaptor protein AP-2 and ubiquitination, thereby increasing translocation of TRKB away from the neuronal surface and directing it towards lysosomal degradation. Accordingly, inhibition of nitric oxide increases TRKB phosphorylation and TRKB-dependent neurite branching in neuronal cultures. In vivo, chronic inhibition of neuronal nitric oxide synthase (nNOS) dramatically reduced TRKB nitration and facilitated TRKB signaling in the visual cortex, and promoted a shift in ocular dominance upon monocular deprivation - an indicator of increased plasticity. Altogether, our data describe and characterize a new molecular brake on plasticity, namely nitration of TRKB receptors.

Remodeling of algal photosystem I through phosphorylation

Photosystem I (PSI) with its associated light-harvesting system is the most important generator of reducing power in photosynthesis. The PSI core complex is highly conserved, whereas peripheral subunits as well as light-harvesting proteins (LHCI) reveal a dynamic plasticity. Moreover, in green alga, PSI-LHCI complexes are found as monomers, dimers, and state transition complexes, where two LHCII trimers are associated. Herein, we show light-dependent phosphorylation of PSI subunits PsaG and PsaH as well as Lhca6. Potential consequences of the dynamic phosphorylation of PsaG and PsaH are structurally analyzed and discussed in regard to the formation of the monomeric, dimeric, and LHCII-associated PSI-LHCI complexes.

Phosphorylation Regulation Mechanism of β2 Integrin for the Binding of Filamin Revealed by Markov State Model

Leukocyte adhesion deficiency-1 (LAD-1) disorder is a severe immunodeficiency syndrome caused by deficiency or mutation of β2 integrin. The phosphorylation on threonine 758 of β2 integrin acts as a molecular switch inhibiting the binding of filamin. However, the switch mechanism of site-specific phosphorylation at the atom level is still poorly understood. To resolve the regulation mechanism, all-atom molecular dynamics simulation and Markov state model were used to study the dynamic regulation pathway of phosphorylation. Wild type system possessed lower binding free energy and fewer number of states than the phosphorylated system. Both systems underwent local disorder-to-order conformation conversion when achieving steady states. To reach steady states, wild type adopted less number of transition paths/shortest path according to the transition path theory than the phosphorylated system. The underlying phosphorylated regulation pathway was from P1 to P0 and then P4 state, and the main driving force should be hydrogen bond and hydrophobic interaction disturbing the secondary structure of phosphorylated states. These studies will shed light on the pathogenesis of LAD-1 disease and lay a foundation for drug development.

Disease-associated patterns of acetylation stabilize tau fibril formation

Assembly of the microtubule-associated protein into tauopathy fibril conformations dictates the pathology of a diversity of diseases. Recent cryogenic Electron Microscopy (cryo-EM) structures have uncovered distinct fibril conformations in different tauopathies but it remains unknown how these structures fold from a single protein sequence. It has been proposed that post-translational modifications may drive tau assembly but no direct mechanism for how modifications drive assembly has emerged. Leveraging established aggregation-regulating tau fragments that are normally inert, we tested the effect of chemical modification of lysines with acetyl groups on tau fragment conversion into amyloid aggregates. We identify specific patterns of acetylation that flank amyloidogenic motifs on the tau fragments that drive rapid fibril assembly. To understand how this pattern of acetylation may drive assembly, we determined a 3.9 Å cryo-EM structure of an amyloid fibril assembled from an acetylated tau fragment. The structure uncovers how lysine acetylation patterns mediate gain-of-function interactions to promote amyloid assembly. Comparison of the structure to an ex vivo tau fibril conformation from Pick's Disease reveals regions of high structural similarity. Finally, we show that our lysine- acetylated sequences exhibit fibril assembly activity in cell-based tau aggregation assays. Our data uncover the dual role of lysine residues in limiting aggregation while their acetylation leads to stabilizing pro-aggregation interactions. Design of tau sequence with specific acetylation patterns may lead to controllable tau aggregation to direct folding of tau into distinct folds.

Redox proteomics and structural analyses provide insightful implications for additional non-catalytic thiol-disulfide motifs in PDIs

Protein disulfide isomerases (PDIs) catalyze redox reactions that reduce, oxidize, or isomerize disulfide bonds and act as chaperones of proteins as they fold. The characteristic features of PDIs are the presence of one or more catalytic thioredoxin (TRX)-like domains harboring typical CXXC catalytic motifs responsible for redox reactions, as well as non-catalytic TRX-like domain. As increasing attention is paid to oxidative post-translational modifications of cysteines (Cys ox-PTMs) with the recognition that they control cellular signaling, strategies to identify sites of Cys ox-PTM by redox proteomics have been optimized. Exploration of an available Cys redoxome dataset supported by modeled structure provided arguments for the existence of an additional non-catalytic thiol-disulfide motif, distinct from those contained in the TRX type patterns, typical of PDIAs. Further structural analysis of PDIA3 and 6 allows us to consider the possibility that this hypothesis could be extended to other members of PDI. These elements invite future studies to decipher the exact role of these non-catalytic thiol-disulfide motifs in the functions of PDIs. Strategies that would allow to validate this hypothesis are discussed.

Umbrella Sampling Simulations Measure Switch Peptide Binding and Hydrophobic Patch Opening Free Energies in Cardiac Troponin

The cardiac troponin (cTn) complex is an important regulatory protein in heart contraction. Upon binding of Ca2+, cTn undergoes a conformational shift that allows the troponin I switch peptide (cTnISP) to be released from the actin filament and bind to the troponin C hydrophobic patch (cTnCHP). Mutations and modifications to this complex can change its sensitivity to Ca2+ and alter the energetics of the transition from the Ca2+-unbound, cTnISP-unbound form to the Ca2+-bound, cTnISP-bound form. We utilized targeted molecular dynamics (TMD) to obtain a trajectory of this transition pathway, followed by umbrella sampling to estimate the free energy associated with the cTnISP-cTnCHP binding and the cTnCHP opening events for wild-type (WT) cTn. We were able to reproduce experimental values for the cTnISP-cTnCHP binding event and obtain cTnCHP opening free energies in agreement with previous computational measurements of smaller cTnC systems. This excellent agreement for WT cTn demonstrated the strength of computational methods in studying the dynamics and energetics of the cTn complex. We then introduced mutations to the cTn complex that cause cardiomyopathy or alter its Ca2+ sensitivity and observed a general decrease in the free energy of opening the cTnCHP. For these same mutations, we observed no general trend in the effect on the cTnISP-cTnCHP binding event. Our method sets the stage for future computational studies on this system that predict the consequences of yet uncharacterized mutations on cTn dynamics and energetics.

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

The 2022 Lady Estelle Wolfson lectureship on neurofilaments

Neurofilament proteins (Nf) have been validated and established as a reliable body fluid biomarker for neurodegenerative pathology. This review covers seven Nf isoforms, Nf light (NfL), two splicing variants of Nf medium (NfM), two splicing variants of Nf heavy (NfH), α -internexin (INA) and peripherin (PRPH). The genetic and epigenetic aspects of Nf are discussed as relevant for neurodegenerative diseases and oncology. The comprehensive list of mutations for all Nf isoforms covers Amyotrophic Lateral Sclerosis, Charcot-Marie Tooth disease, Spinal muscular atrophy, Parkinson Disease and Lewy Body Dementia. Next, emphasis is given to the expanding field of post-translational modifications (PTM) of the Nf amino acid residues. Protein structural aspects are reviewed alongside PTMs causing neurodegenerative pathology and human autoimmunity. Molecular visualisations of NF PTMs, assembly and stoichiometry make use of Alphafold2 modelling. The implications for Nf function on the cellular level and axonal transport are discussed. Neurofilament aggregate formation and proteolytic breakdown are reviewed as relevant for biomarker tests and disease. Likewise, Nf stoichiometry is reviewed with regard to in vitro experiments and as a compensatory mechanism in neurodegeneration. The review of Nf across a spectrum of 87 diseases from all parts of medicine is followed by a critical appraisal of 33 meta-analyses on Nf body fluid levels. The review concludes with considerations for clinical trial design and an outlook for future research.

Anti-colon Cancer Effects of Dendrobium officinale Kimura & Migo Revealed by Network Pharmacology Integrated With Molecular Docking and Metabolomics Studies

Objective

The present study aimed to investigate the potential mechanism of Dendrobium officinale (D. officinale) on colorectal cancer and the relevant targets in the pathway using a network pharmacological approach.

Methods

(1) We identified the major bioactive components of D. officinale by UPLC-ESI-MS/MS and established the in-house library by using the literature mining method. (2) Target prediction was performed by SwissADME and SwissTargetPrediction. (3) A protein–protein interaction (PPI) network and component–target–pathway network (C-T-P network) were constructed. (4) The GO pathways and the KEGG pathway enrichment analysis were carried out by the Metascape database. (5) Molecular docking was performed by AutoDock software. (6) A series of experimental assays including cell proliferation, cell invasion and migration, and TUNEL staining in CRC were performed in CRC cell lines (HT-29, Lovo, SW-620, and HCT-116) to confirm the inhibitory effects of D. officinale.

Results

(1) In total, 396 candidate active components of D. officinale were identified by UPLC-ESI-MS/MS and selected from the database. (2) From OMIM, GeneCards, DrugBank, and TTD databases, 1,666 gene symbols related to CRC were gathered, and (3) 34 overlapping gene symbols related to CRC and drugs were obtained. (4) These results suggested that the anti-CRC components of D. officinale were mainly apigenin, naringenin, caffeic acid, γ-linolenic acid, α-linolenic acid, cis-10-heptadecenoic acid, etc., and the core targets of action were mainly ESR1, EGFR, PTGS2, MMP9, MMP2, PPARG, etc. (5) The proliferation of muscle cells, the regulation of inflammatory response, the response of cells to organic cyclic compounds, and the apoptotic signaling pathway might serve as principal pathways for CRC treatment. (6) The reliability of some important active components and targets was further validated by molecular docking. The molecular docking analysis suggested an important role of apigenin, naringenin, PTGS2, and MMP9 in delivering the pharmacological activity of D. officinale against CRC. (7) These results of the evaluation experiment in vitro suggested that D. officinale had a strong inhibitory effect on CRC cell lines, and it exerted anti-CRC activity by activating CRC cell apoptosis and inhibiting CRC cell migration and invasion.

Conclusion

This study may provide valuable insights into exploring the mechanism of action of D. officinale against CRC.

EZH2 engages TGFβ signaling to promote breast cancer bone metastasis via integrin β1-FAK activation

Bone metastases occur in 50-70% of patients with late-stage breast cancers and effective therapies are needed. The expression of enhancer of zeste homolog 2 (EZH2) is correlated with breast cancer metastasis, but its function in bone metastasis hasn't been well-explored. Here we report that EZH2 promotes osteolytic metastasis of breast cancer through regulating transforming growth factor beta (TGFβ) signaling. EZH2 induces cancer cell proliferation and osteoclast maturation, whereas EZH2 knockdown decreases bone metastasis incidence and outgrowth in vivo. Mechanistically, EZH2 transcriptionally increases ITGB1, which encodes for integrin β1. Integrin β1 activates focal adhesion kinase (FAK), which phosphorylates TGFβ receptor type I (TGFβRI) at tyrosine 182 to enhance its binding to TGFβ receptor type II (TGFβRII), thereby activating TGFβ signaling. Clinically applicable FAK inhibitors but not EZH2 methyltransferase inhibitors effectively inhibit breast cancer bone metastasis in vivo. Overall, we find that the EZH2-integrin β1-FAK axis cooperates with the TGFβ signaling pathway to promote bone metastasis of breast cancer.

Novel Post-translational Modifications in Human Serum Albumin

Aim:

This study aims to identify novel post-translational modifications in human serum albumin by mass spectrometry.

Background:

Serum albumin is the most abundant protein in plasma, has many physiological functions, and is in contact with most of the cells and tissues of the human body. Post-translational modifications (PTMs) may affect functions, stability, and localization of albumin.

Objective:

Identify novel PTMs in human serum albumin by mass spectrometry.

Methods:

Human serum albumin (HSA) was used for tryptic digestion in-solution or in-gel. Mass spectrometry was applied to identify PTMs in HSA. 3-dimensional modeling was applied to explore the potential impact of PTMs on known functions of albumin.

Results:

Here, we report the identification of 61 novel PTMs of human serum albumin. Phosphorylation, glycosylation, nitrosylation, deamidation, methylation, acetylation, palmitoylation, geranylation, and farnesylation are some examples of the identified PTMs. Mass spectrometry was used for the identification of PTMs in a purified HSA and HSA from the human plasma. Threedimensional modeling of albumin with selected PTMs showed the location of these PTMs in the regions involved in albumin interactions with drugs, metals, and fatty acids. The location of PTMs in these regions may modify the binding capacity of albumin.

Conclusion:

This report adds 61 novel PTMs to the catalog of human albumin.

Sirtuin-1 sensitive lysine-136 acetylation drives phase separation and pathological aggregation of TDP-43

Trans-activation response DNA-binding protein of 43  kDa (TDP-43) regulates RNA processing and forms neuropathological aggregates in patients with amyotrophic lateral sclerosis and frontotemporal lobar degeneration. Investigating TDP-43 post-translational modifications, we discovered that K84 acetylation reduced nuclear import whereas K136 acetylation impaired RNA binding and splicing capabilities of TDP-43. Such failure of RNA interaction triggered TDP-43 phase separation mediated by the C-terminal low complexity domain, leading to the formation of insoluble aggregates with pathologically phosphorylated and ubiquitinated TDP-43. Introduction of acetyl-lysine at the identified sites via amber suppression confirmed the results from site-directed mutagenesis. K84-acetylated TDP-43 showed cytoplasmic mislocalization, and the aggregation propensity of K136-acetylated TDP-43 was confirmed. We generated antibodies selective for TDP-43 acetylated at these lysines, and found that sirtuin-1 can potently deacetylate K136-acetylated TDP-43 and reduce its aggregation propensity. Thus, distinct lysine acetylations modulate nuclear import, RNA binding and phase separation of TDP-43, suggesting regulatory mechanisms for TDP-43 pathogenesis.

PEA-15 engages in allosteric interactions using a common scaffold in a phosphorylation-dependent manner

Phosphoprotein enriched in astrocytes, 15 kDa (PEA-15) is a death-effector domain (DED) containing protein involved in regulating mitogen-activated protein kinase and apoptosis pathways. In this molecular dynamics study, we examined how phosphorylation of the PEA-15 C-terminal tail residues, Ser-104 and Ser-116, allosterically mediates conformational changes of the DED and alters the binding specificity from extracellular-regulated kinase (ERK) to Fas-associated death domain (FADD) protein. We delineated that the binding interfaces between the unphosphorylated PEA-15 and ERK2 and between the doubly phosphorylated PEA-15 and FADD are similarly composed of a scaffold that includes both the DED and the C-terminal tail residues of PEA-15. While the unphosphorylated serine residues do not directly interact with ERK2, the phosphorylated Ser-116 engages in strong electrostatic interactions with arginine residues on FADD DED. Upon PEA-15 binding, FADD repositions its death domain (DD) relative to the DED, an essential conformational change to allow the death-inducing signaling complex (DISC) assembly.

Phosphorylation tunes elongation propensity and cohesiveness of INCENP's intrinsically disordered region

The inner centromere protein, INCENP, is crucial for correct chromosome segregation during mitosis. It connects the kinase Aurora B to the inner centromere allowing this kinase to dynamically access its kinetochore targets. However, the function of its central, 440-residue long intrinsically disordered region (IDR) and its multiple phosphorylation sites is unclear. Here, we determined the conformational ensemble of INCENP's IDR, systematically varying the level of phosphorylation, using all-atom and coarse-grain molecular dynamics simulations. Our simulations show that phosphorylation expands INCENP's IDR, both locally and globally, mainly by increasing its overall net charge. The disordered region undergoes critical globule-to-coil conformational transitions and the transition temperature non-monotonically depends on the degree of phosphorylation, with a mildly phosphorylated case of neutral net charge featuring the highest collapse propensity. The IDR transitions from a multitude of globular states, accompanied by several specific internal contacts that reduce INCENP length by loop formation, to weakly interacting and highly extended coiled conformations. Phosphorylation critically shifts the population between these two regimes. It thereby influences cohesiveness and phase behavior of INCENP IDR assemblies, a feature presumably relevant for INCENP's function in the chromosomal passenger complex. Overall, we propose the disordered region of INCENP to act as a phosphorylation-regulated and length-variable component, within the previously defined "dog-leash" model, that thereby regulates how Aurora B reaches its targets for proper chromosome segregation.

FTLD-TDP assemblies seed neoaggregates with subtype-specific features via a prion-like cascade

Morphologically distinct TDP-43 aggregates occur in clinically different FTLD-TDP subtypes, yet the mechanism of their emergence and contribution to clinical heterogeneity are poorly understood. Several lines of evidence suggest that pathological TDP-43 follows a prion-like cascade, but the molecular determinants of this process remain unknown. We use advanced microscopy techniques to compare the seeding properties of pathological FTLD-TDP-A and FTLD-TDP-C aggregates. Upon inoculation of patient-derived aggregates in cells, FTLD-TDP-A seeds amplify in a template-dependent fashion, triggering neoaggregation more efficiently than those extracted from FTLD-TDP-C patients, correlating with the respective disease progression rates. Neoaggregates are sequentially phosphorylated with N-to-C directionality and with subtype-specific timelines. The resulting FTLD-TDP-A neoaggregates are large and contain densely packed fibrils, reminiscent of the pure compacted fibrils present within cytoplasmic inclusions in postmortem brains. In contrast, FTLD-TDP-C dystrophic neurites show less dense fibrils mixed with cellular components, and their respective neoaggregates are small, amorphous protein accumulations. These cellular seeding models replicate aspects of the patient pathological diversity and will be a useful tool in the quest for subtype-specific therapeutics.

BTK operates a phospho-tyrosine switch to regulate NLRP3 inflammasome activity

Activity of the NLRP3 inflammasome, a critical mediator of inflammation, is controlled by accessory proteins, posttranslational modifications, cellular localization, and oligomerization. How these factors relate is unclear. We show that a well-established drug target, Bruton’s tyrosine kinase (BTK), affects several levels of NLRP3 regulation. BTK directly interacts with NLRP3 in immune cells and phosphorylates four conserved tyrosine residues upon inflammasome activation, in vitro and in vivo. Furthermore, BTK promotes NLRP3 relocalization, oligomerization, ASC polymerization, and full inflammasome assembly, probably by charge neutralization, upon modification of a polybasic linker known to direct NLRP3 Golgi association and inflammasome nucleation. As NLRP3 tyrosine modification by BTK also positively regulates IL-1β release, we propose BTK as a multifunctional positive regulator of NLRP3 regulation and BTK phosphorylation of NLRP3 as a novel and therapeutically tractable step in the control of inflammation.

Development of High Affinity Calcitonin Analog Fragments Targeting Extracellular Domains of Calcitonin Family Receptors

The calcitonin and amylin receptors (CTR and AMY receptors) are the drug targets for osteoporosis and diabetes treatment, respectively. Salmon calcitonin (sCT) and pramlintide were developed as peptide drugs that activate these receptors. However, next-generation drugs with improved receptor binding profiles are desirable for more effective pharmacotherapy. The extracellular domain (ECD) of CTR was reported as the critical binding site for the C-terminal half of sCT. For the screening of high-affinity sCT analog fragments, purified CTR ECD was used for fluorescence polarization/anisotropy peptide binding assay. When three mutations (N26D, S29P, and P32HYP) were introduced to the sCT(22–32) fragment, sCT(22–32) affinity for the CTR ECD was increased by 21-fold. CTR was reported to form a complex with receptor activity-modifying protein (RAMP), and the CTR:RAMP complexes function as amylin receptors with increased binding for the peptide hormone amylin. All three types of functional AMY receptor ECDs were prepared and tested for the binding of the mutated sCT(22–32). Interestingly, the mutated sCT(22–32) also retained its high affinity for all three types of the AMY receptor ECDs. In summary, the mutated sCT(22–32) showing high affinity for CTR and AMY receptor ECDs could be considered for developing the next-generation peptide agonists.

Oxygen-dependent changes in binding partners and post-translational modifications regulate the abundance and activity of HIF-1α/2α

Cellular adaptation to low-oxygen environments is mediated in part by the hypoxia-inducible factors (HIFs). Like other transcription factors, the stability and transcriptional activity of HIFs-and consequently, the hypoxic response-are regulated by post-translational modifications (PTMs) and changes in protein-protein interactions. Our current understanding of PTM-mediated regulation of HIFs is primarily based on in vitro protein fragment-based studies typically validated in fragment-expressing cells treated with hypoxia-mimicking compounds. Here, we used immunoprecipitation-based mass spectrometry to characterize the PTMs and binding partners for full-length HIF-1α and HIF-2α under normoxic (21% oxygen) and hypoxic (1% oxygen) conditions. Hypoxia substantially altered the complexity and composition of the HIFα protein interaction networks, particularly for HIF-2α, with the hypoxic networks of both isoforms being enriched for mitochondrial proteins. Moreover, both HIFα isoforms were heavily covalently modified. We identified ~40 PTM sites composed of 13 different types of modification on both HIFα isoforms, including multiple cysteine modifications and an unusual phosphocysteine. More than 80% of the PTMs identified were not previously known and about half exhibited oxygen dependency. We further characterized an evolutionarily conserved phosphorylation of Ser31 in HIF-1α as a regulator of its transcriptional function, and we propose functional roles for Thr406, Thr528, and Ser581 in HIF-2α. These data will help to delineate the different physiological roles of these closely related isoforms in fine-tuning the hypoxic response.

Computational identification and characterization of antigenic properties of Rv3899c of Mycobacterium tuberculosis and its interaction with human leukocyte antigen (HLA)

A rise in drug-resistant tuberculosis (TB) cases demands continued efforts towards the discovery and development of drugs and vaccines. Secretory proteins of Mycobacterium tuberculosis (H37Rv) are frequently studied for their antigenicity and their scope as protein subunit vaccines requires further analysis. In this study, Rv3899c of H37Rv emerges as a potential vaccine candidate on its evaluation by several bioinformatics tools. It is a non-toxic, secretory protein with an 'immunoglobulin-like' fold which does not show similarity with a human protein. Through BlastP and MEME suite analysis, we found Rv3899c homologs in several mycobacterial species and its antigenic score (0.54) to compare well with the known immunogens such as ESAT-6 (0.56) and Rv1860 (0.52). Structural examination of Rv3899c predicted ten antigenic peptides, an accessibility profile of the antigenic determinants constituting B cell epitope-rich regions and a low abundance of antigenic regions (AAR) value. Significantly, STRING analysis showed ESX-2 secretion system proteins and antigenic PE/PPE proteins of H37Rv as the interacting partners of Rv3899c. Further, molecular docking predicted Rv3899c to interact with human leukocyte antigen HLA-DRB1*04:01 through its antigenically conserved motif (RAAEQQRLQRIVDAVARQEPRISWAAGLRDDGTT). Interestingly, the binding affinity was observed to increase on citrullination of its Arg1 residue. Taken together, the computational characterization and predictive information suggest Rv3899c to be a promising TB vaccine candidate, which should be validated experimentally.

Computational Methods Elucidate Consequences of Mutations and Post-translational Modifications on Troponin I Effective Concentration to Troponin C

Ca²⁺ binding to cardiac troponin C (cTnC) causes a conformational shift that exposes a hydrophobic patch (cTnC_HP) for binding of the cTnI switch peptide (cTnI_SP), ultimately resulting in contraction of the heart. The inhibitory peptide (cTnI_IP), attached at the N-terminal end of the cTnI_SP, serves as a tether for the cTnI_SP to the rest of the troponin complex. Due to this tethered nature, the cTnI_SP remains within proximity of the hydrophobic patch region, resulting in the cTnC_HP experiencing an elevated "effective concentration" of the cTnI_SP. Mutations to the cTnI_IP region have been hypothesized to cause disease by affecting the ability of the cTnI_SP to "find" the hydrophobic patch, resulting in alterations to the heart's ability to contract normally. We tested this hypothesis using molecular dynamics (MD) simulations of the troponin complex using a model that contained all three subunits of troponin: C, I, and T. We developed methods that allowed us to quantitatively measure the effective concentration of the cTnI_SP from the simulations. A significant reduction in the cTnI_SP effective concentration was observed when the cTnI_IP was removed from the system, showcasing the importance of a tethered cTnI_SP. Through accelerated MD methods, we proposed the minimum effective concentration of a tethered cTnI_SP to be approximately 21 mM. Modification of the cTnI_IP via PKC-mediated phosphorylation of T143 was shown to significantly increase the estimated effective concentration of cTnI_SP, help the cTnI_SP find the cTnC_HP more effectively, and maintain the relative shape of the cTnI_IP when compared to the native model. All of these data indicate that pT143 may be able to help promote binding of cTnI_SP to the cTnC_HP. We then tested six mutations within the cTnI_IP region that are known cTnC Ca²⁺-sensitizing mutations and have been linked with cardiomyopathy. We did not observe a significant reduction in the effective concentration upon the introduction of these mutations; however, we did observe increased variability in the flexibility and dynamics of the cTnI_IP region when compared to native. Our observations led us to hypothesize that the mechanism by which these cardiomyopathic mutations affect Ca²⁺ sensitivity is by altering the off rate of cTnI_SP from the hydrophobic patch.

PKCε controls the fusion of secretory vesicles in mast cells in a phosphatidic acid-dependent mode

PKCε is highly expressed in mast cells and plays a fundamental role in the antigen-triggered activation of the allergic reaction. Although its regulation by diacylglycerols has been described, its regulation by acidic phospholipids and how this regulation leads to the control of downstream vesicle secretion is barely known. Here, we used structural and evolutionary studies to find the molecular mechanism that explains the selectivity of the C1B domain of PKCε by Phosphatidic Acid (PA). This resided in a collection of Arg residues that form a specific rim on the outer surface of the C1B domain, around the diacylglycerol binding cleft. In RBL-2H3 cells, this basic rim allowed the kinase to respond specifically to phosphatidic acid signals that induced its translocation to the plasma membrane and subsequent activation. Further experiments in cells that overexpress PKCε and a mutant of the PA binding site, showed that PA-dependent PKCε activation increased vesicle degranulation in RBL-2H3 cells, and this correlated with increased SNAP23 phosphorylation. Over-expression of PKCε in these cells also induced an increase in the number of docked vesicles containing SNAP23, when stimulated with PA. This accumulation could be attributed to the stabilizing effect of phosphorylation on the formation of the SNARE complex, which ultimately led to increased release of content in the presence of Ca²⁺ during the fusion process. Therefore, these findings reinforce the importance of PA signaling in the activation of PKCε, which could be an important target to inhibit the exacerbated responses of these cells in the allergic reaction.

Multilevel proteomics reveals host perturbations by SARS-CoV-2 and SARS-CoV

The emergence and global spread of SARS-CoV-2 has resulted in the urgent need for an in-depth understanding of molecular functions of viral proteins and their interactions with the host proteome. Several individual omics studies have extended our knowledge of COVID-19 pathophysiology^1-10. Integration of such datasets to obtain a holistic view of virus-host interactions and to define the pathogenic properties of SARS-CoV-2 is limited by the heterogeneity of the experimental systems. Here we report a concurrent multi-omics study of SARS-CoV-2 and SARS-CoV. Using state-of-the-art proteomics, we profiled the interactomes of both viruses, as well as their influence on the transcriptome, proteome, ubiquitinome and phosphoproteome of a lung-derived human cell line. Projecting these data onto the global network of cellular interactions revealed crosstalk between the perturbations taking place upon infection with SARS-CoV-2 and SARS-CoV at different levels and enabled identification of distinct and common molecular mechanisms of these closely related coronaviruses. The TGF-β pathway, known for its involvement in tissue fibrosis, was specifically dysregulated by SARS-CoV-2 ORF8 and autophagy was specifically dysregulated by SARS-CoV-2 ORF3. The extensive dataset (available at https://covinet.innatelab.org ) highlights many hotspots that could be targeted by existing drugs and may be used to guide rational design of virus- and host-directed therapies, which we exemplify by identifying inhibitors of kinases and matrix metalloproteases with potent antiviral effects against SARS-CoV-2.

Identification of imidazoquinoline derivative (IQD) interacting sites of porcine TLR8 and the underlying species specificity

Toll-like receptor 8 (TLR8), as an important innate immune receptor, can recognize specific ligands, activate intracellular signaling and produce an inflammatory response to kill and eliminate pathogenic microorganisms. Recent studies have resolved the crystal structure of human TLR8 (hTLR8) and two types of ligand binding sites were identified. Among the conserved binding site 1 of hTLR8, the residues interacting with imidazoquinoline derivatives (IQDs) were determined. We previously showed that porcine TLR8 (pTLR8) exhibited species specificity for recognition of the hTLR7 agonist imiquimod (R837). Given the species specificity, the pTLR8 residues interacting with IQDs may be different from hTLR8 counterparts. The present study was aimed to identify the pTLR8 residues interacting with small molecular IQDs. Via molecular docking, the pTLR8 residues interacting with R837 and R848 were predicted. The corresponding mutants were tested for pTLR8 signaling in response to IQDs R837, R848 and CL075, and the results showed that five of nine predicted residues (Y336, K341, K342, F395 and G562) are critical for pTLR8 signaling and these residues are partially different from those reported in hTLR8. Further, we found that the pTLR8 GQKNG motif corresponding to hTLR8 RQSYA exhibited disparity to CL075 stimulation. Our study thus reveals fine TLR8 species specificity which deepens the understanding of TLR8 activation mechanism.

Molecular dynamics shows complex interplay and long-range effects of post-translational modifications in yeast protein interactions

Post-translational modifications (PTMs) play a vital, yet often overlooked role in the living cells through modulation of protein properties, such as localization and affinity towards their interactors, thereby enabling quick adaptation to changing environmental conditions. We have previously benchmarked a computational framework for the prediction of PTMs’ effects on the stability of protein-protein interactions, which has molecular dynamics simulations followed by free energy calculations at its core. In the present work, we apply this framework to publicly available data on Saccharomyces cerevisiae protein structures and PTM sites, identified in both normal and stress conditions. We predict proteome-wide effects of acetylations and phosphorylations on protein-protein interactions and find that acetylations more frequently have locally stabilizing roles in protein interactions, while the opposite is true for phosphorylations. However, the overall impact of PTMs on protein-protein interactions is more complex than a simple sum of local changes caused by the introduction of PTMs and adds to our understanding of PTM cross-talk. We further use the obtained data to calculate the conformational changes brought about by PTMs. Finally, conservation of the analyzed PTM residues in orthologues shows that some predictions for yeast proteins will be mirrored to other organisms, including human. This work, therefore, contributes to our overall understanding of the modulation of the cellular protein interaction networks in yeast and beyond.

Proteins are a diverse set of biological molecules responsible for numerous functions within cells, such as obtaining energy from food or transport of small molecules, and many processes rely on interactions of specific proteins. Moreover, a single protein may acquire different roles depending on cellular requirements and as a response to changes in the environment. A commonly used way to quickly change protein’s function or activity is by introducing small chemical modifications on specific locations within the protein. These modifications can cause the protein to interact in a more or less stable way with other proteins. We have previously developed a computational pipeline for predicting the effect of modifications on interactions of proteins, and in this work we apply it to all yeast proteins with known structures. We find differences in effects on the binding for different types of modifications. Importantly, we demonstrate that the modifications far from the interaction interface also significantly contribute to binding due to their impact on protein’s shape, which is often neglected by other methods. This work contributes to our understanding of the modulation of protein interactions in yeast due to modifications, while our widely applicable method will allow similar investigations in other organisms.

How multisite phosphorylation impacts the conformations of intrinsically disordered proteins

Phosphorylation of intrinsically disordered proteins (IDPs) can produce changes in structural and dynamical properties and thereby mediate critical biological functions. How phosphorylation effects intrinsically disordered proteins has been studied for an increasing number of IDPs, but a systematic understanding is still lacking. Here, we compare the collapse propensity of four disordered proteins, Ash1, the C-terminal domain of RNA polymerase (CTD2’), the cytosolic domain of E-Cadherin, and a fragment of the p130Cas, in unphosphorylated and phosphorylated forms using extensive all-atom molecular dynamics (MD) simulations. We find all proteins to show V-shape changes in their collapse propensity upon multi-site phosphorylation according to their initial net charge: phosphorylation expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. However, force fields including those tailored towards and commonly used for IDPs overestimate these changes. We find quantitative agreement of MD results with SAXS and NMR data for Ash1 and CTD2’ only when attenuating protein electrostatic interactions by using a higher salt concentration (e.g. 350 mM), highlighting the overstabilization of salt bridges in current force fields. We show that phosphorylation of IDPs also has a strong impact on the solvation of the protein, a factor that in addition to the actual collapse or expansion of the IDP should be considered when analyzing SAXS data. Compared to the overall mild change in global IDP dimension, the exposure of active sites can change significantly upon phosphorylation, underlining the large susceptibility of IDP ensembles to regulation through post-translational modifications.

Intrinsically disordered proteins (IDPs) are a class of proteins that lack secondary and tertiary structures and instead explore a broad conformational ensemble. Their functions, from transcriptional regulation to signal transmission, are tightly regulated. IDPs are subject of extensive reversible post-translational modifications (PTMs), such as phosphorylation, methylation and glycosylation. Among these PTMs, phosphorylation is one of the most common and important PTMs. However, the mechanism of how phosphorylation affects the conformations and functions of IDPs remains unclear. To answer this question, we have performed extensive all-atom molecular dynamics simulations for four representative IDPs: Ash1, E-Cadherin, CTD2’ and p130Cas in their unphosphorylated and phosphorylated forms. Our results showed that all IDPs undergo a mild change upon multi-site phosphorylation, which is V-shaped: phosphorylation moderately expands neutral or overall negatively charged IDPs and shrinks positively charged IDPs. More importantly, in two of these IDPs, only two biologically relevant phosphorylation sites suffice to render the adjacent negatively charged active site significantly more exposed to the environment, which implies a higher probability to interact with other binding partners.

Pinpointing cysteine oxidation sites by high-resolution proteomics reveals a mechanism of redox-dependent inhibition of human STING

Protein function is regulated by posttranslational modifications (PTMs), among which reversible oxidation of cysteine residues has emerged as a key regulatory mechanism of cellular responses. Given the redox regulation of virus-host interactions, the identification of oxidized cysteine sites in cells is essential to understand the underlying mechanisms involved. Here, we present a proteome-wide identification of reversibly oxidized cysteine sites in oxidant-treated cells using a maleimide-based bioswitch method coupled to mass spectrometry analysis. We identified 2720 unique oxidized cysteine sites within 1473 proteins with distinct abundances, locations, and functions. Oxidized cysteine sites were found in numerous signaling pathways, many relevant to virus-host interactions. We focused on the oxidation of STING, the central adaptor of the innate immune type I interferon pathway, which is stimulated in response to the detection of cytosolic DNA by cGAS. We demonstrated the reversible oxidation of Cys¹⁴⁸ and Cys²⁰⁶ of STING in cells. Molecular analyses led us to establish a model in which Cys¹⁴⁸ oxidation is constitutive, whereas Cys²⁰⁶ oxidation is inducible by oxidative stress or by the natural ligand of STING, 2'3'-cGAMP. Our data suggest that the oxidation of Cys²⁰⁶ prevented hyperactivation of STING by causing a conformational change associated with the formation of inactive polymers containing intermolecular disulfide bonds. This finding should aid the design of therapies targeting STING that are relevant to autoinflammatory disorders, immunotherapies, and vaccines.

Regulation of Phosphoinositide Levels in the Retina by Protein Tyrosine Phosphatase 1B and Growth Factor Receptor-Bound Protein 14

Protein tyrosine kinases and protein phosphatases play a critical role in cellular regulation. The length of a cellular response depends on the interplay between activating protein kinases and deactivating protein phosphatases. Protein tyrosine phosphatase 1B (PTP1B) and growth factor receptor-bound protein 14 (Grb14) are negative regulators of receptor tyrosine kinases. However, in the retina, we have previously shown that PTP1B inactivates insulin receptor signaling, whereas phosphorylated Grb14 inhibits PTP1B activity. In silico docking of phosphorylated Grb14 and PTP1B indicate critical residues in PTP1B that may mediate the interaction. Phosphoinositides (PIPs) are acidic lipids and minor constituents in the cell that play an important role in cellular processes. Their levels are regulated by growth factor signaling. Using phosphoinositide binding protein probes, we observed increased levels of PI(3)P, PI(4)P, PI(3,4)P₂, PI(4,5)P2, and PI(3,4,5)P₃ in PTP1B knockout mouse retina and decreased levels of these PIPs in Grb14 knockout mouse retina. These observations suggest that the interplay between PTP1B and Grb14 can regulate PIP metabolism.

Conformational interconversion of MLKL and disengagement from RIPK3 precede cell death by necroptosis

Phosphorylation of the MLKL pseudokinase by the RIPK3 kinase leads to MLKL oligomerization, translocation to, and permeabilization of, the plasma membrane to induce necroptotic cell death. The precise choreography of MLKL activation remains incompletely understood. Here, we report Monobodies, synthetic binding proteins, that bind the pseudokinase domain of MLKL within human cells and their crystal structures in complex with the human MLKL pseudokinase domain. While Monobody-32 constitutively binds the MLKL hinge region, Monobody-27 binds MLKL via an epitope that overlaps the RIPK3 binding site and is only exposed after phosphorylated MLKL disengages from RIPK3 following necroptotic stimulation. The crystal structures identified two distinct conformations of the MLKL pseudokinase domain, supporting the idea that a conformational transition accompanies MLKL disengagement from RIPK3. These studies provide further evidence that MLKL undergoes a large conformational change upon activation, and identify MLKL disengagement from RIPK3 as a key regulatory step in the necroptosis pathway.

Insight Into Seeded Tau Fibril Growth From Molecular Dynamics Simulation of the Alzheimer's Disease Protofibril Core

Aggregates of the microtubule associated tau protein are a major constituent of neurofibrillary lesions that define Alzheimer’s disease (AD) pathology. Increasing experimental evidence suggests that the spread of tau neurofibrillary tangles results from a prion-like seeding mechanism in which small oligomeric tau fibrils template the conversion of native, intrinsically disordered, tau proteins into their pathological form. By using atomistic molecular dynamics (MD) simulations, we investigate the stability and dissociation thermodynamics of high-resolution cryo-electron microscopy (cryo-EM) structures of both the AD paired-helical filament (PHF) and straight filament (SF). Non-equilibrium steered MD (SMD) center-of-mass pulling simulations are used to probe the stability of the protofibril structure and identify intermolecular contacts that must be broken before a single tau peptide can dissociate from the protofibril end. Using a combination of exploratory metadynamics and umbrella sampling, we investigate the complete dissociation pathway and compute a free energy profile for the dissociation of a single tau peptide from the fibril end. Different features of the free energy surface between the PHF and SF protofibril result from a different mechanism of tau unfolding. Comparison of wild-type tau PHF and post-translationally modified pSer356 tau shows that phosphorylation at this site changes the dissociation free energy surface of the terminal peptide. These results demonstrate how different protofibril morphologies template the folding of endogenous tau in distinct ways, and how post-translational modification can perturb the folding mechanism.

Stable Isotope Labeling of Amino Acids in Flies (SILAF) Reveals Differential Phosphorylation of Mitochondrial Proteins Upon Loss of OXPHOS Subunits

Drosophila melanogaster has been a workhorse of genetics and cell biology for more than a century. However, proteomic-based methods have been limited due to the complexity and dynamic range of the fly proteome and the lack of efficient labeling methods. Here, we advanced a chemically defined food source into direct stable-isotope labeling of amino acids in flies (SILAF). It allows for the rapid and cost-efficient generation of a large number of larvae or flies, with full incorporation of lysine-[13C6] after six labeling days. SILAF followed by fractionation and enrichment gave proteomic insights at a depth of 7196 proteins and 8451 phosphorylation sites, which substantiated metabolic regulation on enzymatic level. We applied SILAF to quantify the mitochondrial phosphoproteome of an early-stage leucine-rich PPR motif-containing protein (LRPPRC)-knockdown fly model of mitochondrial disease that almost exclusively affects protein levels of the oxidative phosphorylation (OXPHOS) system. While the mitochondrial compartment was hypo-phosphorylated, two conserved phosphosites on OXPHOS subunits NDUFB10 and NDUFA4 were significantly upregulated upon impaired OXPHOS function. The ease and versatility of the method actuate the fruit fly as an appealing model in proteomic and posttranslational modification studies, and it enlarges potential metabolic applications based on heavy amino acid diets.

Phosphorylation-induced changes in the PDZ domain of Dishevelled 3

The PDZ domain of Dishevelled 3 protein belongs to a highly abundant protein recognition motif which typically binds short C-terminal peptides. The affinity of the PDZ towards the peptides could be fine-tuned by a variety of post-translation modifications including phosphorylation. However, how phosphorylations affect the PDZ structure and its interactions with ligands remains elusive. Combining molecular dynamics simulations, NMR titration, and biological experiments, we explored the role of previously reported phosphorylation sites and their mimetics in the Dishevelled PDZ domain. Our observations suggest three major roles for phosphorylations: (1) acting as an on/off PDZ binding switch, (2) allosterically affecting the binding groove, and (3) influencing the secondary binding site. Our simulations indicated that mimetics had similar but weaker effects, and the effects of distinct sites were non-additive. This study provides insight into the Dishevelled regulation by PDZ phosphorylation. Furthermore, the observed effects could be used to elucidate the regulation mechanisms in other PDZ domains.

Molecular dynamics simulations and biochemical characterization of Pf14-3-3 and PfCDPK1 interaction towards its role in growth of human malaria parasite

Scaffold proteins play pivotal role as modulators of cellular processes by operating as multipurpose conformation clamps. 14-3-3 proteins are gold-standard scaffold modules that recognize phosphoSer/Thr (pS/pT) containing conserved motifs, and confer conformational changes leading to modulation of functional parameters of their target proteins. Modulation in functional activity of kinases has been attributed to their interaction with 14-3-3 proteins. Herein, we have annotated and characterized PF3D7_0818200 as 14-3-3 isoform I in Plasmodium falciparum 3D7, and its interaction with one of the key kinases of the parasite, Calcium-Dependent Protein Kinase 1 (CDPK1) by performing various analytical biochemistry and biophysical assays. Molecular dynamics simulation studies indicated that CDPK1 polypeptide sequence (61KLGpS64) behaves as canonical Mode I-type (RXXpS/pT) consensus 14-3-3 binding motif, mediating the interaction. The 14-3-3I/CDPK1 interaction was validated in vitro with ELISA and SPR, which confirmed that the interaction is phosphorylation dependent, with binding affinity constant of 670 ± 3.6 nM. The interaction of 14-3-3I with CDPK1 was validated with well characterized optimal 14-3-3 recognition motifs: Mode I-type ARSHpSYPA and Mode II-type RLYHpSLPA, by simulation studies and ITC. This interaction was found to marginally enhance CDPK1 functional activity. Furthermore, interaction antagonizing peptidomimetics showed growth inhibitory impact on the parasite indicating crucial physiological role of 14-3-3/CDPK1 interaction. Overall, this study characterizes 14-3-3I as a scaffold protein in the malaria parasite and unveils CDPK1 as its previously unidentified target. This sets a precedent for the rational design of 14-3-3 based PPI inhibitors by utilizing 14-3-3 recognition motif peptides, as a potential antimalarial strategy.

Effect of lysine acetylation on the regulation of Trypanosoma brucei glycosomal aldolase activity

Post-translational modifications provide suitable mechanisms for cellular adaptation to environmental changes. Lysine acetylation is one of these modifications and occurs with the addition of an acetyl group to Nε-amino chain of this residue, eliminating its positive charge. Recently, we found distinct acetylation profiles of procyclic and bloodstream forms of Trypanosoma brucei, the agent of African Trypanosomiasis. Interestingly, glycolytic enzymes were more acetylated in the procyclic, which develops in insects and uses oxidative phosphorylation to obtain energy, compared with the bloodstream form, whose main source of energy is glycolysis. Here, we investigated whether acetylation regulates the T. brucei fructose 1,6-bisphosphate aldolase. We found that aldolase activity was reduced in procyclic parasites cultivated in the absence of glucose and partial recovered by in vitro deacetylation. Similarly, acetylation of protein extracts from procyclics cultivated in glucose-rich medium, caused a reduction in the aldolase activity. In addition, aldolase acetylation levels were higher in procyclics cultivated in the absence of glucose compared with those cultivated in the presence of glucose. To further confirm the role of acetylation, lysine residues near the catalytic site were substituted by glutamine in recombinant T. brucei aldolase. These replacements, especially K157, inhibited enzymatic activity, changed the electrostatic surface potential, decrease substrate binding and modify the catalytic pocket structure of the enzyme, as predicted by in silico analysis. Taken together, these data confirm the role of acetylation in regulating the activity of an enzyme from the glycolytic pathway of T. brucei, expanding the factors responsible for regulating important pathways in this parasite.

PyVibMS: a PyMOL plugin for visualizing vibrations in molecules and solids

Visualizing vibrational motions calculated with different ab initio packages requires dedicated post-processing tools. Here, we present a PyMOL plugin called PyVibMS for visualizing the vibrational motions for both molecular and solid systems calculated by mainstream quantum chemical computer programs including Gaussian, Q-Chem, VASP, and CRYSTAL. Benefiting from the continuing development of the PyMOL platform, PyVibMS provides powerful functionalities and user-friendly interface. PyVibMS was written in Python and its open-source nature makes it flexible and sustainable. As an example, the motions of the Konkoli-Cremer local vibrational modes are shown in this work for the first time. PyVibMS is freely available at https://github.com/smutao/PyVibMS . Graphical abstract In this work, a PyMOL plugin named PyVibMS is developed to visualize molecular and lattice vibrations.

A Network Pharmacological Approach to Investigate the Mechanism of Action of Active Ingredients of Epimedii Herba and Their Potential Targets in Treatment of Alzheimer's Disease

Background

Epimedii Herba is a traditional Chinese herbal medicine used to treat central nervous system diseases such as Alzheimer’s disease in China. However, the pharmacological mechanism is unclear. To investigate the mechanisms of Epimedii Herba in the treatment of Alzheimer’s disease, we assessed effective compounds, corresponding targets, and related pathways of Epimedii Herba in the treatment of Alzheimer’s disease based on network pharmacology.

Material/Methods

The active components and targets of Epimedii Herba were obtained through the TCMSP database and the DrugBank database. The DisGeNET database and GeneCards database were used to search for Alzheimer’s disease targets. The common targets of components and disease were obtained by Wayne diagram. Gene ontology (GO) analysis and enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed using the DAVID database. The component-target-pathway interaction network model was constructed using Cytoscape software. Auto Duck Vina software was used for molecular docking to analyze the affinity of the key ingredients and the main targets.

Results

We screened 17 active ingredients and 27 key targets of Epimedii Herba in the treatment of Alzheimer’s disease, which were related to the HIF-1 signaling pathway, TNF signaling pathway, PI3K-Akt signaling pathway, NF-κB signaling pathway, VEGF signaling pathway, and sphingolipid signaling pathway.

Conclusions

Based on network pharmacology, the multi-component, multi-target, and multi-pathway characteristics of Epimedii Herba in the treatment of Alzheimer’s disease were explored. Our results provide new ideas for future pharmacological and experimental research on Epimedii Herba in the treatment of Alzheimer’s disease.

Phosphorylation of the N-terminal domain of ribosomal P-stalk protein uL10 governs its association with the ribosome

The uL10 protein is the main constituent of the ribosomal P-stalk, anchoring the whole stalk to the ribosome through interactions with rRNA. The P-stalk is the core of the GTPase-associated center (GAC), a critical element for ribosome biogenesis and ribosome translational activity. All P-stalk proteins (uL10, P1, and P2) undergo phosphorylation within their C termini. Here, we show that uL10 has multiple phosphorylation sites, mapped also within the N-terminal rRNA-binding domain. Our results reveal that the introduction of a negative charge within the N terminus of uL10 impairs its association with the ribosome. These findings demonstrate that uL10 N-terminal phosphorylation has regulatory potential governing the uL10 interaction with the ribosome and may control the activity of GAC.

Using phosphoproteomics data to understand cellular signaling: a comprehensive guide to bioinformatics resources

Mass spectrometry-based phosphoproteomics is becoming an essential methodology for the study of global cellular signaling. Numerous bioinformatics resources are available to facilitate the translation of phosphopeptide identification and quantification results into novel biological and clinical insights, a critical step in phosphoproteomics data analysis. These resources include knowledge bases of kinases and phosphatases, phosphorylation sites, kinase inhibitors, and sequence variants affecting kinase function, and bioinformatics tools that can predict phosphorylation sites in addition to the kinase that phosphorylates them, infer kinase activity, and predict the effect of mutations on kinase signaling. However, these resources exist in silos and it is challenging to select among multiple resources with similar functions. Therefore, we put together a comprehensive collection of resources related to phosphoproteomics data interpretation, compared the use of tools with similar functions, and assessed the usability from the standpoint of typical biologists or clinicians. Overall, tools could be improved by standardization of enzyme names, flexibility of data input and output format, consistent maintenance, and detailed manuals.

Distinct pseudokinase domain conformations underlie divergent activation mechanisms among vertebrate MLKL orthologues

The MLKL pseudokinase is the terminal effector in the necroptosis cell death pathway. Phosphorylation by its upstream regulator, RIPK3, triggers MLKL's conversion from a dormant cytoplasmic protein into oligomers that translocate to, and permeabilize, the plasma membrane to kill cells. The precise mechanisms underlying these processes are incompletely understood, and were proposed to differ between mouse and human cells. Here, we examine the divergence of activation mechanisms among nine vertebrate MLKL orthologues, revealing remarkable specificity of mouse and human RIPK3 for MLKL orthologues. Pig MLKL can restore necroptotic signaling in human cells; while horse and pig, but not rat, MLKL can reconstitute the mouse pathway. This selectivity can be rationalized from the distinct conformations observed in the crystal structures of horse and rat MLKL pseudokinase domains. These studies identify important differences in necroptotic signaling between species, and suggest that, more broadly, divergent regulatory mechanisms may exist among orthologous pseudoenzymes.

Ten quick tips for homology modeling of high-resolution protein 3D structures

The purpose of this quick guide is to help new modelers who have little or no background in comparative modeling yet are keen to produce high-resolution protein 3D structures for their study by following systematic good modeling practices, using affordable personal computers or online computational resources. Through the available experimental 3D-structure repositories, the modeler should be able to access and use the atomic coordinates for building homology models. We also aim to provide the modeler with a rationale behind making a simple list of atomic coordinates suitable for computational analysis abiding to principles of physics (e.g., molecular mechanics). Keeping that objective in mind, these quick tips cover the process of homology modeling and some postmodeling computations such as molecular docking and molecular dynamics (MD). A brief section was left for modeling nonprotein molecules, and a short case study of homology modeling is discussed.

Posttranslational Modifications Mediate the Structural Diversity of Tauopathy Strains

Tau aggregation into insoluble filaments is the defining pathological hallmark of tauopathies. However, it is not known what controls the formation and templated seeding of strain-specific structures associated with individual tauopathies. Here, we use cryo-electron microscopy (cryo-EM) to determine the structures of tau filaments from corticobasal degeneration (CBD) human brain tissue. Cryo-EM and mass spectrometry of tau filaments from CBD reveal that this conformer is heavily decorated with posttranslational modifications (PTMs), enabling us to map PTMs directly onto the structures. By comparing the structures and PTMs of tau filaments from CBD and Alzheimer's disease, it is found that ubiquitination of tau can mediate inter-protofilament interfaces. We propose a structure-based model in which cross-talk between PTMs influences tau filament structure, contributing to the structural diversity of tauopathy strains. Our approach establishes a framework for further elucidating the relationship between the structures of polymorphic fibrils, including their PTMs, and neurodegenerative disease.