Conformational ensemble of human α-synuclein physiological form predicted by molecular simulations

SPEADI: Accelerated Analysis of IDP-Ion Interactions from MD-Trajectories

The disordered nature of Intrinsically Disordered Proteins (IDPs) makes their structural ensembles particularly susceptible to changes in chemical environmental conditions, often leading to an alteration of their normal functions. A Radial Distribution Function (RDF) is considered a standard method for characterizing the chemical environment surrounding particles during atomistic simulations, commonly averaged over an entire or part of a trajectory. Given their high structural variability, such averaged information might not be reliable for IDPs. We introduce the Time-Resolved Radial Distribution Function (TRRDF), implemented in our open-source Python package SPEADI, which is able to characterize dynamic environments around IDPs. We use SPEADI to characterize the dynamic distribution of ions around the IDPs Alpha-Synuclein (AS) and Humanin (HN) from Molecular Dynamics (MD) simulations, and some of their selected mutants, showing that local ion-residue interactions play an important role in the structures and behaviors of IDPs.

How important is the N-terminal acetylation of alpha-synuclein for its function and aggregation into amyloids?

N-α-acetylation is a frequently occurring post-translational modification in eukaryotic proteins. It has manifold physiological consequences on the regulation and function of several proteins, with emerging studies suggesting that it is a global regulator of stress responses. For decades, in vitro biochemical investigations into the precise role of the intrinsically disordered protein alpha-synuclein (αS) in the etiology of Parkinson's disease (PD) were performed using non-acetylated αS. The N-terminus of α-synuclein is now unequivocally known to be acetylated in vivo, however, there are many aspects of this post-translational modifications that are not understood well. Is N-α-acetylation of αS a constitutive modification akin to most cellular proteins, or is it spatio-temporally regulated? Is N-α-acetylation of αS relevant to the as yet elusive function of αS? How does the N-α-acetylation of αS influence the aggregation of αS into amyloids? Here, we provide an overview of the current knowledge and discuss prevailing hypotheses on the impact of N-α-acetylation of αS on its conformational, oligomeric, and fibrillar states. The extent to which N-α-acetylation of αS is vital for its function, membrane binding, and aggregation into amyloids is also explored here. We further discuss the overall significance of N-α-acetylation of αS for its functional and pathogenic implications in Lewy body formation and synucleinopathies.

α-Synuclein phosphorylation at serine 129 occurs after initial protein deposition and inhibits seeded fibril formation and toxicity

SignificanceConverging evidence points to the build-up of phosphorylated α-synuclein (α-syn) at residue serine 129 (pS129) in Lewy body disease, suggesting its central role in the regulation of α-syn aggregation and neuronal degeneration. However, a comprehensive understanding of the role of α-syn phosphorylation at pS129 in α-synuclenopathies pathogenesis is still lacking. Herein, we study the phosphorylation incidence and its effect on α-syn aggregation propensity and cellular toxicity. Collectively, our data suggest that pS129 occurred subsequent to initial α-syn aggregation, lessened aggregation propensity, and attenuated cytotoxicity through diverse assays. Our findings highlight major implications for a better understanding of the role of a molecular modification on protein aggregation.

Passive Immunization in Alpha-Synuclein Preclinical Animal Models

Alpha-synucleinopathies include Parkinson’s disease, dementia with Lewy bodies, pure autonomic failure and multiple system atrophy. These are all progressive neurodegenerative diseases that are characterized by pathological misfolding and accumulation of the protein alpha-synuclein (αsyn) in neurons, axons or glial cells in the brain, but also in other organs. The abnormal accumulation and propagation of pathogenic αsyn across the autonomic connectome is associated with progressive loss of neurons in the brain and peripheral organs, resulting in motor and non-motor symptoms. To date, no cure is available for synucleinopathies, and therapy is limited to symptomatic treatment of motor and non-motor symptoms upon diagnosis. Recent advances using passive immunization that target different αsyn structures show great potential to block disease progression in rodent studies of synucleinopathies. However, passive immunotherapy in clinical trials has been proven safe but less effective than in preclinical conditions. Here we review current achievements of passive immunotherapy in animal models of synucleinopathies. Furthermore, we propose new research strategies to increase translational outcome in patient studies, (1) by using antibodies against immature conformations of pathogenic αsyn (monomers, post-translationally modified monomers, oligomers and protofibrils) and (2) by focusing treatment on body-first synucleinopathies where damage in the brain is still limited and effective immunization could potentially stop disease progression by blocking the spread of pathogenic αsyn from peripheral organs to the brain.

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.

