For antibody sequence generative modeling, mixture models may be all you need.

MOTIVATION: Antibody therapeutic candidates must exhibit not only tight binding to their target but also good developability properties, especially low risk of immunogenicity. RESULTS: In this work, we fit a simple generative model, SAM, to sixty million human heavy and seventy million human light chains. We show that the probability of a sequence calculated by the model distinguishes human sequences from other species with the same or better accuracy on a variety of benchmark datasets containing >400 million sequences than any other model in the literature, outperforming large language models (LLMs) by large margins. SAM can humanize sequences, generate new sequences, and score sequences for humanness. It is both fast and fully interpretable. Our results highlight the importance of using simple models as baselines for protein engineering tasks. We additionally introduce a new tool for numbering antibody sequences which is orders of magnitude faster than existing tools in the literature. AVAILABILITY AND IMPLEMENTATION: All tools developed in this study are available at https://github.com/Wang-lab-UCSD/AntPack.

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.

Kinetic evidence for multiple aggregation pathways in antibody light chain variable domains.

Aggregation of antibody light chain proteins is associated with the progressive disease light chain amyloidosis. Patient-derived amyloid fibrils are formed from light chain variable domain residues in non-native conformations, highlighting a requirement that light chains unfold from their native structures in order to aggregate. However, mechanistic studies of amyloid formation have primarily focused on the self-assembly of natively unstructured peptides, and the role of native state unfolding is less well understood. Using a well-studied light chain variable domain protein known as WIL, which readily aggregates in vitro under conditions where the native state predominates, we asked how the protein concentration and addition of pre-formed fibril "seeds" alter the kinetics of aggregation. Monitoring aggregation with thioflavin T fluorescence revealed a distinctly non-linear dependence on concentration, with a maximum aggregation rate observed at 8 µM protein. This behavior is consistent with formation of alternate aggregate structures in the early phases of amyloid formation. Addition of N- or C-terminal peptide tags, which did not greatly affect the folding or stability of the protein, altered the concentration dependence of aggregation. Aggregation rates increased in the presence of pre-formed seeds, but this effect did not eliminate the delay before aggregation and became saturated when the proportion of seeds added was greater than 1 in 1600. The complexity of aggregation observed in vitro highlights how multiple species may contribute to amyloid pathology in patients.

abYpap: improvements to the prediction of antibody VH/VL packing using gradient boosted regression.

The Fv region of the antibody (comprising VH and VL domains) is the area responsible for target binding and thus the antibody's specificity. The orientation, or packing, of these two domains relative to each other influences the topography of the Fv region, and therefore can influence the antibody's binding affinity. We present abYpap, an improved method for predicting the packing angle between the VH and VL domains. With the large data set now available, we were able to expand greatly the number of features that could be used compared with our previous work. The machine-learning model was tuned for improved performance using 37 selected residues (previously 13) and also by including the lengths of the most variable 'complementarity determining regions' (CDR-L1, CDR-L2 and CDR-H3). Our method shows large improvements from the previous version, and also against other modeling approaches, when predicting the packing angle.

Nanobodies counteract the toxicity of an amyloidogenic light chain by stabilizing a partially open dimeric conformation.

Light chain amyloidosis (AL) is a systemic disease where fibrillar deposition of misfolded immunoglobulin light chains (LCs) severely affects organ function and results in poor prognosis for patients, especially when heart involvement is severe. Particularly relevant in this context is the cardiotoxicity exerted by still uncharacterized soluble LC species. Here, with the final goal of identifying alternative therapeutic strategies to tackle AL amyloidosis, we produced five llama-derived nanobodies (Nbs) specific against H3, a well-characterized amyloidogenic and cardiotoxic LC from an AL patient with severe cardiac involvement. We found that Nbs are specific and potent agents capable of abolishing H3 soluble toxicity in C. elegans in vivo model. Structural characterization of H3-Nb complexes revealed that the protective effect of Nbs is related to their ability to bind to the H3 VL domain and stabilise an unexpected partially open LC dimer in which the two VL domains no longer interact with each other. Thus, while identifying potent inhibitors of LC soluble toxicity, we also describe the first non-native structure of an amyloidogenic LC that may represent a crucial step in toxicity and aggregation mechanisms.

Pharmacological stabilization of the native state of full-length immunoglobulin light chains to treat light chain amyloidosis.

