A Combined NMR and Molecular Dynamics Study of the Transmembrane Solubility and Diffusion Rate Profile of Dioxygen in Lipid Bilayers

Generating Concentration Gradients across Membranes for Molecular Dynamics Simulations of Periodic Systems

The plasma membrane forms the boundary between a living entity and its environment and acts as a barrier to permeation and flow of substances. Several computational means of calculating permeability have been implemented for molecular dynamics (MD) simulations-based approaches. Except for double bilayer systems, most permeability studies have been performed under equilibrium conditions, in large part due to the challenges associated with creating concentration gradients in simulations utilizing periodic boundary conditions. To enhance the scientific understanding of permeation and complement the existing computational means of characterizing membrane permeability, we developed a non-equilibrium method that enables the generation and maintenance of steady-state gradients in MD simulations. We utilize PBCs advantageously by imposing a directional bias to the motion of permeants so that their crossing of the boundary replenishes the gradient, like a previous study on ions. Under these conditions, a net flow of permeants across membranes may be observed to determine bulk permeability by a direct application of J=PΔc. In the present study, we explore the results of its application to an exemplary O2 and POPC bilayer system, demonstrating accurate and precise permeability measurements. In addition, we illustrate the impact of permeant concentration and the choice of thermostat on the permeability. Moreover, we demonstrate that energetics of permeation can be closely examined by the dissipation of the gradient across the membrane to gain nuanced insights into the thermodynamics of permeability.

Factors determining barrier properties to oxygen transport across model and cell plasma membranes based on EPR spin-label oximetry

This review is motivated by the exciting new area of radiation therapy using a phenomenon termed FLASH in which oxygen is thought to have a central role. Well-established principles of radiation biology and physics suggest that if oxygen has a strong role, it should be the level at the DNA. The key aspect discussed is the rate of oxygen diffusion. If oxygen freely diffuses into cells and rapidly equilibrates, then measurements in the extracellular compartment would enable FLASH to be investigated using existing methodologies that can readily measure oxygen in the extracellular compartment. EPR spin-label oximetry allows evaluation of the oxygen permeability coefficient across lipid bilayer membranes. It is established that simple fluid phase lipid bilayers are not barriers to oxygen transport. However, further investigations indicate that many physical and chemical (compositional) factor can significantly decrease this permeation. In biological cell plasma membranes, the lipid bilayer forms the matrix in which integral membrane proteins are immersed, changing organization and properties of the lipid matrix. To evaluate oxygen permeability coefficients across these complex membranes, oxygen permeation across all membrane domains and components must be considered. In this review, we consider many of the factors that affect (decrease) oxygen permeation across cell plasma membranes. Finally, we address the question, can the plasma membrane of the cell form a barrier to the free diffusion of oxygen into the cell interior? If there is a barrier then this must be considered in the investigations of the role of oxygen in FLASH.

Atomistic simulations modify interpretation of spin-label oximetry data. Part 1: intensified water-lipid interfacial resistances

The role of membrane cholesterol in cellular function and dysfunction has been the subject of much inquiry. A few studies have suggested that cholesterol may slow oxygen diffusive transport, altering membrane physical properties and reducing oxygen permeability. The primary experimental technique used in recent years to study membrane oxygen transport is saturation-recovery electron paramagnetic resonance (EPR) oximetry, using spin-label probes targeted to specific regions of a lipid bilayer. The technique has been used, in particular, to assess the influence of cholesterol on oxygen transport and membrane permeability. The reliability of such EPR recordings at the water-lipid interface near the phospholipid headgroups has been challenged by all-atom molecular dynamics (MD) simulation data that show substantive agreement with spin-label probe measurements throughout much of the bilayer. This work uses further MD simulations, with an updated oxygen model, to determine the location of the maximum resistance to permeation and the rate-limiting barrier to oxygen permeation in 1-palmitoyl,2-oleoylphosphatidylcholine (POPC) and POPC/cholesterol bilayers at 25 and 35°C. The current simulations show a spike of resistance to permeation in the headgroup region that was not detected by EPR but was predicted in early theoretical work by Diamond and Katz. Published experimental nuclear magnetic resonance (NMR) oxygen measurements provide key validation of the MD models and indicate that the positions and relative magnitudes of the phosphatidylcholine resistance peaks are accurate. Consideration of the headgroup-region resistances predicts bilayer permeability coefficients lower than estimated in EPR studies, giving permeabilities lower than the permeability of unstirred water layers of the same thickness. Here, the permeability of POPC at 35°C is estimated to be 13 cm/s, compared with 10 cm/s for POPC/cholesterol and 118 cm/s for simulation water layers of similar thickness. The value for POPC is 12 times lower than estimated from EPR measurements, while the value for POPC/cholesterol is 5 times lower. These findings underscore the value of atomic resolution models for guiding the interpretation of experimental probe-based measurements.

