Mapping oxygen accessibility to ribonuclease a using high-resolution NMR relaxation spectroscopy

Paramagnetic Relaxation Agents for Enhancing Temporal Resolution and Sensitivity in Multinuclear FlowNMR Spectroscopy

Sensitivity in FlowNMR spectroscopy for reaction monitoring often suffers from low levels of pre-magnetisation due to limited residence times of the sample in the magnetic field. While this in-flow effect is tolerable for high sensitivity nuclei such as 1 H and 19 F, it significantly reduces the signal-to-noise ratio in 31 P and 13 C spectra, making FlowNMR impractical for low sensititvity nuclei at low concentrations. Paramagnetic relaxation agents (PRAs), which enhance polarisation and spin-lattice relaxation, could eliminate the adverse in-flow effect and improve the signal-to-noise ratio. Herein, [Co(acac)3 ], [Mn(acac)3 ], [Fe(acac)3 ], [Cr(acac)3 ], [Ni(acac)2 ]3, [Gd(tmhd)3 ] and [Cr(tmhd)3 ] are investigated for their effectiveness in improving signal intensity per unit time in FlowNMR applications under the additional constraint of chemical inertness towards catalytically active transition metal complexes. High-spin Cr(III) acetylacetonates emerged as the most effective compounds, successfully reducing 31 P T1 values four- to five-fold at PRA concentrations as low as 10 mM without causing adverse line broadening. Whereas [Cr(acac)3 ] showed signs of chemical reactivity with a mixture of triphenylphosphine, triphenylphosphine oxide and triphenylphosphate over the course of several hours at 80° C, the bulkier [Cr(tmhd)3 ] was stable and equally effective as a PRA under these conditions. Compatibility with a range of representative transition metal complexes often used in homogeneous catalysis has been investigated, and application of [Cr(tmhd)3 ] in significantly improving 1 H and 31 P{1 H} FlowNMR data quality in a Rh-catalysed hydroformylation reaction has been demonstrated. With the PRA added, 13 C relaxation times were reduced more than six-fold, allowing quantitative reaction monitoring of substrate consumption and product formation by 13 C{1 H} FlowNMR spectroscopy at natural abundance.

Ordered Motions in the Nitric-Oxide Dioxygenase Mechanism of Flavohemoglobin and Assorted Globins with Tightly Coupled Reductases

Nitric-oxide dioxygenases (NODs) activate and combine O₂ with NO to form nitrate. A variety of oxygen-binding hemoglobins with associated partner reductases or electron donors function as enzymatic NODs. Kinetic and structural investigations of the archetypal two-domain microbial flavohemoglobin-NOD have illuminated an allosteric mechanism that employs selective tunnels for O₂ and NO, gates for NO and nitrate, transient O₂ association with ferric heme, and an O₂ and NO-triggered, ferric heme spin crossover-driven, motion-controlled, and dipole-regulated electron-transfer switch. The proposed mechanism facilitates radical-radical coupling of ferric-superoxide with NO to form nitrate while preventing suicidal ferrous-NO formation. Diverse globins display the structural and functional motifs necessary for a similar allosteric NOD mechanism. In silico docking simulations reveal monomeric erythrocyte hemoglobin alpha-chain and beta-chain intrinsically matched and tightly coupled with NADH-cytochrome b₅ oxidoreductase and NADPH-cytochrome P450 oxidoreductase, respectively, forming membrane-bound flavohemoglobin-like mammalian NODs. The neuroprotective neuroglobin manifests a potential NOD role in a close-fitting ternary complex with membrane-bound NADH-cytochrome b₅ oxidoreductase and cytochrome b₅. Cytoglobin interfaces weakly with cytochrome b₅ for O₂ and NO-regulated electron-transfer and coupled NOD activity. The mechanistic model also provides insight into the evolution of O₂ binding cooperativity in hemoglobin and a basis for the discovery of allosteric NOD inhibitors.

