Trifluoperazine regulation of calmodulin binding to Fas: a computational study

Artificial Proteins Designed from G3LEA Contribute to Enhancement of Oxidation Tolerance in E. coli in a Chaperone-like Manner

G3LEA is a family of proteins that exhibit chaperone-like activity when under distinct stress. In previous research, DosH was identified as a G3LEA protein from model extremophile-Deinococcus radiodurans R1 with a crucial core HD domain consisting of eight 11-mer motifs. However, the roles of motifs participating in the process of resistance to stress and their underlying mechanisms remain unclear. Here, eight different proteins with tandem repeats of the same motif were synthesized, named Motif1-8, respectively, whose function and structure were discussed. In this way, the role of each motif in the HD domain can be comprehensively analyzed, which can help in finding possibly crucial amino acid sites. Circular dichroism results showed that all proteins were intrinsically ordered in phosphate buffer, and changed into more α-helical ordered structures with the addition of trifluoroethanol and glycerol. Transformants expressing artificial proteins had significantly higher stress resistance to oxidation, desiccation, salinity and freezing compared with the control group; E. coli with Motif1 and Motif8 had more outstanding performance in particular. Moreover, enzymes and membrane protein protection viability suggested that Motif1 and Motif8 had more positive influences on various molecules, demonstrating a protective role in a chaperone-like manner. Based on these results, the artificial proteins synthesized according to the rule of 11-mer motifs have a similar function to wildtype protein. Regarding the sequence in all motifs, there are more amino acids to produce H bonds and α-helices, and more amino acids to promote interaction between proteins in Motif1 and Motif8; in addition, considering linkers, there are possibly more amino acids forming α-helix and binding substrates in these two proteins, which potentially provides some ideas for us to design potential ideal stress-response elements for synthetic biology. Therefore, the amino acid composition of the 11-mer motif and linker is likely responsible for its biological function.

Molecular and structural basis of interactions of vitamin D3 hydroxyderivatives with aryl hydrocarbon receptor (AhR): An integrated experimental and computational study

To better understand the molecular and structural basis underlying the interaction of vitamin D3 hydroxyderivatives with AhR, molecular simulation was used to probe the binding of 1,20(OH)₂D3, 1,25(OH)₂D3, 20,23(OH)₂D3 and 20(OH)D3 to AhR. qPCR showed that vitamin D3 derivatives stimulate expression of cyp1A1 and cyp1B1 genes that are downstream targets of AhR signaling. These secosteroids stimulated the translocation of the AhR to the nucleus, as measured by flow cytometry and western blotting. Molecular dynamics simulations were used to model the binding of vitamin D3 derivatives to AhR to examine their influence on the structure, conformation and dynamics of the AhR ligand binding domain (LBD). Binding thermodynamics, conformation, secondary structure, dynamical motion and electrostatic potential of AhR were analyzed. The molecular docking scores and binding free energy were all favorable for the binding of D3 derivatives to the AhR. These established ligands and the D3 derivatives are predicted to have different patterns of hydrogen bond formation with the AhR, and varied residue conformational fluctuations and dynamical motion for the LBD. These changes could alter the shape, size and electrostatic potential distribution of the ligand binding pocket, contributing to the different binding affinities of AhR for the natural ligands and D3 derivatives.

Vitamin D3 and its hydroxyderivatives as promising drugs against COVID-19: a computational study

The epidemiologic correlation between the poor prognosis of SARS-CoV-2 infection and vitamin D deficiency has been observed worldwide, however, their molecular mechanisms are not fully understood. In this study, we used combined molecular docking, molecular dynamics simulations and binding free energy analyses to investigate the potentials of vitamin D3 and its hydroxyderivatives as TMPRSS2 inhibitor and to inhibit the SARS-CoV-2 receptor binding domain (RBD) binding to angiotensin-converting enzyme 2 (ACE2), as well as to unveil molecular and structural basis of 1,25(OH)₂D3 capability to inhibit ACE2 and SARS-CoV-2 RBD interactions. The results show that vitamin D3 and its hydroxyderivatives are favorable to bind active site of TMPRSS2 and the binding site(s) between ACE2 and SARS-CoV2-RBD, which indicate that vitamin D3 and its biologically active hydroxyderivatives can serve as TMPRSS2 inhibitor and can inhibit ACE2 binding of SARS-CoV-2 RBD to prevent SARS-CoV-2 entry. Interaction of 1,25(OH)₂D3 with SARS-CoV-2 RBD and ACE2 resulted in the conformation and dynamical motion changes of the binding surfaces between SARS-CoV-2 RBD and ACE2 to interrupt the binding of SARS-CoV-2 RBD with ACE2. The interaction of 1,25(OH)₂D3 with TMPRSS2 also caused the conformational and dynamical motion changes of TMPRSS2, which could affect TMPRSS2 to prime SARS-CoV-2 spike proteins. Our results propose that vitamin D3 and its biologically active hydroxyderivatives are promising drugs or adjuvants in the treatment of COVID-19. Communicated by Ramaswamy H. Sarma.