Updated Evaluation of Cholesterol's Influence on Membrane Oxygen Permeability

There is a surprising gap in knowledge regarding the mechanism of oxygen (O₂) diffusional delivery at the level of tissues and cells. Yet, the effectiveness of tumor radiotherapy, the success of tissue engineering, and healthy metabolism all require ample intracellular oxygen. Tissue-level diffusion takes place in a complex and crowded macromolecular environment. Cholesterol-rich cellular membranes have been thought to reduce oxygen flux. Here, we use atomistic molecular dynamics simulations to update prior estimates of bilayer permeability and related parameters for 1-palmitoyl,2-oleoylphosphatidylcholine (POPC) and POPC/cholesterol bilayers, using a modified O₂ model with improved membrane-water partitioning behavior. This work estimates an oxygen permeability coefficient of 15 ± 1 cm/s for POPC and 11.5 ± 0.4 cm/s for POPC/cholesterol (1:1 molecular ratio) at 37 °C. The permeability of POPC is found to be ~1/3 that of a water layer of similar thickness, and the permeability of POPC/cholesterol is estimated to be 20-30% below that of POPC. Void pathway visualization and free energy data support channeling of oxygen toward the center of cholesterol-incorporating membranes, while partition coefficient data suggest reduced membrane solubility of oxygen due to cholesterol. Further study is needed to understand whether diffusion pathway changes due to cholesterol and other molecular compositional factors influence oxygen availability within tissue.

How does oxygen diffuse from capillaries to tissue mitochondria? Barriers and pathways

Timely delivery of oxygen (O₂ ) to tissue mitochondria is so essential that elaborate circulatory systems have evolved to minimize diffusion distances within tissue. Yet, knowledge is surprisingly limited regarding the diffusion pathway between blood capillaries and tissue mitochondria. An established and growing body of work examines the influence cellular and extracellular structures may have on subcellular oxygen availability. This brief review discusses the physiological and pathophysiological significance of oxygen availability, highlights recent computer modelling studies of transport at the cell-membrane level, and considers alternative diffusion pathways within tissue. Experimental and computer modelling studies suggest that oxygen diffusion may be accelerated by cellular lipids, relative to cytosolic and interstitial fluids. Such acceleration, or 'channelling', would occur due to greatly enhanced oxygen solubility in lipids, especially near the midplane of lipid bilayers. Rapid long-range movement would be promoted by anisotropically enhanced lateral diffusion of oxygen along the midplane and by junctions holding lipid structures in close proximity to one another throughout the tissue. Clarifying the biophysical mechanism of oxygen transport within tissue will shed light on limitations and opportunities in tumour radiotherapy and tissue engineering.

Pathways of Oxygen Diffusion in Cells and Tissues : Hydrophobic Channeling via Networked Lipids

Oxygen delivery to tissue mitochondria relies on simple diffusion in the target cells and tissues. As such, intracellular availability of O₂ in tissue depends on its solubility and diffusivity in complex and heterogeneous macromolecular environments. The path of oxygen diffusion is key to its rate of transfer, especially where pathways of differing favorability are present. Most commonly, aqueous media, such as interstitial fluid and cytoplasm, are assumed to provide the dominant diffusion path. Here, the 'hydrophobic channeling' hypothesis is revisited, and several lines of evidence pointing toward lipid-accelerated oxygen diffusion pathways are discussed. The implications of hydrophobic channeling are considered in light of extended membrane networks in cells and tissues.

A novel tumor and mitochondria dual-targeted photosensitizer showing ultra-efficient photodynamic anticancer activities

We designed and synthesized a simple, but highly effective photosensitizer (G-Mito-Pc), which can precisely target the mitochondria of epidermal growth factor receptor (EGFR)-overexpressing cancer cells, to achieve dual targeting function at both cell and organelle levels in cancer therapy. We further explored the possible molecular mechanism of the enhanced bioactivity of G-Mito-Pc compared to that of the reference photosensitizer using molecular dynamics simulations on their interactions with a physiologically relevant mitochondrial membrane model.

Oxygen distribution in the fluid/gel phases of lipid membranes

The interactions between oxygen and lipid membranes play fundamental roles in basic biological processes (e.g., cellular respiration). Obviously, membrane oxidation is expected to be critically dependent on the distribution and concentration of oxygen in the membrane. Here, we combined theoretical and experimental methods to investigate oxygen partition and distribution in lipid membranes of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) in a temperature range between 298 and 323 K, specifically focusing on the changes caused by the lipid phase and phase transition. Even though oxygen is known to be more concentrated in the center of fluid phase membranes than on the headgroup regions, the distribution profile of oxygen inside gel-phase bilayers remained to be determined. Molecular dynamics simulations now show that the distribution of oxygen inside DPPC bilayers dramatically changes upon crossing the main transition temperature, with oxygen being nearly depleted halfway from the headgroups to the membrane center below the transition temperature. In a parallel approach, singlet oxygen luminescence emission measurements employing the photosensitizer Pheophorbide-a (Pheo) confirmed the differences in oxygen distribution and concentration profiles between gel- and fluid-phase membranes, revealing changes in the microenvironment of the embedded photosensitizer. Our results also reveal that excited triplet state lifetime, as it can be determined from the singlet oxygen luminescence kinetics, is a useful probe to assess oxygen distribution in lipid membranes with distinct lipid compositions.

Diffusion of nitric oxide and oxygen in lipoproteins and membranes studied by pyrene fluorescence quenching

Oxygen and nitric oxide are small hydrophobic molecules that usually need to diffuse a considerable distance to accomplish their biological functions and necessarily need to traverse several lipid membranes. Different methods have been used to study the diffusion of these molecules in membranes and herein we focus in the quenching of fluorescence of pyrenes inserted in the membrane. The pyrene derivatives have long fluorescence lifetimes (around 200 ns) that make them very sensitive to fluorescence quenching by nitric oxide, oxygen and other paramagnetic species. Results show that the apparent diffusion coefficients in membranes are similar to those in water, indicating that diffusion of these molecules in membranes is not considerably limited by the lipids. This high apparent diffusion in membranes is a consequence of both a favorable partition of these molecules in the hydrophobic interior of membranes and a high diffusion coefficient. Altering the composition of the membrane results in slight changes in diffusion, indicating that in most cases the lipid membranes will not hinder the passage of oxygen or nitric oxide. The diffusion of nitric oxide in the lipid core of low density lipoprotein is also very high, supporting its role as an antioxidant. In contrast to the high permeability of membranes to nitric oxide and oxygen, the permeability to other reactive species such as hydrogen peroxide and peroxynitrous acid is nearly five orders of magnitude lower.

Probing the interactions of phthalocyanine-based photosensitizers with model phospholipid bilayer by molecular dynamics simulations

Phthalocyanines (Pc) have received considerable attention in the design of photosensitizers for photodynamic therapy (PDT). It is of great interest to design novel Pc-based photosensitizer with improved biological efficiency, which largely relies on the understanding of the interactions of Pc with cell membranes at a molecular level. Here, via all-atom molecular dynamics simulations, we explored the interaction mechanism between a model phospholipid bilayer and three zinc Pc (ZnPc) molecules with different hydrophobicity in nature. We find that the adsorption and insertion behaviors of ZnPc molecules in the model bilayer are different due to the differing hydrophobicity and interaction patterns with phospholipids. Moreover, our simulations demonstrate that the conjunction of a ZnPc skeleton with a cholesterol moiety may reduce the intrinsic molecular rotation of ZnPc in membranes, presumably leading to an increase of the generation efficiency of reactive oxygen species for PDT. The molecular insights obtained here are likely to help improve the rational design of novel photosensitizers with enhanced cellular uptake and photocytotoxic activity.

Reactive Oxygen and Nitrogen Species at Phospholipid Bilayers: Peroxynitrous Acid and Its Homolysis Products

Peroxynitrite is a powerful and long-lived oxidant generated in vivo. Peroxynitrous acid (ONOOH), its protonated form, may penetrate into phospholipid bilayers and undergo homolytic cleavage to nitrogen dioxide (·NO₂) and hydroxyl radicals (·OH), causing severe nitro-oxidative damage. The membrane environment is thought to influence ONOOH reactions, but the mechanisms remain speculative. Most experimental techniques lack the level of resolution required to keep track of the motion of very reactive species and their interactions with the membrane. Here, we performed molecular dynamics simulations of the permeation, interactions, and dynamics of ONOOH and its homolysis products in the phospholipid membrane environment. We started by developing an ONOOH model that successfully accounted for its conformational equilibria and solvation energies. Membrane permeation of ONOOH was accompanied by conformational changes. ONOOH exhibited a strong tendency to bind to and accumulate at the membrane headgroup region. There, ONOOH homolysis led to ·NO₂ radicals, which in turn partitioned to the membrane interior. About one-third of the ·OH radicals readily escaped to the aqueous phase within 1 ns. However, a significant number of ·OH radicals became trapped at the lipid headgroup region for a longer period. The possible implications for membrane-based nitration and oxidation processes were discussed.

