Next Article in Journal
Phylogeography of Rotavirus G8P[8] Detected in Argentina: Evidence of Transpacific Dissemination
Previous Article in Journal
Isolation, Identification, and Genetic Phylogenetic Analysis of Two Different Genotypes of Bovine Parainfluenza 3 Virus in China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 10
10.3390/v14102222
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Crystal Structures of Epstein–Barr Virus Bcl-2 Homolog BHRF1 Bound to Bid and Puma BH3 Motif Peptides

by
Chathura D. Suraweera
1,†,
Mark G. Hinds
2,* and
Marc Kvansakul
1,*
1
Department of Biochemistry & Chemistry, La Trobe Institute for Molecular Science, La Trobe University, Melbourne, VIC 3086, Australia
2
Bio21 Molecular Science and Biotechnology Institute, The University of Melbourne, Parkville, VIC 3010, Australia
*
Authors to whom correspondence should be addressed.
†
Current address: Monash Biomedical Discovery Institute, Department of Microbiology, Monash University, Melbourne, VIC 3800, Australia.
Viruses 2022, 14(10), 2222; https://doi.org/10.3390/v14102222
Submission received: 29 August 2022 / Revised: 27 September 2022 / Accepted: 4 October 2022 / Published: 9 October 2022
(This article belongs to the Section Human Virology and Viral Diseases)
Browse Figures
Versions Notes

Abstract

:
Apoptosis is a powerful defense mechanism used by multicellular organisms to counteract viral infection. In response to premature host cell suicide, viruses have evolved numerous countermeasures to ensure cell viability to optimize their replication by encoding proteins homologous in structure and function to cellular pro-survival Bcl-2 proteins. Epstein–Barr virus (EBV), a member of the Gammaherpesviridae, encodes the Bcl-2 homolog BHRF1, a potent inhibitor of Bcl-2-mediated apoptosis. BHRF1 acts by directly targeting Bid and Puma, two proapoptotic proteins of the Bcl-2 family. Here, we determined the crystal structures of BHRF1 bound to peptides spanning the Bcl-2 binding motifs (Bcl-2 homology 3 motif, BH3) of Bid and Puma. BHRF1 engages BH3 peptides using the canonical ligand-binding groove of its Bcl-2 fold and maintains a salt bridge between an Arg residue with a conserved Asp residue in the BH3 motif mimicking the canonical ionic interaction seen in host Bcl-2:BH3 motif complexes. Furthermore, both Bid and Puma utilize a fifth binding pocket in the canonical ligand binding groove of BHRF1 to provide an additional hydrophobic interaction distinct from the interactions previously seen with Bak and Bim. These findings provide a structural basis for EBV-mediated suppression of host cell apoptosis and reveal the flexibility of virus encoded Bcl-2 proteins in mimicking key interactions from the endogenous host signaling pathways.

1. Introduction

Epstein–Barr virus (EBV) is a member of the Gammaherpesviridae genus and is the causative agent of infectious mononucleosis (glandular fever) in humans and implicated in lymphoproliferative disorders, gastric, and nasopharyngeal cancer [1]. EBV preferentially infects B-lymphocytes, and it is estimated that greater than 90% of the world adult population has a latent asymptomatic EBV infection [2]. Among the lymphoproliferative diseases associated with EBV are Burkitt’s lymphoma, a quick-onset pediatric B lymphoid-associated carcinoma frequently found in sub-Saharan Africa [3], and Hodgkin’s lymphoma [4]. EBV infection has been closely associated with the development of the autoimmune disease multiple sclerosis (MS) [5,6]. The underlying mechanisms of these diseases is the expression of viral proteins that function to manipulate the host-cell innate immune response and, in particular, inhibit host-cell apoptosis.
Mitochondrial-mediated or intrinsic apoptosis is regulated by the B-cell lymphoma 2 (Bcl-2) family proteins. Activation of proapoptotic members of the Bcl-2 family leads to mitochondrial outer membrane permeabilization (MOMP) and release of cytochrome c and other factors into the cytoplasm that activate the cell-destroying caspase cascade [7,8,9]. The proapoptotic Bcl-2 members are antagonized by their interactions with their pro-survival Bcl-2 counterparts [9]. The mammalian pro-survival Bcl-2 proteins comprise Bcl-2, Bcl-xL, Bcl-w, Mcl-1, A1/Bfl1, and Bcl-B, which function to inhibit intrinsic apoptosis [8]. In contrast, proapoptotic Bcl-2 proteins categorized into two groups: multi-motif apoptosis effectors (Bak, Bax, Bok, and Bid) and BH3-only proteins (Bad, Bid, Bik, Bim, Bmf, Hrk, Noxa, and Puma) [8,10]. The expression levels and interaction between the antagonistic pro-survival and proapoptotic Bcl-2 proteins determines cellular fate [8,9]. The critical role of the Bcl-2 family in life and death decisions in the cell has led to their mimicry by viruses [11].
Numerous large double-stranded DNA viruses, including examples from Poxviridae, Herpesviridae, Iridoviridae, Adenoviridae, and Asfraviridae, have evolved sophisticated immunomodulatory strategies to counter host-cell innate and adaptive immune responses to viral infection. One such strategy is to block host-cell intracellular responses to infection such as apoptosis to viral proliferation, allowing infection to proceed [12,13,14]. This includes expression of direct inhibitors against caspases and mimicry of cellular Bcl-2 proteins to inactivate host proapoptotic proteins and block host-cell intrinsic apoptosis [14]. Viral Bcl-2 (vBcl-2) proteins are highly sequence divergent, but structurally and functionally are well-conserved homologs of cellular pro-survival Bcl-2 proteins [11,13]. Examples from Herpesviridae include those from murine gamma herpesvirus M11 [15], Karposi sarcoma herpes virus (KSHV) KsBcl-2 [16,17], and Epstein–Barr virus (EBV), which unusually bears two Bcl-2 homologs BHRF1 and BALF1 [18,19].
The genome of EBV encodes BHRF1, and elevated expression of BHRF1 leads to resistance to multiple cytotoxic agents [18,20]. During such overexpression of BHRF1, mouse progenitor cells and hematopoietic stem cells have an increased chance of MYC-induced lymphoma development in a mouse model [18]. EBV BHRF1 is a Bcl-2 homolog that potently inhibits host intrinsic apoptosis via direct inactivation of cellular proapoptotic Bcl-2 proteins Bak and Bim [21,22] and selectively interacts with a subset of proapoptotic proteins (Bid, Puma, and Bax) [23] to regulate apoptosis. Here, we report the structural investigation of BHRF1 bound to peptides that span the Bid BH3 and Puma BH3 motif. Our findings provide structural and mechanistic insight into EBV BHRF1-mediated inhibition of apoptosis.

2. Materials and Methods

2.1. Protein Expression and Purification

A synthetic cDNA codon-optimized for Escherichia coli encoding the construct EBV BHRF1ΔC31 (EBV strain B95-8, spanning residues 1–160 with a C-terminal deletion of 31 residues [18]; UniProt KB P03182 and referred to hereafter as BHRF1) was cloned into the bacterial expression vector pET-Duet1 (Genscript) and transformed into E. coli Codon plus (RIL). Protein expression and purification were carried out as previously described [18]. Briefly, BHRF1-expressing E. coli cells were resuspended in 50 mL of lysis buffer A (50 mM Tris, pH 8.5, 300 mM NaCl, and 20 mM imidazole). The cells were lysed using sonication (program 6 (6 s pulse-on, 60 s pulse-off) Model 705 Sonic Dismembrator, Fisher Scientific, Hampton, NH, USA), and the resultant lysate was transferred into SS34 tubes for further centrifugation at 20,000 rpm (JA-25.50 rotor, Beckman Coulter Avanti J-E, Beckman Coulter Australia) for 30 min. The supernatant was loaded onto a HisTrap HP column, 5 mL (GE Healthcare, Little Chalford, UK), equilibrated with buffer A. After sample application, the column was washed with 100 mL of buffer A, and the target protein was eluted with buffer B (50 mM Tris, pH 8.5, 300 mM NaCl, and 300 mM imidazole) and concentrated to a final volume of 5 mL using a centrifugal concentrator with a 3 kDa molecular weight cutoff (Amicon® Ultra 15, Merck Pty. Ltd., Melbourne, Australia). Concentrated BHRF1 protein was subjected to size-exclusion chromatography using a Superdex S75 16/60 column mounted on an ÄKTAExpress system (GE Healthcare) equilibrated in 25 mM HEPES pH 7.5, 150 mM NaCl, 5 mM TCEP (Tris(2-carboxyethyl)phosphine hydrochloride), where it eluted as a single peak in the volume corresponding to the monomeric nature. The final sample purity was estimated to be higher than 95% according to SDS-PAGE analysis.

