SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion

Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 limits endogenous DNA damage, thereby suppressing cGAS-STING-dependent signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a PD-L1 transcriptional regulatory element, thereby promoting PD-L1 expression in cancer cells. SMARCAL1 loss hinders the ability of tumor cells to induce PD-L1 in response to genomic instability, enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a promising target for cancer immunotherapy.

Discriminating physiological from non-physiological interfaces in structures of protein complexes: A community-wide study

Reliably scoring and ranking candidate models of protein complexes and assigning their oligomeric state from the structure of the crystal lattice represent outstanding challenges. A community-wide effort was launched to tackle these challenges. The latest resources on protein complexes and interfaces were exploited to derive a benchmark dataset consisting of 1677 homodimer protein crystal structures, including a balanced mix of physiological and non-physiological complexes. The non-physiological complexes in the benchmark were selected to bury a similar or larger interface area than their physiological counterparts, making it more difficult for scoring functions to differentiate between them. Next, 252 functions for scoring protein-protein interfaces previously developed by 13 groups were collected and evaluated for their ability to discriminate between physiological and non-physiological complexes. A simple consensus score generated using the best performing score of each of the 13 groups, and a cross-validated Random Forest (RF) classifier were created. Both approaches showed excellent performance, with an area under the Receiver Operating Characteristic (ROC) curve of 0.93 and 0.94, respectively, outperforming individual scores developed by different groups. Additionally, AlphaFold2 engines recalled the physiological dimers with significantly higher accuracy than the non-physiological set, lending support to the reliability of our benchmark dataset annotations. Optimizing the combined power of interface scoring functions and evaluating it on challenging benchmark datasets appears to be a promising strategy.

The FANCC-FANCE-FANCF complex is evolutionarily conserved and regulates meiotic recombination

At meiosis, programmed meiotic DNA double-strand breaks are repaired via homologous recombination, resulting in crossovers (COs). From a large excess of DNA double-strand breaks that are formed, only a small proportion gets converted into COs because of active mechanisms that restrict CO formation. The Fanconi anemia (FA) complex proteins AtFANCM, MHF1 and MHF2 were previously identified in a genetic screen as anti-CO factors that function during meiosis in Arabidopsis thaliana. Here, pursuing the same screen, we identify FANCC as a new anti-CO gene. FANCC was previously only identified in mammals because of low primary sequence conservation. We show that FANCC, and its physical interaction with FANCE-FANCF, is conserved from vertebrates to plants. Further, we show that FANCC, together with its subcomplex partners FANCE and FANCF, regulates meiotic recombination. Mutations of any of these three genes partially rescues CO-defective mutants, which is particularly marked in female meiosis. Functional loss of FANCC, FANCE, or FANCF results in synthetic meiotic catastrophe with the pro-CO factor MUS81. This work reveals that FANCC is conserved outside mammals and has an anti-CO role during meiosis together with FANCE and FANCF.

EDS1 complexes are not required for PRR responses and execute TNL-ETI from the nucleus in Nicotiana benthamiana

Heterodimeric complexes incorporating the lipase-like proteins EDS1 with PAD4 or SAG101 are central hubs in plant innate immunity. EDS1 functions encompass signal relay from TIR domain-containing intracellular NLR-type immune receptors (TNLs) towards RPW8-type helper NLRs (RNLs) and, in Arabidopsis thaliana, bolstering of signaling and resistance mediated by cell-surface pattern recognition receptors (PRRs). Increasing evidence points to the activation of EDS1 complexes by small molecule binding. We used CRISPR/Cas-generated mutant lines and agroinfiltration-based complementation assays to interrogate functions of EDS1 complexes in Nicotiana benthamiana. We did not detect impaired PRR signaling in N. benthamiana lines deficient in EDS1 complexes or RNLs. Intriguingly, in assays monitoring functions of SlEDS1-NbEDS1 complexes in N. benthamiana, mutations within the SlEDS1 catalytic triad could abolish or enhance TNL immunity. Furthermore, nuclear EDS1 accumulation was sufficient for N. benthamiana TNL (Roq1) immunity. Reinforcing PRR signaling in Arabidopsis might be a derived function of the TNL/EDS1 immune sector. Although Solanaceae EDS1 functionally depends on catalytic triad residues in some contexts, our data do not support binding of a TNL-derived small molecule in the triad environment. Whether and how nuclear EDS1 activity connects to membrane pore-forming RNLs remains unknown.

Atomic-level evolutionary information improves protein-protein interface scoring

MOTIVATION

The crucial role of protein interactions and the difficulty in characterising them experimentally strongly motivates the development of computational approaches for structural prediction. Even when protein-protein docking samples correct models, current scoring functions struggle to discriminate them from incorrect decoys. The previous incorporation of conservation and coevolution information has shown promise for improving protein-protein scoring. Here, we present a novel strategy to integrate atomic-level evolutionary information into different types of scoring functions to improve their docking discrimination.

RESULTS

