Molecular Dynamics-Assisted Interpretation of Experimentally Determined Intrinsically Disordered Protein Conformational Components: The Case of Human α-Synuclein

The role of solvation on the conformational landscape of α-synuclein

Native ion mobility mass spectrometry has been used extensively to characterize ensembles of intrinsically disordered protein (IDP) conformers, but the extent to which the gaseous measurements provide realistic pictures of the solution conformations for such flexible proteins remains unclear. Therefore, we systematically studied the relationship between the solution and gaseous structural ensembles by measuring electrospray charge state and collision cross section (CCS) distributions for cationic and anionic forms of α-synuclein (αSN), an anionic protein in solution, as well as directly probed gas phase residue to residue distances via ion/ion reactions between gaseous α-synuclein cations and disulfonic acid linkers that form strong electrostatic bonds. We also combined results from in-solution protein crosslinking identified from native tandem mass spectrometry (MS/MS) with an initial αSN ensemble generated computationally by IDPConformerGenerator to generate an experimentally restrained solution ensemble of αSN. CCS distributions were directly calculated for the solution ensembles determined by NMR and compared to predicted gaseous conformers. While charge state and collision cross section distributions are useful for qualitatively describing the relative structural dynamics of proteins and major conformational changes induced by changes to solution states, the predicted and measured gas phase conformers include subpopulations that are significantly different than those expected from completely "freezing" solution conformations and preserving them in the gas phase. However, insights were gained on the various roles of solvent in stabilizing various conformers for extremely dynamic proteins like α-synuclein.

Exploring Intra- and Inter-Regional Interactions in the IDP α-Synuclein Using smFRET and MD Simulations

Theoretical concepts from polymer physics are often used to describe intrinsically disordered proteins (IDPs). However, amino acid interactions within and between regions of the protein can lead to deviations from typical polymer scaling behavior and even to short-lived secondary structures. To investigate the key interactions in the dynamic IDP α-synuclein (αS) at the amino acid level, we conducted single-molecule fluorescence resonance energy transfer (smFRET) experiments and coarse-grained molecular dynamics (CG-MD) simulations. We find excellent agreement between experiments and simulations. Our results show that a physiological salt solution is a good solvent for αS and that the protein is highly dynamic throughout its entire chain, with local intra- and inter-regional interactions leading to deviations from global scaling. Specifically, we observe expansion in the C-terminal region, compaction in the NAC region, and a slightly smaller distance between the C- and N-termini than expected. Our simulations indicate that the compaction in the NAC region results from hydrophobic aliphatic contacts, mostly between valine and alanine residues, and cation-π interactions between lysine and tyrosine. In addition, hydrogen bonds also seem to contribute to the compaction of the NAC region. The expansion of the C-terminal region is due to intraregional electrostatic repulsion and increased chain stiffness from several prolines. Overall, our study demonstrates the effectiveness of combining smFRET experiments with CG-MD simulations to investigate the key interactions in highly dynamic IDPs at the amino acid level.

Clustering Heterogeneous Conformational Ensembles of Intrinsically Disordered Proteins with t-Distributed Stochastic Neighbor Embedding

Intrinsically disordered proteins (IDPs) populate a range of conformations that are best described by a heterogeneous ensemble. Grouping an IDP ensemble into "structurally similar" clusters for visualization, interpretation, and analysis purposes is a much-desired but formidable task, as the conformational space of IDPs is inherently high-dimensional and reduction techniques often result in ambiguous classifications. Here, we employ the t-distributed stochastic neighbor embedding (t-SNE) technique to generate homogeneous clusters of IDP conformations from the full heterogeneous ensemble. We illustrate the utility of t-SNE by clustering conformations of two disordered proteins, Aβ42, and α-synuclein, in their APO states and when bound to small molecule ligands. Our results shed light on ordered substates within disordered ensembles and provide structural and mechanistic insights into binding modes that confer specificity and affinity in IDP ligand binding. t-SNE projections preserve the local neighborhood information, provide interpretable visualizations of the conformational heterogeneity within each ensemble, and enable the quantification of cluster populations and their relative shifts upon ligand binding. Our approach provides a new framework for detailed investigations of the thermodynamics and kinetics of IDP ligand binding and will aid rational drug design for IDPs.

How Cu(II) binding affects structure and dynamics of α-synuclein revealed by molecular dynamics simulations

We report accelerated molecular dynamics simulations of α-Synuclein and its complex with two Cu(II) ions bound to experimentally determined binding sites. Adding two Cu(II) ions, one bound to the N-terminal region and one to the C-terminus, decreases size and flexibility of the peptide while introducing significant new contacts within and between N-terminus and non-Aβ component (NAC). Cu(II) ions also alter the pattern of secondary structure within the peptide, inducing more and longer-lasting elements of secondary structure such as β-strands and hairpins. Free energy surfaces, obtained from reweighting the accelerated molecular dynamics boost potential, further demonstrate the restriction on size and flexibility that results from binding of copper ions.

Native mass spectrometry for the investigation of protein structural (dis)order

A central challenge in structural biology is represented by dynamic and heterogeneous systems, as typically represented by proteins in solution, with the extreme case of intrinsically disordered proteins (IDPs) [1-3]. These proteins lack a specific three-dimensional structure and have poorly organized secondary structure. For these reasons, they escape structural characterization by conventional biophysical methods. The investigation of these systems requires description of conformational ensembles, rather than of unique, defined structures or bundles of largely superimposable structures. Mass spectrometry (MS) has become a central tool in this field, offering a variety of complementary approaches to generate structural information on either folded or disordered proteins [4-6]. Two main categories of methods can be recognized. On one side, conformation-dependent reactions (such as cross-linking, covalent labeling, H/D exchange) are exploited to label molecules in solution, followed by the characterization of the labeling products by denaturing MS [7-11]. On the other side, non-denaturing ("native") MS can be used to directly explore the different conformational components in terms of geometry and structural compactness [12-16]. All these approaches have in common the capability to conjugate protein structure investigation with the peculiar analytical power of MS measurements, offering the possibility of assessing species distributions for folding and binding equilibria and the combination of both. These methods can be combined with characterization of noncovalent complexes [17, 18] and post-translational modifications [19-23]. This review focuses on the application of native MS to protein structure and dynamics investigation, with a general methodological section, followed by examples on specific proteins from our laboratory.