Charge Regulation during Amyloid Formation of α-Synuclein.

Electrostatic interactions play crucial roles in protein function. Measuring pKa value perturbations upon complex formation or self-assembly of e.g. amyloid fibrils gives valuable information about the effect of electrostatic interactions in those processes. Site-specific pKa value determination by solution NMR spectroscopy is challenged by the high molecular weight of amyloid fibrils. Here we report a pH increase during fibril formation of α-synuclein, observed using three complementary experimental methods: pH electrode measurements in water; colorimetric changes of a fluorescent indicator; and chemical shift changes for histidine residues using solution state NMR spectroscopy. A significant pH increase was detected during fibril formation in water, on average by 0.9 pH units from 5.6 to 6.5, showing that protons are taken up during fibril formation. The pH upshift was used to calculate the average change in the apparent pKaave value of the acidic residues, which was found to increase by at least 1.1 unit due to fibril formation. Metropolis Monte Carlo simulations were performed on a comparable system that also showed a proton uptake due to fibril formation. Fibril formation moreover leads to a significant change in proton binding capacitance. Parallel studies of a mutant with five charge deletions in the C-terminal tail revealed a smaller pH increase due to fibril formation, and a smaller change (0.5 units on average) in the apparent pKaave values of the acidic residues. We conclude that the proton uptake during the fibril formation is connected to the high density of acidic residues in the C-terminal tail of α-synuclein.

Introduction

Seven out of the 20 amino acids, Asp, Glu, Lys, Arg, His, Tyr, Cys, as well as the N and C termini, contain ionizable groups, giving rise to the acid/base properties of proteins.[1] In 1917 Sørensen pioneered the research on protein titrations, demonstrating that egg albumin is an ampholyte with both acid and base properties.[2] The ionizable side chain groups (as well as the terminal groups) are the major contributors to the electrostatic interactions of proteins, affecting the structure, solubility, stability, and function, such as the kinetic properties.[3−8]

Proteins are large heteropolymeric polyelectrolytes that can carry many charges simultaneously. The sum of the charged state of the ionizable side chain groups, which in turn depends on interactions and the pKa values of the individual groups, defines the net charge of the protein. Linderstøm-Lang discovered that the protein’s net charge can affect and perturb the pKa values of its titratable groups.[9,10] Since then, pKa perturbations of amino acid residues in proteins have been studied extensively and are frequently found to play important roles in protein function and folding.[11] The quantification of the perturbations of pKa values of individual residues from their intrinsic pKa values thus provide valuable information about the their environment and conformational transitions upon folding and binding.[12]

The pKa values of individual side chain groups can be down-shifted or up-shifted relative to their intrinsic pKa values as measured for model peptides[4,8,13,14] (Table ). The three main factors that contribute to pKa perturbations are Coulomb interactions (charge–charge interactions), charge–dipole interactions, and desolvation.[8,15] pKa values of ionizable groups free in solutions are generally different from the pKa value of the same group buried in a protein interior due to differences in the polarity of the environment (desolvation). The pH dependence of enzyme catalysis can be tailored by modifications of pKa values and altered Coulomb interactions, leading to changes in the electrostatic environment within the active site.[16] The pKa values of ionizable groups within an enzyme’s active sites change as a result of altered electrostatic interactions with other adjacent ionizable groups.[15] Studying pKa perturbations within enzymes’ active sites has given an increased understanding of the role of ionizable residues within enzymatic reactions.[15,17−20]

Table 1

Comparison of the Intrinsic pKa Values of Aspartic Acid and Glutamic Acid to the pKa Values Observed in Small Model Peptides, Folded Proteins, and α-Synuclein Monomer

	Asp	Glu
Intrinsic pKaa	4.0	4.4
pKa in alanine pentapeptideb	3.9	4.3
pKa in random coil model peptidec	3.9 ± 0.1	4.3 ± 0.1
pKa in folded proteinsb	3.5 ± 1.2	4.2 ± 0.9
pKa in full length α-synuclein monomerd	4.0 ± 0.4	4.4 ± 0.3
pKa in N-terminal and NAC region of α-synucleind	3.6 ± 0.05	4.2 ± 0.1
pKa in C-terminal tail of α-synucleind	4.2 ± 0.3	4.6 ± 0.2

From ref (42).

From ref (8).

From ref (43); pKa of the ionizable amino acid residue x of the c-H-Gly-Gly-X-l-Ala-OH peptide.

Average of published values from reference in absence of salt.[30]

Large pKa value shifts are often found for ionizable residues with low solvent accessible surface area (SASA). A study by Pahari and Alexov (2019) showed that the pKa values of Asp with less than 40% SASA range from 0.5 to 9.9. For Glu residues with 30% or less SASA, they found pKa values lower than 2.6 and higher than 6.[15,21] Charged residues on protein surfaces, with greater SASA, are also commonly found to shift, but the degree of perturbations is generally less than 2 units.[15] In the same study Pahari and Alexov found that 68% of the Asp residues have more than 40% SASA, and pKa values range between approximately 2 and 5.3. Similarly, about 75% of the Glu residues have more than 40% SASA and pKa values that range from 2.6 to 6.[21] High pKa values, up to 6.5, are found for Glu residues in a Ca2+-binding site of S100G, also known as calbindin D9k.[22] Perturbations of residues on protein surfaces are generally due to Coulomb interactions (charge–charge interactions).[23] Proteins are evolved for function, which requires a compromise between solubility and stability, which are affected by protein net charge in opposite ways. Electrostatic interaction between groups on protein surfaces[24] as well as within unstructured proteins contribute to the protein stability, where favorable interactions can increase the stability of the protein state, and very commonly, unfavorable electrostatic repulsion decreases the protein stability.[25,26]

The pKa values of individual titratable groups within proteins can be perturbed due to interactions with nearby residues, or due to interactions with neighboring biomolecules or ligands. For example, electrostatic interactions play important roles in the stability and formation of large self-assembling protein structures, such as the tobacco mosaic virus protein and amyloid fibrils.[6,27,13]

α-Synuclein is a 140 residue long and intrinsically disordered protein, for which the formation of amyloid fibrils and their accumulation into Lewy bodies is associated with Parkinson’s disease.[28] The protein sequence is commonly divided into three regions: the N-terminal (amphipathic) region, the central hydrophobic region, also called the NAC region, and the C-terminal (acidic) tail.[29,30] However, the boundaries of the regions may vary depending on the context, as suggested from solved structures.[31] Here, the protein will be divided into the N-terminal tail, the hydrophobic fibril core, and the C-terminal tail (Figure ) based on how the monomeric protein arranges into a fibrillar structure as well as the distribution of charge residues within the sequence.[31,30,29] The last 40 residues of α-synuclein—referred to as the C-terminal tail—are the focus of this project. The tail is highly acidic, consisting of 15 acidic groups (including the C-terminus), and is highly charged, with a net charge of −12 at physiological pH.[32] Previous studies have indicated that the C-terminal tail is not a part of the fibril core and remains unstructured in both the monomeric and fibrillar form.[31,30,33−36] The C-terminal tail is believed to extend from the fibril core like a polymer brush at conditions where it is highly negatively charged.[33]

Figure 1

