A Molecular Dynamics and Quantum Mechanics Analysis of the Effect of DMSO on Enzyme Structure and Dynamics: Subtilisin

Solvent-Modulated Specific Ion Effects: Poly(N-isopropylacrylamide) Brushes in Nonaqueous Electrolytes

Pertinent to cryopreservation as well as energy storage and batteries, nonaqueous electrolytes and their mixtures with water were investigated. In particular, specific ion-induced effects on the modulation of a poly(N-isopropylacrylamide) (PNIPAM) brush were investigated in various dimethyl sulfoxide (DMSO)-water solvent mixtures. Spectroscopic ellipsometry and neutron reflectometry were employed to probe changes in brush swelling and structure, respectively. In water-rich solvents (i.e., pure water and 6 mol % DMSO), PNIPAM undergoes a swollen to collapsed thermotransition with increasing temperature, whereby a forward Hofmeister series was noted; K+ and Li+ electrolytes composed of SCN- and I- salted-in (stabilized) PNIPAM chains, and electrolytes of Cl- and Br- salted-out (destabilized) the polymer. The cation was seen to play a lesser role than that of the anion, merely modulating the magnitude of the anion effect. In 70 mol % DMSO, a collapsed to swollen thermotransition was noted for PNIPAM. Here, concentration-dependent specific ion effects were observed; a forward series was observed in 0.2 mol % electrolytes, whereas increasing the electrolyte concentration to 0.9 mol % led to a series reversal. While no thermotransition was observed in pure DMSO, a solvent-induced specific ion series reversal was noted; SCN- destabilized the brush and Cl- stabilized the brush. Both series reversals are attributed to the delicate balance of interactions between the solvent, solute (ion), and substrate (brush). Namely, the stability of the solvent clusters was hypothesized to drive polymer solvation.

Electrostatic interactions of enzymes in non-aqueous conditions: insights from molecular dynamics simulations

Electrostatic interactions of enzymes and their effects on enzyme activity and stability are poorly understood in non-aqueous conditions. Here, we investigate the contribution of the electrostatic interactions on the stability and activity of enzymes in the non-aqueous environment using molecular dynamics simulations. Lipase was selected as active and lysozyme as inactive model enzymes in non-aqueous media. Hexane was used as a common non-aqueous solvent model. In agreement with the previous experiments, simulations show that lysozyme has more structural instabilities than lipase in hexane. The number of hydrogen bonds and salt bridges of both enzymes is dramatically increased in hexane. In contrast to the other opinions, we show that the increase of the electrostatic interactions in non-aqueous media is not so favorable for enzymatic function and stability. In this condition, the newly formed hydrogen bonds and salt bridges can partially denature the local structure of the enzymes. For lysozyme, the changes in electrostatic interactions occur in all domains including the active site cleft, which leads to enzyme inactivation and destabilization. Interestingly, most of the changes in electrostatic interactions of lipase occur far from the active site regions. Therefore, the active site entrance regions remain functional in hexane. The results of this study reveal how the changes in electrostatic interactions can affect enzyme stability and activity in non-aqueous conditions. Moreover, we show for the first time how some enzymes, such as lipase, remain active in a non-aqueous environment.Communicated by Ramaswamy H. Sarma.

Weak Interactions in Dimethyl Sulfoxide (DMSO)-Tertiary Amide Solutions: The Versatility of DMSO as a Solvent

The structures of equimolar mixtures of the commonly used polar aprotic solvents dimethylformamide (DMF) and dimethylacetamide (DMAc) in dimethyl sulfoxide (DMSO) have been investigated via neutron diffraction augmented by extensive hydrogen/deuterium isotopic substitution. Detailed 3-dimensional structural models of these solutions have been derived from the neutron data via Empirical Potential Structure Refinement (EPSR). The intermolecular center-of-mass (CoM) distributions show that the first coordination shell of the amides comprises ∼13-14 neighbors, of which approximately half are DMSO. In spite of this near ideal coordination shell mixing, the changes to the amide-amide structure are found to be relatively subtle when compared to the pure liquids. Analysis of specific intermolecular atom-atom correlations allows quantitative interpretation of the competition between weak interactions in the solution. We find a hierarchy of formic and methyl C-H···O hydrogen bonds forms the dominant local motifs, with peak positions in the range of 2.5-3.0 Å. We also observe a rich variety of steric and dispersion interactions, including those involving the O═C-N amide π-backbones. This detailed insight into the structural landscape of these important liquids demonstrates the versatility of DMSO as a solvent and the remarkable sensitivity of neutron diffraction, which is critical for understanding weak intermolecular interactions at the nanoscale and thereby tailoring solvent properties to specific applications.

