Investigation of binding phenomenon of NSP3 and p130Cas mutants and their effect on cell signalling

Unraveling the molecular mechanism of novel leukemia mutations on NTRK2 (A203T & R458G) and NTRK3 (E176D & L449F) genes using molecular dynamics simulations approach

Background: NTRK1, NTRK2, and NTRK3 are members of the neurotrophic receptor tyrosine kinases (NTRK) family, which encode TrkA, TrkB, and TrkC receptors, respectively. Hematologic cancers are also linked to point mutations in the NTRK gene's kinase domain. Trk fusions are the most common genetic change associated with oncogenic activity in Trk-driven liquid tumors. Thus, point mutations in NTRK genes may also play a role in tumorigenesis. The structural and functional effect of mutations in Trk-B & Trk-C proteins remains unclear. Methods: In this research, Homology (threading-based approach) modeling and the all-atom molecular dynamics simulations approaches are applied to examine the structural and functional behavior of native and mutant Trk-B and Trk-C proteins at the molecular level. Results: The result of this study reveals how the mutations in Trk-B (A203T & R458G) and Trk-C (E176D & L449F) proteins lost their stability and native conformations. The Trk-B mutant A203T became more flexible than the native protein, whereas the R458G mutation became more rigid than the native conformation of the Trk-B protein. Also, the Trk-C mutations (E176D & L449F) become more rigid compared to the native structure. Conclusions: This structural transition may interrupt the function of Trk-B and Trk-C proteins. Observing the impact of NTRK-2/3 gene alterations at the atomic level could aid in discovering a viable treatment for Trk-related leukemias.

The effect of mutation on neurotoxicity reduction of new chimeric reteplase, a computational study

Background and purpose

Excitotoxicity in nerve cells is a type of neurotoxicity in which excessive stimulation of receptors (such as N-methyl-d-aspartate glutamate receptors (NMDAR)) leads to the influx of high-level calcium ions into cells and finally cell damage or death. This complication can occur after taking some of the plasminogen activators like tissue plasminogen activator and reteplase. The interaction of the kringle2 domain in such plasminogen activator with the amino-terminal domain (ATD) of the NR1 subunit of NMDAR finally leads to excitotoxicity. In this study, we assessed the interaction of two new chimeric reteplase, mutated in the kringle2 domain, with ATD and compared the interaction of wild-type reteplase with ATD, computationally.

Experimental approach

Homology modeling, protein docking, molecular dynamic simulation, and molecular dynamics trajectory analysis were used for the assessment of this interaction.

Findings/Results

The results of the free energy analysis between reteplase and ATD (wild reteplase: -2127.516 ± 0.0, M1-chr: -1761.510 ± 0.0, M2-chr: -521.908 ± 0.0) showed lower interaction of this chimeric reteplase with ATD compared to the wild type.

Conclusion and implications

The decreased interaction between two chimeric reteplase and ATD of NR1 subunit in NMDAR which leads to lower neurotoxicity related to these drugs, can be the start of a way to conduct more tests and if the results confirm this feature, they can be considered potential drugs in acute ischemic stroke treatment.

The Structural Effect of FLT3 Mutations at 835th Position and Their Interaction with Acute Myeloid Leukemia Inhibitors: In Silico Approach

FMS-like tyrosine kinase 3 (FLT3) gene mutations have been found in more than one-third of Acute Myeloid Leukemia (AML) cases. The most common point mutation in FLT3 occurs at the 835th residue (D835A/E/F/G/H/I/N/V/Y), in the activation loop region. The D835 residue is critical in maintaining FLT3 inactive conformation; these mutations might influence the interaction with clinically approved AML inhibitors used to treat the AML. The molecular mechanism of each of these mutations and their interactions with AML inhibitors at the atomic level is still unknown. In this manuscript, we have investigated the structural consequence of native and mutant FLT-3 proteins and their molecular mechanisms at the atomic level, using molecular dynamics simulations (MDS). In addition, we use the molecular docking method to investigate the binding pattern between the FLT-3 protein and AML inhibitors upon mutations. This study apparently elucidates that, due to mutations in the D835, the FLT-3 structure loses its conformation and becomes more flexible compared to the native FLT3 protein. These structural changes are suggested to contribute to the relapse and resistance responses to AML inhibitors. Identifying the effects of FLT3 at the molecular level will aid in developing a personalized therapeutic strategy for treating patients with FLT-3-associated AML.

