Search

BACKGROUND AND AIMS: Primary liver cancers include: Hepatocellular carcinoma (HCC), intrahepatic cholangiocarcinoma (CCA) and combined HCC-CCA tumors (cHCC-CCA). It has been suggested, but not unequivocally proven, that hepatic progenitor cells (HPCs) can contribute to hepatocarcinogenesis. We aimed to determine whether HPCs contribute to HCC, cHCC-CCA or both types of tumors.METHOD: To trace progenitor cells during hepatocarcinogenesis, we generated Mdr2-KO mice that harbor an YFP reporter gene driven by the Foxl1 promoter which is expressed specifically in progenitor cells. These mice (Mdr2-KOFoxl1-CRE;RosaYFP) develop chronic inflammation and HCCs by the age of 14-16 months, followed by cHCC-CCA tumors at the age of 18 months, as we have first observed.RESULTS: In this Mdr2-KOFoxl1-CRE;RosaYFP mouse model, liver progenitor cells are the source of cHCC-CCA tumors, but not the source of HCC. Ablating the progenitors, caused reduction of cHCC-CCA tumors but did not affect HCCs. RNA-seq revealed enrichment of the IL6 signaling pathway in cHCC-CCA tumors compared to HCC tumors. ScRNA-seq analysis revealed that IL6 is expressed from immune and parenchymal cells in senescence, and that IL6 is part of the senescence-associated secretory phenotype (SASP). Administration of anti-IL6 Ab to Mdr2-KOFoxl1-CRE;RosaYFP mice, inhibited the development of cHCC-CCA tumors. By blocking IL6 trans-signaling, cHCC-CCA tumors decreased in number and size, indicating that cHCC-CCA is dependent on IL6 trans-signaling. Furthermore, the administration of a senolytic agent inhibited IL6 and the development of cHCC-CCA tumors.CONCLUSION: Our results demonstrate that cHCC-CCA, but not HCC tumors, originate from HPCs, and that IL6, which derives in part from cells in senescence, plays an important role in this process via IL6 trans-signaling. These findings could enhance new therapeutic approaches for cHCC-CCA liver cancer.LAY SUMMARY: Combined hepatocellular carcinoma - cholangiocarcinoma is the third prevalent liver cancer. We show that the source of this tumor is the liver tissue stem cells and that, this tumor type is dependent on an inflammatory signaling of IL6 and can be inhibited by blocking IL6 signaling or using a senolytic agent.

Universal Minicircle Sequence binding proteins (UMSBPs) are CCHC-type zinc-finger proteins that bind the single-stranded G-rich UMS sequence, conserved at the replication origins of minicircles in the kinetoplast DNA, the mitochondrial genome of kinetoplastids. Trypanosoma brucei UMSBP2 has been recently shown to colocalize with telomeres and to play an essential role in chromosome end protection. Here we report that TbUMSBP2 decondenses in vitro DNA molecules, which were condensed by core histones H2B, H4 or linker histone H1. DNA decondensation is mediated via protein-protein interactions between TbUMSBP2 and these histones, independently of its previously described DNA binding activity. Silencing of the TbUMSBP2 gene resulted in a significant decrease in the disassembly of nucleosomes in T. brucei chromatin, a phenotype that could be reverted, by supplementing the knockdown cells with TbUMSBP2. Transcriptome analysis revealed that silencing of TbUMSBP2 affects the expression of multiple genes in T. brucei, with a most significant effect on the upregulation of the subtelomeric variant surface glycoproteins (VSG) genes, which mediate the antigenic variation in African trypanosomes. These observations suggest that UMSBP2 is a chromatin remodeling protein that functions in the regulation of gene expression and plays a role in the control of antigenic variation in T. brucei.

Table S1: List of primers oligos and siRNAs

Several long noncoding RNAs (lncRNA) are abrogated in cancer but their precise contributions to oncogenesis are still emerging. Here we report that the lncRNA MALAT1 is upregulated in hepatocellular carcinoma and acts as a proto-oncogene through Wnt pathway activation and induction of the oncogenic splicing factor SRSF1. Induction of SRSF1 by MALAT1 modulates SRSF1 splicing targets, enhancing the production of antiapoptotic splicing isoforms and activating the mTOR pathway by modulating the alternative splicing of S6K1. Inhibition of SRSF1 expression or mTOR activity abolishes the oncogenic properties of MALAT1, suggesting that SRSF1 induction and mTOR activation are essential for MALAT1-induced transformation. Our results reveal a mechanism by which lncRNA MALAT1 acts as a proto-oncogene in hepatocellular carcinoma, modulating oncogenic alternative splicing through SRSF1 upregulation. Cancer Res; 77(5); 1155–67. ©2016 AACR.

Supplementary information text: Methods and figure legends

Mesenchymal stem cell-derived extracellular vesicles (MSC-EVs) have been employed in the past decade as therapeutic agents in various diseases, including central nervous system (CNS) disorders. We currently aimed to use MSC-EVs as potential treatment for cerebral small vessel disease (CSVD), a complex disorder with a variety of manifestations. MSC-EVs were intranasally administrated to salt-sensitive hypertension prone SBH/y rats that were DOCA-salt loaded (SBH/y-DS), which we have previously shown is a model of CSVD. MSC-EVs accumulated within brain lesion sites of SBH/y-DS. An in vitro model of an inflammatory environment in the brain demonstrated anti-inflammatory properties of MSC-EVs. Following in vivo MSC-EV treatment, gene set enrichment analysis (GSEA) of SBH/y-DS cortices revealed downregulation of immune system response-related gene sets. In addition, MSC-EVs downregulated gene sets related to apoptosis, wound healing and coagulation, and upregulated gene sets associated with synaptic signaling and cognition. While no specific gene was markedly altered upon treatment, the synergistic effect of all gene alternations was sufficient to increase animal survival and improve the neurological state of affected SBH/y-DS rats. Our data suggest MSC-EVs act as microenvironment modulators, through various molecular pathways. We conclude that MSC-EVs may serve as beneficial therapeutic measure for multifactorial disorders, such as CSVD.

