Computation of Overhauser dynamic nuclear polarization processes reveals fundamental correlation between water dynamics, structure, and solvent restructuring entropy.

Hydration water dynamics, structure, and thermodynamics are crucially important to understand and predict water-mediated properties at molecular interfaces. Yet experimentally and directly quantifying water behavior locally near interfaces at the sub-nanometer scale is challenging, especially at interfaces submerged in biological solutions. Overhauser dynamic nuclear polarization (ODNP) experiments measure equilibrium hydration water dynamics within 8-15 angstroms of a nitroxide spin probe on instantaneous timescales (10 picoseconds to nanoseconds), making ODNP a powerful tool for probing local water dynamics in the vicinity of the spin probe. As with other spectroscopic techniques, concurrent computational analysis is necessary to gain access to detailed molecular level information about the dynamic, structural, and thermodynamic properties of water from experimental ODNP data. We chose a model system that can systematically tune the dynamics of water, a water-glycerol mixture with compositions ranging from 0 to 0.3 mole fraction glycerol. We demonstrate the ability of molecular dynamics (MD) simulations to compute ODNP spectroscopic quantities, and show that translational, rotational, and hydrogen bonding dynamics of hydration water align strongly with spectroscopic ODNP parameters. Moreover, MD simulations show tight correlations between the dynamic properties of water that ODNP captures and the structural and thermodynamic behavior of water. Hence, experimental ODNP readouts of varying water dynamics suggest changes in local structural and thermodynamic hydration water properties.

Evidence for Entropically Controlled Interfacial Hydration in Mesoporous Organosilicas.

At aqueous interfaces, the distribution and dynamics of adsorbates are modulated by the behavior of interfacial water. Hydration of a hydrophobic surface can store entropy via the ordering of interfacial water, which contributes to the Gibbs energy of solute binding. However, there is little experimental evidence for the existence of such entropic reservoirs, and virtually no precedent for their rational design in systems involving extended interfaces. In this study, two series of mesoporous silicas were modified in distinct ways: (1) progressively deeper thermal dehydroxylation, via condensation of surface silanols, and (2) increasing incorporation of nonpolar organic linkers into the silica framework. Both approaches result in decreasing average surface polarity, manifested in a blue-shift in the fluorescence of an adsorbed dye. For the inorganic silicas, hydrogen-bonding of water becomes less extensive as the number of surface silanols decreases. Overhauser dynamic nuclear polarization (ODNP) relaxometry indicates enhanced surface water diffusivity, reflecting a loss of enthalpic hydration. In contrast, organosilicas show a monotonic decrease in surface water diffusivity with decreasing polarity, reflecting enhanced hydrophobic hydration. Molecular dynamics simulations predict increased tetrahedrality of interfacial water for the organosilicas, implying increased ordering near the nm-size organic domains (relative to inorganic silicas, which necessarily lack such domains). These findings validate the prediction that hydrophobic hydration at interfaces is controlled by the microscopic length scale of the hydrophobic regions. They further suggest that the hydration thermodynamics of structurally heterogeneous silica surfaces can be tuned to promote adsorption, which in turn tunes the selectivity in catalytic reactions.

1H-HYSCORE Reveals Structural Details at the Fe(II) Active Site of Taurine:2-Oxoglutarate Dioxygenase.

Proton Hyperfine Sublevel Correlation (1H-HYSCORE) experiments have been used to probe the ligation structure of the Fe(II) active site of taurine:2-oxoglutarate dioxygenase (TauD), a non-heme Fe(II) hydroxylase. To facilitate Electron Paramagnetic Resonance (EPR) experiments, Fe(II) derivatives of the enzyme were studied using nitric oxide as a substitute for molecular oxygen. The addition of NO to the enzyme yields an S = 3/2 {FeNO}7 paramagnetic center characterized by nearly axial EPR spectra with g⊥ = 4 and g|| = 2. Using results from (i) an X-ray crystallographic study of TauD crystallized under anaerobic conditions in the presence of both cosubstrate 2-oxoglutarate and substrate taurine, (ii) a published theoretical description of the {FeNO}7 derivative of this form of the enzyme, and (iii) previous 2H-Electron Spin Echo Envelope Modulation (ESEEM) studies, we were able to assign the proton cross peaks detected in orientation-selected 1H-HYSCORE spectra. Discrete contributions from the protons of two coordinated histidine ligands were resolved. If substrate taurine is absent from the complex, orientation-selective HYSCORE spectra show cross peaks that are less resolved and when combined with information obtained from continuous wave EPR, support an alternate binding scheme for 2-oxoglutarate. HYSCORE studies of TauD in the absence of 2-oxoglutarate show additional 1H cross peaks that can be assigned to two distinct bound water molecules. In addition, 1H and 14N cross peaks that arise from the coordinated histidine side chains show a change in NO coordination for this species. For all of the TauD species, 1H hyperfine couplings and their orientations are sensitive to the detailed electronic structure of the {FeNO}7 center.

