Relating Collective Diffusion, Protein-Protein Interactions and

2 hours ago - To understand and predict the viscosities of highly concentrated monoclonal antibody (mAb) solutions, it is important to characterize th...
0 downloads 0 Views 698KB Size
Subscriber access provided by Nottingham Trent University

Bioengineering

Relating Collective Diffusion, Protein-Protein Interactions and Viscosity of Highly Concentrated Monoclonal Antibodies through Dynamic Light Scattering Barton J Dear, Amjad Chowdhury, Jessica J. Hung, Carl A. Karouta, Kishan Ramachandran, Maria P. Nieto, Logan R. Wilks, Ayush Sharma, Tony Y. Shay, Jason K Cheung, Thomas M. Truskett, and Keith P. Johnston Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.9b03432 • Publication Date (Web): 30 Aug 2019 Downloaded from pubs.acs.org on August 30, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Relating Collective Diffusion, Protein-Protein Interactions and Viscosity of Highly Concentrated Monoclonal Antibodies through Dynamic Light Scattering Barton J. Deara, Amjad Chowdhurya, Jessica J. Hunga, Carl A. Karoutaa, Kishan Ramachandrana, Maria P. Nietoa, Logan R. Wilksa, Ayush Sharmaa, Tony Y. Shaya, Jason K. Cheungc, Thomas M. Trusketta,b and Keith P. Johnston*a a) McKetta Department of Chemical Engineering, The University of Texas at Austin, Austin, TX 78712, United States b) Department of Physics, The University of Texas at Austin, Austin, TX 78712, United States c) Biophysical and Biochemical Characterization, Sterile Formulation Sciences, Merck & Co., Inc., Kenilworth, NJ 07033 USA

*200 E Dean Keeton St. Stop C0400, Austin, TX 78712 Phone: 512 471 4617 Fax: 512 471 7060 email: [email protected]

ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 36

Abstract To understand and predict the viscosities of highly concentrated monoclonal antibody (mAb) solutions, it is important to characterize the protein-protein interactions (PPI) and how they influence the possible formation of protein clusters. Herein, the collective diffusion is measured by dynamic light scattering (DLS) for solutions of three different mAbs, each with a various co-solutes to tune the PPI. The results are combined with measurements of static structure factor, (0), and self-diffusion coefficient, Ds, to understand the behavior of viscosity at high concentration. The small degree of variation in the hydrodynamic factor,

, and solvent protein friction,

, among systems with a

wide range of viscosities suggests solvent-protein interactions have a small influence on viscosity. Measurements of net PPI such as (0) and the diffusion interaction parameter,

, are predictive of high

concentration viscosity across co-solutes for a single mAb, but not for multiple mAbs. In contrast, properties that characterize the presence of clusters in solution, such as the polydispersity index from DLS,

, and the coefficient of protein-protein friction,

, exhibit stronger correlations to

viscosity.

1 ACS Paragon Plus Environment

Page 3 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Introduction Monoclonal antibodies (mAbs) may be used to treat a wide variety of diseases, for example, autoimmune diseases and cancer.1-2 In order to deliver a sufficient dose by subcutaneous administration, highly concentrated (100-300 mg/mL) solutions are required.3-5 Here the separation between mAbs becomes small (1 s for pulsed field gradient nuclear magnetic resonance (PFG-NMR)), the mAb is able to exit the surrounding cage during the course of the measurement leading to more complex dynamics.47-49 The long-time self-diffusion coefficient,

is often related to viscosity through the Generalized Stokes-

Einstein (GSE) equation (1)

= where

and

are the Boltzmann constant and temperature, and

and

are the hydrodynamic

radius and effective viscosity experienced by the tracer, respectively50. If the tracer particle is much smaller than the proteins in solution,

is equivalent to the solvent viscosity,

= 0.

Conversely,

becomes the macroscopic viscosity, , when the tracer particle is much larger than the protein. Thus, tracer particles, ranging from approximately 100 nm to 1 µm, are often used to measure the viscosity of protein solutions.12, 51-53 However, when the sizes of the tracer and mAb are on the same order (for example when the tracer is a labelled mAb),

=0


250 mM) to target a concentration of 250 mM at high concentration as done in previous work.13 After the final concentration was achieved the pH of both solutions was verified (Table 1). The reported SLS (except for mAb4), viscosity and FCS data were previously published as mAb2 and mAb3 SLS and viscosity data can be found in work by Hung et al.,41 mAb4 viscosity data in Dear, et al.,26 and mAb2 FCS data can be found in Hung et al.58 The mAb4 SLS measurements were performed identically to those for mAb2 and mAb3 by Hung et al.41 Dynamic Light Scattering (DLS) The diffusion coefficient of the mAb in each co-solute formulation was measured by DLS using a ZetaPALS zeta potential analyzer (Brookhaven Instruments, Holtsville, NY; P = 660 nm, R = 90°, and thus

= 0.013 nm-1), with some select samples also measured on a Malvern ZetaSizer Nano ZS

(Malvern Panalytical, Westborough, MA; P = 633 nm, R = 173°, and thus samples were sterile-filtered with 0.22 T

= 0.020 nm-1). MAb

PVDF filters (Ultrafree -MC-GV, MilliporeSigma,

Burlington, MA) prior to DLS measurements (to remove dust and particulates) and measured in triplicate, where each run consisted of 4 15-second scans that were averaged together. The recorded values are the average value of all the runs and the error bars indicate the measured standard deviations. The DLS ACFs, ) * were fit using Brookhaven Instruments’ Dynamic Light Scattering software using the quadratic cumulant algorithm77 to obtain

and the polydispersity index, 7

ACS Paragon Plus Environment

, by Eqn. 3

Page 9 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

)(*) = + ,

(-

* *

+

2*

2

2

2

)

(3)

where * is the correlation time, + is the fit baseline amplitude of ) * = 0), * is the fit characteristic mean correlation time with * = (

2 -1

)

, and

2

is the variance of the inverse correlation times fit to

quantify the effect of polydispersity on the diffusion with polydispersity index,

= *2 2.78 In

addition to the quadratic cumulants the ACFs were also fit with a Gaussian distribution model (GDM) (Eqn. 4) similar to what has been previously used to fit FCS data.58, 79 * 2

( )

(4a)

]

(4b)

]

(4c)

150

)(*) = +/1 = 1 1exp - *1

+1 = ) ,

[

/150 ) , 1=1