The blood-brain barrier (BBB) selectively regulates the entry of molecules into the central nervous system (CNS). A crosstalk between brain microvascular endothelial cells (BMECs) and resident CNS cells promotes the acquisition of functional tight junctions (TJs). Retinoic acid (RA), a key signaling molecule during embryonic development, is used to enhance in vitro BBB models' functional barrier properties. However, its physiological relevance and affected pathways are not fully understood. P450 oxidoreductase (POR) regulates the enzymatic activity of microsomal cytochromes. POR-deficient (PORD) patients display impaired steroid homeostasis and cognitive disabilities. Here, we used both patient-specific POR-deficient and CRISPR-Cas9-mediated POR-depleted induced pluripotent stem cell (iPSC)-derived BMECs (iBMECs) to study the role of POR in the acquisition of functional barrier properties. We demonstrate that POR regulates cellular RA homeostasis and that POR deficiency leads to the accumulation of RA within iBMECs, resulting in the impaired acquisition of TJs and, consequently, to dysfunctional development of barrier properties.

Cytochrome c oxidase (COX), a multimeric protein complex, is the final electron acceptor in the mitochondrial electron transfer chain. Primary COX deficiency, caused by mutations in either mitochondrial DNA or nuclear-encoded genes, is a heterogenous group of mitochondrial diseases with a wide range of presentations, ranging from fatal infantile to subtler. We previously reported a patient with primary COX deficiency due to a pathogenic variant in COX4I1 (encoding the common isoform of COX subunit 4, COX4-1), who presented with bone marrow failure, genomic instability, and short stature, mimicking Fanconi anemia (FA). In the present study, we demonstrated that accumulative DNA damage coincided primarily with proliferative cells in the patient’s fibroblasts and in COX4i1 knockdown cells. Expression analysis implicated a reduction in DNA damage response pathways, which was verified by demonstrating impaired recovery from genotoxic insult and decreased DNA repair. The premature senescence of the COX4-1-deficient cells prevented us from undertaking additional studies; nevertheless, taken together, our results indicate replicative stress and impaired nuclear DNA damage response in COX4-1 deficiency. Interestingly, our in vitro findings recapitulated the patient’s presentation and present status.

Sepsis is an excessive, dysregulated immune response to infection that activates inflammatory and coagulation cascades, which may lead to tissue injury, multiple organ dysfunction syndrome and death. Millions of individuals die annually of sepsis. To date, the only treatment available is antibiotics, drainage of the infection source when possible, and organ support in intensive care units. Numerous previous attempts to develop therapeutic treatments, directed at discreet targets of the sepsis cascade, could not cope with the complex pathophysiology of sepsis and failed. Here we describe a novel treatment, based on empty capsids of SV40 (nanocapsids - NCs). Studies in a severe rat sepsis model showed that pre-treatment by NCs led to a dramatic increase in survival, from zero to 75%. Transcript analyses (RNAseq) demonstrated that the NC treatment is a paradigm shift. The NCs affect multiple facets of biological functions. The affected genes are modified with time, adjusting to the recovery processes. The NCs effect on normal control rats was negligible. The study shows that the NCs are capable of coping with diseases with intricate pathophysiology. Further studies are needed to determine whether when applied after sepsis onset, the NCs still improve outcome.

Biological sex is a fundamental trait influencing development, reproduction, pathogenesis, and medical treatment outcomes. Modeling sex differences is challenging because of the masking effect of genetic variability and the hurdle of differentiating chromosomal versus hormonal effects. In this work we developed a cellular model to study sex differences in humans. Somatic cells from a mosaic Klinefelter syndrome patient were reprogrammed to generate isogenic induced pluripotent stem cell (iPSC) lines with different sex chromosome complements: 47,XXY/46,XX/46,XY/45,X0. Transcriptional analysis of the hiPSCs revealed novel and known genes and pathways that are sexually dimorphic in the pluripotent state and during early neural development. Female hiPSCs more closely resembled the naive pluripotent state than their male counterparts. Moreover, the system enabled differentiation between the contributions of X versus Y chromosome to these differences. Taken together, isogenic hiPSCs present a novel platform for studying sex differences in humans and bear potential to promote gender-specific medicine in the future.

Cytochrome- c- oxidase (COX), a multimeric protein complex, is the final electron acceptor in the mitochondrial electron transfer chain. Primary COX deficiency, caused by mutations in either mitochondrial DNA or nuclear-encoded genes, is a heterogenous group of mitochondrial diseases with a wide range of presentations ranging from fatal infantile to subtler. We previously reported a patient with primary COX deficiency due to a pathogenic variant in COX4I1 (encoding the common isoform of COX subunit 4), who presented with bone-marrow failure, genomic instability and short stature, mimicking Fanconi anemia (FA). In the present study, we demonstrated reduced proliferation and premature senescence in this patient’s fibroblasts and in COX4-1 knockdown cells. Accumulative DNA damage coincided primarily with proliferative cells, indicating replicative stress. Expression analysis implicated DNA damage response which was verified by demonstrating impaired recovery from genotoxic insult and decreased DNA repair. Interestingly, our in-vitro findings recapitulate the patient’s presentation and present status. Thus, we suggest that the premature senescence, resulting from accumulative DNA damage in COX4-1 deficiency is a protective mechanism to avoid malignant transformation in a similar manner to what was reported for FA and other “accelerated aging diseases”.