Investigation of nonsynonymous mutations in the spike protein of SARS-CoV-2 and its interaction with the ACE2 receptor by molecular docking and MM/GBSA approach

COVID-19 is an infectious and pathogenic viral disease caused by SARS-CoV-2 that leads to septic shock, coagulation dysfunction, and acute respiratory distress syndrome. The spreading rate of SARS-CoV-2 is higher than MERS-CoV and SARS-CoV. The receptor-binding domain (RBD) of the Spike-protein (S-protein) interacts with the human cells through the host angiotensin-converting enzyme 2 (ACE2) receptor. However, the molecular mechanism of pathological mutations of S-protein is still unclear. In this perspective, we investigated the impact of mutations in the S-protein and their interaction with the ACE2 receptor for SAR-CoV-2 viral infection. We examined the stability of pathological nonsynonymous mutations in the S-protein, and the binding behavior of the ACE2 receptor with the S-protein upon nonsynonymous mutations using the molecular docking and MM_GBSA approaches. Using the extensive bioinformatics pipeline, we screened the destabilizing (L8V, L8W, L18F, Y145H, M153T, F157S, G476S, L611F, A879S, C1247F, and C1254F) and stabilizing (H49Y, S50L, N501Y, D614G, A845V, and P1143L) nonsynonymous mutations in the S-protein. The docking and binding free energy (ddG) scores revealed that the stabilizing nonsynonymous mutations show increased interaction between the S-protein and the ACE2 receptor compared to native and destabilizing S-proteins and that they may have been responsible for the virulent high level. Further, the molecular dynamics simulation (MDS) approach reveals the structural transition of mutants (N501Y and D614G) S-protein. These insights might help researchers to understand the pathological mechanisms of the S-protein and provide clues regarding mutations in viral infection and disease propagation. Further, it helps researchers to develop an efficient treatment approach against this SARS-CoV-2 pandemic.

Effect of novel leukemia mutations (K75E & E222K) on interferon regulatory factor 1 and its interaction with DNA: insights from molecular dynamics simulations and docking studies

Interferon regulatory factor 1 (IRF-1) plays a vital role in cell proliferation and cell differentiation by acting as a tumor suppressor gene and its role is linked to various types of cancers, including leukemia and pre-leukemia myelodysplasia. Mutations in the coding region of the IRF-1 are likely to influence the IRF-1 and its DNA binding affinity. The molecular mechanism of the DNA recognition with the IRF-1 protein upon mutations is still unknown. In this study, we have elucidated the structural and functional behavior of the wild-type and mutant (K75E and E222K) IRF-1 proteins and their corresponding molecular mechanisms with DNA recognition at the molecular level, using molecular dynamics simulations. Furthermore, we also applied the docking approach to examine the binding between the IRF-1 protein and DNA upon mutations. This study evidently explains that, due to mutations, the IRF-1 structure loses its stability and becomes more flexible than the wild-type protein. This structural loss might affect IRF-1-DNA interaction and lead to the inhibition of cancer suppression. Identifying the effects of IRF-1 at the molecular level will be beneficial for designing drugs for IRF-1 associated cancers. These drugs should be designed so that they can help reactivate the IRF-1 function, by increasing the transcriptional activity, to treat leukemia.

Design and production of new chimeric reteplase with enhanced fibrin affinity: a theoretical and experimental study

Plasminogen activators (PAs) are widely used for treatment of disorders caused by clot formation. Fibrin specific PAs are safe drugs from this group because of reducing the incidence of hemorrhage. The newer generation of PAs like tenecteplase, reteplase and desmoteplase were designed with the aim of achieving desirable properties such as improving specificity and affinity to fibrin and increasing half-life. Protein engineering and using of theoretical methods can help to rational and reliable design of new PAs with a set of favorable properties. In the present study, two new chimeric reteplase named M1-chr and M2-chr were designed with the aim of enhancing fibrin affinity also some potential properties include of increasing resistance to plasminogen activator inhibitor-1 and decreasing neurotoxicity. So, finger domain of desmoteplase was added to reteplase as a high fibrin specific domain. Some other point mutations were considering to achieve other mentioned properties. Three dimensional structure of wild-type reteplase and mutants were created by homology modeling and were evaluated by molecular dynamic simulation. Then, mutants docked to fibrin by HADDOCK web tools. Result of theoretical section verified the stability of mutants' structures. Also showed better interaction between M1-chr with fibrin than M2-chr. Wild-type and mutants were produced in bacterial expression system. Experimental assessment showed both mutants have appropriate enzymatic activity also 1.9-fold fibrin binding ability compared to wild-type. Therefore, this study offers new thrombolytic drugs with desirable properties specially enhanced fibrin affinity so they can represent a promising future in cost-effective production of favorable thrombolytic drugs.Communicated by Ramaswamy H. Sarma.