The structural heterogeneity of α-synuclein is governed by several distinct subpopulations with interconversion times slower than milliseconds

α-Synuclein plays an important role in synaptic functions by interacting with synaptic vesicle membrane, while its oligomers and fibrils are associated with several neurodegenerative diseases. The specific monomer structures that promote its membrane binding and self-association remain elusive due to its transient nature as an intrinsically disordered protein. Here, we use inter-dye distance distributions from bulk time-resolved Förster resonance energy transfer as restraints in discrete molecular dynamics simulations to map the conformational space of the α-synuclein monomer. We further confirm the generated conformational ensemble in orthogonal experiments utilizing far-UV circular dichroism and cross-linking mass spectrometry. Single-molecule protein-induced fluorescence enhancement measurements show that within this conformational ensemble, some of the conformations of α-synuclein are surprisingly stable, exhibiting conformational transitions slower than milliseconds. Our comprehensive analysis of the conformational ensemble reveals essential structural properties and potential conformations that promote its various functions in membrane interaction or oligomer and fibril formation.

Amyloid Oligomers: A Joint Experimental/Computational Perspective on Alzheimer's Disease, Parkinson's Disease, Type II Diabetes, and Amyotrophic Lateral Sclerosis

Protein misfolding and aggregation is observed in many amyloidogenic diseases affecting either the central nervous system or a variety of peripheral tissues. Structural and dynamic characterization of all species along the pathways from monomers to fibrils is challenging by experimental and computational means because they involve intrinsically disordered proteins in most diseases. Yet understanding how amyloid species become toxic is the challenge in developing a treatment for these diseases. Here we review what computer, in vitro, in vivo, and pharmacological experiments tell us about the accumulation and deposition of the oligomers of the (Aβ, tau), α-synuclein, IAPP, and superoxide dismutase 1 proteins, which have been the mainstream concept underlying Alzheimer's disease (AD), Parkinson's disease (PD), type II diabetes (T2D), and amyotrophic lateral sclerosis (ALS) research, respectively, for many years.

Albumin Alters the Conformational Ensemble of Amyloid-β by Promiscuous Interactions: Implications for Amyloid Inhibition

Human serum albumin (HSA) is a key endogenous inhibitor of amyloid-β (Αβ) aggregation. In vitro HSA inhibits Aβ fibrillization and targets multiple species along the aggregation pathway including monomers, oligomers, and protofibrils. Amyloid inhibition by HSA has both pathological implications and therapeutic potential, but the underlying molecular mechanism remains elusive. As a first step towards addressing this complex question, we studied the interactions of an Aβ42 monomer with HSA by molecular dynamics simulations. To adequately sample the conformational space, we adapted the replica exchange with solute tempering (REST2) method to selectively heat the Aβ42 peptide in the absence and presence of HSA. Aβ42 binds to multiple sites on HSA with a preference to domain III and adopts various conformations that all differ from the free state. The β-sheet abundances of H14-E22 and A30-M33 regions are significantly reduced by HSA, so are the β-sheet lengths. HSA shifts the conformational ensemble towards more disordered states and alters the β-sheet association patterns. In particular, the frequent association of Q15-V24 and N27-V36 regions into β-hairpin which is critical for aggregation is impeded. HSA primarily interacts with the latter β-region and the N-terminal charged residues. They form promiscuous interactions characterized by salt bridges at the edge of the peptide-protein interface and hydrophobic cores at the center. Consequently, intrapeptide interactions crucial for β-sheet formation are disrupted. Our work builds the bridge between the modification of Aβ conformational ensemble and amyloid inhibition by HSA. It also illustrates the potential of the REST2 method in studying interactions between intrinsically disordered peptides and globular proteins.

Assessing SIRAH's Capability to Simulate Intrinsically Disordered Proteins and Peptides

