Contribution of the surface free energy perturbation to protein-solvent interactions

On the role of excipients in biopharmaceuticals manufacture: Modelling-guided formulation identifies the protective effect of arginine hydrochloride excipient on spray-dried Olipudase alfa recombinant protein

Biopharmaceuticals are labile biomolecules that must be safeguarded to ensure the safety, quality, and efficacy of the product. Batch freeze-drying is an established means of manufacturing solid biopharmaceuticals but alternative technologies such as spray-drying may be more suitable for continuous manufacturing of inhalable biopharmaceuticals. Here we assessed the feasibility of spray-drying Olipudase alfa, a novel parenteral therapeutic enzyme, by evaluating some of its critical quality attributes (CQAs) in a range of excipients, namely, trehalose, arginine (Arg), and arginine hydrochloride (Arg-HCl) in the sucrose/methionine base formulation. The Arg-HCl excipient produced the best gain in CQAs of spray-dried Olipudase with a 63% reduction in reconstitution time and 83% reduction in the optical density of the solution. Molecular dynamics simulations revealed the atomic-scale mechanism of the protein-excipient interactions, substantiating the experimental results. The Arg-HCl effect was explained by the calculated thermal stability and structural order of the protein wherein Arg-HCl acted as a crowding agent to suppress protein aggregation and promote stabilization of Olipudase post-spray-drying. Therefore, by rational selection of appropriate excipients, our experimental and modelling dataset confirms spray-drying is a promising technology for the manufacture of Olipudase and demonstrates the potential to accelerate development of continuous manufacturing of parenteral biopharmaceuticals.

The contrasting roles of co-solvents in protein formulations and food products

Protein aggregation is a major hurdle in developing biopharmaceuticals, in particular protein formulation area, but plays a pivotal role in food products. Co-solvents are used to suppress protein aggregation in pharmaceutical proteins. On the contrary, aggregation is encouraged in the process of food product making. Thus, it is expected that co-solvents play a contrasting role in biopharmaceutical formulation and food products. Here, we show several examples that utilize co-solvents, e.g., salting-out salts, sugars, polyols and divalent cations in promoting protein-protein interactions. The mechanisms of co-solvent effects on protein aggregation and solubility have been studied on aqueous protein solution and applied to develop pharmaceutical formulation based on the acquired scientific knowledge. On the contrary, co-solvents have been used in food industries based on empirical basis. Here, we will review the mechanisms of co-solvent effects on protein-protein interactions that can be applied to both pharmaceutical and food industries and hope to convey knowledge acquired through research on co-solvent interactions in aqueous protein solution and formulation to those involved in food science and provide those involved in protein solution research with the observations on aggregation behavior of food proteins.

Factor defining the effects of tetraalkylammonium chloride on stability, folding, and dynamics of horse cytochrome c

This article describes the molecular mechanism by which tetraalkylammonium chloride ([R4N]Cl: R- = methyl (Me), ethyl (Et), propyl (Pr),butyl (Bu)) modulates the stability, folding, and dynamics of cytochrome c (Cyt c). Analysis of [R4N]Cl effects on thermal/chemical denaturations, millisecond refolding/unfolding kinetics, and slow CO-association kinetics of Cyt c without and with denaturant providing some significant results: (i) [R4N]Cl decreasing the unfolding free energy estimated by thermodynamic and kinetic analysis of thermal/chemical denaturation curves and kinetic chevrons (Log kobs-[GdmCl]) of Cyt c, respectively (ii) hydrophobicity of R-group of [R4N]Cl, preferential inclusion of [R4N]Cl at the protein surface, and destabilizing enthalpic attractive interactions of [Me4N]Cl and steric entropic interactions of [Et4N]Cl,[Pr4N]Cl and [Bu4N]Cl with protein contribute to [R4N]Cl-induced decrease thermodynamic stability of Cyt c (iii) [R4N]Cl exhibits an additive effect with denaturant to decrease thermodynamic stability and refolding rates of Cyt c (iv) low concentrations of [R4N]Cl (≤ 0.5 M) constrain the motional dynamics while the higher concentrations (>0.75 M [R4N]Cl) enhance the structural-fluctuations that denture protein (v) hydrophobicity of R-group of [R4N]Cl alters the [denaturant]-dependent conformational stability, refolding-unfolding kinetics, and CO-association kinetics of Cyt c. Furthermore, the MD simulations depicted that the addition of 1.0 M of [R4N]Cl increased the conformational fluctuations in Cyt c leading to decreased structural stability in the order [Me4N]Cl < [Et4N]Cl < [Pr4N]Cl < [Bu4N]Cl consistent with the experimental results.

Systematic enhancement of protein crystallization efficiency by bulk lysine-to-arginine (KR) substitution

Structural genomics consortia established that protein crystallization is the primary obstacle to structure determination using x-ray crystallography. We previously demonstrated that crystallization propensity is systematically related to primary sequence, and we subsequently performed computational analyses showing that arginine is the most overrepresented amino acid in crystal-packing interfaces in the Protein Data Bank. Given the similar physicochemical characteristics of arginine and lysine, we hypothesized that multiple lysine-to-arginine (KR) substitutions should improve crystallization. To test this hypothesis, we developed software that ranks lysine sites in a target protein based on the redundancy-corrected KR substitution frequency in homologs. This software can be run interactively on the worldwide web at https://www.pxengineering.org/. We demonstrate that three unrelated single-domain proteins can tolerate 5-11 KR substitutions with at most minor destabilization, and, for two of these three proteins, the construct with the largest number of KR substitutions exhibits significantly enhanced crystallization propensity. This approach rapidly produced a 1.9 Å crystal structure of a human protein domain refractory to crystallization with its native sequence. Structures from Bulk KR-substituted domains show the engineered arginine residues frequently make hydrogen-bonds across crystal-packing interfaces. We thus demonstrate that Bulk KR substitution represents a rational and efficient method for probabilistic engineering of protein surface properties to improve crystallization.

