Finding Pathogenic Nssnp's and Their Structural Effect on COPS2 Using Molecular Dynamic Approach
Total Page:16
File Type:pdf, Size:1020Kb
bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.333252; this version posted October 12, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Finding Pathogenic nsSNP's and their structural effect on COPS2 using Molecular Dynamic Approach Ashish Malik1, Kajal Pande1, Abhishek Kumar1, Alekhya Vemula1, Madhuri R1 Vivek Chandramohan2* 1) BiSEP, Bioinformatics and Big data analytics, Department of Biotechnology, Siddaganga Institute of Technology, Tumkur. 2) Department of Biotechnology, Siddaganga Institute of Technology, Tumkur. Corresponding author: Mr. Vivek Chandramohan (M.Sc) Assistant Professor, Backup Coordinator for Biotechnology Skill Enhancement Programme (BiSEP) Department of Biotechnology, Siddaganga Institute of Technology, Tumkur - 572103 bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.333252; this version posted October 12, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Abstract COP9 Signalosome Subunit 2 is a highly conserved multiprotein complex which is involved in the cellular process and developmental process. It is one of the essential components in the COP9 Signalosome Complex (CSN). It is also involved in neuronal differentiation interacting with NIF3L1. The gene involved in neuronal differentiation is negatively regulated due to the transcription co-repressor interaction of NIF3L1 with COPS2. In the present study, we have evaluated the outcome for 90 non- synonymous single nucleotide polymorphisms (nsSNP’s) in COPS2 gene through computational tools. After the analysis, 4 SNP’s (S120C, N144S, Y159H, R173C) were found to be deleterious. The native and mutated structures were prepared using discovery studio and docked to check the interactions with NIF3L1.On the basis of ZDOCK score the top 3 mutations (N144S, Y159H, R173C) were screened out. Further to analyze the effect of amino acid substitution on the molecular structure of protein Molecular Dynamics simulation was carried out. Analysis based on RMSD, RMSF, RG, H-bond showed a significant deviation in the graph, which demonstrated conformation change and instability compared to the wild structure. As it is known mutations in COPS2 gene can disrupt the normal activity of the CSN2 protein which may cause neuronal differentiation. Our results showed N144S, Y159H and R173C mutations are to be more pathogenic and may cause disease KEYWORDS: Homology modeling, MD simulation, protein-protein docking, single nucleotide polymorphism bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.333252; this version posted October 12, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Introduction COP9 Signalosome Subunit 2 (COPS2) is one of the important component of COP9 signalosome complex (CSN) that exists in eukaryotes. There are many synonyms that are in use for COPS2 like, TRIP15, SGN2 and CSN2. CSN is a conserved multi-protein complex (CSN1 to CSN8) that has seemingly evolved in parallel with the ubiquitin- proteasome system [1]. In a eukaryotic organism, CSN complex is associated to be as an essential protein complex that is involved in multiple cellular and developmental processes. Any deletion of CSN subunits can cause serious changes in the gene articulation profile and can cause early lethality of the multicellular organism [2–4]. The CSN complex was initially discovered as a eliminator of photomorphogenesis in Arabidopsis thaliana [1]. The eight mammalian complex subunits are classified according to the molecular weight S1 (57 kDa) to S8 (22 kDa) [5]. COP9 signalosome complex was considered as contaminating component in 26S proteasome preparation from erythrocytes by Seegeret al. and they later separated the two complexes [6]. The complex was cleansed biochemically from pig spleen by [5]. The CSN complex is a crucial regulator of ubiquitin (Ubl) conjugation pathway by mediating the deneddylation of the cullin subunits [7]. The subunits of CSN shows notable sequence homologies to elements of the lid subcomplex of the 26S proteasome which is formed by the assembly of 19S proteasome regulatory complex with the 20S catalytic core [8] and the translation initiation complex eIF3[6]. In eukaryotes, the eIF3 complex is involved in protein synthesis and ribosome biogenesis and it serves as a binding site for assembling other eIFsduring translation initiation [9]. Initially, Trip15/CSN2 was recognized as a thyroid hormone receptor (TR)-interacting protein and acts as a transcriptional corepressor [5]. In yeast two-hybrid screen, Trip15 is reported to interact with the thyroid hormone receptor and the retinoic acid X receptor in a ligand- dependent manner and the N-terminal region of subunit 2 (S2) is identical to Trip15[10]. Even though the mutual interaction between Trip15-nuclear receptor complex and CSN complex is obscure, indicating that the CSN complex assists in several signal transduction pathways [11]. COPS2 is a 51,597 Da protein with a sequence length of 443 amino acids. It also dissipates the interaction between OCT4 and CDK1 by segregating OCT4 and forming a COPS2/CDK1 complex thus hindering the inhibitory effect of OCT4 on CDK1 activation[12]. COPS2 (CSN2) shows transcriptional corepressor activity. The protein-protein interaction occurring inside the CSN complex (21-24) is stabilized by the proteasome COP9 complex initiation factor domain 3 in the C-terminal region of the Trip15/CSN whereas the N-terminal region of CSN2 has been reported to be adequate for the effector function of Trip15/CSN2. Trip15/CSN2 is associated with its binding partner through N terminal (Akiyama et al. 2003). The CS2 domain region (1-275) mediates interaction with NIF3L1, which is a highly conserved gene from bacteria to human. Although it is known that the mouse NIF3L1 gene encodes for a cytoplasmic protein (376 amino acids) and the expression of the protein was detected through mouse embryonic development, still the real-biological function of Nif3l1 is still unknown. In order to study the association of Trip15/CSN2 with Nif3l1, a wet lab study was carried out by Akiyama and co to exhibit that Nif3l1 protein could be translocated into nuclei by Trip15/CSN2 and it is involved in neural differentiation through down-regulation of Oct-¾ transcript which suppresses neurogenic genes. By the result of the study carried out by Akiyama and co, it seems likely that Nif3l1 and Trip15/CSN2 are associated and play an important role in neural differentiation [11]. NIF3L1 is engaged in retinoic acid primed neural differentiation of P19 embryonic carcinoma cells through cooperation with CSN2 [11]. In the complex COPS2 bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.333252; this version posted October 12, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. interacts with COPS1, COPS4, COPS5, COPS6, COPS7A or COPS7B directly[8, 13]. Notably, it also interacts with the ligand binding region of thyroid receptor, CUL1 and CUL2 [8]. Single nucleotide polymorphism is known as the most common type of variants occurring in human genome due to single nucleotide substitution [14, 15]. Non- synonymous SNPs can cause 50% of the disease that is related to inheritance [16]. SNPs can also disrupt the normal function of the protein and denature the structure by changing the stability, folding pattern and ligand binding site [16]. The amino acid changes at key sites within a protein may result in series of conformation changes, including the disintegration of salt bridges, alteration of interaction network, disruption in hydrogen bonds, which in turn may disturb the energy landscape [17]. Protein structure modeling methods have been extensively used for predicting the effects of disease-causing mutation on protein stability and protein-protein interaction [18]. The aim of the present study was to identify the most deleterious nsSNPs within COPS2 gene and their effect on protein structure. Protein-protein docking was done in order to determine the effect of amino acid substitution on the structure of COPS2 which will affect the interaction with Nif3l1. To check the interaction of COPS2 with NIF3L1, protein-protein docking was done. The mutations which are identified to have a positive interrelationship with pathogenicity and have high ZDOCK score were further investigated using MD simulations to determine the changes on its time-dependent physiological affinities and biochemical pathway alterations (Sinha et al. 2018). We utilized RMSD, RMSF, H-Bond, Gyration, in order to examine the motion trajectory and atomic interactions of native and mutant protein structures. Alteration in the protein structure can lead to a change in the binding site of protein which may affect the interacting proteins of COPS2 that may cause disease. bioRxiv preprint doi: https://doi.org/10.1101/2020.10.12.333252; this version posted October 12, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. All rights reserved. No reuse allowed without permission. Materials and Methods SNPs Data Collection The data of human COPS2 gene were obtained from the NCBI database of dbSNP (www.ncbi.nlm.nih.gov/snp),