Characterization of a novel format scFv×VHH single-chain biparatopic antibody against metal binding protein MtsA.

Biparatopic antibodies (bpAbs) are engineered antibodies that bind to multiple different epitopes within the same antigens. bpAbs comprise diverse formats, including fragment-based formats, and choosing the appropriate molecular format for a desired function against a target molecule is a challenging task. Moreover, optimizing the design of constructs requires selecting appropriate antibody modalities and adjusting linker length for individual bpAbs. Therefore, it is crucial to understand the characteristics of bpAbs at the molecular level. In this study, we first obtained single-chain variable fragments and camelid heavy-chain variable domains targeting distinct epitopes of the metal binding protein MtsA and then developed a novel format single-chain bpAb connecting these fragment antibodies with various linkers. The physicochemical properties, binding activities, complex formation states with antigen, and functions of the bpAb were analyzed using multiple approaches. Notably, we found that the assembly state of the complexes was controlled by a linker and that longer linkers tended to form more compact complexes. These observations provide detailed molecular information that should be considered in the design of bpAbs.

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.

Camel-Derived Nanobodies as Potent Inhibitors of New Delhi Metallo-ß-Lactamase-1 Enzyme.

The injudicious usage of antibiotics during infections caused by Gram-negative bacteria leads to the emergence of ß-lactamases. Among them, the NDM-1 enzyme poses a serious threat to human health. Developing new antibiotics or inhibiting ß-lactamases might become essential to reduce and prevent bacterial infections. Nanobodies (Nbs), the smallest antigen-binding single-domain fragments derived from Camelidae heavy-chain-only antibodies, targeting enzymes, are innovative alternatives to develop effective inhibitors. The biopanning of an immune VHH library after phage display has helped to retrieve recombinant antibody fragments with high inhibitory activity against recombinant-NDM-1 enzyme. Nb02NDM-1, Nb12NDM-1, and Nb17NDM-1 behaved as uncompetitive inhibitors against NDM-1 with Ki values in the nM range. Remarkably, IC50 values of 25.0 nM and 8.5 nM were noted for Nb02NDM-1 and Nb17NDM-1, respectively. The promising inhibition of NDM-1 by Nbs highlights their potential application in combating particular Gram-negative infections.

Redirect Tropism of Fowl Adenovirus 4 Vector by Modifying Fiber2 with Variable Domain of Heavy-Chain Antibody.

The variable domain of a heavy-chain antibody (VHH) has the potential to be used to redirect the cell tropism of adenoviral vectors. Here, we attempted to establish platforms to simplify the screening of VHHs for their specific targeting function when being incorporated into the fiber of adenovirus. Both fowl adenovirus 4 (FAdV-4) and simian adenovirus 1 (SAdV-1) have two types of fiber, one of which is dispensable for virus propagation and is a proper site for VHH display. An intermediate plasmid, pMD-FAV4Fs, was constructed as the start plasmid for FAdV-4 fiber2 modification. Foldon from phage T4 fibritin, a trigger for trimerization, was employed to bridge the tail/shaft domain of fiber2 and VHHs against human CD16A, a key membrane marker of natural killer (NK) cells. Through one step of restriction-assembly, the modified fiber2 was transferred to the adenoviral plasmid, which was linearized and transfected to packaging cells. Five FAdV-4 viruses carrying the GFP gene were finally rescued and amplified, with three VHHs being displayed. One recombinant virus, FAdV4FC21-EG, could hardly transduce human 293 or Jurkat cells. In contrast, when it was used at a multiplicity of infection of 1000 viral particles per cell, the transduction efficiency reached 51% or 34% for 293 or Jurkat cells expressing exogenous CD16A. Such a strategy of fiber modification was transplanted to the SAdV-1 vector to construct SAdV1FC28H-EG, which moderately transduced primary human NK cells while the parental virus transduced none. Collectively, we reformed the strategy of integrating VHH to fiber and established novel platforms for screening VHHs to construct adenoviral vectors with a specific tropism.

