Next Article in Journal
Distinct Cargos of Small Extracellular Vesicles Derived from Hypoxic Cells and Their Effect on Cancer Cells
Next Article in Special Issue
Spectroscopic Investigation of the Kinetic Mechanism Involved in the Association of Potyviral VPg with the Host Plant Translation Initiation Factor eIF4E
Previous Article in Journal
Adenosine A2A and A3 Receptors Are Able to Interact with Each Other. A Further Piece in the Puzzle of Adenosine Receptor-Mediated Signaling
Previous Article in Special Issue
Tau and Alpha Synuclein Synergistic Effect in Neurodegenerative Diseases: When the Periphery Is the Core
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 14
10.3390/ijms21145061
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Oscar Palomino-Hernandez
1,2,3,4,†,
Fiamma A. Buratti
5,†,
Pamela S. Sacco
5,
Giulia Rossetti
1,6,7,
Paolo Carloni
1,2,8,* and
Claudio O. Fernandez
5,9,*
1
Computational Biomedicine, Institute for Neuroscience and Medicine (INM-9) and Institute for Advanced Simulations (IAS-5), Forschungszentrum Jülich, 52425 Jülich, Germany
2
Faculty of Mathematics, Computer Science and Natural Sciences, RWTH Aachen, 52425 Aachen, Germany
3
Computation-Based Science and Technology Research Center, The Cyprus Institute, 2121 Nicosia, Cyprus
4
Institute of Life Science, The Hebrew University of Jerusalem, 91904 Jerusalem, Israel
5
Max Planck Laboratory for Structural Biology, Chemistry and Molecular Biophysics of Rosario (MPLbioR, UNR-MPIbpC) and Instituto de Investigaciones para el Descubrimiento de Fármacos de Rosario (IIDEFAR, UNR-CONICET), Universidad Nacional de Rosario, S2002LRK Rosario, Argentina
6
Department of Oncology, Hematology, Oncology, Hemostaseology, and Stem Cell Transplantation University Hospital Aachen, RWTH Aachen University, Pauwelsstraße 30, 52074 Aachen, Germany
7
Jülich Supercomputing Center (JSC), Forschungszentrum Jülich, 52425 Jülich, Germany
8
Institute for Neuroscience and Medicine (INM-11) Forschungszentrum Jülich, 52425 Jülich, Germany
9
Department of NMR-Based Structural Biology, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, D-37077 Göttingen, Germany
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to the manuscript.
Int. J. Mol. Sci. 2020, 21(14), 5061; https://doi.org/10.3390/ijms21145061
Submission received: 31 May 2020 / Revised: 2 July 2020 / Accepted: 13 July 2020 / Published: 17 July 2020
(This article belongs to the Special Issue Intrinsically Disordered Proteins (IDPs): From Physical Chemistry to Pathogenic Mechanisms)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Neurodegenerative disorders such as Parkinson’s and Alzheimer’s are among the human diseases associated with the self-assembly of polypeptides into amyloid structures [1,2]. Toward the design of effective therapeutics to combat these diseases, one of the major unanswered questions of protein aggregation is the propensity of particular primary sequences to aggregate [1,2]. These relevant regions are usually known as aggregation “hot spots”, and their study becomes critical to identify determinants governing protein aggregation and to reveal specific interactions that must be disrupted to prevent amyloid assembly. This knowledge can in turn be used for therapeutic intervention to target the aggregation pathway of these proteins and its associated toxicity [3].
The π-stacking of aromatic residues has been suggested to serve as structural and functional elements for the progression of the self-assembly process [4]. Furthermore, the role of aromatic–aromatic interactions in aggregation was also suggested by an examination of compounds that inhibit fibril formation [5]. Aromatic groups are a frequent feature in these inhibitors [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]; hence, the ability of these molecules to impair amyloid formation might be mediated by π−π interactions with aromatic residues in specific regions of the protein sequence.
In that context, the implication of Tyr residues on the aggregation of the amyloid protein alpha synuclein (αS), a pathological factor in Parkinson’s disease, has been explored by mutagenesis-based approaches [25,26,27,28,29]. In its monomeric, intrinsically disordered state, αS adopts an ensemble of conformations with no rigid secondary structure. However, long-range interactions have been shown to stabilize an aggregation-autoinhibited conformation. The protein contains three Tyr residues in the C-terminal region (Tyr-125, Tyr-133, Tyr-136) and one in the N-terminus (Tyr-39). Cysteine substitution in the latter but not in Tyr-125, Tyr-133, and Tyr-136 positions preferentially enhanced dimer and oligomer formation under oxidative conditions, suggesting a dimerization pathway mediated by this region [25]. On the other hand, substitution of Tyr-39 and Tyr-133 by Ala resulted in substantial inhibition of fibrillation, suggesting that Tyr-39 and the C-terminal Tyr residues might form an aromatic cluster that stabilizes the native state of αS [26]. More recently, the structural elucidation of the inhibitory interaction between αS and phthalocyanine tetrasulfonate (PcTS), a widely used anti-amyloid compound, provided the basis for the rational design of αS variants, indicating that the aromatic side chain of the Tyr residue at position 39 might have a critical role in the fibrillation pathway of αS [5,6]. Phthalocyanines such as PcTS are small molecules inhibiting the formation of fibrillar aggregates for targets such as αS, tau, and prion protein. The aromatic character of these molecules was proposed as a key feature to enable inhibitory interactions with these proteins [5].
Substituting Tyr by Ala [6,26] or Cys [25] not only replaces an aromatic side chain with an aliphatic one, but it also affects considerably the hydrophobicity, the size, and the β-sheet propensity of the side chain. To uncover the role of aromaticity at position 39 for aggregation, here we replaced Tyr-39 with the following residues: (i) Phe, which keeps aromaticity and removes a hydrophilic group such as OH; (ii) Leu, which is similar to Phe in size and hydrophobicity but it is not aromatic; (iii) Ala, which is hydrophobic as Phe and Leu but it is smaller in size. By using a combined approach based on experimental biophysics and replica exchange with solute tempering (REST2) enhanced sampling molecular simulations [30], we analyzed the impact of these mutations on monomeric αS’ structural properties and on the interactions of the protein with PcTS. The structural changes induced by the mutations turned out to be only local, with the monomeric state of αS remaining disordered upon replacement of Tyr-39 by Phe, Leu and Ala. Moreover, the four αS species showed no propensities for specific elements of secondary structure, adopting in all cases typical random-coil conformations. Interestingly, our simulations suggest that the conformations populated by the Y39A and the Y39L variants were more compact than Y39F and wild-type (wt) species. Aromaticity at position 39 resulted in a critical structural determinant for PcTS binding to the N-terminus of αS, an event that mediates the inhibition of this molecule on αS fibril assembly. Overall, these new findings give further support to a role for Tyr-39 in essential inter and intramolecular contacts that might have important repercussions for the function and the dysfunction of αS.

2. Results

2.1. Monomeric αS Remains Disordered after Y39 Mutation