Cosolvent effects on the structure and thermoresponse of a polymer brush: PNIPAM in DMSO–water mixtures

Structural characterisation of thermoresponsive polymer brushes in binary DMSO–water mixtures reveals both LCST and UCST behaviour.

Cryopreservation and Cryobanking of Cells from 100 Coral Species

When coral species become extinct, their genetic resources cannot be recovered. Coral cryobanks can be employed to preserve coral samples and thereby maintain the availability of the samples and increase their potential to be restocked. In this study, we developed a procedure to determine coral species-specific requirements for cryobank freezing through determining suitable cryoprotective agents (CPAs), CPA concentrations, equilibration times, holding durations, viability rates, and cell amounts for banked coral cells, and we established the first ever coral cell cryobank. Coral cells, including supporting and gland cells, epidermal nematocysts, Symbiodiniaceae and symbiotic endoderm cells (SEC) were found from the extracted protocol. Approximately half of the corals from the experimental corals consisted of spindle and cluster cells. Gastrodermal nematocysts were the least common. The overall concentration of Symbiodiniaceae in the coral cells was 8.6%. Freezing using DMSO as a CPA was suitable for approximately half of the corals, and for the other half of species, successful cell cryopreservation was achieved using MeOH and EG. EG and DMSO had similar suitabilities for Acanthastrea, Euphyllia, Favites, Lobophyllia, Pavona, Seriatopora, and Turbinaria, as did EG and MeOH for Acropora, Echinopyllia, and Sinularia and MeOH and DMSO for Platygyra after freezing. At least 14 straws from each species of coral were cryobanked in this study, totaling more than 1884 straws (0.5 mL) with an average concentration of 6.4 × 10⁶ per mL. The results of this study may serve as a framework for cryobanks worldwide and contribute to the long-term conservation of coral reefs.

Unraveling the acid-base characterization and solvent effects on the structural and electronic properties of a bis-bidentate bridging ligand

Understanding the interactions and the solvent effects on the distribution of several species in equilibrium and how it can influence the ¹H-NMR properties, spectroscopy (UV-vis absorption), and the acid-base equilibria can be especially challenging. This is the case of a bis-bidentate bridging ligand bis(2-pyridyl)-benzo-bis(imidazole), where the two pyridyl and four imidazolyl nitrogen atoms can be protonated in different ways, depending on the solvent, generating many isomeric/tautomeric species. Herein, we report a combined theoretical-experimental approach based on a sequential quantum mechanics/molecular mechanics procedure that was successfully applied to describe in detail the acid-base characterization and its effects on the electronic properties of such a molecule in solution. The calculated free-energies allowed the identification of the main species present in solution as a function of the solvent polarity, and its effects on the magnetic shielding of protons (¹H-NMR chemical shifts), the UV-vis absorption spectra, and the acid-base equilibrium constants (pK_as) in aqueous solution. Three acid-base equilibrium constants were experimentally/theoretically determined (pK_a1 = 1.3/1.2, pK_a2 = 2.1/2.2 and pK_a5 = 10.1/11.3) involving mono-deprotonated and mono-protonated cis and trans species. Interestingly, other processes with pK_a3 = 3.7 and pK_a4 = 6.0 were also experimentally determined and assigned to the protonation and deprotonation of dimeric species. The dimerization of the most stable neutral species was investigated by Monte Carlo simulations and its electronic effects were considered for the elucidation of the UV-vis absorption bands, revealing transitions mainly with the charge-transfer characteristic and involving both the monomeric species and the dimeric species. The good matching of the theoretical and experimental results provides an atomistic insight into the solvent effects on the electronic properties of this bis-bidentate bridging ligand.

Anion Effects on the Mixing States of 1-Methyl-3-octylimidazolium Tetrafluoroborate and Bis(trifluoromethylsulfonyl)amide with Methanol, Acetonitrile, and Dimethyl Sulfoxide on the Meso- and Microscopic Scales

