Benchmarking TriadAb using targets from the second antibody modeling assessment.

Computational modeling and design of antibodies has become an integral part of today's research and development in antibody therapeutics. Here we describe the Triad Antibody Homology Modeling (TriadAb) package, a functionality of the Triad protein design platform that predicts the structure of any heavy and light chain sequences of an antibody Fv domain using template-based modeling. To gauge the performance of TriadAb, we benchmarked against the results of the Second Antibody Modeling Assessment (AMA-II). On average, TriadAb produced main-chain carbonyl root-mean-square deviations between models and experimentally determined structures at 1.10 Å, 1.45 Å, 1.41 Å, 3.04 Å, 1.47 Å, 1.27 Å, 1.63 Å in the framework and the six complementarity-determining regions (H1, H2, H3, L1, L2, L3), respectively. The inaugural results are comparable to those reported in AMA-II, corroborating with our internal bench-based experiences that models generated using TriadAb are sufficiently accurate and useful for antibody engineering using the sequence design capabilities provided by Triad.

Quality assessment of VHH models.

Heavy Chain Only Antibodies are specific to Camelid species. Despite the lack of the light chain variable domain, their heavy chain variable domain (VH) domain, named VHH or nanobody, has promising potential applications in research and therapeutic fields. The structural study of VHH is therefore of great interest. Unfortunately, considering the huge amount of sequences that might be produced, only about one thousand of VHH experimental structures are publicly available in the Protein Data Bank, implying that structural model prediction of VHH is a necessary alternative to obtaining 3D information besides its sequence. The present study aims to assess and compare the quality of predictions from different modelling methodologies. Established comparative & homology modelling approaches to recent Deep Learning-based modelling strategies were applied, i.e. Modeller using single or multiple structural templates, ModWeb, SwissModel (with two evaluation schema), RoseTTAfold, AlphaFold 2 and NanoNet. The prediction accuracy was evaluated using RMSD, TM-score, GDT-TS, GDT-HA and Protein Blocks distance metrics. Besides the global structure assessment, we performed specific analyses of Frameworks and CDRs structures. We observed that AlphaFold 2 and especially NanoNet performed better than the other evaluated softwares. Importantly, we performed molecular dynamics simulations of an experimental structure and a NanoNet predicted model of a VHH in order to compare the global structural flexibility and local conformations using Protein Blocks. Despite rather similar structures, substantial differences in dynamical properties were observed, which underlies the complexity of the task of model evaluation.Communicated by Ramaswamy H. Sarma.

Next-CLL, a New Next-Generation Sequencing-Based Method for Assessment of IGHV Gene Mutational Status in Chronic Lymphoid Leukemia.

Current guidelines for patients with chronic lymphocytic leukemia (CLL) recommend mutation status determination of the clonotypic IGHV gene before treatment initiation to guide the choice of first-line therapy. Currently, commercially available next-generation sequencing (NGS) solutions have technical constraints, as they necessitate at least a 2 × 300 bp sequencing, which restricts their use for routine practice. The cost of the commercial kits also represents an important drawback. We present a new method called Next-CLL, a ready-to-use strategy to evaluate IGHV gene mutation status using any NGS device (including 2 × 150 bp sequencers) in routine diagnostic laboratories. The performance of Next-CLL was validated on genomic DNA and cDNA obtained from 80 patients with CLL at diagnosis. Next-CLL identified a productive clone in 100% of cases, whereas PCR with Sanger sequencing led to a 12.5% failure rate. Next-CLL had 100% concordance with the reference technique for IGHV gene identification and allowed assessment of the IGHV mutation status from the leader sequence, following international guidelines. Comparing a large retrospective series of samples, analyzed by using Sanger sequencing (n = 773) or Next-CLL (n = 352), showed no bias in IGHV usage or mutational status, further validating our strategy in the real-life setting. Next-CLL represents a straightforward workflow for IGHV analysis in routine practice to assess clonal architecture and prognosis of patients with CLL.

De novo design and Rosetta-based assessment of high-affinity antibody variable regions (Fv) against the SARS-CoV-2 spike receptor binding domain (RBD).