Inferring B Cell Phylogenies from Paired H and L Chain BCR Sequences with Dowser.

Abs are vital to human immune responses and are composed of genetically variable H and L chains. These structures are initially expressed as BCRs. BCR diversity is shaped through somatic hypermutation and selection during immune responses. This evolutionary process produces B cell clones, cells that descend from a common ancestor but differ by mutations. Phylogenetic trees inferred from BCR sequences can reconstruct the history of mutations within a clone. Until recently, BCR sequencing technologies separated H and L chains, but advancements in single-cell sequencing now pair H and L chains from individual cells. However, it is unclear how these separate genes should be combined to infer B cell phylogenies. In this study, we investigated strategies for using paired H and L chain sequences to build phylogenetic trees. We found that incorporating L chains significantly improved tree accuracy and reproducibility across all methods tested. This improvement was greater than the difference between tree-building methods and persisted even when mixing bulk and single-cell sequencing data. However, we also found that many phylogenetic methods estimated significantly biased branch lengths when some L chains were missing, such as when mixing single-cell and bulk BCR data. This bias was eliminated using maximum likelihood methods with separate branch lengths for H and L chain gene partitions. Thus, we recommend using maximum likelihood methods with separate H and L chain partitions, especially when mixing data types. We implemented these methods in the R package Dowser: https://dowser.readthedocs.io.

Applications and challenges in designing VHH-based bispecific antibodies: leveraging machine learning solutions.

The development of bispecific antibodies that bind at least two different targets relies on bringing together multiple binding domains with different binding properties and biophysical characteristics to produce a drug-like therapeutic. These building blocks play an important role in the overall quality of the molecule and can influence many important aspects from potency and specificity to stability and half-life. Single-domain antibodies, particularly camelid-derived variable heavy domain of heavy chain (VHH) antibodies, are becoming an increasingly popular choice for bispecific construction due to their single-domain modularity, favorable biophysical properties, and potential to work in multiple antibody formats. Here, we review the use of VHH domains as building blocks in the construction of multispecific antibodies and the challenges in creating optimized molecules. In addition to exploring traditional approaches to VHH development, we review the integration of machine learning techniques at various stages of the process. Specifically, the utilization of machine learning for structural prediction, lead identification, lead optimization, and humanization of VHH antibodies.

Recombinant Production and Characterization of VHHs/Nanobodies Targeting Tau to Block Fibrillar Assembly.

Tau protein was extensively studied using nuclear magnetic resonance spectroscopy, providing a powerful way to determine interaction sites between Tau and partner proteins. Here we used this analytical tool to describe the epitopes of Tau-specific VHHs (variable domain of the heavy chain of the heavy chain-only antibodies, aka nanobodies) selected from a synthetic library. An in vitro Tau aggregation assay was subsequently used as a functional screen to check VHH efficacy as aggregation inhibitors. We have observed a correlation between the targeted epitope and the aggregation-inhibition capacity of a series of Tau-specific VHHs.

The immunoglobulin heavy chain super enhancer controls class switch recombination in developing B cells.

Class switch recombination (CSR) plays an important role in adaptive immune response by enabling mature B cells to replace the initial IgM by another antibody class (IgG, IgE or IgA). CSR is preceded by transcription of the IgH constant genes and is controlled by the super-enhancer 3' regulatory region (3'RR) in an activation-specific manner. The 3'RR is composed of four enhancers (hs3a, hs1-2, hs3b and hs4). In mature B cells, 3'RR activity correlates with transcription of its enhancers. CSR can also occur in primary developing B cells though at low frequency, but in contrast to mature B cells, the transcriptional elements that regulate the process in developing B cells are ill-known. In particular, the role of the 3'RR in the control of constant genes' transcription and CSR has not been addressed. Here, by using a mouse line devoid of the 3'RR and a culture system that highly enriches in pro-B cells, we show that the 3'RR activity is indeed required for switch transcription and CSR, though its effect varies in an isotype-specific manner and correlates with transcription of hs4 enhancer only.