Vitamin D and lumisterol derivatives can act on liver X receptors (LXRs)

The interactions of derivatives of lumisterol (L3) and vitamin D3 (D3) with liver X receptors (LXRs) were investigated. Molecular docking using crystal structures of the ligand binding domains (LBDs) of LXRα and β revealed high docking scores for L3 and D3 hydroxymetabolites, similar to those of the natural ligands, predicting good binding to the receptor. RNA sequencing of murine dermal fibroblasts stimulated with D3-hydroxyderivatives revealed LXR as the second nuclear receptor pathway for several D3-hydroxyderivatives, including 1,25(OH)2D3. This was validated by their induction of genes downstream of LXR. L3 and D3-derivatives activated an LXR-response element (LXRE)-driven reporter in CHO cells and human keratinocytes, and by enhanced expression of LXR target genes. L3 and D3 derivatives showed high affinity binding to the LBD of the LXRα and β in LanthaScreen TR-FRET LXRα and β coactivator assays. The majority of metabolites functioned as LXRα/β agonists; however, 1,20,25(OH)3D3, 1,25(OH)2D3, 1,20(OH)2D3 and 25(OH)D3 acted as inverse agonists of LXRα, but as agonists of LXRβ. Molecular dynamics simulations for the selected compounds, including 1,25(OH)2D3, 1,20(OH)2D3, 25(OH)D3, 20(OH)D3, 20(OH)L3 and 20,22(OH)2L3, showed different but overlapping interactions with LXRs. Identification of D3 and L3 derivatives as ligands for LXRs suggests a new mechanism of action for these compounds.

Role of Overexpressed Transcription Factor FOXO1 in Fatal Cardiovascular Septal Defects in Patau Syndrome: Molecular and Therapeutic Strategies

Patau Syndrome (PS), characterized as a lethal disease, allows less than 15% survival over the first year of life. Most deaths owe to brain and heart disorders, more so due to septal defects because of altered gene regulations. We ascertained the cytogenetic basis of PS first, followed by molecular analysis and docking studies. Thirty-seven PS cases were referred from the Department of Pediatrics, King Abdulaziz University Hospital to the Center of Excellence in Genomic Medicine Research, Jeddah during 2008 to 2018. Cytogenetic analyses were performed by standard G-band method and trisomy13 were found in all the PS cases. Studies have suggested that genes of chromosome 13 and other chromosomes are associated with PS. We, therefore, did molecular pathway analysis, gene interaction, and ontology studies to identify their associations. Genomic analysis revealed important chr13 genes such as FOXO1, Col4A1, HMGBB1, FLT1, EFNB2, EDNRB, GAS6, TNFSF1, STARD13, TRPC4, TUBA3C, and TUBA3D, and their regulatory partners on other chromosomes associated with cardiovascular disorders, atrial and ventricular septal defects. There is strong indication of involving FOXO1 (Forkhead Box O1) gene-a strong transcription factor present on chr13, interacting with many septal defects link genes. The study was extended using molecular docking to find a potential drug lead for overexpressed FOXO1 inhibition. The phenothiazine and trifluoperazine showed efficiency to inhibit overexpressed FOXO1 protein, and could be potential drugs for PS/trisomy13 after validation.

Calmodulin antagonist enhances DR5-mediated apoptotic signaling in TRA-8 resistant triple negative breast cancer cells

