Biosynthesis of a clickable pyoverdine via in vivo enzyme engineering of an adenylation domain.

The engineering of non ribosomal peptide synthetases (NRPS) for new substrate specificity is a potent strategy to incorporate non-canonical amino acids into peptide sequences, thereby creating peptide diversity and broadening applications. The non-ribosomal peptide pyoverdine is the primary siderophore produced by Pseudomonas aeruginosa and holds biomedical promise in diagnosis, bio-imaging and antibiotic vectorization. We engineered the adenylation domain of PvdD, the terminal NRPS in pyoverdine biosynthesis, to accept a functionalized amino acid. Guided by molecular modeling, we rationally designed mutants of P. aeruginosa with mutations at two positions in the active site. A single amino acid change results in the successful incorporation of an azido-L-homoalanine leading to the synthesis of a new pyoverdine analog, functionalized with an azide function. We further demonstrated that copper free click chemistry is efficient on the functionalized pyoverdine and that the conjugated siderophore retains the iron chelation properties and its capacity to be recognized and transported by P. aeruginosa. The production of clickable pyoverdine holds substantial biotechnological significance, paving the way for numerous downstream applications.

Cloning of nine glucocorticoid receptor isoforms from the slender African lungfish (Protopterus dolloi).

We wanted to clone the glucocorticoid receptor (GR) from slender African lungfish (Protopterus dolloi) for comparison to the P. dolloi mineralocorticoid receptor (MR), which we had cloned and were characterizing, as well as for comparison to the GRs from humans, elephant shark and zebrafish. However, although sequencing of the genome of the Australian lungfish (Neoceratodus forsteri), as well as, that of the West African lungfish (Protopterus annectens) were reported in the first three months of 2021, we could not retrieve a GR sequence with a BLAST search of GenBank, when we submitted our research for publication in July 2021. Moreover, we were unsuccessful in cloning the GR from slender African lungfish using a cDNA from the ovary of P. dolloi and PCR primers that had successfully cloned a GR from elephant shark, Xenopus and gar GRs. On October 21, 2021 the nucleotide sequence of West African lungfish (P. annectens) GR was deposited in GenBank. We used this GR sequence to construct PCR primers that successfully cloned the GR from the slender spotted lungfish. Here, we report the sequences of nine P. dolloi GR isoforms and explain the basis for the previous failure to clone a GR from slender African lungfish using PCR primers that cloned the GR from elephant shark, Xenopus and gar. Studies are underway to determine corticosteroid activation of these slender African lungfish GRs.

Advances in Vitamin D Receptor Function and Evolution Based on the 3D Structure of the Lamprey Ligand-Binding Domain.

1α,25-dihydroxyvitamin D3 (1,25D3) regulates many physiological processes in vertebrates by binding to the vitamin D receptor (VDR). Phylogenetic analysis indicates that jawless fishes are the most basal vertebrates exhibiting a VDR gene. To elucidate the mechanism driving VDR activation during evolution, we determined the crystal structure of the VDR ligand-binding domain (LBD) complex from the basal vertebratePetromyzon marinus, sea lamprey (lVDR). Comparison of three-dimensional crystal structures of the lVDR-1,25D3 complex with higher vertebrate VDR-1,25D3 structures suggests that 1,25D3 binds to lVDR similarly to human VDR, but with unique features for lVDR around linker regions between H11 and H12 and between H9 and H10. These structural differences may contribute to the marked species differences in transcriptional responses. Furthermore, residue co-evolution analysis of VDR across vertebrates identifies amino acid positions in H9 and the large insertion domain VDR LBD specific as correlated.

Insights into mineralocorticoid receptor homodimerization from a combined molecular modeling and bioinformatics study.

In vertebrates, the mineralocorticoid receptor (MR) is a steroid-activated nuclear receptor (NR) that plays essential roles in water-electrolyte balance and blood pressure homeostasis. It belongs to the group of oxo-steroidian NRs, together with the glucocorticoid (GR), progesterone (PR), and androgen (AR) receptors. Classically, these oxo-steroidian NRs homodimerize and bind to specific genomic sequences to activate gene expression. NRs are multi-domain proteins, and dimerization is mediated by both the DNA (DBD) and ligand binding domains (LBDs), with the latter thought to provide the largest dimerization interface. However, at the structural level, the dimerization of oxo-steroidian receptors LBDs has remained largely a matter of debate and, despite their sequence homology, there is currently no consensus on a common homodimer assembly across the four receptors, that is, GR, PR, AR, and MR. Here, we examined all available MR LBD crystals using different computational methods (protein common interface database, proteins, interfaces, structures and assemblies, protein-protein interaction prediction by structural matching, and evolutionary protein-protein interface classifier, and the molecular mechanics Poisson-Boltzmann surface area method). A consensus is reached by all methods and singles out an interface mediated by helices H9, H10 and the C-terminal F domain as having characteristics of a biologically relevant assembly. Interestingly, a similar assembly was previously identified for GRα, MR closest homolog. Alternative architectures that were proposed for GRα were not observed for MR. These data call for further experimental investigations of oxo-steroid dimer architectures.

Self-Associating Peptides for Modular Bifunctional Conjugation of Tetramer Macromolecules in Living Cells.

Monitoring the assembly of macromolecules to design entities with novel properties can be achieved either chemically creating covalent bonds or by noncovalent connections using appropriate structural motifs. In this report, two self-associating peptides (named K3 and E3) that originate from p53 tetramerization domain were developed as tools for highly specific and noncovalent heterotetramerization of two biomolecules. The pairing/coupling preferences of K3 and E3 were first evaluated by molecular modeling data and confirmed using circular dichroism spectroscopy, size-exclusion chromatography, and biological assays. Regardless of the moieties fused to K3 and E3, these two peptides self-assembled into dimers of dimers to form bivalent heterotetrameric complexes that proved to be extremely stable inside living cells. The benefits of the multivalency in terms of avidity, specificity, and expanded functional activity were strikingly revealed when the proliferating cell nuclear antigen (PCNA), which is essential for DNA replication, was targeted using a heterotetramer presenting both an antibody fragment against PCNA and a specific PCNA binder peptide. In vitro heterotetramerization of these two known PCNA ligands increased their binding efficiencies to PCNA up to 80-fold compared to the best homotetramer counterpart. In cellulo, the heterotetramers were able to efficiently inhibit DNA replication and to trigger cell death. Altogether, we demonstrate that these two biselective self-assembling peptidic domains offer a versatile noncovalent conjugation method that can be easily implemented for protein engineering.

Alternative dimerization interfaces in the glucocorticoid receptor-α ligand binding domain.

BACKGROUND: Nuclear hormone receptors (NRs) constitute a large family of multi-domain ligand-activated transcription factors. Dimerization is essential for their regulation, and both DNA binding domain (DBD) and ligand binding domain (LBD) are implicated in dimerization. Intriguingly, the glucocorticoid receptor-α (GRα) presents a DBD dimeric architecture similar to that of the homologous estrogen receptor-α (ERα), but an atypical dimeric architecture for the LBD. The physiological relevance of the proposed GRα LBD dimer is a subject of debate. METHODS: We analyzed all GRα LBD homodimers observed in crystals using an energetic analysis based on the PISA and on the MM/PBSA methods and a sequence conservation analysis, using the ERα LBD dimer as a reference point. RESULTS: Several dimeric assemblies were observed for GRα LBD. The assembly generally taken to be physiologically relevant showed weak binding free energy and no significant residue conservation at the contact interface, while an alternative homodimer mediated by both helix 9 and C-terminal residues showed significant binding free energy and residue conservation. However, none of the GRα LBD assemblies found in crystals are as stable or conserved as the canonical ERα LBD dimer. GRα C-terminal sequence (F-domain) forms a steric obstacle to the canonical dimer assembly in all available structures. CONCLUSIONS: Our analysis calls for a re-examination of the currently accepted GRα homodimer structure and experimental investigations of the alternative architectures. GENERAL SIGNIFICANCE: This work questions the validity of the currently accepted architecture. This has implications for interpreting physiological data and for therapeutic design pertaining to glucocorticoid research.