Serum immunoglobulin or albumin binding single-domain antibodies that enable tailored half-life extension of biologics in multiple animal species.

Single-domain antibody fragments (sdAbs) can be isolated from heavy-chain-only antibodies that occur in camelids or the heavy chain of conventional antibodies, that also occur in camelids. Therapeutic application of sdAbs is often complicated by their low serum half-life. Fusion to sdAb that bind to long-lived serum proteins albumin or IgG can prolong serum half-life of fusion partners. Such studies mostly focused on human application. For half-life prolongation in multiple animal species novel species cross-reacting sdAb are needed. We here describe the isolation from immunized llamas of sdAbs G6 and G13 that bound IgG of 9-10 species analysed, including horse, dog, cat, and swine, as well as sdAb A12 that bound horse, dog, swine and cat albumin. A12 bound albumin with 13 to 271 nM affinity dependent on the species. G13 affinity was difficult to determine by biolayer interferometry due to low and heterogeneous signals. G13 and G6 compete for the same binding domain on Fab fragments. Furthermore, they both lack the hallmark residues typical of camelid sdAbs derived from heavy-chain antibodies and had sequence characteristics typical of human sdAbs with high solubility and stability. This suggests they are derived from conventional llama antibodies. They most likely bind IgG through pairing with VL domains at the VH-VL interface rather than a paratope involving complementarity determining regions. None of the isolated sdAb interfered with FcRn binding to albumin or IgG, and thus do not prevent endosomal albumin/IgG-sdAb complex recycling. Fusions of albumin-binding sdAb A12 to several tetanus neurotoxin (TeNT) binding sdAbs prolonged the terminal serum half-life in piglets to about 4 days, comparable to authentic swine albumin. However, G13 conferred a much lower half-life of 0.84 days. Similarly, in horse, G13 prolonged half-life to only 1.2 days whereas A12 fused to two TeNT binding domains (T6T16A12) had a half-life of 21 days. The high half-life of T6T16A12, which earlier proved to be a highly potent TeNT antitoxin, further supports its therapeutic value. Furthermore, we have identified several additional sdAbs that enable tailored half-life extension of biologicals in multiple animal species.

Two distinct subpopulations of marginal zone B cells exhibit differential antibody-producing capacities and radioresistance.

Marginal zone (MZ) B cells, which are splenic innate-like B cells that rapidly secrete antibodies (Abs) against blood-borne pathogens, are composed of heterogeneous subpopulations. Here, we showed that MZ B cells can be divided into two distinct subpopulations according to their CD80 expression levels. CD80high MZ B cells exhibited greater Ab-producing, proliferative, and IL-10-secreting capacities than did CD80low MZ B cells. Notably, CD80high MZ B cells survived 2-Gy whole-body irradiation, whereas CD80low MZ B cells were depleted by irradiation and then repleted with one month after irradiation. Depletion of CD80low MZ B cells led to accelerated development of type II collagen (CII)-induced arthritis upon immunization with bovine CII. CD80high MZ B cells exhibited higher expression of genes involved in proliferation, plasma cell differentiation, and the antioxidant response. CD80high MZ B cells expressed more autoreactive B cell receptors (BCRs) that recognized double-stranded DNA or CII, expressed more immunoglobulin heavy chain sequences with shorter complementarity-determining region 3 sequences, and included more clonotypes with no N-nucleotides or with B-1a BCR sequences than CD80low MZ B cells. Adoptive transfer experiments showed that CD21+CD23+ transitional 2 MZ precursors preferentially generated CD80low MZ B cells and that a proportion of CD80low MZ B cells were converted into CD80high MZ B cells; in contrast, CD80high MZ B cells stably remained CD80high MZ B cells. In summary, MZ B cells can be divided into two subpopulations according to their CD80 expression levels, Ab-producing capacity, radioresistance, and autoreactivity, and these findings may suggest a hierarchical composition of MZ B cells with differential stability and BCR specificity.

