Theory and simulations for RNA folding in mixtures of monovalent and divalent cations.

RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion-phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.

Suppression of Electrostatic Mediated Antibody Liquid-Liquid Phase Separation by Charged and Noncharged Preferentially Excluded Excipients.

Isotonic concentrations of inert cosolutes or excipients are routinely used in protein therapeutic formulations to minimize physical instabilities including aggregation, particulation, and precipitation that are often manifested during drug substance/product manufacture and long-term storage. Despite their prevalent use within the biopharmaceutical industry, a more detailed understanding for how excipients modulate the specific protein-protein interactions responsible for these instabilities is still needed so that informed formulation decisions can be made at the earliest stages of development when protein supply and time are limited. In the present report, subisotonic concentrations of the five common formulation excipients, sucrose, proline, sorbitol, glycerol, arginine hydrochloride, and the denaturant urea, were studied for their effect on the room temperature liquid-liquid phase separation of a model monoclonal antibody (mAb-B). Although each excipient lowered the onset temperatures of mAb-B liquid-liquid phase separation to different extents, all six were found to be preferentially excluded from the native state monomer by vapor pressure osmometry, and no apparent correlations to the excipient dependence of mAb-B melting temperatures were observed. These results and those of the effects of solution pH, addition of salt, and impact of a small number of charge mutations were most consistent with a mechanism of local excipient accumulation, to an extent dependent on their type, with the specific residues that mediate mAb-B electrostatic protein-protein interactions. These findings suggest that selection of excipients on the basis of their interaction with the solvent exposed residues of the native state may at times be a more effective strategy for limiting protein-protein interactions at pharmaceutically relevant storage conditions than choosing those that are excluded from the residues of the native state interior.

Machine Learning Models of Antibody-Excipient Preferential Interactions for Use in Computational Formulation Design.

Preferential interactions of formulation excipients govern their impact on the stability properties of proteins in solution. The ability to predict these interactions without the need to perform experiments would enable formulation design to begin early in the development of a new antibody therapeutic. With that in mind, we developed a feature set to numerically describe local regions of an antibody's surface for use in machine learning applications. Then, we used these features to train machine learning models for local antibody-excipient preferential interactions for the excipients sorbitol, sucrose, trehalose, proline, arginine·HCl, and NaCl. Our models had accuracies of up to about 85%. We also used linear (elastic net) models to quantify the contribution of antibody surface features to the preferential interaction coefficients, finding that the carbohydrates and proline tend to have similar important features, while the interactions of arginine·HCl and NaCl are governed by charge features. We present several case studies demonstrating how these machine learning models could be used to predict experimental aggregation and viscosity behavior in solution. Finally, we propose an approach to computational formulation design wherein a panel of excipients may be considered while designing an antibody sequence.

Molecular Computations of Preferential Interaction Coefficients of IgG1 Monoclonal Antibodies with Sorbitol, Sucrose, and Trehalose and the Impact of These Excipients on Aggregation and Viscosity.

Preferential interactions of formulation excipients govern their overall interactions with protein molecules, and molecular dynamics simulations allow for the examination of the interactions at the molecular level. We used molecular dynamics simulations to examine the interactions of sorbitol, sucrose, and trehalose with three different IgG1 antibodies to gain insight into how these excipients impact aggregation and viscosity. We found that sucrose and trehalose reduce aggregation more than sorbitol because of their larger size and their stronger interactions with high-spatial aggregation propensity residues compared to sorbitol. Two of the antibodies had high viscosity in sodium acetate buffer, and for these, we found that sucrose and trehalose tended to have opposite effects on viscosity. The data presented here provide further insight into the mechanisms of interactions of these three carbohydrate excipients with the antibody surface and thus their impact on excipient stabilization of antibody formulations.

Bilateral Effects of Excipients on Protein Stability: Preferential Interaction Type of Excipient and Surface Aromatic Hydrophobicity of Protein.

PURPOSE: Understanding the mechanism of protein-excipient interaction and illuminating the influencing factors on protein stability are key steps in the rational design of protein formulations. The objective of this study was to assess effects of preferential interaction type of excipient and surface aromatic hydrophobicity of protein on protein solution stability. METHODS: The preferential interaction between excipient and aromatic hydrophobic area of protein was investigated by solubility and fluorescence studies of amino acid derivatives in excipient solutions. We examined conformational, colloidal and mechanical stabilities of model proteins with different surface aromatic hydrophobicities, including bovine serum albumin (BSA) and ovalbumin (OVA), and then stability data were visualized by three-index empirical phase diagram. RESULTS: The result showed that preferentially excluded excipients (trehalose, sucrose and sorbitol) protected protein conformation against damage, but they could accelerate mechanical stress-induced aggregation. Preferentially bound excipients (propanediol and arginine) suppressed BSA aggregation, but arginine failed to inhibit OVA aggregation, which might be attributed to the disparate conformational perturbing effects of arginine on aromatic hydrophobic regions of BSA and OVA. CONCLUSIONS: These findings provided strong evidence that excipient possessed bilateral effects, and its application should be determined on different preferential interaction behaviors of excipients with protein, especially with the aromatic hydrophobic region.

