Surface salt bridges contribute to the extreme thermal stability of an FN3-like domain from a thermophilic bacterium.

This study uses differential scanning calorimetry, X-ray crystallography, and molecular dynamics simulations to investigate the structural basis for the high thermal stability (melting temperature 97.5°C) of a FN3-like protein domain from thermophilic bacteria Thermoanaerobacter tengcongensis (FN3tt). FN3tt adopts a typical FN3 fold with a three-stranded beta sheet packing against a four-stranded beta sheet. We identified three solvent exposed arginine residues (R23, R25, and R72), which stabilize the protein through salt bridge interactions with glutamic acid residues on adjacent strands. Alanine mutation of the three arginine residues reduced melting temperature by up to 22°C. Crystal structures of the wild type (WT) and a thermally destabilized (∆Tm -19.7°C) triple mutant (R23L/R25T/R72I) were found to be nearly identical, suggesting that the destabilization is due to interactions of the arginine residues. Molecular dynamics simulations showed that the salt bridge interactions in the WT were stable and provided a dynamical explanation for the cooperativity observed between R23 and R25 based on calorimetry measurements. In addition, folding free energy changes computed using free energy perturbation molecular dynamics simulations showed high correlation with melting temperature changes. This work is another example of surface salt bridges contributing to the enhanced thermal stability of thermophilic proteins. The molecular dynamics simulation methods employed in this study may be broadly useful for in silico surface charge engineering of proteins.

Using Both Plasma and Urine Donor-Derived Cell-Free DNA to Identify Various Renal Allograft Injuries.

BACKGROUND: This study was designed to investigate the association between donor-derived cell-free DNA (dd-cfDNA) and renal allograft injuries. METHODS: This single-center study enrolled 113 adult kidney transplant recipients with kidney biopsies. Plasma and urine dd-cfDNA was detected by target region capture sequencing. RESULTS: Plasma dd-cfDNA fraction was increased in multiple types of injuries, but most significantly in antibody-mediated rejection. Plasma dd-cfDNA fraction in isolated antibody-mediated rejection (1.94%, IQR: 1.15%, 2.33%) was higher than in T cell-mediated rejection (0.55%, IQR: 0.50%, 0.73%, P = 0.002) and negative biopsies (0.58%, IQR: 0.42%, 0.78%, P < 0.001), but lower than in mixed rejection (2.49%, IQR: 1.16%, 4.90%, P = 0.342). Increased urine dd-cfDNA concentration was associated with several types of injury, but most significantly with BK polyomavirus-associated nephropathy. Urine dd-cfDNA concentration in BK polyomavirus-associated nephropathy (12.22 ng/mL, IQR: 6.53 ng/mL, 31.66 ng/mL) was respectively higher than that in T cell-mediated rejection (5.24 ng/mL, IQR: 3.22 ng/mL, 6.99 ng/mL, P = 0.001), borderline change (3.93 ng/mL, IQR: 2.45 ng/mL, 6.30 ng/mL, P < 0.001), and negative biopsies (3.09 ng/mL, IQR: 1.94 ng/mL, 5.05 ng/mL, P < 0.001). Plasma dd-cfDNA fraction was positively associated with glomerulitis (r = 0.365, P < 0.001) and peri-tubular capillaritis (r = 0.344, P < 0.001), while urine dd-cfDNA concentration correlated with tubulitis (r = 0.302, P = 0.002). CONCLUSIONS: Both plasma and urine dd-cfDNA are sensitive markers for renal allograft injuries. The interpretation of a specific disease by dd-cfDNA should be combined with other clinical indicators.

Crystal structure of B-cell co-receptor CD19 in complex with antibody B43 reveals an unexpected fold.

CD19 is a transmembrane protein expressed on malignant B cells, but not in other lineages or other tissues, which makes it an attractive target for monoclonal antibody-mediated immunotherapy. Anti-CD19 antibody B43 was utilized in a bispecific T-cell engager (BiTE) blinatumomab that demonstrated potency for the treatment of relapsed acute lymphoblastic leukemia. To gain insight into the mechanism of action of the antibody, the crystal structure of B43 Fab was determined in complex with CD19 and in the unbound form. The structure revealed the binding epitope, explained the lack of cross-reactivity toward non-human species, and suggested the key-and-lock mechanism of antigen recognition. Most unexpectedly, the structure revealed a unique molecular topology of CD19. Rather than a tandem of c-type immunoglobulin folds predicted from the amino acid sequence, the extracellular domain of CD19 exhibits an elongated ß-sandwich formed by two immunoglobulin folds by swapping their C-terminal halves. This is the first structure of CD19, which has no sequence homologs.