Theoretical Studies on the Effect of Isomerized Aspartic Acid Residues on the Three-Dimensional Structures of Bovine Pancreatic Ribonucleases A

Isomerized aspartic acid (Asp) residues have previously been identified in various aging tissues, and are suspected to contribute to age-related diseases. Asp-residue isomerization occurs nonenzymatically under physiological conditions, resulting in the formation of three types of isomerized Asp (i.e., L-isoAsp, D-Asp, and D-isoAsp) residues. Asp-residue isomerization often accelerates protein aggregation and insolubilization, making structural biology analyses difficult. Recently, Sakaue et al. reported the synthesis of a ribonuclease A (RNase A) in which Asp121 was artificially replaced with different isomerized Asp residues, and experimentally demonstrated that the enzymatic activities of these artificial mutants were completely lost. However, their structural features have not yet been elucidated. In the present study, the three-dimensional (3D) structures of these artificial-mutant RNases A were predicted using molecular dynamics (MD) simulations. The 3D structures of wild-type and artificial-mutant RNases A were converged by 3000-ns MD simulations. Our computational data show that the structures of the active site and the formation frequencies of the appropriate catalytic dyad structures in the artificial-mutant RNases A were quite different from wild-type RNase A. These computational findings may provide an explanation for the experimental data which show that artificial-mutant RNases A lack enzymatic activity. Herein, MD simulations have been used to evaluate the influences of isomerized Asp residues on the 3D structures of proteins.

Allostery in the nitric oxide dioxygenase mechanism of flavohemoglobin

The substrates O₂ and NO cooperatively activate the NO dioxygenase function of Escherichia coli flavohemoglobin. Steady-state and transient kinetic measurements support a structure-based mechanistic model in which O₂ and NO movements and conserved amino acids at the E11, G8, E2, E7, B10, and F7 positions within the globin domain control activation. In the cooperative and allosteric mechanism, O₂ migrates to the catalytic heme site via a long hydrophobic tunnel and displaces LeuE11 away from the ferric iron, which forces open a short tunnel to the catalytic site gated by the ValG8/IleE15 pair and LeuE11. NO permeates this tunnel and leverages upon the gating side chains triggering the CD loop to furl, which moves the E and F-helices and switches an electron transfer gate formed by LysF7, GlnE7, and water. This allows FADH₂ to reduce the ferric iron, which forms the stable ferric–superoxide–TyrB10/GlnE7 complex. This complex reacts with internalized NO with a bimolecular rate constant of 10¹⁰ M⁻¹ s⁻¹ forming nitrate, which migrates to the CD loop and unfurls the spring-like structure. To restart the cycle, LeuE11 toggles back to the ferric iron. Actuating electron transfer with O₂ and NO movements averts irreversible NO poisoning and reductive inactivation of the enzyme. Together, structure snapshots and kinetic constants provide glimpses of intermediate conformational states, time scales for motion, and associated energies.

Nuclear magnetic resonance-based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine-free wild-type T4 lysozyme was investigated using O₂ gas-pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂ -induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C-terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid-search. Two O₂ -accessible sites were identified in internal cavities: One site was consistent with the xenon-binding site in the protein in crystal, and the other site was established to be a novel ligand-binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground-state ensemble transiently develop channels for O₂ association with the internal protein cavities.

Analysis of O2-binding Sites in Proteins Using Gas-Pressure NMR Spectroscopy: Outer Surface Protein A

