Improved serological detection of rheumatoid arthritis: a highly antigenic mimotope of carbonic anhydrase III selected in a murine model by phage display.

INTRODUCTION: Rheumatoid arthritis (RA) is a chronic inflammatory autoimmune disease that affects around 1% of the human population worldwide. RA diagnosis can be difficult as there is no definitive test for its detection. Therefore, the aim of this study was to identify biomarkers that could be used for RA diagnosis.
METHODS: Sera from a collagen-induced arthritis mouse model were used to select potential biomarkers for RA diagnosis by phage display technology. In silico and in vitro analyses were performed to characterize and validate the selected peptides. Samples were classified into three groups: RA; two other immune-mediated rheumatic diseases (systemic lupus erythematosus (SLE) and ankylosing spondylitis (AS)); and healthy controls (HC). Enzyme-linked immunosorbent assay (ELISA) was carried out to determine antibody levels, and diagnostic parameters were determined by constructing receiver operating characteristic curves. Mass spectrometry and Western blot were performed to identify the putative autoantigen that was mimicked by a highly reactive mimotope.
RESULTS: After three rounds of selection, 14 clones were obtained and tested for immunoreactivity analysis against sera from RA and HC groups. The phage-fused peptide with the highest immunoreactivity (M12) was synthesized, and was able to efficiently discriminate RA patients from SLE, AS and HCs (p < 0.0001) by ELISA. The specificity and sensitivity of anti-M12 antibodies for RA diagnosis were 91 % and 84.3 %, respectively. The M12 peptide was identified as one that mimics a predicted antigenic site of the carbonic anhydrase III (CAIII) protein, a ubiquitous biomarker that has been identified in patients with other diseases.
CONCLUSION: M12 is the first peptide associated with the CAIII protein that may be used as an antigen for antibody detection to aid in RA diagnosis with high sensitivity and specificity.

Introduction

Rheumatoid arthritis (RA), the most common inflammatory autoimmune disease, affects 0.8 % of the adult population worldwide [1]. RA diagnosis is largely a clinical one, relying, particularly in the early stages, on the history and examination of the patient, with tests (blood or imaging) sometimes helping to confirm the diagnosis [2]. Serological support to diagnosis has, up to now, been restricted to the determination of rheumatoid factor (RF) and anticitrullinated peptide antibodies (ACPAs), where assays using cyclic citrullinated peptides (CCP) as antigen for ACPA detection have gained wide acceptance [3]. RF presents higher sensitivity as compared with antiCCP antibodies for established disease, with a relatively low specificity. In fact, the RF antibody is not specific for RA due to crossreactivity with many other inflammatory diseases, as well as in elderly healthy individuals [4]. ACPAs are considered a valuable serological biomarker for RA [5] and the diagnostic performance of different generations of CCPs (CCP1, CCP2 and CCP3) have been evaluated in many different studies [6-8]. Differences in cut-off values, specificities and sensitivities exist between the three different generations and also between different assays used for antibody detection. However, antiCCP2 showed better performance characteristics with values of sensitivity ranging from 41 % to 92.2 % and specificity ranging from 65 % to 100 % [9]. At present, the detection of antibodies against CCP2 by enzyme-linked immunosorbent assay (ELISA) is the most widely used assay in studies involving ACPAs worldwide. The combination of RF and antiCCP2 assays demonstrate a positive predictive value close to 100 %, which is much higher than the value of either of the tests alone [10]. The presence of RF and antiCCP has been associated with progressive and destructive disease [11, 12]. Seronegativity in both early and established RA remains a major limitation of these two biomarkers, highlighting the need for new complementary markers that could improve diagnostic sensitivity [13]. Because of the low sensitivity or specificity of the current serological tests, the quest for new efficient auxiliary biomarkers in RA is of clinical relevance.

Animal models of arthritis have contributed to the overall knowledge on RA physiopathology and to the identification of important mediators of inflammation. The collagen-induced arthritis (CIA) mouse model has proven to be a valuable experimental model for inflammatory RA studies [14-17]. After immunization with type II collagen (CII), DBA/1 J mice develop a severe polyarthritis mediated by an autoimmune response that shares many features with human RA [18].

With the goal of identifying new clinically useful biomarkers for RA, we have explored the CIA mouse model and phage display (PD) technology to isolate peptides that can mimic RA autoantigens. PD technology has been widely used by our group and others to screen targeting peptides in drug discovery and biomarker selection, and has been highly effective in discovering peptides with affinities to virtually any target [19-24]. Short peptide sequences selected by PD libraries with high affinity to antibodies, receptors or proteins may present potential applications in diagnostics or therapeutic kits and vaccines [25, 26]. Using the cDNA PD library for autoantigen selection, were recently identified novel autoantibodies in early and seronegative RA patients with sensitivities ranging from 2 % to 29 %, and specificities ranging from 95 % to 100 %. These autoantibodies can be found in 44 % to 67 % of ACPA-negative RA patients [27].

This investigation describes the identification of a short peptide selected by PD technology against sera from CIA mice. This short peptide was characterized by in silico and in vitro strategies, and further tested as a potential biomarker for RA diagnosis in comparison to other rheumatic diseases and healthy controls. Its predicted antigenic epitope target was deduced by mass spectrometry (MS) and Western blot analyses.

