Structure and dynamics of urea/water mixtures investigated by vibrational spectroscopy and molecular dynamics simulation

A generalized van't Hoff relation for the temperature dependence of complex-valued nonlinear spectra

The temperature dependence of spectra can reveal important insights into the structural and dynamical behavior of the system being probed. In the case of linear spectra, this has been exploited to investigate the thermodynamic driving forces governing the spectral response. Indeed, the temperature derivative of a spectrum can be used to obtain effective energetic and entropic profiles as a function of the measured frequency. The former can further be used to predict the temperature-dependent spectrum via a van't Hoff relation. However, these approaches are not directly applicable to nonlinear, complex-valued spectra, such as vibrational sum-frequency generation (SFG) or two-dimensional infrared (2D-IR) photon echo spectra. Here, we show how the energetic and entropic driving forces governing such nonlinear spectra can be determined and used within a generalized van't Hoff relation to predict their temperature dependence. The central idea is to allow the underlying energetic profiles to themselves be complex-valued. We illustrate this approach for 2D-IR spectra of water and SFG spectra of the air-water interface and demonstrate the accuracy of the generalized van't Hoff relationship and its implications for the origin of temperature-dependent spectral changes.

Development of a universal method for vibrational analysis of the terminal alkyne C≡C stretch

The terminal alkyne C≡C stretch has a large Raman scattering cross section in the "silent" region for biomolecules. This has led to many Raman tag and probe studies using this moiety to study biomolecular systems. A computational investigation of these systems is vital to aid in the interpretation of these results. In this work, we develop a method for computing terminal alkyne vibrational frequencies and isotropic transition polarizabilities that can easily and accurately be applied to any terminal alkyne molecule. We apply the discrete variable representation method to a localized version of the C≡C stretch normal mode. The errors of (1) vibrational localization to the terminal alkyne moiety, (2) anharmonic normal mode isolation, and (3) discretization of the Born-Oppenheimer potential energy surface are quantified and found to be generally small and cancel each other. This results in a method with low error compared to other anharmonic vibrational methods like second-order vibrational perturbation theory and to experiments. Several density functionals are tested using the method, and TPSS-D3, an inexpensive nonempirical density functional with dispersion corrections, is found to perform surprisingly well. Diffuse basis functions are found to be important for the accuracy of computed frequencies. Finally, the computation of vibrational properties like isotropic transition polarizabilities and the universality of the localized normal mode for terminal alkynes are demonstrated.

Counteracting Effects of Trimethylamine N-Oxide against Urea in Aqueous Solutions: Insights from Theoretical Two-Dimensional Infrared Spectroscopy

The study of small osmolytes in their aqueous solutions has gained significant attention because of their relevance to structure and thermodynamics of proteins in aqueous media. Special attention has been given to the binary and ternary aqueous solutions of urea and trimethylamine N-oxide (TMAO). Urea is a well-known protein denaturant, while TMAO protects proteins in their native states. Interestingly, TMAO counteracts urea's ability to denature proteins when present in solutions with approximately half of the concentration of urea. Vibrational spectroscopy can improve our understanding of the molecular origin of this counteracting effect because of its sensitivity to local structure and dynamics. We present results of theoretical linear vibrational and two-dimensional infrared (2DIR) spectroscopy of water in the binary and ternary aqueous solutions of TMAO and urea. The 2DIR spectra are calculated using the electronic structure/molecular dynamics approach. The non-Condon effects in spectral transitions are incorporated in the theoretical calculations of 2DIR spectra. It is found that TMAO disrupts the local structure of water, while urea leaves it essentially unaffected. The 2DIR results show that both TMAO and urea slow down the dynamics of spectral diffusion of water. The extent of slowing down is found to be particularly significant for both hydration and bulk water in the presence of TMAO which can be attributed to strong hydrogen bonds between the water and TMAO molecules. The water molecules present in the hydration layer of the solutes in the ternary solutions are found to relax at even slower rates compared to that in their binary solutions in water. The hydrogen bonds between TMAO and urea are found to be not stable. Thus, the counteracting effect of TMAO against urea is seen to take place mainly through water-mediated interactions. Such TMAO-induced effects giving rise to more structured and slower hydrogen-bonded network are successfully captured through 2DIR spectroscopic calculations.

Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy

The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.

