Discrete Molecular Dynamics Simulations Research Articles

Proven Distributed Memory Parallelization of Particle Methods

We provide a mathematically proven parallelization scheme for particle methods on distributed-memory computer systems. Particle methods are a versatile and widely used class of algorithms for computer simulations and numerical predictions in various applications, ranging from continuum fluid dynamics to discrete granular flows and molecular dynamics simulations. Particle methods naturally lend themselves to implementation on parallel computing systems. So far, however, a mathematical proof of correctness and equivalence to sequential implementations was only available for shared-memory parallelism. Here, we leverage a formal definition of the algorithmic class of particle methods to provide a proven parallelization scheme for distributed-memory computers. We prove that thus parallelized particle methods on distributed-memory computers are formally equivalent to their sequential counterpart for a well-defined class of particle methods, and we provide analytical expressions for the speed-up and scalability bounds of this class of algorithms in function of their parameters. The parallelization scheme analyzed here is the basis of many real-world software designs for parallel particle methods. The present analysis is, therefore, of direct relevance to existing and new parallel implementations of particle methods and places them on solid theoretical grounds, rationalizing best practices and providing useful scalability and speedup bounds for benchmarking.

TDP-43 Promotes Amyloid-Beta Toxicity by Delaying Fibril Maturation via Direct Molecular Interaction.

Amyloid-β (Aβ) is a peptide that undergoes self-assembly into amyloid fibrils, which compose the hallmark plaques observed in Alzheimer's disease (AD). TAR DNA-binding protein 43 (TDP-43) is a protein with mislocalization and aggregation implicated in amyotrophic lateral sclerosis and other neurodegenerative diseases. Recent work suggests that TDP-43 may interact with Aβ, inhibiting the formation of amyloid fibrils and worsening AD pathology, but the molecular details of their interaction remain unknown. Using all-atom discrete molecular dynamics simulations, we systematically investigated the direct molecular interaction between Aβ and TDP-43. We found that Aβ monomers were able to bind near the flexible nuclear localization sequence of the N-terminal domain (NTD) of TDP-43, adopting β-sheet rich conformations that were promoted by the interaction. Furthermore, Aβ associated with the nucleic acid binding interface of the tandem RNA recognition motifs of TDP-43 via electrostatic interactions. Using the computational peptide array method, we found the strongest C-terminal domain interaction with Aβ to be within the amyloidogenic core region of TDP-43. With experimental evidence suggesting that the NTD is necessary for inhibiting Aβ fibril growth, we also simulated the NTD with an Aβ40 fibril seed. We found that the NTD was able to strongly bind the elongation surface of the fibril seed via extensive hydrogen bonding and could also diffuse along the lateral surface via electrostatic interactions. Our results suggest that TDP-43 binding to the elongation surface, thereby sterically blocking Aβ monomer addition, is responsible for the experimentally observed inhibition of fibril growth. We conclude that TDP-43 may promote Aβ toxicity by stabilizing the oligomeric state and kinetically delaying fibril maturation.

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.

Experimental and computational investigation of the effect of Hsc70 structural variants on inhibiting amylin aggregation

The misfolding and aggregation of human islet amyloid polypeptide (hIAPP), also known as amylin, have been implicated in the pathogenesis of type 2 diabetes (T2D). Heat shock proteins, specifically, heat shock cognate 70 (Hsc70), are molecular chaperones that protect against hIAPP misfolding and inhibits its aggregation. Nevertheless, there is an incomplete understanding of the mechanistic interactions between Hsc70 domains and hIAPP, thus limiting their potential therapeutic role in diabetes. This study investigates the inhibitory capacities of different Hsc70 variants, aiming to identify the structural determinants that strike a balance between efficacy and cytotoxicity. Our experimental findings demonstrate that the ATPase activity of Hsc70 is not a pivotal factor for inhibiting hIAPP misfolding. We underscore the significance of the C-terminal substrate-binding domain of Hsc70 in inhibiting hIAPP aggregation, emphasizing that the removal of the lid subdomain diminishes the inhibitory effect of Hsc70. Additionally, we employed atomistic discrete molecular dynamics simulations to gain deeper insights into the interaction between Hsc70 variants and hIAPP. Integrating both experimental and computational findings, we propose a mechanism by which Hsc70's interaction with hIAPP monomers disrupts protein-protein connections, primarily by shielding the β-sheet edges of the Hsc70-β-sandwich. The distinctive conformational dynamics of the alpha helices of Hsc70 potentially enhance hIAPP binding by obstructing the exposed edges of the β-sandwich, particularly at the β5-β8 region along the alpha helix interface. This, in turn, inhibits fibril growth, and similar results were observed following hIAPP dimerization. Overall, this study elucidates the structural intricacies of Hsc70 crucial for impeding hIAPP aggregation, improving our understanding of the potential anti-aggregative properties of molecular chaperones in diabetes treatment.

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.

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.