Epitope mapping and structural basis for the recognition of phosphorylated tau by the anti-tau antibody AT8.

Microtubule-associated protein tau becomes abnormally phosphorylated in Alzheimer's disease and other tauopathies and forms aggregates of paired helical filaments (PHF-tau). AT8 is a PHF-tau-specific monoclonal antibody that is a commonly used marker of neuropathology because of its recognition of abnormally phosphorylated tau. Previous reports described the AT8 epitope to include pS202/pT205. Our studies support and extend previous findings by also identifying pS208 as part of the binding epitope. We characterized the phosphoepitope of AT8 through both peptide binding studies and costructures with phosphopeptides. From the cocrystal structure of AT8 Fab with the diphosphorylated (pS202/pT205) peptide, it appeared that an additional phosphorylation at S208 would also be accommodated by AT8. Phosphopeptide binding studies showed that AT8 bound to the triply phosphorylated tau peptide (pS202/pT205/pS208) 30-fold stronger than to the pS202/pT205 peptide, supporting the role of pS208 in AT8 recognition. We also show that the binding kinetics of the triply phosphorylated peptide pS202/pT205/pS208 was remarkably similar to that of PHF-tau. The costructure of AT8 Fab with a pS202/pT205/pS208 peptide shows that the interaction interface involves all six CDRs and tau residues 202-209. All three phosphorylation sites are recognized by AT8, with pT205 acting as the anchor. Crystallization of the Fab/peptide complex under acidic conditions shows that CDR-L2 is prone to unfolding and precludes peptide binding, and may suggest a general instability in the antibody.

Induced conformational change in human IL-4 upon binding of a signal-neutralizing DARPin.

The crystal structure of DARPin 44C12V5 that neutralizes IL-4 signaling has been determined alone and bound to human IL-4. A significant conformational change occurs in the IL-4 upon DARPin binding. The DARPin binds to the face of IL-4 formed by the A and C α-helices. The structure of the DARPin remains virtually unchanged. The conformational changes in IL-4 include a reorientation of the C-helix Trp91 side chain and repositioning of CD-loop residue Leu96. Both side chains move by >9 Å, becoming buried in the central hydrophobic region of the IL-4:DARPin interface. This hydrophobic region is surrounded by a ring of hydrophilic interactions comprised of hydrogen bonds and salt bridges and represents a classical "hotspot." The structures also reveal how the DARPin neutralizes IL-4 signaling. Comparing the IL-4:DARPin complex structure with the structures of IL-4 bound to its receptors (Hage et al., Cell 1999; 97, 271-281; La Porte et al., Cell 2008, 132, 259-272), it is found that the DARPin binds to the same IL-4 face that interacts with the junction of the D1 and D2 domains of the IL-4Rα receptors. Signaling is blocked since IL-4 cannot bind to this receptor, which it must do first before initiating a productive receptor complex with either the IL-13α1 or the γc receptor.

Structural evidence for a constrained conformation of short CDR-L3 in antibodies.

Three Fab structures used as targets in the Antibody Modeling Assessment presented a challenge for modeling CDR-L3 due to deviations from the most typical canonical structure. In all three antibodies CDR-L3 has eight residues, which is one residue shorter than usual, and has a conformation that is rarely observed in crystal structures. We analyzed the sequence and structural determinants of this conformation and found that the "short" CDR-L3 is remarkably rigid and retains the conformation in the interactions with antigens and neighboring CDRs.

C-terminal ß-strand swapping in a consensus-derived fibronectin Type III scaffold.

