The Early Phase of β2m Aggregation: An Integrative Computational Study Framed on the D76N Mutant and the ΔN6 Variant

Dimers of D76N-β2-microglobulin display potent antiamyloid aggregation activity

Self-association of WT β_2-microglobulin (WT-β₂m) into amyloid fibrils is associated with the disorder dialysis related amyloidosis. In the familial variant D76N-β₂m, the single amino acid substitution enhances the aggregation propensity of the protein dramatically and gives rise to a disorder that is independent of renal dysfunction. Numerous biophysical and structural studies on WT- and D76N-β₂m have been performed in order to better understand the structure and dynamics of the native proteins and their different potentials to aggregate into amyloid. However, the structural properties of transient D76N-β₂m oligomers and their role(s) in assembly remained uncharted. Here, we have utilized NMR methods, combined with photo-induced crosslinking, to detect, trap, and structurally characterize transient dimers of D76N-β₂m. We show that the crosslinked D76N-β₂m dimers have different structures from those previously characterized for the on-pathway dimers of ΔN6-β₂m and are unable to assemble into amyloid. Instead, the crosslinked D76N-β₂m dimers are potent inhibitors of amyloid formation, preventing primary nucleation and elongation/secondary nucleation when added in substoichiometric amounts with D76N-β₂m monomers. The results highlight the specificity of early protein-protein interactions in amyloid formation and show how mapping these interfaces can inform new strategies to inhibit amyloid assembly.

Predicting stable binding modes from simulated dimers of the D76N mutant of β 2-microglobulin

•

β2m D76N mutant populates an aggregation-prone monomer (I2) with unstructured termini.

•

MD and MM-PBSA indicate that I2 dimers are stabilized by hydrophobic interactions.

•

The termini regions and BC- and DE-loops are prevalent in the most stable interfaces.

•

The most stable dimer has a limited growth potential without structural rearrangement.

The D76N mutant of the β2m protein is a biologically motivated model system to study protein aggregation. There is strong experimental evidence, supported by molecular simulations, that D76N populates a highly dynamic conformation (which we originally named I2) that exposes aggregation-prone patches as a result of the detachment of the two terminal regions. Here, we use Molecular Dynamics simulations to study the stability of an ensemble of dimers of I2 generated via protein–protein docking. MM-PBSA calculations indicate that within the ensemble of investigated dimers the major contribution to interface stabilization at physiological pH comes from hydrophobic interactions between apolar residues. Our structural analysis also reveals that the interfacial region associated with the most stable binding modes are particularly rich in residues pertaining to both the N- and C-terminus, as well residues from the BC- and DE-loops. On the other hand, the less stable interfaces are stabilized by intermolecular interactions involving residues from the CD- and EF-loops. By focusing on the most stable binding modes, we used a simple geometric rule to propagate the corresponding dimer interfaces. We found that, in the absence of any kind of structural rearrangement occurring at an early stage of the oligomerization pathway, some interfaces drive a self-limited growth process, while others can be propagated indefinitely allowing the formation of long, polymerized chains. In particular, the interfacial region of the most stable binding mode reported here falls in the class of self-limited growth.

The folding space of proteinβ2-microglobulin is modulated by a single disulfide bridge

Protein beta-2-microglobulin (β2m) is classically considered the causative agent of dialysis related amyloidosis, a conformational disorder that affects patients undergoing long-term hemodialysis. The wild type (WT) form, the ΔN6 structural variant, and the D76N mutant have been extensively used as model systems ofβ2m aggregation. In all of them, the native structure is stabilized by a disulfide bridge between the sulphur atoms of the cysteine residues 25 (at B strand) and 80 (at F strand), which has been considered fundamental inβ2m fibrillogenesis. Here, we use extensive discrete molecular dynamics simulations of a full atomistic structure-based model to explore the role of this disulfide bridge as a modulator of the folding space ofβ2m. In particular, by considering different models for the disulfide bridge, we explore the thermodynamics of the folding transition, and the formation of intermediate states that may have the potential to trigger the aggregation cascade. Our results show that the dissulfide bridge affects folding transition and folding thermodynamics of the considered model systems, although to different extents. In particular, when the interaction between the sulphur atoms is stabilized relative to the other intramolecular interactions, or even locked (i.e. permanently established), the WT form populates an intermediate state featuring a well preserved core and two unstructured termini, which was previously detected only for the D76N mutant. The formation of this intermediate state may have important implications in our understanding ofβ2m fibrillogenesis.

When a foreign gene meets its native counterpart: computational biophysics analysis of two PgiC loci in the grass Festuca ovina

