Molecular Mechanism for the Hofmeister Effect Derived from NMR and DSC Measurements on Barnase

Interactions between the protein barnase and co-solutes studied by NMR

Protein solubility and stability depend on the co-solutes present. There is little theoretical basis for selection of suitable co-solutes. Some guidance is provided by the Hofmeister series, an empirical ordering of anions according to their effect on solubility and stability; and by osmolytes, which are small organic molecules produced by cells to allow them to function in stressful environments. Here, NMR titrations of the protein barnase with Hofmeister anions and osmolytes are used to measure and locate binding, and thus to separate binding and bulk solvent effects. We describe a rationalisation of Hofmeister (and inverse Hofmeister) effects, which is similar to the traditional chaotrope/kosmotrope idea but based on solvent fluctuation rather than water withdrawal, and characterise how co-solutes affect protein stability and solubility, based on solvent fluctuations. This provides a coherent explanation for solute effects, and points towards a more rational basis for choice of excipients.

Hydrogen-Deuterium Exchange Dynamics of NISTmAb Measured by Small Angle Neutron Scattering

Understanding protein dynamics and conformational stability holds great significance in biopharmaceutical research. Hydrogen-deuterium exchange (HDX) is a quantitative methodology used to examine these fundamental properties of proteins. HDX involves measuring the exchange of solvent-accessible hydrogens with deuterium, which yields valuable insights into conformational fluctuations and conformational stability. While mass spectrometry is commonly used to measure HDX on the peptide level, we explore a different approach using small-angle neutron scattering (SANS). In this work, SANS is demonstrated as a complementary and noninvasive HDX method (HDX-SANS). By assessing subtle changes in the tertiary and quaternary structure during the exchange process in deuterated buffer, along with the influence of added electrolytes on protein stability, SANS is validated as a complementary HDX technique. The HDX of a model therapeutic antibody, NISTmAb, an IgG1κ, is monitored by HDX-SANS over many hours using several different formulations, including salts from the Hofmeister series of anions, such as sodium perchlorate, sodium thiocyanate, and sodium sulfate. The impact of these formulation conditions on the thermal stability of NISTmAb is probed by differential scanning calorimetry. The more destabilizing salts led to heightened conformational dynamics in mAb solutions even at temperatures significantly below the denaturation point. HDX-SANS is demonstrated as a sensitive and noninvasive technique for quantifying HDX kinetics directly in mAb solution, providing novel information about mAb conformational fluctuations. Therefore, HDX-SANS holds promise as a potential tool for assessing protein stability in formulation.

The protein escape process at the ribosomal exit tunnel has conserved mechanisms across the domains of life

The ribosomal exit tunnel is the primary structure affecting the release of nascent proteins at the ribosome. The ribosomal exit tunnels from different species have elements of conservation and differentiation in structural and physico-chemical properties. In this study, by simulating the elongation and escape processes of nascent proteins at the ribosomal exit tunnels of four different organisms, we show that the escape process has conserved mechanisms across the domains of life. Specifically, it is found that the escape process of proteins follows the diffusion mechanism given by a simple diffusion model, and the median escape time positively correlates with the number of hydrophobic residues and the net charge of a protein for all the exit tunnels considered. These properties hold for 12 distinct proteins considered in two slightly different and improved Gō-like models. It is also found that the differences in physico-chemical properties of the tunnels lead to quantitative differences in the protein escape times. In particular, the relatively strong hydrophobicity of E. coli's tunnel and the unusually high number of negatively charged amino acids on the tunnel's surface of H. marismortui lead to substantially slower escapes of proteins at these tunnels than at those of S. cerevisiae and H. sapiens.

Analytical Workflows to Unlock Predictive Power in Biotherapeutic Developability

PURPOSE

Forming accurate data models that assist the design of developability assays is one area that requires a deep and practical understanding of the problem domain. We aim to incorporate expert knowledge into the model building process by creating new metrics from instrument data and by guiding the choice of input parameters and Machine Learning (ML) techniques.

METHODS

We generated datasets from the biophysical characterisation of 5 monoclonal antibodies (mAbs). We explored combinations of techniques and parameters to uncover the ones that better describe specific molecular liabilities, such as conformational and colloidal instability. We also employed ML algorithms to predict metrics from the dataset.

RESULTS

We found that the combination of Differential Scanning Calorimetry (DSC) and Light Scattering thermal ramps enabled us to identify domain-specific aggregation in mAbs that would be otherwise overlooked by common developability workflows. We also found that the response to different salt concentrations provided information about colloidal stability in agreement with charge distribution models. Finally, we predicted DSC transition temperatures from the dataset, and used the order of importance of different metrics to increase the explainability of the model.

CONCLUSIONS

The new analytical workflows enabled a better description of molecular behaviour and uncovered links between structural properties and molecular liabilities. In the future this new understanding will be coupled with ML algorithms to unlock their predictive power during developability assessment.