The mixing states of two imidazolium-based ionic liquids (ILs) with different anions, 1-methyl-3-octylimidazolium tetrafluoroborate (C₈mimBF₄) and bis(trifluoromethylsulfonyl)amide (C₈mimTFSA), with three molecular liquids (MLs), methanol (MeOH), acetonitrile (AN), and dimethyl sulfoxide (DMSO), have been investigated on both mesoscopic and microscopic scales using small-angle neutron scattering (SANS), infrared (IR), and ¹H and ¹³C nuclear magnetic resonance (NMR) spectroscopy. Additionally, molecular dynamics (MD) simulations have been conducted on the six combinations of ILs and MLs to observe the states of their mixtures on the atomic level. The SANS profiles of the IL-ML mixtures suggested that MeOH molecules only form clusters in both C₈mimBF₄ and C₈mimTFSA, whereas AN and DMSO were homogeneously mixed with ILs on the SANS scale. MeOH clusters are more enhanced in BF₄^--IL than TFSA^--IL. The microscopic interactions among IL cations, anions, and MLs should contribute to the mesoscopic mixing states of the IL-ML mixtures. In fact, the IL cation-anion, cation-ML, anion-ML, and ML-ML interactions observed by IR, NMR, and MD simulations clarified the reasons for the mixing states of the IL-ML binary solutions observed by the SANS experiments. In neat ILs, the imidazolium ring of the IL cation more strongly interacts with BF₄^- than TFSA^- due to the higher charge density of the former. The interaction of anions with the imidazolium ring is more easily loosened on adding MLs to ILs in the order of DMSO > MeOH > AN. It does not significantly depend on the anions. However, the replacement of the anion on the imidazolium ring by an ML depends on the anions; the replacement is more proceeded in the order of MeOH > DMSO > AN in BF₄^--IL, while DMSO > MeOH > AN in TFSA^--IL. On the other hand, the solvation of both anions by MLs is stronger in the order of MeOH > DMSO ≈ AN. Despite the stronger interactions of MeOH with both cations and anions, MeOH molecules are heterogeneously mixed with both ILs to form clusters in the mixtures. Therefore, the self-hydrogen bonding among MeOH molecules most markedly governs the mixing state of the binary solutions among the abovementioned interactions.

The pursuit of resin-dentin bond durability: Simultaneous enhancement of collagen structure and polymer network formation in hybrid layers

OBJECTIVE

Imperfect polymer formation as well as collagen's susceptibility to enzymatic degradation increase the vulnerability of hybrid layers over time. This study investigated the effect of new dimethyl sulfoxide (DMSO)-containing pretreatments on long-term bond strength, hybrid layer quality, monomer conversion and collagen structure.

METHODS

H₃PO₄-etched mid-coronal dentin surfaces from extracted human molars (n = 8) were randomly treated with aqueous and ethanolic DMSO solutions or following the ethanol-wet bonding technique. Dentin bonding was performed with a three-step etch-and-rinse adhesive. Resin-dentin beams (0.8 mm²) were stored in artificial saliva at 37 °C for 24 h and 2.5 years, submitted to microtensile bond strength testing at 0.5 mm/min and semi-quantitative SEM nanoleakage analysis (n = 8). Micro-Raman spectroscopy was used to determine the degree of conversion at different depths in the hybrid layer (n = 6). Changes in the apparent modulus of elasticity of demineralized collagen beams measuring 0.5 × 1.7 × 7 mm (n = 10) and loss of dry mass (n = 10) after 30 days were calculated via three-point bending and precision weighing, respectively.

RESULTS

DMSO-containing pretreatments produced higher bond strengths, which did not change significantly over time presenting lower incidence of water-filled zones. Higher uniformity in monomer conversion across the hybrid layer occurred for all pretreatments. DMSO-induced collagen stiffening was reversible in water, but with lower peptide solubilization.

SIGNIFICANCE

Improved polymer formation and higher stability of the collagen-structure can be attributed to DMSO's unique ability to simultaneously modify both biological and resin components within the hybrid layer. Pretreatments composed of DMSO/ethanol may be a viable-effective alternative to extend the longevity of resin-dentin bonds.

Application of the 3D-RISM-KH molecular solvation theory for DMSO as solvent