Dynamic perspectives into the mechanisms of mutation-induced p53-DNA binding loss and inactivation using active perturbation theory: Structural and molecular insights toward the design of potent reactivators in cancer therapy

The DNA-binding ability of p53 represents the crux of its tumor suppressive activities, which involves transcriptional activation of target genes responsible for apoptosis and cell-cycle arrest. Mutational occurrences within or in close proximity to the DNA-binding surface of p53 have accounted for the loss of direct DNA-binding ability and inactivation implicated in many cases of cancer. Moreover, the design of therapeutic compounds that can restore DNA-binding ability in p53 mutants has been identified as a way forward in curtailing their oncogenic activities. However, there is still the need for more insights into evaluate the perturbations that occur at the DNA-binding interface of mp53 relative to DNA-binding loss, inactivation, and design of potent reactivators, hence the purpose of this study. Therefore, we evaluated p53-structural (R175H) and contact (R273C) mutational effects using tunnel perturbation analysis and other computational tools. We identified significant perturbations in the active tunnels of p53, which resulted in altered geometry and loss, unlike in the wild-type p53. This corroborated with structural, DNA-binding, and interaction network analysis, which showed that loss of flexibility, repulsion of DNA-interactive residues, and instability occurred at the binding interface of both mutants. Also, these mutations altered bonding interactions and network topology at the DNA-binding interface, resulting in the reduction of p53-DNA binding proximity and affinity. Therefore, these findings would aid the structure-based design of novel chemical entities capable of restoring p53-DNA binding and activation.

Bioinformatic Analysis of GJB2 Gene Missense Mutations

Gap junction beta 2 (GJB2) gene is the most commonly mutated connexin gene in patients with autosomal recessive and dominant hearing loss. According to Ensembl (release 74) database, 1347 sequence variations are reported in the GJB2 gene and about 13.5% of them are categorized as missense SNPs or nonsynonymous variant. Because of the high incidence of GJB2 mutations in hearing loss patients, revealing the molecular effect of GJB2 mutations on protein structure may also provide clear point of view regarding the molecular etiology of deafness. Hence, the aim of this study is to analyze structural and functional consequences of all known GJB2 missense variations to the Cx26 protein by applying multiple bioinformatics methods. Two-hundred and eleven nonsynonymous variants were collected from Ensembl release 74, Leiden Open Variation Database (LOVD) and The Human Gene Mutation Database (HGMD). A number of bioinformatic tools were utilized for predicting the effect of GJB2 missense mutations at the sequence, structural, and functional levels. Some of the mutations were found to locate highly conserved regions and have structural and functional properties. Moreover, GJB2 mutations were also found to affect Cx26 protein at the molecular level via loss or gain of disorder, catalytic site, and post-translational modifications, including methylation, glycosylation, and ubiquitination. Findings, presented here, demonstrated the application of bioinformatic algorithms to predict the effects of mutations causing hearing impairment. I expect, this type of analysis will serve as a start point for future experimental evaluation of the GJB2 gene mutations and it will also be helpful in evaluating other deafness-related gene mutations.

Structure and Function of p53-DNA Complexes with Inactivation and Rescue Mutations: A Molecular Dynamics Simulation Study

