Structural Insights into the Interaction between the C-Terminal-Deleted BH3-like Motif Peptide of Hepatitis B Virus X Protein and Bcl-xL.

Hepatitis B virus X protein (HBx) plays a crucial role in the development of hepatocellular carcinoma (HCC) associated with hepatitis B virus (HBV) infection. The full-length HBx protein interacts with Bcl-xL and is involved in the HBV replication and cell death processes. The three hydrophobic residues Trp120, Leu123, and Ile127 of the HBx BH3-like motif are essential for the Bcl-xL-binding. On the other hand, various lengths of C-terminal-truncated HBx mutants are frequently detected in HCC tissues, and these mutants, rather than the full-length HBx, appear to be responsible for HCC development. Notably, the region spanning residues 1-120 of HBx [HBx(1 and 120)] has been strongly associated with an increased risk of HCC development. However, the mode of interaction between HBx(1-120) and Bcl-xL remains unclear. HBx(1-120) possesses only Trp120 among the three hydrophobic residues essential for the Bcl-xL-binding. To elucidate this interaction mode, we employed a C-terminal-deleted HBx BH3-like motif peptide composed of residues 101-120. Here, we present the NMR complex structure of Bcl-xL and HBx(101-120). Our results demonstrate that HBx(101-120) binds to Bcl-xL in a weaker manner. Considering the high expression of Bcl-xL in HCC cells, this weak interaction, in conjunction with the overexpression of Bcl-xL in HCC cells, may potentially contribute to HCC development through the interaction between C-terminal-truncated HBx and Bcl-xL.

The carboxyl-terminal sequence of PUMA binds to both anti-apoptotic proteins and membranes.

Anti-apoptotic proteins such as BCL-XL promote cell survival by sequestering pro-apoptotic BCL-2 family members, an activity that frequently contributes to tumorigenesis. Thus, the development of small-molecule inhibitors for anti-apoptotic proteins, termed BH3-mimetics, is revolutionizing how we treat cancer. BH3 mimetics kill cells by displacing sequestered pro-apoptotic proteins to initiate tumor-cell death. Recent evidence has demonstrated that in live cells the BH3-only proteins PUMA and BIM resist displacement by BH3-mimetics, while others like tBID do not. Analysis of the molecular mechanism by which PUMA resists BH3-mimetic mediated displacement from full-length anti-apoptotic proteins (BCL-XL, BCL-2, BCL-W, and MCL-1) reveals that both the BH3-motif and a novel binding site within the carboxyl-terminal sequence (CTS) of PUMA contribute to binding. Together these sequences bind to anti-apoptotic proteins, which effectively 'double-bolt locks' the proteins to resist BH3-mimetic displacement. The pro-apoptotic protein BIM has also been shown to double-bolt lock to anti-apoptotic proteins however, the novel binding sequence in PUMA is unrelated to that in the CTS of BIM and functions independent of PUMA binding to membranes. Moreover, contrary to previous reports, we find that when exogenously expressed, the CTS of PUMA directs the protein primarily to the endoplasmic reticulum (ER) rather than mitochondria and that residues I175 and P180 within the CTS are required for both ER localization and BH3-mimetic resistance. Understanding how PUMA resists BH3-mimetic displacement will be useful in designing more efficacious small-molecule inhibitors of anti-apoptotic BCL-2 proteins.

Structural Analysis of the Interaction between Bcl-xL and the Noncanonical BH3 Domain of Non-Bcl-2 Family Proteins.

Anti-apoptotic and anti-autophagic Bcl-2 homologues commonly contain a hydrophobic groove in which the BH3 domain is accommodated. The BH3 domain is usually considered a feature of Bcl-2 family members; however, it has also been found in various non-Bcl-2 family proteins. Although interactions among Bcl-2 family members have been extensively investigated and highlighted, those mediated by the BH3 domain of non-Bcl-2 family proteins have not been the focus of substantial research. In this review, the author conducted a structural analysis of Bcl-xL complexed with the BH3 domain of four non-Bcl-2 family proteins, Beclin 1, SOUL, TCTP, and Pxt1, at an atomic level. Although the overall Bcl-xL-binding modes are similar among these proteins, they are characterized by limited sequence conservation of the BH3 consensus motif and differences in residues involved in complex formation. Based on the structural analysis, the author suggests that more "undiscovered" BH3 domain-containing proteins might exist, which have been unidentified due to their limited sequence conservation but can bind to Bcl-2 family proteins and control apoptosis, autophagy, or other biological processes.