Reduced Oxygen Permeability upon Protein Incorporation Within Phospholipid Bilayers

Intracellular oxygenation is key to energy metabolism as well as tumor radiation therapy. Although integral proteins are ubiquitous in membranes, few studies have considered their effects on molecular oxygen permeability. Published experimental work with rhodopsin and bacteriorhodopsin has led to the hypothesis that integral proteins lessen membrane oxygen permeability, as well as the permeability of the lipid region. The current work uses atomistic molecular dynamics simulations to test the influence of an ungated potassium channel protein on the oxygen permeability of palmitoyloleoylphosphatidylcholine (POPC) bilayers with and without cholesterol. Consistent with experiment, whole-membrane oxygen permeability is cut in half upon adding 30 wt% potassium channel protein to POPC, and the apparent permeability of the lipid portion of the membrane decreases by 40%. Unexpectedly, oxygen is found to interact directly with the protein surface, accompanied by a 40% reduction of the apparent whole-membrane diffusion coefficient. Similar effects are seen in systems combining the potassium channel with 1:1 POPC/cholesterol, but the magnitude of permeability reduction is smaller by ~30%. Overall, the simulations indicate that integral proteins can reduce oxygen permeability by altering the diffusional path and the local diffusivity. This effect may be especially important in the protein-dense membranes of mitochondria.

Is the cholesterol bilayer domain a barrier to oxygen transport into the eye lens?

In the eye lens, the oxygen partial pressure is very low and the cholesterol (Chol) content in cell membranes is very high. Disturbance of these quantities results in cataract development. In human lens membranes, both bulk phospholipid-Chol domains and the pure Chol bilayer domains (CBDs) were experimentally detected. It is hypothesized that the CBD constitutes a significant barrier to oxygen transport into the lens. Transmembrane profiles of the oxygen diffusion-concentration product, obtained with electron paramagnetic resonance spin-labeling methods, allow evaluation of the oxygen permeability (P_M) of phospholipid membranes but not the CBD. Molecular dynamics simulation can independently provide components of the product across any bilayer domain, thus allowing evaluation of the P_M across the CBD. Therefore, to test the hypothesis, MD simulation was used. Three bilayers containing palmitoyl-oleoyl-phosphorylcholine (POPC) and Chol were built. The pure Chol bilayer modeled the CBD, the 1:1 POPC-Chol bilayer modeled the bulk membrane in which the CBD is embedded, and the POPC bilayer was a reference. To each model, 200 oxygen molecules were added. After equilibration, the oxygen concentration and diffusion profiles were calculated for each model and multiplied by each other. From the respective product profiles, the P_M of each bilayer was calculated. Favorable comparison with experimental data available only for the POPC and POPC-Chol bilayers validated these bilayer models and allowed the conclusion that oxygen permeation across the CBD is ~10 smaller than across the bulk membrane, supporting the hypothesis that the CBD is a barrier to oxygen transport into the eye lens.

Influence of Cholesterol on the Oxygen Permeability of Membranes: Insight from Atomistic Simulations

Cholesterol is widely known to alter the physical properties and permeability of membranes. Several prior works have implicated cell membrane cholesterol as a barrier to tissue oxygenation, yet a good deal remains to be explained with regard to the mechanism and magnitude of the effect. We use molecular dynamics simulations to provide atomic-resolution insight into the influence of cholesterol on oxygen diffusion across and within the membrane. Our simulations show strong overall agreement with published experimental data, reproducing the shapes of experimental oximetry curves with high accuracy. We calculate the upper-limit transmembrane oxygen permeability of a 1-palmitoyl,2-oleoylphosphatidylcholine phospholipid bilayer to be 52 ± 2 cm/s, close to the permeability of a water layer of the same thickness. With addition of cholesterol, the permeability decreases somewhat, reaching 40 ± 2 cm/s at the near-saturating level of 62.5 mol % cholesterol and 10 ± 2 cm/s in a 100% cholesterol mimic of the experimentally observed noncrystalline cholesterol bilayer domain. These reductions in permeability can only be biologically consequential in contexts where the diffusional path of oxygen is not water dominated. In our simulations, cholesterol reduces the overall solubility of oxygen within the membrane but enhances the oxygen transport parameter (solubility-diffusion product) near the membrane center. Given relatively low barriers to passing from membrane to membrane, our findings support hydrophobic channeling within membranes as a means of cellular and tissue-level oxygen transport. In such a membrane-dominated diffusional scheme, the influence of cholesterol on oxygen permeability is large enough to warrant further attention.