[

ln

*1 2

() *

- 2ln (7)2

ln

*1 2

() *

- 2ln (7)2 = 1

where *1 represents 150 different pre-selected correlation times (logarithmically spaced between 1 and 10000

), which are weighted by a Gaussian distribution with weight

1.

1

was determined by fitting

the average and standard deviation of characteristic correlation times, * and 7, respectively. ) is a normalization constant determined by Eqn. 4c. Representative ACFs with both the cumulants and GDM fits at ~200 and ~10 mg/mL are reported in Fig. S2 of the supporting information for samples of all three mAbs with either 250 mM Arg or 250 mM NaCl, as well as for mAb2 with 1000 mM Im. All DLS values reported are calculated from the quadratic cumulants method, unless specifically stated otherwise. All of the low concentration

values (concentrations up to the point at which

linearly with concentration) were fit to Eqn. 580-82 to obtain

0,

and

, where

began to change nonis the slope of

0

versus concentration at infinite dilution. ( )=

0(1

The values of

(5)

+ 0

and

from Eqn. 1 using

for each mAb solution are reported in Table S1 along with the =

=0

with

=0

9

calculated

being measured by a Cannon-Fenske viscometer.26

Theory 8 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

As derived in previous work54, 63, =

Page 10 of 36

can be calculated by Eqn. 6

(1 - :)2

where

(6)

is the solvent-protein friction per mole protein per mole solvent,

of solvent, : is the volume fraction of protein and

is the molar concentration

is the gas constant. To accurately calculate the

protein volume fraction the hydration layer was included83-85 as follows: :=

; < + =

(7)

?> + ; < + =

where ; is the protein mass and < is the volume of dry protein per mass (set to 0.7407 mL/g for this study),85 is the volume of water per mass (set to 1 mL/g), ?> is the total volume of unbound water, and = is the mass of hydration in terms of the mass of bound water per mass of protein). The value of = was set to 0.7 g/g for this study, which is higher than the 0.4 g/g previously used for BSA,85 because mAbs likely have more bound water than BSA due to their larger surface per volumes. Furthermore, this value produces : similar to those in previous studies.32, 41, 58, 86 Dividing the numerator and denominator of Eqn. 7 by 1 ml yields :=

(8)

< +=

At infinite dilution where : = 0 and (0) = 1, 0

=

=

9

where 9

9

(9)

=0

9

is the solvent-protein friction per mole solvent per mole protein at infinite dilution (with = 0@+),

=

and

9

is the hydrodynamic radius of the particle at infinite dilution.

Combining Eqns. 2, 6, and 9 shows that (0) =

is described by the Stokes-Einstein equation

(1 - :)2

= 0

9

= 0@+

can be directly related to

(1 - :)2

=

0

and :. (10)

Alternatively self-diffusion coefficients are related to both protein squared,

,

and protein-protein friction per mole

, as previously derived54, 63 (11)

+

9 ACS Paragon Plus Environment

Page 11 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

where

is the molar concentration of protein. At infinite dilution Eqn. 11 simplifies to Eqn. 9 as both

and

equal

0.

Assuming that

+

=

@+ at concentrations greater than zero

converts Eqn. 11 into the GSE equation (Eqn. 1). Combining Eqn. 6 and Eqn. 11 yields the following relation between =

and

+

(1 - :)2

=

. (1 - :)2 1

Previous studies have either assumed =

9

and

follows GSE yielding

(12) = 0, yielding the friction ratio, =

= 1,64, 87 or assumed that

=0

.54

The data will be compared to analytical expressions for hard sphere (HS) solutions. HS structure factors, (0) , hydrodynamic factors, (0) , normalized self-diffusion coefficients, and

9 0

, are

calculated from the Carnahan-Starling (Eqn. 13)49, the empirically determined Eqn. 1488 (which is shown to match simulations of (0) 89 - 90 in Fig. S3), and the Van Blaaderen (Eqn. 15)91, respectively. The normalized collective diffusion coefficients of hard spheres,

9 0

, is calculated by Eqn.

16 which comes from combining Eqns. 2, 13 and 14 as follows: (1 - :)4

(0)

= (1 +

(0)

= (1 - :)6.55

(13)

:)2 + :3(: - 4)

(14)

(1 - :)3

9 0

( ) 0

(15)

= 1 + '&0: + 2:2 + 3:3 = (1 - :)2.55[(1 + :)2 + :3(: - 4)]

(16)

Results and Discussion Collective Diffusion We begin with the results of

from DLS at low concentration to approximate the strength of

mAb PPI through kD 6, 9, 13, 15, 24, 34, 92. In Fig. 2, 0

the kD of each mAb solution was determined by fitting

0

versus concentration at low concentration using a weighted linear least squares minimization (further

10 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 36

details can be found in the supporting information). The concentration range for the fit was limited to concentration values less than the point at which the change in

0

began to deviate from linear behavior.

For mAb2 and mAb3, in Figs. 2A and 2B respectively, the mAb diffusion is faster (higher

) for

solutions with Arg and Im (particularly 1000 mM Im) relative to those with NaCl or Lys. This order matches that of relative viscosities,

C,

(Fig. S4 and Table 1) which are lowest for Arg and Im at a given

mAb concentration. The 250 mM Arg and NaCl solutions for each of these mAbs, and the more viscous mAb4, are presented in Fig. 2C. The

values still loosely follow the order of viscosity, as the

larger for NaCl relative to Arg solutions. Furthermore, the

is

values increase from mAb2 to mAb3 to

mAb4 (for relative viscosity ;+D > ;+D( > ;+D )as presented in Table 1. In Fig. 3 the measurements of of the mAb solutions

0

0

are extended to high concentration up to ~250 mg/mL. For all

initially decreases with concentration, but then plateaus around ~125 mg/mL up

to the highest concentrations measured (250 mg/mL), as has previously been seen for mAbs up to 150 mg/mL with high ionic strength.37 This behavior is vastly different from HS behavior (dashed lines in Fig. 3). This indicates that mAb diffusion is strongly affected by interactions beyond those of hard spheres, such as short-range attractions, including dipole interactions, and solvation contributions. As observed for

for a given mAb, higher values of

0

at the plateau (recorded in Table 1) corresponded

closely with the relative viscosities, but the differences between samples are larger at high concentration. It should be noted that high concentration DLS measurements could possibly be affected by multiple scattering93-94 which would cause the reported

0

to be higher than the true value.95-96 The effects of

multiple scattering were addressed by performing a few experiments in the back scattering mode at 173° in the Malvern Zetasizer Nano ZS instrument, where multiple scattering is reduced. As discussed further in the supporting information (Fig. S5), the effect of multiple scattering is relatively small as the measured values of

0

only vary by ~10%. The low degree of multiple scattering is consistent with the

small turbidities (Fig. S6), with 0.98 with the

slopes. Therefore, any trends seen comparing

slope to other properties in Table 1 are essentially the same using 7 slopes (data not shown). Table 1. List of physical properties measured for each mAb solution. Previous studies reported the values for J@K

,41 and J@K

+L

;GH1IG

( ) ;; +

,

.26

Cosolute (mM)

pH

C at 200 mg/ mL

50 NaCl 250 Arg 250 NaCl 250 Im 250 Lys 1000 Im 50 NaCl 250 Arg 250 NaCl 250 Im 250 Lys 250 Arg 250 NaCl

5.38 5.30 5.35 5.92 5.43 6.32 5.30 5.39 5.40 6.14 5.58 6.05 6.08

58 27 59 30 50 13 143 43 156 73 113 70 498

MN

410 135 555 143 500 125 mg/mL)

