Next Article in Journal
The Role of NOX4 in Parkinson’s Disease with Dementia
Previous Article in Journal
Novel Targeted Anti-Tumor Nanoparticles Developed from Folic Acid-Modified 2-Deoxyglucose
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 20
Issue 3
10.3390/ijms20030698
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands

by
Rana Rehan Khalid
1,2,
Abdul Rauf Siddiqi
1,*,
Efstratios Mylonas
3,
Arooma Maryam
1 and
Michael Kokkinidis
2,3,*
1
Department of Biosciences, COMSATS University, Islamabad 45550, Pakistan
2
Department of Biology, University of Crete, 70013 Heraklion, Greece
3
Institute of Molecular Biology and Biotechnology, Foundation for Research and Technology-Hellas (IMBB-FORTH), 70013 Heraklion, Greece
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(3), 698; https://doi.org/10.3390/ijms20030698
Submission received: 21 December 2018 / Revised: 29 January 2019 / Accepted: 3 February 2019 / Published: 6 February 2019
(This article belongs to the Section Macromolecules)
Browse Figures
Versions Notes

Abstract

:
Soluble guanylate cyclase (sGC) regulates numerous physiological processes. The β subunit Heme Nitric Oxide/Oxygen (HNOX) domain makes this protein sensitive to small gaseous ligands. The structural basis of the activation mechanism of sGC under the influence of ligands (NO, O2, CO) is poorly understood. We examine the effect of different ligands on the human sGC HNOX domain. HNOX systems with gaseous ligands were generated and explored using Molecular Dynamics (MD). The distance between heme Fe2+ and histidine in the NO-ligated HNOX (NO-HNOX) system is larger compared to the O2, CO systems. NO-HNOX rapidly adopts the conformation of the five-group metal coordination system. Loops α, β, γ and helix-f exhibit increased mobility and different hydrogen bond networks in NO-HNOX compared to the other systems. The removal of His from the Fe coordination sphere in NO-HNOX is assisted by interaction of the imidazole ring with the surrounding residues which in turn leads to the release of signaling helix-f and activation of the sGC enzyme. Insights into the conformational dynamics of a human sGC HNOX domain, especially for regions which are functionally critical for signal transduction, are valuable in the understanding of cardiovascular diseases.

Graphical Abstract

1. Introduction

There are two types of enzymes with guanylyl cyclase activity; membrane-bound particulate guanylate cyclase (pGC) and cytosolic soluble guanylate cyclase (sGC) [1]. These enzymes are crucial for cGMP production, which regulates numerous physiological and pathophysiological processes [2], such as cardiovascular diseases, platelet aggregation, neurodegeneration, erectile dysfunction, etc. Knowledge of the structure and dynamics of the sGC enzyme is indispensable to understanding the mechanism of function [3]. Many studies have shown that the Heme Nitric Oxide/Oxygen (HNOX) domain of sGC is the starting point of signal transduction. The heme prosthetic group allows this domain to sense the small gaseous ligands CO, O2, and NO [4]. It is the binding of NO to the HNOX domain that stimulates the catalytic domain to cyclize GTP into cGMP, responsible for regulating many downstream pathways. sGC is a 150 kDa macromolecule consisting of subunits α and β [5,6]. The complete structure of sGC has not yet been elucidated, and so therefore, the interdomain interactions of sGC remain poorly understood. Many aspects of the structure and dynamics of the protein and its individual domains still require a detailed and comprehensive analysis.
The crystal structure of the bacterial HNOX domain has been resolved, but essential heme binding pocket residues (Y135, S137, and R139, Y-S-R motif) are also conserved in the human HNOX [7]. In this domain diatomic gaseous molecules such as NO, CO, and O2 bind with the heme group on the so-called distal site, while on the proximal side a coordinate bond forms between heme and histidine 105 (H105) (6c) [8]. The heme group contains propionic acid moieties involved in hydrogen bond interactions with the conserved distal pocket residues of the Y-S-R motif [9]. The reduced heme of sGC favors the binding of diatomic gas molecules, notably CO and NO [10]. Heme is found inside the HNOX domain in two coordination states, i.e., a six coordination (6c) and a five coordination (5c) state. During activation, the gaseous ligand NO binds to heme and changes its state from 6c to 5c. Breaking of the coordinate bond between H105 and iron results in the activation of sGC [11,12,13].
Dealing with metalloproteins like sGC poses significant challenges. Iron (Fe) is a transition metal and can exist in ferrous +2 and ferric +3 states. The Fe2+ coordination state of the HNOX heme favors the binding of small gaseous ligands, such as NO, O2, and CO [14]. The affinity of NO and CO ligands for the heme prosthetic group of sGC is considerably higher than the affinity of O2. Discrimination between these diatomic molecules is controlled by heme distal site residues, whereas residues surrounding the proximal histidine perturb the electron density of the imidazole ring, further contributing to the transfer of charge from histidine to Fe2+. When NO binds to 5c heme, the d orbital electrons of the iron move to an unstable high spin state, with the formation of a transient 6c complex, that converts to 5c with the loss of the histidine-iron interaction [15,16].
In this paper, we elucidate the atomistic dynamics of NO/O2/CO-ligated human HNOX domain in comparison with control apo-HNOX. A comparative modeling approach was used to estimate the structure of the human HNOX domain. The models were then further subjected to detailed dynamics investigation. Stability and hydrogen bond analysis revealed the binding pocket structural transitions in response to different ligands (NO, O2, CO). Furthermore, essential dynamics (ED) were employed to identify noticeable structural variations. The following five complexes were examined; 1. Apo-HNOX (without gaseous ligands), 2. NO-HNOX (5c) 3. NO-HNOX (6c) with the bond between H105 and Fe removed to study the effects of bond-breaking (5c), 4. O2-HNOX (6c) 5. CO-HNOX (6c). Each system is examined for 50 ns.