Cytochrome-c-oxidase (COX) subunit 4 (COX4) plays important roles in the function, assembly and regulation of COX (mitochondrial respiratory complex 4), the terminal electron acceptor of the oxidative phosphorylation (OXPHOS) system. The principal COX4 isoform, COX4-1, is expressed in all tissues, whereas COX4-2 is mainly expressed in the lungs, or under hypoxia and other stress conditions. We have previously described a patient with a COX4-1 defect with a relatively mild presentation compared to other primary COX deficiencies, and hypothesized that this could be the result of a compensatory upregulation of COX4-2. To this end, COX4-1 was downregulated by shRNAs in human foreskin fibroblasts (HFF) and compared to the patient’s cells. COX4-1, COX4-2 and HIF-1α were detected by immunocytochemistry. The mRNA transcripts of both COX4 isoforms and HIF-1 target genes were quantified by RT-qPCR. COX activity and OXPHOS function were measured by enzymatic and oxygen consumption assays, respectively. Pathways were analyzed by CEL-Seq2 and by RT-qPCR. We demonstrated elevated COX4-2 levels in the COX4-1-deficient cells, with a concomitant HIF-1α stabilization, nuclear localization and upregulation of the hypoxia and glycolysis pathways. We suggest that COX4-2 and HIF-1α are upregulated also in normoxia as a compensatory mechanism in COX4-1 deficiency.

ObjectiveCellular senescence limits tumourigenesis by blocking the proliferation of premalignant cells. Additionally, however, senescent cells can exert paracrine effects influencing tumour growth. Senescent cells are present in premalignant pancreatic intraepithelial neoplasia (PanIN) lesions, yet their effects on the disease are poorly characterised. It is currently unknown whether senolytic drugs, aimed at eliminating senescent cells from lesions, could be beneficial in blocking tumour development.DesignTo uncover the functions of senescent cells and their potential contribution to early pancreatic tumourigenesis, we isolated and characterised senescent cells from PanINs formed in a Kras-driven mouse model, and tested the consequences of their targeted elimination through senolytic treatment.ResultsWe found that senescent PanIN cells exert a tumour-promoting effect through expression of a proinflammatory signature that includes high Cox2 levels. Senolytic treatment with the Bcl2-family inhibitor ABT-737 eliminated Cox2-expressing senescent cells, and an intermittent short-duration treatment course dramatically reduced PanIN development and progression to pancreatic ductal adenocarcinoma.ConclusionsThese findings reveal that senescent PanIN cells support tumour growth and progression, and provide a first indication that elimination of senescent cells may be effective as preventive therapy for the progression of precancerous lesions.

Background Mutations in GNE cause a recessive, adult onset myopathy characterized by slowly progressive distal and proximal muscle weakness. Knock-in mice carrying the most frequent mutation in GNE myopathy patients, GneM743T/M743T, usually die few days after birth from severe renal failure, with no muscle phenotype. However, a spontaneous sub-colony remains healthy throughout a normal lifespan without any kidney or muscle pathology. Objective We attempted to decipher the molecular mechanisms behind these phenotypic differences and to determine the mechanisms preventing the kidney and muscles from disease. Methods We analyzed the transcriptome and proteome of kidneys and muscles of sick and healthy GneM743T/M743T mice. Results The sick GneM743T/M743T kidney was characterized by up-regulation of extra-cellular matrix degradation related processes and by down-regulation of oxidative phosphorylation and respiratory electron chain pathway, that was also observed in the asymptomatic muscles. Surprisingly, the healthy kidneys of the GneM743T/M743T mice were characterized by up-regulation of hallmark muscle genes. In addition the asymptomatic muscles of the sick GneM743T/M743T mice showed upregulation of transcription and translation processes. Conclusions Overexpression of muscle physiology genes in healthy GneM743T/M743T mice seems to define the protecting mechanism in these mice. Furthermore, the strong involvement of muscle related genes in kidney may bridge the apparent phenotypic gap between GNE myopathy and the knock-in GneM743T/M743T mouse model and provide new directions in the study of GNE function in health and disease.

Phagocytosis level is an indicator of cancer aggressiveness since both are related to mechanical deformability of the cells., The malignancy potential is correlated with the mechanical deformability of the cancer cells. However, mechanical tests for clinical applications are limited. We present here a Triangular Correlation (TrC) between cell deformability, phagocytic capacity, and cancer aggressiveness, suggesting that phagocytic measurements can be a mechanical surrogate marker of malignancy. The TrC was proved in human prostate cancer cells with different malignancy potential, and in human bladder cancer and melanoma cells that were sorted into subpopulations based solely on their phagocytic capacity. The more phagocytic subpopulations showed elevated aggressiveness ex vivo and in vivo. The uptake potential was preserved, and differences in gene expression and in epigenetic signature were detected. In all cases, enhanced phagocytic and aggressiveness phenotypes were correlated with greater cell deformability and predicted by a computational model. Our multidisciplinary study provides the proof of concept that phagocytic measurements can be applied for cancer diagnostics and precision medicine.

Aging is associated with progressive decline in cell function and with increased damage to macromolecular components. DNA damage, in the form of double-strand breaks (DSBs), increases with age and in turn, contributes to the aging process and age-related diseases. DNA strand breaks triggers a set of highly orchestrated signaling events known as the DNA damage response (DDR), which coordinates DNA repair. However, whether the accumulation of DNA damage with age is a result of decreased repair capacity, remains to be determined. In our study we showed that with age there is a decline in the resolution of foci containing γH2AX and pKAP-1 in diethylnitrosamine (DEN)-treated mouse livers, already evident at a remarkably early age of 6-months. Considerable age-dependent differences in global gene expression profiles in mice livers after exposure to DEN, further affirmed these age related differences in the response to DNA damage. Functional analysis identified p53 as the most overrepresented pathway that is specifically enhanced and prolonged in 6-month-old mice. Collectively, our results demonstrated an early decline in DNA damage repair that precedes ‘old age’, suggesting this may be a driving force contributing to the aging process rather than a phenotypic consequence of old age.