The continued emergence of new SARS-CoV-2 variants has accentuated the growing need for fast and reliable methods for the design of potentially neutralizing antibodies (Abs) to counter immune evasion by the virus. Here, we report on the de novo computational design of high-affinity Ab variable regions (Fv) through the recombination of VDJ genes targeting the most solvent-exposed hACE2-binding residues of the SARS-CoV-2 spike receptor binding domain (RBD) protein using the software tool OptMAVEn-2.0. Subsequently, we carried out computational affinity maturation of the designed variable regions through amino acid substitutions for improved binding with the target epitope. Immunogenicity of designs was restricted by preferring designs that match sequences from a 9-mer library of "human Abs" based on a human string content score. We generated 106 different antibody designs and reported in detail on the top five that trade-off the greatest computational binding affinity for the RBD with human string content scores. We further describe computational evaluation of the top five designs produced by OptMAVEn-2.0 using a Rosetta-based approach. We used Rosetta SnugDock for local docking of the designs to evaluate their potential to bind the spike RBD and performed "forward folding" with DeepAb to assess their potential to fold into the designed structures. Ultimately, our results identified one designed Ab variable region, P1.D1, as a particularly promising candidate for experimental testing. This effort puts forth a computational workflow for the de novo design and evaluation of Abs that can quickly be adapted to target spike epitopes of emerging SARS-CoV-2 variants or other antigenic targets.

Immunoglobulin gene analysis in chronic lymphocytic leukemia in the era of next generation sequencing.

Twenty years after landmark publications, there is a consensus that the somatic hypermutation (SHM) status of the clonotypic immunoglobulin heavy variable (IGHV) gene is an important cornerstone for accurate risk stratification and therapeutic decision-making in patients with chronic lymphocytic leukemia (CLL). The IGHV SHM status has traditionally been determined by conventional Sanger sequencing. However, NGS has heralded a new era in medical diagnostics and immunogenetic analysis is following this trend. There is indeed a growing demand for shifting practice and using NGS for IGHV gene SHM assessment, although it is debatable whether it is always justifiable, at least taking into account financial considerations for laboratories with limited resources. Nevertheless, as this analysis impacts on treatment decisions, standardization of both technical aspects, and data interpretation becomes essential. Also, the need for establishing new recommendations and providing dedicated education and training on NGS-based immunogenetics is greater than ever before. Here we address potential and challenges of NGS-based immunogenetics in CLL. We are convinced that this perspective helps the hematological community to better understand the pros and cons of this new technological development for CLL patient management.

Mass Spectrometry-Based Method Targeting Ig Variable Regions for Assessment of Minimal Residual Disease in Multiple Myeloma.

Multiple myeloma is a systemic malignancy of monoclonal plasma cells that accounts for 10% of hematologic cancers. With development of highly effective therapies for multiple myeloma, minimal residual disease (MRD) assessment has emerged as an important end point for management decisions. Currently, serologic assays lack the sensitivity for MRD assessment, and invasive bone marrow sampling with flow cytometry or molecular methods has emerged as the gold standard. We report a sensitive and robust targeted mass spectrometry proteomics method to detect MRD in serum, without the need of invasive, sequential bone marrow aspirates. The method detects Ig-derived clonotypic tryptic peptides predicted by sequencing the clonal plasma cell Ig genes. A heavy isotope-labeled Ig internal standard is added to patient serum at a known concentration, the Ig is enriched in a light chain type specific manner, and proteins are digested and analyzed by targeted mass spectrometry. Peptides from the constant regions of the λ or κ light chains, Ig heavy chains, and clonotypic peptides unique to the patient monoclonal Igs are targeted. This technique is highly sensitive and specific for the patient-specific monoclonal Igs, even in samples negative by multiparametric flow cytometry. Our method can accurately and precisely detect monoclonal protein in serum of patients treated for myeloma and has broad implications for management of hematologic patients.

Ultra-high-throughput sequencing of the immune receptor repertoire from millions of lymphocytes.

High-throughput sequencing of the variable domains of immune receptors (antibodies and T cell receptors (TCRs)) is of key importance in the understanding of adaptive immune responses in health and disease. However, the sequencing of both immune receptor chains (VH+VL or TCRß/δ+TCRα/γ) at the single-cell level for typical samples containing >10(4) lymphocytes is problematic, because immune receptors comprise two polypeptide chains that are encoded by separate mRNAs. Here we present a technology that allows rapid and low-cost determination of a paired immune receptor repertoire from millions of cells with high precision (>97%). Flow focusing is used to encapsulate single cells in emulsions containing magnetic beads for mRNA capture. The mRNA transcripts are then reverse-transcribed, physically linked to their partners by overlap extension PCR, and interrogated by high-throughput paired-end Illumina sequencing. This protocol describes the construction and operation of the flow-focusing device in detail, as well as the bioinformatics pipeline used to interpret the data. The entire procedure can be performed by a single researcher in under 12 h of effort per sample.

Evaluating the Use of Antibody Variable Region (Fv) Charge as a Risk Assessment Tool for Predicting Typical Cynomolgus Monkey Pharmacokinetics.

