Computing multiple sequence/structure alignments with the T-coffee package

The nsp15 Nuclease as a Good Target to Combat SARS-CoV-2: Mechanism of Action and Its Inactivation with FDA-Approved Drugs

The pandemic caused by SARS-CoV-2 is not over yet, despite all the efforts from the scientific community. Vaccination is a crucial weapon to fight this virus; however, we still urge the development of antivirals to reduce the severity and progression of the COVID-19 disease. For that, a deep understanding of the mechanisms involved in viral replication is necessary. nsp15 is an endoribonuclease critical for the degradation of viral polyuridine sequences that activate host immune sensors. This enzyme is known as one of the major interferon antagonists from SARS-CoV-2. In this work, a biochemical characterization of SARS-CoV-2 nsp15 was performed. We saw that nsp15 is active as a hexamer, and zinc can block its activity. The role of conserved residues from SARS-CoV-2 nsp15 was investigated, and N164 was found to be important for protein hexamerization and to contribute to the specificity to degrade uridines. Several chemical groups that impact the activity of this ribonuclease were also identified. Additionally, FDA-approved drugs with the capacity to inhibit the in vitro activity of nsp15 are reported in this work. This study is of utmost importance by adding highly valuable information that can be used for the development and rational design of therapeutic strategies.

Structure of the Complete Dimeric Human GDAP1 Core Domain Provides Insights into Ligand Binding and Clustering of Disease Mutations

Charcot-Marie-Tooth disease (CMT) is one of the most common inherited neurological disorders. Despite the common involvement of ganglioside-induced differentiation-associated protein 1 (GDAP1) in CMT, the protein structure and function, as well as the pathogenic mechanisms, remain unclear. We determined the crystal structure of the complete human GDAP1 core domain, which shows a novel mode of dimerization within the glutathione S-transferase (GST) family. The long GDAP1-specific insertion forms an extended helix and a flexible loop. GDAP1 is catalytically inactive toward classical GST substrates. Through metabolite screening, we identified a ligand for GDAP1, the fatty acid hexadecanedioic acid, which is relevant for mitochondrial membrane permeability and Ca²⁺ homeostasis. The fatty acid binds to a pocket next to a CMT-linked residue cluster, increases protein stability, and induces changes in protein conformation and oligomerization. The closest homologue of GDAP1, GDAP1L1, is monomeric in its full-length form. Our results highlight the uniqueness of GDAP1 within the GST family and point toward allosteric mechanisms in regulating GDAP1 oligomeric state and function.

Clarinet (CLA-1), a novel active zone protein required for synaptic vesicle clustering and release

Active zone proteins cluster synaptic vesicles at presynaptic terminals and coordinate their release. In forward genetic screens, we isolated a novel Caenorhabditis elegans active zone gene, clarinet (cla-1). cla-1 mutants exhibit defects in synaptic vesicle clustering, active zone structure and synapse number. As a result, they have reduced spontaneous vesicle release and increased synaptic depression. cla-1 mutants show defects in vesicle distribution near the presynaptic dense projection, with fewer undocked vesicles contacting the dense projection and more docked vesicles at the plasma membrane. cla-1 encodes three isoforms containing common C-terminal PDZ and C2 domains with homology to vertebrate active zone proteins Piccolo and RIM. The C-termini of all isoforms localize to the active zone. Specific loss of the ~9000 amino acid long isoform results in vesicle clustering defects and increased synaptic depression. Our data indicate that specific isoforms of clarinet serve distinct functions, regulating synapse development, vesicle clustering and release.

Comparative Protein Structure Modeling Using MODELLER

Comparative protein structure modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and how to use the ModBase database of such models, and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described. © 2016 by John Wiley & Sons, Inc.

A Phylogenomic Framework to Study the Diversity and Evolution of Stramenopiles (=Heterokonts)

