The interconnected network of pathways downstream of the TGFβ, WNT and EGF-families of receptor ligands play an important role in colorectal cancer pathogenesis.We studied and implemented dynamic simulations of multiple downstream pathways and described the section of the signaling network considered as a Molecular Interaction Map (MIM). Our simulations used Ordinary Differential Equations (ODEs), which involved 447 reactants and their interactions.Starting from an initial “physiologic condition”, the model can be adapted to simulate individual pathologic cancer conditions implementing alterations/mutations in relevant onco-proteins. We verified some salient model predictions using the mutated colorectal cancer lines HCT116 and HT29. We measured the amount of MYC and CCND1 mRNAs and AKT and ERK phosphorylated proteins, in response to individual or combination onco-protein inhibitor treatments. Experimental and simulation results were well correlated. Recent independently published results were also predicted by our model.Even in the presence of an approximate and incomplete signaling network information, a predictive dynamic modeling seems already possible. An important long term road seems to be open and can be pursued further, by incremental steps, toward even larger and better parameterized MIMs. Personalized treatment strategies with rational associations of signaling-proteins inhibitors, could become a realistic goal.

No abstract

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).

Current colorectal cancer (CRC) treatment guidelines are primarily based on clinical features, such as cancer stage and grade. However, outcomes may be improved using molecular treatment guidelines. Potentially useful biomarkers include driver mutations and somatically inherited alterations, signaling proteins (their expression levels and (post) translational modifications), mRNAs, micro-RNAs and long noncoding RNAs. Moving to an integrated system is potentially very relevant. To implement such an integrated system: we focus on an important region of the signaling network, immediately above the G1-S restriction point, and discuss the reconstruction of a Molecular Interaction Map and interrogating it with a dynamic mathematical model. Extensive model pretraining achieved satisfactory, validated, performance. The model helps to propose future target combination priorities, and restricts drastically the number of drugs to be finally tested at a cellular, in vivo, and clinical-trial level. Our model allows for the inclusion of the unique molecular profiles of each individual patient's tumor. While existing clinical guidelines are well established, dynamic modeling may be used for future targeted combination therapies, which may progressively become part of clinical practice within the near future.

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.