Position-Dependent Diffusion Tensors in Anisotropic Media from Simulation: Oxygen Transport in and through Membranes

A Bayesian-based methodology is developed to estimate diffusion tensors from molecular dynamics simulations of permeants in anisotropic media, and is applied to oxygen in lipid bilayers. By a separation of variables in the Smoluchowski diffusion equation, the multidimensional diffusion is reduced to coupled one-dimensional diffusion problems that are treated by discretization. The resulting diffusivity profiles characterize the membrane transport dynamics as a function of the position across the membrane, discriminating between diffusion normal and parallel to the membrane. The methodology is first validated with neat water, neat hexadecane, and a hexadecane slab surrounded by water, the latter being a simple model for a lipid membrane. Next, a bilayer consisting of pure 1-palmitoyl 2-oleoylphosphatidylcholine (POPC), and a bilayer mimicking the lipid composition of the inner mitochondrial membrane, including cardiolipin, are investigated. We analyze the detailed time evolution of oxygen molecules, in terms of both normal diffusion through and radial diffusion inside the membrane. Diffusion is fast in the more loosely packed interleaflet region, and anisotropic, with oxygen spreading more rapidly in the membrane plane than normal to it. Visualization of the propagator shows that oxygen enters the membrane rapidly, reaching its thermodynamically favored center in about 1 ns, despite the free energy barrier at the headgroup region. Oxygen transport is quantified by computing the oxygen permeability of the membranes and the average radial diffusivity, which confirm the anisotropy of the diffusion. The position-dependent diffusion constants and free energies are used to construct compartmental models and test assumptions used in estimating permeability, including Overton's rule. In particular, a hexadecane slab surrounded by water is found to be a poor model of oxygen transport in membranes because the relevant energy barriers differ substantially.

Solubility and diffusion of oxygen in phospholipid membranes

The transport of oxygen and other nonelectrolytes across lipid membranes is known to depend on both diffusion and solubility in the bilayer, and to be affected by changes in the physical state and by the lipid composition, especially the content of cholesterol and unsaturated fatty acids. However, it is not known how these factors affect diffusion and solubility separately. Herein we measured the partition coefficient of oxygen in liposome membranes of dilauroyl-, dimiristoyl- and dipalmitoylphosphatidylcholine in buffer at different temperatures using the equilibrium-shift method with electrochemical detection. The apparent diffusion coefficient was measured following the fluorescence quenching of 1-pyrenedodecanoate inserted in the liposome bilayers under the same conditions. The partition coefficient varied with the temperature and the physical state of the membrane, from below 1 in the gel state to above 2.8 in the liquid-crystalline state in DMPC and DPPC membranes. The partition coefficient was directly proportional to the partial molar volume and was then associated to the increase in free-volume in the membrane as a function of temperature. The apparent diffusion coefficients were corrected by the partition coefficients and found to be nearly the same, with a null dependence on viscosity and physical state of the membrane, probably because the pyrene is disturbing the surrounding lipids and thus becoming insensitive to changes in membrane viscosity. Combining our results with those of others, it is apparent that both solubility and diffusion increase when increasing the temperature or when comparing a membrane in the gel to one in the fluid state.

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Reactive oxygen species at phospholipid bilayers: distribution, mobility and permeation

