Molecular and in silico analyses validates pathogenicity of homozygous mutations in the NPR2 gene underlying variable phenotypes of Acromesomelic dysplasia, type Maroteaux

Homozygous and/or heterozygous loss of function mutations in the natriuretic peptide receptor B (NPR2) have been reported in causing acromesomelic dysplasia, type Maroteaux with variable clinical features and idiopathic short stature with nonspecific skeletal deformities. On the other hand, gain of function mutations in the same gene result in overgrowth disorder suggesting that NPR2 and its ligand, natriuretic peptide precursor C (CNP), are the key players of endochondral bone growth. However, the precise mechanism behind phenotypic variability of the NPR2 mutations is not fully understood so far. In the present study, three consanguineous families of Pakistani origin (A, B, C) with variable phenotypes of acromesomelic dysplasia, type Maroteaux were evaluated at clinical and molecular levels. Linkage analysis followed by Sanger sequencing of the NPR2 gene revealed three homozygous mutations including p.(Leu314 Arg), p.(Arg371*), and p.(Arg1032*) in family A, B and C, respectively. In silico structural and functional analyses substantiated that a novel missense mutation [p.(Leu314 Arg)] in family A allosterically affects binding of NPR2 homodimer to its ligand (CNP) which ultimately results in defective guanylate cyclase activity. A nonsense mutation [p.(Arg371*)] in family B entirely removed the transmembrane domain, protein kinase domain and guanylate cyclase domains of the NPR2 resulting in abolishing its guanylate cyclase activity. Another novel mutation [p.(Arg1032*)], found in family C, deteriorated the guanylate cyclase domain of the protein and probably plundered its guanylate cyclase activity. These results suggest that guanylate cyclase activity is the most critical function of the NPR2 and phenotypic severity of the NPR2 mutations is proportional to the reduction in its guanylate cyclase activity.