(0) at 200 mg/mL

0.44 0.57 0.39 0.55 0.41 0.84 0.51 0.61 0.38 0.56 0.39 0.42 0.17

0.21 0.14 0.27 0.18 0.24 0.12 0.16 0.12 0.13 0.15 0.17 0.15 ---

-0.14 ± 0.07 3.59 ± 0.05 -1.43 ± 0.11 0.80 ± 0.07 -1.46 ± 0.85 6.20 ± 0.16 -0.94 ± 0.77 -1.41 ± 0.22 -3.10 ± 0.13 -2.37 ± 0.26 -2.22 ± 0.44 -0.77 ± 0.35 -9.16 ± 1.03

;GH1IG

( ) ;; +

6.56 0.01 5.86 3.04 7.75 250 mM. Higher magnitudes indicate better correlations, while negative values suggest an inverse relationship, such as when attraction is strong

is lower, but viscosity is higher. Correlations are determined from values in Table 1.

24 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

Page 26 of 36

1 2 3 Correlation of all Data 4 J@K +L at / 0 at (0) at C at ;GH1IG 5 ** J@K ** 200 200 125 200 MN 22 (>125 ;; + Slope 6 mg/mL mg/mL*** mg/mL**** mg/mL) mg/mL* 7 C at 200 1.00 0.93 -0.67 -0.82 -0.74 -0.04 0.97 0.76 0.88 1.00 1.00 8 mg/mL 9 1.00 -0.82 -0.87 -0.85 0.09 0.94 0.75 0.92 0.99 0.97 MN 10 1.00 0.87 0.86 -0.10 -0.66 -0.48 -0.78 -0.79 -0.92 11 1.00 0.88 -0.32 -0.69 -0.72 -0.82 -1.00 -0.94 22 12 / 0 13 1.00 -0.55 -0.67 -0.80 -0.78 -0.97 -0.92 (>125 14 mg/mL) 15 (0) at 200 1.00 0.03 0.52 0.16 0.97 1.00 16 mg/mL* 17 ;GH1IG ** 1.00 0.80 0.71 0.93 --18 ;; + 19 J@K ** 1.00 0.51 0.95 --20 1.00 0.35 0.89 Slope 21 at 200 1.00 --22 mg/mL*** 23 J@K +L at 24 1.00 125 25 mg/mL**** 26 *Does not include mAb4 with 250 mM NaCl 27 **Does not include any mAb4 data 28 ***Only includes 3 mAb solutions all with mAb2 (250 mM Arg, 250 mM NaCl, and 1000 mM Im) 29 **** Only includes 3 mAb solutions (mAb2 with 250 mM Arg, mAb4 with 250 mM Arg and mAb4 with 250 mM NaCl) 30 31 As can be seen in Table 3, the parameters that best correlate with viscosity are and the 32 average aggregation number, J@K +L . J@K +L was determined from SAXS structure factor data at 125 33 34 mg/mL that were fit with 12-bead coarse grained molecular dynamics simulations.26 Since both of these 35 36 properties were only measured for three high ionic strength systems (all with mAb2 for ), it is not 37 38 yet fully known if they are better indicators of viscosity than the other measurements. However, because 39 40 they are directly influenced by the presence of clusters in solution there is reason to hypothesize that 41 42 they might be better indicators than , 22, 0 or (0), as many studies suggest that the presence of 43 44 clusters in mAb solutions elevate viscosity.7-10, 26, 39, 41 Recently mAb solution viscosities were related to 45 46 the cluster size distribution.11 Furthermore, this concept has been extended for cluster size distributions 47 from 12-bead model simulations that were fit to SAXS data.26 The is a measure of friction that is 48 49 directly affected by the presence of these clusters,63 and is thus strongly related to viscosity. The J@K +L 50 51 contains direct information about the size of clusters and is influenced not only by the net PPI in 52 53 solution, but also by the presence of specific local anisotropic short ranged attraction.26, 98 Alternatively, 54 55 properties based solely on , including (0), 22, and for the most part 0 (since (0) is mostly 56 57 58 25 59 ACS Paragon Plus Environment 60

( )

( )

Page 27 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

similar for the studied mAbs, Fig. 5) and

, are primarily indicative of net PPI. Computational studies

of mAbs,26 and other systems27, 103 have shown that large differences in cluster formation can occur between systems with the same net value of attraction. Specifically, mAbs with large but relatively weakly attractive regions form far fewer clusters in solution than those with small but strongly attractive regions.26 Furthermore, measures of net PPI can often be correlated strongly with viscosity for a given mAb13, 41 (as seen in Table S4 for correlations of just weakly interacting mAb2 solutions), but less strongly across mAbs (such as in Table 3), especially when nonuniform local attraction is prevalent on particular beads.13, 15, 41 It should also be noted that compared to the other techniques included in Table 3, the measurements of

is somewhat more unwieldy as it requires high concentration

measurements by DLS, SLS and FCS. It should be noted that

could also potentially be measured

by shear rheology given an accurate molecular viscoelasticity model for proteins since viscosity is a measurement of the friction coefficient on the molecular level as suggested by Einstein’s fluctuation dissipation theorem104. After

and J@K

the slope of the linear fit to

+L

the properties that best correlate with high concentration viscosity are and the oligomer ratio from SLS,

;GH1IG

( ) ;; +

. Both of these properties are

correlated more strongly with viscosity than the conventional low concentration properties The

and

22.

is a direct measurement of polydispersity from the presence of clusters in solution which

explains the strong correlation with viscosity. The oligomer ratio and average cluster size from SLS, J@K

, were calculated for the mAb2 and mAb3 systems in a previous study,41 and are determined from

fits of an interacting hard sphere model30-32, 105 (IHS) to the excess Rayleigh ratio, which is directly related to (0),41 across a wide range of mAb concentrations. The IHS model attributes the changes in concentration-dependent scattering to the formation of reversible oligomers, such that the scattering intensity at a given concentration corresponds to weighted contributions from a distribution of monomer, dimer and a characteristic higher order oligomer of fit aggregation number, M. The oligomer ratio is the mass fraction, ,, of the higher order oligomer divided by the sum of mass fractions of the monomer and dimer, while the J@K

is the average mass-based cluster size /1 = '9 9M(,11 . Even though (0) at a given