Tissue-Matched IgH Gene Rearrangement of Circulating Tumor DNA Shows Significant Value in Predicting the Progression of Diffuse Large B Cell Lymphoma.

BACKGROUND: Evidence has demonstrated that monitoring of the variable, diversity, and joining gene segments (VDJ) rearrangement of the immunoglobulin (Ig) genes in the circulating tumor DNA (ctDNA) is of value in predicting the outcomes of diffuse large B cell lymphoma (DLBCL). In this study, we investigated the role of VDJ rearrangement proportion in ctDNA for predicting DLBCL progression. METHODS: Patients diagnosed with newly diagnosed DLBCL were included in this study. The VDJ sequences of IgH were detected using next-generation sequencing (NGS) in formalin-fixed paraffin-embedded tissue and/or peripheral blood. The clonotype of the highest proportion in the peripheral blood was defined as the "dominant circulating clonotype," whilst the clonotype of the highest proportion in matched tissue that is detected in peripheral blood was defined as the "dominant tissue-matched clonotype." The decision tree, a machine learning-based methodology, was used to establish a progression-predicting model through a combination of "dominant tissue-matched clonotype" proportion or "dominant circulating clonotype" proportion, and the clinicopathological information, including age, sex, cell of origin, stage, international prognostic index, lactate dehydrogenase, number of extranodal involvements and ß2-microglobulin. RESULTS: A total of 55 patients with eligible sequencing data were used for prognosis analysis, among which 36 patients had matched tissue samples. The concordance rate of "dominant circulating clonotype" and "dominant tissue-matched clonotype" was 19.44% (7/36). The decision tree model showed that the combination of extranodal involvement event and "dominant circulating clonotype" proportion (≥37%) had a clinical value in predicting the prognosis of DLBCL following combined chemotherapy (sensitivity, 0.63; specificity, 0.81; positive prediction value (PPV), 0.59; negative prediction value, 0.83; kappa value, 0.42). Noticeably, the combination of the "dominant tissue-matched clonotype" and extranodal involvement event showed a higher value in predicting the progression (sensitivity, 0.85; specificity, 0.78; PPV, 0.69; kappa value, 0.64). CONCLUSION: IgH proportion detected in the ctDNA samples traced from tissue samples has a high clinical value in predicting the progression of DLBCL.

Antibody production relies on the tRNA inosine wobble modification to meet biased codon demand.

Antibodies are produced at high rates to provide immunoprotection, which puts pressure on the B cell translational machinery. Here, we identified a pattern of codon usage conserved across antibody genes. One feature thereof is the hyperutilization of codons that lack genome-encoded Watson-Crick transfer RNAs (tRNAs), instead relying on the posttranscriptional tRNA modification inosine (I34), which expands the decoding capacity of specific tRNAs through wobbling. Antibody-secreting cells had increased I34 levels and were more reliant on I34 for protein production than naïve B cells. Furthermore, antibody I34-dependent codon usage may influence B cell passage through regulatory checkpoints. Our work elucidates the interface between the tRNA pool and protein production in the immune system and has implications for the design and selection of antibodies for vaccines and therapeutics.

IGH 3'RR recombination uncovers a non-germinal center imprint and c-MYC-dependent IGH rearrangement in unmutated chronic lymphocytic leukemia.

