Assembly of Stefin B into Polymorphic Oligomers Probed by Discrete Molecular Dynamics

Human stefin B: from its structure, folding, and aggregation to its function in health and disease

Mutations in the gene for human stefin B (cystatin B) cause progressive myoclonic epilepsy type 1 (EPM1), a neurodegenerative disorder. The most common change is dodecamer repeats in the promoter region of the gene, though missense and frameshift mutations also appear. Human stefin B primarily acts as a cysteine cathepsin inhibitor, and it also exhibits alternative functions. It plays a protective role against oxidative stress, likely via reducing mitochondrial damage and thus generating fewer mitochondrial reactive oxygen species (ROS). Accordingly, lack of stefin B results in increased inflammation and NLRP3 inflammasome activation, producing more ROS. The protein is cytosolic but also has an important role in the nucleus, where it prevents cleavage of the N terminal part of histone 3 by inhibiting cathepsins L and B and thus regulates transcription and cell cycle. Furthermore, it has been shown that stefin B is oligomeric in cells and that it has a specific role in the physiology of the synapse and in vesicular transport. On the basis of my research team's data on the structure, folding, and aggregation of stefin B, we have proposed that it might regulate proteostasis, possessing a chaperone-like function. In this review, I synthesize these observations and derive some conclusions on possible sources of EPM1 pathology. The interaction partners of stefin B and other gene mutations leading to EPM1-like pathology are discussed and common pathways are pinpointed.

Soluble State of Villin Headpiece Protein as a Tool in the Assessment of MD Force Fields

Protein self-assembly plays an important role in cellular processes. Whereas molecular dynamics (MD) represents a powerful tool in studying assembly mechanisms, its predictions depend on the accuracy of underlying force fields, which are known to overly promote protein assembly. We here examine villin headpiece domain, HP36, which remains soluble at concentrations amenable to MD studies. The experimental characterization of soluble HP36 at concentrations of 0.05 to 1 mM reveals concentration-independent 90% monomeric and 10% dimeric populations. Extensive all-atom MD simulations at two protein concentrations, 0.9 and 8.5 mM, probe the HP36 dimer population, stability, and kinetics of dimer formation within two MD force fields, Amber ff14SB and CHARMM36m. MD results demonstrate that whereas CHARMM36m captures experimental HP36 monomer populations at the lower concentration, both force fields overly promote HP36 association at the higher concentration. Moreover, contacts stabilizing HP36 dimers are force-field-dependent. CHARMM36m produces consistently higher HP36 monomer populations, lower association rates, and weaker dependence of these quantities on the protein concentration than Amber ff14SB. Nonetheless, the highest monomer populations and dissociation constants are observed when the TIP3P water model in Amber ff14SB is replaced by TIP4P/2005, showcasing the critical role of the water model in addressing the protein solubility problem in MD.

Nucleation of β-rich oligomers and β-barrels in the early aggregation of human islet amyloid polypeptide

The self-assembly of human islet amyloid polypeptide (hIAPP) into β-sheet rich amyloid aggregates is associated with pancreatic β-cell death in type 2 diabetes (T2D). Prior experimental studies of hIAPP aggregation reported the early accumulation of α-helical intermediates before the rapid conversion into β-sheet rich amyloid fibrils, as also corroborated by our experimental characterizations with transmission electron microscopy and Fourier transform infrared spectroscopy. Although increasing evidence suggests that small oligomers populating early hIAPP aggregation play crucial roles in cytotoxicity, structures of these oligomer intermediates and their conformational conversions remain unknown, hindering our understanding of T2D disease mechanism and therapeutic design targeting these early aggregation species. We further applied large-scale discrete molecule dynamics simulations to investigate the oligomerization of full-length hIAPP, employing multiple molecular systems of increasing number of peptides. We found that the oligomerization process was dynamic, involving frequent inter-oligomeric exchanges. On average, oligomers had more α-helices than β-sheets, consistent with ensemble-based experimental measurements. However, in ~4-6% independent simulations, β-rich oligomers expected as the fibrillization intermediates were observed, especially in the pentamer and hexamer simulations. These β-rich oligomers could adopt β-barrel conformations, recently postulated to be the toxic oligomer species but only observed computationally in the aggregates of short amyloid protein fragments. Free-energy analysis revealed high energies of these β-rich oligomers, supporting the nucleated conformational changes of oligomers in amyloid aggregation. β-barrel oligomers of full-length hIAPP with well-defined three-dimensional structures may play an important pathological role in T2D etiology and may be a therapeutic target for the disease.

Graphene quantum dots against human IAPP aggregation and toxicity in vivo

