A conformational fingerprint for amyloidogenic light chains
eLife Assessment
This study addresses an important and long-standing question regarding the molecular mechanism of protein misfolding in Ig light-chain (LC) amyloidosis (AL), a life-threatening condition. By combining advanced techniques, including small-angle X-ray scattering, molecular dynamics simulations, and hydrogen-deuterium exchange mass spectrometry, the authors provide convincing evidence that the ‘H state’ distinguishes amyloidogenic from non-amyloidogenic LCs. These findings not only offer novel insights into LC structural dynamics but also hold promise for guiding therapeutic strategies in AL and will be of particular interest to structural biologists, biophysicists, and many others working on amyloid diseases.
https://doi.org/10.7554/eLife.102002.3.sa0Important: Findings that have theoretical or practical implications beyond a single subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Convincing: Appropriate and validated methodology in line with current state-of-the-art
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
Both immunoglobulin light-chain (LC) amyloidosis (AL) and multiple myeloma (MM) share the overproduction of a clonal LC. However, while LCs in MM remain soluble in circulation, AL LCs misfold into toxic-soluble species and amyloid fibrils that accumulate in organs, leading to distinct clinical manifestations. The significant sequence variability of LCs has hindered the understanding of the mechanisms driving LC aggregation. Nevertheless, emerging biochemical properties, including dimer stability, conformational dynamics, and proteolysis susceptibility, distinguish AL LCs from those in MM under native conditions. This study aimed to identify a2 conformational fingerprint distinguishing AL from MM LCs. Using small-angle X-ray scattering (SAXS) under native conditions, we analyzed four AL and two MM LCs. We observed that AL LCs exhibited a slightly larger radius of gyration and greater deviations from X-ray crystallography-determined or predicted structures, reflecting enhanced conformational dynamics. SAXS data, integrated with molecular dynamics simulations, revealed a conformational ensemble where LCs adopt multiple states, with variable and constant domains either bent or straight. AL LCs displayed a distinct, low-populated, straight conformation (termed H state), which maximized solvent accessibility at the interface between constant and variable domains. Hydrogen-deuterium exchange mass spectrometry experimentally validated this H state. These findings reconcile diverse experimental observations and provide a precise structural target for future drug design efforts.
Introduction
Immunoglobulin light-chain (LC) amyloidosis (AL) is a systemic disease associated with the overproduction and subsequent amyloid aggregation of patient-specific LCs (Merlini et al., 2018; Blancas-Mejia et al., 2018; Cascino et al., 2022; Poshusta et al., 2009). Such aggregation may take place in one or several organs, the heart and kidneys being the most affected ones (Merlini et al., 2018). AL originates from an abnormal proliferation of a plasma cell clone that results in LCs’ overexpression and oversecretion in the bloodstream (Merlini et al., 2018). LCs belonging to both λ and κ isotypes are associated with AL; however, λ-LCs are greatly over-represented in the repertoire of AL patients. Specifically, AL-causing LCs (AL-LCs) most often belong to a specific subset of λ germlines such as IGLV6 (λ6), IGLV1 (λ1), and IGLV3 (λ3) (Perfetti et al., 2012; Kourelis et al., 2017; Comenzo et al., 2001; Absmeier et al., 2023).
λ-LCs are dimeric in solution, with each subunit characterized by two immunoglobulin domains, a constant domain (CL) with a highly conserved sequence and a variable domain (VL) whose extreme sequence variability is the result of genomic recombination and somatic mutations (Bourne et al., 2002; Chiu et al., 2019; Di Noia and Neuberger, 2007; Oberti et al., 2017). VL domains are generally indicated as the key responsible for LC amyloidogenic behavior. The observation that the fibrillar core in most of the structures of ex vivo AL amyloid fibrils consists of VL residues further strengthens this hypothesis (Radamaker et al., 2021a; Radamaker et al., 2021b; Radamaker et al., 2019; Swuec et al., 2019; Puri et al., 2023). However, in a recent cryo-electron microscopy structure, a stretch of residues belonging to the CL domain is also part of the fibrillar core, and mass spectrometry (MS) analysis of several ex vivo fibrils from different patients indicates that amyloids are composed of several LC proteoforms, including full-length LCs (Schulte et al., 2024; Lavatelli et al., 2020; Mazzini et al., 2022; Dasari et al., 2015).
Interestingly, the overproduction of an LC is a necessary but not sufficient condition for the onset of AL. Indeed, the uncontrolled production of a clonal LC is often associated with multiple myeloma (MM), a blood cancer, but only a subset of MM patients develops AL, thus indicating that specific sequence/biophysical properties determine LC amyloidogenicity and AL onset (Oberti et al., 2017; Rottenaicher et al., 2021; Rennella et al., 2019; Klimtchuk et al., 2023). To date, the extreme sequence variability of AL-LCs has prevented the identification of sequence patterns predictive of LC amyloidogenicity; however, it has been reproducibly reported that several biophysical properties correlate with LC aggregation propensity. AL-LCs display a lower thermodynamic and kinetic fold stability compared with non-amyloidogenic LCs found overexpressed in MM patients (hereafter MM-LCs) (Oberti et al., 2017; Mazzini et al., 2022; Dasari et al., 2015; Rottenaicher et al., 2021; Rennella et al., 2019; Klimtchuk et al., 2023).
Previous work on LCs has indicated how differences in conformational dynamics can play a role in the aggregation properties of AL-LCs (Rottenaicher et al., 2021; Rennella et al., 2019; Klimtchuk et al., 2023; Weber et al., 2018; Sun et al., 2023). Oberti et al. compared multiple λ-LCs obtained from either AL patients or MM patients identifying the susceptibility to proteolysis as the best biophysical parameter distinguishing the two sets (Oberti et al., 2017). Weber et al. showed, using a murine-derived κ-LC, how a modification in the linker region can lead to a greater conformational dynamic, an increased susceptibility to proteolysis, as well as an increased in vitro aggregation propensity (Weber et al., 2018). Additionally, AL-LC flexibility and conformational freedom have also been correlated to the proteotoxicity observed in patients affected by cardiac AL and experimentally verified in human cardiac cells and a Caenorhabditis elegans model (Maritan et al., 2020; Broggini et al., 2023). It is noteworthy that the amyloid LCs analyzed in this study were originally purified from patients with cardiac AL.
Here, building on this previous work as well as on our previous experience on β2-microglobulin, another natively folded amyloidogenic protein (Visconti et al., 2019; Sala et al., 2020; Camilloni et al., 2016; Le Marchand et al., 2018; Achour et al., 2020), we aimed at revealing the structural and dynamic determinants of LC amyloidogenicity. With this goal, we investigated the native solution-state dynamics of multiple λ-LCs by combining MD simulations, small-angle X-ray scattering (SAXS), and hydrogen-deuterium exchange mass spectrometry (HDX-MS). Interestingly, we found a unique conformational fingerprint of amyloidogenic LCs corresponding to a low-populated state (defined as H state hereon) characterized by extended linkers, with an accessible VL-CL interface and possible structural rearrangements in the CL-CL interface.
Results
SAXS suggests differences in the conformational dynamics of amyloidogenic and non-amyloidogenic LC
SAXS data were acquired either in bulk or in line with size-exclusion chromatography (SEC) for a set of six LCs previously described (Table 1 and ‘Materials and methods’). Four of these LCs (referred to as H3, H7, and H18 in Oberti et al., 2017 and AL55 in Swuec et al., 2019) were identified in multiple AL patients (AL-LC), while two (referred to as M7 and M10 in Oberti et al., 2017) were identified in MM patients (MM-LC). These LCs cover multiple germlines, with H18 and M7 belonging to the same germline (Table 1). The sequence identity is the largest for H18 and M7 (91.6%), while it is the lowest for AL55 and M7 (75.2%). A table showing the statistics for all pairwise sequence alignments is provided in Supplementary file 1a. For H3, H7, and M7, a crystal structure was previously determined (Oberti et al., 2017), while for H18, AL55, and M10, we obtained a model using either homology modeling (H18 and AL55) or AlphaFold2 (M10). Qualitatively, the SAXS curves in Figure 1 did not reveal any macroscopic deviation of the solution behavior with respect to the crystal or model conformation. For each LC, we compared the experimental and theoretical curves calculated from the LC structures (Table 1) analyzing the residuals and the associated χ2. The analysis indicated a discrepancy between the model conformation and the data in the case of the AL-LCs, which was instead not observed in the case of the MM-LCs. For AL-LC, residuals deviate from normality in the low q region, suggesting some variability in the global size of the system. Additionally, a weak trend distinguishing AL-LC from MM-LCs could be identified in the radius of gyration (Rg) as derived from the Guinier approximation (Table 1). H3, H7, H18, and AL55 have an Rg that is 0.5–0.8 Å larger than M7 and M10, a small but statistically significant difference (p-value<10–5). Overall, the SAXS measurements point to less compact and more structurally heterogeneous AL-LCs compared with more compact and structurally homogeneous MM-LCs.