concentration is only indicative of net PPI, the cluster size distributions are related to (0) (or excess Rayleigh ratio) over a wide concentration range through the IHS model, which is affected by more than simply net PPI.26, 28, 41, 105 For example, for the same net PPI (0) has been shown to change more with 26 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 36

concentration for 12-bead models (representing mAbs) with small strong attraction sites on particular beads than for those with large weak ones.26 Additionally, Fig. S8 specifically demonstrates that the change of (0) with concentration is consistent with the trend in viscosity (mAb4 > mAb3 > mAb2). In terms of practicality both the

and the oligomer ratio only require a single type of measurement

(DLS and SLS, respectively), but determining the oligomer ratio involves more complicated calculations and requires many measurements to accurately determine the curvature of the Rayleigh ratio versus concentration. Conclusions In this work, the dynamic behavior of three mAbs was tuned with co-solute formulations and examined and characterized using DLS up to 250 mg/mL. The

determined at low concentrations with

DLS correlates loosely with the high concentration viscosity across multiple mAbs and co-solute systems, as has been seen previously.13, 15 As the protein concentration increases to ~125 mg/mL, reaches a plateau. Thus,

0

0

does not follow GSE, as it is additionally affected by the osmotic

compressibility (equivalent to (0)), and the solution micro-viscosity, which is intermediate between the solvent and solution viscosities. The plateau values correlate well with viscosities across multiple cosolute systems for a given mAb, but not across multiple mAbs. However, the best predictor of viscosity across multiple mAbs obtained solely from DLS is the

at high concentration, which characterizes

the decay times from diffusion of monomer and small reversible oligomers. The interpretation of the DLS data was more informative upon adding complementary static measurements by SLS,41 or SAXS26 as well as measurements of properties, including

,

,

from FCS58 to obtain a broad set of

. For example, the effects of thermodynamic and hydrodynamic

interactions on the observed diffusive behavior could be delineated. The variations in combining

0

from DLS and

, obtained by

from SLS,41 were small among systems at a given mAb concentration,

as these properties followed HS behavior. The values of the hydrodynamic property

were similar

despite variations in the thermodynamic PPI and widely different viscosities, as seen previously for silica particles with varying amounts of depletion attraction.97 Upon combining from DLS, the

and friction ratio (

-1

from FCS and

) could be determined. These properties were correlated

27 ACS Paragon Plus Environment

Page 29 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

strongly with the aggregation number of clusters as characterized by the J@K

+L

, J@K

,

;GH1IG

( ) ;; +

and

, as well as the viscosity. Thus, the results extend previous work in which viscosities of mAb

solutions were directly related to the presence of clusters, as characterized by structure factors via simulations of SAXS data,26 and even more specifically, the cluster size distributions.11, 26 In contrast, properties related to net PPI strength, such as

0

,

,

22,

and

are only loosely correlated with high

concentration viscosities as they do not capture cluster properties that depend upon the specific strength and location of the attractive regions on the mAbs.26 Acknowledgements The authors gratefully acknowledge Merck Sharp & Dohme Corp., a subsidiary of Merck & Co., Inc., Kenilworth, NJ, USA, Pfizer Inc. and the Welch Foundation (F-1319 and F-1696) for their financial support. The authors are particularly grateful to Dr. Chuck Ekert. The emphasis that Chuck placed on molecular understanding and molecular thermodynamics were very much appreciated for this study by his academic son KPJ and the academic grandchildren. Furthermore, the influence of Chuck on Pablo Debenedetti benefitted this work through TMT. What is it good for? Hopefully more stable proteins in therapeutics and improved protein delivery to aid human health. Supporting Information Comparison of DLS and SLS data for samples prepared using dialysis vs centrifugal diafiltration, representative fits to DLS ACFs, representative turbidities,

,

model comparisons, viscosity data and corresponding model fits, data, and additional correlation values are supplied as

supplemental information. References 1. Reichert, J. M. Therapeutic monoclonal antibodies approved or in review in the European Union or United States. http://www.antibodysociety.org/news/approved_mabs.php (accessed 1/15/2015). 2. Carter, P. J., Potent antibody therapeutics by design. Nature Reviews Immunology 2006, 6 (5), 343-357. 3. Shire, S. J.; Shahrokh, Z.; Liu, J., Challenges in the development of high protein concentration formulations. J. Pharm. Sci. 2004, 93 (6), 1390-1402. 4. Roberts, C. J., Protein aggregation and its impact on product quality. Current Opinion in Biotechnology 2014, 30, 211-217. 5. Wang, W., Advanced protein formulations. Protein Sci. 2015, 24 (7), 1031-9. 6. Yadav, S.; Laue, T. M.; Kalonia, D. S.; Singh, S. N.; Shire, S. J., The influence of charge distribution on selfassociation and viscosity behavior of monoclonal antibody solutions. Mol. Pharmaceutics 2012, 9 (4), 791-802. 28 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 30 of 36