The pharmacokinetic (PK) behavior of monoclonal antibodies in cynomolgus monkeys (cynos) is generally translatable to that in humans. Unfortunately, about 39% of the antibodies evaluated for PKs in cynos have fast nonspecific (or non-target-mediated) clearance (in-house data). An empirical model relating variable region (Fv) charge and hydrophobicity to cyno nonspecific clearance was developed to gauge the risk an antibody would have for fast nonspecific clearance in the monkey. The purpose of this study was to evaluate the predictability of this empirical model on cyno nonspecific clearance with antibodies specifically engineered to have either high or low Fv charge. These amino acid changes were made in the Fv region of two test antibodies, humAb4D5-8 and anti-lymphotoxin α. The humAb4D5-8 has a typical nonspecific clearance in cynos, and by making it more positively charged, the antibody acquires fast nonspecific clearance, and making it less positively charged did not impact its clearance. Anti-lymphotoxin α has fast nonspecific clearance in cynos, and making it more positively charged caused it to clear even faster, whereas making it less positively charged caused it to clear slower and within the typical range. These trends in clearance were also observed in two other preclinical species, mice and rats. The effect of modifying Fv charge on subcutaneous bioavailability was also examined, and in general bioavailability was inversely related to the direction of the Fv charge change. Thus, modifying Fv charge appears to impact antibody PKs, and the changes tended to correlate with those predicted by the empirical model.

Frequencies of SF3B1, NOTCH1, MYD88, BIRC3 and IGHV mutations and TP53 disruptions in Chinese with chronic lymphocytic leukemia: disparities with Europeans.

We studied 307 consecutive Chinese with chronic lymphocytic leukemia (CLL) in diverse disease-stages before and after diverse therapies for mutations in several CLL-related genes. Mutation frequencies were SF3B1, 5%, NOTCH1, 8%, MYD88, 8%, BIRC3, 2%, TP53, 15% and IGHV, 60%. Several of these frequencies differ from those reported in persons of predominately European descent with CLL. Biological and clinical associations were detected including SF3B1 and NOTCH1 mutations with un-mutated IGHV, MYD88 mutations with mutated IGHV, SF3B1 mutations with fludarabine-resistant CLL and NOTCH1 mutation with advanced Binet disease stage and with +12. The NOTCH1 correlation with briefer survival was confirmed in multivariate analyses but the SF3B1 correlation was confounded by concurrent mutations in TP53 and germline IGHV. We show differences in incidence and prognostic impact of mutations in Chinese and CLL compared with persons of predominately European descent with CLL. These data may give insights into the etiology and biology of CLL and suggests different risk stratification models may be needed for different CLL populations.

Quantitative assessment of the robustness of next-generation sequencing of antibody variable gene repertoires from immunized mice.

BACKGROUND: Next-generation sequencing (NGS) of antibody variable regions has emerged as a powerful tool in systems immunology by providing quantitative molecular information on polyclonal humoral immune responses. Reproducible and robust information on antibody repertoires is valuable for basic and applied immunology studies: thus, it is essential to establish the reliability of antibody NGS data. RESULTS: We isolated RNA from antibody-secreting cells (ASCs) from either 1 mouse or a pool of 9 immunized mice in order to simulate both normal and high diversity populations. Next, we prepared three technical replicates of antibody libraries by RT-PCR from each diversity scenario, which were sequenced using the Illumina MiSeq platform resulting in >106 250 bp paired-end reads per replicate. We then assessed the robustness of antibody repertoire data based on clonal identification defined by amino acid sequence of either full-length VDJ region or the complementarity determining region 3 (CDR3). Leveraging modeling approaches adapted from mathematical ecology, we found that in either diversity scenario both CDR3 and VDJ detection nears completeness indicating deep coverage of ASC repertoires. Additionally, we defined reliability thresholds for accurate quantification and ranking of CDR3s and VDJs. Importantly, we show that both factors-(i) replicate sequencing and (ii) sequencing depth-are crucial for robust CDR3 and VDJ detection and ranking. CONCLUSIONS: In summary, we established widely applicable experimental and computational guidelines for robust antibody NGS and analysis, which will help advance systems immunology studies related to the quantitative profiling of antibody responses following infection and vaccination.

Quantification of peptides from immunoglobulin constant and variable regions by LC-MRM MS for assessment of multiple myeloma patients.