Chemical shift assignments of a fusion protein comprising the C-terminal-deleted hepatitis B virus X protein BH3-like motif peptide and Bcl-xL.

Chronic hepatitis B virus (HBV) infection is a major risk factor for the development of liver diseases including fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). HBV has the multifunctional protein, HBV X protein (HBx, 154 residues), which plays key roles in HBV replication and liver disease development. Interaction of HBx through its BH3-like motif with the anti-apoptotic protein Bcl-xL leads to HBV replication and induction of apoptosis, resulting in HCC development. Our previous nuclear magnetic resonance (NMR) study revealed that the HBx BH3-like motif peptide (residues 101-136) binds to the common BH3-binding groove of Bcl-xL. Importantly, a C-terminal-truncated HBx, e.g., residues 1-120 of HBx, is strongly associated with the increased risk of HBV-related HCC development. However, the interaction mode between the C-terminal-truncated HBx and Bcl-xL remains unclear. To elucidate this interaction mode, the C-terminal-deleted HBx BH3-like motif peptide (residues 101-120) was used as a model peptide in this study. To facilitate the NMR analysis, we prepared a fusion protein of HBx (101-120) and Bcl-xL connected with five repeats of the glycine-serine dipeptide as a linker. Here, we report the 1H, 13C, and 15N resonance assignments of the fusion protein. This is the first step for the elucidation of the pathogenesis of liver diseases caused by the interaction between the C-terminal-truncated HBx and Bcl-xL.

Selectivity mechanism of BCL-XL/2 inhibition through in silico investigation.

The BCL-XL protein is among the most important members of the anti-apoptotic subfamily of the BCL-2 protein family, and is currently a promising new target for anti-tumor drug research. However, the BCL-XL/2 proteins have similar structures and functions, which could lead to undesirable side effects because of inhibitors that can bind to both BCL-XL and BCL-2. Therefore, it is crucial to expound on the structural basis of the selective mechanism towards BCL-XL/2 inhibition. In the current study, we employed hybrid computational methods including molecular docking and dynamics simulation, MM/GBSA energy calculation, alanine scanning mutagenesis and Hirshfeld surface analysis to comprehensively reveal the selectivity mechanism towards BCL-XL/2 from multiple perspectives, revealing the significant effects of the BCL-XL residues SER106 and LEU108 as well as the BCL-2 residue ASP103 on the inhibitory selectivity. Overall, our findings provide useful references for the rational design of BCL-XL/2 selective inhibitors with better affinity.

In silico SELEX screening and statistical analysis of newly designed 5mer peptide-aptamers as Bcl-xl inhibitors using the Taguchi method.

Drug development for cancer treatment is a complex process that requires special efforts. Targeting crucial proteins is the most essential purpose of drug design in cancers. Bcl-xl is an anti-apoptotic protein that binds to pro-apoptotic proteins and interrupts their signals. Pro-apoptotic Bcl-xl effectors are short BH3 sequences that form an alpha helix and bind to anti-apoptotic proteins to inhibit their activity. Computational systematic evolution of ligands by exponential enrichment (SELEX) is an exclusive approach for developing peptide aptamers as potential effectors. Here, the amino acids with a high tendency for constructing an alpha-helical structure were selected. Due to the enormous number of pentapeptides, Taguchi method was used to study a selected number of peptides. The binding affinity of the peptides to Bcl-xl was assessed using molecular docking, and after analysis of the obtained results, a final set of optimized peptides was arranged and constructed. For a better comparison, three chemical compounds with approved anti-Bcl-xl activity were selected for comparison with the top-ranked 5mer peptides. The optimized peptides showed considerable binding affinity to Bcl-xl. The molecular dynamics (MD) simulation indicated that the designed peptide (PO5) could create considerable interactions with the BH3 domain of Bcl-xl. The MM/GBSA calculations revealed that these interactions were even stronger than those created by chemical compounds. In silico SELEX is a novel approach to design and evaluate peptide-aptamers. The experimental design improves the SELEX process considerably. Finally, PO5 could be considered a potential inhibitor of Bcl-xl and a potential candidate for cancer treatment.

