Highland games: A benchmarking exercise in predicting biophysical and drug properties of monoclonal antibodies from amino acid sequences.

Biopharmaceutical product and process development do not yet take advantage of predictive computational modeling to nearly the degree seen in industries based on smaller molecules. To assess and advance progress in this area, spirited coopetition (mutually beneficial collaboration between competitors) was successfully used to motivate industrial scientists to develop, share, and compare data and methods which would normally have remained confidential. The first "Highland Games" competition was held in conjunction with the October 2018 Recovery of Biological Products Conference in Ashville, NC, with the goal of benchmarking and assessment of the ability to predict development-related properties of six antibodies from their amino acid sequences alone. Predictions included purification-influencing properties such as isoelectric point and protein A elution pH, and biophysical properties such as stability and viscosity at very high concentrations. Essential contributions were made by a large variety of individuals, including companies which consented to provide antibody amino acid sequences and test materials, volunteers who undertook the preparation and experimental characterization of these materials, and prediction teams who attempted to predict antibody properties from sequence alone. Best practices were identified and shared, and areas in which the community excels at making predictions were identified, as well as areas presenting opportunities for considerable improvement. Predictions of isoelectric point and protein A elution pH were especially good with all-prediction average errors of 0.2 and 1.6 pH unit, respectively, while predictions of some other properties were notably less good. This manuscript presents the events, methods, and results of the competition, and can serve as a tutorial and as a reference for in-house benchmarking by others. Organizations vary in their policies concerning disclosure of methods, but most managements were very cooperative with the Highland Games exercise, and considerable insight into common and best practices is available from the contributed methods. The accumulated data set will serve as a benchmarking tool for further development of in silico prediction tools.

Single Molecule Force Spectroscopy and Molecular Dynamics Simulations as a Combined Platform for Probing Protein Face-Specific Binding.

Biomolecular interactions frequently occur in orientation-specific manner. For example, prior nuclear magnetic resonance spectroscopy experiments in our lab have suggested the presence of a group of strongly binding residues on a particular face of the protein ubiquitin for interactions with Capto MMC multimodal ligands ("Capto" ligands) (Srinivasan, K.; et al. Langmuir 2014, 30 (44), 13205-13216). We present a clear confirmation of those studies by performing single-molecule force spectroscopy (SMFS) measurements of unbinding complemented with molecular dynamics (MD) calculations of the adsorption free energy of ubiquitin in two distinct orientations with self-assembled monolayers (SAMs) functionalized with "Capto" ligands. These orientations were maintained in the SMFS experiments by tethering ubiquitin mutants to SAM surfaces through strategically located cysteines, thus exposing the desired faces of the protein. Analogous orientations were maintained in MD simulations using suitable constraining methods. Remarkably, despite differences between the finer details of experimental and simulation methodologies, they confirm a clear preference for the previously hypothesized binding face of ubiquitin. Furthermore, MD simulations provided significant insights into the mechanism of protein binding onto this multimodal surface. Because SMFS and MD simulations both directly probe protein-surface interactions, this work establishes a key link between experiments and simulations at molecular scale through the determination of protein face-specific binding energetics. Our approach may have direct applications in biophysical systems where face- or orientation-specific interactions are important, such as biomaterials, sensors, and biomanufacturing.

A molecular modeling based method to predict elution behavior and binding patches of proteins in multimodal chromatography.

Multimodal (MM) chromatography provides a powerful means to enhance the selectivity of protein separations by taking advantage of multiple weak interactions that include electrostatic, hydrophobic and van der Waals interactions. In order to increase our understanding of such phenomena, a computationally efficient approach was developed that combines short molecular dynamics simulations and continuum solvent based coarse-grained free energy calculations in order to study the binding of proteins to Self Assembled Monolayers (SAM) presenting MM ligands. Using this method, the free energies of protein-MM SAM binding over a range of incident orientations of the protein can be determined. The resulting free energies were then examined to identify the more "strongly bound" orientations of different proteins with two multimodal surfaces. The overall free energy of protein-MM surface binding was then determined and correlated to retention factors from isocratic chromatography. This correlation, combined with analytical expressions from the literature, was then employed to predict protein gradient elution salt concentrations as well as selectivity reversals with different MM resin systems. Patches on protein surfaces that interacted strongly with MM surfaces were also identified by determining the frequency of heavy atom contacts with the atoms of the MM SAMs. A comparison of these patches to Electrostatic Potential and hydrophobicity maps indicated that while all of these patches contained significant positive charge, only the highest frequency sites also possessed hydrophobicity. The ability to identify key binding patches on proteins may have significant impact on process development for the separation of bioproduct related impurities.