7. Buck, P. M.; Chaudhri, A.; Kumar, S.; Singh, S. K., Highly Viscous Antibody Solutions are a Consequence of Network Formation Caused by Domain-Domain Electrostatic Complementarities: Insights from Coarse-Grained Simulations. Mol. Pharmaceutics 2015, 12 (1), 127-39. 8. Chaudhri, A.; Zarraga, I. E.; Kamerzell, T. J.; Brandt, J. P.; Patapoff, T. W.; Shire, S. J.; Voth, G. A., CoarseGrained Modeling of the Self-Association of Therapeutic Monoclonal Antibodies. J. Phys. Chem. B 2012, 116 (28), 8045-57. 9. Yearley, E. J.; Godfrin, P. D.; Perevozchikova, T.; Zhang, H.; Falus, P.; Porcar, L.; Nagao, M.; Curtis, J. E.; Gawande, P.; Taing, R.; Zarraga, I. E.; Wagner, N. J.; Liu, Y., Observation of Small Cluster Formation in Concentrated Monoclonal Antibody Solutions and its Implications to Solution Viscosity. Biophys. J. 2014, 106 (8), 1763-70. 10. Lilyestrom, W. G.; Yadav, S.; Shire, S. J.; Scherer, T. M., Monoclonal Antibody Self-Association, Cluster Formation, and Rheology at High Concentrations. J. Phys. Chem. B 2013, 117 (21), 6373-84. 11. Kastelic, M.; Dill, K. A.; Kalyuzhnyi, Y. V.; Vlachy, V., Controlling the viscosities of antibody solutions through control of their binding sites. J Mol Liq 2017, 234-242. 12. Schmit, J. D.; He, F.; Mishra, S.; Ketchem, R. R.; Woods, C. E.; Kerwin, B. A., Entanglement model of antibody viscosity. J. Phys. Chem. B 2014, 118 (19), 5044-9. 13. Dear, B. J.; Hung, J. J.; Truskett, T. M.; Johnston, K. P., Contrasting the influence of cationic amino acids on the viscosity and stability of a highly concentrated monoclonal antibody. Pharm Res 2017, 34 (1), 193-207. 14. Pindrus, M.; Shire, S. J.; Kelley, R. F.; Demeule, B.; Wong, R.; Xu, Y.; Yadav, S., Solubility Challenges in High Concentration Monoclonal Antibody Formulations: Relationship with Amino Acid Sequence and Intermolecular Interactions. Mol. Pharmaceutics 2015, 12 (11), 3896-3907. 15. Connolly, B. D.; Petry, C.; Yadav, S.; Demeule, B.; Ciaccio, N.; Moore, J. M.; Shire, S. J.; Gokarn, Y. R., Weak interactions govern the viscosity of concentrated antibody solutions: high-throughput analysis using the diffusion interaction parameter. Biophys. J. 2012, 103 (1), 69-78. 16. Liu, W.; Cellmer, T.; Keerl, D.; Prausnitz, J. M.; Blanch, H. W., Interactions of lysozyme in guanidinium chloride solutions from static and dynamic light-scattering measurements. Biotechnol. Bioeng. 2005, 90 (4), 482490. 17. Zarraga, I. E.; Taing, R.; Zarzar, J.; Luoma, J.; Hsiung, J.; Patel, A.; Lim, F. J., High Shear Rheology and Anisotropy in Concentrated Solutions of Monoclonal Antibodies. J. Pharm. Sci. 2013, 102 (8), 2538-49. 18. Borwankar, A. U.; Dear, B. J.; Twu, A.; Hung, J. J.; Dinin, A. K.; Wilson, B. K.; Yue, J.; Maynard, J. A.; Truskett, T. M.; Johnston, K. P., Viscosity reduction of a concentrated monoclonal antibody with arginine·HCl and arginine·glutamate. Ind. Eng. Chem. Res. 2016, 55 (43), 11225-11234. 19. Hung, J. J.; Borwankar, A. U.; Dear, B. J.; Truskett, T. M.; Johnston, K. P., High concentration tangential flow ultrafiltration of stable monoclonal antibody solutions with low viscosities. J. Membr. Sci. 2016, 508, 113126. 20. Hung, J. J.; Dear, B. J.; Dinin, A. K.; Borwankar, A. U.; Mehta, S. K.; Truskett, T. T.; Johnston, K. P., Improving viscosity and stability of a highly concentrated monoclonal antibody solution with concentrated proline. Pharm. Res. 2018, 35 (7), 133. 21. Guo, Z.; Chen, A.; Nassar, R. A.; Helk, B.; Mueller, C.; Tang, Y.; Gupta, K.; Klibanov, A. M., Structureactivity relationship for hydrophobic salts as viscosity-lowering excipients for concentrated solutions of monoclonal antibodies. Pharm. Res. 2012, 29 (11), 3102-3109. 22. Du, W.; Klibanov, A. M., Hydrophobic Salts Markedly Diminish Viscosity of Concentrated Protein Solutions. Biotechnol. Bioeng. 2011, 108 (3), 632-6. 23. Wang, S.; Zhang, N.; Hu, T.; Dai, W.; Feng, X.; Zhang, X.; Qian, F., Viscosity-Lowering Effect of Amino Acids and Salts on Highly Concentrated Solutions of Two IgG1 Monoclonal Antibodies. Mol. Pharmaceutics 2015, 12 (12), 4478-87. 24. Sorret, L. L.; DeWinter, M. A.; Schwartz, D. K.; Randolph, T. W., Challenges in Predicting Protein-Protein Interactions from Measurements of Molecular Diffusivity. Biophys. J. 2016, 111 (9), 1831-1842. 29 ACS Paragon Plus Environment

Page 31 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

25. Saito, S.; Hasegawa, J.; Kobayashi, N.; Kishi, N.; Uchiyama, S.; Fukui, K., Behavior of Monoclonal Antibodies: Relation Between the Second Virial Coefficient (B (2)) at Low Concentrations and Aggregation Propensity and Viscosity at High Concentrations. Pharm. Res. 2012, 29 (2), 397-410. 26. Dear, B. J.; Bollinger, J. A.; Chowdhury, A.; Hung, J. J.; Wilks, L. R.; Karouta, C. A.; Ramachandran, K.; Shay, T. Y.; Nieto, M. P.; Sharma, A.; Godfrin, P. D.; Truskett, T. M.; Johnston, K. P., X-Ray Scattering and CoarseGrained Simulations for Protein-Protein Interactions of Monoclonal Antibodies. J. Phys. Chem. B 2019, 123, 5274-5290. 27. Quang, L. J.; Sandler, S. I.; Lenhoff, A. M., Anisotropic contributions to protein-protein interactions. J. Chem. Theory Comput. 2014, 10 (2), 835-845. 28. Calero-Rubio, C.; Ghosh, R.; Saluja, A.; Roberts, C. J., Predicting protein-protein interactions of concentrated antibody solutions using dilute solution data and coarse-grained molecular models. J. Pharm. Sci. 2017, 107 (5), 1269-1281. 29. Ghosh, R.; Calero-Rubio, C.; Saluja, A.; Roberts, C. J., Relating Protein-Protein Interactions and Aggregation Rates From Low to High Concentrations. J. Pharm. Sci. 2016, 105 (3), 1086-96. 30. Scherer, T. M.; Liu, J.; Shire, S. J.; Minton, A. I., Intermolecular Interactions of IgG1 Monoclonal Antibodies at High Concentrations Characterized by Light Scattering. J. Phys. Chem. B 2010, 114 (40), 1294812957. 31. Scherer, T. M., Cosolute Effects on the Chemical Potential and Interactions of an IgG1 Monoclonal Antibody at High Concentrations. J. Phys. Chem. B 2013, 117 (8), 2254-66. 32. Scherer, T. M., Role of Cosolute–Protein Interactions in the Dissociation of Monoclonal Antibody Clusters. J. Phys. Chem. B. 2015, 119 (41), 13027-13038. 33. Wang, W.; Lilyestrom, W. G.; Hu, Z. Y.; Scherer, T. M., Cluster Size and Quinary Structure Determine the Rheological Effects of Antibody Self-Association at High Concentrations. J. Phys. Chem. B 2018, 122 (7), 21382154. 34. Mosbaek, C. R.; Konarev, P. V.; Svergun, D. I.; Rischel, C.; Vestergaard, B., High concentration formulation studies of an IgG2 antibody using small angle x-ray scattering. Pharm. Res. 2012, 29 (8), 2225-2235. 35. Corbett, D.; Hebditch, M.; Keeling, R.; Ke, P.; Ekizoglou, S.; Sarangapani, P.; Pathak, J.; Van Der Walle, C. F.; Uddin, S.; Baldock, C.; Avendaño, C.; Curtis, R. A., Coarse-Grained Modeling of Antibodies from Small-Angle Scattering Profiles. J. Phys. Chem. B. 2017, 121 (35), 8276-8290. 36. Castellanos, M. M.; Clark, N. J.; Watson, M. C.; Krueger, S.; McAuley, A.; Curtis, J. E., Role of Molecular Flexibility and Colloidal Descriptions of Proteins in Crowded Environments from Small-Angle Scattering. J. Phys. Chem. B 2016, 120 (49), 12511-12518. 37. Godfrin, P. D.; Zarzar, J.; Zarraga, I. E.; Porcar, L.; Falus, P.; Wagner, N. J.; Liu, Y., Effect of hierarchical cluster formation on the viscosity of concentrated monoclonal antibody formulations studied by neutron scattering. J. Phys. Chem. B. 2016, 120 (2), 278-291. 38. Yearley, E. J.; Zarraga, I. E.; Shire, S. J.; Scherer, T. M.; Gokarn, Y.; Wagner, N. J.; Liu, Y., Small-angle neutron scattering characterization of monoclonal antibody conformations and interactions at high concentrations. Biophys. J. 2013, 105 (3), 720-731. 39. Wang, G.; Varga, Z.; Hofmann, J.; Zarraga, I. E.; Swan, J. W., Structure and Relaxation in Solutions of Monoclonal Antibodies. J. Phys. Chem. B 2018, 122 (11), 2867-2880. 40. Riest, J.; Nagele, G.; Liu, Y.; Wagner, N. J.; Godfrin, P. D., Short-time dynamics of lysozyme solutions with competing short-range attraction and long-range repulsion: Experiment and theory. The Journal of Chemical Physics 2018, 148 (6), 065101. 41. Hung, J. J.; Dear, B. J.; Karouta, C. A.; Godfrin, P. D.; Bollinger, J. A.; Nieto, M. P.; Wilks, L. R.; Shay, T. Y.; Ramachandran, K.; Sharma, A.; Truskett, T. M.; Johnston, K. P., Protein-Protein Interactions of Highly Concentrated Monoclonal Antibody Solutions via Static Light Scattering and Influence on the Viscosity. J. Phys. Chem. B 2019, 123 (4), 739-755.