Sequence and structural properties of α-synuclein fibrils. (A) Distribution of acidic (red) and basic (blue) residues within the α-synuclein sequence. The sequence can be divided into three regions according to how it forms amyloid-fibrils: the fibril core (residues 29–100); the N-terminus (residues 1–28); and the C-terminal tail (101–140) (according to 2NOA.pdb). The C-terminal tail of the 5Q mutant is shown in comparison to the wild-type, where five acidic residues have been mutated to noncharged polar residues, making the tail less acidic. (B) Sequence of α-synuclein. Acidic residues are colored red, and basic residues are colored blue. The positions of the mutations (E114Q+E123Q+E126Q+E137Q+E139Q) of the 5Q mutant are shown below the sequence. (C) Proto-filament structure of full-length α-synuclein (pdb: 2N0A(31)) colored according to the 3 regions. The fibril core, purple; the N-terminal tail, cyan; the C-terminal tail, red. (D) Proto-filament structure of α-synuclein colored according to acidic (Red: Glu, Asp) and basic (Blue: Lys, His) residues. (E) A simplified illustration of α-synuclein fibrils, with the extending tails not being a part of the fibril core.

The mechanism of aggregation of α-synuclein into amyloid fibrils is highly dependent on pH, the process is dominated by surface catalyzed secondary nucleation below pH 6.0 and is much slower at higher pH. The pH dependence is believed to be linked to the high acidity of the C-terminal tail.[37,38] Truncation of the C-terminal tail, as well as charge-screening, has been found to affect the aggregation kinetics of α-synuclein.[37,39−41] Further studies on the effect of the acidic residues in the tail of α-synuclein on the amyloid formation are important in order to gain a deeper understanding of the pH dependence of α-synuclein aggregation and the role of electrostatic interactions in α-synuclein amyloid formation.

The pKa values of the acidic residues of the monomeric α-synuclein have previously been measured using NMR spectroscopy and were found to be higher in the C-terminal tail compared to the first 100 residues (Table ).[30] These pKa value perturbations of the acidic residues in the C-terminal tail could be explained by the high charge density and the close proximity between negative charges within the tail.[37,30]Table gives a comparison of intrinsic pKa values of the acidic residues, Asp and Glu, to their pKa values measured within small model peptides, folded proteins, and α-synuclein monomer.

In this study, the effect of amyloid formation on the pKa values of the acidic residues in the C-terminal tail of α-synuclein was examined. The goal was to compare the pKa values of α-synuclein monomers and fibrils, which could provide a better understanding of the role of electrostatic interactions and the effect of the C-terminal tail on the fibril formation process. The main questions were as follows: Do the pKa values of the acidic residues become upshifted during amyloid formation? If so, is that a result of the high density of acidic residues in the unstructured C-terminal tail? These questions were addressed using three complementary experimental methods to monitor the change in pH upon fibril formation in water as well as through Metropolis Monte Carlo computer simulations. The data for wild-type α-synuclein were compared to a mutant lacking five acidic residues in the C-terminal tail (Figure a,b), to provide a relative measure of the effect of the negative charges in the tail.

Experimental Section

Protein Expression

Escherichia coli (E. coli) BL21* pLysS Ca2+ competent cells were transformed with pET-3a-plasmid containing the gene for wild-type human α-synuclein (not N-acetylated) with E. coli-optimized codons and an ATG start codon (purchased from GenScript, Piscataway, New Jersey; see Supporting Information). The same was done for the mutated α-synuclein used in this study. 0.7 μL of plasmid (100 ng/μL) was carefully mixed with 30 to 40 μL of Ca2+-competent cells and kept on ice for 30 to 60 min, and the sample was heat-pulsed for 45 s at 42 °C and placed on ice for 10 min. The transformed cells were spread on sterile LB agar plates supplemented with chloramphenicol (30 μg/mL) and ampicillin (50 μg/mL) and incubated at 37 °C overnight (ON). The plate was stored at 5 °C for 8 h before a small single colony was selected for inoculation of 50 mL culture (LB medium with chloramphenicol (30 μg/mL) and ampicillin (50 μg/mL)) which was incubated ON at 37 °C with shaking. Next morning, 5 mL overnight culture was transferred to each 500 mL day culture (LB medium with chloramphenicol (30 μg/mL) and ampicillin (50 μg/mL)) in 2.5 L baffled flasks. The day culture was incubated at 37 °C with continuous shaking at 125 rpm. The optical density at 600 nm (OD600 nm) was followed, and when it had reached approximately 0.9–1.0, protein expression was induced with 100 μg/mL isopropyl thio-β-d-galactoside (IPTG). The cells were harvested 4 h after induction, by centrifugation at 6000g for 12 min at 4 °C (JA 8.100 rotor). The supernatant was disposed, and the cells obtained from total of 4 L culture were combined, mixed with 25 mL of water, and frozen. Before harvesting, 1 mL samples were taken from the cultures for testing the expression by SDS-PAGE.

Expression of 15N labeled wild-type α-synuclein was done in M9 minimal medium. In this case, 3 mL of overnight culture in LB medium with 15NH4Cl as the sole nigrogen source, (supplemented with chloramphenicol (30 μL/mL) and ampicillin (50 μL/mL)) were added to 50 mL of middle-day culture in M9 minimal medium supplemented with chloramphenicol (30 μL/mL) and ampicillin (50 μL/mL). When OD600 nm had reached 1.0, then 50 mL of middle-day-culture were transferred to 450 mL of M9 minimal medium (supplemented with chloramphenicol (30 μL/mL) and ampicillin (50 μL/mL)). When D600 nm of the 500 mL culture had reached 0.7–0.9, the expression was induced with 100 μg/mL IPTG. The cells were harvested 4 h after induction as above.

Purification

Cell pellet from 8 L culture was thawed in 100–120 mL of cold buffer A (10 mM Tris/HCl, 1 mM EDTA, pH 7.5) and placed on ice. The volume of the buffer was approximately 5 times the volume of the cell pellet. The pellet was sonicated on ice, using pulse sonication (1 s on, 1 s off) until the mixture remained homogeneous. After sonication the sample was centrifuged for 10 min at 15000g at 4 °C (JA 25.50 rotor). The supernatant was collected (approximately 130 mL). The supernatant was poured into an equal volume of boiling buffer A. In order to precipitate E. coli proteins, the sample was heated to 85 °C. The sample was placed on ice and continuously stirred while cooling down. The sample was centrifuged at 15 000g for 10 min at 4 °C. The supernatant was collected and further purified by ion-exchange chromatography.

The first ion-exchange chromatography step was performed using 100 g of wet diethylaminoethyl (DEAE) cellulose in a column with a diameter of 3.5 cm. Before packing of the column, the resin was washed at least two times with Milli-Q water and two times with buffer A, or until the pH was measured to be around 7.5. All buffers were kept cold, and the purification was performed in a cold room. The column was equilibrated with 100 mL of buffer A. The sample (about 260 mL) was loaded onto the column which was then washed with 100 mL of buffer A. The sample was eluted at a flow rate of 1 mL/min using linear 0–0.5 M NaCl gradient in buffer A (total gradient volume of 1.4 L). SDS-PAGE or agarose gel electrophoresis was used to determine which fractions contained α-synuclein. The fractions containing α-synuclein were pooled, diluted 50:50 with buffer A, and further purified on a DEAE sephacel column (60 g of wet resin in a column with diameter of 2.3 cm). The second ion-exchange chromatography step was performed in the same way as above for the DEAE cellulose column. The absorbance at 280 nm was measured and further analyzed with SDS-PAGE to determine which fraction contained α-synuclein (see appendix). The fractions containing α-synuclein and no detectable impurities were pooled and stored as 1 mL aliquots at −20 °C. The concentration (in the range of 1–3 mg/mL) was determined by absorbance at 280 nm, using an extinction coefficient of ε = 5800 M–1cm–1. All buffers and milli-Q water were filtrated and degassed before use. The corrected mass was confirmed, and the purity was tested using MALDI-TOF and MALDI-TOF-TOF mass spectrometry.

