Search

Lysosomal acid sphingomyelinase (ASM), a critical enzyme in lipid metabolism encoded by the SMPD1 gene, plays a crucial role in sphingomyelin hydrolysis in lysosomes. ASM deficiency leads to acid sphingomyelinase deficiency, a rare genetic disorder with diverse clinical manifestations, and the protein can be found mutated in other diseases. We employed a structure-based framework to comprehensively understand the functional implications of ASM variants, integrating pathogenicity predictions with molecular insights derived from a molecular dynamics simulation in a lysosomal membrane environment. Our analysis, encompassing over 400 variants, establishes a structural atlas of missense variants of lysosomal ASM, associating mechanistic indicators with pathogenic potential. Our study highlights variants that influence structural stability or exert local and long-range effects at functional sites. To validate our predictions, we compared them to available experimental data on residual catalytic activity in 135 ASM variants. Notably, our findings also suggest applications of the resulting data for identifying cases suited for enzyme replacement therapy. This comprehensive approach enhances the understanding of ASM variants and provides valuable insights for potential therapeutic interventions., Competing Interests: Declaration of competing interest None., (Copyright © 2024. Published by Elsevier B.V.)

Lysosomes are pivotal in cellular functions and disease, influencing cancer progression and therapy resistance with Acid Sphingomyelinase (ASM) governing their membrane integrity. Moreover, cation amphiphilic drugs (CADs) are known as ASM inhibitors and have anti-cancer activity, but the structural mechanisms of their interactions with the lysosomal membrane and ASM are poorly explored. Our study, leveraging all-atom explicit solvent molecular dynamics simulations, delves into the interaction of glycosylated ASM with the lysosomal membrane and the effects of CAD representatives, i.e., ebastine, hydroxyebastine and loratadine, on the membrane and ASM. Our results confirm the ASM association to the membrane through the saposin domain, previously only shown with coarse-grained models. Furthermore, we elucidated the role of specific residues and ASM-induced membrane curvature in lipid recruitment and orientation. CADs also interfere with the association of ASM with the membrane at the level of a loop in the catalytic domain engaging in membrane interactions. Our computational approach, applicable to various CADs or membrane compositions, provides insights into ASM and CAD interaction with the membrane, offering a valuable tool for future studies., Competing Interests: None., (© 2024 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations., Competing Interests: Declaration of interests M. Trauelsen, T.M.F., and T.W.S. are co-founders of SOLID Therapeutics, and T.W.S. is a co-founder of Embark Biotech., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Molecular dynamics (MD) simulations are a powerful approach to studying the structure and dynamics of proteins related to health and disease. Advances in the MD field allow modeling proteins with high accuracy. However, modeling metal ions and their interactions with proteins is still challenging. NPL4 is a zinc-binding protein and works as a cofactor for p97 to regulate protein homeostasis. NPL4 is of biomedical importance and has been proposed as the target of disulfiram, a drug recently repurposed for cancer treatment. Experimental studies proposed that the disulfiram metabolites, bis-(diethyldithiocarbamate)‑copper and cupric ions, induce NPL4 misfolding and aggregation. However, the molecular details of their interactions with NPL4 and consequent structural effects are still elusive. Here, biomolecular simulations can help to shed light on the related structural details. To apply MD simulations to NPL4 and its interaction with copper the first important step is identifying a suitable force field to describe the protein in its zinc-bound states. We examined different sets of non-bonded parameters because we want to study the misfolding mechanism and cannot rule out that the zinc may detach from the protein during the process and copper replaces it. We investigated the force-field ability to model the coordination geometry of the metal ions by comparing the results from MD simulations with optimized geometries from quantum mechanics (QM) calculations using model systems of NPL4. Furthermore, we investigated the performance of a force field including bonded parameters to treat copper ions in NPL4 that we obtained based on QM calculations., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field., (© 2023. The Author(s).)

The tumor protein 53 (p53) is involved in transcription-dependent and independent processes. Several p53 variants related to cancer have been found to impact protein stability. Other variants, on the contrary, might have little impact on structural stability and have local or long-range effects on the p53 interactome. Our group previously identified a loop in the DNA binding domain (DBD) of p53 (residues 207-213) which can recruit different interactors. Experimental structures of p53 in complex with other proteins strengthen the importance of this interface for protein-protein interactions. We here characterized with structure-based approaches somatic and germline variants of p53 which could have a marginal effect in terms of stability and act locally or allosterically on the region 207-213 with consequences on the cytosolic functions of this protein. To this goal, we studied 1132 variants in the p53 DBD with structure-based approaches, accounting also for protein dynamics. We focused on variants predicted with marginal effects on structural stability. We then investigated each of these variants for their impact on DNA binding, dimerization of the p53 DBD, and intramolecular contacts with the 207-213 region. Furthermore, we identified variants that could modulate long-range the conformation of the region 207-213 using a coarse-grain model for allostery and all-atom molecular dynamics simulations. Our predictions have been further validated using enhanced sampling methods for 15 variants. The methodologies used in this study could be more broadly applied to other p53 variants or cases where conformational changes of loop regions are essential in the function of disease-related proteins., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Cellular membranes are formed from different lipids in various amounts and proportions depending on the subcellular localization. The lipid composition of membranes is sensitive to changes in the cellular environment, and its alterations are linked to several diseases. Lipids not only form lipid-lipid interactions but also interact with other biomolecules, including proteins. Molecular dynamics (MD) simulations are a powerful tool to study the properties of cellular membranes and membrane-protein interactions on different timescales and resolutions. Over the last few years, software and hardware for biomolecular simulations have been optimized to routinely run long simulations of large and complex biological systems. On the other hand, high-throughput techniques based on lipidomics provide accurate estimates of the composition of cellular membranes at the level of subcellular compartments. Lipidomic data can be analyzed to design biologically relevant models of membranes for MD simulations. Similar applications easily result in a massive amount of simulation data where the bottleneck becomes the analysis of the data. In this context, we developed LipidDyn , a Python-based pipeline to streamline the analyses of MD simulations of membranes of different compositions. Once the simulations are collected, LipidDyn provides average properties and time series for several membrane properties such as area per lipid, thickness, order parameters, diffusion motions, lipid density, and lipid enrichment/depletion. The calculations exploit parallelization, and the pipeline includes graphical outputs in a publication-ready form. We applied LipidDyn to different case studies to illustrate its potential, including membranes from cellular compartments and transmembrane protein domains. LipidDyn is available free of charge under the GNU General Public License from https://github.com/ELELAB/LipidDyn., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.)

The Database of Intrinsically Disordered Proteins (DisProt, URL: https://disprot.org) is the major repository of manually curated annotations of intrinsically disordered proteins and regions from the literature. We report here recent updates of DisProt version 9, including a restyled web interface, refactored Intrinsically Disordered Proteins Ontology (IDPO), improvements in the curation process and significant content growth of around 30%. Higher quality and consistency of annotations is provided by a newly implemented reviewing process and training of curators. The increased curation capacity is fostered by the integration of DisProt with APICURON, a dedicated resource for the proper attribution and recognition of biocuration efforts. Better interoperability is provided through the adoption of the Minimum Information About Disorder (MIADE) standard, an active collaboration with the Gene Ontology (GO) and Evidence and Conclusion Ontology (ECO) consortia and the support of the ELIXIR infrastructure., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Macroautophagy/autophagy is a cellular process to recycle damaged cellular components, and its modulation can be exploited for disease treatments. A key autophagy player is the ubiquitin-like protein MAP1LC3B/LC3B. Mutations and changes in MAP1LC3B expression occur in cancer samples. However, the investigation of the effects of these mutations on MAP1LC3B protein structure is still missing. Despite many LC3B structures that have been solved, a comprehensive study, including dynamics, has not yet been undertaken. To address this knowledge gap, we assessed nine physical models for biomolecular simulations for their capabilities to describe the structural ensemble of MAP1LC3B. With the resulting MAP1LC3B structural ensembles, we characterized the impact of 26 missense mutations from pan-cancer studies with different approaches, and we experimentally validated our prediction for six variants using cellular assays. Our findings shed light on damaging or neutral mutations in MAP1LC3B, providing an atlas of its modifications in cancer. In particular, P32Q mutation was found detrimental for protein stability with a propensity to aggregation. In a broader context, our framework can be applied to assess the pathogenicity of protein mutations or to prioritize variants for experimental studies, allowing to comprehensively account for different aspects that mutational events alter in terms of protein structure and function. Abbreviations : ATG: autophagy-related; Cα: alpha carbon; CG: coarse-grained; CHARMM: Chemistry at Harvard macromolecular mechanics; CONAN: contact analysis; FUNDC1: FUN14 domain containing 1; FYCO1: FYVE and coiled-coil domain containing 1; GABARAP: GABA type A receptor-associated protein; GROMACS: Groningen machine for chemical simulations; HP: hydrophobic pocket; LIR: LC3 interacting region; MAP1LC3B/LC3B microtubule associated protein 1 light chain 3 B; MD: molecular dynamics; OPTN: optineurin; OSF: open software foundation; PE: phosphatidylethanolamine, PLEKHM1: pleckstrin homology domain-containing family M 1; PSN: protein structure network; PTM: post-translational modification; SA: structural alphabet; SLiM: short linear motif; SQSTM1/p62: sequestosome 1; WT: wild-type.

The use of lignocellulosic-based fermentation media will be a necessary part of the transition to a circular bio-economy. These media contain many inhibitors to microbial growth, including acetic acid. Under industrially relevant conditions, acetic acid enters the cell predominantly through passive diffusion across the plasma membrane. The lipid composition of the membrane determines the rate of uptake of acetic acid, and thicker, more rigid membranes impede passive diffusion. We hypothesized that the elongation of glycerophospholipid fatty acids would lead to thicker and more rigid membranes, reducing the influx of acetic acid. Molecular dynamics simulations were used to predict the changes in membrane properties. Heterologous expression of Arabidopsis thaliana genes fatty acid elongase 1 (FAE1) and glycerol-3-phosphate acyltransferase 5 (GPAT5) increased the average fatty acid chain length. However, this did not lead to a reduction in the net uptake rate of acetic acid. Despite successful strain engineering, the net uptake rate of acetic acid did not decrease. We suggest that changes in the relative abundance of certain membrane lipid headgroups could mitigate the effect of longer fatty acid chains, resulting in a higher net uptake rate of acetic acid., (© 2021. The Author(s).)

The role of mitophagy, a process that allows the removal of damaged mitochondria from cells, remains unknown in multiple sclerosis (MS), a disease that is found associated with dysfunctional mitochondria. Here we have qualitatively and quantitatively studied the main players in PINK1-mediated mitophagy in peripheral blood mononuclear cells (PBMCs) of patients with relapsing-remitting MS. We found the variant c.491G>A (rs550510, p.G140E) of NDP52, one of the major mitophagy receptor genes, associated with a MS cohort. Through the characterization of this variant, we discovered that the residue 140 of human NDP52 is a crucial modulator of NDP52/LC3C binding, promoting the formation of autophagosomes in order to drive efficient mitophagy. In addition, we found that in the PBMC population, NDP52 is mainly expressed in B cells and by ensuring efficient mitophagy, it is able to limit the production of the proinflammatory cytokine TNF-α following cell stimulation. In sum, our results contribute to a better understanding of the role of NDP52 in mitophagy and underline, for the first time, a possible role of NDP52 in MS., (© 2021. The Author(s), under exclusive licence to ADMC Associazione Differenziamento e Morte Cellulare.)

Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Lambrughi, Maiani, Aykac Fas, Shaw, Kragelund, Lindorff-Larsen, Teilum, Invernizzi and Papaleo.)

The early secretory pathway and autophagy are two essential and evolutionarily conserved endomembrane processes that are finely interlinked. Although growing evidence suggests that intracellular trafficking is important for autophagosome biogenesis, the molecular regulatory network involved is still not fully defined. In this study, we demonstrate a crucial effect of the COPII vesicle-related protein TFG (Trk-fused gene) on ULK1 puncta number and localization during autophagy induction. This, in turn, affects formation of the isolation membrane, as well as the correct dynamics of association between LC3B and early ATG proteins, leading to the proper formation of both omegasomes and autophagosomes. Consistently, fibroblasts derived from a hereditary spastic paraparesis (HSP) patient carrying mutated TFG (R106C) show defects in both autophagy and ULK1 puncta accumulation. In addition, we demonstrate that TFG activity in autophagy depends on its interaction with the ATG8 protein LC3C through a canonical LIR motif, thereby favouring LC3C-ULK1 binding. Altogether, our results uncover a link between TFG and autophagy and identify TFG as a molecular scaffold linking the early secretion pathway to autophagy., (© 2021 The Authors. Published under the terms of the CC BY NC ND 4.0 license.)

PyInteraph is a software package designed for the analysis of structural communication from conformational ensembles, such as those derived from in silico simulations, under the formalism of protein structure networks. We demonstrate its usage for the calculation and analysis of intramolecular interaction networks derived from three different types of interactions, as well as with a more general protocol based on distances between centers of mass. We use the xPyder PyMOL plug-in to visualize such networks on the three-dimensional structure of the protein. We showcase our protocol on a molecular dynamics trajectory of the Cyclophilin A wild-type enzyme, a well-studied protein in which different allosteric mechanisms have been investigated.

Background: Ammonia lyases are enzymes of industrial and biomedical interest. Knowledge of structure-dynamics-function relationship in ammonia lyases is instrumental for exploiting the potential of these enzymes in industrial or biomedical applications., Methods: We investigated the conformational changes in the proximity of the catalytic pocket of a 3-methylaspartate ammonia lyase (MAL) as a model system. At this scope, we used microsecond all-atom molecular dynamics simulations, analyzed with dimensionality reduction techniques, as well as in terms of contact networks and correlated motions., Results: We identify two regulatory elements in the MAL structure, i.e., the β5-α2 loop and the helix-hairpin-loop subdomain. These regulatory elements undergo conformational changes switching from 'occluded' to 'open' states. The rearrangements are coupled to changes in the accessibility of the active site. The β5-α2 loop and the helix-hairpin-loop subdomain modulate the formation of tunnels from the protein surface to the catalytic site, making the active site more accessible to the substrate when they are in an open state., Conclusions: Our work pinpoints a sequential mechanism, in which the helix-hairpin-loop subdomain of MAL needs to break a subset of intramolecular interactions first to favor the displacement of the β5-α2 loop. The coupled conformational changes of these two elements contribute to modulate the accessibility of the catalytic site., General Significance: Similar molecular mechanisms can have broad relevance in other ammonia lyases with similar regulatory loops. Our results also imply that it is important to account for protein dynamics in the design of variants of ammonia lyases for industrial and biomedical applications., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Autophagy is a conserved and essential intracellular mechanism for the removal of damaged components. Since autophagy deregulation is linked to different kinds of pathologies, it is fundamental to gain knowledge on the fine molecular and structural details related to the core proteins of the autophagy machinery. Among these, the family of human ATG8 proteins plays a central role in recruiting other proteins to the different membrane structures involved in the autophagic pathway. Several experimental structures are available for the members of the ATG8 family alone or in complex with their different biological partners, including disordered regions of proteins containing a short linear motif called LC3 interacting motif. Recently, the first structural details of the interaction of ATG8 proteins with biological membranes came into light. The availability of structural data for human ATG8 proteins has been paving the way for studies on their structure-function-dynamic relationship using biomolecular simulations. Experimental and computational structural biology can help to address several outstanding questions on the mechanism of human ATG8 proteins, including their specificity toward different interactors, their association with membranes, the heterogeneity of their conformational ensemble, and their regulation by post-translational modifications. We here summarize the main results collected so far and discuss the future perspectives within the field and the knowledge gaps. Our review can serve as a roadmap for future structural and dynamics studies of the ATG8 family members in health and disease., (Copyright © 2020 Sora, Kumar, Maiani, Lambrughi, Tiberti and Papaleo.)

The mitochondrial chaperone TRAP1 has been involved in several mitochondrial functions, and modulation of its expression/activity has been suggested to play a role in the metabolic reprogramming distinctive of cancer cells. TRAP1 posttranslational modifications, i.e. phosphorylation, can modify its capability to bind to different client proteins and modulate its oncogenic activity. Recently, it has been also demonstrated that TRAP1 is S-nitrosylated at Cys501, a redox modification associated with its degradation via the proteasome. Here we report molecular dynamics simulations of TRAP1, together with analysis of long-range structural communication, providing a model according to which Cys501 S-nitrosylation induces conformational changes to distal sites in the structure of the protein. The modification is also predicted to alter open and closing motions for the chaperone function. By means of colorimetric assays and site directed mutagenesis aimed at generating C501S variant, we also experimentally confirmed that selective S-nitrosylation of Cys501 decreases ATPase activity of recombinant TRAP1. Coherently, C501S mutant was more active and conferred protection to cell death induced by staurosporine. Overall, our results provide the first in silico, in vitro and cellular evidence of the relevance of Cys501 S-nitrosylation in TRAP1 biology., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The Database of Protein Disorder (DisProt, URL: https://disprot.org) provides manually curated annotations of intrinsically disordered proteins from the literature. Here we report recent developments with DisProt (version 8), including the doubling of protein entries, a new disorder ontology, improvements of the annotation format and a completely new website. The website includes a redesigned graphical interface, a better search engine, a clearer API for programmatic access and a new annotation interface that integrates text mining technologies. The new entry format provides a greater flexibility, simplifies maintenance and allows the capture of more information from the literature. The new disorder ontology has been formalized and made interoperable by adopting the OWL format, as well as its structure and term definitions have been improved. The new annotation interface has made the curation process faster and more effective. We recently showed that new DisProt annotations can be effectively used to train and validate disorder predictors. We believe the growth of DisProt will accelerate, contributing to the improvement of function and disorder predictors and therefore to illuminate the 'dark' proteome., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Apoptosis is an essential defensive mechanism against tumorigenesis. Proteins of the B-cell lymphoma-2 (Bcl-2) family regulate programmed cell death by the mitochondrial apoptosis pathway. In response to intracellular stress, the apoptotic balance is governed by interactions of three distinct subgroups of proteins; the activator/sensitizer BH3 (Bcl-2 homology 3)-only proteins, the pro-survival, and the pro-apoptotic executioner proteins. Changes in expression levels, stability, and functional impairment of pro-survival proteins can lead to an imbalance in tissue homeostasis. Their overexpression or hyperactivation can result in oncogenic effects. Pro-survival Bcl-2 family members carry out their function by binding the BH3 short linear motif of pro-apoptotic proteins in a modular way, creating a complex network of protein-protein interactions. Their dysfunction enables cancer cells to evade cell death. The critical role of Bcl-2 proteins in homeostasis and tumorigenesis, coupled with mounting insight in their structural properties, make them therapeutic targets of interest. A better understanding of gene expression, mutational profile, and molecular mechanisms of pro-survival Bcl-2 proteins in different cancer types, could help to clarify their role in cancer development and may guide advancement in drug discovery. Here, we shed light on the pro-survival Bcl-2 proteins in breast cancer using different bioinformatic approaches, linking -omics with structural data. We analyzed the changes in the expression of the Bcl-2 proteins and their BH3-containing interactors in breast cancer samples. We then studied, at the structural level, a selection of interactions, accounting for effects induced by mutations found in the breast cancer samples. We find two complexes between the up-regulated Bcl2A1 and two down-regulated BH3-only candidates (i.e., Hrk and Nr4a1) as targets associated with reduced apoptosis in breast cancer samples for future experimental validation. Furthermore, we predict L99R, M75R as damaging mutations altering protein stability, and Y120C as a possible allosteric mutation from an exposed surface to the BH3-binding site., Competing Interests: The authors have declared that no competing interests exist.

The centrosome is the master orchestrator of mitotic spindle formation and chromosome segregation in animal cells. Centrosome abnormalities are frequently observed in cancer, but little is known of their origin and about pathways affecting centrosome homeostasis. Here we show that autophagy preserves centrosome organization and stability through selective turnover of centriolar satellite components, a process we termed doryphagy. Autophagy targets the satellite organizer PCM1 by interacting with GABARAPs via a C-terminal LIR motif. Accordingly, autophagy deficiency results in accumulation of large abnormal centriolar satellites and a resultant dysregulation of centrosome composition. These alterations have critical impact on centrosome stability and lead to mitotic centrosome fragmentation and unbalanced chromosome segregation. Our findings identify doryphagy as an important centrosome-regulating pathway and bring mechanistic insights to the link between autophagy dysfunction and chromosomal instability. In addition, we highlight the vital role of centriolar satellites in maintaining centrosome integrity.

HER2/ErbB2 activation turns on transcriptional processes that induce local invasion and lead to systemic metastasis. The early transcriptional changes needed for ErbB2-induced invasion are poorly understood. Here, we link ErbB2 activation to invasion via ErbB2-induced, SUMO-directed phosphorylation of a single serine residue, S27, of the transcription factor myeloid zinc finger-1 (MZF1). Utilizing an antibody against MZF1-pS27, we show that the phosphorylation of S27 correlates significantly (p < 0.0001) with high-level expression of ErbB2 in primary invasive breast tumors. Phosphorylation of MZF1-S27 is an early response to ErbB2 activation and results in increased transcriptional activity of MZF1. It is needed for the ErbB2-induced expression of MZF1 target genes CTSB and PRKCA, and invasion of single-cells from ErbB2-expressing breast cancer spheroids. The phosphorylation of MZF1-S27 is preceded by poly-SUMOylation of K23, which can make S27 accessible to efficient phosphorylation by PAK4. Based on our results, we suggest for an activation mechanism where phosphorylation of MZF1-S27 triggers MZF1 dissociation from its transcriptional repressors such as the CCCTC-binding factor (CTCF). Our findings increase understanding of the regulation of invasive signaling in breast cancer by uncovering a detailed biological mechanism of how ErbB2 activation can rapidly lead to its invasion-promoting target gene expression and invasion.

Myoglobin, encoded by MB, is a small cytoplasmic globular hemoprotein highly expressed in cardiac myocytes and oxidative skeletal myofibers. Myoglobin binds O 2, facilitates its intracellular transport and serves as a controller of nitric oxide and reactive oxygen species. Here, we identify a recurrent c.292C>T (p.His98Tyr) substitution in MB in fourteen members of six European families suffering from an autosomal dominant progressive myopathy with highly characteristic sarcoplasmic inclusions in skeletal and cardiac muscle. Myoglobinopathy manifests in adulthood with proximal and axial weakness that progresses to involve distal muscles and causes respiratory and cardiac failure. Biochemical characterization reveals that the mutant myoglobin has altered O 2 binding, exhibits a faster heme dissociation rate and has a lower reduction potential compared to wild-type myoglobin. Preliminary studies show that mutant myoglobin may result in elevated superoxide levels at the cellular level. These data define a recognizable muscle disease associated with MB mutation.

Several techniques are available to generate conformational ensembles of proteins and other biomolecules either experimentally or computationally. These methods produce a large amount of data that need to be analyzed to identify structure-dynamics-function relationship. In this chapter, we will cover different tools to unveil the information hidden in conformational ensemble data and to guide toward the rationalization of the data. We included routinely used approaches such as dimensionality reduction, as well as new methods inspired by high-order statistics and graph theory.

Congenital myopathies are typically characterised by early onset hypotonia, weakness and hallmark features on biopsy. Despite the rapid pace of gene discovery, ∼50% of patients with a congenital myopathy remain without a genetic diagnosis following screening of known disease genes. We performed exome sequencing on two consanguineous probands diagnosed with a congenital myopathy and muscle biopsy showing selective atrophy/hypotrophy or absence of type II myofibres. We identified variants in the gene (MYL1) encoding the skeletal muscle fast-twitch specific myosin essential light chain (ELC) in both probands. A homozygous essential splice acceptor variant (c.479-2A > G, predicted to result in skipping of exon 5 was identified in Proband 1, and a homozygous missense substitution (c.488T>G, p.(Met163Arg)) was identified in Proband 2. Protein modelling of the p.(Met163Arg) substitution predicted it might impede intermolecular interactions that facilitate binding to the IQ domain of myosin heavy chain, thus likely impacting on the structure and functioning of the myosin motor. MYL1 was markedly reduced in skeletal muscle from both probands, suggesting that the missense substitution likely results in an unstable protein. Knock down of myl1 in zebrafish resulted in abnormal morphology, disrupted muscle structure and impaired touch-evoked escape responses, thus confirming that skeletal muscle fast-twitch specific myosin ELC is critical for myofibre development and function. Our data implicate MYL1 as a crucial protein for adequate skeletal muscle function and that MYL1 deficiency is associated with severe congenital myopathy.

The selective removal of undesired or damaged mitochondria by autophagy, known as mitophagy, is crucial for cellular homoeostasis, and prevents tumour diffusion, neurodegeneration and ageing. The pro-autophagic molecule AMBRA1 (autophagy/beclin-1 regulator-1) has been defined as a novel regulator of mitophagy in both PINK1/PARKIN-dependent and -independent systems. Here, we identified the E3 ubiquitin ligase HUWE1 as a key inducing factor in AMBRA1-mediated mitophagy, a process that takes place independently of the main mitophagy receptors. Furthermore, we show that mitophagy function of AMBRA1 is post-translationally controlled, upon HUWE1 activity, by a positive phosphorylation on its serine 1014. This modification is mediated by the IKKα kinase and induces structural changes in AMBRA1, thus promoting its interaction with LC3/GABARAP (mATG8) proteins and its mitophagic activity. Altogether, these results demonstrate that AMBRA1 regulates mitophagy through a novel pathway, in which HUWE1 and IKKα are key factors, shedding new lights on the regulation of mitochondrial quality control and homoeostasis in mammalian cells.

Lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes that facilitate the degradation of recalcitrant polysaccharides by the oxidative cleavage of glycosidic bonds. They are gaining rapidly increasing attention as key players in biomass conversion, especially for the production of second-generation biofuels. Elucidation of the detailed mechanism of the LPMO reaction is a major step toward the assessment and optimization of LPMO efficacy in industrial biotechnology, paving the way to utilization of sustainable fuel sources. Here, we used density functional theory calculations to study the reaction pathways suggested to date, exploiting a very large active-site model for a fungal AA9 LPMO and using a celloheptaose unit as a substrate mimic. We identify a copper oxyl intermediate as being responsible for H-atom abstraction from the substrate, followed by a rapid, water-assisted hydroxyl rebound, leading to substrate hydroxylation.

Intrinsically disordered proteins (IDPs) are involved in many pivotal cellular processes including phosphorylation and signalling. The structural and functional effects of phosphorylation of IDPs remain poorly understood and difficult to predict. Thus, a need exists to identify motifs that confer phosphorylation-dependent perturbation of the local preferences for forming e.g. helical structures as well as motifs that do not. The disordered distal tail of the Na + /H + exchanger 1 (NHE1) is six-times phosphorylated (S693, S723, S726, S771, T779, S785) by the mitogen activated protein kinase 2 (MAPK1, ERK2). Using NMR spectroscopy, we found that two out of those six phosphorylation sites had a stabilizing effect on transient helices. One of these was further investigated by circular dichroism and NMR spectroscopy as well as by molecular dynamic simulations, which confirmed the stabilizing effect and resulted in the identification of a short linear motif for helix stabilisation: [S/T]-P-{3}-[R/K] where [S/T] is the phosphorylation-site. By analysing IDP and phosphorylation site databases we found that the motif is significantly enriched around known phosphorylation sites, supporting a potential wider-spread role in phosphorylation-mediated regulation of intrinsically disordered proteins. The identification of such motifs is important for understanding the molecular mechanism of cellular signalling, and is crucial for the development of predictors for the structural effect of phosphorylation; a tool of relevance for understanding disease-promoting mutations that for example interfere with signalling for instance through constitutive active and often cancer-promoting signalling., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

Proteins are highly dynamic entities attaining a myriad of different conformations. Protein side chains change their states during dynamics, causing clashes that are propagated at distal sites. A convenient formalism to analyze protein dynamics is based on network theory using Protein Structure Networks (PSNs). Despite their broad applicability, few efforts have been devoted to benchmarking PSN methods and to provide the community with best practices. In many applications, it is convenient to use the centers of mass of the side chains as nodes. It becomes thus critical to evaluate the minimal distance cutoff between the centers of mass which will provide stable network properties. Moreover, when the PSN is derived from a structural ensemble collected by molecular dynamics (MD), the impact of the MD force field has to be evaluated. We selected a dataset of proteins with different fold and size and assessed the two fundamental properties of the PSN, i.e. hubs and connected components. We identified an optimal cutoff of 5 Å that is robust to changes in the force field and the proteins. Our study builds solid foundations for the harmonization and standardization of the PSN approach.

Geometric models of molecular structures are often described as a set of balls, where balls represent individual atoms. The ability to describe and explore the empty space among these balls is important, e.g., in the analysis of the interaction of enzymes with substrates, ligands and solvent molecules. Voronoi diagrams from the field of computational geometry are often used here, because they provide a mathematical description of how the whole space can be divided into regions assigned to individual atoms. This paper introduces a combination of two different types of Voronoi diagrams into a new hybrid Voronoi diagram - one part of this diagram belongs to the additively weighted (aw-Voronoi) diagram and the other to the power diagram. The boundary between them is controlled by a user-defined constant (the probe radius). Both parts are computed by different algorithms, which are already known. The reduced aw-Voronoi diagram is then obtained by removing the power diagram part from the hybrid diagram. Reduced aw-Voronoi diagrams are perfectly tailored for the analysis of dynamic molecular structures, their computation is faster and storage requirements are lower than in the case of complete aw-Voronoi diagrams. Here, we showed their application to key proteins in cancer research such as p53 and ARID proteins as case study. We identified a biologically relevant cavity in p53 structural ensembles generated by molecular dynamics simulations and analyzed its accessibility, attesting the potential of our approach. This method is relevant for cancer research since it permits to depict a dynamical view of cavities and pockets in proteins that could be affected by mutations in the disease. Our approach opens novel prospects for the study of cancer-related proteins by molecular simulations and the identification of novel targets for drug design., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Mutations in neurofibromin, a Ras GTPase-activating protein, lead to the tumor predisposition syndrome neurofibromatosis type 1. Here, we report that cells lacking neurofibromin exhibit enhanced glycolysis and decreased respiration in a Ras/ERK-dependent way. In the mitochondrial matrix of neurofibromin-deficient cells, a fraction of active ERK1/2 associates with succinate dehydrogenase (SDH) and TRAP1, a chaperone that promotes the accumulation of the oncometabolite succinate by inhibiting SDH. ERK1/2 enhances both formation of this multimeric complex and SDH inhibition. ERK1/2 kinase activity is favored by the interaction with TRAP1, and TRAP1 is, in turn, phosphorylated in an ERK1/2-dependent way. TRAP1 silencing or mutagenesis at the serine residues targeted by ERK1/2 abrogates tumorigenicity, a phenotype that is reverted by addition of a cell-permeable succinate analog. Our findings reveal that Ras/ERK signaling controls the metabolic changes orchestrated by TRAP1 that have a key role in tumor growth and are a promising target for anti-neoplastic strategies., (Copyright © 2017 The Author(s). Published by Elsevier Inc. All rights reserved.)

SCAN domains in zinc-finger transcription factors are crucial mediators of protein-protein interactions. Up to 240 SCAN-domain encoding genes have been identified throughout the human genome. These include cancer-related genes, such as the myeloid zinc finger 1 ( MZF1 ), an oncogenic transcription factor involved in the progression of many solid cancers. The mechanisms by which SCAN homo- and heterodimers assemble and how they alter the transcriptional activity of zinc-finger transcription factors in cancer and other diseases remain to be investigated. Here, we provide the first description of the conformational ensemble of the MZF1 SCAN domain cross-validated against NMR experimental data, which are probes of structure and dynamics on different timescales. We investigated the protein-protein interaction network of MZF1 and how it is perturbed in different cancer types by the analyses of high-throughput proteomics and RNASeq data. Collectively, we integrated many computational approaches, ranging from simple empirical energy functions to all-atom microsecond molecular dynamics simulations and network analyses to unravel the effects of cancer-related substitutions in relation to MZF1 structure and interactions.