Several long noncoding RNAs (lncRNA) are abrogated in cancer but their precise contributions to oncogenesis are still emerging. Here we report that the lncRNA MALAT1 is upregulated in hepatocellular carcinoma and acts as a proto-oncogene through Wnt pathway activation and induction of the oncogenic splicing factor SRSF1. Induction of SRSF1 by MALAT1 modulates SRSF1 splicing targets, enhancing the production of antiapoptotic splicing isoforms and activating the mTOR pathway by modulating the alternative splicing of S6K1. Inhibition of SRSF1 expression or mTOR activity abolishes the oncogenic properties of MALAT1, suggesting that SRSF1 induction and mTOR activation are essential for MALAT1-induced transformation. Our results reveal a mechanism by which lncRNA MALAT1 acts as a proto-oncogene in hepatocellular carcinoma, modulating oncogenic alternative splicing through SRSF1 upregulation. Cancer Res; 77(5); 1155–67. ©2016 AACR.

Peptide arrays are a powerful tool for characterizing protein interactions and for identifying and optimizing peptide ligands. Here we demonstrate the use of peptide arrays for performing detailed SAR studies of lead peptide inhibitors of the interaction between the HIV integrase and Rev proteins. The integration of viral DNA into the genome of the host cell, mediated by the viral integrase (IN) enzyme, is a crucial step in the HIV-1 replication cycle. We have recently found that IN activity is regulated by interactions with the HIV-1 Rev protein and identified three lead peptides derived from the Rev-binding interface in IN. Due to their ability to promote dissociation of the Rev–IN complex in HIV infected cells, these peptides caused IN activation leading to multi-integration, genomic instability and specific eradication of such infected cells. Here we explored the mechanism of action of these three IN-derived peptides as the basis for developing improved anti HIV-1 leads. Using peptide array screening, we found that the IN derived peptides bound IN and Rev in a similar pattern. The Rev-binding sites in IN also mediate IN oligomerization while the IN-binding sites in Rev are also involved in Rev oligomerization. A structural homology was found between the oligomerization domains of Rev and IN residues 171–188 and 211–220. We performed SAR studies of the lead inhibitory peptides using a peptide array containing truncated peptides, alanine scan, D-amino acid scan and N-methylated amino acid scan. We screened IN and Rev for binding this array of modified IN-derived peptides. The screening results showed that C-terminal positively charged residues were essential for the interaction of the IN 118–128 and IN 174–188 peptides with both Rev and IN. The peptides could be shortened and modified without loss of binding to IN and Rev. This provides a basis for the future development of shorter peptides with better pharmacological properties that inhibit the Rev–IN interaction. We conclude that peptide arrays are excellent tools to perform detailed SAR binding studies in one short efficient experiment. The SAR study by the peptide array method is a powerful tool for developing improved inhibitors based on a lead peptide sequence.

Numerous genes and molecular pathways are implicated in neurodegenerative proteinopathies, but their inter-relationships are poorly understood. We systematically mapped molecular pathways underlying the toxicity of alpha-synuclein (α-syn), a protein central to Parkinson's disease. Genome-wide screens in yeast identified 332 genes that impact α-syn toxicity. To “humanize” this molecular network, we developed a computational method, TransposeNet. This integrates a Steiner prize-collecting approach with homology assignment through sequence, structure, and interaction topology. TransposeNet linked α-syn to multiple parkinsonism genes and druggable targets through perturbed protein trafficking and ER quality control as well as mRNA metabolism and translation. A calcium signaling hub linked these processes to perturbed mitochondrial quality control and function, metal ion transport, transcriptional regulation, and signal transduction. Parkinsonism gene interaction profiles spatially opposed in the network (ATP13A2/PARK9 and VPS35/PARK17) were highly distinct, and network relationships for specific genes (LRRK2/PARK8, ATXN2, and EIF4G1/PARK18) were confirmed in patient induced pluripotent stem cell (iPSC)-derived neurons. This cross-species platform connected diverse neurodegenerative genes to proteinopathy through specific mechanisms and may facilitate patient stratification for targeted therapy. Keywords: alpha-synuclein; iPS cell; Parkinson’s disease; stem cell; mRNA translation; RNA-binding protein; LRRK2; VPS35; vesicle trafficking; yeast

The HIV-1 Rev and integrase (IN) proteins control important functions in the viral life cycle. We have recently discovered that the interaction between these proteins results in inhibition of IN enzymatic activity. Peptides derived from the Rev and IN binding interfaces have a profound effect on IN catalytic activity: Peptides derived from Rev inhibit IN, while peptides derived from IN stimulate IN activity by inhibiting the Rev-IN interaction. This inhibition leads to multi integration, genomic instability and specific death of virus-infected cells. Here we used protein docking combined with refinement and energy function ranking to suggest a structural model for the Rev-IN complex. Our results indicate that a Rev monomer binds IN at two sites that match our experimental binding data: (1) IN residues 66-80 and 118-128; (2) IN residues 174-188. According to our model, IN binds Rev and its cellular cofactor, lens epithelium derived growth factor (LEDGF), through overlapping interfaces. This supports previous observations that IN is regulated by a tight interplay between Rev and LEDGF. Rev may bind either the IN dimer or tetramer. Accordingly, Rev is suggested to inhibit IN by two possible mechanisms: (i) shifting the oligomerization equilibrium of IN from an active dimer to an inactive tetramer; (ii) displacing LEDGF from IN, resulting in inhibition of IN binding to the viral DNA. Our model is expected to contribute to the development of lead compounds that inhibit the Rev-IN interaction and thus lead to multi-integration of viral cDNA and consequently to apoptosis of HIV-1 infected cells.