Exploration of α/ß/γ-peptidomimetics design for BH3 helical domains.

Systematic incorporation of ring-constrained ß- and γ-amino acid residues into α-helix mimetics engenders stable helical secondary structures. In this paper, functional α/ß/γ-helical peptidomimetics were explored for mimicry of BH3 helical domains, Bim as a pioneering study. The Bim-based α/ß/γ-peptides in an αγααßα-hexad repeat with five helical turns inhibited the interaction between Bak and Bcl-xL with excellent resistance towards proteolytic digestion. Further optimization of the α/ß/γ-backbone strategy will considerably expand the utility of functional α/ß/γ-peptidomimetics, in particular due to its prominent stability against proteolysis.

Artificial intelligence-based identification of octenidine as a Bcl-xL inhibitor.

Apoptosis plays an essential role in maintaining cellular homeostasis and preventing cancer progression. Bcl-xL, an anti-apoptotic protein, is an important modulator of the mitochondrial apoptosis pathway and is a promising target for anticancer therapy. In this study, we identified octenidine as a novel Bcl-xL inhibitor through structural feature-based deep learning and molecular docking from a library of approved drugs. The NMR experiments demonstrated that octenidine binds to the Bcl-2 homology 3 (BH3) domain-binding hydrophobic region that consists of the BH1, BH2, and BH3 domains in Bcl-xL. A structural model of the Bcl-xL/octenidine complex revealed that octenidine binds to Bcl-xL in a similar manner to that of the well-known Bcl-2 family protein antagonist ABT-737. Using the NanoBiT protein-protein interaction system, we confirmed that the interaction between Bcl-xL and Bak-BH3 domains within cells was inhibited by octenidine. Furthermore, octenidine inhibited the proliferation of MCF-7 breast and H1299 lung cancer cells by promoting apoptosis. Taken together, our results shed light on a novel mechanism in which octenidine directly targets anti-apoptotic Bcl-xL to trigger mitochondrial apoptosis in cancer cells.

Impact of Deamidation on the Structure and Function of Antiapoptotic Bcl-xL.

Bcl-xL is an antiapoptotic mitochondrial trans-membrane protein, which is known to play a crucial role in the survival of tumor cells. The deamidation of Bcl-xL is a pivotal switch that regulates its biological function. The potential impact of deamidation on the structure and dynamics of Bcl-xL is directly linked to the intrinsically disordered region (IDR), which is the main site for post-translational modifications (PTMs). In this study, we explored deamidation-induced conformational changes in Bcl-xL to gain insight into its loss of function by performing microsecond-long molecular dynamics (MD) simulations. MD simulation outcomes showed that the IDR motion and interaction patterns have changed notably upon deamidation. Principal component analysis (PCA) demonstrates significant differences between wild-type and deamidated Bcl-xL and suggests that deamidation affects the structure and dynamics of Bcl-xL. The combination of clustering analysis, H-bond analysis, and PCA revealed changes in conformation, interaction, and dynamics upon deamidation. Differences in contact patterns and essential dynamics that lead to a narrowing in the binding groove (BG) are clear indications of deamidation-induced allosteric effects. In line with previous studies, we show that the IDR plays a very important role in the loss of apoptotic functions of Bcl-xL while providing a unique perspective on the underlying mechanism of Bcl-xL deamidation-induced cell death.

Development of a BCL-xL and BCL-2 dual degrader with improved anti-leukemic activity.

