How strongly does trehalose interact with lysozyme in the solid state? Insights from molecular dynamics simulation and inelastic neutron scattering

Essential dynamics of ubiquitin in water and in a natural deep eutectic solvent

Natural deep eutectic solvents (NADESs) comprised of osmolytes are of interest as potential biomolecular (cryo)protectants. However, the way these solvents influence the structure and dynamics of biomolecules as well as the role of water remains poorly understood. We carried out principal component analysis of various secondary structure elements of ubiquitin in water and a betaine : glycerol : water (1 : 2 : ζ; ζ = 0, 1, 2, 5, 10, 20, 45) NADES, from molecular dynamics trajectories, to gain insight into the protein dynamics as it undergoes a transition from a highly viscous anhydrous to an aqueous environment. A crossover of the protein's essential dynamics at ζ ∼ 5, induced by solvent-shell coupled fluctuations, is observed, indicating that ubiquitin might (re)fold in the NADES upon water addition at ζ > ∼5. Further, in contrast to water, the anhydrous NADES preserves ubiquitin's essential modes at high temperatures explaining the protein's seemingly enhanced thermal stability.

Uniform trehalose nanogels for glucagon stabilization

Glucagon is a peptide hormone that acts via receptor-mediated signaling predominantly in the liver to raise glucose levels by hepatic glycogen breakdown or conversion of noncarbohydrate, 3 carbon precursors to glucose by gluconeogenesis. Glucagon is administered to reverse severe hypoglycemia, a clinical complication associated with type 1 diabetes. However, due to low stability and solubility at neutral pH, there are limitations in the current formulations of glucagon. Trehalose methacrylate-based nanoparticles were utilized as the stabilizing and solubilizing moiety in the system reported herein. Glucagon was site-selectively modified to contain a cysteine at amino acid number 24 to covalently attach to the methacrylate-based polymer containing pyridyl disulfide side chains. PEG2000 dithiol was employed as the crosslinker to form uniform nanoparticles. Glucagon nanogels were monitored in Dulbecco's phosphate-buffered saline (DPBS) pH 7.4 at various temperatures to determine its long-term stability in solution. Glucagon nanogels were stable up to at least 5 months by size uniformity when stored at -20 °C and 4 °C, up to 5 days at 25 °C, and less than 12 hours at 37 °C. When glucagon stability was studied by either HPLC or thioflavin T assays, the glucagon was intact for at least 5 months at -20 °C and 4 °C within the nanoparticles at -20 °C and 4 °C and up to 2 days at 25 °C. Additionally, the glucagon nanogels were studied for toxicity and efficacy using various assays in vitro. The findings indicate that the nanogels were nontoxic to fibroblast cells and nonhemolytic to red blood cells. The glucagon in the nanogels was as active as glucagon alone. These results demonstrate the utility of trehalose nanogels towards a glucagon formulation with improved stability and solubility in aqueous solutions, particularly useful for storage at cold temperatures.

Structural Features of the Glassy State and Their Impact on the Solid-State Properties of Organic Molecules in Pharmaceutical Systems

This paper reviews the structure and properties of amorphous active pharmaceutical ingredients (APIs), including small molecules and proteins, in the glassy state (below the glass transition temperature, Tg). Amorphous materials in the neat state and formulated with excipients as miscible amorphous mixtures are included, and the role of absorbed water in affecting glass structure and stability has also been considered. We defined the term "structure" to indicate the way the various molecules in a glass interact with each other and form distinctive molecular arrangements as regions or domains of varying number of molecules, molecular packing, and density. Evidence is presented to suggest that such systems generally exist as heterogeneous structures made up of high-density domains surrounded by a lower density arrangement of molecules, termed the microstructure. It has been shown that the method of preparation and the time frame for handling and storage can give rise to variable glass structures and varying physical properties. Throughout this paper, examples are given of theoretical, computer simulation, and experimental studies which focus on the nature of intermolecular interactions, the size of heterogeneous higher density domains, and the impact of such systems on the relative physical and chemical stability of pharmaceutical systems.

Two populations of protein molecules detected by small-angle neutron and X-ray scattering (SANS and SAXS) in lyophilized protein:lyoprotector (disaccharide) systems

Two protein interaction peaks are observed in pharmaceutically-relevant protein (serum albumin) : disaccharide 1 : 1 and 1 : 3 (w/w) freeze-dried systems for the first time. In samples with a higher disaccharide content, the protein-protein distances are longer for both populations, while the fraction of the protein population with a shorter protein-protein distance is lower. Both factors would favor better stability against aggregation for disaccharide-rich protein formulations. This study provides direct experimental support for a "dilution" hypothesis as a potential stabilization mechanism for freeze-dried protein formulations.

Protecting Proteins from Desiccation Stress Using Molecular Glasses and Gels

Faced with desiccation stress, many organisms deploy strategies to maintain the integrity of their cellular components. Amorphous glassy media composed of small molecular solutes or protein gels present general strategies for protecting against drying. We review these strategies and the proposed molecular mechanisms to explain protein protection in a vitreous matrix under conditions of low hydration. We also describe efforts to exploit similar strategies in technological applications for protecting proteins in dry or highly desiccated states. Finally, we outline open questions and possibilities for future explorations.

An experimentally representative in-silico protocol for dynamical studies of lyophilised and weakly hydrated amorphous proteins

Characterization of biopolymers in both dry and weakly hydrated amorphous states has implications for the pharmaceutical industry since it provides understanding of the effect of lyophilisation on stability and biological activity. Atomistic Molecular Dynamics (MD) simulations probe structural and dynamical features related to system functionality. However, while simulations in homogenous aqueous environments are routine, dehydrated model assemblies are a challenge with systems investigated in-silico needing careful consideration; simulated systems potentially differing markedly despite seemingly negligible changes in procedure. Here we propose an in-silico protocol to model proteins in lyophilised and weakly hydrated amorphous states that is both more experimentally representative and routinely applicable. Since the outputs from MD align directly with those accessed by neutron scattering, the efficacy of the simulation protocol proposed is shown by validating against experimental neutron data for apoferritin and insulin. This work also highlights that without cooperative experimental and simulative data, development of simulative procedures using MD alone would prove most challenging.

Influence of water and trehalose on α- and β-relaxation of freeze-dried lysozyme formulations

Molecular mobility in form of alpha and beta relaxations is considered crucial for characterization of amorphous lyophilizates and reflected in the transition temperatures Tgα and Tgβ. Based on an overview of applied methods to study beta relaxations, Dynamic Mechanical analysis was used to measure Tgα and Tgβ in amorphous freeze-dried samples. Lysozyme and trehalose as well as their mixtures in varying ratios were investigated. Three different residual moisture levels, ranging from roughly 0.5-7 % (w/w), were prepared via equilibration of the freeze-dried samples. Known plasticising effects of water on Tgα were confirmed, also via differential scanning calorimetry. In addition and contrary to expectations, an influence of water on the Tgβ also was observed. On the other hand, an increasing amount of trehalose lowered Tgα but increased Tgβ showing that Tgα and Tgβ are not paired. The findings were interpreted with regard to their underlying molecular mechanisms and a correlation with the known influences of water and trehalose on stability. The results provide encouraging hints for future stability studies of freeze-dried protein formulations, which are urgently needed, not least for reasons of sustainability.

Protein stability in a natural deep eutectic solvent: Preferential hydration or solvent slaving?

Deep eutectic solvents (DESs) emerged as potential alternative solvent media in multiple areas, including biomolecular (cryo)preservation. Herein, we studied the stability of a small protein (ubiquitin) in water and a betaine-glycerol-water (B:G:W) (1:2:ζ; ζ = 0, 1, 2, 5, 10) DES, through molecular dynamics. An AMBER-based model that accurately describes the density and shear viscosity of the DES is proposed. We find that water molecules are largely trapped in the solvent, precluding the formation of a full hydration layer, seemingly opposite to osmolytes' preferential exclusion/preferential hydration mechanism. Although the protein is stable in the DES, structural fluctuations are largely suppressed and only recovered upon sufficient hydration. This is explained by a solvent-slaving mechanism where β-fluctuations are key, with the non-monotonic hydration of some amino acids with the water content providing an explanation to the non-monotonic folding of some proteins in aqueous DESs. A major thermal stability enhancement in the DES is also observed, caused by a similar slowdown of the backbone torsional dynamics. Our results support a kinetic stabilization of the protein in the DES, whereas a possible thermodynamic stabilization does not follow a preferential hydration or water entrapment mechanism.

New insights into the protein stabilizing effects of trehalose by comparing with sucrose

Disaccharides are well known to be efficient stabilizers of proteins, for example in the case of lyophilization or cryopreservation. However, although all disaccharides seem to exhibit bioprotective and stabilizing properties, it is clear that trehalose is generally superior compared to other disaccharides. The aim of this study was to understand this by comparing how the structural and dynamical properties of aqueous trehalose and sucrose solutions influence the protein myoglobin (Mb). The structural studies were based on neutron and X-ray diffraction in combination with empirical potential structure refinement (EPSR) modeling, whereas the dynamical studies were based on quasielastic neutron scattering (QENS) and molecular dynamics (MD) simulations. The results show that the overall differences in the structure and dynamics of the two systems are small, but nevertheless there are some important differences which may explain the superior stabilizing effects of trehalose. It was found that in both systems the protein is preferentially hydrated by water, but that this effect is more pronounced for trehalose, i.e. trehalose forms less hydrogen bonds to the protein surface than sucrose. Furthermore, the rotational motion around dihedrals between the two glucose rings of trehalose is slower than in the case of the dihedrals between the glucose and fructose rings of sucrose. This leads to a less perturbed protein structure in the case of trehalose. The observations indicate that an aqueous environment closest to the protein molecules is beneficial for an efficient bioprotective solution.