: We applied this general strategy to our residue-level statistical potential from InterEvScore and to two atomic-level scores, SOAP-PP and Rosetta interface score (ISC). Including evolutionary information from as few as ten homologous sequences improves the top 10 success rates of individual atomic-level scores SOAP-PP and Rosetta ISC by respectively 6 and 13.5 percentage points, on a large benchmark of 752 docking cases. The best individual homology-enriched score reaches a top 10 success rate of 34.4%. A consensus approach based on the complementarity between different homology-enriched scores further increases the top 10 success rate to 40%.

AVAILABILITY

All data used for benchmarking and scoring results, as well as a Singularity container of the pipeline, are available at http://biodev.cea.fr/interevol/interevdata/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic evidence for the involvement of mismatch repair proteins, PMS2 and MLH3, in a late step of homologous recombination

Homologous recombination (HR) repairs DNA double-strand breaks using intact homologous sequences as template DNA. Broken DNA and intact homologous sequences form joint molecules (JMs), including Holliday junctions (HJs), as HR intermediates. HJs are resolved to form crossover and noncrossover products. A mismatch repair factor, MLH3 endonuclease, produces the majority of crossovers during meiotic HR, but it remains elusive whether mismatch repair factors promote HR in nonmeiotic cells. We disrupted genes encoding the MLH3 and PMS2 endonucleases in the human B cell line, TK6, generating null MLH3^-/- and PMS2^-/- mutant cells. We also inserted point mutations into the endonuclease motif of MLH3 and PMS2 genes, generating endonuclease death MLH3^DN/DN and PMS2^EK/EK cells. MLH3^-/- and MLH3^DN/DN cells showed a very similar phenotype, a 2.5-fold decrease in the frequency of heteroallelic HR-dependent repair of restriction enzyme-induced double-strand breaks. PMS2^-/- and PMS2^EK/EK cells showed a phenotype very similar to that of the MLH3 mutants. These data indicate that MLH3 and PMS2 promote HR as an endonuclease. The MLH3^DN/DN and PMS2^EK/EK mutations had an additive effect on the heteroallelic HR. MLH3^DN/DN/PMS2^EK/EK cells showed normal kinetics of γ-irradiation-induced Rad51 foci but a significant delay in the resolution of Rad51 foci and a 3-fold decrease in the number of cisplatin-induced sister chromatid exchanges. The ectopic expression of the Gen1 HJ re-solvase partially reversed the defective heteroallelic HR of MLH3^DN/DN/PMS2^EK/EK cells. Taken together, we propose that MLH3 and PMS2 promote HR as endonucleases, most likely by processing JMs in mammalian somatic cells.

Dynamics of Ku and bacterial non-homologous end-joining characterized using single DNA molecule analysis

We use single-molecule techniques to characterize the dynamics of prokaryotic DNA repair by non-homologous end-joining (NHEJ), a system comprised only of the dimeric Ku and Ligase D (LigD). The Ku homodimer alone forms a ∼2 s synapsis between blunt DNA ends that is increased to ∼18 s upon addition of LigD, in a manner dependent on the C-terminal arms of Ku. The synapsis lifetime increases drastically for 4 nt complementary DNA overhangs, independently of the C-terminal arms of Ku. These observations are in contrast to human Ku, which is unable to bridge either of the two DNA substrates. We also demonstrate that bacterial Ku binds the DNA ends in a cooperative manner for synapsis initiation and remains stably bound at DNA junctions for several hours after ligation is completed, indicating that a system for removal of the proteins is active in vivo. Together these experiments shed light on the dynamics of bacterial NHEJ in DNA end recognition and processing. We speculate on the evolutionary similarities between bacterial and eukaryotic NHEJ and discuss how an increased understanding of bacterial NHEJ can open the door for future antibiotic therapies targeting this mechanism.

Probing Protein Interaction Networks by Combining MS-Based Proteomics and Structural Data Integration

Protein-protein interactions play a major role in the molecular machinery of life, and various techniques such as AP-MS are dedicated to their identification. However, those techniques return lists of proteins devoid of organizational structure, not detailing which proteins interact with which others. Proposing a hierarchical view of the interactions between the members of the flat list becomes highly tedious for large data sets when done by hand. To help hierarchize this data, we introduce a new bioinformatics protocol that integrates information of the multimeric protein 3D structures available in the Protein Data Bank using remote homology detection, as well as information related to Short Linear Motifs and interaction data from the BioGRID. We illustrate on two unrelated use-cases of different complexity how our approach can be useful to decipher the network of interactions hidden in the list of input proteins, and how it provides added value compared to state-of-the-art resources such as Interactome3D or STRING. Particularly, we show the added value of using homology detection to distinguish between orthologs and paralogs, and to distinguish between core obligate and more facultative interactions. We also demonstrate the potential of considering interactions occurring through Short Linear Motifs.

Design on a Rational Basis of High-Affinity Peptides Inhibiting the Histone Chaperone ASF1