The crystal structures of six different fibronectin Type III consensus-derived Tencon domains, whose solution properties exhibit no, to various degrees of, aggregation according to SEC, have been determined. The structures of the five variants showing aggregation reveal 3D domain swapped dimers. In all five cases, the swapping involves the C-terminal ß-strand resulting in the formation of Tencon dimers in which the target-binding surface is blocked. All of the variants differ in sequence in the FG loop, which is the hinge loop in the ß-strand-swapped dimers. The six tencon variants have between 0 and 5 residues inserted between positions 77 and 78 in the FG loop. Analysis of the structures suggests that a non-glycine residue at position 77 and insertions of <4 residues may destabilize the ß-turn in the FG loop promoting ß-strand swapping. Swapped dimers with an odd number of inserted residues may be less stable, particularly if they contain proline residues, because they cannot form perfect ß-bridges in the FG regions that link the swapped dimers. The Tencon ß-swapped variants with the longest FG sequences are observed to form higher order hexameric or helical oligomeric structures in the crystal correlating well with the aggregation properties of these domains observed in solution. Understanding the structural basis for domain-swapped dimerization and oligomerization will support engineering efforts of the Tencon domain to produce variants with desired biophysical properties.

Structure and specificity of an antibody targeting a proteolytically cleaved IgG hinge.

The functional role of human antihinge (HAH) autoantibodies in normal health and disease remains elusive, but recent evidence supports their role in the host response to IgG cleavage by proteases that are prevalent in certain disorders. Characterization and potential exploitation of these HAH antibodies has been hindered by the absence of monoclonal reagents. 2095-2 is a rabbit monoclonal antibody targeting the IdeS-cleaved hinge of human IgG1. We have determined the crystal structure of the Fab of 2095-2 and its complex with a hinge analog peptide. The antibody is selective for the C-terminally cleaved hinge ending in G236 and this interaction involves an uncommon disulfide in VL CDR3. We probed the importance of the disulfide in VL CDR3 through engineering variants. We identified one variant, QAA, which does not require the disulfide for biological activity or peptide binding. The structure of this variant offers a starting point for further engineering of 2095-2 with the same specificity, but lacking the potential manufacturing liability of an additional disulfide. Proteins 2014; 82:1656-1667. © 2014 Wiley Periodicals, Inc.

Second antibody modeling assessment (AMA-II).

To assess the state of the art in antibody 3D modeling, 11 unpublished high-resolution x-ray Fab crystal structures from diverse species and covering a wide range of antigen-binding site conformations were used as a benchmark to compare Fv models generated by seven structure prediction methodologies. The participants included: Accerlys Inc, Chemical Computer Group (CCG), Schrodinger, Jeff Gray's lab at John Hopkins University, Macromoltek, Astellas Pharma/Osaka University and Prediction of ImmunoGlobulin Structure (PIGS). The sequences of benchmark structures were submitted to the modelers and PIGS, and a set of models were generated for each structure. We provide here an overview of the organization, participants and main results of this second antibody modeling assessment (AMA-II). Also, we compare the results with the first antibody assessment published in this journal (Almagro et al., 2011;79:3050).

Antibody modeling assessment II. Structures and models.

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

N-terminal ß-strand swapping in a consensus-derived alternative scaffold driven by stabilizing hydrophobic interactions.

The crystal structure of an N-terminal ß-strand-swapped consensus-derived tenascin FN3 alternative scaffold has been determined. A comparison with the unswapped structure reveals that the side chain of residue F88 orients differently and packs more tightly with the hydrophobic core of the domain. Dimer formation also results in the burial of a hydrophobic patch on the surface of the domain. Thus, it appears that tighter packing of F88 in the hydrophobic core and burial of surface hydrophobicity provide the driving forces for the N-terminal ß-strand swapping, leading to the formation of a stable compact dimer.

Molecular determinants of ASIC1 modulation by divalent cations.

Acid-sensing ion channels (ASICs) are proton-gated cation channels widely expressed in the nervous system. ASIC gating is modulated by divalent cations as well as small molecules; however, the molecular determinants of gating modulation by divalent cations are not well understood. Previously, we identified two small molecules that bind to ASIC1a at a novel site in the acidic pocket and modulate ASIC1 gating in a manner broadly resembling divalent cations, raising the possibility that these small molecules may help to illuminate the molecular determinants of gating modulation by divalent cations. Here, we examined how these two groups of modulators might interact as well as mutational effects on ASIC1a gating and its modulation by divalent cations. Our results indicate that binding of divalent cations to an acidic pocket site plays a key role in gating modulation of the channel.

Concurrent JCPyV-DNAemia Is Correlated With Poor Graft Outcome in Kidney Transplant Recipients With Polyomavirus-associated Nephropathy.

