Cryo-EM density map fitting driven in-silico structure of human soluble guanylate cyclase (hsGC) reveals functional aspects of inter-domain cross talk upon NO binding

The Biologically Relevant Coordination Chemistry of Iron and Nitric Oxide: Electronic Structure and Reactivity

Nitric oxide (NO) is an important signaling molecule that is involved in a wide range of physiological and pathological events in biology. Metal coordination chemistry, especially with iron, is at the heart of many biological transformations involving NO. A series of heme proteins, nitric oxide synthases (NOS), soluble guanylate cyclase (sGC), and nitrophorins, are responsible for the biosynthesis, sensing, and transport of NO. Alternatively, NO can be generated from nitrite by heme- and copper-containing nitrite reductases (NIRs). The NO-bearing small molecules such as nitrosothiols and dinitrosyl iron complexes (DNICs) can serve as an alternative vehicle for NO storage and transport. Once NO is formed, the rich reaction chemistry of NO leads to a wide variety of biological activities including reduction of NO by heme or non-heme iron-containing NO reductases and protein post-translational modifications by DNICs. Much of our understanding of the reactivity of metal sites in biology with NO and the mechanisms of these transformations has come from the elucidation of the geometric and electronic structures and chemical reactivity of synthetic model systems, in synergy with biochemical and biophysical studies on the relevant proteins themselves. This review focuses on recent advancements from studies on proteins and model complexes that not only have improved our understanding of the biological roles of NO but also have provided foundations for biomedical research and for bio-inspired catalyst design in energy science.

A recursive molecular docking coupled with energy-based pose-rescoring and MD simulations to identify hsGC βH-NOX allosteric modulators for cardiovascular dysfunctions

Modulating the activity of human soluble guanylate cyclase (hsGC) through allosteric regulation of the βH-NOX domain has been considered as an immediate treatment for cardiovascular disorder (CVDs). Currently available βH-NOX domain-specific agonists including cinaciguat are unable to deal with the conundrum raised due to oxidative stress in the case of CVDs and their associated comorbidities. Therefore, the idea of investigating novel compounds for allosteric regulation of hsGC activation has been rekindled to circumvent CVDs. Current study aims to identify novel βH-NOX domain-specific compounds that can selectively turn on sGC functions by modulating the conformational dynamics of the target protein. Through a comprehensive computational drug-discovery approach, we first executed a target-based performance assessment of multiple docking (PLANTS, QVina, LeDock, Vinardo, Smina) scoring functions based on multiple performance metrices. QVina showed the highest capability of selecting true-positive ligands over false positives thus, used to screen 4.8 million ZINC15 compounds against βH-NOX domain. The docked ligands were further probed in terms of contact footprint and pose reassessment through clustering analysis and PLANTS docking, respectively. Subsequently, energy-based AMBER rescoring of top 100 low-energy complexes, per-residue energy decomposition analysis, and ADME-Tox analysis yielded the top three compounds i.e. ZINC000098973660, ZINC001354120371, and ZINC000096022607. The impact of three selected ligands on the internal structural dynamics of the βH-NOX domain was also investigated through molecular dynamics simulations. The study revealed potential electrostatic interactions for better conformational dialogue between βH-NOX domain and allosteric ligands that are critical for the activation of hsGC as compared to the reference compound.

Soluble Guanylate Cyclase Inhibitors Discovered among Natural Compounds

Soluble guanylate cyclase (sGC) is the human receptor of nitric oxide (NO) in numerous kinds of cells and produces the second messenger 3',5'-cyclic guanosine monophosphate (cGMP) upon NO binding to its heme. sGC is involved in many cell signaling pathways both under healthy conditions and under pathological conditions, such as angiogenesis associated with tumor growth. Addressing the selective inhibition of the NO/cGMP pathway is a strategy worthwhile to be investigated for slowing down tumoral angiogenesis or for curing vasoplegia. However, sGC inhibitors are lacking investigation. We have explored a chemical library of various natural compounds and have discovered inhibitors of sGC. The selected compounds were evaluated for their inhibition of purified sGC in vitro and sGC in endothelial cells. Six natural compounds, from various organisms, have IC₅₀ in the range 0.2-1.5 μM for inhibiting the NO-activated synthesis of cGMP by sGC, and selected compounds exhibit a quantified antiangiogenic activity using an endothelial cell line. These sGC inhibitors can be used directly as tools to investigate angiogenesis and cell signaling or as templates for drug design.

Solving the α-helix correspondence problem at medium-resolution Cryo-EM maps through modeling and 3D matching

Cryo-electron microscopy (cryo-EM) has recently emerged as a prominent biophysical method for macromolecular structure determination. Many research efforts have been devoted to produce cryo-EM images, density maps, at near-atomic resolution. Despite many advances in technology, the resolution of the generated density maps may not be sufficiently adequate and informative to directly construct the atomic structure of proteins. At medium-resolution (∼4-10 Å), secondary structure elements (α-helices and β-sheets) are discernible, whereas finding the correspondence of secondary structure elements detected in the density map with those on the sequence remains a challenging problem. In this paper, an automatic framework is proposed to solve α-helix correspondence problem in three-dimensional space. Through modeling of the sequence with the aid of a novel strategy, the α-helix correspondence problem is initially transformed into a complete weighted bipartite graph matching problem. An innovative correlation-based scoring function based on a well-known and robust statistical method is proposed for weighting the graph. Moreover, two local optimization algorithms, which are Greedy and Improved Greedy algorithms, have been presented to find α-helix correspondence. A widely used data set including 16 reconstructed and 4 experimental cryo-EM maps were chosen to verify the accuracy and reliability of the proposed automatic method. The experimental results demonstrate that the automatic method is highly efficient (86.25% accuracy), robust (11.3% error rate), fast (∼1.4 s), and works independently from cryo-EM skeleton.

A new paradigm for gaseous ligand selectivity of hemoproteins highlighted by soluble guanylate cyclase

Nitric oxide (NO), carbon monoxide (CO), and oxygen (O₂) are important physiological messengers whose concentrations vary in a remarkable range, [NO] typically from nM to several μM while [O₂] reaching to hundreds of μM. One of the machineries evolved in living organisms for gas sensing is sensor hemoproteins whose conformational change upon gas binding triggers downstream response cascades. The recently proposed "sliding scale rule" hypothesis provides a general interpretation for gaseous ligand selectivity of hemoproteins, identifying five factors that govern gaseous ligand selectivity. Hemoproteins have intrinsic selectivity for the three gases due to a neutral proximal histidine ligand while proximal strain of heme and distal steric hindrance indiscriminately adjust the affinity of these three gases for heme. On the other hand, multiple-step NO binding and distal hydrogen bond donor(s) specifically enhance affinity for NO and O₂, respectively. The "sliding scale rule" hypothesis provides clear interpretation for dramatic selectivity for NO over O₂ in soluble guanylate cyclase (sGC) which is an important example of sensor hemoproteins and plays vital roles in a wide range of physiological functions. The "sliding scale rule" hypothesis has so far been validated by all experimental data and it may guide future designs for heme-based gas sensors.

Exploring the structural basis of conformational heterogeneity and autoinhibition of human cGMP-specific protein kinase Iα through computational modelling and molecular dynamics simulations

Protein kinase Iα (PKGIα) is a pivotal cyclic guanosine monophosphate (cGMP) signalling protein. Major steps related to the structural plasticity of PKGIα have been inferred but the structural aspects of the auto-inhibition and multidomain tertiary organization of human PKGIα in active and inactive form are not clear. Here we combine computational comparative modelling, protein-protein docking and molecular dynamics (MD) simulations to investigate structural details of the repressed state of the catalytic domain of PKGIα. Exploration of the potential inhibitory conformation of the auto-inhibitory domain (AI) within the catalytic cleft reveals that the pseudo-substrate motif binds with residues of the glycine rich loop and substrate-binding lobe. Dynamic changes as a result of coupling of the catalytic and AI domains are also investigated. The three-dimensional homodimeric models of PKGIα in the active and inactive state indicate that PKGIα in its inactive-state attains a compact globular structure where cyclic nucleotide binding (CNB-A/B) domains are buried, whereas the catalytic domains are inaccessible with their substrate-binding pockets facing the N-terminal of CNB-A. Contrary to this, the active-state model of PKGIα shows an extended conformation where CNB-A/B domains are slightly rearranged and the catalytic domains of homodimer flanking the C-terminal with their substrate binding lobes free to entrap downstream proteins. These findings are consistent with previously reported static images of the multidomain organization of PKGIα. Structural insights pertaining to the conformational heterogeneity and auto-inhibition of PKGIα provided in this study may help to understand the dynamics-driven effective regulation of PKGIα.

Probing the Structural Dynamics of the Catalytic Domain of Human Soluble Guanylate Cyclase

In the nitric oxide (NO) signaling pathway, human soluble guanylate cyclase (hsGC) synthesizes cyclic guanosine monophosphate (cGMP); responsible for the regulation of cGMP-specific protein kinases (PKGs) and phosphodiesterases (PDEs). The crystal structure of the inactive hsGC cyclase dimer is known, but there is still a lack of information regarding the substrate-specific internal motions that are essential for the catalytic mechanism of the hsGC. In the current study, the hsGC cyclase heterodimer complexed with guanosine triphosphate (GTP) and cGMP was subjected to molecular dynamics simulations, to investigate the conformational dynamics that have functional implications on the catalytic activity of hsGC. Results revealed that in the GTP-bound complex of the hsGC heterodimer, helix 1 of subunit α (α:h1) moves slightly inwards and comes close to helix 4 of subunit β (β:h4). This conformational change brings loop 2 of subunit β (β:L2) closer to helix 2 of subunit α (α:h2). Likewise, loop 2 of subunit α (α:L2) comes closer to helix 2 of subunit β (β:h2). These structural events stabilize and lock GTP within the closed pocket for cyclization. In the cGMP-bound complex, α:L2 detaches from β:h2 and establishes interactions with β:L2, which results in the loss of global structure compactness. Furthermore, with the release of pyrophosphate, the interaction between α:h1 and β:L2 weakens, abolishing the tight packing of the binding pocket. This study discusses the conformational changes induced by the binding of GTP and cGMP to the hsGC catalytic domain, valuable in designing new therapeutic strategies for the treatment of cardiovascular diseases.