4.6

CiteScore

2.2

Impact Factor
  • ISSN 1674-8301
  • CN 32-1810/R
Zheyue Wang, Qi Tang, Bende Liu, Wenqing Zhang, Yufeng Chen, Ningfei Ji, Yan Peng, Xiaohui Yang, Daixun Cui, Weiyu Kong, Xiaojun Tang, Tingting Yang, Mingshun Zhang, Xinxia Chang, Jin Zhu, Mao Huang, Zhenqing Feng. A SARS-CoV-2 neutralizing antibody discovery by single cell sequencing and molecular modeling[J]. The Journal of Biomedical Research, 2023, 37(3): 166-178. DOI: 10.7555/JBR.36.20220221
Citation: Zheyue Wang, Qi Tang, Bende Liu, Wenqing Zhang, Yufeng Chen, Ningfei Ji, Yan Peng, Xiaohui Yang, Daixun Cui, Weiyu Kong, Xiaojun Tang, Tingting Yang, Mingshun Zhang, Xinxia Chang, Jin Zhu, Mao Huang, Zhenqing Feng. A SARS-CoV-2 neutralizing antibody discovery by single cell sequencing and molecular modeling[J]. The Journal of Biomedical Research, 2023, 37(3): 166-178. DOI: 10.7555/JBR.36.20220221

A SARS-CoV-2 neutralizing antibody discovery by single cell sequencing and molecular modeling

