PDZome-wide and structural characterization of the PDZ-binding motif of VANGL2.

VANGL2 is a core component of the non-canonical Wnt/Planar Cell Polarity signaling pathway that uses its highly conserved carboxy-terminal type 1 PDZ-binding motif (PBM) to bind a variety of PDZ proteins. In this study, we characterize and quantitatively assess the largest VANGL2 PDZome-binding profile documented so far, using orthogonal methods. The results of our holdup approach support VANGL2 interactions with a large panel of both long-recognized and unprecedented PDZ domains. Truncation and point mutation analyses of the VANGL2 PBM establish that, beyond the strict requirement of the P-0 / V521 and P-2 / T519 amino acids, upstream residues, including E518, Q516 and R514 at, respectively, P-3, P-5 and P-7 further contribute to the robustness of VANGL2 interactions with two distinct PDZ domains, SNX27 and SCRIBBLE-PDZ3. In agreement with these data, incremental amino-terminal deletions of the VANGL2 PBM causes its overall affinity to progressively decline. Moreover, the holdup data establish that the PDZome binding repertoire of VANGL2 starts to diverge significantly with the truncation of E518. A structural analysis of the SYNJ2BP-PDZ/VANGL2 interaction with truncated PBMs identifies a major conformational change in the binding direction of the PBM peptide after the P-2 position. Finally, we report that the PDZome binding profile of VANGL2 is dramatically rearranged upon phosphorylation of S517, T519 and S520. Our crystallographic approach illustrates how SYNJ2BP accommodates a S520-phosphorylated PBM peptide through the ideal positioning of two basic residues, K48 and R86. Altogether our data provides a comprehensive view of the VANGL2 PDZ network and how this network specifically responds to the post-translation modification of distinct PBM residues. These findings should prove useful in guiding future functional and molecular studies of the key PCP component VANGL2.

Lipid exchange in crystal-confined fatty acid binding proteins: X-ray evidence and molecular dynamics explanation.

Fatty acid binding proteins (FABPs) are responsible for the long-chain fatty acids (FAs) transport inside the cell. However, despite the years, since their structure is known and the many studies published, there is no definitive answer about the stages of the lipid entry-exit mechanism. Their structure forms a ß -barrel of 10 anti-parallel strands with a cap in a helix-turn-helix motif, and there is some consensus on the role of the so-called portal region, involving the second α -helix from the cap ( α 2), ß C- ß D, and ß E- ß F turns in FAs exchange. To test the idea of a lid that opens, we performed a soaking experiment on an h-FABP crystal in which the cap is part of the packing contacts, and its movement is strongly restricted. Even in these conditions, we observed the replacement of palmitic acid by 2-Bromohexadecanoic acid (Br-palmitic acid). Our MD simulations reveal a two-step lipid entry process: (i) The travel of the lipid head through the cavity in the order of tens of nanoseconds, and (ii) The accommodation of its hydrophobic tail in hundreds to thousands of nanoseconds. We observed this even in the cases in which the FAs enter the cavity by their tail. During this process, the FAs do not follow a single trajectory, but multiple ones through which they get into the protein cavity. Thanks to the complementary views between experiment and simulation, we can give an approach to a mechanistic view of the exchange process.

Evidence for direct interaction between the oncogenic proteins E6 and E7 of high-risk human papillomavirus (HPV).

Human papillomaviruses (HPVs) are DNA tumor viruses that infect mucosal and cutaneous epithelial cells of more than 20 vertebrates. High-risk HPV causes about 5% of human cancers worldwide, and the viral proteins E6 and E7 promote carcinogenesis by interacting with tumor suppressors and interfering with many cellular pathways. As a consequence, they immortalize cells more efficiently in concert than individually. So far, the networks of E6 and E7 with their respective cellular targets have been studied extensively but independently. However, we hypothesized that E6 and E7 might also interact directly with each other in a novel interaction affecting HPV-related carcinogenesis. Here, we report a direct interaction between E6 and E7 proteins from carcinogenic HPV types 16 and 31. We demonstrated this interaction via cellular assays using two orthogonal methods: coimmunoprecipitation and flow cytometry-based FRET assays. Analytical ultracentrifugation of the recombinant proteins revealed that the stoichiometry of the E6/E7 complex involves two E7 molecules and two E6 molecules. In addition, fluorescence polarization showed that (I) E6 binds to E7 with a similar affinity for HPV16 and HPV31 (in the same micromolar range) and (II) that the binding interface involves the unstructured N-terminal region of E7. The direct interaction of these highly conserved papillomaviral oncoproteins may provide a new perspective for studying HPV-associated carcinogenesis and the overall viral life cycle.

