Structural, functional and molecular dynamics analysis of cathepsin B gene SNPs associated with tropical calcific pancreatitis, a rare disease of tropics

Tropical Calcific Pancreatitis (TCP) is a neglected juvenile form of chronic non-alcoholic pancreatitis. Cathepsin B (CTSB), a lysososmal protease involved in cellular degradation process, is recently been studied as a potential candidate gene in the pathogenesis of TCP. According to cathepsin B hypothesis, mutated CTSB can lead to premature intracellular activation of trypsinogen, which is a key regulatory mechanism in pancreatitis. So far, CTSB mutations have been studied in pancreatitis and neurodegenerative disorders but little is known about the structural and functional effect of variants in CTSB. In this study, we investigated the effect of single nucleotide variants (SNVs) associated with TCP, using molecular dynamics and simulation algorithms. There were two non-synonymous variants in the coding region (L26V and S53G) of CTSB, located in the propeptide region. We tried to predict the effect of these variants on structure and function using multiple algorithms: SIFT, Polyphen2, Panther, SDM sever, i-Mutant2.0 suite, mCSM algorithm and Vadar. Further, using databases like miRdbSNP, PolymiRTS and miRNASNP, two SNPs in 3’UTR region were predicted to affect the miRNA binding sites. Structural mutated models of nsSNP mutants (L26V and S53G) were prepared by MODELLER v9.15 and evaluated using TM-Align, Verify 3D, ProSA and Ramachandran plot. The results showed that the models (L26V and S53G) were of high accuracy. The 3D mutated structures were simulated using GROMACS 5.0 to predict the impact of these SNPs on protein stability. The results from in silico analysis and molecular dynamics simulations suggested that these variants in the propeptide region of cathepsin B could lead to structural and functional changes in the protein. Hence, the structural and functional analysis results have given interim conclusions that these variants can have deleterious effect in TCP and thus should be screen in samples from all TCP patients to decipher its distribution in patient population.