The tumor suppressor protein p53 can lose its function upon DNA-contact mutations (R273C and R273H) in the core DNA-binding domain. The activity can be restored by second-site suppressor or rescue mutations (R273C_T284R, R273H_T284R, and R273H_S240R). In this paper, we elucidate the structural and functional consequence of p53 proteins upon DNA-contact mutations and rescue mutations and the underlying mechanisms at the atomic level by means of molecular dynamics simulations. Furthermore, we also apply the docking approach to investigate the binding phenomena between the p53 protein and DNA upon DNA-contact mutations and rescue mutations. This study clearly illustrates that, due to DNA-contact mutants, the p53 structure loses its stability and becomes more rigid than the native protein. This structural loss might affect the p53-DNA interaction and leads to inhibition of the cancer suppression. Rescue mutants (R273C_T284R, R273H_T284R and R273H_S240R) can restore the functional activity of the p53 protein upon DNA-contact mutations and show a good interaction between the p53 protein and a DNA molecule, which may lead to reactivate the cancer suppression function. Understanding the effects of p53 cancer and rescue mutations at the molecular level will be helpful for designing drugs for p53 associated cancer diseases. These drugs should be designed so that they can help to inhibit the abnormal function of the p53 protein and to reactivate the p53 function (cell apoptosis) to treat human cancer.

Identification and analysis of driver missense mutations using rotation forest with feature selection

Identifying cancer-associated mutations (driver mutations) is critical for understanding the cellular function of cancer genome that leads to activation of oncogenes or inactivation of tumor suppressor genes. Many approaches are proposed which use supervised machine learning techniques for prediction with features obtained by some databases. However, often we do not know which features are important for driver mutations prediction. In this study, we propose a novel feature selection method (called DX) from 126 candidate features' set. In order to obtain the best performance, rotation forest algorithm was adopted to perform the experiment. On the train dataset which was collected from COSMIC and Swiss-Prot databases, we are able to obtain high prediction performance with 88.03% accuracy, 93.9% precision, and 81.35% recall when the 11 top-ranked features were used. Comparison with other various techniques in the TP53, EGFR, and Cosmic2plus datasets shows the generality of our method.

Evidence of colorectal cancer-associated mutation in MCAK: a computational report

Computational prediction of disease-associated non-synonymous polymorphism (nsSNP) has provided a significant platform to filter out the pathological mutations from large pool of SNP datasets at a very low cost input. Several methodologies and complementary protocols have been previously implemented and has provided significant prediction results. Although the previously implicated prediction methods were capable of investigating the most likely deleterious nsSNPs, but due to the lack of genotype-phenotype association analysis, the prediction results lacked in accuracy level. In this work we implemented the computational compilation of protein conformational changes as well as the probable disease-associated phenotypic outcomes. Our result suggested E403K mutation in mitotic centromere-associated kinesin protein as highly damaging and showed strong concordance to the previously observed colorectal cancer mutations aggregation tendency and energy value changes. Moreover, the molecular dynamics simulation results showed major loss in conformation and stability of mutant N-terminal kinesin-like domain structure. The result obtained in this study will provide future prospect of computational approaches in determining the SNPs that may affect the native conformation of protein structure and lead to cancer-associated disorders.

CAS proteins in health and disease: an update

The CAS family of scaffolding proteins has increasingly attracted scrutiny as important for regulation of cancer-associated signaling. BCAR1 (also known as p130Cas), NEDD9 (HEF1, Cas-L), EFS (Sin), and CASS4 (HEPL) are regulated by and mediate cell attachment, growth factor, and chemokine signaling. Altered expression and activity of CAS proteins are now known to promote metastasis and drug resistance in cancer, influence normal development, and contribute to the pathogenesis of heart and pulmonary disease. In this article, we provide an update on recently published studies describing signals regulating and regulated by CAS proteins, and evidence for biological activity of CAS proteins in normal development, cancer, and other pathological conditions.

In-silico screening of cancer associated mutation on PLK1 protein and its structural consequences

The Polo-like kinases (Plks) are a conserved subfamily of serine-threonine protein kinases that have significant roles in cell proliferation. The serine/threonine protein kinases or polo-like kinase 1 (PLK1) exist in centrosome during interphase and is an important regulatory enzyme in cell cycle progression during M phase. Mutations in mammalian PLK1 were found to be over expressed in various human cancers and it is disrupting the binding ability of polo box domain with target peptide. In this analysis we implemented a computational approach to filter the most deleterious and cancer associated mutation on PLK1 protein. We found W414F as the most deleterious and cancer associated by Polyphen 2.0, SIFT, I-mutant 3.0, PANTHER, PhD-SNP, SNP&GO, Mutpred and Dr Cancer tools. Molecular docking and molecular dynamics simulation (MDS) approach was used to investigate the structural and functional behavior of PLK1 protein upon mutation. MDS and docking results showed stability loss in mutant PLK1 protein. Due to mutation, PLK1 protein became more flexible and alters the dynamic property of protein which might affect the interaction with target peptide and leads to cell proliferation. Our study provided a well designed computational methodology to examine the cancer associated nsSNPs and their molecular mechanism. It further helps scientists to develop a drug therapy against PLK1 cancer-associated diseases.

