A multi-scale approach to colorectal cancer: from a biochemical- interaction signaling-network level, to multi-cellular dynamics of malignant transformation. Interplay with mutations and onco-protein inhibitor drugs.

This review article is part of a special Current Cancer Drug Targets issue devoted to colorectal cancer and molecularly targeted treatments. In our paper we made an attempt to connect more basic aspects with preclinical, pharmacological / therapeutic and clinical aspects. Reconstruction of a Molecular Interaction Map (MIM) comprising an important part of the G0 - G1 - S cell cycle transition, was a major component of our review. Such a MIM serves also as a convenient / organized database of a large set of important molecular events. The frequency of mutated / altered signaling-proteins indicates the importance of this signaling-network region. We have considered problems at different scale levels. Our MIM works at a biochemical-interaction level. We have also touched the multi-cellular dynamics of normal and aberrant colon crypts. Until recently, dynamic simulations at a biochemical or multi-cellular scale level were considered as a sort of esoteric approach. We tried to convince the reader, also on the basis of a rapidly growing literature, mostly published in high quality journals, that suspicion towards simulations should dissipate, as the limitations and advantages of their application are better appreciated, opening the door to their permanent adoption in everyday research. What is really required is a more interdisciplinary mentality and an interdisciplinary approach. The prize is a level of understanding going beyond mere intuition.

Dynamic simulations of pathways downstream of ERBB-family, including mutations and treatments: concordance with experimental results.

The pathways downstream of ErbB-family proteins are very important in BC, especially when considering treatment with onco-protein inhibitors. We studied and implemented dynamic simulations of four downstream pathways and described the fragment of the signaling network we evaluated as a Molecular Interaction Map. Our simulations, enacted using Ordinary Differential Equations, involved 242 modified species and complexes, 279 reversible reactions and 111 catalytic reactions. Mutations within a single pathway tended to be mutually exclusive; only inhibitors acting at, or downstream (not upstream), of a given mutation were active. A double alteration along two distinct pathways required the inhibition of both pathways. We started an analysis of sensitivity/robustness of our network, and we systematically introduced several individual fluctuations of total concentrations of independent molecular species. Only very few cases showed significant sensitivity. We transduced the ErbB2 over-expressing BC line, BT474, with the HRAS (V12) mutant, then treated it with ErbB-family and phosphorylated MEK (MEKPP) inhibitors, Lapatinib and U0126, respectively. Experimental and simulation results were highly concordant, showing statistical significance for both pathways and for two respective endpoints, i.e. phosphorylated active forms of ERK and Akt, p one tailed = .0072 and = .0022, respectively. Working with a complex 39 basic species signaling network region, this technology facilitates both comprehension and effective, efficient and accurate modeling and data interpretation. Dynamic network simulations we performed proved to be both practical and valuable for a posteriori comprehension of biological networks and signaling, thereby greatly facilitating handling, and thus complete exploitation, of biological data.

Inhibition of a protein-protein interaction between INI1 and c-Myc by small peptidomimetic molecules inspired by Helix-1 of c-Myc: identification of a new target of potential antineoplastic interest.

c-Myc is a transcription modulator proto-oncogene. When overexpressed, it becomes an important contributor to the multi-hit process of malignant transformation. In two earlier papers in this journal (see refs. 19 , 20) we reported that retro-inverso peptidomimetic molecules inspired by the Helix-1 of c-Myc motif could be sequence-specific antiproliferative agents active in the low micromolar range. We also found that our peptides were not opening the four-alpha-helix Myc:Max bundle. Their antiproliferative activity in cancer cell lines needs the presence of side chains projecting outside of the bundle in the corresponding native H1 motif. This observation suggested interference with an external partner. In this study we investigated the INI1:Myc interaction. INI1 is a subunit of the SWI/SNF complex (component of the enhanceosome surrounding Myc:Max heterodimer). The INI1:Myc interaction was confirmed via pull down, ELISA, and fluorescence anisotropy assays. According to the length of INI1 fragments used, we calculated Kds ranging between 1.3x10(-6) and 4.8x10(-7) M. The three different techniques applied showed that the INI1:Myc interaction was also the target of our retro-inverso peptidomimetic molecules, which seem to bind specifically at INI1. A Myc binding, 21aa INI1 fragment (minimum interacting sequence), could inspire the synthesis of a new class of more selective c-Myc inhibitors.

Sequence specific peptidomimetic molecules inhibitors of a protein-protein interaction at the helix 1 level of c-Myc.

Our work is focused in the broad area of strategies and efforts to inhibit protein-protein interactions. The possible strategies in this field are definitely much more varied than in the case of ATP-pocket inhibitors. In our previous work (10), we reported that a retro-inverso (RI) form of Helix1 (H1) of c-Myc, linked to an RI-internalization sequence arising from the third alpha-helix of Antennapedia (Int) was endowed with an antiproliferative and proapoptotic activity toward the cancer cell lines MCF-7 and HCT-116. The activity apparently was dependent upon the presence of the Myc motif. In this work, by ala-scan mapping of the H1 portion of our molecules with D-aa, we found two amino acids necessary for antiproliferative activity: D-Lys in 4 and D-Arg in 5 (numbers refer to L-forms). In the natural hetero-dimer, these two side chains project to the outside of the four alpha-helix bundle. Moreover, we were able to obtain three peptides more active than the original lead. They strongly reduced cell proliferation and survival (RI-Int-VV-H1-E2A,S6A,F8A; RI-Int-VV-H1-S6A,F8A,R11A; RI-Int-VV-H1-S6A,F8A,Q13A): after 8 days at 10 muM total cell number was approximately 1% of the number of cells initially seeded. In these more potent molecules, the ablated side chains project to the inside in the corresponding natural four alpha-helix bundle. In the present work, we also investigated the behavior of our molecules at the biochemical level. Using both a circular dichroism (CD) and a fluorescence anisotropy approach, we noted that side chains projecting at the interior of the four alpha-helix bundle are needed for inducing the partial unfolding of Myc-H2, without an opening of the leucine zipper. Side chains projecting at the outside are not required for this biochemical effect. However, antiproliferative activity had the opposite requirements: side chains projecting at the outside of the bundle were essential, and, on the contrary, ablation of one side chain at a time projecting at the inside increased rather than decreased biological activity. We conclude that our active molecules probably interfere at the level of a protein-protein interaction between Myc-Max and a third protein of the transcription complex. Finally, CD and nuclear magnetic resonance (NMR) data, plus dynamic simulations, suggest a prevalent random coil conformation of the H1 portion of our molecules, at least in diluted solutions. The introduction of a kink (substitution with proline in positions 5 or 7) led to an important reduction of biological activity. We have also synthesized a longer peptido-mimetic molecule (RI-Int-H1-S6A,F8A-loop-H2) with the intent of obtaining a wider zone of interaction and a stronger interference at the level of the higher-order structure (enhanceosome). RI-Int-H1-S6A,F8A-loop-H2 was less active rather than more active in respect to RI-Int-VV-H1-S6A,F8A, apparently because it has a clear bent to form a beta-sheet (CD and NMR data).