Structural implications of amyloidogenic rare variants Ser282Leu and Gln356Arg identified in h-BRCA1.

Preliminary studies have shown BRCA1 (170-1600) residues to be intrinsically disordered with unknown structural details. However, thousands of clinically reported variants have been identified in this central region of BRCA1. Therefore, we aimed to characterize h-BRCA1(260-553) to assess the structural basis for pathogenicity of two rare missense variants Ser282Leu, Gln356Arg identified from the Indian and Russian populations respectively. Small-angle X-ray scattering analysis revealed WT scores Rg -32 Å, Dmax -93 Å, and Rflex-51% which are partially disordered, whereas Ser282Leu variant displayed a higher degree of disorderedness and Gln356Arg was observed to be aggregated. WT protein also possesses an inherent propensity to undergo a disorder-to-order transition in the presence of cruciform DNA and 2,2,2-Trifluoroethanol (TFE). An increased alpha-helical pattern was observed with increasing concentration of TFE for the Gln356Arg mutant whereas Ser282Leu mutant showed significant differences only at the highest TFE concentration. Furthermore, higher thermal shift was observed for WT-DNA complex compared to the Gln356Arg and Ser282Leu protein-DNA complex. Moreover, mature amyloid-like fibrils were observed with 30 µM thioflavin T (ThT) at 37°C for Ser282Leu and Gln356Arg proteins while the WT protein exists in a protofibril state as observed by TEM. Gln356Arg formed higher-order aggregates with amyloidogenesis over time as monitored by ThT fluorescence. In addition, computational analyses confirmed larger conformational fluctuations for Ser282Leu and Gln356Arg mutants than for the WT. The global structural alterations caused by these variants provide a mechanistic approach for further classification of the variants of uncertain clinical significance in BRCA1 into amyloidogenic variants which may have a significant role in disease pathogenesis.

In silico and structure-based evaluation of deleterious mutations identified in human Chk1, Chk2, and Wee1 protein kinase.

Checkpoint kinases Chk1, Chk2, Wee1 are playing a key role in DNA damage response and genomic integrity. Cancer-associated mutations identified in human Chk1, Chk2, and Wee1 were retrieved to understand the function associated with the mutation and also alterations in the folding pattern. Therefore, an attempt has been made to identify deleterious effect of variants using in silico and structure-based approach. Variants of uncertain significance for Chk1, Chk2, and Wee1 were retrieved from different databases and four prediction servers were employed to predict pathogenicity of mutations. Further, Interpro, I-Mutant 3.0, Consurf, TM-align, and have (y)our protein explained were used for comprehensive study of the deleterious effects of variants. The sequences of Chk1, Chk2, and Wee1 were analyzed using Clustal Omega, and the three-dimensional structures of the proteins were aligned using TM-align. The molecular dynamics simulations were performed to explore the differences in folding pattern between Chk1, Chk2, Wee1 wild-type, and mutant protein and also to evaluate the structural integrity. Thirty-six variants in Chk1, 250 Variants in Chk2, and 29 in Wee1 were categorized as pathogenic using in silico prediction tools. Furthermore, 25 mutations in Chk1, 189 in Chk2, and 14 in Wee1 were highly conserved, possessing deleterious effect and also influencing the protein structure and function. These identified mutations may provide underlying genetic intricacies to serve as potential targets for therapeutic inventions and clinical management.

Crystal structure of clinically reported mutations Gly656Arg, Gly656Glu and Asp751His identified in the kinase domain of EphA7.

Erythropoietin producing hepatocellular (Eph) forms the largest family of receptor tyrosine kinases (RTK). As a family, Eph regulates physiological events such as cell-cell interaction, cell migration, and adhesion. The Kinase domain is the catalytic core of the Eph receptor and is highly conserved sequentially. EphA7 has been recently regarded as a cancer driver gene and comprises several clinically important mutations. Three of the EphA7 mutations Gly656Glu, Gly656Arg, and Asp751His, present in the kinase domain, are predicted to be highly pathogenic. Furthermore, Gly656Glu and Gly656Arg are reported to be hotspot mutations. Considering the importance of mutations, crystals structure of EphA7 Gly656Glu, Gly656Arg, and Asp751His mutants has been explored. Changes in folding pattern and intramolecular interactions were observed in mutant structures. Secondary structural changes were observed in the hinge region of EphA7 Gly656Arg and Asp751His structure, affecting the transition of kinase domain between open and closed conformations. EphA7 Asp751His mutant structure shows a distorted nucleotide-binding groove. Differences were observed in hydrogen bonding and hydrophobic interactions between the catalytic and highly conserved DFG motif in the EphA7 mutants, which may influence the catalytic activity of kinase domain.