AKT kinase pathway: a leading target in cancer research

AKT1, a serine/threonine-protein kinase also known as AKT kinase, is involved in the regulation of various signalling downstream pathways including metabolism, cell proliferation, survival, growth, and angiogenesis. The AKT kinases pathway stands among the most important components of cell proliferation mechanism. Several approaches have been implemented to design an efficient drug molecule to target AKT kinases, although the promising results have not been confirmed. In this paper we have documented the detailed molecular insight of AKT kinase protein and proposed a probable doxorubicin based approach in inhibiting miR-21 based cancer cell proliferation. Moreover, the inhibition of miR-21 activation by raising the FOXO3A concentration seems promising in reducing miR-21 mediated cancer activation in cell. Furthermore, the use of next generation sequencing and computational drug design approaches will greatly assist in designing a potent drug molecule against the associated cancer cases.

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.

In silico screening and molecular dynamics simulation of disease-associated nsSNP in TYRP1 gene and its structural consequences in OCA3

Oculocutaneous albinism type III (OCA3), caused by mutations of TYRP1 gene, is an autosomal recessive disorder characterized by reduced biosynthesis of melanin pigment in the hair, skin, and eyes. The TYRP1 gene encodes a protein called tyrosinase-related protein-1 (Tyrp1). Tyrp1 is involved in maintaining the stability of tyrosinase protein and modulating its catalytic activity in eumelanin synthesis. Tyrp1 is also involved in maintenance of melanosome structure and affects melanocyte proliferation and cell death. In this work we implemented computational analysis to filter the most probable mutation that might be associated with OCA3. We found R326H and R356Q as most deleterious and disease associated by using PolyPhen 2.0, SIFT, PANTHER, I-mutant 3.0, PhD-SNP, SNP&GO, Pmut, and Mutpred tools. To understand the atomic arrangement in 3D space, the native and mutant (R326H and R356Q) structures were modelled. Finally the structural analyses of native and mutant Tyrp1 proteins were investigated using molecular dynamics simulation (MDS) approach. MDS results showed more flexibility in native Tyrp1 structure. Due to mutation in Tyrp1 protein, it became more rigid and might disturb the structural conformation and catalytic function of the structure and might also play a significant role in inducing OCA3. The results obtained from this study would facilitate wet-lab researches to develop a potent drug therapies against OCA3.

Cancer associated E17K mutation causes rapid conformational drift in AKT1 pleckstrin homology (PH) domain

BACKGROUND

AKT1 (v-akt murine thymoma viral oncogene homologue 1) kinase is one of the most frequently activated proliferated and survival pathway of cancer. Recently it has been shown that E17K mutation in the Pleckstrin Homology (PH) domain of AKT1 protein leads to cancer by amplifying the phosphorylation and membrane localization of protein. The mutant has shown resistance to AKT1/2 inhibitor VIII drug molecule. In this study we have demonstrated the detailed structural and molecular consequences associated with the activity regulation of mutant protein.

METHODS

The docking score exhibited significant loss in the interaction affinity to AKT1/2 inhibitor VIII drug molecule. Furthermore, the molecular dynamics simulation studies presented an evidence of rapid conformational drift observed in mutant structure.

RESULTS

There was no stability loss in mutant as compared to native structure and the major cation-π interactions were also shown to be retained. Moreover, the active residues involved in membrane localization of protein exhibited significant rise in NHbonds formation in mutant. The rise in NHbond formation in active residues accounts for the 4-fold increase in the membrane localization potential of protein.

CONCLUSION

The overall result suggested that, although the mutation did not induce any stability loss in structure, the associated pathological consequences might have occurred due to the rapid conformational drifts observed in the mutant AKT1 PH domain.

GENERAL SIGNIFICANCE

The methodology implemented and the results obtained in this work will facilitate in determining the core molecular mechanisms of cancer-associated mutations and in designing their potential drug inhibitors.