Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies

Adsorption of monoclonal antibody fragments at the water-oil interface: A coarse-grained molecular dynamics study

Monoclonal antibodies (mAbs) can undergo structural changes due to interaction with oil-water interfaces during storage. Such changes can lead to aggregation, resulting in a loss of therapeutic efficacy. Therefore, understanding the microscopic mechanism controlling mAb adsorption is crucial to developing strategies that can minimize the impact of interfaces on the therapeutic properties of mAbs. In this study, we used MARTINI coarse-grained molecular dynamics simulations to investigate the adsorption of the Fab and Fc domains of the monoclonal antibody COE3 at the oil-water interface. Our aim was to determine the regions on the protein surface that drive mAb adsorption. We also investigate the role of protein concentration on protein orientation and protrusion to the oil phase. While our structural analyses compare favorably with recent neutron reflectivity measurements, we observe some differences. Unlike the monolayer at the interface predicted by neutron reflectivity experiments, our simulations indicate the presence of a secondary diffused layer near the interface. We also find that under certain conditions, protein-oil interaction can lead to a considerable distortion in the protein structure, resulting in enhanced adsorption behavior.

Exploring the Combined Action of Adding Pertuzumab to Branded Trastuzumab versus Trastuzumab Biosimilars for Treating HER2+ Breast Cancer

The binding activity of various trastuzumab biosimilars versus the branded trastuzumab towards the glycosylated extracellular domain of the human epidermal growth factor receptor 2 (HER2) target in the presence of pertuzumab was investigated. We employed size exclusion chromatography with tetra-detection methodology to simultaneously determine absolute molecular weight, concentration, molecular size, and intrinsic viscosity. All trastuzumab molecules in solution exhibit analogous behavior in their binary action towards HER2 regardless of the order of addition of trastuzumab/pertuzumab. This analogous behavior of all trastuzumab molecules, including biosimilars, highlights the robustness and consistency of their binding activity towards HER2. Furthermore, the addition of HER2 to a mixture of trastuzumab and pertuzumab leads to increased formation of high-order HER2 complexes, up to concentrations of one order of magnitude higher than in the case of sequential addition. The observed increase suggests a potential synergistic effect between these antibodies, which could enhance their therapeutic efficacy in HER2-positive cancers. These findings underscore the importance of understanding the complex interplay between therapeutic antibodies and their target antigens, providing valuable insights for the development of more effective treatment strategies.

Unraveling the Microscopic Mechanism of Molecular Ion Interaction with Monoclonal Antibodies: Impact on Protein Aggregation

Understanding and predicting protein aggregation represents one of the major challenges in accelerating the pharmaceutical development of protein therapeutics. In addition to maintaining the solution pH, buffers influence both monoclonal antibody (mAb) aggregation in solution and the aggregation mechanisms since the latter depend on the protein charge. Molecular-level insight is necessary to understand the relationship between the buffer-mAb interaction and mAb aggregation. Here, we use all-atom molecular dynamics simulations to investigate the interaction of phosphate (Phos) and citrate (Cit) buffer ions with the Fab and Fc domains of mAb COE3. We demonstrate that Phos and Cit ions feature binding mechanisms, with the protein that are very different from those reported previously for histidine (His). These differences are reflected in distinctive ion-protein binding modes and adsorption/desorption kinetics of the buffer molecules from the mAb surface and result in dissimilar effects of these buffer species on mAb aggregation. While His shows significant affinity toward hydrophobic amino acids on the protein surface, Phos and Cit ions preferentially bind to charged amino acids. We also show that Phos and Cit anions provide bridging contacts between basic amino acids in neighboring proteins. The implications of such contacts and their connection to mAb aggregation in therapeutic formulations are discussed.

Mechanistic Insights into the Adsorption of Monoclonal Antibodies at the Water/Vapor Interface

