Oxygen profiles in membranes

Architecture of the pore forming toxin sticholysin I in membranes

Sticholysin I (StI) is a toxin produced by the sea anemone Stichodactyla helianthus and belonging to the actinoporins family. Upon binding to sphingomyelin-containing membranes StI forms oligomeric pores, thereby leading to cell death. According to recent controversial experimental evidences, the pore architecture of actinoporins is a debated topic. Here, we investigated the StI topology in membranes by site-directed spin labeling and electron paramagnetic resonance spectroscopy. The results reveal that StI in membrane exhibits an oligomeric architecture with heterogeneous stoichiometry of predominantly eight or nine protomers, according to the available structural models. The StI topology resembles the conic pore structure reported for the actinoporin fragaceatoxin C. Our data show that StI coexists in two membrane-associated conformations, with the N-terminal segment either attached to the protein core or inserted in the membrane forming the pore. This finding suggests a 'pre-pore' to 'pore' transition determined by a conformational change that detaches the N-terminal segment.

Bacterial denitrifying nitric oxide reductases and aerobic respiratory terminal oxidases use similar delivery pathways for their molecular substrates

The superfamily of heme‑copper oxidoreductases (HCOs) include both NO and O₂ reductases. Nitric oxide reductases (NORs) are bacterial membrane enzymes that catalyze an intermediate step of denitrification by reducing nitric oxide (NO) to nitrous oxide (N₂O). They are structurally similar to heme‑copper oxygen reductases (HCOs), which reduce O₂ to water. The experimentally observed apparent bimolecular rate constant of NO delivery to the deeply buried catalytic site of NORs was previously reported to approach the diffusion-controlled limit (10⁸-10⁹ M^-1 s^-1). Using the crystal structure of cytochrome-c dependent NOR (cNOR) from Pseudomonas aeruginosa, we employed several protocols of molecular dynamics (MD) simulation, which include flooding simulations of NO molecules, implicit ligand sampling and umbrella sampling simulations, to elucidate how NO in solution accesses the catalytic site of this cNOR. The results show that NO partitions into the membrane, enters the enzyme from the lipid bilayer and diffuses to the catalytic site via a hydrophobic tunnel that is resolved in the crystal structures. This is similar to what has been found for O₂ diffusion through the closely related O₂ reductases. The apparent second order rate constant approximated using the simulation data is ~5 × 10⁸ M^-1 s^-1, which is optimized by the dynamics of the amino acid side chains lining in the tunnel. It is concluded that both NO and O₂ reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.

The Synergetic Effects of Combining Structural Biology and EPR Spectroscopy on Membrane Proteins

Protein structures as provided by structural biology such as X-ray crystallography, cryo-electron microscopy and NMR spectroscopy are key elements to understand the function of a protein on the molecular level. Nonetheless, they might be error-prone due to crystallization artifacts or, in particular in case of membrane-imbedded proteins, a mostly artificial environment. In this review, we will introduce different EPR spectroscopy methods as powerful tools to complement and validate structural data gaining insights in the dynamics of proteins and protein complexes such that functional cycles can be derived. We will highlight the use of EPR spectroscopy on membrane-embedded proteins and protein complexes ranging from receptors to secondary active transporters as structural information is still limited in this field and the lipid environment is a particular challenge.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Voltage-sensitive styryl dyes as singlet oxygen targets on the surface of bilayer lipid membrane

Photosensitizers are widely used as photodynamic therapeutic agents killing cancer cells by photooxidation of their components. Development of new effective photosensitive molecules requires profound knowledge of possible targets for reactive oxygen species, especially for its singlet form. Here we studied photooxidation of voltage-sensitive styryl dyes (di-4-ANEPPS, di-8-ANEPPS, RH-421 and RH-237) by singlet oxygen on the surface of bilayer lipid membranes commonly used as cell membrane models. Oxidation was induced by irradiation of a photosensitizer (aluminum phthalocyanine tetrasulfonate) and monitored by the change of dipole potential on the surface of the membrane. We studied the drop of the dipole potential both in the case when the dye molecules were adsorbed on the same side of the lipid bilayer as the photosensitizer (cis-configuration) and in the case when they were adsorbed on the opposite side (trans-configuration). Based on a simple model, we determined the rate of oxidation of the dyes from the kinetics of change of the potential during and after irradiation. This rate is proportional to steady-state concentration of singlet oxygen in the membrane under irradiation. Comparison of the oxidation rates of various dyes reveals that compounds of ANEPPS series are more sensitive to singlet oxygen than RH type dyes, indicating that naphthalene group is primarily responsible for their oxidation.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