The measurement of binding affinities by NMR chemical shift perturbation

We have carried out chemical shift perturbation titrations on three contrasting proteins. The resulting chemical shifts have been analysed to determine the best way to fit the data, and it is concluded that a simultaneous fitting of all raw shift data to a single dissociation constant is both the most accurate and the most precise method. It is shown that the optimal weighting of ¹⁵N chemical shifts to ¹H chemical shifts is protein dependent, but is around the consensus value of 0.14. We show that chemical shift changes of individual residues can be fit to give residue-specific affinities. Residues with affinities significantly stronger than average are found in close contact with the ligand and are suggested to form a rigid contact surface, but only when the binding involves little conformational change. This observation may be of value in analysing binding and conformational change.

Modulation of global stability, ligand binding and catalytic properties of trypsin by anions

Specific salts effect is well-known on stability and solubility of proteins, however, relatively limited knowledge is known regarding the effect on catalytic properties of enzymes. Here, we examined the effect of four sodium anions on thermal stability and catalytic properties of trypsin and binding of the fluorescent probe, p-aminobenzamidine (PAB), to the enzyme. We show that the specific anions effect on trypsin properties agrees with the localization of the anions in the Hofmeister series. Thermal stability of trypsin, T_m, the affinity of the fluorescent probe to the binding site, K_d, and the rate constant, k_cat, of trypsin-catalyzed hydrolysis of the substrate N-benzoyl-L-arginine ethyl ester (BAEE) increase with increasing kosmotropic character of anions in the order: perchlorate<bromide<chloride<sulfate, while the value of Michaelis constant, K_M, decreases. Correlations between the values of T_m, K_d for PAB, k_cat, and K_M for BAEE in the presence of 1 M studied salts suggest interrelation among these parameters of the enzyme. Global stabilization as well as increased rigidity of trypsin is accompanied by strengthening of interaction with fluorescent probe PAB and in accordance with decreasing values of K_M for the substrate BAEE. Strong correlations between parameters characterizing the trypsin properties with the charge densities of anions clearly indicate direct electrostatic interaction as a basis of the specific anion effect on the conformational and functional properties of the enzyme.

When Molecules Meet in Water-Recent Contributions of Supramolecular Chemistry to the Understanding of Molecular Recognition Processes in Water

Molecular recognition processes in water differ from those in organic solvents in that they are mediated to a much greater extent by solvent effects. The hydrophobic effect, for example, causes molecules that only weakly interact in organic solvents to stay together in water. Such water-mediated interactions can be very efficient as demonstrated by many of the synthetic receptors discussed in this review, some of which have substrate affinities matching or even surpassing those of natural binders. However, in spite of considerable success in designing such receptors, not all factors determining their binding properties in water are fully understood. Existing concepts still provide plausible explanations why the reorganization of water molecules often causes receptor-substrate interactions in water to be strongly exothermic rather than entropically favored as predicted by the classical view of the hydrophobic effect.

Hydrophobic and electrostatic interactions modulate protein escape at the ribosomal exit tunnel

After translation, nascent proteins must escape the ribosomal exit tunnel to attain complete folding to their native states. This escape process also frees up the ribosome tunnel for a new translation job. In this study, we investigate the impacts of energetic interactions between the ribosomal exit tunnel and nascent proteins on the protein escape process by molecular dynamics simulations using partially coarse-grained models that incorporate hydrophobic and electrostatic interactions of the ribosome tunnel of Haloarcula marismortui with nascent proteins. We find that, in general, attractive interactions slow down the protein escape process, whereas repulsive interactions speed it up. For the small globular proteins considered, the median escape time correlates with both the number of hydrophobic residues, N_h, and the net charge, Q, of a nascent protein. A correlation coefficient exceeding 0.96 is found for the relation between the median escape time and a combined quantity of N_h + 5.9Q, suggesting that it is ∼6 times more efficient to modulate the escape time by changing the total charge than the number of hydrophobic residues. The estimated median escape times are found in the submillisecond-to-millisecond range, indicating that the escape does not delay the ribosome recycling. For various types of the tunnel model, with and without hydrophobic and electrostatic interactions, the escape time distribution always follows a simple diffusion model that describes the escape process as a downhill drift of a Brownian particle, suggesting that nascent proteins escape along barrier-less pathways at the ribosome tunnel.

Spectroscopic Evidence of the Salt-Induced Conformational Change around the Localized Electric Charges on the Protein Surface of Fibronectin Type III

