Rituraj Purohit

In silico investigation of molecular mechanism of laminopathy caused by a point mutation (R482W) in lamin A/C protein.

Lamin A/C proteins are the major components of a thin proteinaceous filamentous meshwork, the lamina, that underlies the inner nuclear membrane. A few specific mutations in the lamin A/C gene cause a disease with remarkably different clinical features: FPLD, or familial partial lipodystrophy (Dunnigan-type), which mainly affects adipose tissue. Lamin A/C mutant R482W is the key variant that causes FPLD. Biomolecular interaction and molecular dynamics (MD) simulation analysis were performed to understand dynamic behavior of native and mutant structures at atomic level. Mutant lamin A/C (R482W) showed more interaction with its biological partners due to its expansion of interaction surface and flexible nature of binding residues than native lamin A/C. MD simulation clearly indicates that the flexibility of interacting residues of mutant are mainly due to less involvement in formation of inter and intramolecular hydrogen bonds. Our analysis of native and Mutant lamin A/C clearly shows that the structural and functional consequences of the mutation R482W causes FPLD. Because of the pivotal role of lamin A/C in maintaining dynamics of nuclear function, these differences likely contribute to or represent novel mechanisms in laminopathy development.

Role of ELA region in auto-activation of mutant KIT receptor: a molecular dynamics simulation insight.

KIT receptor is the prime target in gastrointestinal stromal tumor (GISTs) therapy. Second generation inhibitor, Sunitinib, binds to an inactivated conformation of KIT receptor and stabilizes it in order to prevent tumor formation. Here, we investigated the dynamic behavior of wild type and mutant D816H KIT receptor, and emphasized the extended A-loop (EAL) region (805–850) by conducting molecular dynamics simulation (∼100 ns). We analyzed different properties such as root mean square cutoff or deviation, root mean square fluctuation, radius of gyration, solvent-accessible surface area, hydrogen bonding network analysis, and essential dynamics. Apart from this, clustering and cross-correlation matrix approach was used to explore the conformational space of the wild type and mutant EAL region of KIT receptor. Molecular dynamics analysis indicated that mutation (D816H) was able to alter intramolecular hydrogen bonding pattern and affected the structural flexibility of EAL region. Moreover, flexible secondary elements, specially, coil and turns were dominated in EAL region of mutant KIT receptor during simulation. This phenomenon increased the movement of EAL region which in turn helped in shifting the equilibrium towards the active kinase conformation. Our atomic investigation of mutant KIT receptor which emphasized on EAL region provided a better insight into the understanding of Sunitinib resistance mechanism of KIT receptor and would help to discover new therapeutics for KIT-based resistant tumor cells in GIST therapy.

Studies on Adaptability of Binding Residues Flap Region of TMC-114 Resistance HIV-1 Protease Mutants

Abstract Drug resistant mutations have severely restricted the success of HIV therapy. These mutations frequently involve the aspartic protease encoded by the virus. Knowledge of the molecular mechanisms underlying the conformational changes of HIV-1 protease mutants may be useful in developing more ffective and longer lasting treatment regimes. The flap regions of the protease are the target of a particular type of mutations occurring far from the active site, which are able to produce significant resistance against the anti-HIV drug TMC-114. We provide insight into the molecular basis of TMC-114 resistance major flap mutations (I50V and I54M) in HIV-1 protease. It reports the shape complementarity and receptor-ligand interaction analysis supported by unrestrained all-atom molecular dynamics simulations of wild and major flap mutants of HIV-1 protease that sample large conformational changes of the flaps and active site binding residues. Both resistant flap mutants showed less atomic interaction toward TMC-114 and more structural deviation compared to wild HIV-protease. It is due to increasing flexibility at TMC-114 binding cavity and deviation of binding residues in 3-D space. Distortion in binding cavity and deviation in binding residues are the result of alteration in hydrogen bonding. Flap region also exhibited similar behaviour due to changes in number of hydrogen bonds during simulations.

Use of Long Term Molecular Dynamics Simulation in Predicting Cancer Associated SNPs

