Tautomeric Effect of Histidine on the Monomeric Structure of Amyloid β

Oct 17, 2016 - Investment flows at J.P. Morgan Healthcare Conference. The J.P. Morgan Healthcare Conference in early January has traditionally been th...
0 downloads 5 Views 925KB Size
Subscriber access provided by UNIV NEW ORLEANS

Article

Tautomeric Effect of Histidine on the Monomeric Structure of Amyloid #-Peptide(1-40) Hu Shi, Baotao Kang, and Jin Yong Lee J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b08685 • Publication Date (Web): 17 Oct 2016 Downloaded from http://pubs.acs.org on October 20, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry B is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Tautomeric Effect of Histidine on the Monomeric Structure of Amyloid β-Peptide(1-40) Hu Shi†, Baotao Kang‡*, and Jin Yong Lee†* †

Department of Chemistry, Sungkyunkwan University, Suwon 440-746, Korea



Department of Chemistry and Chemical Engineering, University of Jinan, Jinan 250022, P. R.

China

ABSTRACT: Histidine state (deprotonated, neutral, and protonated) is considered an important factor influencing the structural properties and aggregation mechanisms in amyloid β-peptides which are associated with the pathogenesis of Alzheimer’s disease. Understanding the structural properties and aggregation mechanisms is a great challenge because two forms (the Nε-H or NδH tautomer) can exist in the free neutral state of histidine. Here, replica exchange molecular dynamics simulation was performed to elucidate the changes of structure and mechanism of aggregation influenced by tautomeric behaviors of histidine in Aβ (1-40). Our results show that sheet dominating conformations can be found in His6(δ)-His13(δ)-His14(δ) (δδδ) isomer with significant antiparallel sheet structure between R5-D7 and L34-G38, as well as between L17-F20 and L34-G38, implying new aggregation mechanism may exist to promote oligomers and/or

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 21

aggregates generation. This work is helpful for understanding the fundamental tautomeric behaviors of neutral histidine in the process of aggregation.

Introduction Alzheimer’s disease (AD), in which protein folding abnormalities is one of the pathogenic causes, is the most common form of dementia affecting more than 35 million people.1 It is pathologically characterized by progressive intracerebral accumulation of amyloid β-peptides (Aβ) and tau protein.2,

3

The major component of plaques, Aβ fragments, are derived from

amyloid precursor protein (APP) by β- and γ-secretase proteolysis.4, 5 Approximately 90 % of the Aβ fragments generated are the 40 residue fragment (Aβ(1-40)).6,

7

Aggregation of the

fragments, including toxic oligomers, plays a crucial role in the pathogenesis of AD.8-11 During Aβ aggregation, the “nucleated growth mechanism” has been accepted to describe the process of development from unfolded monomers to mature fibrils. That is, soluble and natively unfolded monomers first form β-structures, then generate the initial nucleus and proto-fibrils, and finally become mature fibrils.12,

13

Many factors,14-17 such as pH environment, APP mutation, β-/γ-

secretase’s activity, and metal ions, can destroy the frail balance of Aβ production and clearance,18-21 and induce aggregation. Many literatures have focused on the effect of pH and researches have argued that small deviations in physiological pH can sensitively influence the chemical properties and aggregation tendency of Aβ peptides,22 which may be the reason why diseased brains are more acidic than the normal brain.23 The aggregation tendency is significantly high at pH 5 and 6.24 At pH 4 to 7, the β-sheet structure of Aβ(1-42), Aβ(1-28) and Aβ(1-39) is formed in aqueous

ACS Paragon Plus Environment

2

Page 3 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

trifluoroethanol.25 At pH 6, Aβ(1-28) and Aβ(10-42) folding are favorable for aggregation and are driven by hydrophobicity.26 Moreover, aggregation of Aβ(1-42) and shorter peptides is enhanced at pH 3.5 to 6.5.17 After reviewing these literatures, we determined that histidine is the only titratable amino acid with an uncertain condition.27 When in a protonated state, three protonated histidine residues promote the formation of β-sheet in Aβ(1-42) via a reduction in electrostatic repulsion between the two terminal regions,28 implying that the protonation states of histidine influence on peptide properties. However, regarding the unprotonated state, there is still a lack of studies on the effect of tautomeric state in neutral histidine residues in the full-length Aβ(1-40), although in N-AcAβ(13-23)NH2 monomer, only flexible folded structure was obtained in both unprotonated and protonated state of histidine.29. The imidazole side chain exists in two forms, the Nε-H tautomer (NεH, denoted as ε) and the Nδ-H tautomer (NδH, denoted as δ) (Scheme 1), and the population ratio of the two forms (ε/δ) is around 1:0.16.30-32 Tautomeric equilibrium occurs in picosecond.33 However, the population ratio of tautomeric isomers may significantly depend on environmental factors,34 such as histidine position (sequence number) and sidechain-sidechain interactions. Even a slight change in the surrounding physiological environment can perturb the population ratio of tautomeric isomers. Thus, a clear understanding of the tautomeric behavior of histidine residues in Aβ fragments is limited in experiment. However, tautomeric behaviors play an important role in the structural properties and aggregation mechanisms. In this context, we attempted to uncover how tautomeric behaviors effect the aggregation of Aβ by replica exchange molecular dynamics (REMD) simulations. Compared with standard molecular dynamics simulations, REMD is an enhanced sampling method that allows systems of similar potential energies to sample conformations at different temperatures.35

ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 21

