Molecular dynamics study of segment peptides of Bax,Bim, and Mcl-1 BH3 domain of the apoptosis-regulating proteins bound to the anti-apoptotic Mcl-1 protein

Mcl-1 has emerged as a potential therapeutic target in the treatment of several malignancies. Peptides representing BH3 region of pro-apoptotic proteins have been shown to bind the hydrophobic cleft of anti-apoptotic Mcl-1 and this segment is responsible for modulating the apoptotic pathways in living cells. Understanding the molecular basis of protein–peptide interaction is required to develop potent inhibitors specific for Mcl-1. Molecular dynamics simulations were performed for Mcl-1 in complex with three different BH3 peptides derived from Mcl-1, Bax, and Bim. Accordingly, the calculated binding free energies using MM-PBSA method are obtained and comparison with the experimentally determined binding free energies is made. The interactions involving two conserved charged residues (Aspi, and Arg/Lysi-4) and three upstream conserved hydrophobic residues (Leui-5, Ile/Vali-2, and Glyi-1, respectively) of BH3 peptides play the important roles in the structural stability of the complexes. The calculated results exhibit that the interactions of Bim BH3 peptides to Mcl-1 is stronger than the complex with Bax 19BH3 peptides. The hydrophobic residues (position i???9, i???8 and i?+?2) of BH3 peptides can be involved in their inhibitory specificity. The calculated results can be used for designing more effective MCL-1 inhibitors.     

Interaction of a putative BH3 domain of clusterin with anti-apoptotic Bcl-2 family proteins as revealed by NMR spectroscopy

Clusterin (CLU) is a multifunctional glycoprotein that is overexpressed in prostate and breast cancers. Although CLU is known to be involved in the regulation of apoptosis and cell survival, the precise molecular mechanism underlying the pro-apoptotic function of nuclear CLU (nCLU) remains unclear. In this study, we identified a conserved BH3 motif in C-terminal coiled coil (CC2) region of nCLU by sequence analysis and characterized the molecular interaction of the putative nCLU BH3 domain with anti-apoptotic Bcl-2 family proteins by nuclear magnetic resonance (NMR) spectroscopy. The chemical shift perturbation data demonstrated that the nCLU BH3 domain binds to pro-apoptotic BH3 peptide-binding grooves in both Bcl-X_L and Bcl-2. A structural model of the Bcl-X_L/nCLU BH3 peptide complex reveals that the binding mode is remarkably similar to those of other Bcl-X_L/BH3 peptide complexes. In addition, mutational analysis confirmed that Leu323 and Asp328 of nCLU BH3 domain, absolutely conserved in the BH3 motifs of BH3-only protein family, are critical for binding to Bcl-X_L. Taken altogether, our results suggest a molecular basis for the pro-apoptotic function of nCLU by elucidating the residue specific interactions of the BH3 motif in nCLU with anti-apoptotic Bcl-2 family proteins.     

Amphipathic tail-anchoring peptide and Bcl-2 homology domain-3 (BH3) peptides from Bcl-2 family proteins induce apoptosis through different mechanisms

Bcl-2 homology domain-3 (BH3) peptides are potent cancer therapeutic reagents that target regulators of apoptotic cell death in cancer cells. However, their cytotoxic effects are affected by different expression levels of Bcl-2 family proteins. We recently found that the amphipathic tail-anchoring peptide (ATAP) from Bfl-1, a bifunctional Bcl-2 family member, produced strong pro-apoptotic activity by permeabilizing the mitochondrial outer membrane. Here, we test whether the activity of ATAP requires other cellular factors and whether ATAP has an advantage over the BH3 peptides in targeting cancer cells. Confocal microscopic imaging illustrates specific targeting of ATAP to mitochondria, whereas BH3 peptides show diffuse patterns of cytosolic distribution. Although the pro-apoptotic activities of BH3 peptides are largely inhibited by either overexpression of anti-apoptotic Bcl-2 or Bcl-xL or nullification of pro-apoptotic Bax and Bak in cells, the pro-apoptotic function of ATAP is not affected by these cellular factors. Reconstitution of synthetic ATAP into liposomal membranes results in release of fluorescent molecules of the size of cytochrome c from the liposomes, suggesting that the membrane permeabilizing activity of ATAP does not require additional protein factors. Because ATAP can target to the mitochondrial membrane and its pro-apoptotic activity does not depend on the content of Bcl-2 family proteins, it represents a promising candidate for anti-cancer drugs that can potentially overcome the intrinsic apoptosis-resistant nature of cancer cells.     

