1
|
ATP and Tri-Polyphosphate (TPP) Suppress Protein Aggregate Growth by a Supercharging Mechanism. Biomedicines 2021; 9:biomedicines9111646. [PMID: 34829875 PMCID: PMC8616003 DOI: 10.3390/biomedicines9111646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Revised: 10/27/2021] [Accepted: 11/03/2021] [Indexed: 11/17/2022] Open
Abstract
A common strategy to increase aggregation resistance is through rational mutagenesis to supercharge proteins, which leads to high colloidal stability, but often has the undesirable effect of lowering conformational stability. We show this trade-off can be overcome by using small multivalent polyphosphate ions, adenosine triphosphate (ATP) and tripolyphosphate (TPP) as excipients. These ions are equally effective at suppressing aggregation of ovalbumin and bovine serum albumin (BSA) upon thermal stress as monitored by dynamic and static light scattering. Monomer loss kinetic studies, combined with measurements of native state protein–protein interactions and ζ-potentials, indicate the ions reduce aggregate growth by increasing the protein colloidal stability through binding and overcharging the protein. Out of three additional proteins studied, ribonuclease A (RNaseA), α-chymotrypsinogen (α-Cgn), and lysozyme, we only observed a reduction in aggregate growth for RNaseA, although overcharging by the poly-phosphate ions still occurs for lysozyme and α-Cgn. Because the salts do not alter protein conformational stability, using them as excipients could be a promising strategy for stabilizing biopharmaceuticals once the protein structural factors that determine whether multivalent ion binding will increase colloidal stability are better elucidated. Our findings also have biological implications. Recently, it has been proposed that ATP also plays an important role in maintaining intracellular biological condensates and preventing protein aggregation in densely packed cellular environments. We expect electrostatic interactions are a significant factor in determining the stabilizing ability of ATP towards maintaining proteins in non-dispersed states in vivo.
Collapse
|
2
|
Cheng R, Li J, Ríos de Anda I, Taylor TWC, Faers MA, Anderson JLR, Seddon AM, Royall CP. Protein-polymer mixtures in the colloid limit: Aggregation, sedimentation, and crystallization. J Chem Phys 2021; 155:114901. [PMID: 34551522 DOI: 10.1063/5.0052122] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Rui Cheng
- HH Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | - Jingwen Li
- HH Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | | | - Thomas W C Taylor
- HH Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | | | - J L Ross Anderson
- School of Biochemistry, University of Bristol, Bristol BS8 1TD, United Kingdom
| | - Annela M Seddon
- HH Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| | - C Patrick Royall
- HH Wills Physics Laboratory, Tyndall Avenue, Bristol BS8 1TL, United Kingdom
| |
Collapse
|
3
|
Ríos de Anda I, Coutable-Pennarun A, Brasnett C, Whitelam S, Seddon A, Russo J, Anderson JLR, Royall CP. Decorated networks of native proteins: nanomaterials with tunable mesoscopic domain size. SOFT MATTER 2021; 17:6873-6883. [PMID: 34231559 PMCID: PMC8294043 DOI: 10.1039/d0sm02269a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2020] [Accepted: 06/30/2021] [Indexed: 06/13/2023]
Abstract
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.
Collapse
Affiliation(s)
- Ioatzin Ríos de Anda
- H.H. Wills Physics LaboratoryTyndall AvenueBristolBS8 1TLUK
- School of Mathematics, University WalkBristolBS8 1TWUK
| | - Angélique Coutable-Pennarun
- BrisSynBio Synthetic Biology Research Centre, Life Sciences BuildingTyndall AvenueBristolBS8 1TQUK
- School of Biochemistry, University of BristolBristolBS8 1TDUK
| | | | - Stephen Whitelam
- Molecular Foundry, Lawrence Berkeley National LaboratoryBerkeleyCalifornia 94720USA
| | - Annela Seddon
- H.H. Wills Physics LaboratoryTyndall AvenueBristolBS8 1TLUK
- Bristol Centre for Functional Nanomaterials, University of BristolBristolBS8 1TLUK
| | - John Russo
- School of Mathematics, University WalkBristolBS8 1TWUK
- Dipartimento di Fisica and CNR-ISC, Sapienza-Università di RomaPiazzale A. Moro 200185 RomaItaly
| | - J. L. Ross Anderson
- School of Biochemistry, University of BristolBristolBS8 1TDUK
- School of Cellular and Molecular Medicine, University WalkBristolBS8 1TDUK
| | - C. Patrick Royall
- H.H. Wills Physics LaboratoryTyndall AvenueBristolBS8 1TLUK
- Gulliver UMR CNRS 7083, ESPCI Paris, Université PSL75005 ParisFrance
- School of Chemistry, University of BristolCantock's CloseBristolBS8 1TSUK
- Centre for Nanoscience and Quantum InformationTyndall AvenueBristolBS8 1FDUK
| |
Collapse
|
4
|
Lanzaro A, Roche A, Sibanda N, Corbett D, Davis P, Shah M, Pathak JA, Uddin S, van der Walle CF, Yuan XF, Pluen A, Curtis R. Cluster Percolation Causes Shear Thinning Behavior in Concentrated Solutions of Monoclonal Antibodies. Mol Pharm 2021; 18:2669-2682. [PMID: 34121411 DOI: 10.1021/acs.molpharmaceut.1c00198] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
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 ≥103 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.
Collapse
Affiliation(s)
- Alfredo Lanzaro
- Institute for Systems Rheology, Guangzhou University, No. 230 West Outer Ring Road, Guangzhou Higher Education Mega Center, Guangzhou 510006, China
| | - Aisling Roche
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Nicole Sibanda
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Daniel Corbett
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| | - Peter Davis
- Department of Molecular Biology and Biotechnology, University of Sheffield UK, Sheffield S10 2TN, United Kingdom
| | - Maryam Shah
- Division of Pharmacy and Optometry, School of Health Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, M13 9PT, United Kingdom
| | - Jai A Pathak
- Dosage Form Design and Development, Biopharmaceuticals Development, R&D, AstraZeneca, Cambridge CB21 6GH, United Kingdom
| | - Shahid Uddin
- Dosage Form Design and Development, Biopharmaceuticals Development, R&D, AstraZeneca, Cambridge CB21 6GH, United Kingdom
| | - Christopher F van der Walle
- Dosage Form Design and Development, Biopharmaceuticals Development, R&D, AstraZeneca, Cambridge CB21 6GH, United Kingdom
| | - Xue-Feng Yuan
- Institute for Systems Rheology, Guangzhou University, No. 230 West Outer Ring Road, Guangzhou Higher Education Mega Center, Guangzhou 510006, China
| | - Alain Pluen
- Division of Pharmacy and Optometry, School of Health Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, M13 9PT, United Kingdom
| | - Robin Curtis
- School of Chemical Engineering and Analytical Science, Manchester Institute of Biotechnology, University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom
| |
Collapse
|
5
|
Pusara S, Yamin P, Wenzel W, Krstić M, Kozlowska M. A coarse-grained xDLVO model for colloidal protein-protein interactions. Phys Chem Chem Phys 2021; 23:12780-12794. [PMID: 34048523 DOI: 10.1039/d1cp01573g] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
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.
Collapse
Affiliation(s)
- Srdjan Pusara
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
| | - Peyman Yamin
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
| | - Wolfgang Wenzel
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
| | - Marjan Krstić
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany. and Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Mariana Kozlowska
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), Hermann-von-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.
| |
Collapse
|
6
|
Zhou M, Yan Z, Li H, Liu X, Sun P. Application of Affinity-Capture Self-Interaction Nanoparticle Spectroscopy in Predicting Protein Stability, Especially for Co-Formulated Antibodies. Pharm Res 2021; 38:721-732. [PMID: 33754257 DOI: 10.1007/s11095-021-03026-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2020] [Accepted: 03/01/2021] [Indexed: 12/15/2022]
Abstract
PURPOSE From traditional monoclonal antibodies to more and more complex mAb-based formulations, biopharmaceutical faces one challenge after another. To avoid these issues, identification of therapeutic proteins in the initial discovery process that has high stability and low self-interaction would simplify the development of safe and effective antibody therapeutics. METHOD Affinity-capture self-interaction nanoparticle spectroscopy (AC-SINS) is a new prediction method capable of identifying mAbs with different self-association propensity. In this study, 10 formulated monoclonal antibody (mAb) therapeutics include different mAb isotypes and co-formulated antibodies were measured by AC-SINS and some biophysical methods to predict protein stability. The prediction results of all 10 mAbs were compared to their stability data (Δ%monomer and Δ%HMWs) at accelerated (25°C and 40°C) and long-term storage conditions (4°C) as measured by size exclusion chromatography. RESULT AC-SINS method has a good predictive correlation with each mAbs and co-formulated antibodies. There were no physicochemical, intermolecular, or biological interactions that occurred between the two components of co-formulated antibodies which confirmed by Analytical ultracentrifugation (AUC). CONCLUSION Here we discuss the correlation between each method and protein stability, and also use AC-SINS assay to predict the stability of co-formulated antibodies for the first time. This may be an effective way to predict the stability of these complex mAb-based formulations such as co-formulated mAbs.
Collapse
Affiliation(s)
- Meng Zhou
- School of Pharmacy, China Pharmaceutical University, Nanjing, 211198, China
| | - Zhen Yan
- Shanghai Hengrui Pharmaceutical Co., Ltd., Shanghai, 200245, China
| | - Hao Li
- Shanghai Hengrui Pharmaceutical Co., Ltd., Shanghai, 200245, China
| | - Xun Liu
- Shanghai Hengrui Pharmaceutical Co., Ltd., Shanghai, 200245, China.
| | - Piaoyang Sun
- School of Pharmacy, China Pharmaceutical University, Nanjing, 211198, China. .,Shanghai Hengrui Pharmaceutical Co., Ltd., Shanghai, 200245, China.
| |
Collapse
|
7
|
Bakhshandeh A, Frydel D, Levin Y. Charge regulation of colloidal particles in aqueous solutions. Phys Chem Chem Phys 2020; 22:24712-24728. [PMID: 33104140 DOI: 10.1039/d0cp03633a] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
We study the charge regulation of colloidal particles inside aqueous electrolyte solutions. To stabilize a colloidal suspension against precipitation, colloidal particles are synthesized with either acidic or basic groups on their surface. On contact with water, these surface groups undergo proton transfer reactions, resulting in colloidal surface charge. The charge is determined by the condition of local chemical equilibrium between hydronium ions inside the solution and at the colloidal surface. We use a model of Baxter sticky spheres to explicitly calculate the equilibrium dissociation constants and to construct a theory which is able to quantitatively predict the effective charge of colloidal particles with either acidic or basic surface groups. The predictions of the theory for the model are found to be in excellent agreement with the results of Monte Carlo simulations. This theory is further extended to treat colloidal particles with a mixture of both acidic and basic surface groups.
Collapse
Affiliation(s)
- Amin Bakhshandeh
- Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 15051, CEP 91501-970, Porto Alegre, RS, Brazil.
| | | | | |
Collapse
|
8
|
Srivastava R, Chattopadhyaya M, Bandyopadhyay P. Calculation of salt-dependent free energy of binding of β-lactoglobulin homodimer formation and mechanism of dimer formation using molecular dynamics simulation and three-dimensional reference interaction site model (3D-RISM): diffuse salt ions and non-polar interactions between the monomers favor the dimer formation. Phys Chem Chem Phys 2020; 22:2142-2156. [PMID: 31912070 DOI: 10.1039/c9cp05578a] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
There are several important phenomena in chemistry, biology, and physics where molecules (or parts of a molecule) having charges of the same sign come closer together and become stable. DNA condensation, RNA folding, colloid-colloid interactions are some of the examples of this kind. In the current work, we have investigated how β-lactoglobulin, a protein found in milk, in spite of carrying +13 charge, favors the homodimer formation in the presence of salt. We have focussed on calculating the protein-protein binding free energy in the presence of salt and identifying the thermodynamic and microscopic mechanism of the process. Estimation of binding free energy of this salt-dependent process is done by combining molecular dynamics simulation with statistical mechanical theory of three-dimensional reference interaction site model (3D-RISM). Binding free energy is evaluated from the chemical potential of the solutes as opposed to potential of mean force calculation, which gives only a constrained free energy. Our calculated values semi-quantitatively match with the experimental results. By examining the different components of binding free energy, we have found that the role of salt ions (especially of Cl-) is to shift the equilibrium towards the dimer. Non-polar (Lennard-Jones) interactions between the monomers is also favorable to the binding free energy. However, water slightly disfavors the dimer formation. For the microscopic mechanism, heterogeneous of both Na+ and Cl- near the charged residues at the binding interface and change of this charge distribution on dimer formation contribute to the stability. A fine-tuning of enthalpic and entropic effects of salt ions is found to operate at different salt concentrations. Both thermodynamic and microscopic mechanism of dimer formation gives detailed insight into the complex electrostatics of charged protein-protein binding.
Collapse
Affiliation(s)
- Rakesh Srivastava
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India.
| | | | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India.
| |
Collapse
|
9
|
Singh P, Roche A, van der Walle CF, Uddin S, Du J, Warwicker J, Pluen A, Curtis R. Determination of Protein-Protein Interactions in a Mixture of Two Monoclonal Antibodies. Mol Pharm 2019; 16:4775-4786. [PMID: 31613625 DOI: 10.1021/acs.molpharmaceut.9b00430] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The coformulation of monoclonal antibody (mAb) mixtures provides an attractive route to achieving therapeutic efficacy where the targeting of multiple epitopes is necessary. Controlling and predicting the behavior of such mixtures requires elucidating the molecular basis for the self- and cross-protein-protein interactions and how they depend on solution variables. While self-interactions are now beginning to be well understood, systematic studies of cross-interactions between mAbs in solution do not exist. Here, we have used static light scattering to measure the set of self- and cross-osmotic second virial coefficients in a solution containing a mixture of two mAbs, mAbA and mAbB, as a function of ionic strength and pH. mAbB exhibits strong association at a low ionic strength, which is attributed to an electrostatic attraction that is enhanced by the presence of a strong short-ranged attraction of nonelectrostatic origin. Under all solution conditions, the measured cross-interactions are intermediate self-interactions and follow similar patterns of behavior. There is a strong electrostatic attraction at higher pH values, reflecting the behavior of mAbB. Protein-protein interactions become more attractive with an increasing pH due to reducing the overall protein net charges, an effect that is attenuated with an increasing ionic strength due to the screening of electrostatic interactions. Under moderate ionic strength conditions, the reduced cross-virial coefficient, which reflects only the energetic contribution to protein-protein interactions, is given by a geometric average of the corresponding self-coefficients. We show the relationship can be rationalized using a patchy sphere model, where the interaction energy between sites i and j is given by the arithmetic mean of the i-i and j-j interactions. The geometric mean does not necessarily apply to all mAb mixtures and is expected to break down at a lower ionic strength due to the nonadditivity of electrostatic interactions.
Collapse
Affiliation(s)
- Priyanka Singh
- Manchester Pharmacy School , University of Manchester , Manchester M13 9PL , United Kingdom
| | - Aisling Roche
- School of Chemical Engineering and Analytical Science , University of Manchester , Manchester M1 7DN , United Kingdom
| | - Christopher F van der Walle
- School of Chemical Engineering and Analytical Science , University of Manchester , Manchester M1 7DN , United Kingdom.,Dosage Form Design & Development , AstraZeneca , Granta Park , Cambridge CB21 6GH , United Kingdom
| | - Shahid Uddin
- Formulation Sciences CMC , Immunocore , Milton Park , Abingdon OX14 4RW , United Kingdom
| | - Jiali Du
- Dosage Form Design & Development , AstraZeneca , Gaithersburg MD20878 , United States
| | - Jim Warwicker
- School of Chemistry , University of Manchester , Manchester M1 7DN , United Kingdom
| | - Alain Pluen
- Manchester Pharmacy School , University of Manchester , Manchester M13 9PL , United Kingdom
| | - Robin Curtis
- School of Chemical Engineering and Analytical Science , University of Manchester , Manchester M1 7DN , United Kingdom
| |
Collapse
|
10
|
Qin S, Zhou HX. Calculation of Second Virial Coefficients of Atomistic Proteins Using Fast Fourier Transform. J Phys Chem B 2019; 123:8203-8215. [PMID: 31490691 DOI: 10.1021/acs.jpcb.9b06808] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The second virial coefficient, B2, measures a protein solution's deviation from ideal behavior. It is widely used to predict or explain solubility, crystallization condition, aggregation propensity, and critical temperature for liquid-liquid phase separation. B2 is determined by the interaction energy between two protein molecules and, specifically, by the integration of the Mayer f-function in the relative configurational space (translation and rotation) of the two molecules. Simple theoretical models, such as one attributed to Derjaguin, Landau, Verwey, and Overbeek (DLVO), can fit the dependence of B2 on salt concentrations. However, model parameters derived often are physically unrealistic and hardly transferable from protein to protein. Previous B2 calculations incorporating atomistic details were done with limited sampling in the configurational space, due to enormous computational cost. Our FMAP method, based on fast Fourier transform, can considerably accelerate such calculations, and here we adapt it to calculate B2 values for proteins represented at the atomic level in implicit solvent. After tuning of a single parameter in the energy function, FMAPB2 predicts well the B2 values for lysozyme and other proteins over wide ranges of solvent conditions (salt concentration, pH, and temperature). The method is available as a web server at http://pipe.rcc.fsu.edu/fmapb2 .
Collapse
Affiliation(s)
- Sanbo Qin
- Department of Chemistry and Department of Physics , University of Illinois at Chicago , Chicago , Illinois 60607 , United States.,Institute of Molecular Biophysics , Florida State University , Tallahassee , Florida 32306 , United States
| | - Huan-Xiang Zhou
- Department of Chemistry and Department of Physics , University of Illinois at Chicago , Chicago , Illinois 60607 , United States
| |
Collapse
|
11
|
Kim DM, Yao X, Vanam RP, Marlow MS. Measuring the effects of macromolecular crowding on antibody function with biolayer interferometry. MAbs 2019; 11:1319-1330. [PMID: 31401928 PMCID: PMC6748605 DOI: 10.1080/19420862.2019.1647744] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Biotherapeutic proteins are commonly dosed at high concentrations into the blood, which is an inherently complex, crowded solution with substantial protein content. The effects of macromolecular crowding may lead to an appreciable level of non-specific hetero-association in this physiological environment. Therefore, developing a method to characterize the diverse consequences of non-specific interactions between proteins under such non-ideal, crowded conditions, which deviate substantially from those commonly employed for in vitro characterization, is vital to achieving a more complete picture of antibody function in a biological context. In this study, we investigated non-specific interactions between human serum albumin (HSA) and two monoclonal antibodies (mAbs) by static light scattering and determined these interactions are both ionic strength-dependent and mAb-dependent. Using biolayer interferometry (BLI), we assessed the effect of HSA on antigen binding by mAbs, demonstrating that these non-specific interactions have a functional impact on mAb:antigen interactions, particularly at low ionic strength. While this effect is mitigated at physiological ionic strength, our in vitro data support the notion that HSA in the blood may lead to non-specific interactions with mAbs in vivo, with a potential impact on their interactions with antigen. Furthermore, the BLI method offers a high-throughput advantage compared to orthogonal techniques such as analytical ultracentrifugation and is amenable to a greater variety of solution conditions compared to nuclear magnetic resonance spectroscopy. Our study demonstrates that BLI is a viable technology for examining the impact of non-specific interactions on specific biologically relevant interactions, providing a direct method to assess binding events in crowded conditions.
Collapse
Affiliation(s)
- Dorothy M Kim
- Pre-Clinical Development and Protein Chemistry, Regeneron Pharmaceuticals, Inc ., Tarrytown , NY , USA
| | - Xiao Yao
- Pre-Clinical Development and Protein Chemistry, Regeneron Pharmaceuticals, Inc ., Tarrytown , NY , USA
| | - Ram P Vanam
- Pre-Clinical Development and Protein Chemistry, Regeneron Pharmaceuticals, Inc ., Tarrytown , NY , USA
| | - Michael S Marlow
- Pre-Clinical Development and Protein Chemistry, Regeneron Pharmaceuticals, Inc ., Tarrytown , NY , USA.,Biotherapeutics Discovery, Boehringer Ingelheim Pharmaceuticals Inc ., Ridgefield , CT , USA
| |
Collapse
|
12
|
Abstract
The ability of polyvalent anions to influence protein-protein interactions and protein net charge was investigated through solubility and turbidity experiments, determination of osmotic second virial coefficients ( B22), and ζ-potential values for lysozyme solutions. B22 values showed that all anions reduce protein-protein repulsion between positively charged lysozyme molecules, and those anions with higher net valencies are more effective. The polyvalent anions pyrophosphate and tripolyphosphate were observed to induce protein reentrant condensation, which has been previously observed with negatively charged proteins in the presence of trivalent cations. Reentrant condensation is a phenomenon in which low concentrations of polyvalent ions induce protein precipitation, but further increasing polyvalent ion concentration causes the protein precipitate to resolubilize. Interestingly, citrate does not induce lysozyme reentrant condensation despite having a similar charge, size, and shape to pyrophosphate. We observe qualitative differences in protein behavior when compared against negatively charged proteins in solutions of trivalent cations. The polyphosphate ions induce a much stronger protein-protein attraction, which correlates with the occurrence of a liquid-gel transition that replaces the liquid-liquid transition observed with trivalent cations. The results indicate that solutions of polyphosphate ions provide a model system for exploring the link between the protein-phase diagram and model interaction potentials and also highlight the importance that ion-specific effects can have on protein solubility.
Collapse
Affiliation(s)
- Jordan W Bye
- School of Chemical Engineering and Analytical Science , The University of Manchester , Sackville Street , Manchester M13 9PL , U.K
| | - Robin A Curtis
- School of Chemical Engineering and Analytical Science , The University of Manchester , Sackville Street , Manchester M13 9PL , U.K
| |
Collapse
|
13
|
Ferreira GM, Calero-Rubio C, Sathish HA, Remmele RL, Roberts CJ. Electrostatically Mediated Protein-Protein Interactions for Monoclonal Antibodies: A Combined Experimental and Coarse-Grained Molecular Modeling Approach. J Pharm Sci 2019; 108:120-132. [DOI: 10.1016/j.xphs.2018.11.004] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2018] [Revised: 10/27/2018] [Accepted: 11/01/2018] [Indexed: 01/05/2023]
|
14
|
Hong T, Iwashita K, Shiraki K. Viscosity Control of Protein Solution by Small Solutes: A Review. Curr Protein Pept Sci 2018; 19:746-758. [PMID: 29237380 PMCID: PMC6182935 DOI: 10.2174/1389203719666171213114919] [Citation(s) in RCA: 69] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2017] [Revised: 12/01/2017] [Accepted: 12/02/2017] [Indexed: 12/22/2022]
Abstract
Viscosity of protein solution is one of the most troublesome issues for the high-concentration formulation of protein drugs. In this review, we summarize the practical methods that suppress the viscosity of protein solution using small molecular additives. The small amount of salts decreases the viscosity that results from electrostatic repulsion and attraction. The chaotrope suppresses the hydrophobic attraction and cluster formation, which can lower the solution viscosity. Arginine hydrochloride (ArgHCl) also suppresses the solution viscosity due to the hydrophobic and aromatic interactions between protein molecules. The small molecular additives are the simplest resolution of the high viscosity of protein solution as well as understanding of the primary cause in complex phenomena of protein interactions.
Collapse
Affiliation(s)
- Taehun Hong
- Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8573, Japan
| | - Kazuki Iwashita
- Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8573, Japan
| | - Kentaro Shiraki
- Faculty of Pure and Applied Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki305-8573, Japan
| |
Collapse
|
15
|
Safari MS, Byington MC, Conrad JC, Vekilov PG. Polymorphism of Lysozyme Condensates. J Phys Chem B 2017; 121:9091-9101. [DOI: 10.1021/acs.jpcb.7b05425] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Mohammad S. Safari
- Department
of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Road, Houston, Texas 77204-4004, United States
| | - Michael C. Byington
- Department
of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Road, Houston, Texas 77204-4004, United States
| | - Jacinta C. Conrad
- Department
of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Road, Houston, Texas 77204-4004, United States
| | - Peter G. Vekilov
- Department
of Chemical and Biomolecular Engineering, University of Houston, 4726 Calhoun Road, Houston, Texas 77204-4004, United States
- Department
of Chemistry, University of Houston, 3585 Cullen Blvd., Houston, Texas 77204-5003, United States
| |
Collapse
|
16
|
Laber JR, Dear BJ, Martins ML, Jackson DE, DiVenere A, Gollihar JD, Ellington AD, Truskett TM, Johnston KP, Maynard JA. Charge Shielding Prevents Aggregation of Supercharged GFP Variants at High Protein Concentration. Mol Pharm 2017; 14:3269-3280. [PMID: 28870080 DOI: 10.1021/acs.molpharmaceut.7b00322] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding protein stability is central to combatting protein aggregation diseases and developing new protein therapeutics. At the high concentrations often present in biological systems, purified proteins can exhibit undesirable high solution viscosities and poor solubilities mediated by short-range electrostatic and hydrophobic protein-protein interactions. The interplay between protein amino acid sequence, protein structure, and solvent conditions to minimize protein-protein interactions is key to designing well-behaved pharmaceutical proteins. However, theoretical approaches have yet to yield a general framework to address these problems. Here, we analyzed the high concentration behavior of superfolder GFP (sfGFP) and two supercharged sfGFP variants engineered to have formal charges of -18 or +15. Under low cosolute conditions, sfGFP and the -18 variant formed a gel or phase separated at ∼10 mg/mL. Under conditions that screen surface charges, including formulations with high histidine or high NaCl concentrations, all three variants attained concentrations up to 250 mg/mL with moderate viscosities. Moreover, all three variants exhibited very similar viscosity-concentration profiles over this range. This effect was not mimicked by high sugar concentrations that exert excluded-volume effects without shielding charge. Collectively, these data demonstrate that charge shielding neutralizes not only long-range electrostatic interactions but also, surprisingly, short-range electrostatic effects due to surface charge anisotropy. This work shows that supercharged sfGFP behavior under high ionic strength is largely determined by particle geometry, a conclusion that is supported by colloid models and may be applicable to pharmaceutically relevant proteins.
Collapse
Affiliation(s)
- Joshua R Laber
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Barton J Dear
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Matheus L Martins
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Devin E Jackson
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Andrea DiVenere
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Jimmy D Gollihar
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Andrew D Ellington
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Thomas M Truskett
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Keith P Johnston
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| | - Jennifer A Maynard
- Departments of †Chemical Engineering and ‡Molecular Biosciences, University of Texas at Austin , Austin, Texas 78712, United States
| |
Collapse
|
17
|
Corbett D, Hebditch M, Keeling R, Ke P, Ekizoglou S, Sarangapani P, Pathak J, Van Der Walle CF, Uddin S, Baldock C, Avendaño C, Curtis RA. Coarse-Grained Modeling of Antibodies from Small-Angle Scattering Profiles. J Phys Chem B 2017; 121:8276-8290. [DOI: 10.1021/acs.jpcb.7b04621] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Daniel Corbett
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester, M13 9PL, U.K
| | - Max Hebditch
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester, M13 9PL, U.K
| | - Rose Keeling
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester, M13 9PL, U.K
| | - Peng Ke
- Formulation
Sciences, MedImmune Ltd, Aaron Klug Building, Granta Park, Cambridge, CB21 6GH, U.K
| | - Sofia Ekizoglou
- Formulation
Sciences, MedImmune Ltd, Aaron Klug Building, Granta Park, Cambridge, CB21 6GH, U.K
| | - Prasad Sarangapani
- Regeneron Pharmaceuticals, 777
Old Saw Mill River Road, Tarrytown, New York 10591, United States
| | - Jai Pathak
- Vaccine
Research Center, National Institute of Health, 9 West Watkins Mill Road, Suite
250, Gaithersburg, Maryland 20878, United States
| | | | - Shahid Uddin
- Formulation
Sciences, MedImmune Ltd, Aaron Klug Building, Granta Park, Cambridge, CB21 6GH, U.K
| | - Clair Baldock
- Division
of Cell Matrix Biology and Regenerative Medicine, The University of Manchester, Oxford Road, Manchester, M13 9PT, U.K
| | - Carlos Avendaño
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester, M13 9PL, U.K
| | - Robin A. Curtis
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester, M13 9PL, U.K
| |
Collapse
|
18
|
The "Sticky Patch" Model of Crystallization and Modification of Proteins for Enhanced Crystallizability. Methods Mol Biol 2017; 1607:77-115. [PMID: 28573570 DOI: 10.1007/978-1-4939-7000-1_4] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Crystallization of macromolecules has long been perceived as a stochastic process, which cannot be predicted or controlled. This is consistent with another popular notion that the interactions of molecules within the crystal, i.e., crystal contacts, are essentially random and devoid of specific physicochemical features. In contrast, functionally relevant surfaces, such as oligomerization interfaces and specific protein-protein interaction sites, are under evolutionary pressures so their amino acid composition, structure, and topology are distinct. However, current theoretical and experimental studies are significantly changing our understanding of the nature of crystallization. The increasingly popular "sticky patch" model, derived from soft matter physics, describes crystallization as a process driven by interactions between select, specific surface patches, with properties thermodynamically favorable for cohesive interactions. Independent support for this model comes from various sources including structural studies and bioinformatics. Proteins that are recalcitrant to crystallization can be modified for enhanced crystallizability through chemical or mutational modification of their surface to effectively engineer "sticky patches" which would drive crystallization. Here, we discuss the current state of knowledge of the relationship between the microscopic properties of the target macromolecule and its crystallizability, focusing on the "sticky patch" model. We discuss state-of-the-art in silico methods that evaluate the propensity of a given target protein to form crystals based on these relationships, with the objective to design variants with modified molecular surface properties and enhanced crystallization propensity. We illustrate this discussion with specific cases where these approaches allowed to generate crystals suitable for structural analysis.
Collapse
|
19
|
Sarangapani PS, Hudson SD, Jones RL, Douglas JF, Pathak JA. Critical examination of the colloidal particle model of globular proteins. Biophys J 2015; 108:724-37. [PMID: 25650939 DOI: 10.1016/j.bpj.2014.11.3483] [Citation(s) in RCA: 62] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2014] [Revised: 11/01/2014] [Accepted: 11/05/2014] [Indexed: 10/24/2022] Open
Abstract
Recent studies of globular protein solutions have uniformly adopted a colloidal view of proteins as particles, a perspective that neglects the polymeric primary structure of these biological macromolecules, their intrinsic flexibility, and their ability to sample a large configurational space. While the colloidal perspective often serves as a useful idealization in many cases, the macromolecular identity of proteins must reveal itself under thermodynamic conditions in which the native state is no longer stable, such as denaturing solvents and high protein concentrations where macromolecules tend to have screened excluded volume, charge, and hydrodynamic interactions. Under extreme pH conditions, charge repulsion interactions within the protein chain can overcome the attractive hydrogen-bonding interactions, holding it in its native globular state. Conformational changes can therefore be expected to have great significance on the shear viscosity and other rheological properties of protein solutions. These changes are not envisioned in conventional colloidal protein models and we have initiated an investigation of the scattering and rheological properties of model proteins. We initiate this effort by considering bovine serum albumin because it is a globular protein whose solution properties have also been extensively investigated as a function of pH, temperature, ionic strength, and concentration. As we anticipated, near-ultraviolet circular dichroism measurements and intrinsic viscosity measurements clearly indicate that the bovine serum albumin tertiary structure changes as protein concentration and pH are varied. Our findings point to limited validity of the colloidal protein model and to the need for further consideration and quantification of the effects of conformational changes on protein solution viscosity, protein association, and the phase behavior. Small-angle Neutron Scattering measurements have allowed us to assess how these conformational changes influence protein size, shape, and interprotein interaction strength.
Collapse
Affiliation(s)
| | - Steven D Hudson
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland
| | - Ronald L Jones
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland
| | - Jack F Douglas
- Materials Science and Engineering Division, National Institute of Standards and Technology, Gaithersburg, Maryland
| | - Jai A Pathak
- Formulation Sciences Department, MedImmune, Gaithersburg, Maryland.
| |
Collapse
|
20
|
Arzenšek D, Kuzman D, Podgornik R. Hofmeister Effects in Monoclonal Antibody Solution Interactions. J Phys Chem B 2015. [DOI: 10.1021/acs.jpcb.5b02459] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Dejan Arzenšek
- Sandoz Biopharmaceuticals
Mengeš, Lek Pharmaceuticals d.d., Kolodvorska 27, Mengeš SI-1234, Slovenia
- Netica storitve
d.o.o., Reteče 97, Škofja Loka SI-4220, Slovenia
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana SI-1000, Slovenia
| | - Drago Kuzman
- Sandoz Biopharmaceuticals
Mengeš, Lek Pharmaceuticals d.d., Kolodvorska 27, Mengeš SI-1234, Slovenia
| | - Rudolf Podgornik
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, Ljubljana SI-1000, Slovenia
- Department
of Theoretical Physics, J. Stefan Institute, Jamova cesta 39, Ljubljana SI-1000, Slovenia
| |
Collapse
|
21
|
Rakel N, Bauer KC, Galm L, Hubbuch J. From osmotic second virial coefficient (B22 ) to phase behavior of a monoclonal antibody. Biotechnol Prog 2015; 31:438-51. [PMID: 25683855 DOI: 10.1002/btpr.2065] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2014] [Revised: 01/31/2015] [Indexed: 12/14/2022]
Abstract
Antibodies are complex macromolecules and their phase behavior as well as interactions within different solvents and precipitants are still not understood. To shed some light into the processes on a molecular dimension, the occurring self-interactions between antibody molecules were analyzed by means of the osmotic second virial coefficient (B22 ). The determined B22 follows qualitatively the phenomenological Hofmeister series describing the aggregation probability of antibodies for the various solvent compositions. However, a direct correlation between crystallization probability and B22 in form of a crystallization slot does not seem to be feasible for antibodies since the phase behavior is strongly dependent on their anisotropy. Kinetic parameters have to be taken into account due to the molecular size and complexity of the molecules. This is confirmed by a comparison of experimental data with a theoretical phase diagram. On the other hand the solubility is thermodynamically driven and therefore the B22 could be used to establish a universal solubility line for the monoclonal antibody mAb04c and different solvent compositions by using thermodynamic models.
Collapse
Affiliation(s)
- Natalie Rakel
- Section IV: Biomolecular Separation Engineering, Inst. of Engineering in Life Sciences, Karlsruhe Inst. of Technology, Engler-Bunte-Ring 1, Karlsruhe, 76131, Germany; Roche Diagnostics GmbH, Mannheim, Germany
| | | | | | | |
Collapse
|
22
|
Roberts D, Keeling R, Tracka M, van der Walle CF, Uddin S, Warwicker J, Curtis R. Specific Ion and Buffer Effects on Protein–Protein Interactions of a Monoclonal Antibody. Mol Pharm 2014; 12:179-93. [DOI: 10.1021/mp500533c] [Citation(s) in RCA: 109] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- D. Roberts
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, U.K
| | - R. Keeling
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, U.K
| | - M. Tracka
- Formulation
Sciences, MedImmune, Ltd., Aaron Klug Building, Granta Park, Cambridge CB21 6GH, U.K
| | - C. F. van der Walle
- Formulation
Sciences, MedImmune, Ltd., Aaron Klug Building, Granta Park, Cambridge CB21 6GH, U.K
| | - S. Uddin
- Formulation
Sciences, MedImmune, Ltd., Aaron Klug Building, Granta Park, Cambridge CB21 6GH, U.K
| | - J. Warwicker
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, U.K
| | - R. Curtis
- School
of Chemical Engineering and Analytical Science, The University of Manchester, Sackville Street, Manchester M13 9PL, U.K
| |
Collapse
|
23
|
Lapelosa M, Patapoff TW, Zarraga IE. Molecular Simulations of the Pairwise Interaction of Monoclonal Antibodies. J Phys Chem B 2014; 118:13132-41. [DOI: 10.1021/jp508729z] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Mauro Lapelosa
- Department of Late Stage Pharmaceutical Development and ‡Department of Early Stage Pharmaceutical
Development, Genentech Inc., member of Roche, South San Francisco, California 94080, United States
| | - Thomas W. Patapoff
- Department of Late Stage Pharmaceutical Development and ‡Department of Early Stage Pharmaceutical
Development, Genentech Inc., member of Roche, South San Francisco, California 94080, United States
| | - Isidro E. Zarraga
- Department of Late Stage Pharmaceutical Development and ‡Department of Early Stage Pharmaceutical
Development, Genentech Inc., member of Roche, South San Francisco, California 94080, United States
| |
Collapse
|
24
|
Uchiyama S. Liquid formulation for antibody drugs. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2014; 1844:2041-2052. [DOI: 10.1016/j.bbapap.2014.07.016] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/10/2014] [Revised: 07/16/2014] [Accepted: 07/25/2014] [Indexed: 01/21/2023]
|
25
|
Andrews CT, Elcock AH. COFFDROP: A Coarse-Grained Nonbonded Force Field for Proteins Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of Amino Acids. J Chem Theory Comput 2014; 10:5178-5194. [PMID: 25400526 PMCID: PMC4230375 DOI: 10.1021/ct5006328] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2014] [Indexed: 02/06/2023]
Abstract
![]()
We describe the derivation of a set
of bonded and nonbonded coarse-grained
(CG) potential functions for use in implicit-solvent Brownian dynamics
(BD) simulations of proteins derived from all-atom explicit-solvent
molecular dynamics (MD) simulations of amino acids. Bonded potential
functions were derived from 1 μs MD simulations of each of the
20 canonical amino acids, with histidine modeled in both its protonated
and neutral forms; nonbonded potential functions were derived from
1 μs MD simulations of every possible pairing of the amino acids
(231 different systems). The angle and dihedral probability distributions
and radial distribution functions sampled during MD were used to optimize
a set of CG potential functions through use of the iterative Boltzmann
inversion (IBI) method. The optimized set of potential functions—which
we term COFFDROP (COarse-grained Force Field for Dynamic Representation
Of Proteins)—quantitatively reproduced all of the “target”
MD distributions. In a first test of the force field, it was used
to predict the clustering behavior of concentrated amino acid solutions;
the predictions were directly compared with the results of corresponding
all-atom explicit-solvent MD simulations and found to be in excellent
agreement. In a second test, BD simulations of the small protein villin
headpiece were carried out at concentrations that have recently been
studied in all-atom explicit-solvent MD simulations by Petrov and
Zagrovic (PLoS Comput. Biol.2014, 5, e1003638). The anomalously strong intermolecular interactions
seen in the MD study were reproduced in the COFFDROP simulations;
a simple scaling of COFFDROP’s nonbonded parameters, however,
produced results in better accordance with experiment. Overall, our
results suggest that potential functions derived from simulations
of pairwise amino acid interactions might be of quite broad applicability,
with COFFDROP likely to be especially useful for modeling unfolded
or intrinsically disordered proteins.
Collapse
Affiliation(s)
- Casey T Andrews
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| | - Adrian H Elcock
- Department of Biochemistry, University of Iowa , Iowa City, Iowa 52242, United States
| |
Collapse
|
26
|
Kaieda S, Lund M, Plivelic TS, Halle B. Weak self-interactions of globular proteins studied by small-angle X-ray scattering and structure-based modeling. J Phys Chem B 2014; 118:10111-9. [PMID: 25117055 DOI: 10.1021/jp505809v] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We investigate protein-protein interactions in solution by small-angle X-ray scattering (SAXS) and theoretical modeling. The structure factor for solutions of bovine pancreatic trypsin inhibitor (BPTI), myoglobin (Mb), and intestinal fatty acid-binding protein (IFABP) is determined from SAXS measurements at multiple concentrations, from Monte Carlo simulations with a coarse-grained structure-based interaction model, and from analytic approximate solutions of two idealized colloidal interaction models without adjustable parameters. By combining these approaches, we find that the structure factor is essentially determined by hard-core and screened electrostatic interactions. Other soft short-ranged interactions (van der Waals and solvation-related) are either individually insignificant or tend to cancel out. The structure factor is also not significantly affected by charge fluctuations. For Mb and IFABP, with a small net charge and relatively symmetric charge distribution, the structure factor is well described by a hard-sphere model. For BPTI, with a larger net charge, screened electrostatic repulsion is also important, but the asymmetry of the charge distribution reduces the repulsion from that predicted by a charged hard-sphere model with the same net charge. Such charge asymmetry may also amplify the effect of shape asymmetry on the protein-protein potential of mean force.
Collapse
Affiliation(s)
- Shuji Kaieda
- Department of Biophysical Chemistry, Lund University , P.O. Box 124, SE-22100 Lund, Sweden
| | | | | | | |
Collapse
|
27
|
Roberts D, Keeling R, Tracka M, van der Walle CF, Uddin S, Warwicker J, Curtis R. The role of electrostatics in protein-protein interactions of a monoclonal antibody. Mol Pharm 2014; 11:2475-89. [PMID: 24892385 DOI: 10.1021/mp5002334] [Citation(s) in RCA: 113] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Understanding how protein-protein interactions depend on the choice of buffer, salt, ionic strength, and pH is needed to have better control over protein solution behavior. Here, we have characterized the pH and ionic strength dependence of protein-protein interactions in terms of an interaction parameter kD obtained from dynamic light scattering and the osmotic second virial coefficient B22 measured by static light scattering. A simplified protein-protein interaction model based on a Baxter adhesive potential and an electric double layer force is used to separate out the contributions of longer-ranged electrostatic interactions from short-ranged attractive forces. The ionic strength dependence of protein-protein interactions for solutions at pH 6.5 and below can be accurately captured using a Deryaguin-Landau-Verwey-Overbeek (DLVO) potential to describe the double layer forces. In solutions at pH 9, attractive electrostatics occur over the ionic strength range of 5-275 mM. At intermediate pH values (7.25 to 8.5), there is a crossover effect characterized by a nonmonotonic ionic strength dependence of protein-protein interactions, which can be rationalized by the competing effects of long-ranged repulsive double layer forces at low ionic strength and a shorter ranged electrostatic attraction, which dominates above a critical ionic strength. The change of interactions from repulsive to attractive indicates a concomitant change in the angular dependence of protein-protein interaction from isotropic to anisotropic. In the second part of the paper, we show how the Baxter adhesive potential can be used to predict values of kD from fitting to B22 measurements, thus providing a molecular basis for the linear correlation between the two protein-protein interaction parameters.
Collapse
Affiliation(s)
- D Roberts
- School of Chemical Engineering and Analytical Science, The University of Manchester , Sackville Street, Manchester M13 9PL, U.K
| | | | | | | | | | | | | |
Collapse
|
28
|
Mereghetti P, Martinez M, Wade RC. Long range Debye-Hückel correction for computation of grid-based electrostatic forces between biomacromolecules. BMC BIOPHYSICS 2014; 7:4. [PMID: 25045516 PMCID: PMC4082500 DOI: 10.1186/2046-1682-7-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/11/2014] [Accepted: 06/04/2014] [Indexed: 12/01/2022]
Abstract
Background Brownian dynamics (BD) simulations can be used to study very large molecular systems, such as models of the intracellular environment, using atomic-detail structures. Such simulations require strategies to contain the computational costs, especially for the computation of interaction forces and energies. A common approach is to compute interaction forces between macromolecules by precomputing their interaction potentials on three-dimensional discretized grids. For long-range interactions, such as electrostatics, grid-based methods are subject to finite size errors. We describe here the implementation of a Debye-Hückel correction to the grid-based electrostatic potential used in the SDA BD simulation software that was applied to simulate solutions of bovine serum albumin and of hen egg white lysozyme. Results We found that the inclusion of the long-range electrostatic correction increased the accuracy of both the protein-protein interaction profiles and the protein diffusion coefficients at low ionic strength. Conclusions An advantage of this method is the low additional computational cost required to treat long-range electrostatic interactions in large biomacromolecular systems. Moreover, the implementation described here for BD simulations of protein solutions can also be applied in implicit solvent molecular dynamics simulations that make use of gridded interaction potentials.
Collapse
Affiliation(s)
- Paolo Mereghetti
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloß-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany ; Center for Nanotechnology Innovation@NEST, Italian Institute of Technology, Piazza San Silvestro 12, Pisa, Italy
| | - Michael Martinez
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloß-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Rebecca C Wade
- Molecular and Cellular Modeling Group, Heidelberg Institute for Theoretical Studies (HITS), Schloß-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany ; Center for Molecular Biology (ZMBH), University of Heidelberg, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany
| |
Collapse
|
29
|
Quang LJ, Sandler SI, Lenhoff AM. Anisotropic Contributions to Protein–Protein Interactions. J Chem Theory Comput 2014; 10:835-45. [DOI: 10.1021/ct4006695] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Leigh J. Quang
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| | - Stanley I. Sandler
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| | - Abraham M. Lenhoff
- Department of Chemical and
Biomolecular Engineering, University of Delaware, Newark, Delaware 19716, United States of America
| |
Collapse
|
30
|
Hoppe T. A simplified representation of anisotropic charge distributions within proteins. J Chem Phys 2013; 138:174110. [PMID: 23656117 DOI: 10.1063/1.4803099] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Effective coarse-grained representations of protein-protein interaction potentials are vital in the modeling of large scale systems. We develop a method to fit an arbitrary number of effective charges to approximate the electrostatic potential of a protein at a given pH in an ionic solution. We find that the effective charges can reproduce an input potential calculated from a high resolution Poisson-Boltzmann calculation. Since the effective charges used in this model are not constrained to the locations of the original charged groups, the extra degrees of freedom allows us to reproduce the field anisotropy with fewer charges. The fitting procedure uses a number of approximations in the charge magnitudes, initial conditions, and multipoles to speed convergence. The most significant gains are found by fitting the multipole moments of the effective charge potential to the moments of the original field. We show that the Yukawa potential is not only sufficient as a pairwise summation in reproducing the potential, but comes naturally from the linearized expansion of the Poisson-Boltzmann equation. We compute interaction energies and find excellent agreement to the original potential. From the effective charge model we compute the electrostatic contribution to the second virial coefficient.
Collapse
Affiliation(s)
- Travis Hoppe
- National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, Maryland 20892, USA.
| |
Collapse
|
31
|
Stark AC, Andrews CT, Elcock AH. Toward optimized potential functions for protein-protein interactions in aqueous solutions: osmotic second virial coefficient calculations using the MARTINI coarse-grained force field. J Chem Theory Comput 2013; 9. [PMID: 24223529 DOI: 10.1021/ct400008p] [Citation(s) in RCA: 98] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Coarse-grained (CG) simulation methods are now widely used to model the structure and dynamics of large biomolecular systems. One important issue for using such methods - especially with regard to using them to model, for example, intracellular environments - is to demonstrate that they can reproduce experimental data on the thermodynamics of protein-protein interactions in aqueous solutions. To examine this issue, we describe here simulations performed using the popular coarse-grained MARTINI force field, aimed at computing the thermodynamics of lysozyme and chymotrypsinogen self-interactions in aqueous solution. Using molecular dynamics simulations to compute potentials of mean force between a pair of protein molecules, we show that the original parameterization of the MARTINI force field is likely to significantly overestimate the strength of protein-protein interactions to the extent that the computed osmotic second virial coefficients are orders of magnitude more negative than experimental estimates. We then show that a simple down-scaling of the van der Waals parameters that describe the interactions between protein pseudo-atoms can bring the simulated thermodynamics into much closer agreement with experiment. Overall, the work shows that it is feasible to test explicit-solvent CG force fields directly against thermodynamic data for proteins in aqueous solutions, and highlights the potential usefulness of osmotic second virial coefficient measurements for fully parameterizing such force fields.
Collapse
Affiliation(s)
- Austin C Stark
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242
| | | | | |
Collapse
|
32
|
Abstract
AbstractProteins and other biomolecules contain acidic and basic titratable groups that give rise to intricate charge distributions and control electrostatic interactions. ‘Charge regulation’ concerns how the proton equilibria of these sites are perturbed when approached by alien molecular matter such as other proteins, surfaces and membranes, DNA, polyelectrolytes etc. Importantly, this perturbation generates a charge response that leads to attractive intermolecular interactions that can be conveniently described by a single molecular property – the charge capacitance. The capacitance quantifies molecular charge fluctuations, i.e. it is the variance of the mean charge and is an intrinsic property on par with the net charge and the dipole moment. It directly enters the free energy expression for intermolecular interactions and can be obtained experimentally from the derivative of the titration curve or theoretically from simulations. In this review, we focus on the capacitance concept as a predictive parameter for charge regulation and demonstrate how it can be used to estimate the interaction of a protein with other proteins, polyelectrolytes, membranes as well as with ligands.
Collapse
|
33
|
Zhu M, Simone AD, Schenk D, Toth G, Dobson CM, Vendruscolo M. Identification of small-molecule binding pockets in the soluble monomeric form of the Aβ42 peptide. J Chem Phys 2013; 139:035101. [PMID: 23883055 PMCID: PMC5011423 DOI: 10.1063/1.4811831] [Citation(s) in RCA: 65] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The aggregation of intrinsically disordered peptides and proteins is associated with a wide range of highly debilitating neurological and systemic disorders. In this work we explored the potential of a structure-based drug discovery procedure to target one such system, the soluble monomeric form of the Aβ42 peptide. We utilised for this purpose a set of structures of the Aβ42 peptide selected from clusters of conformations within an ensemble generated by molecular dynamics simulations. Using these structures we carried out fragment mapping calculations to identify binding "hot spots" on the monomeric form of the Aβ42 peptide. This procedure provided a set of hot spots with ligand efficiencies comparable to those observed for structured proteins, and clustered into binding pockets. Such binding pockets exhibited a propensity to bind small molecules known to interact with the Aβ42 peptide. Taken together these results provide an initial indication that fragment-based drug discovery may represent a potential therapeutic strategy for diseases associated with the aggregation of intrinsically disordered proteins.
Collapse
Affiliation(s)
- Maximillian Zhu
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, UK
| | - Alfonso De Simone
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, UK
| | - Dale Schenk
- Elan Pharmaceuticals, South San Francisco, CA 94080, USA
| | - Gergely Toth
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, UK
| | | | | |
Collapse
|
34
|
Chinchalikar AJ, Aswal VK, Kohlbrecher J, Wagh AG. Small-angle neutron scattering study of structure and interaction during salt-induced liquid-liquid phase transition in protein solutions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2013; 87:062708. [PMID: 23848716 DOI: 10.1103/physreve.87.062708] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/23/2013] [Indexed: 06/02/2023]
Abstract
The liquid-liquid phase transition (LLPT) in aqueous salt solutions of lysozyme protein has been studied by small-angle neutron scattering. Measurements have been carried out on fixed protein concentration with varying salt concentration approaching LLPT. The data are fitted considering protein interaction by the two Yukawa (2Y) potential which combines short-range attraction and long-range repulsion. We show that LLPT arises because of enhancement of non-DLVO (Derjaguin-Landau-Verwey-Overbeek) short-range attraction without any conformational structural change of the protein. The salt concentration required for LLPT as well as corresponding short-range attraction decreases significantly with increase in protein concentration.
Collapse
Affiliation(s)
- A J Chinchalikar
- Solid State Physics Division, Bhabha Atomic Research Centre, Mumbai-400085, India
| | | | | | | |
Collapse
|
35
|
Brunsteiner M, Flock M, Nidetzky B. Structure based descriptors for the estimation of colloidal interactions and protein aggregation propensities. PLoS One 2013; 8:e59797. [PMID: 23565169 PMCID: PMC3614552 DOI: 10.1371/journal.pone.0059797] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2012] [Accepted: 02/19/2013] [Indexed: 11/19/2022] Open
Abstract
The control of protein aggregation is an important requirement in the development of bio-pharmaceutical formulations. Here a simple protein model is proposed that was used in molecular dynamics simulations to obtain a quantitative assessment of the relative contributions of proteins' net-charges, dipole-moments, and the size of hydrophobic or charged surface patches to their colloidal interactions. The results demonstrate that the strength of these interactions correlate with net-charge and dipole moment. Variation of both these descriptors within ranges typical for globular proteins have a comparable effect. By comparison no clear trends can be observed upon varying the size of hydrophobic or charged patches while keeping the other parameters constant. The results are discussed in the context of experimental literature data on protein aggregation. They provide a clear guide line for the development of improved algorithms for the prediction of aggregation propensities.
Collapse
Affiliation(s)
- Michael Brunsteiner
- Research Center Pharmaceutical Engineering, Graz, Austria
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Graz, Austria
| | - Michaela Flock
- Institute of Inorganic Chemistry, Graz University of Technology, Graz, Austria
| | - Bernd Nidetzky
- Institute of Biotechnology and Biochemical Engineering, Graz University of Technology, Graz, Austria
| |
Collapse
|
36
|
Scherer TM. Cosolute Effects on the Chemical Potential and Interactions of an IgG1 Monoclonal Antibody at High Concentrations. J Phys Chem B 2013; 117:2254-66. [DOI: 10.1021/jp3091717] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Thomas M. Scherer
- Genentech (a Member of the
Roche Group), Late Stage Pharmaceutical Development, 1 DNA Way, South
San Francisco, California 94080, United States
| |
Collapse
|
37
|
Second Osmotic Virial Coefficients and Aggregation of Monoclonal Antibodies by Static Laser Light Scattering. Z PHYS CHEM 2013. [DOI: 10.1524/zpch.2013.0325] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Abstract
The second osmotic virial coefficient and the apparent molar mass of two human and one mouse monoclonal antibodies were measured in different aequeous buffer solutions which also contained sodium chloride or ammonium sulfate, respectively, by static laser light scattering in batch mode. The apparent molar mass indicates aggregation. At a constant pH value of 6.5 the sodium chloride concentration was varied from 0 to 2 M and the ammonium sulfate concentration from 0 to 0.8 M, respectively. A 20 mM sodium-phosphate buffer was used for all experiments. Furthermore the pH value was varied without adding additional salt from 4.5 to 10. The results of the salt dependency are in line with the Hofmeister-series. The results of the pH dependency correspond to the net charge of the molecules.
Collapse
|
38
|
Saito S, Hasegawa J, Kobayashi N, Tomitsuka T, Uchiyama S, Fukui K. Effects of ionic strength and sugars on the aggregation propensity of monoclonal antibodies: influence of colloidal and conformational stabilities. Pharm Res 2013; 30:1263-80. [PMID: 23319172 DOI: 10.1007/s11095-012-0965-4] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2012] [Accepted: 12/13/2012] [Indexed: 10/27/2022]
Abstract
PURPOSE To develop a general strategy for optimizing monoclonal antibody (MAb) formulations. METHODS Colloidal stabilities of four representative MAbs solutions were assessed based on the second virial coefficient (B 2) at 20°C and 40°C, and net charges at different NaCl concentrations, and/or in the presence of sugars. Conformational stabilities were evaluated from the unfolding temperatures. The aggregation propensities were determined at 40°C and after freeze-thawing. The electrostatic potential of antibody surfaces was simulated for the development of rational formulations. RESULTS Similar B 2 values were obtained at 20°C and 40°C, implying little dependence on temperature. B 2 correlated quantitatively with aggregation propensities at 40°C. The net charge partly correlated with colloidal stability. Salts stabilized or destabilized MAbs, depending on repulsive or attractive interactions. Sugars improved the aggregation propensity under freeze-thaw stress through improved conformational stability. Uneven and even distributions of potential surfaces were attributed to attractive and strong repulsive electrostatic interactions. CONCLUSIONS Assessment of colloidal stability at the lowest ionic strength is particularly effective for the development of formulations. If necessary, salts are added to enhance the colloidal stability. Sugars further improved aggregation propensities by enhancing conformational stability. These behaviors are rationally predictable according to the surface potentials of MAbs.
Collapse
Affiliation(s)
- Shuntaro Saito
- Analytical & Quality Evaluation Research Laboratories, Daiichi Sankyo Co., Ltd., 1-12-1, Shinomiya, Hiratsuka-shi, Kanagawa, 254-0014, Japan
| | | | | | | | | | | |
Collapse
|
39
|
Chan HY, Lankevich V, Vekilov PG, Lubchenko V. Anisotropy of the Coulomb interaction between folded proteins: consequences for mesoscopic aggregation of lysozyme. Biophys J 2012; 102:1934-43. [PMID: 22768950 DOI: 10.1016/j.bpj.2012.03.025] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2011] [Revised: 02/16/2012] [Accepted: 03/02/2012] [Indexed: 10/28/2022] Open
Abstract
Toward quantitative description of protein aggregation, we develop a computationally efficient method to evaluate the potential of mean force between two folded protein molecules that allows for complete sampling of their mutual orientation. Our model is valid at moderate ionic strengths and accounts for the actual charge distribution on the surface of the molecules, the dielectric discontinuity at the protein-solvent interface, and the possibility of protonation or deprotonation of surface residues induced by the electric field due to the other protein molecule. We apply the model to the protein lysozyme, whose solutions exhibit both mesoscopic clusters of protein-rich liquid and liquid-liquid separation; the former requires that protein form complexes with typical lifetimes of approximately milliseconds. We find the electrostatic repulsion is typically lower than the prediction of the Derjaguin-Landau-Verwey-Overbeek theory. The Coulomb interaction in the lowest-energy docking configuration is nonrepulsive, despite the high positive charge on the molecules. Typical docking configurations barely involve protonation or deprotonation of surface residues. The obtained potential of mean force between folded lysozyme molecules is consistent with the location of the liquid-liquid coexistence, but produces dimers that are too short-lived for clusters to exist, suggesting lysozyme undergoes conformational changes during cluster formation.
Collapse
Affiliation(s)
- Ho Yin Chan
- Department of Physics, University of Houston, Houston, Texas, USA
| | | | | | | |
Collapse
|
40
|
Chari R, Singh SN, Yadav S, Brems DN, Kalonia DS. Determination of the dipole moments of RNAse SA wild type and a basic mutant. Proteins 2011; 80:1041-52. [PMID: 22213585 DOI: 10.1002/prot.24006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2011] [Revised: 11/18/2011] [Accepted: 11/22/2011] [Indexed: 12/22/2022]
Abstract
In this study, we report the effects of acidic to basic residue point mutations (5K) on the dipole moment of RNAse SA at different pHs. Dipole moments were determined by measuring solution capacitance of the wild type (WT) and the 5K mutant with an impedance analyzer. The dipole moments were then (1) compared with theoretically calculated dipole moments, (2) analyzed to determine the effect of the point mutations, and (3) analyzed for their contribution to overall protein-protein interactions (PPI) in solution as quantitated by experimentally derived second virial coefficients. We determined that experimental and calculated dipoles were in reasonable agreement. Differences are likely due to local motions of residue side chains, which are not accounted for by the calculated dipole. We observed that the proteins' dipole moments increase as the pH is shifted further from their isoelectric points and that the wild-type dipole moments were greater than those of the 5K. This is likely due to an increase in the proportion of one charge (either negative or positive) relative to the other. A greater charge disparity corresponded to a larger dipole moment. Finally, the larger dipole moments of the WT resulted in greater attractive overall PPI for that protein as compared to the 5K.
Collapse
Affiliation(s)
- Ravi Chari
- Pharmaceutics, GPSC, Abbott Bioresearch Center, Worcester, Massachusetts 01605, USA
| | | | | | | | | |
Collapse
|
41
|
Kumar V, Dixit N, Zhou L(L, Fraunhofer W. Impact of short range hydrophobic interactions and long range electrostatic forces on the aggregation kinetics of a monoclonal antibody and a dual-variable domain immunoglobulin at low and high concentrations. Int J Pharm 2011; 421:82-93. [DOI: 10.1016/j.ijpharm.2011.09.017] [Citation(s) in RCA: 61] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2011] [Revised: 08/17/2011] [Accepted: 09/19/2011] [Indexed: 10/17/2022]
|
42
|
Viscosity of concentrated therapeutic protein compositions. Adv Drug Deliv Rev 2011; 63:1107-17. [PMID: 22014592 DOI: 10.1016/j.addr.2011.09.008] [Citation(s) in RCA: 99] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2011] [Accepted: 09/02/2011] [Indexed: 11/22/2022]
Abstract
The use of monoclonal antibodies as therapeutic agents has been increasing steadily over the last decade for the treatment of various conditions. There is often a need to deliver a large dose of the protein, so there is a trend toward developing commercially viable liquid formulations of highly concentrated antibodies. Such concentrated solutions are associated with a number of challenges, including optimization of production processes, plus chemical and physical stability of the final product where solution viscosity becomes a critical quality attribute. Assessment of the rheological characteristics of concentrated compositions is essential as are development strategies to reduce the viscosity. This review covers the state-of-the-art rheology measurement techniques, focusing particularly on concentrated protein solutions. Current understanding of the mechanisms leading to high viscosity and control by formulation parameters is discussed.
Collapse
|
43
|
Mereghetti P, Gabdoulline RR, Wade RC. Brownian dynamics simulation of protein solutions: structural and dynamical properties. Biophys J 2011; 99:3782-91. [PMID: 21112303 DOI: 10.1016/j.bpj.2010.10.035] [Citation(s) in RCA: 76] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2010] [Revised: 10/19/2010] [Accepted: 10/20/2010] [Indexed: 11/17/2022] Open
Abstract
The study of solutions of biomacromolecules provides an important basis for understanding the behavior of many fundamental cellular processes, such as protein folding, self-assembly, biochemical reactions, and signal transduction. Here, we describe a Brownian dynamics simulation procedure and its validation for the study of the dynamic and structural properties of protein solutions. In the model used, the proteins are treated as atomically detailed rigid bodies moving in a continuum solvent. The protein-protein interaction forces are described by the sum of electrostatic interaction, electrostatic desolvation, nonpolar desolvation, and soft-core repulsion terms. The linearized Poisson-Boltzmann equation is solved to compute electrostatic terms. Simulations of homogeneous solutions of three different proteins with varying concentrations, pH, and ionic strength were performed. The results were compared to experimental data and theoretical values in terms of long-time self-diffusion coefficients, second virial coefficients, and structure factors. The results agree with the experimental trends and, in many cases, experimental values are reproduced quantitatively. There are no parameters specific to certain protein types in the interaction model, and hence the model should be applicable to the simulation of the behavior of mixtures of macromolecules in cell-like crowded environments.
Collapse
Affiliation(s)
- Paolo Mereghetti
- Heidelberg Institute for Theoretical Studies, Heidelberg, Germany.
| | | | | |
Collapse
|
44
|
Długosz M, Trylska J. Diffusion in crowded biological environments: applications of Brownian dynamics. BMC BIOPHYSICS 2011; 4:3. [PMID: 21595998 PMCID: PMC3093676 DOI: 10.1186/2046-1682-4-3] [Citation(s) in RCA: 48] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 02/09/2011] [Accepted: 03/02/2011] [Indexed: 01/10/2023]
Abstract
Biochemical reactions in living systems occur in complex, heterogeneous media with total concentrations of macromolecules in the range of 50 - 400 mgml. Molecular species occupy a significant fraction of the immersing medium, up to 40% of volume. Such complex and volume-occupied environments are generally termed 'crowded' and/or 'confined'. In crowded conditions non-specific interactions between macromolecules may hinder diffusion - a major process determining metabolism, transport, and signaling. Also, the crowded media can alter, both qualitatively and quantitatively, the reactions in vivo in comparison with their in vitro counterparts. This review focuses on recent developments in particle-based Brownian dynamics algorithms, their applications to model diffusive transport in crowded systems, and their abilities to reproduce and predict the behavior of macromolecules under in vivo conditions.
Collapse
Affiliation(s)
- Maciej Długosz
- Interdisciplinary Centre for Mathematical and Computational Modeling, University of Warsaw, Żwirki i Wigury 93, 02-089 Warsaw, Poland
| | - Joanna Trylska
- Interdisciplinary Centre for Mathematical and Computational Modeling, University of Warsaw, Żwirki i Wigury 93, 02-089 Warsaw, Poland
| |
Collapse
|
45
|
Abstract
Due to large equilibrium fluctuations of protons at physiological pH, the orthophosphate ion as well as the imidazole group on histidine substantially regulate their charge upon approaching charged interfaces. This implies that these-and comparable-ions function as electostatic "proximity switches" when interacting with lipid membranes, DNA, proteins, etc. Using straightforward statistical thermodynamics as well as mesoscopic computer simulations we quantify the charge regulation mechanism and argue that it is important in a range of biological as well as technical processes.
Collapse
Affiliation(s)
- Mikael Lund
- Department of Theoretical Chemistry, Lund University, POB 124, SE-22100 Lund, Sweden.
| |
Collapse
|
46
|
Abstract
We describe a model for protein crystallization equilibria. The model includes four terms, (1) protein translational entropy opposes crystallization, (2) proteins are attracted to each other by a nonelectrostatic contact free energy favoring crystallization, (3) proteins in the crystal repel each other but, to a greater extent, attract counterions sequestered in the crystal, which favors crystallization, and (4) the translational entropy of the counterions opposes their sequestration into the crystal, opposing crystallization. We treat the electrostatics using the nonlinear Poisson-Boltzmann equation, and we use unit cell information from native protein crystals to determine the boundary conditions. This model predicts the stabilities of protein crystals as functions of temperature, pH, and salt concentrations, in good agreement with the data of Pusey et al. on tetragonal and orthorhombic crystal forms of lysozyme. The experiments show a weak dependence of crystal solubility on pH. According to the model, this is because the entropic cost to neutralize the crystal is compensated by favorable protein-salt interactions. Experiments also show that adding salt stabilizes the crystal. Cohn's empirical law predicts that the logarithm of solubility should be a linear function of salt. The present theory predicts nonlinearity, in better agreement with the experiments. The model shows that the salting out phenomena is not due to more counterion shielding but to lowered counterion translational entropy. Models of this type may help guide faster and better ways to crystallize proteins.
Collapse
Affiliation(s)
- Jeremy D. Schmit
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA 94158
| | - Ken A. Dill
- Department of Pharmaceutical Chemistry, University of California at San Francisco, San Francisco, CA 94158
| |
Collapse
|
47
|
McGuffee SR, Elcock AH. Diffusion, crowding & protein stability in a dynamic molecular model of the bacterial cytoplasm. PLoS Comput Biol 2010; 6:e1000694. [PMID: 20221255 PMCID: PMC2832674 DOI: 10.1371/journal.pcbi.1000694] [Citation(s) in RCA: 524] [Impact Index Per Article: 37.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2009] [Accepted: 01/30/2010] [Indexed: 01/24/2023] Open
Abstract
A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo. The interior of a typical bacterial cell is a highly crowded place in which molecules must jostle and compete with each other in order to carry out their biological functions. The conditions under which such molecules are typically studied in vitro, however, are usually quite different: one or a few different types of molecules are studied as they freely diffuse in a dilute, aqueous solution. There is therefore a significant disconnect between the conditions under which molecules can be most usefully studied and the conditions under which such molecules usually “live”, and developing ways to bridge this gap is likely to be important for properly understanding molecular behavior in vivo. Toward this end, we show in this work that computer simulations can be used to model the interior of bacterial cells at a near atomic level of detail: the rates of diffusion of proteins are matched to known experimental values, and their thermodynamic stabilities are found to be in good agreement with the few measurements that have so far been performed in vivo. While the simulation approach is certainly not free of assumptions, it offers a potentially important complement to experimental techniques and provides a vivid illustration of molecular behavior inside a biological cell that is likely to be of significant educational value.
Collapse
Affiliation(s)
- Sean R. McGuffee
- Department of Biochemistry, University of Iowa, Iowa City, Iowa, United States of America
| | - Adrian H. Elcock
- Department of Biochemistry, University of Iowa, Iowa City, Iowa, United States of America
- * E-mail:
| |
Collapse
|
48
|
Chari R, Jerath K, Badkar AV, Kalonia DS. Long- and short-range electrostatic interactions affect the rheology of highly concentrated antibody solutions. Pharm Res 2009; 26:2607-18. [PMID: 19795191 DOI: 10.1007/s11095-009-9975-2] [Citation(s) in RCA: 113] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2009] [Accepted: 09/14/2009] [Indexed: 01/10/2023]
Abstract
PURPOSE To explain the differences in protein-protein interactions (PPI) of concentrated versus dilute formulations of a model antibody. METHODS High frequency rheological measurements from pH 3.0 to 12.0 quantitated viscoelasticity and PPI at high concentrations. Dynamic light scattering (DLS) characterized PPI in dilute solutions. RESULTS For concentrated solutions at low ionic strength, the storage modulus, a viscosity component and a measure of PPI, is highest at the isoelectric point (pH 9.0) and lowest at pH 5.4. This profile flattens at higher ionic strength but not completely, indicating PPI consist of long-range electrostatics and other short-range attractions. At low concentrations, PPI are near zero at pI but become repulsive as the pH is shifted. Higher salt concentrations completely flatten this profile to zero, indicating that these PPI are mainly electrostatic. CONCLUSIONS This discrepancy occurs because long-range interactions are significant at low concentrations, whereas both long- and short-range interactions are significant at higher concentrations. Computer modeling was used to calculate antibody properties responsible for long- and short-range interactions, i.e. net charge and dipole moment. Charge-charge interactions are repulsive while dipole-dipole interactions are attractive. Their net effect correlated with the storage modulus profile. However, only charge-charge repulsions correlated with PPI determined by DLS.
Collapse
Affiliation(s)
- Ravi Chari
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Connecticut, 69 North Eagleville Rd., U-3092, Storrs, Connecticut 06269, USA
| | | | | | | |
Collapse
|
49
|
Gabdoulline RR, Wade RC. On the contributions of diffusion and thermal activation to electron transfer between Phormidium laminosum plastocyanin and cytochrome f: Brownian dynamics simulations with explicit modeling of nonpolar desolvation interactions and electron transfer events. J Am Chem Soc 2009; 131:9230-8. [PMID: 19518050 DOI: 10.1021/ja809567k] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The factors that determine the extent to which diffusion and thermal activation processes govern electron transfer (ET) between proteins are debated. The process of ET between plastocyanin (PC) and cytochrome f (CytF) from the cyanobacterium Phormidium laminosum was initially thought to be diffusion-controlled but later was found to be under activation control (Schlarb-Ridley, B. G.; et al. Biochemistry 2005, 44, 6232). Here we describe Brownian dynamics simulations of the diffusional association of PC and CytF, from which ET rates were computed using a detailed model of ET events that was applied to all of the generated protein configurations. The proteins were modeled as rigid bodies represented in atomic detail. In addition to electrostatic forces, which were modeled as in our previous simulations of protein-protein association, the proteins interacted by a nonpolar desolvation (hydrophobic) force whose derivation is described here. The simulations yielded close to realistic residence times of transient protein-protein encounter complexes of up to tens of microseconds. The activation barrier for individual ET events derived from the simulations was positive. Whereas the electrostatic interactions between P. laminosum PC and CytF are weak, simulations for a second cyanobacterial PC-CytF pair, that from Nostoc sp. PCC 7119, revealed ET rates influenced by stronger electrostatic interactions. In both cases, the simulations imply significant contributions to ET from both diffusion and thermal activation processes.
Collapse
Affiliation(s)
- Razif R Gabdoulline
- Molecular and Cellular Modeling Group, EML Research gGmbH, Schloss-Wolfsbrunnenweg 33, D-69118 Heidelberg, Germany.
| | | |
Collapse
|
50
|
Persson BA, Lund M. Association and electrostatic steering of alpha-lactalbumin-lysozyme heterodimers. Phys Chem Chem Phys 2009; 11:8879-85. [PMID: 20449034 DOI: 10.1039/b909179c] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The salt and pH dependent association of hen egg white lysozyme with alpha-lactalbumin whey proteins has been studied using molecular level Monte Carlo simulations. A highly uneven charge distribution of alpha-lactalbumin leads to strongly ordered heterodimers that may facilitate the formation of structured, mesoscopic aggregates. This electrostatic steering gives rise to 80% alignment at 5 mM 1 : 1 salt which, due to screening, diminishes to 60% at 100 mM salt. The free energy of interaction minima, dominated by electrostatics, ranges between -9 kT at 1 mM salt to -2 kT at 100 mM (neutral pH). Calculated osmotic second virial cross coefficients indicate complexation in the pH interval 6-10. Multivalent ions are found to effectively destabilize the protein complex and, at constant ionic strength, the order is La(3+) > Ca(2+) > Mg(2+) > Na(+). Upon binding of calcium to alpha-lactalbumin both the interaction and orientational alignment with lysozyme are reduced due to induced changes in the whey protein charge distribution. This potentially explains the experimentally observed absence of supramolecular structuring for the calcium loaded holo alpha-lactalbumin. Where available, good agreement is found with experimental data.
Collapse
Affiliation(s)
- Björn A Persson
- Department of Theoretical Chemistry, Chemical Center, POB 124, S-221 00 Lund, Sweden
| | | |
Collapse
|