Analysis of thermodynamic non-ideality in terms of protein solvation

The molecular basis for hydrodynamic properties of PEGylated human serum albumin

Polyethylene glycol (PEG) conjugation provides a protective modification that enhances the pharmacokinetics and solubility of proteins for therapeutic use. A knowledge of the structural ensemble of these PEGylated proteins is necessary to understand the molecular details that contribute to their hydrodynamic and colligative properties. Because of the large size and dynamic flexibility of pharmaceutically important PEGylated proteins, the determination of structure is challenging. In addition, the hydration of these conjugates that contain large polymers is difficult to determine with traditional methods that identify only first shell hydration water, which does not account for the complete hydrodynamic volume of a macromolecule. Here, we demonstrate that structural ensembles, generated by coarse-grained simulations, can be analyzed with HullRad and used to predict sedimentation coefficients and concentration-dependent hydrodynamic and diffusion nonideality coefficients of PEGylated proteins. A knowledge of these concentration-dependent properties enhances the ability to design and analyze new modified protein therapeutics. HullRad accomplishes this analysis by effectively accounting for the complete hydration of a macromolecule, including that of flexible polymers.

Experimental support for reclassification of the light scattering second virial coefficient from macromolecular solutions as a hydrodynamic parameter

This investigation examines the source of the disparity between experimental values of the light scattering second virial coefficient [Formula: see text] (mL.mol/g2) for proteins and those predicted on the statistical mechanical basis of excluded volume. A much better theoretical description of published results for lysozyme is obtained by considering the experimental parameters to monitor the difference between the thermodynamic excluded volume term and its hydrodynamic counterpart. This involves a combination of parameters quantifying concentration dependence of the translational diffusion coefficient obtained from dynamic light scattering measurements. That finding is shown to account for observations of a strong correlation between [Formula: see text] (mL/g), where M2 is the molar mass (molecular weight) of the macromolecule and the diffusion concentration parameter [Formula: see text] (mL/g). On the grounds that [Formula: see text] is regarded as a hydrodynamic parameter, the same status should be accorded the light scattering second virial coefficient rather than its current incorrect thermodynamic designation as [Formula: see text] (mL.mol/g2), or just B, the osmotic second virial coefficient for protein self-interaction.

Spectroscopic studies on the stability and nucleation-independent fibrillation of partially-unfolded proteins in crowded environment

Fibril formation of globular proteins is driven by attaining an appropriate partially-unfolded conformation. Excluded volume effect exerted by the presence of other macromolecules in the solution, as found in the cellular interior, might affect the conformational state of proteins and alter their fibril formation process. The change in structure, stability and rate of fibril formation of aggregation-prone partially-unfolded states of lysozyme (Lyz) and α-lactalbumin (ALA) in the presence of different sizes of polyethylene glycol (PEG) is examined using spectroscopic methods. Thermal denaturation and far-UV CD studies suggest that Lyz is stabilized by PEGs and the stability increases with increasing concentration of PEGs. However, the stability of ALA depends on the size and concentration of PEG. The change in enthalpy of unfolding indicates the existence of soft-interactions between the proteins and PEG along with excluded volume effect. Fibrillation rate of Lyz is not significantly altered in the presence of lower concentrations of PEGs suggesting that the crowding effect dominates the viscosity-induced retardation of protein association whereas at higher concentrations the rates are reduced. In case of ALA, the rate of fibrillation is drastically reduced; however, there is a marginal increase with the increasing concentration of PEG. The results suggest that the fibril formation is influenced by change in initial conformation of the partially-unfolded states of the proteins and their stability in the presence of the crowding agent. Further, the size and concentration of the crowding agent, and the soft-interaction between the proteins and PEG also affects the fibrillation.

Chromatographic parameter determination for complex biological feedstocks

The application of mechanistic models for chromatography requires accurate model parameters. Especially for complex feedstocks such as a clarified cell harvest, this can still be an obstacle limiting the use of mechanistic models. Another commonly encountered obstacle is a limited amount of sample material and time to determine all needed parameters. Therefore, this study aimed at implementing an approach on a robotic liquid handling system that starts directly with a complex feedstock containing a monoclonal antibody. The approach was tested by comparing independent experimental data sets with predictions generated by the mechanistic model using all parameters determined in this study. An excellent agreement between prediction and experimental data was found verifying the approach. Thus, it can be concluded that RoboColumns with a bed volume of 200 μL can well be used to determine isotherm parameters for predictions of larger scale columns. Overall, this approach offers a new way to determine crucial model input parameters for mechanistic modelling of chromatography for complex biological feedstocks. © 2018 The Authors Biotechnology Progress published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers Biotechnol. Prog., 34:1006–1018, 2018

Foreword to 'Quantitative and analytical relations in biochemistry'-a special issue in honour of Donald J. Winzor's 80th birthday