The structural features of the monomeric wild-type (wt), Y39F, Y39L, and Y39A αS species were investigated by circular dichroism (CD), NMR spectroscopy, along with REST2 simulations [30]. The latter allows one to enhance the sampling of the conformational space of the protein relative to plain molecular dynamics (MD) [30]. As determined for the wild-type protein, the far-UV CD spectra of all single mutants at position 39 were indicative of largely disordered conformations [6], as manifested in prominent negative minima at ~199 nm (Figure 1 Left). The predicted CD spectra of the four simulated variants were consistent with the experimental ones (Figure 1 Right). In particular, it suggests that random coil is the most predominant motif on the wt and the three variants.
The NMR spectra of non-labeled wt αS in Buffer A showed well-resolved clusters of resonances in the 6.0–8.0 ppm (1H) range, comprising the side chains of different aromatic residues: His (aa50), Phe (aa4, aa94), and Tyr (aa39, aa125, aa133, aa136) (Figure 2) [6,31]. The replacement of tyrosine at position 39 by phenylalanine, leucine, and alanine residues in the 1D 1H NMR spectra of the mutant species was confirmed by the lack of signals corresponding to the Tyr-39 residue (Figure 2B–D) and the detection of new signals corresponding to the presence of an additional phenylalanine residue in the Y39F variant (Figure 2B).
The overlaid 1H-15N heteronuclear single quantum correlation (HSQC) spectra of a 50 µM sample of uniformly 15N-labeled wt αS and its Y39L variant recorded in Buffer A at 15 °C is shown in Figure 3A. The resonances were well resolved and sharp, with a poor dispersion of chemical shifts, reflecting the disordered nature of the backbone (Figure S1). As shown in Figure 3B, changes induced by the mutations were only local and were restricted to the immediate vicinity of the Tyr residue. The absence of noticeable changes in the structural properties indicated that the monomeric state of the protein variants remained disordered upon replacement of Tyr-39 by Phe, Leu, and Ala.

2.2. Structural Propensities of the Y39 αS Variants

We also evaluated the different protein variants in terms of their conformational properties. 3JHN-Hα couplings are reliable quantitative reporters of the time-averaged distribution of the backbone torsion angles φ and are frequently used to probe the propensity of intrinsically disordered proteins (IDP) to sample different regions of conformational space [32]. Therefore, we measured residue-specific 3JHN-Hα couplings in all of these proteins. As shown in Figure 3C, the values measured for the wt and Y39F, Y39L, and Y39A αS proteins were typical of random-coil conformations and essentially indistinguishable among the mutant species studied.
The determination of the hydrodynamic properties of a macromolecule was extensively applied in the study of conformational changes accompanying processes such as molecular association and folding [33,34]. We then employed pulse field gradient-NMR (PFG-NMR) to measure the hydrodynamic properties of wt αS and its mutant variants at position 39. The values of the hydrodynamic radius (RH) determined for the protein in its native state (RH = 31.0 ± 0.4 Å) were consistent with previous determinations [31,35], whereas no changes were detected in the Y39F variant (RH = 30.9 ± 0.4 Å). Compared with the parameters measured for the Y39A (RH = 28.0 ± 0.4 Å) and the Y39L (RH = 28.4 ± 0.4 Å) mutants, our results suggest that the latter mutations might affect the size of native αS or cause some collapse to a more compact species. The small percentage deviation across the RH values of the four αS species allowed us to detect differences in the compactness of the protein. This was shown by the following facts: (i) the measured RH values of αS in its native state (31.9 Å) and in the presence of 8 M urea (35.0 Å) reflected the decreased persistence of the residual long-range interactions in αS, as demonstrated conclusively by paramagnetic relaxation enhancement measurements [36]. Thus, the change of ~10% (3.1 Å) in RH was significant and consistent with the transition of αS from its native state towards an ensemble of more extended, unfolded conformations. (ii) The measured RH value of αS in dilute solution (26.6 ± 0.5 Å) changed to 22.5 ± 0.6 Å in 1 M glucose. The decrease of ~15% (4 Å) was consistent with the fact that the protein adopts more compact states under crowding conditions [37,38].
In addition to the previous results, we expressed the compactness of the protein in terms of “compaction factors”, which relate the measured RH with empirical estimations of minimum and maximum values expected for the length of the polypeptide chain for each conformational ensemble of a protein [36,39,40]. The values ranged from 1.0 for a natively folded protein to 0.0 for a pure random coil. The calculated compaction factors turned out to be 0.30 for wt and for Y39F αS, consistently with previously published results [36]. The compaction factors were instead 0.50 both for Y39L αS and Y39A αS. Such an increase of ~65% in compaction factor values was significantly indicative of a higher degree of compaction for Y39A and Y39L relative to wt αS. We concluded that the Y39A and the Y39L protein variants adopted an ensemble of more compact conformations in solution than the Y39F and the wt proteins.
The REST2 simulations were validated through comparison with experimental data. The predicted chemical shifts N, Cα, and Cβ atoms of wt correlated well with the corresponding experimental values, while the C atoms had more outliers, making the correlation only fair (Figure 4). The wt protein exhibited partial (10%) α-helical propensities within its first ten residues (Figure S2), consistently with the experiment [41]. The REST2 simulations suggested that the four αS species shared similar levels of secondary structure. Indeed, the content of secondary structure did not change significantly on passing from the wt to the three mutants investigated here (Table 1). However, the latter exhibited a different local organization of the secondary structure’s motifs (Figure S2). These results are in agreement with the above reported 3JHN-Hα profiles, indicating that changes induced by the mutations were only local and were restricted to the immediate vicinity of the position-39 residue.
The compactness of the proteins was further investigated through a projection of the radius of gyration (RG) and the end-to-end distance (EE) as reaction coordinates, as performed for other IDP analyses [42]. Figure 5 shows that wt and Y39F shared a similar compactness profile. They populated mostly regions of the conformational space with relatively large RG and EE, therefore having preference for more expanded structures. In contrast, Y39L and Y39A mostly populated regions with progressively lower values of RG and EE (additional information on the minima is reported in Table S1). At the qualitative level, this suggested the following trend in compactness: Y39A > Y39L > Y39F ~ wt, which is consistent with the experimental results (vide supra).
We next analyzed the proteins intramolecular hydrogen bonds (HB) and hydrophobic contacts across the simulations of the four species.
Overall, Y39L and Y39F αS exhibited more HB contacts than the wt, with significant changes in the non-amyloid-β component (NAC) region (Figure 6 and Table S2). Instead, the number of salt-bridges is similar across the four proteins (data not shown). As far as the position 39 was concerned, we found that Y39 side chain in wt was mostly solvent exposed, with an average of ~1.7 HBs with the solvent during the simulation. It formed intramolecular HB contacts (with K21/G25/T44 backbone units) for only 2% of the simulated time. The backbone unit of the residue across the four species was also mostly solvated. It formed HB contacts with the backbone units of neighboring residues for only less than 10% of the simulated time. We concluded that the content of intramolecular HB contacts was rather small, albeit it varied across the four species (Figure 6).
The overall number of hydrophobic contacts increased on passing from wt to Y39F and Y39L αS, while it decreased on passing from wt to Y39A αS (Figure 6 and Table S3). However, while the first three species formed only intra-region contacts, Y39A αS also displayed inter-region contacts. The residue in position 39 formed significant interactions with its neighboring residues; it interacted with (i) L38 for almost the entire simulation time across the four species; (ii) V40 for ~50 to ~70% of simulated time on passing from wt/Y39A to Y39F and Y39L αS; (iii) both G47 and V48 for ~55% of the simulated time in Y39A αS, ~70% in Y39L αS/wt; ~90% in Y39F αS.
We concluded that the intramolecular contacts of residue in position 39 across the four species involved the N-term residues, and they were basically hydrophobic in nature. Indeed, as pointed out above, the HB contacts of the residue in position 39 are formed mostly with the solvent across the four species. These contacts may have played a role for the observed increase of compactness on passing from wt and Y39F αS to the other two mutants (Figure 5). On one hand, Y39L αS showed a higher number of hydrophobic contacts than Y39F αS and wt (Figure 6). This might have contributed to the formation of a more compact structure. On the other hand, as already mentioned above, Y39A αS was the only species out of the four investigated here forming significant inter-region residue hydrophobic contacts. The latter may have been important for increasing the compactness of the protein. The difference between Y39L and Y39A variants may have been caused, at least in part, by the different conformational properties of the two residues [43].