Comparison of small-angle X-ray scattering (SAXS) data and single light-chain (LC) structures.
Kratky plots of for amyloidosis (AL) and multiple myeloma (MM) LCs. The experimental (orange) and theoretical (black) curves (top panels) and the associated residuals (bottom panels) indicate that AL-LC solution behavior deviates from reference structures more than MM-LC. SAXS was measured as follows: H3 measured in bulk (Hamburg), 3.4 mg/ml; H7 measured in bulk (Hamburg), 3.4 mg/ml; H18 measured by size-exclusion chromatography (SEC)-SAXS (ESRF [European Synchrotron radiation facility]) with the injection concentration of at 2.8 mg/ml; AL55 measured in bulk (ESRF), 2.6 mg/ml; M7 measured in bulk (Hamburg), 3.6 mg/ml; and M10 measured by SEC-SAXS with the injection concentration of 6.7 mg/ml (ESRF). Theoretical SAXS curves were calculated using crysol (Manalastas-Cantos et al., 2021). Log-log plots are shown in Figure 1—figure supplement 1.
LC systems studied in this work.
The table includes information about the germline, phenotype, method to obtain structure, the SAXS curves, and the radius of gyration derived from the SAXS data for all the model proteins studied in this work.
LC | Germline | Phenotype | Structure | SAXS χ2 q < 0.5 Å (q < 0.3) | Rg (SAXS) (nm) |
---|---|---|---|---|---|
H3 | IGLV1-44*01 | AL | 5MTL | 1.6 (1.9) | 2.57 ± 0.02 |
H7 | IGLV1-51*01 | AL | 5MUH | 2.8 (4.0) | 2.56 ± 0.02 |
H18 | IGLV3-19*01 | AL | Homology | 1.6 (1.9) | 2.56 ± 0.01 |
AL55 | IGLV6-57*02 | AL | Homology | 5.1 (7.8) | 2.58 ± 0.04 |
M7 | IGLV3-19*01 | MM | 5MVG | 1.2 (1.2) | 2.50 ± 0.02 |
M10 | IGLV2-14*03 | MM | AF2 | 1.2 (1.2) | 2.51 ± 0.01 |
-
AL, amyloidosis; LC, light chain; MM, multiple myeloma; SAXS, small-angle X-ray scattering.
MD simulations reveal a conformational fingerprint for amyloidogenic light chains
To investigate the conformational dynamics of the six LCs, we performed metadynamics metainference (M&M) MD simulations employing the SAXS curves (q < 0.3 Å) as restraints (‘Materials and methods’; Bonomi et al., 2016b; Bonomi et al., 2016a; Paissoni et al., 2020; Paissoni and Camilloni, 2021). Metainference is a Bayesian framework that allows the integration of experimental knowledge on-the-fly in MD simulations, improving the latter while accounting for the uncertainty in the data and their interpretation. Metadynamics is an enhanced sampling technique that is able to speed up the sampling of the conformational space of complex systems. The combination of SAXS and MD simulations has been shown to be effective for multidomain proteins, as well as for intrinsically disordered proteins (Saad et al., 2021; Ahmed et al., 2021; Thomasen and Lindorff-Larsen, 2022).
For each LC, we performed two independent M&M simulations coupled by the SAXS restraint, accumulating around 120–180 µs of MD per protein (‘Materials and methods’ and Table 2). The resulting conformational ensembles resulted in a generally improved agreement with the SAXS data employed as restraints (Table 2 and Figure 2A). To investigate the differences in the LC local flexibility, we first analyzed the root mean square fluctuations (RMSFs) for the CL and VL separately, averaging over the chains and the replicates (Figure 2B). The RMSF indicates comparable flexibility in most of the regions, with differences localized in the termini and in some loops. The VLs of amyloidogenic LCs are generally more flexible than the MM ones, but this may be associated with the lengths of their complementarity-determining regions (CDRs). Indeed, M10 has the longest and most flexible CDR1 and the shortest and least flexible CDR3. Unexpectedly, there are some differences also in the CL domains. Here, in Figure 2B, the AL-LCs are always more flexible in at least one region even if these differences are relatively small. Overall, the RMSF does not provide a clear indication to differentiate AL and MM-LCs. To provide a global description of the dynamics of the six LC systems, we then introduced two collective variables, namely the elbow angle, describing the relative orientation of VL and CL dimers, and the distance between the VL and CL dimers center of mass, illustrated in Figure 2C.