The molecular solvation theory in the form of the Three-Dimensional Reference Interaction Site Model (3D-RISM) with Kovalenko-Hirata (KH) closure relation is benchmarked for use with dimethyl sulfoxide (DMSO) as solvent for (bio)-chemical simulation within the framework of integral equation formalism. Several force field parameters have been tested to correctly reproduce solvation free energy in DMSO, ion solvation in DMSO, and DMSO coordination prediction. Our findings establish a united atom (UA) type parameterization as the best model of DMSO for use in 3D-RISM-KH theory based calculations.

Directed sortase A evolution for efficient site-specific bioconjugations in organic co-solvents

Directed sortase A evolution yielded the variants R159G and D165Q/D186G/K196V with increased resistance (2.2-fold) and catalytic efficiency (6.3-fold) in 45% (v/v) dimethylsulfoxide. Interestingly, D165Q/D186G/K196V also showed an up to 4.7-fold increased activity for the conjugation of hydrophobic peptides/amines in co-solvents. MD simulations revealed that conformational mobilities are important for the gained resistance.

Molecular dynamics study of unfolding of lysozyme in water and its mixtures with dimethyl sulfoxide

All-atom explicit solvent molecular dynamics was used to study the process of unfolding of hen egg white lysozyme in water and mixtures of water with dimethyl sulfoxide at different compositions. We have determined the kinetic parameters of unfolding at a constant temperature 450K. For each run, the time of disruption of the tertiary structure of lysozyme t_u was defined as the moment when a certain structural criterion computed from the trajectory reaches its critical value. A good agreement is observed between the results obtained using several different criteria. The secondary structure according to DSSP calculations is found to be partially unfolded to the moment of disruption of tertiary structure, but some of its elements keep for a long time after that. The values of t_u averaged over ten 30ns-long trajectories for each solvent composition are shown to decrease very rapidly with addition of dimethyl sulfoxide, and rather small amounts of dimethyl sulfoxide are found to change the pathway of unfolding. In pure water, despite the loss of tertiary contacts and disruption of secondary structure elements, the protein preserves its compact globular state at least over 130ns of simulation, while even at 5mol percents of dimethyl sulfoxide it loses its compactness within 30ns. The proposed methodology is a generally applicable tool to quantify the rate of protein unfolding in simulation studies.

Scrutiny of electrostatic-driven conformational ordering of polypeptide chains in DMSO: a study with a model oligopeptide

The molecular mechanism of DMSO-induced stabilisation of β-sheets is attributed to the combination of polar electrostatic interactions among side chains, and backbone desolvation through bulky side chains which promotes backbone hydrogen bonding.

DMSO enhanced conformational switch of an interfacial enzyme

Interfacial proteins function in unique heterogeneous solvent environments, such as water-oil interfaces. One important example is microbial lipase, which is activated in an oil-water emulsion phase and has many important enzymatic functions. A unique aprotic dipolar organic solvent, dimethyl sulfoxide (DMSO), has been shown to increase the activity of lipases, but the mechanism behind this enhancement is still unknown. Here, all-atom molecular dynamics simulations of lipase in a binary solution were performed to examine the effects of DMSO on the dynamics of the gating mechanism. The amphiphilic α5 region of the lipase was a focal point for the analysis, since the structural ordering of α5 has been shown to be important for gating under other perturbations. Compared to the closed-gorge ensemble in an aqueous environment, the conformational ensemble shifts towards open-gorge structures in the presence of DMSO solvents. Increased width of the access channel is particularly prevalent in 45% and 60% DMSO concentrations (w/w). As the amount of DMSO increases, the α5 region of the lipase becomes more α-helical, as we previously observed in studies that address water-oil interfacial and high pressure activation. We believe that the structural ordering of α5 plays an essential role on gating and lipase activity.

Excitation Dynamics in Hetero-bichromophoric Calixarene Systems

In this work, the dynamics of electronic energy transfer (EET) in bichromophoric donor-acceptor systems, obtained by functionalizing a calix[4]arene scaffold with two dyes, was experimentally and theoretically characterized. The investigated compounds are highly versatile, due to the possibility of linking the dye molecules to the cone or partial cone structure of the calix[4]arene, which directs the two active units to the same or opposite side of the scaffold, respectively. The dynamics and efficiency of the EET process between the donor and acceptor units was investigated and discussed through a combined experimental and theoretical approach, involving ultrafast pump-probe spectroscopy and density functional theory based characterization of the energetic and spectroscopic properties of the system. Our results suggest that the external medium strongly determines the particular conformation adopted by the bichromophores, with a direct effect on the extent of excitonic coupling between the dyes and hence on the dynamics of the EET process itself.