2.3. Aromaticity at Position 39 and PcTS Binding

We then studied the mutants Y39F, Y39L, and Y39A of αS, aiming to determine the role of these residues in directing the binding of the anti-amyloid compound PcTS to the N-terminus (Figure 7). As previously shown, complete loss of PcTS binding to the 35–41 region was observed when the Tyr residue in position 39 was replaced by Ala [6]. Replacing the aromatic side-chain by leucine in the Y39L mutant was also sufficient to impair PcTS binding. Conversely, removing the aromatic side chain in position 39 and replacing it with Phe did not alter PcTS interaction, as revealed by the binding features of PcTS at this site. Altogether, these results demonstrate unequivocally that aromaticity at position 39 also had a critical role for PcTS binding to the N-terminus of αS and that specific aromatic interactions with the Y39 residue provided a central mechanistic basis for the inhibitory process of PcTS on αS fibrillation.

3. Discussion

In this work, we investigated the importance of aromaticity at position 39 in αS through the analysis of the impact of Y39F, Y39L, and Y39A mutations on the structural properties of its monomeric state and its interaction with the anti-amyloid agent PcTS.
By using biophysical techniques, we showed that removing aromatic functionality at position 39 of monomeric αS does not affect the disordered nature or the secondary structure propensities of the protein. Added to that, our experimental and simulated results indicate that replacement of Tyr-39 by Ala and Leu leads to protein variants that populate more compact conformations. These αS mutants show a modified number and type of contacts compared to wt and Y39F species. Assuming more expanded conformations for the αS species containing aromatic residues at position 39 and based on the experimental evidences from previous studies performed on the Y39A protein variant [6,26], one could then predict the following trend of protein fibril assembly propensities: wt ~ Y39F > Y39A, Y39L.
On the other hand, contrasting with the subtle changes induced by mutations on the protein structure, removing aromaticity at position 39 not only affects the compactness of the protein but also impacts strongly on the interaction of αS with PcTS. The inhibitory mechanism exerted by phthalocyanines on amyloid assembly is proposed to be a direct consequence of their interaction with target proteins [5]. These molecular interactions are influenced strongly by the highly aromatic character of the cyclic tetrapyrrole ring system, which contributes importantly to the ability of these molecules to bind strongly and selectively to a protein via π−π interactions with aromatic residues. In our previous studies, we showed that the basis for the inhibitory effect of this compound on the amyloid assembly of αS relies on its binding to the Y39 site at the N-terminus sequence of αS [6]. Added to that, in the current study, we performed the structural characterization of the PcTS interaction with both the Y39A and the Y39L variants of the protein, concluding that this interaction is primarily driven by π−π stacking interactions between the aromatic side-chain of Tyr-39 and the aromatic ring system of phthalocyanines. These findings provide clear evidence that aromatic interactions are required for the formation of the αS–PCTS complex, a key event in the molecular pathway that leads to the inhibition of αS fibril formation.
Indeed, the fact that residue 39 was found to be solvent exposed not only in the monomer but also in transient species formed during the early stages of αS fibrillation [44] implies that Y39 could also be targeted directly by PcTS in early aggregate αS species, a fact that would be consistent with the evidences reporting that PcTS molecules appear mostly incorporated into the αS aggregates. In that direction, the crystal structures of fibrillar αS reveal that Tyr-39 is involved in π−π self-stacking interactions (Figure 8), with ring distances around 4.86 ± 0.12 Å, a very small inter-ring angle, and a χ2 dihedral angle of −80 degrees, indicating that Tyr-39 π-stacking might play a role for the aggregation of monomers or oligomers in the earliest stages of αS fibril formation. However, this is not relevant for later stages of the process, such as the dimerization of fibrils, for which the NAC region is known to play a major role [45,46,47]. Indeed, the core region of amyloid fibrils of αS was shown to begin somewhere in the range of residues 31–39 [48,49], precisely in the close vicinity of the primary site targeted by phthalocyanines.
It is noteworthy that interactions governing the protein–compound reactions seem to be of the same nature of those responsible for protein–protein association reactions. Hence, perturbation of intra- and intermolecular interactions on passing from wt αS to the Y39A and the Y39L variants might also affect both the π–π stacking formation and the aggregation process. Our results reported here for the first time open the door toward experimental studies in the interface between molecular biophysics and cellular biology aimed to investigate the relative contributions of π-stacking versus hydrophobic packing in the aggregation of αS.

4. Materials and Methods

4.1. Proteins and Reagents

Unlabeled and 15N isotopically enriched N-terminally acetylated αS and its Y39F, Y39L, and Y39A variants were obtained by co-transforming Escherichia coli BL21 with the plasmid harboring the corresponding protein gene and a second one that encodes for the components of yeast NatB acetylase complex [50]. Both plasmids carried different antibiotic resistance, namely ampicillin and chloramphenicol, to select the doubly transformed E. coli colonies. Purification was carried out as previously reported [51] with the exception that both antibiotics were included in the growth flasks to avoid plasmid purge during growth and expression. The final purity of the protein samples was determined by SDS-PAGE. Purified protein samples were dissolved in 20 mM 2-(N-morpholino)ethanesulfonic acid (MES) buffer supplemented with 100 mM NaCl at pH 6.5 (Buffer A). Protein concentrations were determined spectrophotometrically by measuring absorption at 274 nm and using an epsilon value of 5600 M−1·cm−1.

4.2. NMR Experiments