Internal cavities in proteins produce conformational fluctuations and enable the binding of small ligands. Here, we report a NMR analysis of O₂-binding sites by O₂-induced paramagnetic relaxation enhancements (PREs) on amide groups of proteins in solution. Outer surface protein A contains a nonglobular single-layer β-sheet that connects the N- and C-terminal globular domains. Several cavities have been observed in both domains of the crystallized protein structure. The receptor-binding sites are occluded and line the largest cavity of the C-terminal domain. We observed significant O₂-induced PREs for amide protons located around the largest cavity and at the central β-sheet. We suggested three potential O₂-accessible sites in the protein based on the 1/r⁶ distance dependence of the PRE. Two sites were in or close to the largest cavity and the third site was in the surface crevice of the central β-sheet. These results provide, to our knowledge, the first evidence of ligand binding to the surface crevice and cavity of the protein in solution. Because O₂ generally binds more specifically to hydrophobic rather than hydrophilic cavities within a protein, the results also indicated that the receptor-binding sites lining the largest cavity were in the hydrophobic environment in the ground-state conformation. Molecular dynamics simulations permitted the visualization of the rotational and translational motions of O₂ within the largest cavity, egress of O₂ from the cavity, and ingress of O₂ in the surface crevice of the β-sheet. These molecular dynamics simulation results qualitatively explained the O₂-induced changes in NMR observations. Exploring cavities that are sufficiently dynamic to enable access by small molecules can be a useful strategy for the design of stable proteins and their ligands.

Detecting O2 binding sites in protein cavities

Internal cavities are important elements in protein structure, dynamics, stability and function. Here we use NMR spectroscopy to investigate the binding of molecular oxygen (O₂) to cavities in a well-studied model for ligand binding, the L99A mutant of T4 lysozyme. On increasing the O₂ concentration to 8.9 mM, changes in ¹H, ¹⁵N, and ¹³C chemical shifts and signal broadening were observed specifically for backbone amide and side chain methyl groups located around the two hydrophobic cavities of the protein. O₂-induced longitudinal relaxation enhancements for amide and methyl protons could be adequately accounted for by paramagnetic dipolar relaxation. These data provide the first experimental demonstration that O₂ binds specifically to the hydrophobic, and not the hydrophilic cavities, in a protein. Molecular dynamics simulations visualized the rotational and translational motions of O₂ in the cavities, as well as the binding and egress of O₂, suggesting that the channel consisting of helices D, E, G, H, and J could be the potential gateway for ligand binding to the protein. Due to strong paramagnetic relaxation effects, O₂ gas-pressure NMR measurements can detect hydrophobic cavities when populated to as little as 1%, and thereby provide a general and highly sensitive method for detecting oxygen binding in proteins.

Cosolute paramagnetic relaxation enhancements detect transient conformations of human uracil DNA glycosylase (hUNG)

The human DNA repair enzyme uracil DNA glycosylase (hUNG) locates and excises rare uracil bases that arise in DNA from cytosine deamination or through dUTP incorporation by DNA polymerases. Previous NMR studies of hUNG have revealed millisecond time scale dynamic transitions in the enzyme-nonspecific DNA complex, but not the free enzyme, that were ascribed to a reversible clamping motion of the enzyme as it scans along short regions of duplex DNA in its search for uracil. Here we further probe the properties of the nonspecific DNA binding surface of {(2)H(12)C}{(15)N}-labeled hUNG using a neutral chelate of a paramagnetic Gd(3+) cosolute (Gd(HP-DO3A)). Overall, the measured paramagnetic relaxation enhancements (PREs) on R(2) of the backbone amide protons for free hUNG and its DNA complex were in good agreement with those calculated based on their relative exposure observed in the crystal structures of both enzyme forms. However, the calculated PREs systematically underestimated the experimental PREs by large amounts in discrete regions implicated in DNA recognition and catalysis: active site loops involved in DNA recognition (268-274, 246-250), the uracil binding pocket (143-148, 169-170), a transient extrahelical base binding site (214-216), and a remote hinge region (129-132) implicated in dynamic clamping. These reactive hot spots were not correlated with structural, hydrophobic, or solvent exchange properties that might be common to these regions, leaving the possibility that the effects arise from dynamic sampling of exposed conformations that are distinct from the static structures. Consistent with this suggestion, the above regions have been previously shown to be flexible based on relaxation dispersion measurements and course-grained normal-mode analysis. A model is suggested where the intrinsic dynamic properties of these regions allows sampling of transient conformations where the backbone amide groups have greater average exposure to the cosolute as compared to the static structures. We conclude that PREs derived from the paramagnetic cosolute reveal dynamic hot spots in hUNG and that these regions are highly correlated with substrate binding and recognition.