Reactive oxygen species (ROS) are involved in biochemical processes such as redox signaling, aging, carcinogenesis and neurodegeneration. Although biomembranes are targets for reactive oxygen species attack, little is known about the role of their specific interactions. Here, molecular dynamics simulations were employed to determine the distribution, mobility and residence times of various reactive oxygen species at the membrane-water interface. Simulations showed that molecular oxygen (O2) accumulated at the membrane interior. The applicability of this result to singlet oxygen ((1)O2) was discussed. Conversely, superoxide (O2(-)) radicals and hydrogen peroxide (H2O2) remained at the aqueous phase. Both hydroxyl (HO) and hydroperoxyl (HO2) radicals were able to penetrate deep into the lipid headgroups region. Due to membrane fluidity and disorder, these radicals had access to potential peroxidation sites along the lipid hydrocarbon chains, without having to overcome the permeation free energy barrier. Strikingly, HO2 radicals were an order of magnitude more concentrated in the headgroups region than in water, implying a large shift in the acid-base equilibrium between HO2 and O2(-). In comparison with O2, both HO and HO2 radicals had lower lateral mobility at the membrane. Simulations revealed that there were intermittent interruptions in the H-bond network around the HO radicals at the headgroups region. This effect is expected to be unfavorable for the H-transfer mechanism involved in HO diffusion. The implications for lipid peroxidation and for the effectiveness of membrane antioxidants were evaluated.

Topology and immersion depth of an integral membrane protein by paramagnetic rates from dissolved oxygen

In studies of membrane proteins, knowledge of protein topology can provide useful insight into both structure and function. In this work, we present a solution NMR method for the measurement the tilt angle and average immersion depth of alpha helices in membrane proteins, from analysis of the paramagnetic relaxation rate enhancements arising from dissolved oxygen. No modification to the micelle or protein is necessary, and the topology of both transmembrane and amphipathic helices are readily determined. We apply this method to the measure the topology of a monomeric mutant of phospholamban (AFA-PLN), a 52-residue membrane protein containing both an amphipathic and a transmembrane alpha helix. In dodecylphosphocholine micelles, the amphipathic helix of AFA-PLN was found to have a tilt angle of 87° ± 1° and an average immersion depth of 13.2 Å. The transmembrane helix was found to have an average immersion depth of 5.4 Å, indicating residues 41 and 42 are closest to the micelle centre. The resolution of paramagnetic relaxation rate enhancements from dissolved oxygen compares favourably to those from Ni (II), a hydrophilic paramagnetic species.

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)).

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.

An NMR method for the determination of protein binding interfaces using TEMPOL-induced chemical shift perturbations

BACKGROUND

The determination of protein-protein interfaces is of crucial importance to understand protein function and to guide the design of compounds. To identify protein-protein interface by NMR spectroscopy, 13C NMR paramagnetic shifts induced by freely diffusing 4-hydroxy-2, 2, 6, 6-tetramethyl-piperidine-1-oxyl (TEMPOL) are promising, because TEMPOL affects distinct 13C NMR chemical shifts of the solvent accessible nuclei belonging to proteins of interest, while 13C nuclei within the interior of the proteins may be distinguished by a lack of such shifts.

METHOD

We measured the 13C NMR paramagnetic shifts induced by TEMPOL by recording 13C-(13)C TOCSY spectra for ubiquitin in the free state and the complex state with yeast ubiquitin hydrolase1 (YUH1).

RESULTS

Upon complexation of ubiquitin with YUH1, 13C NMR paramagnetic shifts associated with the protein binding interface were reduced by 0.05 ppm or more. The identified interfacial atoms agreed with the prior X-ray crystallographic data.

CONCLUSIONS

The TEMPOL-induced 13C chemical shift perturbation is useful to determine precise protein-protein interfaces.

GENERAL SIGNIFICANCE

The present method is a useful method to determine protein-protein interface by NMR, because it has advantages in easy sample preparations, simple data analyses, and wide applicabilities.

Generation of superoxide and singlet oxygen from alpha-tocopherolquinone and analogues

Three potential routes to generation of reactive oxygen species (ROS) from alpha-tocopherolquinone (alpha-TQ) have been identified. The quinone of the water-soluble vitamin E analogue Trolox C (Trol-Q) is reduced by hydrated electron and isopropanol alpha-hydroxyalkyl radical, and the resulting semiquinone reacts with molecular oxygen to form superoxide with a second order rate constant of 1.3 x 10(8) dm(3)/mol/s, illustrating the potential for redox cycling. Illumination (UV-A, 355 nm) of the quinone of 2,2,5,7,8-pentamethyl-6-hydroxychromanol (PMHC-Q) leads to a reactive short-lived (ca. 10(- 6) s) triplet state, able to oxidise tryptophan with a second order rate constant greater than 10(9) dm(3)/mol/s. The triplet states of these quinones sensitize singlet oxygen formation with quantum yields of about 0.8. Such potentially damaging reactions of alpha-TQ may in part account for the recent findings that high levels of dietary vitamin E supplementation lack any beneficial effect and may lead to slightly enhanced levels of overall mortality.