Direct calculation of the temperature dependence of 2D-IR spectra: Urea in water

A method for directly calculating the temperature derivative of two-dimensional infrared (2D-IR) spectra from simulations at a single temperature is presented. The approach is demonstrated by application to the OD stretching spectrum of isotopically dilute aqueous (HOD in H₂O) solutions of urea as a function of concentration. Urea is an important osmolyte because of its ability to denature proteins, which has motivated significant interest in its effect on the structure and dynamics of water. The present results show that the temperature dependence of both the linear IR and 2D-IR spectra, which report on the underlying energetic driving forces, is more sensitive to urea concentration than the spectra themselves. Additional physical insight is provided by calculation of the contributions to the temperature derivative from different interactions, e.g., water-water, water-urea, and urea-urea, present in the system. Finally, it is demonstrated how 2D-IR spectra at other temperatures can be obtained from only room temperature simulations.

Does Urea Preferentially Interact with Amide Moieties or Nonpolar Sidechains? A Question Answered Through a Judicious Selection of Model Systems

The transfer model suggests that urea unfolds proteins mainly by increasing the solubility of the amide backbone, probably through urea-induced increase in hydrogen bonding. Other studies suggest that urea addition increases the magnitude of solvent-solute van der Waals interactions, which increases the solubility of nonpolar sidechains. More recent analyses hypothesize that urea has a similar effect in increasing the solubility of backbone and sidechain groups. In this work, we compare the effects of urea addition on the solvation of amides and alkyl groups. At first, we study the effects of urea addition upon solvent hydrogen bonding acidity and basicity through the perturbation in the fluorescence spectrum of probes 1-AN and 1-DMAN. Our results demonstrate that the solvent's hydrogen bonding properties are minimally affected by urea addition. Subsequently, we show that urea addition does not perturb the intra-molecular hydrogen bonding in salicylic acid significantly. Finally, we investigate how urea preferentially interacts with amide and alkyl groups moieties in water by comparing the effects of urea addition upon the solubility of acetaminophen and 4-tertbutylphenol. We show that urea affects amide and t-butyl solubility (lowers the transfer free energy of both amide (backbone) and alkyl (sidechain) groups) in a similar fashion. In other words, preferential interaction of urea with both moieties contributes to protein denaturation.

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.

Investigating the binding mechanism of topiramate with bovine serum albumin using spectroscopic and computational methods

Various spectroscopic techniques involving fluorescence spectroscopy, circular dichroism (CD), and computational approaches were used to elucidate the molecular aspects of interaction between the antiepileptic drug topiramate and the multifunctional transport protein bovine serum albumin (BSA) under physiological conditions. Topiramate quenched BSA fluorescence in a static quenching mode, according to the Stern-Volmer quenching constant (K_sv ) data derived from fluorescence spectroscopy for the topiramate-BSA complex. The binding constant was also used to calculate the binding affinity for the topiramate-BSA interaction. Fluorescence and circular dichroism experiments demonstrate that the protein's tertiary structure is affected by the microenvironmental alterations generated by topiramate binding to BSA. To establish the exact binding site, interacting residues, and interaction forces involved in the binding of topiramate to BSA, molecular modeling and simulation approaches were used. According to the MMPBSA calculations, the average binding energy between topiramate and BSA is -421.05 kJ/mol. Topiramate was discovered to have substantial interactions with BSA, changing the structural dynamic and Gibbs free energy landscape patterns.

Water-Induced Restructuring of the Surface of a Deep Eutectic Solvent

We study the molecular-scale structure of the surface of Reline, a DES made from urea and choline chloride, using heterodyne-detected vibrational sum frequency generation (HD-VSFG). Reline absorbs water when exposed to the ambient atmosphere, and following structure-specific changes at the Reline/air interface is crucial and difficult. For Reline (dry, 0 wt %, w/w, water) we observe vibrational signatures of both urea and choline ions at the surface. Upon increase of the water content, there is a gradual depletion of urea from the surface, an enhanced alignment, and an enrichment of the surface with choline cations, indicating surface speciation of ChCl. Above 40% w/w water content, choline cations abruptly deplete from the surface, as evidenced by the decrease of the vibrational signal of the -CH2- groups of choline and the rapid rise of a water signal. Above 60% w/w water content, the surface spectrum of aqueous Reline becomes indistinguishable from that of neat water.