Nitroxides: Versatile Reporters and Reactants

Stable nitroxyl radicals are important tools in chemistry, biophysics, biology, and materials science. Their stability and the sensitivity of their EPR spectra to the local environment make them valuable molecular probes. This review seeks to give an overview of the developments in the field of nitroxide spin probes and their various applications, with the main focus on the pH-sensitive imidazoline nitroxide family.

Dioxygen transmembrane distributions and partitioning thermodynamics in lipid bilayers and micelles

Cellular respiration, mediated by the passive diffusion of oxygen across lipid membranes, is key to many basic cellular processes. In this work, we report the detailed distribution of oxygen across lipid bilayers and examine the thermodynamics of oxygen partitioning via NMR studies of lipids in a small unilamellar vesicle (SUV) morphology. Dissolved oxygen gives rise to paramagnetic chemical shift perturbations and relaxation rate enhancements, both of which report on local oxygen concentration. From SUVs containing the phospholipid sn-2-perdeuterio-1-myristelaidoyl, 2-myristoyl-sn-glycero-3-phosphocholine (MLMPC), an analogue of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), we deduced the complete trans-bilayer oxygen distribution by measuring (13)C paramagnetic chemical shifts perturbations for 18 different sites on MLMPC arising from oxygen at a partial pressure of 30 bar. The overall oxygen solubility at 45 °C spans a factor of 7 between the bulk water (23.7 mM) and the bilayer center (170 mM) and is lowest in the vicinity of the phosphocholine headgroup, suggesting that oxygen diffusion across the glycerol backbone should be the rate-limiting step in diffusion-mediated passive transport of oxygen across the lipid bilayer. Lowering of the temperature from 45 to 25 °C gave rise to a slight decrease of the oxygen solubility within the hydrocarbon interior of the membrane. An analysis of the temperature dependence of the oxygen solubility profile, as measured by (1)H paramagnetic relaxation rate enhancements, reveals that oxygen partitioning into the bilayer is entropically favored (ΔS° = 54 ± 3 J K(-1) mol(-1)) and must overcome an enthalpic barrier (ΔH° = 12.0 ± 0.9 kJ mol(-1)).

Water-proton-spin-lattice-relaxation dispersion of paramagnetic protein solutions

The paramagnetic contributions to water-proton-spin-lattice relaxation rate constants in protein systems spin-labeled with nitroxide radicals were re-examined. As noted by others, the strength of the dipolar coupling between water protons and the protein-bound nitroxide radical often appears to be larger than physically reasonable when the relaxation is assumed to be controlled by 3-dimensional diffusive processes in the vicinity of the spin label. We examine the effects of the surface in biasing the diffusive exploration of the radical region and derive a relaxation model that incorporates 2-dimensional dynamics at the interfacial layer. However, we find that the local 2-dimensional dynamics changes the shape of the relaxation dispersion profile but does not necessarily reproduce the low-field relaxation efficiency found by experiment. We examine the contributions of long-range dipolar couplings between the paramagnetic center and protein-bound-water molecules and find that the contributions from these several long range couplings may be competitive with translational contributions because the correlation time for global rotation of the protein is approximately 1000 times longer than that for the diffusive motion of water at the interfacial region. As a result the electron-proton dipolar coupling to rare protein-bound-water-molecule protons may be significant for protein systems that accommodate long-lived-water molecules. Although the estimate of local diffusion coefficients is not seriously compromised because it derives from the Larmor frequency dependence, these several contributions confound efforts to fit relaxation data quantitatively with unique models.

A solution NMR approach to the measurement of amphiphile immersion depth and orientation in membrane model systems

Oxygen and Ni(II) are ideal paramagnetic species for NMR studies of immersion depth since they establish prominent concentration gradients across the membrane-water interface of either bilayers or micelles. Corresponding gradients of paramagnetic shifts and relaxation rates are observed by NMR for membrane embedded amphiphiles. Specifically, upon dissolution of oxygen at a partial pressure of 20 bar or more, (13)C NMR spectra of membrane embedded amphiphiles reveal chemical shift perturbations which depend sensitively on average immersion depth in the membrane. Similarly, depth-dependent enhancements of spin-lattice relaxation rates can be detected by (1)H NMR. Generally, such paramagnetic effects depend both on steric or accessibility factors and on the local concentration of the paramagnet. The steric terms can be factored out by combining paramagnetic rates arising from O(2) and Ni, in the form of a ratio, R(1P)(O(2))/R(1P)(Ni). The natural logarithm of this ratio is shown to depend linearly on immersion depth in a micelle. The analysis is verified using molecular dynamics simulations of dodecylphosphocholine in a detergent micelle, while thorough consideration of the paramagnetic rate data also allows for the determination of the orientation of imipramine in the micelle. Thus, a complete picture of topology arises from this approach which is readily applicable to studies of drugs and amphiphiles in fast-tumbling bicelles, small unilamellar vesicles, and micelles.

Measuring the dynamic surface accessibility of RNA with the small paramagnetic molecule TEMPOL

The surface accessibility of macromolecules plays a key role in modulating molecular recognition events. RNA is a complex and dynamic molecule involved in many aspects of gene expression. However, there are few experimental methods available to measure the accessible surface of RNA. Here, we investigate the accessible surface of RNA using NMR and the small paramagnetic molecule TEMPOL. We investigated two RNAs with known structures, one that is extremely stable and one that is dynamic. For helical regions, the TEMPOL probing data correlate well with the predicted RNA surface, and the method is able to distinguish subtle variations in atom depths, such as the relative accessibility of pyrimidine versus purine aromatic carbon atoms. Dynamic motions are also detected by TEMPOL probing, and the method accurately reports a previously characterized pH-dependent conformational transition involving formation of a protonated C-A pair and base flipping. Some loop regions are observed to exhibit anomalously high accessibility, reflective of motions that are not evident within the ensemble of NMR structures. We conclude that TEMPOL probing can provide valuable insights into the surface accessibility and dynamics of RNA, and can also be used as an independent means of validating RNA structure and dynamics in solution.

Oxygen accessibility to ribonuclease a: quantitative interpretation of nuclear spin relaxation induced by a freely diffusing paramagnet

The nuclear spin relaxation induced by a freely diffusing paramagnetic center provides a direct measure of intermolecular accessibility. A number of factors are involved in a quantitative interpretation of relaxation data including excluded volume effects, solvation differences, and the details of the electron spin relaxation in the paramagnetic center. In the case where the electron relaxation time is short compared with correlation times describing the electron-nuclear coupling, the nuclear spin relaxation rates may be related to the effective local concentration of the paramagnetic center at different locations about the solute of interest. The local concentrations may in turn be related to differences in the local free energies of interaction between the diffusing paramagnet and the cosolute. We demonstrate this approach for the case of ribonuclease A and deduce surface free energy differences for a large number of protein proton sites. We find that the oxygen accessibility is poorly represented by hard-sphere models such as computed solvent or steric accessibility. There is a distribution of local intermolecular interactions with a width of the order of RT that dominates the report of the intermolecular exploration of the protein by this simple solute.

MD and NMR studies of alpha-bungarotoxin surface accessibility

Protein surface accessibility represents a dimension of structural biology which has not been discussed in details so far, in spite of its fundamental role in controlling the molecular recognition process. In the present report the surface accessibility of alpha-bungarotoxin, a small and well characterized protein, has been investigated by analyzing its interaction with solvent and paramagnetic molecules in an integrated way. The presence of strong hydration sites, identified by a combined analysis of MD simulation and NMR results, seems to prevent the access of Gd(III)DTPA-BMA to the protein surface. On the contrary, the limited hydration of the alpha-bungarotoxin active site favors frequent encounters between the paramagnetic probe and the protein in the latter region. All the data obtained here for alpha-bungarotoxin suggest that shape and stability of the solvation shell control its surface accessibility and, hence, intermolecular interactions in a way which could be common to many other proteins.

The Reaction of Hydrogen Atoms with Methionine Residues: A Model of Reductive Radical Stress Causing Tandem Protein-Lipid Damage

The occurrence of tandem damage, due to reductive radical stress involving proteins and lipids, is shown by using a biomimetic model. It is made of unsaturated lipid vesicle suspensions in phosphate buffer in the presence of methionine, either as a single amino acid or as part of a protein such as RNase A, which contains four methionine residues. The radical process starts with the formation of H(.) atoms by reaction of solvated electrons with dihydrogen phosphate anions, which selectively attack the thioether function of methionine. The modification of methionine to alpha-aminobutyric acid is accompanied by the formation of thiyl radicals, which in turn cause the isomerization of the cis fatty acid residues to the trans isomers. The relationship between methionine modification and lipid damage and some details of the reductive radical stress obtained by proteomic analysis of irradiated RNase A are presented.

Spin relaxation measurements of electrostatic bias in intermolecular exploration

We utilize the paramagnetic contribution to proton spin-lattice relaxation rate constants induced by freely diffusing charged paramagnetic centers to investigate the effect of charge on the intermolecular exploration of a protein by the small molecule. The proton NMR spectrum provided 255 resolved resonances that report how the explorer molecule local concentration varies with position on the surface. The measurements integrate over local dielectric constant variations, and, in principle, provide an experimental characterization of the surface free energy sampling biases introduced by the charge distribution on the protein. The experimental results for ribonuclease A obtained using positive, neutral, and negatively charged small nitroxide radicals are qualitatively similar to those expected from electrostatic calculations. However, while systematic electrostatic trends are apparent, the three different combinations of the data sets do not yield internally consistent values for the electrostatic contribution to the intermolecular free energy. We attribute this failure to the weakness of the electrostatic sampling bias for charged nitroxides in water and local variations in effective translational diffusion constant at the water-protein interface, which enters the nuclear spin relaxation equations for the nitroxide-proton dipolar coupling.

An NMR Method for the Determination of Protein-Binding Interfaces Using Dioxygen-Induced Spin−Lattice Relaxation Enhancement

Using oxygen as a paramagnetic probe, researchers can routinely study topologies and protein-binding interfaces by NMR. The paramagnetic contribution to the amide (1)H spin-lattice relaxation rates (R(1)(P)) have been studied for uniformly (2)H,(15)N-labeled FB protein, a 60-residue three-helix bundle, constituting the B domain of protein A. Through TROSY versions of inversion-recovery experiments, R(1)(P) could be determined. R(1)(P) was then measured in the presence of a stoichiometric equivalent of an unlabeled Fc fragment of immunoglobulin (Ig) G, and the ratio of R(1)(P) of the FB-Fc complex to that of free FB [i.e., R(1)(P)(complex)/R(1)(P)(free)] was determined for each observable residue. Regions of helix I and helix II, which were previously known to interact with Fc, were readily identified as belonging to the binding interface by their characteristically reduced values of R(1)(P)(complex)/R(1)(P)(free). The method of comparing oxygen-induced spin-lattice relaxation rates of free protein and protein-protein complexes, to detect binding interfaces, offers greater sensitivity than chemical shift perturbation, while it is not necessary to heavily deuterate the labeled protein, as is the case in cross saturation experiments.