NMR spectra were recorded on a Bruker 600 MHz HD Avance III spectrometer equipped with a cryogenically cooled triple resonance inverse (TCI) 1H (13C/15N) probe. One-dimensional 1D 1H-NMR experiments were acquired at 15 °C on 100 µM unlabeled αS samples dissolved in Buffer A. Two-dimensional 2D 1H-15N heteronuclear single quantum correlation (HSQC) experiments were performed with pulsed-field gradient enhanced pulse sequences on 50 µM 15N-labeled protein samples dissolved in Buffer A at 15 °C [52]. Aggregation did not occur under these low temperature conditions and absence of stirring. Mean weighted chemical-shift displacements (MW 1H-15N ΔCS) were calculated as [(Δδ1H)2 + (Δδ15N)2/25]1/2 [52]. Three-bond HN-Hα coupling constants (3JHN-Hα) were obtained from the ratio between the intensities of the diagonal peaks and the cross-peaks in the HNHA experiment, which are three-dimensional experiments designed to accurately determine three-bond HN-Hα J-coupling constants [53]. Three-bond HN-Hα coupling constants (3JHN-Hα) are sensitive to the torsion angle φ populated by each residue in the protein sequence and thus report on secondary structure content. This coupling fell in the range 3.0–6.0 Hz for an α-helix and 8.0–11.0 Hz for a β-sheet structure. For a random-coil, a weighted average of these values was observed, which typically ranged between 6.0 and 8.0 Hz for most residues [54,55]. Pulse field gradient-NMR experiments were acquired at 15 °C on 100 µM unlabeled αS samples dissolved in D2O and containing dioxane as an internal radius standard (2.12 Å) and viscosity probe. A series of 20 one-dimensional spectra were collected as a function of gradient amplitude. The gradient strength was shifted from 1.69 to 33.72 G cm−1 in a linear manner. For the mapping experiments with the PcTS compound, 1H-15N HSQC amide cross-peaks affected during titrations with the ligand molecule were identified by comparing their intensities (I) with those of the same cross-peaks in the data set of free protein (I0) [6]. The I/I0 ratios of non-overlapping cross-peaks belonging to residues in the N-terminal region were plotted as a function of the protein sequence to obtain the intensity profiles. Acquisition and processing of NMR spectra were performed using TOPSPIN 7.0 (Bruker Biospin, https://www.bruker.com/nc.html Bruker Scientific LLC. Billerica, MA 01821. USA). The 2D spectra analysis and visualization were performed with CCPN (Collaborative Computational Project for NMR).

4.3. CD Spectroscopy

αS samples were diluted 10-fold in Buffer A, and CD spectra were recorded on a JASCO J-530 spectropolarimeter (JASCO, https://jascoinc.com, Mary’s Court Easton, MD 21601. USA).

4.4. Molecular Simulations

The calculations were based on our predicted conformational ensemble of wt monomeric, human N-terminus acetylated alpha synuclein [56]. Our prediction was shown to be fully consistent with available experimental data [56]. We first clustered the conformational ensemble according to the gromos algorithm proposed by Daura et al. [57] and implemented in the software GROMACS 2016.4 (http://www.gromacs.org) [58]. Then, we selected the most extended structure among the cluster representatives. The Y39F, the Y39L, and the Y39A mutations were created using the Schrödinger software [59]. The wt and the three variants were solvated in a truncated dodecahedral box insuring at least 1.0 nm of solvation shell and neutralized with Na+ and Cl− and ions until achieving an excess salt concentration of 15 mM. The topology and the coordinate files were prepared using GROMACS 2016.4 patched with the software PLUMED 2.4.0 (https://www.plumed.org) [60].
The force field for the water, the counterions, and the protein was a99SB-disp [61]. Simulations based on the latter have been shown to reproduce a variety of structural properties of non-folded proteins [61]. The Particle Mesh Ewald method [62] was used to treat the long-range interactions with a real space cutoff of 12 Å. The same cutoff was used for the van der Waals interactions. The LINear Constraint Solver (LINCS) algorithm [63] was used to constrain all bond lengths. Constant pressure and temperature conditions were considered coupling the systems to a Nosé–Hoover thermostat [64,65] at 288 K and a Parrinello–Rahman barostat [66] at 1 atm. A timestep of 2 fs was used for the simulations.
After 50 ns of plain MD simulations, the wt and the three variants underwent 25 ns of Hamiltonian replica exchange with solute tempering (REST2) simulations [30]. The same computational setup as above was applied. For an adequate sampling, 26 replicas were distributed in a range from 288 K to 533 K, achieving an exchange probability around 0.25 during the simulation. The range of temperatures was chosen for promoting the conformational sampling without increasing significantly the number of residues outside the allowed Ramachandran areas, similar to previous works [56].
A variety of properties were calculated for the last 20 ns of the replica with the lowest temperature, either using the full ensemble of structures (items (i)–(iv)) or five representative conformations spanning ca. 65% of the conformational space, identified with a cluster analysis using the gmx cluster and the gromos algorithm [57] in GROMACS 2016.4 (items (v)–(vii)): (i) the radius of gyration (RG) and end-to end distance (EE) with the gmx polystat code from GROMACS 2016.4 [58]; (ii) the inter-residue contacts with the gmx select code from GROMACS 2016.4; (iii) the H-bonds (HB) and the salt-bridges (SB) with the gmx hbond code from GROMACS 2016.4 or VMD [67]. A hydrogen bond was considered to be formed when acceptor and donor heavy atoms were at distance of 3.5 Å or lower, and the angle formed by acceptor, hydrogen atoms, and donor was of 30 degrees or less, while, for a salt bridge, the distance between the nitrogen atoms of basic residues and the oxygen atoms of acid residues was set to no more than 3.2 Å; (iv) the hydrophobic contacts analysis, the distance between heavy atom side-chains, was computed with gmx distance and gmx mdmat codes present in GROMACS 2016.4. The following residues were considered for the analysis: Ala, Leu, Met, Val, Phe, Pro, Tyr, and Ile. A cut-off of 6.0 Å was set. We considered contacts occurring for 15% or more of the simulation time. Analysis with different choices of these conditions turned out not to differ from our (qualitative) conclusions regarding the hydrophobic interactions in Section 2.2 (data not shown); (v) the circular dichroism (CD) spectrum with PDB2CD (https://pdb2cd.cryst.bbk.ac.uk) [68]; (vi) the chemical shifts of the C, Cα, Cβ, and N atoms, with SHIFTX2 1.10 (http://www.shiftx2.ca) [69]; (vii) the secondary structure elements with DSSP (https://swift.cmbi.umcn.nl/gv/dssp/) [70].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/14/5061/s1. Figure S1. Overlaid contour plots of 1H-15N HSQC spectra of 50 μM wt (black) and Y39L αS (red), showing all amide resonances detected. Most-affected residues located in the vicinity of the mutated site are labeled. Figure S2. Secondary structure content for individual residues across the four species investigated here during the simulation. Table S1. Characteristics of the RG and EE distributions across the four variants. The centers depicted in Figure 6 are described as the respective joint distributions. Additionally, the values of the mode and average of the marginal distributions (the univariate, independent ones) are depicted below. The conclusions of the study were based on the joint distribution of EE and RG. Table S2. Intramolecular HB content for the four variants within the structural ensemble. Table S3. Intramolecular hydrophobic contacts for the four variants within the structural ensemble.

Author Contributions

Conceptualization, C.O.F., G.R., and P.C.; methodology, all.; software, O.P.-H., F.A.B., and P.S.S.; validation, all.; formal analysis, O.P.-H., F.A.B., and P.S.S.; investigation, all.; resources, C.O.F. and P.C.; data curation, O.P.-H., F.A.B., and P.S.S.; writing—original draft preparation, all.; writing—review and editing, all; visualization, O.P.-H., F.A.B., and P.S.S.; supervision, C.O.F., G.R., and P.C.; project administration, C.O.F., G.R., and P.C.; funding acquisition, C.O.F., G.R., and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

O.P. was funded by the Marie Skłodowska-Curie Innovative Training Network grant 765048, “STIMULATE”.

Acknowledgments

C.O.F. thanks Universidad Nacional de Rosario (UNR) and ANPCyT- FONCyT (PICT 2014-3704 and PICT 2017-4665) for financial support. C.O.F. thanks MPIbpC (P10390) for support. P.S.S. thanks UNR for fellowships. F.A.B. thanks CONICET for fellowship. O. P-H. has been supported by the European Union’s Horizon 2020 research and innovation programme under Grant Agreement No 765048. The authors gratefully acknowledge the computing time granted by the JARA Vergabegremium and provided on the JARA Partition part of the supercomputer CLAIX at RWTH Aachen University, with project ID: jara0209.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Soto, C. Unfolding the role of protein misfolding in neurodegenerative diseases. Nat. Rev. Neurosci. 2003, 4, 49–60. [Google Scholar] [CrossRef] [PubMed]
  2. Chiti, F.; Dobson, C.M. Protein misfolding, functional amyloid, and human disease. Annu. Rev. Biochem. 2006, 75, 333–366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Liu, T.; Bitan, G. Modulating self-assembly of amyloidogenic proteins as a therapeutic approach for neurodegenerative diseases: Strategies and mechanisms. ChemMedChem 2012, 7, 359–374. [Google Scholar] [CrossRef] [PubMed]
  4. Gazit, E. A possible role for pi-stacking in the self-assembly of amyloid fibrils. FASEB J. 2002, 16, 77–83. [Google Scholar] [CrossRef]
  5. Valiente-Gabioud, A.A.; Miotto, M.C.; Chesta, M.E.; Lombardo, V.; Binolfi, A.; Fernandez, C.O. Phthalocyanines as molecular scaffolds to block disease-associated protein aggregation. Acc. Chem. Res. 2016, 49, 801–808. [Google Scholar] [CrossRef] [Green Version]
  6. Lamberto, G.R.; Binolfi, A.; Orcellet, M.L.; Bertoncini, C.W.; Zweckstetter, M.; Griesinger, C.; Fernandez, C.O. Structural and mechanistic basis behind the inhibitory interaction of PcTS on alpha-synuclein amyloid fibril formation. Proc. Natl. Acad. Sci. USA 2009, 106, 21057–21062. [Google Scholar] [CrossRef] [Green Version]
  7. Bulic, B.; Pickhardt, M.; Khlistunova, I.; Biernat, J.; Mandelkow, E.M.; Mandelkow, E.; Waldmann, H. Rhodanine-based tau aggregation inhibitors in cell models of tauopathy. Angew. Chem. Int. Ed. 2007, 46, 9215–9219. [Google Scholar] [CrossRef]
  8. Schenk, D.; Basi, G.S.; Pangalos, M.N. Treatment strategies targeting amyloid β-protein. Cold Spring Harb. Perspect. Med. 2012, 2, a006387. [Google Scholar] [CrossRef] [Green Version]
  9. Masuda, M.; Suzuki, N.; Taniguchi, S.; Oikawa, T.; Nonaka, T.; Iwatsubo, T.; Hisanaga, S.-I.; Goedert, M.; Hasegawa, M. Small molecule inhibitors of alpha-synuclein filament assembly. Biochemistry 2006, 45, 6085–6094. [Google Scholar] [CrossRef]
  10. Caughey, B.; Caughey, W.S.; Kocisko, D.A.; Lee, K.S.; Silveira, J.R.; Morrey, J.D. Prions and transmissible spongiform encephalopathy (TSE) chemotherapeutics: A common mechanism for anti-TSE compounds? Acc. Chem. Res. 2006, 39, 646–653. [Google Scholar] [CrossRef]
  11. Ehrnhoefer, D.E.; Bieschke, J.; Boeddrich, A.; Herbst, M.; Masino, L.; Lurz, R.; Engemann, S.; Pastore, A.; Wanker, E.E.; Ehrnhoefer, D.E. EGCG redirects amyloidogenic polypeptides into unstructured, off-pathway oligomers. Nat. Struct. Mol. Biol. 2008, 15, 558–566. [Google Scholar] [CrossRef] [PubMed]
  12. Caughey, W.S.; Priola, S.A.; Kocisko, D.A.; Raymond, L.D.; Ward, A.; Caughey, B. Cyclic tetrapyrrrole sulfonation, metals, and oligomerization in antiprion activity. Antimicrob. Agents Chemother. 2007, 51, 3887–3894. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wagner, J.; Ryazanov, S.; Leonov, A.; Levin, J.; Shi, S.; Schmidt, F.; Prix, C.; Pan-Montojo, F.; Bertsch, U.; Mitteregger-Kretzschmar, G. Anle138b: A novel oligomer modulator for disease-modifying therapy of neurodegenerative diseases such as prion and Parkinson’s disease. Acta Neuropathol. 2013, 125, 795–813. [Google Scholar] [CrossRef] [Green Version]
  14. Levin, J.; Schmidt, F.; Boehm, C.; Prix, C.; Bötzel, K.; Ryazanov, S.; Leonov, A.; Griesinger, C.; Giese, A. The oligomer modulator anle138b inhibits disease progression in a Parkinson mouse model even with treatment started after disease onset. Acta Neuropathol. 2014, 127, 779–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Scherzer-Attali, R.; Shaltiel-Karyo, R.; Adalist, Y.H.; Segal, D.; Gazit, E. Generic inhibition of amyloidogenic proteins by two naphthoqui-none-tryptophan hybrid molecules. Proteins 2012, 80, 1962–1973. [Google Scholar] [CrossRef] [PubMed]
  16. Frydman-Marom, A.; Shaltiel-Karyo, R.; Moshe, S.; Gazit, E. The generic amyloid formation inhibition effect of a designed small aromatic β-breaking peptide. Amyloid 2011, 18, 119–127. [Google Scholar] [CrossRef]
  17. Chakraborty, R.; Sahoo, S.; Halder, N.; Rath, H.; Chattopadhyay, K. Conformational-Switch Based Strategy Triggered by [18] π Heteroannulenes toward Reduction of Alpha Synuclein Oligomer Toxicity. ACS Chem. Neurosci. 2019, 10, 573–587. [Google Scholar] [CrossRef]
  18. Lizarraga, F.G.; Socías, S.B.; Avila, C.L.; Torres-Bugeau, C.M.; Barbosa, L.R.; Binolfi, A.; Sepúlveda-Díaz, J.E.; Del Bel, E.; Fernández, C.O.; Papy-Garcia, D.; et al. Repurposing doxycycline for synucleinopathies: Remodelling of α-synuclein oligomers towards non-toxic parallel beta-sheet structured species. Sci. Rep. 2017, 7, 41755. [Google Scholar] [CrossRef] [Green Version]
  19. Pujols, J.; Díaz, S.P.; Lázaro, D.F.; Peccati, F.; Pinheiro, F.; González, D.; Carija, A.; Navarro, S.; Giménez, M.C.; García, J.; et al. Small molecule inhibits α-synuclein aggregation, disrupts amyloid fibrils, and prevents degeneration of dopaminergic neurons. Proc. Natl. Acad. Sci. USA 2018, 115, 10481–10486. [Google Scholar] [CrossRef] [Green Version]
  20. Valdinocci, D.; Grant, G.D.; Dickson, T.; Pountney, D. Epothilone D inhibits microglia-mediated spread of alpha-synuclein aggregates. Mol. Cell Neurosci. 2018, 89, 80–94. [Google Scholar] [CrossRef] [Green Version]
  21. Schwab, K.; Frahm, S.; Horsley, D.; Rickard, J.E.; Melis, V.; Goatman, E.A.; Magbagbeolu, M.; Douglas, M.; Leith, M.G.; Baddeley, T.C.; et al. A protein aggregation inhibitor, leuco-methylthioninium bis(hydromethanesulfonate), decreases α-synuclein inclusions in a transgenic mouse model of synucleinopathy. Front. Mol. Neurosci. 2017, 10, 447. [Google Scholar] [CrossRef] [PubMed]
  22. Palazzi, L.; Bruzzone, E.; Bisello, G.; Leri, M.; Stefani, M.; Bucciantini, M.; De Laureto, P.P. Oleuropein aglycone stabilizes the monomeric α-synuclein and favours the growth of non-toxic aggregates. Sci. Rep. 2018, 8, 8337. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Jha, N.N.; Kumar, R.; Panigrahi, R.; Navalkar, A.; Ghosh, D.; Sahay, S.; Mondal, M.; Kumar, A.; Maji, S.K. Comparison of α-synuclein fibril inhibition by four different amyloid inhibitors. ACS Chem. Neurosci. 2017, 8, 2722–2733. [Google Scholar] [CrossRef] [PubMed]
  24. Reiner, A.M.; Schmidt, F.; Ryazanov, S.; Leonov, A.; Weckbecker, D.; Deeg, A.A.; Griesinger, C.; Giese, A.; Zinth, W. Photophysics of diphenyl-pyrazole compounds in solutions and α-synuclein aggregates. Biochim. Biophys. Acta Gen. Subj. 2018, 1862, 800–807. [Google Scholar] [CrossRef] [Green Version]
  25. Zhou, W.; Freed, C.R. Tyrosine-to-cysteine modification of human α-synuclein enhances protein aggregation and cellular toxicity. J. Biol. Chem. 2004, 279, 10128–10135. [Google Scholar] [CrossRef] [Green Version]
  26. Ulrih, N.P.; Barry, C.H.; Fink, A.L. Impact of Tyr to Ala mutations on α-synuclein fibrillation and structural properties. Biochim. Biophys. Acta 2008, 1782, 581–585. [Google Scholar] [CrossRef] [Green Version]
  27. Kaylor, J.; Bodner, N.; Edridge, S.; Yamin, G.; Hong, D.P.; Fink, A.L. Characterization of oligomeric intermediates in alpha-synuclein fibrillation: FRET studies of Y125W/Y133F/Y136F alpha-synuclein. J. Mol. Biol. 2005, 353, 357–372. [Google Scholar] [CrossRef]
  28. Zhu, M.; Rajamani, S.; Kaylor, J.; Han, S.; Zhou, F.; Fink, A.L. The flavonoid baicalein inhibits fibrillation of alpha-synuclein and disaggregates existing fibrils. J. Biol. Chem. 2004, 279, 26846–26857. [Google Scholar] [CrossRef] [Green Version]
  29. Burai, R.; Ait-Bouziad, N.; Chiki, A.; Lashuel, H.A. Elucidating the Role of Site-Specific Nitration of α-Synuclein in the Pathogenesis of Parkinson’s Disease via Protein Semisynthesis and Mutagenesis. J. Am. Chem. Soc. 2015, 137, 5041–5052. [Google Scholar] [CrossRef]
  30. Wang, L.; Friesner, R.A.; Berne, B.J. Replica exchange with solute scaling: A more efficient version of replica exchange with solute tempering (REST2). J. Phys. Chem. B 2011, 115, 9431–9438. [Google Scholar] [CrossRef] [Green Version]
  31. Lamberto, G.R.; Torres-Monserrat, V.; Bertoncini, C.W.; Salvatella, X.; Zweckstetter, M.; Griesinger, C.; Fernandez, C.O. Towards the discovery of effective polycyclic inhibitors of alpha-Synuclein amyloid assembly. J. Biol. Chem. 2011, 286, 32036–32044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Gentile, I.; Garro, H.A.; Delgado Ocaña, S.; Gonzalez, N.; Strohäker, T.; Schibich, D.; Quintanar, L.; Sambrotta, L.; Zweckstetter, M.; Griesinger, C.; et al. Interaction of Cu(I) with the Met-X3-Met motif of alpha-synuclein: Binding ligands, affinity and structural features. Metallomics 2018, 10, 1383–1389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Jones, J.A.; Wilkins, D.K.; Smith, L.J.; Dobson, C.M. Characterization of protein unfolding by NMR diffusion measurements. J. Biomol. NMR 1997, 10, 199–203. [Google Scholar] [CrossRef]
  34. Altieri, A.S.; Hinton, D.P.; Byrd, R.A. Association of Biomolecular Systems via Pulsed Field Gradient NMR Self-Diffusion Measurements. J. Am. Chem. Soc. 1995, 117, 7566–7567. [Google Scholar] [CrossRef]
  35. Binolfi, A.; Rasia, R.M.; Bertoncini, C.W.; Ceolin, M.; Zweckstetter, M.; Griesinger, C.; Jovin, T.M.; Fernández, C.O. Interaction of α-Synuclein with Divalent Metal Ions Reveals Key Differences: A Link between Structure, Binding Specificity and Fibrillation Enhancement. J. Am. Chem. Soc. 2006, 128, 9893–9901. [Google Scholar] [CrossRef]
  36. Bertoncini, C.W.; Rasia, R.M.; Lamberto, G.R.; Binolfi, A.; Zweckstetter, M.; Griesinger, C.; Fernandez, C.O. Structural characterization of the intrinsically unfolded protein beta-synuclein, a natural negative regulator of alpha-synuclein aggregation. J. Mol. Biol. 2007, 372, 708–722. [Google Scholar] [CrossRef]
  37. Morar, A.S.; Olteanu, A.; Young, G.B.; Pielak, G.J. Solvent-induced collapse of alpha-synuclein and acid-denatured cytochrome c. Protein Sci. 2001, 10, 2195–2199. [Google Scholar] [CrossRef]
  38. Davis-Searles, P.R.; Saunders, A.J.; Erie, D.A.; Winzor, D.J.; Pielak, G.J. Interpreting the effects of small uncharged solutes on protein-folding equilibria. Annu. Rev. Biophys. Biomol. Struct. 2001, 30, 271–306. [Google Scholar] [CrossRef] [Green Version]
  39. Wilkins, D.K.; Grimshaw, S.B.; Receveur, V.; Dobson, C.M.; Jones, J.A.; Smith, L.J. Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques. Biochemistry 1999, 38, 16424–16431. [Google Scholar] [CrossRef]
  40. Allison, J.R.; Rivers, R.C.; Christodoulou, J.C.; Vendruscolo, M.; Dobson, C.M. A relationship between the transient structure in the monomeric state and the aggregation propensities of α-synuclein and β-synuclein. Biochemistry 2014, 53, 7170–7183. [Google Scholar] [CrossRef]
  41. Maltsev, A.S.; Ying, J.; Bax, A. Impact of N-terminal acetylation of alpha-synuclein on its random coil and lipid binding properties. Biochemistry 2012, 51, 5004–5013. [Google Scholar] [CrossRef] [PubMed]
  42. Wise-Scira, O.; Dunn, A.; Aloglu, A.K.; Sakallioglu, I.T.; Coskuner, O. Structures of the E46K mutant-type alpha-synuclein protein and impact of E46K mutation on the structures of the wild-type alpha-synuclein protein. ACS Chem. Neurosci. 2013, 4, 498–508. [Google Scholar] [CrossRef] [Green Version]
  43. Hollingsworth, S.A.; Karplus, P.A. A fresh look at the Ramachandran plot and the occurrence of standard structures in proteins. Biomol. Concepts 2010, 1, 271–283. [Google Scholar] [CrossRef] [PubMed]
  44. Dusa, A.; Kaylor, J.; Edridge, S.; Bodner, N.; Hong, N.-P.; Fink, A.L. Characterization of oligomers during α-synuclein aggregation using intrinsic tryptophan fluorescence. Biochemistry 2006, 45, 2752–2760. [Google Scholar] [CrossRef]
  45. Guerrero-Ferreira, R.; Taylor, N.M.I.; Mona, D.; Ringler, P.; Lauer, M.E.; Riek, R.; Britschgi, M.; Stahlberg, H. Cryo-EM structure of alpha-synuclein fibrils. Elife 2018, 7, e36402. [Google Scholar] [CrossRef] [PubMed]
  46. Li, Y.; Zhao, C.; Luo, F.; Liu, Z.; Gui, X.; Luo, Z.; Zhang, X.; Li, D.; Liu, C.; Li, X. Amyloid fibril structure of alpha-synuclein determined by cryo-electron microscopy. Cell Res. 2018, 28, 897–903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Meade, R.M.; Fairlie, D.P.; Mason, J.M. Alpha-synuclein structure and Parkinson’s disease—Lessons and emerging principles. Mol. Neurodegener 2019, 14, 29. [Google Scholar] [CrossRef] [Green Version]
  48. Der-Sarkissian, A.; Jao, C.C.; Chen, J.; Langen, R. Structural organization of alpha-synuclein fibrils studied by site-directed spin labeling. J. Biol. Chem. 2003, 278, 37530–37535. [Google Scholar] [CrossRef] [Green Version]
  49. Heise, H.; Hoyer, W.; Becker, S.; Andronesi, O.C.; Riedel, D.; Baldus, M. Molecular-level secondary structure, polymorphism, and dynamics of full-length alpha-synuclein fibrils studied by solid-state NMR. Proc. Natl. Acad. Sci. USA 2005, 102, 15871–15876. [Google Scholar] [CrossRef] [Green Version]
  50. Johnson, M.; Geeves, M.A.; Mulvihill, D.P. Production of amino-terminally acetylated recombinant proteins in E. coli. Methods Mol. Biol. 2013, 981, 193–200. [Google Scholar] [CrossRef] [Green Version]
  51. Hoyer, W.; Cherny, D.; Subramaniam, V.; Jovin, T.M. Impact of the acidic C-terminal region comprising amino acids 109-140 on alpha-synuclein aggregation in vitro. Biochemistry 2004, 43, 16233–16242. [Google Scholar] [CrossRef]
  52. Cavanagh, J.; Fairbrother, W.; Palmer, A., III; Skelton, N. Protein NMR Spectroscopy: Principles and Practice; Academic Press: Cambridge, MA, USA, 1995; ISBN 9780080515298. [Google Scholar]
  53. Vuister, G.W.; Bax, A. Quantitative J correlation: A new approach for measuring homonuclear three-bond J(HNHA) coupling constants in 15N-enriched proteins. J. Am. Chem. Soc. 1993, 115, 7772–7777. [Google Scholar] [CrossRef]
  54. Hill, R.B.; Flanagan, J.M.; Prestegard, J.H. 1H and 15N magnetic resonance assignments, secondary structure, and tertiary fold of Escherichia coli DnaJ(1-78). Biochemistry 1995, 34, 5587–5596. [Google Scholar] [CrossRef] [PubMed]
  55. Serrano, L. Comparison between the phi distribution of the amino acids in the protein database and NMR data indicates that amino acids have various phi propensities in the random coil conformation. J. Mol. Biol. 1995, 254, 322–333. [Google Scholar] [CrossRef] [PubMed]
  56. Rossetti, G.; Musiani, F.; Abad, E.; DiBenedetto, D.; Mouhib, H.; Fernández, C.O.; Carloni, P. Conformational ensemble of human alpha-synuclein physiological form predicted by molecular simulations. Phys. Chem. Chem. Phys. 2016, 18, 5702–5706. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Daura, X.; Gademann, K.; Jaun, B.; Seebach, D.; Van Gunsteren, W.F.; Mark, A.E. Peptide Folding: When Simulation Meets Experiment. Angew. Chem. Int. Ed. 1999, 38, 236–240. [Google Scholar] [CrossRef]
  58. Páll, S.; Abraham, M.; Kutzner, C.; Hess, B.; Lindahl, E. Tackling exascale software challenges in molecular dynamics simulations with GROMACS. In Solving Software Challenges for Exascale; Markidis, S., Laure, E., Eds.; Springer International Publishing: Cham, Switzerland, 2015; pp. 3–27. [Google Scholar]
  59. Schrödinger. Schrödinger Release 2019-3: Maestro. 2019: New York, NY. Available online: https://www.schrodinger.com (accessed on 15 July 2020).
  60. Tribello, G.A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED 2: New feathers for an old bird. Comput. Phys. Commun. 2014, 185, 604–613. [Google Scholar] [CrossRef] [Green Version]
  61. Robustelli, P.; Piana, S.; Shaw, D.E. Developing a molecular dynamics force field for both folded and disordered protein states. Proc. Natl. Acad. Sci. USA 2018, 115, E4758–E4766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089. [Google Scholar] [CrossRef] [Green Version]
  63. Hess, B.; Bekker, H.; Berendsen, H.J.; Fraaije, J.G. LINCS: A linear constraint solver for molecular simulations. J. Comp. Chem. 1998, 18, 9. [Google Scholar] [CrossRef]
  64. Nosé, S. A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 2006, 52, 255–268. [Google Scholar] [CrossRef]
  65. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A Gen. Phys. 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [Green Version]
  66. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Chem. Phys. 1998, 52, 7182–7190. [Google Scholar] [CrossRef]
  67. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  68. Mavridis, L.; Janes, R.W. PDB2CD: A web-based application for the generation of circular dichroism spectra from protein atomic coordinates. Bioinformatics 2017, 33, 56–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Han, B.; Liu, Y.; Ginzinger, S.W.; Wishart, D.S. SHIFTX2: Significantly improved protein chemical shift prediction. J. Biomol. NMR 2011, 50, 43–57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Kabsch, W.; Sander, C. Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22, 2577–2637. [Google Scholar] [CrossRef]
Ijms 21 05061 g001 550
Figure 1. The position-39 αS variants are largely disordered in the monomeric state. The far-UV circular dichroism (CD) spectra of wt (black), Y39F (blue), Y39A (red), and Y39L (green) αS monomers show spectral features that are characteristic of random coil conformations (Left). In comparison, the predicted CD from the molecular ensemble for the four variants is also shown (Right).
Figure 1. The position-39 αS variants are largely disordered in the monomeric state. The far-UV circular dichroism (CD) spectra of wt (black), Y39F (blue), Y39A (red), and Y39L (green) αS monomers show spectral features that are characteristic of random coil conformations (Left). In comparison, the predicted CD from the molecular ensemble for the four variants is also shown (Right).
Ijms 21 05061 g001
Ijms 21 05061 g002 550
Figure 2. 1H NMR of aromatic side chains of the position-39 αS proteins. Spectra were registered at 15 °C in Buffer A of samples containing 100 µM wt (A), Y39F (B), Y39A (C), and Y39L αS (D). Asterisks indicate the lack of peaks corresponding to the Tyr-39 residue. Circles indicate the detection of new signals assigned to the Phe residue in position 39.
Figure 2. 1H NMR of aromatic side chains of the position-39 αS proteins. Spectra were registered at 15 °C in Buffer A of samples containing 100 µM wt (A), Y39F (B), Y39A (C), and Y39L αS (D). Asterisks indicate the lack of peaks corresponding to the Tyr-39 residue. Circles indicate the detection of new signals assigned to the Phe residue in position 39.
Ijms 21 05061 g002
Ijms 21 05061 g003 550
Figure 3. Changes induced by the mutations are only local and restricted to the immediate vicinity of the position-39 residue. (A) Overlaid contour plots of 1H-15N heteronuclear single quantum correlation (HSQC) spectra of 50 µM wt (black) and Y39L αS (red). Most affected residues located in the vicinity of the mutated site are labeled. (B) Differences in the mean weighted chemical shifts (MW 1H-15N ΔCS) displacements between wt and Y39F, Y39A, and Y39L αS mutants. (C) Panel shows the 3JHN-Hα profiles measured for the wt and Y39F, Y39A, and Y39L αS proteins. In all cases, experiments were recorded at 15 °C using 15N isotopically enriched protein samples (50 µM) dissolved in Buffer A.
Figure 3. Changes induced by the mutations are only local and restricted to the immediate vicinity of the position-39 residue. (A) Overlaid contour plots of 1H-15N heteronuclear single quantum correlation (HSQC) spectra of 50 µM wt (black) and Y39L αS (red). Most affected residues located in the vicinity of the mutated site are labeled. (B) Differences in the mean weighted chemical shifts (MW 1H-15N ΔCS) displacements between wt and Y39F, Y39A, and Y39L αS mutants. (C) Panel shows the 3JHN-Hα profiles measured for the wt and Y39F, Y39A, and Y39L αS proteins. In all cases, experiments were recorded at 15 °C using 15N isotopically enriched protein samples (50 µM) dissolved in Buffer A.
Ijms 21 05061 g003
Ijms 21 05061 g004 550
Figure 4. Predicted versus experimental chemical shifts of wt αS. Scatterplots showing the comparison between experimental chemical shifts of C, C, Cβ, and N atoms and those calculated from our replica exchange with solute tempering (REST2) simulations using the SHIFTX2 code (see Methods for details).
Figure 4. Predicted versus experimental chemical shifts of wt αS. Scatterplots showing the comparison between experimental chemical shifts of C, C, Cβ, and N atoms and those calculated from our replica exchange with solute tempering (REST2) simulations using the SHIFTX2 code (see Methods for details).
Ijms 21 05061 g004
Ijms 21 05061 g005 550
Figure 5. Conformational preferences observed for wt αS and its Y39A, Y39F, and Y39L mutants. Density maps of the radius of gyration (RG) and end-to-end (EE) distance reaction coordinates for the simulated structural ensembles.
Figure 5. Conformational preferences observed for wt αS and its Y39A, Y39F, and Y39L mutants. Density maps of the radius of gyration (RG) and end-to-end (EE) distance reaction coordinates for the simulated structural ensembles.
Ijms 21 05061 g005
Ijms 21 05061 g006 550
Figure 6. Hydrogen bonds (HB) and hydrophobic contacts (HC) of the four species investigated here. First panel from the left: distribution of HB contacts, with average values of 26.5, 29.6, 32.4, and 29.6 for wt, Y39F, Y39L, and Y39A αS, respectively. Second panel: HB contacts in each protein region. An HB contact was assumed to be present if persistent for 8% of the simulation or longer. Third panel: Distribution of HC contacts, with average values at 53.0, 54.8, 53.8, and 52.7 for wt, Y39F, Y39L, and Y39A αS, respectively. Fourth panel: HC contacts in each protein region. The HC contacts were assumed to be present if persistent for 15% of the simulation or longer. In the second and the fourth panels, error bars are reported. They were computed from the standard deviations of the distributions in the first and the third panels, respectively.
Figure 6. Hydrogen bonds (HB) and hydrophobic contacts (HC) of the four species investigated here. First panel from the left: distribution of HB contacts, with average values of 26.5, 29.6, 32.4, and 29.6 for wt, Y39F, Y39L, and Y39A αS, respectively. Second panel: HB contacts in each protein region. An HB contact was assumed to be present if persistent for 8% of the simulation or longer. Third panel: Distribution of HC contacts, with average values at 53.0, 54.8, 53.8, and 52.7 for wt, Y39F, Y39L, and Y39A αS, respectively. Fourth panel: HC contacts in each protein region. The HC contacts were assumed to be present if persistent for 15% of the simulation or longer. In the second and the fourth panels, error bars are reported. They were computed from the standard deviations of the distributions in the first and the third panels, respectively.
Ijms 21 05061 g006
Ijms 21 05061 g007 550
Figure 7. Structural characterization of the interaction between phthalocyanine tetrasulfonate and the Y39 αS mutants monitored by NMR spectroscopy. I/I0 profiles of backbone amide groups of 50 µM wt, Y39F, Y39L, and Y39A αS proteins in the presence of 50 µM phthalocyanine tetrasulfonate (PcTS). Binding sites for the αS-PcTS complexes at the N-terminal region of the protein are shown. 1H–15N HSQC spectra were recorded at 15 °C using 15N isotopically enriched αS samples (50 μM) dissolved in Buffer A.
Figure 7. Structural characterization of the interaction between phthalocyanine tetrasulfonate and the Y39 αS mutants monitored by NMR spectroscopy. I/I0 profiles of backbone amide groups of 50 µM wt, Y39F, Y39L, and Y39A αS proteins in the presence of 50 µM phthalocyanine tetrasulfonate (PcTS). Binding sites for the αS-PcTS complexes at the N-terminal region of the protein are shown. 1H–15N HSQC spectra were recorded at 15 °C using 15N isotopically enriched αS samples (50 μM) dissolved in Buffer A.
Ijms 21 05061 g007
Ijms 21 05061 g008 550
Figure 8. Structural features of fibrillated wt αS proteins. Boxplots showing the distribution of values for ring-to-ring distances (left) and χ2 dihedral angles (right) for fibrillated αS proteins in the PDB for Tyr-39.
Figure 8. Structural features of fibrillated wt αS proteins. Boxplots showing the distribution of values for ring-to-ring distances (left) and χ2 dihedral angles (right) for fibrillated αS proteins in the PDB for Tyr-39.
Ijms 21 05061 g008
Table
Table 1. Number of residues forming different secondary structure elements, as predicted by our REST2 simulation. The averages values and the corresponding standard deviations are reported.
Table 1. Number of residues forming different secondary structure elements, as predicted by our REST2 simulation. The averages values and the corresponding standard deviations are reported.
Secondary Structure wtY39FY39LY39A
Coil72 ± 667 ± 659 ± 668 ± 6
Bend30 ± 531 ± 432 ± 525 ± 5
Turn19 ± 523 ± 625 ± 524 ± 6
α-Helix11 ± 48 ± 511 ± 411 ± 4
3-Helix5 ± 47 ± 47 ± 46 ± 4
β-Sheet2 ± 41 ± 23 ± 32 ± 2
β-Bridge1 ± 23 ± 24 ± 24 ± 2

Share and Cite

MDPI and ACS Style

Palomino-Hernandez, O.; Buratti, F.A.; Sacco, P.S.; Rossetti, G.; Carloni, P.; Fernandez, C.O. Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly. Int. J. Mol. Sci. 2020, 21, 5061. https://doi.org/10.3390/ijms21145061

AMA Style

Palomino-Hernandez O, Buratti FA, Sacco PS, Rossetti G, Carloni P, Fernandez CO. Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly. International Journal of Molecular Sciences. 2020; 21(14):5061. https://doi.org/10.3390/ijms21145061

Chicago/Turabian Style

Palomino-Hernandez, Oscar, Fiamma A. Buratti, Pamela S. Sacco, Giulia Rossetti, Paolo Carloni, and Claudio O. Fernandez. 2020. "Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly" International Journal of Molecular Sciences 21, no. 14: 5061. https://doi.org/10.3390/ijms21145061

APA Style

Palomino-Hernandez, O., Buratti, F. A., Sacco, P. S., Rossetti, G., Carloni, P., & Fernandez, C. O. (2020). Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly. International Journal of Molecular Sciences, 21(14), 5061. https://doi.org/10.3390/ijms21145061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop