Emergence of MUC1 in Mammals for Adaptation of Barrier Epithelia

The mucin 1 (MUC1) gene was discovered based on its overexpression in human breast cancers. Subsequent work demonstrated that MUC1 is aberrantly expressed in cancers originating from other diverse organs, including skin and immune cells. These findings supported a role for MUC1 in the adaptation of barrier tissues to infection and environmental stress. Of fundamental importance for this evolutionary adaptation was inclusion of a SEA domain, which catalyzes autoproteolysis of the MUC1 protein and formation of a non-covalent heterodimeric complex. The resulting MUC1 heterodimer is poised at the apical cell membrane to respond to loss of homeostasis. Disruption of the complex releases the MUC1 N-terminal (MUC1-N) subunit into a protective mucous gel. Conversely, the transmembrane C-terminal (MUC1-C) subunit activates a program of lineage plasticity, epigenetic reprogramming and repair. This MUC1-C-activated program apparently evolved for barrier tissues to mount self-regulating proliferative, inflammatory and remodeling responses associated with wound healing. Emerging evidence indicates that MUC1-C underpins inflammatory adaptation of tissue stem cells and immune cells in the barrier niche. This review focuses on how prolonged activation of MUC1-C by chronic inflammation in these niches promotes the cancer stem cell (CSC) state by establishing auto-inductive nodes that drive self-renewal and tumorigenicity.

Computational Alanine Scanning Reveals Common Features of TCR/pMHC Recognition in HLA-DQ8-Associated Celiac Disease

In HLA-DQ8-associated celiac disease, Gliadin-γ1 or Gliadin-α1 peptide is presented to the cell surface and recognized by several types of T-cell receptor (TCR), but it is still unclear how the TCR, peptide, and the major histocompatibility complex (MHC) act together to trigger celiac disease. For now, most of the analysis is based on static crystal structures. And the detailed information about these structures based on energetic interaction is still lacking. Here, we took four types of celiac disease-related MHC-peptide-TCR structures from three patients to perform computational alanine scanning calculations using the molecular mechanics generalized born surface area (MM/GBSA) approach combined with a recently developed interaction entropy (IE) method to identify the key residues on TCR, peptide, and MHC. Our study aims to shed some light on the interaction mechanism of this complex protein interaction system. Based on detailed computational analysis and mutational calculations, important binding interactions in these triple-interaction complexes are analyzed, and critical residues responsible for TCR/pMHC recognition pattern in HLA-DQ8-associated celiac disease are presented. These detailed analysis and computational result should help shed light on our understanding of the celiac disease and the development of the medical treatment.

Molecular basis of SMAC-XIAP binding and the effect of electrostatic polarization

X-chromosome-linked inhibitor of apoptosis (XIAP) inhibits cell apoptosis. Overexpression of XIAP is widely found in human cancers. Second mitochondria-derived activator of caspase (SMAC) protein inhibits XIAP through binding with Baculovirus Inhibitor of apoptosis protein Repeat (BIR) 3 or BIR2 domain of XIAP. In this study, molecular dynamics (MD) simulations and the alanine scanning calculations by MM-GBSA_IE method were used to investigate the protein-peptide interaction between BIR3 and BIR2 domains of XIAP and SMAC peptide. Energetic contribution of each binding residue is calculated and hotspots on both XIAP and SMAC were identified using computational alanine scanning with interaction entropy method. We found that electrostatic polarization is important in stabilizing the protein-protein complex structure in MD simulation. By using polarized protein-specific charges, much better agreement with experimental result is obtained for calculated binding free energies compared to those using standard (nonpolarizable) AMBER force field. In particular, excellent correlation between calculated binding free energies in alanine scanning with mutational experimental data was obtained for BIR3/SMAC binding.Communicated by Ramaswamy H. Sarma.

Exploring the mechanism how Marburg virus VP35 recognizes and binds dsRNA by molecular dynamics simulations and free energy calculations