PURPOSE: Quantitative MS assays for Igs are compared with existing clinical methods in samples from patients with plasma cell dyscrasias, for example, multiple myeloma (MM). EXPERIMENTAL DESIGN: Using LC-MS/MS data, Ig constant region peptides, and transitions were selected for LC-MRM MS. Quantitative assays were used to assess Igs in serum from 83 patients. RNA sequencing and peptide-based LC-MRM are used to define peptides for quantification of the disease-specific Ig. RESULTS: LC-MRM assays quantify serum levels of Igs and their isoforms (IgG1-4, IgA1-2, IgM, IgD, and IgE, as well as kappa (κ) and lambda (λ) light chains). LC-MRM quantification has been applied to single samples from a patient cohort and a longitudinal study of an IgE patient undergoing treatment, to enable comparison with existing clinical methods. Proof-of-concept data for defining and monitoring variable region peptides are provided using the H929 MM cell line and two MM patients. CONCLUSIONS AND CLINICAL RELEVANCE: LC-MRM assays targeting constant region peptides determine the type and isoform of the involved Ig and quantify its expression; the LC-MRM approach has improved sensitivity compared with the current clinical method, but slightly higher inter-assay variability. Detection of variable region peptides is a promising way to improve Ig quantification, which could produce a dramatic increase in sensitivity over existing methods, and could further complement current clinical techniques.

Assessment of fully automated antibody homology modeling protocols in molecular operating environment.

The success of antibody-based drugs has led to an increased demand for predictive computational tools to assist antibody engineering efforts surrounding the six hypervariable loop regions making up the antigen binding site. Accurate computational modeling of isolated protein loop regions can be quite difficult; consequently, modeling an antigen binding site that includes six loops is particularly challenging. In this work, we present a method for automatic modeling of the FV region of an immunoglobulin based upon the use of a precompiled antibody x-ray structure database, which serves as a source of framework and hypervariable region structural templates that are grafted together. We applied this method (on common desktop hardware) to the Second Antibody Modeling Assessment (AMA-II) target structures as well as an experimental specialized CDR-H3 loop modeling method. The results of the computational structure predictions will be presented and discussed.

Antibody modeling assessment II. Structures and models.

To assess the state-of-the-art in antibody structure modeling, a blinded study was conducted. Eleven unpublished Fab crystal structures were used as a benchmark to compare Fv models generated by seven structure prediction methodologies. In the first round, each participant submitted three non-ranked complete Fv models for each target. In the second round, CDR-H3 modeling was performed in the context of the correct environment provided by the crystal structures with CDR-H3 removed. In this report we describe the reference structures and present our assessment of the models. Some of the essential sources of errors in the predictions were traced to the selection of the structure template, both in terms of the CDR canonical structures and VL/VH packing. On top of this, the errors present in the Protein Data Bank structures were sometimes propagated in the current models, which emphasized the need for the curated structural database devoid of errors. Modeling non-canonical structures, including CDR-H3, remains the biggest challenge for antibody structure prediction.

Assessment of the effect of triton X-114 on the physicochemical properties of an antibody fragment.

The effect of Triton X-114 on the physicochemical properties of a single-chain antibody fragment (scFv) has been studied. According to the far UV circular dichroism spectroscopy, the secondary structure of the recombinant antibody was not significantly affected by the presence of Triton. From the antibody tertiary structure analysis, it was found that the surfactant could be located around the tryptophan molecules accessible to the solvent, diminishing the polarity of its environment but maintaining most of the protein structure integrity. However, in certain conditions of high temperature and high concentration of denaturant molecules, the presence of TX could compromise the antibody fragment stability. These results represent a previous step in designing scFv purification protocols and should be considered prior to developing scFv liquid-liquid extraction procedures.

Monoclonal antibody binding-site diversity assessment with a cell-based clustering assay.

The diversity of a panel of antibodies that target a specific antigen can be established in various assay formats. In conventional epitope binning assays purified antibodies are tested in a pairwise manner: two antibodies that compete with each other for binding to an antigen are grouped into the same cluster or bin, while they are assigned to two different clusters when they do not compete. Here we present a high through put assay that enables grouping of crude hybridoma supernatants without a need for antibody purification. In addition, the assay does not require recombinant protein, because it is conducted on cells that express the antigen of interest. Hence, one can use the antibody-clustering assay for cell surface proteins that are not amenable to purification. Heavy chain variable region (VH) sequencing shows that VH composition within clusters is conserved. Finally, the assay is in good agreement with a conventional epitope binning assay with purified antigen.

Development of scoring functions for antibody sequence assessment and optimization.