Thermodynamics and solvent linkage of macromolecule-ligand interactions.

Binding involves two steps, desolvation and association. While water is ubiquitous and occurs at high concentration, it is typically ignored. In vitro experiments typically use infinite dilution conditions, while in vivo, the concentration of water is decreased due to the presence of high concentrations of molecules in the cellular milieu. This review discusses isothermal titration calorimetry approaches that address the role of water in binding. For example, use of D2O allows the contribution of solvent reorganization to the enthalpy component to be assessed. Further, the addition of osmolytes will decrease the water activity of a solution and allow effects on Ka to be determined. In most cases, binding becomes tighter in the presence of osmolytes as the desolvation penalty associated with binding is minimized. In other cases, the osmolytes prefer to interact with the ligand or protein, and if their removal is more difficult than shedding water, then binding can be weakened. These complicating layers can be discerned by different slopes in ln(Ka) vs osmolality plots and by differential scanning calorimetry in the presence of the osmolyte.

Thermodynamics of modulation of interaction of α-helix inducer 2, 2, 2-trifluoroethanol with lysozyme in presence of cationic, anionic and non-ionic surfactants.

The interactions of anionic sodium dodecyl sulphate (SDS), cationic cetyltrimethylammonium bromide (CTAB) and nonionic triton X-100 (TX-100) surfactants with lysozyme at pH = 2.4 have been studied individually as well as in combination with 2,2,2-trifluoroetanol (TFE). Urea has also been used in combination with surfactants. By using these combinations, efforts have been made to obtain partially folded conformations of the protein in the presence of surfactants and effect of α-helix inducer 2,2,2-trifluoroethanol on these intermediate states. Thermodynamic analysis of all these interactions has been done employing a combination of UV-visible, fluorescence and circular dichroism spectroscopies. The results have been correlated with each other and characterized qualitatively as well as quantitatively. At lower concentration of surfactant, the thermodynamic parameters indicated the destabilizing effect of SDS, stabilizing effect of CTAB and unappreciable destabilizing impact of TX-100 on lysozyme. The enhancement in destabilization effect or reduction in stabilization effect of surfactants on lysozyme in the presence of TFE and urea has also been indicated.Communicated by Ramaswamy H. Sarma.

Next generation multimodal chromatography resins via an iterative mapping approach: Chemical diversity, high-throughput screening, and chromatographic modelling.

Multimodal chromatography resins are becoming a key tool in the purification of biomolecules. The main objective of this research was the establishment of an iterative framework for the rapid development of new multimodal resins to provide novel selectivity for the future purification challenges. A large chemically diverse virtual library of 100 multimodal Capto™ MMC ligand analogues was created, and a broad array of chemical descriptors were calculated for each ligand in silico. Principal component analysis (PCA) was used to map the chemical diversity and guide selection of ligands for synthesis and coupling to the Capto ImpRes agarose base matrix. Twelve new ligands were prepared in two groups: 'group one' consist of L00-L07 and 'group two' consist of L08-L12. These ligands are diverse in the influence of varied secondary interactions such as hydrophobic interactions, H-bonding, etc. Additional resin prototypes were also prepared to look at the chromatographic impact of ligand density variation. High-throughput plate-based studies were performed for parallel resin screening for batch-binding of six model proteins at different chromatographic binding pH and sodium chloride concentration conditions. Principal component analysis of the binding data provided a chromatographic diversity map leading to the identification of ligands with improved binding. Further, the new ligands have improved separation resolution between a monoclonal antibody (mAb1) and product related impurities, a Fab fragment and high molecular weight (HMW) aggregates, using linear salt gradient elutions. To quantify the importance of secondary interactions, analysis of the retention factor of mAb1 on the ligands at various isocratic conditions lead to estimations of (a) the total number of water molecules and counter salt ions released during adsorption, and (b) hydrophobic contact area (HCA). The iterative mapping approach of chemical and chromatography diversity maps described in the paper proves to be a promising method for identifying new chromatography ligands for biopharmaceutical purification challenges.

The Role of Cosolvent-Water Interactions in Effects of the Media on Functionality of Enzymes: A Case Study of Photobacterium leiognathi Luciferase.

