Mobility of Xe atoms within the oxygen diffusion channel of cytochrome ba(3) oxidase

High-pressure crystallography shows noble gas intervention into protein-lipid interaction and suggests a model for anaesthetic action

In this work we examine how small hydrophobic molecules such as inert gases interact with membrane proteins (MPs) at a molecular level. High pressure atmospheres of argon and krypton were used to produce noble gas derivatives of crystals of three well studied MPs (two different proton pumps and a sodium light-driven ion pump). The structures obtained using X-ray crystallography showed that the vast majority of argon and krypton binding sites were located on the outer hydrophobic surface of the MPs – a surface usually accommodating hydrophobic chains of annular lipids (which are known structural and functional determinants for MPs). In conformity with these results, supplementary in silico molecular dynamics (MD) analysis predicted even greater numbers of argon and krypton binding positions on MP surface within the bilayer. These results indicate a potential importance of such interactions, particularly as related to the phenomenon of noble gas-induced anaesthesia.

Noble gases are known to interact with proteins and can be good anaesthetics in hyperbaric conditions. This study identifies argon and krypton binding sites on membrane proteins and proposes as a hypothesis that noble gases, by altering protein/lipid contacts, may affect protein function.

Comparative study of the effects of high hydrostatic pressure per se and high argon pressure on urate oxidase ligand stabilization

The stability of the tetrameric enzyme urate oxidase in complex with excess of 8-azaxanthine was investigated either under high hydrostatic pressure per se or under a high pressure of argon. The active site is located at the interface of two subunits, and the catalytic activity is directly related to the integrity of the tetramer. This study demonstrates that applying pressure to a protein–ligand complex drives the thermodynamic equilibrium towards ligand saturation of the complex, revealing a new binding site. A transient dimeric intermediate that occurs during the pressure-induced dissociation process was characterized under argon pressure and excited substates of the enzyme that occur during the catalytic cycle can be trapped by pressure. Comparison of the different structures under pressure infers an allosteric role of the internal hydrophobic cavity in which argon is bound, since this cavity provides the necessary flexibility for the active site to function.

Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling

Heme-copper oxidases (HCO), nitric oxide reductases (NOR), and sulfite reductases (SiR) catalyze the multi-electron and multi-proton reductions of O2, NO, and SO32-, respectively. Each of these reactions is important to drive cellular energy production through respiratory metabolism and HCO, NOR, and SiR evolved to contain heteronuclear active sites containing heme/copper, heme/nonheme iron, and heme-[4Fe-4S] centers, respectively. The complexity of the structures and reactions of these native enzymes, along with their large sizes and/or membrane associations, make it challenging to fully understand the crucial structural features responsible for the catalytic properties of these active sites. In this review, we summarize progress that has been made to better understand these heteronuclear metalloenzymes at the molecular level though study of the native enzymes along with insights gained from biomimetic models comprising either small molecules or proteins. Further understanding the reaction selectivity of these enzymes is discussed through comparisons of their similar heteronuclear active sites, and we offer outlook for further investigations.

Mechanisms of Energy Transduction by Charge Translocating Membrane Proteins

Life relies on the constant exchange of different forms of energy, i.e., on energy transduction. Therefore, organisms have evolved in a way to be able to harvest the energy made available by external sources (such as light or chemical compounds) and convert these into biological useable energy forms, such as the transmembrane difference of electrochemical potential (Δμ̃). Membrane proteins contribute to the establishment of Δμ̃ by coupling exergonic catalytic reactions to the translocation of charges (electrons/ions) across the membrane. Irrespectively of the energy source and consequent type of reaction, all charge-translocating proteins follow two molecular coupling mechanisms: direct- or indirect-coupling, depending on whether the translocated charge is involved in the driving reaction. In this review, we explore these two coupling mechanisms by thoroughly examining the different types of charge-translocating membrane proteins. For each protein, we analyze the respective reaction thermodynamics, electron transfer/catalytic processes, charge-translocating pathways, and ion/substrate stoichiometries.

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.

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.

Cytochrome aa3 Oxygen Reductase Utilizes the Tunnel Observed in the Crystal Structures To Deliver O2 for Catalysis

