Unraveling paralog-specific Notch signaling through thermodynamics of ternary complex formation and transcriptional activation of chimeric receptors.

Notch signaling in humans is mediated by four paralogous receptors that share conserved architectures and possess overlapping, yet non-redundant functions. The receptors share a canonical activation pathway wherein upon extracellular ligand binding, the Notch intracellular domain (NICD) is cleaved from the membrane and translocates to the nucleus where its N-terminal RBP-j-associated molecule (RAM) region and ankyrin repeat (ANK) domain bind transcription factor CSL and recruit co-activator Mastermind-like-1 (MAML1) to activate transcription. However, different paralogs can lead to distinct outcomes. To better understand paralog-specific differences in Notch signaling, we performed a thermodynamic analysis of the Notch transcriptional activation complexes for all four Notch paralogs using isothermal titration calorimetry. Using chimeric constructs, we find that the RAM region is the primary determinant of stability of binary RAMANK:CSL complexes, and that the ANK regions are largely the determinants of MAML1 binding to pre-formed RAMANK:CSL complexes. Free energies of these binding reactions (ΔGRA and ΔGMAML) vary among the four Notch paralogs, although variations for Notch2, 3, and 4 offset in the free energy of the ternary complex (ΔGTC, where ΔGTC = ΔGRA + ΔGMAML). To probe how these affinity differences affect Notch signaling, we performed transcriptional activation assays with the paralogous and chimeric NICDs, and analyzed the results with an independent multiplicative model that quantifies contributions of the paralogous RAM, ANK, and C-terminal regions (CTR) to activation. This analysis shows that transcription activation correlates with ΔGTC, but that activation is further modified by CTR identity in a paralog-specific way.

The importance of input sequence set to consensus-derived proteins and their relationship to reconstructed ancestral proteins.

A protein sequence encodes its energy landscape-all the accessible conformations, energetics, and dynamics. The evolutionary relationship between sequence and landscape can be probed phylogenetically by compiling a multiple sequence alignment of homologous sequences and generating common ancestors via Ancestral Sequence Reconstruction or a consensus protein containing the most common amino acid at each position. Both ancestral and consensus proteins are often more stable than their extant homologs-questioning the differences between them and suggesting that both approaches serve as general methods to engineer thermostability. We used the Ribonuclease H family to compare these approaches and evaluate how the evolutionary relationship of the input sequences affects the properties of the resulting consensus protein. While the consensus protein derived from our full Ribonuclease H sequence alignment is structured and active, it neither shows properties of a well-folded protein nor has enhanced stability. In contrast, the consensus protein derived from a phylogenetically-restricted set of sequences is significantly more stable and cooperatively folded, suggesting that cooperativity may be encoded by different mechanisms in separate clades and lost when too many diverse clades are combined to generate a consensus protein. To explore this, we compared pairwise covariance scores using a Potts formalism as well as higher-order sequence correlations using singular value decomposition (SVD). We find the SVD coordinates of a stable consensus sequence are close to coordinates of the analogous ancestor sequence and its descendants, whereas the unstable consensus sequences are outliers in SVD space.