F1000Research F1000Research 2046-1402 F1000 Research Limited London, UK 10.12688/f1000research.13612.2 Research Article Articles Crystal structures of a llama VHH antibody BCD090-M2 targeting human ErbB3 receptor

[version 2; peer review: 2 approved]

 Eliseev Igor E. Conceptualization Investigation Methodology Visualization Writing – Original Draft Preparation Writing – Review & Editing https://orcid.org/0000-0002-3344-2513 a 1 Yudenko Anna N. Conceptualization Investigation Methodology Writing – Original Draft Preparation 1 Vysochinskaya Vera V. Investigation 1 Svirina Anna A. Investigation 1 Evstratyeva Anna V. Investigation Resources 2 Drozhzhachih Maria S. Investigation Resources 2 Krendeleva Elena A. Investigation Resources 2 Vladimirova Anna K. Investigation Resources 2 Nemankin Timofey A. Investigation Resources 2 Ekimova Viktoria M. Investigation Resources 2 Ulitin Andrey B. Investigation Resources 2 Lomovskaya Maria I. Investigation 2 Yakovlev Pavel A. Investigation https://orcid.org/0000-0002-6585-9508 2 Bukatin Anton S. Investigation 1 Knyazev Nickolay A. Investigation 1 Moiseenko Fedor V. Investigation 1 Chakchir Oleg B. Funding Acquisition Project Administration Supervision 1 St. Petersburg National Research Academic University RAS, St. Petersburg, 194021, Russian Federation CJSC Biocad, St. Petersburg, 198515, Russian Federation eliseev@spbau.ru

No competing interests were disclosed.

 4 7 2018 2018 7 57 2 7 2018 Copyright: © 2018 Eliseev IE et al. 2018 This is an open access article distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background: The ability of ErbB3 receptor to functionally complement ErbB1-2 and induce tumor resistance to their inhibitors makes it a unique target in cancer therapy by monoclonal antibodies. Here we report the expression, purification and structural analysis of a new anti-ErbB3 single-chain antibody.

Methods: The VHH fragment of the antibody was expressed in E. coli SHuffle cells as a SUMO fusion, cleaved by TEV protease and purified to homogeneity. Binding to the extracellular domain of ErbB3 was studied by surface plasmon resonance. For structural studies, the antibody was crystallized by hanging-drop vapor diffusion in two different forms.

Results: We developed a robust and efficient system for recombinant expression of single-domain antibodies. The purified antibody was functional and bound ErbB3 with K D =15±1 nM. The crystal structures of the VHH antibody in space groups C2 and P1 were solved by molecular replacement at 1.6 and 1.9 A˚ resolution. The high-quality electron density maps allowed us to build precise atomic models of the antibody and the putative paratope. Surprisingly, the CDR H2 existed in multiple distant conformations in different crystal forms, while the more complex CDR H3 had a low structural variability. The structures were deposited under PDB entry codes 6EZW and 6F0D.

Conclusions: Our results may facilitate further mechanistic studies of ErbB3 inhibition by single-chain antibodies. Besides, the solved structures will contribute to datasets required to develop new computational methods for antibody modeling and design.

 cancer therapeutic antibodies receptor tyrosine kinase HER3 single-domain antibody nanobody crystal structure Ministry of Education and Science of the Russian Federation RFMEFI57716X0217 CJSC Biocad This work is funded by the Ministry of Education and Science of the Russian Federation (contract 14.577.21.0217, unique identifier RFMEFI57716X0217) and co-funded by CJSC Biocad. Amendments from Version 1

This version includes new analysis of the antibody binding to the extracellular domain of ErbB3. Based on the referee report by Prof. Ferguson, we revised the experimental SPR protocol. Mainly, we eliminated harsh regeneration step, and greatly increased the length of the injection, so that the binding reaches equilibrium. We updated the Figure 2 and the associated dataset 1 with new SPR data.

 Introduction

Receptor tyrosine kinases ErbB1-4 (HER1-4) receive inputs from growth factors and transmit signals to the cell nucleus, thus regulating key cellular processes such as growth, differentiation, migration, and apoptosis 1 . Aberrations of ErbB signaling, caused by mutations or receptor overexpression, are associated with the development of a wide variety of cancers. The essential role of ErbB receptors in tumor development makes them a unique target in cancer therapy by monoclonal antibodies 2 . Therapeutic antibodies often act on the first stage of signal transduction by inhibiting ligand binding or receptor dimerization.

The first two members of the family, ErbB1 (EGFR, HER1) and ErbB2 (HER2/neu), were early recognized as promising drug targets because of their frequent overexpression in tumors. The examples of successful application of anti-ErbB antibodies in cancer treatment include cetuximab (anti-EGFR) in head and neck cancer therapy and trastuzumab (anti-HER2) in breast cancer treatment. The role of the third member, ErbB3, has long been underestimated because it lacks intrinsic tyrosine kinase activity. However, its ability to form functional dimers with ErbB1-2 and to confer resistance to their inhibitors makes ErbB3 an important drug target 3 . Particularly, it was shown that inhibition of ErbB2 with lapatinib caused transcriptional up-regulation of ErbB3, which was then phosphorylated by residual ErbB2 kinase activity thus limiting antitumor effect 4 . A comprehensive clinical study revealed that ErbB3 overexpression was a significant marker of reduced survival in patients with breast cancer 5 .