Stramenopiles or heterokonts constitute one of the most speciose and diverse clades of protists. It includes ecologically important algae (such as diatoms or large multicellular brown seaweeds), as well as heterotrophic (e.g., bicosoecids, MAST groups) and parasitic (e.g., Blastocystis, oomycetes) species. Despite their evolutionary and ecological relevance, deep phylogenetic relationships among stramenopile groups, inferred mostly from small-subunit rDNA phylogenies, remain unresolved, especially for the heterotrophic taxa. Taking advantage of recently released stramenopile transcriptome and genome sequences, as well as data from the genomic assembly of the MAST-3 species Incisomonas marina generated in our laboratory, we have carried out the first extensive phylogenomic analysis of stramenopiles, including representatives of most major lineages. Our analyses, based on a large data set of 339 widely distributed proteins, strongly support a root of stramenopiles lying between two clades, Bigyra and Gyrista (Pseudofungi plus Ochrophyta). Additionally, our analyses challenge the Phaeista-Khakista dichotomy of photosynthetic stramenopiles (ochrophytes) as two groups previously considered to be part of the Phaeista (Pelagophyceae and Dictyochophyceae), branch with strong support with the Khakista (Bolidophyceae and Diatomeae). We propose a new classification of ochrophytes within the two groups Chrysista and Diatomista to reflect the new phylogenomic results. Our stramenopile phylogeny provides a robust phylogenetic framework to investigate the evolution and diversification of this group of ecologically relevant protists.

Comparative Protein Structure Modeling Using MODELLER

Comparative protein structure modeling predicts the three-dimensional structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and how to use the ModBase database of such models, and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described. © 2016 by John Wiley & Sons, Inc.

Inference of Episodic Changes in Natural Selection Acting on Protein Coding Sequences via CODEML

This unit provides protocols for using the CODEML program from the PAML package to make inferences about episodic natural selection in protein-coding sequences. The protocols cover inference tasks such as maximum likelihood estimation of selection intensity, testing the hypothesis of episodic positive selection, and identifying sites with a history of episodic evolution. We provide protocols for using the rich set of models implemented in CODEML to assess robustness, and for using bootstrapping to assess if the requirements for reliable statistical inference have been met. An example dataset is used to illustrate how the protocols are used with real protein-coding sequences. The workflow of this design, through automation, is readily extendable to a larger-scale evolutionary survey. © 2016 by John Wiley & Sons, Inc.

Evolution of Klk4 and enamel maturation in eutherians

Kallikrein-related peptidase 4 (KLK4) is a secreted serine protease that degrades residual enamel proteins to facilitate their removal by ameloblasts, which increases mineralization and hardens the enamel. Mutations in human KLK4 cause hypomaturation amelogenesis imperfecta. Enamel formed by Klk4 null mice is normal in thickness and prism structure, but the enamel layer retains proteins, is hypomineralized, and undergoes rapid attrition following tooth eruption. We searched multiple databases, retrieved Klk4 and Klk5 from various mammalian genomes, and identified Klk4 in 46 boreoeutherian genomes. In non-Boreoeutheria, Klk4 was detected in only one afrotherian genome (as a pseudogene), and not in the other six afrotherian, two xenarthran, or three marsupial genomes. In contrast, Klk5 was detected in both marsupial and eutherian mammals. Our phylogenetic and mutation rate analyses support the hypothesis that Klk4 arose from Klk5 by gene duplication near the divergence of Afrotheria, Xenarthra and Boreoeutheria, and that functionally-differentiated Klk4 survived only in Boreoeutheria. Afrotherian mammals share the feature of delayed dental eruption relative to boreoeutherian mammals. KLK4 shortens the time required for enamel maturation and could have alleviated negative selection following mutations that resulted in thicker enamel or earlier tooth eruption, without reducing enamel hardness or causing dental attrition.

Protein purification using PDZ affinity chromatography

PDZ domains function in nature as protein-binding domains within scaffold and membrane-associated proteins. They comprise approximately 90 residues and undergo specific, high-affinity interactions with complementary C-terminal peptide sequences, other PDZ domains, and/or phospholipids. We have previously shown that the specific, strong interactions of PDZ domains with their ligands make them well suited for use in affinity chromatography. This unit provides protocols for the PDZ affinity chromatography procedure that are applicable for the purification of proteins that contain PDZ domains or PDZ domain-binding ligands, either naturally or introduced by genetic engineering. We detail the preparation of affinity resins composed of PDZ domains or PDZ domain peptide ligands coupled to solid supports. These resins can be used to purify proteins containing endogenous or genetically introduced PDZ domains or ligands, eluting the proteins with free PDZ domain peptide ligands.

Comparative Modeling of Drug Target Proteins☆

In this perspective, we begin by describing the comparative protein structure modeling technique and the accuracy of the corresponding models. We then discuss the significant role that comparative prediction plays in drug discovery. We focus on virtual ligand screening against comparative models and illustrate the state-of-the-art by a number of specific examples.

