Truncation of the constant domain drives amyloid formation by immunoglobulin light chains.

AL amyloidosis is a life-threatening disease caused by deposition of immunoglobulin light chains. While the mechanisms underlying light chains amyloidogenesis in vivo remain unclear, several studies have highlighted the role that tissue environment and structural amyloidogenicity of individual light chains have in the disease pathogenesis. AL natural deposits contain both full-length light chains and fragments encompassing the variable domain (VL) as well as different length segments of the constant region (CL), thus highlighting the relevance that proteolysis may have in the fibrillogenesis pathway. Here, we investigate the role of major truncated species of the disease-associated AL55 light chain that were previously identified in natural deposits. Specifically, we study structure, molecular dynamics, thermal stability, and capacity to form fibrils of a fragment containing both the VL and part of the CL (133-AL55), in comparison with the full-length protein and its variable domain alone, under shear stress and physiological conditions. Whereas the full-length light chain forms exclusively amorphous aggregates, both fragments generate fibrils, although, with different kinetics, aggregate structure, and interplay with the unfragmented protein. More specifically, the VL-CL 133-AL55 fragment entirely converts into amyloid fibrils microscopically and spectroscopically similar to their ex vivo counterpart and increases the amorphous aggregation of full-length AL55. Overall, our data support the idea that light chain structure and proteolysis are both relevant for amyloidogenesis in vivo and provide a novel biocompatible model of light chain fibrillogenesis suitable for future mechanistic studies.

Clinical ApoA-IV amyloid is associated with fibrillogenic signal sequence.

Apolipoprotein A-IV amyloidosis is an uncommon form of the disease normally resulting in renal and cardiac dysfunction. ApoA-IV amyloidosis was identified in 16 patients attending the National Amyloidosis Centre and in eight clinical samples received for histology review. Unexpectedly, proteomics identified the presence of ApoA-IV signal sequence residues (p.18-43 to p.20-43) in 16/24 trypsin-digested amyloid deposits but in only 1/266 non-ApoA-IV amyloid samples examined. These additional signal residues were also detected in the cardiac sample from the Swedish patient in which ApoA-IV amyloid was first described, and in plasma from a single cardiac ApoA-IV amyloidosis patient. The most common signal-containing peptide observed in ApoA-IV amyloid, p.20-43, and to a far lesser extent the N-terminal peptide, p.21-43, were fibrillogenic in vitro at physiological pH, generating Congo red-positive fibrils. The addition of a single signal-derived alanine residue to the N-terminus has resulted in markedly increased fibrillogenesis. If this effect translates to the mature circulating protein in vivo, then the presence of signal may result in preferential deposition as amyloid, perhaps acting as seed for the main circulating native form of the protein; it may also influence other ApoA-IV-associated pathologies. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd. on behalf of The Pathological Society of Great Britain and Ireland.

Conversion of the Native N-Terminal Domain of TDP-43 into a Monomeric Alternative Fold with Lower Aggregation Propensity.

TAR DNA-binding protein 43 (TDP-43) forms intraneuronal cytoplasmic inclusions associated with amyotrophic lateral sclerosis and ubiquitin-positive frontotemporal lobar degeneration. Its N-terminal domain (NTD) can dimerise/oligomerise with the head-to-tail arrangement, which is essential for function but also favours liquid-liquid phase separation and inclusion formation of full-length TDP-43. Using various biophysical approaches, we identified an alternative conformational state of NTD in the presence of Sulfobetaine 3-10 (SB3-10), with higher content of α-helical structure and tryptophan solvent exposure. NMR shows a highly mobile structure, with partially folded regions and ß-sheet content decrease, with a concomitant increase of α-helical structure. It is monomeric and reverts to native oligomeric NTD upon SB3-10 dilution. The equilibrium GdnHCl-induced denaturation shows a cooperative folding and a somewhat lower conformational stability. When the aggregation processes were compared with and without pre-incubation with SB3-10, but at the identical final SB3-10 concentration, a slower aggregation was found in the former case, despite the reversible attainment of the native conformation in both cases. This was attributed to protein monomerization and oligomeric seeds disruption by the conditions promoting the alternative conformation. Overall, the results show a high plasticity of TDP-43 NTD and identify strategies to monomerise TDP-43 NTD for methodological and biomedical applications.

