1
|
Yuan G, Salipante PF, Hudson SD, Gillilan RE, Huang Q, Hatch HW, Shen VK, Grishaev AV, Pabit S, Upadhya R, Adhikari S, Panchal J, Blanco MA, Liu Y. Flow Activation Energy of High-Concentration Monoclonal Antibody Solutions and Protein-Protein Interactions Influenced by NaCl and Sucrose. Mol Pharm 2024; 21:4553-4564. [PMID: 39163212 DOI: 10.1021/acs.molpharmaceut.4c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/22/2024]
Abstract
The solution viscosity and protein-protein interactions (PPIs) as a function of temperature (4-40 °C) were measured at a series of protein concentrations for a monoclonal antibody (mAb) with different formulation conditions, which include NaCl and sucrose. The flow activation energy (Eη) was extracted from the temperature dependence of solution viscosity using the Arrhenius equation. PPIs were quantified via the protein diffusion interaction parameter (kD) measured by dynamic light scattering, together with the osmotic second virial coefficient and the structure factor obtained through small-angle X-ray scattering. Both viscosity and PPIs were found to vary with the formulation conditions. Adding NaCl introduces an attractive interaction but leads to a significant reduction in the viscosity. However, adding sucrose enhances an overall repulsive effect and leads to a slight decrease in viscosity. Thus, the averaged (attractive or repulsive) PPI information is not a good indicator of viscosity at high protein concentrations for the mAb studied here. Instead, a correlation based on the temperature dependence of viscosity (i.e., Eη) and the temperature sensitivity in PPIs was observed for this specific mAb. When kD is more sensitive to the temperature variation, it corresponds to a larger value of Eη and thus a higher viscosity in concentrated protein solutions. When kD is less sensitive to temperature change, it corresponds to a smaller value of Eη and thus a lower viscosity at high protein concentrations. Rather than the absolute value of PPIs at a given temperature, our results show that the temperature sensitivity of PPIs may be a more useful metric for predicting issues with high viscosity of concentrated solutions. In addition, we also demonstrate that caution is required in choosing a proper protein concentration range to extract kD. In some excipient conditions studied here, the appropriate protein concentration range needs to be less than 4 mg/mL, remarkably lower than the typical concentration range used in the literature.
Collapse
Affiliation(s)
- Guangcui Yuan
- Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Paul F Salipante
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Steven D Hudson
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Richard E Gillilan
- Center for High-Energy X-ray Sciences at CHESS, Cornell University, Ithaca, New York 14853, United States
| | - Qingqiu Huang
- Center for High-Energy X-ray Sciences at CHESS, Cornell University, Ithaca, New York 14853, United States
| | - Harold W Hatch
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Vincent K Shen
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Alexander V Grishaev
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Suzette Pabit
- Analytical Enabling Capabilities, Merck & Co., Inc., Rahway, New Jersey 07065, United States
| | - Rahul Upadhya
- Analytical Enabling Capabilities, Merck & Co., Inc., Rahway, New Jersey 07065, United States
| | - Sudeep Adhikari
- Analytical Enabling Capabilities, Merck & Co., Inc., Rahway, New Jersey 07065, United States
| | - Jainik Panchal
- Sterile and Specialty Products, Merck & Co., Inc., Kenilworth, New Jersey 07033, United States
| | - Marco A Blanco
- Discovery Pharmaceutical Sciences, Merck & Co., Inc., West Point, Pennsylvania 19486, United States
| | - Yun Liu
- Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States
| |
Collapse
|
2
|
Polimeni M, Zaccarelli E, Gulotta A, Lund M, Stradner A, Schurtenberger P. A multi-scale numerical approach to study monoclonal antibodies in solution. APL Bioeng 2024; 8:016111. [PMID: 38425712 PMCID: PMC10902793 DOI: 10.1063/5.0186642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 01/30/2024] [Indexed: 03/02/2024] Open
Abstract
Developing efficient and robust computational models is essential to improve our understanding of protein solution behavior. This becomes particularly important to tackle the high-concentration regime. In this context, the main challenge is to put forward coarse-grained descriptions able to reduce the level of detail, while retaining key features and relevant information. In this work, we develop an efficient strategy that can be used to investigate and gain insight into monoclonal antibody solutions under different conditions. We use a multi-scale numerical approach, which connects information obtained at all-atom and amino-acid levels to bead models. The latter has the advantage of reproducing the properties of interest while being computationally much faster. Indeed, these models allow us to perform many-protein simulations with a large number of molecules. We can, thus, explore conditions not easily accessible with more detailed descriptions, perform effective comparisons with experimental data up to very high protein concentrations, and efficiently investigate protein-protein interactions and their role in phase behavior and protein self-assembly. Here, a particular emphasis is given to the effects of charges at different ionic strengths.
Collapse
Affiliation(s)
- Marco Polimeni
- Division of Physical Chemistry, Lund University, Lund, Sweden
| | - Emanuela Zaccarelli
- Institute for Complex Systems, National Research Council (ISC-CNR), Piazzale Aldo Moro 5, 00185 Rome, Italy
| | | | - Mikael Lund
- Division of Computational Chemistry, Lund University, Lund, Sweden
| | - Anna Stradner
- Division of Physical Chemistry, Lund University, Lund, Sweden
| | | |
Collapse
|
3
|
Lou H, Zhang Y, Kuczera K, Hageman MJ, Schöneich C. Molecular Dynamics Simulation of an Iron(III) Binding Site on the Fc Domain of IgG1 Relevant for Visible Light-Induced Protein Fragmentation. Mol Pharm 2024; 21:501-512. [PMID: 38128475 DOI: 10.1021/acs.molpharmaceut.3c00612] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Molecular dynamics simulations were employed to investigate the interaction between Fe(III) and an iron-binding site composed of THR259, ASP252, and GLU261 on the Fc domain of an IgG1. The goal was to provide microscopic mechanistic information for the photochemical, iron-dependent site-specific oxidative fragmentation of IgG1 at THR259 reported in our previous paper. The distance between Fe(III) and residues of interest as well as the binding pocket size was examined for both protonated and deprotonated THR259. The Fe(III) binding free energy (ΔG) was estimated by using an umbrella sampling approach. The pKa shift of the THR259 hydroxyl group caused by the presence of nearby Fe(III) was estimated based on a thermodynamic cycle. The simulation results show that Fe(III) resides inside the proposed binding pocket and profoundly changes the pocket configuration. The ΔG values indicate that the pocket possesses a strong binding affinity for Fe(III). Furthermore, Fe(III) profoundly lowers the pKa value of the THR259 hydroxyl group by 5.4 pKa units.
Collapse
Affiliation(s)
- Hao Lou
- Department of Pharmaceutical Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
- Biopharmaceutical Innovation and Optimization Center, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yilue Zhang
- Department of Pharmaceutical Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| | - Krzysztof Kuczera
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
- Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66045, United States
| | - Michael J Hageman
- Department of Pharmaceutical Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
- Biopharmaceutical Innovation and Optimization Center, University of Kansas, Lawrence, Kansas 66047, United States
| | - Christian Schöneich
- Department of Pharmaceutical Chemistry, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
4
|
Brudar S, Breydo L, Chung E, Dill KA, Ehterami N, Phadnis K, Senapati S, Shameem M, Tang X, Tayyab M, Hribar-Lee B. Antibody association in solution: cluster distributions and mechanisms. MAbs 2024; 16:2339582. [PMID: 38666507 PMCID: PMC11057677 DOI: 10.1080/19420862.2024.2339582] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Accepted: 04/02/2024] [Indexed: 05/01/2024] Open
Abstract
Understanding factors that affect the clustering and association of antibodies molecules in solution is critical to their development as therapeutics. For 19 different monoclonal antibody (mAb) solutions, we measured the viscosities, the second virial coefficients, the Kirkwood-Buff integrals, and the cluster distributions of the antibody molecules as functions of protein concentration. Solutions were modeled using the statistical-physics Wertheim liquid-solution theory, representing antibodies as Y-shaped molecular structures of seven beads each. We found that high-viscosity solutions result from more antibody molecules per cluster. Multi-body properties such as viscosity are well predicted experimentally by the 2-body Kirkwood-Buff quantity, G22, but not by the second virial coefficient, B22, and well-predicted theoretically from the Wertheim protein-protein sticking energy. Weakly interacting antibodies are rate-limited by nucleation; strongly interacting ones by propagation. This approach gives a way to relate micro to macro properties of solutions of associating proteins.
Collapse
Affiliation(s)
- Sandi Brudar
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
| | - Leonid Breydo
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Elisha Chung
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
- Department of Chemistry and Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY, USA
| | - Nasim Ehterami
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Ketan Phadnis
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Samir Senapati
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Mohammed Shameem
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Xiaolin Tang
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Muhammmad Tayyab
- Formulation Development Group, Regeneron Pharmaceuticals, Tarrytown, NY, USA
| | - Barbara Hribar-Lee
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Ljubljana, Slovenia
| |
Collapse
|
5
|
Vlachy V, Kalyuzhnyi YV, Hribar-Lee B, Dill KA. Protein Association in Solution: Statistical Mechanical Modeling. Biomolecules 2023; 13:1703. [PMID: 38136574 PMCID: PMC10742237 DOI: 10.3390/biom13121703] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 11/16/2023] [Accepted: 11/20/2023] [Indexed: 12/24/2023] Open
Abstract
Protein molecules associate in solution, often in clusters beyond pairwise, leading to liquid phase separations and high viscosities. It is often impractical to study these multi-protein systems by atomistic computer simulations, particularly in multi-component solvents. Instead, their forces and states can be studied by liquid state statistical mechanics. However, past such approaches, such as the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory, were limited to modeling proteins as spheres, and contained no microscopic structure-property relations. Recently, this limitation has been partly overcome by bringing the powerful Wertheim theory of associating molecules to bear on protein association equilibria. Here, we review these developments.
Collapse
Affiliation(s)
- Vojko Vlachy
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000 Ljubljana, Slovenia;
| | | | - Barbara Hribar-Lee
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, 1000 Ljubljana, Slovenia;
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, New York, NY 11794, USA;
- Department of Chemistry, Physics and Astronomy, Stony Brook University, New York, NY 11790, USA
| |
Collapse
|
6
|
Chowdhury AA, Manohar N, Lanzaro A, Kimball WD, Witek MA, Woldeyes MA, Majumdar R, Qian KK, Xu S, Gillilan RE, Huang Q, Truskett TM, Johnston KP. Characterizing Protein-Protein Interactions and Viscosity of a Monoclonal Antibody from Low to High Concentration Using Small-Angle X-ray Scattering and Molecular Dynamics Simulations. Mol Pharm 2023; 20:5563-5578. [PMID: 37782765 DOI: 10.1021/acs.molpharmaceut.3c00484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Understanding protein-protein interactions and formation of reversible oligomers (clusters) in concentrated monoclonal antibody (mAb) solutions is necessary for designing stable, low viscosity (η) concentrated formulations for processing and subcutaneous injection. Here we characterize the strength (K) of short-range anisotropic attractions (SRA) for 75-200 mg/mL mAb2 solutions at different pH and cosolute conditions by analyzing structure factors (Seff(q)) from small-angle X-ray scattering (SAXS) using coarse-grained molecular dynamics simulations. Best fit simulations additionally provide cluster size distributions, fractal dimensions, cluster occluded volume, and mAb coordination numbers. These equilibrium properties are utilized in a model to account for increases in viscosity caused by occluded volume in the clusters (packing effects) and dissipation of stress across lubricated fractal clusters. Seff(q) is highly sensitive to K at 75 mg/mL where mAbs can mutually align to form SRA contacts but becomes less sensitive at 200 mg/mL as steric repulsion due to packing becomes dominant. In contrast, η at 200 mg/mL is highly sensitive to SRA and the average cluster size from SAXS/simulation, which is observed to track the cluster relaxation time from shear thinning. By analyzing the distribution of sub-bead hot spots on the 3D mAb surface, we identify a strongly attractive hydrophobic patch in the complementarity determining region (CDR) at pH 4.5 that contributes to the high K and consequently large cluster sizes and high η. Adding NaCl screens electrostatic interactions and increases the impact of hydrophobic attraction on cluster size and raises η, whereas nonspecific binding of Arg attenuates all SRA, reducing η. The hydrophobic patch is absent at higher pH values, leading to smaller K, smaller clusters, and lower η. This work constitutes a first attempt to use SAXS and CG modeling to link both structural and rheological properties of concentrated mAb solutions to the energetics of specific hydrophobic patches on mAb surfaces. As such, our work opens an avenue for future research, including the possibility of designing coarse-grained models with physically meaningful interacting hot spots.
Collapse
Affiliation(s)
- Amjad A Chowdhury
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Neha Manohar
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Alfredo Lanzaro
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - William D Kimball
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Marta A Witek
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | | | - Ranajoy Majumdar
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | - Ken K Qian
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | - Shifeng Xu
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Richard E Gillilan
- Center for High Energy X-ray Sciences at CHESS, Cornell University, Ithaca, New York 14853, United States
| | - Qingqiu Huang
- Center for High Energy X-ray Sciences at CHESS, Cornell University, Ithaca, New York 14853, United States
| | - Thomas M Truskett
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Department of Physics, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Keith P Johnston
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
7
|
Xu AY, Blanco MA, Castellanos MM, Meuse CW, Mattison K, Karageorgos I, Hatch HW, Shen VK, Curtis JE. Role of Domain-Domain Interactions on the Self-Association and Physical Stability of Monoclonal Antibodies: Effect of pH and Salt. J Phys Chem B 2023; 127:8344-8357. [PMID: 37751332 PMCID: PMC10561141 DOI: 10.1021/acs.jpcb.3c03928] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Revised: 09/11/2023] [Indexed: 09/28/2023]
Abstract
Monoclonal antibodies (mAbs) make up a major class of biotherapeutics with a wide range of clinical applications. Their physical stability can be affected by various environmental factors. For instance, an acidic pH can be encountered during different stages of the mAb manufacturing process, including purification and storage. Therefore, understanding the behavior of flexible mAb molecules in acidic solution environments will benefit the development of stable mAb products. This study used small-angle X-ray scattering (SAXS) and complementary biophysical characterization techniques to investigate the conformational flexibility and protein-protein interactions (PPI) of a model mAb molecule under near-neutral and acidic conditions. The study also characterized the interactions between Fab and Fc fragments under the same buffer conditions to identify domain-domain interactions. The results suggest that solution pH significantly influences mAb flexibility and thus could help mAbs remain physically stable by maximizing local electrostatic repulsions when mAbs become crowded in solution. Under acidic buffer conditions, both Fab and Fc contribute to the repulsive PPI observed among the full mAb at a low ionic strength. However, as ionic strength increases, hydrophobic interactions lead to the self-association of Fc fragments and, subsequently, could affect the aggregation state of the mAb.
Collapse
Affiliation(s)
- Amy Y. Xu
- Department
of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Marco A. Blanco
- Discovery
Pharmaceutical Sciences, Merck Research
Laboratories, Merck & Co., Inc, West Point, Pennsylvania 19486, United States
| | - Maria Monica Castellanos
- Institute
for Bioscience and Biotechnology Research, University of Maryland, Rockville, Maryland 20850, United States
- NIST
Center for Neutron Research, National Institute
of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Curtis W. Meuse
- Institute
for Bioscience and Biotechnology Research, University of Maryland, Rockville, Maryland 20850, United States
- Biomolecular
Measurement Division, National Institute
of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Kevin Mattison
- Malvern
Panalytical, Westborough, Massachusetts 01581, United States
| | - Ioannis Karageorgos
- Institute
for Bioscience and Biotechnology Research, University of Maryland, Rockville, Maryland 20850, United States
- Biomolecular
Measurement Division, National Institute
of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Harold W. Hatch
- Chemical
Sciences Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Vincent K. Shen
- Chemical
Sciences Division, Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, United States
| | - Joseph E. Curtis
- NIST
Center for Neutron Research, National Institute
of Standards and Technology, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
8
|
Mosca I, Pounot K, Beck C, Colin L, Matsarskaia O, Grapentin C, Seydel T, Schreiber F. Biophysical Determinants for the Viscosity of Concentrated Monoclonal Antibody Solutions. Mol Pharm 2023; 20:4698-4713. [PMID: 37549226 DOI: 10.1021/acs.molpharmaceut.3c00440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/09/2023]
Abstract
Monoclonal antibodies (mAbs) are particularly relevant for therapeutics due to their high specificity and versatility, and mAb-based drugs are hence used to treat numerous diseases. The increased patient compliance of self-administration motivates the formulation of products for subcutaneous (SC) administration. The associated challenge is to formulate highly concentrated antibody solutions to achieve a significant therapeutic effect, while limiting their viscosity and preserving their physicochemical stability. Protein-protein interactions (PPIs) are in fact the root cause of several potential problems concerning the stability, manufacturability, and delivery of a drug product. The understanding of macroscopic viscosity requires an in-depth knowledge on protein diffusion, PPIs, and self-association/aggregation. Here, we study the self-diffusion of different mAbs of the IgG1 subtype in aqueous solution as a function of the concentration and temperature by quasi-elastic neutron scattering (QENS). QENS allows us to probe the short-time self-diffusion of the molecules and therefore to determine the hydrodynamic mAb cluster size and to gain information on the internal mAb dynamics. Small-angle neutron scattering (SANS) is jointly employed to probe structural details and to understand the nature and intensity of PPIs. Complementary information is provided by molecular dynamics (MD) simulations and viscometry, thus obtaining a comprehensive picture of mAb diffusion.
Collapse
Affiliation(s)
- Ilaria Mosca
- Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, Tübingen 72076, Germany
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | - Kévin Pounot
- Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, Tübingen 72076, Germany
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | - Christian Beck
- Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, Tübingen 72076, Germany
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | - Louise Colin
- Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, Tübingen 72076, Germany
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | - Olga Matsarskaia
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | | | - Tilo Seydel
- Institut Max von Laue - Paul Langevin, 71 Av. des Martyrs, Grenoble 38042, France
| | - Frank Schreiber
- Institut für Angewandte Physik, Universität Tübingen, Auf der Morgenstelle 10, Tübingen 72076, Germany
| |
Collapse
|
9
|
Chowdhury AA, Manohar N, Witek MA, Woldeyes MA, Majumdar R, Qian KK, Kimball WD, Xu S, Lanzaro A, Truskett TM, Johnston KP. Subclass Effects on Self-Association and Viscosity of Monoclonal Antibodies at High Concentrations. Mol Pharm 2023. [PMID: 37191356 DOI: 10.1021/acs.molpharmaceut.3c00023] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
The effects of a subclass of monoclonal antibodies (mAbs) on protein-protein interactions, formation of reversible oligomers (clusters), and viscosity (η) are not well understood at high concentrations. Herein, we quantify a short-range anisotropic attraction between the complementarity-determining region (CDR) and CH3 domains (KCDR-CH3) for vedolizumab IgG1, IgG2, or IgG4 subclasses by fitting small-angle X-ray scattering (SAXS) structure factor Seff(q) data with an extensive library of 12-bead coarse-grained (CG) molecular dynamics simulations. The KCDR-CH3 bead attraction strength was isolated from the strength of long-range electrostatic repulsion for the full mAb, which was determined from the theoretical net charge and a scaling parameter ψ to account for solvent accessibility and ion pairing. At low ionic strength (IS), the strongest short-range attraction (KCDR-CH3) and consequently the largest clusters and highest η were observed with IgG1, the subclass with the most positively charged CH3 domain. Furthermore, the trend in KCDR-CH3 with the subclass followed the electrostatic interaction energy between the CDR and CH3 regions calculated with the BioLuminate software using the 3D mAb structure and molecular interaction potentials. Whereas the equilibrium cluster size distributions and fractal dimensions were determined from fits of SAXS with the MD simulations, the degree of cluster rigidity under flow was estimated from the experimental η with a phenomenological model. For the systems with the largest clusters, especially IgG1, the inefficient packing of mAbs in the clusters played the largest role in increasing η, whereas for other systems, the relative contribution from stress produced by the clusters was more significant. The ability to relate η to short-range attraction from SAXS measurements at high concentrations and to theoretical characterization of electrostatic patches on the 3D surface is not only of fundamental interest but also of practical value for mAb discovery, processing, formulation, and subcutaneous delivery.
Collapse
Affiliation(s)
- Amjad A Chowdhury
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Neha Manohar
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Marta A Witek
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | | | - Ranajoy Majumdar
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | - Ken K Qian
- Eli Lilly and Company, Indianapolis, Indiana 46225, United States
| | - William D Kimball
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Shifeng Xu
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Alfredo Lanzaro
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Thomas M Truskett
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Department of Physics, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Keith P Johnston
- McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|