Methods

Study subjects

For this study, we used serum samples from patients who fulfilled the RA diagnostic criteria according to the 2010 classification [28]. The demographic and laboratory characteristics of all studied subjects are presented in Table 1. Well-characterized serum samples used in this study were requested from Biobanco-IMM, Lisbon Academic Medical Center, Lisbon, Portugal. Samples from healthy subjects who did not have any arthritis symptoms were used as the healthy control (HC) group. To provide data on assay specificity, samples from patients with other rheumatic diseases (nonRA) were also analyzed. The subjects enrolled in this study were classified into three groups: RA, 172 patients (151 females and 21 males, mean age 53.7 ± 11.2 years); HC, 113 subjects (79 females and 34 males, mean age 58.8 ± 8.7 years); and nonRA, 32 patients with other immune-mediated rheumatic diseases that consisted of 19 systemic lupus erythematosus (SLE) and 13 ankylosing spondylitis (AS) patients. For RA patients, data recorded were age at disease diagnosis, disease duration, medication status, presence of RF and antiCCP, erythrocyte sedimentation rate (ESR), levels of C-reactive protein (CRP), and tender and swollen joint counts obtained at the time of blood sample collection. Blood samples were allowed to clot and then centrifuged at 4,000 rpm for 10 minutes. All sera samples were filtrated with a microcell filter (0.22 μm) to eliminate red blood cell fragments and bacteria, and then frozen at −80 °C immediately until analysis. The parts of the study which involved human subjects was approved by the ethics committees of the Centro Hospitalar Lisboa Norte, Hospital de Santa Maria and the Hospital Garcia de Orta, Lisbon, Portugal, and conducted in accordance with the regulations governing clinical trials such as the Declaration of Helsinki (2008). An informed consent form was signed by all subjects enrolled in this study before any protocol procedure was carried out.

Table 1

Demographic and laboratory characteristics of the studied population

	RA	HC	SLE	AS
Age (years; mean ± SD)	53.7 ± 11.2	58.8 ± 8.7	55.7 ± 5.7	49.5 ± 7.5
Gender (female/male)	151/21	79/34	16/3	4/9
Caucasian (n (%))	146 (85 %)	104 (92 %)	18 (95 %)	11 (85 %)
Disease duration (years; mean ± SD)	8.1 ± 8.4		9.2 ± 6.5	12 ± 9.1
Rheumatoid factor positive (n)	102 (59.3 %)		NA	NA
CRP (mg/dl; mean ± SD)	2.4 ± 1.1		3.1 ± 0.9	2.7 ± 0.8
ESR (mm/hour; mean ± SD)	32 ± 26		33.4 ± 19.5	NA
AntiCCP positive (n)	75 (43.6 %)		NA	NA
Biologic therapy (n in treatment)	62 (36 %)		4 (21 %)	3 (23 %)
NSAIDs (n in treatment)	127 (74 %)		9 (47 %)	5 (38 %)

CRP value ≥0.8 mg/dl is considered elevated. ESR value ≥15 mm/hour is considered elevated. antiCCP anticyclic citrullinated peptide antibodies, AS ankylosing spondylitis, CRP C-reactive protein, ESR erythrocyte sedimentation rate, HC healthy control, NA not available, NSAIDs nonsteroidal anti-inflammatory drugs, RA rheumatoid arthritis, SD standard deviation, SLE systemic lupus erythematosus

CIA mice

A total of twelve male DBA/1 J mice (8–10 weeks old) were housed at the animal facility of the Federal University of Uberlândia (Uberlândia-MG, Brazil) where they fed ad libitum. This study was carried out in strict accordance with the recommendations contained in terms for the use of animals in research and teaching of the Federal University of Uberlândia (Uberlândia, MG, Brazil), in compliance with the National Guidelines, as set forth by the Institutional Animal Care (Law number 11.794, 2008). The study protocol was approved by the Animal Ethics Committee from the Federal University of Uberlândia (UFU) under the number CEUA/UFU 059/10.

Arthritis induction and assessment of CIA mice

Arthritis was induced in the isogenic murine model DBA/1 J (males; 8–10 weeks old) as described previously with some modifications [18]. Briefly, a commercial bovine CII (Becton Sigma-Aldrich, St. Louis, MO, USA) emulsified in complete Freund’s adjuvant (CFA) (Sigma-Aldrich) was used for mice immunization. DBA/1 J mice (n = 9) were immunized by intradermal injections at the base of the tail with 100 μg bovine CII. For the control group, mice (n = 3) were immunized with phosphate-buffered saline (PBS) in CFA. Mice were boosted at day 30 with 50 μg bovine CII in CFA. After immunization, mice were examined twice a week when paw edema was measured using manual calipers and arthritis signs were visually scored by evaluating joint inflammation, using an established scoring system from 0 to 4 [29], where 0 = no evidence of erythema and swelling, 1 = erythema and mild swelling confined to the tarsals or ankle joints, 2 = erythema and mild swelling extending from the ankle to the tarsals, 3 = erythema and moderate swelling extending from the ankle to metatarsal joints, 4 = erythema and severe swelling encompass the ankle, foot and digits, or ankylosis of the limb. Around 60 days after immunization, mice presenting acute arthritis (score 4) were sacrificed under anesthesia. Blood (for sera extraction) and inflamed joints were then extracted for peptide selection through PD technology.