Anti-silencing function 1 (ASF1) is a conserved H3-H4 histone chaperone involved in histone dynamics during replication, transcription, and DNA repair. Overexpressed in proliferating tissues including many tumors, ASF1 has emerged as a promising therapeutic target. Here, we combine structural, computational, and biochemical approaches to design peptides that inhibit the ASF1-histone interaction. Starting from the structure of the human ASF1-histone complex, we developed a rational design strategy combining epitope tethering and optimization of interface contacts to identify a potent peptide inhibitor with a dissociation constant of 3 nM. When introduced into cultured cells, the inhibitors impair cell proliferation, perturb cell-cycle progression, and reduce cell migration and invasion in a manner commensurate with their affinity for ASF1. Finally, we find that direct injection of the most potent ASF1 peptide inhibitor in mouse allografts reduces tumor growth. Our results open new avenues to use ASF1 inhibitors as promising leads for cancer therapy.

Structural and functional characterization of the PDZ domain of the human phosphatase PTPN3 and its interaction with the human papillomavirus E6 oncoprotein

The human protein tyrosine phosphatase non-receptor type 3 (PTPN3) is a PDZ (PSD-95/Dlg/ZO-1) domain-containing phosphatase with a tumor-suppressive or a tumor-promoting role in many cancers. Interestingly, the high-risk genital human papillomavirus (HPV) types 16 and 18 target the PDZ domain of PTPN3. The presence of a PDZ binding motif (PBM) on E6 confers interaction with a number of different cellular PDZ domain-containing proteins and is a marker of high oncogenic potential. Here, we report the molecular basis of interaction between the PDZ domain of PTPN3 and the PBM of the HPV E6 protein. We combined biophysical, NMR and X-ray experiments to investigate the structural and functional properties of the PDZ domain of PTPN3. We showed that the C-terminal sequences from viral proteins encompassing a PBM interact with PTPN3-PDZ with similar affinities to the endogenous PTPN3 ligand MAP kinase p38γ. PBM binding stabilizes the PDZ domain of PTPN3. We solved the X-ray structure of the PDZ domain of PTPN3 in complex with the PBM of the HPV E6 protein. The crystal structure and the NMR chemical shift mapping of the PTPN3-PDZ/peptide complex allowed us to pinpoint the main structural determinants of recognition of the C-terminal sequence of the E6 protein and the long-range perturbations induced upon PBM binding.

The trimeric coiled-coil HSBP1 protein promotes WASH complex assembly at centrosomes

The Arp2/3 complex generates branched actin networks that exert pushing forces onto different cellular membranes. WASH complexes activate Arp2/3 complexes at the surface of endosomes and thereby fission transport intermediates containing endocytosed receptors, such as α5β1 integrins. How WASH complexes are assembled in the cell is unknown. Here, we identify the small coiled-coil protein HSBP1 as a factor that specifically promotes the assembly of a ternary complex composed of CCDC53, WASH, and FAM21 by dissociating the CCDC53 homotrimeric precursor. HSBP1 operates at the centrosome, which concentrates the building blocks. HSBP1 depletion in human cancer cell lines and in Dictyostelium amoebae phenocopies WASH depletion, suggesting a critical role of the ternary WASH complex for WASH functions. HSBP1 is required for the development of focal adhesions and of cell polarity. These defects impair the migration and invasion of tumor cells. Overexpression of HSBP1 in breast tumors is associated with increased levels of WASH complexes and with poor prognosis for patients.

Disentangling the molecular determinants for Cenp-F localization to nuclear pores and kinetochores

Cenp-F is a multifaceted protein implicated in cancer and developmental pathologies. The Cenp-F C-terminal region contains overlapping binding sites for numerous proteins that contribute to its functions throughout the cell cycle. Here, we focus on the nuclear pore protein Nup133 that interacts with Cenp-F both at nuclear pores in prophase and at kinetochores in mitosis, and on the kinase Bub1, known to contribute to Cenp-F targeting to kinetochores. By combining in silico structural modeling and yeast two-hybrid assays, we generate an interaction model between a conserved helix within the Nup133 β-propeller and a short leucine zipper-containing dimeric segment of Cenp-F. We thereby create mutants affecting the Nup133/Cenp-F interface and show that they prevent Cenp-F localization to the nuclear envelope, but not to kinetochores. Conversely, a point mutation within an adjacent leucine zipper affecting the kinetochore targeting of Cenp-F KT-core domain impairs its interaction with Bub1, but not with Nup133, identifying Bub1 as the direct KT-core binding partner of Cenp-F. Finally, we show that Cenp-E redundantly contributes together with Bub1 to the recruitment of Cenp-F to kinetochores.

FIGL1 and its novel partner FLIP form a conserved complex that regulates homologous recombination

Homologous recombination is central to repair DNA double-strand breaks, either accidently arising in mitotic cells or in a programed manner at meiosis. Crossovers resulting from the repair of meiotic breaks are essential for proper chromosome segregation and increase genetic diversity of the progeny. However, mechanisms regulating crossover formation remain elusive. Here, we identified through genetic and protein-protein interaction screens FIDGETIN-LIKE-1 INTERACTING PROTEIN (FLIP) as a new partner of the previously characterized anti-crossover factor FIDGETIN-LIKE-1 (FIGL1) in Arabidopsis thaliana. We showed that FLIP limits meiotic crossover together with FIGL1. Further, FLIP and FIGL1 form a protein complex conserved from Arabidopsis to human. FIGL1 interacts with the recombinases RAD51 and DMC1, the enzymes that catalyze the DNA strand exchange step of homologous recombination. Arabidopsis flip mutants recapitulate the figl1 phenotype, with enhanced meiotic recombination associated with change in counts of DMC1 and RAD51 foci. Our data thus suggests that FLIP and FIGL1 form a conserved complex that regulates the crucial step of strand invasion in homologous recombination.

Homologous recombination is a DNA repair mechanism that is essential to preserve the integrity of genetic information and thus to prevent cancer formation. Homologous recombination is also used during sexual reproduction to generate genetic diversity in the offspring by shuffling parental chromosomes. Here, we identified a novel protein complex (FLIP-FIGL1) that regulates homologous recombination and is conserved from plants to mammals. This suggests the existence of a novel mode of regulation at a central step of homologous recombination.

Protein-Protein Docking Using Evolutionary Information

The structural modeling of protein complexes by docking simulations has been attracting increasing interest with the rise of proteomics and of the number of experimentally identified binary interactions. Structures of unbound partners, either modeled or experimentally determined, can be used as input to sample as extensively as possible all putative binding modes and single out the most plausible ones. At the scoring step, evolutionary information contained in the joint multiple sequence alignments of both partners can provide key insights to recognize correct interfaces. Here, we describe a computational protocol based on the InterEvDock web server to exploit coevolution constraints in protein-protein docking methods. We provide methodology guidelines to prepare the input protein structures and generate improved alignments. We also explain how to extract and use the information returned by the server through the analysis of two representative examples.

Direct interaction between exocyst and Wave complexes promotes cell protrusions and motility

Coordination between membrane trafficking and actin polymerization is fundamental in cell migration, but a dynamic view of the underlying molecular mechanisms is still missing. The Rac1 GTPase controls actin polymerization at protrusions by interacting with its effector, the Wave regulatory complex (WRC). The exocyst complex, which functions in polarized exocytosis, has been involved in the regulation of cell motility. Here, we show a physical and functional connection between exocyst and WRC. Purified components of exocyst and WRC directly associate in vitro, and interactions interfaces are identified. The exocyst-WRC interaction is confirmed in cells by co-immunoprecipitation and is shown to occur independently of the Arp2/3 complex. Disruption of the exocyst-WRC interaction leads to impaired migration. By using time-lapse microscopy coupled to image correlation analysis, we visualized the trafficking of the WRC towards the front of the cell in nascent protrusions. The exocyst is necessary for WRC recruitment at the leading edge and for resulting cell edge movements. This direct link between the exocyst and WRC provides a new mechanistic insight into the spatio-temporal regulation of cell migration.

HHalign-Kbest: exploring sub-optimal alignments for remote homology comparative modeling

MOTIVATION

The HHsearch algorithm, implementing a hidden Markov model (HMM)-HMM alignment method, has shown excellent alignment performance in the so-called twilight zone (target-template sequence identity with ∼20%). However, an optimal alignment by HHsearch may contain small to large errors, leading to poor structure prediction if these errors are located in important structural elements.

RESULTS

HHalign-Kbest server runs a full pipeline, from the generation of suboptimal HMM-HMM alignments to the evaluation of the best structural models. In the HHsearch framework, it implements a novel algorithm capable of generating k-best HMM-HMM suboptimal alignments rather than only the optimal one. For large proteins, a directed acyclic graph-based implementation reduces drastically the memory usage. Improved alignments were systematically generated among the top k suboptimal alignments. To recognize them, corresponding structural models were systematically generated and evaluated with Qmean score. The method was benchmarked over 420 targets from the SCOP30 database. In the range of HHsearch probability of 20-99%, average quality of the models (TM-score) raised by 4.1-16.3% and 8.0-21.0% considering the top 1 and top 10 best models, respectively.

AVAILABILITY AND IMPLEMENTATION

http://bioserv.rpbs.univ-paris-diderot.fr/services/HHalign-Kbest/ (source code and server).

CONTACT

guerois@cea.fr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

New Insights into the Coupling between Microtubule Depolymerization and ATP Hydrolysis by Kinesin-13 Protein Kif2C

