Hidden kinetic traps in multidomain folding highlight the presence of a misfolded but functionally competent intermediate.

Although more than 75% of the proteome is composed of multidomain proteins, current knowledge of protein folding is based primarily on studies of isolated domains. In this work, we describe the folding mechanism of a multidomain tandem construct comprising two distinct covalently bound PDZ domains belonging to a protein called Whirlin, a scaffolding protein of the hearing apparatus. In particular, via a synergy between NMR and kinetic experiments, we demonstrate the presence of a misfolded intermediate that competes with productive folding. In agreement with the view that tandem domain swapping is a potential source of transient misfolding, we demonstrate that such a kinetic trap retains native-like functional activity, as shown by the preserved ability to bind its physiological ligand. Thus, despite the general knowledge that protein misfolding is intimately associated with dysfunction and diseases, we provide a direct example of a functionally competent misfolded state. Remarkably, a bioinformatics analysis of the amino acidic sequence of Whirlin from different species suggests that the tendency to perform tandem domain swapping between PDZ1 and PDZ2 is highly conserved, as demonstrated by their unexpectedly high sequence identity. On the basis of these observations, we discuss on a possible physiological role of such misfolded intermediate.

Ultrarare heterozygous pathogenic variants of genes causing dominant forms of early-onset deafness underlie severe presbycusis.

Presbycusis, or age-related hearing loss (ARHL), is a major public health issue. About half the phenotypic variance has been attributed to genetic factors. Here, we assessed the contribution to presbycusis of ultrarare pathogenic variants, considered indicative of Mendelian forms. We focused on severe presbycusis without environmental or comorbidity risk factors and studied multiplex family age-related hearing loss (mARHL) and simplex/sporadic age-related hearing loss (sARHL) cases and controls with normal hearing by whole-exome sequencing. Ultrarare variants (allele frequency [AF] < 0.0001) of 35 genes responsible for autosomal dominant early-onset forms of deafness, predicted to be pathogenic, were detected in 25.7% of mARHL and 22.7% of sARHL cases vs. 7.5% of controls (P = 0.001); half were previously unknown (AF < 0.000002). MYO6, MYO7A, PTPRQ, and TECTA variants were present in 8.9% of ARHL cases but less than 1% of controls. Evidence for a causal role of variants in presbycusis was provided by pathogenicity prediction programs, documented haploinsufficiency, three-dimensional structure/function analyses, cell biology experiments, and reported early effects. We also established Tmc1N321I/+ mice, carrying the TMC1:p.(Asn327Ile) variant detected in an mARHL case, as a mouse model for a monogenic form of presbycusis. Deafness gene variants can thus result in a continuum of auditory phenotypes. Our findings demonstrate that the genetics of presbycusis is shaped by not only well-studied polygenic risk factors of small effect size revealed by common variants but also, ultrarare variants likely resulting in monogenic forms, thereby paving the way for treatment with emerging inner ear gene therapy.

Phylogenetic analysis of Harmonin homology domains.

BACKGROUND: Harmonin Homogy Domains (HHD) are recently identified orphan domains of about 70 residues folded in a compact five alpha-helix bundle that proved to be versatile in terms of function, allowing for direct binding to a partner as well as regulating the affinity and specificity of adjacent domains for their own targets. Adding their small size and rather simple fold, HHDs appear as convenient modules to regulate protein-protein interactions in various biological contexts. Surprisingly, only nine HHDs have been detected in six proteins, mainly expressed in sensory neurons. RESULTS: Here, we built a profile Hidden Markov Model to screen the entire UniProtKB for new HHD-containing proteins. Every hit was manually annotated, using a clustering approach, confirming that only a few proteins contain HHDs. We report the phylogenetic coverage of each protein and build a phylogenetic tree to trace the evolution of HHDs. We suggest that a HHD ancestor is shared with Paired Amphipathic Helices (PAH) domains, a four-helix bundle partially sharing fold and functional properties. We characterized amino-acid sequences of the various HHDs using pairwise BLASTP scoring coupled with community clustering and manually assessed sequence features among each individual family. These sequence features were analyzed using reported structures as well as homology models to highlight structural motifs underlying HHDs fold. We show that functional divergence is carried out by subtle differences in sequences that automatized approaches failed to detect. CONCLUSIONS: We provide the first HHD databases, including sequences and conservation, phylogenic trees and a list of HHD variants found in the auditory system, which are available for the community. This case study highlights surprising phylogenetic properties found in orphan domains and will assist further studies of HHDs. We unveil the implication of HHDs in their various binding interfaces using conservation across families and a new protein-protein surface predictor. Finally, we discussed the functional consequences of three identified pathogenic HHD variants involved in Hoyeraal-Hreidarsson syndrome and of three newly reported pathogenic variants identified in patients suffering from Usher Syndrome.