Total protein extraction from inflamed joints of CIA mice

Inflamed joints were removed from CIA mice, macerated with liquid nitrogen and suspended in extraction buffer (20 mM Tris–HCl pH 7.2, 10 mM EDTA, 2 mM EGTA, 250 mM sucrose, 1 mM DTT, 1 mM Benzamidine, 1 mM PMSF). The resulting material was transferred to a microtube, and centrifuged at 20,000 × g for 30 minutes at 4 °C. The supernatant was collected and the concentration of the extracted proteins was determined by the Bradford method [30]. As many inflammatory cytokines indicative of RA are expressed in the cartilage, ligament, pannus, articular capsule, and so forth, of CIA mice [31], we used these proteins to dissociate the selected phage clones from target antibodies by competitive elution.

Purification and isolation of immunoglobulin G from sera

Immunoglobulin G (IgG) purification and isolation from sera of mice with acute arthritis and naïve mice was performed with Dynabeads® Protein G ( Invitrogen, Carlsbad, CA, USA), following the manufacturer’s instructions. Briefly, a total of 50 μl magnetic beads were washed twice with TBS-T 0.1 % (Tris-buffered saline: 50 mM Tris–HCl, 150 mM NaCl, pH 7.5 plus 0.1 % Tween 20) and then incubated with 100 μl pooled serum for 1 hour at room temperature. After binding, the bead–IgG complex was blocked with TBS plus 5 % bovine serum albumin (BSA) at 37 °C for 1 hour, washed three times with TBS-T 0.1 % and then resuspended in 200 μl TBS.

Mimotope selection through phage display

For mimotope (phages expressing peptides on their surface) selection a PhD-12mer phage display peptide library kit (New England Biolabs, Beverly, MA, USA) was screened against IgG purified from CIA mice. This is a combinatorial library of random dodecapeptides fused to the N-terminus of the minor coat protein (pIII) of M13 phages. The library consists of 2.7 × 109 diverse sequences that were amplified once to yield about 50 copies of each peptide sequence.

Based on scientific evidence showing that CIA mice and RA patients share several pathological features, circulating autoantibodies to common targets [29], and also because arthritis can be much easier to monitor in experimental animal models than in humans due to the complexity of the symptoms, we decided to perform the mimotope selection against IgG purified from CIA mice.

The strategy adopted for the mimotope selection consisted of a subtractive step to remove nonspecific phages by pre-incubating the phage peptide library with IgG purified from naïve mice. In order to remove nonspecific phages, a volume of 7 × 108 beads/ml coupled with IgG purified from naïve mice serum was incubated with 1 × 1011 phage particles from the PhD-12 library in 200 μl TBS-T 0.1 % solution for 30 minutes at room temperature. After paramagnetic bead precipitation using the Magnetic Particle Concentrator (Dynal MPC™; Invitrogen), unbound phages were collected and incubated with beads coupled with IgG purified from CIA mice presenting signs of acute arthritis, following incubation for 30 minutes at room temperature. The unbound phages were discarded this time by washing ten times with TBS-T 0.1 %. For competitive elution, the bound phages were incubated for 30 minutes with 10 μg of the total protein extracted from inflamed joints of CIA mice. The eluted phages were amplified in Escherichia coli ER2738 strain (New England Biolabs) and purified by PEG-NaCl precipitation. After each of the three rounds of selection, individual bacterial colonies containing amplified phage clones were grown in a microtiter plate and titrated as described elsewhere [32].

DNA extraction and sequencing

Phage DNA was isolated from 1 ml overnight cultures by precipitation with 1/6 volume PEG/NaI (20 % w/w, polyethylene glycol 8000) and iodide buffer (10 mm Tris–HCl (pH 8.0), 1 mm EDTA, and 4 m NaI). Phage DNA was precipitated with absolute ethanol, followed by a wash with 70 % ethanol, and resuspended in 20 μl Milli-Q water. Electrophoresis was performed on 0.8 % agarose gel stained with ethidium bromide solution in order to verify DNA quality. Sequencing reactions were carried out using the DyEnamic ET Dye Terminator Cycle Sequencing Kit (GE Healthcare, Pittsburg, PA, USA), with the primer −96 M13 (5’-OH CCC TCA TAG TTA GCG TAA CG-3’) following the manufacturer’s instructions, and detection was performed in a MegaBace 1000 Genetic Analyzer (Amersham Biosciences, Little Chalfont, UK) automatic capillary sequencer.

Bioinformatic analysis

A tool that can be found on the Sequence Manipulation Suite collection of JavaScript programs [33] was utilized to obtain the reverse complementary sequences of the DNA extracted from phages. Amino acid sequences were deducted by ExPASy Proteomics and Sequence Analysis tool [34, 35]. For sequence similarity analysis, multiple sequence alignment was performed by ClustalW2 online server [36]. Three-dimensional structure prediction was performed using The Pepitope Server [37, 38], and Immune Epitope Database and Analysis Resource [39]. The PyMOL (available at [40]) plugin was used for showing the peptide surface. Antigenicity prediction was carried out using Kolaskar and Tongaonkar antigenicity scale [41, 42].