Quantitative fragmentomics allow affinity mapping of interactomes.

Human protein networks have been widely explored but most binding affinities remain unknown, hindering quantitative interactome-function studies. Yet interactomes rely on minimal interacting fragments displaying quantifiable affinities. Here, we measure the affinities of 65,000 interactions involving PDZ domains and their target PDZ-binding motifs (PBM) within a human interactome region particularly relevant for viral infection and cancer. We calculate interactomic distances, identify hot spots for viral interference, generate binding profiles and specificity logos, and explain selected cases by crystallographic studies. Mass spectrometry experiments on cell extracts and literature surveys show that quantitative fragmentomics effectively complements protein interactomics by providing affinities and completeness of coverage, putting a full human interactome affinity survey within reach. Finally, we show that interactome hijacking by the viral PBM of human papillomavirus E6 oncoprotein substantially impacts the host cell proteome beyond immediate E6 binders, illustrating the complex system-wide relationship between interactome and function.

A scalable strategy to solve structures of PDZ domains and their complexes.

The human PDZome represents one of the largest globular domain families in the human proteome, with 266 instances. These globular domains typically interact with C-terminal peptide motifs found in thousands of human proteins. Despite previous efforts, not all PDZ domains have experimentally solved structures and most of their complexes remain to be solved. Here, a simple and cost-effective strategy is proposed for the crystallization of PDZ domains and their complexes. A human annexin A2 fusion tag was used as a crystallization chaperone and the structures of nine PDZ domains were solved, including five domains that had not yet been solved. Finally, these novel experimental structures were compared with AlphaFold predictions and it is speculated how predictions and experimental methods could cooperate in order to investigate the structural landscapes of entire domain families and interactomes.

High-Risk Mucosal Human Papillomavirus 16 (HPV16) E6 Protein and Cutaneous HPV5 and HPV8 E6 Proteins Employ Distinct Strategies To Interfere with Interferon Regulatory Factor 3-Mediated Beta Interferon Expression.

Persistent infection with some mucosal α-genus human papillomaviruses (HPVs; the most prevalent one being HPV16) can induce cervical carcinoma, anogenital cancers, and a subset of head and neck squamous cell carcinoma (HNSCC). Cutaneous ß-genus HPVs (such as HPV5 and HPV8) associate with skin lesions that can progress into squamous cell carcinoma with sun exposure in Epidermodysplasia verruciformis patients and immunosuppressed patients. Here, we analyzed mechanisms used by E6 proteins from the α- and ß-genus to inhibit the interferon-ß (IFNB1) response. HPV16 E6 mediates this effect by a strong direct interaction with interferon regulatory factor 3 (IRF3). The binding site of E6 was localized within a flexible linker between the DNA-binding domain and the IRF-activation domain of IRF3 containing an LxxLL motif. The crystallographic structure of the complex between HPV16 E6 and the LxxLL motif of IRF3 was solved and compared with the structure of HPV16 E6 interacting with the LxxLL motif of the ubiquitin ligase E6AP. In contrast, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3-binding domain (IBiD) of the CREB-binding protein (CBP), a key transcriptional coactivator in IRF3-mediated IFN-ß expression. IMPORTANCE Persistent HPV infections can be associated with the development of several cancers. The ability to persist depends on the ability of the virus to escape the host immune system. The type I interferon (IFN) system is the first-line antiviral defense strategy. HPVs carry early proteins that can block the activation of IFN-I. Among mucosal α-genus HPV types, the HPV16 E6 protein has a remarkable property to strongly interact with the transcription factor IRF3. Instead, cutaneous HPV5 and HPV8 E6 proteins bind to the IRF3 cofactor CBP. These results highlight the versatility of E6 proteins to interact with different cellular targets. The interaction between the HPV16 E6 protein and IRF3 might contribute to the higher prevalence of HPV16 than that of other high-risk mucosal HPV types in HPV-associated cancers.

Structure of High-Risk Papillomavirus 31 E6 Oncogenic Protein and Characterization of E6/E6AP/p53 Complex Formation.

The degradation of p53 is a hallmark of high-risk human papillomaviruses (HPVs) of the alpha genus and HPV-related carcinogenicity. The oncoprotein E6 forms a ternary complex with the E3 ubiquitin ligase E6-associated protein (E6AP) and tumor suppressor protein p53 targeting p53 for ubiquitination. The extent of p53 degradation by different E6 proteins varies greatly, even for the closely related HPV16 and HPV31. HPV16 E6 and HPV31 E6 display high sequence identity (∼67%). We report here, for the first time, the structure of HPV31 E6 bound to the LxxLL motif of E6AP. HPV16 E6 and HPV31 E6 are structurally very similar, in agreement with the high sequence conservation. Both E6 proteins bind E6AP and degrade p53. However, the binding affinities of 31 E6 to the LxxLL motif of E6AP and p53, respectively, are reduced 2-fold and 5.4-fold compared to 16 E6. The affinity of E6-E6AP-p53 ternary complex formation parallels the efficacy of the subsequent reaction, namely, degradation of p53. Therefore, closely related E6 proteins addressing the same cellular targets may still diverge in their binding efficiencies, possibly explaining their different phenotypic or pathological impacts.IMPORTANCE Variations of carcinogenicity of human papillomaviruses are related to variations of the E6 and E7 interactome. While different HPV species and genera are known to target distinct host proteins, the fine differences between E6 and E7 of closely related HPVs, supposed to target the same cellular protein pools, remain to be addressed. We compare the oncogenic E6 proteins of the closely related high-risk HPV31 and HPV16 with regard to their structure and their efficiency of ternary complex formation with their cellular targets p53 and E6AP, which results in p53 degradation. We solved the crystal structure of 31 E6 bound to the E6AP LxxLL motif. HPV16 E6 and 31 E6 structures are highly similar, but a few sequence variations lead to different protein contacts within the ternary complex and, as quantified here, an overall lower binding affinity of 31 E6 than 16 E6. These results align with the observed lower p53 degradation potential of 31 E6.

Dual Specificity PDZ- and 14-3-3-Binding Motifs: A Structural and Interactomics Study.

Protein-protein interaction motifs are often alterable by post-translational modifications. For example, 19% of predicted human PDZ domain-binding motifs (PBMs) have been experimentally proven to be phosphorylated, and up to 82% are theoretically phosphorylatable. Phosphorylation of PBMs may drastically rewire their interactomes, by altering their affinities for PDZ domains and 14-3-3 proteins. The effect of phosphorylation is often analyzed by performing "phosphomimetic" mutations. Here, we focused on the PBMs of HPV16-E6 viral oncoprotein and human RSK1 kinase. We measured the binding affinities of native, phosphorylated, and phosphomimetic variants of both PBMs toward the 266 human PDZ domains. We co-crystallized all the motif variants with a selected PDZ domain to characterize the structural consequence of the different modifications. Finally, we elucidated the structural basis of PBM capture by 14-3-3 proteins. This study provides novel atomic and interactomic insights into phosphorylatable dual specificity motifs and the differential effects of phosphorylation and phosphomimetic approaches.

Conformational editing of intrinsically disordered protein by α-methylation.

Intrinsically disordered proteins (IDPs) constitute a large portion of "Dark Proteome" - difficult to characterize or yet to be discovered protein structures. Here we used conformationally constrained α-methylated amino acids to bias the conformational ensemble in the free unstructured activation domain of transcriptional coactivator ACTR. Different sites and patterns of substitutions were enabled by chemical protein synthesis and led to distinct populations of α-helices. A specific substitution pattern resulted in a substantially higher binding affinity to nuclear coactivator binding domain (NCBD) of CREB-binding protein, a natural binding partner of ACTR. The first X-ray structure of the modified ACTR domain - NCBD complex visualized a unique conformation of ACTR and confirmed that the key α-methylated amino acids are localized within α-helices in the bound state. This study demonstrates a strategy for characterization of individual conformational states of IDPs.

A novel small-molecule inhibitor of the human papillomavirus E6-p53 interaction that reactivates p53 function and blocks cancer cells growth.