How Sugars Protect Dry Protein Structure

Extremotolerant organisms and industry exploit sugars as desiccation protectants, with trehalose being widely used by both. How sugars, in general, and the hydrolytically stable sugar trehalose, in particular, protect proteins is poorly understood, which hinders the rational design of new excipients and implementation of novel formulations for preserving lifesaving protein drugs and industrial enzymes. We employed liquid-observed vapor exchange nuclear magnetic resonance (LOVE NMR), differential scanning calorimetry (DSC), and thermal gravimetric analysis (TGA) to show how trehalose and other sugars protect two model proteins: the B1 domain of streptococcal protein G (GB1) and truncated barley chymotrypsin inhibitor 2 (CI2). Residues with intramolecular H-bonds are most protected. The LOVE NMR and DSC data indicate that vitrification may be protective. Combining LOVE NMR and TGA data shows that water retention is not important. Our data suggest that sugars protect protein structure as they dry by strengthening intraprotein H-bonds and water replacement and that trehalose is the stress-tolerance sugar of choice because of its covalent stability.

A Long Journey into the Investigation of the Structure–Dynamics–Function Paradigm in Proteins through the Activities of the Palermo Biophysics Group

An overview of the biophysics activity at the Department of Physics and Chemistry Emilio Segrè of the University of Palermo is given. For forty years, the focus of the research has been on the protein structure–dynamics–function paradigm, with the aim of understanding the molecular basis of the relevant mechanisms and the key role of solvent. At least three research lines are identified; the main results obtained in collaboration with other groups in Italy and abroad are presented. This review is dedicated to the memory of Professors Massimo Ugo Palma, Maria Beatrice Palma Vittorelli, and Lorenzo Cordone, which were the founders of the Palermo School of Biophysics. We all have been, directly or indirectly, their pupils; we miss their enthusiasm for scientific research, their deep physical insights, their suggestions, their strict but always constructive criticisms, and, most of all, their friendship. This paper is dedicated also to the memory of Prof. Hans Frauenfelder, whose pioneering works on nonexponential rebinding kinetics, protein substates, and energy landscape have inspired a large part of our work in the field of protein dynamics.

Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme

Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.

Reversible processes in collagen dehydration: A molecular dynamics study

Collagen dehydration is an unavoidable damaging process that causes the lack of fibers' physical properties and it is usually irreversible. However, the identification of low hydration conditions that permit a recovering of initial collagen features after a rehydration treatment is particularly of interest. Monitoring structural changes by means of MD simulations, we investigated the hydration-dehydration-rehydration cycle of two microfibril models built on different fragments of the sequence of rat tail collagen type I. The microfibrils have different hydropathic features, to investigate the influence of amino acid composition on the whole process. We showed that with low hydration at a level corresponding to the first shell, microfibril gains in compactness and tubularity. Crucially, some water molecules remain trapped inside the fibrils, allowing, by rehydrating, a recovery of the initial collagen structural features. Water rearranges in cluster around the protein, and its first layer is more anchored to the surface. However, these changes in distribution and mobility in low hydration conditions get back with rehydration.

Unraveling the Influence of Osmolytes on Water Hydrogen-Bond Network: From Local Structure to Graph Theory Analysis

Water structure in aqueous osmolyte solutions, deduced from the slight alteration in the water-water radial distribution function, the decrease in water-water hydrogen bonding, and tetrahedral ordering based only on the orientation of nearest water molecules derived from the molecular dynamics simulations, appears to have been perturbed. A careful analysis, however, reveals that the hydrogen bonding and the tetrahedral ordering around a water molecule in binary solutions remain intact as in neat water when the contribution of osmolyte-water interactions is appropriately incorporated. Furthermore, the distribution of the water binding energies and the water excess chemical potential of solvation in solutions are also pretty much the same as in neat water. Osmolytes are, therefore, well integrated into the hydrogen-bond network of water. Indeed, osmolytes tend to preferentially hydrogen bond with water molecules and their interaction energies are strongly correlated to their hydrogen-bonding capability. The graph network analysis, further, illustrates that osmolytes act as hubs in the percolated hydrogen-bond network of solutions. The degree of hydrogen bonding of osmolytes predominantly determines all of the network properties. Osmolytes like ethanol that form fewer hydrogen bonds than a water molecule disrupt the water hydrogen-bond network, while other osmolytes that form more hydrogen bonds effectively increase the connectivity among water molecules. Our observation of minimal variation in the local structure and the vitality of osmolyte-water hydrogen bonds on the solution network properties clearly imply that the direct interaction between protein and osmolytes is solely responsible for the protein stability. Further, the relevance of hydrogen bonds on solution properties suggests that the hydrogen-bonding interaction among protein, water, and osmolyte could be the key determinant of the protein conformation in solutions.