Tex19 and Sectm1 concordant molecular phylogenies support co-evolution of both eutherian-specific genes.

BACKGROUND: Transposable elements (TE) have attracted much attention since they shape the genome and contribute to species evolution. Organisms have evolved mechanisms to control TE activity. Testis expressed 19 (Tex19) represses TE expression in mouse testis and placenta. In the human and mouse genomes, Tex19 and Secreted and transmembrane 1 (Sectm1) are neighbors but are not homologs. Sectm1 is involved in immunity and its molecular phylogeny is unknown. METHODS: Using multiple alignments of complete protein sequences (MACS), we inferred Tex19 and Sectm1 molecular phylogenies. Protein conserved regions were identified and folds were predicted. Finally, expression patterns were studied across tissues and species using RNA-seq public data and RT-PCR. RESULTS: We present 2 high quality alignments of 58 Tex19 and 58 Sectm1 protein sequences from 48 organisms. First, both genes are eutherian-specific, i.e., exclusively present in mammals except monotremes (platypus) and marsupials. Second, Tex19 and Sectm1 have both duplicated in Sciurognathi and Bovidae while they have remained as single copy genes in all further placental mammals. Phylogenetic concordance between both genes was significant (p-value < 0.05) and supported co-evolution and functional relationship. At the protein level, Tex19 exhibits 3 conserved regions and 4 invariant cysteines. In particular, a CXXC motif is present in the N-terminal conserved region. Sectm1 exhibits 2 invariant cysteines and an Ig-like domain. Strikingly, Tex19 C-terminal conserved region was lost in Haplorrhini primates while a Sectm1 C-terminal extra domain was acquired. Finally, we have determined that Tex19 and Sectm1 expression levels anti-correlate across the testis of several primates (ρ = -0.72) which supports anti-regulation. CONCLUSIONS: Tex19 and Sectm1 co-evolution and anti-regulated expressions support a strong functional relationship between both genes. Since Tex19 operates a control on TE and Sectm1 plays a role in immunity, Tex19 might suppress an immune response directed against cells that show TE activity in eutherian reproductive tissues.

Increased frequency of single base substitutions in a population of transcripts expressed in cancer cells.

BACKGROUND: Single Base Substitutions (SBS) that alter transcripts expressed in cancer originate from somatic mutations. However, recent studies report SBS in transcripts that are not supported by the genomic DNA of tumor cells. METHODS: We used sequence based whole genome expression profiling, namely Long-SAGE (L-SAGE) and Tag-seq (a combination of L-SAGE and deep sequencing), and computational methods to identify transcripts with greater SBS frequencies in cancer. Millions of tags produced by 40 healthy and 47 cancer L-SAGE experiments were compared to 1,959 Reference Tags (RT), i.e. tags matching the human genome exactly once. Similarly, tens of millions of tags produced by 7 healthy and 8 cancer Tag-seq experiments were compared to 8,572 RT. For each transcript, SBS frequencies in healthy and cancer cells were statistically tested for equality. RESULTS: In the L-SAGE and Tag-seq experiments, 372 and 4,289 transcripts respectively, showed greater SBS frequencies in cancer. Increased SBS frequencies could not be attributed to known Single Nucleotide Polymorphisms (SNP), catalogued somatic mutations or RNA-editing enzymes. Hypothesizing that Single Tags (ST), i.e. tags sequenced only once, were indicators of SBS, we observed that ST proportions were heterogeneously distributed across Embryonic Stem Cells (ESC), healthy differentiated and cancer cells. ESC had the lowest ST proportions, whereas cancer cells had the greatest. Finally, in a series of experiments carried out on a single patient at 1 healthy and 3 consecutive tumor stages, we could show that SBS frequencies increased during cancer progression. CONCLUSION: If the mechanisms generating the base substitutions could be known, increased SBS frequency in transcripts would be a new useful biomarker of cancer. With the reduction of sequencing cost, sequence based whole genome expression profiling could be used to characterize increased SBS frequency in patient's tumor and aid diagnostic.