Monoclonal antibodies (mAbs) are active components of therapeutic formulations that interact with the water-vapor interface during manufacturing, storage, and administration. Surface adsorption has been demonstrated to mediate antibody aggregation, which leads to a loss of therapeutic efficacy. Controlling mAb adsorption at interfaces requires a deep understanding of the microscopic processes that lead to adsorption and identification of the protein regions that drive mAb surface activity. Here, we report all-atom molecular dynamics (MD) simulations of the adsorption behavior of a full IgG1-type antibody at the water/vapor interface. We demonstrate that small local changes in the protein structure play a crucial role in promoting adsorption. Also, interfacial adsorption triggers structural changes in the antibody, potentially contributing to the further enhancement of surface activity. Moreover, we identify key amino acid sequences that determine the adsorption of antibodies at the water-air interface and outline strategies to control the surface activity of these important therapeutic proteins.

Structure and interaction of therapeutic proteins in solution: a combined simulation and experimental study

The aggregation of therapeutic proteins in solution has attracted significant interest, driving efforts to understand the relationship between microscopic structural changes and protein-protein interactions determining aggregation processes in solution. Additionally, there is substantial interest in being able to predict aggregation based on protein structure as part of molecular developability assessments. Molecular Dynamics provides theoretical tools to complement experimental studies and to interrogate and identify the microscopic mechanisms determining aggregation. Here we perform all-atom MD simulations to study the structure and inter-protein interaction of the Fab and Fc fragments of the monoclonal antibody (mAb) COE3. We unravel the role of ion-protein interactions in building the ionic double layer and determining effective inter-protein interaction. Further, we demonstrate, using various state-of-the-art force fields (charmm, gromos, amber, opls/aa), that the protein solvation, ionic structure and protein-protein interaction depend significantly on the force field parameters. We perform SANS and Static Light Scattering experiments to assess the accuracy of the different forcefields. Comparison of the simulated and experimental results reveal significant differences in the forcefields' performance, particularly in their ability to predict the protein size in solution and inter-protein interactions quantified through the second virial coefficients. In addition, the performance of the forcefields is correlated with the protein hydration structure.

Investigating the Orientation of an Interfacially Adsorbed Monoclonal Antibody and Its Fragments Using Neutron Reflection

Interfacial adsorption is a molecular process occurring during the production, purification, transport, and storage of antibodies, with a direct impact on their structural stability and subsequent implications on their bioactivities. While the average conformational orientation of an adsorbed protein can be readily determined, its associated structures are more complex to characterize. Neutron reflection has been used in this work to investigate the conformational orientations of the monoclonal antibody COE-3 and its Fab and Fc fragments at the oil/water and air/water interfaces. Rigid body rotation modeling was found to be suitable for globular and relatively rigid proteins such as the Fab and Fc fragments but less so for relatively flexible proteins such as full COE-3. Fab and Fc fragments adopted a 'flat-on' orientation at the air/water interface, minimizing the thickness of the protein layer, but they adopted a substantially tilted orientation at the oil/water interface with increased layer thickness. In contrast, COE-3 was found to adsorb in tilted orientations at both interfaces, with one fragment protruding into the solution. This work demonstrates that rigid-body modeling can provide additional insights into protein layers at various interfaces relevant to bioprocess engineering.

Stability Convergence in Antibody Coformulations

Combined administration of antibody therapeutics has proven to be beneficial for patients with cancer or infectious diseases. As a result, there is a growing trend toward multiple antibodies premixed into a single product form and delivered to patients as a fixed-dose coformulation. However, combining antibodies into a single coformulation could be challenging as proteins have the potential to interact and alter their stability and degradation profiles in the mixture, compared to that in isolation. We show that in two specific antibody-antibody coformulations, the more stable antibody component increased the stability of the less stable component, which in return destabilized the more stable component, hence exhibiting an overall convergence of stability in the coformulation.

Understanding the Stabilizing Effect of Histidine on mAb Aggregation: A Molecular Dynamics Study