2. Results and Discussion

2.1. Homology Modeling

The HNOX domain harbors a heme prosthetic group which can sense diatomic gaseous molecules such as CO, NO, and O2. The human HNOX N-terminal domain (1–186 a.a) has 34% sequence identity with the nostoc bacterial template. Modeling results showed that 94.6% of the residues are found in favored regions of the Ramachandran plot and 97.8% in allowed regions. After energy minimization through YASARA (Figure 1) 98.4% of residues occupied favored regions and 99.5% allowed regions. The global structure quality score [17] of the predicted model is statistically significant (TM = 0.83), which makes the model reliable for further structural interpretation.

2.2. QM Calculations

The Metal Center Parameter Builder (MCPB.py) was used in conjunction with Gaussian09 to parametrize the active site (His-heme-NO/O2/CO) using models of differing sizes to balance speed and accuracy [18] (Figure 2). Smaller models were used to estimate the iron associated bond and angle parameters, while the larger models were used for partial charge calculations. The accuracy of the optimized small model metal-associated bonds and angles was evaluated with Chimera [19]. The coordination states of aforesaid systems remain stable during optimization and force constant calculations (Figure 2c,f,i). The same level of theory and the same basis set (B3LYP/6-31G*) was applied to all complexes to ensure local minima were found. The force constants of the active site atoms were estimated with the Seminario method [20] and all bonded and non-bonded type parameters were within acceptable range. The histidine imidazole ring adopts a stagger conformation forming a 45° angle with the plane of heme opposite to pyrrolic nitrogen atoms. This conformation is favorable to the binding of gaseous ligand such as O2. Furthermore, the distance between the histidine imidazole ring and iron is high, in the NO-ligated system, in contrast to the CO and O2 ones. Moreover, O2- and NO-bonded complexes have bent conformations; however, the CO-ligated complex adopts a linear configuration with a 179° angle (Figure 2c,f,i).

2.3. MD Stability Analysis

CPPTRAJ [21] was used to evaluate the conformational stability of the HNOX systems over the course of the simulation. The root-mean-square deviation (RMSD) of the five systems; apo-HNOX, 6c NO/CO/O2− HNOX and 5c NO-HNOX, was analyzed for all backbone atoms, and the RMSD results indicate that the largest variations were observed in the NO/CO-HNOX systems in contrast to the O2 and apo system that did not experience noteworthy RMSD variations (Figure 3). In the 6c NO-HNOX system the Fe-H105 bond length increased abruptly at the time interval of 4–6 ns. This is possibly the point at which the activation of the HNOX domain occurs with the loss of the Fe-H105 interaction. Both NO-bonded HNOX complexes exhibit an average RMSD of 2.7 Å in the time interval of 6–25 ns while from 25 ns to 50 ns the average RMSD of 5c NO-HNOX was further increased to 4.4 Å which can be attributed to the release of signaling helix-f. Of the other systems, the CO-ligated one underwent the highest RMSD fluctuations with an average RMSD of 3.4 Å. The loops α, β and the flanking helix-f were the regions which contributed the most to the RMSD increase.
Root mean square fluctuation (RMSF) stability analysis helps to investigate the mobility of residues under the influence of different ligands over the course of the entire simulation. The peaks represent areas of increased residual flexibility (Figure 4). The results indicate that O2 binding to the HNOX domain has minor effects on the residual fluctuations, similar to apo-HNOX. In NO-bonded complexes regions 105–118 (helix-f, loop α) and 124-129 (loop β) are significantly more mobile with a maximum RMSF of 4.2 Å for loop α. The 5c NO-HNOX system exhibited a similar behavior with a slightly smaller peak of 3.6 Å. In contrast, the highest peak of the CO-ligated system at 5.1 Å is observed at the loop β region. These observations suggest that the aforementioned sites possibly contribute to the release of heme from H105.
Another conformational stability indicator, the radius of gyration (Rg), provides a measure of the compactness of the protein by monitoring the effective size of protein. Rg is defined as the mass-weighted root mean square distance of a group of atoms from their common center of mass (Figure 5). The Rg-analysis results are consistent with previous stability findings. The O2-ligated 6c HNOX radius of gyration remains stable at ~16.9 Å. The 5c NO-HNOX system has the highest average Rg at 17.4 Å. The broadest range of configurational space variation was observed for the 6c NO-HNOX system with an Rg below 17 Å at the beginning of the simulation, stabilizing at ~17.3 Å after 6ns. Finally, the average Rg of apo and CO-ligated systems was ~17 Å.

2.4. Hydrogen Bond Occurrence

