The Indian hedgehog protein belongs to the Hedgehog family of signaling molecules, a group of secreted morphogens that is conserved across vertebrates and includes Sonic hedgehog, Indian hedgehog, and Desert hedgehog. While Sonic hedgehog is the most prevalent member and governs limb growth, neural patterning, and organ development, Indian hedgehog is the star of the skeletal show. It regulates chondrocyte differentiation, ossification, and skeletal morphogenesis, and it does so partly through its interaction with parathyroid hormone-related peptide, a molecule that controls the proliferation and maturation of cartilage cells in developing fetal bones. Mice engineered to lack Indian hedgehog display profound abnormalities in bone formation, and in humans, mutations in the IHH gene have already been linked to brachydactyly type A-1, a condition marked by shortened fingers and toes, and to acrocapitofemoral dysplasia, characterized by short limbs, a relatively large head, and a narrow thorax. One well-known example is the Asn100Asp substitution, which is associated with brachydactyly type A1.

The starting point for the new analysis was the dbSNP database, the National Center for Biotechnology Information’s comprehensive catalog of human genetic variation. The researchers identified a total of 436 single nucleotide polymorphisms within the IHH gene, of which 280 were classified as non-synonymous, meaning that the single-letter DNA change actually alters the amino acid sequence of the resulting protein. Non-synonymous variants are the most consequential class of common genetic variation because they can directly reshape a protein’s structure and function, although many turn out to be harmless. The challenge, and the central aim of the study, was to separate the potentially disease-causing variants from the benign majority using computational methods that are far faster and cheaper than laboratory experiments.

To accomplish this triage, the team deployed a battery of seven independent prediction algorithms: PANTHER, FATHMM, SIFT, PolyPhen, Mutation Assessor, SNPs&GO, and Provean. Each tool interrogates a mutation from a different angle. FATHMM relies on a hidden Markov model and flags variants with scores above 0.4 as potentially pathogenic. PolyPhen evaluates sequence homology and assigns probability scores, treating values above 0.5 as deleterious. SIFT weighs sequence conservation and amino acid properties, classifying mutations with scores below 0.05 as harmful, while Provean applies a cutoff of minus 2.5. The researchers adopted a deliberately stringent consensus rule: only variants predicted to be damaging by at least five of the seven tools advanced to the next stage. This approach, previously applied to genes such as HBA1, BRCA1, and BRCA2, narrowed the field from 280 non-synonymous variants to 116 candidates with a high likelihood of functional impact.

Structure and stability came next. A protein’s three-dimensional fold is essential to its function, and mutations that destabilize that fold can trigger misfolding and loss of activity. The team used two complementary tools, I-Mutant 3.0 and DynaMUT, to estimate how each candidate mutation would change the protein’s free energy of stability. I-Mutant 3.0, which is built on a support vector machine with roughly 80 percent accuracy, classifies mutations as neutral, destabilizing, or stabilizing based on predicted changes in free energy, while DynaMUT combines graph-based signatures with normal mode dynamics to reach a consensus verdict. The verdict was striking: 103 to 104 of the candidate mutations were predicted to destabilize the protein. Even more telling, 53 of the 116 deleterious variants fell within functionally significant regions of the protein, particularly the stretch between amino acids 28 and 202, which is involved in metal binding, cholesterol transfer, lipidation, glycosylation, and autocleavage, with calcium and zinc ions contributing to structural integrity.

To visualize what these mutations actually do to the protein’s architecture, the researchers retrieved the native structure of the N-terminal signaling domain of Indian hedgehog from the Protein Data Bank, using the structure with the identifier 3K7G. This domain is the business end of the protein, the part that binds the PTCH1 receptor on target cells, so mutations here are especially likely to disrupt signaling. Mutant models were generated by homology modeling with SWISS-MODEL, achieving more than 90 percent sequence identity with the native protein, and were validated with PROCHECK, which assesses stereochemical quality through Ramachandran plots. Most mutant models placed more than 90 percent of residues in favored regions with G-factor values above 0.5, confirming that the models were reliable enough for further comparison. Superimposing the mutant structures onto the native one with PyMOL revealed root-mean-square deviation values exceeding 1 angstrom for many variants, indicating significant structural deviation, while energy minimization analyses showed altered energy profiles that suggested the mutations could substantially reshape the protein’s interaction landscape.