To Keep or Not to Keep? The Question of Crystallographic Waters for Enzyme Simulations in Organic Solvent

The use of enzymes in non-aqueous solvents expands the use of biocatalysts to hydrophobic substrates, with the ability to tune selectivity of reactions through solvent selection. Non-aqueous enzymology also allows for fundamental studies on the role of water and other solvents in enzyme structure, dynamics, and function. Molecular dynamics simulations serve as a powerful tool in this area, providing detailed atomic information about the effect of solvents on enzyme properties. However, a common protocol for non-aqueous enzyme simulations does not exist. If you want to simulate enzymes in non-aqueous solutions, how many and which crystallographic waters do you keep? In the present work, this question is addressed by determining which crystallographic water molecules lead most quickly to an equilibrated protein structure. Five different methods of selecting and keeping crystallographic waters are used in order to discover which crystallographic waters lead the protein structure to reach an equilibrated structure more rapidly in organic solutions. It is found that buried waters contribute most to rapid equilibration in organic solvent, with slow-diffusing waters giving similar results.

(1)H NMR spectra of alcohols in hydrogen bonding solvents: DFT/GIAO calculations of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in hydrogen bonding solvents have been computed on the basis of density functional theory by applying the gauge-including atomic orbital method to geometry-optimized alcohol/solvent complexes. The OH proton shifts and hydrogen bond distances for methanol or ethanol complexed with pyridine depend very much on the functional employed and very little on the basis set, provided it is sufficiently large to give the correct quasi-linear hydrogen bond geometry. The CH proton shifts are insensitive to both the functional and the basis set. NMR shifts for all protons in several alcohol/pyridine complexes are calculated at the Perdew, Burke and Ernzerhof PBE0/cc-pVTZ//PBE0/6-311 + G(d,p) level in the gas phase. The results correlate with the shifts for the pyridine-complexed alcohols, determined by analysing data from the NMR titration of alcohols against pyridine. More pragmatically, computed shifts for a wider range of alcohols correlate with experimental shifts in neat pyridine. Shifts for alcohols in dimethylsulfoxide, based on the corresponding complexes in the gas phase, correlate well with the experimental values, but the overall root mean square difference is high (0.23 ppm), shifts for the OH, CHOH and other CH protons being systematically overestimated, by averages of 0.42, 0.21 and 0.06 ppm, respectively. If the computed shifts are corrected accordingly, a very good correlation is obtained with a gradient of 1.00 ± 0.01, an intercept of 0.00 ± 0.02 ppm and a root mean square difference of 0.09 ppm. This is a modest improvement on the result of applying the CHARGE programme to a slightly different set of alcohols. Some alcohol complexes with acetone and acetonitrile were investigated both in the gas phase and in a continuum of the relevant solvent.

Solvation structure and transport properties of alkali cations in dimethyl sulfoxide under exogenous static electric fields

A combination of molecular dynamics simulations and pulsed field gradient nuclear magnetic resonance spectroscopy is used to investigate the role of exogenous electric fields on the solvation structure and dynamics of alkali ions in dimethyl sulfoxide (DMSO) and as a function of temperature. Good agreement was obtained, for select alkali ions in the absence of an electric field, between calculated and experimentally determined diffusion coefficients normalized to that of pure DMSO. Our results indicate that temperatures of up to 400 K and external electric fields of up to 1 V nm(-1) have minimal effects on the solvation structure of the smaller alkali cations (Li(+) and Na(+)) due to their relatively strong ion-solvent interactions, whereas the solvation structures of the larger alkali cations (K(+), Rb(+), and Cs(+)) are significantly affected. In addition, although the DMSO exchange dynamics in the first solvation shell differ markedly for the two groups, the drift velocities and mobilities are not significantly affected by the nature of the alkali ion. Overall, although exogenous electric fields induce a drift displacement, their presence does not significantly affect the random diffusive displacement of the alkali ions in DMSO. System temperature is found to have generally a stronger influence on dynamical properties, such as the DMSO exchange dynamics and the ion mobilities, than the presence of electric fields.