Hydrogen bond analysis was performed to elucidate the binding modes of 6c O2/NO/CO-HNOX and 5c NO-HNOX systems (Figure 6). Although the number of O2-HNOX binding pocket residues that take part in electrostatic interactions with heme is smaller, they show a more robust hydrogen bond network compared to the NO and CO complexes. The binding modes pattern of apo-HNOX (not shown) is very similar to the O2-HNOX, except for the Y2 and R116 electrostatic interactions. The hydrogen bond fractional occurrence is high because the O2-HNOX (and apo-HNOX) system does not undergo as drastic conformational transitions as the NO/CO-bonded systems. Residues Y135, R137, R139, R116, and H105 are the major contributors of hydrogen bonds with the heme propionic acid moiety in all four systems. Weaker hydrogen bond interactions were also found for M1, Y49, A117, E138, and L142 in the different systems (Table 1). H105 hydrogen donor residues such as A109, L108, and acceptor residues such as 101L, 102A, and 118P may be involved in the perturbation of the imidazole ring electron density leading to the formation of the 5c heme-NO complex (Table 2). Earlier studies suggest that the interaction of H105 with the surrounding residues drive the imidazole ring to a high spin state, while the slightly negative charge of NO favors the low spin state, resulting in stretching of the bond between iron (Fe2+) and H105 [22,23,24]. Furthermore, interaction of R116 with the heme carboxylic groups (O1A, O2A, O1D, O2D) is stronger in 5c NO-HNOX and 6c CO-HNOX complexes, compared to the apo, 6c O2-HNOX and 6c NO-HNOX complexes (Table 1, Figure 6). The liberation of the heme cofactor from H105 might be contributing to increasing R116 hydrogen bonding with the heme carboxylic groups, compared to 6c NO-HNOX. This interaction possibly enhances the structural fluctuations of loop α, further contributing to sGC activation. Previous studies have shown that residue R116 resides on the periphery of HNOX, facing towards the adjacent PAS domain [25], implicating it in cyclase activity. Moreover, the apo-HNOX and O2-HNOX complexes heme propionic moiety bears considerable hydrogen bond interactions with Y2, Y112, and Y135 enhancing the overall compactness of the system (Table 1, Figure 6). It is also discerned that heme carboxyl group in CO-HNOX showed more robust interaction with R116 amongst all four systems, although there is no H105-Fe bond cleavage observed in CO-HNOX system (Table 1).
To visualize the entirety of the electrostatic interactions of heme with the diatomic ligands (NO, CO, O2), all systems were treated with the Ligplot tool [26]. The H105-Fe bond undergoes a bond stretching of 3.23 Å in the NO-ligated system (Figure 7) contributing to the release of signaling helix-f while establishing a robust interaction with R116. In contrast, the bond length remained less than 2.5 Å in the CO− and O2-ligated systems.

2.5. Essential Dynamics

Essential dynamics was performed to identify the overall normalized pattern of motion in all four systems. The principal component analysis facilitates further elucidation of the sGC activation mechanism. The diagonalization of the covariance matrix contributes to getting more meaningful configurational space, revealing functionally critical structural transitions [27]. The backbone atoms of apo, CO/NO/O2-HNOX complexes were considered to perform essential dynamics. Interestingly, the first two principal components correspond to 85% of the overall positional transitions in all systems. The motion of the first two dominant modes was examined through the VMD normal mode analysis plugin and the porcupine plot was used to show the magnitude and direction of selected eigenvectors [28].
Critical regions of the human HNOX domain that contribute to high variance during essential dynamic analysis are shown in Figure 8. The regions include helix-f and loop α and β in the NO-ligated HNOX complexes, although in the case of the CO-ligated complex the helix-f displacement was much less pronounced. Interestingly, in the control simulation (apo-HNOX), helix-f and loop α and β regions experienced trivial structural transitions, in contrast with other systems. These results are compatible with the stability analysis. Moreover, the length of the Fe-H105 bond strongly fluctuates in the 6c NO-HNOX system. Experimental data [11] has shown that the Fe-H105 bond breaking event upon NO binding to heme, during the sGC activation, contributes to the liberation of signaling helix-f. The complex after the bond breaking event is described by the 5c NO-HNOX system. The 5c NO-HNOX system undergoes significant structural transitions almost similar to 6c NO-HNOX. In contrast, the O2-ligated HNOX did not experience any significant change and the system maintained its compactness throughout the simulation (Figure 9). The structural transitions of helix-f and loops α and β are in agreement with earlier experimental and theoretical studies which suggest they play a crucial role in the activation of the soluble guanylate cyclase enzyme, after binding of the nitric oxide axial ligand [29,30]. The strong interaction of residues S137, R139 in loop γ and R116 in loop α with the heme moiety possibly contribute to the overall dynamic changes associated with sGC activation.

3. Materials and Methods

3.1. Comparative Modeling of Human sGC HNOX Domain

To predict the structure of the HNOX domain, Protein BLAST [31] was used to find suitable protein structure templates. The 186 N-terminal residues that encode the HNOX domain of the sGC β subunit (NCBI Reference Sequence: NP_000848.1) [32], were subjected to comparative modeling. The structure of the nostoc bacterial sGC HNOX domain (PDB ID: 2O09) [6] was used as a template for the human sGC HNOX domain, with which it shares 34% identity. Homology modeling was done with Modeller9.16 [33], and the best model was selected based on the z-DOPE, estimated RMSD, estimated overlap scores and the quality of the Ramachandran plot [34]. Steric clashes were eliminated through energy minimization [35], and the structure was validated with MolProbity [36]. Furthermore, the builder module of pymol [37] was used to ligate NO, O2 and CO molecules to the heme group of the HNOX domain.

3.2. QM/MM Calculations

MCPB.py, a module of the Amber16 program package [38], which deals explicitly with metalloproteins [18], was employed to parametrize the H105-heme-NO/O2/CO active sites, acting as a bridge between Quantum mechanics (QM) and MD simulations. MCPB.py generates two models of which the smaller one is used for bond and angle parameters calculation and the larger one for charge calculation. The hybrid functional method B3LYP [39] with the 6-31G* basis set were applied for the optimization of the smaller model (heme-Fe-NO/O2/CO) while the force constant calculation was performed with the Seminario method [20]. This theoretical method uses a Cartesian Hessian matrix to calculate harmonic bond and angles. Atomic charges were estimated from the large model with the Merz-Kollman Restrained Electro Static Potential (RESP) method [40]. All calculations were performed using Gaussian09 [41], with the help of MCPB.py.

3.3. MD Simulations