Kinesin-13 proteins depolymerize microtubules in an ATP hydrolysis-dependent manner. The coupling between these two activities remains unclear. Here, we first studied the role of the kinesin-13 subfamily-specific loop 2 and of the KVD motif at the tip of this loop. Shortening the loop, the lysine/glutamate interchange and the additional Val to Ser substitution all led to Kif2C mutants with decreased microtubule-stimulated ATPase and impaired depolymerization capability. We rationalized these results based on a structural model of the Kif2C-ATP-tubulin complex derived from the recently determined structures of kinesin-1 bound to tubulin. In this model, upon microtubule binding Kif2C undergoes a conformational change governed in part by the interaction of the KVD motif with the tubulin interdimer interface. Second, we mutated to an alanine the conserved glutamate residue of the switch 2 nucleotide binding motif. This mutation blocks motile kinesins in a post-conformational change state and inhibits ATP hydrolysis. This Kif2C mutant still depolymerized microtubules and yielded complexes of one Kif2C with two tubulin heterodimers. These results demonstrate that the structural change of Kif2C-ATP upon binding to microtubule ends is sufficient for tubulin release, whereas ATP hydrolysis is not required. Overall, our data suggest that the conformation reached by kinesin-13s upon tubulin binding is similar to that of tubulin-bound, ATP-bound, motile kinesins but that this conformation is adapted to microtubule depolymerization.

Structural insight into how the human helicase subunit MCM2 may act as a histone chaperone together with ASF1 at the replication fork

MCM2 is a subunit of the replicative helicase machinery shown to interact with histones H3 and H4 during the replication process through its N-terminal domain. During replication, this interaction has been proposed to assist disassembly and assembly of nucleosomes on DNA. However, how this interaction participates in crosstalk with histone chaperones at the replication fork remains to be elucidated. Here, we solved the crystal structure of the ternary complex between the histone-binding domain of Mcm2 and the histones H3-H4 at 2.9 Å resolution. Histones H3 and H4 assemble as a tetramer in the crystal structure, but MCM2 interacts only with a single molecule of H3-H4. The latter interaction exploits binding surfaces that contact either DNA or H2B when H3-H4 dimers are incorporated in the nucleosome core particle. Upon binding of the ternary complex with the histone chaperone ASF1, the histone tetramer dissociates and both MCM2 and ASF1 interact simultaneously with the histones forming a 1:1:1:1 heteromeric complex. Thermodynamic analysis of the quaternary complex together with structural modeling support that ASF1 and MCM2 could form a chaperoning module for histones H3 and H4 protecting them from promiscuous interactions. This suggests an additional function for MCM2 outside its helicase function as a proper histone chaperone connected to the replication pathway.

Structural basis for signaling by exclusive EDS1 heteromeric complexes with SAG101 or PAD4 in plant innate immunity

Biotrophic plant pathogens encounter a postinfection basal resistance layer controlled by the lipase-like protein enhanced disease susceptibility 1 (EDS1) and its sequence-related interaction partners, senescence-associated gene 101 (SAG101) and phytoalexin deficient 4 (PAD4). Maintainance of separate EDS1 family member clades through angiosperm evolution suggests distinct functional attributes. We report the Arabidopsis EDS1-SAG101 heterodimer crystal structure with juxtaposed N-terminal α/β hydrolase and C-terminal α-helical EP domains aligned via a large conserved interface. Mutational analysis of the EDS1-SAG101 heterodimer and a derived EDS1-PAD4 structural model shows that EDS1 signals within mutually exclusive heterocomplexes. Although there is evolutionary conservation of α/β hydrolase topology in all three proteins, a noncatalytic resistance mechanism is indicated. Instead, the respective N-terminal domains appear to facilitate binding of the essential EP domains to create novel interaction surfaces on the heterodimer. Transitions between distinct functional EDS1 heterodimers might explain the central importance and versatility of this regulatory node in plant immunity.

InterEvScore: a novel coarse-grained interface scoring function using a multi-body statistical potential coupled to evolution

MOTIVATION

Structural prediction of protein interactions currently remains a challenging but fundamental goal. In particular, progress in scoring functions is critical for the efficient discrimination of near-native interfaces among large sets of decoys. Many functions have been developed using knowledge-based potentials, but few make use of multi-body interactions or evolutionary information, although multi-residue interactions are crucial for protein-protein binding and protein interfaces undergo significant selection pressure to maintain their interactions.

RESULTS

This article presents InterEvScore, a novel scoring function using a coarse-grained statistical potential including two- and three-body interactions, which provides each residue with the opportunity to contribute in its most favorable local structural environment. Combination of this potential with evolutionary information considerably improves scoring results on the 54 test cases from the widely used protein docking benchmark for which evolutionary information can be collected. We analyze how our way to include evolutionary information gradually increases the discriminative power of InterEvScore. Comparison with several previously published scoring functions (ZDOCK, ZRANK and SPIDER) shows the significant progress brought by InterEvScore.

AVAILABILITY

http://biodev.cea.fr/interevol/interevscore

CONTACT

guerois@cea.fr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Perturbation of Arabidopsis amino acid metabolism causes incompatibility with the adapted biotrophic pathogen Hyaloperonospora arabidopsidis