"Water Association" Band in Saccharide Amorphous Matrices: Role of Residual Water on Bioprotection

Saccharides protect biostructures against adverse environmental conditions mainly by preventing large scale motions leading to unfolding. The efficiency of this molecular mechanism, which is higher in trehalose with respect to other sugars, strongly depends on hydration and sugar/protein ratio. Here we report an Infrared Spectroscopy study on dry amorphous matrices of the disaccharides trehalose, maltose, sucrose and lactose, and the trisaccharide raffinose. Samples with and without embedded protein (Myoglobin) are investigated at different sugar/protein ratios, and compared. To inspect matrix properties we analyse the Water Association Band (WAB), and carefully decompose it into sub-bands, since their relative population has been shown to effectively probe water structure and dynamics in different matrices. In this work the analysis is extended to investigate the structure of protein-sugar-water samples, for the first time. Results show that several classes of water molecules can be identified in the protein and sugar environment and that their relative population is dependent on the type of sugar and, most important, on the sugar/protein ratio. This gives relevant information on how the molecular interplay between residual waters, sugar and protein molecules affect the biopreserving properties of saccharides matrices.

The biology of tardigrade disordered proteins in extreme stress tolerance

Abstract

Disordered proteins have long been known to help mediate tolerance to different abiotic stresses including freezing, osmotic stress, high temperatures, and desiccation in a diverse set of organisms. Recently, three novel families of intrinsically disordered proteins were identified in tardigrades, microscopic animals capable of surviving a battery of environmental extremes. These three families include the Cytoplasmic-, Secreted-, and Mitochondrial- Abundant Heat Soluble (CAHS, SAHS, and MAHS) proteins, which are collectively termed Tardigrade Disordered Proteins (TDPs). At the level of sequence conservation TDPs are unique to tardigrades, and beyond their high degree of disorder the CAHS, SAHS, and MAHS families do not resemble one another. All three families are either highly expressed constitutively, or significantly enriched in response to desiccation. In vivo, ex vivo, and in vitro experiments indicate functional roles for members of each TDP family in mitigating cellular perturbations induced by various abiotic stresses. What is currently lacking is a comprehensive and holistic understanding of the fundamental mechanisms by which TDPs function, and the properties of TDPs that allow them to function via those mechanisms. A quantitative and systematic approach is needed to identify precisely what cellular damage TDPs work to prevent, what sequence features are important for these functions, and how those sequence features contribute to the underlying mechanisms of protection. Such an approach will inform us not only about these fascinating proteins, but will also provide insights into how the sequence of a disordered protein can dictate its functional, structural, and dynamic properties.

Water confined in solutions of biological relevance

In many relevant situations, water is not in its bulk form but instead attached to some substrates or filling some cavities. We shall call water in the latter environment confined water as opposed to bulk water. It is known that the confined water is essential for the stability and the function of biological macromolecules. In this paper, we provide a review of the experimental and computational advances over the past decades concerning the understanding of the structure and dynamics of water confined in aqueous solutions of biological relevance. Examples involving water in solution of organic solutes (cryoprotectants such as dimethylsulfoxide (DMSO), sugars such as trehalose) are provided.

New insight into the molecular mechanism of the trehalose effect on urate oxidase stability