Immunoglobulin light chain amyloidosis (AL) is a cancer of plasma cells that secrete unstable full-length immunoglobulin light chains. These light chains misfold and aggregate, often with aberrant endoproteolysis, leading to organ toxicity. AL is currently treated by pharmacological elimination of the clonal plasma cells. Since it remains difficult to completely kill these cells in the majority of patients, we seek a complementary drug that inhibits light chain aggregation, which should diminish organ toxicity. We discovered a small-molecule binding site on full-length immunoglobulin light chains by structurally characterizing hit stabilizers emerging from a high-throughput screen seeking small molecules that protect full-length light chains from conformational excursion-linked endoproteolysis. The x-ray crystallographic characterization of 7 structurally distinct hit native-state stabilizers provided a structure-based blueprint, reviewed herein, to design more potent stabilizers. This approach enabled us to transform hits with micromolar affinity into stabilizers with nanomolar dissociation constants that potently prevent light chain aggregation.

Mechanistic insights into the aggregation pathway of the patient-derived immunoglobulin light chain variable domain protein FOR005.

Systemic antibody light chain (AL) amyloidosis is characterized by deposition of amyloid fibrils. Prior to fibril formation, soluble oligomeric AL protein has a direct cytotoxic effect on cardiomyocytes. We focus on the patient derived λ-III AL variable domain FOR005 which is mutated at five positions with respect to the closest germline protein. Using solution-state NMR spectroscopy, we follow the individual steps involved in protein misfolding from the native to the amyloid fibril state. Unfavorable mutations in the complementary determining regions introduce a strain in the native protein structure which yields partial unfolding. Driven by electrostatic interactions, the protein converts into a high molecular weight, oligomeric, molten globule. The high local concentration of aggregation prone regions in the oligomer finally catalyzes the conversion into fibrils. The topology is determined by balanced electrostatic interactions in the fibril core implying a 180° rotational switch of the beta-sheets around the conserved disulfide bond.

Unfolding and Aggregation Pathways of Variable Domains from Immunoglobulin Light Chains.

Light chain amyloidosis is the most common form of systemic amyloidosis. This disease is caused by the formation and deposition of amyloid fibers made from immunoglobulin light chains. Environmental conditions such as pH and temperature can affect protein structure and induce the development of these fibers. Several studies have shed light on the native state, stability, dynamics, and final amyloid state of these proteins; however, the initiation process and the fibril formation pathway remain poorly understood structurally and kinetically. To study this, we analyzed the unfolding and aggregation process of the 6aJL2 protein under acidic conditions, with temperature changes, and upon mutation, using biophysical and computational techniques. Our results suggest that the differences in amyloidogenicity displayed by 6aJL2 under these conditions are caused by traversing different aggregation pathways, including unfolded intermediates and the formation of oligomers.

Widespread amyloidogenicity potential of multiple myeloma patient-derived immunoglobulin light chains.

BACKGROUND: In a range of human disorders such as multiple myeloma (MM), immunoglobulin light chains (IgLCs) can be produced at very high concentrations. This can lead to pathological aggregation and deposition of IgLCs in different tissues, which in turn leads to severe and potentially fatal organ damage. However, IgLCs can also be highly soluble and non-toxic. It is generally thought that the cause for this differential solubility behaviour is solely found within the IgLC amino acid sequences, and a variety of individual sequence-related biophysical properties (e.g. thermal stability, dimerisation) have been proposed in different studies as major determinants of the aggregation in vivo. Here, we investigate biophysical properties underlying IgLC amyloidogenicity. RESULTS: We introduce a novel and systematic workflow, Thermodynamic and Aggregation Fingerprinting (ThAgg-Fip), for detailed biophysical characterisation, and apply it to nine different MM patient-derived IgLCs. Our set of pathogenic IgLCs spans the entire range of values in those parameters previously proposed to define in vivo amyloidogenicity; however, none actually forms amyloid in patients. Even more surprisingly, we were able to show that all our IgLCs are able to form amyloid fibrils readily in vitro under the influence of proteolytic cleavage by co-purified cathepsins. CONCLUSIONS: We show that (I) in vivo aggregation behaviour is unlikely to be mechanistically linked to any single biophysical or biochemical parameter and (II) amyloidogenic potential is widespread in IgLC sequences and is not confined to those sequences that form amyloid fibrils in patients. Our findings suggest that protein sequence, environmental conditions and presence and action of proteases all determine the ability of light chains to form amyloid fibrils in patients.