Genetic ablation of pregnancy zone protein promotes breast cancer progression by activating TGF-ß/SMAD signaling.

PURPOSE: Pregnancy zone protein (PZP) is best known as protease inhibitor and its concentration in human blood plasma increases dramatically during pregnancy. Recent investigation revealed a role of PZP inactivating germ-line mutation in breast cancer predisposition, and therefore we designed a study to evaluate functional involvement of this protein in tumor pathogenesis. METHODS: PZP knockout cells were generated utilizing the CRISPR-Cas9 approach in MCF7 and T47D (breast cancer) cell lines, and colony formation, cell proliferation, and migration assays carried out. TGF-ß and SMAD expression studies were performed using qRT-PCR and Western blot. PZP expression in tumor vs normal tissue was compared using meta-analyses of data records of breast cancer patients (n = 1211) included in the TCGA consortium registry as well as in independent cohorts of hormone receptor-positive (n = 118) and triple-negative breast cancer (TNBC) patients (n = 116). RESULTS: We demonstrated that genetic ablation of PZP efficiently inhibits tamoxifen-induced apoptosis and enhances cell proliferation, migration, and colony-forming capacity. We found a significant increase in survival fraction of CRISPR/Cas9-mediated PZP knockout clones compared to wild-type counterpart after tamoxifen treatment (p < 0.05). The PZP knockout significantly promoted breast cancer cell migration (p < 0.01) in vitro. We observed high expression of TGF-ß2 ligand, TGF-ß- receptor 2, and upregulation of phosphorylated regulatory-SMADs (pSMAD2 and pSMAD3) activating the pro-survival function of TGF-ß/SMAD signaling in PZP knockout clones. Meta-analyses of data records of breast cancer patients indicated that low PZP expression is associated with poor overall survival at 6 years (51.7% vs 62.9% in low vs high expressers, respectively; p = 0.026). We also observed a significantly lower PZP mRNA expression in TNBC as compared with hormone receptor-positive tumors (p = 0.019). CONCLUSION: Taken together, our results suggest that genetic ablation of PZP results in tumor progression and low expression of PZP is associated with poor survival of breast cancer patients.

Exome sequencing study of Russian breast cancer patients suggests a predisposing role for USP39.

PURPOSE: Germline variants in known breast cancer (BC) predisposing genes explain less than half of hereditary BC cases. This study aimed to identify missing genetic determinants of BC. METHODS: Whole exome sequencing (WES) of lymphocyte DNA was performed for 49 Russian patients with clinical signs of genetic BC predisposition, who lacked Slavic founder mutations in BRCA1, BRCA2, CHEK2, and NBS1 genes. RESULTS: Bioinformatic analysis of WES data was allowed to compile a list of 229 candidate mutations. 79 of these mutations were subjected to a three-stage case-control analysis. The initial two stages, which involved up to 797 high-risk BC patients, 1504 consecutive BC cases, and 1081 healthy women, indicated a potentially BC-predisposing role for 6 candidates, i.e., USP39 c.*208G > C, PZP p.Arg680Ter, LEPREL1 p.Pro636Ser, SLIT3 p.Arg154Cys, CREB3 p.Lys157Glu, and ING1 p.Pro319Leu. USP39 c.*208G > C was strongly associated with triple-negative breast tumors (p = 0.0001). In the third replication stage, we genotyped the truncating variant of PZP (rs145240281) and the potential splice variant of USP39 (rs112653307) in three independent cohorts of Russian, Byelorussian, and German ancestry, comprising a total of 3216 cases and 2525 controls. The data obtained for USP39 rs112653307 supported the association identified in the initial stages (the combined OR 1.72, p = 0.035). CONCLUSIONS: This study suggests the role of a rare splicing variant in BC susceptibility. USP39 encodes an ubiquitin-specific peptidase that regulates cancer-relevant tumor suppressors including CHEK2. Further epidemiological and functional studies involving these gene variants are warranted.