Immunoreactivity of the selected mimotopes by phage-ELISA

For immunoreactivity measurements of selected mimotopes against sera from RA and HC groups, a phage-ELISA assay was performed. A 96-well Maxisorp™ microtiter plate (NUNC, New York, NY, USA) was coated in triplicate with anti-M13 monoclonal antibody (Amersham Biosciences) diluted (1:100) in carbonate buffer (0.1 M NaHCO3, pH 8.6) overnight at 4 °C. The plate was washed once with TBS-T 0.5 % and blocked for 1 hour at 37 °C with 5 % BSA diluted in TBS. Additionally, the plate was washed twice and incubated with culture supernatant containing amplified phage particles (~ 1011 pfu/μl) for 1 hour at 37 °C. The plate was washed three times followed by incubation with serum pools from the RA and HC groups diluted (1:100) in TBS-T 0.5 % plus 5 % BSA for 1 hour at 37 °C. The plate was washed three times more with TBS-T 0.5 % followed by incubation with HRP-conjugated rabbit anti-human IgG (Roche Applied Science, Indianapolis, IN, USA.) diluted (1:5,000) in TBS-T 0.5 % plus 5 % BSA for 1 hour at 37 °C. The ELISA plate was washed three times, revealed with OPD SigmaFast™ (Sigma-Aldrich) and read at 492 nm. The reactivity obtained by the wild-type M13 phage without displaying any peptide (performed for each sample tested) was used for data adjustment, where the final optical density (OD) values obtained for each mimotope were adjusted by subtracting the corresponding OD values obtained by the wild-type M13 phage.

Purification of human anti-M12 antibody and M12 mimotope by immunoprecipitation

For purification of IgG antibodies that bind to the M12 mimotope, 1 × 1011 of M12 phage particles were covalently bound to Dynabeads® (Invitrogen). Thereafter, a volume of 50 μl of the solution containing the bead–M12 complex was separately incubated with 100 μl pooled sera from RA patients and HCs for 1 hour at room temperature under shaking for IgG:M12 binding. After incubation, the complex was washed ten times with TBS-T 0.1 % and the unbound nonspecific IgG present in the supernatant was discarded. The IgG that bound with high affinity to the M12 phage clone were eluted with 100 μl elution buffer (0.2 M Glycine-HCl, pH 2.2 and BSA 1 mg/ml) after incubation for 10 minutes at room temperature, followed by neutralization with 15 μl 1 M Tris–HCl (pH 9.1). The eluted IgG was incubated with 50 μl magnetic beads, as previous described. Subsequently, the solution containing the bead–IgG complex was incubated with 1 μg total protein extracted from inflamed joints of CIA mice for 1 hour at room temperature, under shaking condition for protein binding. The unbound proteins were discarded by washing ten times with TBS-T 0.1 %, followed by elution of the bound proteins with 100 μl 0.8 M acetic acid (pH 2.0). Proteins eluted were dried and submitted for MS analysis.

Western blot analysis

Total proteins extracted from inflamed joints of CIA mice (1 μg) were separated by 10 % SDS-PAGE and transferred to a nitrocellulose membrane (GE Healthcare). The membrane was blocked for 1 hour with 3 % BSA in PBS, and then rinsed three times in washing buffer containing PBS-T 0.1 %. Thereafter, the membrane was incubated overnight at 4 °C with IgG purified from sera of RA patients and HCs diluted at 100 μg/ml in PBS plus 3 % BSA, washed three times and incubated for 1 hour with rabbit anti-human IgG conjugated with peroxidase (Roche Applied Science) diluted 1:5,000 in PBS plus 3 % BSA. The subsequent washing steps and detection procedures were performed according to the ECL Plus manual (GE Healthcare).

Mass spectrometry

Protein digestion: dried proteins were suspended in 30 μl 0.2 % RapiGest and vortexed for 5 minutes. Protein digestion was carried out as described elsewhere with a few modifications [43]. Briefly, 10 μl 50 mM ammonium bicarbonate was added to the protein suspension to a final volume of 40 μl. Protein samples were denatured with 25 μl 0.2 % (w:v) RapiGest SF for 15 minutes at 80 °C, reduced with 2.5 μl 100 mM dithiothreitol at 60 °C for 30 minutes, alkylated with 2.5 μl 300 mM iodoacetamide at room temperature, and enzymatically digested at 37 °C overnight with trypsin at a 1:100 (w/w) enzyme to protein ratio. Then, 10 μl 5 % trifluoroacetic acid (TFA) was added to the digestion mixture to hydrolyze the RapiGest, and samples were incubated at 37 °C for 90 minutes. The tryptic peptide solution was then centrifuged at 14,000 rpm for 30 minutes at 6 °C, and the pH of the supernatant was adjusted to 2.6 by the addition of 10 μl 3 % acetonitrile 0.1 % formic acid.