Histidine, a widely used buffer in monoclonal antibody (mAb) formulations, is known to reduce antibody aggregation. While experimental studies suggest a nonelectrostatic, nonstructural (relating to secondary structure preservation) origin of the phenomenon, the underlying microscopic mechanism behind the histidine action is still unknown. Understanding this mechanism will help evaluate and predict the stabilizing effect of this buffer under different experimental conditions and for different mAbs. We have used all-atom molecular dynamics simulations and contact-based free energy calculations to investigate molecular-level interactions between the histidine buffer and mAbs, which lead to the observed stability of therapeutic formulations in the presence of histidine. We reformulate the Spatial Aggregation Propensity index by including the buffer–protein interactions. The buffer adsorption on the protein surface leads to lower exposure of the hydrophobic regions to water. Our analysis indicates that the mechanism behind the stabilizing action of histidine is connected to the shielding of the solvent-exposed hydrophobic regions on the protein surface by the buffer molecules.

Open Circuit Potential as a Tool for the Assessment of Binding Kinetics and Reagentless Protein Quantitation

A microfluidic open circuit potential label-free protein assay was developed for the reagentless quantification of C-reactive protein (CRP), a model protein target, and further utilized to assess target-receptor binding kinetics. Generated sensors have very high baseline stabilities (<1% change in 100 min) and high levels of selectivity in complex media. Real-time assays are fast (<20 min), of high sensitivity (1 ng/mL limit of detection for CRP in serum), and resolve kinetic and thermodynamic characteristics that correlate well with those resolved optically. The assay shows excellent correlation with an enzyme-linked immunosorbent assay analysis of patient samples. The methodology has value in potentially underpinning a low-cost, rapid, and sensitive single-step biomarker quantification.

Towards an improved prediction of concentrated antibody solution viscosity using the Huggins coefficient

The viscosity of a monoclonal antibody solution must be monitored and controlled as it can adversely affect product processing, packaging and administration. Engineering low viscosity mAb formulations is challenging as prohibitive amounts of material are required for concentrated solution analysis, and it is difficult to predict viscosity from parameters obtained through low-volume, high-throughput measurements such as the interaction parameter, k_D, and the second osmotic virial coefficient, B₂₂. As a measure encompassing the effect of intermolecular interactions on dilute solution viscosity, the Huggins coefficient, k_h, is a promising candidate as a parameter measureable at low concentrations, but indicative of concentrated solution viscosity. In this study, a differential viscometry technique is developed to measure the intrinsic viscosity, [η], and the Huggins coefficient, k_h, of protein solutions. To understand the effect of colloidal protein-protein interactions on the viscosity of concentrated protein formulations, the viscometric parameters are compared to k_D and B₂₂ of two mAbs, tuning the contributions of repulsive and attractive forces to the net protein-protein interaction by adjusting solution pH and ionic strength. We find a strong correlation between the concentrated protein solution viscosity and the k_h but this was not observed for the k_D or the b₂₂, which have been previously used as indicators of high concentration viscosity. Trends observed in [η] and k_h values as a function of pH and ionic strength are rationalised in terms of protein-protein interactions.

Protein-polymer mixtures in the colloid limit: Aggregation, sedimentation, and crystallization

While proteins have been treated as particles with a spherically symmetric interaction, of course in reality, the situation is rather more complex. A simple step toward higher complexity is to treat the proteins as non-spherical particles and that is the approach we pursue here. We investigate the phase behavior of the enhanced green fluorescent protein (eGFP) under the addition of a non-adsorbing polymer, polyethylene glycol. From small angle x-ray scattering, we infer that the eGFP undergoes dimerization and we treat the dimers as spherocylinders with aspect ratio L/D - 1 = 1.05. Despite the complex nature of the proteins, we find that the phase behavior is similar to that of hard spherocylinders with an ideal polymer depletant, exhibiting aggregation and, in a small region of the phase diagram, crystallization. By comparing our measurements of the onset of aggregation with predictions for hard colloids and ideal polymers [S. V. Savenko and M. Dijkstra, J. Chem. Phys. 124, 234902 (2006) and Lo Verso et al., Phys. Rev. E 73, 061407 (2006)], we find good agreement, which suggests that the behavior of the eGFP is consistent with that of hard spherocylinders and ideal polymers.

Decorated networks of native proteins: nanomaterials with tunable mesoscopic domain size