Spectroscopic Investigation of Phosphorus Mineralization as Affected by the Calcite-Water Interfacial Chemistry

The mineralization and bioavailability of phytic acid, the predominant organic phosphorus (OP) species in many soils, have generally been rendered limited due to its interaction with soil minerals. In particularly calcareous and neutral to slightly alkaline soils, phytic acid is known to actively react with calcite, although how this interaction affects phytic acid mineralization is still unknown. This study, therefore, investigated the mechanisms regarding how the calcite-water interface influences phytic acid mineralization by phytase, at pHs 6 and 8 using in situ spectroscopic techniques including solution nuclear magnetic resonance and attenuated total reflection Fourier transform infrared spectroscopy. The findings indicated a pH-specific effect of the calcite-water interface. Inhibited phytase activity and thus impaired phytic acid mineralization were induced by calcite at pH 6, while the opposite effect was observed at pH 8. How the interaction between phytic acid and calcite and between phytase and calcite differed between the two pH values contributed to the pH-specific effect. The results demonstrate the importance of soil pH, enzyme-, and OP-clay mineral interactions in controlling the mineralization and transformation of OP and, consequently, the release of phosphate in soils. The findings can also provide implications for the management of calcite-rich and limed soils.

Mechanistic Insight into Poly-Reactivity of Immune Antibodies upon Acid Denaturation or Arginine Mutation in Antigen-Binding Regions

The poly-reactivity of antibodies is defined as their binding to specific antigens as well as to related proteins and also to unrelated targets. Poly-reactivity can occur in individual molecules of natural serum antibodies, likely due to their conformation flexibility, and, for therapeutic antibodies, it plays a critical role in their clinical development. On the one hand, it can enhance their binding to target antigens and cognate receptors, but, on the other hand, it may lead to a loss of antibody function by binding to off-target proteins. Notably, poly-reactivity has been observed in antibodies subjected to treatments with dissociating, destabilizing or denaturing agents, in particular acidic pH, a common step in the therapeutic antibody production process involving the elution of Protein-A bound antibodies and viral clearance using low pH buffers. Additionally, poly-reactivity can emerge during the affinity maturation in the immune system, such as the germinal center. This review delves into the underlying potential causes of poly-reactivity, highlighting the importance of conformational flexibility, which can be further augmented by the acid denaturation of antibodies and the introduction of arginine mutations into the complementary regions of antibody-variable domains. The focus is placed on a particular antibody's acid conformation, meticulously characterized through circular dichroism, differential scanning calorimetry, and sedimentation velocity analyses. By gaining a deeper understanding of these mechanisms, we aim to shed light on the complexities of antibody poly-reactivity and its implications for therapeutic applications.

Elucidating the mechanisms of additive effects at high concentrations on hydrophobic interaction chromatography

Hydrophobic interaction chromatography (HIC) is a commonly used chromatography technique for purifying proteins. It utilizes salting-out salts to facilitate the binding of native proteins to weakly hydrophobic ligands. There have been three proposed mechanisms for the promoting effects of salting-out salts, which include the dehydration of proteins by salts, cavity theory, and salt exclusion. To evaluate the above three mechanisms, an HIC study was conducted on Phenyl Sepharose using four different additives. These additives included a salting-out salt (NH₄)₂SO₄, sodium phosphate that increases the surface tension of water, a salting-in salt MgCl₂, and an amphiphilic protein-precipitant polyethylene glycol (PEG). Results indicated that the first two salts resulted in protein binding, while MgCl₂ and PEG led to flow-through. These findings were then used to interpret the three proposed mechanisms, which showed that MgCl₂ and PEG deviated from the dehydration mechanism, and MgCl₂ also deviated from the cavity theory. The observed effects of these additives on HIC were reasonably well explained for the first time by their interactions with proteins.

Agarose native gel electrophoresis analysis of thermal aggregation controlled by Hofmeister series

The effects of salting-in and salting-out salts defined by Hofmeister series on the solution state of bovine serum albumin (BSA) in 50 mM Tris-HCl buffer at pH 7.4 before and after thermal unfolding at 80 °C for 5 min were examined using agarose native gel electrophoresis and mass photometry. Gel electrophoresis showed that salting-in MgCl₂, CaCl₂ and NaSCN resulted in formation of intermediate structures of BSA upon heating on native gel, while heating in buffer alone resulted in aggregated bands. Mass photometry showed large loss of monomer and oligomers when heated in this buffer, but retaining these structures in the presence of 1 M MgCl₂ and NaSCN. To our surprise, salting-out MgSO₄ also showed a similar effect on gel electrophoresis and mass photometry. Salting-out NaCl and (NH₄)₂SO₄ resulted in smearing and aggregated bands, which were supported by mass photometry. Aggregation-suppressive ArgHCl also showed oligomer aggregates upon gel electrophoresis and mass photometry.

Non-Affinity Purification of Antibodies

Currently, purification of antibodies is mainly carried out using a platform technology composed primarily of Protein A chromatography as a capture step, regardless of the scale. However, Protein A chromatography has a number of drawbacks, which are summarized in this review. As an alternative, we propose a simple small-scale purification protocol without Protein A that uses novel agarose native gel electrophoresis and protein extraction. For large-scale antibody purification, we suggest mixed-mode chromatography that can in part mimic the properties of Protein A resin, focusing on 4-Mercapto-ethyl-pyridine (MEP) column chromatography.