Aβ monomers transiently sample oligomer and fibril-like configurations: ensemble characterization using a combined MD/NMR approach

Amyloid β (Aβ) peptides are a primary component of fibrils and oligomers implicated in the etiology of Alzheimer's disease (AD). However, the intrinsic flexibility of these peptides has frustrated efforts to investigate the secondary and tertiary structure of Aβ monomers, whose conformational landscapes directly contribute to the kinetics and thermodynamics of Aβ aggregation. In this work, de novo replica exchange molecular dynamics (REMD) simulations on the microseconds-per-replica timescale are used to characterize the structural ensembles of Aβ42, Aβ40, and M35-oxidized Aβ42, three physiologically relevant isoforms with substantially different aggregation properties. J-coupling data calculated from the REMD trajectories were compared to corresponding NMR-derived values acquired through two different pulse sequences, revealing that all simulations converge on the order of hundreds of nanoseconds-per-replica toward ensembles that yield good agreement with experiment. Though all three Aβ species adopt highly heterogeneous ensembles, these are considerably more structured compared to simulations on shorter timescales. Prominent in the C-terminus are antiparallel β-hairpins between L17-A21, A30-L36, and V39-I41, similar to oligomer and fibril intrapeptide models that expose these hydrophobic side chains to solvent and may serve as hotspots for self-association. Compared to reduced Aβ42, the absence of a second β-hairpin in Aβ40 and the sampling of alternate β topologies by M35-oxidized Aβ42 may explain the reduced aggregation rates of these forms. A persistent V24-K28 bend motif, observed in all three species, is stabilized by buried backbone to side-chain hydrogen bonds with D23 and a cross-region salt bridge between E22 and K28, highlighting the role of the familial AD-linked E22 and D23 residues in Aβ monomer folding. These characterizations help illustrate the conformational landscapes of Aβ monomers at atomic resolution and provide insight into the early stages of Aβ aggregation pathways.

Simplicity within the complexity: bilateral impact of DMSO on the functional and unfolding patterns of α-chymotrypsin

New understanding of the fundamental links between protein stability, conformational flexibility and function, can be gained through synergic studies on their catalytic and folding/unfolding properties under the influence of stabilizing/destabilizing additives. We explored an impact of dimethyl sulfoxide (DMSO), the moderate effector of multilateral action, on the kinetic (functional) and thermodynamic (thermal unfolding) patterns of a hydrolytic enzyme, α-chymotrypsin (α-CT), over a wide range of additive concentrations, 0-70% (v/v). Both the calorimetric and kinetic data exhibited rich behavior pointing to the complex interplay of global/local stability (and flexibility) patterns. The complex action of DMSO is explained through the negative and positive preferential solvation motifs that prevail for the extreme opposite, native-like and unfolded states, respectively, implying essential stabilization of compact domains by enhancement of interfacial water networks and destabilization of a flexible active site by direct binding of DMSO to the unoccupied specific positions intended for elongated polypeptide substrates.

Calcium-ion-induced stabilization of the protease from Bacillus cereus WQ9-2 in aqueous hydrophilic solvents: effect of calcium ion binding on the hydration shell and intramolecular interactions

The neutral protease WQ from Bacillus cereus is stable in various aqueous organic mixtures, with the exception of those containing acetonitrile (ACN) and dimethylformamide (DMF). The stability of the enzyme in aqueous hydrophilic solvents was dramatically enhanced with the addition of calcium ions, with the degree of improvement in the half-life relative to different solutions ranging from fourfold to more than 70-fold. Studies of the kinetic constants showed that calcium ions induced slight conformational changes in the active site of the enzyme in aqueous ACN. We investigated the molecular mechanisms underlying this stabilizing effect by employing a combination of biophysical techniques and molecular dynamics simulation. In aqueous ACN, the intrinsic fluorescence and circular dichroism analysis demonstrated that the addition of calcium ions induced a relatively compact conformation and maintained both the native-like microenvironment near the tryptophan residues and the secondary structure. Alternatively, homology modeling confirmed the location of four calcium-ion-binding sites in the enzyme, and molecular dynamics simulation revealed that three other calcium ions were bound to the surface of the enzyme. Calcium ions, known as a type of kosmotrope, can strongly bond with water molecules, thus aiding in the formation of the regional hydration shell required for the maintenance of enzyme activity. In addition, the introduction of calcium ions resulted in the formation of additional ionic interactions, providing propitious means for protein stabilization. Thus, the stronger intramolecular interactions were also expected to contribute partially to the enhanced stability of the enzyme in an aqueous organic solvent.

