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

Programming xenon diffusion in maltose-binding protein

Protein interiors contain void space that can bind small gas molecules. Determination of gas pathways and kinetics in proteins has been an intriguing and challenging task. Here, we combined computational methods and the hyperpolarized xenon-129 chemical exchange saturation transfer (hyper-CEST) NMR technique to investigate xenon (Xe) exchange kinetics in maltose-binding protein (MBP). A salt bridge ∼9 Å from the Xe-binding site formed upon maltose binding and slowed the Xe exchange rate, leading to a hyper-CEST 129Xe signal from maltose-bound MBP. Xe dissociation occurred faster than dissociation of the salt bridge, as shown by 13C NMR spectroscopy and variable-B1 hyper-CEST experiments. "Xe flooding" molecular dynamics simulations identified a surface hydrophobic site, V23, that has good Xe binding affinity. Mutations at this site confirmed its role as a secondary exchange pathway in modulating Xe diffusion. This shows the possibility for site-specifically controlling xenon protein-solvent exchange. Analysis of the available MBP structures suggests a biological role of MBP's large hydrophobic cavity to accommodate structural changes associated with ligand binding and protein-protein interactions.

Small-Molecule Tunnels in Metalloenzymes Viewed as Extensions of the Active Site

Rigorous substrate selectivity is a hallmark of enzyme catalysis. This selectivity is generally ascribed to a thermodynamically favorable process of substrate binding to the enzyme active site based upon complementary physiochemical characteristics, which allows both acquisition and orientation. However, this chemical selectivity is more difficult to rationalize for diminutive molecules that possess too narrow a range of physical characteristics to allow either precise positioning or discrimination between a substrate and an inhibitor. Foremost among these small molecules are dissolved gases such as H₂, N₂, O₂, CO, CO₂, NO, N₂O, NH₃, and CH₄ so often encountered in metalloenzyme catalysis. Nevertheless, metalloenzymes have evolved to metabolize these small-molecule substrates with high selectivity and efficiency.The soluble methane monooxygenase enzyme (sMMO) acts upon two of these small molecules, O₂ and CH₄, to generate methanol as part of the C1 metabolic pathway of methanotrophic organisms. sMMO is capable of oxidizing many alternative hydrocarbon substrates. Remarkably, however, it will preferentially oxidize methane, the substrate with the fewest discriminating physical characteristics and the strongest C-H bond. Early studies led us to broadly attribute this specificity to the formation of a "molecular sieve" in which a methane- and oxygen-sized tunnel provides a size-selective route from bulk solvent to the completely buried sMMO active site. Indeed, recent cryogenic and serial femtosecond ambient temperature crystallographic studies have revealed such a route in sMMO. A detailed study of the sMMO tunnel considered here in the context of small-molecule tunnels identified in other metalloenzymes reveals three discrete characteristics that contribute to substrate selectivity and positioning beyond that which can be provided by the active site itself. Moreover, the dynamic nature of many tunnels allows an exquisite coordination of substrate binding and reaction phases of the catalytic cycle. Here we differentiate between the highly selective molecular tunnel, which allows only the one-dimensional transit of small molecules, and the larger, less-selective channels found in typical enzymes. Methods are described to identify and characterize tunnels as well as to differentiate them from channels. In metalloenzymes which metabolize dissolved gases, we posit that the contribution of tunnels is so great that they should be considered to be extensions of the active site itself. A full understanding of catalysis by these enzymes requires an appreciation of the roles played by tunnels. Such an understanding will also facilitate the use of the enzymes or their synthetic mimics in industrial or pharmaceutical applications.

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.

Menthol Binding to the Human α4β2 Nicotinic Acetylcholine Receptor Facilitated by Its Strong Partitioning in the Membrane

We utilize various computational methodologies to study menthol's interaction with multiple organic phases, a lipid bilayer, and the human α4β2 nicotinic acetylcholine receptor (nAChR), the most abundant nAChR in the brain. First, force field parameters developed for menthol are validated in alchemical free energy perturbation simulations to calculate solvation free energies of menthol in water, dodecane, and octanol and compare the results against experimental data. Next, umbrella sampling is used to construct the free energy profile of menthol permeation across a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The results from a flooding simulation designed to study the water-membrane partitioning of menthol in a POPC lipid bilayer are used to determine the penetration depth and the preferred orientation of menthol in the bilayer. Finally, employing both docking and flooding simulations, menthol is shown to bind to different sites on the human α4β2 nAChR. The most likely binding mode of menthol to a desensitized membrane-embedded α4β2 nAChR is identified to be via a membrane-mediated pathway in which menthol binds to the sites at the lipid-protein interface after partitioning in the membrane. A rare but distinct binding mode in which menthol binds to the extracellular opening of receptor's ion permeation pore is also reported.

Dimeric structures of quinol-dependent nitric oxide reductases (qNORs) revealed by cryo-electron microscopy

Quinol-dependent nitric oxide reductases (qNORs) are membrane-integrated, iron-containing enzymes of the denitrification pathway, which catalyze the reduction of nitric oxide (NO) to the major ozone destroying gas nitrous oxide (N₂O). Cryo-electron microscopy structures of active qNOR from Alcaligenes xylosoxidans and an activity-enhancing mutant have been determined to be at local resolutions of 3.7 and 3.2 Å, respectively. They unexpectedly reveal a dimeric conformation (also confirmed for qNOR from Neisseria meningitidis) and define the active-site configuration, with a clear water channel from the cytoplasm. Structure-based mutagenesis has identified key residues involved in proton transport and substrate delivery to the active site of qNORs. The proton supply direction differs from cytochrome c-dependent NOR (cNOR), where water molecules from the cytoplasm serve as a proton source similar to those from cytochrome c oxidase.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

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.

Functional importance of Glutamate-445 and Glutamate-99 in proton-coupled electron transfer during oxygen reduction by cytochrome bd from Escherichia coli

The recent X-ray structure of the cytochrome bd respiratory oxygen reductase showed that two of the three heme components, heme d and heme b₅₉₅, have glutamic acid as an axial ligand. No other native heme proteins are known to have glutamic acid axial ligands. In this work, site-directed mutagenesis is used to probe the roles of these glutamic acids, E445 and E99 in the E. coli enzyme. It is concluded that neither glutamate is a strong ligand to the heme Fe and they are not the major determinates of heme binding to the protein. Although very important, neither glutamate is absolutely essential for catalytic function. The close interactions between the three hemes in cyt bd result in highly cooperative properties. For example, mutation of E445, which is near heme d, has its greatest effects on the properties of heme b₅₉₅ and heme b₅₅₈. It is concluded that 1) O₂ binds to the hydrophilic side of heme d and displaces E445; 2) E445 forms a salt bridge with R448 within the O₂ binding pocket, and both residues play a role to stabilize oxygenated states of heme d during catalysis; 3) E445 and E99 are each protonated accompanying electron transfer to heme d and heme b₅₉₅, respectively; 4) All protons used to generate water within the heme d active site come from the cytoplasm and are delivered through a channel that must include internal water molecules to assist proton transfer: [cytoplasm] → E107 → E99 (heme b₅₉₅) → E445 (heme d) → oxygenated heme d.

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.

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.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.