Ab Initio Molecular Dynamics Study of Aqueous Solutions of Magnesium and Calcium Nitrates: Hydration Shell Structure, Dynamics and Vibrational Echo Spectroscopy

Ab initio molecular dynamics simulations are performed to study the hydration shell structure, dynamics, and vibrational echo spectroscopy of aqueous Mg(NO₃)₂ and Ca(NO₃)₂ solutions. The hydration shell structure is probed through calculations of various ion-ion and ion-water radial and spatial distribution functions. On the dynamical side, calculations have been made for the hydrogen bond dynamics of hydration shells and also residence dynamics and lifetimes of water in different solvation environments. Subsequently, we looked at the dynamics of frequency fluctuations of OD modes of heavy water in different hydration environments. Specifically, the temporal decay of spectral observables of two-dimensional infrared (2DIR) spectroscopy, three pulse echo peak shift (3PEPS) measurements and also of time correlations of frequency fluctuations are calculated to investigate the dynamics of vibrational spectral diffusion of water in different hydration environments in these solutions. The OD stretch frequencies of water molecules in the vicinity of both divalent cations are found to be red-shifted and also fluctuating at a slower rate than other water molecules present in the solutions. The Mg²⁺ ions are found to be strongly hydrated which can be linked to their lower tendency to form contact ion-pairs and essentially no water exchange between the cationic hydration shells and bulk during the time scale of the current simulations. The stronger hydration of Mg²⁺ ions make their hydration shells structurally and dynamically more rigid and make the dynamics of hydrogen bonds and vibrational spectral diffusion, as revealed through spectral observables of 2DIR and 3PEPS slower than that for the Ca²⁺ ions. The structural and spectral dynamics of water molecules outside the cationic solvation shells in the Mg(NO₃)₂ solution are also found to be relatively slower than that of the Ca(NO₃)₂ solution and pure water which show the effects of stronger electric fields of Mg²⁺ ions extending beyond their first hydration shells. Also, water molecules in the hydration shells of the NO₃^- ions are found to relax at a slower rate in the Mg(NO₃)₂ solution which manifests the effect countercations have on anionic hydration shells for divalent metal nitrate solutions.

Effects of Stearyl Alcohol Monolayer on the Structure, Dynamics and Vibrational Sum Frequency Generation Spectroscopy of Interfacial Water

The structure, dynamics and vibrational spectroscopy of water surface covered by a monolayer of stearyl alcohol (STA) are investigated by means of molecular dynamics simulations and vibrational sum frequency generation...

Direct Observation of the Orientation of Urea Molecules at Charged Interfaces

Dissolving urea into water induces special solvation properties that play a crucial role in many biological processes. Here we probe the properties of urea molecules at charged aqueous interfaces using heterodyne-detected vibrational sum-frequency generation (HD-VSFG) spectroscopy. We find that at the neat water/air interface urea molecules do not yield a significant sum-frequency generation signal. However, upon the addition of ionic surfactants, we observe two vibrational bands at 1660 and 1590 cm^-1 in the HD-VSFG spectrum, assigned to mixed bands of the C═O stretch and NH₂ bend vibrations of urea. The orientation of the urea molecules depends on the sign of the charge localized at surface and closely follows the orientation of the neighboring water molecules. We demonstrate that urea is an excellent probe of the local electric field at aqueous interfaces.

Photophysical Properties of BADAN Revealed in the Study of GGBP Structural Transitions

The fluorescent dye BADAN (6-bromoacetyl-2-dimetylaminonaphtalene) is widely used in various fields of life sciences, however, the photophysical properties of BADAN are not fully understood. The study of the spectral properties of BADAN attached to a number of mutant forms of GGBP, as well as changes in its spectral characteristics during structural changes in proteins, allowed to shed light on the photophysical properties of BADAN. It was shown that spectral properties of BADAN are determined by at least one non-fluorescent and two fluorescent isomers with overlapping absorbing bands. It was found that BADAN fluorescence is determined by the unsolvated "PICT" (planar intramolecular charge transfer state) and solvated "TICT" (twisted intramolecular charge transfer state) excited states. While "TICT" state can be formed both as a result of the "PICT" state solvation and as a result of light absorption by the solvated ground state of the dye. BADAN fluorescence linked to GGBP/H152C apoform is quenched by Trp 183, but this effect is inhibited by glucose intercalation. New details of the changes in the spectral characteristics of BADAN during the unfolding of the protein apo and holoforms have been obtained.