Taurine promotes B-cell activation by interaction with the VH /VL framework regions of B-cell receptor.

Taurine (Tau) is a special sulphur-containing amino acid and has been widely used as a dietary supplement. Although Tau exists in lymphocytes in large quantities, the physiological significance of Tau to modulate human immunity is unknown. In the present study, we first found that Tau regulates the B-cell receptor (BCR)-mediated signal transduction and induces the B cells activation. The IgG production of mice after ovalbumin immunization was also increased by Tau administration. Moreover, the isothermal titration calorimetry and surface plasmon resonance analysis have shown that Tau specifically bound to the IgG2a-BCR. The Tau could bind to IgG F(ab')2 regions via fluorescence spectroscopy analysis. In the molecular docking analysis, Tau bound to the framework regions (FRs) of variable region of the heavy chains (VH ) and in the light chains (VL ) of IgG2a-BCR. Our results suggested that Tau could improve the activation of B cells by interaction with the VH /VL FRs of BCR.

Functional antibodies exhibit light chain coherence.

The vertebrate adaptive immune system modifies the genome of individual B cells to encode antibodies that bind particular antigens1. In most mammals, antibodies are composed of heavy and light chains that are generated sequentially by recombination of V, D (for heavy chains), J and C gene segments. Each chain contains three complementarity-determining regions (CDR1-CDR3), which contribute to antigen specificity. Certain heavy and light chains are preferred for particular antigens2-22. Here we consider pairs of B cells that share the same heavy chain V gene and CDRH3 amino acid sequence and were isolated from different donors, also known as public clonotypes23,24. We show that for naive antibodies (those not yet adapted to antigens), the probability that they use the same light chain V gene is around 10%, whereas for memory (functional) antibodies, it is around 80%, even if only one cell per clonotype is used. This property of functional antibodies is a phenomenon that we call light chain coherence. We also observe this phenomenon when similar heavy chains recur within a donor. Thus, although naive antibodies seem to recur by chance, the recurrence of functional antibodies reveals surprising constraint and determinism in the processes of V(D)J recombination and immune selection. For most functional antibodies, the heavy chain determines the light chain.

A method to obtain purified free light chain monomers and dimers from urine samples of patients with multiple myeloma.

Antibody light chains are synthesized in excess by plasma cells, and this excess can be secreted into biological fluids as dimers or monomers in various proportions. Structural differences between monomers or dimers of free light chains (FLC) can affect their biological functions and possibly their pathogenicity. They also may exhibit differential immune reactivity, perhaps explaining discrepant quantifications when measured by different immunoreagents. Having purified FLC monomers and dimers available can be useful for studying their properties. Here we propose a simple preparatory procedure to purify FLC monomers and dimers from urine samples of patients with plasma cell disorders. Two representative urine samples containing lambda or kappa FLC were loaded into a nonreducing sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). The gel strips containing separate monomers and dimers were excised, electroeluted, and the FLC recovered. The FLC were recovered from SDS-PAGE gel in sufficient amounts to be quantified by UV and two automated nephelometric assays immunochemical. The procedure was found to be simple, reproducible, and with a high yield, thus offering the opportunity to compare different assays. Not all urine samples are suitable for this procedure, but this approach allows for the purification of FLC monomers and dimers from many selected urine samples which maintain their oligomeric organization.

Molecular Mechanism of Pathogenesis and Treatment Strategies for AL Amyloidosis.

In amyloid light-chain (AL) amyloidosis, small B-cell clones (mostly plasma cell clones) present in the bone marrow proliferate and secrete unstable monoclonal free light chains (FLCs), which form amyloid fibrils that deposit in the interstitial tissue, resulting in organ injury and dysfunction. AL amyloidosis progresses much faster than other types of amyloidosis, with a slight delay in diagnosis leading to a marked exacerbation of cardiomyopathy. In some cases, the resulting heart failure is so severe that chemotherapy cannot be administered, and death sometimes occurs within a few months. To date, many clinical studies have focused on therapeutics, especially chemotherapy, to treat this disease. Because it is necessary to promptly lower FLC, the causative protein of amyloid, to achieve a hematological response, various anticancer agents targeting neoplastic plasma cells are used for the treatment of this disease. In addition, many basic studies using human specimens to elucidate the pathophysiology of AL have been conducted. Gene mutations associated with AL, the characteristics of amyloidogenic LC, and the structural specificity of amyloid fibrils have been clarified. Regarding the mechanism of cellular and tissue damage, the mass effect due to amyloid deposition, as well as the toxicity of pre-fibrillar LC, is gradually being elucidated. This review outlines the pathogenesis and treatment strategies for AL amyloidosis with respect to its molecular mechanisms.

CDR loop interactions can determine heavy and light chain pairing preferences in bispecific antibodies.

As the current biotherapeutic market is dominated by antibodies, the design of different antibody formats, like bispecific antibodies, is critical to the advancement of the field. In contrast to monovalent antibodies, which consist of two identical antigen-binding sites, bispecific antibodies can target two different epitopes by containing two different antigen-binding sites. Thus, the rise of new formats as successful therapeutics has reignited the interest in advancing and facilitating the efficient production of bispecific antibodies. Here, we investigate the influence of point mutations in the antigen-binding site, the paratope, on heavy and light chain pairing preferences by using molecular dynamics simulations. In agreement with experiments, we find that specific residues in the antibody variable domain (Fv), i.e., the complementarity-determining region (CDR) L3 and H3 loops, determine heavy and light chain pairing preferences. Excitingly, we observe substantial population shifts in CDR-H3 and CDR-L3 loop conformations in solution accompanied by a decrease in bispecific IgG yield. These conformational changes in the CDR3 loops induced by point mutations also influence all other CDR loop conformations and consequentially result in different CDR loop states in solution. However, besides their effect on the obtained CDR loop ensembles, point mutations also lead to distinct interaction patterns in the VH-VL interface. By comparing the interaction patterns among all investigated variants, we observe specific contacts in the interface that drive heavy and light chain pairing. Thus, these findings have broad implications in the field of antibody engineering and design because they provide a mechanistic understanding of antibody interfaces, by identifying critical factors driving the pairing preferences, and thus can help to advance the design of bispecific antibodies.

Amyloidogenic immunoglobulin light chain kinetic stabilizers comprising a simple urea linker module reveal a novel binding sub-site.

In immunoglobulin light chain (LC) amyloidosis, the misfolding, or misfolding and misassembly of LC a protein or fragments thereof resulting from aberrant endoproteolysis, causes organ damage to patients. A small molecule "kinetic stabilizer" drug could slow or stop these processes and improve prognosis. We previously identified coumarin-based kinetic stabilizers of LCs that can be divided into four components, including a "linker module" and "distal substructure". Our prior studies focused on characterizing carbamate, hydantoin, and spirocyclic urea linker modules, which bind in a solvent-exposed site at the VL-VL domain interface of the LC dimer. Here, we report structure-activity relationship data on 7-diethylamino coumarin-based kinetic stabilizers. This substructure occupies the previously characterized "anchor cavity" and the "aromatic slit". The potencies of amide and urea linker modules terminating in a variety of distal substructures attached at the 3-position of this coumarin ring were assessed. Surprisingly, crystallographic data on a 7-diethylamino coumarin-based kinetic stabilizer reveals that the urea linker module and distal substructure attached at the 3-position bind a solvent-exposed region of the full-length LC dimer distinct from previously characterized sites. Our results further elaborate the small-molecule binding surface of LCs that could be occupied by potent and selective LC kinetic stabilizers.

Cu(II) Binding Increases the Soluble Toxicity of Amyloidogenic Light Chains.

Light chain amyloidosis (AL) is caused by the aberrant overproduction of immunoglobulin light chains (LCs). The resulting abnormally high LC concentrations in blood lead to deposit formation in the heart and other target organs. Organ damage is caused not only by the accumulation of bulky amyloid deposits, but extensive clinical data indicate that circulating soluble LCs also exert cardiotoxic effects. The nematode C. elegans has been validated to recapitulate LC soluble toxicity in vivo, and in such a model a role for copper ions in increasing LC soluble toxicity has been reported. Here, we applied microscale thermophoresis, isothermal calorimetry and thermal melting to demonstrate the specific binding of Cu2+ to the variable domain of amyloidogenic H7 with a sub-micromolar affinity. Histidine residues present in the LC sequence are not involved in the binding, and yet their mutation to Ala reduces the soluble toxicity of H7. Copper ions bind to and destabilize the variable domains and induce a limited stabilization in this domain. In summary, the data reported here, elucidate the biochemical bases of the Cu2+-induced toxicity; moreover, they also show that copper binding is just one of the several biochemical traits contributing to LC soluble in vivo toxicity.