Effect of guanidine and arginine on protein-ligand interactions in multimodal cation-exchange chromatography.

The addition of fluid phase modifiers provides significant opportunities for increasing the selectivity of multimodal chromatography. In order to optimize this selectivity, it is important to understand the fundamental interactions between proteins and these modifiers. To this end, molecular dynamics (MD) simulations were first performed to study the interactions of guanidine and arginine with three proteins. The simulation results showed that both guanidine and arginine interacted primarily with the negatively charged regions on the proteins and that these regions could be readily predicted using electrostatic potential maps. Protein surface characterization was then carried out using computationally efficient coarse-grained techniques for a broader set of proteins which exhibited interesting chromatographic retention behavior upon the addition of these modifiers. It was shown that proteins exhibiting an increased retention in the presence of guanidine possessed hydrophobic regions adjacent to negatively charged regions on their surfaces. In contrast, proteins which exhibited a decreased binding in the presence of guanidine did not have hydrophobic regions adjacent to negatively charged patches. These results indicated that the effect of guanidine could be described as a combination of competitive binding, charge neutralization and increased hydrophobic interactions for certain proteins. In contrast, arginine resulted in a significant decrease in protein retention times primarily due to competition for the resin and steric effects, with minimal accompanying increase in hydrophobic interactions. The approach presented in this paper which employs MD simulations to guide the application of coarse-grained approaches is expected to be extremely useful for methods development in downstream bioprocesses. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 33:435-447, 2017.

The effect of geometrical presentation of multimodal cation-exchange ligands on selective recognition of hydrophobic regions on protein surfaces.

The effects of spatial organization of hydrophobic and charged moieties on multimodal (MM) cation-exchange ligands were examined by studying protein retention behavior on two commercial chromatographic media, Capto™ MMC and Nuvia™ cPrime™. Proteins with extended regions of surface-exposed aliphatic residues were found to have enhanced retention on the Capto MMC system as compared to the Nuvia cPrime resin. The results further indicated that while the Nuvia cPrime ligand had a strong preference for interactions with aromatic groups, the Capto MMC ligand appeared to interact with both aliphatic and aromatic clusters on the protein surfaces. These observations were formalized into a new set of protein surface property descriptors, which quantified the local distribution of electrostatic and hydrophobic potentials as well as distinguishing between aromatic and aliphatic properties. Using these descriptors, high-performing quantitative structure-activity relationship (QSAR) models (R(2)>0.88) were generated for both the Capto MMC and Nuvia cPrime datasets at pH 5 and pH 6. Descriptors of electrostatic properties were generally common across the four models; however both Capto MMC models included descriptors that quantified regions of aliphatic-based hydrophobicity in addition to aromatic descriptors. Retention was generally reduced by lowering the ligand densities on both MM resins. Notably, elution order was largely unaffected by the change in surface density, but smaller and more aliphatic proteins tended to be more affected by this drop in ligand density. This suggests that modulating the exposure, shape and density of the hydrophobic moieties in multimodal chromatographic systems can alter the preference for surface exposed aliphatic or aromatic residues, thus providing an additional dimension for modulating the selectivity of MM protein separation systems.

Interactions of Multimodal Ligands with Proteins: Insights into Selectivity Using Molecular Dynamics Simulations.