Patients with triple negative breast cancer (TNBC) have no successful "targeted" treatment modality, which represents a priority for novel therapy strategies. Upregulated death receptor 5 (DR5) expression levels in breast cancer cells compared to normal cells enable TRA-8, a DR5 specific agonistic antibody, to specifically target malignant cells for apoptosis without inducing normal hepatocyte apoptosis. Drug resistance is a common obstacle in TRAIL-based therapy for TNBC. Calmodulin (CaM) is overexpressed in breast cancer. In this study, we characterized the novel function of CaM antagonist in enhancing TRA-8 induced cytotoxicity in TRA-8 resistant TNBC cells and its underlying molecular mechanisms. Results demonstrated that CaM antagonist(s) enhanced TRA-8 induced cytotoxicity in a concentration and time-dependent manner for TRA-8 resistant TNBC cells. CaM directly bound to DR5 in a Ca²⁺ dependent manner, and CaM siRNA promoted DR5 recruitment of FADD and caspase-8 for DISC formation and TRA-8 activated caspase cleavage for apoptosis in TRA-8 resistant TNBC cells. CaM antagonist, trifluoperazine, enhanced TRA-8 activated DR5 oligomerization, DR5-mediated DISC formation, and TRA-8 activated caspase cleavage for apoptosis, and decreased anti-apoptotic pERK, pAKT, XIAP, and cIAP-1 expression in TRA-8 resistant TNBC cells. These results suggest that CaM could be a key regulator to mediate DR5-mediated apoptotic signaling, and suggests a potential strategy for using CaM antagonists to overcome drug resistance of TRAIL-based therapy for TRA-8 resistant TNBC.

Calmodulin Binding to Death Receptor 5-mediated Death-Inducing Signaling Complex in Breast Cancer Cells

Activation of death receptor-5 (DR5) leads to the formation of death-inducing signaling complex (DISC) for apoptotic signaling. TRA-8, a DR5 specific agonistic antibody, has demonstrated significant cytotoxic activity in vitro and in vivo without inducing hepatotoxicity. Calmodulin (CaM) that is overexpressed in breast cancer plays a critical role in regulating DR5-mediated apoptosis. However, the mechanism of CaM in regulating DR5-mediated apoptotic signaling remains unknown. In this study, we characterized CaM binding to DR5-mediated DISC for apoptosis in TRA-8 sensitive breast cancer cell lines using co-immunoprecipitation, fluorescence microscopic imaging, caspase signaling analysis, and cell viability assay. Results show that upon DR5 activation, CaM was recruited into DR5-mediated DISC in a calcium dependent manner. CaM antagonist, trifluoperazine (TFP), inhibited CaM recruitment into the DISC and attenuated DISC formation. DR5 oligomerization is critical for DISC formation for apoptosis. TFP decreased TRA-8 activated DR5 oligomerization, which was consistent with TFP's effect on DR5-mediated DISC formation. TFP and Ca²⁺ chelator, EGTA, impeded TRA-8-activated caspase-dependent apoptotic signaling, and TFP decreased TRA-8-induced cell cytotoxicity. These results demonstrated CaM binding to DR5-mediated DISC in a calcium dependent manner and may identify CaM as a key regulator of DR5-mediated DISC formation for apoptosis in breast cancer. J. Cell. Biochem. 118: 2285-2294, 2017. © 2017 Wiley Periodicals, Inc.

Structural Insight for Roles of DR5 Death Domain Mutations on Oligomerization of DR5 Death Domain-FADD Complex in the Death-Inducing Signaling Complex Formation: A Computational Study

Death receptor 5 (DR5)-induced apoptosis that prioritizes the death of tumor cells has been proposed as one of the promising cancer therapies. In this process, oligomerized DR5 death domain (DD) binding to Fas-associated death domain (FADD) leads to FADD activating caspase-8, which marks the formation of the death-inducing signaling complex (DISC) that initiates apoptosis. DR5 DD mutations found in cancer cells have been suggested to play an important pathological role, the mechanism through which those mutants prevent the DR5-activated DISC formation is not clear yet. This study sought to provide structural and molecular insight for the roles of four selected DR5 DD mutations (E355K, E367K, K415N, and L363F) in the oligomerization of DR5 DD-FADD complex during the DISC formation. Results from the molecular dynamics simulations show that the simulated mutants induce conformational, dynamical motions and interactions changes in the DR5 DD-FADD tetramer complex, including changes in a protein's backbone flexibility, less exposure of FADD DED's caspase-8 binding site, reduced H-bonding and hydrophobic contacts at the DR5 DD-FADD DD binding, altered distribution of the electrostatic potentials and correlated motions of residues, and reduced binding affinity of DR5 DD binding to FADD. This study provides structural and molecular insight for the influence of DR5 DD mutations on oligomerization of DR5 DD-FADD complex, which is expected to foster understanding of the DR5 DD mutants' resistance mechanism against DR5-activated DISC formation.