The STIL protein participates in mitosis and malignant transformation by regulating centrosomal duplication. Using biophysical methods we studied the structure and interactions of STIL. We revealed that its central domain is intrinsically disordered and mediates protein-protein interactions of STIL. The intrinsic disorder may provide STIL with the conformational flexibility required for its multitude binding.

Several peptides that specifically bind the HIV-1 integrase (IN) and either inhibit or stimulate its enzymatic activity were developed in our laboratories. Kinetic studies using 3'-end processing and strand-transfer assays were performed to study the mode of action of these peptides. The effects of the various peptides on the interaction between IN and its substrate DNA were also studied by fluorescence anisotropy. On the basis of our results, we divided these IN-interacting peptides into three groups: (a) IN-inhibitory peptides, whose binding to IN decrease its affinity for the substrate DNA - these peptides increased the K(m) of the IN-DNA interaction, and were thus inhibitory; (b) peptides that slightly increased the K(m) of the IN-DNA interaction, but in addition modified the V(max) and K(cat) values of the IN, and thus stimulated or inhibited IN activity, respectively; and (c) peptides that bound IN but had no effect on its enzymatic activity. We elucidated the approximate binding sites of the peptides in the structure of IN, providing structural insights into their mechanism of action. The IN-stimulating peptide bound IN in several specific sites that did not bind any of the inhibitory peptides. This may account for its unique activity.

We used computational methods to study the interaction between two key proteins in apoptosis regulation: the transcription factor NF-κ-B (NFκB) and the proapoptotic protein ASPP2. The C-terminus of ASPP2 contains ankyrin repeats and SH3 domains (ASPP2ANK-SH3) that mediate interactions with numerous apoptosis-related proteins, including the p65 subunit of NFκB (NFκBp65). Using peptide-based methods, we have recently identified the interaction sites between NFκBp65 and ASPP2ANK-SH3 (Rotem et al., J Biol Chem 283, 18990–18999). Here we conducted a computational study of protein docking and molecular dynamics to obtain a structural model of the complex between the full length proteins and propose a mechanism for the interaction. We found that ASPP2ANK-SH3 binds two sites in NFκBp65, at residues 236–253 and 293–313 that contain the nuclear localization signal (NLS). These sites also mediate the binding of NFκB to its natural inhibitor IκB, which also contains ankyrin repeats. Alignment of the ankyrin repeats of ASPP2ANK-SH3 and IκB revealed that both proteins share highly similar interfaces at their binding sites to NFκB. Protein docking of ASPP2ANK-SH3 and NFκBp65, as well as molecular dynamics simulations of the proteins, provided structural models of the complex that are energetically similar to the NFκB-IκB determined structure. Our results show that ASPP2ANK-SH3 binds NFκBp65 in a similar manner to its natural inhibitor IκB, suggesting a possible novel role for ASPP2 as an NFκB inhibitor. Proteins 2009. © 2009 Wiley-Liss, Inc.

We have characterized the molecular basis of the interaction between ASPP2 and Bcl-2, which are key proteins in the apoptotic pathway. The C-terminal ankyrin repeats and SH3 domain of ASPP2 (ASPP2 Ank-SH3 ) mediate its interactions with the antiapoptotic protein Bcl-2. We used biophysical and computational methods to identify the interaction sites of Bcl-2 and its homologues with ASPP2. Using peptide array screening, we found that ASPP2 Ank-SH3 binds two homologous sites in all three Bcl proteins tested: ( i ) the conserved BH4 motif, and ( ii ) a binding site for proapoptotic regulators. Quantitative binding studies revealed that binding of ASPP2 Ank-SH3 to the Bcl-2 family members is selective at two levels: ( i ) interaction with Bcl-2-derived peptides is the tightest compared to peptides from the other family members, and ( ii ) within Bcl-2, binding of ASPP2 Ank-SH3 to the BH4 domain is tightest. Sequence alignment of the ASPP2-binding peptides combined with binding studies of mutated peptides revealed that two nonconserved positions where only Bcl-2 contains positively charged residues account for its tighter binding. The experimental binding results served as a basis for docking analysis, by which we modeled the complexes of ASPP2 Ank-SH3 with the full-length Bcl proteins. Using peptide arrays and quantitative binding studies, we found that Bcl-2 binds three loops in ASPP2 Ank-SH3 with similar affinity, in agreement with our predicted model. Based on our results, we propose a mechanism in which ASPP2 induces apoptosis by inhibiting functional sites of the antiapoptotic Bcl-2 proteins.

The unique properties of fullerenes have raised the interest of using them for biomedical applications. Within this framework, the interactions of fullerenes with proteins have been an exciting research target, yet little is known about how native proteins can bind fullerenes, and what is the nature of these interactions. Moreover, though some proteins have been shown to interact with fullerenes, up to date, no crystal structure of such complexes was obtained. Here we report docking studies aimed at examining the interactions of fullerene in two forms (C60 nonsubstituted fullerene and carboxyfullerene) with four proteins that are known to bind fullerene derivatives: HIV protease, fullerene-specific antibody, human serum albumin, and bovine serum albumin. Our work provides docking models with detailed binding pockets information, which closely match available experimental data. We further compare the predicted binding sites using a novel multiple binding site alignment method. A high similarity between the physicochemical properties and surface geometry was found for fullerene's binding sites of HIV protease and the human and bovine serum albumins.

