Fourier transform infrared spectroscopic and theoretical study of water interactions with glycine and its N-methylated derivatives

Specific Proton-Donor Properties of Glycine Betaine. Metric Parameters and Enthalpy of Noncovalent Interactions in its Dimer, Water Complexes and Crystalline Hydrate

Trimethylglycine (glycine betaine, GB) is an important organic osmolyte that accumulates in various plant species in response to environmental stresses and has significant potential as a bioactive agent with low environmental impact. It is assumed that the hydration of GB is playing an important role in the protective mechanism. The hydration and aggregation properties of GB have not yet been studied in detail at the atomistic level. In this work, noncovalent interactions in the GB dimer and its complexes with water and crystalline monohydrate are studied. Depending on the object, periodic and non-periodic DFT calculations are used. Particular attention is paid to the metric parameters and enthalpies of intermolecular hydrogen bonds. The identification of noncovalent interactions is carried out by means of the Bader analysis of periodic or non-periodic electron density. The enthalpy of hydrogen bonds is estimated using the Rosenberg formula (PCCP 2 (2000) 2699). The specific proton donor properties of glycine betaine are due to its ability to form intermolecular C-H∙∙∙O bonds with the oxygen atom of a water molecule or the carboxylate group of a neighboring GB. The enthalpy of these bonds can be significantly greater than 10 kJ/mol. The water molecule that forms a hydrogen bond with the carboxylate group of GB also interacts with its CH groups through lone pairs of electrons. The C-H∙∙∙O bonds contribute up to 40% of the total entropy of the GB-water interaction, which is about 45 kJ/mol. The possibility of identifying C-H∙∙∙O bonds by the proton nuclear magnetic resonance method is discussed.

Interaction of glycine with Li+ in the (H2O)n (n = 0-8) clusters

CONTEXT

We investigated the interaction between glycine and Li⁺ in water environment based on the Gly·Li⁺(H₂O)_n (n = 0-8) cluster. Our study shows that for Gly·Li⁺, Li⁺ binds to both carbonyl oxygen and amino nitrogen to form a bidentate structure, and the first three water molecules preferentially interact with Li⁺. For n = 0-5, the complexes of Gly·Li⁺(H₂O)_n exist in neutral form, and when the water number reached 6, the complex can coexist in neutral and zwitterionic form, then zwitterionic structures are dominant for n = 7, 8. The analyses by RDG, AIM, and ESP in conjunction with the calculated interaction energies show that the interaction between Li⁺ and Gly decreases gradually with the water molecules involved successively from n = 1 to 6 and then increases for n = 7-8. Additionally, the infrared spectra of Gly·Li⁺(H₂O)_n (n = 0-8) are also calculated.

METHODS

The initial structures were optimized using Gaussian 09 program package in B3LYP-D3 (BJ)/6-311G(d, p) method, and the frequency was calculated with 6-311 + G(2d, p) basis set. GaussView5.0.9 was used to view simulation infrared spectra. The noncovalent interaction method (NCl), energy decomposition (EDA), atoms in molecules (AIM) analysis, and electrostatic potential (ESP) analyses were conducted using Multiwfn software to gain a deeper understanding of the interaction properties of Gly, Li⁺, and water.

Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations

We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.

Mechanism of Osmolyte Stabilization-Destabilization of Proteins: Experimental Evidence

In this work, we investigated the influence of stabilizing (N,N,N-trimethylglycine) and destabilizing (urea) osmolytes on the hydration spheres of biomacromolecules in folded forms (trpzip-1 peptide and hen egg white lysozyme—hewl) and unfolded protein models (glycine—GLY and N-methylglycine—NMG) by means of infrared spectroscopy. GLY and NMG were clearly limited as minimal models for unfolded proteins and should be treated with caution. We isolated the spectral share of water changed simultaneously by the biomacromolecule/model molecule and the osmolyte, which allowed us to provide unambiguous experimental arguments for the mechanism of stabilization/destabilization of proteins by osmolytes. In the case of both types of osmolytes, the decisive factor determining the equilibrium folded/unfolded state of protein was the enthalpy effect exerted on the hydration spheres of proteins in both forms. In the case of stabilizing osmolytes, enthalpy was also favored by entropy, as the unfolded state of a protein was more entropically destabilized than the folded state.

Hydration of the Zwitterionic and Protonated Forms of Glycine Betaine Probed by Soft X-ray Emission Spectroscopy Coupled with Chemometrics