Antibody development is still associated with substantial risks and difficulties as single mutations can radically change molecule properties like thermodynamic stability, solubility or viscosity. Since antibody generation methodologies cannot select and optimize for molecule properties which are important for biotechnological applications, careful sequence analysis and optimization is necessary to develop antibodies that fulfil the ambitious requirements of future drugs. While efforts to grab the physical principles of undesired molecule properties from the very bottom are becoming increasingly powerful, the wealth of publically available antibody sequences provides an alternative way to develop early assessment strategies for antibodies using a statistical approach which is the objective of this paper. Here, publically available sequences were used to develop heuristic potentials for the framework regions of heavy and light chains of antibodies of human and murine origin. The potentials take into account position dependent probabilities of individual amino acids but also conditional probabilities which are inevitable for sequence assessment and optimization. It is shown that the potentials derived from human sequences clearly distinguish between human sequences and sequences from mice and, hence, can be used as a measure of humaness which compares a given sequence with the phenotypic pool of human sequences instead of comparing sequence identities to germline genes. Following this line, it is demonstrated that, using the developed potentials, humanization of an antibody can be described as a simple mathematical optimization problem and that the in-silico generated framework variants closely resemble native sequences in terms of predicted immunogenicity.

Clinical characteristics, response to therapy, and survival of African American patients diagnosed with chronic lymphocytic leukemia: joint experience of the MD Anderson Cancer Center and Duke University Medical Center.

BACKGROUND: Little is known regarding racial disparities in characteristics and outcomes among patients with chronic lymphocytic leukemia (CLL). METHODS: The characteristics and outcomes of untreated African American (AA) patients with CLL (n = 84) were analyzed and compared with a reference nonblack (NB) patient population (n = 1571). RESULTS: At the time of presentation, AA patients had lower median hemoglobin levels (12.9 g/dL vs 13.7 g/dL), higher ß2 microglobulin levels (2.7 mg/dL vs 2.4 mg/dL), greater frequency of constitutional symptoms (27% vs 10%), unmutated immunoglobulin heavy-chain variable region (IGHV) mutation status (65% vs 47%), ζ-chain-associated protein kinase 70 (ZAP70) expression (58% vs 32%), and deletion of chromosome 17p or chromosome 11q (28% vs 17%; P ≤ 02 for each comparison). Fifty-one percent of AA patients and 39% of NB patients required first-line therapy and 91% and 88%, respectively, received chemoimmunotherapy. Overall response rates to treatment were 85% for AA patients and 94% for NB patients (P = .06); and the complete response rates were 56% and 58%, respectively (P = .87). The median survival of AA patients was shorter compared with that of NB patients (event-free survival: 36 months vs 61 months; P = .007; overall survival: 152 months vs not reached; P = .0001). AA race was an independent predictor of shorter event-free and overall survival in multivariable regression models. CONCLUSIONS: The current results indicated that AA patients with CLL have more unfavorable prognostic characteristics and shorter survival compared with their NB counterparts.

Accuracy and coverage assessment of Oryctolagus cuniculus (rabbit) genes encoding immunoglobulins in the whole genome sequence assembly (OryCun2.0) and localization of the IGH locus to chromosome 20.

We report on the analyses of genes encoding immunoglobulin heavy and light chains in the rabbit 6.51× whole genome assembly. This OryCun2.0 assembly confirms previous mapping of the duplicated IGK1 and IGK2 loci to chromosome 2 and the IGL lambda light chain locus to chromosome 21. The most frequently rearranged and expressed IGHV1 that is closest to IG DH and IGHJ genes encodes rabbit VHa allotypes. The partially inbred Thorbecke strain rabbit used for whole-genome sequencing was homozygous at the IGK but heterozygous with the IGHV1a1 allele in one of 79 IGHV-containing unplaced scaffolds and IGHV1a2, IGHM, IGHG, and IGHE sequences in another. Some IGKV, IGLV, and IGHA genes are also in other unplaced scaffolds. By fluorescence in situ hybridization, we assigned the previously unmapped IGH locus to the q-telomeric region of rabbit chromosome 20. An approximately 3-Mb segment of human chromosome 14 including IGH genes predicted to map to this telomeric region based on synteny analysis could not be located on assembled chromosome 20. Unplaced scaffold chrUn0053 contains some of the genes that comparative mapping predicts to be missing. We identified discrepancies between previous targeted studies and the OryCun2.0 assembly and some new BAC clones with IGH sequences that can guide other studies to further sequence and improve the OryCun2.0 assembly. Complete knowledge of gene sequences encoding variable regions of rabbit heavy, kappa, and lambda chains will lead to better understanding of how and why rabbits produce antibodies of high specificity and affinity through gene conversion and somatic hypermutation.