A synthetic peptide corresponding to a region of the human pericentriolar material 1 (PCM‐1) protein binds β‐amyloid (Aβ1‐42) oligomers

We have recently reported that a ~19‐kDa polypeptide, rPK‐4, is a protein kinase Cs inhibitor that is 89% homologous to the 1171–1323 amino acid region of the 228‐kDa human pericentriolar material‐1 (PCM‐1) protein (Chakravarthy et al. 2012). We have now discovered that rPK‐4 binds oligomeric amyloid‐β peptide (Aβ)_1‐42 with high affinity. Most importantly, a PCM‐1‐selective antibody co‐precipitated Aβ and amyloid β precursor protein (AβPP) from cerebral cortices and hippocampi from AD (Alzheimer's disease) transgenic mice that produce human AβPP and Aβ_1‐42, suggesting that PCM‐1 may interact with amyloid precursor protein/Aβ in vivo. We have identified rPK‐4′s Aβ‐binding domain using a set of overlapping synthetic peptides. We have found with ELISA, dot‐blot, and polyacrylamide gel electrophoresis techniques that a ~ 5 kDa synthetic peptide, amyloid binding peptide (ABP)‐p4‐5 binds Aβ_1‐42 at nM levels. Most importantly, ABP‐p4‐5, like rPK‐4, appears to preferentially bind Aβ_1‐42 oligomers, believed to be the toxic AD‐drivers. As expected from these observations, ABP‐p4‐5 prevented Aβ_1‐42 from killing human SH‐SY5Y neuroblastoma cells via apoptosis. These findings indicate that ABP‐p4‐5 is a possible candidate therapeutic for AD.     

Promiscuous and Selective: How Intrinsically Disordered BH3 Proteins Interact with Their Pro-survival Partner MCL-1

The BCL-2 family of proteins plays a central role in regulating cell survival and apoptosis. Disordered BH3-only proteins bind promiscuously to a number of different BCL-2 proteins, with binding affinities that vary by orders of magnitude. Here we investigate the basis for these differences in affinity. We show that eight different disordered BH3 proteins all bind to their BCL-2 partner (MCL-1) very rapidly, and that the differences in sequences result in different dissociation rates. Similarly, mutation of the binding surface of MCL-1 generally affects association kinetics in the same way for all BH3 peptides but has significantly different effects on the dissociation rates. Importantly, we infer that the evolution of homologous, competing interacting partners has resulted in complexes with significantly different lifetimes.     

3D-QSAR and molecular recognition of Klebsiella pneumoniae NDM-1 inhibitors

New Delhi metallo-β-lactamase-1 (NDM-1) as a target for the development of anti-superbug agents, plays an important role in the resistance of β-lactam antibiotics and has received worldwide attention. Sulfhydryl propionic acid derivatives can effectively inhibit the catalytic activity of NDM-1, but the quantitative structure–activity relationship (QSAR) and inhibitor-target recognition mechanism both remain unclear. In this work, CoMFA and CoMSIA models of sulfhydryl propionic acid inhibitors with high predictive ability were obtained, from which the effect of different substituents on the inhibitory activity against NDM-1 were revealed at the molecular level. Then, two 120-nanosecond comparative molecular dynamics (MD) simulations for NDM-1 enzyme and NDM-1-inhibitor complex systems were performed to study the recognition and inhibition mechanism of sulfhydryl propionic acid derivatives. The binding of inhibitors alters the entire H-bond network of the NDM-1 system accompanied by the formation of strong interactions with I35, W93, H120, H122, D124, H189 and H250, further weakens the recognition of NDM-1 by its endogenic substrates. Finally, the results of free energy landscape and conformation cluster analyses show that NDM-1 underwent a significant conformational change after the association with sulfhydryl propionic acid inhibitors. Our findings can provide theoretical support and help for anti-superbug agents design based on the structures of NDM-1 and sulfhydryl propionic acid derivatives.     

α4 phosphoprotein interacts with EDD E3 ubiquitin ligase and poly(A)‐binding protein

Mammalian α4 phosphoprotein, the homolog of yeast Tap42, is a component of the mammalian target‐of‐rapamycin (mTOR) pathway that regulates ribogenesis, the initiation of translation, and cell‐cycle progression. α4 is known to interact with the catalytic subunit of protein phosphatase 2A (PP2Ac) and to regulate PP2A activity. Using α4 as bait in yeast two‐hybrid screening of a human K562 erythroleukemia cDNA library, EDD (E3 isolated by differential display) E3 ubiquitin ligase was identified as a new protein partner of α4. EDD is the mammalian ortholog of Drosophila hyperplastic discs gene (hyd) that controls cell proliferation during development. The EDD protein contains a PABC domain that is present in poly(A)‐binding protein (PABP), suggesting that PABP may also interact with α4. PABP recruits translation factors to the poly(A)‐tails of mRNAs. In the present study, immunoprecipitation/immunoblotting (IP/IB) analyses showed a physical interaction between α4 and EDD in rat Nb2 T‐lymphoma and human MCF‐7 breast cancer cell lines. α4 also interacted with PABP in Nb2, MCF‐7 and the human Jurkat T‐leukemic and K562 myeloma cell lines. COS‐1 cells, transfected with Flag‐tagged‐pSG5‐EDD, gave a (Flag)‐EDD–α4 immunocomplex. Furthermore, deletion mutants of α4 were constructed to determine the binding site for EDD. IP/IB analysis showed that EDD bound to the C‐terminal region of α4, independent of the α4‐PP2Ac binding site. Therefore, in addition to PP2Ac, α4 interacts with EDD and PABP, suggesting its involvement in multiple steps in the mTOR pathway that leads to translation initiation and cell‐cycle progression. J. Cell. Biochem. 110: 1123–1129, 2010. Published 2010 Wiley‐Liss, Inc.     

A computational study of binding between 3-(4-fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide and tubulin

3-(4-Fluorophenyl)-N-((4-fluorophenyl)sulphonyl)acrylamide (FFSA) is a potential tubulin polymerisation inhibitor. In this article, a theoretical study of the binding between FFSA and tubulin in colchicine site was carried out by molecular docking, molecular dynamics (MD) simulation and binding free energy calculations. The docking calculations preliminarily indicate that there are three possible binding modes 1, 2 and 3; MD simulations and binding free energy calculations identify that binding mode 2 is the most favourable, with the lowest binding free energy of ? 29.54 kcal/mol. Moreover, our valuable results for the binding are as follows: the inhibitor FFSA is suitably located at the colchicine site of tubulin, where it not only interacts with residues Leu248β, Lys254β, Leu255β, Lys352β, Met259β and Val181a by hydrophilic interaction, but also interacts with Val181α and Thr179α by hydrogen bond interaction. These two factors are both essential for FFSA strongly binding to tubulin. These theoretical results help understanding the action mechanism and designing new compounds with higher affinity to tubulin.     

Structural conservation of residues in BH1 and BH2 domains of Bcl-2 family proteins

The sequence of Bcl-2 homology domains, BH1 and BH2, is known to be conserved among anti- and pro-apoptotic members of Bcl-2 family proteins. But structural conservation of these domains with respect to functionally active residues playing role in heterodimerization-mediated regulation of apoptosis has never been elucidated. Here, we have suggested the formation of an active site by structurally conserved residues in BH1 (glycine, arginine) and BH2 (tryptophan) domains of Bcl-2 family members, which also accounts for the functional effect of known mutations in BH1 (G145A, G145E) and BH2 (W188A) domains of Bcl-2.     

Computational design of small molecular modulators of protein–protein interactions with a novel thermodynamic cycle: Allosteric inhibitors of HIV‐1 integrase

Targeting protein–protein interactions for therapeutic development involves designing small molecules to either disrupt or enhance a known PPI. For this purpose, it is necessary to compute reliably the effect of chemical modifications of small molecules on the protein–protein association free energy. Here we present results obtained using a novel thermodynamic free energy cycle, for the rational design of allosteric inhibitors of HIV‐1 integrase (ALLINI) that act specifically in the early stage of the infection cycle. The new compounds can serve as molecular probes to dissect the multifunctional mechanisms of ALLINIs to inform the discovery of new allosteric inhibitors. The free energy protocol developed here can be more broadly applied to study quantitatively the effects of small molecules on modulating the strengths of protein–protein interactions.     

BNIP3 heterodimerizes with Bcl-2/Bcl-X(L) and induces cell death independent of a Bcl-2 homology 3 (BH3) domain at both mitochondrial and nonmitochondrial sites

BNIP3 (formerly NIP3) is a pro-apoptotic, mitochondrial protein classified in the Bcl-2 family based on limited sequence homology to the Bcl-2 homology 3 (BH3) domain and COOH-terminal transmembrane (TM) domain. BNIP3 expressed in yeast and mammalian cells interacts with survival promoting proteins Bcl-2, Bcl-X(L), and CED-9. Typically, the BH3 domain of pro-apoptotic Bcl-2 homologues mediates Bcl-2/Bcl-X(L) heterodimerization and confers pro-apoptotic activity. Deletion mapping of BNIP3 excluded its BH3-like domain and identified the NH(2) terminus (residues 1-49) and TM domain as critical for Bcl-2 heterodimerization, and either region was sufficient for Bcl-X(L) interaction. Additionally, the removal of the BH3-like domain in BNIP3 did not diminish its killing activity. The TM domain of BNIP3 is critical for homodimerization, pro-apoptotic function, and mitochondrial targeting. Several TM domain mutants were found to disrupt SDS-resistant BNIP3 homodimerization but did not interfere with its killing activity or mitochondrial localization. Substitution of the BNIP3 TM domain with that of cytochrome b(5) directed protein expression to nonmitochondrial sites and still promoted apoptosis and heterodimerization with Bcl-2 and Bcl-X(L). We propose that BNIP3 represents a subfamily of Bcl-2-related proteins that functions without a typical BH3 domain to regulate apoptosis from both mitochondrial and nonmitochondrial sites by selective Bcl-2/Bcl-X(L) interactions.     

Direct Interaction of Bax and Bak Proteins with Bcl-2 Homology Domain 3 (BH3)-only Proteins in Living Cells Revealed by Fluorescence Complementation

The key event in the mitochondrial pathway of apoptosis is the activation of Bax and Bak by BH3-only proteins through a molecular mechanism that is still a matter of debate. Here we studied interactions among anti- and proapoptotic proteins of the Bcl-2 family in living cells by using bimolecular fluorescence complementation analysis. Our results indicate that the antiapoptotic proteins Mcl-1 and Bcl-x_L bind preferably to the BH3-only proteins Bim, PUMA, and Noxa but can also bind to Bak and Bax. We also found a direct interaction between Bim, PUMA, or Noxa with either Bax or Bak during apoptosis induction. In HeLa cells, interaction of Bim with Bax occurs in cytosol, and then Bim-Bax complexes translocate to mitochondria. Complexes of either PUMA or Noxa with Bax or Bak were always detected at mitochondria. Overexpression of Bcl-x_L or Mcl-1 delayed Bim/Bax translocation to mitochondria. These results reveal the ability of main BH3-only proteins to directly activate Bax and Bak in living cells and suggest that a complex network of interactions regulate the function of Bcl-2 family members during apoptosis.     

The beta-isoform of the BRCA2 and CDKN1A(p21)-interacting protein (BCCIP) stabilizes nuclear RPL23/uL14

BRCA2 and CDKN1A(p21,CIP1)-interacting protein (BCCIP) is an evolutionary conserved protein implicated in maintenance of genome stability and cell cycle progression. Two isoforms of BCCIP with distinct C-terminal domains exist in humans. We show that mammalian BCCIPβ, but not BCCIPα, forms a ternary complex with the ribosomal protein RPL23/uL14 and the pre-60S trans-acting factor eIF6. Complex formation is dependent on an intact C-terminal domain of BCCIPβ. Depletion of BCCIPβ reduces the pool of free RPL23, and decreases eIF6 levels in nucleoli. Overexpression of BCCIPβ leads to nucleoplasmic accumulation of extra-ribosomal RPL23 and stabilizes overexpressed RPL23, suggesting that BCCIPβ functions as nuclear chaperone for RPL23.     

Yeast two-hybrid screens imply that GGNBP1, GGNBP2 and OAZ3 are potential interaction partners of testicular germ cell-specific protein GGN1

Gametogenetin (Ggn) is a testicular germ cell-specific gene specifically expressed from late pachytene spermatocytes through round spermatids. The function of gametogenetin protein 1 (GGN1) remains unknown. Here, we used the yeast two-hybrid approach to look for more GGN1 interacting proteins. We found that gametogenetin binding protein 1 (GGNBP1), gametogenetin binding protein 2 (GGNBP2) and ornithine decarboxylase antizyme 3 (OAZ3) were potential GGN1 interaction partners. We determined the regions mediating the interactions and further showed the interactions between the proteins in mammalian cells by colocalization and coimmunoprecipitation experiments. Our work suggested that GGN1, GGNBP1, GGNBP2 and OAZ3 could be involved in a common process associated with spermatogenesis.     

Conformational behavior of temporin A and temporin L in aqueous solution: a computational/experimental study

Molecular dynamics (MD) simulations and circular dichroism (CD) experiments were carried out on aqueous temporin A and L, two short peptides belonging to an interesting class of natural substances known to be active mainly against Gram-positive/negative bacteria and fungi. Experimental results indicate the higher propensity of temporin L, with respect to temporin A, in forming alpha-helical structures. These results were revisited by long-timescale MD simulations, in which their alpha-helical propensity was investigated in the absence of trifluoroethanol. Results clearly show the higher stability of alpha-helix conformations in temporin L; moreover, an interestingly strong mechanical analogy emerges since both temporins show the same residue interval (from 7 to 10) as the most energetically accessible for alpha-helix formation. Such studies provide some intriguing structural and mechanical evidence that may help in better understanding and rationalizing the conformational behaviour of temporins in water solution and, ultimately, the inner principles of their microbial targets selectivity and mechanism of action at the level of cell membranes.     

Molecular basis of the Tfs1/Ira2 interaction: a combined protein engineering and molecular modelling study

Tfs1p and Ylr179cp are yeast proteins belonging to the PEBP family. Tfs1p, but not Ylr179cp, has been shown to interact with and inhibit Ira2p, a GTPase-activating protein of Ras. Tfs1p has been shown to be a specific inhibitor of the CPY protease and the 3D structure of the complex has been resolved. To shed light on the molecular determinants of Tfs1p involved in the Tfs1/Ira2 interaction, the 3D structure of Ylr179cp has been modelled and compared to that of Tfs1p. Tfs1p point mutants and Tfs1 hybrid proteins combining regions of Tfs1p and Ylr179cp were also designed and their function was tested. Results, interpreted from a structural point of view, show that the accessibility of the surface pocket of Tfs1p, its N-terminal region and the specific electrostatic properties of a large surface region containing these two elements, play a crucial role in this interaction.