Urate oxidase (EC 1.7.3.3) is a key enzyme in the purine metabolism which is applied in the treatment of gout and also, as a diagnostic reagent for uric acid detection. In the current study, the trehalose (TRE) effects as an additive on the structural stability and function of uricase were investigated. For recombinant expression of UOX in E. coli BL21 cells, firstly the coding sequence was subcloned into the pET-28a vector and after induction with IPTG, the recombinant UOX was purified by affinity chromatography using a Ni-NTA agarose column. To specify the trehalose effects on the urate oxidase (UOX) structure, optimum pH, optimum temperature, kinetic and thermodynamic parameters and also, the intrinsic fluorescence of UOX in the absence and presence of trehalose were examined. The UOX half-life is 24.32 min at 40 °C, whereas the UOX-TRE has a higher half-life (32.09 min) at this temperature. Generally, our findings confirm that trehalose has a protective effect on the enzyme structure. Optimum pH and temperature were 9 and 25 °C, respectively for both the naked and treated enzymes and their activity retained 42.18 and 64.80%, respectively after 48 h of incubation at room temperature. Also, theoretical results indicate that the random coil of the enzyme was converted to α-helix and β-sheet in the presence of trehalose which may preserve the integrity of the active site of the enzyme and increased the enzymatic activity. The MD simulation results indicated greater stability of the uricase structure in the presence of trehalose.Communicated by Ramaswamy H. Sarma.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Mechanism of Trehalose-Induced Protein Stabilization from Neutron Scattering and Modeling

The sugar molecule trehalose has been proven to be an excellent stabilizing cosolute for the preservation of biological materials. However, the stabilizing mechanism of trehalose has been much debated during the previous decades, and it is still not fully understood, partly because it has not been completely established how trehalose molecules structure around proteins. Here, we present a molecular model of a protein-water-trehalose system, based on neutron scattering results obtained from neutron diffraction, quasielastic neutron scattering, and different computer modeling techniques. The structural data clearly show how the proteins are preferentially hydrated, and analysis of the dynamical properties show that the protein residues are slowed down because of reduced dynamics of the protein hydration shell, rather than because of direct trehalose-protein interactions. These findings, thereby, strongly support previous models related to the preferential hydration model and contradict other models based on water replacement at the protein surface. Furthermore, the results are important for understanding the specific role of trehalose in biological stabilization and, more generally, for providing a likely mechanism of how cosolutes affect the dynamics of proteins.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

Bioprotection Can Be Tuned with a Proper Protein/Saccharide Ratio: The Case of Solid Amorphous Matrices

Saccharides, and in particular trehalose, are well known for their high efficiency in protecting biostructures against adverse environmental conditions. The protein dynamics is known to be highly inhibited in a low-water trehalose host medium, the inhibition being markedly dependent on the amount of residual water. Besides hydration, the protein/sugar ratio is expected to affect the properties of saccharide amorphous matrices. In this work, we report an infrared spectroscopy study in dry amorphous matrices of various sugars (the disaccharides trehalose, maltose, sucrose, and lactose, and the trisaccharide raffinose) containing myoglobin, at different protein/sugar ratios. We analyze the stretching band of the bound CO molecule and the water association band. Such bands have already been successfully exploited for the simultaneous study of thermal evolution of a matrix and embedded protein. The results show a high dependence of protein and matrix signals on the protein/sugar ratio, the system behavior evolving from situations where (i) the protein slaves the matrix to (ii) protein ↔ matrix coupling/uncoupling, then to (iii) the matrix slaving the protein, with increasing sugar concentration. This supports a mutual protein ↔ matrix structural and dynamic influence in low hydrated systems, indicating that the protein/solvent master and slave paradigm does not strictly hold, but the mutual relationship depends on the relative concentrations. Furthermore, for each sugar, an optimal protein/sugar concentration ratio can be identified, which maximizes the protein preservation; under such a condition, the water content is minimal.

Freezing of Aqueous Solutions and Chemical Stability of Amorphous Pharmaceuticals: Water Clusters Hypothesis

Molecular mobility has been traditionally invoked to explain physical and chemical stability of diverse pharmaceutical systems. Although the molecular mobility concept has been credited with creating a scientific basis for stabilization of amorphous pharmaceuticals and biopharmaceuticals, it has become increasingly clear that this approach represents only a partial description of the underlying fundamental principles. An additional mechanism is proposed herein to address 2 key questions: (1) the existence of unfrozen water (i.e., partial or complete freezing inhibition) in aqueous solutions at subzero temperatures and (2) the role of water in the chemical stability of amorphous pharmaceuticals. These apparently distant phenomena are linked via the concept of water clusters. In particular, freezing inhibition is associated with the confinement of water clusters in a solidified matrix of an amorphous solute, with nanoscaled water clusters being observed in aqueous glasses using wide-angle neutron scattering. The chemical instability is suggested to be directly related to the catalysis of proton transfer by water clusters, considering that proton transfer is the key elementary reaction in many chemical processes, including such common reactions as hydrolysis and deamidation.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

Molecular Packing, Hydrogen Bonding, and Fast Dynamics in Lysozyme/Trehalose/Glycerol and Trehalose/Glycerol Glasses at Low Hydration