Domain tethering impacts dimerization and DNA-mediated allostery in the human transcription factor FoxP1.

Transcription factors are multidomain proteins with specific DNA binding and regulatory domains. In the human FoxP subfamily (FoxP1, FoxP2, FoxP3, and FoxP4) of transcription factors, a 90 residue-long disordered region links a Leucine Zipper (ZIP)-known to form coiled-coil dimers-and a Forkhead (FKH) domain-known to form domain swapping dimers. We used replica exchange discrete molecular dynamics simulations, single-molecule fluorescence experiments, and other biophysical tools to understand how domain tethering in FoxP1 impacts dimerization at ZIP and FKH domains and how DNA binding allosterically regulates their dimerization. We found that domain tethering promotes FoxP1 dimerization but inhibits a FKH domain-swapped structure. Furthermore, our findings indicate that the linker mediates the mutual organization and dynamics of ZIP and FKH domains, forming closed and open states with and without interdomain contacts, thus highlighting the role of the linkers in multidomain proteins. Finally, we found that DNA allosterically promotes structural changes that decrease the dimerization propensity of FoxP1. We postulate that, upon DNA binding, the interdomain linker plays a crucial role in the gene regulatory function of FoxP1.

Differential Binding and Conformational Dynamics of Tau Microtubule-Binding Repeats with a Preformed Amyloid-β Fibril Seed.

Both senile plaques formed by amyloid-β (Aβ) and neurofibrillary tangles (NFTs) comprised of tau are pathological hallmarks of Alzheimer's disease (AD). The accumulation of NFTs better correlates with the loss of cognitive function than senile plaques, but NFTs are rarely observed without the presence of senile plaques. Hence, cross-seeding of tau by preformed Aβ amyloid fibril seeds has been proposed to drive the aggregation of tau and exacerbate AD progression, but the molecular mechanism remains unknown. Here, we first identified cross-interaction hotspots between Aβ and tau using atomistic discrete molecular dynamics simulations (DMD) and confirmed the critical role of the four microtubule-binding repeats of tau (R1-R4) in the cross-interaction with Aβ. We further investigated the binding structure and dynamics of each tau repeat with a preformed Aβ fibril seed. Specifically, R1 and R3 preferred to bind the Aβ fibril lateral surface instead of the elongation end. In contrast, R2 and R4 had higher binding propensities to the fibril elongation end than the lateral surface, enhancing β-sheet content by forming hydrogen bonds with the exposed hydrogen bond donors and acceptors. Together, our results suggest that the four repeats play distinct roles in driving the binding of tau to different surfaces of an Aβ fibril seed. Binding of tau to the lateral surface of Aβ fibril can increase the local concentration, while the binding to the elongation surface promotes β-sheet formation, both of which reduce the free energy barrier for tau aggregation nucleation and subsequent fibrillization.

Structural and dynamic insights into α-synuclein dimer conformations

Parkinson disease is associated with the aggregation of the protein α-synuclein. While α-synuclein can exist in multiple oligomeric states, the dimer has been a subject of extensive debates. Here, using an array of biophysical approaches, we demonstrate that α-synuclein invitro exhibits primarily a monomer-dimer equilibrium in nanomolar concentrations and up to a few micromolars. We then use spatial information from hetero-isotopic cross-linking mass spectrometry experiments as restrains in discrete molecular dynamics simulations to obtain the ensemble structure of dimeric species. Out of eight structural sub-populations of dimers, we identify one that is compact, stable, abundant, and exhibits partially exposed β-sheet structures. This compact dimer is the only one where the hydroxyls of tyrosine 39 are in proximity that may promote dityrosine covalent linkage upon hydroxyl radicalization, which is implicated in α-synuclein amyloid fibrils. We propose that this α-synuclein dimer features etiological relevance to Parkinson disease.