Identical Water Dynamics in Acrylamide Hydrogels, Polymers, and Monomers in Solution: Ultrafast IR Spectroscopy and Molecular Dynamics Simulations

The dynamics and structure of water in polyacrylamide hydrogels (PAAm-HG), polyacrylamide, and acrylamide solutions are investigated using ultrafast infrared experiments on the OD stretch of dilute HOD/H₂O and molecular dynamics simulations. The amide moiety of the monomer/polymers interacts strongly with water through hydrogen bonding (H-bonding). The FT-IR spectra of the three systems indicate that the range of H-bond strengths is relatively unchanged from bulk water. Vibrational population relaxation measurements show that the amide/water H-bonds are somewhat weaker but fall within the range of water/water H-bond strengths. A previous study of water dynamics in PAAm-HG suggested that the slowing observed was due to increasing confinement with concentration. Here, for the same concentrations of the amide moiety, the experimental results demonstrate that the reorientational dynamics (infrared pump-probe experiments) and structural dynamics (two-dimensional infrared spectroscopy) are identical in the three acrylamide systems studied. Molecular dynamics simulations of the water orientational relaxation in aqueous solutions of the acrylamide monomer, trimer, and pentamer are in good agreement with the experimental results and are essentially chain length independent. The simulations show that there is a slower, low-amplitude (<7%) decay component not accessible by the experiments. The simulations examine the dynamics and structure of water H-bonded to acrylamide, in the first solvent shell, and beyond for acrylamide monomers and short chains. The experiments and simulations show that the slowing of water dynamics in PAAm-HG is not caused by confinement in the polymer network but rather by interactions with individual acrylamide moieties.

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.

Remdesivir Strongly Binds to RNA-Dependent RNA Polymerase, Membrane Protein, and Main Protease of SARS-CoV-2: Indication From Molecular Modeling and Simulations

Development of new drugs is a time-taking and expensive process. Comprehensive efforts are being made globally toward the search of therapeutics against SARS-CoV-2. Several drugs such as remdesivir, favipiravir, ritonavir, and lopinavir have been included in the treatment regimen and shown effective results in several cases. Among the existing broad-spectrum antiviral drugs, remdesivir is found to be more effective against SARS-CoV-2. Remdesivir has broad-spectrum antiviral action against many single-stranded RNA viruses including pathogenic SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). In this study, we proposed that remdesivir strongly binds to membrane protein (Mprotein), RNA-dependent RNA polymerase (RDRP), and main protease (Mprotease) of SARS-CoV-2. It might show antiviral activity by inhibiting more than one target. It has been found that remdesivir binds to Mprotease, Mprotein, and RDRP with -7.8, -7.4, and -7.1 kcal/mol, respectively. The structure dynamics study suggested that binding of remdesivir leads to unfolding of RDRP. It has been found that strong binding of remdesivir to Mprotein leads to decrease in structural deviations and gyrations. Additionally, the average solvent-accessible surface area of Mprotein decreases from 127.17 to 112.12 nm², respectively. Furthermore, the eigenvalues and the trace of the covariance matrix were found to be low in case of Mprotease-remdesivir, Mprotein-remdesivir, and RDRP-remdesivir. Binding of remdesivir to Mprotease, Mprotein, and RDRP reduces the average motions in protein due to its strong binding. The MMPBSA calculations also suggested that remdesivir has strong binding affinity with Mprotein, Mprotease, and RDRP. The detailed analysis suggested that remdesivir has more than one target of SARS-CoV-2.

Should deep eutectic solvents be treated as a mixture of two components or as a pseudo-component?