BACKGROUND: Co-infection of JC polyomavirus (JCPyV) and BK polyomavirus (BKPyV) is uncommon in kidney transplant recipients, and the prognosis is unclear. This study aimed to investigate the effect of concurrent JCPyV-DNAemia on graft outcomes in BKPyV-infected kidney transplant recipients with polyomavirus-associated nephropathy (PyVAN). METHODS: A total of 140 kidney transplant recipients with BKPyV replication and PyVAN, 122 without concurrent JCPyV-DNAemia and 18 with JCPyV-DNAemia were included in the analysis. Least absolute shrinkage and selection operator regression analysis and multivariate Cox regression analysis were used to identify prognostic factors for graft survival. A nomogram for predicting graft survival was created and evaluated. RESULTS: The median tubulitis score in the JCPyV-DNAemia-positive group was higher than in JCPyV-DNAemia-negative group ( P  = 0.048). At last follow-up, the graft loss rate in the JCPyV-DNAemia-positive group was higher than in the JCPyV-DNAemia-negative group (50% versus 25.4%; P  = 0.031). Kaplan-Meier analysis showed that the graft survival rate in the JCPyV-DNAemia-positive group was lower than in the JCPyV-DNAemia-negative group ( P  = 0.003). Least absolute shrinkage and selection operator regression and multivariate Cox regression analysis demonstrated that concurrent JCPyV-DNAemia was an independent risk factor for graft survival (hazard ratio = 4.808; 95% confidence interval: 2.096-11.03; P  < 0.001). The nomogram displayed favorable discrimination (C-index = 0.839), concordance, and clinical applicability in predicting graft survival. CONCLUSIONS: Concurrent JCPyV-DNAemia is associated with a worse graft outcome in BKPyV-infected kidney transplant recipients with PyVAN.

Upregulation of BMP1 through ncRNAs correlates with adverse outcomes and immune infiltration in clear cell renal cell carcinoma.

BACKGROUND: Renal cell carcinoma (RCC) accounts for approximately 2-3% of all adult malignancies. Clear cell renal cell carcinoma (ccRCC), which comprises 70-80% of all RCC cases, is the most common histological subtype. METHODS: ccRCC transcriptome data and clinical information were downloaded from the TCGA database. We used the TCGA and GEPIA databases to analyze relative expression of BMP1 in various types of human cancer. GEPIA was used to perform survival analysis for BMP1 in various cancer types. Upstream binding miRNAs of BMP1 were obtained through several important target gene prediction tools. StarBase was used to predict candidate miRNAs that may bind to BMP1 and candidate lncRNAs that may bind to hsa-miR-532-3p. We analyzed the association between expression of BMP1 and immune cell infiltration levels in ccRCC using the TIMER website. The relationship between BMP1 expression levels and immune checkpoint expression levels was also investigated. RESULTS: BMP1 was upregulated in GBM, HNSC, KIRC, KIRP and STAD and downregulated in KICH and PRAD. Combined with OS and DFS, BMP1 can be used as a biomarker for poor prognosis among patients with KIRC. Through expression analysis, survival analysis and correlation analysis, LINC00685, SLC16A1-AS1, PVT1, VPS9D1-AS1, SNHG15 and the CCDC18-AS1/hsa-miR-532-3p/BMP1 axis were established as the most potential upstream ncRNA-related pathways of BMP1 in ccRCC. Furthermore, we found that BMP1 levels correlated significantly positively with tumor immune cell infiltration, biomarkers of immune cells, and immune checkpoint expression. CONCLUSION: Our results demonstrate that ncRNA-mediated high expression of BMP1 is associated with poor prognosis and tumor immune infiltration in ccRCC.

Multivalent human antibody-centyrin fusion protein to prevent and treat Staphylococcus aureus infections.

Treating and preventing infections by antimicrobial-resistant bacterial pathogens is a worldwide problem. Pathogens such as Staphylococcus aureus produce an array of virulence determinants, making it difficult to identify single targets for the development of vaccines or monoclonal therapies. We described a human-derived anti-S. aureus monoclonal antibody (mAb)-centyrin fusion protein ("mAbtyrin") that simultaneously targets multiple bacterial adhesins, resists proteolysis by bacterial protease GluV8, avoids Fc engagement by S. aureus IgG-binding proteins SpA and Sbi, and neutralizes pore-forming leukocidins via fusion with anti-toxin centyrins, while maintaining Fc- and complement-mediated functions. Compared with the parental mAb, mAbtyrin protected human phagocytes and boosted phagocyte-mediated killing. The mAbtyrin also reduced pathology, reduced bacterial burden, and protected from different types of infections in preclinical animal models. Finally, mAbtyrin synergized with vancomycin, enhancing pathogen clearance in an animal model of bacteremia. Altogether, these data establish the potential of multivalent mAbs for treating and preventing S. aureus diseases.