Effect of prolonged exposure to organic solvents on the active site environment of subtilisin Carlsberg

The potential of enzyme catalysis as a tool for organic synthesis is nowadays indisputable, as is the fact that organic solvents affect an enzyme's activity, selectivity and stability. Moreover, it was recently realized that an enzyme's initial activity is substantially decreased after prolonged exposure to organic media, an effect that further hampers their potential as catalysts for organic synthesis. Regrettably, the mechanistic reasons for these effects are still debatable. In the present study we have made an attempt to explain the reasons behind the partial loss of enzyme activity on prolonged exposure to organic solvents. Fluorescence spectroscopic studies of the serine protease subtilisin Carlsberg chemically modified with polyethylene glycol (PEG-SC) and inhibited with a Dancyl fluorophore, and dissolved in two organic solvents (acetonitrile and 1,4-dioxane) indicate that when the enzyme is initially introduced into these solvents, the active site environment is similar to that in water; however prolonged exposure to the organic medium causes this environment to resemble that of the solvent in which the enzyme is dissolved. Furthermore, kinetic studies show a reduction on both V(max) and K(M) as a result of prolonged exposure to the solvents. One interpretation of these results is that during this prolonged exposure to organic solvents the active-site fluorescent label inhibitor adopts a different binding conformation. Extrapolating this to an enzymatic reaction we argue that substrates bind in a less catalytically favorable conformation after the enzyme has been exposed to organic media for several hours.

Ligand-Exchange Processes on Solvated Lithium Cations: DMSO and Water/DMSO Mixtures

Solutions of LiClO(4) in solvent mixtures consisting of dimethylsulfoxide (DMSO) and water, or DMSO and gamma-butyrolactone, were studied by (7)Li NMR spectroscopy (for complexation by cryptands in gamma-butyrolactone as a solvent, see: E. Pasgreta, R. Puchta, M. Galle, N. J. R. van Eikema Hommes, A. Zahl, R. van Eldik, J. Incl. Phen., 2007, 58, 81-88). Chemical shifts indicate that the Li(+) ion is coordinated by four DMSO molecules. In the binary solvent mixture of water and DMSO, no selective solvation is detected, thus indicating that on increasing the water content of the solvent mixture, DMSO is gradually displaced by water in the coordination sphere of Li(+). The ligand-exchange mechanism of Li(+) ions solvated by DMSO and water/DMSO mixtures was studied using DFT calculations. Ligand exchange on [Li(DMSO)(4)](+) was found to follow a limiting associative (A) mechanism. The displacement of coordinated H(2)O by DMSO in [Li(H(2)O)(4)](+) follows an associative interchange mechanism. The suggested mechanisms are discussed in reference to available experimental and theoretical data.

Modeling hydration mechanisms of enzymes in nonpolar and polar organic solvents

A comprehensive study of the hydration mechanism of an enzyme in nonaqueous media was done using molecular dynamics simulations in five organic solvents with different polarities, namely, hexane, 3-pentanone, diisopropyl ether, ethanol, and acetonitrile. In these solvents, the serine protease cutinase from Fusarium solani pisi was increasingly hydrated with 12 different hydration levels ranging from 5% to 100% (w/w) (weight of water/weight of protein). The ability of organic solvents to 'strip off' water from the enzyme surface was clearly dependent on the nature of the organic solvent. The rmsd of the enzyme from the crystal structure was shown to be lower at specific hydration levels, depending on the organic solvent used. It was also shown that organic solvents determine the structure and dynamics of water at the enzyme surface. Nonpolar solvents enhance the formation of large clusters of water that are tightly bound to the enzyme, whereas water in polar organic solvents is fragmented in small clusters loosely bound to the enzyme surface. Ions seem to play an important role in the stabilization of exposed charged residues, mainly at low hydration levels. A common feature is found for the preferential localization of water molecules at particular regions of the enzyme surface in all organic solvents: water seems to be localized at equivalent regions of the enzyme surface independently of the organic solvent employed.