An efficient method to predict protein thermostability in alanine mutation

The relationship between protein sequence and its thermodynamic stability is a critical aspect of computational protein design. In this work, we present a new theoretical method to calculate the free energy change (ΔΔG) resulting from a single-point amino acid mutation to alanine in a protein sequence. The method is derived based on physical interactions and is very efficient in estimating the free energy changes caused by a series of alanine mutations from just a single molecular dynamics (MD) trajectory. Numerical calculations are carried out on a total of 547 alanine mutations in 19 diverse proteins whose experimental results are available. The comparison between the experimental ΔΔG_exp and the calculated values shows a generally good correlation with a correlation coefficient of 0.67. Both the advantages and limitations of this method are discussed. This method provides an efficient and valuable tool for protein design and engineering.

Trehalose and its applications in the food industry

Trehalose is a nonreducing disaccharide composed of two glucose molecules linked by α, α-1,1-glycosidic bond. It is present in a wide variety of organisms, including bacteria, fungi, insects, plants, and invertebrate animals. Trehalose has distinct physical and chemical properties that have been investigated for their biological importance in a range of prokaryotic and eukaryotic species. Emerging research on trehalose has identified untapped opportunities for its application in the food, medical, pharmaceutical, and cosmetics industries. This review summarizes the chemical and biological properties of trehalose, its occurrence and metabolism in living organisms, its protective role in molecule stabilization, and natural and commercial production methods. Utilization of trehalose in the food industry, in particular how it stabilizes protein, fat, carbohydrate, and volatile compounds, is also discussed in depth. Challenges and opportunities of its application in specific applications (e.g., diagnostics, bioprocessing, ingredient technology) are described. We conclude with a discussion on the potential of leveraging the unique molecular properties of trehalose in molecular stabilization for improving the safety, quality, and sustainability of our food systems.

The protein-stabilizing effects of TMAO in aqueous and non-aqueous conditions

We present a new water-dependent molecular mechanism for the widely-used protein stabilizing osmolyte, trimethylamine N-oxide (TMAO), whose mode of action has remained controversial. Classical interpretations, such as osmolyte exclusion from the vicinity of protein, cannot adequately explain the behavior of this osmolyte and were challenged by recent data showing the direct interactions of TMAO with proteins, mainly via hydrophobic binding. Solvent effect theories also fail to propose a straightforward mechanism. To explore the role of water and the hydrophobic association, we disabled osmolyte-protein hydrophobic interactions by replacing water with hexane and using lipase enzyme as an anhydrous-stable protein. Biocatalysis experiments showed that under this non-aqueous condition, TMAO does not act as a stabilizer, but strongly deactivates the enzyme. Molecular dynamics (MD) simulations reveal that TMAO accumulates near the enzyme and makes many hydrogen bonds with it, like denaturing osmolytes. Some TMAO molecules even reach the active site and interact strongly with the catalystic traid. In aqueous solvent, the enzyme functions well: the extent of TMAO interactions is reduced and can be divided into both polar and non-polar terms. Structural analysis shows that in water, some TMAO molecules bind to the enzyme surface like a surfactant. We show that these interactions limit water-protein hydrogen bonds and unfavorable water-hydrophobic surface contacts. Moreover, a more hydrophobic environment is formed in the solvation layer, which reduces water dynamics and subsequently, rigidifies the backbone in aqueous solution. We show that osmolyte amphiphilicity and protein surface heterogeneity can address the weaknesses of exclusion and solvent effect theories about the TMAO mechanism.

A simple theory for interfacial properties of dilute solutions

Recent studies suggest that cosolute mixtures may exert significant non-additive effects upon protein stability. The corresponding liquid–vapor interfaces may provide useful insight into these non-additive effects. Accordingly, in this work, we relate the interfacial properties of dilute multicomponent solutions to the interactions between solutes. We first derive a simple model for the surface excess of solutes in terms of thermodynamic observables. We then develop a lattice-based statistical mechanical perturbation theory to derive these observables from microscopic interactions. Rather than adopting a random mixing approximation, this dilute solution theory (DST) exactly treats solute–solute interactions to lowest order in perturbation theory. Although it cannot treat concentrated solutions, Monte Carlo (MC) simulations demonstrate that DST describes the interactions in dilute solutions with much greater accuracy than regular solution theory. Importantly, DST emphasizes a fundamental distinction between the “intrinsic” and “effective” preferences of solutes for interfaces. DST predicts that three classes of solutes can be distinguished by their intrinsic preference for interfaces. While the surface preference of strong depletants is relatively insensitive to interactions, the surface preference of strong surfactants can be modulated by interactions at the interface. Moreover, DST predicts that the surface preference of weak depletants and weak surfactants can be qualitatively inverted by interactions in the bulk. We also demonstrate that DST can be extended to treat surface polarization effects and to model experimental data. MC simulations validate the accuracy of DST predictions for lattice systems that correspond to molar concentrations.

Exploring the structural basis of conformational alterations of myoglobin in the presence of spermine through computational modeling, molecular dynamics simulations, and spectroscopy methods

Spermine as polyamines can have interaction with the myoglobin (Mb). The intent of this pondering to evaluate the impact of spermine on Mb properties, for example, the structure and thermal stability. For this analysis, the following approaches are employed. Thermodynamics, molecular dynamics (MD), and docking and the use of other spectroscopic procedures. The results of fluorescence spectroscopy and docking showed that binding spermine to Mb was spontaneous. Spermine quenched the fluorescence of Mb through the static quenching process. The thermal stability of Mb was incremented when the concentration of spermine increased. The CD spectra showed Mb's secondary structure shift with a rise in β-sheet and a decrease in α-helicity Mb's in spermine presence. Molecular docking and MD simulation outcomes demonstrate that electrostatic forces show a critical function in stabilizing of this complex, which is in conforming to spectroscopic results.Communicated by Ramaswamy H. Sarma.