The most revealing experiments were computational docking studies that asked how the mutations affect the handshake between Indian hedgehog and its receptor, PTCH1. Under normal conditions, PTCH1 suppresses a protein called Smoothened; when a Hedgehog ligand binds PTCH1, the brake is released and the signaling cascade proceeds. Using the PatchDock server, which exploits shape complementarity principles, and refining selected solutions with FIREDOCK, the team compared the binding behavior of native and mutant proteins. The results were nuanced. Several mutations, including G69S, L166M, G137D/E142K, and A184V/A184D, produced markedly higher geometric shape complementarity scores than the wild-type protein, with A184V and A184D showing the greatest potential to enhance Ihh-PTCH1 binding and possibly hyperactivate the downstream pathway, a scenario consistent with diseases of inappropriate bone growth. Meanwhile, mutations C28Y, Y49C, and E68G displayed strongly negative atomic contact energies, indicating favorable binding energetics, whereas R66H and R32P showed reduced values that could stabilize the complex in unusual ways. Intriguingly, some variants such as G29R and G202C improved structural complementarity while simultaneously worsening desolvation energy, a reminder that shape and energetics do not always move in lockstep. The overall pattern echoes the E95K mutation already associated with brachydactyly type A1, lending plausibility to the idea that altered receptor binding is a genuine disease mechanism.

The study also zoomed out to place Indian hedgehog in its cellular context. Using the STRING database, the researchers constructed a protein-protein interaction network containing Ihh and ten partner proteins, connected by 41 edges with a high average local clustering coefficient of 0.89 and a PPI enrichment p-value of 3.62e-13, indicating statistically significant connectivity. The network includes PTCH1 and PTCH2, the receptors that negatively regulate the pathway; Smoothened, the signal transmitter; HHIP, a negative regulator that restricts Hedgehog diffusion; DISP, which helps release Hedgehog from producing cells; the co-receptors Cdo and Boc; Hedgehog acyltransferase, which performs the palmitoylation essential for signaling activity; GAS1, a positive regulator; and PTHLH, the parathyroid hormone-related peptide that governs chondrocyte behavior at the growth plate. Mutations in many of these partners are independently associated with skeletal malformations, including Gorlin syndrome, and the authors argue that structural changes in Ihh could ripple through this entire network, disrupting the crosstalk between Hedgehog, Wnt, BMP, and Notch signaling that coordinates skeletal formation and repair.

The practical implications extend toward the clinic. Because the deleterious mutations appear to alter receptor binding and pathway activity, the authors suggest that the Ihh protein could serve as a potential drug target, and that restoring signaling balance with small molecules or biologics might one day help patients with IHH-related disorders. Alternative strategies could target downstream effectors of the Wnt or BMP pathways, activate SMAD proteins, or inhibit negative regulators such as HHIP. Biomarkers derived from gene expression signatures or protein interaction profiles altered by IHH mutations, including changes in PTHLH expression or GLI transcription factor activity, could support earlier diagnosis and personalized therapy. The authors are careful to note that these are computational predictions, and that structural and biochemical experiments will be needed to validate the functional consequences of the flagged variants. Even so, by distilling 436 raw variants down to a prioritized set of 53 mutations concentrated in the receptor-binding domain, the study provides exactly the kind of focused roadmap that experimental geneticists and drug developers need to translate a database of single-letter changes into a deeper understanding of skeletal disease.

Subject of Research: Bioinformatics prediction of pathogenic non-synonymous variants in the IHH gene linked to bone development abnormalities

Article Title: Predicting pathogenic variants in the IHH gene associated with bone development abnormalities: a bioinformatics approach

Article References: Sunkar, S., Pavankumar, G., Hariharan, V., Namrata, K., Nachiyar, C. V., & Athista, M. (2025). Predicting pathogenic variants in the IHH gene associated with bone development abnormalities: a bioinformatics approach. Discover Biotechnology, 2(1), Article 8. https://doi.org/10.1007/s44340-025-00016-z

Image Credits: AI Generated

DOI: 10.1007/s44340-025-00016-z

Keywords: IHH gene, Indian hedgehog, nsSNPs, bone development, brachydactyly, acrocapitofemoral dysplasia, Hedgehog signaling, PTCH1, protein stability, molecular docking, bioinformatics, skeletal dysplasia