Soft X-ray absorption and emission spectra of glycine betaine (GB) have been measured at the O K-edge in neutral and strongly acidic solutions. The absorption spectra of the neutral solutions have a resonance peak at 532.6 eV, assigned to the transition to the π* orbital, whereas in the acidic solutions, the peak is shifted by -0.3 eV. The emission spectra taken as a function of the GB concentration have been analyzed by means of a modified classical least-squares regression method to obtain the hydration number of the solute. The analysis is successful when the emission spectra have been acquired at the energy of a slightly detuned resonance, giving 28 and 24 as the minimum values for the zwitterionic and protonated GB, respectively. The number of 28 accords with the reported values for the number of water molecules in the first hydration layer of the zwitterion and is greater than that obtained by other experimental techniques. The obtained numbers are used to discuss the hydration structure of GB with the aid of ab initio molecular orbital calculations. The hydration structure of the protonated form of GB is explored for the first time.

Distinctive behavior and two-dimensional vibrational dynamics of water molecules inside glycine solvation shell

We present a first principles molecular dynamics study of a deuterated aqueous solution of a single glycine moiety to explore the structure, dynamics, and two-dimensional infrared spectra of water molecules found in the solvation shell of glycine.

Understanding the function of water during the gelation of globular proteins by temperature-dependent near infrared spectroscopy

Water plays an indispensable role in the gelation of proteins, but its function still remains unclear. In this work, the variation of water species with the structural changes of globular proteins was investigated using temperature-dependent near infrared (NIR) spectroscopy. Ovalbumin (OVA) was used as a model protein, which forms a gel-like structure as the temperature increases through three phases, i.e., phase I (native), phase II (molten globule state), and phase III (gel state). The structural change and the content variation of different water species in the three phases of gelation were analyzed by two-dimensional correlation NIR spectroscopy and Gaussian fitting. A decrease in the water species with two hydrogen bonds (S2) was found and the change follows the same phases as OVA. In the first two phases, the change occurs after those of other water species but in the third phase, the change is faster than that of free water species. The result indicates that in the native and molten globule states, S2 is located in the hydration shell of OVA to maintain the stability of the protein structure, and then in the gel state, high temperature weakens the hydrogen bonding of S2 and leads to the destruction of the hydration shell, making OVA clusters form a gel structure.

Phonon Spectrometric Evaluation of the Solute-Solvent Interface in Solutions of Glycine and Its N-Methylated Derivatives

From the perspective of O:H-O bond cooperativity, we analyzed the solute capability of transiting the O:H-O bond from the mode of ordinary water to the hydration state and its consequence on the solution viscosity and surface stress. Phonon spectrometric results suggest that glycine and its N-methyl derivatives strongly affect the surrounding solvent molecules through H ↔ H repulsion and dipolar polarization. The H ↔ H interproton repulsion disrupts the surface stress, and the polarization enhances the solution viscosity.

Taurine as a water structure breaker and protein stabilizer

The enhancing effect on the water structure has been confirmed for most of the osmolytes exhibiting both stabilizing and destabilizing properties in regard to proteins. The presented work concerns osmolytes, which should be classified as "structure breaking" solutes: taurine and N,N,N-trimethyltaurine (TMT). Here, we combine FTIR spectroscopy, DSC calorimetry and DFT calculations to gain an insight into the interactions between osmolytes and two proteins: lysozyme and ubiquitin. Despite high structural similarity, both osmolytes exert different influence on protein stability: taurine is a stabilizer, TMT is a denaturant. We show also that taurine amino group interacts directly with the side chains of proteins, whereas TMT does not interact with proteins at all. Although two solutes weaken on average the structure of the surrounding water, their hydration spheres are different. Taurine is surrounded by two populations of water molecules: bonded with weak H-bonds around sulfonate group, and strongly bonded around amino group. The strong hydrogen-bonded network of water molecules around the amino group of taurine further improves properties of enhanced protein hydration sphere and stabilizes the native protein form. Direct interactions of this group with surface side chains provide a proper orientation of taurine and prevents the [Formula: see text] group from negative influence. The weakened [Formula: see text] hydration sphere of TMT breaks up the hydrogen-bonded network of water around the protein and destabilizes it. However, TMT at low concentration stabilize both proteins to a small extent. This effect can be attributed to an actual osmophobic effect which is overcome if the concentration increases.

Molecular basis of the osmolyte effect on protein stability: a lesson from the mechanical unfolding of lysozyme