Comparative study of the stabilities of synthetic in vitro and natural ex vivo transthyretin amyloid fibrils.

Systemic amyloidosis caused by extracellular deposition of insoluble fibrils derived from the pathological aggregation of circulating proteins, such as transthyretin, is a severe and usually fatal condition. Elucidation of the molecular pathogenic mechanism of the disease and discovery of effective therapies still represents a challenging medical issue. The in vitro preparation of amyloid fibrils that exhibit structural and biochemical properties closely similar to those of natural fibrils is central to improving our understanding of the biophysical basis of amyloid formation in vivo and may offer an important tool for drug discovery. Here, we compared the morphology and thermodynamic stability of natural transthyretin fibrils with those of fibrils generated in vitro either using the common acidification procedure or primed by limited selective cleavage by plasmin. The free energies for fibril formation were -12.36, -8.10, and -10.61 kcal mol-1, respectively. The fibrils generated via plasmin cleavage were more stable than those prepared at low pH and were thermodynamically and morphologically similar to natural fibrils extracted from human amyloidotic tissue. Determination of thermodynamic stability is an important tool that is complementary to other methods of structural comparison between ex vivo fibrils and fibrils generated in vitro Our finding that fibrils created via an in vitro amyloidogenic pathway are structurally similar to ex vivo human amyloid fibrils does not necessarily establish that the fibrillogenic pathway is the same for both, but it narrows the current knowledge gap between in vitro models and in vivo pathophysiology.

Insights on peptide topology in the computational design of protein ligands: the example of lysozyme binding peptides.

Herein, we compared the ability of linear and cyclic peptides generated in silico to target different protein sites: internal pockets and solvent-exposed sites. We selected human lysozyme (HuL) as a model target protein combined with the computational evolution of linear and cyclic peptides. The sequence evolution of these peptides was based on the PARCE algorithm. The generated peptides were screened based on their aqueous solubility and HuL binding affinity. The latter was evaluated by means of scoring functions and atomistic molecular dynamics (MD) trajectories in water, which allowed prediction of the structural features of the protein-peptide complexes. The computational results demonstrated that cyclic peptides constitute the optimal choice for solvent exposed sites, while both linear and cyclic peptides are capable of targeting the HuL pocket effectively. The most promising binders found in silico were investigated experimentally by surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and electrospray ionization mass spectrometry (ESI-MS) techniques. All tested peptides displayed dissociation constants in the micromolar range, as assessed by SPR; however, both NMR and ESI-MS suggested multiple binding modes, at least for the pocket binding peptides. A detailed NMR analysis confirmed that both linear and cyclic pocket peptides correctly target the binding site they were designed for.

Plasminogen activation triggers transthyretin amyloidogenesis in vitro.

Systemic amyloidosis is a usually fatal disease caused by extracellular accumulation of abnormal protein fibers, amyloid fibrils, derived by misfolding and aggregation of soluble globular plasma protein precursors. Both WT and genetic variants of the normal plasma protein transthyretin (TTR) form amyloid, but neither the misfolding leading to fibrillogenesis nor the anatomical localization of TTR amyloid deposition are understood. We have previously shown that, under physiological conditions, trypsin cleaves human TTR in a mechano-enzymatic mechanism that generates abundant amyloid fibrils in vitro In sharp contrast, the widely used in vitro model of denaturation and aggregation of TTR by prolonged exposure to pH 4.0 yields almost no clearly defined amyloid fibrils. However, the exclusive duodenal location of trypsin means that this enzyme cannot contribute to systemic extracellular TTR amyloid deposition in vivo Here, we therefore conducted a bioinformatics search for systemically active tryptic proteases with appropriate tissue distribution, which unexpectedly identified plasmin as the leading candidate. We confirmed that plasmin, just as trypsin, selectively cleaves human TTR between residues 48 and 49 under physiological conditions in vitro Truncated and full-length protomers are then released from the native homotetramer and rapidly aggregate into abundant fibrils indistinguishable from ex vivo TTR amyloid. Our findings suggest that physiological fibrinolysis is likely to play a critical role in TTR amyloid formation in vivo Identification of this surprising intersection between two hitherto unrelated pathways opens new avenues for elucidating the mechanisms of TTR amyloidosis, for seeking susceptibility risk factors, and for therapeutic innovation.