2.2. Crystallization and Structure Determination

Crystals for BHRF1–BidBH3 and BHRF1–PumaBH3 complexes were obtained by mixing BHRF1 with human Bid (34-mer, 76-SESQEDIIRNIARHLAQVGDSMDRSIPPGLVNGL-109): Puma (26-mer, 130-EEQWAREIGAQLRRMADDLNAQYERR-155) peptide using a 1:1.5 molar ratio as described [24] and concentrated using a centrifugal concentrator with 3 kDa molecular weight cutoff (Amicon® Ultra 0.5, Merck Pty. Ltd., Melbourne, Australia) to 5 mg/mL. Concentrated protein was immediately used for crystallization trials. Initial high-throughput sparse matrix screening was performed using 96-well sitting drop trays (Swissic, Neuheim, Switzerland).
BHRF1–BidBH3 crystals were grown using the sitting drop vapor diffusion method at 20 °C in 0.1 M sodium acetate, pH 4.6, 8% (w/v) PEG 4000, or 0.4 M ammonium phosphate monobasic. The crystals were flash-cooled at −173 °C in mother liquor supplemented with 20% ethylene glycol. The BHRF1–BidBH3 complex formed single cuboidal crystals belonging to space group P6522 with a = 94.20 Å, b = 94.20 Å, c = 455.58 Å, α = 90.0°, β = 90.0°, and γ = 120.0° in the hexagonal crystal system.
Diffraction data were acquired at the Australian Synchrotron MX2 beamline [25] using an Eiger 16M detector with an oscillation range of 0.1° per frame and a wavelength of 0.9537 Å. The diffraction data were integrated using XDS [26] and scaled using AIMLESS [27]. A molecular replacement solution was obtained using PHASER as implemented in PHENIX [28], using the previously solved structure of BHRF1–BimBH3 (PDB ID 2WH6) [18] as the initial search model. BHRF1–BidBH3 crystals contained five molecules of BHRF1 and five BidBH3 peptides in the asymmetric unit, with a 53.7% solvent content and final TFZ and LLG values of 47.3 and 3316.2, respectively. The final model of the BHRF1–BidBH3 complex was built manually over several cycles using Coot [29] and refined using PHENIX [30] with final Rwork/Rfree values of 0.231/0.274, yielding a model where 99% of residues were in the favored region of the Ramachandran plot with no outliers.
BHRF1–PumaBH3 crystals were obtained in 0.4 M ammonium monobasic phosphate. The crystals were flash-cooled at −173 °C in mother liquor supplemented with 20% ethylene glycol as cryoprotectant. The complex formed single cubic crystals belonging to space group P3221 with a = 62.77 Å, b = 62.77 Å, c = 92.60 Å, α = 90.00°, β = 90.00°, and γ = 120.00° in the trigonal crystal system. Diffraction data collection, integration, and scaling were performed as described above. The molecular replacement was carried out using PHASER [28] as described above. BHRF1–PumaBH3 crystals contained a single molecule of BHRF1 and a single PumaBH3 peptide, with a 47.93% solvent content and final TFZ and LLG values of 34.5 and 4338.35, respectively. The final model of BHRF1–PumaBH3 was built manually over several cycles using Coot [29] and refined using PHENIX [30] with final Rwork/Rfree values of 21.7/25.2, yielding a model where 99% of residues were in the favored region of the Ramachandran plot with no outliers.
The X-ray crystallographic data collection and refinement statistics for the structures determined are presented in Table 1. Sequence and structural analysis was performed using the DALI server [31], and sequence statistics are presented in Table 2. In the absence of an experimental structure of F72W BHRF1, AlphaFold (version 2.1.0) was used to model the BHRF1 mutant that selectively binds Puma [20], using the Google Colab online server. The structure of F72W BHRF1(1–160) was calculated by providing the sequence with the standard parameters provided [32]. The model produced was built with the per-residue accuracy of the structure as determined by the predicted local distance difference test (pLDDT) at the >90% confidence level over the entire sequence, and it was further confirmed by superimposition with experimentally determined structure of BHRF1–Puma.
All images for BHRF1–BidBH3 and BHRF1–PumaBH3 complexes and superimpositions were generated using PyMOL molecular graphic system version 1.8.6.0 (Schrödinger, LLC, New York, NY, USA). The two structures were deposited at the Protein Data Bank with PDB accession codes 7P33 and 7P9W for the BidBH3 and PumaBH3 complexes, respectively. All raw images were deposited at the SBGridDB [33] using their PDB accession codes. All software was accessed through the SBGrid suite [34].

3. Results