Cytochrome aa₃ is the terminal respiratory enzyme of all eukaryotes and many bacteria and archaea, reducing O₂ to water and harnessing the free energy from the reaction to generate the transmembrane electrochemical potential. The diffusion of O₂ to the heme-copper catalytic site, which is buried deep inside the enzyme, is the initiation step of the reaction chemistry. Our previous molecular dynamics (MD) study with cytochrome ba₃, a homologous enzyme of cytochrome aa₃ in Thermus thermophilus, demonstrated that O₂ diffuses from the lipid bilayer to its reduction site through a 25 Å long tunnel inferred by Xe binding sites detected by X-ray crystallography [Mahinthichaichan, P., Gennis, R., and Tajkhorshid, E. (2016) Biochemistry 55, 1265-1278]. Although a similar tunnel is observed in cytochrome aa₃, this putative pathway appears partially occluded between the entrances and the reduction site. Also, the experimentally determined second-order rate constant for O₂ delivery in cytochrome aa₃ (∼10⁸ M^-1 s^-1) is 10 times slower than that in cytochrome ba₃ (∼10⁹ M^-1 s^-1). A question to be addressed is whether cytochrome aa₃ utilizes this X-ray-inferred tunnel as the primary pathway for O₂ delivery. Using complementary computational methods, including multiple independent flooding MD simulations and implicit ligand sampling calculations, we probe the O₂ delivery pathways in cytochrome aa₃ of Rhodobacter sphaeroides. All of the O₂ molecules that arrived in the reduction site during the simulations were found to diffuse through the X-ray-observed tunnel, despite its apparent constriction, supporting its role as the main O₂ delivery pathway in cytochrome aa₃. The rate constant for O₂ delivery in cytochrome aa₃, approximated using the simulation results, is 10 times slower than in cytochrome ba₃, in agreement with the experimentally determined rate constants.

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N

In heme proteins, the efficient transport of ligands such as NO or O₂ to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it.

Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.

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.

Exploring membrane respiratory chains

Acquisition of energy is central to life. In addition to the synthesis of ATP, organisms need energy for the establishment and maintenance of a transmembrane difference in electrochemical potential, in order to import and export metabolites or to their motility. The membrane potential is established by a variety of membrane bound respiratory complexes. In this work we explored the diversity of membrane respiratory chains and the presence of the different enzyme complexes in the several phyla of life. We performed taxonomic profiles of the several membrane bound respiratory proteins and complexes evaluating the presence of their respective coding genes in all species deposited in KEGG database. We evaluated 26 quinone reductases, 5 quinol:electron carriers oxidoreductases and 18 terminal electron acceptor reductases. We further included in the analyses enzymes performing redox or decarboxylation driven ion translocation, ATP synthase and transhydrogenase and we also investigated the electron carriers that perform functional connection between the membrane complexes, quinones or soluble proteins. Our results bring a novel, broad and integrated perspective of membrane bound respiratory complexes and thus of the several energetic metabolisms of living systems. This article is part of a Special Issue entitled 'EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2-6, 2016', edited by Prof. Paolo Bernardi.

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.

Oxygen diffusion pathways in a cofactor-independent dioxygenase

A combination of experimental and computational studies reveals the main O₂ diffusion pathways, providing insight into how cofactor-independent oxygenases control stereospecificity and prevent oxidative inactivation.

Molecular oxygen plays an important role in a wide variety of enzymatic reactions. Through recent research efforts combining computational and experimental methods a new view of O₂ diffusion is emerging, where specific channels guide O₂ to the active site. The focus of this work is DpgC, a cofactor-independent oxygenase. Molecular dynamics simulations, together with mutagenesis experiments and xenon-binding data, reveal that O₂ reaches the active site of this enzyme using three main pathways and four different access points. These pathways connect a series of dynamic hydrophobic pockets, concentrating O₂ at a specific face of the enzyme substrate. Extensive molecular dynamics simulations provide information about which pathways are more frequently used. This data is consistent with the results of kinetic measurements on mutants and is difficult to obtain using computational cavity-location methods. Taken together, our results reveal that although DpgC is rare in its ability of activating O₂ in the absence of cofactors or metals, the way O₂ reaches the active site is similar to that reported for other O₂-using proteins: multiple access channels are available, and the architecture of the pathway network can provide regio- and stereoselectivity. Our results point to the existence of common themes in O₂ access that are conserved among very different types of proteins.

CO-dynamics in the active site of cytochrome c oxidase

The transfer of CO from heme a3 to the Cu(B) site in Cytochrome c oxidase (CcO) after photolysis is studied using molecular dynamics simulations using an explicitly reactive, parametrized potential energy surface based on density functional theory calculations. After photodissociation from the heme-Fe, the CO ligand rebinds to the Cu(B) site on the sub-picosecond time scale. Depending on the simulation protocol the characteristic time ranges from 260 fs to 380 fs which compares with an estimated 450 fs from experiment based on the analysis of the spectral changes as a function of time delay after the photodissociating pulse. Following photoexcitation ≈90% of the ligands are found to rebind to either the Cu(B) (major component, 85%) or the heme-Fe (minor component, 2%) whereas about 10% remain in an unbound state. The infrared spectra of unbound CO in the active site is broad and featureless and no appreciable shift relative to gas-phase CO is found, which is in contrast to the situation in myoglobin. These observations explain why experimentally, unbound CO in the binuclear site of CcO has not been found as yet.

Substrate pathways in the nitrogenase MoFe protein by experimental identification of small molecule binding sites

In the nitrogenase molybdenum-iron (MoFe) protein, we have identified five potential substrate access pathways from the protein surface to the FeMo-cofactor (the active site) or the P-cluster using experimental structures of Xe pressurized into MoFe protein crystals from Azotobacter vinelandii and Clostridium pasteurianum. Additionally, all published structures of the MoFe protein, including those from Klebsiella pneumoniae, were analyzed for the presence of nonwater, small molecules bound to the protein interior. Each pathway is based on identification of plausible routes from buried small molecule binding sites to both the protein surface and a metallocluster. Of these five pathways, two have been previously suggested as substrate access pathways. While the small molecule binding sites are not conserved among the three species of MoFe protein, residues lining the pathways are generally conserved, indicating that the proposed pathways may be accessible in all three species. These observations imply that there is unlikely a unique pathway utilized for substrate access from the protein surface to the active site; however, there may be preferred pathways such as those described here.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

Conserved glycine 232 in the ligand channel of ba3 cytochrome oxidase from Thermus thermophilus

Knowing how the protein environment modulates ligand pathways and redox centers in the respiratory heme-copper oxidases is fundamental for understanding the relationship between the structure and function of these enzymes. In this study, we investigated the reactions of O2 and NO with the fully reduced G232V mutant of ba3 cytochrome c oxidase from Thermus thermophilus (Tt ba3) in which a conserved glycine residue in the O2 channel of the enzyme was replaced with a bulkier valine residue. Previous studies of the homologous mutant of Rhodobacter sphaeroides aa3 cytochrome c oxidase suggested that the valine completely blocked the access of O2 to the active site [Salomonsson, L., et al. (2004) Proc. Natl. Acad. Sci. U.S.A. 101, 11617-11621]. Using photolabile O2 and NO carriers, we find by using time-resolved optical absorption spectroscopy that the rates of O2 and NO binding are not significantly affected in the Tt ba3 G232V mutant. Classical molecular dynamics simulations of diffusion of O2 to the active site in the wild-type enzyme and G232V mutant show that the insertion of the larger valine residue in place of the glycine appears to open up other O2 and NO exit/entrance pathways that allow these ligands unhindered access to the active site, thus compensating for the larger valine residue.

Density functional study for the bridged dinuclear center based on a high-resolution X-ray crystal structure of ba3 cytochrome c oxidase from Thermus thermophilus

Strong electron density for a peroxide type dioxygen species bridging the Fea3 and CuB dinuclear center (DNC) was observed in the high-resolution (1.8 Å) X-ray crystal structures (PDB entries 3S8G and 3S8F) of ba3 cytochrome c oxidase (CcO) from Thermus thermophilus. The crystals represent the as-isolated X-ray photoreduced CcO structures. The bridging peroxide was proposed to arise from the recombination of two radiation-produced HO(•) radicals formed either very near to or even in the space between the two metals of the DNC. It is unclear whether this peroxide species is in the O2(2-), O2(•)(-), HO2(-), or the H2O2 form and what is the detailed electronic structure and binding geometry including the DNC. In order to answer what form of this dioxygen species was observed in the DNC of the 1.8 Å X-ray CcO crystal structure (3S8G), we have applied broken-symmetry density functional theory (BS-DFT) geometric and energetic calculations (using OLYP potential) on large DNC cluster models with different Fea3-CuB oxidation and spin states and with O2(2-), O2(•)(-), HO2(-), or H2O2 in the bridging position. By comparing the DFT optimized geometries with the X-ray crystal structure (3S8G), we propose that the bridging peroxide is HO2(-). The X-ray crystal structure is likely to represent the superposition of the Fea3(2+)-(HO2(-))-CuB(+) DNC's in different states (Fe(2+) in low spin (LS), intermediate spin (IS), or high spin (HS)) with the majority species having the proton of the HO2(-) residing on the oxygen atom (O1) which is closer to the Fea3(2+) site in the Fea3(2+)-(HO-O)(-)-CuB(+) conformation. Our calculations show that the side chain of Tyr237 is likely trapped in the deprotonated Tyr237(-) anion form in the 3S8G X-ray crystal structure.

