Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91

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.

Comparison of force fields to study the zinc-finger containing protein NPL4, a target for disulfiram in cancer therapy

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.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins

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.

Directed Evolution of (R)-2-Hydroxyglutarate Dehydrogenase Improves 2-Oxoadipate Reduction by 2 Orders of Magnitude

Pathway engineering is commonly employed to improve the production of various metabolites but may incur in bottlenecks due to the low catalytic activity of a particular reaction step. The reduction of 2-oxoadipate to (R)-2-hydroxyadipate is a key reaction in metabolic pathways that exploit 2-oxoadipate conversion via α-reduction to produce adipic acid, an industrially important platform chemical. Here, we engineered (R)-2-hydroxyglutarate dehydrogenase from Acidaminococcus fermentans (Hgdh) with the aim of improving 2-oxoadipate reduction. Using a combination of computational analysis, saturation mutagenesis, and random mutagenesis, three mutant variants with a 100-fold higher catalytic efficiency were obtained. As revealed by rational analysis of the mutations found in the variants, this improvement could be ascribed to a general synergistic effect where mutation A206V played a key role since it boosted the enzyme’s activity by 4.8-fold. The Hgdh variants with increased activity toward 2-oxoadipate generated within this study pave the way for the bio-based production of adipic acid.

Unraveling membrane properties at the organelle-level with LipidDyn

•

Lipidomics of organelles could be used to design models for molecular simulations.

•

The bottleneck is the analysis and rationalization of the data from simulations.

•

LipidDyn is an automated pipeline to streamline the analyses of lipid bilayers.

•

LipidDyn allows to collect analysis in a non-time-consuming and reproducible manner.

•

We applied LipidDyn to different case studies to illustrate its potential.

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.

DisProt in 2022: improved quality and accessibility of protein intrinsic disorder annotation

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 conformational and mutational landscape of the ubiquitin-like marker for autophagosome formation in cancer

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.

Characterization of a natural variant of human NDP52 and its functional consequences on mitophagy

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.

Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs

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.

TFG binds LC3C to regulate ULK1 localization and autophagosome formation

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.

Structure and Dynamics in the ATG8 Family From Experimental to Computational Techniques

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.

DisProt: intrinsic protein disorder annotation in 2020

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.

Structure-function investigation of 3-methylaspartate ammonia lyase reveals substrate molecular determinants for the deamination reaction

The enzymatic reactions leading to the deamination of β-lysine, lysine, or 2-aminoadipic acid are of great interest for the metabolic conversion of lysine to adipic acid. Enzymes able to carry out these reactions are not known, however ammonia lyases (EC 4.3.1.-) perform deamination on a wide range of substrates. We have studied 3-methylaspartate ammonia lyase (MAL, EC 4.3.1.2) as a potential candidate for protein engineering to enable deamination towards β-lysine, that we have shown to be a competitive inhibitor of MAL. We have characterized MAL activity, binding and inhibition properties on six different compounds that would allow to define the molecular determinants necessary for MAL to deaminate our substrate of interest. Docking calculations showed that β-lysine as well as the other compounds investigated could fit spatially into MAL catalytic pocket, although they probably are weak or very transient binders and we identified molecular determinants involved in the binding of the substrate. The hydrophobic interactions formed by the methyl group of 3-methylaspartic acid, together with the presence of the amino group on carbon 2, play an essential role in the appropriate binding of the substrate. The results showed that β-lysine is able to fit and bind in MAL catalytic pocket and can be potentially converted from inhibitor to substrate of MAL upon enzyme engineering. The characterization of the binding and inhibition properties of the substrates tested here provide the foundation for future and more extensive studies on engineering MAL that could lead to a MAL variant able to catalyse this challenging deamination reaction.

Conformational gating in ammonia lyases

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.

S-nitrosylation affects TRAP1 structure and ATPase activity and modulates cell response to apoptotic stimuli

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.

Alterations of the interactome of Bcl-2 proteins in breast cancer at the transcriptional, mutational and structural level

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.

Selective autophagy maintains centrosome integrity and accurate mitosis by turnover of centriolar satellites

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.

Bi-allelic mutations in MYL1 cause a severe congenital myopathy

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.

Release of transcriptional repression via ErbB2-induced, SUMO-directed phosphorylation of myeloid zinc finger-1 serine 27 activates lysosome redistribution and invasion

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.

Analyzing Biomolecular Ensembles

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.

HUWE1 E3 ligase promotes PINK1/PARKIN-independent mitophagy by regulating AMBRA1 activation via IKKα

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.

Mitophagy is crucial for mitochondrial quality control and maintenance of cellular homeostasis. Here the authors identify an E3 ubiquitin ligase, HUWE1, that collaborates with LC3-interacting protein AMBRA1 to induce mitochondrial clearance.

Absence of Neurofibromin Induces an Oncogenic Metabolic Switch via Mitochondrial ERK-Mediated Phosphorylation of the Chaperone TRAP1

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.

Catalytic Mechanism of Fungal Lytic Polysaccharide Monooxygenases Investigated by First-Principles Calculations

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.

An optimal distance cutoff for contact-based Protein Structure Networks using side-chain centers of mass

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.

A phosphorylation-motif for tuneable helix stabilisation in intrinsically disordered proteins - Lessons from the sodium proton exchanger 1 (NHE1)

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.

Hybrid Voronoi diagrams, their computation and reduction for applications in computational biochemistry

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.

The Mutational Landscape of the Oncogenic MZF1 SCAN Domain in Cancer

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.

DNA-binding protects p53 from interactions with cofactors involved in transcription-independent functions

Binding-induced conformational changes of a protein at regions distant from the binding site may play crucial roles in protein function and regulation. The p53 tumour suppressor is an example of such an allosterically regulated protein. Little is known, however, about how DNA binding can affect distal sites for transcription factors. Furthermore, the molecular details of how a local perturbation is transmitted through a protein structure are generally elusive and occur on timescales hard to explore by simulations. Thus, we employed state-of-the-art enhanced sampling atomistic simulations to unveil DNA-induced effects on p53 structure and dynamics that modulate the recruitment of cofactors and the impact of phosphorylation at Ser215. We show that DNA interaction promotes a conformational change in a region 3 nm away from the DNA binding site. Specifically, binding to DNA increases the population of an occluded minor state at this distal site by more than 4-fold, whereas phosphorylation traps the protein in its major state. In the minor conformation, the interface of p53 that binds biological partners related to p53 transcription-independent functions is not accessible. Significantly, our study reveals a mechanism of DNA-mediated protection of p53 from interactions with partners involved in the p53 transcription-independent signalling. This also suggests that conformational dynamics is tightly related to p53 signalling.

The dynamics of the β-propeller domain in Kelch protein KLHL40 changes upon nemaline myopathy-associated mutation

The nemaline myopathy-associated E528K mutation in the KLHL40 alters the communication between the Kelch propeller blades.

E2 superfamily of ubiquitin-conjugating enzymes: constitutively active or activated through phosphorylation in the catalytic cleft

Protein phosphorylation is a modification that offers a dynamic and reversible mechanism to regulate the majority of cellular processes. Numerous diseases are associated with aberrant regulation of phosphorylation-induced switches. Phosphorylation is emerging as a mechanism to modulate ubiquitination by regulating key enzymes in this pathway. The molecular mechanisms underpinning how phosphorylation regulates ubiquitinating enzymes, however, are elusive. Here, we show the high conservation of a functional site in E2 ubiquitin-conjugating enzymes. In catalytically active E2s, this site contains aspartate or a phosphorylatable serine and we refer to it as the conserved E2 serine/aspartate (CES/D) site. Molecular simulations of substrate-bound and -unbound forms of wild type, mutant and phosphorylated E2s, provide atomistic insight into the role of the CES/D residue for optimal E2 activity. Both the size and charge of the side group at the site play a central role in aligning the substrate lysine toward E2 catalytic cysteine to control ubiquitination efficiency. The CES/D site contributes to the fingerprint of the E2 superfamily. We propose that E2 enzymes can be divided into constitutively active or regulated families. E2s characterized by an aspartate at the CES/D site signify constitutively active E2s, whereas those containing a serine can be regulated by phosphorylation.

PyInteraph: a framework for the analysis of interaction networks in structural ensembles of proteins