Pulsed EPR characterization of HIV-1 protease conformational sampling and inhibitor-induced population shifts.

The conformational landscape of HIV-1 protease (PR) can be experimentally characterized by pulsed-EPR double electron-electron resonance (DEER). For this characterization, nitroxide spin labels are attached to an engineered cysteine residue in the flap region of HIV-1 PR. DEER distance measurements from spin-labels contained within each flap of the homodimer provide a detailed description of the conformational sampling of apo-enzyme as well as induced conformational shifts as a function of inhibitor binding. The distance distribution profiles are further interpreted in terms of a conformational ensemble scheme that consists of four unique states termed "curled/tucked", "closed", "semi-open" and "wide-open" conformations. Reported here are the DEER results for a drug-resistant variant clinical isolate sequence, V6, in the presence of FDA approved protease inhibitors (PIs) as well as a non-hydrolyzable substrate mimic, CaP2. Results are interpreted in the context of the current understanding of the relationship between conformational sampling, drug resistance, and kinetic efficiency of HIV-1PR as derived from previous DEER and kinetic data for a series of HIV-1PR constructs that contain drug-pressure selected mutations or natural polymorphisms. Specifically, these collective results support the notion that inhibitor-induced closure of the flaps correlates with inhibitor efficiency and drug resistance. This body of work also suggests DEER as a tool for studying conformational sampling in flexible enzymes as it relates to function.

Continuous wave W- and D-band EPR spectroscopy offer "sweet-spots" for characterizing conformational changes and dynamics in intrinsically disordered proteins.

Site-directed spin labeling (SDSL) electron paramagnetic resonance (EPR) spectroscopy is a powerful tool for characterizing conformational sampling and dynamics in biological macromolecules. Here we demonstrate that nitroxide spectra collected at frequencies higher than X-band (∼9.5 GHz) have sensitivity to the timescale of motion sampled by highly dynamic intrinsically disordered proteins (IDPs). The 68 amino acid protein IA3, was spin-labeled at two distinct sites and a comparison of X-band, Q-band (35 GHz) and W-band (95 GHz) spectra are shown for this protein as it undergoes the helical transition chemically induced by tri-fluoroethanol. Experimental spectra at W-band showed pronounced line shape dispersion corresponding to a change in correlation time from ∼0.3 ns (unstructured) to ∼0.6 ns (α-helical) as indicated by comparison with simulations. Experimental and simulated spectra at X- and Q-bands showed minimal dispersion over this range, illustrating the utility of SDSL EPR at higher frequencies for characterizing structural transitions and dynamics in IDPs.

Nanoscale water-polymer interactions tune macroscopic diffusivity of water in aqueous poly(ethylene oxide) solutions.