30 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 32 of 36

42. Bucciarelli, S.; Myung, J. S.; Farago, B.; Das, S.; Vliegenthart, G. A.; Holderer, O.; Winkler, R. G.; Schurtenberger, P.; Gompper, G.; Stradner, A., Dramatic influence of patchy attractions on short-time protein diffusion under crowded conditions. Science Advances 2016, 2 (12). 43. Roosen-Runge, F.; Hennig, M.; Zhang, F.; Jacobs, R. M.; Sztucki, M.; Schober, H.; Seydel, T.; Schreiber, F., Protein self-diffusion in crowded solutions. Proceedings of the National Academy of Sciences of the United States of America 2011, 108 (29), 11815-20. 44. Grimaldo, M.; Roosen-Runge, F.; Zhang, F.; Seydel, T.; Schreiber, F., Diffusion and dynamics of gammaglobulin in crowded aqueous solutions. J. Phys. Chem. B 2014, 118 (25), 7203-9. 45. Grimaldo, M.; Roosen-Runge, F.; Hennig, M.; Zanini, F.; Zhang, F.; Zamponi, M.; Jalarvo, N.; Schreiber, F.; Seydel, T., Salt-Induced Universal Slowing Down of the Short-Time Self-Diffusion of a Globular Protein in Aqueous Solution. J Phys Chem Lett 2015, 6 (13), 2577-82. 46. Bucciarelli, S.; Casal-Dujat, L.; De Michele, C.; Sciortino, F.; Dhont, J.; Bergenholtz, J.; Farago, B.; Schurtenberger, P.; Stradner, A., Unusual dynamics of concentration fluctuations in solutions of weakly attractive globular proteins. J Phys Chem Lett 2015, 6 (22), 4470-4. 47. Ochab-Marcinek, A.; Holyst, R., Scale-dependent diffusion of spheres in solutions of flexible and rigid polymers: mean square displacement and autocorrelation function for FCS and DLS measurements. Soft Matter 2011, 7 (16), 7366-7374. 48. Ochab-Marcinek, A.; Wieczorek, S. A.; Ziebacz, N.; Holyst, R., The effect of depletion layer on diffusion of nanoparticles in solutions of flexible and polydisperse polymers. Soft Matter 2012, 8 (43), 11173-11179. 49. Banchio, A. J.; Nagele, G., Short-time transport properties in dense suspensions: from neutral to chargestabilized colloidal spheres. The Journal of Chemical Physics 2008, 128 (10), 104903. 50. Kholodenko, A. L.; Douglas, J. F., Generalized Stokes-Einstein equation for spherical particle suspensions. Phys. Rev. E 1995, 51 (2), 1081-1090. 51. He, F.; Becker, G. W.; Litowski, J. R.; Narhi, L. O.; Brems, D. N.; Razinkov, V. I., High-throughput dynamic light scattering method for measuring viscosity of concentrated protein solutions. Anal. Biochem. 2010, 399 (1), 141-3. 52. He, F.; Woods, C.; Litowski, J.; Roschen, L.; Gadgil, H.; Razinkov, V.; Kerwin, B., Effect of sugar molecules on the viscosity of high concentration monoclonal antibody solutions. Pharm. Res. 2011, 28 (7), 1552-1560. 53. Josephson, L. L.; Furst, E. M.; Galush, W. J., Particle tracking microrheology of protein solutions. Journal of Rheology 2016, 60 (4), 531-540. 54. Nicoud, L.; Jagielski, J.; Pfister, D.; Lazzari, S.; Massant, J.; Lattuada, M.; Morbidelli, M., Kinetics of Monoclonal Antibody Aggregation from Dilute toward Concentrated Conditions. J. Phys. Chem. B 2016, 120 (13), 3267-80. 55. Kalwarczyk, T.; Sozanski, K.; Ochab-Marcinek, A.; Szymanski, J.; Tabaka, M.; Hou, S.; Holyst, R., Motion of nanoprobes in complex liquids within the framework of the length-scale dependent viscosity model. Adv Colloid Interface Sci 2015, 223, 55-63. 56. Lavalette, D.; Hink, M. A.; Tourbez, M.; Tétreau, C.; Visser, A. J., Proteins as micro viscosimeters: Brownian motion revisited. European Biophysics Journal 2006, 35 (6), 517-522. 57. Zorrilla, S.; Hink, M. A.; Visser, A. J.; Lillo, M. P., Translational and rotational motions of proteins in a protein crowded environment. Biophys. Chem. 2007, 125 (2-3), 298-305. 58. Hung, J. J.; Zeno, W. F.; Chowdhury, A. A.; Dear, B. J.; Ramachandran, K.; Nieto, M. P.; Shay, T. Y.; Karouta, C. A.; Hayden, C. C.; Cheung, J. K.; Truskett, T. M.; Stachowiak, J. C.; Johnston, K. P., Self-diffusion of a highly concentrated monoclonal antibody by fluorescence correlation spectroscopy: insight into protein–protein interactions and self-association. Soft Matter 2019, 15 (33), 6660-6676. 59. Wang, Y. Q.; Li, C. G.; Pielak, G. J., Effects of Proteins on Protein Diffusion. Journal of the American Chemical Society 2010, 132 (27), 9392-9397. 60. Roos, M.; Link, S.; Balbach, J.; Krushelnitsky, A.; Saalwachter, K., NMR-detected brownian dynamics of alphaB-crystallin over a wide range of concentrations. Biophys. J. 2015, 108 (1), 98-106. 31 ACS Paragon Plus Environment