Computational prediction of cancer associated SNPs from the large pool of SNP dataset is now being used as a tool for detecting the probable oncogenes, which are further examined in the wet lab experiments. The lack in prediction accuracy has been a major hurdle in relying on the computational results obtained by implementing multiple tools, platforms and algorithms for cancer associated SNP prediction. Our result obtained from the initial computational compilations suggests the strong chance of Aurora-A G325W mutation (rs11539196) to cause hepatocellular carcinoma. The implementation of molecular dynamics simulation (MDS) approaches has significantly aided in raising the prediction accuracy of these results, but measuring the difference in the convergence time of mutant protein structures has been a challenging task while setting the simulation timescale. The convergence time of most of the protein structures may vary from 10 ns to 100 ns or more, depending upon its size. Thus, in this work we have implemented 200 ns of MDS to aid the final results obtained from computational SNP prediction technique. The MDS results have significantly explained the atomic alteration related with the mutant protein and are useful in elaborating the change in structural conformations coupled with the computationally predicted cancer associated mutation. With further advancements in the computational techniques, it will become much easier to predict such mutations with higher accuracy level.

Structural basis for the resilience of Darunavir (TMC114) resistance major flap mutations of HIV-1 protease

To understand the origin of the apparent low sensitivity to mutations exhibited by Darunavir, the binding energetics of this inhibitor to the HIV-1 protease was studied. Our research indicates that the observed effectiveness of Darunavir against the wild type HIV-1 protease is due to an extremely high affinity towards the wild-type and a relatively mild effect to the I50V and I54M mutations is due to low affinity towards the inhibitor. Good affinity of Darunavir accounts for the additive effects of well accommodation at binding site, good ligand-receptor electrostatic and van der waals energy while, the low susceptibility to I50V and I54M can be rationalized in terms of flexibility in the binding site residues that do not permit drug accommodation to the binding site distortions created by the mutation. The major flap mutations I50V and I54M lower the binding affinity of Darunavir by altering the position of binding site residues in 3D space. It decreases the electrostatic and van der waals interaction energy and further reduction in total receptor-ligand interaction energy. The results summarized in this paper emphasize the importance of shape complementarity and protein flexibility analysis of binding residual interactions in drug design. These data together with an interaction energy and flexibility analysis have established rigorous guidelines for the design of new and more powerful inhibitors. The principles learned from the HIV-1 protease can be applied to other design problems.

Computational screening and molecular dynamics simulation of disease associated nsSNPs in CENP-E.

Aneuploidy and chromosomal instability (CIN) are hallmarks of most solid tumors. Mutations in centroemere proteins have been observed in promoting aneuploidy and tumorigenesis. Recent studies reported that Centromere-associated protein-E (CENP-E) is involved in inducing cancers. In this study we investigated the pathogenic effect of 132 nsSNPs reported in CENP-E using computational platform. Y63H point mutation found to be associated with cancer using SIFT, Polyphen, PhD-SNP, MutPred, CanPredict and Dr. Cancer tools. Further we investigated the binding affinity of ATP molecule to the CENP-E motor domain. Complementarity scores obtained from docking studies showed significant loss in ATP binding affinity of mutant structure. Molecular dynamics simulation was carried to examine the structural consequences of Y63H mutation. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (R(g)), solvent accessibility surface area (SASA), energy value, hydrogen bond (NH Bond), eigenvector projection, trace of covariance matrix and atom density analysis results showed notable loss in stability for mutant structure. Y63H mutation was also shown to disrupt the native conformation of ATP binding region in CENP-E motor domain. Docking studies for remaining 18 mutations at 63rd residue position as well as other two computationally predicted disease associated mutations S22L and P69S were also carried to investigate their affect on ATP binding affinity of CENP-E motor domain. Our study provided a promising computational methodology to study the tumorigenic consequences of nsSNPs that have not been characterized and clear clue to the wet lab scientist.

Computational investigation of pathogenic nsSNPs in CEP63 protein