The preparation and characterization of the stable human serum albumin (HSA)-C3 isomer of tris-malonic acid [C60]fullerene complex is reported. Other than the anti-fullerene antibody, a stable protein-fullerene complex with a native protein has never been observed. This study may provide valuable answers to the growing concern regarding the effects of carbonaceous nanomaterials on human health on one hand and, on the other, may lead to the development of novel antioxidant therapeutic agents, radiopharmaceuticals, and components for bioelectronic devices.

We present a novel method for multiple alignment of protein structures and detection of structural motifs. To date, only a few methods are available for addressing this task. Most of them are based on a series of pairwise comparisons. In contrast, MASS (Multiple Alignment by Secondary Structures) considers all the given structures at the same time. Exploiting the secondary structure representation aids in filtering out noisy results and in making the method highly efficient and robust. MASS disregards the sequence order of the secondary structure elements. Thus, it can find non-sequential and even non-topological structural motifs. An important novel feature of MASS is subset alignment detection: It does not require that all the input molecules be aligned. Rather, MASS is capable of detecting structural motifs shared only by a subset of the molecules. Given its high efficiency and capability of detecting subset alignments, MASS is suitable for a broad range of challenging applications: It can handle large-scale protein ensembles (on the order of tens) that may be heterogeneous, noisy, topologically unrelated and contain structures of low resolution. Keywords: Structural Bioinformatics, Subset alignment, Non-sequential alignment, Non-topological motif, Supersecondary structural motif Availability: MASS program and more information about it are available on http://bioinfo3d.cs.tau.ac.il/MASS Contact: oranit@post.tau.ac.il *To whom correspondence should be addressed.

Current data suggest that globular domains may form amyloids via different mechanisms. Nevertheless, there are indications that the initiation of the process takes place invariably in the less stable segments of a protein domain. We have studied the sequence and structural conservation of beta(2)-microglobulin that deposits into fibrils in dialysis-related amyloidosis. The dataset includes 51 high-resolution non-redundant structures of the antibody constant domain-like proteins (C1) and 132 related sequences. We describe a set of 30 conserved residues. Among them, 23 are conserved structurally, 16 are conserved sequentially and nine are conserved both sequentially and structurally. Strands A (12-18), G (91-95) and D (45-55) are the less conserved and stable segments of the domain, while strands B (22-28), C (36-41), E (62-70) and F (78-83) are the conserved and stable segments. We find that the conserved residues form a cluster with a network of interactions. The observed pattern of conservation is consistent with experimental data including H/D exchange, urea denaturation and limited proteolysis that suggest that strands A and G do not participate in the amyloid fibril. Additionally, the low conservation of strand D is consistent with the observation that this strand may acquire different conformations as seen in crystal structures of bound and isolated beta(2)-microglobulin. We used a docking technique to suggest a model for a fibril via stacking of beta(2)-microglobulin monomers. Our analysis suggests that the favored monomer building block for fibril elongation is the conformation of the isolated beta(2)-microglobulin, without the beta-bulge on strand D and without strands A and G participating in the fibril beta-sheet structure. This monomer retains all the conserved residues and their network of interactions, increasing the likelihood of its existence in solution. The inter-strand interaction between the two (monomer) building blocks forms a new continuous beta-sheet such that addition of monomers results in a fibril model that has the characteristic cross-beta structure.

The mechanism through which globular proteins transform into amyloid fibrils is still not understood. Here we analyze the structure and sequence conservation to assess the differential stability of segments from two structurally related protein families: the amyloidogenic gelsolin-like and its structurally related cofilin-like. The two families belong to the actin depolymerizing proteins, with a central β-sheet stacked between 2 and 4 α-helices. Although sequentially remote, the two families share regions of high and low conservation and stability. Our results show a highly conserved hydrophobic and aromatic cluster, located at a central buried β-hairpin. The geometry of the aromatic residues with respect to each other is strictly conserved, suggesting involvement in strand registering and β-sheet stabilization. Consistent with experiment, we find a region of weak conservation and stability at one of the exposed β-strands (strand B in the gelsolin-like family). This region was recently found to be affected by a point mutation-mediated destabilization of the human gelsolin domain 2, which facilitates the first proteolytic event in the formation of the amyloidogenic fragment. Thus, both experimental and computational conservation analyses suggest that this unstable region may constitute a first step in amyloid formation. Our analysis uses a recently developed multiple-structure comparison algorithm in which molecules are aligned simultaneously. Proteins 2003;51:266–282. © 2003 Wiley-Liss, Inc.

Protein–protein interactions (PPIs) govern all aspects of cell function and, as such, are a major target for research and therapeutic intervention. A major rate-limiting step in PPI research is the expression and purification of full-length proteins. The use of peptides to study PPIs significantly facilitates the structural and biophysical characterization of PPIs as well as the effort to develop drugs to control PPIs. Here we describe examples for the use of peptides to study PPI and some of the important experimental methods that are used in the field. Peptides have proved to be excellent tools to study PPIs and have been contributing both for understanding mechanisms of PPIs as well as for drug design for PPI modulation.

Several peptides that specifically bind the HIV-1 integrase (IN) and either inhibit or stimulate its enzymatic activity were developed in our laboratories. Kinetic studies using 3'-end processing and strand-transfer assays were performed to study the mode of action of these peptides. The effects of the various peptides on the interaction between IN and its substrate DNA were also studied by fluorescence anisotropy. On the basis of our results, we divided these IN-interacting peptides into three groups: (a) IN-inhibitory peptides, whose binding to IN decrease its affinity for the substrate DNA - these peptides increased the K(m) of the IN-DNA interaction, and were thus inhibitory; (b) peptides that slightly increased the K(m) of the IN-DNA interaction, but in addition modified the V(max) and K(cat) values of the IN, and thus stimulated or inhibited IN activity, respectively; and (c) peptides that bound IN but had no effect on its enzymatic activity. We elucidated the approximate binding sites of the peptides in the structure of IN, providing structural insights into their mechanism of action. The IN-stimulating peptide bound IN in several specific sites that did not bind any of the inhibitory peptides. This may account for its unique activity.

