Structural Perturbation of Monomers Determines the Amyloid Aggregation Propensity of Calcitonin Variants

Computational Investigation of Coaggregation and Cross-Seeding between Aβ and hIAPP Underpinning the Cross-Talk in Alzheimer's Disease and Type 2 Diabetes

The coexistence of amyloid-β (Aβ) and human islet amyloid polypeptide (hIAPP) in the brain and pancreas is associated with an increased risk of Alzheimer's disease (AD) and type 2 diabetes (T2D) due to their coaggregation and cross-seeding. Despite this, the molecular mechanisms underlying their interaction remain elusive. Here, we systematically investigated the cross-talk between Aβ and hIAPP using atomistic discrete molecular dynamics (DMD) simulations. Our results revealed that the amyloidogenic core regions of both Aβ (Aβ10-21 and Aβ30-41) and hIAPP (hIAPP8-20 and hIAPP22-29), driving their self-aggregation, also exhibited a strong tendency for cross-interaction. This propensity led to the formation of β-sheet-rich heterocomplexes, including potentially toxic β-barrel oligomers. The formation of Aβ and hIAPP heteroaggregates did not impede the recruitment of additional peptides to grow into larger aggregates. Our cross-seeding simulations demonstrated that both Aβ and hIAPP fibrils could mutually act as seeds, assisting each other's monomers in converting into β-sheets at the exposed fibril elongation ends. The amyloidogenic core regions of Aβ and hIAPP, in both oligomeric and fibrillar states, exhibited the ability to recruit isolated peptides, thereby extending the β-sheet edges, with limited sensitivity to the amino acid sequence. These findings suggest that targeting these regions by capping them with amyloid-resistant peptide drugs may hold potential as a therapeutic approach for addressing AD, T2D, and their copathologies.

Molecular Insights into the Effects of F16L and F19L Substitutions on the Conformation and Aggregation Dynamics of Human Calcitonin

Human calcitonin (hCT) regulates calcium-phosphorus metabolism, but its amyloid aggregation disrupts physiological activity, increases thyroid carcinoma risk, and hampers its clinical use for bone-related diseases like osteoporosis and Paget's disease. Improving hCT with targeted modifications to mitigate amyloid formation while maintaining its function holds promise as a strategy. Understanding how each residue in hCT's amyloidogenic core affects its structure and aggregation dynamics is crucial for designing effective analogues. Mutants F16L-hCT and F19L-hCT, where Phe residues in the core are replaced with Leu as in nonamyloidogenic salmon calcitonin, showed different aggregation kinetics. However, the molecular effects of these substitutions in hCT are still unclear. Here, we systematically investigated the folding and self-assembly conformational dynamics of hCT, F16L-hCT, and F19L-hCT through multiple long-time scale independent atomistic discrete molecular dynamics (DMD) simulations. Our results indicated that the hCT monomer primarily assumed unstructured conformations with dynamic helices around residues 4-12 and 14-21. During self-assembly, the amyloidogenic core of hCT14-21 converted from dynamic helices to β-sheets. However, substituting F16L did not induce significant conformational changes, as F16L-hCT exhibited characteristics similar to those of wild-type hCT in both monomeric and oligomeric states. In contrast, F19L-hCT exhibited substantially more helices and fewer β-sheets than did hCT, irrespective of their monomers or oligomers. The substitution of F19L significantly enhanced the stability of the helical conformation for hCT14-21, thereby suppressing the helix-to-β-sheet conformational conversion. Overall, our findings elucidate the molecular mechanisms underlying hCT aggregation and the effects of F16L and F19L substitutions on the conformational dynamics of hCT, highlighting the critical role of F19 as an important target in the design of amyloid-resistant hCT analogs for future clinical applications.

Molecular Insights into the Differential Effects of Acetylation on the Aggregation of Tau Microtubule-Binding Repeats