Liquid chromatography/MS analysis: separation of peptides was performed using nanoliquid chromatography employing reverse phase. Peptides were injected into nanoLC through a nanoACQUITY system (Waters, Manchester, UK). Samples were first trapped in a Symmetry C18 5 μm, 180 mm × 20 mm column (Waters) with 0.1 % TFA in 3 % acetonitrile, then peptides were eluted from the trap column to an HSS T3 1.8 μm, 75 μm × 15 cm analytical column (Waters; mobile phase A, water with 0.1 % formic acid and B, 0.1 % formic acid in acetonitrile). Mass spectrometric acquisition was achieved in a Synapt MS Q-TOF mass spectrometer equipped with a NanoLockSpray source in the positive ion mode (Waters). For all measurements, the mass spectrometer was operated in the ‘V’ mode with a typical resolving power of at least 12,500. Data-independent scanning (MSE) experiments were performed by switching between low (3 eV) and elevated collision energies (15–50 eV) applied to the trap ‘T-wave’ cell filled with argon. Scan times of 0.8 s were used for low- and high-energy scans from m/z 50 to 2,000.

Protein identification and database analysis: protein identification was performed in Global Server v.2.5 (PLGS, Waters) with mouse UniProtKB Complete Proteome database. Up to one maximum missed cleavage by trypsin was allowed; fixed modification by carbamidomethylation (cysteine) and variable modifications by acetyl N-terminal and oxidation (methionine) were considered. The precursor and fragment ion mass error tolerances were adjusted to 10 and 20 ppm, respectively (default values). The criteria used for a positive protein match were at least three fragment ions per peptide, seven fragment ions per protein, and at least one peptide per protein hit. A false-positive discovery rate was allowed up to 4 %.

Peptide design and synthesis

After bioinformatic analysis of the selected clones, the M12 peptide sequence was designed and chemically synthesized by GenScript USA Inc. (Piscataway, NJ, USA). The peptide was constructed with 17 residues (CNVNSKSPVERITGGGS), with amidation of the C-terminal and conjugation of the BSA to cysteine at the N-terminus, a design strategy to increase sensitivity and decrease crossreactivity on ELISA assay [44].

Anti-M12 antibody detection by ELISA

A specific ELISA test was carried out to determine the M12 synthetic peptide reactivity to sera from patients with RA, SLE, AS and HCs. Ninety-six-well Maxisorp™ microtiter plates (NUNC) were coated with synthetic M12 (1.5 μg/ml) in carbonate buffer (0.1 M NaHCO3, pH 8.6), and incubated overnight at 4 °C. After blocking with 3 % BSA in PBS at 37 °C for 1 hour, 100 μl/well of sera from the different groups were diluted in blocking buffer (1:100) and incubated for 1 hour at 37 °C under gentle agitation. After four washes with PBS-T 0.05 %, HRP-conjugated rabbit anti-human IgG (Roche Applied Science) diluted (1:5,000) in the blocking buffer was incubated for 1 hour at 37 °C under gentle agitation. The ELISA plate was washed three times with PBS-T 0.05 %, revealed with OPD SigmaFast™ (Sigma-Aldrich) and read at 492 nm. All samples were tested in triplicate. The optimum point of reaction for anti-M12 antibody detection was determined using the receiver operating characteristic (ROC) curve, where a cut-off point was determined as the value of the parameter corresponding to the highest possible sensitivity without losing specificity. To calculate the ROC curve, sensitivity and specificity, we considered the control groups as a single group. Each serum sample was tested without M12 as negative control. The final OD values obtained for each RA, HC, SLE or AS samples were adjusted by subtracting the corresponding OD value obtained by the negative control. After data adjustment, OD values obtained for each sample from all groups were divided by the cut-off value for data normalization. The values obtained are expressed as reactivity index (RI), where samples presenting RI ≥1 were considered positives.

Statistical analysis

Unpaired t test with Welch’s correction was used to evaluate the differences in sera reactivity in ELISA assays among groups for phage clones and the synthetic peptide. Sensitivity and specificity parameters were calculated based on the ROC curve analysis, and Fisher’s exact test was used for categorical data. To estimate the positive predictive accuracy, the area under the curve (AUC) was also determined. Pearson’s correlation was used for analysis among variables. All statistical analyses were performed using GraphPad Prism 5.0 software (GraphPad Software, Inc., San Diego, CA). p values less than 0.05 were considered statistically significant.

Results

Arthritis induction

Arthritis induction was efficient in CIA mice. A total of 80 % of male DBA/1 J mice developed acute arthritis around 60 days after CII immunization. The severity and incidence of arthritis were assessed as described in Methods. No manifestation of arthritis was observed in mice in the control group treated with PBS in CFA. A schematic workflow illustrating the steps employed in this study is shown in Fig. 1.

Fig. 1

Schematic workflow illustrating the steps involved in the study. a Selection step, and b validation step. Arthritis was induced in collagen-induced arthritis (CIA) mice by immunization with type II collagen (CII). Immunoglobulin G (IgG) was purified by protein G beads from serum of mice presenting with acute arthritis and naïve mice. A phage display library was used to select mimotopes against purified IgG from CIA mice (three cycles of biopanning), with a subtractive step against IgG from naïve mice. IgG-bound phages from CIA mice were competitively eluted by incubation with total proteins extracted from inflamed joints. Total proteins extracted from inflamed joints were also used for anti-M12 antibody recognition in Western blot. Immunoreactivity of each selected mimotope was tested by phage exyme-linked immunosorbent assay (ELISA). After DNA extraction and in silico and in vitro analysis, the most reactive mimotope (M12) was identified as a peptide that mimics a predicted antigenic site of the human carbonic anhydrase III protein. Validation of the M12 synthetic peptide as a possible rheumatoid arthritis (RA) autoantigen was carried out by ELISA assay

