Non-heme iron and manganese species with terminal oxo ligands are proposed to be key intermediates in a variety of biological and synthetic systems; however, the stabilization of these types of complexes has proven difficult because of the tendency to form oxo-bridged complexes. Described herein are the design, isolation, and properties for a series of mononuclear Fe(III) and Mn(III) complexes with terminal oxo or hydroxo ligands. Isolation of the complexes was facilitated by the tripodal ligand tris[(N'-tert-butylureaylato)-N-ethyl]aminato ([H(3)1](3-)), which creates a protective hydrogen bond cavity around the M(III)-O(H) units (M(III) = Fe and Mn). The M(III)-O(H) complexes are prepared by the activation of dioxygen and deprotonation of water. In addition, the M(III)-O(H) complexes can be synthesized using oxygen atom transfer reagents such as N-oxides and hydroxylamines. The [Fe(III)H(3)1(O)](2-) complex also can be made using sulfoxides. These findings support the proposal of a high valent M(IV)-oxo species as an intermediate during dioxygen cleavage. Isotopic labeling studies show that oxo ligands in the [M(III)H(3)1(O)](2-) complexes come directly from the cleavage of dioxygen: for [Fe(III)H(3)1(O)](2-) the nu(Fe-(16)O) = 671 cm(-1), which shifts 26 cm(-1) in [Fe(III)H(3)1((18)O)](2-) (nu(Fe-(18)O) = 645 cm(-1)); a nu(Mn-(16)O) = 700 cm(-1) was observed for [Mn(III)H(3)1((16)O)](2-), which shifts to 672 cm(-1) in the Mn-(18)O isotopomer. X-ray diffraction studies show that the Fe-O distance is 1.813(3) A in [Fe(III)H(3)1(O)](2-), while a longer bond is found in [Fe(III)H(3)1(OH)](-) (Fe-O at 1.926(2) A); a similar trend was found for the Mn(III)-O(H) complexes, where a Mn-O distance of 1.771(5) A is observed for [Mn(III)H(3)1(O)](2-) and 1.873(2) A for [Mn(III)H(3)1(OH)](-). Strong intramolecular hydrogen bonds between the urea NH groups of [H(3)1](3-) and the oxo and oxygen of the hydroxo ligand are observed in all the complexes. These findings, along with density functional theory calculations, indicate that a single sigma-bond exists between the M(III) centers and the oxo ligands, and additional interactions to the oxo ligands arise from intramolecular H-bonds, which illustrates that noncovalent interactions may replace pi-bonds in stabilizing oxometal complexes.

The mechanism of the OH bond reorientation in liquid methanol and ethanol is examined. It is found that the extended jump model, recently developed for water, describes the OH reorientation in these liquids. The slower reorientational dynamics in these alcohols compared to water can be explained by two key factors. The alkyl groups on the alcohol molecules exclude potential partners for hydrogen bonding exchanges, an effect that grows with the size of the alkyl chain. This increases the importance of the reorientation of intact hydrogen bonds, which also slows with increasing size of the alcohol and becomes the dominant reorientation pathway.

The low-frequency collective vibrational modes in proteins as well as the protein-water interface have been suggested as dominant factors controlling the efficiency of biochemical reactions and biological energy transport. It is thus crucial to uncover the mystery of the hydration structure and dynamics as well as their coupling to collective motions of proteins in aqueous solutions. Here, we report dielectric properties of aqueous bovine serum albumin protein solutions as a model system using an extremely sensitive dielectric spectrometer with frequencies spanning from megahertz to terahertz. The dielectric relaxation spectra reveal several polarization mechanisms at the molecular level with different time constants and dielectric strengths, reflecting the complexity of protein-water interactions. Combining the effective-medium approximation and molecular dynamics simulations, we have determined collective vibrational modes at terahertz frequencies and the number of water molecules in the tightly bound and loosely bound hydration layers. High-precision measurements of the number of hydration water molecules indicate that the dynamical influence of proteins extends beyond the first solvation layer, to around 7 Å distance from the protein surface, with the largest slowdown arising from water molecules directly hydrogen-bonded to the protein. Our results reveal critical information of protein dynamics and protein-water interfaces, which determine biochemical functions and reactivity of proteins.

The enzyme Candida Antarctica lipase B (CALB) serves here as a model for understanding connections among hydration layer dynamics, solvation shell structure, and protein surface structure. The structure and dynamics of water molecules in the hydration layer were characterized for regions of the CALB surface, divided around each α-helix, β-sheet, and loop structure. Heterogeneous hydration dynamics were observed around the surface of the enzyme, in line with spectroscopic observations of other proteins. Regional differences in the structure of the biomolecular hydration layer were found to be concomitant with variations in dynamics. In particular, it was seen that regions of higher density exhibit faster water dynamics. This is analogous to the behavior of bulk water, where dynamics (diffusion coefficients) are connected to water structure (density and tetrahedrality) by excess (or pair) entropy, detailed in the Rosenfeld scaling relationship. Additionally, effects of protein surface topology and hydrophobicity on water structure and dynamics were evaluated using multiregression analysis, showing that topology has a somewhat larger effect on hydration layer structure-dynamics. Concave and hydrophobic protein surfaces favor a less dense and more tetrahedral solvation layer, akin to a more ice-like structure, with slower dynamics. Results show that pairwise entropies of local hydration layers, calculated from regional radial distribution functions, scale logarithmically with local hydration dynamics. Thus, the Rosenfeld relationship describes the heterogeneous structure-dynamics of the hydration layer around the enzyme CALB. These findings raise the question of whether this may be a general principle for understanding the structure-dynamics of biomolecular solvation.

Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.

Gigahertz-to-terahertz spectroscopy of macromolecules in aqueous environments provides an important approach for identifying their global and transient molecular structures, as well as directly assessing hydrogen-bonding. We report dielectric properties of zwitterionic dodecylphosphocholine (DPC) micelles in aqueous solutions over a wide frequency range, from 50 MHz to 1.12 THz. The dielectric relaxation spectra reveal different polarization mechanisms at the molecular level, reflecting the complexity of DPC micelle-water interactions. We have made a deconvolution of the spectra into different components and combined them with the effective-medium approximation to separate delicate processes of micelles in water. Our measurements demonstrate reorientational motion of the DPC surfactant head groups within the micelles, and two levels of hydration water shells, including tightly and loosely bound hydration water layers. From the dielectric strength of bulk water in DPC solutions, we found that the number of waters in hydration shells is approximately constant at 950 ± 45 water molecules per micelle in DPC concentrations up to 400 mM, and it decreases after that. At terahertz frequencies, employing the effective-medium approximation, we estimate that each DPC micelle is surrounded by a tightly bound layer of 310 ± 45 water molecules that behave as if they are an integral part of the micelle. Combined with molecular dynamics simulations, we determine that tightly bound waters are directly hydrogen-bonded to oxygens of DPC, while loosely bound waters reside within 4 Å of micellar atoms. The dielectric response of DPC micelles at terahertz frequencies yields, for the first time, experimental information regarding the largest scale, lowest frequency collective motions in micelles. DPC micelles are a relatively simple biologically relevant system, and this work paves the way for more insight into future studies of hydration and dynamics of biomolecular systems with gigahertz-to-terahertz spectroscopy.

Fluorine NMR spectroscopy is a powerful tool for studying biomolecular structure, dynamics, and ligand binding, yet the origins of 19F chemical shifts are not well understood. Herein, we use electronic structure calculations to describe the changes in 19F chemical shifts of 2F- and 4F-histidine/(5-methyl)-imidazole upon acid titration. While the protonation of the 2F species results in a deshielded chemical shift, protonation of the 4F results in an opposite, shielded chemical shift. The deshielding of 2F-histidine/(5-methyl)-imidazole upon protonation can be rationalized by concomitant decreases in charge density on fluorine and a reduced dipole moment. These correlations do not hold for 4F-histidine/(5-methyl)-imidazole, however. Molecular orbital calculations reveal that for the 4F species, there are no lone pair electrons on the fluorine until protonation. Analysis of a series of 4F-imidazole analogues, all with delocalized fluorine electron density, indicates that the deshielding of 19F chemical shifts through substituent effects correlates with increased C-F bond polarity. In summary, the delocalization of fluorine electrons in the neutral 4F species, with gain of a lone pair upon protonation may help explain the difficulty in developing a predictive framework for fluorine chemical shifts. Ideas debated by chemists over 40 years ago, regarding fluorine's complex electronic effects, are shown to have relevance for understanding and predicting fluorine NMR spectra.

In this paper, the role of entropy and internal energy in determining the free energy surfaces of a model dye molecule dissolved in nanoconfined solvent is investigated by Monte Carlo simulations. Three solvents (CH3I, CH3CN, and CH3OH) confined in smooth-walled, spherical nanocavities with radii of 10 and 15 Å are considered. The results provide insight into the driving forces that determine the differences in position in the nanocavity between the ground and the excited states of the solute, differences that have previously been shown to be important in time-dependent fluorescence in these systems [J. Chem. Phys. 2004, 120, 8125]. The results indicate that entropic contributions are central in controlling the shape of the free energy surfaces of the ground and excited-state solute. The effects of solvent molecule shape and charge are explored, and the possible molecular-level origins of the entropic effects are discussed.