In the last years, a growing interest has been gathering around the ability of Molecular Dynamics (MD) to provide insight into the paths of long-range structural communication in biomolecules. The knowledge of the mechanisms related to structural communication helps in the rationalization in atomistic details of the effects induced by mutations, ligand binding, and the intrinsic dynamics of proteins. We here present PyInteraph, a tool for the analysis of structural ensembles inspired by graph theory. PyInteraph is a software suite designed to analyze MD and structural ensembles with attention to binary interactions between residues, such as hydrogen bonds, salt bridges, and hydrophobic interactions. PyInteraph also allows the different classes of intra- and intermolecular interactions to be represented, combined or alone, in the form of interaction graphs, along with performing network analysis on the resulting interaction graphs. The program also integrates the network description with a knowledge-based force field to estimate the interaction energies between side chains in the protein. It can be used alone or together with the recently developed xPyder PyMOL plugin through an xPyder-compatible format. The software capabilities and associated protocols are here illustrated by biologically relevant cases of study. The program is available free of charge as Open Source software via the GPL v3 license at http://linux.btbs.unimib.it/pyinteraph/.

A comparative study of Whi5 and retinoblastoma proteins: from sequence and structure analysis to intracellular networks

Cell growth and proliferation require a complex series of tight-regulated and well-orchestrated events. Accordingly, proteins governing such events are evolutionary conserved, even among distant organisms. By contrast, it is more singular the case of “core functions” exerted by functional analogous proteins that are not homologous and do not share any kind of structural similarity. This is the case of proteins regulating the G1/S transition in higher eukaryotes–i.e., the retinoblastoma (Rb) tumor suppressor Rb—and budding yeast, i.e., Whi5. The interaction landscape of Rb and Whi5 is quite large, with more than one hundred proteins interacting either genetically or physically with each protein. The Whi5 interactome has been used to construct a concept map of Whi5 function and regulation. Comparison of physical and genetic interactors of Rb and Whi5 allows highlighting a significant core of conserved, common functionalities associated with the interactors indicating that structure and function of the network—rather than individual proteins—are conserved during evolution. A combined bioinformatics and biochemical approach has shown that the whole Whi5 protein is highly disordered, except for a small region containing the protein family signature. The comparison with Whi5 homologs from Saccharomycetales has prompted the hypothesis of a modular organization of structural disorder, with most evolutionary conserved regions alternating with highly variable ones. The finding of a consensus sequence points to the conservation of a specific phosphorylation rhythm along with two disordered sequence motifs, probably acting as phosphorylation-dependent seeds in Whi5 folding/unfolding. Thus, the widely disordered Whi5 appears to act as a hierarchical, “date hub” that has evolutionary assayed an original way of modular organization before being supplanted by the globular, multi-domain structured Rb, more suitable to cover the role of a “party hub”.

Intramolecular interactions stabilizing compact conformations of the intrinsically disordered kinase-inhibitor domain of Sic1: a molecular dynamics investigation

Cyclin-dependent kinase inhibitors (CKIs) are key regulatory proteins of the eukaryotic cell cycle, which modulate cyclin-dependent kinase (Cdk) activity. CKIs perform their inhibitory effect by the formation of ternary complexes with a target kinase and its cognate cyclin. These regulators generally belong to the class of intrinsically disordered proteins (IDPs), which lack a well-defined and organized three-dimensional (3D) structure in their free state, undergoing folding upon binding to specific partners. Unbound IDPs are not merely random-coil structures, but can present intrinsically folded structural units (IFSUs) and collapsed conformations. These structural features can be relevant to protein function in vivo. The yeast CKI Sic1 is a 284-amino acid IDP that binds to Cdk1 in complex with the Clb5,6 cyclins, preventing phosphorylation of G1 substrates and, therefore, entrance to the S phase. Sic1 degradation, triggered by multiple phosphorylation events, promotes cell-cycle progression. Previous experimental studies pointed out a propensity of Sic1 and its isolated domains to populate both extended and compact conformations. The present contribution provides models for compact conformations of the Sic1 kinase-inhibitory domain (KID) by all-atom molecular dynamics (MD) simulations in explicit solvent and in the absence of interactors. The results are integrated by spectroscopic and spectrometric data. Helical IFSUs are identified, along with networks of intramolecular interactions. The results identify a group of putative hub residues and networks of electrostatic interactions, which are likely to be involved in the stabilization of the globular states.