Selection of mimotopes by phage display

The enrichment of phages recovered after each round of biopanning was determined by the output to input ratio. The increase from 8 × 104 in the first round to 4 × 105 after three rounds of affinity selection showed a clear enrichment of phage particles (Table 2). These data indicate a successful affinity selection of phages that specifically recognized IgG present in sera of DBA1/J mice with acute arthritis. A total of 37 randomly selected mimotopes were obtained after three rounds of biopanning using a phage displayed 12-mer random peptide library. From the 37 mimotopes selected, 14 presented different sequences. Alignment analysis revealed some consensus sequence between the mimotopes selected (Fig. 2a), indicating that these motifs were positively selected during the biopanning. The selected mimotopes showed different reactivities when tested against human sera by phage-ELISA, and all of them were able to discriminate RA patients from HCs. The reactivity values were very similar for all tested mimotopes, except for M12, which showed the highest values of absorbance (Fig. 2b) and the highest difference compared to the other clones (p < 0.05). Due to the higher reactivity compared to the other mimotopes, the common motif shared with another five clones, and the highest ratio of RA:HC, we focused on M12 peptide for further investigation. However, the other mimotopes selected will be explored in future studies.

Table 2

Enrichment of phage for each round of selection from phage display peptide library

	Number of phage particles
Round	Input (cpu)	Output (cpu)	Ratio (output/input)
1st	1 × 1011	8 × 104	8 × 10−7
2nd	1 × 1011	9 × 103	9 × 10−8
3rd	1 × 1011	4 × 105	4 × 10−6

cpu phage units

Fig. 2

Perfomance of the mimotopes selected by phage display. a Multiple sequence alignment with all selected mimotopes showing their consensus sequences and frequency. b Reactivity obtained though the interaction of the mimotopes selected and pooled sera from rheumatoid arthritis (RA) patients and health controls (HC). aFrequency is defined as the ratio of the number of phage clones expressing a common peptide sequence to that of the total phage clones obtained in the biopanning

SDS-PAGE fractionated proteins from inflamed joints of CIA mice were transferred to a nitrocellulose membrane. Western blot revealed that purified IgG from RA patients against the M12 mimotope recognized a protein with a molecular weight of approximately 29 kDa (Fig. 3).

Fig. 3

Western blot analyses. a The Ponceau staining shows 1 μg of total proteins extracted from inflamed joints of CIA mice (§) separated by 10 % SDS-PAGE. b The same membrane stained with Ponceau was probed with anti-M12 antibodies purified from pooled sera of RA patients (†) and HCs (‡). Anti-M12 antibodies displayed strong reactivity with a protein extracted from inflamed joints of CIA mice presenting with mass of approximately 29 kDa (arrow). Carbonic anhydrase III (CAIII) was identified as the protein target that is mimicked by the M12 mimotope. The membrane was cropped to allow differential incubation with anti-M12 antibodies purified from RA patients and HCs

Antigen target and epitope prediction

MS analysis identified 15 putative proteins associated with anti-M12 IgG, where 14 proteins presented few peptide matches, with very low scores and sequence coverage and none of them presented a molecular mass close to 29 kDa. Therefore, based on the molecular weight of the protein target by anti-M12 antibodies observed in the Western blot results (29 kDa), the alignment of putative proteins with the M12 peptide sequence, and the sequencing coverage by MS, we have identified the carbonic anhydrase III (CAIII) [Uniprot/Swiss-Prot:P16015] as the protein target that the M12 peptide mimics. The CAIII protein showed 49.62 % of sequencing coverage, including the peptide sequence, which was aligned by conserved and semi-conserved amino acid (aa) residues between position 24–35 (KGDNQSPIELHT).

The multiple sequence alignment revealed several homologous sequences between M12 peptide and the CAIII protein sequences from mouse and human. Nine (75 %) conserved or semi-conserved aa residues of the M12 peptide sequence matched a domain of the CAIII protein (Fig. 4a). A three-dimensional structural alignment was performed to predict the putative epitope site of the M12 peptide in the CAIII protein structure [PDB:3UYN], which confirmed its surface exposure and corroborated to the possible antibody-binding region in the external sequences of the predicted protein (Fig. 4b). Five (41.6 %) aa residues (SPVET) of the M12 peptide was matched to an antigenic region of the human CAIII protein (Fig. 5).

Fig. 4

Peptide localization. a Multiple alignments of M12 peptide and deduced amino acid sequences of carbonic anhydrase III from Homo sapiens and Mus musculus. Conserved (stars) and semi-conserved residues (gray) are highlighted. Valine (V) amino acid is star marked because it matches a position between M12 peptide and human CAIII protein. b Model of the three-dimensional structure predicted in the PyMol server software for the human carbonic anhydrase III with M12 peptide localization

Fig. 5