Deep eutectic solvents (DESs) and dilutions thereof (mainly in H₂O but also in many other non-aqueous solvents and co-solvent mixtures) have recently attracted great attention. It is well known that DES dilutions exhibit deviations from ideality. Interestingly, the treatment of DES as a mixture of two components or a pseudo-component is by no means trivial when determining deviations in density and, mainly, in viscosity. Herein, we studied aqueous dilutions of one of the most widely studied DES, this is, that composed of choline chloride and urea in a 1:2 molar ratio (e.g., ChCl2U). Using density and viscosity data reported in previous works, we calculated the excess molar volumes (V^E) and excess viscosities (ln η^E) considering ChCl2U as either a mixture of two components or a pseudo-component, that is, taking the DES molecular weight as M_ChCl2U = f_ChClM_ChCl + f_UM_U = 86.58 g mol^-1 (with f_ChCl = 1/3 and f_U = 2/3) or as M^* _ChCl2U = M_ChCl + 2 M_U = 259.74 g mol^-1. We found that neither the sign of V^E and V^E* nor their evolution with temperature was influenced by the use of either M_ChCl2U or M^* _ChCl2U, and only the absolute magnitude of the deviation and the DES content (in wt. %) at which the minimum appears exhibited some differences. However, ln η^E and ln η^E* exhibited opposite signs, negative and positive, respectively. The odd achievement of negative ln η^E in aqueous dilutions of ChCl2U characterized by the formation of HB networks suggest the treatment of ChCl2U as a pseudo-component as more appropriate. Moreover, the role played by the presence of U in the evolution of ln η^E* with temperature was also discussed.

Structure and Interactions in Perfluorooctanoate Micellar Solutions Revealed by Small-Angle Neutron Scattering and Molecular Dynamics Simulations Studies: Effect of Urea

The self-assembly of surfactants in aqueous solution can be modulated by the presence of additives including urea, which is a well-known protein denaturant and also present in physiological fluids and agricultural runoff. This study addresses the effects of urea on the structure of micelles formed in water by the fluorinated surfactant perfluoro-n-octanoic acid ammonium salt (PFOA). Analysis of small-angle neutron scattering (SANS) experiments and atomistic molecular dynamics (MD) simulations provide consensus strong evidence for the direct mechanism of urea action on micellization: urea helps solvate the hydrophobic micelle core by localizing at the surface of the core in the place of some water molecules. Consequently, urea decreases electrostatic interactions at the micelle shell, changes the micelle shape from prolate ellipsoid to sphere, and decreases the number of surfactant molecules associating in a micelle. These findings inform the interactions and behavior of surface active per- and polyfluoroalkyl substances (PFAS) released in the aqueous environment and biota.

Pyrazoline analogs as potential anticancer agents and their apoptosis, molecular docking, MD simulation, DNA binding and antioxidant studies

N-formyl pyrazoline derivatives (3a-3l) were designed and synthesized via Michael addition reaction through cyclization of chalcones with hydrazine hydrate in presence of formic acid. The structural elucidation of N-formyl pyrazoline derivatives was carried out by various spectroscopic techniques such as ¹H, ¹³C NMR, FT-IR, UV-visible spectroscopy, mass spectrometry and elemental analysis. Anticancer activity of the pyrazoline derivatives (3a-3l) was evaluated against human lung cancer (A549), fibrosarcoma cell lines (HT1080) and human primary normal lung cells (HFL-1) by MTT assay. The results of anticancer activity showed that potent analogs 3b and 3d exhibited promising activity against A549 (IC₅₀ = 12.47 ± 1.08 and 14.46 ± 2.76 µM) and HT1080 (IC₅₀ = 11.40 ± 0.66 and 23.74 ± 13.30 µM) but low toxic against the HFL-1 (IC₅₀ = 116.47 ± 43.38 and 152.36 ± 22.18 µM). The anticancer activity of potent derivatives (3b and 3d) against A549 cancer cell line was further confirmed by flow cytometry based approach. DNA binding interactions of the pyrazoline derivatives 3b and 3d have been carried out with calf thymus DNA (Ct-DNA) using absorption, fluorescence and viscosity measurements, circular dichroism and cyclic voltammetry. Antioxidant potential of N-formyl pyrazoline derivatives (3a-3l) has been also estimated through DPPH (2,2-diphenyl-1-picrylhydrazyl) free radical and H₂O_2. Results revealed that all the compounds exhibited significant antioxidant activity. In silico molecular modelling and ADMET properties of pyrazoline derivatives were also studied using PyRx software against topoisomerase II receptor with PDB ID: 1ZXM to explore their best hits. MD simulation of 3b and 3d was also carried out with topoisomerase II for structure-function correlation in a protein. HuTopoII inhibitory activity of the analogs (3a-3l) was examined by relaxation assay at varying concentrations 100-1000 µM.

Solvatochromism in urea/water and urea-derivative/water solutions