Analysis of proteins by agarose native gel electrophoresis in the presence of solvent additives

Solvent additives, including NaCl, arginine hydrochloride (ArgHCl), glycine and sucrose, are used to enhance protein stability or reduce protein aggregation. Here, we studied the effects of these additives on proteins using agarose native gel electrophoresis. Since these additives are used at relatively high concentration, we first confirmed that they do not interfere with the performance of the native gel electrophoresis. Agarose native gel electrophoresis showed that aggregation of bovine serum albumin (BSA) induced by heating was slightly reduced by NaCl and ArgHCl. On the contrary, glycine and sucrose had marginal effects. ArgHCl and NaCl promoted heat aggregation of monoclonal antibody (mAb), while glycine and sucrose stabilized the native mAb. Arginine methyl ester inhibited heat aggregation of lysozyme and, to a much lesser extent, BSA. These results show that agarose native gel electrophoresis can be used to analyze the effects of solvent additives on proteins subjected to heat stresses. SYPRO Orange that stains only unfolded proteins confirmed unfolded structures of soluble aggregates.

Small-to-large length scale transition of TMAO interaction with hydrophobic solutes

We report the effect of trimethylamine N-oxide (TMAO) on the solvation of nonpolar solutes in water studied with molecular dynamics (MD) simulations and free-energy calculations. The simulation data indicate the occurrence of a length scale crossover in the TMAO interaction with repulsive Weeks-Chandler-Andersen (WCA) solutes: while TMAO is depleted from the hydration shell of a small WCA solute (methane) and increases the free-energy cost of solute-cavity formation, it preferentially binds to a large WCA solute (α-helical polyalanine), reducing the free-energy cost of solute-cavity formation via a surfactant-like mechanism. Significantly, we show that this surfactant-like behaviour of TMAO reinforces the solvent-mediated attraction between large WCA solutes by means of an entropic force linked to the interfacial accumulation of TMAO. Specifically, this entropic force arises from the natural tendency of adsorbed TMAO molecules to mix back into the bulk. It therefore favours solute-solute contact states that minimise the surface area exposed to the solvent and have a small overall number of TMAO molecules adsorbed. In contrast to the well-known depletion force, its effect is compensated by enthalpic solute-solvent interactions. Correspondingly, the hydrophobic association free energy of the large α-helical solutes passes through a minimum at low TMAO concentration when cohesive solute-solvent van der Waals interactions are considered. The observations reported herein are reminiscent to cosolvent effects on hydrophobic polymer coil-globule collapse free energies (Bharadwaj et al., Commun. Chem. 2020, 3, 165) and may be of general significance in systems whose properties are determined by hydrophobic self-assembly.

Classification of protein solubilizing additives by fluorescence assay

Aromatic interaction plays a crucial role in controlling protein interaction by additives. Here we investigated the interaction of protein salting-in (solubilizing) additives with tryptophan (Trp), tyrosine (Tyr), indole, and proteins based on their fluorescence spectra. Five salting-in additives, i.e., arginine (Arg), urea, guanidine (Gdn), ethylene glycol (EG), and magnesium chloride (MgCl₂), showed different effects on the fluorescence properties of Trp and Tyr. Arg significantly reduced fluorescence intensity of Trp and Tyr, as was the case for glycine to a lesser extent. MgCl₂ and calcium chloride (CaCl₂) showed little effect on the aromatic fluorescence spectra. Gdn also showed little effect on the aromatic fluorescence spectra of Trp and Tyr even at high concentrations. EG increased the aromatic fluorescence intensity of Trp and Tyr with blue-shifted emission wavelength. Urea enhanced fluorescence of Trp and Tyr without altering emission wavelength. These results indicate that the protein solubilizing additives interact with aromatic groups differently.

Insight into the protein salting-in mechanism of arginine, magnesium chloride and ethylene glycol: Solvent interaction with aromatic solutes

Key factors in the salting-in effects on proteins of additives are their interactions with aromatic groups. We studied the interaction of four aromatic solutes, benzyl alcohol (BA), phenol, 4-hydroxybenzyl alcohol (4-HBA) and methyl gallate (MG), with different salting-in additives, arginine hydrochloride (ArgHCl), magnesium chloride (MgCl₂), ethylene glycol (EG), and guanidine hydrochloride (GdnHCl) using solubility measurements. We used sodium chloride (NaCl) as a control. MgCl₂ decreased the solubility of the four aromatic solutes with weak solute dependence. In contrast, ArgHCl, GdnHCl, and EG increased the solubility of four aromatic solutes with a similar solute dependence. Their salting-in effects were weaker on BA and 4-HBA and stronger on phenol and MG. These results indicate that attached groups alter the aromatic properties, affecting the interactions between the benzene ring and these three additives. More importantly, the observed results demonstrate that the salting-in mechanism is different between MgCl₂, EG and ArgHCl, which should play a role in their effects on protein solubility.

A tribute to Dr. Serge N. Timasheff, our mentor

Dr. Serge N. Timasheff, our mentor and friend, passed away in 2019. This article is a collection of tributes from his postdoctoral fellows, friends, and daughter, who all have been associated with or influenced by him or his research. Dr. Timasheff is a pioneer of research on thermodynamic linkage between ligand interaction and macromolecular reaction. We all learned a great deal from Dr. Timasheff, not only about science but also about life.

Analysis of the binding modes of the first- and second-generation antiandrogens with respect to F876L mutation

Androgen receptor (AR) is an important target for the treatment of prostate cancer, and mutations in the AR have an important impact on the resistance of existing drugs. In this work, we performed molecular dynamics simulations of the existing marketed antiandrogens flutamide, nilutamide, bicalutamide, enzalutamide, apalutamide, darolutamide, and its main metabolite ORM15341 in complex with the wild-type and F876L mutant AR. We calculated the residue-specific binding free energy contribution of the wild-type and mutant ARs with the AS-IE method and analyzed the hotspot residues and the binding free energy contributions of specific residues before and after the mutation. In addition, we analyzed the total binding obtained by adding residue binding energy contributions and compared the results with experimental values. The obtained residue-specific binding information should be very helpful in understanding the mechanism of drug resistance with respect to specific mutations and in the design of new generation drugs against possible new mutations.

Osmoprotectant Coated Thermostable Gold Nanoparticles Efficiently Restrict Temperature-Induced Amyloid Aggregation of Insulin

Naturally occurring osmoprotectants are known to prevent aggregation of proteins under various stress factors including extreme pH and elevated temperature conditions. Here, we synthesized gold nanoparticles coated with selected osmolytes (proline, hydroxyproline, and glycine) and examined their effect on temperature-induced amyloid-formation of insulin hormone. These uniform, thermostable, and hemocompatible gold nanoparticles were capable of inhibiting both spontaneous and seed-induced amyloid aggregation of insulin. Both quenching and docking experiments suggest a direct interaction between the osmoprotectant-coated nanoparticles and aggregation-prone hydrophobic stretches of insulin. Circular-dichroism results confirmed the retention of insulin's native structure in the presence of these nanoparticles. Unlike the indirect solvent-mediated effect of free osmolytes, the inhibition effect of osmolyte-coated gold nanoparticles was observed to be mediated through their direct interaction with insulin. The results signify the protection of the exposed aggregation-prone domains of insulin from temperature-induced self-assembly through osmoprotectant-coated nanoparticles, and such effect may inspire the development of osmolyte-based antiamyloid nanoformulations.

Hydrophobic interactions control the self-assembly of DNA and cellulose

Desoxyribosenucleic acid, DNA, and cellulose molecules self-assemble in aqueous systems. This aggregation is the basis of the important functions of these biological macromolecules. Both DNA and cellulose have significant polar and nonpolar parts and there is a delicate balance between hydrophilic and hydrophobic interactions. The hydrophilic interactions related to net charges have been thoroughly studied and are well understood. On the other hand, the detailed roles of hydrogen bonding and hydrophobic interactions have remained controversial. It is found that the contributions of hydrophobic interactions in driving important processes, like the double-helix formation of DNA and the aqueous dissolution of cellulose, are dominating whereas the net contribution from hydrogen bonding is small. In reviewing the roles of different interactions for DNA and cellulose it is useful to compare with the self-assembly features of surfactants, the simplest case of amphiphilic molecules. Pertinent information on the amphiphilic character of cellulose and DNA can be obtained from the association with surfactants, as well as on modifying the hydrophobic interactions by additives.

Computational analysis of binding free energies, hotspots and the binding mechanism of Bcl-xL/Bcl-2 binding to Bad/Bax

Hierarchical clustering tree of residues providing contributions to system binding based on the binding free energy of specific residues for (A) Bcl-xL systems (B) Bcl-2 systems.

Enhancing the crystallisation of insulin using amino acids as soft-templates to control nucleation

Amino acid as soft templates in promoting nucleation of insulin.

Production of Lysozyme-PLGA-Loaded Microparticles for Controlled Release Using Hot-Melt Extrusion

Biopharmaceuticals are usually administered intravenously with frequent dosing regimens which may decrease patient compliance. Controlled-release formulations allow to reduce the frequency of injections while providing a constant dosing of the biopharmaceutical over extended periods. These formulations are typically produced by emulsions, requiring high amounts of organic solvents and have limited productivity. Hot-melt extrusion (HME) is an alternative technology to produce controlled drug delivery systems. It is a continuous solvent-free process, leading to a small ecological footprint and higher productivity. However, it may induce thermolabile compounds' degradation. In this work, the impact of the formulation and extrusion temperature on lysozyme's bioactivity and release profile of poly(lactic-co-glycolic acid) (PLGA)-based extended release formulations were evaluated using a design-of-experiments (DoE) approach. The lysozyme-loaded PLGA microparticles were produced by HME followed by milling. It was observed that the in vitro release (IVR) profile was mainly affected by the drug load; higher drug load led to higher burst and total lysozyme release after 14 days. HME temperature seemed to decrease lysozyme's activity although this correlation was not statistically significant (p value = 0.0490). Adding polyethylene glycol 400 (PEG 400) as a plasticizer to the formulation had no significant impact on the lysozyme release profile. The burst release was effectively mitigated with the inclusion of a washing step. Washing the microparticles with water reduced the burst release by 80% whereas washing them with a poly(vinyl alcohol) (PVA) aqueous solution eliminated it. In conclusion, HME is demonstrated to be suitable in producing controlled-release microparticles of small biopharmaceuticals. Graphical abstract.

A method for efficient calculation of thermal stability of proteins upon point mutations

A method for efficient prediction of the relative stability of a protein due to a single amino acid point mutation is presented. In this approach, we calculate the free energy change due to an arbitrary point mutation of a protein from a single MD trajectory of the wild type protein. The method is tested on 27 diverse protein systems with a total of 853 mutations and the calculated relative free energies show a generally good correlation with the experimental values (a correlation coefficient of 0.63). Comparison with the free energy perturbation (FEP) method and the recently developed machine learning methods on two different benchmark data sets shows that the current method is computationally efficient and also numerically reliable for predicting the changes in thermostability upon an arbitrary point mutation of a protein. A discussion is provided on how to further improve the accuracy of the method for the prediction of thermostability of proteins.

An accurate free energy estimator: based on MM/PBSA combined with interaction entropy for protein-ligand binding affinity

The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method is constantly used to calculate the binding free energy of protein-ligand complexes, and has been shown to effectively balance computational cost against accuracy. The relative binding affinities obtained by the MM/PBSA approach are acceptable, while it usually overestimates the absolute binding free energy. This paper proposes four free energy estimators based on the MM/PBSA for enthalpy change combined with interaction entropy (IE) for entropy change using different weights for individual energy terms. The ΔG_{PBSA_IE} method is determined to be an optimal estimator based on its performance in terms of the correlation between experimental and theoretical values and error estimations. This approach is optimized using high-quality experimental values from a training set containing 84 protein-ligand systems, and the coefficients for the sum of electrostatic energy and polar solvation free energy, van der Waals (vdW) energy, non-polar solvation energy and entropy change are obtained by multivariate linear fitting to the corresponding experimental values. A comparison between the traditional MM/PBSA method and this method shows that the correlation coefficient is improved from 0.46 to 0.72 and the slope of the regression line increases from 0.10 to 1.00. More importantly, the mean absolute error (MAE) is significantly reduced from 22.52 to 1.59 kcal mol^-1. Furthermore, the numerical stability of this method is validated on a test set with a similar correlation coefficient, slope and MAE to those of the training set. Based on the above advantages, the ΔG_{PBSA_IE} method can be a powerful tool for a reliable and accurate estimation of binding free energy and plays a significant role in a detailed energetic investigation of protein-ligand interaction.

Molecular Mechanism of Selective Binding of NMS-P118 to PARP-1 and PARP-2: A Computational Perspective

The nuclear protein poly (ADP-ribose) polymerase-1 (PARP-1) inhibitors have been proven effective to potentiate both chemotherapeutic agents and radiotherapy. However, a major problem of most current PARP inhibitors is their lack of selectivity for PARP-1 and its closest isoform PARP-2. NMS-P118 is a highly selective PARP inhibitor that binds PARP-1 stronger than PARP-2 and has many advantages such as excellent pharmacokinetic profiles. In this study, molecular dynamics (MD) simulations of NMS-P118 in complex with PARP-1 and PARP-2 were performed to understand the molecular mechanism of its selectivity. Alanine scanning together with free energy calculation using MM/GBSA and interaction entropy reveal key residues that are responsible for the selectivity. Although the conformation of the binding pockets and NMS-P118 are very similar in PARP-1 and PARP-2, most of the hot-spot residues in PARP-1 have stronger binding free energy than the corresponding residues in PARP-2. Detailed analysis of the binding energy shows that the 4′4-difluorocyclohexyl ring on NMS-P118 form favorable hydrophobic interaction with Y889 in PARP-1. In addition, the H862 residue in PARP-1 has stronger binding free energy than H428 in PARP-2, which is due to shorter distance and stronger hydrogen bonds. Moreover, the negatively charged E763 residue in PARP-1 forms stronger electrostatic interaction energy with the positively charged NMS-P118 than the Q332 residue in PARP-2. These results rationalize the selectivity of NMS-P118 and may be useful for designing novel selective PARP inhibitors.

Improving the performance of the MM/PBSA and MM/GBSA methods in recognizing the native structure of the Bcl-2 family using the interaction entropy method

Correct discrimination of native structure plays an important role in drug design. IE method significantly improves the performance of MM/PB(GB)SA method in discriminating native and decoy structures in protein–ligand/protein systems of Bcl-2 family.

Understanding conformational diversity of heat shock protein 90 (HSP90) and binding features of inhibitors to HSP90 via molecular dynamics simulations

Heat shock protein 90 (HSP90) is a promising target for treatment of cancer, and inhibitor bindings can generate efficient suppression on tumor in multiple ways. In this work, 140-ns molecular dynamics simulations were performed on six systems. Principal component analysis was subsequently carried out to explore the conformational diversity of HSP90. The results suggest that inhibitor bindings induce large conformational changes of HSP90, which tends to enlarge the volume of the binding pocket to facilitate the entrance of inhibitors. Hierarchical clustering analyses, the calculation of the energy contribution of each atom, and the analyses of hydrogen-bonding interactions were performed. The results indicate that 20 residues in group A of the hierarchical tree are responsible for major contributions, and van der Waals interactions as well as hydrogen-bonding interactions between important residues in HSP90 and key regions of inhibitors are the main force for promoting inhibitor bindings. We expect that this work can provide useful theoretical information for development of efficient inhibitors targeting HSP90.

Hofmeister Order of Anions on Protein Stability Originates from Lifshitz-van der Waals Dispersion Interaction with the Protein Phase

The mechanism underpinning the Hofmeister order of anions on protein stability and other physical and biological processes has been a mystery since its discovery in 1888. In the present study, we investigated electrostatic and Lifshitz-van der Waals (L-vdW) dispersion (electrodynamic) interactions between Hofmeister salts and four monomeric globular proteins. It is shown that structure-stabilizing salts exerted positive L-vdW pressure, whereas structure-destabilizing salts exerted negative L-vdW pressure on proteins. The relative order of the L-vdW pressure followed the Hofmeister series and it overshadowed the electrostatic pressure at high salt concentrations. The net change in the thermal denaturation temperature (ΔT_d) of proteins in 0.8 M Hofmeister salt solutions followed a linear relationship (r² > 0.8) with the net electrodynamic pressure regardless of the physicochemical differences between proteins. This study also revealed that segregation of anions into structure stabilizers and destabilizers depended on the dielectric susceptibility of the anion in the ultraviolet region: ions having absorbtion spectrum in the ultraviolet region (e.g., Cl^-, Br^-, I^-, and SCN^-) exerted a negative electrodynamic pressure, whereas those with absorbtion spectrum only in the infrared region, for example, SO₄^2-, exerted a positive electrodynamic pressure. The lack of ultraviolet absorption of SO₄^2- ions was because of quenching of ultraviolet radiation by water at below 170 nm.