The purpose of this special issue is to honour Professor Donald J. Winzor's long career as a researcher and scientific mentor, and to celebrate the milestone of his 80th birthday. Throughout his career, Don has been renowned for his development of clever approximations to difficult quantitative relations governing a range of biophysical measurements. The theme of this special issue, 'Quantitative and analytical relations in biochemistry', was chosen to reflect this aspect of Don's scientific approach.

Concentration dependence of translational diffusion coefficients for globular proteins

This investigation examines published results of traditional diffusion experiments on ovalbumin and bovine serum albumin to determine the extent to which assumed concentration independence of the translational diffusion coefficient is a reasonable approximation in the analysis of boundary spreading in sedimentation velocity experiments on proteins.

Reassessment of the size of the supermolecular state of Dishevelled-3

Amendment of the interpretation of recently published size-exclusion chromatography data for Dishevelled-3 on Superdex 200 and Sephacryl S-400 has led to an increase in the estimated size of the supermolecular state from 2000 to 35 000 kDa, a value that essentially duplicates the redetermined and reported estimates obtained by fluorescence correlation spectroscopy on live cells. The earlier discrepancy between the sizes of the extensively aggregated form of this scaffold protein in vivo and in vitro is thereby eliminated.

Reactions of dipolar bio-molecules in nano-capsules--example of folding-unfolding process

The confinement of chemical reactions in nano-capsules can lead to a dramatic effect on the equilibrium constant of these latter. Indeed, capillary effects due to the curvature and surface energy of nano-capsules can alter in a noticeable way the evolution of reactions occurring within. Nano-encapsulation of bio-materials has attracted lately wide interest from the scientific community because of the great potential of its applications in biomedical areas and targeted therapies. The present paper focuses one's attention on alterations of conformation mechanisms due to extremely confining and interacting solvated dipolar macromolecules at their isoelectric point. As a specific example studied here, the folding-unfolding reaction of proteins (particularly RNase A and creatine kinase CK) is drastically changed when encapsulated in solid inorganic hollow nano-capsules. The effects demonstrated in this work can be extended to a wide variety of nano-encapsulation situations. The design and sizing of nano-capsules can even make use of the effects shown in the present study to achieve better and more effective encapsulation.

Allowance for the effect of protein charge in the characterization of nonideal solute self-association by sedimentation equilibrium

This theoretical investigation explores the use of statistical-mechanical approaches to characterize the reversible tetramerization of a protein monomer with the size and charge characteristics of serum albumin under conditions where consideration of nearest-neighbor interactions suffices to describe effects of thermodynamic non-ideality. Such analysis of simulated sedimentation equilibrium distributions points to the adequacy of both the scaled particle theory and potential-of-mean-force methods for determining the self-association constant. Although the latter method usually entails the assignment of a magnitude to monomer net charge, this requirement can be obviated to some extent by repeating the analysis for a range of monomer charges and identifying the most appropriate value as that associated with a minimum in the sum-of-squares-of-residuals (SSR) of the best-fit descriptions of the sedimentation equilibrium distribution. Reasonable estimates of the association constant are usually obtained from corresponding analyses of the same sedimentation equilibrium distributions with activity coefficients obtained by scaled particle theory, an approach which also involves the identification of parameters on the basis of a minimum in SSR. However, the value of monomer charge determined must be regarded as a curve-fitting parameter rather than a true measure of monomer charge. Similar qualifications are shown to prevail in the scaled particle theory approach, which also involves the identification of parameters (the effective monomer volume and the polymer/monomer volume ratio) on the basis of a minimum in SSR. We therefore recommend discontinuation of the practice whereby quite precise distinction between modes of self-association has been attempted on the grounds of the physical credibility of the magnitudes of these additional curve-fitting parameters.

The acacia gum arabinogalactan fraction is a thin oblate ellipsoid: a new model based on small-angle neutron scattering and ab initio calculation

Acacia gum is a branched complex polysaccharide whose main chain consists of 1,3-linked beta-D-galactopyranosyl units. Acacia gum is defined as a heteropolysaccharide since it contains approximately 2% of a polypeptide. The major molecular fraction (F1) accounting for approximately 88% of the total acacia gum mass is an arabinogalactan peptide with a weight-average molecular weight of 2.86 x 10(5) g/mol. The molecular structure of F1 is actually unknown. From small angle neutron scattering experiments in charge screening conditions, F1 appeared to be a dispersion of two-dimensional structures with a radius of gyration of approximately 6.5 nm and an inner dense branched structure. Inverse Fourier transform of F1 scattering form factor revealed a disk-like morphology with a diameter of approximately 20 nm and a thickness below 2 nm. Ab initio calculations on the pair distance distribution function produced a porous oblate ellipsoid particle with a central intricated "network". Both transmission electron microscopy and atomic force microscopy confirm the thin disk model and structural dimensions. The model proposed is a breakthrough in the field of arabinogalactan-protein-type macromolecules. In particular, concerning the site of biosynthesis of these macromolecules, the structural dimensions found in this study would be in agreement with a phloem-mediated long-distance transport. In addition, the structure of F1 could also explain the low viscosity of acacia gum solutions, and its ability to self-assemble and to interact with proteins.