[START domain-containing proteins: a review of their role in lipid transport and exchange]. / Les protéines à domaine START, des trafiquants intracellulaires de lipides.

Fifteen START domain-containing proteins exist in mammals. On the basis of their structural homology, this family is divided into several sub-families consisting mainly of non-vesicular intracellular lipid carriers. With the exception of the Thioesterase-START subfamily, the other subfamilies are represented among invertebrates. The START domain is always located in the C-terminus of the protein. It is a module of about 210 residues that binds lipids, including sterols. Cholesterol, 25-hydroxycholesterol, phosphatidylcholine, phosphatidylethanolamine and ceramides are ligands for STARD1/STARD3-6, STARD5, STARD2/STARD10, STARD10 and STARD11, respectively. The lipids or sterols bound by the remaining 7 START proteins are unknown. The START domain can be regarded as a lipid-exchange and/or a lipid-sensing domain. The START domain consists in a deep lipid-binding pocket--that shields the hydrophic ligand from the external aqueous environment--covered by a lid formed by a C-terminal alpha helix. Within the same subgroup, such as the sterols-carriers subgroup, different START domains have similar biochemical properties; however, their expression profile and their subcellular localization distinguish them and are critical for their different biological functions. START proteins act in a variety of distinct physiological processes, such as lipid transfer between intracellular compartments, lipid metabolism and modulation of signaling events. Mutation or misexpression of START proteins is linked to pathological processes, including genetic disorders, autoimmune diseases and cancers.

Combinatorial signalling controls Neurogenin2 expression at the onset of spinal neurogenesis.

A central issue during embryonic development is to define how different signals cooperate in generating unique cell types. To address this issue, we focused on the function and the regulation of the proneural gene Neurogenin2 (Neurog2) during early mouse spinal neurogenesis. We showed that Neurog2 is first expressed in cells within the neural plate anterior to the node from the 5 somite-stage. The analysis of Neurog2 mutants established a role for this gene in triggering neural differentiation during spinal cord elongation. We identified a 798 base pair enhancer element (Neurog2-798) upstream of the Neurog2 coding sequence that directs the early caudal expression of Neurog2. Embryo culture experiments showed that Retinoic Acid (RA), Sonic hedgehog (Shh) and Fibroblast Growth Factor signals act in concert on this enhancer to control the spatial and temporal induction of Neurog2. We further demonstrated by transgenesis that two RA response elements and a Gli binding site within the Neurog2-798 element are absolutely required for its activity, strongly suggesting that the regulation of Neurog2 early expression by RA and Shh signals is direct. Our data thus support a model where signal integration at the level of a single enhancer constitutes a key mechanism to control the onset of neurogenesis.

ADCK3, an ancestral kinase, is mutated in a form of recessive ataxia associated with coenzyme Q10 deficiency.

Muscle coenzyme Q(10) (CoQ(10) or ubiquinone) deficiency has been identified in more than 20 patients with presumed autosomal-recessive ataxia. However, mutations in genes required for CoQ(10) biosynthetic pathway have been identified only in patients with infantile-onset multisystemic diseases or isolated nephropathy. Our SNP-based genome-wide scan in a large consanguineous family revealed a locus for autosomal-recessive ataxia at chromosome 1q41. The causative mutation is a homozygous splice-site mutation in the aarF-domain-containing kinase 3 gene (ADCK3). Five additional mutations in ADCK3 were found in three patients with sporadic ataxia, including one known to have CoQ(10) deficiency in muscle. All of the patients have childhood-onset cerebellar ataxia with slow progression, and three of six have mildly elevated lactate levels. ADCK3 is a mitochondrial protein homologous to the yeast COQ8 and the bacterial UbiB proteins, which are required for CoQ biosynthesis. Three out of four patients tested showed a low endogenous pool of CoQ(10) in their fibroblasts or lymphoblasts, and two out of three patients showed impaired ubiquinone synthesis, strongly suggesting that ADCK3 is also involved in CoQ(10) biosynthesis. The deleterious nature of the three identified missense changes was confirmed by the introduction of them at the corresponding positions of the yeast COQ8 gene. Finally, a phylogenetic analysis shows that ADCK3 belongs to the family of atypical kinases, which includes phosphoinositide and choline kinases, suggesting that ADCK3 plays an indirect regulatory role in ubiquinone biosynthesis possibly as part of a feedback loop that regulates ATP production.

Tex19, a mammalian-specific protein with a restricted expression in pluripotent stem cells and germ line.

Although the properties of embryonic stem (ES) cells make these cells very attractive in the field of replacement therapy, the molecular mechanisms involved in the maintenance of their pluripotency are not fully characterized. Starting from the observation that most pluripotent markers are also expressed by spermatogonia stem cells, we identified Tex19 as a new potential pluripotency marker. We show that Tex19 is a mammalian-specific protein duplicated in mouse and rat, renamed Tex19.1 and Tex19.2, whereas only one form is found in human. In mouse, both forms are localized on chromosome 11 and transcribed in opposite directions. Tex19 proteins are well conserved, showing two highly conserved domains that do not present any similarity with any other known domains. We show that Tex19.2 is specifically detected in the male somatic gonad lineage, whereas Tex19.1 expression is very similar to that of Oct4. Transcripts are maternally inherited, and expression starts as soon as the early embryo and later is limited to the germ line. Tex19.1 transcripts were also detected in mouse pluripotent stem cells, and expression of Tex19.1, like that of Oct4, decreases after murine embryonic stem and germ cell differentiation. Human TEX19 was more closely related to murine Tex19.1 and was also detected in adult testis and in undifferentiated ES cells. By immunofluorescence, we found that Tex19.1 protein localizes to the nucleus of mouse ES and inner cell mass cells. All these results suggest that Tex19.1, as well as human TEX19, could be a new factor involved in the maintenance of self-renewal or pluripotency of stem cells.

SAGETTARIUS: a program to reduce the number of tags mapped to multiple transcripts and to plan SAGE sequencing stages.

SAGE (Serial Analysis of Gene Expression) experiments generate short nucleotide sequences called 'tags' which are assumed to map unambiguously to their original transcripts (1 tag to 1 transcript mapping). Nevertheless, many tags are generated that do not map to any transcript or map to multiple transcripts. Current bioinformatics resources, such as SAGEmap and TAGmapper, have focused on reducing the number of unmapped tags. Here, we describe SAGETTARIUS, a new high-throughput program that performs successive precise Nla3 and Sau3A tag to transcript mapping, based on specifically designed Virtual Tag (VT) libraries. First, SAGETTARIUS decreases the number of tags mapped to multiple transcripts. Among the various mapping resources compared, SAGETTARIUS performed the best in this respect by decreasing up to 11% the number of multiply mapped tags. Second, SAGETTARIUS allows the establishment of a guideline for SAGE experiment sequencing efforts through efficient mapping of the CRT (Cytoplasmic Ribosomal protein Transcripts)-specific tags. Using all publicly available human and mouse Nla3 SAGE experiments, we show that sequencing 100,000 tags is sufficient to map almost all CRT-specific tags and that four sequencing stages can be identified when carrying out a human or mouse SAGE project. SAGETTARIUS is web interfaced and freely accessible to academic users.

PromAn: an integrated knowledge-based web server dedicated to promoter analysis.