Page 33 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

61. Roos, M.; Ott, M.; Hofmann, M.; Link, S.; Rossler, E.; Balbach, J.; Krushelnitsky, A.; Saalwachter, K., Coupling and Decoupling of Rotational and Translational Diffusion of Proteins under Crowding Conditions. Journal of the American Chemical Society 2016, 138 (32), 10365-72. 62. Balbo, J.; Mereghetti, P.; Herten, D. P.; Wade, R. C., The shape of protein crowders is a major determinant of protein diffusion. Biophys. J. 2013, 104 (7), 1576-84. 63. Vink, H., Mutual diffusion and self-diffusion in the frictional formalism of non-equilibrium thermodynamics. Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases 1985, 81 (7), 1725-1730. 64. Le Bon, C.; Nicolai, T.; Kuil, M. E.; Hollander, J. G., Self-Diffusion and Cooperative Diffusion of Globular Proteins in Solution. J. Phys. Chem. B. 1999, 103 (46), 10294-10299. 65. Dharmaraj, V. L.; Godfrin, P. D.; Liu, Y.; Hudson, S. D., Rheology of clustering protein solutions. Biomicrofluidics 2016, 10 (4), 043509. 66. Soraruf, D.; Roosen-Runge, F.; Grimaldo, M.; Zanini, F.; Schweins, R.; Seydel, T.; Zhang, F.; Roth, R.; Oettel, M.; Schreiber, F., Protein cluster formation in aqueous solution in the presence of multivalent metal ions-a light scattering study. Soft Matter 2014, 10 (6), 894-902. 67. Pan, W.; Galkin, O.; Filobelo, L.; Nagel, R. L.; Vekilov, P. G., Metastable mesoscopic clusters in solutions of sickle-cell hemoglobin. Biophys. J. 2007, 92 (1), 267-77. 68. Li, Y.; Lubchenko, V.; Vekilov, P. G., The use of dynamic light scattering and brownian microscopy to characterize protein aggregation. Rev. Sci. Instrum. 2011, 82 (5), 053106. 69. Vorontsova, M. A.; Chan, H. Y.; Lubchenko, V.; Vekilov, P. G., Lack of Dependence of the Sizes of the Mesoscopic Protein Clusters on Electrostatics. Biophys. J. 2015, 109 (9), 1959-68. 70. Castellanos, Maria M.; Pathak, Jai A.; Leach, W.; Bishop, Steven M.; Colby, Ralph H., Explaining the nonNewtonian character of aggregating monoclonal antibody aolutions using small-angle neutron scattering. Biophys. J. 2014, 107 (2), 469-476. 71. Neergaard, M. S.; Kalonia, D. S.; Parshad, H.; Nielsen, A. D.; Moller, E. H.; van de Weert, M., Viscosity of high concentration protein formulations of monoclonal antibodies of the IgG1 and IgG4 subclass - prediction of viscosity through protein-protein interaction measurements. Eur J Pharm Sci 2013, 49 (3), 400-10. 72. Borwankar, A. U.; Dinin, A. K.; Laber, J. R.; Twu, A.; Wilson, B. K.; Maynard, J. A.; Truskett, T. M.; Johnston, K. P., Tunable equilibrium nanocluster dispersions at high protein concentrations. Soft Matter 2013, 9 (6), 1766. 73. Johnston, K. P.; Maynard, J. A.; Truskett, T. M.; Borwankar, A.; Miller, M. A.; Wilson, B.; Dinin, A. K.; Khan, T. A.; Kaczorowski, K. J., Concentrated dispersions of equilibrium protein nanoclusters that reversibly dissociate into active monomers. ACS Nano 2012, 6 (2), 1357-1369. 74. Dear, B. J.; Hung, J. J.; Truskett, T. M.; Johnston, K. P., Contrasting the Influence of Cationic Amino Acids on the Viscosity and Stability of a Highly Concentrated Monoclonal Antibody. Pharm. Res. 2017, 34 (1), 193-207. 75. Dear, B. J.; Hung, J. J.; Laber, J. R.; Wilks, L. R.; Sharma, A.; Truskett, T. M.; Johnston, K. P., Enhancing Stability and Reducing Viscosity of a Monoclonal Antibody With Cosolutes by Weakening Protein-Protein Interactions. J. Pharm. Sci. 2019, 108 (8), 2517-2526. 76. Calero-Rubio, C.; Saluja, A.; Roberts, C. J., Coarse-grained antibody models for "weak" protein-protein interactions from low to high concentrations. J. Phys. Chem. B 2016, 120 (27), 6592-6605. 77. Hassan, P. A.; Rana, S.; Verma, G., Making Sense of Brownian Motion: Colloid Characterization by Dynamic Light Scattering. Langmuir 2014. 78. Hassan, P. A.; Rana, S.; Verma, G., Making sense of Brownian motion: colloid characterization by dynamic light scattering. Langmuir : the ACS journal of surfaces and colloids 2015, 31 (1), 3-12. 79. Pal, N.; Dev Verma, S.; Singh, M. K.; Sen, S., Fluorescence correlation spectroscopy: an efficient tool for measuring size, size-distribution and polydispersity of microemulsion droplets in solution. Anal Chem 2011, 83 (20), 7736-44.

32 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 34 of 36