Chronic lymphocytic leukemia (CLL) is an incurable indolent non-Hodgkin lymphoma characterized by tumor B cells that weakly express a B-cell receptor. The mutational status of the variable region (IGHV) within the immunoglobulin heavy chain (IGH) locus is an important prognosis indicator and raises the question of the CLL cell of origin. Mutated IGHV gene CLL are genetically imprinted by activation-induced cytidine deaminase (AID). AID is also required for IGH rearrangements: class switch recombination and recombination between switch Mu (Sµ) and the 3' regulatory region (3'RR) (Sµ-3'RRrec). The great majority of CLL B cells being unswitched led us to examine IGH rearrangement blockade in CLL. Our results separated CLL into two groups on the basis of Sµ-3'RRrec counts per sample: Sµ-3'RRrecHigh cases (mostly unmutated CLL) and Sµ-3'RRrecLow cases (mostly mutated CLL), but not based on the class switch recombination junction counts. Sµ-3'RRrec appeared to be ongoing in Sµ-3'RRrecHigh CLL cells and comparison of Sµ-3'RRrec junction structural features pointed to different B-cell origins for both groups. In accordance with IGHV mutational status and PIM1 mutation rate, Sµ-3'RRrecHigh CLL harbor a non-germinal center experienced B-cell imprint while Sµ-3'RRrecLow CLL are from AID-experienced B cells from a secondary lymphoid organ. In addition to the proposals already made concerning the CLL cell of origin, our study highlights that analysis of IGH recombinatory activity can identify CLL cases from different origins. Finally, on-going Sµ-3'RRrec in Sµ-3'RRrecHigh cells appeared to presumably be the consequence of high c-MYC expression, as c-MYC overexpression potentiated IGH rearrangements and Sµ-3'RRrec, even in the absence of AID for the latter.

BCL3 rearrangements in B-cell lymphoid neoplasms occur in two breakpoint clusters associated with different diseases.

The t(14;19)(q32;q13) often juxtaposes BCL3 with immunoglobulin heavy chain (IGH) resulting in overexpression of the gene. In contrast to other oncogenic translocations, BCL3 rearrangement (BCL3-R) has been associated with a broad spectrum of lymphoid neoplasms. Here we report an integrative whole-genome sequence, transcriptomic, and DNA methylation analysis of 13 lymphoid neoplasms with BCL3-R. The resolution of the breakpoints at single base-pair revealed that they occur in two clusters at 5' (n=9) and 3' (n=4) regions of BCL3 associated with two different biological and clinical entities. Both breakpoints were mediated by aberrant class switch recombination of the IGH locus. However, the 5' breakpoints (upstream) juxtaposed BCL3 next to an IGH enhancer leading to overexpression of the gene whereas the 3' breakpoints (downstream) positioned BCL3 outside the influence of the IGH and were not associated with its expression. Upstream BCL3-R tumors had unmutated IGHV, trisomy 12, and mutated genes frequently seen in chronic lymphocytic leukemia (CLL) but had an atypical CLL morphology, immunophenotype, DNA methylome, and expression profile that differ from conventional CLL. In contrast, downstream BCL3-R neoplasms were atypical splenic or nodal marginal zone lymphomas (MZL) with mutated IGHV, complex karyotypes and mutated genes typical of MZL. Two of the latter four tumors transformed to a large B-cell lymphoma. We designed a novel fluorescence in situ hybridization assay that recognizes the two different breakpoints and validated these findings in 17 independent tumors. Overall, upstream or downstream breakpoints of BCL3-R are mainly associated with two subtypes of lymphoid neoplasms with different (epi)genomic, expression, and clinicopathological features resembling atypical CLL and MZL, respectively.

Presence of leukemic clone-specific immunoglobulin heavy chain rearrangements in neonatal blood spots of children with B-cell precursor acute lymphoblastic leukemia.

INTRODUCTION: Childhood B-cell precursor acute lymphoblastic leukemia (BCP-ALL) can be traced back to birth using leukemic clone-specific immunoglobulin heavy chain (IGH) rearrangements, implying prenatal origin of this disease. METHODS: We retrospectively analyzed neonatal blood spots (Guthrie cards) of 24 patients with childhood BCP-ALL aged 1-9.6 years (median 3.1 years) for the presence of clonotypic IGH rearrangements identified in diagnostic bone marrow samples. Based on the sequences of IGH rearrangements, 2 patient-specific primers were designed for each patient and used in semi-nested polymerase chain reaction for the detection of preleukemic clones at birth. RESULTS: Clonotypic IGH rearrangements were detected in neonatal blood spots of 54.2% of patients (13/24). In two cases with double IGH rearrangements detected at diagnosis, only one rearrangement was present at birth, while in the third case both leukemic rearrangements were detected in neonatal blood. Guthrie card-positive findings were significantly more frequent in children ≤5 years of age than in older children (p = 0.011). Regarding patients' characteristics at birth and at diagnosis, Guthrie card-positivity was not associated with sex, birth weight and mother's age, as well as with white blood cell count, percentage of bone marrow blasts, immunophenotype and the presence of ETV6/RUNX1 and TCF3/PBX1 fusion genes at diagnosis. CONCLUSION: Our study confirms that a large proportion of childhood BCP-ALL originates in utero, regardless of the molecular subtype defined by chromosomal aberrations. The observed trend toward younger age at diagnosis in Guthrie card-positive versus Guthrie card-negative patients implies that the age at diagnosis depends on the presence of preleukemic clone at birth, as well as on the timing of postnatal transforming genetic events.

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.

Enteric pharmacokinetics of monomeric and multimeric camelid nanobody single-domain antibodies.

Single-domain antibodies (sdAbs) derived from Camelidae heavy-chain-only antibodies (also called nanobodies or VHHs) have advantages over conventional antibodies in terms of their small size and stability to pH and temperature extremes, their ability to express well in microbial hosts, and to be functionally multimerized for enhanced properties. For these reasons, VHHs are showing promise as enteric disease therapeutics, yet little is known as to their pharmacokinetics (PK) within the digestive tract. To improve understanding of enteric VHH PK, we investigated the functional and structural stability of monomeric and multimeric camelid VHH-agents following in vitro incubation with intestinal extracts (chyme) from rabbits and pigs or fecal extracts from human sources, and in vivo in rabbits. The results showed that unstructured domains such as epitopic tags and flexible spacers composed of different amino acid sequences were rapidly degraded by enteric proteases while the functional core VHHs were much more stable to these treatments. Individual VHHs were widely variable in their functional stability to GI tract proteases. Some VHH-based agents which neutralize enteric Shiga toxin Stx2 displayed a functional stability to chyme incubations comparable to that of Stx2-neutralizing IgG and IgA mAbs, thus indicating that selected nanobodies can approach the functional stability of conventional immunoglobulins. Enteric PK data obtained from in vitro incubation studies were consistent with similar incubations performed in vivo in rabbit surgical gut loops. These findings have broad implications for enteric use of VHH-based agents, particularly VHH fusion proteins.

Minimal residual disease detection by next-generation sequencing of different immunoglobulin gene rearrangements in pediatric B-ALL.

While the prognostic role of immunoglobulin heavy chain locus (IGH) rearrangement in minimal residual disease (MRD) in pediatric B-acute lymphoblastic leukemia (B-ALL) has been reported, the contribution of light chain loci (IGK/IGL) remains elusive. This study is to evaluate the prognosis of IGH and IGK/IGL rearrangement-based MRD detected by next-generation sequencing in B-ALL at the end of induction (EOI) and end of consolidation (EOC). IGK/IGL rearrangements identify 5.5% of patients without trackable IGH clones. Concordance rates for IGH and IGK/IGL are 79.9% (cutoff 0.01%) at EOI and 81.0% (cutoff 0.0001%) at EOC, respectively. Patients with NGS-MRD < 0.01% at EOI or <0.0001% at EOC present excellent outcome, with 3-year event-free survival rates higher than 95%. IGH-MRD is prognostic at EOI/EOC, while IGK-MRD at EOI/EOC and IGL-MRD at EOI are not. At EOI, NGS identifies 26.2% of higher risk patients whose MRD < 0.01% by flow cytometry. However, analyzing IGK/IGL along with IGH fails to identify additional higher risk patients both at EOI and at EOC. In conclusion, IGH is crucial for MRD monitoring while IGK and IGL have relatively limited value.