The challenges posed by intrinsically disordered proteins (IDPs) to atomistic and coarse-grained (CG) simulations are boosting efforts to develop and reparametrize current force fields. An assessment of the dynamical behavior of IDPs' and unstructured peptides with the CG SIRAH force field suggests that the current version achieves a fair description of IDPs' conformational flexibility. Moreover, we found a remarkable capability to capture the effect of point mutations in loosely structured peptides.

Liquid-Liquid Phase Separation by Intrinsically Disordered Protein Regions of Viruses: Roles in Viral Life Cycle and Control of Virus-Host Interactions

Intrinsically disordered proteins (IDPs) are unable to adopt a unique 3D structure under physiological conditions and thus exist as highly dynamic conformational ensembles. IDPs are ubiquitous and widely spread in the protein realm. In the last decade, compelling experimental evidence has been gathered, pointing to the ability of IDPs and intrinsically disordered regions (IDRs) to undergo liquid-liquid phase separation (LLPS), a phenomenon driving the formation of membrane-less organelles (MLOs). These biological condensates play a critical role in the spatio-temporal organization of the cell, where they exert a multitude of key biological functions, ranging from transcriptional regulation and silencing to control of signal transduction networks. After introducing IDPs and LLPS, we herein survey available data on LLPS by IDPs/IDRs of viral origin and discuss their functional implications. We distinguish LLPS associated with viral replication and trafficking of viral components, from the LLPS-mediated interference of viruses with host cell functions. We discuss emerging evidence on the ability of plant virus proteins to interfere with the regulation of MLOs of the host and propose that bacteriophages can interfere with bacterial LLPS, as well. We conclude by discussing how LLPS could be targeted to treat phase separation-associated diseases, including viral infections.

Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly

Recent studies suggest that Tyr-39 might play a critical role for both the normal function and the pathological dysfunction of α-synuclein (αS), an intrinsically disordered protein involved in Parkinson’s disease. We perform here a comparative analysis between the structural features of human αS and its Y39A, Y39F, and Y39L variants. By the combined application of site-directed mutagenesis, biophysical techniques, and enhanced sampling molecular simulations, we show that removing aromatic functionality at position 39 of monomeric αS leads to protein variants populating more compact conformations, conserving its disordered nature and secondary structure propensities. Contrasting with the subtle changes induced by mutations on the protein structure, removing aromaticity at position 39 impacts strongly on the interaction of αS with the potent amyloid inhibitor phthalocyanine tetrasulfonate (PcTS). Our findings further support the role of Tyr-39 in forming essential inter and intramolecular contacts that might have important repercussions for the function and the dysfunction of αS.

Single-vesicle imaging quantifies calcium's regulation of nanoscale vesicle clustering mediated by α-synuclein

Although numerous studies have shown that the protein α-synuclein (α-Syn) plays a central role in Parkinson's disease, dementia with Lewy bodies, and other neurodegenerative diseases, the protein's physiological function remains poorly understood. Furthermore, despite recent reports suggesting that, under the influence of Ca²⁺, α-Syn can interact with synaptic vesicles, the mechanisms underlying that interaction are far from clear. Thus, we used single-vesicle imaging to quantify the extent to which Ca²⁺ regulates nanoscale vesicle clustering mediated by α-Syn. Our results revealed not only that vesicle clustering required α-Syn to bind to anionic lipid vesicles, but also that different concentrations of Ca²⁺ exerted different effects on how α-Syn induced vesicle clustering. In particular, low concentrations of Ca²⁺ inhibited vesicle clustering by blocking the electrostatic interaction between the lipid membrane and the N terminus of α-Syn, whereas high concentrations promoted vesicle clustering, possibly due to the electrostatic interaction between Ca²⁺ and the negatively charged lipids that is independent of α-Syn. Taken together, our results provide critical insights into α-Syn's physiological function, and how Ca²⁺ regulates vesicle clustering mediated by α-Syn.

N-Terminal Proteoforms in Human Disease