Scheme 1. Tautomeric forms of histidine

The Aβ(1-40) monomer has three histidine residues located at position 6, 13, and 14 of the amino acid sequence. Each histidine has two tautomeric isomers (ε and δ), so there are eight isomers in the Aβ(1-40) monomer. The eight tautomeric isomers are denoted as (εεδ) in cases when the histidine residues are ε, ε, and δ isomers at amino acid positions 6, 13, and 14, respectively. Likewise, other isomers are denoted as (εεε), (εδε), (εδδ), (δεε), (δεδ), (δδε), and (δδδ) according to tautomeric forms. These isomers were investigated to obtain insight into the structural and aggregation features of Aβ(1-40). In order to determine if there is an influence of different force fields, compared with the ff99SB non-polarization force field in our REMD simulations, we employed the ff02 polarization force field for further standard MD simulation in (δδδ) and (εεε) isomers. Highly consistent ensembles (Figure S2) indicate that there is little dependence on the choice of force field. Furthermore, we added a 300 ns (δδε) simulation using (δδδ) secondary structure assignment as the initial structure. Spontaneous structural changes also supported our claim that configuration of the histidine residues is critical to the structural features (Figure S2).

ACS Paragon Plus Environment

4

Page 5 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Simulation Details The crystal structure of the full-length Aβ(1-40) (PDB entry 1BA4), taken from the Protein Data Bank (http://www.rcsb.org/), was used in the present study. The eight tautomeric isomers were prepared by residual substitution with the Nε-H or Nδ-H tautomer at position number of 6, 13, and 14 of the amino acid sequence. The ff99SB non-polarization force field36 and implicit solvent generalized Born model were employed to perform the REMD simulations with the Amber 12/Amber14 package37. The AMBER ff99SB force field has been frequently used to simulate the proteins misfolding,38-40 which presents a careful reparametrization of the backbone torsion terms in ff99 and achieves a much better balance of four basic secondary structure elements.41 The SHAKE algorithm was applied to constrain the bonds including hydrogen. Particle mesh Ewald was used to describe the long-range electrostatics. The cutoff distance has an infinite value for non-bonded interactions. The temperature was controlled by the Langevin thermostat with coupling constants of 1.0 ps. The time step is 2.0 fs. Each simulation consisted of 10 replicas with the temperature range from 310 to 510 K. Exchanges between replicas were attempted every 4 ps. In the current simulations, more than 700 ns of data (each replica) were collected for each isomer. Our REMD results show that the acceptance ratio is ~22 %. Only the converged trajectories at 310 K were used to analyze the data. In addition, structural changes were exhibited using Visual Molecular Dynamics (VMD).42 We calculated the Pearson correlation coefficient (PCC) and Root-mean-square deviation (RMSD) (Table 1 in the main text) between our simulations and experimental data obtained from references (Ref. 55-57 in the main text). Our calculated 3JHNHA-coupling constants were derived from each converged trajectory using the Gromacs 5.043 analysis module. Coefficient values for the Karplus equation were 6.51, -1.76, and 1.60.44

ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

Results and Discussion

0.5 helix sheet

Probability

0.4 0.3 0.2 0.1 0.0 εεε

εεδ

εδε

εδδ

δεε

δεδ

δδε

δδδ

1.0 εεε εεδ εδε εδδ δεε δεδ δδε δδδ

0.8

Para

0.6 0.4 0.2 0.0 0

10

20

30

40

30

40

30

40

Residue 1.0 0.8

Anti

0.6 0.4 0.2 0.0 0

10

20

Residue 1.0 0.8

α-Helix

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 21

0.6 0.4 0.2 0.0 0

10

20

Residue

Figure 1. Probability distributions of averaged sheet and α-helix content in each isomer (A), and residual probability distributions of parallel β-sheet (B), antiparallel β-sheet (C), and α-Helix (D).

ACS Paragon Plus Environment

6

Page 7 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

After careful evaluation of convergence (Figure S1 and Table S1), the probability distributions of helix and sheet formation (Figure 1A) are remarkably different for eight isomers (details in Table S2), indicating that diverse neutral histidine tautomeric states (ε or δ) at sequence 6, 13, and 14 can dramatically influence the structural properties of Aβ(1-40). In particular, as the crucial secondary structural characteristics, the β-sheet content of (δδδ) (22.4 %) is much higher than that of other tautomeric states. As opposed to N-AcAβ(13-23)NH2 monomer that assumes a flexible folded structure,29 the (εεε) isomer exhibits ~49 % parallel β-bridge content in H6 and H14 (Figure 1B). In the (εδε) isomer, parallel β-bridge can be formed between Q15 and G38 with a probability of ~41 % (Figure 1B), while antiparallel β-bridge was generated between F20 and D23 with content of ~24 % (Figure 1C). Although the (εεε) and (εδε) isomers show significant chances to form β-bridge structures, these structures cannot further lead to formation of β-sheet. In contrast, antiparallel β-sheet can be found in the (δδδ) isomer, arising mostly between R5-D7 (content: 87 ~ 88 %) and L34-G38 (content: 37 ~ 90 %), and between L17-F20 (content: 39 ~ 88 %) and L34-G38 (Figure 1C). In the (δδδ) isomer, meanwhile, only E22-G33 exhibit marginal chances (< 29 %) to be involved in α-helical structures (Figure 1D). The observation of the dominant antiparallel β-sheet structure implies that the formation of the (δδδ) isomer may play a crucial role in Aβ oligomerization and/or aggregation because of the fatal structural transition from an individual soluble structure to a regular sheet structure.12, 45, 46

ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 21

Figure 2. Most abundant five conformational states and their probability distributions obtained by hierarchical agglomerative approach with cutoff of 5 Å via best-fit coordinate Root-mean-

ACS Paragon Plus Environment

8

Page 9 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

square deviation (RMSD) using alpha C atoms. Blue band represents L17-A21 region; green band represents sheet structure; blue sphere represents the N-terminus.

Based on cluster analysis, the (εεε) (δεδ), and (δδε) isomers were found to be unfolded or have a random propensity with significant structural diversity, which appear to be in a non-toxic state (Figure 2). Although a dominant structure was observed for the (εδε) isomer, the structure is rich in turn and coil structures rather than sheet or helix structures. Besides, owing to low sheet contents and high helix contents, the (εεδ), (εδδ), and (δεε) isomers also may not be toxic. It has been proposed that the toxicity of Aβ results from conformational transitions from α-helix or random coil to β-sheet structures.47 In contrast, the (δδδ) isomer is a β-sheet dominating conformer with 88.0 % of population based on cluster analysis (Figure 2). It is worth mentioning that the β-sheet structures occur in the N-terminus, the central L17-A21 region (recognition unit), and the C-terminal region. In particular, the β-sheet between the L17-A21 region and C-terminal region was reported as the intrinsic property in aggregation deposits. In this regard, the (δδδ) isomer could promote oligomers formation, and such oligomers may function as toxicity agents in the pathology of AD.11, 48, 49

ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 21

Figure 3. Most dominating conformer of (δδδ) isomer (A) analyzed by cluster analysis with same approach (describe in Figure 2) but using all atoms. Insert figure (B) shows hydrogen bond networks between strand regions. Contact map (C) of (δδδ) isomer based on Alpha C interaction (upper triangle) and residual sidechain interaction (lower triangle) analysis. In contact map, X and Y axis represent residue number index, color (Z axis) represents distance between X residual

ACS Paragon Plus Environment

10

Page 11 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

index and Y residual index. Region (a) and (b) show antiparallel β-sheet structures. Region (c) figure out sidechain interactions between R5, E11, H14, and Q15.

In order to understand the structural details of (δδδ) isomer, cluster analysis was employed using the same approach (above) but based on all atoms. The probability of the most abundant five conformational states were 54.1%, 11.7%, 8.1%, 6.4%, and 2.2%, respectively (Figure S3). The most dominant conformer is shown in Figure 3A. There is a tendency to adapt the antiparallel β-sheet between L17-V18 and V36-G38 (Figure 3C, upper triangle, (a)), which is in good agreement with the previous findings that the β-sheet primarily occurs at the CHC and Cterminus.50, 51 Meanwhile, another substantial antiparallel β-sheet character occurs between R5D7 and V36-G38 (Figure 3C, upper triangle, (b)). The same tendency was found between current study and previously reported studies that β-content can occur in N-terminus, especially, in the first six residues.52, 53 Ball et al. also reported the formation of two β-sheets in three β-strand structures from MD simulations: one β-sheet is between CHC and G9-H13, and the other one between CHC and residues M35-G37.54 To compare the present study with this previous study, (εεε) isomer was selected because ε isomer is the default histidine state in simulations. After analyzing 2384 frames of the (εεε) isomer (Figure S4), anti-parallel sheet was obtained in residues G9-H13, CHC, and residue S26-V39, which is consistent with the above findings. Additionally, in the (δδδ) isomer, exposure of the recognition unit (L17-A21 region) and charged residues implies that the (δδδ) isomer has ideal conditions as a potential seed to promote oligomers and/or aggregates generation. Moreover, the most commonly populated long-distance backbone hydrogen bonds between the CHC and the C-terminal regions reveal that L17/G37, V18/G18, and F19/M35 (Figure 3C, upper triangle, (a)) play an important role in the stability of

ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 21

the antiparallel β-sheet structure. Meanwhile, F4/V39, R5/G38, H6/G37, D7/V36, and S8/M35 (Figure 3C upper triangle, (b)) are crucial for the interaction between the N-terminal and Cterminal regions. In sidechain hydrogen bond networks (Figure 3B, and Figure 3C lower triangle, (c)), H14/E11, R5/E11, R5/V40, Q15/R5, and Q15/V40 hydrogen bonds were obtained. Compared with His6(δ) and His13(δ), we can see that His14(δ) plays an important role in the structural stability of (δδδ) via a hydrogen bond in H14/E11. In order to uncover the role of δ in (δδδ), (δδε), (δεδ), and (εδδ), the hydrogen bond population was collected in Table S3. In (δδδ), 44.90 % of the hydrogen bond population was obtained in E11/H14, and there was no interaction with H6 and H13. Compared with (δδδ), (δδε) lost the E11/H14 interaction, while interactions in Y10/H14, H14/L17, and H14/K16 were found due to substitution of ε for δ in residue 14. Interestingly, the interactions of L17 and K16 with H13 were detected in (δεδ) and (εδδ) with population range from 14.05 % to 16.10 %. Such backbone hydrogen bonds would obstruct the interaction between the central L17-A21 region and C-terminal region. This explains why (δδδ) shows antiparallel β-sheet feature while (δδε), (δεδ), and (εδδ) do not.

Table 1. Pearson correlation coefficient (PCC) and RMSD between our REMD calculated (Table S4) and experimental 3JHNHA-coupling constants.

εεε

with Ref. 55

with Ref. 56

with Ref. 57

PCC

RMSD(Hz)

PCC

RMSD(Hz)

PCC

RMSD(Hz)

0.44

0.97

0.44

0.97

0.54

0.72

ACS Paragon Plus Environment

12

Page 13 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

εεδ

-0.06

1.41

-0.09

1.21

0.40

0.82

εδε

0.16

1.38

-0.02

1.47

0.78

0.71

εδδ

0.01

1.24

0.11

1.08

0.62

0.70

δεε

0.13

1.59

0.17

1.56

0.58

0.67

δεδ

0.40

0.84

0.44

0.79

0.53

0.75

δδε

0.50

0.77

0.65

0.66

0.39

0.89

δδδ

0.30

1.07

0.28

0.97

0.48

0.80

PCC is derived from averaged 3JHNHA-coupling constant between simulated and experimental data. RMSD is standard deviation by calculating averaged deviation between simulated and experimental data.

Although it is too difficult to distinguish specific isomers in experiment, the ensemble of specific tautomeric state or mixed states should be in agreement with previous experimental data. In this regard, we calculated the PCC and RMSD (Table 1) between our simulations and experimental data obtained from references55-57. Compared with the experimental data in Ref. 55, maximum PCC of 0.50 with RMSD of 0.77 was located at the (δδε) isomer. When compared with Ref. 56, PCC changed into 0.65 with RMSD of 0.66, implying these experimental data are in agreement with the (δδε) isomer. Compared with Ref. 57, the PCC is 0.78 with RMSD of 0.71 for the Aβ(25-35) fragment, so we would conclude that the (εδε) isomer was preferred in Ref. 57. Significant combinations of the tautomeric behaviors of histidine could contribute to the

ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 21

diversity in experimental results which may be derived from different experimental conditions and environments. Up to now, the aggregation mechanisms influenced by the tautomeric state of histidine in Aβ(1-40) can be speculated out. The populations of Aβ(1-40) and Aβ(1-42) in human central nerve system are almost equivalent in both AD patients and normal humans. However, the rate of clearance is significantly reduced for AD patients.58 Considering that the tautomeric effect is a fundamental behavior in the free neutral histidine state, the (δδδ) isomer with rich sheet character would contribute to the pathway that produces Aβ oligomers and/or aggregates. In the nonpathogenic pathway, the production and clearance of Aβ pieces (monomers, oligomers, or aggregates) have a perfect balance for maintaining a normal healthy concentration of Aβ in the brain.59 However, in the pathogenic pathway, a weakened clearance (Apolipoprotein E4 inheritance, faulty Aβ degradation, etc60) destroys the frail balance, and induces the histopathological properties of AD. The present study proposes that the frail balance can be shifted by developing small molecules to inhibit the (δδδ) isomer.

Conclusion In summary, REMD simulations were carried out to investigate the structural properties and aggregation mechanisms of Aβ(1-40) under different tautomeric behaviors of histidine in the free neutral state. We found that different isomers display distinctly different compositions of secondary structures. Among eight isomers, the (εεε), (δεδ), and (δδε) isomers are structurally disordered in the absence of any dominant conformational features. The (εδε) isomer does have a dominant conformation that however is composed mainly of turn and coil structures. The (εεδ),

ACS Paragon Plus Environment

14

Page 15 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(εδδ), and (δεε) isomers exhibit also dominant conformations with regular secondary structures but these structures are largely helical. Only the (δδδ) isomer could assume dominant conformations enriched with β-sheet structures. Secondary structure analysis, cluster analysis, and contact map suggested that these conformations are featured with antiparallel β-sheet between R5-D7 (content: 87 ~ 88 %) and L34-G38 (content: 37 ~ 90 %), as well as those between L17-F20 (content: 39 ~ 88 %) and L34-G38. Furthermore, we found that H14 plays an more important role in structural stability in the (δδδ) isomer as compared with H6 and H13. The interactions of each histidine with L17 and K16 play a critical role to different structural features of (δδε), (δεδ), and (εδδ) in comparison with (δδδ), since backbone hydrogen bonds would obstruct interactions between central L17-A21 region and C-terminal region. Together, our study argues for the tautomeric states of neutral histidine can influence the aggregation of Aβ(1-40), which would be useful for understanding the pathology of AD and design of drugs for therapeutic treatment. Our study opens a new avenue to the study of the effects of tautomerism of histidine in protein misfolding diseases.

ASSOCIATED CONTENT Supporting Information. Convergence details, secondary structure distributions in each isomer, structure properties of (δδδ) isomer in ff02 polarization force field, Hydrogen bond interaction, 3

JHNHA-coupling constant and standard deviation in each isomer.

AUTHOR INFORMATION Corresponding Author

ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 21

* Jin Yong Lee, Email: [email protected], Tel: +82-031-299-4560 * Baotao Kang, Email: [email protected], Tel: +86-0531-827605961 Notes The authors declare no competing financial interests.

ACKNOWLEDGMENT This work was supported by National Research Foundation (NRF) grants funded by the Ministry of Education, Science and Technology, subjected to the project EDISON (Education-research Integration through Simulation on the Net, Grant No.: 2012M3C1A6035359). KISTI supercomputing center through the strategic support program for the supercomputing application research [No. KSC-2014-C2-056].

ABBREVIATIONS AD, Alzheimer’s disease; Aβ peptides, amyloid β-peptides; REMD, replica exchange molecular dynamics; RMSD, Root-mean-square deviation; PCC, Pearson correlation coefficient.

REFERENCES 1. Goedert, M.; Spillantini, M. G. A century of Alzheimer's disease. Science 2006, 314, 777-781. 2. Hardy, J.; Selkoe, D. J. The amyloid hypothesis of Alzheimer's disease: progress and problems on the road to therapeutics. Science 2002, 297, 353-356.

ACS Paragon Plus Environment

16

Page 17 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3. Alonso, A.; Zaidi, T.; Novak, M.; Grundke-Iqbal, I.; Iqbal, K. Hyperphosphorylation induces self-assembly of tau into tangles of paired helical filaments/straight filaments. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 6923-6928. 4. Masters, C. L.; Simms, G.; Weinman, N. A.; Multhaup, G.; McDonald, B. L.; Beyreuther, K. Amyloid plaque core protein in Alzheimer disease and down syndrome. Proc. Natl. Acad. Sci. U.S.A. 1985, 82, 4245-4249. 5. Kang, J.; Lemaire, H. G.; Unterbeck, A.; Salbaum, J. M.; Masters, C. L.; Grzeschik, K. H.; Multhaup, G.; Beyreuther, K.; Muller-Hill, B. The precursor of Alzheimer's disease amyloid A4 protein resembles a cell-surface receptor. Nature 1987, 325, 733-736. 6. LaFerla, F. M.; Green, K. N.; Oddo, S. Intracellular amyloid-[beta] in Alzheimer's disease. Nat. Rev. Neurosci. 2007, 8, 499-509. 7. Lovell, M. A. A potential role for alterations of zinc and zinc transport proteins in the progression of Alzheimer’s disease. J. Alzheimer's Dis. 2009, 16, 471-483. 8. Jakob-Roetne, R.; Jacobsen, H. Alzheimer's disease: from pathology to therapeutic approaches. Angew. Chem. Int. Ed. Engl. 2009, 48, 3030-3059. 9. Hamley, I. W. The amyloid beta peptide: a chemist's perspective. Role in Alzheimer's and fibrillization. Chem. Rev. 2012, 112, 5147-5192. 10. Kepp, K. P. Bioinorganic chemistry of Alzheimer's disease. Chem. Rev. 2012, 112, 51935239. 11. Lambert, M. P.; Barlow, A. K.; Chromy, B. A.; Edwards, C.; Freed, R.; Liosatos, M.; Morgan, T. E.; Rozovsky, I.; Trommer, B.; Viola, K. L.; et al. Diffusible, nonfibrillar ligands derived from Aβ1–42 are potent central nervous system neurotoxins. Proc. Nat. Acad. Sci. U.S.A. 1998, 95, 6448-6453. 12. Wilson, M. R.; Yerbury, J. J.; Poon, S. Potential roles of abundant extracellular chaperones in the control of amyloid formation and toxicity. Mol. Biosyst. 2008, 4, 42-52. 13. Savelieff, M. G.; Lee, S.; Liu, Y.; Lim, M. H. Untangling amyloid-beta, tau, and metals in Alzheimer's disease. ACS Chem. Biol. 2013, 8, 856-865. 14. Zhang, C. E.; Wei, W.; Liu, Y. H.; Peng, J. H.; Tian, Q.; Liu, G. P.; Zhang, Y.; Wang, J. Z. Hyperhomocysteinemia increases beta-amyloid by enhancing expression of gamma-secretase and phosphorylation of amyloid precursor protein in rat brain. Am. J. Pathol. 2009, 174, 14811491. 15. Shi, H.; Kang, B.; Lee, J. Y. Zn2+ effect on structure and residual hydrophobicity of amyloid β-peptide monomers. J. Phys. Chem. B 2014, 118, 10355-10361. 16. Zhang, R.; Ai, H.; Zhu, X.; Li, Q. Molecular simulations of human and mouse Aβ1–16 at different pH values: structural characteristics toward understanding Cu2+-coordinated amyloid beta spheres. ChemPhysChem 2016, 17, 1656-1668. 17. Burdick, D.; Soreghan, B.; Kwon, M.; Kosmoski, J.; Knauer, M.; Henschen, A.; Yates, J.; Cotman, C.; Glabe, C. Assembly and aggregation properties of synthetic Alzheimer's A4/beta amyloid peptide analogs. J. Biol. Chem. 1992, 267, 546-554. 18. Haass, C.; Selkoe, D. J. Soluble protein oligomers in neurodegeneration: lessons from the Alzheimer's amyloid beta-peptide. Nat. Rev. Mol. Cell Biol. 2007, 8, 101-112. 19. Crouch, P. J.; Tew, D. J.; Du, T.; Nguyen, D. N.; Caragounis, A.; Filiz, G.; Blake, R. E.; Trounce, I. A.; Soon, C. P.; Laughton, K.; et al. Restored degradation of the Alzheimer's amyloid-beta peptide by targeting amyloid formation. J. Neurochem. 2009, 108, 1198-1207. 20. Blennow, K.; de Leon, M. J.; Zetterberg, H. Alzheimer's disease. Lancet 2006, 368, 387403.

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 21

21. Mattson, M. P. Pathways towards and away from Alzheimer's disease. Nature 2004, 430, 631-639. 22. Fezoui, Y.; Hartley, D. M.; Harper, J. D.; Khurana, R.; Walsh, D. M.; Condron, M. M.; Selkoe, D. J.; Lansbury, P. T.; Fink, A. L.; Teplow, D. B. An improved method of preparing the amyloid beta-protein for fibrillogenesis and neurotoxicity experiments. Amyloid 2000, 7, 166178. 23. Yates, C. M.; Butterworth, J.; Tennant, M. C.; Gordon, A. Enzyme activities in relation to pH and lactate in postmortem brain in Alzheimer-type and other dementias. J. Neurochem. 1990, 55, 1624-1630. 24. Snyder, S. W.; Ladror, U. S.; Wade, W. S.; Wang, G. T.; Barrett, L. W.; Matayoshi, E. D.; Huffaker, H. J.; Krafft, G. A.; Holzman, T. F. Amyloid-beta aggregation: selective inhibition of aggregation in mixtures of amyloid with different chain lengths. Biophys. J. 1994, 67, 12161228. 25. Barrow, C. J.; Zagorski, M. G. Solution structures of beta peptide and its constituent fragments: relation to amyloid deposition. Science 1991, 253, 179-182. 26. Khandogin, J.; Brooks, C. L. Linking folding with aggregation in Alzheimer's β-amyloid peptides. Proc. Nat. Acad. Sci. U.S.A. 2007, 104, 16880-16885. 27. Ma, K.; Clancy, E. L.; Zhang, Y.; Ray, D. G.; Wollenberg, K.; Zagorski, M. G. Residuespecific pKa measurements of the β-peptide and mechanism of pH-induced amyloid formation. J. Am. Chem. Soc. 1999, 121, 8698-8706. 28. Olubiyi, O. O.; Strodel, B. Structures of the amyloid beta-peptides Abeta1-40 and Abeta1-42 as influenced by pH and a D-peptide. J. Phys. Chem. B 2012, 116, 3280-3291. 29. Opare, S. K. A.; Petoyan, A.; Mehrazma, B.; Rauk, A. Molecular dynamics study of the monomers and dimers of N-AcAβ(13–23)NH2: on the effect of pH on the aggregation of the amyloid beta peptide of Alzheimer’s disease. Can. J. Chem. 2015, 93, 1-9. 30. Henry, B.; Tekely, P.; Delpuech, J. J. pH and pK determinations by high-resolution solidstate 13C NMR: acid-base and tautomeric equilibria of lyophilized L-histidine. J. Am. Chem. Soc. 2002, 124, 2025-2034. 31. Ivanov, I.; Klein, M. L. Deprotonation of a histidine residue in aqueous solution using constrained ab initio molecular dynamics. J. Am. Chem. Soc. 2002, 124, 13380-13381. 32. Huang, Z.; Lin, Z.; Song, C. Protonation processes and electronic spectra of histidine and related ions. J. Phys. Chem. A 2007, 111, 4340-4352. 33. Ghosh, A.; Tucker, M. J.; Gai, F. 2D IR spectroscopy of histidine: probing side-chain structure and dynamics via backbone amide vibrations. J. Phys. Chem. B 2014, 118, 7799-7805. 34. Vila, J. A.; Arnautova, Y. A.; Vorobjev, Y.; Scheraga, H. A. Assessing the fractions of tautomeric forms of the imidazole ring of histidine in proteins as a function of pH. Proc. Natl. Acad. Sci. U.S.A. 2011, 108, 5602-5607. 35. Sugita, Y.; Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 1999, 314, 141-151. 36. Salomon-Ferrer, R.; Götz, A. W.; Poole, D.; Le Grand, S.; Walker, R. C. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. explicit solvent particle mesh ewald. J. Chem. Theory. Comput. 2013, 9, 3878-3888. 37. Case, D. A.; Darden, T. A.; Cheatham, T. E. I.; Simmerling, C. L.; Wang, J.; Duke, R. E.; Luo, R.; Walker, R. C.; Zhang, W.; Merz, K. M.; et al. AMBER12 2012, University of California, San Francisco.

ACS Paragon Plus Environment

18

Page 19 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

38. Qiao, Q.; Bowman, G. R.; Huang, X. Dynamics of an intrinsically disordered protein reveal metastable conformations that potentially seed aggregation. J. Am. Chem. Soc. 2013, 135, 16092-16101. 39. Fawzi, N. L.; Phillips, A. H.; Ruscio, J. Z.; Doucleff, M.; Wemmer, D. E.; Head-Gordon, T. Structure and dynamics of the Aβ(21–30) peptide from the interplay of NMR experiments and molecular simulations. J. Am. Chem. Soc. 2008, 130, 6145-6158. 40. Ning, L.; Pan, D.; Zhang, Y.; Wang, S.; Liu, H.; Yao, X. Effects of the pathogenic mutation A117V and the protective mutation H111S on the folding and aggregation of PrP106126: insights from replica exchange molecular dynamics simulations. PloS one 2015, 10, e0125899. 41. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins: Structure, Function, and Bioinformatics 2006, 65, 712-725. 42. Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graph. 1996, 14, 33-38, 27-28. 43. Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19-25. 44. Vuister, G. W.; Bax, A. Quantitative J correlation: a new approach for measuring homonuclear three-bond J(HNH.alpha.) coupling constants in 15N-enriched proteins. J. Am. Chem. Soc. 1993, 115, 7772-7777. 45. Chiti, F.; Dobson, C. M. Protein misfolding, functional amyloid, and human disease. Annu. Rev. Biochem. 2006, 75, 333-366. 46. Harper, J. D.; Lieber, C. M.; Lansbury, P. T. Atomic force microscopic imaging of seeded fibril formation and fibril branching by the Alzheimer's disease amyloid-β protein. Chem Biol. 1997, 4, 951-959. 47. Serpell, L. C. Alzheimer's amyloid fibrils: structure and assembly. Biochim. Biophys. Acta 2000, 1502, 16-30. 48. Haupt, C.; Leppert, J.; Rönicke, R.; Meinhardt, J.; Yadav, J. K.; Ramachandran, R.; Ohlenschläger, O.; Reymann, K. G.; Görlach, M.; Fändrich, M. Structural basis of β-amyloiddependent synaptic dysfunctions. Angew. Chem. Int. Ed. 2012, 51, 1576-1579. 49. Hoyer, W.; Grönwall, C.; Jonsson, A.; Ståhl, S.; Härd, T. Stabilization of a β-hairpin in monomeric Alzheimer's amyloid-β peptide inhibits amyloid formation. Proc. Nat. Acad. Sci. U.S.A. 2008, 105, 5099-5104. 50. Sgourakis, N. G.; Merced-Serrano, M.; Boutsidis, C.; Drineas, P.; Du, Z.; Wang, C.; Garcia, A. E. Atomic-level characterization of the ensemble of the Aβ(1–42) monomer in water using unbiased molecular dynamics simulations and spectral algorithms. J. Mol. Biol. 2011, 405, 570-583. 51. Sgourakis, N. G.; Yan, Y.; McCallum, S.; Wang, C.; Garcia, A. E. The Alzheimer’s peptides Aβ40 and 42 adopt distinct conformations in water: A combined MD / NMR study. J. Mol. Biol. 2007, 368, 1448-1457. 52. Ball, K. A.; Phillips, A. H.; Nerenberg, P. S.; Fawzi, N. L.; Wemmer, D. E.; HeadGordon, T. Homogeneous and heterogeneous tertiary structure ensembles of amyloid-β peptides. Biochemistry 2011, 50, 7612-7628.

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 21

53. Lin, Y.-S.; Bowman, G. R.; Beauchamp, K. A.; Pande, V. S. Investigating how peptide length and a pathogenic mutation modify the structural ensemble of amyloid beta monomer. Biophys. J. 2012, 102, 315-324. 54. Ball, K. A.; Phillips, A. H.; Wemmer, D. E.; Head-Gordon, T. Differences in β-strand populations of monomeric Aβ40 and Aβ42. Biophys. J. 2013, 104, 2714-2724. 55. Rosenman, D. J.; Connors, C. R.; Chen, W.; Wang, C.; García, A. E. Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach. J. Mol. Biol. 2013, 425, 3338-3359. 56. Yan, Y.; McCallum, S. A.; Wang, C. M35 oxidation induces Aβ40-like structural and dynamical changes in Aβ42. J. Am. Chem. Soc. 2008, 130, 5394-5395. 57. Tauro, S.; Coutinho, E.; Srivastava, S. A rare occurrence of a β-turn in an amyloid βA4 peptide. Magn. Reson. Chem. 2002, 40, 211-218. 58. Mawuenyega, K. G.; Sigurdson, W.; Ovod, V.; Munsell, L.; Kasten, T.; Morris, J. C.; Yarasheski, K. E.; Bateman, R. J. Decreased clearance of CNS β-amyloid in Alzheimer’s disease. Science 2010, 330, 1774-1774. 59. Bateman, R. J.; Munsell, L. Y.; Morris, J. C.; Swarm, R.; Yarasheski, K. E.; Holtzman, D. M. Quantifying CNS protein production and clearance rates in humans using in vivo stable isotope labeling, immunoprecipitation, and tandem mass spectrometry. Nat. Med. 2006, 12, 856861. 60. Selkoe, D. J.; Hardy, J. The amyloid hypothesis of Alzheimer's disease at 25 years. EMBO Mol. Med. 2016, 8, 595-608.

ACS Paragon Plus Environment

20

Page 21 of 21

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC Graphic

ACS Paragon Plus Environment

21