Influence of ABCB-1, ERCC-1 and ERCC-2 gene polymorphisms on response to capecitabine and oxaliplatin (CAPOX) treatment in colorectal cancer (CRC) patients of South India.

WHAT IS KNOWN AND OBJECTIVE: High interindividual response variability was reported with capecitabine and oxaliplatin (CAPOX) regimen in colorectal cancer (CRC). The single nucleotide polymorphisms (SNPs) of the genes related to drug efflux transport (ABCB1) and DNA repair (ERCC) could result in altered tumour response. Hence, this study was designed to assess the influence of ABCB1, ERCC-1 and ERCC-2 gene polymorphisms on tumour response to CAPOX treatment in CRC patients of South Indian origin. PATIENTS AND METHODS: A total of 145 newly diagnosed CRC patients were included in the final analysis. Response to CAPOX treatment in the adjuvant setting was assessed in terms of disease-free survival rate (DFSR) and overall survival rate (OSR) at 3 years, whereas in the palliative setting, the response was assessed as progression-free survival rate (PFSR) and OSR at 3 years. Five millilitres of the venous blood sample was collected from each patient for genomic DNA extraction by the manual phenol-chloroform method. Genotyping and allelic discrimination analysis were done using real-time PCR (RT-PCR). RESULTS AND DISCUSSION: With ABCB1 gene polymorphism rs1045642 (A > G), patients with AG/GG genotype showed better DFSR [P value = .02, OR = 2 (CI: 1.5-3)] and PFSR [P value = .02, OR = 1.6 (CI: 1.1-2.5)] when compared to AA genotype in the adjuvant and palliative settings, respectively. Similarly with rs1128503 (A > G) polymorphism, patients with AG/GG genotype were found to have better DFSR [P value = .02, OR = 1.9 (CI: 1.3-3)] and PFSR [P value = .01, OR = 2 (CI: 1.1-3.7)] when compared to AA genotype. However, we did not find any association between CAPOX response and ABCB1 gene polymorphisms in a binary logistic regression when non-genetic predictors were considered for analysis. We did not find any association with ERCC1 (rs11615 A > G) and ERCC2 (rs13181 T > G) gene polymorphisms with respect to CAPOX response in either of the treatment settings. WHAT IS NEW AND CONCLUSION: The response to CAPOX treatment was found to be influenced by the ABCB1 gene variants (rs1128503 and rs1045642), thereby strengthening their predictive role. No association was found between ERCC1 (rs11615 A > G), ERCC2 (rs13181 T > G) gene polymorphisms and tumour response to CAPOX treatment in CRC patients of South Indian origin.

Identification and characterization of GRIP domain Golgin PpImh1 from Pichia pastoris.

Budding yeast Pichia pastoris has highly advanced secretory pathways resembling mammalian systems, an advantage that makes it a suitable model system to study vesicular trafficking. Golgins are large Golgi-resident proteins, primarily reported to play role in cargo vesicle capture, but details of such mechanisms are yet to be deciphered. Golgins that localize to the Golgi via their GRIP domain, a C-terminal Golgi anchoring domain, are known as GRIP domain Golgins. In this present study, we have identified and functionally characterized a homologue of one such GRIP domain Golgin protein, Imh1, from the budding yeast P. pastoris. We have demonstrated that the GRIP domain present at the C-terminal of P. pastoris Imh1 (PpImh1) functions as its Golgi-targeting sequence. Using a combination of yeast two-hybrid analysis, dynamic light scattering and electron microscopy, we have shown that PpImh1 can self-associate and form a homodimer. Analysis of purified recombinant PpImh1 by CD spectroscopy indicates the presence of an 85% α-helical structure, a characteristic of high-content α-helical coiled-coil sequences normally present in other Golgin family proteins. Two-hybrid analysis indicated self-interaction between C-terminal fragments, yet N-terminal fragments do not mediate any such form of self-interaction, suggesting that PpImh1 may form a parallel dimer. Electron microscopy data indicates that PpImh1 forms extended rod-like homo-dimeric molecules with splayed N-terminal end which can act as a tether for capturing vesicles. Our study provides the first evidence in support of the dimeric Y-shaped structure for any Golgin in the budding yeast.

