Adenylyl cyclases (ACs) are membrane bound enzymes that catalyze the production of cAMP from ATP in response to the activation by G-protein Gαs. Different isoforms of ACs are ubiquitously expressed in different tissues involved in regulatory mechanisms in response to specific stimulants. There are 9 AC isoforms present in humans, with AC5 and AC6 proposed to play a vital role in cardiac functions. The activity of AC6 is sensitive to nitric oxide, such that nitrosylation of the protein might regulate its function. However, the information on structural determinants of nitrosylation in ACs and how they interact with Gαs is limited. Here we used homology modeling to build a molecular model of human AC6 bound to Gαs. Based on this 3D model, we predict the nitrosylation amenable cysteines, and identify potential novel ligands of AC6 using virtual ligand screening. Our model suggests Cys1004 in AC6 (subunit C2) and Cys174 in Gαs present at the AC-Gαs interface as the possible residues that might undergo reversible nitrosylation. Docking analysis predicted novel ligands of AC6 that include forskolin-based compounds and its derivatives. Further work involving site-directed mutagenesis of the predicted residues will allow manipulation of AC activity using novel ligands, and crucial insights on the role of nitrosylation of these proteins in pathophysiological conditions.