Natural and artificial proteins with designer properties and functionalities offer unparalleled opportunity for functional nanoarchitectures formed through self-assembly. However, to exploit this potential we need to design the system such that assembly results in desired architecture forms while avoiding denaturation and therefore retaining protein functionality. Here we address this challenge with a model system of fluorescent proteins. By manipulating self-assembly using techniques inspired by soft matter where interactions between the components are controlled to yield the desired structure, we have developed a methodology to assemble networks of proteins of one species which we can decorate with another, whose coverage we can tune. Consequently, the interfaces between domains of each component can also be tuned, with potential applications for example in energy - or electron - transfer. Our model system of eGFP and mCherry with tuneable interactions reveals control over domain sizes in the resulting networks.

Cluster Percolation Causes Shear Thinning Behavior in Concentrated Solutions of Monoclonal Antibodies

High-concentration (>100 g/L) solutions of monoclonal antibodies (mAbs) are typically characterized by anomalously large solution viscosity and shear thinning behavior for strain rates ≥10³ s^-1. Here, the link between protein-protein interactions (PPIs) and the rheology of concentrated solutions of COE-03 and COE-19 mAbs is studied by means of static and dynamic light scattering and microfluidic rheometry. By comparing the experimental data with predictions based on the Baxter sticky hard-sphere model, we surprisingly find a connection between the observed shear thinning and the predicted percolation threshold. The longest shear relaxation time of mAbs was much larger than that of model sticky hard spheres within the same region of the phase diagram, which is attributed to the anisotropy of the mAb PPIs. Our results suggest that not only the strength but also the patchiness of short-range attractive PPIs should be explicitly accounted for by theoretical approaches aimed at predicting the shear rate-dependent viscosity of dense mAb solutions.

A coarse-grained xDLVO model for colloidal protein-protein interactions

Colloidal protein-protein interactions (PPIs) of attractive and repulsive nature modulate the solubility of proteins, their aggregation, precipitation and crystallization. Such interactions are very important for many biotechnological processes, but are complex and hard to control, therefore, difficult to be understood in terms of measurements alone. In diluted protein solutions, PPIs can be estimated from the osmotic second virial coefficient, B22, which has been calculated using different methods and levels of theory. The most popular approach is based on the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory and its extended versions, i.e. xDLVO. Despite much efforts, these models are not fully quantitative and must be fitted to experiments, which limits their predictive value. Here, we report an extended xDLVO-CG model, which extends existing models by a coarse-grained representation of proteins and the inclusion of an additional ion-protein dispersion interaction term. We demonstrate for four proteins, i.e. lysozyme (LYZ), subtilisin (Subs), bovine serum albumin (BSA) and immunoglobulin (IgG1), that semi-quantitative agreement with experimental values without the need to fit to experimental B22 values. While most likely not the final step in the nearly hundred years of research in PPIs, xDLVO-CG is a step towards predictive PPIs calculations that are transferable to different proteins.

Biophysical Characterization of Binary Therapeutic Monoclonal Antibody Mixtures

Coformulations containing two therapeutic monoclonal antibodies (mAbs) could offer various benefits like enhanced therapeutic efficacy and better patient compliance. However, there are very few published studies on coformulations and binary mixtures of mAbs. It remains unclear to what extent mAbs with different physicochemical properties can be combined in solution without detrimental effects on protein stability. Here, we present a study including six model mAbs of the IgG1 subclass that are commercially available. In silico and biophysical characterization shows that the proteins have different physicochemical properties. Thus, their combinations represent various scenarios for coformulation development. We prepared all possible binary mixtures of the six mAbs and determined several biophysical parameters that are assessed during early-stage protein drug product development. The measured biophysical parameters are indicative of the conformational protein stability (inflection points of the thermal protein unfolding transitions) and the colloidal protein stability (aggregation onset temperatures and interaction parameter k_D from dynamic light scattering). Remarkably, all 15 binary mAb mixtures do not exhibit biophysical parameters that indicate inferior conformational or colloidal stability compared to the least stable mAb in the mixture. Our findings suggest that the coformulation of some therapeutic monoclonal antibodies of the IgG1 subclass could be possible in a straightforward way as severe detrimental effects on the stability of these proteins in binary mixtures were not observed.