Water and glycerol are well-known to facilitate the structural relaxation of amorphous protein matrices. However, several studies evidenced that they may also limit fast (∼picosecond-nanosecond, ps-ns) and small-amplitude (∼Å) motions of proteins, which govern their stability in freeze-dried sugar mixtures. To determine how they interact with proteins and sugars in glassy matrices and, thereby, modulate their fast dynamics, we performed molecular dynamics (MD) simulations of lysozyme/trehalose/glycerol (LTG) and trehalose/glycerol (TG) mixtures at low glycerol and water concentrations. Upon addition of glycerol and/or water, the glass transition temperature, T_g, of LTG and TG mixtures decreases, the molecular packing of glasses is improved, and the mean-square displacements (MSDs) of lysozyme and trehalose either decrease or increase, depending on the time scale and on the temperature considered. A detailed analysis of the hydrogen bonds (HBs) formed between species reveals that water and glycerol may antiplasticize the fast dynamics of lysozyme and trehalose by increasing the total number and/or the strength of the HBs they form in glassy matrices.

The water association band as a marker of hydrogen bonds in trehalose amorphous matrices

The water association band is a suitable marker of residual water behavior in bioprotective trehalose matrices.

Investigating Structure and Dynamics of Proteins in Amorphous Phases Using Neutron Scattering

In order to increase shelf life and minimize aggregation during storage, many biotherapeutic drugs are formulated and stored as either frozen solutions or lyophilized powders. However, characterizing amorphous solids can be challenging with the commonly available set of biophysical measurements used for proteins in liquid solutions. Therefore, some questions remain regarding the structure of the active pharmaceutical ingredient during freezing and drying of the drug product and the molecular role of excipients. Neutron scattering is a powerful technique to study structure and dynamics of a variety of systems in both solid and liquid phases. Moreover, neutron scattering experiments can generally be correlated with theory and molecular simulations to analyze experimental data. In this article, we focus on the use of neutron techniques to address problems of biotechnological interest. We describe the use of small-angle neutron scattering to study the solution structure of biological molecules and the packing arrangement in amorphous phases, that is, frozen glasses and freeze-dried protein powders. In addition, we discuss the use of neutron spectroscopy to measure the dynamics of glassy systems at different time and length scales. Overall, we expect that the present article will guide and prompt the use of neutron scattering to provide unique insights on many of the outstanding questions in biotechnology.

Structural and thermodynamic aspects of plasticization and antiplasticization in glassy encapsulation and biostabilization matrices

The structural and thermodynamic properties of glassy carbohydrate matrices for the encapsulation and biostabilization of sensitive bioactive compounds, such as pharmaceutically active proteins and oxidation-sensitive compounds, are reviewed in the context of the plasticization and antiplasticization of glassy carbohydrates of intermediate and high molecular weight by low molecular weight diluents. Plasticization and antiplasticization may be monitored either by dynamic measures or by structural and thermodynamic features of the glassy matrices. Specifically, it is demonstrated that the decrease in size of the molecular free volume holes with increasing diluent content, as determined by positron annihilation lifetime spectroscopy (PALS), is related to the antiplasticization of glassy carbohydrate matrices, resulting in increased barrier properties of the glassy matrix. As far as could be ascertained from the available data, the regimes as identified by PALS map on those detected by neutron scattering and dielectric spectroscopy for glassy matrices consisting of trehalose and the diluent glycerol. The review is concluded by a survey of the published results on the stability of bioactive compounds encapsulated in carbohydrate glasses and an overview of outstanding questions.

Size and molecular flexibility of sugars determine the storage stability of freeze-dried proteins

Protein-based biopharmaceuticals are generally produced as aqueous solutions and stored refrigerated to obtain sufficient shelf life. Alternatively, proteins may be freeze-dried in the presence of sugars to allow storage stability at ambient conditions for prolonged periods. However, to act as a stabilizer, these sugars should remain in the glassy state during storage. This requires a sufficiently high glass transition temperature (Tg). Furthermore, the sugars should be able to replace the hydrogen bonds between the protein and water during drying. Frequently used disaccharides are characterized by a relatively low Tg, rendering them sensitive to plasticizing effects of residual water, which strongly reduces the Tg values of the formulation. Larger sugars generally have higher Tgs, but it is assumed that these sugars are limited in their ability to interact with the protein due to steric hindrance. In this paper, the size and molecular flexibility of sugars was related to their ability to stabilize proteins. Four diverse proteins varying in size from 6 kDa to 540 kDa were freeze-dried in the presence of different sugars varying in size and molecular flexibility. Subsequently, the different samples were subjected to an accelerated stability test. Using protein specific assays and intrinsic fluorescence, stability of the proteins was monitored. It was found that the smallest sugar (disaccharide trehalose) best preserved the proteins, but also that the Tg of the formulations was only just high enough to maintain sufficient vitrification. When trehalose-based formulations are exposed to high relative humidities, water uptake by the product reduces the Tgs too much. In that respect, sugars with higher Tgs are desired. Addition of polysaccharide dextran 70 kDa to trehalose greatly increased the Tg of the formulation. Moreover, this combination also improved the stability of the proteins compared to dextran only formulations. The molecularly flexible oligosaccharide inulin 4 kDa provided better stabilization than the similarly sized but molecularly rigid oligosaccharide dextran 6 kDa. In conclusion, the results of this study indicate that size and molecular flexibility of sugars affect their ability to stabilize proteins. As long as they maintain vitrified, smaller and molecularly more flexible sugars are less affected by steric hindrance and thus better capable at stabilizing proteins.