Reliance of biotrophic pathogens on living plant tissues to propagate implies strong interdependence between host metabolism and nutrient uptake by the pathogen. However, factors determining host suitability and establishment of infection are largely unknown. We describe a loss-of-inhibition allele of ASPARTATE KINASE2 and a loss-of-function allele of DIHYDRODIPICOLINATE SYNTHASE2 identified in a screen for Arabidopsis thaliana mutants with increased resistance to the obligate biotrophic oomycete Hyaloperonospora arabidopsidis (Hpa). Through different molecular mechanisms, these mutations perturb amino acid homeostasis leading to overaccumulation of the Asp-derived amino acids Met, Thr, and Ile. Although detrimental for the plant, the mutations do not cause defense activation, and both mutants retain full susceptibility to the adapted obligate biotrophic fungus Golovinomyces orontii (Go). Chemical treatments mimicking the mutants' metabolic state identified Thr as the amino acid suppressing Hpa but not Go colonization. We conclude that perturbations in amino acid homeostasis render the mutant plants unsuitable as an infection substrate for Hpa. This may be explained by deployment of the same amino acid biosynthetic pathways by oomycetes and plants. Our data show that the plant host metabolic state can, in specific ways, influence the ability of adapted biotrophic strains to cause disease.

Delineation of the Xrcc4-interacting region in the globular head domain of cernunnos/XLF

In mammals, the majority of DNA double-strand breaks are processed by the nonhomologous end-joining (NHEJ) pathway, composed of seven factors: Ku70, Ku80, DNA-PKcs, Artemis, Xrcc4 (X4), DNA-ligase IV (L4), and Cernunnos/XLF. Cernunnos is part of the ligation complex, constituted by X4 and L4. To improve our knowledge on the structure and function of Cernunnos, we performed a systematic mutagenesis study on positions selected from an analysis of the recent three-dimensional structures of this factor. Ten of 27 screened mutants were nonfunctional in several DNA repair assays. Outside amino acids critical for the expression and stability of Cernunnos, we identified three amino acids (Arg(64), Leu(65), and Leu(115)) essential for the interaction with X4 and the proper function of Cernunnos. Docking the crystal structures of the two factors further validated this probable interaction surface of Cernunnos with X4.

PCSK9 dominant negative mutant results in increased LDL catabolic rate and familial hypobetalipoproteinemia

OBJECTIVE

Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a central player in the regulation of cholesterol homeostasis, increasing the low-density lipoprotein (LDL) receptor degradation. Our study aimed at exploring the pathogenic consequences in vivo and in vitro of a PCSK9 prodomain mutation found in a family with hypobetalipoproteinemia (FHBL).

METHODS AND RESULTS

A white 49-year-old diabetic man had profound FBHL (LDLC: 16 mg/dL) whereas his daughter and sister displayed a milder phenotype (LDLC 44 mg/dL and 57 mg/dL, respectively), all otherwise healthy with a normal liver function. A monoallelic PCSK9 double-mutant R104C/V114A cosegregated with FBHL, with no mutation found at other FHBL-causing loci. A dose-effect was also found in FBHL relatives for plasma APOB and PCSK9 (very-low to undetectable in proband, approximately 50% decreased in sister and daughter) and LDL catabolic rate (256% and 88% increased in proband and daughter). Transient transfection in hepatocytes showed severely impaired processing and secretion of the double mutant which acted as a dominant negative over secretion of wild-type PCSK9.

CONCLUSIONS

These results show that heterozygous PCSK9 missense mutations may associate with profound hypobetalipoproteinemia and constitute the first direct evidence in human that decrease of plasma LDLC concentrations associated to PCSK9 LOF mutations are attributable to an increased clearance rate of LDL.

Sequence search methods and scoring functions for the design of protein structures

This chapter focuses on the methods developed for the automatic or semiautomatic design of protein structures. We present several algorithms for the exploration of the sequence space and scoring of the designed models. There are now several successful designs that have been achieved using these approaches such as the stabilization of a protein fold, the stabilization of a protein-protein complex interface, and the optimization of a protein function. A rapid presentation of the methodologies is followed by a detailed analysis of two test case studies. The first one deals with the redesign of a protein hydrophobic core and the second one with the stabilization of a protein structure through mutations at the surface. The different approaches are compared and the consistency of the predictions with the experimental data are discussed. All the programs tested in these protocols are freely available through the internet and may be applied to a wide range of design issues.

Energy estimation in protein design

The progress achieved by several groups in the field of computational protein design shows that successful design methods include two major features: efficient algorithms to deal with the combinatorial exploration of sequence space and optimal energy functions to rank sequences according to their fitness for the given fold.

Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations

We have developed a computer algorithm, FOLDEF (for FOLD-X energy function), to provide a fast and quantitative estimation of the importance of the interactions contributing to the stability of proteins and protein complexes. The predictive power of FOLDEF was tested on a very large set of point mutants (1088 mutants) spanning most of the structural environments found in proteins. FOLDEF uses a full atomic description of the structure of the proteins. The different energy terms taken into account in FOLDEF have been weighted using empirical data obtained from protein engineering experiments. First, we considered a training database of 339 mutants in nine different proteins and optimised the set of parameters and weighting factors that best accounted for the changes in stability of the mutants. The predictive power of the method was then tested using a blind test mutant database of 667 mutants, as well as a database of 82 protein-protein complex mutants. The global correlation obtained for 95 % of the entire mutant database (1030 mutants) is 0.83 with a standard deviation of 0.81 kcal mol(-1) and a slope of 0.76. The present energy function uses a minimum of computational resources and can therefore easily be used in protein design algorithms, and in the field of protein structure and folding pathways prediction where one requires a fast and accurate energy function. FOLDEF is available via a web-interface at http://fold-x.embl-heidelberg.de

Protein design based on folding models

The basic rules governing the folding of small, single-domain proteins are being discovered. New algorithms that can predict the major features of the folding process give the opportunity to design and optimise protein folding in a rational way. Recent experimental works suggest that sequence-specific features should be integrated in folding models to improve their performance.

15N NMR relaxation as a probe for helical intrinsic propensity: the case of the unfolded D2 domain of annexin I

The isolated D2 domain of annexin I is unable to adopt a tertiary fold but exhibits both native and non-native residual structures. It thus constitutes an attractive model for the investigation of dynamics of partially folded states in the context of protein folding and stability. 15N relaxation parameters of the D2 domain have been acquired at three different magnetic fields, 500, 600 and 800 MHz. This enables the estimation of the contribution of conformational exchange to the relaxation parameters on the micro- to millisecond time scale, thus providing a suitable data set for the description of motions on the pico- and nanosecond time scale. The analysis of the seven spectral densities obtained (J(0), J(50 MHz), J(60 MHz), J(80 MHz), <J(500 MHz)>, <J(600 MHz)>, <J(800 MHz)>) provides complementary and meaningful results on the conformational features of the D2 domain structure previously depicted by chemical shift and NOE data. Especially, residual helix segments exhibit distinct dynamical behaviors that are related to their intrinsic helical propensity. Beside the spectral density analysis, a series of models derived from the Lipari and Szabo model-free approach are investigated. Two models containing three parameters are able to reproduce equally well the experimental data within experimental errors but provide different values of order parameters and correlation times. The inability to find a unique model to describe the data emphasizes the difficulty to use and interpret the model-free parameters in the case of partially or fully unfolded proteins consisting of a wide range of interconverting conformers.

The SH3-fold family: experimental evidence and prediction of variations in the folding pathways

To investigate the relationships between protein topology, amino acid sequence and folding mechanisms, the folding transition state of the Sso7d protein has been characterised both experimentally and theoretically. Although Sso7d protein has a similar topology to that of the SH3 domains, the structure of its transition state is different from that of alpha-spectrin and src SH3 domains previously studied. The folding algorithm, Fold-X, including an energy function with specific sequence features, accounts for these differences and reproduces with a good agreement the set of experimental phi(double dagger-U) values obtained for the three proteins. Our analysis shows that taking into account sequence features underlying protein topology is critical for an accurate prediction of the folding process.

Dodecylphosphocholine micelles as a membrane-like environment: new results from NMR relaxation and paramagnetic relaxation enhancement analysis

To further examine to what extent a dodecyl-phosphocholine (DPC) micelle mimics a phosphatidylcholine bilayer environment, we performed 13C, 2H, and 31P NMR relaxation measurements. Our data show that the dynamic behavior of DPC phosphocholine groups at low temperature (12 degrees C) corresponds to that of a phosphatidylcholine interface at high temperature (51 degrees C). In the presence of helical peptides, a PMP1 fragment, or an annexin fragment, the DPC local dynamics are not affected whereas the DPC aggregation number is increased to match an appropriate area/volume ratio for accommodating the bound peptides. We also show that quantitative measurements of paramagnetic relaxation enhancements induced by small amounts of spin-labeled phospholipids on peptide proton signals provide a meaningful insight on the location of both PMP1 and annexin fragments in DPC micelles. The paramagnetic contributions to the relaxation were extracted from intra-residue cross-peaks of NOESY spectra for both peptides. The location of each peptide in the micelles was found consistent with the corresponding relaxation data. As illustrated by the study of the PMP1 fragment, paramagnetic relaxation data also allow us to supply the missing medium-range NOEs and therefore to complete a standard conformational analysis of peptides in micelles.

A conformational equilibrium in a protein fragment caused by two consecutive capping boxes: 1H-, 13C-NMR, and mutational analysis