The effect of salt on the electrostatic interaction of a protein is an important issue, because addition of salt affects protein stability and association/aggregation. Although adding salt is a generally recognized strategy to improve protein stability, this improvement does not necessarily occur. The lack of an effect upon the addition of salt was previously confirmed for the tenth fibronectin type III domain from human fibronectin (FN3) by thermal stability analysis. However, the detailed molecular mechanism is unknown. In the present study, by employing the negatively charged carboxyl triad on the surface of FN3 as a case study, the molecular mechanism of the inefficient NaCl effect on protein stability was experimentally addressed using spectroscopic methods. Complementary analysis using Raman spectroscopy and 8-anilino-1-naphthalenesulfonic acid fluorescence revealed the three-phase behavior of the salt-protein interaction between NaCl and FN3 over a wide salt concentration range from 100 mM to 4.0 M, suggesting that the Na⁺-specific binding to the negatively charged carboxyl triad causes a local conformational change around the binding site with an accompanying structural change in the overall protein, which contributes to the protein's structural destabilization. This spectroscopic evidence clarifies the molecular understanding of the inefficiency of salt to improve protein stability. The findings will inform the optimization of formulation conditions.

Anion binding to ubiquitin and its relevance to the Hofmeister effects

Although the non-covalent interactions between proteins and salts contributing to the Hofmeister effects have been generally mapped, there are many questions regarding the specifics of these interactions. We report here studies involving the small protein ubiquitin and salts of polarizable anions. These studies reveal a complex interplay between the reverse Hofmeister effect at low pH, the salting-in Hofmeister effect at higher pH, and six anion binding sites in ubiquitin at the root of these phenomena. These sites are all located at protuberances of preorganized secondary structure, and although stronger at low pH, are still apparent when ubiquitin possesses no net charge. These results demonstrate the traceability of these Hofmeister phenomena and suggest new strategies for understanding the supramolecular properties of proteins.

Binding of excipients is a poor predictor for aggregation kinetics of biopharmaceutical proteins

One of the major challenges in formulation development of biopharmaceuticals is improving long-term storage stability, which is often achieved by addition of excipients to the final formulation. Finding the optimal excipient for a given protein is usually done using a trial-and-error approach, due to the lack of general understanding of how excipients work for a particular protein. Previously, preferential interactions (binding or exclusion) of excipients with proteins were postulated as a mechanism explaining diversity in the stabilisation effects. Weak preferential binding is however difficult to quantify experimentally, and the question remains whether the formulation process should seek excipients which preferentially bind with proteins, or not. Here, we apply solution NMR spectroscopy to comprehensively evaluate protein-excipient interactions between therapeutically relevant proteins and commonly used excipients. Additionally, we evaluate the effect of excipients on thermal and colloidal protein stability, on aggregation kinetics and protein storage stability at elevated temperatures. We show that there is a weak negative correlation between the strength of protein-excipient interactions and effect on enhancing protein thermal stability. We found that the overall protein-excipient binding per se can be a poor criterion for choosing excipients enhancing formulation stability. Experiments on a diverse set of excipients and test proteins reveal that while excipients affect all of the different aspects of protein stability, the effects are very much protein specific, and care must be taken to avoid apparent generalisations if a smaller dataset is being used.

Salt Effects on Hydrophobic Solvation: Is the Observed Salt Specificity the Result of Excluded Volume Effects or Water Mediated Ion-Hydrophobe Association?

The solubility of hydrophobic molecules in water is sensitive to salt addition in an ion-specific manner. Such "salting-out" and "salting-in" properties have been shown to be a major contributor to the measured ion-specific Hofmeister effects that are observed in many biophysical phenomena. Various theoretical models have suggested a number of disparate mechanisms for salting-out (salting-in) of hydrophobic moieties, the most popular of which include preferential interaction, water-mediated association, and electrostriction models. However, a complete molecular level description of this ion-specificity is not yet available. This work investigates the ion-specific nature of hydrophobic solvation by studying how sodium and chloride salts affect the thermodynamics of 1,2-hexanediol micellization. The results of this study are analyzed in terms of scaled-particle theory and we show that salt addition can affect hydrophobic solvation in two modalities: salt addition changes the cavitation free energy; salt addition also influences the solvent-solute interaction energy by changing the hydration of the hydrophobic solute. These two effects are salt specific in nature and we suggest that for small hydrophobic solutes these effects are the main cause of salt-specific Hofmeister effects on their solubility.

Quantification of 24 circulating endocannabinoids, endocannabinoid-related compounds, and their phospholipid precursors in human plasma by UHPLC-MS/MS

Endocannabinoids and endocannabinoid-related compounds (ERCs) are involved in many physiological processes. They are released on demand from phosphoinositide and N-acylphosphatidyl ethanolamine (NAPE) precursors and comprise 2-monoacylglycerols (2-MGs) and FA ethanolamides (FEAs). Despite the abundance of advanced quantitative methods, however, their determined concentrations in blood plasma are inconsistent because 2-MGs and FEAs undergo artifactual de novo formation, chemical isomerization, and degradation during sample collection and storage. For a comprehensive survey of these compounds in blood and plasma, we have developed and validated an ultra-HPLC-MS/MS method to quantify 24 endocannabinoids, ERCs, and their phospholipid precursors. Immediate acidification of EDTA-blood to pH 5.8 blocked artifactual FEA formation for at least 4 h on ice. The 2-MGs were stabilized after plasma harvest with 0.5 M potassium thiocyanate at pH 4.7. FEA and MG plasma concentrations in six healthy volunteers ranged between 0.04-3.48 and 0.63-6.18 ng/ml, respectively. Interestingly, only 1-5% of circulating FEAs were present in their free form, while the majority was bound to NAPEs. Similarly, 97% of 2-arachidonoylglycerol (2-AG) was bound to a potential phosphoinositide pool. The herein-described stabilization and extraction methods may now be used to reliably and comprehensively quantify endocannabinoids, ERCs, and their phospholipid precursors in clinical studies.

Science and art of protein formulation development

Protein pharmaceuticals have become a significant class of marketed drug products and are expected to grow steadily over the next decade. Development of a commercial protein product is, however, a rather complex process. A critical step in this process is formulation development, enabling the final product configuration. A number of challenges still exist in the formulation development process. This review is intended to discuss these challenges, to illustrate the basic formulation development processes, and to compare the options and strategies in practical formulation development.

Controlling Phase Separation of Lysozyme with Polyvalent Anions

The ability of polyvalent anions to influence protein-protein interactions and protein net charge was investigated through solubility and turbidity experiments, determination of osmotic second virial coefficients ( B₂₂), and ζ-potential values for lysozyme solutions.  B₂₂ values showed that all anions reduce protein-protein repulsion between positively charged lysozyme molecules, and those anions with higher net valencies are more effective. The polyvalent anions pyrophosphate and tripolyphosphate were observed to induce protein reentrant condensation, which has been previously observed with negatively charged proteins in the presence of trivalent cations. Reentrant condensation is a phenomenon in which low concentrations of polyvalent ions induce protein precipitation, but further increasing polyvalent ion concentration causes the protein precipitate to resolubilize. Interestingly, citrate does not induce lysozyme reentrant condensation despite having a similar charge, size, and shape to pyrophosphate. We observe qualitative differences in protein behavior when compared against negatively charged proteins in solutions of trivalent cations. The polyphosphate ions induce a much stronger protein-protein attraction, which correlates with the occurrence of a liquid-gel transition that replaces the liquid-liquid transition observed with trivalent cations. The results indicate that solutions of polyphosphate ions provide a model system for exploring the link between the protein-phase diagram and model interaction potentials and also highlight the importance that ion-specific effects can have on protein solubility.

The effect of Hofmeister anions on water structure at protein surfaces

To understand the effects of specific ions on protein-water interactions and the thermodynamic stability of proteins in salt solutions, we use a molecular dynamics (MD) simulation to examine the water structure, orientational distribution, and dynamics near the surface of ubiquitin. In particular, we consider NaCl, NaBF₄, NaSCN, and NaClO₄ salt solutions containing ubiquitin, where the anions of the latter three salts are well-known chaotropic ions in the Hofmiester anion series. The number of hydrogen bonds (H-bonds) per water molecule is found to decrease significantly at the ubiquitin-water interface, indicating a significant disruption of the water H-bonding network. The distribution of the water H-bond numbers near the protein surface is modulated by dissolved ions, and the extent of the ion effect on the H-bonding network structure follows the order of the Hofmeister anion series, while there are no specific ion effects on water properties at distances larger than 5 Å from the protein surface. From detailed analyses of the surface area, volume, and root-mean-square deviation (RMSD) of ubiquitin, we show that changes in the properties of the protein could originate from the disruption of the water H-bond network induced by ions with a higher affinity for the protein surface instead of direct protein residue-ion interactions. An interesting observation made here is that the orientational distribution of water molecules at the protein-water interface is close to random, but there is a slight preference for interfacial water molecules with a straddle structure within 2.5 Å of the protein surface, where one of the two OH groups points away from the protein surface and the other points toward the surface. In addition, comparing the MD simulation results for ubiquitin solutions with dissolved NaSCN and KSCN, we show that Na⁺ affects the water H-bonding structure at the protein surface more than K⁺. It is clear that the H-bonding network structure of water more than one water layer away from the protein surface is not distinguishably different from that of neat water. We thus anticipate that the present work will provide insights into the scale of specific ion effects on the H-bonding structure and orientational distribution of water in the vicinity of protein surfaces in aqueous solutions.

NMR spectroscopic studies of a TAT-derived model peptide in imidazolium-based ILs: influence on chemical shifts and the cis/trans equilibrium state

The impact of ionic liquids on the chemical shifts and the cis/trans equilibrium state of a model peptide was systematically investigated by NMR spectroscopy.