PROteolysis-TArgeting Chimeras (PROTACs) have emerged as an innovative drug development platform. However, most PROTACs have been generated empirically because many determinants of PROTAC specificity and activity remain elusive. Through computational modelling of the entire NEDD8-VHL Cullin RING E3 ubiquitin ligase (CRLVHL)/PROTAC/BCL-xL/UbcH5B(E2)-Ub/RBX1 complex, we find that this complex can only ubiquitinate the lysines in a defined band region on BCL-xL. Using this approach to guide our development of a series of ABT263-derived and VHL-recruiting PROTACs, we generate a potent BCL-xL and BCL-2 (BCL-xL/2) dual degrader with significantly improved antitumor activity against BCL-xL/2-dependent leukemia cells. Our results provide experimental evidence that the accessibility of lysines on a target protein plays an important role in determining the selectivity and potency of a PROTAC in inducing protein degradation, which may serve as a conceptual framework to guide the future development of PROTACs.

Computational Evaluation of Abrogation of HBx-Bcl-xL Complex with High-Affinity Carbon Nanotubes (Fullerene) to Halt the Hepatitis B Virus Replication.

Hepatitis B virus (HBV) is the world's most prevalent chronic viral infection. More than 350 million individuals are chronic carriers of the virus, with an estimated 2 billion infected persons. For instance, the role of HBx protein in attachment and infection is very obvious and consequently deemed as an important druggable target. Targeting the interface and discovering novel drugs greatly advanced the field of therapeutics development. Therefore, in the current study, HBx to Bcl-xL is abrogated on high-affinity carbon nanotubes using computational structural biology tools. Our analysis revealed that among the total 62 carbon fullerenes, only 13 compounds exhibited inhibitory activity against HBx, which was further confirmed through IFD-based rescoring. Structural dynamics investigation revealed stable binding, compactness, and hydrogen bonds reprogramming. Moreover, the binding free energy calculation results revealed that the top hits1-4 possess the total binding energy of -54.36 kcal/mol (hit1), -50.81 kcal/mol (hit2), -47.09 kcal/mol (hit3), and -45.59 kcal/mol for hit4. In addition, the predicted KD values and bioactivity scores further validated the inhibitory potential of these top hits. The identified compounds need further in vitro and in vivo validation to aid the treatment process of HBV.

Difference in the binding mechanisms of ABT-263/43b with Bcl-xL/Bcl-2: computational perspective on the accurate binding free energy analysis.

B-cell lymphoma/leukemia gene-2(Bcl-2) protein family known for regulating cell cycle arrest and subsequent cell death is highly expressed in a variety of cancers. Among them, the Bcl-xL and Bcl-2 are two essential proteins in the Bcl-2 family. In the present work, the differences in binding modes as between the two proteins and two ligands ABT-263/43b were investigated and compared. And the computational alanine scanning combined with the recently developed interaction entropy (AS-IE) method was employed for predicting their binding free energies and finding those amino acids that were more critical during the binding process. The result showed that the binding free energy calculated by the AS-IE method was more in line with experimental values than the molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method. Besides, no significant difference was found between Bcl-xL and ABT-263/43b in the binding free energy, which Bcl-xL showed slightly weaker binding free energy to 43b because of the fewer number of key residues with interactions. Nonetheless, compared with the Bcl-2 and 43b complex, the Bcl-2 and ABT-263 system had greater number of key residues interacting with ABT-263, in particular, contribute favorably, resulting in a stronger binding ability for the Bcl-2 and ABT-263 systems. The van der Waals and hydrogen bond contributions were significant in the four protein-ligand complexes. Overall, Tyr108 was found to be the common key residues in the Bcl-xL-ligand complex, while Tyr105, Glu100, and Glu143 were established as the common key residue in the Bcl-2-ligand systems. We hope that the predicted hot spot residues and their energy distributions can guide the design of peptide and small-molecule drugs targeting Bcl-xL and Bcl-2.

Novel Synthesized N-Ethyl-Piperazinyl-Amides of C2-Substituted Oleanonic and Ursonic Acids Exhibit Cytotoxic Effects through Apoptotic Cell Death Regulation.