Structure-based optimization of a PDZ-binding motif within a viral peptide stimulates neurite outgrowth.

Protection of neuronal homeostasis is a major goal in the management of neurodegenerative diseases. Microtubule-associated Ser/Thr kinase 2 (MAST2) inhibits neurite outgrowth, and its inhibition therefore represents a potential therapeutic strategy. We previously reported that a viral protein (G-protein from rabies virus) capable of interfering with protein-protein interactions between the PDZ domain of MAST2 and the C-terminal moieties of its cellular partners counteracts MAST2-mediated suppression of neurite outgrowth. Here, we designed peptides derived from the native viral protein to increase the affinity of these peptides for the MAST2-PDZ domain. Our strategy involved modifying the length and flexibility of the noninteracting sequence linking the two subsites anchoring the peptide to the PDZ domain. Three peptides, Neurovita1 (NV1), NV2, and NV3, were selected, and we found that they all had increased affinities for the MAST2-PDZ domain, with Kd values decreasing from 1300 to 60 nm, while target selectivity was maintained. A parallel biological assay evaluating neurite extension and branching in cell cultures revealed that the NV peptides gradually improved neural activity, with the efficacies of these peptides for stimulating neurite outgrowth mirroring their affinities for MAST2-PDZ. We also show that NVs can be delivered into the cytoplasm of neurons as a gene or peptide. In summary, our findings indicate that virus-derived peptides targeted to MAST2-PDZ stimulate neurite outgrowth in several neuron types, opening up promising avenues for potentially using NVs in the management of neurodegenerative diseases.

Quantifying domain-ligand affinities and specificities by high-throughput holdup assay.

Many protein interactions are mediated by small linear motifs interacting specifically with defined families of globular domains. Quantifying the specificity of a motif requires measuring and comparing its binding affinities to all its putative target domains. To this end, we developed the high-throughput holdup assay, a chromatographic approach that can measure up to 1,000 domain-motif equilibrium binding affinities per day. After benchmarking the approach on 210 PDZ-peptide pairs with known affinities, we determined the affinities of two viral PDZ-binding motifs derived from human papillomavirus E6 oncoproteins for 209 PDZ domains covering 79% of the human 'PDZome'. We obtained sharply sequence-dependent binding profiles that quantitatively describe the PDZome recognition specificity of each motif. This approach, applicable to many categories of domain-ligand interactions, has wide potential for quantifying the specificities of interactomes.

Molecular Basis of the Interaction of the Human Protein Tyrosine Phosphatase Non-receptor Type 4 (PTPN4) with the Mitogen-activated Protein Kinase p38γ.

The human protein tyrosine phosphatase non-receptor type 4 (PTPN4) prevents cell death induction in neuroblastoma and glioblastoma cell lines in a PDZ·PDZ binding motifs-dependent manner, but the cellular partners of PTPN4 involved in cell protection are unknown. Here, we described the mitogen-activated protein kinase p38γ as a cellular partner of PTPN4. The main contribution to the p38γ·PTPN4 complex formation is the tight interaction between the C terminus of p38γ and the PDZ domain of PTPN4. We solved the crystal structure of the PDZ domain of PTPN4 bound to the p38γ C terminus. We identified the molecular basis of recognition of the C-terminal sequence of p38γ that displays the highest affinity among all endogenous partners of PTPN4. We showed that the p38γ C terminus is also an efficient inducer of cell death after its intracellular delivery. In addition to recruiting the kinase, the binding of the C-terminal sequence of p38γ to PTPN4 abolishes the catalytic autoinhibition of PTPN4 and thus activates the phosphatase, which can efficiently dephosphorylate the activation loop of p38γ. We presume that the p38γ·PTPN4 interaction promotes cellular signaling, preventing cell death induction.

Quantitative and dynamic analysis of PTEN phosphorylation by NMR.

The dual lipid and protein phosphatase PTEN is a tumor suppressor controlling key biological processes, such as cell growth, proliferation and neuro-survival. Its activity and intracellular trafficking is finely regulated notably by multi-site phosphorylation of its C-terminal tail. The reversible and highly dynamic character of these regulatory events confers a temporal dimension to the cell for triggering crucial decisions. In this review, we describe how a recently developed time-resolved NMR spectroscopy approach unveils the dynamic establishment of the phosphorylation events of PTEN C-terminal tail controlled by CK2 and GSK3ß kinases. Two cascades of reactions have been identified, in vitro and in extracts of human neuroblastoma cells. They are triggered independently on two nearby clusters of sites (S380-S385 and S361-S370) and occur on different timescales. In each cascade, the reactions follow an ordered model with a distributive kinetic mechanism. The vision of these cascades as two delay timers activating distinct or time-delayed regulatory responses gives a temporal dimension on PTEN regulation and is discussed in relation to the known functional roles of each cluster.

Deciphering the unconventional peptide binding to the PDZ domain of MAST2.

Phosphatase and tensin homologue (PTEN) and microtubule-associated serine threonine kinase 2 (MAST2) are key negative regulators of survival pathways in neuronal cells. The two proteins interact via the PDZ (PSD-95, Dlg1, Zo-1) domain of MAST2 (MAST2-PDZ). During infection by rabies virus, the viral glycoprotein competes with PTEN for interaction with MAST2-PDZ and promotes neuronal survival. The C-terminal PDZ-binding motifs (PBMs) of the two proteins bind similarly to MAST2-PDZ through an unconventional network of connectivity involving two anchor points. Combining stopped-flow fluorescence, analytical ultracentrifugation (AUC), microcalorimetry and NMR, we document the kinetics of interaction between endogenous and viral ligands to MAST2-PDZ as well as the dynamic and structural effects of these interactions. Viral and PTEN peptide interactions to MAST2-PDZ occur via a unique kinetic step which involves both canonical C-terminal PBM binding and N-terminal anchoring. Indirect effects induced by the PBM binding include modifications to the structure and dynamics of the PDZ dimerization surface which prevent MAST2-PDZ auto-association. Such an energetic communication between binding sites and distal surfaces in PDZ domains provides interesting clues for protein regulation overall.

Retinal Vascular Changes in Heart Failure with Preserved Ejection Fraction Using Optical Coherence Tomography Angiography.

Background: Systemic microvascular regression and dysfunction are considered important underlying mechanisms in heart failure with preserved ejection fraction (HFpEF), but retinal changes are unknown. Methods: This prospective study aimed to investigate whether retinal microvascular and structural parameters assessed using optical coherence tomography angiography (OCT-A) differ between patients with HFpEF and control individuals (i.e., capillary vessel density, thickness of retina layers). We also aimed to assess the associations of retinal parameters with clinical and echocardiographic parameters in HFpEF. HFpEF patients, but not controls, underwent echocardiography. Macula-centered 6 × 6 mm volume scans were computed of both eyes. Results: Twenty-two HFpEF patients and 24 controls without known HFpEF were evaluated, with an age of 74 [68-80] vs. 68 [58-77] years (p = 0.027), and 73% vs. 42% females (p = 0.034), respectively. HFpEF patients showed vascular degeneration compared to controls, depicted by lower macular vessel density (p < 0.001) and macular ganglion cell-inner plexiform layer thickness (p = 0.025), and a trend towards lower total retinal volume (p = 0.050) on OCT-A. In HFpEF, a lower total retinal volume was associated with markers of diastolic dysfunction (septal e', septal and average E/e': R2 = 0.38, 0.36, 0.25, respectively; all p < 0.05), even after adjustment for age, sex, diabetes mellitus, or atrial fibrillation. Conclusions: Patients with HFpEF showed clear levels of retinal vascular changes compared to control individuals, and retinal alterations appeared to be associated with markers of more severe diastolic dysfunction in HFpEF. OCT-A may therefore be a promising technique for monitoring systemic microvascular regression and cardiac diastolic dysfunction.

Role of the iron axial ligands of heme carrier HasA in heme uptake and release.

The hemophore protein HasA from Serratia marcescens cycles between two states as follows: the heme-bound holoprotein, which functions as a carrier of the metal cofactor toward the membrane receptor HasR, and the heme-free apoprotein fishing for new porphyrin to be taken up after the heme has been delivered to HasR. Holo- and apo-forms differ for the conformation of the two loops L1 and L2, which provide the axial ligands of the iron through His(32) and Tyr(75), respectively. In the apo-form, loop L1 protrudes toward the solvent far away from loop L2; in the holoprotein, closing of the loops on the heme occurs upon establishment of the two axial coordination bonds. We have established that the two variants obtained via single point mutations of either axial ligand (namely H32A and Y75A) are both in the closed conformation. The presence of the heme and one out of two axial ligands is sufficient to establish a link between L1 and L2, thanks to the presence of coordinating solvent molecules. The latter are stabilized in the iron coordination environment by H-bond interactions with surrounding protein residues. The presence of such a water molecule in both variants is revealed here through a set of different spectroscopic techniques. Previous studies had shown that heme release and uptake processes occur via intermediate states characterized by a Tyr(75)-iron-bound form with open conformation of loop L1. Here, we demonstrate that these states do not naturally occur in the free protein but can only be driven by the interaction with the partner proteins.

Interactions of the protein tyrosine phosphatase PTPN3 with viral and cellular partners through its PDZ domain: insights into structural determinants and phosphatase activity.

The human protein tyrosine phosphatase non-receptor type 3 (PTPN3) is a phosphatase containing a PDZ (PSD-95/Dlg/ZO-1) domain that has been found to play both tumor-suppressive and tumor-promoting roles in various cancers, despite limited knowledge of its cellular partners and signaling functions. Notably, the high-risk genital human papillomavirus (HPV) types 16 and 18 and the hepatitis B virus (HBV) target the PDZ domain of PTPN3 through PDZ-binding motifs (PBMs) in their E6 and HBc proteins respectively. This study focuses on the interactions between the PTPN3 PDZ domain (PTPN3-PDZ) and PBMs of viral and cellular protein partners. We solved the X-ray structures of complexes between PTPN3-PDZ and PBMs of E6 of HPV18 and the tumor necrosis factor-alpha converting enzyme (TACE). We provide new insights into key structural determinants of PBM recognition by PTPN3 by screening the selectivity of PTPN3-PDZ recognition of PBMs, and by comparing the PDZome binding profiles of PTPN3-recognized PBMs and the interactome of PTPN3-PDZ. The PDZ domain of PTPN3 was known to auto-inhibit the protein's phosphatase activity. We discovered that the linker connecting the PDZ and phosphatase domains is involved in this inhibition, and that the binding of PBMs does not impact this catalytic regulation. Overall, the study sheds light on the interactions and structural determinants of PTPN3 with its cellular and viral partners, as well as on the inhibitory role of its PDZ domain on its phosphatase activity.

New mutation in the ß1 propeller domain of LRP4 responsible for congenital myasthenic syndrome associated with Cenani-Lenz syndrome.

Congenital myasthenic syndromes (CMS) are a clinically and genetically heterogeneous group of rare diseases due to mutations in neuromuscular junction (NMJ) protein-coding genes. Until now, many mutations encoding postsynaptic proteins as Agrin, MuSK and LRP4 have been identified as responsible for increasingly complex CMS phenotypes. The majority of mutations identified in LRP4 gene causes bone diseases including CLS and sclerosteosis-2 and rare cases of CMS with mutations in LRP4 gene has been described so far. In the French cohort of CMS patients, we identified a novel LRP4 homozygous missense mutation (c.1820A > G; p.Thy607Cys) within the ß1 propeller domain in a patient presenting CMS symptoms, including muscle weakness, fluctuating fatigability and a decrement in compound muscle action potential in spinal accessory nerves, associated with congenital agenesis of the hands and feet and renal malformation. Mechanistic expression studies show a significant decrease of AChR aggregation in cultured patient myotubes, as well as altered in vitro binding of agrin and Wnt11 ligands to the mutated ß1 propeller domain of LRP4 explaining the dual phenotype characterized clinically and electoneuromyographically in the patient. These results expand the LRP4 mutations spectrum associated with a previously undescribed clinical association involving impaired neuromuscular transmission and limb deformities and highlighting the critical role of a yet poorly described domain of LRP4 at the NMJ. This study raises the question of the frequency of this rare neuromuscular form and the future diagnosis and management of these cases.

Neuroinvasion and anosmia are independent phenomena upon infection with SARS-CoV-2 and its variants.

Anosmia was identified as a hallmark of COVID-19 early in the pandemic, however, with the emergence of variants of concern, the clinical profile induced by SARS-CoV-2 infection has changed, with anosmia being less frequent. Here, we assessed the clinical, olfactory and neuroinflammatory conditions of golden hamsters infected with the original Wuhan SARS-CoV-2 strain, its isogenic ORF7-deletion mutant and three variants: Gamma, Delta, and Omicron/BA.1. We show that infected animals develop a variant-dependent clinical disease including anosmia, and that the ORF7 of SARS-CoV-2 contributes to the induction of olfactory dysfunction. Conversely, all SARS-CoV-2 variants are neuroinvasive, regardless of the clinical presentation they induce. Taken together, this confirms that neuroinvasion and anosmia are independent phenomena upon SARS-CoV-2 infection. Using newly generated nanoluciferase-expressing SARS-CoV-2, we validate the olfactory pathway as a major entry point into the brain in vivo and demonstrate in vitro that SARS-CoV-2 travels retrogradely and anterogradely along axons in microfluidic neuron-epithelial networks.

Ordered phosphorylation events in two independent cascades of the PTEN C-tail revealed by NMR.

PTEN phosphatase is a tumor suppressor controlling notably cell growth, proliferation and survival. The multisite phosphorylation of the PTEN C-terminal tail regulates PTEN activity and intracellular trafficking. The dynamical nature of such regulatory events represents a crucial dimension for timing cellular decisions. Here we show that NMR spectroscopy allows reporting on the order and kinetics of clustered multisite phosphorylation events. We first unambiguously identify in vitro seven bona fide sites modified by CK2 and GSK3ß kinases and two new sites on the PTEN C-terminal tail. Then, monitoring the formation of transient intermediate phosphorylated states, we determine the sequence of these reactions and calculate their apparent rate constants. Finally, we assess the dynamic formation of these phosphorylation events induced by endogenous kinases directly in extracts of human neuroblastoma cells. Taken together, our data indicate that two cascades of events controlled by CK2 and GSK3ß occur independently on two clusters of sites (S380-S385 and S361-S370) and that in each cluster the reactions follow an ordered model with a distributive kinetic mechanism. Besides emphasizing the ability of NMR to quantitatively and dynamically follow post-translational modifications, these results bring a temporal dimension on the establishment of PTEN phosphorylation cascades.

Deciphering the Molecular Interaction Between the Adhesion G Protein-Coupled Receptor ADGRV1 and its PDZ-Containing Regulator PDZD7.

Hearing relies on the transduction of sound-evoked vibrations into electrical signals, occurring in the stereocilia bundle of inner ear hair cells. The G protein-coupled receptor (GPCR) ADGRV1 and the multi-PDZ protein PDZD7 play a critical role in the formation and function of stereocilia through their scaffolding and signaling properties. During hair cell development, the GPCR activity of ADGRV1 is specifically inhibited by PDZD7 through an unknown mechanism. Here, we describe the key interactions mediated by the two N-terminal PDZ domains of PDZD7 and the cytoplasmic domain of ADGRV1. Both PDZ domains can bind to the C-terminal PDZ binding motif (PBM) of ADGRV1 with the critical contribution of atypical C-terminal ß extensions. The two PDZ domains form a supramodule in solution, stabilized upon PBM binding. Interestingly, we showed that the stability and binding properties of the PDZ tandem are affected by two deafness-causing mutations located in the binding grooves of PDZD7 PDZ domains.

Interactions of Severe Acute Respiratory Syndrome Coronavirus 2 Protein E With Cell Junctions and Polarity PSD-95/Dlg/ZO-1-Containing Proteins.

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.

PDZ-Containing Proteins Targeted by the ACE2 Receptor.

Angiotensin-converting enzyme 2 (ACE2) is a main receptor for SARS-CoV-2 entry to the host cell. Indeed, the first step in viral entry is the binding of the viral trimeric spike (S) protein to ACE2. Abundantly present in human epithelial cells of many organs, ACE2 is also expressed in the human brain. ACE2 is a type I membrane protein with an extracellular N-terminal peptidase domain and a C-terminal collectrin-like domain that ends with a single transmembrane helix and an intracellular 44-residue segment. This C-terminal segment contains a PDZ-binding motif (PBM) targeting protein-interacting domains called PSD-95/Dlg/ZO-1 (PDZ). Here, we identified the human PDZ specificity profile of the ACE2 PBM using the high-throughput holdup assay and measuring the binding intensities of the PBM of ACE2 against the full human PDZome. We discovered 14 human PDZ binders of ACE2 showing significant binding with dissociation constants' values ranging from 3 to 81 µM. NHERF, SHANK, and SNX27 proteins found in this study are involved in protein trafficking. The PDZ/PBM interactions with ACE2 could play a role in ACE2 internalization and recycling that could be of benefit for the virus entry. Interestingly, most of the ACE2 partners we identified are expressed in neuronal cells, such as SHANK and MAST families, and modifications of the interactions between ACE2 and these neuronal proteins may be involved in the neurological symptoms of COVID-19.

Automatic Bayesian Weighting for SAXS Data.

Small-angle X-ray scattering (SAXS) experiments are important in structural biology because they are solution methods, and do not require crystallization of protein complexes. Structure determination from SAXS data, however, poses some difficulties. Computation of a SAXS profile from a protein model is expensive in CPU time. Hence, rather than directly refining against the data, most computational methods generate a large number of conformers and then filter the structures based on how well they satisfy the SAXS data. To address this issue in an efficient manner, we propose here a Bayesian model for SAXS data and use it to directly drive a Monte Carlo simulation. We show that the automatic weighting of SAXS data is the key to finding optimal structures efficiently. Another key problem with obtaining structures from SAXS data is that proteins are often flexible and the data represents an average over a structural ensemble. To address this issue, we first characterize the stability of the best model with extensive molecular dynamics simulations. We analyse the resulting trajectories further to characterize a dynamic structural ensemble satisfying the SAXS data. The combination of methods is applied to a tandem of domains from the protein PTPN4, which are connected by an unstructured linker. We show that the SAXS data contain information that supports and extends other experimental findings. We also show that the conformation obtained by the Bayesian analysis is stable, but that a minor conformation is present. We propose a mechanism in which the linker may maintain PTPN4 in an inhibited enzymatic state.

PDZ Sample Quality Assessment by Biochemical and Biophysical Characterizations.

PDZ domains are small globular domains involved in protein-protein interactions. They participate in a wide range of critical cellular processes. These domains, very abundant in the human proteome, are widely studied by high-throughput interactomics approaches and by biophysical and structural methods. However, the quality of the results is strongly related to the optimal folding and solubility of the domains. We provide here a detailed description of protocols for a strict quality assessment of the PDZ constructs. We describe appropriate experimental approaches that have been selected to overcome the small size of such domains to check the purity, identity, homogeneity, stability, and folding of samples.

Host PDZ-containing proteins targeted by SARS-CoV-2.

Small linear motifs targeting protein interacting domains called PSD-95/Dlg/ZO-1 (PDZ) have been identified at the C terminus of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins E, 3a, and N. Using a high-throughput approach of affinity-profiling against the full human PDZome, we identified sixteen human PDZ binders of SARS-CoV-2 proteins E, 3A, and N showing significant interactions with dissociation constants values ranging from 3 to 82 µm. Six of them (TJP1, PTPN13, HTRA1, PARD3, MLLT4, LNX2) are also recognized by SARS-CoV while three (NHERF1, MAST2, RADIL) are specific to SARS-CoV-2 E protein. Most of these SARS-CoV-2 protein partners are involved in cellular junctions/polarity and could be also linked to evasion mechanisms of the immune responses during viral infection. Among the binders of the SARS-CoV-2 proteins E, 3a, or N, seven significantly affect viral replication under knock down gene expression in infected cells. This PDZ profiling identifying human proteins potentially targeted by SARS-CoV-2 can help to understand the multifactorial severity of COVID19 and to conceive effective anti-coronaviral agents for therapeutic purposes.