Fundamental niche prediction of the pathogenic yeasts Cryptococcus neoformans and Cryptococcus gattii in Europe.

Fundamental niche prediction of Cryptococcus neoformans and Cryptococcus gattii in Europe is an important tool to understand where these pathogenic yeasts have a high probability to survive in the environment and therefore to identify the areas with high risk of infection. In this study, occurrence data for C. neoformans and C. gattii were compared by MaxEnt software with several bioclimatic conditions as well as with soil characteristics and land use. The results showed that C. gattii distribution can be predicted with high probability along the Mediterranean coast. The analysis of variables showed that its distribution is limited by low temperatures during the coldest season, and by heavy precipitations in the driest season. C. neoformans var. grubii is able to colonize the same areas of C. gattii but is more tolerant to cold winter temperatures and summer precipitations. In contrast, the C. neoformans var. neoformans map was completely different. The best conditions for its survival were displayed in sub-continental areas and not along the Mediterranean coasts. In conclusion, we produced for the first time detailed prediction maps of the species and varieties of the C. neoformans and C. gattii species complex in Europe and Mediterranean area.

Multilocus sequence typing analysis reveals that Cryptococcus neoformans var. neoformans is a recombinant population.

Cryptococcus neoformans var. neoformans (serotype D) represents about 30% of the clinical isolates in Europe and is present less frequently in the other continents. It is the prevalent etiological agent in primary cutaneous cryptococcosis as well as in cryptococcal skin lesions of disseminated cryptococcosis. Very little is known about the genotypic diversity of this Cryptococcus subtype. The aim of this study was to investigate the genotypic diversity among a set of clinical and environmental C. neoformans var. neoformans isolates and to evaluate the relationship between genotypes, geographical origin and clinical manifestations. A total of 83 globally collected C. neoformans var. neoformans isolates from Italy, Germany, France, Belgium, Denmark, Greece, Turkey, Thailand, Japan, Colombia, and the USA, recovered from different sources (primary and secondary cutaneous cryptococcosis, disseminated cryptococcosis, the environment, and animals), were included in the study. All isolates were confirmed to belong to genotype VNIV by molecular typing and they were further investigated by MLST analysis. Maximum likelihood phylogenetic as well as network analysis strongly suggested the existence of a recombinant rather than a clonal population structure. Geographical origin and source of isolation were not correlated with a specific MLST genotype. The comparison with a set of outgroup C. neoformans var. grubii isolates provided clear evidence that the two varieties have different population structures.

Role of MERIT40 in stabilization of BRCA1 complex: a protein-protein interaction study.

MERIT40 is a novel associate of the BRCA1-complex, thus play an essential role in DNA damage repair mechanism. It is the least implicit protein and its structural and functional aspects of regulating the stability of BRCA1-MERIT40 complex remain equivocal. Analysis of protein-protein interactions between BRCA1 and its cellular binding partners like ABRAXAS, RAP80 and MERIT40 would help to understand the role of protein complex integrity in DNA repair mechanism. The recombinant proteins were purified and their structural aspects were elucidated by spectroscopic methods. Interaction analysis was carried out to determine binding partners of MERIT40. MERIT40 showed interaction with bridging molecule, called ABRAXAS, thus generate a scaffold among various members which further stabilizes the entire complex. It acts as an adapter molecule by interacting with BRCA1-BRCT in non-phosphorylation dependent manner. The feature enlighten on structural and interaction profile of BRCA1-complex member to elucidate their role in complex stability and DNA repair process.

Tetrameric ZBRK1 DNA binding domain has affinity towards cognate DNA in absence of zinc ions.

Zinc finger transcription regulatory proteins play crucial roles in cell-cycle regulation, DNA damage response and tumor genesis. Human ZBRK1 is a zinc-finger transcription repressor protein, which recognizes double helical DNA containing consensus sequences of 5'GGGXXXCAGXXXTTT3'. In the present study, we have purified recombinant DNA binding domain of ZBRK1, and studied binding with zinc ions and DNA, using biophysical techniques. The elution profile of the purified protein suggests that this ZBRK1 forms a homotetramer in solution. Dissociation and pull down assays also suggest that this domain forms a higher order oligomer. The ZBRK1-DNA binding domain acquires higher stability in the presence of zinc ions and DNA. The secondary structure of the ZBRK1-DNA complex is found to be significantly altered from the standard B-DNA conformation.