More Information
  • Corresponding author:

    Jin Zhu, Huadong Medical Institute of Biotechniques, 293 Zhongshan Road, Nanjing, Jiangsu 210028, China. Tel: +86-25-84514223, E-mail: zhujin1968@njmu.edu.cn

    Mao Huang, Department of Respiratory and Critical Care Medicine, the First Affiliated Hospital of Nanjing Medical University, 300 Guangzhou Road, Nanjing, Jiangsu 210029, China. Tel: +86-25-83718836, E-mail: hm6114@163.com

    Zhenqing Feng, National Health Commission Key Laboratory of Antibody Technique, Jiangsu Province Engineering Research Center of Antibody Drug, Department of Pathology, Nanjing Medical University, 101 Longmian Avenue, Nanjing, Jiangsu 211166, China. Tel: +86-25-86869411, E-mail: fengzhenqing@njmu.edu.cn

  • Received Date: October 04, 2022
  • Revised Date: December 05, 2022
  • Accepted Date: December 08, 2022
  • Available Online: December 11, 2022
  • Although vaccines have been developed, mutations of SARS-CoV-2, especially the dominant B.1.617.2 (delta) and B.1.529 (omicron) strains with more than 30 mutations on their spike protein, have caused a significant decline in prophylaxis, calling for the need for drug improvement. Antibodies are drugs preferentially used in infectious diseases and are easy to get from immunized organisms. The current study combined molecular modeling and single memory B cell sequencing to assess candidate sequences before experiments, providing a strategy for the fabrication of SARS-CoV-2 neutralizing antibodies. A total of 128 sequences were obtained after sequencing 196 memory B cells, and 42 sequences were left after merging extremely similar ones and discarding incomplete ones, followed by homology modeling of the antibody variable region. Thirteen candidate sequences were expressed, of which three were tested positive for receptor binding domain recognition but only one was confirmed as having broad neutralization against several SARS-CoV-2 variants. The current study successfully obtained a SARS-CoV-2 antibody with broad neutralizing abilities and provided a strategy for antibody development in emerging infectious diseases using single memory B cell BCR sequencing and computer assistance in antibody fabrication.
  • The SARS-CoV-2 pandemic broke out in late 2019 and resulted in six million deaths and 628 million infections by the end of October 2022, according to the World Health Organization (WHO). Sharing 82% of its sequence with the SARS virus[1], SARS-CoV-2 is an RNA virus that invades host cells by binding to angiotensin converting enzyme 2 (ACE2) with its receptor binding domain (RBD) within its spike (S) protein[2]. New variants of SARS-CoV-2 have emerged one by one: from the first variant B.1.1.7 to the latest B.1.1.529 (omicron), which had 15 mutations within RBD and were able to escape from neutralizing antibodies[3]. The three-year pandemic saw a dominant shift from SARS-CoV-2 wild type (WT) strain to the variant B.1.617.2 (delta) that was again followed by the omicron, resulting in the need for faster production and broader reactivity drugs, as current vaccines are unable to protect people from infection[4].

    Neutralizing antibodies have become the premier drug developed for infectious diseases for their accessibility from immunized organisms. For SARS-CoV-2, antibodies targeting RBD are more likely to prevent virus entry by blocking the ACE2 receptor or inactivating the spike protein that is an important site for antibody development[5].

    The commonly used antibody production methods are hybridoma and phage display. Hybridoma involves immunized animals and humanization, while phage display obtains antibodies with randomly matched heavy and light chain pairs[6]. Nowadays, single cell sequencing has accelerated the speed of obtaining antibody sequences by seeking certain memory B cell receptor (BCR) genes directly from infected patients without the introduction of animals or phages[79]. Thus, antibodies obtained in this way are fully human in origin with paired heavy and light chains. Conversely, along with the shortened time consumption, the selection of sequence candidates becomes harder. A memory B cell gives a theoretical antibody sequence, and numerous results can be obtained.

    Molecular modeling (MM) has been used in the development of drug design since the 1950s, which explains and predicts the 3D structure or behavior of molecules[10]. By docking with virtual 3D structures generated by computer algorithms, MM could calculate interaction energy to predict how possible a reaction would occur[11]. To date, several specialized software or programs have been invented: Rosetta, Schrödinger, Discovery Studio, GROMACS, Autodock, etc. MM together with single cell sequencing will fasten the speed of antibody discovery since real world experiments, such as enzyme-linked immunosorbent assay (ELISA), may be abridged. Furthermore, MM is the basis of an AI-based antibody design, which has become a promising method, especially after the emergence of the deep learning algorithm AlphaFold[12].

    In this article, single memory B cell BCR sequencing was employed to obtain antibody sequences from a convalescence patient. MM and clustering were performed to select candidate sequences that were further characterized by ELISA, bio-layerinterference (BLI), and pseudovirus neutralization assays. An antibody with broad neutralizing abilities against several SARS-CoV-2 variants, including omicron, was obtained, providing a strategy for antibody development of the emerging infectious disease.

    Twenty-seven donors, who were tested positive for SARS-CoV-2 infection, were treated and recovered three months after the onset between January 24th and April 1st, 2020, and one healthy donor was also recruited. The relevant details of each donor are summarized in Supplementary Table 1 (available online). Eight milliliters of blood were collected from each donor. Sera and peripheral blood mononuclear cells (PBMCs) were isolated and stored at −80 ℃. All procedures were approved by the Ethics Committee of the First People's Hospital of Jiangxia District, Wuhan, China. The study using human blood abided by the principles of the Declaration of Helsinki.

    293T-ACE2 cells and pseudovirus of SARS-CoV-2 WT, B.1.1.7, B.1.351, P.1, B.1.617.2, and B.1.1.529 were provided by Nanjing Vazyme Biotech Co., Ltd. (Nanjing, Jiangsu, China). 293T-ACE2 cells were cultured at 37 ℃ and 5% CO2 with 10% fetal bovine serum in DMEM. HEK 293 freestyle cells were cultured at 37 ℃ and 8% CO2 with 120 rpm in FreeStyle 293 expression medium (Cat. #12338018, Gibco, Carlsbad, CA, USA) for antibody expression.

    Spike protein (Cat. #40589-V08B1, SinoBiological, Beijing, China) and nuclear protein (Cat. #40588-V08B, SinoBiological) were coated in 100 ng/well respectively. Five microliters of serum were added with dilution from 1:1000 to 1:12800 followed by a secondary antibody (1:2000 dilution; Cat. #A0293, Sigma-Aldrich, St. Louis, MO, USA). TMB (Cat. #ab171523, Abcam, Cambridge, UK) was added, and the absorbency was measured at a wavelength of 450 nm.

    PBMCs were isolated from sample W25 and resuspended by 40 mL PBS at the amount of 2.24×108 cells followed by B cell enrichment[13] (Cat. #17954, StemCell Technologies, Vancouver, BC, Canada). B cells (8.2×106) were obtained, incubated with spike and RBD proteins (Cat. Nos. 40589-V08B1 [SinoBiological], 40592-V08H [SinoBiological], T80302 [GenScript, Nanjing, Jiangsu, China], and DRA36 [Novoprotein, Suzhou, Jiangsu, China]), stained with CD19-APC-cy7, IgG-BV421 and IgM-BV605 antibodies (Cat. Nos. 561742, 562581 and 562977, Beckman Coulter, Bria, CA, USA) and sorted into a 96 well plate by flow cytometry.

    Single memory B cell BCR sequencing was performed by GEXSCOPE sequencing technology[14] (Cat. #4183121, Singleron, Nanjing, Jiangsu, China) with mRNA capture and cDNA reverse-transcription. Two PCR primers were designed for the CH1 region of heavy and light antibody chains by two rounds of semi-nest PCR specific for BCR enrichment products. Specific PCR primers for the library were designed by PCR Handle sequences on the primers of the second semi-nest PCR. Libraries were built, and transcriptome was sequenced followed by BCR mapping and V(D)J gene matching.

    Fragment of variable region (Fv) models was created by Antibody Modeling Cascade in Discovery Studio 2019 (DS19), which used an optimized template for homologous modeling[15]. The models were then evaluated by the Verify Protein (Profiles-3D) and Ramachandran plot for reliability. Docking was done by inputting the modeled antibodies and RBD structure (7BWJ in PDB)[16] in the ZDOCK program. Two thousand poses were obtained after the ZDOCK program, and then the most suitable one was inputted into the RDOCK program to refine the antigen-antibody complex. The interaction energy of the configuration with the best score in RDOCK was calculated by the Analyze Trajectory program[17].

    The antibody sequences were clustered by Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/) and adjusted by V gene distribution and amino acid compositions. A phylogenetic tree of the antibody sequences was generated by ChiPlot (https://xiaochi.chifei3d.com/static/xiaochiPlot/src/index.html).

    DNA sequences coding the antibody Fv region were synthesized onto the expressing plasmids that contain the constant region of VH and VL by Genscript. VH and VL plasmids were co-transfected into HEK 293 freestyle cells with PEI (Cat. #AC04L092, Life-iLab, Shanghai, China), and the supernatant was purified by the AKTA Pure 25 system with a Hitrap protein A column (Cat. #29048576, Cytiva, Washington, DC, USA).

    RBD (Cat. #CQ201-00, Vazyme, China) was coated at 2 μg/mL. One hundred microliters of expressed antibody was added at the beginning concentration of 2 μg/mL and a half successively for the following wells. Secondary antibody (1:10 000 dilution; Cat. #SA00001-17, Protentech, Chicago, IL, USA) was added followed by TMB and absorbency measurement at a wavelength of 450 nm.

    For sera, the competitive ELISA assay was employed with the Add&Read SARS-CoV-2 S protein RBD/ACE2 Kit (DD2202, Vazyme, China). Four microliters of seurm samples with dilutions from 1:10 to 1:101 920 were added, followed by 4 μL Tag2-ACE2 and Tag1-RBD. The absorbance ratio at 665/620 nm was measured after the involvement of anti-Tag1-Eu and anti-Tag2-A2 antibodies.

    For antibodies, the competitive ELISA assay was the same as the RBD binding assay, except for the replacement of 100 μL of expressed antibody by the mixture of 50 μL ACE2 (2 μg/mL, CG205-00, Vazyme) and 50 μL of expressed antibody.

    The affinity was measured by BLItz (Octet N1, Sartorius, Göttingen, Niedersachsen, Germany) with Protein G sensors (Cat. #18-5082, Sartorius). The protein G sensors were first soaked in PBS (with 0.02% Tween and 0.1% BSA) for a baseline, loaded with antibodies for 120 s and 30 s in PBS for another baseline. Then, the sensors were soaked in RBD with different concentrations for an association curve and in PBS for a disassociation curve. Kon, Koff, and Kd were calculated after all runs were finished.

    Eleven continuous peptides of RBD (Supplementary Table 2, available online) were synthesized by GL Biochem, Shanghai. The assay was the same as the ELISA binding assay except that peptides were coated at a concentration of 10 μg/mL and that the antibody was diluted 10 times successively after a beginning concentration of 10 μg/mL.

    The pseudovirus (Cat. #DD1402, Vazyme) was diluted to 1×104–2×104 TCID50/mL and mixed with antibodies of different concentrations. Fifty microliters of 293T-ACE2 cells were added and co-cultured with the virus. Relative light unit (RLU) was measured after 48 h of culture. For the sera, the samples were heated for 30 min at 56 ℃ before the assay.

    Serum antibody titers, RBD binding, and RBD-ACE2 inhibition of 27 samples were measured once, while serum neutralization was tested in three replicates and displayed as the mean±SD. Characterization of the antibodies expressed was assessed by ELISA and pseudovirus neutralizing assay in three replicates and displayed as the mean±SD. Affinity was measured with six different concentrations, and Kon, Koff, and Kd were calculated by Fortebio analysis software of BLItz. Properties of antibody amino acid composition were analyzed according to clusters using one-way ANOVA in graphpad prism 9.0. Half maximal effective concentration (EC50) and half maximal inhibitory concentration (IC50) were calculated using nonlinear regression, four parameters in graphpad prism 9.0.

    Blood samples from 27 donors were tested for antibody titer. The titers of antibodies against S protein and nucleoprotein (N) varied a lot: S antibodies and N antibodies could be detected at the diluting ratio of 1:8000 and 1:64000, respectively (Supplementary Fig. 1, available online). However, antibody titer showed personal deviation at the diluting ratio of 1:2000 (Fig. 1A). ELISA assay found that samples W2, W6, W12, W15, W21, W22, and W25 exhibited a strong binding ability to S protein (Fig. 1B), as well as a strong inhibition to RBD-ACE2 interactions (Fig. 1C). Some of these samples were chosen for pseudovirus neutralization assay, and the neutralizing ability was detected (Fig. 1D and Supplementary Table 1). Among these samples, W25 had the highest S antibody titer of 1:64000, its relative EC50 was 1.423, its inhibition IC50 was 5.028 and its neutralization IC50 was 0.025, and thus W25 was selected for single memory B cell sequencing.

    Figure  1.  Convalescence serum was able to neutralize SARS-CoV-2.
    Twenty-seven blood samples were collected from donors who had recovered from SARS-CoV-2 infection. Before sequencing memory B cells, all the sera were detected for antigen recognition. A: Two virus antigens, spike (S) and nucleoprotein (N), were tested in enzyme-linked immunosorbent assay (ELISA) with sera dilution ratio of 1:2000 (P<0.001, n=27). B: Gradient dilution was performed to get a response curve and relative EC50 for further identification of spike protein reactivity. C: The inhibition of RBD binding to its ligand ACE2 was validated for the sera in competitive ELISAs with gradient dilution and the relative IC50 was calculated. D: The neutralization or protection was examined in pseudovirus neutralizing assay. The pseudovirus was loaded with luciferase genes, and RLU was positively related to viral survival. Samples with strong recognition, high inhibition, and neutralization were in colors while others in black. EC50 and IC50 were summarized in Supplementary Table 1. For the binding and inhibition, n=27 and for the neutralization, n=10. Data are shown as mean±SD. Student's t-test was used for statistical analysis. EC50 and IC50 were calculated using nonlinear regression, four parameters in GraphPad Prism 9.0. OD450: absorbency at 450 nm; OD665/620: fluorescence ratio at 665/620 nm; RLU: relative light unit; [Ab]: antibody concentration; EC50: half maximal effective concentration; IC50: half maximal inhibitory concentration.

    PBMCs of W25 were separated, followed by fluorescence activated cell sorting. Spike protein with His-tag was used to pick out S-specific IgG expressing B cells (Fig. 2A). The results showed that the radio of His positive cells was significantly different between convalescence and healthy donors (Fig. 2B). A total of 198 single B cells were obtained and underwent single memory B cell BCR sequencing, followed by transcribing mRNA into cDNA with additional two rounds of semi-nest PCR using specific primers of both the 5′ and 3′ ends (Fig. 2C). A total of 128 sequences were obtained after BCR mapping and packaging from sequence libraries. Forty-two sequences were left by combining highly similar ones and discarding incomplete ones after alignment with the international ImMunoGeneTics information system (IMGT) database (Fig. 2D).

    Figure  2.  Antibody sequences were generated by single memory B cell BCR sequencing.
    A: Single B cells were sorted from peripheral blood mononuclear cells (PBMCs) of sample W25 by first discarding dead cells (left panel), gating for IgG subset (middle panel), and picking out S-His positive cells (right panel). B: The radio of His positive cells in the healthy donor (right panel, 0.031%) was lower than that in sample W25 (left panel, 0.058%). C: Antibody sequences originated from separated single B cells by reverse transcription and two rounds of semi-nest PCR with specific primers. D: Forty-two sequences were left for molecular modeling analysis after PBMCs isolation (2.24×108 cells), B cell separation (8.2×106 cells), single B cells expressing IgG antibodies against S protein (198 cells) and BCR sequence libraries (128 sequences). HC: healthy control.

    The interaction energy given by MM was used to evaluate the potential for antibody neutralization. Here, the ACE2-RBD complex resolved by X-ray[18] was calculated for the interaction energy, resulting in −57.00 kcal/mol that was taken as a standard. Since heavy chain complementarity determining region 3 (HCDR3) played the most important role in an antigen-antibody interaction, it was assumed that only antibodies whose HCDR3-RBD interaction energy was lower or near −57.00 kcal/mol were able to neutralize SARS-CoV-2 by preventing the combination of RBD and ACE2 on host cells.

    Homology models of Fv regions were first built for the 42 antibodies and then qualified by the Verify Protein (Profiles-3D) program and Ramachandran plot (Supplementary Fig. 2, available online). The results (Supplementary Table 3, available online) showed that all the models were theoretically reliable for all scores given by the Profiles-3D program at a higher-than-expected score, and the percentage of amino acids within the credibility interval was higher than 94% in the Ramachandran plot.

    Modeled structures and RBD-antibody complex (7BWJ in PDB)[16] were run in the ZDOCK program, followed by a refinement of the RDOCK program (Fig. 3A). The best pose in the RDOCK program was calculated for interaction energy, and only two antibodies had the HCDR3-RBD interaction energy near −57.00 kcal/mol: −53.09 kcal/mol for H52 and −56.49 kcal/mol for M20 (Fig. 3B). Besides HCDR3-RBD, the Fv-RBD interaction energy was also calculated, and H52 had the lowest energy of −132.96 kcal/mol.

    Figure  3.  Molecular modeling predicted potential SARS-CoV-2 neutralizing antibodies.
    Using molecular modeling programs, schematic plots of antibodies and antigen-antibody complex were created and displayed. A: Fragment of the variable region (Fv) model of the H52 antibody (upper left) was generated by antibody modeling cascade. Receptor binding domain (RBD) structure[16] (upper right) was used as antigen input and H52-RBD complex was created by RDOCK program (lower). B: Interaction energy was calculated for Fv-RBD and HCDR3-RBD. The interaction energy of RBD-ACE2 calculated by the program was –57.00 kcal/mol and was taken as a standard. Only two antibodies, M20 and H52 had HCDR3-RBD energy near –57.00 kcal/mol (M20: –56.49 kcal/mol; H52: –53.09 kcal/mol) and H52 had the lowest Fv-RBD energy (–132.96 kcal/mol). LCDR3: light chain complementarity determining region 3; HCDR3: heavy chain complementarity determining region 3; RBD: receptor binding domain.

    To cover these sequences as much as possible, some of which were selected as representatives to be expressed and validated with further considerations other than the interaction energy outlined above. The 42 antibodies were clustered into six clusters by their amino acid sequences and optimized with V gene types and MM interaction energy. As CDR3 was of great importance in antigen recognition, the length and amino acid compositions were all analyzed (Fig. 4). Among the six clusters, averages of CDR3 length, alkaline, polar and acidic amino acids were calculated. Those with parameters higher than the averages were taken as "elite" antibodies, while those lower than average were taken as control groups. The six clusters had some deviation in their parameters: cluster three had the lowest interaction energy and the fewest polar amino acids (Supplementary Fig. 3, available online). Taking all these factors together, seven "elite" and six control antibodies were selected as candidates for the verification (Supplementary Table 4, available online).

    Figure  4.  Clustering determined candidate sequences to be expressed.
    Six clusters were generated by amino acid similarities. Other parameters, such as interaction energy and amino acid composition, were shown in the heat map. Interaction energies of different parts were listed in columns as a fragment of variable region (Fv), heavy chain (VH), light chain (VL), HCDR3, and LCDR3. Candidate antibodies were selected in each cluster with the consideration of the parameters. Those with lower HCDR3 interaction, longer CDR3, and more polar or alkaline/acidic amino acids were regarded as elite (red star) and those with contrary features were control groups (blue circle). Details are shown in Supplementary Table 4, available online. LCDR3: light chain complementarity determining region 3; HCDR3: heavy chain complementarity determining region 3.

    Seven elite and six control antibodies were expressed and purified. RBD recognition was first tested by ELISA, which found that only three antibodies M20, H45 and H52 were positive, while others showed no binding (Fig. 5A). The binding ability of M20 and H45 to RBD significantly declined in competitive ELISA with ACE2 as an inhibitor (Fig. 5B). As the absorbancy at 450 nm was lower than 0.5 when the concentration of H45 was 1 μg/mL in competitive ELISA, H45 was considered to have no neutralization potential. BLI was used to measure the affinity of the three antibodies. In line with ELISA, H52 showed the best affinity with 10–7 mol and 10–5 for both M20 and H45, respectively (Fig. 5C). Compared with others and the two, M20 and H52 that had neutralization potential, there were longer light chain complementarity determining region 3 (LCDR3) and more serine in HCDR3 (Fig. 5D). The properties of these groups were compared using one-way ANOVA in GraphPad Prism 9.0.

    Figure  5.  Candidate antibodies were verified in real experiments.
    Thirteen candidate antibodies that belonged to six clusters were expressed and purified. A: Antigen recognition was firstly tested by enzyme-linked immunosorbent assay (ELISA). Antibodies in red were "elite" ones while those in black were controls in the clustering analysis. B: The three antibodies, H52, M20, and H45, which were positive in antigen recognition, were further evaluated in competitive ELISA where ACE2 was involved. C: Affinity was measured using bio-layer interference (BLI) with gradient receptor binding domain concentrations and Kon, Koff, Kd were calculated. D: The two antibodies, H52 and M20, which exhibited better binding and inhibition, were analyzed in properties of HCDR3 and LCDR3 (*P<0.05, **P<0.01, n=13). Data are shown as mean±SD. For ELISA, three replicates were included. In BLI, six concentrations were employed, and Kd was calculated according to association and disassociation curve in BLI software BLItz Pro-1.2.0.49. Statistic analysis of antibody properties was performed using one-way ANOVA in GraphPad Prism 9.0. OD450: absorbency at 450 nm; [Ab]: antibody concentration.

    RBD was divided into 11 peptides (Supplementary Table 2), named P1–P11. Among these peptides, P4, P6, and P11 were confirmed to be the top binding areas by ELISA (Fig. 6A ). Additional three peptides (P4-1, P4-2 and P4-3, Supplementary Table 2) were synthesized with the replacement of certain amino acid to alanine. The binding between P4-1/3 to H52 significantly decreased, indicating that P4 plays an important role in the RBD-H52 interaction (Fig. 6B). Mutations in RBD in some important SARS-CoV-2 variants were illustrated, and no mutated sites were found in all three peptides (Fig. 6C). The structure of RBD showed that although the three peptides were not adjacent linearly, they were close spatially. Most mutation sites were detected at the interface of RBD-ACE2, which was opposite to peptides P4, P6, and P11 (Fig. 6D).

    Figure  6.  Antibody H52 bound to conserved regions of the receptor binding domain.
    To study the epitope of antibody H52, peptides of receptor binding domain (RBD) were synthesized with 20 amino acids for each (named P1 to P10) and 23 for the last one (named P11). A: enzyme-linked immunosorbent assay (ELISA) was performed with an antibody gradient concentration from 2 μg/mL. The area under the curve (AUC) was calculated and P4, P6, and P11 showed higher binding ability. B: Three peptides, named P4-1 to P4-3 were synthesized by replacement of a single amino acid to alanine and ELISA was performed again. AUC showed deceased binding ability. C: The mutated amino acids within RBD were highlighted with yellow background for B.1.1.7, B.1.351, P.1, B.1.617.2 (delta), C.37, and B.1.1.529 (omicron). No mutation was found in P4, P6 or P11. D: In a stereo schematic diagram, P4 (red), P6 (blue), and P11 (black) lie in the opposite side, while mutation sites (purple) were focused on RBD-ACE2 interface. AUC was calculated in GraphPad Prism 9.0 by three replicates in ELISA and was shown in a single number. ACE2: angiotensin converting enzyme 2.

    RBD recognition and ACE2 inhibition were employed, in which H52 was able to bind to and inhibit several SARS-CoV-2 variants with different effects (Fig. 7A and B). Neutralizing ability of H52 was confirmed by the pseudovirus neutralization assay, which showed an IC50 of 2.66 μg/mL (Fig. 7C). Besides, different protection was discovered against these variants, and the results showed that H52 was a potent neutralizing antibody against variant P.1 with an IC50 of 0.17 μg/mL, but had an unsatisfying performance against B.1.617.2 variant with an IC50 of 10.54 μg/mL (Fig. 7C).

    Figure  7.  Antibody H52 was capable of neutralizing several SARS-CoV-2 variants.
    The reactive spectrum of antibody H52 was validated. A: The RBD proteins of SARS-CoV-2 variants B.1.1.7, B.1.351, B.1.617.2 (delta), P.1, C.37, and B.1.529 (omicron) were used as antigens. Antibody H52 was diluted with a gradient in enzyme-linked immunosorbent assays (ELISAs). B: The inhibition of RBD-ACE2 interaction was detected in competitive ELISA. Ratio of the mutated RBD and ACE2 was 2:1. C: Pseudovirus neutralizing assays were employed to evaluate antibody protection against SARS-CoV-2 and its variants B.1.1.7, B.1351, P.1, B.1.617.2, and B.1.1.529. Data are shown as mean±SD. Each experiment includes three replicates. IC50 was calculated using nonlinear regression, four parameters in GraphPad Prism 9.0. RBD: receptor binding domain; ACE2: angiotensin converting enzyme 2; OD450: absorbency at 450 nm; RLU: relative light unit; [Ab]: antibody concentration; IC50: half maximal inhibitory concentration.

    Neutralizing antibodies are important drugs against infectious diseases. For the emerging variants of SARS-CoV-2, more strategies for antibody production are needed.

    In the current study, 27 donors who recovered from SARS-CoV-2 infection were recruited. As antibodies underwent a class switch from IgM to IgG with an enhanced humoral immunity[19], the time for sampling was three months after the elimination of the virus to ensure that IgG made up the main subset.

    Single cell sequencing accelerates antibody discovery, and thus became the most popular method in the SARS-CoV-2 antibody discovery[2023]. However, a considerable number of cells is needed. The 10× Genomics and BD Rhapsody platform requires 500 and 2000 cells at least, respectively[24]. The current study obtained H52 antibody from 198 single memory B cells by GEXSCOPE BCR sequencing with specific primers, which encoded the leading and constant sequence, covered a broad read and yielded a higher rate of output sequences.

    The V gene distribution had some features of the 42 antibodies: VH3 occupied the most of all VH genes, which was in line with the common type of B cell repertoire[25], while kappa and lambda light chains were equal in amount, which was different from normal, when kappa is nearly twice the frequency of lambda[26]. VK3-20, VK3-15, VK3-11, VK1-5, VK2-30 and VK1-39-01 are the major six V genes of Kappa chain[27] and the dominant one in the current study is VK1-39-01. In addition, besides VK1-39-01, the bias for VH3-30, VH3-53 and VL3-21 was detected in SARS-CoV-2 infected patients[28], which was inconsistent with our antibody H52 (VL1-40, VH1-18) and M20 (VL2-23, VH3-21) due to a lack of samples.

    MM predicted potential antibodies for neutralization, especially H52 that was characterized by broad reactivities to many SARS-CoV-2 variants, indicating the desirability of this method for sequence analysis. HCDR3 is in the core of antigen-antibody interaction, and we postulate that only if the interaction energy of HCDR3-antigen is lower or close to that of the antigen-receptor will it be able to prevent their reaction. Moreover, the Fv is where the antibody can recognize its target antigen, so we postulate that the lower the energy is, the more powerful the antibody is, when the reaction of antigen to the receptor can be prevented by the antibody. These assumptions were tested in the current study, and antibody H52, which met the threshold and had the lowest Fv-antigen interaction energy, was the most efficient at neutralizing SARS-CoV-2.

    Although MM has been used in predicting antigen-antibody interaction recently[29], its limitations should not be ignored. The major concern about MM is its accuracy— that it does not match experimental results exactly[30], which is caused by a variety of factors: inadequate templates for modeling[31], conformational change in the actual reaction[32] and inaccuracies in the dock algorithm that calculates intricate macromolecules systems[33]. In the curretn study, the MM prediction was based on blocking of ACE2-RBD interactions, which was only one of neutralizing mechanisms. Antibodies that employed other mechanisms, such as steric hindrance or conformational disruption, may not be selected in this way, indicating the existence of false negative errors. As a newly emerged method involved in multi-subjects, MM still needs to be improved for better performance.

    Similar sequences have similar structures and possess similar functions. For antibodies, similar sequences may recognize a certain kind of epitope, which makes clustering a reasonable way to narrow down candidates from numerous sequences. Another advantage of the clustering analysis is that if one representative of a cluster behaves excellently, this cluster will be of value for a detailed study, by which better ones would be discovered or sequence information for functional improvement would be found. In the current study, seven "elite" antibodies from six clusters were tested in ELISA, but only three of them recognized RBD, the rest might recognize non-RBD regions (NTD, or S2 subunit) or even resulted from mis-captured B cells. These clusters do have some but not statistically significant deviaion, which illustrates the complexity of antibody sequences.

    The length of HCDR3 and the framework pairing are regarded as factors for antibody function[3435]. In the six clusters, we selected control sequences with factors higher or lower than the average to investigate key elements in SARS-CoV-2 antibodies. The two antibodies, H52 and M20, showed a better inhibition, which was in line with MM prediction, and exhibited some common characteristics: longer LCDR3 and more serine in HCDR3, which might provide some reference for SARS-CoV-2 antibody selection.

    SARS-CoV-2 invades host cells using RBD. Studies revealed that within its spike protein and structural amino acids, A475, N487, E484, Y453 and N501 play important roles in the interaction between RBD, virus and ACE2 in host cells[18]. Several variants of concern (VOCs), such as B.1.1.7, B.1.351, P.1, B.1.617.2, and B.1.1.529 as proposed by the WHO, possess mutations around those crucial amino acids mentioned above, enhancing virus spreading and evasion[3639]. Antibody H52 here revealed a strong binding to P4 (379–398), P6 (419–438), and P11 (519–541). These areas were away from the RBD-ACE2 interface and had few strands or loops around, providing possibilities for direct antibody contact. The three peptides are not adjacent linearly but close spatially, indicating a conformational epitope of antibody H52. More importantly, the peptides have no mutations among VOCs, suggesting the conservativeness of these areas and that the evolution of SARS-CoV-2 virus is mainly to strengthen RBD-ACE2 interaction. Furthermore, the latest variants of SARS-CoV-2, including omicron and its lineage BA.1, BA.2, BA.3[40] and their recombinant omicron XE, XD and XF, retained conservativeness on peptide P4, P6, and P11. The antibody H52 may retain broad neutralization for these variants and remain to be tested in future experiments.

    The neutralization ability of H52 is not satisfactory enough to fight against the virus, which may just be the result of indirect inhibition by steric hindrance. On the other hand, an indirect inhibition gives rise to a reduced chance of escape, or perhaps the limited neutralization may be enhanced by adding other potent antibodies for cocktail therapies. Moreover, the conserved region P4, P6 and P11 may be used as the target for more powerful antibodies or vaccine production rather than the RBD-ACE2 interface that mutates fast.

    With the emergence of SARS-CoV-2 mutations, several neutralizing antibodies have displayed declined protection, which is called antibody resistance[41]. It is of great significance to develop more powerful antibodies, but this is hard for traditional methods in a limited time. However, single cell sequencing generates numerous antibody sequences, and MM helps select interested ones. For the virus variants, the mutated RBD protein can be uploaded as a template to select antibodies with a high reactivity in silico and thus decrease the time spent, which is helpful to avoid antibody resistance.

    In the current study, we employed single memory B cell BCR sequencing from a donor with a high S protein antibody titer and neutralization. Forty-two sequences from sequence libraries were selected for MM and clustering analysis. Thirteen antibodies were expressed, three of which were tested positive for RBD recognition, and one exhibited a broad neutralization. These results have provided a strategy for screening antibodies from numerous sequences, which is able to expedite the antibody research period. Additionally, we propose that peptides P4, P6, and P11 are conserved regions of SARS-CoV-2 RBD, and may be a target for better antibody and vaccine design.

    This study was supported by the Jiangsu Provincial Key Research and Development Program (Grant No. BE2020616), the National Key R&D Program of China (Grant No. 2018YFC1200603), the National Science and Technology Major Project (Grant No. 2019SWAQ05-5-4) and Jiangsu Key Lab of Cancer Biomarkers, Prevention and Treatment, Collaborative Innovation Center for Cancer Personalized Medicine, Nanjing Medical University.

    We thank those donors whose sera were used in the study. We thank the First People's Hospital of Jiangxia District, Wuhan and Nanjing Vazyme Biotech Co., Ltd for their help on the study.

    CLC number: R392, Document code: A

    The authors reported no conflict of interests.

    △These authors contributed equally to this work.

  • [1]
    Chan JFW, Kok KH, Zhu Z, et al. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan[J]. Emerg Microbes Infect, 2020, 9(1): 221–236. doi: 10.1080/22221751.2020.1719902
    [2]
    Wang Q, Zhang Y, Wu L, et al. Structural and functional basis of SARS-CoV-2 entry by using human ACE2[J]. Cell, 2020, 181(4): 894–904.e9. doi: 10.1016/j.cell.2020.03.045
    [3]
    He X, Hong W, Pan X, et al. SARS-CoV-2 Omicron variant: characteristics and prevention[J]. MedComm, 2021, 2(4): 838–845. doi: 10.1002/mco2.110
    [4]
    Liu C, Ginn HM, Dejnirattisai W, et al. Reduced neutralization of SARS-CoV-2 B. 1.617 by vaccine and convalescent serum[J]. Cell, 2021, 184(16): 4220–4236.e13. doi: 10.1016/j.cell.2021.06.020
    [5]
    Corti D, Purcell LA, Snell G, et al. Tackling COVID-19 with neutralizing monoclonal antibodies[J]. Cell, 2021, 184(12): 3086–3108. doi: 10.1016/j.cell.2021.05.005
    [6]
    Kuhn P, Fühner V, Unkauf T, et al. Recombinant antibodies for diagnostics and therapy against pathogens and toxins generated by phage display[J]. Proteomics Clin Appl, 2016, 10(9–10): 922–948. doi: 10.1002/prca.201600002
    [7]
    Guthmiller JJ, Dugan HL, Neu KE, et al. An efficient method to generate monoclonal antibodies from human B cells[J]. Methods Mol Biol, 2019, 1904: 109–145. https://pubmed.ncbi.nlm.nih.gov/30539468/
    [8]
    Cao Y, Su B, Guo X, et al. Potent neutralizing antibodies against SARS-CoV-2 identified by high-throughput single-cell sequencing of convalescent patients' B cells[J]. Cell, 2020, 182(1): 73–84.e16. doi: 10.1016/j.cell.2020.05.025
    [9]
    Shi R, Shan C, Duan X, et al. A human neutralizing antibody targets the receptor-binding site of SARS-CoV-2[J]. Nature, 2020, 584(7819): 120–124. doi: 10.1038/s41586-020-2381-y
    [10]
    Broughton HB. Molecular modeling[J]. Curr Opin Chem Biol, 1997, 1(3): 392–398. doi: 10.1016/S1367-5931(97)80079-8
    [11]
    Ferreira LG, Dos Santos RN, Oliva G, et al. Molecular docking and structure-based drug design strategies[J]. Molecules, 2015, 20(7): 13384–13421. doi: 10.3390/molecules200713384
    [12]
    Shan S, Luo S, Yang Z, et al. Deep learning guided optimization of human antibody against SARS-CoV-2 variants with broad neutralization[J]. Proc Natl Acad Sci U S A, 2022, 119(11): e2122954119. doi: 10.1073/pnas.2122954119
    [13]
    Kong L, Ju B, Chen Y, et al. Key gp120 glycans pose roadblocks to the rapid development of VRC01-class antibodies in an HIV-1-infected Chinese donor[J]. Immunity, 2016, 44(4): 939–950. doi: 10.1016/j.immuni.2016.03.006
    [14]
    Zheng H, Chen Y, Li J, et al. Longitudinal analyses reveal distinct immune response landscapes in lung and intestinal tissues from SARS-CoV-2-infected rhesus macaques[J]. Cell Rep, 2022, 39(8): 110864. doi: 10.1016/j.celrep.2022.110864
    [15]
    Chen Q, Xian D, Xu W, et al. Affinity improvement of the fully human anti-TSLP recombinant antibody[J]. Mol Med Rep, 2019, 21(2): 759–767. doi: 10.3892/mmr.2019.10880
    [16]
    Ju B, Zhang Q, Ge J, et al. Human neutralizing antibodies elicited by SARS-CoV-2 infection[J]. Nature, 2020, 584(7819): 115–119. doi: 10.1038/s41586-020-2380-z
    [17]
    Zhang Q, Li R, Peng W, et al. Identification of the active constituents and significant pathways of Guizhi-Shaoyao-Zhimu decoction for the treatment of diabetes mellitus based on molecular docking and network pharmacology[J]. Comb Chem High Throughput Screen, 2019, 22(9): 584–598. https://pubmed.ncbi.nlm.nih.gov/31642770/
    [18]
    Lan J, Ge J, Yu J, et al. Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor[J]. Nature, 2020, 581(7807): 215–220. doi: 10.1038/s41586-020-2180-5
    [19]
    Yu K, Lieber MR. Current insights into the mechanism of mammalian immunoglobulin class switch recombination[J]. Crit Rev Biochem Mol Biol, 2019, 54(4): 333–351. doi: 10.1080/10409238.2019.1659227
    [20]
    Sun Y, Ho M. Emerging antibody-based therapeutics against SARS-CoV-2 during the global pandemic[J]. Antib Ther, 2020, 3(4): 246–256. doi: 10.1093/ABT/TBAA025
    [21]
    Wen W, Su W, Tang H, et al. Immune cell profiling of COVID-19 patients in the recovery stage by single-cell sequencing[J]. Cell Discov, 2020, 6: 31. https://pubmed.ncbi.nlm.nih.gov/32377375/
    [22]
    Dugan HL, Stamper CT, Li L, et al. Profiling B cell immunodominance after SARS-CoV-2 infection reveals antibody evolution to non-neutralizing viral targets[J]. Immunity, 2021, 54(6): 1290–1303.e7. doi: 10.1016/j.immuni.2021.05.001
    [23]
    Copin R, Baum A, Wloga E, et al. The monoclonal antibody combination REGEN-COV protects against SARS-CoV-2 mutational escape in preclinical and human studies[J]. Cell, 2021, 184(15): 3949–3961.e11. doi: 10.1016/j.cell.2021.06.002
    [24]
    Gao C, Zhang M, Chen L. The comparison of two single-cell sequencing platforms: BD rhapsody and 10x genomics chromium[J]. Curr Genomics, 2020, 21(8): 602–609. doi: 10.2174/1389202921999200625220812
    [25]
    Guo Y, Chen K, Kwong PD, et al. cAb-Rep: a database of curated antibody repertoires for exploring antibody diversity and predicting antibody prevalence[J]. Front Immunol, 2019, 10: 2365. doi: 10.3389/fimmu.2019.02365
    [26]
    Collins AM, Watson CT. Immunoglobulin light chain gene rearrangements, receptor editing and the development of a self-tolerant antibody repertoire[J]. Front Immunol, 2018, 9: 2249. doi: 10.3389/fimmu.2018.02249
    [27]
    Jackson KJL, Wang Y, Gaeta BA, et al. Divergent human populations show extensive shared IGK rearrangements in peripheral blood B cells[J]. Immunogenetics, 2012, 64(1): 3–14. doi: 10.1007/s00251-011-0559-z
    [28]
    Gaebler C, Wang Z, Lorenzi JCC, et al. Evolution of antibody immunity to SARS-CoV-2[J]. Nature, 2021, 591(7851): 639–644. doi: 10.1038/s41586-021-03207-w
    [29]
    Teraguchi S, Saputri DS, Llamas-Covarrubias MA, et al. Methods for sequence and structural analysis of B and T cell receptor repertoires[J]. Comput Struct Biotechnol J, 2020, 18: 2000–2011. doi: 10.1016/j.csbj.2020.07.008
    [30]
    Singh A, Dauzhenka T, Kundrotas PJ, et al. Application of docking methodologies to modeled proteins[J]. Proteins, 2020, 88(9): 1180–1188. doi: 10.1002/prot.25889
    [31]
    Long X, Jeliazkov JR, Gray JJ. Non-H3 CDR template selection in antibody modeling through machine learning[J]. PEERJ, 2019, 7: e6179. doi: 10.7717/peerj.6179
    [32]
    Vakser IA. Protein-protein docking: from interaction to interactome[J]. Biophys J, 2014, 107(8): 1785–1793. doi: 10.1016/j.bpj.2014.08.033
    [33]
    Ambrosetti F, Jiménez-García B, Roel-Touris J, et al. Modeling antibody-antigen complexes by information-driven docking[J]. Structure, 2020, 28(1): 119–129.e2. doi: 10.1016/j.str.2019.10.011
    [34]
    Sankar K, Hoi KH, Hötzel I. Dynamics of heavy chain junctional length biases in antibody repertoires[J]. Commun Biol, 2020, 3(1): 207. doi: 10.1038/s42003-020-0931-3
    [35]
    Venet S, Ravn U, Buatois V, et al. Transferring the characteristics of naturally occurring and biased antibody repertoires to human antibody libraries by trapping CDRH3 sequences[J]. PLoS One, 2012, 7(8): e43471. doi: 10.1371/journal.pone.0043471
    [36]
    Galloway SE, Paul P, MacCannell DR, et al. Emergence of SARS-CoV-2 B. 1.1. 7 lineage - United States, December 29, 2020-January 12, 2021[J]. MMWR Morb Mortal Wkly Rep, 2021, 70(3): 95–99. doi: 10.15585/mmwr.mm7003e2
    [37]
    Wibmer CK, Ayres F, Hermanus T, et al. SARS-CoV-2 501Y. V2 escapes neutralization by South African COVID-19 donor plasma[J]. Nat Med, 2021, 27(4): 622–625. doi: 10.1038/s41591-021-01285-x
    [38]
    Faria NR, Mellan TA, Whittaker C, et al. Genomics and epidemiology of the P. 1 SARS-CoV-2 lineage in Manaus, Brazil[J]. Science, 2021, 372(6544): 815–821. doi: 10.1126/science.abh2644
    [39]
    Jhun H, Park HY, Hisham Y, et al. SARS-CoV-2 Delta (B. 1.617. 2) variant: a unique T478K mutation in receptor binding motif (RBM) of spike gene[J]. Immune Netw, 2021, 21(5): e32. doi: 10.4110/in.2021.21.e32
    [40]
    Chen J, Wei G. Omicron BA. 2 (B. 1.1. 529.2): high potential to becoming the next dominating variant[J]. Res Sq, 2022, doi: 10.21203/rs.3.rs-1362445/v1.
    [41]
    Beeraka NM, Sukocheva OA, Lukina E, et al. Development of antibody resistance in emerging mutant strains of SARS CoV-2: impediment for COVID-19 vaccines[J]. Rev Med Virol, 2022, 32(5): e2346. doi: 10.1002/RMV.2346
  • Related Articles

    [1]Ya Xie, Zifeng Xu, Yumin Zhang, Yisheng Li, Pengyu Du, Chun Wang. First-episode psychiatric disorder risk from SARS-CoV-2 infection: A clinical analysis with Chinese psychiatric inpatients[J]. The Journal of Biomedical Research, 2025, 39(1): 50-60. DOI: 10.7555/JBR.38.20240005
    [2]Ziyu Zhao, Yi Zhou, Jinxing Guan, Yan Yan, Jing Zhao, Zhihang Peng, Feng Chen, Yang Zhao, Fang Shao. The relationship between compartment models and their stochastic counterparts: A comparative study with examples of the COVID-19 epidemic modeling[J]. The Journal of Biomedical Research, 2024, 38(2): 175-188. DOI: 10.7555/JBR.37.20230137
    [3]Qianqian Chen, Chunmei Hu, Wei Lu, Tianxing Hang, Yan Shao, Cheng Chen, Yanli Wang, Nan Li, Linling Jin, Wei Wu, Hong Wang, Xiaoning Zeng, Weiping Xie. Characteristics of alveolar macrophages in bronchioalveolar lavage fluids from active tuberculosis patients identified by single-cell RNA sequencing[J]. The Journal of Biomedical Research, 2022, 36(3): 167-180. DOI: 10.7555/JBR.36.20220007
    [4]Misevic Gradimir. Single-cell omics analyses with single molecular detection: challenges and perspectives[J]. The Journal of Biomedical Research, 2021, 35(4): 264-276. DOI: 10.7555/JBR.35.20210026
    [5]Liu Zhi, Wang Jianwei, Ge Yiyue, Xu Yuyu, Guo Mengchen, Mi Kai, Xu Rui, Pei Yang, Zhang Qiankun, Luan Xiaoting, Hu Zhibin, Chi Ying, Liu Xingyin. SARS-CoV-2 encoded microRNAs are involved in the process of virus infection and host immune response[J]. The Journal of Biomedical Research, 2021, 35(3): 216-227. DOI: 10.7555/JBR.35.20200154
    [6]Chen Wei, Hu Zhiliang, Yi Changhua, Chi Yun, Xiong Qingfang, Tan Chee Wah, Yi Yongxiang, Wang Lin-Fa. An unusual COVID-19 case with over four months of viral shedding in the presence of low neutralizing antibodies: a case report[J]. The Journal of Biomedical Research, 2020, 34(6): 470-474. DOI: 10.7555/JBR.34.20200099
    [7]Guan Jinxing, Wei Yongyue, Zhao Yang, Chen Feng. Modeling the transmission dynamics of COVID-19 epidemic: a systematic review[J]. The Journal of Biomedical Research, 2020, 34(6): 422-430. DOI: 10.7555/JBR.34.20200119
    [8]Liu Shuying, Lu Shan. Antibody responses in COVID-19 patients[J]. The Journal of Biomedical Research, 2020, 34(6): 410-415. DOI: 10.7555/JBR.34.20200134
    [9]Yuan Huang, Peng Quan, Yongwei Wang, Dongsheng Zhang, Mingwan Zhang, Rui Li, Nan Jiang. Host-guest interaction of β-cyclodextrin with isomeric ursolic acid and oleanolic acid: physicochemical characterization and molecular modeling study[J]. The Journal of Biomedical Research, 2017, 31(5): 395-407. DOI: 10.7555/JBR.31.20160073
    [10]Daniel G Rosen, Zhihong Zhang, Weiwei Shan, Jinsong Liu. Morphological and molecular basis of ovarian serous carcinoma[J]. The Journal of Biomedical Research, 2010, 24(4): 257-263. DOI: 10.1016/S1674-8301(10)60036-X
  • Other Related Supplements

  • Cited by

    Periodical cited type(1)

    1. Hickerson BT, Khalenkov AM, Xie T, et al. Interchangeability of the Assays Used to Assess the Activity of Anti-SARS-CoV-2 Monoclonal Antibodies. Viruses, 2023, 15(8): 1698. DOI:10.3390/v15081698
    1. Hickerson BT, Khalenkov AM, Xie T, et al. Interchangeability of the Assays Used to Assess the Activity of Anti-SARS-CoV-2 Monoclonal Antibodies. Viruses, 2023, 15(8): 1698. DOI:10.3390/v15081698

    Other cited types(0)

Catalog

    Figures(7)

    Article Metrics

    Article views (962) PDF downloads (181) Cited by(1)
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return