A series of novel hybrid chalcone N-ethyl-piperazinyl amide derivatives of oleanonic and ursonic acids were synthesized, and their cytotoxic potential was evaluated in vitro against the NCI-60 cancer cell line panel. Compounds 4 and 6 exhibited the highest overall anticancer activity, with GI50 values in some cases reaching nanomolar values. Thus, the two compounds were further assessed in detail in order to identify a possible apoptosis- and antiangiogenic-based mechanism of action induced by the assessed compounds. DAPI staining revealed that both compounds induced nuclei condensation and overall cell morphological changes consistent with apoptotic cell death. rtPCR analysis showed that up-regulation of pro-apoptotic Bak gene combined with the down-regulation of the pro-survival Bcl-XL and Bcl-2 genes caused altered ratios between the pro-apoptotic and anti-apoptotic proteins' levels, leading to overall induced apoptosis. Molecular docking analysis revealed that both compounds exhibited high scores for Bcl-XL inhibition, suggesting that compounds may induce apoptotic cell death through targeted anti-apoptotic protein inhibition, as well. Ex vivo determinations showed that both compounds did not significantly alter the angiogenesis process on the tested cell lines.

High-resolution analysis of the conformational transition of pro-apoptotic Bak at the lipid membrane.

Permeabilization of the outer mitochondrial membrane by pore-forming Bcl2 proteins is a crucial step for the induction of apoptosis. Despite a large set of data suggesting global conformational changes within pro-apoptotic Bak during pore formation, high-resolution structural details in a membrane environment remain sparse. Here, we used NMR and HDX-MS (Hydrogen deuterium exchange mass spectrometry) in lipid nanodiscs to gain important high-resolution structural insights into the conformational changes of Bak at the membrane that are dependent on a direct activation by BH3-only proteins. Furthermore, we determined the first high-resolution structure of the Bak transmembrane helix. Upon activation, α-helix 1 in the soluble domain of Bak dissociates from the protein and adopts an unfolded and dynamic potentially membrane-bound state. In line with this finding, comparative protein folding experiments with Bak and anti-apoptotic BclxL suggest that α-helix 1 in Bak is a metastable structural element contributing to its pro-apoptotic features. Consequently, mutagenesis experiments aimed at stabilizing α-helix 1 yielded Bak variants with delayed pore-forming activity. These insights will contribute to a better mechanistic understanding of Bak-mediated membrane permeabilization.

Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL.

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Ligand-Receptor Interactions and Drug Design.

In silico rational drug design is one of the major pylons in the drug discovery process. Drugs usually act on specific targets such as proteins, DNA, and lipid bilayers. Thus, molecular docking is an essential part of the rational drug design process. Molecular docking uses specific algorithms and scoring functions to reveal the strength of the interaction of the ligand to its target. AutoDock is a molecular docking suite that offers a variety of algorithms to tackle specific problems. These algorithms include Monte Carlo Simulated Annealing (SA), a Genetic Algorithm (GA), and a hybrid local search GA, also known as the Lamarckian Genetic Algorithm (LGA). This chapter aims to acquaint the reader with the docking process using AutoDockTools (GUI of AutoDock). Furthermore, herein is described the docking process of calf thymus DNA with three metal complexes, as a potential metallo-therapeutics as also the docking process of the plant flavonoid quercetin to the antiapoptotic protein BcL-xL.

Computational analysis of binding free energies, hotspots and the binding mechanism of Bcl-xL/Bcl-2 binding to Bad/Bax.

The anti-apoptotic proteins B-cell lymphoma-extra large (Bcl-xL) and B-cell lymphoma/leukemia-2 (Bcl-2) are members of the Bcl-2 protein family, and they play important roles in regulating apoptosis and cell cycle retardation. However, the binding mechanisms of Bcl-xL/Bcl-2 with their associated agonists, including Bcl-2-associated death promoter (Bad) and Bcl-2-associated X protein (Bax), are not well understood. In the present study, the recently developed interaction entropy approach was employed for the calculation of entropic contribution, and the computational alanine scanning method was used to identify the hot spot in the protein-protein interactions between Bcl-xL/Bcl-2 and Bad/Bax. The calculated binding free energies and their ranks for the four systems were in good agreement with the experimental results. Computational analysis shows that there are more hot-spot residues in the Bcl-xL/Bad complex than that in the Bcl-xL/Bax complex, leading to a stronger binding affinity in the former. It is interesting to find that the reason for the stronger binding affinity of Bcl-2 to Bad than to Bax is different for the Bcl-xL system. Although there are more hot-spot residues in the Bcl-2/Bax system than in the Bcl-2/Bad complex, there are also more negatively contributing residues in the Bcl-2/Bax. Our study identified Arg104, Tyr105, Leu116, and Leu134 to be the common key residues in the Bcl-xL complexes, and Arg107, Tyr108, Phe112, Gln118, Leu137, Arg146, and Tyr202 are common key residues in the Bcl-2 complexes. These results would provide valuable information for the design of potent inhibitors of Bcl-xL/Bcl-2.