It was previously reported that the EBV Bcl-2 homolog BHRF1 has a restricted BH3 interaction profile toward peptides spanning the BH3 motif (26-mer) of proapoptotic Bcl-2 proteins and binds only human Bak, Bax, Bim, Bid, and Puma BH3-motif peptides with nanomolar to sub-micromolar affinities [18]. To understand the molecular basis of the interaction between BHRF1 and Bid BH3 and Puma BH3 in modulating the inhibition of apoptosis, we determined the crystal structures of BHRF1 in complex with BidBH3 (Figure 1a) and PumaBH3 (Figure 1b). The structure of the BHRF1–BidBH3 complex was refined to a resolution of 2.78 Å, whereas the BHRF1–Puma BH3 complex was refined to 2.0 Å (Table 1) resolution. Both crystal structures of BHRF1 BH3 complexes were solved by molecular replacement using the previously solved structure of BHRF1–BimBH3 (PDB ID 2WH6) [18] as an initial search model. For the BHRF1–Bid BH3 complex, a clear and continuous electron density map was observed for BHRF1 residues 3–155 and BidBH3 motif residues 78–107 (Figure 2a). Despite sharing very weak sequence identity with its cellular counterparts Bcl-2 (21.7%), Bcl-xL (17.9%), Bcl-w (18.4%), Mcl-1 (9.6%), Bfl-1/A1 (17.8%), and Bcl-B (21.5%), BHRF1 adopts the conserved alpha helical Bcl-2 fold comprising eight α-helices (Figure 1a,b), with helices α2–α5 forming the canonical hydrophobic ligand-binding groove which provides the interaction site for the BH3-motif peptide (Figure 1a,b). The crystal structure of BHRF1 revealed a globular monomeric configuration that has been observed in other Bcl-2 homologs including human Bcl-xL (Figure 1c) and herpesvirus Bcl-2 homologs including KSHV KsBcl-2 (Figure 1d) [17], murine gamma herpes virus Bcl-2 homolog M11 [15], several poxviruses including sheep poxvirus SPPV14 [35,36], myxomavirus M11L [37], tanapoxvirus TANV16L [38] (Figure 1e) (which adopts both monomeric and dimeric configurations), fowl poxvirus FPV039 [39], canary poxvirus CNP058 [40], and African swine fever virus (ASFV) A179L [41,42]. A structural investigation and similarity analysis was performed using DALI server [31], which identified KSHV KsBcl-2 (PDB ID 1K3K, 7QTX) [17] as the closest viral Bcl-2 homolog of the BHRF1–Bid and BHRF1–Puma complexes with an RMSD of 2.8 Å over 116 Cα atoms and human Bcl-xL (PDB ID 2M04) (Figure 1c) [43] as the closest cellular Bcl-2 homolog with an RMSD of 2.5 Å over 133 Cα atoms where there was no significant difference in RMSD value observed for searches of the BH3-motif complexes. Bid BH3 is bound to BHRF1 via the canonical hydrophobic ligand-binding groove (Figure 1a). The four canonical BH3 motif-defining residues from Bid, I86, L90, V93, and M97, engage four hydrophobic pockets of the canonical BHRF1 ligand-binding groove with an additional hydrophobic pocket occupied by a fifth hydrophobic residue, I83 (Figure 2c). Furthermore, these hydrophobic interactions are supplemented with two ionic interactions (salt bridges) between D95Bid and R100BHRF1, the hallmark of the pro-survival Bcl-2 interaction between a conserved arginine in the BH1 motif and aspartate of the BH3 motif, which is found in almost all pro-survival Bcl-2–BH3 complexes [9] and R88Bid and E89BHRF1, as well as two hydrogen bonding interactions between D98Bid and G99BHRF1, and between R84Bid and E89BHRF1.
The BHRF1–Puma BH3 complex showed a clear and continuous electron density map for BHRF1 residues 2–94 and 96–156, as well as Puma BH3 motif peptide residues 131–151 (Figure 2b). Similar to the BHRF1–Bid BH3 complex, BHRF1–Puma BH3 adopted a globular monomeric topology (Figure 1b). The four conserved hydrophobic residues of the Puma BH3 motif, I137, L141 M144, and L148, protrude into four hydrophobic pockets in the canonical ligand-binding groove, with an additional hydrophobic pocket occupied by a fifth hydrophobic residue W133 (Figure 2d). Furthermore, a conserved ionic interaction was also found in the BHRF1–Puma BH3 complex between R100BHRF1 and D146Puma. These salt bridges were supplemented by additional hydrogen bonds between the carboxyl group of E89BHRF1 and the NδH side-chain group of R135Puma, and between the G99BHRF1 main-chain amide group and the side-chain carboxyl and amide group of N149Puma (Figure 2d).
In a search for BHRF1 variants with restricted ligand selectivity, Fitzsimmons et al. [20] investigated F72W BHRF1 and showed that it selectively bound Puma but much more weaky than the native sequence, and that it did not selectively bind Bid (Figure 3a). The lowered affinity of Puma is consistent with biological experiments showing that F72W BHRF1 does not protect cells from Puma-dependent apoptotic pathways [20]. In order to examine the possible molecular reasons for the alteration in the BH3 binding specificity for the F72W BHRF1 mutant, we generated an AlphaFold model [32]. The F72W BHRF1(1–160) AlphaFold model superimposition with the experimentally determined structure of BHRF1–PumaBH3 resulted in an RMSD value of 0.48 Å, confirming that the AlphaFold-generated structure was nearly identical to the experimental structures over the backbone atoms. Pocket 5 of F72W BHRF1 is shown in Figure 3b, and the main differences are small changes in the pocket including the presence of a hydrogen bond between T76 of the α3–α4 connecting loop and W72, in addition to a conformational change in the α3–α4 loop. A combination of altered steric and electrostatic interactions is likely responsible for the reduction in the ability of pocket 5 to accommodate a residue from the BH3 peptide. Figure 3b shows pocket 5 for F72W BHRF1, and Figure 3c,d show pocket 5 for F72W BHRF1 and BHRF1–BidBH3, respectively, colored by hydrophobicity. Figure 3e,f show the same region in the liganded forms of BHRF1 with BidBH3 and PumaBH3, respectively.

4. Discussion