The collection of chemically different protein variants, or proteoforms, by far exceeds the number of protein-coding genes in the human genome. Major contributors are alternative splicing and protein modifications. In this review, we focus on those proteoforms that differ at their N termini with a molecular link to disease. We describe the main underlying mechanisms that give rise to such N-terminal proteoforms, these being splicing, initiation of protein translation, and protein modifications. Given their role in several human diseases, it is becoming increasingly clear that several of these N-terminal proteoforms may have potential as therapeutic interventions and/or for diagnosing and prognosing their associated disease.

Enhanced Molecular Dynamics Simulations of Intrinsically Disordered Proteins

Molecular dynamics simulations represent a powerful tool to gain insights into structural and dynamical features of biomolecular systems. Nevertheless, their recognized limitation in terms of achievable timescales becomes particularly severe when dealing with slow processes. In such cases, the employment of enhanced sampling methods, which allow accelerating the characterization of rare events in a timeframe consistent with conventional computational resources, results as crucial. In particular, such advanced techniques have proven highly valuable in the context of protein folding and, specifically, to explore the conformational ensemble spanned by intrinsically disordered proteins (IDPs). Here, we describe how to set up molecular dynamics simulations with one of these enhanced sampling approaches (namely, Parallel Tempering Metadynamics in the Well-Tempered Ensemble) using the N_TAIL peptide as a test case.

Depicting Conformational Ensembles of α-Synuclein by Single Molecule Force Spectroscopy and Native Mass Spectroscopy

Description of heterogeneous molecular ensembles, such as intrinsically disordered proteins, represents a challenge in structural biology and an urgent question posed by biochemistry to interpret many physiologically important, regulatory mechanisms. Single-molecule techniques can provide a unique contribution to this field. This work applies single molecule force spectroscopy to probe conformational properties of α-synuclein in solution and its conformational changes induced by ligand binding. The goal is to compare data from such an approach with those obtained by native mass spectrometry. These two orthogonal, biophysical methods are found to deliver a complex picture, in which monomeric α-synuclein in solution spontaneously populates compact and partially compacted states, which are differently stabilized by binding to aggregation inhibitors, such as dopamine and epigallocatechin-3-gallate. Analyses by circular dichroism and Fourier-transform infrared spectroscopy show that these transitions do not involve formation of secondary structure. This comparative analysis provides support to structural interpretation of charge-state distributions obtained by native mass spectrometry and helps, in turn, defining the conformational components detected by single molecule force spectroscopy.

An arsenal of methods for the experimental characterization of intrinsically disordered proteins - How to choose and combine them?

In this review, we detail the most common experimental approaches to assess and characterize protein intrinsic structural disorder, with the notable exception of NMR and EPR spectroscopy, two ideally suited approaches that will be described in depth in two other reviews within this special issue. We discuss the advantages, the limitations, as well as the caveats of the various methods. We also describe less common and more demanding approaches that enable achieving further insights into the conformational properties of IDPs. Finally, we present recent developments that have enabled assessment of structural disorder in living cells, and discuss the currently available methods to model IDPs as conformational ensembles.

Investigating targets for neuropharmacological intervention by molecular dynamics simulations

Medical research has identified over 500 brain disorders. Among these, there are still only very few neuropathologies whose causes are fully understood and, consequently, very few drugs whose mechanism of action is known. No FDA drug has been identified for major neurodegenerative diseases, such as Alzheimer's and Parkinson's. We still lack effective treatments and strategies for modulating progression or even early neurodegenerative disease onset diagnostic tools. A great support toward the highly needed identification of neuroactive drugs comes from computer simulation methods and, in particular, from molecular dynamics (MD). This provides insight into structure-function relationship of a target and predicts structure, dynamics and energetics of ligand/target complexes under biologically relevant conditions like temperature and physiological saline concentration. Here, we present examples of the predictive power of MD for neuroactive ligands/target complexes. This brief survey from our own research shows the usefulness of partnerships between academia and industry, and from joint efforts between experimental and theoretical groups.

Conformational ensemble of native α-synuclein in solution as determined by short-distance crosslinking constraint-guided discrete molecular dynamics simulations

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain β-structure elements and hairpins.