Mutations at proximal cysteine residues in PML impair ATO binding by destabilizing the RBCC domain.

Acute promyelocytic leukemia (APL) is characterized by the fusion gene promyelocytic leukemia-retinoic acid receptor-alpha (PML-RARA) and is conventionally treated with arsenic trioxide (ATO). ATO binds directly to the RING finger, B-box, coiled-coil (RBCC) domain of PML and initiates degradation of the fusion oncoprotein PML-RARA. However, the mutational hotspot at C212-S220 disrupts ATO binding, leading to drug resistance in APL. Therefore, structural consequences of these point mutations in PML that remain uncertain require comprehensive analysis. In this study, we investigated the structure-based ensemble properties of the promyelocytic leukemia-RING-B-box-coiled-coil (PML-RBCC) domains and ATO-resistant mutations. Oligomeric studies reveal that PML-RBCC wild-type and mutants C212R, S214L, A216T, L217F, and S220G predominantly form tetramers, whereas mutants C213R, A216V, L218P, and D219H tend to form dimers. The stability of the dimeric mutants was lower, exhibiting a melting temperature (Tm) reduction of 30 °C compared with the tetrameric mutants and wild-type PML protein. Furthermore, the exposed surface of the C213R mutation rendered it more prone to protease digestion than that of the C212R mutation. The spectroscopic analysis highlighted ATO-induced structural alterations in S214L, A216V, and D219H mutants, in contrast to C213R, L217F, and L218P mutations. Moreover, the computational analysis revealed that the ATO-resistant mutations C213R, A216V, L217F, and L218P caused changes in the size, shape, and flexibility of the PML-RBCC wild-type protein. The mutations C213R, A216V, L217F, and L218P destabilize the wild-type protein structure due to the adaptation of distinct conformational changes. In addition, these mutations disrupt several hydrogen bonds, including interactions involving C212, C213, and C215, which are essential for ATO binding. The local and global structural features induced by these mutations provide mechanistic insight into ATO resistance and APL pathogenesis.

Variant of uncertain significance Arg866Cys enhances disorderedness of h-BRCA1 (759-1064) region.

High structural flexibility has been reported in the central region of BRCA1, which hinders the structural and functional evaluations of mutations identified in the domain. Additionally, the need to categorize variants of unknown significance (VUS) has increased due to the growth in the number of variants reported in clinical settings. Therefore, unraveling the disease-causing mechanism of VUS identified in different functional domains of BRCA1 is still challenging. The current study uses a multidisciplinary approach to assess the structural impact of BRCA1 Arg866Cys mutation discovered in the central domain of BRCA1. The structural alterations have been characterized using Circular-Dichroism spectroscopy, nano-DSF, and molecular-dynamics simulations. BRCA1 Arg866Cys mutant demonstrated more flexibility and lesser affinity to DNA than the wild-type protein. The BRCA1(759-1064) wild-type protein was shown to be a ßII-rich protein with an induced D-O transition in the presence of DNA and 2,2,2-Trifluoroethanol (TFE). The protein's alpha-helical composition did not significantly change in the presence of TFE, besides an increase in ß-turns and loops. Under Transmission Electron Microscopes (TEM), amyloid-like fibrils structure was detected for Arg866Cys mutant whereas the wild-type protein showed amorphous aggregates. An increased ThT fluorescence indicated ß-rich composition and aggregation-prone behaviour for BRCA1 wild-type protein, while the fluorescence intensity was significantly quenched in the Arg866Cys mutant. Furthermore, increased conformational flexibility in the Arg866Cys variant was observed by principal component analysis. This work aims to comprehend the inherently disordered region of BRCA1 as well as the impact of missense mutations on folding patterns and binding to DNA for functional aspects.

Preliminary crystallographic studies of BRCA1 BRCT-ABRAXAS complex.