Viruses have evolved countermeasures to hijack host-cell suicide to ensure cellular persistence and optimize viral replication [11,13,14]. Many viruses regulate intrinsic apoptosis through expressing homologs of the Bcl-2 family [11], including the oncogenic herpesvirus EBV. BHRF1, one of the two Bcl-2 homologs in the EBV genome, is known to stimulate tumorigenesis by inactivating the cellular proapoptotic proteins Bak, Bax, Bim, Bid, and Puma. Furthermore, this inhibition drives resistance to chemotherapy [18,20]. Here, we determined the crystal structures of EBV BHRF1 in complex with the BH3 motif peptides of human proapoptotic proteins Bid and Puma, two key BH3-only proteins in Bcl-2-regulated apoptosis. In contrast to Puma, Bid is activated by proteolytic cleavage by caspase-8 during death receptor-mediated apoptosis that results in a fragment, truncated Bid (tBid), bearing the BH3 motif that interacts with pro-survival Bcl-2 proteins [44]; however, like Bax and Bak, it can also elicit MOMP [45] and, therefore, links cell surface receptor-initiated apoptosis (extrinsic apoptosis) with intrinsic apoptosis. Bid interacts with five cellular pro-survival proteins (Bcl-2, Bcl-xL, Bcl-w, Mcl-1, and A1) to neutralize their activity but not with Bcl-B [46,47], and it additionally interacts with most viral Bcl-2 proteins with high affinity [23]. Although the primary protein sequence of BHRF1 significantly varies from its cellular Bcl-2 homologous counterparts (Table 2), it is highly conserved at a structural level (Figure 4). A detailed understanding of Bcl-2 family protein function in apoptosis regulation is not only crucial for identifying their biological function, but is important in the development of new therapeutic targets directed against this family [47,48,49]. Indeed, as a major mediator of programmed cell death, there is substantial interest in resolving the role of Bcl-2 family proteins at a molecular level with the aim of targeting them for their role in cancer progression [9,50].
Structural analysis of Bid and PumaBH3 motif interactions with BHRF1 presented here showed that binding occurred in a similar manner to what has previously been reported for viral and mammalian pro-survival Bcl-2–BH3 complexes [23]. A comparison of the BHRF1–Bid BH3 complex with KSHV KsBcl-2–Bid BH3 [17] or VARV F1L–Bid BH3 [51] shows that their mode of BH3 interaction is nearly identical, despite having significantly different affinities (110 nM, 20 nM, and 3220 nM toward Bid BH3, respectively). A similar trend can be observed in BHRF1–Puma BH3 and KsBcl-2–Puma BH3 [17] complexes, where BHRF1 shows nearly threefold reduced affinity for Puma BH3. Superimposition of the BHRF1 complex structures BHRF1–BidBH3 and BHRF1–Puma BH3 gave an RMSD value of 0.5 Å over the entire BHRF1 backbone, indicating the close structural similarity of the complexes. Pairwise analysis showed that the closest structural neighbors of the BHRF1–Bid and BHRF1–Puma complexes were KsBcl-2 from the herpesvirus KSHV, while the nearest mammalian homolog was human Bcl-xL. Although the Bcl-2 fold was highly similar (Figure 4), the BH3 binding selectivity differed (Figure 3a) [18].
A comparison between BHRF1 interactions with BidBH3 or PumaBH3 identified three conserved BHRF1 residues, E89, G99, and R100, in the binding groove which formed polar interactions with the BH3 motif peptides Bid or Puma (Figure 2 and Figure 4). Among these, E89 and R100 formed ionic interactions, and R100 made the highly conserved canonical salt bridge interaction with the conserved aspartic acid residue from the BH3 motif. This salt bridge is a hallmark interaction between pro-survival Bcl-2 protein and the BH3 motif of proapoptotic Bcl-2 protein [9], conserved across evolution [52,53], also reported for the BHRF1–BimBH3 and BHRF1–BakBH3 complexes [18]. The BH1 motif is located at the N-terminus of α5, and, although it generally conforms to the sequence XXGR, where XX are most commonly NW, but may vary, with the GR dipeptide very strongly conserved. In BHRF1, this BH1 region is SLGR, and the glycine and arginine residues make polar contacts with both Bid BH3 and Puma BH3 peptides. Site-directed mutagenesis studies in the BH1 motif of human pro-survival Bcl-2 proteins abolished BH3 binding [47], and similar behavior was reported with BHRF1. A charge reversal mutant in BHRF1 in the BH1 region, BHRF1 R100D, lost its affinity for Bid, Bak, and Bax and significantly lowered the affinity for BimBH3 (~15-fold reduction) and PumaBH3 (~6-fold reduction) [20]. These data suggest that residue R100 of BHRF1 is a crucial mediator for BH3 binding. Human Bcl-xL and KSHV KsBcl2 have conserved BH1 NWGR motifs and show comparatively high affinities for both Bid and Puma [17,47]. In contrast, the structurally equivalent region to the NWGR motif in VARV F1L is replaced by the sequence LGVR, and the canonical salt bridge interaction is absent [51]. These observations highlight the importance of the BH1 equivalent region for determining interactions between Bcl-2 family proteins and provide a structural basis for the selectivity variation observed for proapoptotic BH3 ligands.
Structural analysis of the BHRF1 complexes showed the importance of a fifth hydrophobic pocket for BHRF1 BH3 ligand binding. Figure 4 presents a comparison of BidBH3 (Figure 4a–d) and Puma BH3 (Figure 4e–h) complexes. Examination of the BHRF1 and Bcl-xL complexes liganded with Bid BH3 (Figure 4a,b) revealed that Bid BH3 utilized the canonical four conserved hydrophobic residues of the BH3 motif peptides to interact with four hydrophobic pockets of the Bcl-2 binding groove. In contrast to Bcl-xL–Bid BH3 (Figure 4b) and the previously determined structures of BHRF1–Bim and Bak BH3 motif complexes [18], Bid BH3 engages a fifth hydrophobic pocket of BHRF1 (Figure 3 and Figure 4a). Residue I83Bid occupies a pocket formed by BHRF1 residues F72, T76, and V79 (Figure 3e), all residues highly conserved or conservatively substituted in BHRF1 variants. The use of a fifth hydrophobic pocket has previously been observed in several other viral Bcl-2–Bid BH3 complexes such as KsBcl-2 (Figure 4c) and the poxvirus variola virus F1L (Figure 4d) [17,51]. The fifth hydrophobic pocket of BHRF1 is utilized by Puma BH3, where the structurally equivalent residue to Bid I83 is A134 or W133 (Figure 3f and Figure 4e). The binding mode for Bcl-xL–Puma BH3 (Figure 4f) differs from that of the BHRF1 complex. The fifth pocket has previously been observed in the KSHV KsBcl-2 Puma complex (Figure 4g) and tanapox 16L, showing the variation found within the exact interactions that the viral Bcl-2 proteins can form in binding the same ligand (Figure 4h). Binding studies have emphasized the importance of the fifth pocket to BHRF1 to ligand selectivity and affinity. Fitzsimmons et al. showed that F72W BHRF1 lost binding to almost all BH3 peptides except for PumaBH3, which showed a ~4-fold reduction in affinity (Figure 3a) [20]. The predicted changes rendered in pocket 5 by the presence of the Trp residue include changes in backbone α3–α4 loop conformation, as well as an altered conformation of Thr76, and these facets change the hydrophobicity of the pocket (Figure 3c,d). Collectively, these changes have the effect of disfavoring Bid binding but only reducing Puma BH3 affinity (Figure 3a). An alternative fifth hydrophobic pocket has been observed for the orfvirus Bcl-2 protein Puma complex, ORFV125–Puma BH3, using residue Y152Puma [54], as well as for Bax BH3 complexes where the residue Bax BH3 M74 uses a fifth pocket in complexes with Mcl-1, Bcl-xL [55], Bcl-2 [56], ORFV125 [57], and SPPV14 [36]. Collectively, these data suggest that virus-encoded pro-survival Bcl-2 proteins make extensive use of an additional hydrophobic pocket to bind select host pro-death Bcl-2 proteins, which has implications for selectivity and affinity for host pro-death Bcl-2 proteins by their viral Bcl-2 antagonists. Conversely, obstruction of one of the four canonical hydrophobic pockets is employed by the ancient pro-survival Bcl-2 protein BHP2 from a sponge [58]. Clearly, the flexibility within the Bcl-2 fold enables multiple avenues to achieve BH3 interactor selectivity. Intriguingly, the most promiscuous host pro-death Bcl-2 protein Bim has not been shown to use a fifth hydrophobic pocket whilst typically displaying high affinities, but it can be bound with high selectivity by several virus-encoded Bcl-2 proteins including vaccinia virus F1L [59] and grouper iridovirus GIV66 [60], as well as by human Bcl-b [46].
Apoptosis induction by viral infection has been shown to be dependent on Puma for apoptosis via Bax/Bak activation [61]; viral infection leads to caspase-8 activation and, hence, Bid cleavage and activation [14]. Thus, inhibiting both Puma and Bid can play significant roles in suppressing virus induced apoptosis. Understanding the precise atomic details of the BHRF1 Puma and Bid complexes will potentially allow us to target the interactions in a molecular design process. The exploitation of the fifth hydrophobic pocket of BHRF1 presents a unique opportunity for designing specific inhibitors. Inhibition of BHRF1 has been shown to induce apoptosis in EBV-infected cells [62], holding out the prospect that small-molecule inhibitors of BHRF1 may be of utility in inducing apoptosis in EBV-infected cells. Our findings may provide important targets for the design of peptidomimetic drugs to mimic these interactions, and the structures we present here provide the mechanistic basis that could be potentially exploited for targeting EBV infections.

Author Contributions

C.D.S., acquisition, analysis, and interpretation of data, and drafting and revising the article; M.G.H., conceptualization and design, analysis and interpretation of data, and drafting and revising the article; M.K., conceptualization and design, acquisition, analysis, and interpretation of data, and drafting and revising the article. All authors read and agreed to the published version of the manuscript.

Funding