As the population ages, neurodegenerative diseases such as Parkinson’s disease will become an increasing problem in many countries. Aggregation of the protein α-synuclein is the primary cause of Parkinson’s disease, but there is still a dearth of structural information pertaining to the native, non-aggregating form of this protein. A better understanding the structural state of the native protein may prove useful for the design of new therapeutics to combat this disease. In order to obtain more structural information on this protein, we have recently modelled the native α-synuclein protein. These models were generated using a novel approach which combines protein crosslinking and discrete molecular dynamics simulations. We have found that the α-synuclein protein can adopt several shapes, all with a similar topology, resembling a three fingered closed claw. A region of the protein important for aggregation was found to be protected from the surrounding biological environment in these conformations, and the stabilization of these structures may be a fruitful avenue for future drug research into mitigating the cause and effect of Parkinson’s disease.

N-Terminal Acetylation Preserves α-Synuclein from Oligomerization by Blocking Intermolecular Hydrogen Bonds

The abnormal aggregation of α-synuclein (α-Syn) is closely associated with Parkinson's disease. Different post-translational modifications of α-Syn have been identified and contribute distinctly in α-Syn aggregation and cytotoxicity. Recently, α-Syn was reported to be N-terminally acetylated in cells, yet the functional implication of this modification, especially in α-Syn oligomerization, remains unclear. By using a solid-state nanopore system, we found that N-terminal acetylation can significantly decrease α-Syn oligomerization. Replica-exchange molecular dynamics simulations further revealed that addition of an acetyl group at the N-terminus disrupts intermolecular hydrogen bonds, which slows down the initial α-Syn oligomerization. Our finding highlights the essential role of N-terminal acetylation of α-Syn in preserving its native conformation against pathological aggregation.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Environmental and genetic factors support the dissociation between α-synuclein aggregation and toxicity

Synucleinopathies are a group of progressive disorders characterized by the abnormal aggregation and accumulation of α-synuclein (aSyn), an abundant neuronal protein that can adopt different conformations and biological properties. Recently, aSyn pathology was shown to spread between neurons in a prion-like manner. Proteins like aSyn that exhibit self-propagating capacity appear to be able to adopt different stable conformational states, known as protein strains, which can be modulated both by environmental and by protein-intrinsic factors. Here, we analyzed these factors and found that the unique combination of the neurodegeneration-related metal copper and the pathological H50Q aSyn mutation induces a significant alteration in the aggregation properties of aSyn. We compared the aggregation of WT and H50Q aSyn with and without copper, and assessed the effects of the resultant protein species when applied to primary neuronal cultures. The presence of copper induces the formation of structurally different and less-damaging aSyn aggregates. Interestingly, these aggregates exhibit a stronger capacity to induce aSyn inclusion formation in recipient cells, which demonstrates that the structural features of aSyn species determine their effect in neuronal cells and supports a lack of correlation between toxicity and inclusion formation. In total, our study provides strong support in favor of the hypothesis that protein aggregation is not a primary cause of cytotoxicity.

The Impact of N-terminal Acetylation of α-Synuclein on Phospholipid Membrane Binding and Fibril Structure

Human α-synuclein (αS) has been shown to be N terminally acetylated in its physiological state. This modification is proposed to modulate the function and aggregation of αS into amyloid fibrils. Using bacterially expressed acetylated-αS (NTAc-αS) and endogenous αS (Endo-αS) from human erythrocytes, we show that N-terminal acetylation has little impact on αS binding to anionic membranes and thus likely not relevant for regulating membrane affinity. N-terminal acetylation does have an effect on αS aggregation, resulting in a narrower distribution of the aggregation lag times and rates. 2D-IR spectra show that acetylation changes the secondary structure of αS in fibrils. This difference may arise from the slightly higher helical propensity of acetylated-αS in solution leading to a more homogenous fibril population with different fibril structure than non-acetylated αS. We speculate that N-terminal acetylation imposes conformational restraints on N-terminal residues in αS, thus predisposing αS toward specific interactions with other binding partners or alternatively decrease nonspecific interactions.

Computational scrutiny of the effect of N-terminal proline and residue stereochemistry in the nucleation of α-helix fold

N-Terminal l- to d-residue mutation nucleate helical fold in Ac–^DAla–^LAla₃–NHMe (Ib, m2), Ac–^DPro–^LAla₃–NHMe (IIb, m1), and Ac–^DPro–^LPro–^LAla₂–NHMe (IIIb, m2) peptides.