Antigenicity predictions using the Kolaskar-Tongaonkar algorithm. Threshold (0.993); average (0.993); minimum antigenicity (0.898); maximum antigenicity (1.097). Regions with antigenic propensity scale above 0.993 are antigenic. Window size and center position were 7 and 4, respectively. *Localization of M12 peptide in the CAIII protein (aa 24–35). Matched aa residues are presented in red.

Evaluation of the M12 peptide as a possible biomarker for RA diagnosis

The M12 synthetic peptide was tested by ELISA with individual sera obtained from 172 RA patients, 113 HCs, 19 SLE patients and 13 AS patients to evaluate its potential as a diagnostic. The cut-off point was determined as the value of the parameter corresponding to the highest sensitivity without lowering the specificity. Of 172 RA patients, 145 (84.3 %) were positive to anti-M12 antibodies. The synthetic molecule was able to efficiently discriminate sera from RA patients and HCs (p < 0.0001), SLE (p < 0.0001) and AS (p < 0.0001) (Fig. 6a). The ROC curve analysis constructed based on the control groups was significant (p < 0.0001; AUC = 0.946). Based on the determined cut-off point, M12 synthetic peptide presented a specificity of 91 % and sensitivity of 84.3 % (Fig. 6b). For humoral immune response evaluation, 55 (32 %) serum samples from RA patients were simultaneously positive for RF, antiCCP and anti-M12 antibodies, while 45 (26.2 %) samples were positive only for anti-M12 and 12 (7 %) were negative for the three antibodies (Table 3).

Fig. 6

Anti-M12 antibody detection by ELISA. Detection of anti-M12 antibodies in sera from rheumatoid arthritis (RA) patients (n = 172), healthy controls (HC) (n = 113), systemic lupus erythematosus (SLE) patients (n = 19) and ankylosing spondylitis (AS) patients (n = 13) and the respectively receiver operating characteristic (ROC) curve. a Sera from the four groups individually tested for their ability to bind to synthetic M12 peptide. Horizontal line = cut-off value; error bars show mean and standard deviation. b ROC curve constructed based on the control groups. The area under the curve (AUC), sensitivity (Se), specificity (Sp) and corresponding p-value are indicated inside the graph

Table 3

Antibodies positivity in the rheumatoid arthritis population

Antibody detection			Positivity
Rheumatoid factor	AntiCCP	Anti-M12	n	%
+	−	−	8	4.6
−	+	−	0	0
−	−	+	45	26.2
+	+	−	7	4
+	−	+	32	18.6
−	+	+	13	7.6
+	+	+	55	32
−	−	−	12	7
Total			172	100

+ positive, − negative, CCP cyclic citrullinated peptide

Correlation of anti-M12 humoral response with clinical variables

The weakly associated variables with anti-M12 levels from Gaussian populations (Pearson) analysis were disease duration (r = 0.1753) and the use of biological therapy (r = −0.2934) (Table 4).

Table 4

Clinical variables associated with detection of anti-M12 antibodies in patients with theumatoid arthritis

Variables	Anti-M12
	Pearson r	p	R squared
Ethnicity	†	†	†
Gender	†	†	†
Age	†	†	†
Disease duration	0.1753	<0.05	0.03072
Rheumatoid factor	†	†	†
ESR	†	†	†
CRP	†	†	†
AntiCCP	†	†	†
Tender joints	†	†	†
Swollen joints	†	†	†
Biologic therapy	−0.2934	<0.005	0.08616
NSAIDs	†	†	†

Correlations were performed from Gaussian populations (Pearson) with a confidence interval of 95 %. †Values not statistically significant. CCP cyclic citrullinated peptide, CRP C-reactive protein, ESR erythrocyte sedimentation rate, NSAID nonsteroidal anti-inflammatory drug

Discussion

In order to identify new biomarkers for the improvement of RA diagnosis, we have used the CIA mouse model and PD technologies to select peptides that could mimic RA autoantigens. This approach allowed us to discover a new peptide with 12 amino acids (M12) that encompasses an antigenic domain of the CAIII protein. The use of the synthetic M12 peptide as an autoantigen in ELISA confirmed the presence of specific anti-M12 antibodies in sera of patients with RA with high specificity and sensitivity.

We have used the DBA1/J mouse as a model for inflammatory RA, where arthritis can be induced in susceptible mouse strains by immunization with CII [45]. DBA/1 J mice are CIA models of particular interest, since this strain is not known to have any immunologic aberrations or other pathologic defects [46]. CIA mice are widely used as an animal model for RA studies [14–17, 47–49]. We have chosen this experimental model due to the fact that genomic similarities between mice and humans are quite high, reaching 98 % [50]; it is much easier to monitor arthritis in mice than in humans, and also it is a unique experimental opportunity to increase our understanding of human arthritis. We confirmed the efficiency of the arthritis induction in our CIA mice strain, when a total of 80 % of DBA/1 J developed signs of arthritis around 60 days after CII immunization. We have used a subtractive PD selection against IgG present in sera of CIA mice that aims to select new RA biomarkers. We have chosen this technology due to its successful applications in biomarker discovery, especially when short random peptides are selected against specific target molecules, such as those selected against IgG present in sera for antigen discovery [51, 52].