Water penetration profile at the protein-lipid interface in Na,K-ATPase membranes

The affinity of ionized fatty acids for the Na,K-ATPase is used to determine the transmembrane profile of water penetration at the protein-lipid interface. The standardized intensity of the electron spin echo envelope modulation (ESEEM) from (2)H-hyperfine interaction with D2O is determined for stearic acid, n-SASL, spin-labeled systematically at the C-n atoms throughout the chain. In both native Na,K-ATPase membranes from shark salt gland and bilayers of the extracted membrane lipids, the D2O-ESEEM intensities of fully charged n-SASL decrease progressively with position down the fatty acid chain toward the terminal methyl group. Whereas the D2O intensities decrease sharply at the n = 9 position in the lipid bilayers, a much broader transition region in the range n = 6 to 10 is found with Na,K-ATPase membranes. Correction for the bilayer population in the membranes yields the intrinsic D2O-intensity profile at the protein-lipid interface. For positions at either end of the chains, the D2O concentrations at the protein interface are greater than in the lipid bilayer, and the positional profile is much broader. This reveals the higher polarity, and consequently higher intramembrane water concentration, at the protein-lipid interface. In particular, there is a significant water concentration adjacent to the protein at the membrane midplane, unlike the situation in the bilayer regions of this cholesterol-rich membrane. Experiments with protonated fatty acid and phosphatidylcholine spin labels, both of which have a considerably lower affinity for the Na,K-ATPase, confirm these results.

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

Membrane docking of the C2 domain from protein kinase Cα as seen by polarized ATR-IR. The role of PIP₂

We have used attenuated total internal reflection infrared spectroscopy (ATR-IR) spectroscopy to study the association of the C2 domain from protein kinase Cα (PKCα) with different phospholipid membranes, so as to characterise the mode of membrane docking and its modulation by the second-messenger lipid PIP₂. In parallel, we have also examined the membrane interaction of the C2 domain from cytosolic phospholipase A₂. PIP₂ did not induce significant changes in secondary structure of the membrane-bound PKCα-C2 domain, nor did binding of the PKCα-C2 domain change the dichroic ratios of the lipid chains, whereas the C2 domain from phospholipase A₂ did perturb the lipid chain orientation. Measurements of the dichroic ratios for the amide I and amide II protein bands were combined so as to distinguish the tilt of the β-sheets from that of the β-strands within the sheet. When associated with POPC/POPS membranes, the β-sandwich of the PKCα-C2 domain is inclined at an angle α=35° to the membrane normal, i.e., is oriented more nearly perpendicular than parallel to the membrane. In the process of membrane docking, the tilt angle increases to α=44° in the presence of PIP₂, indicating that the β-sandwich comes closer to the membrane surface, so confirming the importance of this lipid in determining docking of the C2 domain and consequent activation of PKCα.

Spin-Label EPR for Determining Polarity and Proticity in Biomolecular Assemblies: Transmembrane Profiles

Hyperfine couplings and g-values of nitroxyl spin labels are sensitive to polarity and hydrogen bonding in the environment probed. The dependences of these electronic paramagnetic resonance (EPR) properties on environmental dielectric permittivity and proticity are reviewed. Calibrations are given, in terms of the Block-Walker reaction field and local proton donor concentration, for the nitroxides that are commonly used in spin labeling of lipids and proteins. Applications to studies of the transverse polarity profiles in lipid bilayers, which constitute the permeability barrier of biological membranes, are reviewed. Emphasis is given to parallels with the permeation profiles of oxygen and nitric oxide that are determined from spin-label relaxation enhancements by using nonlinear continuous-wave EPR and saturation recovery EPR, and with permeation profiles of D(2)O that are determined by using (2)H electron spin echo envelope modulation spectroscopy.