The BRCA1 holoenzyme complex plays an important role in DNA damage repair. ABRAXAS is a newly discovered component of this complex and its C-terminal region directly binds to the BRCA1 BRCT domain. Single crystals of the BRCA1 BRCT-ABRAXAS complex grown by co-crystallization belonged to space group P4(1)2(1)2, with unit-cell parameters a = b = 187.18, c = 85.31 Å. Diffraction data were collected on the BM-14 beamline at the ESRF. Molecular-replacement calculations using Phaser led to three molecules in the asymmetric unit and a high solvent content of 76%.

In silico and structure-based assessment to classify VUS identified in the α-helical domain of BRCA2.

Breast cancer type 2 susceptibility (BRCA2) protein plays a crucial role in DNA double-strand breaks repair mechanism by homologous recombination. Pathogenic mutations in the BRCA2 gene confer an increased risk of hereditary breast and ovarian cancer (HBOC). Different missense mutations are identified from a larger cohort of patient populations in the BRCA2. However, most missense mutations are classified as 'Variants of Uncertain Significance' (VUS) due to a lack of data from structural, functional, and clinical assessments. Therefore, this study focused on assessing VUS identified in the α-helical domain of h-BRCA2 using different in silico tools and structure-based molecular dynamics simulation. A total of 286 identified VUS were evaluated using Align-GVGD, PROVEAN and PANTHER servers and 18 variants were predicted to be pathogenic. Further, out of 18 variants analyzed using the ConSurf server, 16 variants were found to be evolutionary conserved. These 16 conserved variants were submitted to PremPS and Dynamut server to assess the effect of the mutation at the protein structure level; 12 mutations were predicted to have a destabilizing effect on the native protein structure. Finally, molecular dynamics simulations revealed 5 variants BRCA2 Cys2646Tyr, Asp2665Val, Trp2619Arg, Trp2619Ser and Tyr2660Cys can alter the folding pattern and need further validation using in vitro, structural and in vivo studies to classify as pathogenic.Communicated by Ramaswamy H. Sarma.

In silico structural analysis of secretory clusterin to assess pathogenicity of mutations identified in the evolutionarily conserved regions.

Clusterin (CLU) is a secreted glycoprotein, heterodimeric in nature, and is expressed in a wide variety of tissues and body fluids such as serum and plasma. CLU has also been known to be a promising biomarker for cell death, malignancy, cancer progression, and resistance development. However, the lack of a CLU crystal structure obstructs understanding the possible role of reported mutations on the structure, and the subsequent effects on downstream signaling pathways and cancer progression. Considering the importance of crystal structure, a model structure of the pre-secretory isoform of CLU was built to predict the effect of mutations at the molecular level. Ab initio model was built using RaptorX, and loop refinement and energy minimization were carried out with ModLoop, ModRefiner, and GalaxyWeb servers. The cancer associated mutational spectra of CLU was retrieved from the cBioPortal server and 117 unique missense mutations were identified. Evolutionarily conserved regions and pathogenicity of mutations identified in CLU were analyzed using ConSurf and Rhapsody, respectively. Furthermore, sequence and structure-based mutational analysis were carried out with iSTABLE, DynaMut and PremPS servers. Molecular dynamics simulations were carried out with GROMACS for 50 ns to determine the stability of the wild type and mutant protein structures. A dynamically stable model structure of pre-secretory CLU (psCLU) which has high concurrence with the sequence based secondary structure predictions has been explored. Changes in the intra-atomic interactions and folding pattern between wild type and mutant structures were observed. To our conclusion, eleven mutations with the highest structural and functional significance have been predicted to have pathogenic and deleterious effects.Communicated by Ramaswamy H. Sarma.

In Silico-Based Structural Evaluation to Categorize the Pathogenicity of Mutations Identified in the RAD Class of Proteins.

