A Functional Investigation of Antibody Fc-FcRn Variant Binding Guided by In Silico Free Energy Perturbation Methods

Accurate calculation of energy changes upon mutation is a key requirement for the effective use of computational methods in protein design. In this study, we applied free energy perturbation (FEP) calculations to predict the effects of mutations on the binding free energy between the immunoglobulin G (IgG) antibody fragment-crystallizable (Fc) region and the neonatal Fc receptor (FcRn), an interaction that is primarily responsible for antibody half-life. We assembled an extensive experimental dataset of Fc-FcRn binding affinities for wild-type (wt) and mutant complexes, including values from literature and from newly measured results. Starting from a crystal structure of the M252Y/S254T/T256E ("YTE") Fc variant bound to FcRn, we prepared all-atom models of human IgG1-subtype wt and YTE variant Fc-FcRn complexes, adding explicit hydrogens and assigning protonation states for key ionizable residues. Initial results using standard FEP protocols to compute relative binding free energies were promising but exhibited multiple outliers. By accounting for coupling effects for FEP mutations near key histidine residues, we improved the results for several outliers, suggesting such coupling as an important approach for pH-sensitive systems. Further, upon determining new crystal structures of wt Fc and three Fc variants at multiple pH values, we observed subtle conformational changes in unbound Fc; by accounting for these conformational changes in FEP calculations, we additionally improved agreement with experiment. The detailed structural and energetic analyses of the Fc-FcRn system we present here thus provide a practical energy-calculation framework to enable rational in silico design of novel Fc variants. SIGNIFICANCE: Antibody mutations that increase half-life are of high medical importance; unfortunately, these have been difficult to predict computationally. In this study of the Fc-FcRn complex, which controls antibody half-life via differential affinity at different pHs, we demonstrate a successful computational approach using free energy perturbation (FEP) calculations. Our prediction of accurate binding energies across a wide range of cases speaks to the power of the FEP methodology in navigating the free energy landscapes of dynamic molecular complexes. Furthermore, we show that accurate Fc-FcRn affinity calculations required careful consideration of conformational flexibility between bound and unbound states, contributing to our functional understanding of a system that will be important for future rational antibody-design efforts.