This research was funded by the Australian Research Council (Fellowship FT130101349 to MK) and La Trobe University (scholarship to CDS).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data supporting the findings of this manuscript are available from the corresponding authors upon reasonable request. Coordinate files were deposited at the Protein Data Bank (https://www.rcsb.org/, accessed on 13 March 2021) using accession codes 7P33 and 7P9W for the BHRF1–BidBH3 and BHRF1–PumaBH3 complexes, respectively. The raw X-ray diffraction data were deposited at the SBGrid Data Bank [33] (https://data.sbgrid.org/data/, accessed on 1 January 2020) using their PDB accession codes.

Acknowledgments

The authors thank staff at the MX beamlines at the Australian Synchrotron for help with X-ray data collection. The authors also thank the ACRF for their support of the Eiger MX detector at the Australian Synchrotron MX2 beamline and the Comprehensive Proteomics Platform at La Trobe University for core instrument support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Young, L.S.; Yap, L.F.; Murray, P.G. Epstein-Barr virus: More than 50 years old and still providing surprises. Nat. Rev. Cancer 2016, 16, 789–802. [Google Scholar] [CrossRef] [PubMed]
  2. Fugl, A.; Andersen, C.L. Epstein-Barr virus and its association with disease—A review of relevance to general practice. BMC Fam. Pract. 2019, 20, 62. [Google Scholar] [CrossRef] [PubMed]
  3. Coghill, A.E.; Proietti, C.; Liu, Z.; Krause, L.; Bethony, J.; Prokunina-Olsson, L.; Obajemu, A.; Nkrumah, F.; Biggar, R.J.; Bhatia, K.; et al. The Association between the Comprehensive Epstein–Barr Virus Serologic Profile and Endemic Burkitt Lymphoma. Cancer Epidemiol. Biomark. Prev. 2020, 29, 57–62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mundo, L.; Del Porro, L.; Granai, M.; Siciliano, M.C.; Mancini, V.; Santi, R.; Marcar, L.; Vrzalikova, K.; Vergoni, F.; Di Stefano, G.; et al. Frequent traces of EBV infection in Hodgkin and non-Hodgkin lymphomas classified as EBV-negative by routine methods: Expanding the landscape of EBV-related lymphomas. Mod. Pathol. 2020, 33, 2407–2421. [Google Scholar] [CrossRef] [PubMed]
  5. Lanz, T.V.; Brewer, R.C.; Ho, P.P.; Moon, J.S.; Jude, K.M.; Fernandez, D.; Fernandes, R.A.; Gomez, A.M.; Nadj, G.S.; Bartley, C.M.; et al. Clonally expanded B cells in multiple sclerosis bind EBV EBNA1 and GlialCAM. Nature 2022, 603, 321–327. [Google Scholar] [CrossRef]
  6. Bjornevik, K.; Cortese, M.; Healy, B.C.; Kuhle, J.; Mina, M.J.; Leng, Y.; Elledge, S.J.; Niebuhr, D.W.; Scher, A.I.; Munger, K.L.; et al. Longitudinal analysis reveals high prevalence of Epstein-Barr virus associated with multiple sclerosis. Science 2022, 375, 296–301. [Google Scholar] [CrossRef]
  7. Green, D.R. The Mitochondrial Pathway of Apoptosis: Part I: MOMP and Beyond. Cold Spring Harb. Perspect. Biol. 2022, 14, a041038. [Google Scholar] [CrossRef]
  8. Green, D.R. The Mitochondrial Pathway of Apoptosis Part II: The BCL-2 Protein Family. Cold Spring Harb. Perspect. Biol. 2022, 14, a041046. [Google Scholar] [CrossRef]
  9. Banjara, S.; Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. The Bcl-2 Family: Ancient Origins, Conserved Structures, and Divergent Mechanisms. Biomolecules 2020, 10, 128. [Google Scholar] [CrossRef] [Green Version]
  10. Ke, F.S.; Holloway, S.; Uren, R.T.; Wong, A.W.; Little, M.H.; Kluck, R.M.; Voss, A.K.; Strasser, A. The BCL-2 family member BID plays a role during embryonic development in addition to its BH3-only protein function by acting in parallel to BAX, BAK and BOK. EMBO J. 2022, 41, e110300. [Google Scholar] [CrossRef]
  11. Kvansakul, M.; Caria, S.; Hinds, M.G. The Bcl-2 Family in Host-Virus Interactions. Viruses 2017, 9, 290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Youle, R.J.; Strasser, A. The BCL-2 protein family: Opposing activities that mediate cell death. Nat. Rev. Mol. Cell Biol. 2008, 9, 47–59. [Google Scholar] [CrossRef] [PubMed]
  13. Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. Poxviral Strategies to Overcome Host Cell Apoptosis. Pathogens 2021, 10, 6. [Google Scholar] [CrossRef] [PubMed]
  14. Galluzzi, L.; Brenner, C.; Morselli, E.; Touat, Z.; Kroemer, G. Viral control of mitochondrial apoptosis. PLoS Pathog. 2008, 4, e1000018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Sinha, S.; Colbert, C.L.; Becker, N.; Wei, Y.; Levine, B. Molecular basis of the regulation of Beclin 1-dependent autophagy by the gamma-herpesvirus 68 Bcl-2 homolog M11. Autophagy 2008, 4, 989–997. [Google Scholar] [CrossRef] [PubMed]
  16. Sarid, R.; Sato, T.; Bohenzky, R.A.; Russo, J.J.; Chang, Y. Kaposi’s sarcoma-associated herpesvirus encodes a functional Bcl-2 homologue. Nat. Med. 1997, 3, 293–298. [Google Scholar] [CrossRef] [PubMed]
  17. Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. Structural Insight into KsBcl-2 Mediated Apoptosis Inhibition by Kaposi Sarcoma Associated Herpes Virus. Viruses 2022, 14, 738. [Google Scholar] [CrossRef] [PubMed]
  18. Kvansakul, M.; Wei, A.H.; Fletcher, J.I.; Willis, S.N.; Chen, L.; Roberts, A.W.; Huang, D.C.; Colman, P.M. Structural basis for apoptosis inhibition by Epstein-Barr virus BHRF1. PLoS Pathog. 2010, 6, e1001236. [Google Scholar] [CrossRef] [Green Version]
  19. Fitzsimmons, L.; Kelly, G.L. EBV and Apoptosis: The Viral Master Regulator of Cell Fate? Viruses 2017, 9, 339. [Google Scholar] [CrossRef] [Green Version]
  20. Fitzsimmons, L.; Cartlidge, R.; Chang, C.; Sejic, N.; Galbraith, L.C.A.; Suraweera, C.D.; Croom-Carter, D.; Dewson, G.; Tierney, R.J.; Bell, A.I.; et al. EBV BCL-2 homologue BHRF1 drives chemoresistance and lymphomagenesis by inhibiting multiple cellular pro-apoptotic proteins. Cell Death Differ. 2020, 27, 1554–1568. [Google Scholar] [CrossRef]
  21. Desbien, A.L.; Kappler, J.W.; Marrack, P. The Epstein-Barr virus Bcl-2 homolog, BHRF1, blocks apoptosis by binding to a limited amount of Bim. Proc. Natl. Acad. Sci. USA 2009, 106, 5663–5668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Flanagan, A.M.; Letai, A. BH3 domains define selective inhibitory interactions with BHRF-1 and KSHV BCL-2. Cell Death Differ. 2008, 15, 580–588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Kvansakul, M.; Hinds, M.G. Structural biology of the Bcl-2 family and its mimicry by viral proteins. Cell Death Dis. 2013, 4, e909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kvansakul, M.; Czabotar, P.E. Preparing Samples for Crystallization of Bcl-2 Family Complexes. Methods Mol. Biol. 2016, 1419, 213–229. [Google Scholar] [CrossRef] [PubMed]
  25. Aragão, D.; Aishima, J.; Cherukuvada, H.; Clarken, R.; Clift, M.; Cowieson, N.P.; Ericsson, D.J.; Gee, C.L.; Macedo, S.; Mudie, N.; et al. MX2: A high-flux undulator microfocus beamline serving both the chemical and macromolecular crystallography communities at the Australian Synchrotron. J. Synchrotron Radiat. 2018, 25, 885–891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kabsch, W. XDS. Acta Cryst. D Biol. Cryst. 2010, 66, 125–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Evans, P. Scaling and assessment of data quality. Acta Crystallogr. Sect. D Biol. Crystallogr. 2006, 62, 72–82. [Google Scholar] [CrossRef]
  28. McCoy, A.J. Solving structures of protein complexes by molecular replacement with Phaser. Acta Crystallogr. Sect. D Biol. Crystallogr. 2007, 63, 32–41. [Google Scholar] [CrossRef] [Green Version]
  29. Emsley, P.; Lohkamp, B.; Scott, W.G.; Cowtan, K. Features and development of Coot. Acta Crystallogr. Sect. D Biol. Crystallogr. 2010, 66, 486–501. [Google Scholar] [CrossRef] [Green Version]
  30. Afonine, P.V.; Grosse-Kunstleve, R.W.; Echols, N.; Headd, J.J.; Moriarty, N.W.; Mustyakimov, M.; Terwilliger, T.C.; Urzhumtsev, A.; Zwart, P.H.; Adams, P.D. Towards automated crystallographic structure refinement with phenix.refine. Acta Crystallogr. D Biol. Crystallogr. 2012, 68, 352–367. [Google Scholar] [CrossRef]
  31. Holm, L.; Rosenström, P. Dali server: Conservation mapping in 3D. Nucleic Acids Res. 2010, 38, W545–W549. [Google Scholar] [CrossRef] [PubMed]
  32. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef] [PubMed]
  33. Meyer, P.A.; Socias, S.; Key, J.; Ransey, E.; Tjon, E.C.; Buschiazzo, A.; Lei, M.; Botka, C.; Withrow, J.; Neau, D.; et al. Data publication with the structural biology data grid supports live analysis. Nat. Commun. 2016, 7, 10882. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Morin, A.; Eisenbraun, B.; Key, J.; Sanschagrin, P.C.; Timony, M.A.; Ottaviano, M.; Sliz, P. Collaboration gets the most out of software. Elife 2013, 2, e01456. [Google Scholar] [CrossRef] [PubMed]
  35. Okamoto, T.; Campbell, S.; Mehta, N.; Thibault, J.; Colman, P.M.; Barry, M.; Huang, D.C.; Kvansakul, M. Sheeppox virus SPPV14 encodes a Bcl-2-like cell death inhibitor that counters a distinct set of mammalian proapoptotic proteins. J. Virol. 2012, 86, 11501–11511. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Suraweera, C.D.; Burton, D.R.; Hinds, M.G.; Kvansakul, M. Crystal structures of the sheeppox virus encoded inhibitor of apoptosis SPPV14 bound to the proapoptotic BH3 peptides Hrk and Bax. FEBS Lett. 2020, 594, 2016–2026. [Google Scholar] [CrossRef]
  37. Kvansakul, M.; van Delft, M.F.; Lee, E.F.; Gulbis, J.M.; Fairlie, W.D.; Huang, D.C.; Colman, P.M. A structural viral mimic of prosurvival Bcl-2: A pivotal role for sequestering proapoptotic Bax and Bak. Mol. Cell 2007, 25, 933–942. [Google Scholar] [CrossRef] [PubMed]
  38. Suraweera, C.D.; Anasir, M.I.; Chugh, S.; Javorsky, A.; Impey, R.E.; Hasan Zadeh, M.; Soares da Costa, T.P.; Hinds, M.G.; Kvansakul, M. Structural insight into tanapoxvirus-mediated inhibition of apoptosis. FEBS J. 2020, 287, 3733–3750. [Google Scholar] [CrossRef] [PubMed]
  39. Anasir, M.I.; Caria, S.; Skinner, M.A.; Kvansakul, M. Structural basis of apoptosis inhibition by the fowlpox virus protein FPV039. J. Biol. Chem. 2017, 292, 9010–9021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Anasir, M.I.; Baxter, A.A.; Poon, I.K.H.; Hulett, M.D.; Kvansakul, M. Structural and Functional Insight into Canarypox Virus CNP058 Mediated Regulation of Apoptosis. Viruses 2017, 9, 305. [Google Scholar] [CrossRef]
  41. Banjara, S.; Caria, S.; Dixon, L.K.; Hinds, M.G.; Kvansakul, M. Structural Insight into African Swine Fever Virus A179L-Mediated Inhibition of Apoptosis. J. Virol. 2017, 91, e02228-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Banjara, S.; Shimmon, G.L.; Dixon, L.K.; Netherton, C.L.; Hinds, M.G.; Kvansakul, M. Crystal Structure of African Swine Fever Virus A179L with the Autophagy Regulator Beclin. Viruses 2019, 11, 789. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Follis, A.V.; Chipuk, J.E.; Fisher, J.C.; Yun, M.K.; Grace, C.R.; Nourse, A.; Baran, K.; Ou, L.; Min, L.; White, S.W.; et al. PUMA binding induces partial unfolding within BCL-xL to disrupt p53 binding and promote apoptosis. Nat. Chem. Biol. 2013, 9, 163–168. [Google Scholar] [CrossRef] [Green Version]
  44. Billen, L.P.; Shamas-Din, A.; Andrews, D.W. Bid: A Bax-like BH3 protein. Oncogene 2008, 27, S93–S104. [Google Scholar] [CrossRef] [Green Version]
  45. Flores-Romero, H.; Hohorst, L.; John, M.; Albert, M.C.; King, L.E.; Beckmann, L.; Szabo, T.; Hertlein, V.; Luo, X.; Villunger, A.; et al. BCL-2-family protein tBID can act as a BAX-like effector of apoptosis. EMBO J. 2022, 41, e108690. [Google Scholar] [CrossRef]
  46. Rautureau, G.J.P.; Yabal, M.; Yang, H.; Huang, D.C.S.; Kvansakul, M.; Hinds, M.G. The restricted binding repertoire of Bcl-B leaves Bim as the universal BH3-only prosurvival Bcl-2 protein antagonist. Cell Death Dis. 2012, 3, e443. [Google Scholar] [CrossRef] [Green Version]
  47. Kvansakul, M.; Hinds, M.G. The Bcl-2 family: Structures, interactions and targets for drug discovery. Apoptosis 2015, 20, 136–150. [Google Scholar] [CrossRef]
  48. Roberts, A.W. Therapeutic development and current uses of BCL-2 inhibition. Hematol. Am. Soc. Hematol. Educ. Program. 2020, 2020, 1–9. [Google Scholar] [CrossRef]
  49. Cory, S.; Roberts, A.W.; Colman, P.M.; Adams, J.M. Targeting BCL-2-like Proteins to Kill Cancer Cells. Trends Cancer 2016, 2, 443–460. [Google Scholar] [CrossRef]
  50. Fuchs, Y.; Steller, H. Live to die another way: Modes of programmed cell death and the signals emanating from dying cells. Nat. Rev. Mol. Cell Biol. 2015, 16, 329–344. [Google Scholar] [CrossRef]
  51. Marshall, B.; Puthalakath, H.; Caria, S.; Chugh, S.; Doerflinger, M.; Colman, P.M.; Kvansakul, M. Variola virus F1L is a Bcl-2-like protein that unlike its vaccinia virus counterpart inhibits apoptosis independent of Bim. Cell Death Dis. 2015, 6, e1680. [Google Scholar] [CrossRef]
  52. Popgeorgiev, N.; Sa, J.D.; Jabbour, L.; Banjara, S.; Nguyen, T.T.M.; Akhavan, E.S.A.; Gadet, R.; Ralchev, N.; Manon, S.; Hinds, M.G.; et al. Ancient and conserved functional interplay between Bcl-2 family proteins in the mitochondrial pathway of apoptosis. Sci. Adv. 2020, 6, eabc4149. [Google Scholar] [CrossRef]
  53. Suraweera, C.D.; Banjara, S.; Hinds, M.G.; Kvansakul, M. Metazoans and Intrinsic Apoptosis: An Evolutionary Analysis of the Bcl-2 Family. Int. J. Mol. Sci. 2022, 23, 3691. [Google Scholar] [CrossRef] [PubMed]
  54. Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. Structural Investigation of Orf Virus Bcl-2 Homolog ORFV125 Interactions with BH3-Motifs from BH3-Only Proteins Puma and Hrk. Viruses 2021, 13, 1374. [Google Scholar] [CrossRef] [PubMed]
  55. Czabotar, P.E.; Lee, E.F.; Thompson, G.V.; Wardak, A.Z.; Fairlie, W.D.; Colman, P.M. Mutation to Bax beyond the BH3 domain disrupts interactions with pro-survival proteins and promotes apoptosis. J. Biol. Chem. 2011, 286, 7123–7131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Ku, B.; Liang, C.; Jung, J.U.; Oh, B.H. Evidence that inhibition of BAX activation by BCL-2 involves its tight and preferential interaction with the BH3 domain of BAX. Cell Res. 2011, 21, 627–641. [Google Scholar] [CrossRef]
  57. Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. Crystal structures of ORFV125 provide insight into orf virus-mediated inhibition of apoptosis. Biochem. J. 2020, 477, 4527–4541. [Google Scholar] [CrossRef]
  58. Caria, S.; Hinds, M.G.; Kvansakul, M. Structural insight into an evolutionarily ancient programmed cell death regulator—The crystal structure of marine sponge BHP2 bound to LB-Bak-2. Cell Death Dis. 2017, 8, e2543. [Google Scholar] [CrossRef]
  59. Campbell, S.; Thibault, J.; Mehta, N.; Colman, P.M.; Barry, M.; Kvansakul, M. Structural Insight into BH3 Domain Binding of Vaccinia Virus Antiapoptotic F1L. J. Virol. 2014, 88, 8667. [Google Scholar] [CrossRef] [Green Version]
  60. Banjara, S.; Mao, J.; Ryan, T.M.; Caria, S.; Kvansakul, M. Grouper iridovirus GIV66 is a Bcl-2 protein that inhibits apoptosis by exclusively sequestering Bim. J. Biol. Chem. 2018, 293, 5464–5477. [Google Scholar] [CrossRef]
  61. Papaianni, E.; El Maadidi, S.; Schejtman, A.; Neumann, S.; Maurer, U.; Marino-Merlo, F.; Mastino, A.; Borner, C. Phylogenetically Distant Viruses Use the Same BH3-Only Protein Puma to Trigger Bax/Bak-Dependent Apoptosis of Infected Mouse and Human Cells. PLoS ONE 2015, 10, e0126645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Procko, E.; Berguig, G.Y.; Shen, B.W.; Song, Y.; Frayo, S.; Convertine, A.J.; Margineantu, D.; Booth, G.; Correia, B.E.; Cheng, Y.; et al. A computationally designed inhibitor of an Epstein-Barr viral Bcl-2 protein induces apoptosis in infected cells. Cell 2014, 157, 1644–1656. [Google Scholar] [CrossRef] [PubMed]