Despite prophylactic vaccination campaigns, human papillomavirus (HPV)-induced cancers still represent a major medical issue for global population, thus specific anti-HPV drugs are needed. Since the ability of HPV E6 oncoprotein to promote p53 degradation is linked to tumor progression, E6 has been proposed as an ideal target for cancer treatment. Using the crystal structure of the E6/E6AP/p53 complex, we performed an in silico screening of small-molecule libraries against a highly conserved alpha-helix in the N-terminal domain of E6 involved in the E6-p53 interaction. We discovered a compound able to inhibit the E6-mediated degradation of p53 through disruption of E6-p53 binding both in vitro and in cells. This compound could restore p53 intracellular levels and transcriptional activity, reduce the viability and proliferation of HPV-positive cancer cells, and block 3D cervospheres formation. Mechanistic studies revealed that the compound anti-tumor activity mainly relies on induction of cell cycle arrest and senescence. Our data demonstrate that the disruption of the direct E6-p53 interaction can be obtained with a small-molecule compound leading to specific antitumoral activity in HPV-positive cancer cells and thus represents a new approach for anti-HPV drug development.

Discovery of N-Substituted 3-Amino-4-(3-boronopropyl)pyrrolidine-3-carboxylic Acids as Highly Potent Third-Generation Inhibitors of Human Arginase I and II.

Recent efforts to identify new highly potent arginase inhibitors have resulted in the discovery of a novel family of (3R,4S)-3-amino-4-(3-boronopropyl)pyrrolidine-3-carboxylic acid analogues with up to a 1000-fold increase in potency relative to the current standards, 2-amino-6-boronohexanoic acid (ABH) and N-hydroxy-nor-l-arginine (nor-NOHA). The lead candidate, with an N-2-amino-3-phenylpropyl substituent (NED-3238), example 43, inhibits arginase I and II with IC50 values of 1.3 and 8.1 nM, respectively. Herein, we report the design, synthesis, and structure-activity relationships for this novel series of inhibitors, along with X-ray crystallographic data for selected examples bound to human arginase II.

Structural Basis of Outstanding Multivalent Effects in Jack Bean α-Mannosidase Inhibition.

Multivalent design of glycosidase inhibitors is a promising strategy for the treatment of diseases involving enzymatic hydrolysis of glycosidic bonds in carbohydrates. An essential prerequisite for successful applications is the atomic-level understanding of how outstanding binding enhancement occurs with multivalent inhibitors. Herein we report the first high-resolution crystal structures of the Jack bean α-mannosidase (JBα-man) in apo and inhibited states. The three-dimensional structure of JBα-man in complex with the multimeric cyclopeptoid-based inhibitor displaying the largest binding enhancements reported so far provides decisive insight into the molecular mechanisms underlying multivalent effects in glycosidase inhibition.

Design, synthesis, structure-activity relationships and X-ray structural studies of novel 1-oxopyrimido[4,5-c]quinoline-2-acetic acid derivatives as selective and potent inhibitors of human aldose reductase.

Human aldose reductase (AKR1B1, AR) is a key enzyme of the polyol pathway, catalyzing the reduction of glucose to sorbitol at high glucose concentrations, as those found in diabetic condition. Indeed, AKR1B1 overexpression is related to diabetes secondary complications and, in some cases, with cancer. For many years, research has been focused on finding new AKR1B1 inhibitors (ARIs) to overcome these diseases. Despite the efforts, most of the new drug candidates failed because of their poor pharmacokinetic properties and/or unacceptable side effects. Here we report the synthesis of a series of 1-oxopyrimido[4,5-c]quinoline-2-acetic acid derivatives as novel ARIs. IC50 assays and X-ray crystallographic studies proved that these compounds are promising hits for further drug development, with high potency and selectivity against AKR1B1. Based on the determined X-ray structures with hit-to-lead compounds, we designed and synthesized a second series that yielded lead compound 68 (Kiappvs. AKR1B1 = 73 nM). These compounds are related to the previously reported 2-aminopyrimido[4,5-c]quinolin-1(2H)-ones, which exhibit antimitotic activity. Regardless of their similarity, the 2-amino compounds are unable to inhibit AKR1B1 while the 2-acetic acid derivatives are not cytotoxic against fibrosarcoma HT-1080 cells. Thus, the replacement of the amino group by an acetic acid moiety changes their biological activity, improving their potency as ARIs.

Targeting Acidic Mammalian chitinase Is Effective in Animal Model of Asthma.

This article highlights our work toward the identification of a potent, selective, and efficacious acidic mammalian chitinase (AMCase) inhibitor. Rational design, guided by X-ray analysis of several inhibitors bound to human chitotriosidase (hCHIT1), led to the identification of compound 7f as a highly potent AMCase inhibitor (IC50 values of 14 and 19 nM against human and mouse enzyme, respectively) and selective (>150× against mCHIT1) with very good PK properties. This compound dosed once daily at 30 mg/kg po showed significant anti-inflammatory efficacy in HDM-induced allergic airway inflammation in mice, reducing inflammatory cell influx in the BALF and total IgE concentration in plasma, which correlated with decrease of chitinolytic activity. Therapeutic efficacy of compound 7f in the clinically relevant aeroallergen-induced acute asthma model in mice provides a rationale for developing AMCase inhibitor for the treatment of asthma.

Structural variability of E. coli thioredoxin captured in the crystal structures of single-point mutants.

Thioredoxin is a ubiquitous small protein that catalyzes redox reactions of protein thiols. Additionally, thioredoxin from E. coli (EcTRX) is a widely-used model for structure-function studies. In a previous paper, we characterized several single-point mutants of the C-terminal helix (CTH) that alter global stability of EcTRX. However, spectroscopic signatures and enzymatic activity for some of these mutants were found essentially unaffected. A comprehensive structural characterization at the atomic level of these near-invariant mutants can provide detailed information about structural variability of EcTRX. We address this point through the determination of the crystal structures of four point-mutants, whose mutations occurs within or near the CTH, namely L94A, E101G, N106A and L107A. These structures are mostly unaffected compared with the wild-type variant. Notably, the E101G mutant presents a large region with two alternative traces for the backbone of the same chain. It represents a significant shift in backbone positions. Enzymatic activity measurements and conformational dynamics studies monitored by NMR and molecular dynamic simulations show that E101G mutation results in a small effect in the structural features of the protein. We hypothesize that these alternative conformations represent samples of the native-state ensemble of EcTRX, specifically the magnitude and location of conformational heterogeneity.

Structural basis for the inhibition of AKR1B10 by the C3 brominated TTNPB derivative UVI2008.

UVI2008, a retinoic acid receptor (RAR) ß/γ agonist originated from C3 bromine addition to the parent RAR pan-agonist 4-[(E)-2-(5,6,7,8-tetrahydro-5,5,8,8-tetramethyl-2-naphthalenyl)-1-propenyl]benzoic acid (TTNPB), is also a selective inhibitor of aldo-keto reductase family member 1B10 (AKR1B10). Thus, it might become a lead drug for the design of compounds targeting both activities, as an AKR1B10 inhibitor and RAR agonist, which could constitute a novel therapeutic approach against cancer and skin-related diseases. Herein, the X-ray structure of the methylated Lys125Arg/Val301Leu AKR1B10 (i.e. AKME2MU) holoenzyme in complex with UVI2008 was determined at 1.5 Å resolution, providing an explanation for UVI2008 selectivity against AKR1B10 (IC50 = 6.1 µM) over the closely related aldose reductase (AR, IC50 = 70 µM). The carboxylic acid group of UVI2008 is located in the anion-binding pocket, at hydrogen-bond distance of catalytically important residues Tyr49 and His111. The inhibitor bromine atom can only fit in the wider active site of AKR1B10, mainly because of the native Trp112 side-chain orientation, not possible in AR. In AKR1B10, Trp112 native conformation, and thus UVI2008 binding, is facilitated through interaction with Gln114. IC50 analysis of the corresponding Thr113Gln mutant in AR confirmed this hypothesis. The elucidation of the binding mode of UVI2008 to AKR1B10, along with the previous studies on the retinoid specificity of AKR1B10 and on the stilbene retinoid scaffold conforming UVI2008, could indeed be used to foster the drug design of bifunctional antiproliferative compounds.

IDD388 Polyhalogenated Derivatives as Probes for an Improved Structure-Based Selectivity of AKR1B10 Inhibitors.