Aggregation of tau protein into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The microtubule-binding (MTB) domain of tau, containing either three or four repeats with sequence similarities, plays an important role in determining tau's aggregation. Previous studies have reported that abnormal acetylation of lysine residues displays a distinct effect on the formation of pathological tau aggregates. However, the underlying molecular mechanism remains mostly elusive. In this study, we performed extensive replica exchange molecular dynamics (REMD) simulations of 144 μs in total to systematically investigate the dimerization of four tau MTB repeats and explore the impacts of Lys280 (K280) or Lys321 (K321) acetylation on the conformational ensembles of the R2 or R3 dimer. Our results show that R3 is the most prone to aggregation among the four repeats, followed by R2 and R4, while R1 displays the weakest aggregation propensity with a disordered structure. Acetylation of K280 could promote the aggregation of R2 peptides by increasing the formation of β-sheet structures and strengthening the interchain interaction. However, K321 acetylation decreases the β-sheet content of the R3 dimer, reduces the ability of R3 peptides to form long β-strands, and promotes the stable helix structure formation. The salt bridge and Y310-Y310 π-π stacking interactions of the R3 dimer are greatly weakened by K321 acetylation, resulting in the inhibition of dimerization. This study uncovers the structural ensembles of tau MTB repeats and provides mechanistic insights into the influences of acetylation on tau aggregation, which may deepen the understanding of the pathogenesis of tauopathies.

Unraveling the underlying mechanisms of reduced amyloidogenic properties in human calcitonin via double mutations

The therapeutic efficacy of peptide-based drugs is commonly hampered by the intrinsic propensity to aggregation. A notable example is human calcitonin (hCT), a peptide hormone comprising 32 amino acids, which is synthesized and secreted by thyroid gland parafollicular cells (C cells). This hormone plays a vital role in regulating blood calcium levels and upholding bone integrity. Despite its physiological importance, utilizing hCT as a drug is hampered by its inclination to form amyloid. To address this limitation, an alternative is provided by salmon calcitonin (sCT), which possesses a lower aggregation propensity. Although sharing the same disulfide bond at the N terminus as hCT, sCT differs from hCT at a total of 16 amino acid positions. However, due to the dissimilarity in sequences, using sCT as a clinical replacement occasionally results in adverse side effects in patients. Earlier investigations have highlighted the significant roles of Tyr-12 and Asn-17 in inducing the formation of amyloid fibrils. By introducing double mutations at these sites, the ability to hinder aggregation can be significantly augmented. This study delves into the oligomerization and helical structure formation of the hCT double mutant (Y12LN17H hCT, noted as DM hCT), as well as two single mutants (Y12L and N17H), aiming to elucidate the mechanism behind hCT fibrillization. In addition, computational prediction tools were employed again to identify potential substitutes. Although the results yielded were not entirely satisfactory, a comparison between the newly examined and previously found hCT double mutants provides insights into the reduced aggregation propensity of the latter. This research endeavor holds the promise of informing the design of more effective therapeutic peptide drugs in the future.

Computational insights into the cross-talk between medin and Aβ: implications for age-related vascular risk factors in Alzheimer's disease

The aggregation of medin forming aortic medial amyloid is linked to arterial wall degeneration and cerebrovascular dysfunction. Elevated levels of arteriolar medin are correlated with an increased presence of vascular amyloid-β (Aβ) aggregates, a hallmark of Alzheimer's disease (AD) and vascular dementia. The cross-interaction between medin and Aβ results in the formation of heterologous fibrils through co-aggregation and cross-seeding processes both in vitro and in vivo. However, a comprehensive molecular understanding of the cross-interaction between medin and Aβ-two intrinsically disordered proteins-is critically lacking. Here, we employed atomistic discrete molecular dynamics simulations to systematically investigate the self-association, co-aggregation and also the phenomenon of cross-seeding between these two proteins. Our results demonstrated that both Aβ and medin were aggregation prone and their mixture tended to form β-sheet-rich hetero-aggregates. The formation of Aβ-medin hetero-aggregates did not hinder Aβ and medin from recruiting additional Aβ and medin peptides to grow into larger β-sheet-rich aggregates. The β-barrel oligomer intermediates observed in the self-aggregations of Aβ and medin were also present during their co-aggregation. In cross-seeding simulations, preformed Aβ fibrils could recruit isolated medin monomers to form elongated β-sheets. Overall, our comprehensive simulations suggested that the cross-interaction between Aβ and medin may contribute to their pathological aggregation, given the inherent amyloidogenic tendencies of both medin and Aβ. Targeting medin, therefore, could offer a novel therapeutic approach to preserving brain function during aging and AD by improving vascular health.

Deciphering the influence of Y12L and N17H substitutions on the conformation and oligomerization of human calcitonin

The abnormal aggregation of human calcitonin (hCT) hormone peptides impairs their physiological function, leading to harmful immune responses and cytotoxicity, which limits their clinical utility. Interestingly, a representative hCT analog incorporating Y12L and N17H substitutions (DM-hCT) has shown reduced aggregation tendencies while maintaining bioactivity. But the molecular mechanism of Y12L and N17H substitutions on the conformational dynamics of hCT remains unclear. Here, we systematically investigated the folding and self-assembly dynamics of hCT and DM-hCT using atomistic discrete molecular dynamics (DMD) simulations. Our findings revealed that hCT monomers predominantly adopted unstructured conformations with dynamic helices. Oligomerization of hCT resulted in the formation of β-sheet-rich aggregates and β-barrel intermediates. The Y12L and N17H substitutions enhanced helical conformations and suppressed β-sheet formation in both monomers and oligomers. These substitutions stabilized the dynamic helices and disrupted aromatic interactions responsible for β-sheet formation at residue 12. Notably, DM-hCT assemblies still exhibited β-sheets in phenylalanine-rich and C-terminal hydrophobic regions, suggesting that future optimizations should focus on these areas. Our simulations provide insights into the molecular mechanisms underlying hCT aggregation and the amyloid-resistant effects of Y12L and N17H substitutions. These findings have valuable implications for the development of clinical hCT analogs.

Co-aggregation of α-synuclein with amyloid-β stabilizes β-sheet-rich oligomers and enhances the formation of β-barrels

Alzheimer's disease (AD) and Parkinson's disease (PD) are the two most common neurodegenerative diseases with markedly different pathological features of β-amyloid (Aβ) plaques and α-synuclein (αS) Lewy bodies (LBs), respectively. However, clinical overlaps in symptoms and pathologies between AD and PD are commonly observed caused by the cross-interaction between Aβ and αS. To uncover the molecular mechanisms behind their overlapping symptoms and pathologies, we computationally investigated the impact of αS on an Aβ monomer and dimerization using atomistic discrete molecular dynamics simulations (DMD). Our results revealed that αS could directly interact with Aβ monomers and dimers, thus forming β-sheet-rich oligomers, including potentially toxic β-barrel intermediates. The binding hotspot involved the second half of the N-terminal domain and NAC region in αS, along with residues 10-21 and 31-42 in Aβ. In their hetero-complex, the binding hotspot primarily assumed a β-sheet core buried inside, which was dynamically shielded by the highly charged, amyloid-resistant C-terminus of αS. Because the amyloid prion region was the same as the binding hotspot being buried, their fibrillization may be delayed, causing the toxic oligomers to increase. This study sheds light on the intricate relationship between Aβ and αS and provides insights into the overlapping pathology of AD and PD.

Unveiling Medin Folding and Dimerization Dynamics and Conformations via Atomistic Discrete Molecular Dynamics Simulations

Medin is a principal component of localized amyloid found in the vasculature of individuals over 50 years old. Its amyloid aggregation has been linked to endothelial dysfunction and vascular inflammation, contributing to the pathogenesis of various vascular diseases. Despite its significance, the structures of the medin monomer, oligomer, and fibril remain elusive, and the dynamic processes of medin aggregation are not fully understood. In this study, we comprehensively investigated the medin folding and dimerization dynamics and conformations using atomistic discrete molecular dynamics simulations. Our simulation results suggested that the folding initiation of the medin involved the formation of β-sheets around medin30-41 and medin42-50, with subsequent capping of other segments to their β-sheet edges. Medin monomers typically consisted of three or four β-strands, along with a dynamic N-terminal helix. Two isolated medin peptides readily aggregated into a β-sheet-rich dimer, displaying a strong aggregation propensity. Dimerization of medin not only enhanced the β-sheet conformations but also led to the formation of β-barrel oligomers. The aggregation tendencies of medin1-18 and medin19-29 were relatively weak. However, the segments of medin30-41 and medin42-50 played a crucial role as they primarily formed a β-sheet core and facilitated medin1-18 and medin19-29 to form intra- and interpeptide β-sheets. The findings highlight the critical role of the medin30-41 and medin42-50 regions in stabilizing the monomer structure and driving the medin amyloid aggregation. These regions could potentially serve as promising targets for designing antiamyloid inhibitors against amyloid aggregation of medin. Additionally, our study provides a full picture of the monomer conformations and dimerization dynamics for medin, which will help better understand the pathology of medin aggregation.

Dissecting the Self-assembly Dynamics of Imperfect Repeats in α-Synuclein

The pathological aggregation of α-synuclein (αS) into amyloid fibrils is the hallmark of Parkinson's disease (PD). The self-assembly and membrane interactions of αS are mainly governed by the seven imperfect 11-residue repeats of the XKTKEGVXXXX motif around residues 1-95. However, the particular role of each repeat in αS fibrillization remains unclear. To answer this question, we studied the aggregation dynamics of each repeat with up to 10 peptides in silico by conducting multiple independent micro-second atomistic discrete molecular dynamics simulations. Our simulations revealed that only repeats R3 and R6 readily self-assembled into β-sheet-rich oligomers, while the other repeats remained as unstructured monomers with weak self-assembly and β-sheet propensities. The self-assembly process of R3 featured frequent conformational changes with β-sheet formation mainly in the non-conserved hydrophobic tail, whereas R6 spontaneously self-assembled into extended and stable cross-β structures. These results of seven repeats are consistent with their structures and organization in recently solved αS fibrils. As the primary amyloidogenic core, R6 was buried inside the central cross-β core of all αS fibrils, attracting the hydrophobic tails of adjacent R4, R5, and R7 repeats forming β-sheets around R6 in the core. Further away from R6 in the sequence but with a moderate amyloid aggregation propensity, the R3 tail could serve as a secondary amyloidogenic core and form independent β-sheets in the fibril. Overall, our results demonstrate the critical role of R3 and R6 repeats in αS amyloid aggregation and suggest their potential as targets for the peptide-based and small-molecule amyloid inhibitors.

Evolution of large Aβ16-22 aggregates at atomic details and potential of mean force associated to peptide unbinding and fragmentation events

Atomic characterization of large nonfibrillar aggregates of amyloid polypeptides cannot be determined by experimental means. Starting from β-rich aggregates of Y and elongated topologies predicted by coarse-grained simulations and consisting of more than 100 Aβ16-22 peptides, we performed atomistic molecular dynamics (MD), replica exchange with solute scaling (REST2), and umbrella sampling simulations using the CHARMM36m force field in explicit solvent. Here, we explored the dynamics within 3 μs, the free energy landscape, and the potential of mean force associated with either the unbinding of one single peptide in different configurations within the aggregate or fragmentation events of a large number of peptides. Within the time scale of MD and REST2, we find that the aggregates experience slow global conformational plasticity, and remain essentially random coil though we observe slow beta-strand structuring with a dominance of antiparallel beta-sheets over parallel beta-sheets. Enhanced REST2 simulation is able to capture fragmentation events, and the free energy of fragmentation of a large block of peptides is found to be similar to the free energy associated with fibril depolymerization by one chain for longer Aβ sequences.