This new data stimulated the development of anti-ErbB3 antibodies, which are at various stages of clinical trials 6 . The rational antibody design requires knowledge of molecular mechanism of ErbB3 inhibition. Recently, several structures of ErbB3-antibody complexes were solved 7– 9 . Surprisingly, these structures showed that antibodies target entirely different epitopes on the receptor: extracellular domain I 7 , domains II and IV 8 , or domain III 9 .

In Russia, anti-ErbB3 antibodies are developed by BIOCAD biotechnology company. The phage display selection of antibodies from immunized llamas and subsequent sequencing allowed the identification of several anti-ErbB3 single-chain antibodies. As a part of our ongoing effort to elucidate the molecular mechanism of ErbB3 inhibition and ultimately open up a possibility of therapeutic application of these antibodies, we study their thermodynamic stability 10 , functional properties, and structure. In this work, we describe the expression, purification, crystallization and structural analysis of the variable fragment of an antibody BCD090-M2, which demonstrated an affinity to the extracellular domain of ErbB3 in preliminary experiments.

 Methods Plasmid construction

Gene fragment encoding the VHH fragment of the antibody was cloned into pSolSUMO expression vector (Lucigen), following the manufacturer’s recommendations. Briefly, the fragment was PCR amplified using the primers: forward 5’-aatctgtacttccagggtcaggtgcagctggtgcag-3’, reverse 5’-gtggcggccgctctattatgaggagacggtgaccgt-3’, with the first 18 nucleotides in both primers matching the ends of linearized pSolSUMO vector. Following the amplification, the fragment was mixed with linearized pSolSUMO vector and used to transform chemically competent E. cloni 10G cells (Lucigen). The resultant plasmid pSolSUMO-BCD090-M2, encoding SUMO with N-terminal hexahistidine tag fused to BCD090-M2 through a TEV recognition site, was sequenced and used for further protein expression in E. coli SHuffle T7 Express cells (NEB).

 Protein expression and purification

Chemically competent E. coli SHuffle T7 Express cells (NEB) were transformed by pSolSUMO-BCD090-M2, and single colonies were used to start small-scale overnight cultures. Then 2–4 l bacterial cultures were inoculated by 1:100 volume of overnight culture and grown in 2xYT supplemented with 50 µg/ml kanamycin at 37°C. At OD 0.6–0.8, protein expression was induced by the addition of L-Rhamnose to a final concentration 5 mM, temperature was lowered to 27–29°C and cells were grown overnight for additional 14–15 h. Cells were then harvested by centrifugation at 10000g (5 min), resuspended in IMAC buffer (50 mM Na 2HPO 4 pH 8.0, 0.3 M NaCl, 5 mM Imidazole) with 1 mM PMSF, 0.5 mM EDTA as protease inhibitors and lysed by ultrasonication. Cell debris were pelleted by centrifugation at 40000g (20 min), and the cell extract supernatant was filtered through 0.22 µm membrane. The solution was loaded on 1 ml IMAC column cOmplete (Roche) at 0.5–1 ml/min, the column was washed with 20–40 column volumes of IMAC buffer, and then the protein was eluted by IMAC buffer with 0.3 M Imidazole.

After elution from IMAC column, sample purity was usually higher than 90% as judged by SDS-PAGE, and the protein was cleaved by TEV protease. Sample was first dialyzed against TEV buffer (30 mM Tris pH 8.0, 0.5 mM EDTA, 1 mM DTT) overnight at 4°C, then mixed with TEV protease at 1:40 to 1:80 enzyme to substrate ratio and cleaved for 4 h at 25°C with mild agitation. Histidine-tagged SUMO was then removed by three repeats of negative IMAC chromatography in batch mode. Sample in TEV buffer with 50 mM NaCl and 5 mM Imidazole was mixed with 200 µl Ni-NTA agarose resin (Qiagen) and incubated for 30 mins with mild agitation, the agarose beads were pelleted by a short centrifugation, and the supernatant was taken and used for the next round of SUMO depletion.

Finally, the cleaved VHH antibody BCD090-M2 was purified by an additional polishing step of high-resolution cation exchange chromatography on a MonoS 5/50 GL column (GE Healthcare). The protein was dialyzed against IEX buffer (20 mM Na Acetate pH 6.0) overnight at 4°C, loaded on a pre-equilibrated column, and eluted by IEX buffer with 0–0.5 M NaCl gradient over 20 column volumes. The peak fractions analyzed by SDS-PAGE were pooled, dialyzed overnight against Sample buffer (20 mM HEPES pH 7.5, 50 mM NaCl), and concentrated on a 10 kDa Amicon centrifuge concentrators (Millipore). Protein concentration was measured spectrophotometrically with the parameters ε=27055 M −1cm −1, MW=13955 Da calculated from the amino acid sequence with the ProtParam tool 11 .