The conformational properties of an 18 residues peptide spanning the entire sequence, L1KTPA5QFDAD10ELRAA15MKG, of the first helix (A-helix) of domain 2 of annexin I, were thoroughly investigated. This fragment exhibits several singular features, and in particular, two successive potential capping boxes, T3xxQ6 and D8xxE11. The former corresponds to the native hydrogen bond network stabilizing the alpha helix N-terminus in the protein; the latter is a non-native capping box able to break the helix at residue D8, and is observed in the domain 2 partially folded state. Using 2D-NMR techniques, we showed that two main populations of conformers coexist in aqueous solution. The first corresponds to a single helix extending from T3 to K17. The second corresponds to a broken helix at residue Ds. Four mutants, T3A, F7A, D8A, and E11A, were designed to further analyze the role of key amino acids in the equilibrium between the two ensembles of conformers. The sensitivity of NMR parameters to account for the variations in the populations of conformers was evaluated for each peptide. Our data show the delta13Calpha chemical shift to be the most relevant parameter. We used it to estimate the population ratio in the various peptides between the two main ensembles of conformers, the full helix and the broken helix. For the WT, E11A, and F7A peptides, these ratios are respectively 35/65, 60/40, 60/40. Our results were compared to the data obtained from helix/coil transition algorithms.

Exploring the folding pathways of annexin I, a multidomain protein. I. non-native structures stabilize the partially folded state of the isolated domain 2 of annexin I

Proteins of the annexin family constitute very attractive models because of their four approximately 70 residue domains, D1 to D4, exhibiting an identical topology comprising five helix segments with only a limited sequence homology of approximately 30%. We focus on the isolated D2 domain, which is only partially folded. A detailed analysis of this equilibrium partially folded state in aqueous solution and micellar solution using 15N-1H multidimensional NMR is presented. Comparison of the residual structure of the entire domain with that of shorter fragments indicates the presence of long-range transient hydrophobic interactions that slightly stabilize the secondary structure elements. The unfolded domain tends to behave as a four-helix, rather than as a five-helix domain. The ensemble of residual structures comprises: (i) a set of native structures consisting of three regions with large helix populations, in rather sharp correspondence with A, B and E helices, and a small helix population in the second part of the C helix; (ii) a set of non-native local structures corresponding to turn-like structures stabilized by several side-chain to side-chain interactions and helix-disruptive side-chains to backbone interactions. Remarkably, residues involved in these local non-native interactions are also involved, in the native structure, in structurally important non-local interactions. During the folding process of annexin I, the local non-native interactions have to switch to native long-range interactions. This structural switch reveals the existence of a sequence-encoded regulation of the folding pathways and kinetics, and emphasizes the key role of the non-native local structures in this regulation.

Exploring the folding pathways of annexin I, a multidomain protein. II. Hierarchy in domain folding propensities may govern the folding process

In the context of exploring the relationship between sequence and folding pathways, the multi-domain proteins of the annexin family constitute very attractive models. They are constituted of four approximately 70-residue domains, named D1 to D4, with identical topologies but only limited sequence homology of approximately 30%. The domains are organized in a pseudochiral circular arrangement. Here, we report on the folding propensity of the D1 domain of annexin I obtained from overexpression in Escherichia coli. Unlike the D2 domain, which is only partially folded, the isolated D1 domain exhibits autonomous refolding in pure aqueous solution. Similarly, the D3 domain and D2-D3 module were obtained from expression in E. coli but were found to be largely unfolded. No conclusion could be drawn for the D4 domain because it was not possible to extract it from the bacterial inclusion bodies. The data allow us to propose a plausible scenario for the annexin I folding. This working model states that firstly the D1 domain folds, and the D2 and D3 domains remain partly unfolded, facilitating the docking of the D4 domain to the D1 domain. In a second step, the D1 and D4 domains dock, and D4 may fold if already not folded. The final step starts with the stabilization of the D1-D4 module. This stabilization is crucial for allowing the non-native local interactions inside the still partially unfolded D2 domain to switch to the native long-range interactions involving D4. This switch allows the complete folding of D2 and D3. The model proposes a sequential and hierarchical process for the folding of annexin I and emphasizes the role of both native framework and non-native structures in the process.

Folding properties of an annexin I domain: a 1H-15N NMR and CD study

The annexin fold consists of four 70-residue domains with markedly homologous sequences and nearly identical structures. Each domain contains five helices designated A to E. Domain 2 of annexin I was obtained by chemical synthesis including ten specifically labeled residues and studied by 1H-15N NMR and circular dichroism (CD). In pure aqueous solution this annexin domain presents, at most, 25% of residual helix secondary structure compared to 75%-85% for the native helix content and thus does not constitute an autonomous folding unit. Dodecylphosphocholine (DPC) micelles were used to provide the annexin domain with non-specific hydrophobic interactions. The structuring effect of micelles was thoroughly investigated by CD and 1H-15N NMR. Most, but not all, of the native helix secondary structure was recovered at DPC saturation. NMR data made it possible to determine the intrinsic helix propensity hierarchy of the different helix segments of the domain: A approximately B approximately E > C, D. This hierarchy is remarkably well correlated with the location of the helices in the native protein since A, B, and E helices are those in contact with the remaining parts of the protein. This result tends to support the view that, for large proteins like annexins (35 kDa), high intrinsic secondary structure propensities, at least helix propensity, in selected protein segments is necessary for a correct folding process. As a consequence this also indicates that important information concerning the folding pathway is encoded in the protein sequence.