Why a diaminopyrrolic tripodal receptor binds mannosides in acetonitrile but not in water?

Intermolecular interactions involving carbohydrates and their natural receptors play important roles in several biological processes. The development of synthetic receptors is very useful to study these recognition processes. Recently, it was synthetized a diaminopyrrolic tripodal receptor that is selective for mannosides, which are obtained from mannose, a sugar with significant relevance in living systems. However, this receptor is significantly more active in acetonitrile than in water. In this work, we performed several molecular dynamics and constant-pH molecular dynamics simulations in acetonitrile and water to evaluate the conformational space of the receptor and to understand the molecular detail of the receptor–mannoside interaction. The protonation states sampled by the receptor show that the positive charges are always as distant as possible in order to avoid large intramolecular repulsions. Moreover, the conformational space of the receptor is very similar in water above pH 4.0 and in acetonitrile. From the simulations with the mannoside, we observe that the interactions are more specific in acetonitrile (mainly hydrogen bonds) than in water (mainly hydrophobic). Our results suggest that the readiness of the receptor to bind mannoside is not significantly affected in water (above pH 4.0). Probably, the hydrogen bond network that is formed in acetonitrile (which is weaker in water) is the main reason for the higher activity in this solvent. This work also presents a new implementation of the stochastic titration constant-pH molecular dynamics method to a synthetic receptor of sugars and attests its ability to describe the protonation/conformation coupling in these molecules.

The importance of individual protein molecule dynamics in developing and assessing solid state protein preparations

Processing protein solutions into the solid state is a common approach for generating stable amorphous protein mixtures that are suitable for long-term storage. Great care is typically given to protecting the protein native structure during the various drying steps that render it into the amorphous solid state. However, many studies illustrate that chemical and physical degradations still occur in spite of this amorphous material having good glassy properties and it being stored at temperatures below its glass transition temperature (Tg). Because of these persistent issues and recent biophysical studies that have refined the debate ascribing meaning to the molecular dynamical transition temperature and Tg of protein molecules, we provide an updated discussion on the impact of assessing and managing localized, individual protein molecule nondiffusive motions in the context of proteins being prepared into bulk amorphous mixtures. Our aim is to bridge the pharmaceutical studies addressing bulk amorphous preparations and their glassy behavior, with the biophysical studies historically focused on the nondiffusive internal protein dynamics and a protein's activity, along with their combined efforts in assessing the impact of solvent hydrogen-bonding networks on local stability. We also provide recommendations for future research efforts in solid-state formulation approaches.

Structural and dynamical characteristics of trehalose and sucrose matrices at different hydration levels as probed by FTIR and high-field EPR

Some organisms can survive complete dehydration and high temperatures by adopting an anhydrobiotic state in which the intracellular medium contains large amounts of disaccharides, particularly trehalose and sucrose. Trehalose is most effective also in protecting isolated in vitro biostructures. In an attempt to clarify the molecular mechanisms of disaccharide bioprotection, we compared the structure and dynamics of sucrose and trehalose matrices at different hydration levels by means of high-field W-band EPR and FTIR spectroscopy. The hydration state of the samples was characterized by FTIR spectroscopy and the structural organization was probed by EPR using a nitroxide radical dissolved in the respective matrices. Analysis of the EPR spectra showed that the structure and dynamics of the dehydrated matrices as well as their evolution upon re-hydration differ substantially between trehalose and sucrose. The dehydrated trehalose matrix is homogeneous in terms of distribution of the residual water and spin-probe molecules. In contrast, dehydrated sucrose forms a heterogeneous matrix. It is comprised of sucrose polycrystalline clusters and several bulk water domains. The amorphous form was found only in 30% (volume) of the sucrose matrix. Re-hydration leads to a structural homogenization of the sucrose matrix, whilst in the trehalose matrix several domains develop differing in the local water/radical content and radical mobility. The molecular model of the matrices provides an explanation for the different protein-matrix dynamical coupling observed in dried ternary sucrose and trehalose matrices, and accounts for the superior efficacy of trehalose as a bioprotectant. Furthermore, for bacterial photosynthetic reaction centers it is shown that at low water content the protein-matrix coupling is modulated by the sugar/protein molar ratio in sucrose matrices only. This effect is suggested to be related to the preference for sucrose, rather than trehalose, as a bioprotective disaccharide in some anhydrobiotic organisms.

NMR study on the interaction of trehalose with lactose and its effect on the hydrogen bond interaction in lactose

Trehalose, a well-known stress-protector of biomolecules, has been investigated for its effect on the mobility, hydration and hydrogen bond interaction of lactose using diffusion-ordered NMR spectroscopy and NMR of hydroxy protons. In ternary mixtures of trehalose, lactose and water, the two sugars have the same rate of diffusion. The chemical shifts, temperature coefficients, vicinal coupling constants and ROE of the hydroxy protons in trehalose, lactose and sucrose were measured for the disaccharides alone in water/acetone-d₆ solutions as well as in mixtures. The data indicated that addition of trehalose did not change significantly the strength of the hydrogen bond interaction between GlcOH3 and GalO5' in lactose. Small upfield shifts were however measured for all hydroxy protons when the sugar concentration was increased. The chemical shift of the GlcOH3 signal in lactose showed less change, attributed to the spatial proximity to GalO5'. Chemical exchange between hydroxy protons of lactose and trehalose was observed in the ROESY NMR spectra. Similar effects were observed with sucrose indicating no specific effect of trehalose at the concentrations investigated (73 to 763 mg/mL) and suggesting that it is the concentration of hydroxy groups more than the type of sugars which is guiding intermolecular interactions.

Proteins in amorphous saccharide matrices: structural and dynamical insights on bioprotection

Bioprotection by sugars, and in particular trehalose peculiarity, is a relevant topic due to the implications in several fields. The underlying mechanisms are not yet clearly elucidated, and remain the focus of current investigations. Here we revisit data obtained at our lab on binary sugar/water and ternary protein/sugar/water systems, in wide ranges of water content and temperature, in the light of the current literature. The data here discussed come from complementary techniques (Infrared Spectroscopy, Molecular Dynamics simulations, Small Angle X-ray Scattering and Calorimetry), which provided a consistent description of the bioprotection by sugars from the atomistic to the macroscopic level. We present a picture, which suggests that protein bioprotection can be explained in terms of a strong coupling of the biomolecule surface to the matrix via extended hydrogen-bond networks, whose properties are defined by all components of the systems, and are strongly dependent on water content. Furthermore, the data show how carbohydrates having similar hydrogen-bonding capabilities exhibit different efficiency in preserving biostructures.

Microscopic mechanism of protein cryopreservation in an aqueous solution with trehalose

In order to investigate the cryoprotective mechanism of trehalose on proteins, we use molecular dynamics computer simulations to study the microscopic dynamics of water upon cooling in an aqueous solution of lysozyme and trehalose. We find that the presence of trehalose causes global retardation of the dynamics of water. Comparing aqueous solutions of lysozyme with/without trehalose, we observe that the dynamics of water in the hydration layers close to the protein is dramatically slower when trehalose is present in the system. We also analyze the structure of water and trehalose around the lysozyme and find that the trehalose molecules form a cage surrounding the protein that contains very slow water molecules. We conclude that the transient cage of trehalose molecules that entraps and slows the water molecules prevents the crystallisation of protein hydration water upon cooling.

The effect of complex solvents on the structure and dynamics of protein solutions: The case of Lysozyme in trehalose/water mixtures

We present a Molecular Dynamics simulation study of the effect of trehalose concentration on the structure and dynamics of individual proteins immersed in trehalose/water mixtures. Hen egg-white Lysozyme is used in this study and trehalose concentrations of 0%, 10%, 20%, 30% and 100% by weight are explored. Surprisingly, we have found that changes in trehalose concentration do not change the global structural characteristics of the protein as measured by standard quantities like the mean square deviation, radius of gyration, solvent accessible surface area, inertia tensor and asphericity. Only in the limit of pure trehalose these metrics change significantly. Specifically, we found that the protein is compressed by 2% when immersed in pure trehalose. At the amino acid level there is noticeable rearrangement of the surface residues due to the change in polarity of the surrounding environment with the addition of trehalose. From a dynamic perspective, our computation of the Incoherent Intermediate Scattering Function shows that the protein slows down with increasing trehalose concentration; however, this slowdown is not monotonic. Finally, we also report in-depth results for the hydration layer around the protein including its structure, hydrogen-bonding characteristics and dynamic behavior at different length scales.