Five different complexes, apo-HNOX, 6c O2/NO/CO-HNOX and 5c NO-HNOX (where the Fe-His bond is broken) were explored with MD simulations. The calculated active site (H105-heme-NO/O2/CO) parameters were combined with the ff14SB force field to fully parametrize the protein systems [42]. The Amber16 tleap module was used to neutralize the system and add the missing hydrogen atoms. The protonation state of H105 was set manually depending on the presence of a bond between the imidazole ring and iron. The systems were immersed in orthogonal boxes with TIP3P-model water molecules [43]. All systems were simulated enabling periodic boundary conditions, and the Particle Mesh Ewald (PME) [44] method was employed to calculate long-range electrostatic interactions with a cutoff of 10 Å. Initially, the systems were minimized for 10,000 steps. During the subsequent annealing, the Langevin thermostat was operated to control the temperature from 0 to 300 K at constant volume (NVT) for 200 ps. The SHAKE algorithm was employed to constrain all bonds involving hydrogen atoms [45]. All systems were investigated for 30 ns at 300 K and 1 atm constant pressure (NPT). Trajectory snapshots were obtained every 2 fs. Basic root mean square deviation (RMSD), root mean square fluctuation (RMSF) and radius of gyration (Rg) analysis was done to evaluate the systems’ stability throughout the simulation. Furthermore, hydrogen bond analysis and essential dynamics [46] were carried out with the CPPTRAJ module [20]. Visualization of the trajectories was achieved with VMD [47], while the plots were drawn with the help of Xmgrace [48] and Gnuplot [49].

3.4. Hydrogen Bond Occupancy

The analysis of hydrogen bonds is indispensable to studying the binding modes of complexes. For our systems, the hydrogen bonds of heme with the distal binding pocket residues of the HNOX domain, Y2, Y135, S137, E138, L142, R139, and the proximal binding pocket residues H105, R116, P118 were analyzed with CPPTRAJ. Additionally, solute-solute, solute-solvent, and solute-solvent-solute bridges hydrogen bond analysis of the entire HNOX systems was also carried out. The lifetime of hydrogen bond occupancies is indicative of the overall binding mode occupancies throughout the simulation and was visualized with Gnuplot.

3.5. Principal Component Analysis (PCA)

Essential dynamics were carried out to evaluate and interpret the multivariate data after detecting all correlated variable clusters. Furthermore, these correlated variables were converted to a small number of more meaningful uncorrelated variables [46]. The CPPTRAJ module of Amber16 was used to strip the solvent and ions from the trajectories before PCA analysis. The trajectories snapshots were then aligned over the average minimized structures. PCA was applied to all backbone atoms of the HNOX-O2/NO/CO and apo-HNOX systems. Diagonalization of covariance matrices assists in calculating the principal components of motion [50]. The trajectories were projected over eigenvectors corresponding to the first three largest eigenvalues of the correlation matrix. The most prominent fluctuations of different modes were investigated with the NMWiz VMD plugin [51].

4. Conclusions

Soluble guanylate cyclase is a regulatory protein involved in a multitude of physiological processes in humans such as platelet aggregation, vasodilation, and neurotransmission. Deciphering the ways in which it is regulated and regulates other proteins requires detailed structural characterization. The HNOX domain is the starting point in the signal transduction mechanism of this heterodimeric, multidomain protein. In the current study, we analyzed the effect of the physiologically relevant diatomic gaseous ligands CO/NO/O2 on the dynamics of the human HNOX domain, by employing a wide range of computational techniques. Residues in helix-f and loops α and β exhibited the most variable behavior, upon binding of the different gaseous ligands, highlighting the functional importance of these regions. These regions are located on the periphery of the inter-domain interface (HNOX-PAS) suggesting that they are involved in sGC activation. Among all complexes, 5c NO-HNOX experienced the highest positional displacement in helix-f and loop α, potentially facilitating the communication with adjacent sGC domains. The 5c NO-HNOX and 6c CO-HNOX residue R116 exhibits robust hydrogen bond interactions with the heme moiety, which could be critical to the communication with the neighboring domain. Interestingly, these interactions were either weak or absent in 6c NO-HNOX and 6c O2-HNOX. This work highlights important residues for sGC activation under physiological conditions that were not identified in our previous work [52] on the binding of sGC drugs. It is thus possible that a new class of compounds targeting the Fe-heme interface can be explored as potential sGC modulators. Exploring the HNOX binding modes is an important step in understanding the sGC activation mechanism and ultimately revealing regions with potential as drug targets for various sGC-associated diseases.

Author Contributions

M.K., and A.R.S. conceived and designed the experiments; R.R.K. performed the experiments and wrote the paper; R.R.K., E.M., and A.M. analyzed the data; A.R.S., E.M. and M.K. revised and improved the paper.

Funding

This research was partially funded by Erasmus Mundus Action 2, Strand I 552850-EM-1-2014-1-IT-ERA MUNDUS-EMA21 and Higher Education Commission (HEC) of Pakistan, grant number NRPU 4050. This work was co-financed by Greece and the European Union in the framework of the project “NATIONAL RESEARCH INFRASTRUCTURES ON INTEGRATED STRUCTURAL BIOLOGY, DRUG SCREENING EFFORTS AND DRUG TARGET FUNCTIONAL CHARACTERIZATION—NSPIRED” from the EPAnEK 2014-2020 action, operational programme “Competitiveness-Entrepreneurship-Innovation”.

Acknowledgments