Structural Characterizations of the Fas Receptor and the Fas-Associated Protein with Death Domain Interactions

The Fas receptor is a representative death receptor, and the Fas-associated protein with death domain (FADD) is a crucial adapter protein needed to support the Fas receptor's activity. The Fas-FADD interactions constitute an important signaling pathway that ultimately induces apoptosis or programmed cell death in biological systems. The interactions responsible for this cell-death process are governed by the binding process of the Fas ligand to the Fas, followed by the caspase cascade activation. Using a computational approach, the present communication explores certain essential structural aspects of the Fas-FADD death domains and their interfacial interactions.

Molecular insight into the effect of lipid bilayer environments on thrombospondin-1 and calreticulin interactions

Thrombospondin-1 (TSP1) binding to cell surface calreticulin (CRT) stimulates the association of CRT with low-density lipoprotein (LDL) receptor-related protein (LRP1) to signal focal adhesion disassembly and engagement of cellular activities. A recent study demonstrated that membrane rafts are necessary for TSP1-mediated focal adhesion disassembly, but the molecular role of membrane rafts in mediating TSP1-CRT-LRP1 signaling is unknown. In this study, we investigated the effect of lipid bilayer environments on TSP1 and CRT interactions via atomically detailed molecular dynamics simulations. Results showed that the microscopic structural properties of lipid molecules and mesoscopic mechanical properties and electrostatic potential of the bilayer were significantly different between a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer and a raftlike lipid bilayer [a POPC/cholesterol (CHOL) raftlike lipid bilayer or a POPC/CHOL/sphingomyelin (SM) raftlike lipid bilayer], and the difference was enhanced by SM lipids in a raftlike lipid bilayer. These bilayer property differences affect the interactions of CRT with the bilayer, further influencing CRT conformation and TSP1-CRT interactions. A raftlike lipid bilayer stabilized CRT conformation as compared to a POPC bilayer environment. TSP1 binding to CRT resulted in a conformation for the CRT N-domain more "open" than that of the CRT P-domain in a raftlike lipid bilayer environment, which could facilitate binding of CRT to LRP1 to engage downstream signaling. The open conformational changes of CRT by binding to TSP1 in a raftlike lipid bilayer were enhanced by SM lipids in a lipid bilayer. The direct interactions of both the N- and P-domains of CRT with the bilayer contribute to the more open conformation of CRT in the TSP1-CRT complex on a raftlike lipid bilayer as compared to that on a POPC bilayer. The interactions of CRT or the TSP1-CRT complex with the lipid bilayer also caused CHOL molecules and/or lipids to be more coordinated and to aggregate into patchlike regions in the raftlike lipid bilayers. The lipid and CHOL molecule coordination and aggregation could in turn affect the interactions of CRT with the membrane raft, thereby altering TSP1-CRT interactions and CRT conformational changes that potentially regulate its interactions with LRP1. This study provides molecular insights into the role of lipid bilayer environments in TSP1-CRT interactions and in the CRT conformational changes that are predicted to facilitate binding of CRT to LRP1 to engage downstream signaling events.

Structural insights into the mechanism of calmodulin binding to death receptors

The death receptors Fas, p75NTRand DR6 are key components of extrinsically activated apoptosis. Characterization of how they interact with the adaptors is crucial in order to unravel the signalling mechanisms. However, the exact conformation that their intracellular death domain adopts upon binding downstream partners remains unclear. One model suggests that it adopts a typical compact fold, whilst a second model proposed an open conformation. Calmodulin (CaM), a major calcium sensor, has previously been reported to be one of the Fas adaptors that modulate apoptosis. This work reports that CaM also binds directly to the death domains of p75NTRand DR6, indicating that it serves as a common modulator of the death receptors. Two crystal structures of CaM in complexes with the corresponding binding regions of Fas and p75NTRare also reported. Interestingly, the precise CaM-binding sites were mapped to different regions: helix 1 in Fas and helix 5 in p75NTRand DR6. A novel 1–11 motif for CaM binding was observed in p75NTR. Modelling the complexes of CaM with full-length receptors reveals that the opening of the death domains would be essential in order to expose their binding sites for CaM. These results may facilitate understanding of the diverse functional repertoire of death receptors and CaM and provide further insights necessary for the design of potential therapeutic peptide agents.