Zinc-Epigallocatechin-3-gallate Network-Coated Nanocomposites against the Pathogenesis of Amyloid-Beta.

The aggregation of amyloid beta (Aβ) is a hallmark of Alzheimer's disease (AD), a major cause of dementia and an unmet challenge in modern medicine. In this study, we constructed a biocompatible metal-phenolic network (MPN) comprising a polyphenol epigallocatechin gallate (EGCG) scaffold coordinated by physiological Zn(II). Upon adsorption onto gold nanoparticles, the MPN@AuNP nanoconstruct elicited a remarkable potency against the amyloid aggregation and toxicity of Aβ in vitro. The superior performance of MPN@AuNP over EGCG@AuNP was attributed to the porosity and hence larger surface area of the MPN in comparison with that of EGCG alone. The atomic detail of Zn(II)-EGCG coordination was unraveled by density functional theory calculations and the structure and dynamics of Aβ aggregation modulated by the MPN were further examined by discrete molecular dynamics simulations. As MPN@AuNP also displayed a robust capacity to cross a blood-brain barrier model through the paracellular pathway, and given the EGCG's function as an anti-amyloidosis and antioxidation agent, this MPN-based strategy may find application in regulating the broad AD pathology beyond protein aggregation inhibition.

Direct observation of single-molecular hIAPP fibril elongation reveal molecular mechanisms of morphological dependence and asymmetric growth.

Rapid growth of amyloid fibrils via seeded conformational conversion of monomers, is a critical step of fibrillization and important for the prion-like transmission of amyloid diseases. Amyloid fibrils are often polymorphic with relative populations morphology-dependent, but the molecular mechanisms of fibril elongation and corresponding morphological dependence remain poorly understood. Here, we applied replica exchange all-atom discrete molecular dynamics simulations to study the conformational conversion of hIAPP, an amyloid peptide associated with type-2 diabetes, catalyzed by two representative fibrils with distinct morphologies where hIAPP adopted S-shaped and extended conformations correspondingly. For both cases, the incorporation of monomers into the preformed fibrils with parallel in-register alignments was observed in our unbiased simulations. The fibril elongation free energy landscape derived from simulations is characterized by multiple basins corresponding to different intermediate states. Fibril morphology affected monomer binding at fibril elongation and lateral surfaces, and also the conformational conversion dynamics at the elongation surface, resulting into different elongation rates and thus relative populations for different fibril morphologies as observed experimentally. We also investigated the dynamics of hIAPP monomers bound to the fibril elongation surfaces at opposite ends. We observed distinct free energy landscapes in agreement with experimentally-observed asymmetric fibril growth, which can be attributed to subtle physicochemical difference of two ends in terms of contacting surfaces and patterns of available hydrogen bond donors and acceptors. Together, our results offered molecular insights to seeded conformational conformation of fibril elongation, morphological dependence, asymmetry fibril growth, and could shed light on amyloid strains with distinct disease progression.

Differential binding dynamics of four microtubule-binding repeats of tau with an amyloid-β fibril seed.

Both senile plaques formed by amyloid-beta (Aβ) and neurofibrillary tangles comprised of tau are hallmarks of Alzheimer's disease (AD). The accumulation of neurofibrillary tangles better correlates with the loss of cognitive function than senile plaques, but neurofibrillary tangles are rarely observed without the presence of senile plaques. Hence, cross-seeding of tau by preformed Aβ amyloid fibril seeds has been proposed to drive the aggregation of tau and exacerbate AD progression, but the molecular mechanism remains poorly understood. Herein, we identified cross-interaction hotspots between Aβ and tau using atomistic discrete molecular dynamics simulations (DMD) and confirmed the critical role of the four microtubule-binding repeats of tau (R1-R4) in the cross-interaction with Aβ. We further investigated the binding structure and dynamics of each tau repeat with a preformed Aβ fibril seed. Specifically, R1 and R3 preferred to bind the lateral instead of elongation surfaces of the fibril, preserving similar secondary structures as in the solution. In contrast, R2 and R4 had higher binding propensities to the elongation than lateral surfaces, adopting β-sheet structures by forming hydrogen bonds with the exposed hydrogen bond donors and acceptors. Together, our results suggest that four MTBD repeats play distinct roles in driving the binding of tau to two different surfaces of an Aβ fibril seed. The binding of tau to the lateral surface of Aβ fibril can increase the local concentration, while the binding in the elongation surface promotes β-sheet formation, both of which reduce the free energy barrier for tau aggregation nucleation and subsequent fibrillization. Our computational study offers molecular insight into the experimentally observed cross-seeding of Tau aggregation by pre-formed Aβ fibril seeds and helps better understand the pathological synergy of Aβ and tau aggregation in the etiology of AD.

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