This work reports the experimental measurements of solvent acidity (SA), basicity (SB), and solvent dipolarity and polarizability (SPP) for water solutions with urea (U) and its molecular derivatives, monomethyl-urea (MU), 1,3-dimethyl-urea (DMU) and tetramethyl-urea (TMU). These solvatochromic parameters are applied to understanding the variation of indexes of refraction and densities and other physico-chemical properties reported for these solutions. These properties are well correlated to the SA, SB, and SPP solvent parameters of these solutions. As a result, from the characterization of the physico-chemical properties, one can infer that urea and its molecular derivatives are mainly modifiers in the structure of liquid water. The solvatochromic parameters indicate the possible existence of different mechanisms in the denaturation process of proteins in these urea/water solutions.

Two-Dimensional Infrared Spectroscopy of Aqueous Solutions of Metal Nitrates: Slowdown of Spectral Diffusion in the Presence of Divalent Cations

The hydrogen-bonded network of water can be affected both structurally and dynamically by the presence of ions. In the present study, we have considered three aqueous solutions of metal nitrates to investigate the effects of divalent cations (Mg²⁺ and Ca²⁺), compared to that of monovalent Na⁺ ions, on hydrogen-bond fluctuations and vibrational spectral diffusion through calculations of linear and two-dimensional infrared spectra of these solutions at room temperature. We have employed the methods of molecular dynamics simulations using effective polarizable models of ions combined with quantum mechanical calculations of transition variables and statistical mechanical calculations of spectral response functions of vibrational spectroscopy. Divalent cations are found to have much stronger and longer-ranged effects on the structure and dynamics of the hydrogen-bonded network than that induced by the monovalent sodium ions. The blue shifts in the calculated linear spectra are found to follow the Hofmeister trend for the cations. The 2D-IR spectral lineshape and intensity corresponding to three-pulse echo peak shift (3PEPS) experiments are calculated. The timescales of these nonlinear spectral responses and also frequency-time correlations show significant slowing down of spectral diffusion for solutions containing divalent Mg²⁺ and Ca²⁺ ions compared to the corresponding dynamics of the solution containing Na⁺ ions. Unlike NaNO₃ solution, the relaxation of frequency and dipole orientational fluctuations of anion-bound water in Mg(NO₃)₂ and Ca(NO₃)₂ solutions are found to be somewhat slower than bulk water, which can be attributed to the presence of divalent cations whose effects go beyond their first solvation shells. This is also seen in the dynamics of bulk water in these solutions which is found to be notably slower for the solutions containing divalent cations than that in the NaNO₃ solution. Unlike Mg²⁺ and Ca²⁺ ions, no specific cationic effect is observed for the Na⁺ ions.

Mechanistic insights into the urea-induced denaturation of human sphingosine kinase 1

Sphingosine kinase 1 (SphK1) plays a significant role in various cellular processes, including cell proliferation, apoptosis, and angiogenesis. SphK1 is considered as an attractive target for drug development owing to its connection with several diseases, including cancer. In the current work, the urea-induced unfolding of SphK1 was performed at pH 8.0 and 25 °C using CD and fluorescence spectroscopy. SphK1 follows a biphasic unfolding transition (N ⇌ I ⇌ D) with an intermediate (I) state populated around 4.0 M urea concentration. The circular dichroism ([θ]₂₂₂) and fluorescence emission spectra (λ_max) of SphK1 with increasing concentrations of urea were analyzed to calculate Gibbs free energy (ΔG⁰) for both the transitions (N ⇌ I and I ⇌ D). A significant overlap of both the transitions obtained by two spectroscopic properties ([θ]₂₂₂ and λ_max) was observed, indicating that both N ⇌ I and I ⇌ D transition follow two-step equilibrium unfolding pattern. Also, we performed 100 ns molecular dynamics (MD) simulations to get atomistic insights into the structural changes in SphK1 with increasing urea concentrations. Our results showed a consistent pattern of the SphK1 unfolding with increasing urea concentrations. Together, spectroscopic and MD simulation findings provide deep insights into the unfolding mechanism and conformational features of SphK1.

Mechanistic insights into the urea-induced denaturation of a non-seleno thiol specific antioxidant human peroxiredoxin 6