For affinity measurements, we produced an extracellular domain (residues 21–643) of the human ErbB3 receptor using a pEE vector with a CMV promoter carrying the ErbB3 gene fragment followed by a hexahistidine tag and a FLAG-tag. CHO-T-HC cells were transfected with PEI and grown one day in HyCell TransFx-C media (GE Healthcare) at 37°C. On the day 2, the temperature was lowered to 32°C, and cells were grown for an additional 8 days. Then the cells were harvested by sterile filtration through Opticap XL capsule filters (Millipore), the clarified culture fluid was supplemented with 1 mM NiCl 2 and 10 mM Imidazole, and loaded on a HisTrap HP (GE Healthcare) column equilibrated with IMAC buffer. The column was washed with 10 volumes of IMAC buffer, and the protein was eluted with 0.3 M Imidazole. After elution from IMAC column, the protein was further purified by size-exclusion chromatography on HiLoad 16/600 Superdex 200pg (GE Healthcare), dialyzed against PBS and concentrated on a 10 kDa Amicon centrifuge concentrators (Millipore).

 Affinity measurement

Interaction of the recombinantly expressed VHH antibody BCD090-M2 with the extracellular domain of the ErbB3 receptor was studied by surface plasmon resonance technique using a Biacore T200 instrument (GE Healthcare). The purified extracellular domain was diluted in 10 mM Na Acetate buffer pH 4.5 to a final concentration of 100 µg/ml and immobilized on CM5 chip via amine coupling with EDC/NHS, following the manufacturer’s recommendations. The BCD090-M2 stock solution was serially diluted in HBS buffer supplemented with BSA (10 mM HEPES pH 7.4, 0.15 M NaCl, 50 µg/ml BSA) to prepare concentration series from 1 µM down to 0.49 nM. Each sample was injected to a cell with immobilized receptor and a reference cell at 10 µl/min flow rate and association/dissociation time of 30 min. All measurements were performed at 37°C. Between the samples the chip surface was equilibrated with HBS for 10 min without an additional regeneration step. The sensograms were reference-subtracted and analyzed in Biacore T200 evaluation software. The equilibrium dissociation constant was obtained by fitting the response measured at 5 s before the end of association phase as a function of analyte concentration. The reported value was a mean calculated from three experiments.

 Crystallization