Nonequivalence of second virial coefficients from sedimentation equilibrium and static light scattering studies of protein solutions

Experimental data for ovalbumin and lysozyme are presented to highlight the nonequivalence of second virial coefficients obtained for proteins by sedimentation equilibrium and light scattering. Theoretical considerations confirm that the quantity deduced from sedimentation equilibrium distributions is B(22), the osmotic second virial coefficient describing thermodynamic nonideality arising solely from protein self-interaction. On the other hand, the virial coefficient determined by light scattering is shown to reflect the combined contributions of protein-protein and protein-buffer interactions to thermodynamic nonideality of the protein solution. Misidentification of the light scattering parameter as B(22) accounts for published reports of negative osmotic second virial coefficients as indicators of conditions conducive to protein crystal growth. Finally, textbook assertions about the equivalence of second virial coefficients obtained by sedimentation equilibrium and light scattering reflect the restriction of consideration to single-solute systems. Although sedimentation equilibrium distributions for buffered protein solutions are, indeed, amenable to interpretation in such terms, the same situation does not apply to light scattering measurements because buffer constituents cannot be regarded as part of the solvent: instead they must be treated as non-scattering cosolutes.

Characterization of weak protein dimerization by direct analysis of sedimentation equilibrium distributions: the INVEQ approach

Closer scrutiny has been accorded a recently reported procedure for characterizing weak protein dimerization by sedimentation equilibrium (INVEQ) in which the equilibrium distribution is analyzed as a dependence of radial distance on solute concentration rather than of solute concentration on radial distance. By demonstrating theoretically that the fundamental parameter derived from the analysis is simply the difference between the dimerization constant and the osmotic second virial coefficient for monomer-monomer interaction, this investigation refutes the original claim that independent estimates of these two parameters can be obtained by nonlinear curve fitting of the sedimentation equilibrium distribution. This criticism also applies to conventional analyses of sedimentation distributions by the commonly employed Beckman Origin and NONLIN software. Numerically simulated distributions are then analyzed to demonstrate limitations of the procedure and also to indicate a means of improving the reliability of the returned estimate of the dimerization constant. These features are illustrated by applying the original and revised analytical procedures to a sedimentation equilibrium distribution for alpha-chymotrypsin (pH 4.0, I 0.05 M).

Negative second virial coefficients as predictors of protein crystal growth: evidence from sedimentation equilibrium studies that refutes the designation of those light scattering parameters as osmotic virial coefficients

The effects of ammonium sulphate concentration on the osmotic second virial coefficient (BAA/MA) for equine serum albumin (pH 5.6, 20 degrees C) have been examined by sedimentation equilibrium. After an initial steep decrease with increasing ammonium sulphate concentration, BAA/MA assumes an essentially concentration-independent magnitude of 8-9 ml/g. Such behaviour conforms with the statistical-mechanical prediction that a sufficient increase in ionic strength should effectively eliminate the contributions of charge interactions to BAA/MA but have no effect on the covolume contribution (8.4 ml/g for serum albumin). A similar situation is shown to apply to published sedimentation equilibrium data for lysozyme (pH 4.5). Although termed osmotic second virial coefficients and designated as such (B22), the negative values obtained in published light scattering studies of both systems have been described incorrectly because of the concomitant inclusion of the protein-salt contribution to thermodynamic nonideality of the protein. Those negative values are still valid predictors of conditions conducive to crystal growth inasmuch as they do reflect situations in which there is net attraction between protein molecules. However, the source of attraction responsible for the negative virial coefficient stems from the protein-salt rather than the protein-protein contribution, which is necessarily positive.

Molecular crowding effects of linear polymers in protein solutions

Measurement of protein-polymer second virial coefficients (BAP) by sedimentation equilibrium studies of carbonic anhydrase and cytochrome c in the presence of dextrans (T10-T80) has revealed an inverse dependence of BAP upon dextran molecular mass that conforms well with the behaviour predicted for the excluded-volume interaction between a spherical protein solute A and a random-flight representation of the polymeric cosolute P. That model of the protein-polymer interaction is also shown to provide a reasonable description of published gel chromatographic and equilibrium dialysis data on the effect of polymer molecular mass on BAP for human serum albumin in the presence of polyethylene glycols, a contrary finding from analysis of albumin solubility measurements being rejected on theoretical grounds. Inverse dependence upon polymer chainlength is also the predicted excluded-volume effect on the strength of several types of macromolecular equilibria-protein isomerization, protein dimerization, and 1:1 complex formation between dissimilar protein reactants. It is therefore concluded that published experimental observations of the reverse dependence, preferential reaction enhancement within DNA replication complexes by larger polyethylene glycols, must reflect the consequences of cosolute chemical interactions that outweigh those of thermodynamic nonideality arising from excluded-volume effects.

Effect of water potential on sol-gel transition and intermolecular interaction of gelatin near the transition temperature

The sol-gel transition of gelatin, measured by thermal analysis and viscosity measurement, was analyzed in terms of the change in hydration state of polymer molecules. A new thermodynamic model was proposed in which the effect of water potential is explicitly taken into account for the evaluation of the free energy change in the sol-gel transition process. Because of the large number of water molecules involved and the small free energy change in the transition process, the contribution of water activity, a(W), was proved to be not negligible in the sol-gel transition process in solutions containing such low-molecular cosolutes as sugars, glycerol, urea, and formamide. The gel-stabilization effect of sugars and glycerol was linear with a(W), which seemed consistent with the contribution of water potential in the proposed model. The different stabilization effect among sugars and glycerol was explained by the difference in solvent ordering, which affects hydrophobic interaction among protein molecules. The gel-destabilization effect of urea and formamide could be explained only by the direct binding of them to protein molecules through hydrogen bonding. On the contrary, the polymer-polymer interaction, measured by the viscosity analysis, in polyethyleneglycol and dextran solutions was not sensitive to the change in a(W), suggesting that no substantial change in hydration state with a(W) occurred in these polymer solutions.

Analytical exclusion chromatography

This review summarizes the development of exclusion chromatography, also termed gel filtration, molecular-sieve chromatography and gel permeation chromatography, for the quantitative characterization of solutes and solute interactions. As well as affording a means of determining molecular mass and molecular mass distribution, the technique offers a convenient way of characterizing solute self-association and solute-ligand interactions in terms of reaction stoichiometry and equilibrium constant. The availability of molecular-sieve media with different selective porosities ensures that very little restriction is imposed on the size of solute amenable to study. Furthermore, access to a diverse array of assay procedures for monitoring the column eluate endows analytical exclusion chromatography with far greater flexibility than other techniques from the viewpoint of solute concentration range that can be examined. In addition to its widely recognized prowess as a means of solute separation and purification, exclusion chromatography thus also possesses considerable potential for investigating the functional roles of the purified solutes.

Allowance for thermodynamic non-ideality in the characterization of protein self-association by frontal exclusion chromatography: hemoglobin revisited

This investigation re-examines theoretical aspects of the allowance for effects of thermodynamic non-ideality on the characterization of protein self-association by frontal exclusion chromatography, and thereby provides methods of analysis with greater thermodynamic rigor than those used previously. Their application is illustrated by reappraisal of published exclusion chromatography data for hemoglobin on the controlled-pore-glass matrix CPG-120. The equilibrium constant of 100/M that is obtained for dimerization of the alpha(2)beta(2) species by this means is also deduced from re-examination of published studies of concentrated hemoglobin solutions by osmotic pressure and sedimentation equilibrium methods.

Thermodynamic non-ideality as an alternative source of the effect of sucrose on the thrombin-catalyzed hydrolysis of peptide p-nitroanilide substrates

The inhibitory effect of sucrose on the kinetics of thrombin-catalyzed hydrolysis of the chromogenic substrate S-2238 (D-phenylalanyl-pipecolyl-arginoyl-p-nitroanilide) is re-examined as a possible consequence of thermodynamic non-ideality-an inhibition originally attributed to the increased viscosity of reaction mixtures. However, those published results may also be rationalized in terms of the suppression of a substrate-induced isomerization of thrombin to a slightly more expanded (or more asymmetric) transition state prior to the irreversible kinetic steps that lead to substrate hydrolysis. This reinterpretation of the kinetic results solely in terms of molecular crowding does not signify the lack of an effect of viscosity on any reaction step(s) subject to diffusion control. Instead, it highlights the need for development of analytical procedures that can accommodate the concomitant operation of thermodynamic non-ideality and viscosity effects.

The hydration problem in solution biophysics: an introduction

The background and purpose of the British Biophysical Society Discussion meeting, The Hydration Problem in Solution Biophysics, held at the University of Glasgow, 12 September 2000, is described, particularly in relation to previous meetings in this field. Whereas a study of the nature and dynamic properties of water associated with a molecule is an important topic by itself, the collection of papers based on this meeting focus mainly on its affect in interpreting biophysical data in terms of macromolecular shape in a solution environment, particularly under dilute and very dilute systems. The techniques considered are largely hydrodynamically or thermodynamically based and supplemented by molecular modeling strategies; but in the context of how these could be used in conjunction with techniques like X-ray crystallography, NMR and neutron scattering.