Characterization of calmodulin-Fas death domain interaction: an integrated experimental and computational study

The Fas death receptor-activated death-inducing signaling complex (DISC) regulates apoptosis in many normal and cancer cells. Qualitative biochemical experiments demonstrate that calmodulin (CaM) binds to the death domain of Fas. The interaction between CaM and Fas regulates Fas-mediated DISC formation. A quantitative understanding of the interaction between CaM and Fas is important for the optimal design of antagonists for CaM or Fas to regulate the CaM-Fas interaction, thus modulating Fas-mediated DISC formation and apoptosis. The V254N mutation of the Fas death domain (Fas DD) is analogous to an identified mutant allele of Fas in lpr-cg mice that have a deficiency in Fas-mediated apoptosis. In this study, the interactions of CaM with the Fas DD wild type (Fas DD WT) and with the Fas DD V254N mutant were characterized using isothermal titration calorimetry (ITC), circular dichroism spectroscopy (CD), and molecular dynamics (MD) simulations. ITC results reveal an endothermic binding characteristic and an entropy-driven interaction of CaM with Fas DD WT or with Fas DD V254N. The Fas DD V254N mutation decreased the association constant (Ka) for CaM-Fas DD binding from (1.79 ± 0.20) × 10(6) to (0.88 ± 0.14) × 10(6) M(-1) and slightly increased a standard state Gibbs free energy (ΔG°) for CaM-Fas DD binding from -8.87 ± 0.07 to -8.43 ± 0.10 kcal/mol. CD secondary structure analysis and MD simulation results did not show significant secondary structural changes of the Fas DD caused by the V254N mutation. The conformational and dynamical motion analyses, the analyses of hydrogen bond formation within the CaM binding region, the contact numbers of each residue, and the electrostatic potential for the CaM binding region based on MD simulations demonstrated changes caused by the Fas DD V254N mutation. These changes caused by the Fas DD V254N mutation could affect the van der Waals interactions and electrostatic interactions between CaM and Fas DD, thereby affecting CaM-Fas DD interactions. Results from this study characterize CaM-Fas DD interactions in a quantitative way, providing structural and thermodynamic evidence of the role of the Fas DD V254N mutation in the CaM-Fas DD interaction. Furthermore, the results could help to identify novel strategies for regulating CaM-Fas DD interactions and Fas DD conformation and thus to modulate Fas-mediated DISC formation and thus Fas-mediated apoptosis.

Trifluoperazine, an antipsychotic agent, inhibits cancer stem cell growth and overcomes drug resistance of lung cancer

RATIONALE

Cancer stem cell (CSC) theory has drawn much attention, with evidence supporting the contribution of stem cells to tumor initiation, relapse, and therapy resistance.

OBJECTIVES

To screen drugs that target CSCs to improve the current treatment outcome and overcome drug resistance in patients with lung cancer.

METHODS

We used publicly available embryonic stem cell and CSC-associated gene signatures to query the Connectivity Map for potential drugs that can, at least in part, reverse the gene expression profile of CSCs. High scores were noted for several phenothiazine-like antipsychotic drugs, including trifluoperazine. We then treated lung CSCs with different EGFR mutation status with trifluoperazine to examine its anti-CSC properties. Lung CSCs resistant to epidermal growth factor receptor-tyrosine kinase inhibitor or cisplatin were treated with trifluoperazine plus gefitinib or trifluoperazine plus cisplatin. Animal models were used for in vivo validation of the anti-CSC effect and synergistic effect of trifluoperazine with gefitinib.

MEASUREMENTS AND MAIN RESULTS

We demonstrated that trifluoperazine inhibited CSC tumor spheroid formation and down-regulated the expression of CSC markers (CD44/CD133). Trifluoperazine inhibited Wnt/β-catenin signaling in gefitinib-resistant lung cancer spheroids. The combination of trifluoperazine with either gefitinib or cisplatin overcame drug resistance in lung CSCs. Trifluoperazine inhibited the tumor growth and enhanced the inhibitory activity of gefitinib in lung cancer metastatic and orthotopic CSC animal models.

CONCLUSIONS

Using in silico drug screening by Connectivity Map followed by empirical validations, we repurposed an existing phenothiazine-like antipsychotic drug, trifluoperazine, as a potential anti-CSC agent that could overcome epidermal growth factor receptor-tyrosine kinase inhibitor and chemotherapy resistance.