A complex heterogeneous intracellular environment seems to affect enzymatic catalysis by changing the mobility of biomolecules, their stability, and their conformational states, as well as by facilitating or hindering continuously occurring interactions. The evaluation and description of the influence of the cytoplasmic matrix components on enzymatic activity are problems that remain unsolved. In this work, we aimed to determine the mechanisms of action of two-component media with cosolvents of various molecular sizes on the complex multi-stage bioluminescent reaction catalyzed by bacterial luciferase. Kinetic and structural effects of ethylene glycol, glycerol, sorbitol, glucose, sucrose, dextran, and polyethylene glycol on bacterial luciferase were studied using stopped-flow and fluorescence spectroscopy techniques and molecular dynamics simulations. We have found that diffusion limitations in the presence of cosolvents promote the stabilization of flavin substrate and peroxyflavin intermediate of the reaction, but do not provide any advantages in bioluminescence quantum yield, because substrate binding is slowed down as well. The catalytic constant of bacterial luciferase has been found to be viscosity-independent and correlated with parameters of water-cosolvent interactions (Norrish constant, van der Waals interaction energies). Crowding agents, in contrast to low-molecular-weight cosolvents, had little effect on peroxyflavin intermediate decay and enzyme catalytic constant. We attributed specific kinetic effects to the preferential interaction of the cosolvents with enzyme surface and their penetration into the active site.

Functioning of a Fluorescein pH-Probe in Aqueous Media: Impact of Temperature and Viscosity.

In this work, we considered the influence of viscogenic agents (glycerol, sucrose) as well as the temperature on the fluorescent characteristics of fluorescein at pH 6.5 in order to describe the acid-base status of local environment in terms of a spectrally detectable dianion-anion equilibrium. The protolytic equilibrium of fluorescein was found to depend on the solvent viscosity in a complex way. Whereas in the presence of sucrose the ratiometric signal of fluorescein (I488/I435) remains rather unchanged, the addition of glycerol (up to 40% w/w) results in the increase of the signal (up to 19%), that can be attributed to the different mechanisms of cosolvents effects on dye molecules in the ground state. Molecular dynamics of the dye in the presence of glycerol and sucrose revealed that the cosolvents preferentially interact with fluorescein monoanion and dianion, displacing water molecules from the local environment which in turn reduces the average number of the hydrogen bonds between xanthene ring of the dye and water molecules. The ratiometric signal demonstrates linear growth with the temperature in the range of 10-80 °C regardless of the presence of viscogenic agents. A linear correlation between the temperature sensitivity of the ratiometric signal and the change in the molar enthalpy of the proton dissociation reaction in buffer and viscous media was determined.

Molecular computations of preferential interactions of proline, arginine.HCl, and NaCl with IgG1 antibodies and their impact on aggregation and viscosity.

Preferential interactions of excipients with the antibody surface govern their effect on the stability of antibodies in solution. We probed the preferential interactions of proline, arginine.HCl (Arg.HCl), and NaCl with three therapeutically relevant IgG1 antibodies via experiment and simulation. With simulations, we examined how excipients interacted with different types of surface patches in the variable region (Fv). For example, proline interacted most strongly with aromatic surfaces, Arg.HCl was included near negative residues, and NaCl was excluded from negative residues and certain hydrophobic regions. The differences in interaction of different excipients with the same surface patch on an antibody may be responsible for variations in the antibody's aggregation, viscosity, and self-association behaviors in each excipient. Proline reduced self-association for all three antibodies and reduced aggregation for the antibody with an association-limited aggregation mechanism. The effects of Arg.HCl and NaCl on aggregation and viscosity were highly dependent on the surface charge distribution and the extent of exclusion from highly hydrophobic patches. At pH 5.5, both tended to increase the aggregation of an antibody with a strongly positive charge on the Fv, while only NaCl reduced the aggregation of the antibody with a large negative charge patch on the Fv. Arg.HCl reduced the viscosities of antibodies with either a hydrophobicity-driven mechanism or a charge-driven mechanism. Analysis of this data presents a framework for understanding how amino acid and ionic excipients interact with different protein surfaces, and how these interactions translate to the observed stability behavior.

Protein Solvent Interaction: Transition of Protein-solvent Interaction Concept from Basic Research into Solvent Manipulation of Chromatography.

