H-NOX domains display different tunnel systems for ligand migration

Molecular dynamics simulations suggest Thiosemicarbazones can bind p53 cancer mutant R175H

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC - p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi‑sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.

Gas Sensing by Bacterial H-NOX Proteins: An MD Study

Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O₂ diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O₂ > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.

Bacterial Haemoprotein Sensors of NO: H-NOX and NosP

Low concentrations of nitric oxide (NO) modulate varied behaviours in bacteria including biofilm dispersal and quorum sensing-dependent light production. H-NOX (haem-nitric oxide/oxygen binding) is a haem-bound protein domain that has been shown to be involved in mediating these bacterial responses to NO in several organisms. However, many bacteria that respond to nanomolar concentrations of NO do not contain an annotated H-NOX domain. Nitric oxide sensing protein (NosP), a newly discovered bacterial NO-sensing haemoprotein, may fill this role. The focus of this review is to discuss structure, ligand binding, and activation of H-NOX proteins, as well as to discuss the early evidence for NO sensing and regulation by NosP domains. Further, these findings are connected to the regulation of bacterial biofilm phenotypes and symbiotic relationships.

Regulation of Neuronal Oxygen Responses in C. elegans Is Mediated through Interactions between Globin 5 and the H-NOX Domains of Soluble Guanylate Cyclases

Soluble guanylate cyclases (sGCs) are gas-binding proteins that control diverse physiological processes such as vasodilation, platelet aggregation, and synaptic plasticity. In the nematode Caenorhabditis elegans, a complex of sGCs, GCY-35 and GCY-36, functions in oxygen (O2) sensing. Previous studies suggested that the neuroglobin GLB-5 genetically interacts with GCY-35, and that the inhibitory effect of GLB-5 on GCY-35 function is necessary for fast recovery from prolonged hypoxia. In this study, we identified mutations in gcy-35 and gcy-36 that impact fast recovery and other phenotypes associated with GLB-5, without undermining sGC activity. These mutations, heb1 and heb3, change conserved amino acid residues in the regulatory H-NOX domains of GCY-35 and GCY-36, respectively, and appear to suppress GLB-5 activity by different mechanisms. Moreover, we observed that short exposure to 35% O2 desensitized the neurons responsible for ambient O2 sensing and that this phenomenon does not occur in heb1 animals. These observations may implicate sGCs in neuronal desensitization mechanisms far beyond the specific case of O2 sensing in nematodes. The conservation of functionally important regions of sGCs is supported by examining site-directed mutants of GCY-35, which suggested that similar regions in the H-NOX domains of O2 and NO-sensing sGCs are important for heme/gas interactions. Overall, our studies provide novel insights into sGC activity and regulation, and implicate similar structural determinants in the control of both O2 and NO sensors. Significance statement: Soluble guanylate cyclases (sGCs) control essential and diverse physiological processes, including memory processing. We used Caenorhabditis elegans to explore how a neuroglobin inhibits a complex of oxygen-sensing sGCs, identifying sGC mutants that resist inhibition. Resistance appears to arise by two different mechanisms: increased basal sGC activity or disruption of an interaction with neuroglobin. Our findings demonstrate that the inhibition of sGCs by neuroglobin is essential for rapid adaptation to either low or high oxygen levels, and that similar structural regions are key for regulating both oxygen and nitric oxide sensors. Based on our structural and functional analyses, we present the hypothesis that neuroglobin-sGC interactions may be generally important for adaptation processes, including those in organisms with more complex neurological functions.

Entry of water into the distal heme pocket of soluble guanylate cyclase β1 H-NOX domain alters the ligated CO structure: a resonance Raman and in silico simulation study

Water accessing into the heme pocket and alters the structures of CO–sGC (heme), exhibiting two different v_Fe–CO stretching modes.

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.

Theoretical investigation on the diatomic ligand migration process and ligand binding properties in non-O2-binding H-NOX domain

The Nostoc sp (Ns) H-NOX (heme-nitric oxide or OXygen-binding) domain shares 35% sequence identity with soluble guanylate cyclase (sGC) and exhibits similar ligand binding property with the sGC. Previously, our molecular dynamic (MD) simulation work identified that there exists a Y-shaped tunnel system hosted in the Ns H-NOX interior, which servers for ligand migration. The tunnels were then confirmed by Winter et al. [PNAS 2011;108(43):E 881-889] recently using x-ray crystallography with xenon pressured conditions. In this work, to further investigate how the protein matrix of Ns H-NOX modulates the ligand migration process and how the distal residue composition affects the ligand binding prosperities, the free energy profiles for nitric oxide (NO), carbon monooxide (CO), and O2 migration are explored using the steered MDs simulation and the ligand binding energies are calculated using QM/MM schemes. The potential of mean force profiles suggest that the longer branch of the tunnel would be the most favorable route for NO migration and a second NO trapping site other than the distal heme pocket along this route in the Ns H-NOX was identified. On the contrary, CO and O2 would prefer to diffuse via the shorter branch of the tunnel. The QM/MM (quantum mechanics/molecular mechanics) calculations suggest that the hydrophobic distal pocket of Ns H-NOX would provide an approximately vacuum environment and the ligand discrimination would be determined by the intrinsic binding properties of the diatomic gas ligand to the heme group.