Centrosomes are central regulators of mitosis that are often amplified in cancer cells. Centrosomal protein of 63kDa (CEP63) is a centrosomal protein that has an effective role in mitotic spindle assembly and cell cycle regulation. Genetic alterations in CEP63 coding gene has been widely studied for inducing aneuploidy and solid tumors in humans. The nonsynonymous single nucleotide polymorphisms (nsSNPs) are a genetic variant resulting in amino acid substitution and are reported in a wide range of human diseases. Here we report one new SNP (rs112926188) in a CEP63 coding region that can potentially disrupt the structure and basic functionality of a CEP63 protein. We used extensive functional and structural level analyses of an available SNP in a CEP63 coding gene. Furthermore the disease-association analysis was carried out to examine the possible pathogenic variant among the available dataset. To understand atomic arrangement in 3D space, native and pathogenic mutant structures were modeled. Molecular dynamics simulations were performed to understand structural consequences of prioritized deleterious mutation. Our analysis showed that rs112926188 allele substituting proline at the 61st residue position (L61P) produced more flexibility in 3D space. Moreover the flexible nature of mutant L61P was validated by a hydrogen bond network. This nature of mutant L61P CEP63 may restrict the recruitment of essential centrosomal proteins to their respective location and may play an active role in inducing aneuploidy.

A novel computational and structural analysis of nsSNPs in CFTR gene

Single Nucleotide Polymorphisms (SNPs) are being intensively studied to understand the biological basis of complex traits and diseases. The Genetics of human phenotype variation could be understood by knowing the functions of SNPs. In this study using computational methods, we analyzed the genetic variations that can alter the expression and function of the CFTR gene responsible candidate for causing cystic fibrosis. We applied an evolutionary perspective to screen the SNPs using a sequence homology-based SIFT tool, which suggested that 17 nsSNPs (44%) were found to be deleterious. The structure-based approach PolyPhen server suggested that 26 nsSNPS (66%) may disrupt protein function and structure. The PupaSuite tool predicted the phenotypic effect of SNPs on the structure and function of the affected protein. Structure analysis was carried out with the major mutation that occurred in the native protein coded by CFTR gene, and which is at amino acid position F508C for nsSNP with id (rs1800093). The amino acid residues in the native and mutant modeled protein were further analyzed for solvent accessibility, secondary structure and stabilizing residues to check the stability of the proteins. The SNPs were further subjected to iHAP analysis to identify htSNPs, and we report potential candidates for future studies on CFTR mutations.

Mutational analysis of TYR gene and its structural consequences in OCA1A.

Oculocutaneous albinism type 1A (OCA1A) is the most severe form of albinism characterized by a complete lack of melanin production throughout life and is caused by mutations in the TYR gene. TYR gene codes tyrosinase protein to its relation with melanin formation by knowing the function of these SNPs. Based on the computational approaches, we have analyzed the genetic variations that could change the functional behaviour by altering the structural arrangement in TYR protein which is responsible for OCA1A. Consequences of mutation on TYR structure were observed by analyzing the flexibility behaviour of native and mutant tyrosinase protein. Mutations T373K, N371Y, M370T and P313R were suggested as high deleterious effect on TYR protein and it is responsible for OCA1A which were also endorsed with previous in vivo experimental studies. Based on the quantitative assessment and flexibility analysis of OCA1A variants, T373K showed the most deleterious effect. Our analysis determines that certain mutations can affect the dynamic properties of protein and can lead to disease conditions. This study provides a significant insight into the underlying molecular mechanism involved in albinism associated with OCA1A.

Mutational Analysis of Oculocutaneous Albinism: A Compact Review

Oculocutaneous albinism (OCA) is an autosomal recessive disorder caused by either complete lack of or a reduction of melanin biosynthesis in the melanocytes. The OCA1A is the most severe type with a complete lack of melanin production throughout life, while the milder forms OCA1B, OCA2, OCA3, and OCA4 show some pigment accumulation over time. Mutations in TYR, OCA2, TYRP1, and SLC45A2 are mainly responsible for causing oculocutaneous albinism. Recently, two new genes SLC24A5 and C10orf11 are identified that are responsible to cause OCA6 and OCA7, respectively. Also a locus has been mapped to the human chromosome 4q24 region which is responsible for genetic cause of OCA5. In this paper, we summarized the clinical and molecular features of OCA genes. Further, we reviewed the screening of pathological mutations of OCA genes and its molecular mechanism of the protein upon mutation by in silico approach. We also reviewed TYR (T373K, N371Y, M370T, and P313R), OCA2 (R305W), TYRP1 (R326H and R356Q) mutations and their structural consequences at molecular level. It is observed that the pathological genetic mutations and their structural and functional significance of OCA genes will aid in development of personalized medicine for albinism patients.