Support by the High Performance Computing facility of the University of Crete DataCenter is acknowledged. We thank Alexis Molfetas for carefully reading the manuscript and helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pyriochou, A.; Papapetropoulos, A. Soluble guanylyl cyclase: More secrets revealed. Cell Signal. 2005. [CrossRef] [PubMed]
  2. Beste, K.Y.; Spangler, C.M.; Burhenne, H.; Koch, K.W.; Shen, Y.; Tang, W.J.; Kaever, V.; Seifert, R. Nucleotidyl Cyclase Activity of Particulate Guanylyl Cyclase A: Comparison with Particulate Guanylyl Cyclases E and F, Soluble Guanylyl Cyclase and Bacterial Adenylyl Cyclases Cyaa and Edema Factor. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed]
  3. Campbell, M.G.; Underbakke, E.S.; Potter, C.S.; Carragher, B.; Marletta, M.A. Single-particle EM reveals the higher-order domain architecture of soluble guanylate cyclase. Proc. Natl. Acad. Sci. USA 2014, 111, 2960–2965. [Google Scholar] [CrossRef] [Green Version]
  4. Jain, R.; Chan, M.K. Mechanisms of ligand discrimination by heme proteins. J. Biol. Inorg. Chem. 2003, 8, 1–11. [Google Scholar] [CrossRef] [PubMed]
  5. Friebe, A.; Koesling, D. Regulation of nitric oxide-sensitive guanylyl cyclase. Circ. Res. 2003, 93, 96–105. [Google Scholar] [CrossRef] [PubMed]
  6. Ma, X.; Sayed, N.; Beuve, A.; Van Den Akker, F. NO and CO differentially activate soluble guanylyl cyclase via a heme pivot-bend mechanism. EMBO J. 2007, 26, 578–588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Pellicena, P.; Karow, D.S.; Boon, E.M.; Marletta, M.A.; Kuriyan, J. Crystal structure of an oxygen-binding heme domain related to soluble guanylate cyclases. Proc. Natl. Acad. Sci. USA 2004, 101, 12854–12859. [Google Scholar] [CrossRef] [Green Version]
  8. Bikiel, D.E.; Boechi, L.; Capece, L.; Crespo, A.; De Biase, P.M.; Di Lella, S.; González Lebrero, M.C.; Martí, M.a.; Nadra, A.D.; Perissinotti, L.L.; et al. Modeling heme proteins using atomistic simulations. Phys. Chem. Chem. Phys. 2006, 8, 5611–5628. [Google Scholar] [CrossRef]
  9. Uchida, T.; Sato, E.; Sato, A.; Sagami, I.; Shimizu, T.; Kitagawa, T. CO-dependent activity-controlling mechanism of heme-containing CO-sensor protein, neuronal PAS domain protein 2. J. Biol. Chem. 2005, 280, 21358–21368. [Google Scholar] [CrossRef]
  10. Derbyshire, E.R.; Marletta, M.A. Biochemistry of Soluble Guanylate Cyclase. In cGMP: Generators, Effectors and Therapeutic Implications; Handbook of Experimental Pharmacology; Springer: Berlin/Heidelberg, Germany, 2009; Volume 191, pp. 17–32. [Google Scholar]
  11. Enemark, J.H.; Feltham, R.D. Principles of structure, bonding, and reactivity for metal nitrosyl complexes. Coord. Chem. Rev. 1974, 13, 339–406. [Google Scholar] [CrossRef]
  12. Karow, D.S.; Pan, D.; Tran, R.; Pellicena, P.; Presley, A.; Mathies, R.A.; Marletta, M.A. Spectroscopic characterization of the soluble guanylate cyclase-like heme domains from Vibrio cholerae and Thermoanaerobacter tengcongensis. Biochemistry 2004, 43, 10203–10211. [Google Scholar] [CrossRef] [PubMed]
  13. Ballou, D.P.; Zhao, Y.; Brandish, P.E.; Marletta, M.A. Revisiting the kinetics of nitric oxide (NO) binding to soluble guanylate cyclase: The simple NO-binding model is incorrect. Proc. Natl. Acad. Sci. USA 2002, 99, 12097–12101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Frenking, G.; Fröhlich, N. The Nature of the Bonding in Transition-Metal Compounds. Chem. Rev. 2000, 100, 717–774. [Google Scholar] [CrossRef] [PubMed]
  15. Rovira, C.; Kunc, K.; Hutter, J.; Ballone, P.; Parrinello, M. Equilibrium geometries and electronic structure of iron-porphyrin complexes: A density functional study. J. Phys. Chem. A 1997, 101, 8914–8925. [Google Scholar] [CrossRef]
  16. Ellison, M.K.; Schulz, C.E.; Scheidt, W.R. Nitrosyliron(III) porphyrinates: Porphyrin core conformation and FeNO geometry. Any correlation? J. Am. Chem. Soc. 2002, 124, 13833–13841. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Skolnick, J. Scoring function for automated assessment of protein structure template quality. Proteins Struct. Funct. Genet. 2004, 57, 702–710. [Google Scholar] [CrossRef] [PubMed]
  18. Li, P.; Merz, K.M. MCPB.py: A Python Based Metal Center Parameter Builder. J. Chem. Inf. Model. 2016, 56, 599–604. [Google Scholar] [CrossRef]
  19. Goddard, T.D.; Huang, C.C.; Ferrin, T.E. Visualizing density maps with UCSF Chimera. J. Struct. Biol. 2007, 157, 281–287. [Google Scholar] [CrossRef] [PubMed]
  20. Seminario, J.M. Calculation of intramolecular force fields from second-derivative tensors. Int. J. Quantum Chem. 1996, 60, 1271–1277. [Google Scholar] [CrossRef]
  21. Roe, D.R.; Cheatham, T.E. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef]
  22. Schelvis, J.P.M.; Seibold, S.A.; Cerda, J.F.; Garavito, R.M.; Babcock, G.T. Interaction of Nitric Oxide with Prostaglandin Endoperoxide H Synthase-1: Implications for Fe−His Bond Cleavage in Heme Proteins. J. Phys. Chem. B 2000, 104, 10844–10850. [Google Scholar] [CrossRef]
  23. Martí, M.A.; Capece, L.; Crespo, A.; Doctorovich, F.; Estrin, D.A. Nitric oxide interaction with cytochrome c’ and its relevance to guanylate cyclase. Why does the iron histidine bond break? J. Am. Chem. Soc. 2005, 127, 7721–7728. [Google Scholar] [CrossRef] [PubMed]
  24. Zhao, Y.; Brandish, P.; Ballou, D.; Marletta, M. A Molecular Basis for Nitric Oxide Sensing by Soluble Guanylate Cyclase. Proc. Natl. Acad. Sci. USA 1999, 96, 14753–14758. [Google Scholar] [CrossRef] [PubMed]
  25. Underbakke, E.S.; Iavarone, A.T.; Marletta, M.A. Higher-order interactions bridge the nitric oxide receptor and catalytic domains of soluble guanylate cyclase. Proc. Natl. Acad. Sci. USA 2013, 110, 6777–6782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Laskowski, R.A.; Swindells, M.B. LigPlot+: Multiple ligand-protein interaction diagrams for drug discovery. J. Chem. Inf. Model. 2011. [Google Scholar] [CrossRef] [PubMed]
  27. Teeter, M.M.; Case, D.A. Harmonic and quasiharmonic descriptions of crambin. J. Phys. Chem. 1990, 94, 8091–8097. [Google Scholar] [CrossRef]
  28. Tai, K.; Shen, T.; Henchman, R.H.; Bourne, Y.; Marchot, P.; McCammon, J.A. Mechanism of acetylcholinesterase inhibition by fasciculin: A 5-ns molecular dynamics simulation. J. Am. Chem. Soc. 2002, 124, 6153–6161. [Google Scholar] [CrossRef]
  29. Stone, J.R.; Marletta, M.A. Spectral and kinetic studies on the activation of soluble guanylate cyclase by nitric oxide. Biochemistry 1996, 35, 1093–1099. [Google Scholar] [CrossRef]
  30. Capece, L.; Estrin, D.A.; Marti, M.A. Dynamical characterization of the heme NO oxygen binding (HNOX) domain. Insight into soluble guanylate cyclase allosteric transition. Biochemistry 2008, 47, 9416–9427. [Google Scholar] [CrossRef]
  31. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990. [Google Scholar] [CrossRef]
  32. Giuili, G.; Scholl, U.; Bulle, F.; Guellaën, G. Molecular cloning of the cDNAs coding for the two subunits of soluble guanylyl cyclase from human brain. FEBS Lett. 1992, 304, 83–88. [Google Scholar] [CrossRef] [Green Version]
  33. Webb, B.; Sali, A. Comparative protein structure modeling using MODELLER. Curr. Protoc. Bioinform. 2014, 2014, 5.6.1–5.6.32. [Google Scholar] [CrossRef]
  34. Lovell, S.C.; Davis, I.W.; Arendall, W.B.; de Bakker, P.I.W.; Word, J.M.; Prisant, M.G.; Richardson, J.S.; Richardson, D.C. Structure validation by Calpha geometry: Phi, psi and Cbeta deviation. Proteins 2003, 50, 437–450. [Google Scholar] [CrossRef]
  35. Krieger, E.; Koraimann, G.; Vriend, G. Increasing the precision of comparative models with YASARA NOVA—A self-parameterizing force field. Proteins Struct. Funct. Genet. 2002, 47, 393–402. [Google Scholar] [CrossRef] [PubMed]
  36. Chen, V.B.; Arendall, W.B.; Headd, J.J.; Keedy, D.A.; Immormino, R.M.; Kapral, G.J.; Murray, L.W.; Richardson, J.S.; Richardson, D.C. MolProbity: All-atom structure validation for macromolecular crystallography. Acta Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 12–21. [Google Scholar] [CrossRef] [PubMed]
  37. DeLano, W.L. The PyMOL Molecular Graphics System. San Carlos, CA, USA. 2002. Available online: http://www.pymol.org (accessed on 2 February 2019).
  38. Case, D.A.; Betz, R.M.; Botello-Smith, W.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A., III; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; et al. AMBER 2016, University of California, San Francisco; University of California: San Francisco, CA, USA, 2016. [Google Scholar]
  39. Tirado-Rives, J.; Jorgensen, W.L. Performance of B3LYP density functional methods for a large set of organic molecules. J. Chem. Theory Comput. 2008, 4, 297–306. [Google Scholar] [CrossRef]
  40. Chakravorty, D.K.; Wang, B.; Lee, C.W.; Giedroc, D.P.; Merz, K.M. Simulations of allosteric motions in the zinc sensor CzrA. J. Am. Chem. Soc. 2012, 134, 3367–3376. [Google Scholar] [CrossRef] [PubMed]
  41. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09; Revision D.01; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar] [CrossRef]
  42. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef]
  43. Mark, P.; Nilsson, L. Structure and dynamics of the TIP3P, SPC, and SPC/E water models at 298 K. J. Phys. Chem. A 2001, 105, 9954–9960. [Google Scholar] [CrossRef]
  44. Darden, T.; Perera, L.; Li, L.; Lee, P. New tricks for modelers from the crystallography toolkit: The particle mesh Ewald algorithm and its use in nucleic acid simulations. Structure 1999, 7. [Google Scholar] [CrossRef]
  45. Gonnet, P. P-SHAKE: A quadratically convergent SHAKE in O (n2). J. Comput. Phys. 2007, 220, 740–750. [Google Scholar] [CrossRef]
  46. Amadei, A.; Linssen, A.B.M.; Berendsen, H.J.C. Essential dynamics of proteins. Proteins Struct. Funct. Bioinform. 1993, 17, 412–425. [Google Scholar] [CrossRef] [PubMed]
  47. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  48. Turner, P.J. XMGRACE; Version 5.1.19; Center for Coastal and Land-Margin Research, Oregon Graduate Institute of Science and Technology: Beaverton, OR, USA, 2005. [Google Scholar]
  49. Williams, T.; Kelley, C. Gnuplot 4.4: An Interactive Plotting Program. 2010. Available online: http://gnuplot.info (accessed on 2 February 2019).
  50. Van Aalten, D.M.; Findlay, J.B.; Amadei, A.; Berendsen, H.J. Essential dynamics of the cellular retinol-binding protein—Evidence for ligand-induced conformational changes. Protein Eng. 1995, 8, 1129–1135. [Google Scholar] [CrossRef] [PubMed]
  51. Bakan, A.; Meireles, L.M.; Bahar, I. ProDy: Protein dynamics inferred from theory and experiments. Bioinformatics 2011, 27, 1575–1577. [Google Scholar] [CrossRef] [PubMed]
  52. Khalid, R.R.; Tahir, M.; Maryam, A. Comparative studies of dynamics effects of BAY60-2770 and BAY58-2667 Binding with Human and Bacterial H-NOX domains. Molecules 2018, 23, 2141. [Google Scholar] [CrossRef] [PubMed]
Ijms 20 00698 g001 550
Figure 1. Predicted model of the 6c coordination state of the Human HNOX domain where the binding pocket heme is ligated to NO (6c NO-HNOX). The Y135-S137-R139 (Y-S-R) motif that bears strong interactions with the heme propionic moiety is also labeled (yellow). Functionally critical regions are also labeled as helix-f (red), loop α, β, and γ (yellow). A simple chemical scheme of heme is shown in the right panel. A detailed scheme of the heme interactions for each system are shown in Figure 7.
Figure 1. Predicted model of the 6c coordination state of the Human HNOX domain where the binding pocket heme is ligated to NO (6c NO-HNOX). The Y135-S137-R139 (Y-S-R) motif that bears strong interactions with the heme propionic moiety is also labeled (yellow). Functionally critical regions are also labeled as helix-f (red), loop α, β, and γ (yellow). A simple chemical scheme of heme is shown in the right panel. A detailed scheme of the heme interactions for each system are shown in Figure 7.
Ijms 20 00698 g001
Ijms 20 00698 g002 550
Figure 2. HNOX active site models used in the parametrization procedure. (a,b), (d,e), and (g,h) represent the small and large model systems of active sites His-heme-NO, His-heme-O2 and His-heme-CO, respectively. Panels (c,f,i) show the optimized active site of the smaller models (a,d,g) that was further used for the parameterization of His-heme-NO/O2/CO.
Figure 2. HNOX active site models used in the parametrization procedure. (a,b), (d,e), and (g,h) represent the small and large model systems of active sites His-heme-NO, His-heme-O2 and His-heme-CO, respectively. Panels (c,f,i) show the optimized active site of the smaller models (a,d,g) that was further used for the parameterization of His-heme-NO/O2/CO.
Ijms 20 00698 g002
Ijms 20 00698 g003 550
Figure 3. Time-dependent RMSD plot of the backbone atoms of apo, 5c NO-HNOX, 6c NO-HNOX, and 6c O2/CO-HNOX. Magenta, black, red, green and blue represent the RMSD plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX and 6c CO-HNOX, respectively.
Figure 3. Time-dependent RMSD plot of the backbone atoms of apo, 5c NO-HNOX, 6c NO-HNOX, and 6c O2/CO-HNOX. Magenta, black, red, green and blue represent the RMSD plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX and 6c CO-HNOX, respectively.
Ijms 20 00698 g003
Ijms 20 00698 g004 550
Figure 4. The plot depicts the residual fluctuations of human HNOX domain (1–186 a.a.) under the influence of gaseous ligands NO, O2 and CO. The x-axis shows the residue number, as well as the secondary structure elements of the H-NOX domain; helices are represented as red springs (HA, HB, HC, HD, HE, HF, HG), loops as green lines (α, β, γ, δ), and strands as yellow arrows (S1, S2, S3, S4). The y-axis represents the residual mobility around its mean position. Magenta, black, red, green, and blue represent the RMSD plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX and 6c CO-HNOX, respectively.
Figure 4. The plot depicts the residual fluctuations of human HNOX domain (1–186 a.a.) under the influence of gaseous ligands NO, O2 and CO. The x-axis shows the residue number, as well as the secondary structure elements of the H-NOX domain; helices are represented as red springs (HA, HB, HC, HD, HE, HF, HG), loops as green lines (α, β, γ, δ), and strands as yellow arrows (S1, S2, S3, S4). The y-axis represents the residual mobility around its mean position. Magenta, black, red, green, and blue represent the RMSD plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX and 6c CO-HNOX, respectively.
Ijms 20 00698 g004
Ijms 20 00698 g005 550
Figure 5. Time-dependent plot of the radius of gyration (Rg) for apo, 6c NO/CO/O2-HNOX and 5c NO-HNOX systems. Magenta, black, red, green and blue represent the plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX, and 6c CO-HNOX, respectively.
Figure 5. Time-dependent plot of the radius of gyration (Rg) for apo, 6c NO/CO/O2-HNOX and 5c NO-HNOX systems. Magenta, black, red, green and blue represent the plots of apo, 5c NO-HNOX, 6c NO-HNOX, 6c O2-HNOX, and 6c CO-HNOX, respectively.
Ijms 20 00698 g005
Ijms 20 00698 g006 550
Figure 6. Binding modes of (a) 5c NO-HNOX, (b) 6c NO-HNOX, (c) O2-HNOX and (d) CO-HNOX complexes with heme carboxyl groups. The x-axis represents the time (ns) while the y-axis shows the heme (HM) propionic moiety (O1A, O2A, O1D, O2D) and critical binding pocket residues atoms interacting by hydrogen bonds. Yellow color illustrates the presence of hydrogen bond while black color represents absence.
Figure 6. Binding modes of (a) 5c NO-HNOX, (b) 6c NO-HNOX, (c) O2-HNOX and (d) CO-HNOX complexes with heme carboxyl groups. The x-axis represents the time (ns) while the y-axis shows the heme (HM) propionic moiety (O1A, O2A, O1D, O2D) and critical binding pocket residues atoms interacting by hydrogen bonds. Yellow color illustrates the presence of hydrogen bond while black color represents absence.
Ijms 20 00698 g006
Ijms 20 00698 g007 550
Figure 7. (a–d) represents the non-bonded interactions of 5c NO-HNOX, 6c NO-HNOX, O2-HNOX and CO-HNOX systems respectively. The propionic moiety of heme (O1A, O2A, O1D, O2D) prominently interacts with functionally important residues such as Y112, Y135, S137, R139, R116 and Y2 throughout the simulation.
Figure 7. (a–d) represents the non-bonded interactions of 5c NO-HNOX, 6c NO-HNOX, O2-HNOX and CO-HNOX systems respectively. The propionic moiety of heme (O1A, O2A, O1D, O2D) prominently interacts with functionally important residues such as Y112, Y135, S137, R139, R116 and Y2 throughout the simulation.
Ijms 20 00698 g007
Ijms 20 00698 g008 550
Figure 8. Porcupine plot illustrating the highest principal component variance of the first (PC-1) and second (PC-2) modes of motion. Highly fluctuating regions are highlighted in red. The length of the cyan “needles” represents the degree of mobility.
Figure 8. Porcupine plot illustrating the highest principal component variance of the first (PC-1) and second (PC-2) modes of motion. Highly fluctuating regions are highlighted in red. The length of the cyan “needles” represents the degree of mobility.
Ijms 20 00698 g008
Ijms 20 00698 g009 550
Figure 9. Superposition of HNOX-NO (yellow) and O2 (magenta) systems. The heme-O2 ligand is shown as transparent spheres while the heme-NO ligand is shown as sticks. Regions assumed to be important for signal transduction (helix-f, loops α, β and γ) experienced noticeable structural changes after ligation of NO to HNOX.
Figure 9. Superposition of HNOX-NO (yellow) and O2 (magenta) systems. The heme-O2 ligand is shown as transparent spheres while the heme-NO ligand is shown as sticks. Regions assumed to be important for signal transduction (helix-f, loops α, β and γ) experienced noticeable structural changes after ligation of NO to HNOX.
Ijms 20 00698 g009
Table
Table 1. Comparison of the apo-HNOX and HNOX-NO/O2/CO-ligated complexes (5c, 6c heme coordination states) hydrogen bond interactions of the heme propionic moiety with binding pocket residues.
Table 1. Comparison of the apo-HNOX and HNOX-NO/O2/CO-ligated complexes (5c, 6c heme coordination states) hydrogen bond interactions of the heme propionic moiety with binding pocket residues.
ComplexesAcceptorDonorPercent of Frames of 30 ns Simulation Showing the Hydrogen BondAverage DistanceAverage Angle
Apo-HNOXHM-187@OR139@N96.862.79155.98
HM-187@OS137@OH94.792.72156.46
HM-187@OY135@OH87.412.70158.41
HM-187@OY2@N44.392.84153.42
HM-187@OR116@N28.942.86159.21
HM-187@OY112@OH2.062.76164.44
5c HNOX-NOHM-187@OR116@N54.752.92145.26
HM-187@OS137@OH35.372.81150.52
HM-187@OR139@N30.232.90149.28
HM-187@OL142@N4.742.90155.25
HM-187@OY49@OH0.822.85147.41
HM-187@OY135@OH0.182.84148.21
HM-187@OE138@N0.072.89141.79
6c HNOX-NOHM-187@OR139@N96.622.89145.54
HM-187@OS137@OH72.482.85150.51
HM-187@OY135@OH40.072.85150.70
HM-187@OR116@N9.332.91146.68
HM-187@OA117@N3.262.88142.18
HM-187@OE138@N1.302.93155.31
6c HNOX-O2HM-187@OR139@N93.062.82152.35
HM-187@OS137@OH91.082.73156.89
HM-187@OY135@OH87.662.69160.04
HM-187@OY2@N17.572.82153.32
HM-187@OR116@N8.722.85156.24
HM-187@OY112@OH1.342.73165.64
6c HNOX-COHM-187@OR116@N90.822.79155.25
HM-187@OR139@N89.552.74154.96
HM-187@OS137@OH82.142.64165.79
HM-187@OY135@OH8.112.65165.85
HM-187@OY112@OH3.662.56159.99
HM-187@OE138@N1.462.85158.29
HM-187@OM1@N1.192.78143.47
HM-187@OY2@N0.342.85156.35
Table
Table 2. Hydrogen bond occupancy of the histidine imidazole ring in the 5c NO-HNOX system.
Table 2. Hydrogen bond occupancy of the histidine imidazole ring in the 5c NO-HNOX system.
AcceptorDonorFramesFractionDistance ÅAngle
L101@OH105@N637800.82832.7872163.1623
H105@OA109@N391300.50822.8682158.5987
H105@OL108@N47460.06162.9051150.2395
P118@OH105@N27490.03572.8107155.5122
D102@OH105@N47220.06142.8373151.0958

Share and Cite

MDPI and ACS Style

Khalid, R.R.; Siddiqi, A.R.; Mylonas, E.; Maryam, A.; Kokkinidis, M. Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands. Int. J. Mol. Sci. 2019, 20, 698. https://doi.org/10.3390/ijms20030698

AMA Style

Khalid RR, Siddiqi AR, Mylonas E, Maryam A, Kokkinidis M. Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands. International Journal of Molecular Sciences. 2019; 20(3):698. https://doi.org/10.3390/ijms20030698

Chicago/Turabian Style

Khalid, Rana Rehan, Abdul Rauf Siddiqi, Efstratios Mylonas, Arooma Maryam, and Michael Kokkinidis. 2019. "Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands" International Journal of Molecular Sciences 20, no. 3: 698. https://doi.org/10.3390/ijms20030698

APA Style

Khalid, R. R., Siddiqi, A. R., Mylonas, E., Maryam, A., & Kokkinidis, M. (2019). Dynamic Characterization of the Human Heme Nitric Oxide/Oxygen (HNOX) Domain under the Influence of Diatomic Gaseous Ligands. International Journal of Molecular Sciences, 20(3), 698. https://doi.org/10.3390/ijms20030698

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop