| Literature DB >> 35283834 |
Yanlei Zhu1, Flavio Alvarez1, Nicolas Wolff1, Ariel Mechaly2, Sébastien Brûlé3, Benoit Neitthoffer4, Sandrine Etienne-Manneville4, Ahmed Haouz2, Batiste Boëda4, Célia Caillet-Saguy1.
Abstract
The C-terminus of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) protein E contains a PBM (PDZ-binding motif) targeting PDZ (PSD-95/Dlg/ZO-1) domains, which is identical to the PBM of SARS-CoV. The latter is involved in the pathogenicity of the virus. Recently, we identified 10 human PDZ-containing proteins showing significant interactions with SARS-CoV-2 protein E PBM. We selected several of them involved in cellular junctions and cell polarity (TJP1, PARD3, MLLT4, and LNX2) and MPP5/PALS1 previously shown to interact with SARS-CoV E PBM. Targeting cellular junctions and polarity components is a common strategy by viruses to hijack cell machinery to their advantage. In this study, we showed that these host PDZ domains TJP1, PARD3, MLLT4, LNX2, and MPP5/PALS1 interact in a PBM-dependent manner in vitro and colocalize with the full-length E protein in cellulo, sequestrating the PDZ domains to the Golgi compartment. We solved three crystal structures of complexes between human LNX2, MLLT4, and MPP5 PDZs and SARS-CoV-2 E PBM highlighting its binding preferences for several cellular targets. Finally, we showed different affinities for the PDZ domains with the original SARS-CoV-2 C-terminal sequence containing the PBM and the one of the beta variant that contains a mutation close to the PBM. The acquired mutations in the E protein localized near the PBM might have important effects both on the structure and the ion-channel activity of the E protein and on the host machinery targeted by the variants during the infection.Entities:
Keywords: PDZ-binding motif; PDZ-containing protein; SARS-CoV-2; cell junctions and polarity; envelope protein; host–pathogen interactions; protein–protein interaction
Year: 2022 PMID: 35283834 PMCID: PMC8909127 DOI: 10.3389/fmicb.2022.829094
Source DB: PubMed Journal: Front Microbiol ISSN: 1664-302X Impact factor: 5.640
FIGURE 1Schematic domain representation and sequences of the protein E constructs. (A) Domain delimitations and sequence of the full-length protein E of SARS-CoV-2. NTD, TM, and CTD correspond to the N-terminal domain, the transmembrane domain, and the C-terminal domain respectively. (B) Sequence of the full-length protein E of SARS-CoV-2 beta variant (B.1.351). (C) Sequences of the SARS-CoV-2 WT and the beta variant E C-terminal encompassing the PSD-95/Dlg/ZO-1-binding motif (PBM). The PBMs are underlined. The mutation is highlighted in red.
Selected human PSD-95/Dlg/ZO-1 (PDZ) domains targeted by the SARS-CoV-2 E protein PSD-95/Dlg/ZO-1-binding motif (PBM).
| PDZ domain | Uniprot | Delimitations | PDB structure/alphafold model | α 2 residue facing p-2 |
| ZO-1 PDZ2 | Q07157 | 184-264 | 2JWE | Leu |
| LNX2 PDZ2 | Q8N448 | 334-426 | 7QCT | Pro |
| PARD3 PDZ3 | Q8TEW0 | 582-685 | Alphafold prediction ( | Asn |
| MLLT4 | P55196 | 1002-1095 | 7QCR | Gln |
| MPP5 | Q8N3R9 | 238-336 | 7QCS | Val |
PDZ name, Uniprot code, delimitations of the constructs used in this work (residue numbers correspond to the numbering in the full-length human protein), PDB of the structure model used, and the residue of PDZ helix α2 facing the PBM residue at position -2.
FIGURE 2PDZ domains interact with the SARS-CoV-2 full-length E protein in a PBM-dependent manner. (A) Schematic representation of the GFP-E constructs used in this study. Construct 1 corresponds to the GFP alone, construct 2 and 3 to the full-length WT and mutated E protein respectively, and construct 4 to the last 12 C-terminal residues of E protein. TM and PBM correspond to the transmembrane domain and to the PDZ-binding motif, respectively. (B) Input fraction of GFP-tagged viral E gene construct expressed in HEK293 cells and used for pull-down assay. Samples were analyzed by immunoblotting using anti-GFP antibody. (C) Gluthatione S-transferase (GST) pull-down binding results with immobilized GST or GST-PDZ domains (GST-PDZ ZO1, GST-PDZ LNX2, GST-PDZ MLLT4, GST-PDZ MPP5, and GST-PDZ PARD3) used as affinity resin and incubated with the HEK293 cell lysates. GST-tagged proteins were stained by Ponceau S (top panels). The bound fraction was analyzed by immunoblotting using an anti-GFP antibody (bottom panels). Note that blots represent discontinuous panels from the gels when black line delimitation is present. The images have been cropped to frame the relevant region.
FIGURE 3X-ray structures of the PDZ domains of MLLT4, MPP5, and LNX2 bound to the SARS-CoV-2 protein E PBM. (A–C) The asymmetric unit of the PDZ domains of MLLT4, MPP5, and LNX2, respectively, bound to the SARS-CoV-2 protein E PBM shown as red sticks. (A) Selected interchain polar contacts related to the swapped dimer between the fragments Lys 1014–Gly 1017 of each chain are shown in orange, and the associated residues are shown as sticks. (D–F) Detailed views of the PDZ domains bound to SARS-CoV-2 protein E PBM. Important residues are labeled and shown as sticks. Intermolecular H-bonds and polar contacts are reported as orange dashed lines.
X-ray data collection and refinement statistics.
| LNX2-PDZ/SARS-CoV-2-E-PBM | MLLT4-PDZ/SARS-CoV-2-E-PBM | MPP5-PDZ/SARS-CoV-2-E-PBM | |
|
| 2.5 M sodium chloride 0.1 M sodium acetate pH 4.5 0.2 M lithium sulfate | 1 M tri-sodium citrate, 0.1 M TRIS-HCl pH 7.0 0.2 M NaCl | 0.5 M lithium chloride 1.6 M ammonium sulfate |
| Beamline | Proxima-1 | Proxima-1 | Proxima-2A |
| Resolution range | 48.97–3.197 (3.312–3.197) | 37.81–2.281 (2.362–2.281) | 44.38–2.804 (2.904–2.804) |
| Space group | I 4 3 2 | C 1 2 1 | C 2 2 21 |
| Unit cell | 183.23 183.23 183.23 90 90 90 | 76.16 43.49 56.07 90 96.88 90 | 88.75 131.55 70.93 90 90 90 |
| Total reflections | 708,906 (68,100) | 57,314 (4,867) | 44,758 (4,462) |
| Unique reflections | 9,009 (861) | 8,373 (746) | 8,040 (799) |
| Multiplicity | 78.7 (77.6) | 6.8 (6.5) | 5.6 (5.6) |
| Completeness (%) | 99.69 (97.95) | 98.90 (89.23) | 75.85 (78.61) |
| Mean I/sigma(I) | 20.81 (1.39) | 6.59 (1.16) | 9.33 (1.31) |
| Wilson B-factor | 123.76 | 44.58 | 67.85 |
| R-merge | 0.2679 (3.061) | 0.2185 (1.454) | 0.1238 (1.321) |
| R-meas | 0.2697 (3.082) | 0.2367 (1.577) | 0.1365 (1.455) |
| R-pim | 0.03029 (0.3451) | 0.08997 (0.6014) | 0.05584 (0.5944) |
| CC1/2 | 0.999 (0.459) | 0.988 (0.506) | 0.999 (0.689) |
| CC* | 1 (0.793) | 0.997 (0.82) | 1 (0.903) |
| Reflections used in refinement | 8,992 (860) | 8,373 (746) | 8,036 (794) |
| Reflections used for R-free | 451 (44) | 419 (37) | 399 (39) |
| R-work | 0.2281 (0.3597) | 0.2208 (0.3409) | 0.2633 (0.4147) |
| R-free | 0.2395 (0.3810) | 0.2523 (0.4239) | 0.3049 (0.4974) |
| C(work) | 0.908 (0.646) | 0.916 (0.670) | 0.936 (0.687) |
| CC(free) | 0.905 (0.593) | 0.911 (0.399) | 0.891 (0.744) |
| Number of non-hydrogen atoms | 1,483 | 1,410 | 2,204 |
| Macromolecules | 1,460 | 1,331 | 2,204 |
| Ligands | 0 | 5 | 0 |
| Solvent | 23 | 74 | 0 |
| Protein residues | 190 | 183 | 287 |
| RMS(bonds) | 0.012 | 0.012 | 0.013 |
| RMS(angles) | 1.53 | 1.56 | 1.66 |
| Ramachandran favored (%) | 96.15 | 97.66 | 96.39 |
| Ramachandran allowed (%) | 3.30 | 1.75 | 3.25 |
| Ramachandran outliers (%) | 0.55 | 0.58 | 0.36 |
| Rotamer outliers (%) | 11.69 | 7.64 | 12.24 |
| Clashscore | 10.15 | 6.51 | 5.61 |
| Average B-factor | 130.45 | 51.33 | 89.21 |
| Macromolecules | 130.97 | 51.32 | 89.21 |
| Ligands | 128.47 | ||
| Solvent | 97.36 | 46.31 | |
| PDB entry ID | 7QCT | 7QCR | 7QCS |
Statistics for the highest-resolution shell are shown in parentheses.
FIGURE 4The viral E protein recruits the PDZ domain of ZO1 to the Golgi apparatus. (A) Hela cell cotransfected with encoding constructs GFP-PDZ2 of ZO1 and ALFA-E (red) displays Golgi recruitment of both proteins as indicated by GM130 staining (Magenta). (B) Cotransfection of GFP-ZO1 with ALFA-E PBM mutant (red) does not lead to the recruitment of GFP-ZO1 to the Golgi apparatus. (C) Hela cell cotransfected with GFP-PDZ1 of SCRIB and ALFA-E (red) does not show recruitment of the proteins to the Golgi apparatus (Magenta). White arrowheads indicate Golgi apparatus position. Nuclei are stained with DAPI. Bars correspond to 10 μm.
Kd values between the E protein PBM of SARS-CoV-2 WT and of the P71L mutant for the selected PDZ domains.
| Kd (μ M) | ||
|
| ||
| SARS-CoV-2 WT | SARS-CoV-2 P71L | |
| ZO1 PDZ2 | 15 (±11) | 133 (±28) |
| MPP5 PDZ | 30 (±24) | 22 (±2) |
| LNX2 PDZ2 | 289 (±116) | 47 (±20) |
| PARD3 PDZ3 | 341 (±57) | No interaction detected |
| MLLT4 PDZ | 569 (±129) | No interaction detected |
The data are representative of two independent experiments and error bars correspond to the standard deviation.
Hydrodynamic parameters of ZO-1 PDZ2 and MLLT4 PDZ derived from the analysis of analytical ultracentrifugation.
| ZO-1 1 mg/mL | MLLT4-PBM 1 mg/mL | |||||
| Sedimentation coefficient (S) | 2.1 | 2.8 | 3.7 | 0.5 | 1.4 | 2.3 |
| Frictional ratio f/f0 (average) | 1.4 | 1.2 | ||||
| Molecular weight estimated by analytical ultracentrifugation (AUC) (kDa) | 26 | 40 | 60 | 2.8 | 12 | 26 |
| Peak proportion% | 49 | 49 | 2 | 25 | 65 | 10 |