Peroxiredoxin 6 (Prdx6) is a unique enzyme among mammalian peroxiredoxins as it lacks resolving cysteine. It is found to be involved in number of different diseases including tumours and its expression level is highest in lungs as compared to other organs. It has been found that Prdx6 plays a significant role different metabolic diseases, ocular damage, neurodegeneration and male infertility. It is a bifunctional protein having phospholipase A₂ and peroxidase (also has the ability to reduce phospholipid hydroperoxides) activities. In order to complete the peroxidise reaction cycle it requires glutathione catalyzed by glutathione S-transferase. Equilibrium unfolding and conformational stability of Prdx6 was studied by using urea as a chemical denaturant to understand the changes it goes under cellular stress conditions. Three different spectroscopic methods were employed to monitor urea-induced denaturation. From the results obtained, it was found that the urea denaturation of Prdx6 follows a variable two state process due to non-coincidence of the normalized transition curves obtained from different optical probes. The different denaturation curves were normalized and thermodynamic parameters, ΔG_D^o, Gibbs free energy change related to the urea-induced denaturation, midpoint of denaturation (C_m), and m = (δΔG_D / [urea]) were obtained. The structural information of Prdx6 were further analysed by several parameters obtained by 100 ns MD simulation. The results of MD simulation clearly favour the outcome of spectroscopic studies.

Hydrophobic pattern of alkylated ureas markedly affects water rotation and hydrogen bond dynamics in aqueous solution

Alkylated ureas are frequently used amphiphiles to mediate biomolecule water interactions, yet their hydrophobic substitution pattern critically affects their function. These differences can be traced back to their hydration, which is poorly understood. Here, we investigate subtle effects of the hydrophobic pattern of ureas on hydration dynamics using a combination of linear and non-linear infrared spectroscopies on the OD stretching vibration of HDO. Isomeric 1,3-dimethylurea (1,3-DMU), 1,1-dimethylurea (1,1-DMU) and 1-ethylurea (1-EU) exhibit very similar and rather weak modulation of the water hydrogen-bond strength distribution. Yet, only 1,3-DMU and 1,1-DMU enhance the hydrogen-bond heterogeneity and slow-down its fluctuation dynamics. In turn, rotational dynamics of water molecules, which is dominated by hydrogen bond switches, is significantly impeded in the presence of 1,3-DMU and only weakly by 1,1-DMU and 1-EU. These marked differences can be explained by both excluded volume effects in hydration and self-aggregation, which may be the key to their biotechnological function.

Vibrational echo spectroscopy of aqueous sodium bromide solutions from first principles simulations

A theoretical study of the time-dependent vibrational echo spectroscopy of sodium bromide solutions in deuterated water at two different concentrations of 0.5 and 5.0 M and at temperatures of 300 and 350 K is presented using the method of ab initio molecular dynamics simulations. The instantaneous fluctuations in frequencies of local OD stretch modes are calculated using time-series analysis of the simulated trajectories. The third-order polarization and intensities of three pulse photon-echo are calculated from ab initio simulations. The timescales of vibrational spectral diffusion are determined from the frequency time correlation functions (FTCF) and short-time slope of three pulse photon echo (S3PE) calculated within the second-order cumulant and Condon approximations. It is found that under ambient conditions, the rate of vibrational spectral diffusion becomes slower with increase in ionic concentration. Decay of S3PE calculated for different systems give timescales, which are in close agreement with those of FTCF and also with the results of experimental time-dependent vibrational spectroscopic experiments. © 2019 Wiley Periodicals, Inc.

Permeability of P and K-nutrient through polystyrene membrane from aqueous solutions of urea + KH2PO4

With the polymer-coated fertilizer as background, the permeability of P- and K-nutrient through a representative polymer membrane-polystyrene membrane were investigated by measuring their permeability in the solutions of KH2PO4-water and urea-KH2PO4-water at nominal temperature of 298 K using the Ussing chamber method. To analyze and interpret the variation of permeability with solute concentration, the solubility of permeate in polymer membrane were determined experimentally and the permeate diffusion coefficient were assessed by the measurements of density and apparent molar volume of the aqueous fertilizer solutions. An interesting “increase-decrease” trend for the permeability of both phosphorous (P)-nutrient, and potassium (K)-nutrient fertilizer with permeate concentration was observed, in which the increases in permeability at low concentrations of permeate could be attributed to the increase in solubility of KH2PO4 in polymer while the decreases in permeability at high concentrations was due to the decrease in diffusion coefficient of permeate in polymer membrane. Finally, the release kinetics of these nutrients from a PS-coated urea-KH2PO4 compound fertilizer granule was predicted using the Shaviv’s model along with the permeability data of P- and K-nutrient generated.