Molecular replacement: tricks and treats

Molecular replacement is the method of choice for X-ray crystallographic structure determination provided that suitable structural homologues are available in the PDB. Presently, there are ~80,000 structures in the PDB (8074 were deposited in the year 2012 alone), of which ~70% have been solved by molecular replacement. For successful molecular replacement the model must cover at least 50% of the total structure and the Cα r.m.s.d. between the core model and the structure to be solved must be less than 2 Å. Here, an approach originally implemented in the CaspR server (http://www.igs.cnrs-mrs.fr/Caspr2/index.cgi) based on homology modelling to search for a molecular-replacement solution is discussed. How the use of as much information as possible from different sources can improve the model(s) is briefly described. The combination of structural information with distantly related sequences is crucial to optimize the multiple alignment that will define the boundaries of the core domains. PDB clusters (sequences with ≥30% identical residues) can also provide information on the eventual changes in conformation and will help to explore the relative orientations assumed by protein subdomains. Normal-mode analysis can also help in generating series of conformational models in the search for a molecular-replacement solution. Of course, finding a correct solution is only the first step and the accuracy of the identified solution is as important as the data quality to proceed through refinement. Here, some possible reasons for failure are discussed and solutions are proposed using a set of successful examples.

Fife, a Drosophila Piccolo-RIM homolog, promotes active zone organization and neurotransmitter release

Neuronal communication depends on the precisely orchestrated release of neurotransmitter at specialized sites called active zones (AZs). A small number of scaffolding and cytoskeletal proteins comprising the cytomatrix of the active zone (CAZ) are thought to organize the architecture and functional properties of AZs. The majority of CAZ proteins are evolutionarily conserved, underscoring the fundamental similarities in neurotransmission at all synapses. However, core CAZ proteins Piccolo and Bassoon have long been believed exclusive to vertebrates, raising intriguing questions about the conservation of the molecular mechanisms that regulate presynaptic properties. Here, we present the identification of a piccolo-rim-related gene in invertebrates, together with molecular phylogenetic analyses that indicate the encoded proteins may represent Piccolo orthologs. In accordance, we find that the Drosophila homolog, Fife, is neuronal and localizes to presynaptic AZs. To investigate the in vivo function of Fife, we generated a deletion of the fife locus. We find that evoked neurotransmitter release is substantially decreased in fife mutants and loss of fife results in motor deficits. Through morphological analysis of fife synapses, we identify underlying AZ abnormalities including pervasive presynaptic membrane detachments and reduced synaptic vesicle clustering. Our data demonstrate the conservation of a Piccolo-related protein in invertebrates and identify critical roles for Fife in regulating AZ structure and function. These findings suggest the CAZ is more conserved than previously thought, and open the door to a more complete understanding of how CAZ proteins regulate presynaptic structure and function through genetic studies in simpler model systems.

SwissRegulon, a database of genome-wide annotations of regulatory sites: recent updates

Identification of genomic regulatory elements is essential for understanding the dynamics of cellular processes. This task has been substantially facilitated by the availability of genome sequences for many species and high-throughput data of transcripts and transcription factor (TF) binding. However, rigorous computational methods are necessary to derive accurate genome-wide annotations of regulatory sites from such data. SwissRegulon (http://swissregulon.unibas.ch) is a database containing genome-wide annotations of regulatory motifs, promoters and TF binding sites (TFBSs) in promoter regions across model organisms. Its binding site predictions were obtained with rigorous Bayesian probabilistic methods that operate on orthologous regions from related genomes, and use explicit evolutionary models to assess the evidence of purifying selection on each site. New in the current version of SwissRegulon is a curated collection of 190 mammalian regulatory motifs associated with ∼340 TFs, and TFBS annotations across a curated set of ∼35 000 promoters in both human and mouse. Predictions of TFBSs for Saccharomyces cerevisiae have also been significantly extended and now cover 158 of yeast’s ∼180 TFs. All data are accessible through both an easily navigable genome browser with search functions, and as flat files that can be downloaded for further analysis.

Structure of an MmyB-like regulator from C. aurantiacus, member of a new transcription factor family linked to antibiotic metabolism in actinomycetes

Actinomycetes are important bacterial sources of antibiotics and other secondary metabolites. Many antibiotic gene clusters are controlled by pathway-specific activators that act in response to growth conditions. Here we present the crystal structure of an MmyB-like transcription regulator MltR (PDB code 3pxp) (Caur_2278) from Chloroflexus aurantiacus, in complex with a fatty acid (myristic acid). MltR is a distant homolog of the methylenomycin activator MmyB and consists of an Xre-type N-terminal DNA-binding domain and a C-terminal ligand-binding module that is related to the Per-Arnt-Sim (PAS) domain. This structure has enabled identification of a new family of bacterial transcription factors that are distributed predominantly in actinomycetes. Bioinformatics analysis of MltR and other characterized family members suggest that they are likely associated with antibiotic and fatty acid metabolism in actinomycetes. Streptomyces coelicolor SCO4944 is a candidate as an ancestral member of the family. Its ortholog in S. griseus, SGR_6891, is induced by A-factor, a γ-butyrolactone that controls antibiotic production and development, and is adjacent to the A-factor synthase gen, afsA. The location of mltR/mmyB homologs, in particular those adjacent to less well-studied antibiotic-related genes, makes them interesting genetic markers for identifying new antibiotic genes. A model for signal-triggered DNA-binding by MltR is proposed.

Programmed ribosomal frameshifting in the expression of the regulator of intestinal stem cell proliferation, adenomatous polyposis coli (APC)

A programmed ribosomal frameshift (PRF) in the decoding of APC (adenomatous polyposis coli) mRNA has been identified and characterized in Caenorhabditis worms, Drosophila and mosquitoes. The frameshift product lacks the C-terminal approximately one-third of the product of standard decoding and instead has a short sequence encoded by the -1 frame which is just 13 residues in C. elegans, but is 125 in D. melanogaster. The frameshift site is A_AA.A_AA.C in Caenorhabditids, fruit flies and the mosquitoes studied while a variant A_AA.A_AA.A is found in some other nematodes. The predicted secondary RNA structure of the downstream stimulators varies considerably in the species studied. In the twelve sequenced Drosophila genomes, it is a long stem with a four-way junction in its loop. In the five sequenced Caenorhabditis species, it is a short RNA pseudoknot with an additional stem in loop 1. The efficiency of frameshifting varies significantly, depending on the particular stimulator within the frameshift cassette, when tested with reporter constructs in rabbit reticulocyte lysates. Phylogenetic analysis of the distribution of APC programmed ribosomal frameshifting cassettes suggests it has an ancient origin and raises questions about a possibility of synthesis of alternative protein products during expression of APC in other organisms such as humans. The origin of APC as a PRF candidate emerged from a prior study of evolutionary signatures derived from comparative analysis of the 12 fly genomes. Three other proposed PRF candidates (Xbp1, CG32736, CG14047) with switches in conservation of reading frames are likely explained by mechanisms other than PRF.

A pilot study of bacterial genes with disrupted ORFs reveals a surprising profusion of protein sequence recoding mediated by ribosomal frameshifting and transcriptional realignment

Bacterial genome annotations contain a number of coding sequences (CDSs) that, in spite of reading frame disruptions, encode a single continuous polypeptide. Such disruptions have different origins: sequencing errors, frameshift, or stop codon mutations, as well as instances of utilization of nontriplet decoding. We have extracted over 1,000 CDSs with annotated disruptions and found that about 75% of them can be clustered into 64 groups based on sequence similarity. Analysis of the clusters revealed deep phylogenetic conservation of open reading frame organization as well as the presence of conserved sequence patterns that indicate likely utilization of the nonstandard decoding mechanisms: programmed ribosomal frameshifting (PRF) and programmed transcriptional realignment (PTR). Further enrichment of these clusters with additional homologous nucleotide sequences revealed over 6,000 candidate genes utilizing PRF or PTR. Analysis of the patterns of conservation apparently associated with nontriplet decoding revealed the presence of both previously characterized frameshift-prone sequences and a few novel ones. Since the starting point of our analysis was a set of genes with already annotated disruptions, it is highly plausible that in this study, we have identified only a fraction of all bacterial genes that utilize PRF or PTR. In addition to the identification of a large number of recoded genes, a surprising observation is that nearly half of them are expressed via PTR-a mechanism that, in contrast to PRF, has not yet received substantial attention.

Key enzymes of the retinoid (visual) cycle in vertebrate retina

A major goal in vision research over the past few decades has been to understand the molecular details of retinoid processing within the retinoid (visual) cycle. This includes the consequences of side reactions that result from delayed all-trans-retinal clearance and condensation with phospholipids that characterize a variety of serious retinal diseases. Knowledge of the basic retinoid biochemistry involved in these diseases is essential for development of effective therapeutics. Photoisomerization of the 11-cis-retinal chromophore of rhodopsin triggers a complex set of metabolic transformations collectively termed phototransduction that ultimately lead to light perception. Continuity of vision depends on continuous conversion of all-trans-retinal back to the 11-cis-retinal isomer. This process takes place in a series of reactions known as the retinoid cycle, which occur in photoreceptor and RPE cells. All-trans-retinal, the initial substrate of this cycle, is a chemically reactive aldehyde that can form toxic conjugates with proteins and lipids. Therefore, much experimental effort has been devoted to elucidate molecular mechanisms of the retinoid cycle and all-trans-retinal-mediated retinal degeneration, resulting in delineation of many key steps involved in regenerating 11-cis-retinal. Three particularly important reactions are catalyzed by enzymes broadly classified as acyltransferases, short-chain dehydrogenases/reductases and carotenoid/retinoid isomerases/oxygenases. This article is part of a Special Issue entitled: Retinoid and Lipid Metabolism.

A multiple-template approach to protein threading

Most threading methods predict the structure of a protein using only a single template. Due to the increasing number of solved structures, a protein without solved structure is very likely to have more than one similar template structures. Therefore, a natural question to ask is if we can improve modeling accuracy using multiple templates. This article describes a new multiple-template threading method to answer this question. At the heart of this multiple-template threading method is a novel probabilistic-consistency algorithm that can accurately align a single protein sequence simultaneously to multiple templates. Experimental results indicate that our multiple-template method can improve pairwise sequence-template alignment accuracy and generate models with better quality than single-template models even if they are built from the best single templates (P-value <10(-6)) while many popular multiple sequence/structure alignment tools fail to do so. The underlying reason is that our probabilistic-consistency algorithm can generate accurate multiple sequence/template alignments. In another word, without an accurate multiple sequence/template alignment, the modeling accuracy cannot be improved by simply using multiple templates to increase alignment coverage. Blindly tested on the CASP9 targets with more than one good template structures, our method outperforms all other CASP9 servers except two (Zhang-Server and QUARK of the same group). Our probabilistic-consistency algorithm can possibly be extended to align multiple protein/RNA sequences and structures.

Phylogenetic and functional characterization of the hAT transposon superfamily

Transposons are found in virtually all organisms and play fundamental roles in genome evolution. They can also acquire new functions in the host organism and some have been developed as incisive genetic tools for transformation and mutagenesis. The hAT transposon superfamily contains members from the plant and animal kingdoms, some of which are active when introduced into new host organisms. We have identified two new active hAT transposons, AeBuster1, from the mosquito Aedes aegypti and TcBuster from the red flour beetle Tribolium castaneum. Activity of both transposons is illustrated by excision and transposition assays performed in Drosophila melanogaster and Ae. aegypti and by in vitro strand transfer assays. These two active insect transposons are more closely related to the Buster sequences identified in humans than they are to the previously identified active hAT transposons, Ac, Tam3, Tol2, hobo, and Hermes. We therefore reexamined the structural and functional relationships of hAT and hAT-like transposase sequences extracted from genome databases and found that the hAT superfamily is divided into at least two families. This division is supported by a difference in target-site selections generated by active transposons of each family. We name these families the Ac and Buster families after the first identified transposon or transposon-like sequence in each. We find that the recently discovered SPIN transposons of mammals are located within the family of Buster elements.

Comparative protein structure modeling using Modeller

Functional characterization of a protein sequence is one of the most frequent problems in biology. This task is usually facilitated by accurate three-dimensional (3-D) structure of the studied protein. In the absence of an experimentally determined structure, comparative or homology modeling can sometimes provide a useful 3-D model for a protein that is related to at least one known protein structure. Comparative modeling predicts the 3-D structure of a given protein sequence (target) based primarily on its alignment to one or more proteins of known structure (templates). The prediction process consists of fold assignment, target-template alignment, model building, and model evaluation. This unit describes how to calculate comparative models using the program MODELLER and discusses all four steps of comparative modeling, frequently observed errors, and some applications. Modeling lactate dehydrogenase from Trichomonas vaginalis (TvLDH) is described as an example. The download and installation of the MODELLER software is also described.