Protein aggregation suppressor arginine as an effective mouth cleaning agent

We have tested here whether or not arginine, a well-known aggregation suppressor, is effective in removing bacterial cells, which may present a potential risk of accidental pneumonia infection in aged individuals, from the oral mucosal membranes. This is based on the ability of arginine to suppress protein-protein interaction and surface adsorption and increase the solubility of organic compounds. Twelve student volunteers were subjected to mouthwashes with saline, citrate buffer (pH 3.5), arginine (pH 3.5) and a commercial Listerine. Insignificant effects were observed with saline and citrate buffer, whereas arginine and Listerine mouthwashes led to significant reduction of bacterial cells from the dorsal side of the volunteer's tongue. Arginine also appeared to disrupt biofilms present in the mouth.

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.

Effect of double mutations T790M/L858R on conformation and drug-resistant mechanism of epidermal growth factor receptor explored by molecular dynamics simulations

Epidermal growth factor receptor (EGFR) is one of the most promising targets for the treatment of cancers. Double mutations T790M/L858R lead to different degrees of drug resistance toward inhibitors. In this study, molecular dynamics (MD) simulations followed by principal component analysis are performed to study the conformational changes of EGFR induced by T790M/L858R. The results suggest that T790M/L858R cause obvious disturbance of the structural stability of EGFR. Molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) and residue-based free energy decomposition methods are integrated to explore the drug-resistant mechanism of T790M/L858R toward inhibitors. The results show that the decrease in van der Waals interactions of inhibitors with the mutated EFGR relative to the wild-type (WT) one is the main force inducing drug resistance of T790M/L858R toward inhibitors TAK-285, while drug resistance toward W2P and HKI-272 is dominated by the decrease in van der Waals interactions and the increase in polar interactions. We expect that the information obtained from this study can aid rational design of effective drugs to relieve drug resistance of EGFR induced by T790M/L858R.

The impact of interior dielectric constant and entropic change on HIV-1 complex binding free energy prediction

At present, the calculated binding free energy obtained using the molecular mechanics/Poisson-Boltzmann (Generalized-Born) surface area (MM/PB(GB)SA) method is overestimated due to the lack of knowledge of suitable interior dielectric constants in the simulation on the interaction of Human Immunodeficiency Virus (HIV-1) protease systems with inhibitors. Therefore, the impact of different values of the interior dielectric constant and the entropic contribution when using the MM/PB(GB)SA method to calculate the binding free energy was systemically evaluated. Our results show that the use of higher interior dielectric constants (1.4-2.0) can clearly improve the predictive accuracy of the MM/PBSA and MM/GBSA methods, and computational errors are significantly reduced by including the effects of electronic polarization and using a new highly efficient interaction entropy (IE) method to calculate the entropic contribution. The suitable range for the interior dielectric constant is 1.4-1.6 for the MM/PBSA method; within this range, the correlation coefficient fluctuates around 0.84, and the mean absolute error fluctuates around 2 kcal/mol. Similarly, an interior dielectric constant of 1.8-2.0 produces a correlation coefficient of approximately 0.76 when using the MM/GBSA method. In addition, the entropic contribution of each individual residue was further calculated using the IE method to predict hot-spot residues, and the detailed binding mechanisms underlying the interactions of the HIV-1 protease, its inhibitors, and bridging water molecules were investigated. In this study, the use of a higher interior dielectric constant and the IE method can improve the calculation accuracy of the HIV-1 system.

Feasible stabilization of chondroitinase abc enables reduced astrogliosis in a chronic model of spinal cord injury

AIMS

Usually, spinal cord injury (SCI) develops into a glial scar containing extracellular matrix molecules including chondroitin sulfate proteoglycans (CSPGs). Chondroitinase ABC (ChABC), from Proteus vulgaris degrading the glycosaminoglycan (GAG) side chains of CSPGs, offers the opportunity to improve the final outcome of SCI. However, ChABC usage is limited by its thermal instability, requiring protein structure modifications, consecutive injections at the lesion site, or implantation of infusion pumps.

METHODS

Aiming at more feasible strategy to preserve ChABC catalytic activity, we assessed various stabilizing agents in different solutions and demonstrated, via a spectrophotometric protocol, that the 2.5 mol/L Sucrose solution best stabilized ChABC as far as 14 days in vitro.

RESULTS

ChABC activity was improved in both stabilizing and diluted solutions at +37°C, that is, mimicking their usage in vivo. We also verified the safety of the proposed aqueous sucrose solution in terms of viability/cytotoxicity of mouse neural stem cells (NSCs) in both proliferating and differentiating conditions in vitro. Furthermore, we showed that a single intraspinal treatment with ChABC and sucrose reduced reactive gliosis at the injury site in chronic contusive SCI in rats and slightly enhanced their locomotor recovery.

CONCLUSION

Usage of aqueous sucrose solutions may be a feasible strategy, in combination with rehabilitation, to ameliorate ChABC-based treatments to promote the regeneration of central nervous system injuries.

Modulation of the Formation of Aβ- and Sup35NM-Based Amyloids by Complex Interplay of Specific and Nonspecific Ion Effects

In vitro formation of highly ordered protein aggregates, amyloids, is influenced by the presence of ions. Here, we have studied the effect of anions on amyloid fibril formation by two different amyloidogenic proteins, human amyloid beta-42 (Aβ42), associated with Alzheimer disease and produced recombinantly with an N-terminal methionine (Met-Aβ42), and histidine-tagged NM fragment of Sup35 protein (Sup35NM-His6), a yeast release factor controlling protein-based inheritance, at pH values above and below their isoelectric points. We demonstrate here that pH plays a critical role in determining the effect of ions on the aggregation of Met-Aβ42 and Sup35NM-His6. Further, the electrophoretic mobilities of Met-Aβ42 and Sup35NM-His6 were measured in the presence of different anions at pH above and below the isoelectric points to understand how anions interact with these proteins when they bear a net positive or negative charge. We find that although ion-protein interactions generally follow expectations based on the anion positions within the Hofmeister series, there are qualitative differences in the aggregation behavior of Met-Aβ42 and Sup35NM-His6. These differences arise from a competition between nonspecific charge neutralization and screening effects and specific ion adsorption and can be explained by the different biochemical and biophysical properties of Met-Aβ42 and Sup35NM-His6.

L-arginine modified magnetic nanoparticles: green synthesis and characterization

In recent years, there has been considerable interest in Arg which is a unique, nontoxic, and biocompatible biomolecule since it can be utilized as an agent for the functionalization and subsequent stabilization of MNPs against oxidation and aggregation, during or after a synthesis procedure. Our studies demonstrate that Arg has great impacts on MNPs with the decreasing size of the particle. Also, saturation magnetization and electrostatic interactions of RMNPs have a direct impact on biological molecules such as proteins and nucleic acids. By controlling the concentration of Arg, it is possible to accurately control the above-mentioned characteristics, which are useful tools for applications such as connecting to antibodies, catalysis, drug loading, and modification of MNP stability. In the current study, three RMNPs with different Arg densities, i.e. 0.42, 1.62, and 2.29 μg per mg were successfully synthesized through a simple co-precipitation method (named p _0.5, p ₁, and p _1.5, respectively) and verified by colorimetric determination. Also, the as-synthesized RMNP powders were characterized by XRD, SEM/EDAX, FTIR, VSM, and zeta potential analysis. The presence of a magnetic core was proved by XRD, FTIR, and EDAX. Colorimetric analysis showed the existence of Arg in the synthesized samples. According to the zeta potential and VSM results, increasing the cap of Arg on the MNP surface leads to an increase in the surface charge and decrease in the magnetization of the RMNPs, respectively.

Glutamate Induced Thermal Equilibrium Intermediate and Counteracting Effect on Chemical Denaturation of Proteins

When organisms are subjected to stress conditions, one of their adaptive responses is accumulation of small organic molecules called osmolytes. These osmolytes affect the structure and stability of the biological macromolecules including proteins. The present study examines the effect of a negatively charged amino acid osmolyte, glutamate (Glu), on two model proteins, ribonuclease A (RNase A) and α-lactalbumin (α-LA), which have positive and negative surface charges at pH 7, respectively. These proteins follow two-state unfolding transitions during both heat and chemical induced denaturation processes. The addition of Glu stabilizes the proteins against temperature and induces an early equilibrium intermediate during unfolding. The stability is found to be enthalpy-driven, and the free energy of stabilization is more for α-LA compared to RNase A. The decrease in the partial molar volume and compressibility of both of the proteins in the presence of Glu suggests that the proteins attain a more compact state through surface hydration which could provide a more stable conformation. This is also supported by molecule dynamic simulation studies which demonstrate that the water density around the proteins is increased upon the addition of Glu. Further, the intermediates could be completely destabilized by lower concentrations (∼0.5 M) of guanidinium chloride and salt. However, urea subverts the Glu-induced intermediate formed by α-LA, whereas it only slightly destabilizes in the case of RNase A which has a positive surface charge and could possess charge-charge interactions with Glu. This suggests that, apart from hydration, columbic interactions might also contribute to the stability of the intermediate. Gdm-induced denaturation of RNase A and α-LA in the absence and the presence of Glu at different temperatures was carried out. These results also show the Glu-induced stabilization of both of the proteins; however, all of the unfolding transitions followed two-state transitions during chemical denaturation. The extent of stability exerted by Glu is higher for RNase A at higher temperature, whereas it provides more stability for α-LA at lower temperature. Thus, the experiments indicate that Glu induces a thermal equilibrium intermediate and increases the thermodynamic stability of proteins irrespective of their surface charges. The extent of stability varies between the proteins in a temperature-dependent manner.

Trypsin-Ligand binding affinities calculated using an effective interaction entropy method under polarized force field

Molecular dynamics (MD) simulation in the explicit water is performed to study the interaction mechanism of trypsin-ligand binding under the AMBER force field and polarized protein-specific charge (PPC) force field combined the new developed highly efficient interaction entropy (IE) method for calculation of entropy change. And the detailed analysis and comparison of the results of MD simulation for two trypsin-ligand systems show that the root-mean-square deviation (RMSD) of backbone atoms, B-factor, intra-protein and protein-ligand hydrogen bonds are more stable under PPC force field than AMBER force field. Our results demonstrate that the IE method is superior than the traditional normal mode (Nmode) method in the calculation of entropy change and the calculated binding free energy under the PPC force field combined with the IE method is more close to the experimental value than other three combinations (AMBER-Nmode, AMBER-IE and PPC-Nmode). And three critical hydrogen bonds between trypsin and ligand are broken under AMBER force field. However, they are well preserved under PPC force field. Detailed binding interactions of ligands with trypsin are further analyzed. The present work demonstrates that the polarized force field combined the highly efficient IE method is critical in MD simulation and free energy calculation.