Osmolytes are a class of small organic molecules that shift the protein folding equilibrium. For this reason, they are accumulated by organisms under environmental stress and find applications in biotechnology where proteins need to be stabilized or dissolved. However, despite years of research, debate continues over the exact mechanisms underpinning the stabilizing and denaturing effect of osmolytes. Here, we simulated the mechanical denaturation of lysozyme in different solvent conditions to study the molecular mechanism by which two biologically relevant osmolytes, denaturing (urea) and stabilizing (betaine), affect the folding equilibrium. We found that urea interacts favorably with all types of residues via both hydrogen bonds and dispersion forces, and therefore accumulates in a diffuse solvation shell around the protein. This not only provides an enthalpic stabilization of the unfolded state, but also weakens the hydrophobic effect, as hydrophobic forces promote the association of urea with nonpolar residues, facilitating the unfolding. In contrast, we observed that betaine is excluded from the protein backbone and nonpolar side chains, but is accumulated near the basic residues, yielding a nonuniform distribution of betaine molecules at the protein surface. Spatially resolved solvent–protein interaction energies further suggested that betaine behaves in a ligand- rather than solvent-like manner and its exclusion from the protein surface arises mostly from the scarcity of favorable binding sites. Finally, we found that, in the presence of betaine, the reduced ability of water molecules to solvate the protein results in an additional enthalpic contribution to the betaine-induced stabilization.

Convergent Mechanistic Features between the Structurally Diverse N- and O-Methyltransferases: Glycine N-Methyltransferase and Catechol O-Methyltransferase

Although an enormous and still growing number of biologically diverse methyltransferases have been reported and identified, a comprehensive understanding of the enzymatic methyl transfer mechanism is still lacking. Glycine N-methyltransferase (GNMT), a member of the family that acts on small metabolites as the substrate, catalyzes methyl transfer from S-adenosyl-l-methionine (AdoMet) to glycine to form S-adenosyl-l-homocysteine and sarcosine. We report primary carbon (¹²C/¹⁴C) and secondary (¹H₃/³H₃) kinetic isotope effects at the transferred methyl group, together with ¹H₃/³H₃ binding isotope effects for wild-type GNMT and a series of Tyr21 mutants. The data implicate a compaction effect in the methyl transfer step that is conferred by the protein structure. Furthermore, a remarkable similarity of properties is observed between GNMT and catechol O-methyltransferase, despite significant differences between these enzymes with regard to their active site structures and catalyzed reactions. We attribute these results to a catalytically relevant reduction in the methyl donor–acceptor distance that is dependent on a tyrosine side chain positioned behind the methyl-bearing sulfur of AdoMet.

Study of hydration of sarcosine, formation of its zwitterion and their different oligomers in aqueous media: a Raman spectroscopic and theoretical study

Raman spectra of the biologically important molecule sarcosine (SAR) (C3H7NO2) were studied experimentally in aqueous solution at different concentrations. These spectra were also calculated theoretically using density functional theory (DFT) at the B3LYP/6-311++G(d,p) level. Further, all the observed normal modes were assigned through potential energy distribution (PED). Geometry optimization of SAR produced its three conformers with slightly different energies. The lowest energy conformer of SAR was selected for a systematic solvation study wherein different numbers of water molecules (nW, n=1-9) were placed near it. In the SAR-9W complex, the SAR molecule is located almost at the center of the cage of 9 water molecules. Geometries of different oligomers of SAR (dimer, trimer, tetramer and pentamer) were also optimized in aqueous media taking the input structures from crystallographic data and using the polarizable continuum (PCM). Proton transfer required for the formation of the zwitterionic form of SAR was found to occur when the number of water molecules in the first hydration shell was six or more.

Chemometric method of spectra analysis leading to isolation of lysozyme and CtDNA spectra affected by osmolytes

In this paper we present a chemometric method of analysis leading to isolation of Fourier transform infrared (FT-IR) spectra of biomacromolecules (HEW lysozyme, ctDNA) affected by osmolytes (trimethylamine-N-oxide and N,N,N-trimethylglycine, respectively) in aqueous solutions. The method is based on the difference spectra method primarily used to characterize the structure of solvent affected by solute. The cyclical usage of factor analysis allows precise information to be obtained on the shape of "affected spectra" of analyzed biomacromolecules. "Affected spectra" of selected biomacromolecules give valuable information on their structure in the presence of the osmolytes in solution, as well as on the level of perturbation in dependence of osmolyte concentration. The method also gives a possibility of insight into the mechanism of interaction in presented types of systems. It can be easily adapted to various chemical and biochemical problems where vibrational or ultraviolet-visible (UV-Vis) spectroscopy is used.

Computational investigation on microsolvation of the osmolyte glycine betaine [GB (H(2)O)(1-7)]

The preferential interactions of glycine betaine (GB) with solvent components and the effect of solvent on its stability have been examined. In particular, the microsolvation of organic osmolyte and widely important osmoprotectant in nature as glycine betaine has been reported by using M06 method. A number of configurations (b(X) (a-z)) of the clusters for one to seven water molecules (× = 1-7) have been considered for the microsolvation. Structures of stable conformers are obtained and denoted as b1a, b2a, b3a, b4a, b5a, b6a and b7a. It is observed from the interaction energy difference (∆E) that only seven water molecules can be accommodated in the first solvation shell to stabilize GB. It is also observed that the calculated relative energy using M06 is in close agreement with calculations at the MP2 level of theory.