The BCD090-M2 crystallization conditions were screened using the commercial sparse-matrix screens Classics I and II, AmSO 4 (Qiagen), Clear Strategy I and II, Morpheus (Molecular Dimensions). Crystallization experiments were set up by sitting-drop vapor diffusion method in 96-well plates at 19°C. Each crystallization drop consisted of 100 nl protein solution at a concentration of 17 mg/ml in Sample buffer (20 mM HEPES pH 7.5, 50 mM NaCl) and 100 nl reservoir solution. The screening revealed the two classes of promising conditions, one with salts of carboxylic acids and PEG (Morpheus #73, #76), and the other with divalent cations and PEG (Classics II #64). The crystallization experiments with the identified conditions were reproduced in 24-well Linbro plates by hanging-drop vapor diffusion method with 2 µl drop volume, 1:1 ratio, and 0.5 ml reservoir volume at 20°C. In the case of first crystallization condition (Morpheus #73), well-formed crystals of 0.2–0.3 mm size appeared after 3–4 days in the hanging-drop experiments. Preliminary X-ray experiments showed diffraction up to 1.6 Å, therefore no further optimization was attempted; the crystals for data collection were obtained using the original precipitant solution #73 from Morpheus screen. In the case of the second crystallization condition (Classics II #64), optimization experiments were made to increase crystal size and improve morphology. It appeared that among divalent cations only Cd 2+ was essential for crystallization, the optimized reservoir solution had the following composition: 0.1 M MES pH 6.5, 12% PEG 3350, 5 mM CdSO 4. Large crystals up to 0.7 mm usually appeared after 4–5 days and diffracted below 2.0 Å.

 Data collection and processing

The crystals grown in the first crystallization condition (Morpheus #73) were mounted in loops, cryoprotected in the mother liquor with 25% glycerol, and flash cooled in cold nitrogen gas stream. The crystals grown with Cd 2+ deteriorated upon soaking in different cryoprotectant solutions, and so were mounted in thin-walled quartz capillaries (Hampton Research) for room-temperature data collection. All diffraction data were collected on a Kappa Apex II diffractometer (Bruker AXS) using CuK α radiation. The datasets were integrated with SAINT V8.18C and scaled with SADABS v. 2008/1 software 12 . The crystal grown in the first condition diffracted to 1.6 Å, the unit cell parameters were a=65.76 Å, b=38.93 Å, c=47.48 Å, α= γ=90°, β =102.24°, the space group C2 was determined with XPREP v. 2008/2 12 . Notably, the crystal grown in the second condition with Cd 2+ appeared triclinic with unit cell parameters a=35.77 Å, b=41.53 Å, c=46.49 Å, α=89.99°, β =67.92°, γ=76.06° and two copies of the VHH antibody in the asymmetric unit, and diffracted to a slightly lower resolution of 1.9 Å. The details of data collection and processing are summarized in Table 1.

Data collection and refinement statistics.

Statistics for the highest-resolution shell are shown in parentheses.

Dataset name PDB entry code 6EZW 6F0D
Diffraction source Kappa Apex II (Bruker AXS) Kappa Apex II (Bruker AXS)
Detector Apex II CCD Apex II CCD
Temperature, K 100.0 293.15
Distance, mm 38 38
Image width, ° 1.0 1.0
Images 599 2613
Wavelength, Å 1.54184 1.54184
Resolution range, Å 33.3–1.598 (1.656–1.598) 32.34–1.9 (1.968–1.9)
Space group C 1 2 1 P 1
Unit cell: a b c, Å, α β γ, ◦ 65.76      38.93 47.48           90 102.24 90 35.77      41.53 46.49      89.99 67.92 76.06
Total reflections 80881 (3994) 246903 (9578)
    unique 15511 (1365) 18859 (1885)
Multiplicity 5.2 (2.9) 13.1 (5.1)
Completeness, % 96.50 (81.41) 96.75 (91.94)
Mean I/ σ(I) 13.89 (2.06) 14.86 (2.26)
Wilson B-factor 9.36 15.28
R-merge 0.073 (0.446) 0.144 (0.806)
R-meas 0.081 (0.539) 0.15 (0.906)
R-pim 0.033 (0.296) 0.040 (0.406)
CC1/2 0.998 (0.794) 0.997 (0.412)
CC* 0.999 (0.941) 0.999 (0.764)
Reflections used in refinement 15157 (1266) 18248 (1733)
Reflections used for R-free 1517 (126) 1829 (174)
R-work 0.182 (0.28) 0.174 (0.257)
R-free 0.213 (0.288) 0.218 (0.322)
CC(work) 0.964 (0.810) 0.971 (0.828)
CC(free) 0.945 (0.778) 0.951 (0.674)
Number of non-hydrogen atoms 1171 2119
    protein 981 1962
    ligands 0 2
    solvent 190 155
Protein residues 128 256
RMS (bonds) 0.006 0.014
RMS (angles) 0.82 1.53
Ramachandran
    favored, % 98.41 97.22
    allowed, % 1.59 2.78
    outliers, % 0.00 0.00
Rotamer outliers, % 0.78 1.96
Clashscore 8.89 5.76
Avg. B-factor 14.36 19.86
    protein 12.20 19.29
    ligands – 15.44
    solvent 25.49 27.19
Results and discussion

The llama VHH antibody BCD090-M2 was expressed in soluble form in the cytoplasm of E. coli SHuffle T7 Express cells as a SUMO fusion. The SHuffle strain 13 has deletions of the genes trxB and gor, and constitutively expresses a chromosomal copy of the disulfide bond isomerase DsbC to promote formation of correct disulfide bonds in recombinant proteins. The system proved very efficient for single-chain antibody expression, the typical protein yield in our experiments was 50 mg of a fusion protein per liter of a bacterial culture after the first IMAC step. For the IMAC we used new chromatography media cOmplete (Roche), which withstands high EDTA and DTT concentrations and demonstrates a different binding strength and specificity compared to traditional Ni-NTA resins. Particularly, our histidine-tagged fusion protein eluted at 30 mM Imidazole concentration in gradient elution experiments, and addition of only 5 mM Imidazole to cell extract and wash buffer efficiently suppressed non-specific binding, resulting in greater than 90% purity after the first chromatography step. The protein was then cleaved by TEV protease to obtain untagged antibody BCD090-M2 with nearly native N-terminus, differing from the original sequence by a single additional glycine residue left from the TEV recognition site. The extent of cleavage was monitored by SDS-PAGE and typically was more than 70% at 1:80 enzyme ratio and more than 85% at 1:40 ratio, as shown in Figure 1 (left panel).

Expression and purification of BCD090-M2 VHH antibody.

The antibody was expressed in the cytoplasm of E. coli SHuffle cells as His 6-SUMO fusion, purified by IMAC and cleaved with TEV protease at 1:80 or 1:40 enzyme:substrate ratio. The bottom band marked with a red arrow corresponds to BCD090-M2, the middle band marked with a blue arrow is a histidine-tagged SUMO, and a top band is an uncleaved protein. After the cleavage, BCD090-M2 was further purified by negative IMAC and cation-exchange chromatography on MonoS, the SDS-PAGE analysis of the purified antibody is presented in the right panel.

The bottom band with apparent molecular mass 14 kDa corresponds to processed BCD090-M2, the middle band corresponds to His 6-SUMO (13 kDa) and migrates anomalously slow probably due to positively charged histidine tag, and the top band is the uncleaved protein. The processed BCD090-M2 was separated from His 6-SUMO and intact fusion protein by negative IMAC chromatography in batch regime, and then polished by an additional step of cation exchange chromatography on a MonoS column. After elution from MonoS, the VHH antibody was almost pure as judged by SDS-PAGE shown in Figure 1 (right panel).

The surface plasmon resonance experiment confirmed that the recombinant VHH antibody was functional and efficiently bound to the immobilized receptor. The representative experimental binding sensograms are shown in Figure 2A, and the processed data fitted with an equilibrium binding model is shown in Figure 2B. The experiments showed little variation and yielded mean equilibrium dissociation constant for monovalent binding K D =15±1 nM. Although the monovalent affinity is a fundamental characteristic of antibody-antigen interaction, the avidity of a full-length bivalent antibody can be much higher. Typically, the enhancement of avidity due to a bivalent interaction is 3–4 orders of magnitude and depends strongly on the surface concentration of antigen 14 .

Binding of BCD090-M2 VHH antibody to the extracellular domain of human ErbB3 receptor analyzed by surface plasmon resonance.

The purified extracellular domain (residues 21-643) was immobilized on CM5 chip surface. A) binding sensograms measured for different concentrations of BCD090-M2; B) fitting of the measured response as a function of concentration of free BCD090-M2 in solution with a simple bimolecular equilibrium binding model and averaging over three experiments gives K D =15±1 nM.