Temperature dependence of the ultrafast vibrational echo spectroscopy of OD modes in liquid water from first principles simulations

The temperature dependence of the vibrational spectral diffusion of OD modes in liquid water is investigated through calculations of vibrational echo spectral observables from first principles molecular dynamics.

Contrasting Effects of Guanidinium Chloride and Urea on the Activity and Unfolding of Lysozyme

Cosolvents play an important role in regulating the stability and function of proteins present in the cell. We studied the role of cosolvents, urea and guanidinium chloride (GdmCl), which act as protein denaturants, in the catalytic activity and structural stability of the protein lysozyme using activity measurements, spectroscopy, and molecular dynamics simulations. We find that the activity of lysozyme increases on the addition of urea, whereas it decreases sharply on the addition of GdmCl. At low GdmCl concentrations ([GdmCl] < 4 M), the activity of lysozyme decreases, even though there is no significant perturbation in the structure of the lysozyme folded state. We find that this is due to the strong interaction of the Gdm⁺ ion with the residues Asp52 and Glu35, which are present in the lysozyme catalytic site. In contrast, urea interacts with Trp63 present in the loop region present near the active site of lysozyme, inducing minor conformational changes in lysozyme, which can increase the activity of lysozyme. At higher denaturant concentrations, experiments show that GdmCl completely denatures the protein, whereas the folded state is stable in the presence of urea. We further show that GdmCl denatures lysozyme with the disulfide bonds intact in the protein, whereas urea denatures the protein only when the disulfide bonds are broken using reducing agents.

Noncovalent Interactions between Trimethylamine N-Oxide (TMAO), Urea, and Water

Trimethylamine N-oxide (TMAO) and urea are two important osmolytes with their main significance to the biophysical field being in how they uniquely interact with proteins. Urea is a strong protein destabilizing agent, whereas TMAO is known to counteract urea's deleterious effects. The exact mechanisms by which TMAO stabilizes and urea destabilizes folded proteins continue to be debated in the literature. Although recent evidence has suggested that urea binds directly to amino acid side chains to make protein folding less thermodynamically favored, it has also been suggested that urea acts indirectly to denature proteins by destabilizing the surrounding hydrogen bonding water networks. Here, we elucidate the molecular level mechanism of TMAO's unique ability to counteract urea's destabilizing nature by comparing Raman spectroscopic frequency shifts to the results of electronic structure calculations of microsolvated molecular clusters. Experimental and computational data suggest that the addition of TMAO into an aqueous solution of urea induces blue shifts in urea's H-N-H symmetric bending modes, which is evidence for direct interactions between the two cosolvents.

Unravelling the unfolding mechanism of human integrin linked kinase by GdmCl-induced denaturation

Integrin-linked kinase (ILK) is a ubiquitously expressed Ser/Thr kinase which plays significant role in the cell-matrix interactions and growth factor signalling. In this study, guanidinium chloride (GdmCl)-induced unfolding of kinase domain of ILK (ILK_193-446) was carried out at pH 7.5 and 25 °C. Eventually, denaturation curves of mean residue ellipticity at 222 nm ([θ]₂₂₂) and fluorescence emission spectrum were analysed to estimate stability parameters. The optical properties maximum emission (λ_max) and difference absorption coefficient at 292 nm (Δε₂₉₂) were analysed. The denaturation curve was measured only in the GdmCl molar concentration ranging 3.0-4.2 M because protein was aggregating below 3.0 M of GdmCl concentrations. The denaturation process of ILK_193-446 was found as reversible at [GdmCl] ≥ 3.0 M. Moreover, a coincidence of normalized denaturation curves of optical properties ([θ]₂₂₂, Δε₂₉₂ and λ_max) suggesting that GdmCl-induced denaturation of ILK_193-446 is a two-state process. In addition, 100 ns molecular dynamics simulations were performed to see the effects of GdmCl on the structure and stability of ILK_193-446. Both the spectroscopic and molecular dynamics approaches provided clear insights into the stability and conformational properties of ILK_193-446.