Previously, we have reviewed in this journal (Arakawa, T., Kita, Y., Curr. Protein Pept. Sci., 15, 608-620, 2014) the interaction of arginine with proteins and various applications of this solvent additive in the area of protein formulations and downstream processes. In this special issue, we expand the concept of protein-solvent interaction into the analysis of the effects of solvent additives on various column chromatography, including mixed-mode chromatography. Earlier in our research, we have studied the interactions of such a variety of solvent additives as sugars, salts, amino acids, polymers and organic solvents with a variety of proteins, which resulted in mechanistic understanding on their protein stabilization and precipitation effects, the latter known as Hofmeister series. While such a study was then a pure academic research, rapid development of genetic engineering technologies and resultant biotechnologies made it a valuable knowledge in fully utilizing solvent additives in manipulation of protein solution, including column chromatography.

Excluded Cosolvent in Chromatography.

The concept of cosolvent exclusion was developed by a group of Timasheff's laboratory in 1970-1990 and is currently used widely to explain the effects of a variety of cosolvents on the stability and solubility of macromolecules. Not surprisingly, these concepts have had substantial influence in the fields of formulation, protein folding and unfolding, but they have perhaps more surprisingly found their way into the field of chromatography. A variety of excluded cosolvents have been used to enhance binding and resolution of proteins and other macromolecules in ion exchange, hydroxyapatite, affinity, and hydrophobic interaction chromatography. These cosolvents include salting-out salts, amino acids and polymers, and frequently polyethylene glycol (PEG). A new mode of chromatography, termed "steric exclusion chromatography," was recently introduced. It employs hydroxylated solid phase surfaces. Steric exclusion of the PEG stabilizes the association of macromolecules with the solid phase. Elution is achieved by reducing the PEG concentration. Magnetic particles are also used in this chromatography. This review summarizes the concepts of preferential cosolvent exclusion and its applications in column chromatography.

Isotonic concentrations of excipients control the dimerization rate of a therapeutic immunoglobulin G1 antibody during refrigerated storage based on their rank order of native-state interaction.

Inert co-solutes, or excipients, are often included in protein biologic formulations to adjust the tonicity of liquid dosage forms intended for subcutaneous delivery. Despite the low concentration of their use, many of these excipients alter protein-protein interactions such as dimerization and aggregation rates of high concentration monoclonal antibody (mAb) therapeutics to varying extents during long-term refrigerated clinical storage, challenging the formulation scientist to make informed excipient selections at the earliest stages of development when protein supply and time are often limited. The objectives of this study were to better understand how isotonic concentrations of excipients influence the dimerization rates of a model mAb stored at refrigerated and room temperatures and explore protein sparing biophysical methods capable of predicting this dependence. Despite their prevalence of use in the biopharmaceutical industry, methods for assessing conformational stability such differential scanning calorimetry and isothermal equilibrium unfolding showed little predictive power and we highlight some of the assumptions and technical challenges of their use with mAbs. Conversely, measures of colloidal stability of the native-state such as preferential interaction coefficients measured by vapor pressure osmometry and solubility assessed by polyethylene-glycol induced precipitation correlated reasonably well with the mAb dimerization data and are most consistent with the excipients tested minimizing dimerization by interacting favorably with the residues comprising the protein-protein association interface.

The role of the concentration scale in the definition of transfer free energies.

The Gibbs free energy of transferring a solute at infinite dilution between two solvents quantifies differences in solute-solvent interactions - if the transfer takes place at constant molarity of the solute. Yet, many calculation formulae and measuring instructions that are commonly used to quantify solute-solvent interactions correspond to transfer processes in which not the molarity of the solute but its concentration measured in another concentration scale is constant. Here, we demonstrate that in this case, not only the change in solute-solvent interactions is quantified but also the entropic effect of a volume change during the transfer. Consequently, the "phenomenon" which is known as "concentration-scale dependence" of transfer free energies is simply explained by a volume-entropy effect. Our explanations are of high importance for the study of cosolvent effects on protein stability.

Lysozyme adsorption onto a cation-exchanger: mechanism of interaction study based on the analysis of retention chromatographic data.

In this study, based on the analysis of retention chromatographic data, we examined the adsorption of lysozyme onto carboxymethyl cellulose. Lysozyme retention data was collected at pH 5 and pH 8. The sodium chloride (NaCl) concentration in the mobile phase ranged from 300mM to 500mM and the temperature for this study varied from 288K to 308K. The retention measurements generated from these experimental conditions were analyzed with the Van't Hoff method, the preferential interaction model and the stoichiometric displacement model. Endothermic heats-of-adsorption and increases in entropy were observed under certain experimental conditions. These data suggest the presence of entropic driving forces such as the release of water and/or possibly structural changes in lysozyme molecules adsorbed to the surface of carboxymethyl cellulose. The modest observed exergonic adsorption ΔG° and the preferential interaction analysis corroborate the presence of water-release for this study. Additional analysis with the stoichiometric displacement model method revealed negligible changes in the structure of lysozyme molecules in contact with the surface of carboxymethyl cellulose.