Computational prediction and in vitro analysis of potential physiological ligands of the bile acid binding site in cytochrome c oxidase

A conserved bile acid site has been crystallographically defined in the membrane domain of mammalian and Rhodobacter sphaeroides cytochrome c oxidase (RsCcO). Diverse amphipathic ligands were shown previously to bind to this site and affect the electron transfer equilibrium between heme a and a3 cofactors by blocking the K proton uptake path. Current studies identify physiologically relevant ligands for the bile acid site using a novel three-pronged computational approach: ROCS comparison of ligand shape and electrostatics, SimSite3D comparison of ligand binding site features, and SLIDE screening of potential ligands by docking. Identified candidate ligands include steroids, nicotinamides, flavins, nucleotides, retinoic acid, and thyroid hormones, which are predicted to make key protein contacts with the residues involved in bile acid binding. In vitro oxygen consumption and ligand competition assays on RsCcO wildtype and its Glu101Ala mutant support regulatory activity and specificity of some of these ligands. An ATP analog and GDP inhibit RsCcO under low substrate conditions, while fusidic acid, cholesteryl hemisuccinate, retinoic acid, and T3 thyroid hormone are more potent inhibitors under both high and low substrate conditions. The sigmoidal kinetics of RsCcO inhibition in the presence of certain nucleotides is reminiscent of previously reported ATP inhibition of mammalian CcO, suggesting regulation involving the conserved core subunits of both mammalian and bacterial oxidases. Ligand binding to the bile acid site is noncompetitive with respect to cytochrome c and appears to arrest CcO in a semioxidized state with some resemblance to the "resting" state of the enzyme.

Molecular dynamics simulations reveal highly permeable oxygen exit channels shared with water uptake channels in photosystem II

Photosystem II (PSII) catalyzes the oxidation of water in the conversion of light energy into chemical energy in photosynthesis. Water delivery and oxygen removal from the oxygen evolving complex (OEC), buried deep within PSII, are critical requirements to facilitate the reaction and minimize reactive oxygen damage. It has often been assumed that water and oxygen travel through separate channels within PSII, as demonstrated in cytochrome c oxidase. This study describes all-atom molecular dynamics simulations of PSII designed to investigate channels by fully characterizing the distribution and permeation of both water and oxygen. Interestingly, most channels found in PSII were permeable to both oxygen and water, however individual channels exhibited different energetic barriers for the two solutes. Several routes for oxygen diffusion within PSII with low energy permeation barriers were found, ensuring its fast removal from the OEC. In contrast, all routes for water showed significant energy barriers, corresponding to a much slower permeation rate for water through PSII. Two major factors were responsible for this selectivity: (1) hydrogen bonds between water and channel amino acids, and (2) steric restraints. Our results reveal the presence of a shared network of channels in PSII optimized to both facilitate the quick removal of oxygen and effectively restrict the water supply to the OEC to help stabilize and protect it from small water soluble inhibitors.

Ligand access to the active site in Thermus thermophilus ba(3) and bovine heart aa(3) cytochrome oxidases

Knowledge of the structure and dynamics of the ligand channel(s) in heme-copper oxidases is critical for understanding how the protein environment modulates the functions of these enzymes. Using photolabile NO and O(2) carriers, we recently found that NO and O(2) binding in Thermus thermophilus (Tt) ba(3) is ~10 times faster than in the bovine enzyme, indicating that inherent structural differences affect ligand access in these enzymes. Using X-ray crystallography, time-resolved optical absorption measurements, and theoretical calculations, we investigated ligand access in wild-type Tt ba(3) and the mutants, Y133W, T231F, and Y133W/T231F, in which tyrosine and threonine in the O(2) channel of Tt ba(3) are replaced by the corresponding bulkier tryptophan and phenylalanine, respectively, present in the aa(3) enzymes. NO binding in Y133W and Y133W/T231F was found to be 5 times slower than in wild-type ba(3) and the T231F mutant. The results show that the Tt ba(3) Y133W mutation and the bovine W126 residue physically impede NO access to the binuclear center. In the bovine enzyme, there is a hydrophobic "way station", which may further slow ligand access to the active site. Classical simulations of diffusion of Xe to the active sites in ba(3) and bovine aa(3) show conformational freedom of the bovine F238 and the F231 side chain of the Tt ba(3) Y133W/T231F mutant, with both residues rotating out of the ligand channel, resulting in no effect on ligand access in either enzyme.