PDB2ENTROPY and PDB2TRENT: Conformational and Translational-Rotational Entropy from Molecular Ensembles.

Entropy calculation is an important step in the postprocessing of molecular dynamics trajectories or predictive models. In recent years the nearest neighbor method has emerged as a powerful method to deal in a flexible way with the dimensionality of the problem. Here we provide two programs, PBD2ENTROPY and PDB2TRENT that compute the conformational and translational-rotational entropy, respectively, based on the nearest neighbor method. PDB2ENTROPY takes in input two files containing the following: (1) conformational ensembles of the same molecule(s) in PDB format and (2) definitions of torsion angles (a default file is provided where additional user definitions can be easily implemented). PDB2TRENT takes in a file containing samples of the complexed molecules, a string specifying atoms providing the reference framework to superimpose samples, and a string specifying atoms used to compute rotation and translation of one molecule with respect to the other. The C programs and sample demonstration data are available on the GitHub repository (URL: http://github.com/federico-fogolari/pdb2entropy and http://github.com/federico-fogolari/pdb2trent ).

Automation of peak-tracking analysis of stepwise perturbed NMR spectra.

We describe a new algorithmic approach able to automatically pick and track the NMR resonances of a large number of 2D NMR spectra acquired during a stepwise variation of a physical parameter. The method has been named Trace in Track (TINT), referring to the idea that a gaussian decomposition traces peaks within the tracks recognised through 3D mathematical morphology. It is capable of determining the evolution of the chemical shifts, intensity and linewidths of each tracked peak.The performances obtained in term of track reconstruction and correct assignment on realistic synthetic spectra were high above 90% when a noise level similar to that of experimental data were considered. TINT was applied successfully to several protein systems during a temperature ramp in isotope exchange experiments. A comparison with a state-of-the-art algorithm showed promising results for great numbers of spectra and low signal to noise ratios, when the graduality of the perturbation is appropriate. TINT can be applied to different kinds of high throughput chemical shift mapping experiments, with quasi-continuous variations, in which a quantitative automated recognition is crucial.

Computational design of cyclic peptides for the customized oriented immobilization of globular proteins.

The oriented immobilization of proteins, key for the development of novel responsive biomaterials, relies on the availability of effective probes. These are generally provided by standard approaches based on in vivo maturation and in vitro selection of antibodies and/or aptamers. These techniques can suffer technical problems when a non-immunogenic epitope needs to be targeted. Here we propose a strategy to circumvent this issue by in silico design. In our method molecular binders, in the form of cyclic peptides, are computationally evolved by stochastically exploring their sequence and structure space to identify high-affinity peptides for a chosen epitope of a target globular protein: here a solvent-exposed site of ß2-microglobulin (ß2m). Designed sequences were screened by explicit solvent molecular dynamics simulations (MD) followed by experimental validation. Five candidates gave dose-response surface plasmon resonance signals with dissociation constants in the micromolar range. One of them was further analyzed by means of isothermal titration calorimetry, nuclear magnetic resonance, and 250 ns of MD. Atomic-force microscopy imaging showed that this peptide is able to immobilize ß2m on a gold surface. In short, we have shown by a variety of experimental techniques that it is possible to capture a protein through an epitope of choice by computational design.

Accuracy assessment of the linear Poisson-Boltzmann equation and reparametrization of the OBC generalized Born model for nucleic acids and nucleic acid-protein complexes.

The generalized Born model in the Onufriev, Bashford, and Case (Onufriev et al., Proteins: Struct Funct Genet 2004, 55, 383) implementation has emerged as one of the best compromises between accuracy and speed of computation. For simulations of nucleic acids, however, a number of issues should be addressed: (1) the generalized Born model is based on a linear model and the linearization of the reference Poisson-Boltmann equation may be questioned for highly charged systems as nucleic acids; (2) although much attention has been given to potentials, solvation forces could be much less sensitive to linearization than the potentials; and (3) the accuracy of the Onufriev-Bashford-Case (OBC) model for nucleic acids depends on fine tuning of parameters. Here, we show that the linearization of the Poisson Boltzmann equation has mild effects on computed forces, and that with optimal choice of the OBC model parameters, solvation forces, essential for molecular dynamics simulations, agree well with those computed using the reference Poisson-Boltzmann model.

Monitoring the interaction between ß2-microglobulin and the molecular chaperone αB-crystallin by NMR and mass spectrometry: αB-crystallin dissociates ß2-microglobulin oligomers.

The interaction at neutral pH between wild-type and a variant form (R3A) of the amyloid fibril-forming protein ß2-microglobulin (ß2m) and the molecular chaperone αB-crystallin was investigated by thioflavin T fluorescence, NMR spectroscopy, and mass spectrometry. Fibril formation of R3Aß2m was potently prevented by αB-crystallin. αB-crystallin also prevented the unfolding and nonfibrillar aggregation of R3Aß2m. From analysis of the NMR spectra collected at various R3Aß2m to αB-crystallin molar subunit ratios, it is concluded that the structured ß-sheet core and the apical loops of R3Aß2m interact in a nonspecific manner with the αB-crystallin. Complementary information was derived from NMR diffusion coefficient measurements of wild-type ß2m at a 100-fold concentration excess with respect to αB-crystallin. Mass spectrometry acquired in the native state showed that the onset of wild-type ß2m oligomerization was effectively reduced by αB-crystallin. Furthermore, and most importantly, αB-crystallin reversibly dissociated ß2m oligomers formed spontaneously in aged samples. These results, coupled with our previous studies, highlight the potent effectiveness of αB-crystallin in preventing ß2m aggregation at the various stages of its aggregation pathway. Our findings are highly relevant to the emerging view that molecular chaperone action is intimately involved in the prevention of in vivo amyloid fibril formation.

Structure, folding dynamics, and amyloidogenesis of D76N ß2-microglobulin: roles of shear flow, hydrophobic surfaces, and α-crystallin.

Systemic amyloidosis is a fatal disease caused by misfolding of native globular proteins, which then aggregate extracellularly as insoluble fibrils, damaging the structure and function of affected organs. The formation of amyloid fibrils in vivo is poorly understood. We recently identified the first naturally occurring structural variant, D76N, of human ß2-microglobulin (ß2m), the ubiquitous light chain of class I major histocompatibility antigens, as the amyloid fibril protein in a family with a new phenotype of late onset fatal hereditary systemic amyloidosis. Here we show that, uniquely, D76N ß2m readily forms amyloid fibrils in vitro under physiological extracellular conditions. The globular native fold transition to the fibrillar state is primed by exposure to a hydrophobic-hydrophilic interface under physiological intensity shear flow. Wild type ß2m is recruited by the variant into amyloid fibrils in vitro but is absent from amyloid deposited in vivo. This may be because, as we show here, such recruitment is inhibited by chaperone activity. Our results suggest general mechanistic principles of in vivo amyloid fibrillogenesis by globular proteins, a previously obscure process. Elucidation of this crucial causative event in clinical amyloidosis should also help to explain the hitherto mysterious timing and location of amyloid deposition.

Degradation versus fibrillogenesis, two alternative pathways modulated by seeds and glycosaminoglycans.

The mechanism that converts native human transthyretin into amyloid fibrils in vivo is still a debated and controversial issue. Commonly, non-physiological conditions of pH, temperature, or organic solvents are used in in vitro models of fibrillogenesis of globular proteins. Transthyretin amyloid formation can be achieved under physiological conditions through a mechano-enzymatic mechanism involving specific serine proteases such as trypsin or plasmin. Here, we investigate S52P and L111M transthyretin variants, both causing a severe form of systemic amyloidosis mostly targeting the heart at a relatively young age with heterogeneous phenotype among patients. Our studies on thermodynamics show that both proteins are significantly less stable than other amyloidogenic variants. However, despite a similar thermodynamic stability, L111M variant seems to have enhanced susceptibility to cleavage and a lower tendency to form fibrils than S52P in the presence of specific proteases and biomechanical forces. Heparin strongly enhances the fibrillogenic capacity of L111M transthyretin, but has no effect on the S52P variant. Fibrillar seeds similarly affect the fibrillogenesis of both proteins, with a stronger effect on the L111M variant. According to our model of mechano-enzymatic fibrillogenesis, both full-length and truncated monomers, released after the first cleavage, can enter into fibrillogenesis or degradation pathways. Our findings show that the kinetics of the two processes can be affected by several factors, such as intrinsic amyloidogenicity due to the specific mutations, environmental factors including heparin and fibrillar seeds that significantly accelerate the fibrillogenic pathway.

Single-shot NMR measurement of protein unfolding landscapes.

The transient unfolding events from the native state of a protein towards higher energy states can be closely investigated by studying the process of hydrogen exchange. Here, we present BLUU-Tramp (Biophysics Laboratory University of Udine-Temperature ramp), a new method to measure the rates for the exchange process and the underlying equilibrium thermodynamic parameters, using just a single sample preparation, in a single experiment that lasts some 20 to 60h depending on the protein thermal stability, to record hundreds of points over a virtually continuous temperature window. The method is suitable also in presence of other proteins in the sample, if only the target protein is (15)N-labelled. This allows the complete thermodynamic description of the unfolding landscape at an atomic level in the presence of small or macromolecular ligands or cosolutes, or in physiological environments. The method was successfully tested with human ubiquitin. Then the unfolding thermodynamic parameters were satisfactorily determined for the amyloidogenic protein ß(2)-microglobulin, in aqueous buffer and in synovial liquid, that is the natural medium of amyloid deposition in joints.

Carnosine inhibits Aß(42) aggregation by perturbing the H-bond network in and around the central hydrophobic cluster.

Aggregation of the amyloid-ß peptide (Aß) into fibrillar structures is a hallmark of Alzheimer's disease. Thus, preventing self-assembly of the Aß peptide is an attractive therapeutic strategy. Here, we used experimental techniques and atomistic simulations to investigate the influence of carnosine, a dipeptide naturally occurring in the brain, on Aß aggregation. Scanning force microscopy, circular dichroism and thioflavin T fluorescence experiments showed that carnosine does not modify the conformational features of Aß42 but nonetheless inhibits amyloid growth. Molecular dynamics (MD) simulations indicated that carnosine interacts transiently with monomeric Aß42 by salt bridges with charged side chains, and van der Waals contacts with residues in and around the central hydrophobic cluster ((17)LVFFA(21)). NMR experiments on the nonaggregative fragment Aß12-28 did not evidence specific intermolecular interactions between the peptide and carnosine, in agreement with MD simulations. However, a close inspection of the spectra revealed that carnosine interferes with the local propensity of the peptide to form backbone hydrogen bonds close to the central hydrophobic cluster (residues E22, S26 and N27). Finally, MD simulations of aggregation-prone Aß heptapeptide segments show that carnosine reduces the propensity to form intermolecular backbone hydrogen bonds in the region 18-24. Taken together, the experimental and simulation results (cumulative MD sampling of 0.2 ms) suggest that, despite the inability of carnosine to form stable contacts with Aß, it might block the pathway toward toxic aggregates by perturbing the hydrogen bond network near residues with key roles in fibrillogenesis.

Bluues server: electrostatic properties of wild-type and mutated protein structures.

MOTIVATION: Electrostatic calculations are an important tool for deciphering many functional mechanisms in proteins. Generalized Born (GB) models offer a fast and convenient computational approximation over other implicit solvent-based electrostatic models. Here we present a novel GB-based web server, using the program Bluues, to calculate numerous electrostatic features including pKa-values and surface potentials. The output is organized allowing both experts and beginners to rapidly sift the data. A novel feature of the Bluues server is that it explicitly allows to find electrostatic differences between wild-type and mutant structures. AVAILABILITY: The Bluues server, examples and extensive help files are available for non-commercial use at URL: http://protein.bio.unipd.it/bluues/.

A differential equation for the Generalized Born radii.

The Generalized Born (GB) model offers a convenient way of representing electrostatics in complex macromolecules like proteins or nucleic acids. The computation of atomic GB radii is currently performed by different non-local approaches involving volume or surface integrals. Here we obtain a non-linear second-order partial differential equation for the Generalized Born radius, which may be solved using local iterative algorithms. The equation is derived under the assumption that the usual GB approximation to the reaction field obeys Laplace's equation. The equation admits as particular solutions the correct GB radii for the sphere and the plane. The tests performed on a set of 55 different proteins show an overall agreement with other reference GB models and "perfect" Poisson-Boltzmann based values.

Generalized Born forces: surface integral formulation.

Generalized Born (GB) models offer a convenient alternative to Poisson-Boltzmann based models. In the last decade, the GB radii computed based on the exact results obtained for a charge embedded in a conducting sphere have proven to be accurate also for the complex molecular shapes of proteins. The surface integral formulation of the theory has been much less explored than the volume integral formulation. In this work, we provide the exact equations for the GB solvation forces in the surface integral formulation, which are non-trivial due to the non-negligible dependence of GB radii on atomic positions and due to the discontinuity in the derivative of the solvent accessible surface point positions with respect to atomic positions. The equations derived here provide a useful reference for developing faster approximations.

Pathological self-aggregation of ß(2)-microglobulin: a challenge for protein biophysics.

The pathological aggregation of b(2)-microglobulin (b2m) is examined starting from the relevance of some structural aspects of the protein. The systemic deposition of b2m fibrils has been ascribed to several factors, but no conclusive evidence emerged so far. The characterization of b2m aggregates by direct investigation through electron microscopy, atomic force microscopy, solid state NMR and other solid state techniques provides important structural and morphological information on the assembly, but no clues about the mechanism of the aggregation process. The most relevant mechanistic hypotheses are critically reviewed. In addition to the mechanisms exclusively based on structural features, also the recently reported prion-like conversion is analyzed and shown to hardly comply with some established conditions of the fibrillogenic process. An alternative mechanism is recalled that does not require rare events and involves only the full-length protein in proximity of collagen, i.e. the environment that physiologically supports deposition.

Reduction of conformational mobility and aggregation in W60G ß2-microglobulin: assessment by 15N NMR relaxation.

The amyloid pathology associated with long-term haemodialysis is due to the deposition of ß2-microglobulin, the non-polymorphic light chain of class I major histocompatibility complex, that accumulates at bone joints into amyloid fibrils. Several lines of evidence show the relevance of the tryptophan residue at position 60 for the fibrillogenic transition of the protein. A comparative (15)N NMR relaxation analysis is presented for wild-type human ß2-microglobulin and W60G ß2-microglobulin, i.e. the mutant with a glycyne replacing the natural tryptophan residue at position 60. The experimental data, collected at 11.4 T and 310 K, were analyzed by means of the reduced spectral density approach. Molecular dynamics (MD) simulations and corresponding thermodynamic integration, together with hydrodynamic calculations were performed to support data interpretation. The analysis results for the mutant protein are consistent with a reduced aggregation with respect to the wild-type counterpart, as a consequence of an increased conformational rigidity probed by either NMR relaxation and MD simulations. Although dynamics in solution is other than fibrillar competence, the assessed properties of the mutant protein can be related with its reduced ability of forming fibrils when seeded in 20% trifluoroethanol.