The ASPP proteins are apoptosis regulators: ASPP1 and ASPP2 promote, while iASPP inhibits, apoptosis. The mechanism by which these different outcomes are achieved is still unknown. The C-terminal ankyrin repeats and SH3 domain (ANK-SH3) mediate the interactions of the ASPP proteins with major apoptosis regulators such as p53, Bcl-2, and NFκB. The structure of the complex between ASPP2ANK-SH3 and the core domain of p53 (p53CD) was previously determined. We have recently characterized the individual interactions of ASPP2ANK-SH3 with Bcl-2 and NFκB, as well as a regulatory intramolecular interaction with the proline rich domain of ASPP2. Here we compared the ASPP interactions at two levels: ASPP2ANK-SH3 with different proteins, and different ASPP family members with each protein partner. We found that the binding sites of ASPP2 to p53CD, Bcl-2, and NFκB are different, yet lie on the same face of ASPP2ANK-SH3. The intramolecular binding site to the proline rich domain overlaps the three intermolecular binding sites. To reveal the basis of functional diversity in the ASPP family, we compared their protein-binding domains. A subset of surface-exposed residues differentiates ASPP1 and ASPP2 from iASPP: ASPP1/2 are more negatively charged in specific residues that contact positively charged residues of p53CD, Bcl-2, and NFκB. We also found a gain of positive charge at the non-protein binding face of ASPP1/2, suggesting a role in electrostatic direction towards the negatively charged protein binding face. The electrostatic differences in binding interfaces between the ASPP proteins may be one of the causes for their different function. Copyright © 2010 John Wiley & Sons, Ltd.

We used computational methods to study the interaction between two key proteins in apoptosis regulation: the transcription factor NF-kappa-B (NFkappaB) and the proapoptotic protein ASPP2. The C-terminus of ASPP2 contains ankyrin repeats and SH3 domains (ASPP2(ANK-SH3)) that mediate interactions with numerous apoptosis-related proteins, including the p65 subunit of NFkappaB (NFkappaB(p65)). Using peptide-based methods, we have recently identified the interaction sites between NFkappaB(p65) and ASPP2(ANK-SH3) (Rotem et al., J Biol Chem 283, 18990-18999). Here we conducted a computational study of protein docking and molecular dynamics to obtain a structural model of the complex between the full length proteins and propose a mechanism for the interaction. We found that ASPP2(ANK-SH3) binds two sites in NFkappaB(p65), at residues 236-253 and 293-313 that contain the nuclear localization signal (NLS). These sites also mediate the binding of NFkappaB to its natural inhibitor IkappaB, which also contains ankyrin repeats. Alignment of the ankyrin repeats of ASPP2(ANK-SH3) and IkappaB revealed that both proteins share highly similar interfaces at their binding sites to NFkappaB. Protein docking of ASPP2(ANK-SH3) and NFkappaB(p65), as well as molecular dynamics simulations of the proteins, provided structural models of the complex that are energetically similar to the NFkappaB-IkappaB determined structure. Our results show that ASPP2(ANK-SH3) binds NFkappaB(p65) in a similar manner to its natural inhibitor IkappaB, suggesting a possible novel role for ASPP2 as an NFkappaB inhibitor.

The human immunodeficiency virus type-1 (HIV-1) Vif protein neutralizes the cellular defense mechanism against the virus. The C-terminal domain of Vif (CTD, residues 141-192) mediates many of its interactions. Full-length Vif is difficult to purify in large amounts, hence the only available structure of Vif is of residues 140-155 within the ElonginBC complex. Other structural information, derived from modeling and indirect experiments, indicates that the Vif CTD may be unstructured. Here, we chemically synthesized the Vif CTD using pseudo-proline-building blocks, studied its solution structure in the unbound state using biophysical techniques and found that it is unstructured under physiological conditions. The circular dichroism (CD) spectrum of Vif CTD showed a pattern of random coil with residual helical structure. The (15)N-HSQC nuclear magnetic resonance (NMR) spectrum was characteristic of natively unfolded peptides. Vif CTD eluted from an analytical gel filtration column earlier than expected, indicating an extended conformation. Disorder predictions found the CTD to be unstructured, in agreement with our experimental results. CD experiments showed that Vif CTD underwent a conformational change upon interacting with membrane-mimicking DPC micelles, but not upon binding to a peptide derived from its binding region in ElonginC. Our results provide direct evidence for the unfolded structure of the free Vif CTD and indicate that it may gain structure upon binding its natural ligands.

ASPP2 is a pro-apoptotic protein that stimulates the p53-mediated apoptotic response. The C terminus of ASPP2 contains ankyrin (Ank) repeats and a SH3 domain, which mediate its interactions with numerous partner proteins such as p53, NFkappaB, and Bcl-2. It also contains a proline-rich domain (ASPP2 Pro), whose structure and function are unclear. Here we used biophysical and biochemical methods to study the structure and the interactions of ASPP2 Pro, to gain insight into its biological role. We show, using biophysical and computational methods, that the ASPP2 Pro domain is natively unfolded. We found that the ASPP2 Pro domain interacts with the ASPP2 Ank-SH3 domains, and mapped the interaction sites in both domains. Using a combination of peptide array screening, biophysical and biochemical techniques, we found that ASPP2 Ank-SH3, but not ASPP2 Pro, mediates interactions of ASPP2 with peptides derived from its partner proteins. ASPP2 Pro-Ank-SH3 bound a peptide derived from its partner protein NFkappaB weaker than ASPP2 Ank-SH3 bound this peptide. This suggested that the presence of the proline-rich domain inhibited the interactions mediated by the Ank-SH3 domains. Furthermore, a peptide from ASPP2 Pro competed with a peptide derived from NFkappaB on binding to ASPP2 Ank-SH3. Based on our results, we propose a model in which the interaction between the ASPP2 domains regulates the intermolecular interactions of ASPP2 with its partner proteins.