Biophysical Characterization of Pro-apoptotic BimBH3 Peptides Reveals an Unexpected Capacity for Self-Association.

Bcl-2 proteins orchestrate the mitochondrial pathway of apoptosis, pivotal for cell death. Yet, the structural details of the conformational changes of pro- and antiapoptotic proteins and their interactions remain unclear. Pulse dipolar spectroscopy (double electron-electron resonance [DEER], also known as PELDOR) in combination with spin-labeled apoptotic Bcl-2 proteins unveils conformational changes and interactions of each protein player via detection of intra- and inter-protein distances. Here, we present the synthesis and characterization of pro-apoptotic BimBH3 peptides of different lengths carrying cysteines for labeling with nitroxide or gadolinium spin probes. We show by DEER that the length of the peptides modulates their homo-interactions in the absence of other Bcl-2 proteins and solve by X-ray crystallography the structure of a BimBH3 tetramer, revealing the molecular details of the inter-peptide interactions. Finally, we prove that using orthogonal labels and three-channel DEER we can disentangle the Bim-Bim, Bcl-xL-Bcl-xL, and Bim-Bcl-xL interactions in a simplified interactome.

Chemical Translational Biology: Redefining Druggability of Protein-Protein Interactions.

Chemical biologists use chemical tools to answer biological questions. The translational application of these principles has led to an explosion in the discovery and druggability of new protein targets, including protein-protein interactions (PPIs). Proteins tend to interact with other macromolecules using relatively large and featureless binding surfaces, which has hampered traditional drug discovery efforts, particularly for interactions with weaker affinity. In this article, I discuss several emerging strategies for targeting PPIs, including computational and structural methods and novel screening approaches. In particular, I focus on hijacking intrinsic protein allosteric pathways for the discovery and design of small-molecule and peptide ligands.

Computational Design of BH3-Mimetic Peptide Inhibitors That Can Bind Specifically to Mcl-1 or Bcl-XL: Role of Non-Hot Spot Residues.

Interactions between pro- and anti-apoptotic Bcl-2 proteins decide the fate of the cell. The BH3 domain of pro-apoptotic Bcl-2 proteins interacts with the exposed hydrophobic groove of their anti-apoptotic counterparts. Through their design and development, BH3 mimetics that target the hydrophobic groove of specific anti-apoptotic Bcl-2 proteins have the potential to become anticancer drugs. We have developed a novel computational method for designing sequences with BH3 domain features that can bind specifically to anti-apoptotic Mcl-1 or Bcl-XL. In this method, we retained the four highly conserved hydrophobic and aspartic residues of wild-type BH3 sequences and randomly substituted all other positions to generate a large number of BH3-like sequences. We modeled 20000 complex structures with Mcl-1 or Bcl-XL using the BH3-like sequences derived from five wild-type pro-apoptotic BH3 peptides. Peptide-protein interaction energies calculated from these models for each set of BH3-like sequences resulted in negatively skewed extreme value distributions. The selected BH3-like sequences from the extreme negative tail regions have highly favorable interaction energies with Mcl-1 or Bcl-XL. They are enriched in acidic and basic residues when they bind to Mcl-1 and Bcl-XL, respectively. With the charged residues often away from the binding interface, the overall electric field generated by the charged residues results in strong long-range electrostatic interaction energies between the peptide and the protein giving rise to high specificity. Cell viability studies of representative BH3-like peptides further validated the predicted specificity. This study has revealed the importance of non-hot spot residues in BH3-mimetic peptides in providing specificity to a particular anti-apoptotic protein.