"Stapling" scFv for multispecific biotherapeutics of superior properties.

Single-chain fragment variable (scFv) domains play an important role in antibody-based therapeutic modalities, such as bispecifics, multispecifics and chimeric antigen receptor T cells or natural killer cells. However, scFv domains exhibit lower stability and increased risk of aggregation due to transient dissociation ("breathing") and inter-molecular reassociation of the two domains (VL and VH). We designed a novel strategy, referred to as stapling, that introduces two disulfide bonds between the scFv linker and the two variable domains to minimize scFv breathing. We named the resulting molecules stapled scFv (spFv). Stapling increased thermal stability (Tm) by an average of 10°C. In multiple scFv/spFv multispecifics, the spFv molecules display significantly improved stability, minimal aggregation and superior product quality. These spFv multispecifics retain binding affinity and functionality. Our stapling design was compatible with all antibody variable regions we evaluated and may be widely applicable to stabilize scFv molecules for designing biotherapeutics with superior biophysical properties.

Optimizing nucleotide sequence ensembles for combinatorial protein libraries using a genetic algorithm.

Protein libraries are essential to the field of protein engineering. Increasingly, probabilistic protein design is being used to synthesize combinatorial protein libraries, which allow the protein engineer to explore a vast space of amino acid sequences, while at the same time placing restrictions on the amino acid distributions. To this end, if site-specific amino acid probabilities are input as the target, then the codon nucleotide distributions that match this target distribution can be used to generate a partially randomized gene library. However, it turns out to be a highly nontrivial computational task to find the codon nucleotide distributions that exactly matches a given target distribution of amino acids. We first showed that for any given target distribution an exact solution may not exist at all. Formulated as a constrained optimization problem, we then developed a genetic algorithm-based approach to find codon nucleotide distributions that match as closely as possible to the target amino acid distribution. As compared with the previous gradient descent method on various objective functions, the new method consistently gave more optimized distributions as measured by the relative entropy between the calculated and the target distributions. To simulate the actual lab solutions, new objective functions were designed to allow for two separate sets of codons in seeking a better match to the target amino acid distribution.

Endoplasmic Reticulum Stress Mediates Renal Tubular Vacuolation in BK Polyomavirus-Associated Nephropathy.

Objective: This study aimed to explore the molecular mechanism of cytoplasmic vacuolation caused by BK polyomavirus (BKPyV) and thus search for potential target for drug repurposing. Methods: Morphological features of BK polyomavirus-associated nephropathy (BKPyVAN) were studied under light and electron microscopes. Microarray datasets GSE75693, GSE47199, and GSE72925 were integrated by ComBat, and differentially expressed genes (DEGs) were analyzed using limma. Furthermore, the endoplasmic reticulum (ER)-related genes obtained from GenCLiP 2.0 were intersected with DEGs. GO and KEGG enrichment pathways were performed with intersection genes by R package clusterProfiler. The single-cell RNA sequencing (scRNA-seq) from a BKPyVAN recipient was analyzed with a dataset (GSE140989) downloaded from Gene Expression Omnibus (GEO) as control for gene set variation analysis (GSVA). Immunohistochemistry and electron microscopy of kidney sections from drug-induced ERS mouse models were performed to explore the association of ERS and renal tubular vacuolation. Protein-protein interaction (PPI) network of the intersection genes was constructed to identify hub target. AutoDock was used to screen Food and Drug Administration (FDA)-approved drugs that potentially targeted hub gene. Results: Light and electron microscopes exhibited obvious intranuclear inclusions, vacuoles, and virus particles in BKPyV-infected renal tubular cells. Transcriptome analysis revealed 629 DEGs between samples of BKPyVAN and stable transplanted kidneys, of which 16 were ER-associated genes. GO analysis with the intersection genes illustrated that ERS-related pathways were significantly involved, and KEGG analysis showed a prominent enrichment of MAPK, Toll-like receptor, and chemokine signaling pathways. GSVA analysis of the proximal tubule revealed similar pathways enrichment. An electron microscope image of the kidney from ERS mouse models showed an obvious renal tubular vacuolation with prominent activation of ERS markers verified by immunohistochemistry. Furthermore, DDIT3 was identified as the hub gene based on PPI analysis, and ZINCOOOOO1531009 (Risedronate) was indicated to be a potential drug for DDIT3. Conclusion: ERS was involved in renal tubular cytoplasmic vacuolation in BKPyVAN recipients. Risedronate was screened as a potential drug for BKPyVAN by targeting DDIT3.

