Abstract
High-resolution homology models are useful in structure-based protein engineering applications, especially when a crystallographic structure is unavailable. Here, we report the development and implementation of RosettaAntibody, a protocol for homology modeling of antibody variable regions. The protocol combines comparative modeling of canonical complementarity determining region (CDR) loop conformations and de novo loop modeling of CDR H3 conformation with simultaneous optimization of V(L)-V(H) rigid-body orientation and CDR backbone and side-chain conformations. The protocol was tested on a benchmark of 54 antibody crystal structures. The median root mean square deviation (rmsd) of the antigen binding pocket comprised of all the CDR residues was 1.5 A with 80% of the targets having an rmsd lower than 2.0 A. The median backbone heavy atom global rmsd of the CDR H3 loop prediction was 1.6, 1.9, 2.4, 3.1, and 6.0 A for very short (4-6 residues), short (7-9), medium (10-11), long (12-14) and very long (17-22) loops, respectively. When the set of ten top-scoring antibody homology models are used in local ensemble docking to antigen, a moderate-to-high accuracy docking prediction was achieved in seven of fifteen targets. This success in computational docking with high-resolution homology models is encouraging, but challenges still remain in modeling antibody structures for sequences with long H3 loops. This first large-scale antibody-antigen docking study using homology models reveals the level of "functional accuracy" of these structural models toward protein engineering applications.