The expression and purification of the α-synuclein 5Q mutant was performed in the same way as for the wild-type, except for the sample being eluted using a linear 0–0.35 M NaCl gradient in buffer instead of 0–0.5 M, in both the first and second ion-exchange chromatography step.

Preparation of Samples

Prior to each experiment, monomeric α-synuclein was isolated by size exclusion chromatography (SEC) in the desired buffer. When buffer free protein was needed, the sample was subsequently loaded onto a desalting column in order to exchange the buffer for milli-Q water in order to get α-synuclein in water.

Size-Exclusion Chromatography

To further concentrate the sample, the fractions containing α-synuclein were lyophilized and dissolved in a smaller volume of 6 M guanidinium hydrochloride. Prior to loading the sample onto the SEC column, the sample was incubated at room temperature for at least 1 h to ensure that the protein was properly dissolved. The sample was injected onto a Superdex 75 Increase 10/300 GL (GE Healthcare) column using a fast protein liquid chromatography (FPLC) system (Bio-RAD, BIOLOgic Duo Flow, USA). The sample was eluted at 0.7 mL/min in the desired buffer. The center of the monomeric peak (1–1.5 mL) was collected into low-binding tubes (Genuine Axygen Quality). The absorbance at 280 nm was used to follow the elution of the monomeric α-synuclein. The concentration was determined from the absorbance at 280 nm using an extinction coefficient of ε = 5800 M–1 cm–1. All buffer solutions were prepared fresh, degassed, and filtered before usage.

Desalting Column

To obtain α-synuclein in pure water, without any buffer, a HiTrap Desalting column (GE Healtcare) was used to separate α-synuclein from buffer. The sample (1 mL sample collected from precious SEC chromatography step) was injected onto the desalting column and eluted at 1 mL/min in degassed and filtrated Milli-Q water using a FPLC system. The center of the protein peak was collected into low-binding tubes and used for further experiments. The concentration was determined from the absorbance at 280 nm using an extinction coefficient of ε = 5800 M–1 cm–1.

Fibril Formation Kinetics

Monomeric α-synuclein was prepared the same day as the experiment started and kept on ice until the beginning of the experiment. The α-synuclein fibrils were formed by incubating 100 μL of monomeric α-synuclein in a nontreated 96 well polystyrene black plate (3631 Costar) under quiescent conditions. The plate had a transparent bottom and was sealed with a transparent SealPlate film in order to avoid evaporation. The fibril formation was indirectly measured by supplementing the samples with 20 μM Thioflavin T (ThT), and the fluorescence was monitored over time using a Fluostar plate reader (BMG Labtech, Offenburg, Germany) at 37 °C without shaking using excitation and emission wavelengths of 440 and 480 nm, respectively.

Confirming the Presence of Amyloid Fibrils

Far-UV Circular Dichroism (CD) Spectroscopy

The Far-UV CD spectra were recorded between 260 and 190 nm at 20 °C in a 0.1 mm quartz cuvette using a Jasco J-815 CD spectrometer. The spectra were acquired with a scanning speed of 50 nm/min, continuous scanning mode, digital integration time per data point (D.I.T) of 8 s, and sensitivity set to standard. The data were averaged over three accumulations.

ANS Fluorescence

The fluorescence of the optical probe, 8-Anilino-1-napthalenesulfonate (ANS), is highly sensitive to the environment and can easily increase about 20 times in a nonpolar or rigid environment.[44] Binding of ANS to hydrophobic patches of amyloid fibrils results in a significant increase in fluorescence intensity.[45] For validation of the presence of fibrils, ANS (10 μM) was supplemented to aliquots of the samples of interest, before and after fibril formation. The samples were excited at 395 nm, and the emission spectra were recorded from 410 nm to 600 (PerkinElmer luminescence Spectrometer LS-50B, UK).

Detection of pH Changes

Detection of pH Change Using pH Meter

After isolation of monomeric samples in pure water using a desalting column, the samples were diluted to 20 or 70 μM α-synuclein in filtrated Milli-Q water. The pH of an aliquot from the sample was measured with an Orion star A211 pH meter and an Orion ROSS Combination pH Micro Electrode. The rest of the sample was incubated at 37 °C in a 2 mL low binding-tube (Genuine Axygen Quality), with continuous stirring at 700 rpm using a micro stirring bar (8 mm × 1.5 mm, Polytetrafluoroethylene (PTFE) coated, strong Alnico V magnetic core, round smooth surface). The presence of fibrils was confirmed using CD spectroscopy and/or ANS fluorescence using a withdrawn aliquot. The pH of the sample containing amyloid fibrils was measured using the same pH electrode.

Continuous pH Measurement during Fibril Formation Using Probe Drum

The pH was measured continuously using an Orion star A211 pH meter connected to an Orion ROSS Combination pH Micro Electrode and a Probe Drum instrument (Probation Laboratories Sweden AB). A 1 mL aliquot of 70 μM α-synuclein prepared in water was placed in a Quartz SUPRASIL (HellmaAnalytics) cuvette with a path length of 10 mm. A micro stirring bar (8 mm × 1.5 mm, PTFE coated) was added to the cuvette, and the sample was stirred at 237 rpm (setting nr. 4) at 37 °C and pH was measured continuously.

In some of these experiments, the pH, static light scattering, and ThT fluorescence were measured simultaneously using Probe Drum. The aggregation and pH measurements were performed the same way as explained above, except for including ThT in the sample (20 μM). Spectra were recorded from 300 to 730 nm at each time point. The ThT fluorescence was measured using a UV-light source of 392 nm, and the average fluorescence signal between 475 and 485 nm was plotted against time. The static light scattering was measured at 90° angle using a laser of 635 nm, and the average signal between 633 and 638 nm was plotted against time.

pH Detected Using pH Indicator

Monomeric α-synuclein was isolated in a weak buffer system consisting of 0.5 mM 2-(N-morpholino)ethanesulfonic acid (MES), 12.5 μM resazurin, and 0.02% NaN3 at pH 5.5. Part of the monomeric sample was stored at 5 °C to prevent fibril formation; this was done for later color comparison between monomeric and fibrillar sample. The other part of the sample was incubated in 2 mL low-binding tubes at 37 °C with stirring at 900 rpm (IKA IKAMAG REO Magnetic Stirrer) with a micro stirring bar (8 mm × 1.5 mm, PTFE coated). After 1 day of incubation, the color of the incubated sample was compared, at room temperature, to the color of the monomeric sample. In addition, an aliquot of the samples was taken before and after the incubation at 37 °C, and the pH was measured using a pH electrode. Also far-UV CD (260–190 nm) spectroscopy was used to verify whether the samples were monomeric or fibrillar by checking what secondary structural component was dominant in the sample.

pH Change Detected Using NMR Spectroscopy by Using Histidine as a pH Sensor