Duplicative horizontal gene transfer may bring two previously separated homologous genes together, which may raise questions about the interplay between the gene products. One such gene pair is the “native” PgiC1 and “foreign” PgiC2 in the perennial grass Festuca ovina. Both PgiC1 and PgiC2 encode cytosolic phosphoglucose isomerase, a dimeric enzyme whose proper binding is functionally essential. Here, we use biophysical simulations to explore the inter-monomer binding of the two homodimers and the heterodimer that can be produced by PgiC1 and PgiC2 in F. ovina. Using simulated native-state ensembles, we examine the structural properties and binding tightness of the dimers. In addition, we investigate their ability to withstand dissociation when pulled by a force. Our results suggest that the inter-monomer binding is tighter in the PgiC2 than the PgiC1 homodimer, which could explain the more frequent occurrence of the foreign PgiC2 homodimer in dry habitats. We further find that the PgiC1 and PgiC2 monomers are compatible with heterodimer formation; the computed binding tightness is comparable to that of the PgiC1 homodimer. Enhanced homodimer stability and capability of heterodimer formation with PgiC1 are properties of PgiC2 that may contribute to the retaining of the otherwise redundant PgiC2 gene.

Edge Strand Dissociation and Conformational Changes in Transthyretin under Amyloidogenic Conditions

During amyloidogenesis, proteins undergo conformational changes that allow them to aggregate and assemble into insoluble, fibrillar structures. Soluble oligomers that form during this process typically contain 2-24 monomeric subunits and are cytotoxic. Before the formation of these soluble oligomers, monomeric species first adopt aggregation-competent conformations. Knowledge of the structures of these intermediate states is invaluable to the development of molecular strategies to arrest pathological amyloid aggregation. However, the highly dynamic and interconverting nature of amyloidogenic species limits biophysical characterization of their structures during amyloidogenesis. Here, we use molecular dynamics simulations to probe conformations sampled by monomeric transthyretin under amyloidogenic conditions. We show that certain β-strands in transthyretin tend to unfold and sample nonnative conformations and that the edge strands in one β-sheet (the DAGH sheet) are particularly susceptible to conformational changes in the monomeric state. We also find that changes in the tertiary structure of transthyretin can be associated with disruptions to the secondary structure. We evaluated the conformations produced by molecular dynamics by calculating how well molecular-dynamics-derived structures reproduced NMR-derived interatomic distances. Finally, we leverage our computational results to produce experimentally testable hypotheses that may aid experimental explorations of pathological conformations of transthyretin.

The Early Phase of β2-Microglobulin Aggregation: Perspectives From Molecular Simulations

Protein β2-microglobulin is the causing agent of two amyloidosis, dialysis related amyloidosis (DRA), affecting the bones and cartilages of individuals with chronic renal failure undergoing long-term hemodialysis, and a systemic amyloidosis, found in one French family, which impairs visceral organs. The protein’s small size and its biomedical significance attracted the attention of theoretical scientists, and there are now several studies addressing its aggregation mechanism in the context of molecular simulations. Here, we review the early phase of β2-microglobulin aggregation, by focusing on the identification and structural characterization of monomers with the ability to trigger aggregation, and initial small oligomers (dimers, tetramers, hexamers etc.) formed in the so-called nucleation phase. We focus our analysis on results from molecular simulations and integrate our views with those coming from in vitro experiments to provide a broader perspective of this interesting field of research. We also outline directions for future computer simulation studies.

The role of the IT-state in D76N β2-microglobulin amyloid assembly: A crucial intermediate or an innocuous bystander?

The D76N variant of human β₂-microglobulin (β₂m) is the causative agent of a hereditary amyloid disease. Interestingly, D76N-associated amyloidosis has a distinctive pathology compared with aggregation of WT-β₂m, which occurs in dialysis-related amyloidosis. A folding intermediate of WT-β₂m, known as the I_T-state, which contains a nonnative trans Pro-32, has been shown to be a key precursor of WT-β₂m aggregation in vitro. However, how a single amino acid substitution enhances the rate of aggregation of D76N-β₂m and gives rise to a different amyloid disease remained unclear. Using real-time refolding experiments monitored by CD and NMR, we show that the folding mechanisms of WT- and D76N-β₂m are conserved in that both proteins fold slowly via an I_T-state that has similar structural properties. Surprisingly, however, direct measurement of the equilibrium population of I_T using NMR showed no evidence for an increased population of the I_T-state for D76N-β₂m, ruling out previous models suggesting that this could explain its enhanced aggregation propensity. Producing a kinetically trapped analog of I_T by deleting the N-terminal six amino acids increases the aggregation rate of WT-β₂m but slows aggregation of D76N-β₂m, supporting the view that although the folding mechanisms of the two proteins are conserved, their aggregation mechanisms differ. The results exclude the I_T-state as the origin of the rapid aggregation of D76N-β₂m, suggesting that other nonnative states must cause its high aggregation rate. The results highlight how a single substitution at a solvent-exposed site can affect the mechanism of aggregation and the resulting disease.