Human calcitonin (hCT) is a polypeptide hormone that participates in calcium-phosphorus metabolism. Irreversible aggregation of 32-amino acid hCT into β-sheet-rich amyloid fibrils impairs physiological activity and increases the risk of medullary carcinoma of the thyroid. Amyloid-resistant hCT derivatives substituting critical amyloidogenic residues are of particular interest for clinical applications as therapeutic drugs against bone-related diseases. Uncovering the aggregation mechanism of hCT at the molecular level, therefore, is important for the design of amyloid-resistant hCT analogues. Here, we investigated the aggregation dynamics of hCT, non-amyloidogenic salmon calcitonin (sCT), and two hCT analogues with reduced aggregation tendency─TL-hCT and phCT─using long timescale discrete molecular dynamics simulations. Our results showed that hCT monomers mainly adopted unstructured conformations with dynamically formed helices around the central region. hCT self-assembled into helix-rich oligomers first, followed by a conformational conversion into β-sheet-rich oligomers with β-sheets formed by residues 10-30 and stabilized by aromatic and hydrophobic interactions. Our simulations confirmed that TL-hCT and phCT oligomers featured more helices and fewer β-sheets than hCT. Substitution of central aromatic residues with leucine in TL-hCT and replacing C-terminal hydrophobic residue with hydrophilic amino acid in phCT only locally suppressed β-sheet propensities in the central region and C-terminus, respectively. Having mutations in both central and C-terminal regions, sCT monomers and dynamically formed oligomers predominantly adopted helices, confirming that both central aromatic and C-terminal hydrophobic residues played important roles in the fibrillization of hCT. We also observed the formation of β-barrel intermediates, postulated as the toxic oligomers in amyloidosis, for hCT but not for sCT. Our computational study depicts a complete picture of the aggregation dynamics of hCT and the effects of mutations. The design of next-generation amyloid-resistant hCT analogues should consider the impact on both amyloidogenic regions and also take into account the amplification of transient β-sheet population in monomers upon aggregation.

Molecular Insights into the Misfolding and Dimerization Dynamics of the Full-Length α-Synuclein from Atomistic Discrete Molecular Dynamics Simulations.

The misfolding and pathological aggregation of α-synuclein forming insoluble amyloid deposits is associated with Parkinson's disease, the second most common neurodegenerative disease in the world population. Characterizing the self-assembly mechanism of α-synuclein is critical for discovering treatments against synucleinopathies. The intrinsically disordered property, high degrees of freedom, and macroscopic timescales of conformational conversion make its characterization extremely challenging in vitro and in silico. Here, we systematically investigated the dynamics of monomer misfolding and dimerization of the full-length α-synuclein using atomistic discrete molecular dynamics simulations. Our results suggested that both α-synuclein monomers and dimers mainly adopted unstructured formations with partial helices around the N-terminus (residues 8-32) and various β-sheets spanning the residues 35-56 (N-terminal tail) and residues 61-95 (NAC region). The C-terminus mostly assumed an unstructured formation wrapping around the lateral surface and the elongation edge of the β-sheet core formed by an N-terminal tail and NAC regions. Dimerization enhanced the β-sheet formation along with a decrease in the unstructured content. The inter-peptide β-sheets were mainly formed by the N-terminal tail and NACore (residues 68-78) regions, suggesting that these two regions played critical roles in the amyloid aggregation of α-synuclein. Interactions of the C-terminus with the N-terminal tail and the NAC region were significantly suppressed in the α-synuclein dimer, indicating that the interaction of the C-terminus with the N-terminal tail and NAC regions could prevent α-synuclein aggregation. These results on the structural ensembles and early aggregation dynamics of α-synuclein will help understand the nucleation and fibrillization of α-synuclein.

Activation of Disulfide Redox Switch in REDD1 Promotes Oxidative Stress Under Hyperglycemic Conditions.