The pH change during fibril formation was detected by comparing the chemical shifts of the δ2 proton in the imidazole ring of histidine before and after the fibrillation of α-synuclein. 70 μM α-synuclein was isolated in 99.9% D2O, without any buffer by using a two-step preparation, including SEC on a Superdex 75 column followed by SEC on a G25 desalting column. The sample was split into two parts, one to be kept monomeric and one to become fibrillar. The monomeric sample was stored on ice until the measurements started. The fibrillar sample was prepared by incubating the sample at 37 °C, with stirring at 700 rpm (IKA IKAMAG REO Magnetic Stirrer) with a micro stirring bar (0.8 cm × 0.3 cm, PTFE coated) for 2 days. The presence of fibrils was confirmed using far-UV CD spectroscopy. Solution state 1D 1H (proton) NMR spectra were acquired for the monomeric and fibrillar samples in 100% D2O at 300 K using an Agilent VNMRS DirectDrive spectrometer, operating at a 1H frequency of 499.9 MHz, equipped with a triple resonance probe. The 1H NMR spectra were acquired using a single-90°-pulse experiment. The water signal was suppressed by presaturation. A total of 256 transients of 16 k data points spanning a spectral width of 16 ppm were collected.

Measurement of Monomer Concentration after Fibril Formation

Wild-type and 5Q mutant were prepared in water as explained in sections and 2.3.2. The samples were diluted to 20 μM. An aliquot of the sample was stored on ice until further analysis with SDS-PAGE (10–20% Tris-Tricine gel). The samples of both variants were split into three tubes each, containing a 500 μL sample. The samples were incubated at 37 °C in 2 mL low-binding tubes (Genuine Axygen Quality), with continuous stirring at 700 rpm using a micro stirring bar (8 mm × 1.5 mm) for 2 days. The presence of fibrils was confirmed using CD spectroscopy, of a withdrawn aliquot. Two-step separation of fibrils and soluble monomers was performed as follows: 300 μL of each sample were centrifuged for 10 min at 14600g and room temperature in a 1.5 mL low-binding tube. Next, 150 μL of the supernatant were removed from the top of the sample and transferred to another 1.5 mL low-binding tube. The supernatant was centrifuged again for 25 min at 14 600g at room temperature. Thereafter, 50 μL were taken of the supernatant for further analysis with SDS-PAGE. CD spectroscopy was used to verify the separation of monomers and fibrils, by taking spectra of both supernatant and resuspended pellet.

Monomer (that had been stored on ice) and the supernatants of wild-type and 5Q mutant were loaded onto a polyacrylamide gel. The concentration of the monomer in equilibrium with fibrils was calculated by analyzing the intensity of the gel bands, using the program ImageJ. The intensity of each band was compared to that of a 6 μM monomeric sample.

Cryogenic Transmission Electron Microscopy (Cryo-TEM)

The fibrils were formed in a 96-well plate as previously explained. The protein concentration was 35 μM. The samples were supplemented with 20 μM ThT, and its fluorescence was recorded to ensure that the sample contained fibrils and had reached the final plateau before collection at 55 h. The fibrils tend to stick to the bottom of the well, and therefore the sample was collected by carefully scraping the bottom of the well, using an Eppendorf tip. Before the sample was frozen, it was carefully pipetted up and down. All samples were treated in the same way. To ensure a stable temperature and avoiding loss of solution during sample preparation, a controlled environment vitrification system was used. Samples were prepared as thin liquid films (<300 nm thick) on glow-discharged treated lacey carbon film coated copper grids and plunged into liquid ethane at −180 °C. In this way the original microstructures are preserved as we can avoid component segmentation and rearrangement in addition to water crystallization as the samples are vitrified. Samples were stored under liquid N2 until measured and then transferred using an Oxford CT3500 cryoholder and its workstation into the electron microscope (Philips CM120 BioTWIN Cryo) equipped with a postcolumn energy filter (Gatan GIF100). An acceleration voltage of 120 kV was used, and images were recorded digitally with a CCD camera under low electron dose conditions.

Diffusion Ordered NMR Spectroscopy

Diffusion ordered NMR spectroscopy (DOSY) was used to characterize the size of the NMR detected molecular species. Here the signal intensity of the NAC region of the protein monomer was followed as a function of gradient strength in both the absence and presence of fibrils. Measurements were done in 100% D2O at 300 K on a Bruker Avance III HD 800 spectrometer (Bruker Biospin, Rheinstetten, Germany), operating at a 1H frequency of 799.9 MHz, equipped with a 5 mm cold probe using the ledbpgp2s1d sequence with stimulated echo and bipolar gradient pulses.[46] Spectra were recorded at four gradient strengths with a diffusion delay of 100 ms and a gradient pulse duration of 4 ms. A total of 32 transients of 32 k data points spanning a spectral width of 16 ppm were collected. The intensities were integrated over the regions 7.29–6.68, 4.59–3.70, 3.10–2.11, 2.00–1.28, and 0.90–0.72 ppm and normalized to the intensity at the lowest gradient strength.

Two-Dimensional Heteronuclear NMR Spectroscopy

2D 15N–1H heteronuclear single quantum coherence (HSQC) spectra were recorded in 5% D2O/95% H2O at 300 K on the same Bruker Avance III HD 800 spectrometer, using the hsqcetgpsi2 sequence. A total of 2048 × 128 data points spanning 14 × 30 ppm in the 1H and 15N dimensions were collected.

Molecular Modeling: Metropolis Monte Carlo Simulations

To study the ionization states of α-synuclein monomers and fibrils, constant pH Metropolis Monte Carlo simulations were performed using a coarse-grained protein model where protonation states are allowed to fluctuate according to their chemical environment. The model has previously been used to study α-synuclein[41] and presented in detail elsewhere (https://dx.doi.org/ 10/mrv). The use of constant pH in the simulations explicitly accounts for any coupling between titrating groups. All residues of the protein are included in the simulations. Amino acid residues are treated as single beads either connected with harmonic bonds (N-term 1–28, C-term 101–140) or kept fixed (residues 29–100) according to their mass center positions from the fibril structure found in PDB 2N0A. Amino acid beads interact through a combined Lennard–Jones and screened Coulomb potential where the latter uses a Bjerrum length of 7 Å and a screening length corresponding to 0.5 mM 1:1 salt. Residues participating in acid/base equilibria are subject to Monte Carlo charge swap moves mimicking (de)protonation and associated with a trial energy contribution of Δu = kBT(pH – pKa0) + Δuelec, where kB is Boltzmann’s constant, pKa0 is the intrinsic acid dissociation constant of the amino acid model compound, and Δuelec is the change in electrostatic energy. In addition to swap moves, configurational space in the NVT ensemble is sampled using atomic translations, chain pivot, and crankshaft moves at 300 K. Average protonation states are sampled over one million Monte Carlo sweeps where, in each sweep, N residue move attempts are carried out. Both a single monomer and a fibril (represented by ten planes, PDB: 2N0A) of α-synuclein are studied. The 2N0A structure was chosen as it includes all the 140 residues of α-synuclein. For the monomer, a single chain is simulated with all beads connected by harmonic springs. For the fibril, the middle part is kept as described above and connected with flexible N-terminal and C-terminal ends (see above and Figure a for a snapshot). All simulations were performed using the Faunus software,[47] and an electronic notebook (Jupyter Notebook) for reproducing all simulation results is provided. This allows the reader to investigate alternative fibril structures, add a finite free monomer concentration, or explore a change in any model parameter. For more information, see https://doi.org/10.5281/zenodo.4729860.

Figure 7

Results from Monte Carlo simulations of α-synuclein monomer and fibril (10 planes of a fibril as found in PDB: 2N0A) states. (A) Snapshot showing the rigid middle part of the fibril, flanked by flexible N-termini (blue) and C-termini (red). (B) Partial charge of carboxyl groups of glutamic acid and aspartic acid in the fibril (magenta) and in the monomer (green) as a function of pH. The dashed, black line shows the nonperturbed titration curve given by the model pKa0 value. (C) Differences in residue partial charges between α-synuclein in a monomer and in a fibril at constant pH 6.5. Shifts (charge in fibrillar minus charge in monomeric state) are indicated by the colorbars and letters correspond to the amino acid sequence with a rigid middle part surrounded by flexible (shown in “wavy” text) N-terminal and C-terminal parts. (D) Same as (C), but for the 5Q-mutant.

Results

Formation of α-Synuclein Fibrils in Water

The investigation of proton uptake by α-synuclein during fibril formation required that the reaction was followed in water, i.e. in the absence of buffer components except the protein. The kinetic behavior and the morphology of α-synuclein in water were therefore investigated and compared to fibrils formed in 10 mM MES/NaOH, pH 5.5. The fibril formation was followed by supplementing control samples with 20 μM ThT and monitoring the increase in ThT fluorescence (Figure ). The pH measurements below were for protein only in water without ThT.

Figure 2

Ultrastructure of α-synuclein fibrils formed in pure water and 10 mM MES/NaOH, pH 5.5. Fibrils were formed by incubating the monomeric α-synuclein in a nontreated polystyrene plate at quiescent conditions. Samples were supplemented with 20 μM ThT and the fluorescence intensity measured at different time points. The protein concentration was 35 μM. (A) Raw data. Red: 10 mM MES/NaOH (pH 5.5), t1/2 = 13.4 h, two replicates. Blue: in pure water, t1/2 = 38.0 h, four replicates. The bold dots show the average of the replicates. Individual replicates are shown with thinner dots. The t1/2 was calculated by fitting a sigmoidal curve to the average data. (B) Normalized data, same method and data as those in A. Cryo-TEM images of the α-synuclein fibril samples collected after reaching the plateau in ThT fluorescence. (C and D) α-Synuclein fibrils formed in water. (E and F) α-Synuclein fibrils formed in 10 mM MES, pH 5.5.

The aggregation kinetics (Figure a,b), indicate that amyloid fibrils form more slowly in pure water, with less steep transition and a t1/2 of 38.0 h, than in 10 mM MES/NaOH at pH 5.5 where t1/2 is 13.2 h.

The samples were harvested at the end of the aggregation reaction or after 56 h (Figure ). Cryogenic transmission electron microscopy (Cryo-TEM) was used to evaluate if there were any ultrastructural differences between the fibrils formed in water and in 10 mM MES/NaOH pH 5.5. A clear difference is seen in the morphology of the fibrils formed in water and in buffer (Figure c–f). The fibrils formed in water are generally longer, thinner, and more dispersed over the grid than those formed in buffer. Some of the fibrils are very long, up to few micrometres. Short and thin fibrils are commonly close to the mature fibrils. The fibrils in water are more separated with longer distances between the fibrils and the filaments than the fibrils in the pH 5.5 buffer, which are more compact and tangled, forming big dense aggregates. This is presumably a result of stronger electrostatic repulsion between the protofilaments in water compared to buffer where salt screening becomes significant.

Detection of pH Increase during Fibril Formation by pH Electrode

Using pure water as a solvent makes it possible to measure the effect of amyloid fibril formation on the pH of the sample and subsequently use the pH to calculate the average shift in pKa values. In water without buffer, any changes in pH of the sample depends only on changes in pKa values that occur upon aggregation. Alteration of the pKa values will change the number of free protons in solution, i.e. the pH of the sample. Therefore, a recorded change in pH within the sample (closed system) can be related to changes in pKa values during fibril formation. The pH of the sample was measured using a pH electrode before and after the formation of amyloid fibrils. An increase of as much as 0.9 pH units (from pH 5.6 to 6.5) was detected for 20 μM α-synuclein in water (Figure a). From this pH change, it is possible to calculate the change in the apparent pKaave values of α-synuclein during fibril formation (see below). The effect on pH is amplified at higher protein concentration (70 μM) (Table S1).

Figure 3

Change in pH during amyloid formation of wild-type α-synuclein and an α-synuclein 5Q mutant. (A) pH change during amyloid formation detected by pH electrode. The proteins (wild-type α-synuclein and the 5Q mutant) were isolated in pure water, and pH was measured before and after fibril formation. On average, pH changed from pH 5.6 to pH 6.5 for the wild-type protein and on average from pH 6.4 to 6.6 for the 5Q mutants. (B) pH change detected using the pH sensitive dye, resazurin. Monomeric α-synuclein in weak buffer containing resazurin, at pH 5.8, resulted in purple color. The fibrillar sample gave rise to blue color. For comparison, the pH was measured to be 6.4 using a pH electrode. (C) pH change detected using NMR spectroscopy using the δ2 proton of histidine in monomeric α-synuclein as a pH sensor. The red and blue traces show the spectrum before and after fibril formation, respectively. The large decrease in chemical shift detected for the monomer in the fibrillar sample is indicative of a pH increase.

To investigate if the pH increase is related to the high density of acidic residues in the C-terminal tail, the experiment was repeated for a mutant in which five acidic (glutamic acid) residues in the C-terminal tail had been replaced with noncharged polar (glutamine) residues (see Figure ), denoted the 5Q mutant. The increase in pH during fibrillation was only 0.2 pH units for the 5Q mutant (from pH 6.4 to 6.6; Figure a) at 20 μM and slightly higher at 70 μM (Table S1).

The secondary structure of the protein was studied using far-UV CD spectroscopy before and after fibril formation to verify the conversion of monomeric protein (mostly random coil structure) to amyloid fibrils (mostly β-sheet structure) (Figure c). As can be seen from Figure c, the fibrillar sample gave a CD spectrum that had a negative peak at 217 nm, a typical characteristic of a sample containing β-sheet structure, which indicates the conversion of the monomers to fibrils. Aliquots of the same samples were also supplemented with the fluorescent dye, ANS, and the fluorescence spectra were recorded. The fluorescence intensity of the fibrillar sample increased significantly compared to monomers (Figure a and 4b). This confirms the presence of amyloid fibrils in the fibrillar sample.

Figure 4

Verification of the presence of fibrils by optical spectroscopy. (A) ANS Fluorescence spectra before and after fibril formation of wild-type α-synuclein in water with the pH change measured using a pH electrode. (B) ANS Fluorescence spectra before and after fibril formation of the α-synuclein 5Q mutant in water with the pH change measured using a pH electrode. (C) CD-spectra before and after fibril formation of wild-type and 5Q mutant α-synuclein in water with pH changes measured using a pH electrode. (D) CD-spectra of wild-type α-synuclein before and after fibril formation in a weak buffer containing resazurin. (E) CD-spectra of α-synuclein in 100% D2O before and after fibril formation as monitored using NMR spectroscopy.

Change in pH Detected with pH Indicator

The pH change during fibril formation was also monitored using a complementary method that does not involve a pH electrode. α-Synuclein was allowed to form fibrils in a weak buffer system, constituting 0.5 mM MES and 12.5 μM resazurin, 0.02% NaN3 at pH 5.8. Resazurin is a pH indicator that shows a color change from purple to blue from pH 5.5 to 6.5. The starting pH of the sample was pH 5.8, giving the sample a purple color. The sample was then allowed to fibrillate, and the color of the fibril containing sample was compared to the monomeric sample (Figure b).

A change in the color from purple to blue was observed as fibrils were formed, indicating an increase in pH. For comparison, the pH of this sample was also measured to 6.4 with a pH electrode after fibril formation, in accordance with the color change. As before, CD-spectra were recorded before and after the fibril formation, to verify the conversion from monomeric to fibrillar sample (Figure d). Absorbance spectra were also recorded for comparison between the monomeric and fibrillar samples with indicator, but could not be used for quantification of the pH shift because of extensive baseline elevation for the fibrillar sample due to light scattering. Therefore, the indicator data serve as a qualitative confirmation of the upshift in pH.

Increase in pH Detected with NMR Spectroscopy

In addition to using a pH sensitive dye, solution NMR spectroscopy was used as a second complementary method to report on the change in pH during fibril formation. Here, the chemical shift of Hδ2 of the imidazole ring in His50 of the monomer was used to monitor the pH.[48−50] The pKa value of histidine within the monomeric α-synuclein in 20 mM sodium phosphate buffer, pH 7, had previously been measured to be 6.78.[30] At pKa ± 1, the buffer capacity of the imidazole side chain is within 33% of its maximum, and at pKa ± 1.5, it is within 10% of its maximum. The pH of monomeric and fibrillar samples in water (pH 5.6 to 6.5) falls within the buffer region of the histidine residue, making it possible to utilize it as a pH sensor. All the signals obtained in this experiment come from monomers, as the fibrils are too big to be observed by liquid state NMR (see in more detail in next paragraph), and therefore we can assume that any change in protonation of His50 is due to the change in pH of the solution and not due to pKa perturbations resulting from conformation changes. A sample containing 70 μM α-synuclein in 99.9% D2O was analyzed by solution state 1D 1H (proton) NMR before and after amyloid fibril formation. The peak corresponding to the Hδ2 of the imidazole ring in His50 shifted from 8.55 to 8.38 ppm (Figure c). CD-spectra were recorded before and after the fibril formation, to verify the formation of fibrils (Figure e).

2D Heteronuclear and Diffusion Ordered NMR Spectroscopy

2D 15N–1H HSQC NMR spectra showing resolved signals from amides in nearly all amino acid residues were run on both the monomeric and fibrillar sample (Figure a). The chemical shift differences between the two samples were very limited except for the region around His50.

Figure 5

(A) 2D 15N-1H HSQC NMR spectra of 15N labeled α-synuclein in 5% D2O/95% H2O. (B) Relative intensities from diffusion ordered (DOSY) NMR spectra. Monomer signals in the absence of fibril are shown in red (diamonds in B), and in the presence of fibrils in blue (triangles in B). In B, the signal intensity is calculated relative to the intensity at a gradient strength of 0.05 of max strength.

Diffusion ordered NMR spectroscopy (DOSY) was used to distinguish between signals originating from monomers (faster diffusing species) and fibrils (slower diffusing species). The highest gradient strength used was strong enough to suppress the signals in the purely monomeric sample (Figure b). The signals in the fibril containing sample were only slightly less suppressed, indicating that they originate from monomeric species in fast exchange with a slower diffusing species such as fibrils rather than from the fibril itself. We can say that chemical shift change in the His50 originates from the unstructured monomers in both the monomeric and the fibrillar sample and is presumably due to a change in pH, causing a change in the protonation state of the His50.

Continuous Measurement of pH during Fibril Formation

The pH of a solution containing freshly purified monomeric α-synuclein was measured continuously during formation of amyloid fibrils using a Probe Drum instrument equipped with a pH electrode. A small pH drop was detected at the start of the aggregation (first 3 h), followed by a gradual increase in pH over more than 40 h (Figure ). Other samples, supplemented by ThT, allowed simultaneous measurement in the same sample of pH and ThT fluorescence as well as light scattering. These data show that the onset of the pH increase coincided with the onset of ThT fluorescence or scattering (Figure S1).

Figure 6

Time evolution of pH followed with pH electrode using Probe Drum. 70 μM sample was stirred at 237 rpm at 37 °C.

Calculations of the Average pKa Upshift

As mentioned above, an increase in pH from 5.6 to 6.5 was observed for 20 μM α-synuclein in water. The experiments were performed in a closed system with a defined number of protons which makes it possible to calculate the change in the average apparent pKa value (ΔpKaave) from the change in pH during the fibril formation (ΔpH). In order to do so the degree of protonation, Q, was calculated for each individual titratable acidic group in α-synuclein (including the C-terminus), according toThis was done by assuming that the pKa values, and therefore the dissociation constants (Ka) of the acidic groups of the monomeric α-synuclein, were the same as published by Croke et al.[30] The average protonation state of the acidic groups Qave was calculated from the protonation state, Q, of individual groups.

By knowing the number of titratable acidic groups (# acidic groups) at this pH range, the Qave, and the concentration (C) of protein molecules, the number of protons bound to the acidic groups can be calculated according toThe number of protons free in solution (free[H+ ]mono) in the monomeric sample was then calculated from the pH value of the sample:The total number of protons in the monomeric sample (total[H+]mono) is equal to the sum of the bound and the free protons:The total proton concentration in the monomeric sample (total[H+]mono) and the fibrillar sample (total[H+]fib) was assumed to be the same, as the pH change occurs within the same sample, a closed system:From now on we will consider the fibrillar sample. The concentration of protons free in solution in the fibrillar sample (free[H+]fib) was calculated from the pH measured in the fibrillar sample:The total concentration of protons bound to the acidic groups (bound[H+]fibtot) within the fibrillar sample was then calculated according toWe have observed that, during the aggregation of α-synuclein, the monomers are not fully consumed into fibrils, in agreement with Gray et al. (2015).[51] Therefore, the remaining monomer concentration (Cfibmono) after fibril formation in water was measured (see SI, Figures S2 to S4)) and was taken into account in the calculations. The monomeric Qave was assumed to be the same in both samples. The monomeric Qave and the monomeric concentration in the fibrillar sample (Cfibmono) was used to calculate the number of protons bound to monomers in the fibrillar sample (bound[H+]fibmono).From that the number of protons bound to fibrils within the fibrillar sample was calculated (bound[H+]fibfib)To calculate the average Ka for the fibrillar sample, the concentration of acidic protein side-chain groups in the fibrils that were not protonated (free[P]fibfib) was calculated according toFrom that the proton equilibrium constant Ka for the fibrillar sample could be obtained:The Ka was then used to get the apparent pKaave value of the acidic groups within the fibrillar sample:For comparison, the apparent pKaave was also calculated for the 5Q mutant. Here, Qave was calculated using the pKa value of the remaining acidic residues, i.e. five groups less than for the wild-type protein. The concentration of free protons in solution was then calculated from the measured pH values of the 5Q mutant sample, before and after fibril formation.

The calculated apparent pKaave values of the fibrillar sample and the ΔpKaave are given in Table . The ΔpKaave for the wild-type sample was 1.1 while for the 5Q mutant, having 5 less acidic groups in the tail, it was 0.5. This suggests that this increase in pKa values during amyloid formation of α-synuclein is related to the high density of acidic residues in the C-terminal tail.

Table 2

Apparent pKaave Values of the Acidic Residues of α-Synuclein Calculated for Both the Monomers and Fibrilsa

	pKa monomers	pKa fibrils	ΔpKaave
Wild-type	4.3	5.4	1.1
5Q mutant	4.2	4.7	0.5

Based on this, the upshift in the apparent pKaave values, ΔpKaave, was calculated.

The contribution of autoionization of water was calculated and was found to have no pronounced effect on the calculated pKaave, with a 0.0004 units higher ΔpKaave value for the wild-type α-synuclein, i.e. beyond the precision of our measurements.

The number of protons taken up per monomer during fibril formation were calculated to be about 2 protons per monomer in nonbuffered closed system with a defined number of protons (our experimental system), which partially neutralized the protein.

Molecular Modeling of Protonation States

α-Synuclein Fibril

The ionization states of α-synuclein monomer and fibril (represented by ten planes, PDB: 2N0A) were studied using constant pH Metropolis Monte Carlo simulations; see Figure . Differences in residue partial charges in a monomeric and fibril state were calculated for both the wild-type α-synuclein and the α-synuclein 5Q mutant at constant pH 6.5 (Figure ) and constant pH 5.5 (Figure S5).

α-Synuclein Monomers

In addition to experiments, the ionization state of α-synuclein monomer can be estimated using Metropolis Monte Carlo Computer simulations. To validate the peptide interaction model described earlier, we first studied the stoichiometric pKa-values of aspartate and glutamate side chains in an isolated α-synuclein monomer. These have previously been measured using NMR,[30] and in Figure we show how the experiment and simulation correlate.

Figure 8

Correlation between measured[30] and simulated stoichiometric pKa-values of aspartates and glutamates in α-synuclein monomers. The Metropolis Monte Carlo simulations were performed at 15 mM salt and the pKa-value was defined as the pH where the degree of ionization is 0.5. Values marked in red indicate that the sites are found in the C-terminal tail. Only a single peptide chain is present in the simulation cell, corresponding to infinite protein dilution, while the ionic strength was set to 15 mM.

Discussion

Proton Uptake during Amyloid Fibril Formation

In this work, the effect of aggregation on the pKa values of α-synuclein was investigated and an increase in pKa was observed by measuring the change in pH during fibril formation of α-synuclein in water, using three complementary methods. A pH increase from 5.6 to 6.5 for 20 μM α-synuclein was quantified using a pH electrode. This implies a proton uptake, i.e. an increased affinity for protons, indicating an upshift in the apparent pKa values during fibril formation. Based on those results, the change in the apparent average pKa value (ΔpKaave) of the titratable acidic residues in α-synuclein was calculated to be ΔpKaave = 1.1, i.e. an increase in average pKa values from 4.3 to 5.4.

The proton uptake during fibril formation in a nonbuffered closed system was found to be about 2 protons per monomer, in good agreement with the simulations performed at pH 6.5, which also showed an uptake of about 2 protons per monomer during fibril formation.

In a closed system, the protons taken up by the fibril must come from somewhere. This can only in part be explained by the reduction in proton concentration in the water, and indeed the protein itself serves as an internal buffer. Protons released by the monomeric fraction are taken up by the fibrillar fraction. As fibril formation proceeds, the monomers that remain in solution (having lower apparent pKaave values compared to fibrils) release protons, which are in turn taken up by the monomers that become a part of a fibril (having higher apparent pKaave values of the acidic residues). Therefore, the charge state of α-synuclein in the monomeric state and fibrillar state differ in the fibrillar sample. After fibril formation, the net charge of free monomers has increased while the net charge of monomers in the fibrils have decreased. This lowers the electrostatic repulsion within the fibril.

It should be noted that the monomeric pKa values used in all the calculations are obtained from Croke et al. for 250 μM α-synuclein in buffer.[30] However, buffer and protein concentrations affect the ion binding affinities of charged proteins due to electrostatic salt- and self-screening.[52,53] The protein concentration in our samples, 20 μM, was more than 10 times lower than in the case of Croke et al., and the samples contained no other buffer than the protein itself. Also, from our studies, we see that the pH value of pure monomer in water is 5.6, indicating that the monomeric pKa values might be closer to this value than to 4.3.[30] In order to account for this uncertainty, the pKa upshift was thus recalculated using several values for the pKa values of monomer in the absence of fibrils. The result showed that the average upshift in pKa value would not be much different if the pKa values of the monomers would be higher (see SI, Figure S6).

Result of High Density of Acidic Residues within the Unstructured C-Terminal Tail?

The pH change during fibril formation was measured in the same way for an α-synuclein 5Q mutant containing five fewer acidic residues within the C-terminal tail, making the tail less acidic and less charged. This was done in order to gain better insight into the effect of the high charge density in the C-terminal tail, to investigate if the upshift in pKaave could be explained by the close proximity of acidic residues of the unstructured C-terminal tails extending from the fibril core. The pH increase during fibrillation was significantly lower for the 5Q mutant, giving an average pH increase from 6.4 to 6.6 for a sample of 20 μM α-synuclein. The pKaave was calculated to increase from 4.2 to 4.7, or a ΔpKaave of 0.5 was found. These results demonstrate that decreased acidity and charge density of the C-terminal tail results in a smaller pH increase during fibril formation. This implies that the pH increases during fibril formation (proton uptake), and therefore, the upshift of pKa values during amyloid formation of α-synuclein is linked to the high density of the acidic residues in the C-terminal tail. Those results were in good agreement with the results obtained from Monte Carlo simulations, also showing significantly less effect on protonation states during fibril formation of the 5Q mutant.

Electrostatically induced pKa Shifts

Our results are line with a previous study by Urry et al. (1993) on upshifted pKa values. They demonstrated a correlation between the decrease in density of acidic residues and the pKa decrease, and suggested that the pKa upshifts are primarily electrostatically induced and not induced due to desolvation.[11] Upshifts in pKa values due to the high density of acidic residues on protein surfaces of halophilic proteins has been reported.[54] Likewise, Kesvatera et al. found an asymmetry in pKa values in the EF-hand protein S100G with highly upshifted pKa values in the highly acidic loop regions and unperturbed or down-shifted values at the opposite end of the protein.[22] These observations are in accordance with our system where the extent of the upshifts was related to the charge density within the tails. After submission of this work, a study appeared with chemical shift changes in the C-terminal tail of α-synuclein upon fibril formation measured by solid-state NMR indicating a similar upshift in pKa values.[55]

Upshifts in pKa values due to changes in charge–charge interactions can be understood through Coulomb law, by calculating the Coulomb force between two adjacent charges:where F stands for the electrostatic force or the Coulomb force between two charges (q1, q2) separated by the distance r. When C-terminal tails come together at the fibrillar surface, the distance between the charges of different tails decreases, increasing the Coulomb force between tails. The increase in electrostatic force between charges is counteracted by upshift in pKa values, making the tails proportionally more neutral (lower q1 and q2) and therefore decreasing the net force between adjacent tails.

Molecular Simulations

Figure shows that the simulated results for monomeric α-synuclein are in fair agreement with NMR experiments,[30] with similar upshifts and downshifts for aspartic and glutamic acids; the encircled red points correspond to residues in the C-terminal region. While deviations are expected due to the simplified interaction scheme, the data show that general trends in the observed pKa shifts can be explained by an approximate electric potential due to the surrounding chemical environment.

Simulation data are in agreement with experimental data, and both show that protons are taken up during fibril formation of α-synuclein. Figure a shows a simulation snapshot of an α-synuclein fibril where the rigid middle is surrounded by flexible N- and C-termini. Significant differences are observed between the protonation states of carboxyl groups in the monomeric state and the fibril (Figure b). In the α-synuclein monomer, the pH titration curves are approximately symmetric while highly stretched in the fibril where full protonation is never reached even at the lowest pH. Fibril formation thus leads to a significant change in proton binding capacitance. This is because the densely packed fibril is cationic at low pH (see Figure S7), and the carboxyl groups are thus exposed to a strong positive potential, opposing protonation. The opposite is seen at high pH where the electrostatic environment is negative and protonation is favored compared to ideal titration. Figure b highlights the difficulty in describing electric perturbations by a simple upshift or downshift of the thermodynamic pKa-value. This is formally seen in the definition of the stoichiometric or apparent acid dissociation constant, Ka* = KaΓ(pH), where Γ is an activity coefficient product reflecting the local chemical environment of the site. Since Γ in nontrivial ways varies with pH, conformation, and other solution conditions, so does Ka*. Thus, the pKa “shift” is in reality a function of pH whereby the resulting titration curve can deviate far from a merely translated Henderson–Hasselbalch behavior. Still, from pH 6 and up, the titratable sites in the C-terminal region are on average predicted to be more protonated in the fibrillar state than in the monomeric state which of course agrees with the qualitative notion of a pKa upshift during fibrillation.

It is instructive to trace the changes in the average protonation state of individual residues when embedding a monomer into a fibril. This is shown in Figure c where the FASTA letter background color shows the charge increment (blue) or decrement (red). At pH 6.5 the largest differences are observed in the C-terminal with all carboxyl groups being less negatively charged than in the monomer. This is due to the strong negative electric potential exerted by the surrounding C-terminal ends, causing proton uptake. In agreement with our experimental observations, the charge shifts are strongly subdued in the 5Q mutant (Figure d) and only marginal charge shifts between α-synuclein monomers and fibrils are observed. It is important to note that the simulations are performed at constant pH, i.e. fully buffered, while in the experiments proton uptake upon fibril formation is monitored through a pH change.

It has been argued that an elimination of transient interactions[32,56] between the N- and C-terminus of wild-type monomers upon formation of fibrils could be an alternative explanation to the pKa upshift. However, the results presented in Figure c and 7d show asymmetric perturbations with very small changes in the N-terminus. It should be emphasized that all residues are included in the simulations and we do not add restraints to the monomer conformation; it is thus allowed to take any transient conformation governed by the sequence. We see a clear change in the charge states in the C-terminus when monomer forms fibrils, and that the effect is strongly subdued for the 5Q mutant (see Figure c and d). The combination of our experimental and theoretical work strongly supports that the close stacking of the C-termini has a major role in the pKa upshift during fibril formation.

Driving Force of Fibril Formation

The thermodynamic stability of proteins is dependent on the ionization state of the different amino acid residues. As seen from Coulomb’s law, the pKa values of the residues within α-synuclein, and therefore the charge density of the C-terminal tail, can affect the thermodynamic stability of the fibrils. A decrease in electrostatic repulsion between adjacent tails, due to the pK upshifts, results in the fibril formation being thermodynamically more favorable.

According to literature, the hydrophobic effect and hydrogen bonding are two main intermolecular interactions stabilizing amyloid fibrils.[37,57−60] Similar to formation of globular protein structures, the fibril core is stabilized through the hydrophobic effect and hydrogen bonding within the β-sheets of the amyloid structure,[57] although hydrogen bonding may be more or less iso-energetic over monomer and fibril as long as all potential hydrogen bonding groups are engaged in hydrogen bonds with water molecules or other protein groups. This is reflected in the known amyloid structures, where the fibrillar core has been found to profoundly constitute hydrophobic and polar, noncharged residues (Figure ).[31,61,62] Highly charged C-terminal tails and extensive electrostatic repulsion between adjacent monomers will presumably oppose the hydrophobic effect, making fibril formation less favorable. Without upshifts in the pKa values and a higher degree of neutralization of the tails on the fibrillar surface, the favorable interactions between adjacent monomers would presumably be overridden by the unfavorable electrostatic repulsion.

Increase in pKa Values and the Preventive Role of the Tail

Previous studies of α-synuclein have shown that truncation of the tail,[39,40,63,64] charge-shielding,[40] and neutralization[36,37] affects the aggregation propensity. It has been suggested that the tail has a preventive role against aggregation.[39] Protective characteristics of negative charges has been inferred for other systems; for example, the high number of negatively charged residues on the surface of halophilic proteins might prevent aggregation through electrostatic repulsion.[54]

It is interesting to speculate on the implications of the higher degree of neutralization of the tails within the fibrils, which makes the surface less charged at mildly acidic pH. The pKa upshift and consequently the higher degree of neutralization presumably modulates the interaction between the fibrillar surface and the approaching monomers with implications for surface-catalyzed secondary nucleation.

Comparison to Other Systems

Perturbations of pKa values are commonly found in enzyme active sites and within protein cores. Quantifications of pKa upshifts are frequently reported in the literature. However, to our knowledge, the upshift in pKa values during amyloid fibril formation has not been reported. Proton uptake or release during fibrillation has been observed for lysozyme and glucagon in a previous study.[65] This leads to neutralization of the protein net charge and indicates an upshift in pKa values within the fibrillar state in these cases.

In addition, pKa perturbations have been observed during other self-assembly reactions. The pKa value of a C-terminal carboxyl acid group was found to upshift during fiber and hydrogel formation of Fmoc-diphenylalanine.[66,67] Additionally, upshifts in pKa values of fatty acids were found to be a result of a decrease in intermolecular distance upon membrane self-assembly and dependent on the chain length, which modulates the intermolecular interactions and thereby the distance between the acidic headgroups at the interface.[68] These additional examples demonstrate that pKa value perturbations are likely to occur within all self-assembly processes of charged substances with ionizable groups to counteract the increase in local concentration of negative charges at surfaces upon a decrease in distance between molecules.

Conclusion

A significant upshift in pKa values of at least 1.1 units on average accompanies amyloid fibril formation of α-synuclein. We demonstrated, measured, and calculated the upshift using three complementary experimental methods as well performing Metropolis Monte Carlo simulations on a comparable system. A significant change in proton binding capacitance upon fibril formation is due to the accumulation of titrating groups on the fibril surface. The significantly smaller upshift for an α-synuclein 5Q mutant having five less charged residues in the C-terminal tail implies that the magnitude of the upshift in pKa is related to the density of acidic residues in the C-terminal tail. The observed pKa perturbations during fibril formation provide further insights into the importance and the role of electrostatic interactions in the fibril formation of α-synuclein. The results have implications for other amyloid forming proteins and self-assembling systems in general, for which the methods reported here could provide quantitative information on pKa perturbations and the role of electrostatic interaction in self-assembly processes.