The development of biocompatible nanomaterials has become a new frontier in the detection, treatment and prevention of human amyloid diseases. Here we demonstrated the use of graphene quantum dots (GQDs) as a potent inhibitor against the in vivo aggregation and toxicity of human islet amyloid polypeptide (IAPP), a hallmark of type 2 diabetes. GQDs initiated contact with IAPP through electrostatic and hydrophobic interactions as well as hydrogen bonding, which subsequently drove the peptide fibrillization off-pathway to eliminate the toxic intermediates. Such interactions, probed in vitro by a thioflavin T kinetic assay, fluorescence quenching, circular dichroism spectroscopy, a cell viability assay and in silico by discrete molecular dynamics simulations, translated to a significant recovery of embryonic zebrafish from the damage elicited by IAPP in vivo, as indicated by improved hatching as well as alleviated reactive oxygen species production, abnormality and mortality of the organism. This study points to the potential of using zero-dimensional nanomaterials for in vivo mitigation of a range of amyloidosis.

Distinct oligomerization and fibrillization dynamics of amyloid core sequences of amyloid-beta and islet amyloid polypeptide

A direct observation of amyloid aggregation from isolated peptides to cross-β fibrils is crucial for understanding the nucleation-dependence process, but the corresponding macroscopic timescales impose a major computational challenge. Using rapid all-atom discrete molecular dynamics simulations, we capture the oligomerization and fibrillization dynamics of the amyloid core sequences of amyloid-β (Aβ) in Alzheimer's disease and islet amyloid polypeptide (IAPP) in type-2 diabetes, namely Aβ16-22 and IAPP22-28. Both peptides and their mixture spontaneously assemble into cross-β aggregates in silico, but follow distinct pathways. Aβ16-22 is highly aggregation-prone with a funneled free energy basin toward multi-layer β-sheet aggregates. IAPP22-28, on the other hand, features the accumulation of unstructured oligomers before the nucleation of β-sheets and growth into double-layer β-sheet aggregates. In the presence of Aβ16-22, the aggregation of IAPP22-28 is promoted by forming co-aggregated multi-layer β-sheets. Our study offers a detailed molecular insight to the long-postulated oligomerization-nucleation process in the amyloid aggregations.

Amino acid substitutions [K16A] and [K28A] distinctly affect amyloid β-protein oligomerization

Amyloid β-protein (A β) assembles into oligomers that play a seminal role in Alzheimer's disease (AD), a leading cause of dementia among the elderly. Despite undisputed importance of A β oligomers, their structure and the basis of their toxicity remain elusive. Previous experimental studies revealed that the [K16A] substitution strongly inhibits toxicity of the two predominant A β alloforms in the brain, A β 40 and A β 42, whereas the [K28A] substitution exerts only a moderate effect. Here, folding and oligomerization of [A16]A β 40, [A28]A β 40, [A16]A β 42, and [A28]A β 42 are examined by discrete molecular dynamics (DMD) combined with a four-bead implicit solvent force field, DMD4B-HYDRA, and compared to A β 40 and A β 42 oligomer formation. Our results show that both substitutions promote A β 40 and A β 42 oligomerization and that structural changes to oligomers are substitution- and alloform-specific. The [K28A] substitution increases solvent-accessible surface area of hydrophobic residues and the intrapeptide N-to-C terminal distance within oligomers more than the [K16A] substitution. The [K16A] substitution decreases the overall β-strand content, whereas the [K28A] substitution exerts only a modest change. Substitution-specific tertiary and quaternary structure changes indicate that the [K16A] substitution induces formation of more compact oligomers than the [K28A] substitution. If the structure-function paradigm applies to A β oligomers, then the observed substitution-specific structural changes in A β 40 and A β 42 oligomers are critical for understanding the structural basis of A β oligomer toxicity and correct identification of therapeutic targets against AD.

Stabilization of native amyloid β-protein oligomers by Copper and Hydrogen peroxide Induced Cross-linking of Unmodified Proteins (CHICUP)

Oligomeric assemblies are postulated to be proximate neurotoxic species in human diseases associated with aberrant protein aggregation. Their heterogeneous and transient nature makes their structural characterization difficult. Size distributions of oligomers of several amyloidogenic proteins, including amyloid β-protein (Aβ) relevant to Alzheimer's disease (AD), have been previously characterized in vitro by photo-induced cross-linking of unmodified proteins (PICUP) followed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). Due to non-physiological conditions associated with the PICUP chemistry, Aβ oligomers cross-linked by PICUP may not be representative of in vivo conditions. Here, we examine an alternative Copper and Hydrogen peroxide Induced Cross-linking of Unmodified Proteins (CHICUP), which utilizes naturally occurring divalent copper ions and hydrogen peroxide and does not require photo activation. Our results demonstrate that CHICUP and PICUP applied to the two predominant Aβ alloforms, Aβ40 and Aβ42, result in similar oligomer size distributions. Thioflavin T fluorescence data and atomic force microscopy images demonstrate that both CHICUP and PICUP stabilize Aβ oligomers and attenuate fibril formation. Relative to noncross-linked peptides, CHICUP-treated Aβ40 and Aβ42 cause prolonged disruption to biomimetic lipid vesicles. CHICUP-stabilized Aβ oligomers link the amyloid cascade, metal, and oxidative stress hypotheses of AD into a more comprehensive understanding of the molecular basis of AD pathology. Because copper and hydrogen peroxide are elevated in the AD brain, CHICUP-stabilized Aβ oligomers are biologically relevant and should be further explored as a new therapeutic target.