Human enzyme aldo-keto reductase family member 1B10 (AKR1B10) has evolved as a tumor marker and promising antineoplastic target. It shares high structural similarity with the diabetes target enzyme aldose reductase (AR). Starting from the potent AR inhibitor IDD388, we have synthesized a series of derivatives bearing the same halophenoxyacetic acid moiety with an increasing number of bromine (Br) atoms on its aryl moiety. Next, by means of IC50 measurements, X-ray crystallography, WaterMap analysis, and advanced binding free energy calculations with a quantum-mechanical (QM) approach, we have studied their structure-activity relationship (SAR) against both enzymes. The introduction of Br substituents decreases AR inhibition potency but improves it in the case of AKR1B10. Indeed, the Br atoms in ortho position may impede these drugs to fit into the AR prototypical specificity pocket. For AKR1B10, the smaller aryl moieties of MK181 and IDD388 can bind into the external loop A subpocket. Instead, the bulkier MK184, MK319, and MK204 open an inner specificity pocket in AKR1B10 characterized by a π-π stacking interaction of their aryl moieties and Trp112 side chain in the native conformation (not possible in AR). Among the three compounds, only MK204 can make a strong halogen bond with the protein (-4.4 kcal/mol, using QM calculations), while presenting the lowest desolvation cost among all the series, translated into the most selective and inhibitory potency AKR1B10 (IC50 = 80 nM). Overall, SAR of these IDD388 polyhalogenated derivatives have unveiled several distinctive AKR1B10 features (shape, flexibility, hydration) that can be exploited to design novel types of AKR1B10 selective drugs.

X-Ray Crystal Structure of the Full Length Human Chitotriosidase (CHIT1) Reveals Features of Its Chitin Binding Domain.

Chitinases are enzymes that catalyze the hydrolysis of chitin. Human chitotriosidase (CHIT1) is one of the two active human chitinases, involved in the innate immune response and highly expressed in a variety of diseases. CHIT1 is composed of a catalytic domain linked by a hinge to its chitin binding domain (ChBD). This latter domain belongs to the carbohydrate-binding module family 14 (CBM14 family) and facilitates binding to chitin. So far, the available crystal structures of the human chitinase CHIT1 and the Acidic Mammalian Chitinase (AMCase) comprise only their catalytic domain. Here, we report a crystallization strategy combining cross-seeding and micro-seeding cycles which allowed us to obtain the first crystal structure of the full length CHIT1 (CHIT1-FL) at 1.95 Å resolution. The CHIT1 chitin binding domain (ChBDCHIT1) structure shows a distorted ß-sandwich 3D fold, typical of CBM14 family members. Accordingly, ChBDCHIT1 presents six conserved cysteine residues forming three disulfide bridges and several exposed aromatic residues that probably are involved in chitin binding, including the highly conserved Trp465 in a surface- exposed conformation. Furthermore, ChBDCHIT1 presents a positively charged surface which may be involved in electrostatic interactions. Our data highlight the strong structural conservation of CBM14 family members and uncover the structural similarity between the human ChBDCHIT1, tachycitin and house mite dust allergens. Overall, our new CHIT1-FL structure, determined with an adapted crystallization approach, is one of the few complete bi-modular chitinase structures available and reveals the structural features of a human CBM14 domain.

The role of the C-terminal region on the oligomeric state and enzymatic activity of Trypanosoma cruzi hypoxanthine phosphoribosyl transferase.

Hypoxanthine phosphoribosyl transferase from Trypanosoma cruzi (TcHPRT) is a critical enzyme for the survival of the parasite. This work demonstrates that the full-length form in solution adopts a stable and enzymatically active tetrameric form, exhibiting large inter-subunit surfaces. Although this protein irreversibly aggregates during unfolding, oligomerization is reversible and can be modulated by low concentrations of urea. When the C-terminal region, which is predicted as a disordered stretch, is excised by proteolysis, TcHPRT adopts a dimeric state, suggesting that the C-terminal region acts as a main guide for the quaternary arrangement. These results are in agreement with X-ray crystallographic data presented in this work. On the other hand, the C-terminal region exhibits a modulatory role on the enzyme, as attested by the enhanced activity observed for the dimeric form. Bisphosphonates act as substrate-mimetics, uncovering long-range communications among the active sites. All in all, this work contributes to establish new ways applicable to the design of novel inhibitors that could eventually result in new drugs against parasitic diseases.