80. Yadav, S.; Sreedhara, A.; Kanai, S.; Liu, J.; Lien, S.; Lowman, H.; Kalonia, D.; Shire, S., Establishing a Link Between Amino Acid Sequences and Self-Associating and Viscoelastic Behavior of Two Closely Related Monoclonal Antibodies. Pharm. Res. 2011, 28 (7), 1750-1764. 81. Yadav, S.; Liu, J.; Shire, S. J.; Kalonia, D. S., Specific Interactions in High Concentration Antibody Solutions Resulting in High Viscosity. J Pharm Sci 2010, 99 (3), 1152-1168. 82. Connolly, B D.; Petry, C.; Yadav, S.; Demeule, B.; Ciaccio, N.; Moore, J M. R.; Shire, S J.; Gokarn, Y R., Weak Interactions Govern the Viscosity of Concentrated Antibody Solutions: High-Throughput Analysis Using the Diffusion Interaction Parameter. Biophys. J. 2012, 103 (1), 69-78. 83. Squire, P. G.; Himmel, M. E., Hydrodynamics and protein hydration. Arch. Biochem. Biophys. 1979, 196 (1), 165-177. 84. Perkins, S. J., Protein volumes and hydration effects. European Journal of Biochemistry 1986, 157 (1), 169-180. 85. Sarangapani, P. S.; Hudson, S. D.; Migler, K. B.; Pathak, J. A., The Limitations of an Exclusively Colloidal View of Protein Solution Hydrodynamics and Rheology. Biophys. J. 2013, 105 (10), 2418-26. 86. Godfrin, P. D.; Zarraga, I. E.; Zarzar, J.; Porcar, L.; Falus, P.; Wagner, N. J.; Liu, Y., Effect of Hierarchical Cluster Formation on the Viscosity of Concentrated Monoclonal Antibody Formulations Studied by Neutron Scattering. J. Phys. Chem. B. 2016, 120 (2), 278-291. 87. Zettl, U.; Hoffmann, S. T.; Koberling, F.; Krausch, G.; Enderlein, J.; Harnau, L.; Ballauff, M., Self-Diffusion and Cooperative Diffusion in Semidilute Polymer Solutions As Measured by Fluorescence Correlation Spectroscopy. Macromolecules 2009, 42 (24), 9537-9547. 88. Margaret, R.; Annette, F.-T., Colloidal dispersions. Edited by W. B. Russel, D. A. Saville & W. R. Schowalter, Cambridge University Press, Cambridge, UK, 1989, xvii + 506 pp., price: £60.00. ISBN 0 521 34188 4. Journal of Chemical Technology & Biotechnology 1992, 54 (2), 201-202. 89. Ladd, A. J. C., Hydrodynamic transport coefficients of random dispersions of hard spheres. The Journal of Chemical Physics 1990, 93 (5), 3484-3494. 90. Segrè, P. N.; Behrend, O. P.; Pusey, P. N., Short-time Brownian motion in colloidal suspensions: Experiment and simulation. Phys. Rev. E 1995, 52 (5), 5070-5083. 91. van Blaaderen, A.; Peetermans, J.; Maret, G.; Dhont, J. K. G., 5 +1 " 61! 66 " of spherical colloidal particles measured with fluorescence recovery after photobleaching. The Journal of Chemical Physics 1992, 96 (6), 4591-4603. 92. Yadav, S.; Shire, S. J.; Kalonia, D. S., Factors affecting the viscosity in high concentration solutions of different monoclonal antibodies. J. Pharm. Sci. 2010, 99 (12), 4812-4829. 93. Pusey, P., Pusey, P. Suppression of multiple scattering by photon cross-correlation techniques. Curr. Opin. Colloid Interface Sci. 4, 177-185. 1999; Vol. 4, p 177-185. 94. Stieber, F.; Richtering, W., Fiber-Optic-Dynamic-Light-Scattering and Two-Color-Cross-Correlation Studies of Turbid, Concentrated, Sterically Stabilized Polystyrene Latex. Langmuir 1995, 11 (12), 4724-4727. 95. Thomas, J. C.; Tjin, S. C., Fiber optic dynamic light scattering (FODLS) from moderately concentrated suspensions. Journal of Colloid and Interface Science 1989, 129 (1), 15-31. 96. Dhont, J. K. G.; de Kruif, C. G.; Vrij, A., Light scattering in colloidal dispersions: Effects of multiple scattering. Journal of Colloid and Interface Science 1985, 105 (2), 539-551. 97. Ramakrishnan, S.; Shah, S. A.; Ruggeri, L.; Chen, Y. L.; Schweizer, K. S.; Zukoski, C. F., Collective diffusion in colloid-polymer suspensions: relative role of thermodynamics and hydrodynamics. Langmuir : the ACS journal of surfaces and colloids 2009, 25 (18), 10507-14. 98. Woldeyes, M. A.; Qi, W.; Razinkov, V. I.; Furst, E. M.; Roberts, C. J., How Well Do Low- and HighConcentration Protein Interactions Predict Solution Viscosities of Monoclonal Antibodies? J. Pharm. Sci. 2019, 108 (1), 142-154. 99. Varga, Z.; Wang, G.; Swan, J., The hydrodynamics of colloidal gelation. Soft Matter 2015, 11 (46), 90099019. 33 ACS Paragon Plus Environment

Page 35 of 36 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

100. Pathak, Jai A.; Sologuren, Rumi R.; Narwal, R., Do Clustering Monoclonal Antibody Solutions Really Have a Concentration Dependence of Viscosity? Biophys. J. 2013, 104 (4), 913-923. 101. Holyst, R.; Bielejewska, A.; Szymanski, J.; Wilk, A.; Patkowski, A.; Gapinski, J.; Zywocinski, A.; Kalwarczyk, T.; Kalwarczyk, E.; Tabaka, M.; Ziebacz, N.; Wieczorek, S. A., Scaling form of viscosity at all length-scales in poly(ethylene glycol) solutions studied by fluorescence correlation spectroscopy and capillary electrophoresis. Phys Chem Chem Phys 2009, 11 (40), 9025-32. 102. Zarraga, I. E.; Taing, R.; Zarzar, J.; Luoma, J.; Hsiung, J.; Patel, A.; Lim, F. J., High shear rheology and anisotropy in concentrated solutions of monoclonal antibodies. J. Pharm. Sci. 2013, 102 (8), 2538-2549. 103. Smallenburg, F.; Sciortino, F., Liquids more stable than crystals in particles with limited valence and flexible bonds. Nature Physics 2013, 9 (9), 554-558. 104. Marconi, U. M. B.; Puglisi, A.; Rondoni, L.; Vulpiani, A., Fluctuation–dissipation: Response theory in statistical physics. Physics Reports 2008, 461 (4), 111-195. 105. Minton, A. P., The effective hard particle model provides a simple, robust, and broadly applicable description of nonideal behavior in concentrated solutions of bovine serum albumin and other nonassociating proteins. J. Pharm. Sci. 2007, 96 (12), 3466-3469.

34 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Paragon Plus Environment

Page 36 of 36