Observed Antibody Space: A diverse database of cleaned, annotated, and translated unpaired and paired antibody sequences.

The antibody repertoires of individuals and groups have been used to explore disease states, understand vaccine responses, and drive therapeutic development. The arrival of B-cell receptor repertoire sequencing has enabled researchers to get a snapshot of these antibody repertoires, and as more data are generated, increasingly in-depth studies are possible. However, most publicly available data only exist as raw FASTQ files, making the data hard to access, process, and compare. The Observed Antibody Space (OAS) database was created in 2018 to offer clean, annotated, and translated repertoire data. In this paper, we describe an update to OAS that has been driven by the increasing volume of data and the appearance of paired (VH/VL) sequence data. OAS is now accessible via a new web server, with standardized search parameters and a new sequence-based search option. The new database provides both nucleotides and amino acids for every sequence, with additional sequence annotations to make the data Minimal Information about Adaptive Immune Receptor Repertoire compliant, and comments on potential problems with the sequence. OAS now contains 25 new studies, including severe acute respiratory syndrome coronavirus 2 data and paired sequencing data. The new database is accessible at http://opig.stats.ox.ac.uk/webapps/oas/, and all data are freely available for download.

Dynamic protein structures in normal function and pathologic misfolding in systemic amyloidosis.

Dynamic and disordered regions in native proteins are often critical for their function, particularly in ligand binding and signaling. In certain proteins, however, such regions can contribute to misfolding and pathologic deposition as amyloid fibrils in vivo. For example, dynamic and disordered regions can promote amyloid formation by destabilizing the native structure, by directly triggering the aggregation, by promoting protein condensation, or by acting as sites of early proteolytic cleavage that favor a release of aggregation-prone fragments or facilitate fibril maturation. At the same time, enhanced dynamics in the native protein state accelerates proteolytic degradation that counteracts amyloid accumulation in vivo. Therefore, the functional need for dynamic protein regions must be balanced against their inherently labile nature. How exactly this balance is achieved and how is it shifted upon amyloidogenic mutations or post-translational modifications? To illustrate possible scenarios, here we review the beneficial and pathologic roles of dynamic and disordered regions in the native states of three families of human plasma proteins that form amyloid precursors in systemic amyloidoses: immunoglobulin light chain, apolipoproteins, and serum amyloid A. Analysis of structure, stability and local dynamics of these diverse proteins and their amyloidogenic variants exemplifies how disordered/dynamic regions can provide a functional advantage as well as an Achilles heel in pathologic amyloid formation.

A Conservative Point Mutation in a Dynamic Antigen-binding Loop of Human Immunoglobulin λ6 Light Chain Promotes Pathologic Amyloid Formation.

Immunoglobulin light chain (LC) amyloidosis (AL) is a life-threatening human disease wherein free mono-clonal LCs deposit in vital organs. To determine what makes some LCs amyloidogenic, we explored patient-based amyloidogenic and non-amyloidogenic recombinant LCs from the λ6 subtype prevalent in AL. Hydrogen-deuterium exchange mass spectrometry, structural stability, proteolysis, and amyloid growth studies revealed that the antigen-binding CDR1 loop is the least protected part in the variable domain of λ6 LC, particularly in the AL variant. N32T substitution in CRD1 is identified as a driver of amyloid formation. Substitution N32T increased the amyloidogenic propensity of CDR1 loop, decreased its protection in the native structure, and accelerated amyloid growth in the context of other AL substitutions. The destabilizing effects of N32T propagated across the molecule increasing its dynamics in regions ∼30 Å away from the substitution site. Such striking long-range effects of a conservative point substitution in a dynamic surface loop may be relevant to Ig function. Comparison of patient-derived and engineered proteins showed that N32T interactions with other substitution sites must contribute to amyloidosis. The results suggest that CDR1 is critical in amyloid formation by other λ6 LCs.