Dynamic change of heme environment in soluble guanylate cyclase and complexation of NO-independent drug agents with H-NOX domain

Soluble guanylate cyclase is a heterodimer receptor that functions in several signal transduction pathways. Conversion of guanosine 5'-triphosphate to 3',5'-cyclic monophosphate second messenger at the catalytic domain is regulated by the changes at heme nitric oxide/oxygen domain of the β-subunit. To better understand conformational changes at heme site that may impact on activities of catalytic domain, three soluble guanylate cyclase homolog proteins with heme at Fe-His state were investigated, and their dynamic behaviors were monitored in both unliganded (apo) and complex with heme. As a result of dynamic conformational changes, Lys110, Asp45, Arg135, and Glu41 were found interacting with the site gate, which may interfere with transportation of small molecules in and out of the heme site. An alternative binding site adjacent to that of heme was identified. Binding affinity of several nitric oxide-independent activators and heme-dependent stimulators was examined, and their binding modes in the heme site and in the alternative binding site in the human soluble guanylate cyclase enzyme were computationally simulated. The calculated binding energies were used as criteria to filter results of virtual high-throughput screenings based on FlexX ligand-docking algorithm and absorption, distribution, metabolism, excretion, and toxicity properties on databases of available drugs. The identified drugs from virtual high-throughput screening have been suggested for experimental investigations, based on which they may either be directly repurposed or require structural modifications for better physico-chemical and pharmacological properties.

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

We use a form of "freeze-trap, kinetic crystallography" to explore the migration of Xe atoms away from the dinuclear heme a(3)/Cu(B) center in Thermus thermophilus cytochrome ba(3) oxidase. This enzyme is a member of the heme-copper oxidase superfamily and is thus crucial for dioxygen-dependent life. The mechanisms involved in the migration of oxygen, water, electrons, and protons into and/or out of the specialized channels of the heme-copper oxidases are generally not well understood. Pressurization of crystals with Xe gas previously revealed a O(2) diffusion channel in cytochrome ba(3) oxidase that is continuous, Y-shaped, 18-20 Å in length and comprised of hydrophobic residues, connecting the protein surface within the bilayer to the a(3)-Cu(B) center in the active site. To understand movement of gas molecules within the O(2) channel, we performed crystallographic analysis of 19 Xe laden crystals freeze-trapped in liquid nitrogen at selected times between 0 and 480 s while undergoing outgassing at room temperature. Variation in Xe crystallographic occupancy at five discrete sites as a function of time leads to a kinetic model revealing relative degrees of mobility of Xe atoms within the channel. Xe egress occurs primarily through the channel formed by the Xe1 → Xe5 → Xe3 → Xe4 sites, suggesting that ingress of O(2) is likely to occur by the reverse of this process. The channel itself appears not to undergo significant structural changes during Xe migration, thereby indicating a passive role in this important physiological function.

Carbon Dioxide Migration Pathways in Proteins

Some of the most important biological processes, such as carbon fixation, are dependent on protein-gas interactions. The motion of CO2 through the enzyme phosphoenolpyruvate carboxykinase was investigated using extensive all-atom molecular dynamics simulations. Three discrete migration pathways were located, suggesting the protein directs the movement of CO2. The chemical nature of these pathways is discussed, as are their biotechnological ramifications.

Structure and regulation of soluble guanylate cyclase

Nitric oxide (NO) is an essential signaling molecule in biological systems. In mammals, the diatomic gas is critical to the cyclic guanosine monophosphate (cGMP) pathway as it functions as the primary activator of soluble guanylate cyclase (sGC). NO is synthesized from l-arginine and oxygen (O(2)) by the enzyme nitric oxide synthase (NOS). Once produced, NO rapidly diffuses across cell membranes and binds to the heme cofactor of sGC. sGC forms a stable complex with NO and carbon monoxide (CO), but not with O(2). The binding of NO to sGC leads to significant increases in cGMP levels. The second messenger then directly modulates phosphodiesterases (PDEs), ion-gated channels, or cGMP-dependent protein kinases to regulate physiological functions, including vasodilation, platelet aggregation, and neurotransmission. Many studies are focused on elucidating the molecular mechanism of sGC activation and deactivation with a goal of therapeutic intervention in diseases involving the NO/cGMP-signaling pathway. This review summarizes the current understanding of sGC structure and regulation as well as recent developments in NO signaling.

Tunnels modulate ligand flux in a heme nitric oxide/oxygen binding (H-NOX) domain

Interior topological features, such as pockets and channels, have evolved in proteins to regulate biological functions by facilitating the diffusion of biomolecules. Decades of research using the globins as model heme proteins have clearly highlighted the importance of gas pockets around the heme in controlling the capture and release of O(2). However, much less is known about how ligand migration contributes to the diverse functions of other heme protein scaffolds. Heme nitric oxide/oxygen binding (H-NOX) domains are a conserved family of gas-sensing heme proteins with a divergent fold that are critical to numerous signaling pathways. Utilizing X-ray crystallography with xenon, a tunnel network has been shown to serve as a molecular pathway for ligand diffusion. Structure-guided mutagenesis results show that the tunnels have unexpected effects on gas-sensing properties in H-NOX domains. The findings provide insights on how the flux of biomolecules through protein scaffolds modulates protein chemistry.