The protein was successfully crystallized in two different forms: in space group C2 with a single copy of BCD090-M2 in the asymmetric unit and in P1 with two molecules and two cadmium ions in the unit cell. After data collection and processing, all further data analysis procedures, including phasing, model building, and refinement, were conducted in Phenix software suite v. 1.11.1_2575 15 . To solve the structures by molecular replacement, we selected a set of single-domain antibodies from PDB with the highest homology with BCD090-M2. We then processed the search models with Sculptor 16 to delete residues that were not aligned with the target and to prune sidechains by the Schwarzenbacher algorithm 17 , and performed molecular replacement using Phaser v. 2.7.16 18 . The best molecular replacement solutions for both datasets were obtained with a nanobody targeting complement receptor Vsig4, PDB:5IMK 19 . Then we used phenix.autobuild 20 to automatically build the framework regions of the BCD090-M2 and Coot v. 0.8.6.1 21 to manually fit the missing CDRs into the experimental electron density. The structures were refined using phenix.refine 22 (see Table 1 for refinement statistics) and deposited to the Protein Data Bank under entry codes 6EZW and 6F0D. All figures were generated using PyMOL.

The overall structure of BCD090-M2 crystallized in space group C2 is shown in Figure 3. The framework regions are drawn as grey ribbons, and the CDRs are colored in orange (H1), blue (H2), and red (H3). The CDRs were defined according to Kabat 23 as shown in Figure 4.

Ribbon diagram of the crystal structure of BCD090-M2 in space group C2 at 1.6 Å resolution (PDB:6EZW).

Framework regions are grey and CDRs are colored.

 Amino acid sequence of BCD090-M2.

The numbering is the same as in the PDB files 6EZW and 6F0D. The CDR regions according to Kabat definition are indicated with color and frames. In the case of CDR H1, the AbM definition (dashed frame, residues 27-36) was used.

For the CDR H1, we also used AbM definition, a combination of Kabat and Chotia 24 definitions, which better matches the loop in protein structure. The high resolution of the dataset and good quality of the electron density maps allowed us to build a precise atomic model of the antibody. We then used a PyIgClassify 25 to analyze the structures of BCD090-M2 CDRs using a comprehensive database of antibody CDR loop conformations. The CDR H1 and H2 belonged to two large common clusters, H1-13-1 and H2-10-2 respectively. The CDR H3 was not assigned to any known structure cluster, probably due to size (18 residues) and complex structure of the loop, which has a cis-proline and four aromatic residues. It is consistent with an observation that CDR H3 loops are very diverse in structure and a few clusters of significant size are present in the database.

The structure of the BCD090-M2 dimer crystallized in space group P1 with cadmium ions is shown in Figure 5. Two cadmium ions in the unit cell are pictured as green spheres, and two symmetry-related ions from the neighbor unit cell are shown in dots. Each cadmium ion is bound by residues Asp100, Glu114, and Asp116 belonging to CDR H3, and by N-terminal glycine residue of an antibody molecule from neighbor unit cell. Thus, intermolecular interaction through cadmium ions effectively define the lattice of a crystal. This finding is consistent with a previous general observation that cadmium can induce the formation of protein crystals or improve their quality 26 .

Ribbon diagram of the crystal structure of BCD090-M2 in space group P1 at 1.9 Å resolution (PDB:6F0D).

Asymmetric unit contains two copies of the VHH antibody and two cadmium ions which are pictured as green spheres, two symmetry-related ions from the neighbor unit cell are shown in dots. Framework regions are grey and CDRs are colored.

Finally, we analyzed the structural variability of the overall antibody fold and the putative paratope by comparing the three solved structures 6EZW, 6F0D (A), and 6F0D (B). We first performed the structural alignment based on C α atoms of the framework regions, which is shown in Figure 6A.

Structural comparison of BCD090-M2 crystallized in different space groups.

A) the three structures, 6EZW, 6F0D chain A, and 6F0D chain B, were superimposed, the structural alignment was based on C α atoms of framework regions. CDRs H1 and H3 showed little structural variation, while CDR H2 adopted different conformations in 6EZW and 6F0D, which are shown in blue and green color. B) Electron density map (F-obs, Phi-model, 1.0 σ) for CDR H2 (residues 51-67) of 6EZW. The quality of electron density maps allowed us to unambiguously trace the loop and place amino-acid side chains.

The overall fold of the antibody had a very low structural variability, the RMSD between framework C α atoms in 6EZW and 6F0D was 0.19 Å. The conformational mobility of the short CDR H1 was also low, with C α RMSD between 6EZW and 6F0D 0.28 and 0.31 Å for chains A and B, respectively. Surprisingly, the most complex CDR H3 had a rigid conformation, which was not significantly altered by binding of a cadmium ion and intermolecular interaction in the asymmetric unit of 6F0D. The C α RMSD for CDR H3 between two structures was 0.61 (chain A) and 0.57 Å (chain B). The largest structural variation was observed in the loop H2, as seen in Figure 6A. The C α RMSD was 1.81 and 1.77 Å for chains A and B, respectively, and the largest difference was observed in the region Arg53-Leu57. Despite the probable flexibility of this loop, in both crystals it was well-resolved in the electron density maps, as shown in Figure 6B, allowing us to unambiguously place all residues. As a result of the observed structural change, the CDRs H2 in 6F0D were attributed to a different from 6EZW distant structural clusters H2-10-6 (chain A) or H2-10-7 (chain B) by PyIgClassify 25 .

http://dx.doi.org/10.5256/f1000research.13612.d209631

Conclusions

In conclusion, here were present expression and purification, functional and structural analysis of a new single-domain llama antibody against human ErbB3. We crystallized the antibody in two different forms and solved high-resolution structures, giving an insight to the organization of the putative ErbB3 binding paratope. We believe this data may facilitate further studies of mechanisms of ErbB3 inhibition by single-chain antibodies. Besides, the solved structures will contribute to datasets that are required to develop new robust computational methods for antibody modeling and design.

 Data availability

The data referenced by this article are under copyright with the following copyright statement: Copyright: � 2018 Eliseev IE et al.

Data associated with the article are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CC0 1.0 Public domain dedication).

The atomic coordinates and structure factors can be accessed under PDB codes 6EZW and 6F0D.

Dataset 1: Uncropped gel from Figure 1 and raw output data from Biacore software. DOI, 10.5256/f1000research.13612.d209631 27

The plasmids and recombinant proteins used in this study are available from Igor Eliseev (corresponding author) upon request.

 Acknowledgements

We acknowledge technical support by the SPC facility at EMBL Hamburg. We also thank Ms. Maria V. Mitkevich for remarkable administrative support.

 Citri A Yarden Y : EGF-ERBB signalling: towards the systems level. Nat Rev Mol Cell Biol. 2006;7(7):505–516. 16829981 10.1038/nrm1962 Polanovski OL Lebedenko EN Deyev SM : ERBB oncogene proteins as targets for monoclonal antibodies. Biochemistry (Mosc). 2012;77(3):227–245. 22803941 10.1134/S0006297912030029 Kol A Terwisscha van Scheltinga AG Timmer-Bosscha H : HER3, serious partner in crime: therapeutic approaches and potential biomarkers for effect of HER3-targeting. Pharmacol Ther. 2014;143(1):1–11. 24513440 10.1016/j.pharmthera.2014.01.005 Garrett JT Olivares MG Rinehart C : Transcriptional and posttranslational up-regulation of HER3 (ErbB3) compensates for inhibition of the HER2 tyrosine kinase. Proc Natl Acad Sci U S A. 2011;108(12):5021–5026. 21385943 10.1073/pnas.1016140108 3064360 Chiu CG Masoudi H Leung S : HER-3 overexpression is prognostic of reduced breast cancer survival: a study of 4046 patients. Ann Surg. 2010;251(6):1107–1116. 20485140 10.1097/SLA.0b013e3181dbb77e Zhang N Chang Y Rios A : HER3/ErbB3, an emerging cancer therapeutic target. Acta Biochim Biophys Sin (Shanghai). 2016;48(1):39–48. 26496898 10.1093/abbs/gmv103 Mirschberger C Schiller CB Schräml M : RG7116, a therapeutic antibody that binds the inactive HER3 receptor and is optimized for immune effector activation. Cancer Res. 2013;73(16):5183–5194. 23780344 10.1158/0008-5472.CAN-13-0099 Garner AP Bialucha CU Sprague ER : An antibody that locks HER3 in the inactive conformation inhibits tumor growth driven by HER2 or neuregulin. Cancer Res. 2013;73(19):6024–6035. 23928993 10.1158/0008-5472.CAN-13-1198 3924863 Lee S Greenlee EB Amick JR : Inhibition of ErbB3 by a monoclonal antibody that locks the extracellular domain in an inactive configuration. Proc Natl Acad Sci U S A. 2015;112(43):13225–13230. 26460020 10.1073/pnas.1518361112 4629334 Eliseev IE Yudenko AN Besedina NA : Thermodynamic analysis of the conformational stability of a single-domain therapeutic antibody. Tech Phys Lett. 2017;43(12):1088–1091. 10.1134/S1063785017120045 Gasteiger E Hoogland C Gattiker A : Protein identification and analysis tools on the ExPASy server.In JM Walker, editor, The Proteomics Protocols Handbook. Humana Press,2005;571–607. 10.1385/1-59259-890-0:571 Bruker: SAINT, SADABS, XPREP. Bruker AXS Inc., Madison, Wisconsin, USA, 2008–2012. Lobstein J Emrich CA Jeans C : SHuffle, a novel Escherichia coli protein expression strain capable of correctly folding disulfide bonded proteins in its cytoplasm. Microb Cell Fact. 2012;11(1):56. 22569138 10.1186/1475-2859-11-56 3526497 Kaufman EN Jain RK : Effect of bivalent interaction upon apparent antibody affinity: experimental confirmation of theory using fluorescence photobleaching and implications for antibody binding assays. Cancer Res. 1992;52(15):4157–4167. 1638531 Adams PD Afonine PV Bunkóczi G : PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 2):213–221. 20124702 10.1107/S0907444909052925 2815670 Bunkóczi G Read RJ : Improvement of molecular-replacement models with Sculptor. Acta Crystallogr D Biol Crystallogr. 2011;67(Pt 4):303–312. 21460448 10.1107/S0907444910051218 3069745 Schwarzenbacher R Godzik A Grzechnik SK : The importance of alignment accuracy for molecular replacement. Acta Crystallogr D Biol Crystallogr. 2004;60(Pt 7):1229–1236. 15213384 10.1107/S0907444904010145 McCoy AJ Grosse-Kunstleve RW Adams PD : Phaser crystallographic software. J Appl Crystallogr. 2007;40(Pt 4):658–674. 19461840 10.1107/S0021889807021206 2483472 Wen Y Ouyang Z Schoonooghe S : Structural evaluation of a nanobody targeting complement receptor Vsig4 and its cross reactivity. Immunobiology. 2017;222(6):807–813. 27889311 10.1016/j.imbio.2016.11.008 Terwilliger TC Grosse-Kunstleve RW Afonine PV : Iterative model building, structure refinement and density modification with the PHENIX AutoBuild wizard. Acta Crystallogr D Biol Crystallogr. 2008;64(Pt 1):61–69. 18094468 10.1107/S090744490705024X 2394820 Emsley P Lohkamp B Scott WG : Features and development of Coot. Acta Crystallogr D Biol Crystallogr. 2010;66(Pt 4):486–501. 20383002 10.1107/S0907444910007493 2852313 Afonine PV Grosse-Kunstleve RW Echols N : Towards automated crystallographic structure refinement with phenix.refine. Acta Crystallogr D Biol Crystallogr. 2012;68(Pt 4):352–367. 22505256 10.1107/S0907444912001308 3322595 Kabat EA Te Wu T Foeller C : Sequences of proteins of immunological interest. NIH Publication No. 91-3242, 5th edition,1991. Reference Source Al-Lazikani B Lesk AM Chothia C : Standard conformations for the canonical structures of immunoglobulins. J Mol Biol. 1997;273(4):927–948. 9367782 10.1006/jmbi.1997.1354 Adolf-Bryfogle J Xu Q North B : PyIgClassify: a database of antibody CDR structural classifications. Nucleic Acids Res. 2015;43(Database issue):D432–D438. 25392411 10.1093/nar/gku1106 4383924 Trakhanov S Kreimer DI Parkin S : Cadmium-induced crystallization of proteins: II. Crystallization of the Salmonella typhimurium histidine-binding protein in complex with L-histidine, L-arginine, or L-lysine. Protein Sci. 1998;7(3):600–604. 9541391 10.1002/pro.5560070308 2143943 Eliseev IE Yudenko AN Vysochinskaya VV : Dataset 1 in: Crystal structures of a llama VHH antibody BCD090-M2 targeting human ErbB3 receptor. F1000Research. 2018. Data Source 10.5256/f1000research.16763.r35737 Reviewer response for version 2 Ferguson Kathryn M. 1 Referee Department of Physiology, Perelman School of Medicine, Philadelphia, PA, USA

Competing interests: No competing interests were disclosed.

 15 8 2018 Copyright: © 2018 Ferguson KM 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

The SPR data in this revised manuscript are very much improved.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

 10.5256/f1000research.14786.r32141 Reviewer response for version 1 Ferguson Kathryn M. 1 Referee Department of Physiology, Perelman School of Medicine, Philadelphia, PA, USA

Competing interests: No competing interests were disclosed.

 10 4 2018 Copyright: © 2018 Ferguson KM 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve-with-reservations

The authors purify a VHH fragment corresponding to an camelid single-chain antibody previously selected from a Llama derived phage antibody library. They measure binding of this VHH (BCD090-M2) to the extracellular region of ErbB3 that was immobilized to a CM5 biosenor chip, using Surface plasmon resonance (SPR). An apparent K D of 1μM is reported. The authors solve the X-ray crystal structure of the VHH in two different crystal forms and observe differences in the conformation of CDR H2. 

The SPR sensograms do not come to a plateau during each injection, which is most obvious at the higher concentrations, so it is not a saturation binding response for each concentration. The Response data that are plotted in Figure 2B for each concentration of VHH do not correspond to the response values 5 seconds before the end of the injection in the sensograms plotted in Figure 2A. It looks more like the values 50 seconds before end of injection. There is something wrong here. Perhaps the data in A and B are from two different experiments? No error is given on the K D value so it is not clear if these data were repeated. I also worry about use of regeneration using glycine at pH 2. It is highly likely that this will denature the immobilized ErbB3. These data are thus questionable. No information on binding for the intact antibody is referenced or provided. The VHH is clearly folded. The question is whether it really binds to ErbB3. 

The structure determination is technically sound.

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Partly

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Partly

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.

Eliseev Igor St. Petersburg National Research Academic University RAS, Russian Federation

Competing interests: No competing interests were disclosed.

 3 7 2018

Dear Prof. Ferguson,

We appreciate your comments on the SPR experiment, which appeared both insightful and very helpful for us.

We first addressed the regeneration issue. The initial choice of regeneration buffer (glycine-HCl pH 2.0) was made because it efficiently eluted several anti-ErbB3 VHH fragments we were working with. We probed the receptor with a different antibody and saw no significant loss of binding after treatment with glycine buffer. However, this alternative antibody targeted a different epitope on the receptor, while the epitope of BCD090-M2 appeared very sensitive to low pH, as you suggested. We tested various regeneration conditions, and although their effect on the receptor was not as detrimental as that of glycine, we still saw some memory effects. Therefore, we finally adopted a protocol with a long dissociation phase (30 min) followed by a 10 min equilibration with running buffer, which was sufficient for complete elution of the antibody.

Then we increased the length of the injection to 30 min to ensure that all binding curves reach equilibrium at the end of association phase. The measured binding response as a function of antibody concentration was fitted with simple 1:1 binding model by a nonlinear least-squares method which gave equilibrium constant K D=14 nM. We repeated the experiment twice and saw a remarkable reproducibility; the mean K D was 15±1 nM. Finally, we updated the Figures 2A and 2B to present new experimental data.

 10.5256/f1000research.14786.r30112 Reviewer response for version 1 Hu Shi 1 Referee Department of Biophysics, College of Basic Medical Sciences, Second Military Medical University, Shanghai, China

Competing interests: No competing interests were disclosed.

 5 3 2018 Copyright: © 2018 Hu S 2018 This is an open access peer review report distributed under the terms of the Creative Commons Attribution Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. recommendation approve

For the immunoglobulin-γ (IgG) antibodies, the structure usually is assembled from two identical heavy (H)-chain and two identical light (L)-chain polypeptides is well established and highly conserved in mammals. The L chain of these immunoglobulins comprises two domains, whereas the H chain folds into four domains. Currently, most of therapeutic antibodies used in clinical development are belong to such type. One exception to this conventional mammalian IgG structure is found in sera of Camelidae. Unlike the conventional heterotetrameric antibodies these sera possess special IgG antibodies. IgG antibodies, known as heavy-chain antibodies (HCAbs), are devoid of the L chain polypeptide and are unique because they lack the first constant domain (CH1). At its N-terminal region, the H chain of the homodimeric protein contains a dedicated variable domain, referred to as VHH, which serves to associate with its cognate antigen. The VHH in an HCAb is the structural and functional equivalent of the Fab fragment (antigen-binding fragment) of conventional antibodies.

In this report, Igor et al expressed a recombinant single-domain antibody (VHH form llamas) targeting ErbB3 Protein crystallography study was used for determine the structure of VHH in different crystal forms. It is a well done crystallography report providing readers more data for the structure of VHH, and I only have a few minor suggestions:

 Do the authors have any information or data about the epitope of BCD090-M2?

Single-domain antibodies from llamas are new generation of therapeutical antibody drugs, can authors make another 3-D Alignment figure to compare BCD090-M2 with human antibodies?

Is the work clearly and accurately presented and does it cite the current literature?

Yes

If applicable, is the statistical analysis and its interpretation appropriate?

Not applicable

Are all the source data underlying the results available to ensure full reproducibility?

Yes

Is the study design appropriate and is the work technically sound?

Yes

Are the conclusions drawn adequately supported by the results?

Yes

Are sufficient details of methods and analysis provided to allow replication by others?

Yes

Reviewer Expertise:

NA

I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard.

Eliseev Igor St. Petersburg National Research Academic University RAS, Russian Federation

Competing interests: No competing interests were disclosed.

 3 7 2018

Dear Dr. Shi,

We thank you for reading and commenting on our manuscript. Let us respond to your comments:

Initially, we hypothesized that the antibody BCD090-M2 binds to the extracellular domain III of ErbB3 (ECDIII). To test this hypothesis, we have conducted pairwise SPR binding experiments with intact ErbB3 ECD and another anti-ErbB3 antibody which targets ECDIII but saw simultaneous binding of two antibodies. Then we expressed and purified ECDIII (amino acids 329-532 of ErbB3), and also observed no interaction with BCD090-M2 in direct binding experiments. We are now performing additional experiments to determine the epitope of BCD090-M2 on the ErbB3 extracellular domain and elucidate the binding mechanism.

The alignment of the framework regions of our llama antibody and heavy chains of human antibodies is indeed very good. However, the CDRs of our antibody does not align well with human antibodies, e.g., those anti-ErbB3 antibodies we cited in the manuscript (target ECD I, III, II and IV). It is interesting whether llama and human antibodies recognize the same molecular features, and we certainly plan to compare our llama antibody with human anti-ErbB3 antibodies after identification of its epitope.