RAD genes, known as double-strand break repair proteins, play a major role in maintaining the genomic integrity of a cell by carrying out essential DNA repair functions via double-strand break repair pathways. Mutations in the RAD class of proteins show high susceptibility to breast and ovarian cancers; however, adequate research on the mutations identified in these genes has not been extensively reported for their deleterious effects. Changes in the folding pattern of RAD proteins play an important role in protein-protein interactions and also functions. Missense mutations identified from four cancer databases, cBioPortal, COSMIC, ClinVar, and gnomAD, cause aberrant conformations, which may lead to faulty DNA repair mechanisms. It is therefore necessary to evaluate the effects of pathogenic mutations of RAD proteins and their subsequent role in breast and ovarian cancers. In this study, we have used eight computational prediction servers to analyze pathogenic mutations and understand their effects on the protein structure and function. A total of 5122 missense mutations were identified from four different cancer databases, of which 1165 were predicted to be pathogenic using at least five pathogenicity prediction servers. These mutations were characterized as high-risk mutations based on their location in the conserved domains and subsequently subjected to structural stability characterization. The mutations included in the present study were selected from clinically relevant mutants in breast cancer pedigrees. Comparative folding patterns and intra-atomic interaction results showed alterations in the structural behavior of RAD proteins, specifically RAD51C triggered by mutations G125V and L138F and RAD51D triggered by mutations S207L and E233G.

In-silico and structure-based assessment to evaluate pathogenicity of missense mutations associated with non-small cell lung cancer identified in the Eph-ephrin class of proteins.

Ephs belong to the largest family of receptor tyrosine kinase and are highly conserved both sequentially and structurally. The structural organization of Eph is similar to other receptor tyrosine kinases; constituting the extracellular ligand binding domain, a fibronectin domain followed by intracellular juxtamembrane kinase, and SAM domain. Eph binds to respective ephrin ligand, through the ligand binding domain and forms a tetrameric complex to activate the kinase domain. Eph-ephrin regulates many downstream pathways that lead to physiological events such as cell migration, proliferation, and growth. Therefore, considering the importance of Eph-ephrin class of protein in tumorigenesis, 7,620 clinically reported missense mutations belonging to the class of variables of unknown significance were retrieved from cBioPortal and evaluated for pathogenicity. Thirty-two mutations predicted to be pathogenic using SIFT, Polyphen-2, PROVEAN, SNPs&GO, PMut, iSTABLE, and PremPS in-silico tools were found located either in critical functional regions or encompassing interactions at the binding interface of Eph-ephrin. However, seven were reported in nonsmall cell lung cancer (NSCLC). Considering the relevance of receptor tyrosine kinases and Eph in NSCLC, these seven mutations were assessed for change in the folding pattern using molecular dynamic simulation. Structural alterations, stability, flexibility, compactness, and solvent-exposed area was observed in EphA3 Trp790Cys, EphA7 Leu749Phe, EphB1 Gly685Cys, EphB4 Val748Ala, and Ephrin A2 Trp112Cys. Hence, it can be concluded that the evaluated mutations have potential to alter the folding pattern and thus can be further validated by in-vitro, structural and in-vivo studies for clinical management.

Evaluation of Pathogenicity and Structural Alterations for the Mutations Identified in the Conserved Region of the C-Terminal Kinase Domain of Human-Ribosomal S6 Kinase 1.

Human-ribosomal s6 kinase 1 (h-RSK1) is an effector kinase of the Ras/MAPK signaling pathway, which is involved in the regulation of the cell cycle, proliferation, and survival. RSKs comprise two functionally distinct kinase domains at the N-terminal (NTKD) and C-terminal (CTKD) separated by a linker region. The mutations in RSK1 may have the potential to provide an extra benefit to the cancer cell to proliferate, migrate, and survive. The present study focuses on evaluating the structural basis for the missense mutations identified at the C-terminal kinase domain of human-RSK1. A total of 139 mutations reported on RSK1 were retrieved from cBioPortal, where 62 were located at the CTKD region. Furthermore, 10 missense mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, Arg726Gln, His533Asn, Pro613Leu, Ser720Cys, Arg725Gln, and Ser732Phe were predicted to be deleterious using in silico tools. To our observation, these mutations are located in the evolutionarily conserved region of RSK1 and shown to alter the inter- and intramolecular interactions and also the conformational stability of RSK1-CTKD. The molecular dynamics (MD) simulation study further revealed that the five mutations Arg434Pro, Thr701Met, Ala704Thr, Arg725Trp, and Arg726Gln showed maximum structural alterations in RSK1-CTKD. Thus, based on the in silico and MD simulation analysis, it can be concluded that the reported mutations may serve as potential candidates for further functional studies.