The highly reactive peptide M12 was able to efficiently discriminate RA patients from HCs as well as from patients affected by SLE and AS, two other immune-mediated rheumatic diseases. In fact, the M12 peptide was not the most enriched mimotope after the rounds of phage amplification; however, the most frequently selected mimotope will not necessarily be the most reactive – selection is independent, which means that increasing one does not eliminate the effect of the other [53]. Many events during viral infection may contribute to increasing the amplification of specific phages. The biological reasons for abundances and growth advantages of some mimotopes have been discussed in many reports and include events like the phages ability to bind to bacterial pili, and interference with packing or infection [54-56]. The M12 peptide sequence was further characterized by reverse binding using its specific IgG antibodies to capture and identify its antigen target by MS, which led us to the CAIII protein. The CAIII protein is a member of a multigene family that encodes carbonic anhydrase isozymes. It is an abundant muscle protein and plays an important role in facilitated CO2 diffusion and diverse processes involving H+ and -HCO3 transport [57]. The CAIII gene encodes a protein of 260 aa and a molecular weight of approximately 29 kDa, which is well conserved throughout humans and mice. This protein has already been identified as an autoantigen expressed in the synovial membrane of RA patients [58] and its circulating autoantibodies have been found in many diseases [59]. The Kolaskar and Tongaonkar antigenicity scale is a semi-empirical method that makes use of physicochemical properties of aa residues and their frequencies of occurrence in experimentally known segmental epitopes, and is used to predict antigenic determinants on proteins. Using this parameter, this method that can predict antigenic determinants with 75 % accuracy [41] has been used in many studies [60-62]. Amino acid residues of the M12 peptide matched conservative and semi-conservative residues in a predicted antigenic and exposed site of a putative epitope of the human CAIII protein, which probably gave this peptide the ability to be recognized by specific IgG during the mimotope selection. This hypothesis is supported by the fact that semi-conserved aa residues present similar physicochemical properties with the original residue allowing antibody binding [63]. Of interest, a synthetic peptide between 24–54 aa residues of the human CAIII, that encompasses the M12 epitope region, has been commercially used to generate antiCAIII antibodies in rabbits (Product ID: ABIN391954). This fact probably reinforces the importance of the CAIII protein in RA and the potential use of this peptide as an antigen for antibody detection to aid in diagnosis.

A wide range of serum biomarkers has been assessed to improve diagnosis and prognosis of RA. However, only the RF and antiCCP antibodies have gained wide acceptance [10]. The use of recombinant CAIII protein as an antigen has confirmed the presence of specific antiCAIII antibodies in RA sera [58, 59]; however, its sensitivity considering the entire CAIII protein was only 17 %. In fact, when compared with HCs, the specificity of circulating autoantibodies was high (100 %), but reached only medium to low specificity when compared with other autoimmune diseases, such as SLE (67 %) [58]. The crossreactivity and lack of sensitivity of the whole CAIII protein suggests the presence of multiple epitopes sharing common regions with other proteins, as well as the presence of immunodominant epitopes that surpass the response to the M12 critical epitope, which may have generated either false positive results or an insufficient reactivity. On the other hand, the use of the M12 peptide conferred more sensitivity and specificity in the detection of RA, since it was able to discriminate RA patients from HCs, and patients with AS and SLE, with high accuracy.

Several studies have measured the three different generations of antiCCP antibodies in RA patients in comparison with control groups. These antibodies present specificities in the range of 65 % to 100 % (mean of 93.4 %) and sensitivities in the range of 42 % to 92.2 % (mean of 64.3 %) [9, 64–71]. RF testing presents specificities in the range 85 % to 90.29 % and sensitivities in the range 46.26 % to 90 % [72-74]. The diagnostic parameters for the M12 peptide for RA diagnosis presented specificity of 91 % and sensitivity of 84.3 %, which is very close to the mean of values reported for antiCCP and RF antibodies, the two most currently used serum biomarkers for RA. On the other hand, anti-M12 antibodies were detected in 26.2 % of samples that are negative for antiCCP and RF. We speculate that the high percentage of positive samples for anti-M12 antibodies might be due to the different tests used to cut-off calculation, or even the patients disease activity at the time of sample collection or individual response to therapy. In this context, high levels of antiCCP and RF antibodies in RA patients have already been associated with an insufficient response to therapy [75, 76], while anti-M12 antibodies were weakly associated with therapeutic response in our RA cohort, which may help to explain the percentage of positive samples only for M12. It is widely accepted that the identification of RA at early stages and consequently the implementation of effective treatment strategies can significantly improve patient prognosis [77]. Anti-M12 antibodies showed a weak association with disease duration, suggesting its potential use as a biomarker in different stages during the process of RA development.

Conclusions

We have selected and identified a peptide that is capable of detecting specific circulating IgG in the serum of RA patients with high specificity and sensitivity. Although antibodies against synthetic M12 peptide were detected in patients with early and established RA, its potential use in the diagnosis at different stages of the disease remains to be studied during follow-up in a well characterized cohort of patients with early RA, and also its specificity as a serum biomarker should be further studied in a much larger cohort of patients with other inflammatory diseases. This peptide mimics a predicted antigenic region of the human carbonic anhydrase III by linear sequence analysis and could be used as an antigen for detection of specific RA autoantibodies.