Fundamental understanding of protein-ligand interactions is important to the development of efficient bioseparations in multimodal chromatography. Here we employ molecular dynamics (MD) simulations to investigate the interactions of three different proteins--ubiquitin, cytochrome C, and α-chymotrypsinogen A, sampling a range of charge from +1e to +9e--with two multimodal chromatographic ligands containing similar chemical moieties--aromatic, carboxyl, and amide--in different structural arrangements. We use a spherical harmonic expansion to analyze ligand and individual moiety density profiles around the proteins. We find that the Capto MMC ligand, which contains an additional aliphatic group, displays stronger interactions than Nuvia CPrime ligand with all three proteins. Studying the ligand densities at the moiety level suggests that hydrophobic interactions play a major role in determining the locations of high ligand densities. Finally, the greater structural flexibility of the Capto MMC ligand compared to that of the Nuvia cPrime ligand allows for stronger structural complementarity and enables stronger hydrophobic interactions. These subtle and not-so-subtle differences in binding affinities and modalities for multimodal ligands can result in significantly different binding behavior towards proteins with important implications for bioprocessing.

Investigation of protein selectivity in multimodal chromatography using in silico designed Fab fragment variants.

In this study, a unique set of antibody Fab fragments was designed in silico and produced to examine the relationship between protein surface properties and selectivity in multimodal chromatographic systems. We hypothesized that multimodal ligands containing both hydrophobic and charged moieties would interact strongly with protein surface regions where charged groups and hydrophobic patches were in close spatial proximity. Protein surface property characterization tools were employed to identify the potential multimodal ligand binding regions on the Fab fragment of a humanized antibody and to evaluate the impact of mutations on surface charge and hydrophobicity. Twenty Fab variants were generated by site-directed mutagenesis, recombinant expression, and affinity purification. Column gradient experiments were carried out with the Fab variants in multimodal, cation-exchange, and hydrophobic interaction chromatographic systems. The results clearly indicated that selectivity in the multimodal system was different from the other chromatographic modes examined. Column retention data for the reduced charge Fab variants identified a binding site comprising light chain CDR1 as the main electrostatic interaction site for the multimodal and cation-exchange ligands. Furthermore, the multimodal ligand binding was enhanced by additional hydrophobic contributions as evident from the results obtained with hydrophobic Fab variants. The use of in silico protein surface property analyses combined with molecular biology techniques, protein expression, and chromatographic evaluations represents a previously undescribed and powerful approach for investigating multimodal selectivity with complex biomolecules.

Investigation into the molecular and thermodynamic basis of protein interactions in multimodal chromatography using functionalized nanoparticles.

Although multimodal chromatography offers significant potential for bioseparations, there is a lack of molecular level understanding of the nature of protein binding in these systems. In this study a nanoparticle system is employed that can simulate a chromatographic resin surface while also being amenable to isothermal titration calorimetry (ITC) and solution NMR. ITC and NMR titration experiments are carried out with (15)N-labeled ubiquitin to investigate the interactions of ubiquitin with nanoparticles functionalized with two industrially important multimodal ligands. The ITC results suggest that binding to both multimodal ligand surfaces is entropically driven over a range of temperatures and that this is due primarily to the release of surface bound waters. In order to reveal structural details of the interaction process, binding-induced chemical shift changes obtained from the NMR experiments are employed to obtain dissociation constants of individual amino acid residues on the protein surface. The residue level information obtained from NMR is then used to identify a preferred binding face on ubiquitin for interaction to both multimodal ligand surfaces. In addition, electrostatic potential and spatial aggregation propensity maps are used to determine important protein surface property data that are shown to correlate well with the molecular level information obtained from NMR. Importantly, the data demonstrate that the cluster of interacting residues on the protein surface act co-operatively to give rise to multimodal binding affinities several orders of magnitude greater than those obtained previously for interactions with free solution ligands. The use of NMR and ITC to study protein interactions with functionalized nanoparticles offers a new tool for obtaining important molecular and thermodynamic insights into protein affinity in multimodal chromatographic systems.

Application of a spherical harmonics expansion approach for calculating ligand density distributions around proteins.

Protein-ligand interactions are central to many biological applications, including molecular recognition, protein formulations, and bioseparations. Complex, multisite ligands can have affinities for different locations on a protein's surface, depending on the chemical and topographical complementarity. We employ an approach based on the spherical harmonic expansion to calculate spatially resolved three-dimensional atomic density profiles of water and ligands in the vicinity of macromolecules. To illustrate the approach, we first study the hydration of model C180 buckyball solutes, with nonspherical patterns of hydrophobicity/-philicity on their surface. We extend the approach to calculate density profiles of increasingly complex ligands and their constituent groups around a protein (ubiquitin) in aqueous solution. Analysis of density profiles provides information about the binding face of the protein and the preferred orientations of ligands on the binding surface. Our results highlight that the spherical harmonic expansion based approach is easy to implement and efficient for calculation and visualization of three-dimensional density profiles around spherically nonsymmetric and topographically and chemically complex solutes.