Protein folding and protein binding are similar processes. In both, structural units combinatorially associate with each other. In the case of folding, we mostly handle relatively small units, building blocks or domains, that are covalently linked. In the case of multi-molecular binding, the subunits are relatively large and are associated only by non-covalent bonds. Experimentally, the difficulty in the determination of the structures of such large assemblies increases with the complex size and the number of components it contains. Computationally, the prediction of the structures of multi-molecular complexes has largely not been addressed, probably owing to the magnitude of the combinatorial complexity of the problem. Current docking algorithms mostly target prediction of pairwise interactions. Here our goal is to predict the structures of multi-unit associations, whether these are chain-connected as in protein folding, or separate disjoint molecules in the assemblies. We assume that the structures of the single units are known, either through experimental determination or modeling. Our aim is to combinatorially assemble these units to predict their structure. To address this problem we have developed CombDock. CombDock is a combinatorial docking algorithm for the structural units assembly problem. Below, we briefly describe the algorithm and present examples of its various applications to folding and to multi-molecular assemblies. To test the robustness of the algorithm, we use inaccurate models of the structural units, derived either from crystal structures of unbound molecules or from modeling of the target sequences. The algorithm has been able to predict near-native arrangements of the input structural units in almost all of the cases, suggesting that a combinatorial approach can overcome the imperfect shape complementarity caused by the inaccuracy of the models. In addition, we further show that through a combinatorial docking strategy it is possible to enhance the predictions of pairwise interactions involved in a multi-molecular assembly.

Obtaining atomic resolution structural models of amyloid fibrils is currently impossible, yet crucial for our understanding of the amyloid mechanism. Different pathways in the transformation of a native globular domain to an amyloid fibril invariably involve domain destabilization. Hence, locating the unstable segments of a domain is important for understanding its amyloidogenic transformation and possibly control it. Since relative conservation is suggested to relate to local stability [H. Benyamini, K. Gunasekaran, H. Wolfson, R. Nussinov, Conservation and amyloid formation: a study of the gelsolin-like family, Proteins 51 (2003) 266–282. [24]], we performed an extensive, sequence and structure conservation analysis of the β 2 -microglobulin (β 2 -m) domain. Our dataset include 51 high resolution structures belonging to the “C1 set domain” family and 132 clustered PSI-BLAST search results. Segments of the β 2 -m domain corresponding to strands A (residues 12–18), D (45–55) and G (91–95) were found to be less conserved and stable, while the central strands B (residues 22–28), C (36–41), E (62–70) and F (78–83) were found conserved and stable. Our findings are supported by accumulating observations from various experimental methods, including urea denaturation, limited proteolysis, H/D exchange and structure determination by both NMR and X-ray crystallography. We used our conservation findings together with experimental literature information to suggest a structural model for the polymerized unit of β 2 -m. Pairwise protein docking and subsequent monomer stacking in the same manner suggest a fibril model consistent with the cross-β structure.

The majority of proteins function when associated in multimolecular assemblies. Yet, prediction of the structures of multimolecular complexes has largely not been addressed, probably due to the magnitude of the combinatorial complexity of the problem. Docking applications have traditionally been used to predict pairwise interactions between molecules. We have developed an algorithm that extends the application of docking to multimolecular assemblies. We apply it to predict quaternary structures of both oligomers and multi-protein complexes. The algorithm predicted well a near-native arrangement of the input subunits for all cases in our data set, where the number of the subunits of the different target complexes varied from three to ten. In order to simulate a more realistic scenario, unbound cases were tested. In these cases the input conformations of the subunits are either unbound conformations of the subunits or a model obtained by a homology modeling technique. The successful predictions of the unbound cases, where the input conformations of the subunits are different from their conformations within the target complex, suggest that the algorithm is robust. We expect that this type of algorithm should be particularly useful to predict the structures of large macromolecular assemblies, which are difficult to solve by experimental structure determination.

We present MASS (Multiple Alignment by Secondary Structures), a novel highly efficient method for structural alignment of multiple protein molecules and detection of common structural motifs. MASS is based on a two-level alignment, using both secondary structure and atomic representation. Utilizing secondary structure information aids in filtering out noisy solutions and achieves efficiency and robustness. Currently, only a few methods are available for addressing the multiple structural alignment task. In addition to using secondary structure information, the advantage of MASS as compared to these methods is that it is a combination of several important characteristics: (1) While most existing methods are based on series of pairwise comparisons, and thus might miss optimal global solutions, MASS is truly multiple, considering all the molecules simultaneously; (2) MASS is sequence order-independent and thus capable of detecting nontopological structural motifs; (3) MASS is able to detect not only structural motifs, shared by all input molecules, but also motifs shared only by subsets of the molecules. Here, we show the application of MASS to various protein ensembles. We demonstrate its ability to handle a large number (order of tens) of molecules, to detect nontopological motifs and to find biologically meaningful alignments within nonpredefined subsets of the input. In particular, we show how by using conserved structural motifs, one can guide protein–protein docking, which is a notoriously difficult problem. MASS is freely available at http://bioinfo3d.cs.tau.ac.il/MASS/.