The stress response protein regulated in development and DNA damage response 1 (REDD1) has been implicated in visual deficits in patients with diabetes. The aim here was to investigate the mechanism responsible for the increase in retinal REDD1 protein content that is observed with diabetes. We found that REDD1 protein expression was increased in the retina of streptozotocin-induced diabetic mice in the absence of a change in REDD1 mRNA abundance or ribosome association. Oral antioxidant supplementation reduced retinal oxidative stress and suppressed REDD1 protein expression in the retina of diabetic mice. In human retinal Müller cell cultures, hyperglycemic conditions increased oxidative stress, enhanced REDD1 expression, and inhibited REDD1 degradation independently of the proteasome. Hyperglycemic conditions promoted a redox-sensitive cross-strand disulfide bond in REDD1 at C150/C157 that was required for reduced REDD1 degradation. Discrete molecular dynamics simulations of REDD1 structure revealed allosteric regulation of a degron upon formation of the disulfide bond that disrupted lysosomal proteolysis of REDD1. REDD1 acetylation at K129 was required for REDD1 recognition by the cytosolic chaperone HSC70 and degradation by chaperone-mediated autophagy. Disruption of REDD1 allostery upon C150/C157 disulfide bond formation prevented the suppressive effect of hyperglycemic conditions on REDD1 degradation and reduced oxidative stress in cells exposed to hyperglycemic conditions. The results reveal redox regulation of REDD1 and demonstrate the role of a REDD1 disulfide switch in development of oxidative stress.

Substoichiometric Inhibition of Insulin against IAPP Aggregation Is Attenuated by the Incompletely Processed N-Terminus of proIAPP.

Substoichiometric aggregation inhibition of human islet amyloid polypeptide (IAPP), the hallmark of type 2 diabetes impacting millions of people, is crucial for developing clinic therapies, yet it remains challenging given that many candidate inhibitors require high doses. Intriguingly, insulin, the key regulatory polypeptide on blood glucose levels that are cosynthesized, costored, and cosecreted with IAPP by pancreatic β cells, has been identified as a potent inhibitor that can suppress IAPP amyloid aggregation at substoichiometric concentrations. Here, we computationally investigated the molecular mechanisms of the substoichiometric inhibition of insulin against the aggregation of IAPP and the incompletely processed IAPP (proIAPP) using discrete molecular dynamics simulations. Our results suggest that the amyloid aggregations of both IAPP and proIAPP might be disrupted by insulin through its binding with the shared amyloidogenic core sequences. However, the N-terminus of proIAPP competed with the amyloidogenic core sequences for the insulin interactions, resulting in attenuated inhibition by insulin. Moreover, insulin preferred to bind the elongation surfaces of IAPP seeds with fibril-like structure, with a stronger affinity than that of IAPP monomers. The capping of elongation surfaces by a small amount of insulin sterically prohibited the seed growth via monomer addition, achieving the substoichiometric inhibition. Together, our computational results provided molecular insights for the substoichiometric inhibition of insulin against IAPP aggregation, also the weakened effect on proIAPP. The uncovered substoichiometric inhibition by capping the elongation of amyloid seeds or fibrils may guide the rational designs of new potent inhibitors effective at low doses.

Multifunctions of α‐Synuclein Explained by Its Dynamic Heterogeneous Conformations with a Hierarchy of Transition Times

α‐Synuclein plays essential roles in synaptic vesicle homeostasis and neurotransmitter release through the interaction with membrane, while α‐synuclein aggregates have been associated with several neurodegenerative diseases, especially Parkinson’s disease. The different roles of α‐synuclein could be attributed to the heterogenous conformations adopted by this intrinsically disordered protein. To demystify the physiological functions of α‐synuclein and understand its pathogenic mechanism of Parkinson’s disease, characterizing the monomeric structural ensemble and identifying the aggregation‐prone and non‐aggregation‐prone structures are of particular importance. Here, we use inter‐residue distance distributions derived from time‐resolved FRET experiments as constraints to guide discrete molecular dynamics simulations of α‐synuclein monomer. We explore the conformational space of α‐synuclein and verify the generated conformational ensemble by additional experiments including far‐UV circular dichroism spectrum and cross‐linking mass spectrometry. We find that some conformational states of α‐synuclein are surprisingly stable displaying dynamic transitions less than milliseconds. A comprehensive analysis of the conformational ensemble uncovers important structural features and potential conformations that are critical to stabilize the monomeric state or induce different oligomerization pathways.