Light-chain small-angle X-ray scattering (SAXS)-driven molecular dynamics simulations.
(A) Kratky plots and associated residuals (bottom panels) comparing experimental (orange), and theoretical (black) curves obtained by averaging over the metainference ensemble for H3, H7, H18, AL55, M7, and M10, respectively. Theoretical SAXS curves were calculated using crysol (Manalastas-Cantos et al., 2021). (B) Residue-wise root mean square fluctuations (RMSFs) obtained by averaging the two metainference replicates and the two equivalent domains for the six systems studied. The top panel shows data for the variable domains, while the bottom panel shows data for the constant domain. Residues are reported using Chothia and Lesk numbering (Al-Lazikani et al., 1997). (C) Schematic representation of two global collective variables used to compare the conformational dynamics of the different systems, namely the distance between the center of mass of the variable domain (VL) and constant domain (CL) dimers and the angle describing the bending of the two domain dimers.
Metainference simulations performed in this work for the six systems.
For each metainference simulation, the simulation time per replica with the number of replicas is reported; the χ2 of the resulting conformational ensemble with the experimental SAXS curve, the range q < 0.3 Å is the one used as restraint in the simulation; the average radius of gyration with error estimated by block averaging.
LC code | Simulation | Length per replica (ns) (# replicas) | SAXS χ2 q < 0.5 Å (q < 0.3) | Average Rg (nm) |
---|---|---|---|---|
H3 | M&M 1 | 1530 (60) | 1.2 (1.1) | 2.56 ± 0.02 |
M&M 2 | 1520 (60) | 1.2 (1.1) | 2.56 ± 0.02 | |
H7 | M&M 1 | 1627 (60) | 1.1 (1.2) | 2.54 ± 0.04 |
M&M 2 | 1545 (60) | 1.1 (1.2) | 2.54 ± 0.05 | |
H18 | M&M 1 | 1643 (60) | 1.4 (1.6) | 2.53 ± 0.03 |
M&M 2 | 1529 (60) | 1.4 (1.7) | 2.54 ± 0.03 | |
AL55 | M&M 1 | 1545 (60) | 2.7 (3.0) | 2.59 ± 0.02 |
M&M 2 | 1591 (60) | 2.5 (2.7) | 2.58 ± 0.05 | |
M7 | M&M 1 | 1623 (60) | 1.2 (1.2) | 2.52 ± 0.03 |
M&M 2 | 1530 (60) | 1.1 (1.2) | 2.52 ± 0.03 | |
M10 | M&M 1 | 987 (60) | 1.1 (1.1) | 2.53 ± 0.07 |
M&M 2 | 995 (60) | 1.1 (1.1) | 2.54 ± 0.05 |
-
LC, light chain; M&M, metadynamics metainference; SAXS, small-angle X-ray scattering.
In Figure 3, we report the free energy surfaces (FESs) obtained from the processing of the two replicates of each LC as a function of the elbow angle and the CL-VL distance calculated from their center of mass. The visual inspection of the FES indicates converged simulation: in all cases, the replicates explore a comparable free-energy landscape with comparable features. All six LC FESs share common features: a relatively continuous low free energy region along the diagonal, spanning configurations where the CL and VL are bent and close to each other (state LB), and configurations where the CL and VL domains are straight and at relative distance between 3.4 and 4.1 nm (state LS). A subset of LCs, namely H18, M7, and AL55, display conformations where the domains are straight in line (elbow angle >2.5 rad) and in close vicinity, with a relative distance between the center of mass of <3.4 nm (state G). Of note, H18 and M7 belong to the same germline, letting us speculate that this state G may be germline specific. Most importantly, only the AL-LCs display configurations with CL and VL straight in line but well separated at relative distances >4.1 nm; this state H seems to be a fingerprint specific for AL-LCs. A set of configurations exemplifying the four states is reported in Figure 3. The estimates of the populations for the four states LB, LS, G, and H are reported in Table 3. The quantitative analysis indicates that, within the statistical significance of the simulations, states LB and LS represent in all cases most of the conformational space. In the case of H18, AL55, and M7, the compact state G is also significantly populated (10–34%). The state H, associated with amyloidogenic LCs, is populated between 5 and 10% in H3, H7, H18, and AL55 and <1% in M7 and M10.

Free energy surfaces (FESs) for the six light-chain systems under study by metadynamics metainference molecular dynamics simulations.
For each system, the simulations are performed in duplicate. The x-axis represents the elbow angle, indicating the relative bending of the constant and variable domains (in radians), while the y-axis represents the distance in nm between the center of mass of the constant domain (CL) and variable domain (VL) dimers. The free energy is shown with color and isolines every 2kBT corresponding to 5.16 kJ/mol. On each FES are represented four regions (green, LB state; red, LS state; blue, H state; and black, G state), highlighting the main conformational states. For each region, a representative structure is reported in a rectangle of the same color.
Populations of the four states shown in Figure 3 resulting from the two independent metadynamics metainference simulations performed for each of the six LCs.
The population of the H state, which is supposed to be a fingerprint specific for AL-LCs, is in bold.
% | H3 | H7 | H18 | AL55 | M7 | M10 |
---|---|---|---|---|---|---|
LB | 62.0 ± 0.4 | 72.3 ± 2.5 | 22.9 ± 2.8 | 48.3 ± 0.1 | 46.8 ± 0.1 | 48.4 ± 0.3 |
LS | 33.0 ± 0.2 | 15.2 ± 2.7 | 38.4 ± 2.1 | 32.5 ± 2.0 | 35.0 ± 0.1 | 49.1 ± 1.5 |
G | 0.2 ± 0.1 | 0.8 ± 0.4 | 33.8 ± 3.8 | 10.5 ± 1.4 | 17.6 ± 0.1 | 1.8 ± 0.6 |
H | 4.8 ± 0.5 | 11.7 ± 0.5 | 5.0 ± 1.0 | 8.7 ± 0.5 | 0.6 ± 0.2 | 0.8 ± 0.6 |
-
AL, amyloidosis; LC, light chain.
To identify additional differences between the conformations observed in state H and the rest of the conformational space, we focused our attention on the VL-VL and CL-CL dimerization interfaces. In Figure 4, we showed the free energy as a function of the distance between the CL domains versus the distance between the VL domains for each of the four states for one of the two simulations performed on H3; the same analysis for all other simulations is shown in Figure 4—figure supplements 1–6. From the comparison of the FESs, it is clear that only in the conformations corresponding to the state H do the CL-CL dimers display an alternative configuration. In the case of H3, the CL-CL domains in the H state are characterized by a shift toward the configurations characterized by a larger distance; the same is observed in the case of H18 and AL55, while in the case of H7, the H state is characterized by a smaller distance between the CL domains.

Free energy surfaces for the four substates identified in Figure 3 in the case of the first H3 metainference simulation.
The x-axis shows the distance between the centers of mass of the constant domains, while the y-axis shows the distance between the centers of mass of the variable domains. The free energy is shown with color and isolines every 2kBT corresponding to 5.16 kJ/mol. See also Figure 4—figure supplements 1–6 for the same analysis on the other simulations.
Our conformational ensembles allowed us to hypothesize a conformational fingerprint for AL proteins, namely the presence of a weakly but significantly populated state (H) characterized by a more extended quaternary structure, with VL and CL dimers well separated, and with perturbed CL-CL interfaces.
HDX independently validates the amyloidogenic LC conformational fingerprint
To gain further molecular insight into how the dynamics of the tertiary and quaternary structures can be differentiated in AL- and MM-LCs, HDX-MS was performed on our set of proteins. HDX-MS probes the protein dynamics by monitoring the hydrogen-to-deuterium uptake over time and the obtained data well complement structural, biophysical, and computational data. Four LCs from our set (H3, H7, AL55, and M10) yielded good peptide sequence coverages of 98.6, 92.5, 98.6, and 99.1%, respectively, with a redundancy of >4.0 (Figure 5—figure supplements 1–4 and Supplementary file 1b). H18 and M7 were not included in this analysis due to their poor sequence coverage and were not further investigated. Due to sequence heterogeneity among our proteins, the common peptide analyses were not performed. Therefore, the relative deuterium exchange at different HDX time points of 0.5–240 min with respect to zero exchange time was used to compare the dynamics of individual proteins. The average uptake at different time points for selected regions including residues 34–50 and 152–180 is included in Figure 5—figure supplement 5. The individual peptide mapping at all HDX time points is given in Figure 5—figure supplement 6. The deuterium uptakes at 30 min HDX time showed the most pronounced differences between different proteins, which were chosen to illustrate the key structural features in the main figure panel (Figure 5 and Figure 5—videos 1–4).

Hydrogen-deuterium exchange mass spectrometry analysis.
(A) The top panel represents the simplified presentation of the primary structure of a light chain, including variable domain (VL) and constant domain (CL). The location of β-strands according to Chothia and Lesk (Al-Lazikani et al., 1997). (B) The structural mapping and butterfly plot of relative deuterium uptake (Da) of H3. Chain A of H3 structure is colored in a gradient of blue–white–red for an uptake of 0–30% at an exchange time of 30 min. Chain B is shown in gray (right-hand panel). The butterfly plot showing relative deuterium uptake at all time points from 0.5 to 240 min on a gradient of light to dark blue (left-hand panel). (C–E) are the figures corresponding to H7, AL55, and M10, respectively, with the same color coding as in (B). The overall sequence coverage for all proteins was >90%, with a redundancy of >4.0. The VL-VL domains interface peptides covering amino acid residues 34–50 are labeled in orange, while the CL-CL interface region containing 152–180 amino acids is labeled in magenta. Collectively, they form VL-CL interface, which is important to define the H state.
HDX-MS analysis revealed subtle structural dynamics of the individual proteins. The most significant difference between the AL and MM-LCs is observed for residues 34–50, which are part of both the VL-VL dimerization interface and, more importantly in the context of this work, the CL-VL interface. These residues show significantly higher deuterium uptakes in all H-proteins, with H3 being the highest, implying that AL-LCs dimeric interfaces (VL-VL and CL-VL) are more dynamic and hence significantly more destabilized than in M10 (Figure 5 and Figure 5—figure supplement 5). The highly dynamic VL-VL interface of H3 also correlates well with its open VL-VL interface in a crystal structure (PDB: 8P89), which houses two nanobodies interacting with each VL in a dimeric structure (Broggini et al., 2023). On the other hand, residues 54–70, which are not part of either interface, show a higher deuterium uptake and hence more dynamics in M10 protein, which may be a result of redistribution of dynamics in the regions away from the rigid VL-VL interface to stabilize the overall structure as also observed for other proteins previously (Puri and Hsu, 2022; Ko et al., 2019; Figure 5 and Figure 5—figure supplement 5). In contrast, the VL-CL hinge regions (residues 110–120) show homogeneously high flexibility in all the proteins except H7, indicating their higher accessible surface areas (Figure 5). As expected, the CL domains from residues 170–200 also show a similar pattern of average deuterium uptakes and hence flexibility in AL- and MM-LCs, with a minor difference contributed by the rigid VL-CL interface containing residues 165–180 (Figure 5—figure supplements 1 and 4B). This region shows significantly less deuterium uptake (rigid) in both AL and MM proteins compared with other peptides in the CL domain (Figure 5—figure supplements 1–4B). However, comparing the average uptake for this region (residue 152–180) between AL and MM proteins shows that H3 and AL55 have higher uptake than M10. In contrast, H7 is an exception with the lowest deuterium uptake in this region (Figure 5—figure supplement 5, panel 152–180). These data are particularly interesting in the light of our simulations. The dimeric conformations identified in the H state (Figure 3) are characterized by higher accessibility for the CL-VL interface, which is in agreement with the increased accessibility for the regions 34–50 on the VL and 159–180 in the CL observed in the HDX-MS analysis. Notably, the H state of H7 is the only one in which the CL-CL interface is remarkably compact (see Figure 4—figure supplement 2), consistent with the lower HDX for residues 152–180 observed in H7. Overall, the HDX-MS data provide an independent validation of the H state predicted from our conformational ensembles.
Discussion
Understanding the molecular determinants of AL has been hampered by its high-sequence variability in contrast to its highly conserved 3D structure (Absmeier et al., 2023). In this work, building on our previous studies highlighting susceptibility to proteolysis as a property that can discriminate between AL-LC and MM-LC, as well as the role of conformational dynamics in protein aggregation, we characterized LC conformational dynamics under the assumption that AL-LC proteins, despite their sequence diversity, may share a property that emerges at the level of their dynamics. We combined SAXS measurements with MD simulations under the integrative framework of metainference to generate conformational ensembles representing the native state conformational dynamics of four AL-LC and two MM-LC. While SAXS alone already indicated possible differences, its combination with MD allowed us to observe a possible low-populated state, which we refer to as state H, characterized by well-separated VL and CL dimers and a perturbed CL-CL interface, which is significantly populated in AL-LCs while only marginally populated in MM-LCs. Notably, high-energy states associated with amyloidogenic proteins have been previously identified in the case of SH3 (Neudecker et al., 2012) and β2m (Le Marchand et al., 2018). Here, HDX measurements allowed us to independently validate this state by observing increased accessibility in CL-VL interface regions. Importantly, given the limited size of our LC set, we observed in Peterle et al., 2021 that a comparison of λ-6 LCs by HDX-MS revealed differences between AL and MM LCs in a region overlapping with the one reported here (residues 35–50). Specifically, this region exhibited increased deuterium uptake in AL compared with MM proteins. A T32N substitution was shown to mitigate aggregation propensity and reduce deuterium uptake in the AL protein, while the reverse substitution, N32T, applied to the germline sequence, increased it (Peterle et al., 2021). Finally, the H state observed here is comparable to that previously observed for an aggregation-prone mutation in a monomeric κ LC (Weber et al., 2018).
Having established a conformational fingerprint for AL-LC proteins, it would be tempting to identify possible mutations that could be associated with the presence of the H state. Comparing the sequences and structures of M7 and H18, both of which belong to the IGLV3-19*01 germline, we can identify a single mutation, A40G, that could easily be associated with the appearance of the H state in H18. This mutation is located in the 37–43 loop, which H/D exchange showed to be more accessible in our three AL-LCs than in our MM-LC (see Figure 5 and Figure 5—figure supplement 5), and it breaks a hydrophobic interaction with the methyl group of T165, as observed in the crystal structure of M7 (PDB: 5MVG and Figure 5—figure supplement 6), potentially making T165 more accessible in H18 than in M7 (see Figure 5 and Figure 5—figure supplement 5). Comparing the H18 and M7 sequences with the germline reference sequence, we see that position 40 in IGLV3-19*01 is a glycine (see Figure 6—figure supplement 1). This would suggest the intriguing interpretation that the G40A mutation in M7 may increase the interdomain stability compared with the germline sequence, making it less susceptible to aggregation. However, it should also be noted that, while this framework position is a glycine in H3, H7, and AL55, it is also a glycine in M10. Previous research has often focused on identifying, on a case-by-case basis, the key mutations that may be considered responsible for the emergence of the aggregation propensity under the assumption that such aggregation propensity should not be present in germline sequences, but this assumption may be misleading given the observation that few germlines are strongly over-represented in AL, suggesting that these starting germline sequences may be inherently more aggregation-prone than the germline genes that are absent or rarely found in AL patients. More generally, by comparing our AL-LC sequences with their germline references (Figure 6—figure supplements 2–7), we observe that all mutations fall exclusively in the variable domain, allowing us to exclude for these systems a direct role for residues in the linker region, as observed in Weber et al., 2018 and Nokwe et al., 2015, or in the constant domain, as observed in Rottenaicher et al., 2023. Many mutations fall in the CDR regions, as expected, but others are found in the framework regions, both near the dimerization interface and in other regions of the protein. Somatic mutations in the dimerization interfaces have been identified previously as possibly responsible for toxicity of AL-LC (Garofalo et al., 2021). Regarding mutations in the CDRs, it has been suggested that AL-LC proteins may exhibit frustrated CDR2 and CDR3 loops, with few key residues populating the left-hand alpha helix or other high-energy conformations (Pradhan et al., 2023), resulting in the destabilization of the VL. In Figure 6—figure supplements 2–7, we have analyzed the Ramachandran plot obtained from our conformational ensembles, focusing only on those residues that most populate the left-hand alpha helix region, which are marked with a red circle symbol and whose Ramachandran is reported. Our data indicate the presence of residues populating left-hand alpha regions in the Ramachandran plot, but these are also found in the case of MM-LCs, so our simulations do not allow us to confirm or exclude this mechanism in our set of protein systems.
In conclusion, our study provides a novel, complementary, perspective on the determinants of the misfolding propensity of AL-LCs that we schematize in Figure 6. The identification of a high-energy state, with perturbed CL dimerization interfaces, extended linkers, and accessible regions in both the VL-CL and VL-VL interfaces, may be the common feature interplaying with specific properties shown by previous work, including the direct or indirect destabilization of both the VL-VL and CL-CL dimerization interfaces (Rottenaicher et al., 2021; Rennella et al., 2019; Rottenaicher et al., 2023; Garofalo et al., 2021; Pradhan et al., 2023; Peterson et al., 2010; Kazman et al., 2020). Our conformational fingerprint is also consistent with the observation that protein stability does not fully correlate with the tendency to aggregate, whereas susceptibility to proteolysis and conformational dynamics may better capture the differences between AL-LC and MM-LC. In this context, our data allow us to rationally suggest that targeting the constant domain region at the CL-VL interface, which is more labile in the H state, may be a novel strategy to search for molecules against LC aggregation in AL.

Schematic representation summarizing our findings in the context of previous work on the biophysical properties of amyloidogenic light chains (LCs).
We propose that the H state is the conformational fingerprint distinguishing amyloidosis (AL) LCs from other LCs, which, together with other features, contributes to the amyloidogenicity of AL LCs.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Strain, strain background (Escherichia coli) | DH5α | NEB 5-alpha | Cat# C2987H | Chemical competent cells |
Strain, strain background (E. coli) | BL21(DE3) | NEB | Cat# C2527H | Chemical competent cells |
Peptide, recombinant protein | H3 | Oberti et al., 2017 | NA | |
Peptide, recombinant protein | H7 | Oberti et al., 2017 | NA | |
Peptide, recombinant protein | AL55 | Puri et al., 2025 | NA | |
Peptide, recombinant protein | H18 | Oberti et al., 2017 | NA | |
Peptide, recombinant protein | M7 | Oberti et al., 2017 | NA | |
Peptide, recombinant protein | M10 | Oberti et al., 2017 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-H3 | Oberti et al., 2017 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-H7 | Oberti et al., 2017 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-AL55 | Puri et al., 2025 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-H18 | Oberti et al., 2017 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-M7 | Oberti et al., 2017 | NA | |
Recombinant DNA reagent (plasmid) | pET21b(+)-M10 | Oberti et al., 2017 | NA | |
Software, algorithm | PLGS | Waters | RRID:SCR_016664 | Section ‘Hydrogen-deuterium mass exchange spectrometry’ http://www.waters.com/ |
Software, algorithm | DynamX | Waters | Section ‘Hydrogen-deuterium mass exchange spectrometry’ http://www.waters.com/ | |
Software, algorithm | PyMOL | Schrodinger LLC | RRID:SCR_000305 | Section ‘Hydrogen-deuterium mass exchange spectrometry’ http://www.pymol.org/ |
Software, algorithm | GraphPad Prism | GraphPad | RRID:SCR_002798 | Section ‘Hydrogen-deuterium mass exchange spectrometry’ http://www.graphpad.com/ |
Software, algorithm | ATSAS 3 package | Manalastas-Cantos et al., 2021 | RRID:SCR_015648 | Section ‘Small-angle X-ray scattering’ http://www.embl-hamburg.de/biosaxs/atsas-online/ |
Software, algorithm | GROMACS 2019 | Abraham et al., 2015 | RRID:SCR_014565 | Section ‘Molecular dynamics simulations’ https://www.gromacs.org/ |
Software, algorithm | PLUMED2 v2.9 | Tribello et al., 2014 | RRID:SCR_021952 | Section ‘Molecular dynamics simulations’ https://www.plumed.org/ |
Other | Hi prep QFF- 20 ml column | Cytiva | Product code 28936543 | Section ‘Protein production and purification’ |
Other | Superdex 200 10/300 increase | Cytiva | Cat# 28990946 | Section ‘Protein production and purification’ |
Other | Superdex 75 10/300 GL | Cytiva | Cat# 17517401 | Section ‘Protein production and purification’ |
Other | Immobilized pepsin digestion column | Waters | (Waters Enzymate BEH Pepsin, 2.1 × 30 mm) | Section ‘Hydrogen-deuterium mass exchange spectrometry’ |
Other | Isopropyl-β-d-thiogalactopyranoside (IPTG) | Himedia | RM2578 | Section ‘Protein production and purification’ |
Other | Guanidinium chloride | Sigma-Aldrich | G3272-10KG | Section ‘Protein production and purification’ |
Other | D2O | Sigma-Aldrich | CAS 7789-20-0 | Section ‘Hydrogen-deuterium mass exchange spectrometry’ |
LC production and purification
Request a detailed protocolRecombinant AL- (H3, H7, H18, and AL5) and M- (M7, M10) proteins were produced and purified from the host Escherichia coli strain BL21(DE3). First, the competent BL21(DE3) cells were transformed with plasmid pET21(b+), which contains genes encoding H3, H7, H18, AL55, M7, and M10 proteins. The transformed cells were selected for each plasmid by growing them on LB agar plates containing the antibiotic ampicillin at a final concentration of 100 µg/ml. For overexpression of protein, one colony was picked from each plate and grown overnight in 20 ml of LB broth containing ampicillin at a final concentration of 100 µg/ml. The overnight-grown cells were then used to inoculate a secondary culture in 1 l of LB broth. The cells were grown until the turbidity (OD600nm) reached between 0.6 and 0.8 and protein expression was subsequently induced by adding 0.5 mM isopropyl-β-d-thiogalactopyranoside for 4 h. The bacterial cells containing overexpressed LCs were then harvested using a Backman Coulter centrifuge at 6000 rpm for 20 min at 4°C. All the proteins were overexpressed as inclusion bodies. For protein purification, the inclusion bodies were isolated by cell lysis induced by sonication. The purification of inclusion bodies was performed by washing them with buffer containing 10 mM Tris (pH 8) and 1% triton X 100. The purified inclusion bodies were unfolded with buffer containing 6.0M guanidinium chloride (GdnHCl) for 4 h at 4°C. The unfolded LCs were then refolded in a buffer containing reduced and oxidized glutathione to assist in disulfide bond formation. The refolded proteins were subjected to anion exchange and SEC steps for final purification. The level of protein purity was checked on 12% sodium dodecyl sulfate-polyacrylamide gel electrophoresis gels. The final protein concentration was measured using molecular weight and extinction coefficient of individual proteins. The purified proteins were stored at –20°C for further use. Size exclusion coupled multiangle light scattering confirmed the dimeric assembly of purified proteins (Figure 1—figure supplement 2), which were used for SAXS and HDX-MS experiments mentioned below.
Small-angle X-ray scattering
Request a detailed protocolFor SAXS analysis, H3 was diluted to 3.4 mg/ml, H7 was diluted to 3.4 mg/ml, H18 was diluted to 2.8 mg/ml, AL55 was diluted to 2.6 mg/ml, M7 was diluted to 3.6 mg/ml, in 20 mM Tris-HCl, 150 mM NaCl, pH 8. H3, H7, and M7 batch data were collected at the P12 BioSAXS beamline of the EMBL Hamburg Synchrotron (Blanchet et al., 2015), while AL55 batch data and H18 and M10 SEC data were collected at the BM29 BioSAXS beamline of the ESRF, Grenoble (Pernot et al., 2010). For SEC-SAXS, H18 and M10 were injected into a superdex 200 increase 10/300 GL column previously equilibrated in 20 mM Tris-HCl, 150 mM NaCl, pH 8, at a concentration of 2.8 mg/ml and 6.7 mg/ml, respectively (see also Figure 1—figure supplement 3). SAXS data were processed using programs PRIMUS and GNOM within the ATSAS package (Manalastas-Cantos et al., 2021). Data are deposited in the SASBDB (Valentini et al., 2015) and available with accession codes SASDVL4, SASDVM4, SASDVN4, SASDVK4, SASDVP4, and SASDVQ4.
Molecular dynamics simulations
Request a detailed protocolThe available crystallographic structures of H3, H7, and M7 (PDB: 5mtl, 5muh, and 5mvg, respectively; Oberti et al., 2017) were used as starting conformations, using Modeller to add missing residues (Webb and Sali, 2016). H18 and AL55 were modeled by homology modeling using SwissModel (Waterhouse et al., 2018), while M10 was modeled using AF2 (Jumper et al., 2021). Simulations were performed using GROMACS 2019 (Abraham et al., 2015) and the PLUMED2 software (Tribello et al., 2014; Bonomi and Camilloni, 2017; Bonomi, 2019), using AMBER-DES force field and TIP4P-D water (Piana et al., 2020; Piana et al., 2015). During in-vacuum minimization, RMSD (Root mean square deviation)-restraints were imposed to enhance the symmetry between the two constant and the two variable domains. The systems were solvated in a periodic dodecahedron box, initially 1.2 nm larger than the protein in each direction, neutralized with Na and Cl ions to reach a salt concentration of 10 mM, then minimized and equilibrated at the temperature of 310 K and pressure of 1 atm using the Berendsen thermostat and barostat. Two independent 900-ns-long plain MD simulations were run to generate reliable and independent starting conformations for the M&M simulations (Bonomi et al., 2016b; Bonomi et al., 2016a; Löhr et al., 2017). Also, 30 conformations were extracted from each simulation and duplicated by inverting the two chains to obtain 60 starting conformations symmetrically distributed with respect to chain inversion.
M&M production simulations were run in duplicate using 60 replicas, each replica evolved for ~1 μs (Table 2). Simulations were performed in the NPT ensemble maintaining the temperature at 310 K with the Bussi thermostat (Bussi et al., 2007) and the pressure of 1 atm with the Parrinello–Rahman barostat Parrinello and Rahman, 1981; the electrostatic was treated using the particle mesh Ewald scheme with a short-range cut‐off of 0.9 nm, and van der Waals interaction cut‐off was set to 0.9 nm. To reduce the computational cost, the hydrogen mass repartitioning scheme was used (Hopkins et al., 2015): the mass of heavy atoms was repartitioned into the bonded hydrogen atoms using the heavyh flag in the pdb2gmx tool, and the LINCS algorithm was used to constraint all bonds, allowing to use a time step of 5 fs. In these simulations, Parallel Bias Metadynamics (Pfaendtner and Bonomi, 2015) was used to enhance the sampling, combined with well-tempered metadynamics and the multiple-walker scheme, where Gaussians with an initial height of 1.0 kJ/mol were deposited every 0.5 ps using a bias factor of 10. Five CVs(Collective Variable) were biased, including combinations of phi/psi dihedral angles of the linker regions (i.e., residues connecting variable and constant domains) in the two chains, combinations of chi dihedral angles of the linker regions in the two chains, combination of inter-domain contacts between the variable and the constant domains. The width of the Gaussians was 0.07, 0.12, and 120 for the combination of phi/psi, of chi dihedral angles and combination of contacts, respectively. Metainference was used to include SAXS restraints, using the hySAXS hybrid approach described in Paissoni et al., 2020, Paissoni and Camilloni, 2021, and Ballabio et al., 2023. A set of 13 representative SAXS intensities at different scattering angles, ranging between 0.015 Å−1 and 0.25 Å−1 and equally spaced, was used as restraints. These intensities were extracted from experimental data after performing regularization with the Distance Distribution tool of Primus, based on Gnom (Manalastas-Cantos et al., 2021). Metainference was applied every five steps using a single Gaussian noise per data point and sampling a scaling factor between experimental and calculated SAXS intensities with a flat prior between 0.5 and 1.5. The aggregate sampling from the 60 replicas was reweighted using the final metadynamics bias to obtain a conformational ensemble where each conformation has an associated statistical weight (Branduardi et al., 2012). Convergence and error estimates were assessed by the inspection of the two replicated M&M run. SAXS data were then recalculated using crysol (Manalastas-Cantos et al., 2021). All relevant simulation data are available on Zenodo (cf. Dataset S1 in Supplementary file 1).
Hydrogen-deuterium mass exchange spectrometry
Request a detailed protocolSYNAPT G2-HDMS system (Waters Corporation, USA) equipped with a LEAP robotic liquid handler was used to perform HDX-MS measurements in a fully automated mode as described previously (Puri and Hsu, 2022; Ko et al., 2019; Masson et al., 2019; Puri et al., 2022). The data collection was carried out by a 20-fold dilution of H3, H7, AL55, and M10 proteins (100 µM) with the labeling buffer 1× phosphate-buffered saline prepared in D2O (pD 7.4) to trigger HDX for 0, 0.5, 1,10, 30, 120, and 240 min at 25°C in technical triplicates. Each reaction was quenched by mixing the labeled protein with quench buffer (50 mM sodium phosphate, 250 mM TCEP(tris(2-carboxyethyl)phosphine), 3.0 M GdnHCl [pH 2]) in a 1:1 ratio at 0°C. Online digestion was then performed using an immobilized pepsin digestion column (Waters Enzymate BEH Pepsin, 2.1 × 30 mm). The digested peptides were trapped using a C18 trapping column (Acquity BEH VanGuard 1.7 µm, 2.1 × 5.0 mm) and separated by a linear acetonitrile gradient of 5–40%. Protein Lynx Global Server (PLGS) and DynamX (Waters Corporation) were used to identify the individual peptides, and subsequently, data processing using parameters: maximum peptide length of 25, the minimum intensity of 1000; minimum ion per amino acid of 0.1; maximum MH+ error of 5 ppm; and a file threshold of 3. A reference molecule [(Glu1)-fibrinopeptide B human (CAS# 103213-49-6, Merck, USA)] was used to lock mass with an expected molecular weight of 785.8426 Da. The obtained peptide coverage of H3, H7, AL55, and M10 was 98.6, 92.5, 98.6, and 99.1%, respectively, with a redundancy of >4.0. Backbone amide groups exhibited a relative deuterium uptake (with respect to the zero exchange time data) of 0–30% within 4 h of exchange time. The relative deuterium uptake data at different HDX times were then used to generate heat maps for each amino acid. The obtained data were mapped on individual protein structures in a gradient of blue–white–red showing 0–30% of uptake, respectively (DeLano, 2025). Red denotes the dynamic peptides, while peptide colors in blue are rigid. The common peptide analysis between different model LCs was not performed in this analysis as they all generate different peptides due to both sequence heterogeneity in LCs, especially in the VL domain, and non-specific cleavage of pepsin enzyme used to generate peptides after deuterium exchanged for MS analysis. Therefore, all interpretation was done on individual relative exchange data.
Data availability
Small-angle x-ray scattering data have been deposited in SASBDB under accession codes SASDVL4, SASDVM4, SASDVN4, SASDVK4, SASDVP4, and SASDVQ4. All molecular dynamics simulation trajectories and associated statistical weights are available on Zenodo.
-
ZenodoSimulations data for the paper "A conformational fingerprint for amyloidogenic light chains".https://doi.org/10.5281/zenodo.12731283
References
-
Antibodies gone bad - the molecular mechanism of light chain amyloidosisThe FEBS Journal 290:1398–1419.https://doi.org/10.1111/febs.16390
-
Refinement of α-synuclein ensembles against SAXS data: comparison of force fields and methodsFrontiers in Molecular Biosciences 8:654333.https://doi.org/10.3389/fmolb.2021.654333
-
Standard conformations for the canonical structures of immunoglobulinsJournal of Molecular Biology 273:927–948.https://doi.org/10.1006/jmbi.1997.1354
-
Accurate and efficient SAXS/SANS implementation including solvation layer effects suitable for molecular simulationsJournal of Chemical Theory and Computation 19:8401–8413.https://doi.org/10.1021/acs.jctc.3c00864
-
Immunoglobulin light chain amyloid aggregationChemical Communications 54:10664–10674.https://doi.org/10.1039/c8cc04396e
-
Versatile sample environments and automation for biological solution X-ray scattering experiments at the P12 beamline (PETRA III, DESY)Journal of Applied Crystallography 48:431–443.https://doi.org/10.1107/S160057671500254X
-
Metainference: a Bayesian inference method for heterogeneous systemsScience Advances 2:e1501177.https://doi.org/10.1126/sciadv.1501177
-
Integrative structural and dynamical biology with PLUMED-ISDBBioinformatics 33:3999–4000.https://doi.org/10.1093/bioinformatics/btx529
-
Three-dimensional structure of an immunoglobulin light-chain dimer with amyloidogenic propertiesActa Crystallographica. Section D, Biological Crystallography 58:815–823.https://doi.org/10.1107/s0907444902004183
-
Metadynamics with adaptive gaussiansJournal of Chemical Theory and Computation 8:2247–2254.https://doi.org/10.1021/ct3002464
-
Canonical sampling through velocity rescalingThe Journal of Chemical Physics 126:014101.https://doi.org/10.1063/1.2408420
-
Single-molecule real-time sequencing of the M protein: Toward personalized medicine in monoclonal gammopathiesAmerican Journal of Hematology 97:E389–E392.https://doi.org/10.1002/ajh.26684
-
Proteomic detection of immunoglobulin light chain variable region peptides from amyloidosis patient biopsiesJournal of Proteome Research 14:1957–1967.https://doi.org/10.1021/acs.jproteome.5b00015
-
Molecular mechanisms of antibody somatic hypermutationAnnual Review of Biochemistry 76:1–22.https://doi.org/10.1146/annurev.biochem.76.061705.090740
-
Long-time-step molecular dynamics through hydrogen mass repartitioningJournal of Chemical Theory and Computation 11:1864–1874.https://doi.org/10.1021/ct5010406
-
Mass spectrometry characterization of light chain fragmentation sites in cardiac AL amyloidosis: insights into the timing of proteolysisThe Journal of Biological Chemistry 295:16572–16584.https://doi.org/10.1074/jbc.RA120.013461
-
Metadynamic metainference: convergence towards force field independent structural ensembles of a disordered peptideThe Journal of Chemical Physics 146:165102.https://doi.org/10.1063/1.4981211
-
ATSAS 3.0: expanded functionality and new tools for small-angle scattering data analysisJournal of Applied Crystallography 54:343–355.https://doi.org/10.1107/S1600576720013412
-
Inherent biophysical properties modulate the toxicity of soluble amyloidogenic light chainsJournal of Molecular Biology 432:845–860.https://doi.org/10.1016/j.jmb.2019.12.015
-
Systemic immunoglobulin light chain amyloidosisNature Reviews. Disease Primers 4:38.https://doi.org/10.1038/s41572-018-0034-3
-
The antibody light-chain linker is important for domain stability and amyloid formationJournal of Molecular Biology 427:3572–3586.https://doi.org/10.1016/j.jmb.2015.09.012
-
Determination of protein structural ensembles by hybrid-resolution SAXS restrained molecular dynamicsJournal of Chemical Theory and Computation 16:2825–2834.https://doi.org/10.1021/acs.jctc.9b01181
-
How to determine accurate conformational ensembles by metadynamics metainference: a chignolin study caseFrontiers in Molecular Biosciences 8:694130.https://doi.org/10.3389/fmolb.2021.694130
-
Polymorphic transitions in single crystals: a new molecular dynamics methodJournal of Applied Physics 52:7182–7190.https://doi.org/10.1063/1.328693
-
Efficient sampling of high-dimensional free-energy landscapes with parallel bias metadynamicsJournal of Chemical Theory and Computation 11:5062–5067.https://doi.org/10.1021/acs.jctc.5b00846
-
Water dispersion interactions strongly influence simulated structural properties of disordered protein statesThe Journal of Physical Chemistry. B 119:5113–5123.https://doi.org/10.1021/jp508971m
-
Development of a force field for the simulation of single-chain proteins and protein-protein complexesJournal of Chemical Theory and Computation 16:2494–2507.https://doi.org/10.1021/acs.jctc.9b00251
-
Impacts of cancer-associated mutations on the structure-activity relationship of BAP1Journal of Molecular Biology 434:167553.https://doi.org/10.1016/j.jmb.2022.167553
-
Oxidation of catalytic cysteine of human deubiquitinase BAP1 triggers misfolding and aggregation in addition to functional lossBiochemical and Biophysical Research Communications 599:57–62.https://doi.org/10.1016/j.bbrc.2022.02.029
-
The critical role of the variable domain in driving proteotoxicity and aggregation in full-length light chainsJournal of Molecular Biology 437:168958.https://doi.org/10.1016/j.jmb.2025.168958
-
Molecular mechanism of amyloidogenic mutations in hypervariable regions of antibody light chainsThe Journal of Biological Chemistry 296:100334.https://doi.org/10.1016/j.jbc.2021.100334
-
High conformational flexibility of the E2F1/DP1/DNA ComplexJournal of Molecular Biology 433:167119.https://doi.org/10.1016/j.jmb.2021.167119
-
Role of conformational dynamics in pathogenic protein aggregationCurrent Opinion in Chemical Biology 73:102280.https://doi.org/10.1016/j.cbpa.2023.102280
-
Conformational ensembles of intrinsically disordered proteins and flexible multidomain proteinsBiochemical Society Transactions 50:541–554.https://doi.org/10.1042/BST20210499
-
PLUMED 2: New feathers for an old birdComputer Physics Communications 185:604–613.https://doi.org/10.1016/j.cpc.2013.09.018
-
SASBDB, a repository for biological small-angle scattering dataNucleic Acids Research 43:D357–D363.https://doi.org/10.1093/nar/gku1047
-
Investigating the molecular basis of the aggregation propensity of the pathological D76N mutant of beta-2 microglobulin: role of the denatured stateInternational Journal of Molecular Sciences 20:396.https://doi.org/10.3390/ijms20020396
-
SWISS-MODEL: homology modelling of protein structures and complexesNucleic Acids Research 46:W296–W303.https://doi.org/10.1093/nar/gky427
-
The antibody light-chain linker regulates domain orientation and amyloidogenicityJournal of Molecular Biology 430:4925–4940.https://doi.org/10.1016/j.jmb.2018.10.024
Article and author information
Author details
Funding
Ministero dell'Università e della Ricerca (PRIN 20207XLJB2)
- Stefano Ricagno
Fondazione Cariplo (Telethon GJC23044)
- Stefano Ricagno
Ministero della Salute (#GR-2018-12368387)
- Stefano Ricagno
Fondazione AIRC per la ricerca sul cancro ETS (IG 2024 ID 30307)
- Stefano Ricagno
Academia Sinica (AS-CDA-109- L08)
- Shang-Te Danny Hsu
National Science and Technology Council (113-2123-M-001-010-)
- Shang-Te Danny Hsu
National Science and Technology Council (110-2113-M-001-050-MY3)
- Shang-Te Danny Hsu
National Science and Technology Council (113-2811-M-001-110)
- Manoj K Sriramoju
Academia Sinica (AS-IV-114-L04)
- Shang-Te Danny Hsu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We acknowledge Martin A Schroer and Dmitri Svergun at EMBL Hamburg, and Sonia Longhi at AFMB Marseille, for discussion and support on the SAXS data acquisition and analysis. We acknowledge PRACE for awarding us access to Piz Daint at CSCS, Switzerland. We acknowledge CINECA for an award under the ISCRA initiative for the availability of high-performance computing resources and support. SP acknowledges Fondazione Veronesi for a postdoctoral fellowship. We also acknowledge Dr. Min-Feng Karen Hsu and Mr. Yong-Sheng Wang for the protein quality controls, and Dr. Shu-Yu Lin and Mr. Ming-Jie Tsai of the Academia Sinica Common Mass Spectrometry Facilities (AS-CFII-111-209), funded by the Academia Sinica Core Facility and Innovative Instrument Project, for supporting the HDX-MS experiments.
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.102002. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Paissoni, Puri et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 419
- views
-
- 11
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Neuroscience
The buildup of knot-like RNA structures in brain cells may be the key to understanding how uncontrolled protein aggregation drives Alzheimer’s disease.
-
- Biochemistry and Chemical Biology
- Neuroscience
As the world population ages, new molecular targets in aging and Alzheimer’s disease (AD) are needed to combat the expected influx of new AD cases. Until now, the role of RNA structure in aging and neurodegeneration has largely remained unexplored. In this study, we examined human hippocampal postmortem tissue for the formation of RNA G-quadruplexes (rG4s) in aging and AD. We found that rG4 immunostaining strongly increased in the hippocampus with both age and with AD severity. We further found that neurons with the accumulation of phospho-tau immunostaining contained rG4s, rG4 structure can drive tau aggregation, and rG4 staining density depended on APOE genotype in the human tissue examined. Combined with previous studies showing the dependence of rG4 structure on stress and the extreme power of rG4s at oligomerizing proteins, we propose a model of neurodegeneration in which chronic rG4 formation is linked to proteostasis collapse. These morphological findings suggest that further investigation of RNA structure in neurodegeneration is a critical avenue for future treatments and diagnoses.