Abstract
The human epidermal growth factor receptor (EGFR/ErbB) family consists of four members (ErbB1-4) and belongs to the superfamily of receptor tyrosine kinases (RTKs). The ErbB family members participate in multiple cellular pathways and are the key players in several cancers (brain, breast, lung etc.). Activation of these family members depends on their extracellular domains forming back-to-back hetero/homo dimers. Moreover, dimers are glycosylated, which is a crucial post-translational modification that affects the conformation and function of the protein. Here, molecular modeling and molecular docking are used to comprehensively investigate the dimerization mechanism in glycosylated back-to-back active dimer formation in the entire ErbB receptors for the first time. Results showed that 21 out of 37 clusters of active back-to-back dimers formed by all family members are through heterodimerization. Including; ErbB1-ErbB3/ErbB4, ErbB2-ErbB3/ErbB4 and ErbB3-ErbB4. Ranking ErbB2-ErbB3 as the most stabilized back-to-back dimeric construct. While glycan arrangements favor both homo/hetero dimerization at the dimeric interfaces, it promotes heterodimerization by stabilizing and packing the ligand binding sites of EGFR and ErbB2 respectively. These findings pave the path to future heterodimeric interface/glycan targeting rational anti-cancer drug designs for ErbB receptors.