PromAn is a modular web-based tool dedicated to promoter analysis that integrates distinct complementary databases, methods and programs. PromAn provides automatic analysis of a genomic region with minimal prior knowledge of the genomic sequence. Prediction programs and experimental databases are combined to locate the transcription start site (TSS) and the promoter region within a large genomic input sequence. Transcription factor binding sites (TFBSs) can be predicted using several public databases and user-defined motifs. Also, a phylogenetic footprinting strategy, combining multiple alignment of large genomic sequences and assignment of various scores reflecting the evolutionary selection pressure, allows for evaluation and ranking of TFBS predictions. PromAn results can be displayed in an interactive graphical user interface, PromAnGUI. It integrates all of this information to highlight active promoter regions, to identify among the huge number of TFBS predictions those which are the most likely to be potentially functional and to facilitate user refined analysis. Such an integrative approach is essential in the face of a growing number of tools dedicated to promoter analysis in order to propose hypotheses to direct further experimental validations. PromAn is publicly available at http://bips.u-strasbg.fr/PromAn.

vALId: validation of protein sequence quality based on multiple alignment data.

The validation of sequences is essential to perform accurate phylogeny and structure/function analysis. However among the thousands of protein sequences available in the public databases, most have been predicted in silico and have not systematically undergone a quality verification. It has recently become evident that they often contain sequence errors. To address the problem of automatic protein quality control, we have developed vALId, an interactive web interfaced software. Taking advantage of high quality multiple alignments of complete protein sequences (MACS), vALId first warns about the presence of suspicious insertions, deletions (indels) and divergent segments, and second, proposes corrections based on transcripts and genome contigs. In a first evaluation test, hundreds of indels and divergent segments were randomly generated in a manually refined MACS. The sensitivity (Sn) and specificity (Sp) of indel detection were excellent (0.96) while the mean Sn(0.49) and Sp(0.56) of divergent segment delineation depended on the percent identity between sequence neighbors. In a second test, 6195 sequences in 100 MACS corresponding to different functional and structural protein families were analyzed. 65% of the sequences were in silico predictions and 44% of eukaryote predicted proteins were partially incorrect with at least one suspicious indel or divergent segment.

PipeAlign: A new toolkit for protein family analysis.

PipeAlign is a protein family analysis tool integrating a five step process ranging from the search for sequence homologues in protein and 3D structure databases to the definition of the hierarchical relationships within and between subfamilies. The complete, automatic pipeline takes a single sequence or a set of sequences as input and constructs a high-quality, validated MACS (multiple alignment of complete sequences) in which sequences are clustered into potential functional subgroups. For the more experienced user, the PipeAlign server also provides numerous options to run only a part of the analysis, with the possibility to modify the default parameters of each software module. For example, the user can choose to enter an existing multiple sequence alignment for refinement, validation and subsequent clustering of the sequences. The aim is to provide an interactive workbench for the validation, integration and presentation of a protein family, not only at the sequence level, but also at the structural and functional levels. PipeAlign is available at http://igbmc.u-strasbg.fr/PipeAlign/.

M13 endopeptidases: New conserved motifs correlated with structure, and simultaneous phylogenetic occurrence of PHEX and the bony fish.

M13 endopeptidase alignments have focused mainly on mammalian sequences and on the active site region defining the catalytic sequence signatures. Aligning all available M13 from bacteria to human on a full-length basis, we have performed a sequence analysis. This enabled us to highlight the origin and function of the M13 PHEX subtype family endopeptidase (phosphate regulating gene with homologies to endopeptidases on the X chromosome). New evolutionary conserved regions in both prokaryotes and eukaryotes have been detected and eukaryotic-specific regions clearly delineated. Using the recently solved neprilysin structure, we have observed that all new motifs, except one, localize in the spatial vicinity of the previously reported catalytic signatures. Interestingly, a highly hydrophobic pocket containing three newly reported motifs is centered by the C-terminal tryptophan residue. Extensive M13 searches in complete and in progress higher eukaryotic genomes have lead to the identification of Danio rerio as the simplest organism having PHEX. Finally, the human PHEX substrate, the parathyroid hormone-related peptide, PTHrP(107-139), is absent in bony fish: this suggests the existence of further PHEX substrates common to both bony fishes and higher vertebrates.