The separation and anti-fouling performance of water purification membranes is governed by both macroscopic and molecular-scale water properties near polymer surfaces. However, even for poly(ethylene oxide) (PEO) - ubiquitously used in membrane materials - there is little understanding of whether or how the molecular structure of water near PEO surfaces affects macroscopic water diffusion. Here, we probe both time-averaged bulk and local water dynamics in dilute and concentrated PEO solutions using a unique combination of experimental and simulation tools. Pulsed-Field Gradient NMR and Overhauser Dynamic Nuclear Polarization (ODNP) capture water dynamics across micrometer length scales in sub-seconds to sub-nanometers in tens of picoseconds, respectively. We find that classical models, such as the Stokes-Einstein and Mackie-Meares relations, cannot capture water diffusion across a wide range of PEO concentrations, but that free volume theory can. Our study shows that PEO concentration affects macroscopic water diffusion by enhancing the water structure and altering free volume. ODNP experiments reveal that water diffusivity near PEO is slower than in the bulk in dilute solutions, previously not recognized by macroscopic transport measurements, but the two populations converge above the polymer overlap concentration. Molecular dynamics simulations reveal that the reduction in water diffusivity occurs with enhanced tetrahedral structuring near PEO. Broadly, we find that PEO does not simply behave like a physical obstruction but directly modifies water's structural and dynamic properties. Thus, even in simple PEO solutions, molecular scale structuring and the impact of polymer interfaces is essential to capturing water diffusion, an observation with important implications for water transport through structurally complex membrane materials.

Spin labeling and Double Electron-Electron Resonance (DEER) to Deconstruct Conformational Ensembles of HIV Protease.

An understanding of macromolecular conformational equilibrium in biological systems is oftentimes essential to understand function, dysfunction, and disease. For the past few years, our lab has been utilizing site-directed spin labeling (SDSL), coupled with electron paramagnetic resonance (EPR) spectroscopy, to characterize the conformational ensemble and ligand-induced conformational shifts of HIV-1 protease (HIV-1PR). The biomedical importance of characterizing the fractional occupancy of states within the conformational ensemble critically impacts our hypothesis of a conformational selection mechanism of drug-resistance evolution in HIV-1PR. The purpose of the following chapter is to give a timeline perspective of our SDSL EPR approach to characterizing conformational sampling of HIV-1PR. We provide detailed instructions for the procedure utilized in analyzing distance profiles for HIV-1PR obtained from pulsed electron-electron double resonance (PELDOR). Specifically, we employ a version of PELDOR known as double electron-electron resonance (DEER). Data are processed with the software package "DeerAnalysis" (http://www.epr.ethz.ch/software), which implements Tikhonov regularization (TKR), to generate a distance profile from electron spin-echo amplitude modulations. We assign meaning to resultant distance profiles based upon a conformational sampling model, which is described herein. The TKR distance profiles are reconstructed with a linear combination of Gaussian functions, which is then statistically analyzed. In general, DEER has proven powerful for observing structural ensembles in proteins and, more recently, nucleic acids. Our goal is to present our advances in order to aid readers in similar applications.

Measuring the orientation of taurine in the active site of the non-heme Fe(II)/α-ketoglutarate-dependent taurine hydroxylase (TauD) using electron spin echo envelope modulation (ESEEM) spectroscopy.

The position and orientation of taurine near the non-heme Fe(II) center of the α-ketoglutarate (α-KG)-dependent taurine hydroxylase (TauD) was measured using Electron Spin Echo Envelope Modulation (ESEEM) spectroscopy. TauD solutions containing Fe(II), α-KG, and natural abundance taurine or specifically deuterated taurine were prepared anaerobically and treated with nitric oxide (NO) to make an S = 3/2 {FeNO}(7) complex that is suitable for robust analysis with EPR spectroscopy. Using ratios of ESEEM spectra collected for TauD samples having natural abundance taurine or deuterated taurine, (1)H and (14)N modulations were filtered out of the spectra and interactions with specific deuterons on taurine could be studied separately. The Hamiltonian parameters used to calculate the amplitudes and line shapes of frequency spectra containing isolated deuterium ESEEM were obtained with global optimization algorithms. Additional statistical analysis was performed to validate the interpretation of the optimized parameters. The strongest (2)H hyperfine coupling was to a deuteron on the C1 position of taurine and was characterized by an effective dipolar distance of 3.90 ± 0.25 Å from the {FeNO}(7) paramagnetic center. The principal axes of this C1-(2)H hyperfine coupling and nuclear quadrupole interaction tensors were found to make angles of 26 ± 5 and 52 ± 17°, respectively, with the principal axis of the {FeNO}(7) zero-field splitting tensor. These results are discussed within the context of the orientation of substrate taurine prior to the initiation of hydrogen abstraction.