Filoviruses often cause terrible infectious disease which has not been successfully dealt with pharmacologically. All filoviruses encode a unique protein termed VP35 which can mask doubled-stranded RNA to deactivate interferon. The interface of VP35-dsRNA would be a feasible target for structure-based antiviral agent design. To explore the essence of VP35-dsRNA interaction, molecular dynamics simulation combined with MM-GBSA calculations were performed on Marburg virus VP35-dsRNA complex and several mutational complexes. The energetic analysis indicates that nonpolar interactions provide the main driving force for the binding process. Although the intermolecular electrostatic interactions play important roles in VP35-dsRNA interaction, the whole polar interactions are unfavorable for binding which result in a low binding affinity. Compared with wild type VP35, the studied mutants F228A, R271A, and K298A have obviously reduced binding free energies with dsRNA reflecting in the reduction of polar or nonpolar interactions. The results also indicate that the loss of binding affinity for one dsRNA strand would abolish the total binding affinity. Three important residues Arg271, Arg294, and Lys298 which makes the largest contribution for binding in VP35 lose their binding affinity significantly in mutants. The uncovering of VP35-dsRNA recognition mechanism will provide some insights for development of antiviral drug.

Probing Origin of Binding Difference of inhibitors to MDM2 and MDMX by Polarizable Molecular Dynamics Simulation and QM/MM-GBSA Calculation

Binding abilities of current inhibitors to MDMX are weaker than to MDM2. Polarizable molecular dynamics simulations (MD) followed by Quantum mechanics/molecular mechanics generalized Born surface area (QM//MM-GBSA) calculations were performed to investigate the binding difference of inhibitors to MDM2 and MDMX. The predicted binding free energies not only agree well with the experimental results, but also show that the decrease in van der Walls interactions of inhibitors with MDMX relative to MDM2 is a main factor of weaker bindings of inhibitors to MDMX. The analyses of dihedral angles based on MD trajectories suggest that the closed conformation formed by the residues M53 and Y99 in MDMX leads to a potential steric clash with inhibitors and prevents inhibitors from arriving in the deep of MDMX binding cleft, which reduces the van der Waals contacts of inhibitors with M53, V92, P95 and L98. The calculated results using the residue-based free energy decomposition method further prove that the interaction strength of inhibitors with M53, V92, P95 and L98 from MDMX are obviously reduced compared to MDM2. We expect that this study can provide significant theoretical guidance for designs of potent dual inhibitors to block the p53-MDM2/MDMX interactions.

Free energy of binding of coiled-coil complexes with different electrostatic environments: the influence of force field polarisation and capping

Coiled-coils are well known protein-protein interaction motifs, with the leucine zipper region of activator protein-1 (AP-1) consisting of the c-Jun and c-Fos proteins being a typical example. Molecular dynamics (MD) simulations using the MM/GBSA method have been used to predict the free energy of interaction of these proteins. The influence of force field polarisation and capping on the predicted free energy of binding of complexes with different electrostatic environments (net charge) were investigated. Although both force field polarisation and peptide capping are important for the prediction of the absolute free energy of binding, peptide capping has the largest influence on the predicted free energy of binding. Polarisable simulations appear better suited to determine structural properties of the complexes of these proteins while non-polarisable simulations seem to give better predictions of the associated free energies of binding.

MMGBSA as a tool to understand the binding affinities of filamin-peptide interactions

Filamins (FLN) are large dimeric proteins that cross-link actin and work as important scaffolds in human cells. FLNs consist of an N-terminal actin-binding domain followed by 24 immunoglobulin-like domains (FLN1-24). FLN domains are divided into four subgroups based on their amino acid sequences. One of these subgroups, including domains 4, 9, 12, 17, 19, 21, and 23, shares a similar ligand-binding site between the β strands C and D. Several proteins, such as integrins β2 and β7, glycoprotein Ibα (GPIbα), and migfilin, have been shown to bind to this site. Here, we computationally estimated the binding free energies of filamin A (FLNa) subunits with bound peptides using the molecular mechanics-generalized Born surface area (MMGBSA) method. The obtained computational results correlated well with the experimental data, and they ranked efficiently both the binding of one ligand to all used FLNa-domains and the binding of all used ligands to FLNa21. Furthermore, the steered molecular dynamics (SMD) simulations pinpointed the binding hot spots for these complexes. These results demonstrate that molecular dynamics combined with free energy calculations are applicable to estimating the energetics of protein-protein interactions and can be used to direct the development of novel FLN function modulators.