Discovery and pharmacological characterization of cetrelimab (JNJ-63723283), an anti-programmed cell death protein-1 (PD-1) antibody, in human cancer models.

PURPOSE: Preclinical characterization of cetrelimab (JNJ-63723283), a fully humanized immunoglobulin G4 kappa monoclonal antibody targeting programmed cell death protein-1 (PD-1), in human cancer models. METHODS: Cetrelimab was generated by phage panning against human and cynomolgus monkey (cyno) PD-1 extracellular domains (ECDs) and affinity maturation. Binding to primate and rodent PD-1 ECDs, transfected and endogenous cell-surface PD-1, and inhibition of ligand binding were measured. In vitro activity was evaluated using cytomegalovirus recall, mixed lymphocyte reaction, staphylococcal enterotoxin B stimulation, and Jurkat-PD-1 nuclear factor of activated T cell reporter assays. In vivo activity was assessed using human PD-1 knock-in mice implanted with MC38 tumors and a lung patient-derived xenograft (PDX) model (LG1306) using CD34 cord-blood-humanized NSG mice. Pharmacodynamics, toxicokinetics, and safety were assessed in cynos following single and/or repeat intravenous dosing. RESULTS: Cetrelimab showed high affinity binding to human (1.72 nM) and cyno (0.90 nM) PD-1 and blocked binding of programmed death-ligand 1 (PD-L1; inhibitory concentration [IC] 111.7 ng/mL) and PD-L2 (IC 138.6 ng/mL). Cetrelimab dose-dependently increased T cell-mediated cytokine production and stimulated cytokine expression. Cetrelimab 10 mg/kg reduced mean MC38 tumor volume in PD-1 knock-in mice at Day 21 (P < 0.0001) versus control. In a PDX lung model, 10 mg/kg cetrelimab (every 5 days for six cycles) increased frequency of peripheral T cells and reduced (P < 0.05) mean tumor volume versus control. Activity was consistent with that of established PD-1 inhibitors. Cetrelimab dosing was well tolerated in cynos and mean drug exposure increase was dose-dependent. CONCLUSION: Cetrelimab potently inhibits PD-1 in vitro and in vivo, supporting its clinical evaluation.

Antibody modeling assessment.

A blinded study to assess the state of the art in three-dimensional structure modeling of the variable region (Fv) of antibodies was conducted. Nine unpublished high-resolution x-ray Fab crystal structures covering a wide range of antigen-binding site conformations were used as benchmark to compare Fv models generated by four structure prediction methodologies. The methodologies included two homology modeling strategies independently developed by CCG (Chemical Computer Group) and Accerlys Inc, and two fully automated antibody modeling servers: PIGS (Prediction of ImmunoGlobulin Structure), based on the canonical structure model, and Rosetta Antibody Modeling, based on homology modeling and Rosetta structure prediction methodology. The benchmark structure sequences were submitted to Accelrys and CCG and a set of models for each of the nine antibody structures were generated. PIGS and Rosetta models were obtained using the default parameters of the servers. In most cases, we found good agreement between the models and x-ray structures. The average rmsd (root mean square deviation) values calculated over the backbone atoms between the models and structures were fairly consistent, around 1.2 Å. Average rmsd values of the framework and hypervariable loops with canonical structures (L1, L2, L3, H1, and H2) were close to 1.0 Å. H3 prediction yielded rmsd values around 3.0 Å for most of the models. Quality assessment of the models and the relative strengths and weaknesses of the methods are discussed. We hope this initiative will serve as a model of scientific partnership and look forward to future antibody modeling assessments.