Effects of urea on selectivity and protein-ligand interactions in multimodal cation exchange chromatography.

Nuclear magnetic resonance (NMR) and molecular dynamics (MD) simulations were employed in concert with chromatography to provide insight into the effect of urea on protein-ligand interactions in multimodal (MM) chromatography. Chromatographic experiments with a protein library in ion exchange (IEX) and MM systems indicated that, while urea had a significant effect on protein retention and selectivity for a range of proteins in MM systems, the effects were much less pronounced in IEX. NMR titration experiments carried out with a multimodal ligand, and isotopically enriched human ubiquitin indicated that, while the ligand binding face of ubiquitin remained largely intact in the presence of urea, the strength of binding was decreased. MD simulations were carried out to provide further insight into the effect of urea on MM ligand binding. These results indicated that, while the overall ligand binding face of ubiquitin remained the same, there was a reduction in the occupancy of the MM ligand interaction region along with subtle changes in the residues involved in these interactions. This work demonstrates the effectiveness of urea in enhancing selectivity in MM chromatographic systems and also provides an in-depth analysis of how MM ligand-protein interactions are altered in the presence of this fluid phase modifier.

Probing multimodal ligand binding regions on ubiquitin using nuclear magnetic resonance, chromatography, and molecular dynamics simulations.

Site-directed mutagenesis, nuclear magnetic resonance (NMR) chemical shift perturbation experiments, and molecular dynamics (MD) simulations are employed in concert with chromatographic experiments to provide insight into protein-ligand interactions in multimodal chromatographic systems. In previous studies, a preferred binding region was identified on the surface of the protein ubiquitin for binding with a multimodal ligand. In this study, site-directed mutagenesis is used to enable direct NMR evaluation of the mutant protein as compared to the wild type. It is found that reversing the charge of a key residue (K6E) in the proposed preferred binding region results in substantial decreases in the magnitude of the ligand-induced NMR chemical shift perturbations relative to those detected for the wild type protein, particularly for residues located within the preferred binding region. These NMR results also indicate a decrease in ligand affinity, consistent with the weaker chromatographic retention observed for the mutant as compared to the wild type on a multimodal cation exchange resin. MD simulation results provide additional insight at a molecular level and demonstrate that many residues located within the preferred binding region exhibit weaker binding interactions due to the mutation. The analysis suggests that multimodal ligand binding consists of initial localization of the ligand by long-ranged electrostatic interactions followed by multiple short-ranged synergistic interactions to attain high affinities of the ligand to specific residues.

Mobile phase modifier effects in multimodal cation exchange chromatography.

This study examines protein adsorption behavior and the effects of mobile phase modifiers in multimodal chromatographic systems. Chromatography results with a diverse protein library indicate that multimodal and ion exchange resins have markedly different protein binding behavior and selectivity. NMR results corroborate the stronger binding observed for the multimodal system and provide insight into the structural basis for the observed binding behavior. Protein-binding affinity and selectivity in multimodal and ion exchange systems are then examined using a variety of mobile phase modifiers. Arginine and guanidine are found to have dramatic effects on protein adsorption, yielding changes in selectivity in both chromatographic systems. While sodium caprylate leads to slightly weaker chromatographic retention for most proteins, certain proteins exhibit significant losses in retention in both systems. The presence of a competitive binding mechanism between the multimodal ligand and sodium caprylate for binding to ubiquitin is confirmed using STD NMR. Polyol mobile phase modifiers are shown to result in increased retention for weakly bound proteins and decreased retention for strongly bound proteins, indicating that the overall retention behavior is determined by a balance between changes in electrostatic and hydrophobic interactions. This work provides an improved understanding of protein adsorption and mobile phase modifier effects in multimodal chromatographic systems and sets the stage for future work to develop more selective protein separation systems.