Spin labeling EPR

Site-directed spin labeling in combination with electron paramagnetic resonance spectroscopy has emerged as an efficient tool to elucidate the structure and conformational dynamics of biomolecules under native-like conditions. This article summarizes the basics as well as recent progress of site-directed spin labeling. Continuous wave EPR spectra analyses and pulse EPR techniques are reviewed with special emphasis on applications to the sensory rhodopsin-transducer complex mediating the photophobic response of the halophilic archaeum Natronomonas pharaonis and the photosynthetic reaction center from Rhodobacter sphaeroides R26.

Reaction fields in the environment of fluorescent probes: polarity profiles in membranes

Fluorescent probes in biological systems are sensitive to environmental polarity by virtue of their response to the reaction field created by polarization of the dielectric medium. Classically, fluorophore solvatochromism is analyzed in terms of the Lippert equation and later variants, all of which rely upon the original reaction field of Onsager. A recent survey of the solvent dependence of EPR spin-label probes, which are responsive solely to the reaction field in the ground state without the complication of excited states, shows that the reaction field of Block and Walker performs best in describing the polarity dependence. In this model, the step-function transition to the bulk dielectric medium used by Onsager is replaced by a graded transition. Analysis of the Stokes shifts for representative fluorescent membrane probes, such as PRODAN, DANSYL, and anthroyl fatty acid, reveals that, of several different reaction fields (including that of Onsager), the Block-Walker model best describes the dependence on solvent dielectric constant and refractive index for the different probes simultaneously. This is after full allowance is made for all contributions involving polarizability of the fluorophore, a point that is frequently neglected or treated incorrectly in studies using biological fluorescent probes. By using the full range of polar and apolar solvents, it is then possible to establish a common reference for the polarity dependence of different fluorophores and to relate this also to the polarity dependence of biologically relevant spin-label EPR probes. An important application is calibration of the transmembrane polarity profile recorded by fluorescent probes in terms of the high-resolution profile obtained from site-specifically spin-labeled lipid chains.

Intramembrane water associated with TOAC spin-labeled alamethicin: electron spin-echo envelope modulation by D2O

Alamethicin is a 20-residue, hydrophobic, helical peptide, which forms voltage-sensitive ion channels in lipid membranes. The helicogenic, nitroxyl amino acid TOAC was substituted isosterically for Aib at residue positions 1, 8, or 16 in a F50/5 alamethicin analog to enable EPR studies. Electron spin-echo envelope modulation (ESEEM) spectroscopy was used to investigate the water exposure of TOAC-alamethicin introduced into membranes of saturated or unsaturated diacyl phosphatidylcholines that were dispersed in D2O. Echo-detected EPR spectra were used to assess the degree of assembly of the peptide in the membrane, via the instantaneous diffusion from intermolecular spin-spin interactions. The profile of residue exposure to water differs between membranes of saturated and unsaturated lipids. In monounsaturated dioleoyl phosphatidylcholine, D2O-ESEEM intensities decrease from TOAC(1) to TOAC(8) and TOAC(16) but not uniformly. This is consistent with a transmembrane orientation for the protoassembled state, in which TOAC(16) is located in the bilayer leaflet opposite to that of TOAC(1) and TOAC(8). Relative to the monomer in fluid bilayers, assembled alamethicin is disposed asymmetrically about the bilayer midplane. In saturated dimyristoyl phosphatidylcholine, the D2O-ESEEM intensity is greatest for TOAC(8), indicating a more superficial location for alamethicin, which correlates with the difference in orientation between gel- and fluid-phase membranes found by conventional EPR of TOAC-alamethicin in aligned phosphatidylcholine bilayers. Increasing alamethicin/lipid ratio in saturated phosphatidylcholine shifts the profile of water exposure toward that with unsaturated lipid, consistent with proposals of a critical concentration for switching between the two different membrane-associated states.

Conformational changes and protein stability of the pro-apoptotic protein Bax

Pro-apoptotic Bax is a soluble and monomeric protein under normal physiological conditions. Upon its activation substantial structural rearrangements occur: The protein inserts into the mitochondrial outer membrane and forms higher molecular weight oligomers. Subsequently, the cells can undergo apoptosis. In our studies, we focused on the structural rearrangements of Bax during oligomerization and on the protein stability. Both protein conformations exhibit high stability against thermal denaturation, chemically induced unfolding and proteolytic processing. The oligomeric protein is stable up to 90 degrees C as well as in solutions of 8 M urea or 6 M guanidinium hydrochloride. Helix 9 appears accessible in the monomer but hidden in the oligomer assessed by proteolysis. Tryptophan fluorescence indicates that the environment of the C-terminal protein half becomes more apolar upon oligomerization, whereas the loop region between helices 1 and 2 gets solvent exposed.

Spin-diffusion couples proton relaxation rates for proteins in exchange with a membrane interface

Changes in nuclear spin-lattice relaxation rates that are induced by a freely diffusing paramagnetic relaxation agent are examined for a protein in solution and compared to the case where the protein binds to a membrane. In the solution case, the intramolecular cross-relaxation rates are modest and large differences are observed in the oxygen induced protein-proton relaxation rates. In the case where a dynamic equilibrium between solution and membrane-bound environments is established, the intramolecular (1)H cross-relaxation rates for the protein protons increase dramatically because of the slow reorientational motion in the membrane-bound environment. As a consequence, all protein protons relax with nearly the same spin-lattice relaxation rate constants when bound to the membrane, and site specific relaxation effects of the diffusing paramagnet are suppressed. Slowly reorienting sites or rotationally immobilized sites sampled by observable molecules in vivo will demonstrate similar relaxation leveling effects.

Crystallographic studies of Xe and Kr binding within the large internal cavity of cytochrome ba3 from Thermus thermophilus: structural analysis and role of oxygen transport channels in the heme-Cu oxidases

Cytochrome ba3 is a cytochrome c oxidase from the plasma membrane of Thermus thermophilus and is the preferred terminal enzyme of cellular respiration at low dioxygen tensions. Using cytochrome ba 3 crystals pressurized at varying conditions under Xe or Kr gas, and X-ray data for six crystals, we identify the relative affinities of Xe and Kr atoms for as many as seven distinct binding sites. These sites track a continuous, Y-shaped channel, 18-20 A in length, lined by hydrophobic residues, which leads from the surface of the protein where two entrance holes, representing the top of the Y, connect the bilayer to the a3-CuB center at the base of the Y. Considering the increased affinity of O2 for hydrophobic environments, the hydrophobic nature of the channel, its orientation within the bilayer, its connection to the active site, its uniform diameter, its virtually complete occupation by Xe, and its isomorphous presence in the native enzyme, we infer that the channel is a diffusion pathway for O2 into the dinuclear center of cytochrome ba3. These observations provide a basis for analyzing similar channels in other oxidases of known structure, and these structures are discussed in terms of mechanisms of O2 transport in biological systems, details of CO binding to and egress from the dinuclear center, the bifurcation of the oxygen-in and water-out pathways, and the possible role of the oxygen channel in aerobic thermophily.

Solvent dependence of the rotational diffusion of TOAC-spin-labeled alamethicin

Three derivatives of the hydrophobic, channel-forming peptaibiotic alamethicin (F50/5) have been synthesized, the original Aib residue at position 1, 8, or 16 being replaced with the spin-labeled amino acid TOAC (=2,2,6,6-tetramethylpiperidin-1-oxyl-4-amino-4-carboxylic acid). Electron-paramagnetic-resonance (EPR) spectroscopy was used to characterize the rotational diffusion of these compounds in five isotropic solvents of differing viscosity and polarity, including MeOH, EtOH, PrOH, i-PrOH, and hexanol (HxOH). In MeOH, the labeled alamethicins were found to rotate anisotropically as a monomer (axial ratio a/b=3). In aliphatic alcohols of increasing viscosity (and hydrophobicity), the rotational correlation times progressively increased. Even in HxOH, the (fivefold) increase in correlation time was no greater than the increase in viscosity. We conclude that TOAC-labeled alamethicins remain monomeric in these solvents of relatively high polarity.