Viruses 14 02222 g001 550
Figure 1. EBV BHRF1 binds BH3 motif peptides of proapoptotic Bcl-2 proteins using the canonical ligand-binding groove. (a) Ribbon representation of the crystal structure of BHRF1–Bid BH3 (PDB ID 7P33) complex and its 180° rotation. BHRF1 is shown in salmon, while BidBH3 is shown in orange, with BHRF1 helices labeled α1–α8. The view is into the hydrophobic binding groove formed by helices α2–α5. (b) Ribbon representation of the crystal structure of BHRF1–Puma BH3 (PDB ID 7P9W) complex. BHRF1 is shown as the salmon cartoon, while PumaBH3 is shown in purple. The view is as in (a). (c) Cartoon of Bcl-xL–BidBH3 complex. Human Bcl-xL is shown in light blue, while Bid is shown in orange. (d) Cartoon of KSHV KsBcl-2–Puma BH3 complex. KsBcl-2 is shown as an olive-colored cartoon, while Puma is shown in purple. I Cartoon of tanapoxvirus TANV16L–Puma BH3 complex. TANV16L is shown as a red-colored cartoon, while Puma is shown in purple. The views in (c–e) are as in (a).
Figure 1. EBV BHRF1 binds BH3 motif peptides of proapoptotic Bcl-2 proteins using the canonical ligand-binding groove. (a) Ribbon representation of the crystal structure of BHRF1–Bid BH3 (PDB ID 7P33) complex and its 180° rotation. BHRF1 is shown in salmon, while BidBH3 is shown in orange, with BHRF1 helices labeled α1–α8. The view is into the hydrophobic binding groove formed by helices α2–α5. (b) Ribbon representation of the crystal structure of BHRF1–Puma BH3 (PDB ID 7P9W) complex. BHRF1 is shown as the salmon cartoon, while PumaBH3 is shown in purple. The view is as in (a). (c) Cartoon of Bcl-xL–BidBH3 complex. Human Bcl-xL is shown in light blue, while Bid is shown in orange. (d) Cartoon of KSHV KsBcl-2–Puma BH3 complex. KsBcl-2 is shown as an olive-colored cartoon, while Puma is shown in purple. I Cartoon of tanapoxvirus TANV16L–Puma BH3 complex. TANV16L is shown as a red-colored cartoon, while Puma is shown in purple. The views in (c–e) are as in (a).
Viruses 14 02222 g001
Viruses 14 02222 g002 550
Figure 2. Electron density maps of BHRF1–Bid BH3 and Puma BH3 complexes and detailed view of the BHRF1–Bid BH3 and Puma BH3 interfaces. (a) The 2Fo-Fc electron density maps of BHRF1–Bid BH3 complex interface contoured at 1.0 σ. BHRF1 is shown as salmon ribbon, while Bid BH3 is shown as orange-colored sticks. (b) The 2Fo-Fc electron density maps of BHRF1–PumaBH3 complex interface contoured at 1.0 σ. BHRF1 is shown as a salmon ribbon, while Puma BH3 is shown as purple-colored sticks. Electron density distribution is shown in gray. (c) The BHRF1 backbone, floor of the binding groove, and surface are shown in salmon, yellow, and gray, respectively, whilst the Bid BH3 ribbon is shown in orange. The five hydrophobic residues of Bid BH3 (I83, I86, L90, V93, and M97) are protruding into the binding groove, and conserved ionic interactions formed by BHRF1 R100 and Bid BH3 D95 are labeled, as well as all other residues involved in ionic interactions and hydrogen bonds between protein and peptide. (d) BHRF1 is shown as in (c), whilst the Puma BH3 ribbon is shown in purple. The five hydrophobic residues of Puma BH3 (W133, I137, L141 M144, and L148) are protruding into the binding groove, and the conserved salt bridge interaction formed by BHRF1 R100 and Puma BH3 D146 is labeled, as well as other residues involved in additional ionic interactions and hydrogen bonds. Ionic interactions and hydrogen bonds are shown as dotted black lines, and the participating residues are labeled. The canonical hydrophobic binding pockets h1–h4 are indicated by black circles, while pocket h5 is indicated by a red circle. The molecular orientations depicted in (a–d) are identical to those in Figure 1.
Figure 2. Electron density maps of BHRF1–Bid BH3 and Puma BH3 complexes and detailed view of the BHRF1–Bid BH3 and Puma BH3 interfaces. (a) The 2Fo-Fc electron density maps of BHRF1–Bid BH3 complex interface contoured at 1.0 σ. BHRF1 is shown as salmon ribbon, while Bid BH3 is shown as orange-colored sticks. (b) The 2Fo-Fc electron density maps of BHRF1–PumaBH3 complex interface contoured at 1.0 σ. BHRF1 is shown as a salmon ribbon, while Puma BH3 is shown as purple-colored sticks. Electron density distribution is shown in gray. (c) The BHRF1 backbone, floor of the binding groove, and surface are shown in salmon, yellow, and gray, respectively, whilst the Bid BH3 ribbon is shown in orange. The five hydrophobic residues of Bid BH3 (I83, I86, L90, V93, and M97) are protruding into the binding groove, and conserved ionic interactions formed by BHRF1 R100 and Bid BH3 D95 are labeled, as well as all other residues involved in ionic interactions and hydrogen bonds between protein and peptide. (d) BHRF1 is shown as in (c), whilst the Puma BH3 ribbon is shown in purple. The five hydrophobic residues of Puma BH3 (W133, I137, L141 M144, and L148) are protruding into the binding groove, and the conserved salt bridge interaction formed by BHRF1 R100 and Puma BH3 D146 is labeled, as well as other residues involved in additional ionic interactions and hydrogen bonds. Ionic interactions and hydrogen bonds are shown as dotted black lines, and the participating residues are labeled. The canonical hydrophobic binding pockets h1–h4 are indicated by black circles, while pocket h5 is indicated by a red circle. The molecular orientations depicted in (a–d) are identical to those in Figure 1.
Viruses 14 02222 g002
Viruses 14 02222 g003 550
Figure 3. Interaction profiles and detailed view of the fifth hydrophobic pocket of BHRF1. (a) Summary of BH3 interaction profiles of BHRF1, BHRF1 R100D, and BHRF1 F72W mutant and its closely related Bcl-2 homologs, human Bcl-xL, KSHV KsBcl-2, VARV F1L, and tanapoxvirus TANV16L. (b) Enlarged view of the additional hydrophobic pocket (fifth pocket) in the AlphaFold structure of F72W BHRF1 (c) View of the fifth pocket in BHRF1–BidBH3 with the ligand removed and colored according to residue hydrophobicity; selected residues of BHRF1 are labeled. (d) An identical view of F72W BHRF1 as in (c). The surface, hydrophobic pockets of the binding groove, and key residues are indicated. In (c,d), the surface is colored according to residue hydrophobicity: red: F, V, L, I, M; orange: W, Y, C, A, G, P, S, T; pink: N, E, H; gray: Q, R, K, D. (e) View of the fifth pocket in BHRF1–BidBH3. The surface of the BHRF1 is shown in gray with side-chains of the key hydrophobic residues located in the pocket shown as sticks. The BHRF1 backbone (salmon) and BidBH3 backbone (orange) are shown as tubes. Selected residues of Bid and BHRF1 are labeled. (f) Similar to (e), a view of the fifth hydrophobic pocket of BHRF1–Puma BH3 (purple).
Figure 3. Interaction profiles and detailed view of the fifth hydrophobic pocket of BHRF1. (a) Summary of BH3 interaction profiles of BHRF1, BHRF1 R100D, and BHRF1 F72W mutant and its closely related Bcl-2 homologs, human Bcl-xL, KSHV KsBcl-2, VARV F1L, and tanapoxvirus TANV16L. (b) Enlarged view of the additional hydrophobic pocket (fifth pocket) in the AlphaFold structure of F72W BHRF1 (c) View of the fifth pocket in BHRF1–BidBH3 with the ligand removed and colored according to residue hydrophobicity; selected residues of BHRF1 are labeled. (d) An identical view of F72W BHRF1 as in (c). The surface, hydrophobic pockets of the binding groove, and key residues are indicated. In (c,d), the surface is colored according to residue hydrophobicity: red: F, V, L, I, M; orange: W, Y, C, A, G, P, S, T; pink: N, E, H; gray: Q, R, K, D. (e) View of the fifth pocket in BHRF1–BidBH3. The surface of the BHRF1 is shown in gray with side-chains of the key hydrophobic residues located in the pocket shown as sticks. The BHRF1 backbone (salmon) and BidBH3 backbone (orange) are shown as tubes. Selected residues of Bid and BHRF1 are labeled. (f) Similar to (e), a view of the fifth hydrophobic pocket of BHRF1–Puma BH3 (purple).
Viruses 14 02222 g003
Viruses 14 02222 g004 550
Figure 4. Comparison of BHRF1–Bid BH3 and Puma BH3 interaction sites with other closely related Bcl-2–BH3 complexes. (a) The BHRF1 surface and BH3 interaction sites are shown in gray and red respectively, whilst the Bid BH3 ribbon is shown in orange. (b) Human Bcl-xL–Bid BH3; (c) KSHV KsBcl-2–Bid BH3; (d) VARV F1L–Bid BH3; (e) BHRF1–Puma BH3; (f) human Bcl-xL–Puma BH3; (g) KSHV KsBcl-2–Puma BH3; (h) tanapoxvirus TANV16L–Puma BH3. The interaction formed by Bcl-2–BH3 residues and hydrophobic residues in the ligand-binding groove are labeled. Only residues located in the binding pockets are shown. Polar residues close to the BH3 ligand on the Bcl-2 protein are shown in red. The color schemes for the BH3 peptide and views in (b–h) are as in Figure 1.
Figure 4. Comparison of BHRF1–Bid BH3 and Puma BH3 interaction sites with other closely related Bcl-2–BH3 complexes. (a) The BHRF1 surface and BH3 interaction sites are shown in gray and red respectively, whilst the Bid BH3 ribbon is shown in orange. (b) Human Bcl-xL–Bid BH3; (c) KSHV KsBcl-2–Bid BH3; (d) VARV F1L–Bid BH3; (e) BHRF1–Puma BH3; (f) human Bcl-xL–Puma BH3; (g) KSHV KsBcl-2–Puma BH3; (h) tanapoxvirus TANV16L–Puma BH3. The interaction formed by Bcl-2–BH3 residues and hydrophobic residues in the ligand-binding groove are labeled. Only residues located in the binding pockets are shown. Polar residues close to the BH3 ligand on the Bcl-2 protein are shown in red. The color schemes for the BH3 peptide and views in (b–h) are as in Figure 1.
Viruses 14 02222 g004
Table
Table 1. X-ray crystallographic data collection and refinement statistics. Values in parentheses are for the highest resolution shell.
Table 1. X-ray crystallographic data collection and refinement statistics. Values in parentheses are for the highest resolution shell.
BHRF1–Bid BH3
(PDB ID 7P33)		BHRF1–Puma BH3
(PDB ID 7P9W)
Data collection
Space groupP6522P3221
Cell dimensions
a, b, c (Å)94.20, 94.20, 455.5862.77, 62.77, 92.60
α, β, γ (°)90, 90, 12090, 90,120
Wavelength (Å)0.95370.9537
Resolution (Å)46.7–2.78 (2.88–2.78) 35.25–2.00 (2.07–2.0)
Rsym or Rmerge0.04 (0.38)0.05 (0.51)
I/σI7.86 (1.14)6.76 (1.15)
Completeness (%)98.8 (90.9)99.7 (99.9)
CC1/20.99 (0.66)0.99 (0.57)
Redundancy2.0 (2.0)2.0 (2.0)
Refinement
Resolution (Å)46.7–2.78 (2.88–2.78)35.25–2.00 (2.07–2.0)
No. reflections32,68814,737
Rwork/Rfree0.231/0.2740.217/0.252
Clashscore1.205.18
No. atoms
Protein71921422
Ligand/ion957
Water6254
B-factors
Protein62.6842.03
Ligand/ion104.7851.63
Water51.2949.48
R.m.s. deviations
Bond lengths (Å)0.0020.005
Bond angles (°)0.400.73
Table
Table 2. Percentage sequence identity of BHRF1 with its mammalian counterparts.
Table 2. Percentage sequence identity of BHRF1 with its mammalian counterparts.
Bcl-2Bcl-xLBcl-wMcl-1Bfl-1/A1Bcl-B
BHRF121.7%17.9%18.4%9.6%17.821.5%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Suraweera, C.D.; Hinds, M.G.; Kvansakul, M. Crystal Structures of Epstein–Barr Virus Bcl-2 Homolog BHRF1 Bound to Bid and Puma BH3 Motif Peptides. Viruses 2022, 14, 2222. https://doi.org/10.3390/v14102222

AMA Style

Suraweera CD, Hinds MG, Kvansakul M. Crystal Structures of Epstein–Barr Virus Bcl-2 Homolog BHRF1 Bound to Bid and Puma BH3 Motif Peptides. Viruses. 2022; 14(10):2222. https://doi.org/10.3390/v14102222

Chicago/Turabian Style

Suraweera, Chathura D., Mark G. Hinds, and Marc Kvansakul. 2022. "Crystal Structures of Epstein–Barr Virus Bcl-2 Homolog BHRF1 Bound to Bid and Puma BH3 Motif Peptides" Viruses 14, no. 10: 2222. https://doi.org/10.3390/v14102222

APA Style

Suraweera, C. D., Hinds, M. G., & Kvansakul, M. (2022). Crystal Structures of Epstein–Barr Virus Bcl-2 Homolog BHRF1 Bound to Bid and Puma BH3 Motif Peptides. Viruses, 14(10), 2222. https://doi.org/10.3390/v14102222

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop