Soluble guanylate cyclase (sGC), the main target of nitric oxide (NO), has been proven to have a significant role in coronary artery disease, pulmonary hypertension, erectile dysfunction, and myocardial infarction. One of its agonists, BAY 41-2272 (Riociguat), has been recently approved for treatment of pulmonary arterial hypertension (PHA), while some others are in clinical phases of development. However, the location of the binding sites for the two known types of agonists, heme-dependent stimulators and heme-independent activators, is a matter of debate, particularly for the first group where both a location on the regulatory (H-NOX) and on the catalytic domain have been suggested by different authors. Here, we address its potential location on the catalytic domain, the unique well characterized at the structural level, by an "in silico" approach. Homology models of the catalytic domain of sGC in "inactive" or "active" conformations were constructed using the structure of previously described crystals of the catalytic domains of "inactive" sGCs (2WZ1, 3ET6) and of "active" adenylate cyclase (1CJU). Each model was submitted to six independent molecular dynamics simulations of about 1 μs. Docking of YC-1, a classic heme-dependent stimulator, to all frames of representative trajectories of "inactive" and "active" conformations, followed by calculation of absolute binding free energies with the linear interaction energy (LIE) method, revealed a potential high-affinity binding site on the "active" structure. The site, located between the pseudo-symmetric and the catalytic site just over the loop β2 -β3 , does not overlap with the forskolin binding site on adenylate cyclases. Proteins 2016; 84:1534-1548. © 2016 Wiley Periodicals, Inc.