Searching protein structure databases with DaliLite v.3

Mechanistic insights into the regulation of plant phosphate homeostasis by the rice SPX2 - PHR2 complex

Phosphate (Pi) starvation response (PHR) transcription factors play key roles in plant Pi homeostasis maintenance. They are negatively regulated by stand-alone SPX proteins, cellular receptors for inositol pyrophosphate (PP-InsP) nutrient messengers. How PP-InsP-bound SPX interacts with PHRs is poorly understood. Here, we report crystal structures of the rice SPX2/InsP₆/PHR2 complex and of the PHR2 DNA binding (MYB) domain in complex with target DNA at resolutions of 3.1 Å and 2.7 Å, respectively. In the SPX2/InsP₆/PHR2 complex, the signalling-active SPX2 assembles into a domain-swapped dimer conformation and binds two copies of PHR2, targeting both its coiled-coil (CC) oligomerisation domain and MYB domain. Our results reveal that the SPX2 senses PP-InsPs to inactivate PHR2 by establishing severe steric clashes with the PHR2 MYB domain, preventing DNA binding, and by disrupting oligomerisation of the PHR2 CC domain, attenuating promoter binding. Our findings rationalize how PP-InsPs activate SPX receptor proteins to target PHR family transcription factors.

SPX receptors regulate plant phosphate response via PHR transcription factors. Here, based on crystal structure analysis of rice PHR2 complexes, the authors propose that SPX2 regulates PHR2 by preventing DNA binding and oligomerisation of the PHR2 CC domain.

Fission yeast Duf89 and Duf8901 are cobalt/nickel-dependent phosphatase-pyrophosphatases that act via a covalent aspartyl-phosphate intermediate

Domain of Unknown Function 89 (DUF89) proteins are metal-dependent phosphohydrolases. Exemplary DUF89 enzymes differ in their metal and phosphosubstrate preferences. Here, we interrogated the activities and structures of two DUF89 paralogs from fission yeast-Duf89 and Duf8901. We find that Duf89 and Duf8901 are cobalt/nickel-dependent phosphohydrolases adept at hydrolyzing p-nitrophenylphosphate and PP_i. Crystal structures of metal-free Duf89 and Co²⁺-bound Duf8901 disclosed two enzyme conformations that differed with respect to the position of a three-helix module, which is either oriented away from the active site in Duf89 or forms a lid over the active site in Duf8901. Lid closure results in a 16 Å movement of Duf8901 Asp195, vis-à-vis Asp199 in Duf89, that brings Asp195 into contact with an octahedrally coordinated cobalt. Reaction of Duf8901 with BeCl₂ and NaF in the presence of divalent cations Co²⁺, Ni²⁺, or Zn²⁺ generated covalent Duf8901-(Asp248)-beryllium trifluoride (BeF₃)•Co²⁺, Duf8901-(Asp248)-BeF₃•Ni²⁺, or Duf8901-(Asp248)-BeF₃•Zn²⁺ adducts, the structures of which suggest a two-step catalytic mechanism via formation and hydrolysis of an enzyme-(aspartyl)-phosphate intermediate. Alanine mutations of Duf8901 Asp248, Asn249, Lys401, Asp286, and Asp195 that interact with BeF₃•Co²⁺ squelched p-nitrophenylphosphatase activity. A 1.8 Å structure of a Duf8901-(Asp248)-AlF₄-OH₂•Co²⁺ transition-state mimetic suggests an associative mechanism in which Asp195 and Asp363 orient and activate the water nucleophile. Whereas deletion of the duf89 gene elicited a phenotype in which expression of phosphate homeostasis gene pho1 was derepressed, deleting duf8901 did not, thereby hinting that the DUF89 paralogs have distinct functional repertoires in vivo.

Computational Study on the Dynamics of Mycobacterium Tuberculosis RNA Polymerase Assembly

Gene regulation is an intricate phenomenon involving precise function of many macromolecular complexes. Molecular basis of this phenomenon is highly complex and cannot be fully understood using a single technique. Computational approaches can play a crucial role in overall understanding of functional and mechanistic features of a protein or an assembly. Large amounts of structural data pertaining to these complexes are publicly available. In this project, we took advantage of the availability of the structural information to unravel functional intricacies of Mycobacterium tuberculosis RNA polymerase upon interaction with RbpA. In this article, we discuss how the knowledge on protein structure and dynamics can be exploited to study function using various computational tools and resources. Overall, this article provides an overview of various computational methods which can be efficiently used to understand the role of any protein. We hope especially the nonexperts in the field could benefit from our article.

Structural basis of target DNA recognition by CRISPR-Cas12k for RNA-guided DNA transposition

The type V-K CRISPR-Cas system, featured by Cas12k effector with a naturally inactivated RuvC domain and associated with Tn7-like transposon for RNA-guided DNA transposition, is a promising tool for precise DNA insertion. To reveal the mechanism underlying target DNA recognition, we determined a cryoelectron microscopy (cryo-EM) structure of Cas12k from cyanobacteria Scytonema hofmanni in complex with a single guide RNA (sgRNA) and a double-stranded target DNA. Coupled with mutagenesis and in vitro DNA transposition assay, our results revealed mechanisms for the recognition of the GGTT protospacer adjacent motif (PAM) sequence and the structural elements of Cas12k critical for RNA-guided DNA transposition. These structural and mechanistic insights should aid in the development of type V-K CRISPR-transposon systems as tools for genome editing.

The extracellular region of bovine milk butyrophilin exhibits closer structural similarity to human myelin oligodendrocyte glycoprotein than to immunological BTN family receptors

Bovine butyrophilin (BTN1A1) is an abundant type I transmembrane glycoprotein exposed on the surface of milk fat globules. We have solved the crystal structure of its extracellular region via multiple wavelength anomalous dispersion after incorporation of selenomethionine into the bacterially produced protein. The butyrophilin ectodomain exhibits two subdomains with immunoglobulin fold, each comprising a β-sandwich with a central disulfide bridge as well as one N-linked glycosylation. The fifth Cys residue at position 193 is unpaired and prone to forming disulfide crosslinks. The apparent lack of a ligand-binding site or receptor activity suggests a function predominantly as hydrophilic coat protein to prevent coagulation of the milk fat droplets. While there is less structural resemblance to members of the human butyrophilin family such as BTN3A, which play a role as immune receptors, the N-terminal bovine butyrophilin subdomain shows surprising similarity to the human myelin oligodendrocyte glycoprotein, a protein exposed on the surface of myelin sheaths. Thus, our study lends structural support to earlier hypotheses of a correlation between the consumption of cow milk and prevalence of neurological autoimmune diseases and may offer guidance for the breeding of cattle strains that express modified butyrophilin showing less immunological cross-reactivity.

TMEM106B in humans and Vac7 and Tag1 in yeast are predicted to be lipid transfer proteins

TMEM106B is an integral membrane protein of late endosomes and lysosomes involved in neuronal function, its overexpression being associated with familial frontotemporal lobar degeneration, and point mutation linked to hypomyelination. It has also been identified in multiple screens for host proteins required for productive SARS-CoV-2 infection. Because standard approaches to understand TMEM106B at the sequence level find no homology to other proteins, it has remained a protein of unknown function. Here, the standard tool PSI-BLAST was used in a nonstandard way to show that the lumenal portion of TMEM106B is a member of the late embryogenesis abundant-2 (LEA-2) domain superfamily. More sensitive tools (HMMER, HHpred, and trRosetta) extended this to predict LEA-2 domains in two yeast proteins. One is Vac7, a regulator of PI(3,5)P₂ production in the degradative vacuole, equivalent to the lysosome, which has a LEA-2 domain in its lumenal domain. The other is Tag1, another vacuolar protein, which signals to terminate autophagy and has three LEA-2 domains in its lumenal domain. Further analysis of LEA-2 structures indicated that LEA-2 domains have a long, conserved lipid-binding groove. This implies that TMEM106B, Vac7, and Tag1 may all be lipid transfer proteins in the lumen of late endocytic organelles.

Rational engineering of a human GFP-like protein scaffold for humanized targeted nanomedicines

Green fluorescent protein (GFP) is a widely used scaffold for protein-based targeted nanomedicines because of its high biocompatibility, biological neutrality and outstanding structural stability. However, being immunogenicity a major concern in the development of drug carriers, the use of exogenous proteins such as GFP in clinics might be inadequate. Here we report a human nidogen-derived protein (HSNBT), rationally designed to mimic the structural and functional properties of GFP as a scaffold for nanomedicine. For that, a GFP-like β-barrel, containing the G2 domain of the human nidogen, has been rationally engineered to obtain a biologically neutral protein that self-assembles as 10nm-nanoparticles. This scaffold is the basis of a humanized nanoconjugate, where GFP, from the well-characterized protein T22-GFP-H6, has been substituted by the nidogen-derived GFP-like HSNBT protein. The resulting construct T22-HSNBT-H6, is a humanized CXCR4-targeted nanoparticle that selectively delivers conjugated genotoxic Floxuridine into cancer CXCR4+ cells. Indeed, the administration of T22-HSNBT-H6-FdU in a CXCR4-overexpressing colorectal cancer mouse model results in an even more efficient selective antitumoral effect than that shown by its GFP-counterpart, in absence of systemic toxicity. Therefore, the newly developed GFP-like protein scaffold appears as an ideal candidate for the development of humanized protein nanomaterials and successfully supports the tumor-targeted nanoscale drug T22-HSNBT-H6-FdU. STATEMENT OF SIGNIFICANCE: Targeted nanomedicine seeks for humanized and biologically neutral protein carriers as alternative of widely used but immunogenic exogenous protein scaffolds such as green fluorescent protein (GFP). This work reports for the first time the rational engineering of a human homolog of the GFP based in the human nidogen (named HSNBT) that shows full potential to be used in humanized protein-based targeted nanomedicines. This has been demonstrated in T22-HSNBT-H6-FdU, a humanized CXCR4-targeted protein nanoconjugate able to selectively deliver its genotoxic load into cancer cells.

Completeness and Consistency in Structural Domain Classifications

Domain classifications are a useful resource for computational analysis of the protein structure, but elements of their composition are often opaque to potential users. We perform a comparative analysis of our classification ECOD against the SCOPe, SCOP2, and CATH domain classifications with respect to their constituent domain boundaries and hierarchal organization. The coverage of these domain classifications with respect to ECOD and to the PDB was assessed by structure and by sequence. We also conducted domain pair analysis to determine broad differences in hierarchy between domains shared by ECOD and other classifications. Finally, we present domains from the major facilitator superfamily (MFS) of transporter proteins and provide evidence that supports their split into domains and for multiple conformations within these families. We find that the ECOD and CATH provide the most extensive structural coverage of the PDB. ECOD and SCOPe have the most consistent domain boundary conditions, whereas CATH and SCOP2 both differ significantly.

Structural Basis for a Dual Function ATP Grasp Ligase That Installs Single and Bicyclic ω-Ester Macrocycles in a New Multicore RiPP Natural Product

Among the ribosomally synthesized and post-translationally modified peptide (RiPP) natural products, "graspetides" (formerly known as microviridins) contain macrocyclic esters and amides that are formed by ATP-grasp ligase tailoring enzymes using the side chains of Asp/Glu as acceptors and Thr/Ser/Lys as donors. Graspetides exhibit diverse patterns of macrocylization and connectivities exemplified by microviridins, that have a caged tricyclic core, and thuringin and plesiocin that feature a "hairpin topology" with cross-strand ω-ester bonds. Here, we characterize chryseoviridin, a new type of multicore RiPP encoded by Chryseobacterium gregarium DS19109 (Phylum Bacteroidetes) and solve a 2.44 Å resolution crystal structure of a quaternary complex consisting of the ATP-grasp ligase CdnC bound to ADP, a conserved leader peptide and a peptide substrate. HRMS/MS analyses show that chryseoviridin contains four consecutive five- or six-residue macrocycles ending with a microviridin-like core. The crystal structure captures respective subunits of the CdnC homodimer in the apo or substrate-bound state revealing a large conformational change in the B-domain upon substrate binding. A docked model of ATP places the γ-phosphate group within 2.8 Å of the Asp acceptor residue. The orientation of the bound substrate is consistent with a model in which macrocyclization occurs in the N- to C-terminal direction for core peptides containing multiple Thr/Ser-to-Asp macrocycles. Using systematically varied sequences, we validate this model and identify two- or three-amino acid templating elements that flank the macrolactone and are required for enzyme activity in vitro. This work reveals the structural basis for ω-ester bond formation in RiPP biosynthesis.

eCOMPASS: evaluative comparison of multiple protein alignments by statistical score

Motivation

Detecting subtle biologically relevant patterns in protein sequences often requires the construction of a large and accurate multiple sequence alignment (MSA). Methods for constructing MSAs are usually evaluated using benchmark alignments, which, however, typically contain very few sequences and are therefore inappropriate when dealing with large numbers of proteins.

Results

eCOMPASS addresses this problem using a statistical measure of relative alignment quality based on direct coupling analysis (DCA): to maintain protein structural integrity over evolutionary time, substitutions at one residue position typically result in compensating substitutions at other positions. eCOMPASS computes the statistical significance of the congruence between high scoring directly coupled pairs and 3D contacts in corresponding structures, which depends upon properly aligned homologous residues. We illustrate eCOMPASS using both simulated and real MSAs.

Availability and implementation

The eCOMPASS executable, C++ open source code and input data sets are available at https://www.igs.umaryland.edu/labs/neuwald/software/compass

Supplementary information

Supplementary data are available at Bioinformatics online.

The structure of a major surface antigen SAG19 from Eimeria tenella unifies the Eimeria SAG family

In infections by apicomplexan parasites including Plasmodium, Toxoplasma gondii, and Eimeria, host interactions are mediated by proteins including families of membrane-anchored cysteine-rich surface antigens (SAGs) and SAG-related sequences (SRS). Eimeria tenella causes caecal coccidiosis in chickens and has a SAG family with over 80 members making up 1% of the proteome. We have solved the structure of a representative E. tenella SAG, EtSAG19, revealing that, despite a low level of sequence similarity, the entire Eimeria SAG family is unified by its three-layer αβα fold which is related to that of the CAP superfamily. Furthermore, sequence comparisons show that the Eimeria SAG fold is conserved in surface antigens of the human coccidial parasite Cyclospora cayetanensis but this fold is unrelated to that of the SAGs/SRS proteins expressed in other apicomplexans including Plasmodium species and the cyst-forming coccidia Toxoplasma gondii, Neospora caninum and Besnoitia besnoiti. However, despite having very different structures, Consurf analysis showed that Eimeria SAG and Toxoplasma SRS families each exhibit marked hotspots of sequence hypervariability that map to their surfaces distal to the membrane anchor. This suggests that the primary and convergent purpose of the different structures is to provide a platform onto which sequence variability can be imposed.

NMR structure and dynamics of Q4DY78, a conserved kinetoplasid-specific protein from Trypanosoma cruzi

The 106-residue protein Q4DY78 (UniProt accession number) from Trypanosoma cruzi is highly conserved in the related kinetoplastid pathogens Trypanosoma brucei and Leishmania major. Given the essentiality of its orthologue in T. brucei, the high sequence conservation with other trypanosomatid proteins, and the low sequence similarity with mammalian proteins, Q4DY78 is an attractive protein for structural characterization. Here, we solved the structure of Q4DY78 by solution NMR and evaluated its backbone dynamics. Q4DY78 is composed of five α -helices and a small, two-stranded antiparallel β-sheet. The backbone RMSD is 0.22 ± 0.05 Å for the representative ensemble of the 20 lowest-energy structures. Q4DY78 is overall rigid, except for N-terminal residues (V⁸ to I¹⁰), residues at loop 4 (K⁵⁷ to G⁶⁵) and residues at the C-terminus (F⁸⁹ to F¹¹²). Q4DY78 has a short motif FPCAP that could potentially mediate interactions with the host cytoskeleton via interaction with EVH1 (Drosophila Enabled (Ena)/Vasodilator-stimulated phosphoprotein (VASP) homology 1) domains. Albeit Q4DY78 lacks calcium-binding motifs, its fold resembles that of eukaryotic calcium-binding proteins such as calcitracin, calmodulin, and polcacin Bet V4. We characterized this novel protein with a calcium binding fold without the capacity to bind calcium.

Structural and functional study of Legionella pneumophila effector RavA

Legionella pneumophila is an intracellular pathogen that causes Legionnaire's disease in humans. This bacterium can be found in freshwater environments as a free-living organism, but it is also an intracellular parasite of protozoa. Human infection occurs when inhaled aerosolized pathogen comes into contact with the alveolar mucosa and replicates in alveolar macrophages. Legionella enters the host cell by phagocytosis and redirects the Legionella-containing phagosomes from the phagocytic maturation pathway. These nascent phagosomes fuse with ER-derived secretory vesicles and membranes forming the Legionella-containing vacuole. Legionella subverts many host cellular processes by secreting over 300 effector proteins into the host cell via the Dot/Icm type IV secretion system. The cellular function for many Dot/Icm effectors is still unknown. Here, we present a structural and functional study of L. pneumophila effector RavA (Lpg0008). Structural analysis revealed that the RavA consists of four ~85 residue long α-helical domains with similar folds, which show only a low level of structural similarity to other protein domains. The ~90 residues long C-terminal segment is predicted to be natively unfolded. We show that during L. pneumophila infection of human cells, RavA localizes to the Golgi apparatus and to the plasma membrane. The same localization is observed when RavA is expressed in human cells. The localization signal resides within the C-terminal sequence C⁴⁰⁹ WTSFCGLF⁴¹⁷ . Yeast-two-hybrid screen using RavA as bait identified RAB11A as a potential binding partner. RavA is present in L. pneumophila strains but only distant homologs are found in other Legionella species, where the number of repeats varies.

Structural basis for a bacterial Pip system plant effector recognition protein

A number of plant-associated proteobacteria have LuxR family transcription factors that we refer to as PipR subfamily members. PipR proteins play roles in interactions between bacteria and their plant hosts, and some are important for bacterial virulence of plants. We identified an ethanolamine derivative, N-(2-hydroxyethyl)-2-(2-hydroxyethylamino) acetamide (HEHEAA), as a potent effector of PipR-mediated gene regulation in the plant endophyte Pseudomonas GM79. HEHEAA-dependent PipR activity requires an ATP-binding cassette-type active transport system, and the periplasmic substrate-binding protein (SBP) of that system binds HEHEAA. To begin to understand the molecular basis of PipR system responses to plant factors we crystallized a HEHEAA-responsive SBP in the free- and HEHEAA-bound forms. The SBP, which is similar to peptide-binding SBPs, was in a closed conformation. A narrow cavity at the interface of its two lobes is wide enough to bind HEHEAA, but it cannot accommodate peptides with side chains. The polar atoms of HEHEAA are recognized by hydrogen-bonding interactions, and additional SBP residues contribute to the binding site. This binding mode was confirmed by a structure-based mutational analysis. We also show that a closely related SBP from the plant pathogen Pseudomonas syringae pv tomato DC3000 does not recognize HEHEAA. However, a single amino acid substitution in the presumed effector-binding pocket of the P. syringae SBP converted it to a weak HEHEAA-binding protein. The P. syringae PipR depends on a plant effector for activity, and our findings imply that different PipR-associated SBPs bind different effectors.

Molecular Cloning, Purification and Characterization of Mce1R of Mycobacterium tuberculosis

The mce1 operon of Mycobacterium tuberculosis, important for lipid metabolism/transport, host cell invasion, modulation of host immune response and pathogenicity, is under the transcriptional control of Mce1R. Hence characterizing Mce1R is an important step for novel anti-tuberculosis drug discovery. The present study reports functional and in silico characterization of Mce1R. In this work, we have computationally modeled the structure of Mce1R and have validated the structure by computational and experimental methods. Mce1R has been shown to harbor the canonical VanR-like structure with a flexible N-terminal 'arm', carrying conserved positively charged residues, most likely involved in the operator DNA binding. The mce1R gene has been cloned, expressed, purified and its DNA-binding activity has been measured in vitro. The K_d value for Mce1R-operator DNA interaction has been determined to be 0.35 ± 0.02 µM which implies that Mce1R binds to DNA with moderate affinity compared to the other FCD family of regulators. So far, this is the first report for measuring the DNA-binding affinity of any VanR-type protein. Despite significant sequence similarity at the N-terminal domain, the wHTH motif of Mce1R exhibits poor conservancy of amino acid residues, critical for DNA-binding, thus results in moderate DNA-binding affinity. The N-terminal DNA-binding domain is structurally dynamic while the C-terminal domain showed significant stability and such profile of structural dynamics is most likely to be preserved in the structural orthologs of Mce1R. In addition to this, a cavity has been detected in the C-terminal domain of Mce1R which contains a few conserved residues. Comparison with other FCD family of regulators suggests that most of the conserved residues might be critical for binding to specific ligand. The max pKd value and drug score for the cavity are estimated to be 9.04 and 109 respectively suggesting that the cavity represents a suitable target site for novel anti-tuberculosis drug discovery approaches.

Jumbo Phages: A Comparative Genomic Overview of Core Functions and Adaptions for Biological Conflicts

Jumbo phages have attracted much attention by virtue of their extraordinary genome size and unusual aspects of biology. By performing a comparative genomics analysis of 224 jumbo phages, we suggest an objective inclusion criterion based on genome size distributions and present a synthetic overview of their manifold adaptations across major biological systems. By means of clustering and principal component analysis of the phyletic patterns of conserved genes, all known jumbo phages can be classified into three higher-order groups, which include both myoviral and siphoviral morphologies indicating multiple independent origins from smaller predecessors. Our study uncovers several under-appreciated or unreported aspects of the DNA replication, recombination, transcription and virion maturation systems. Leveraging sensitive sequence analysis methods, we identify novel protein-modifying enzymes that might help hijack the host-machinery. Focusing on host–virus conflicts, we detect strategies used to counter different wings of the bacterial immune system, such as cyclic nucleotide- and NAD⁺-dependent effector-activation, and prevention of superinfection during pseudolysogeny. We reconstruct the RNA-repair systems of jumbo phages that counter the consequences of RNA-targeting host effectors. These findings also suggest that several jumbo phage proteins provide a snapshot of the systems found in ancient replicons preceding the last universal ancestor of cellular life.

Classification and phylogeny for the annotation of novel eukaryotic GNAT acetyltransferases

The enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily count more than 870 000 members through all kingdoms of life and share the same structural fold. GNAT enzymes transfer an acyl moiety from acyl coenzyme A to a wide range of substrates including aminoglycosides, serotonin, glucosamine-6-phosphate, protein N-termini and lysine residues of histones and other proteins. The GNAT subtype of protein N-terminal acetyltransferases (NATs) alone targets a majority of all eukaryotic proteins stressing the omnipresence of the GNAT enzymes. Despite the highly conserved GNAT fold, sequence similarity is quite low between members of this superfamily even when substrates are similar. Furthermore, this superfamily is phylogenetically not well characterized. Thus functional annotation based on sequence similarity is unreliable and strongly hampered for thousands of GNAT members that remain biochemically uncharacterized. Here we used sequence similarity networks to map the sequence space and propose a new classification for eukaryotic GNAT acetyltransferases. Using the new classification, we built a phylogenetic tree, representing the entire GNAT acetyltransferase superfamily. Our results show that protein NATs have evolved more than once on the GNAT acetylation scaffold. We use our classification to predict the function of uncharacterized sequences and verify by in vitro protein assays that two fungal genes encode NAT enzymes targeting specific protein N-terminal sequences, showing that even slight changes on the GNAT fold can lead to change in substrate specificity. In addition to providing a new map of the relationship between eukaryotic acetyltransferases the classification proposed constitutes a tool to improve functional annotation of GNAT acetyltransferases.

Enzymes of the GCN5-related N-acetyltransferase (GNAT) superfamily transfer an acetyl group from one molecule to another. This reaction is called acetylation and is one of the most common reactions inside the cell. The GNAT superfamily counts more than 870 000 members through all kingdoms of life. Despite sharing the same fold the GNAT superfamily is very diverse in terms of amino acid sequence and substrates. The eight N-terminal acetyltransferases (NatA, NatB, etc.. to NatH) are a GNAT subtype which acetylates the free amine group of polypeptide chains. This modification is called N-terminal acetylation and is one of the most abundant protein modifications in eukaryotic cells. This subtype is also characterized by a high sequence diversity even though they share the same substrate. In addition, the phylogeny of the superfamily is not characterized. This hampers functional annotation based on sequence similarity, and discovery of novel NATs. In this work we set out to solve the problem of the classification of eukaryotic GCN5-related acetyltransferases and report the first classification framework of the superfamily. This framework can be used as a tool for annotation of all GCN5-related acetyltransferases. As an example of what can be achieved we report in this paper the computational prediction and in vitro verification of the function of two previously uncharacterized N-terminal acetyltransferases. We also report the first acetyltransferase phylogenetic tree of the GCN5 superfamily. It indicates that N-terminal acetyltransferases do not constitute one homogeneous protein family, but that the ability to bind and acetylate protein N-termini had evolved more than once on the same acetylation scaffold. We also show that even small changes in key positions can lead to altered enzyme specificity.

Experimentally based structural model of Yih1 provides insight into its function in controlling the key translational regulator Gcn2

Yeast impact homolog 1 (Yih1), or IMPACT in mammals, is part of a conserved regulatory module controlling the activity of General Control Nonderepressible 2 (Gcn2), a protein kinase that regulates protein synthesis. Yih1/IMPACT is implicated not only in many essential cellular processes, such as neuronal development, immune system regulation and the cell cycle, but also in cancer. Gcn2 must bind to Gcn1 in order to impair the initiation of protein translation. Yih1 hinders this key Gcn1-Gcn2 interaction by binding to Gcn1, thus preventing Gcn2-mediated inhibition of protein synthesis. Here, we solved the structures of the two domains of Saccharomyces cerevisiae Yih1 separately using Nuclear Magnetic Resonance and determined the relative positions of the two domains using a range of biophysical methods. Our findings support a compact structural model of Yih1 in which the residues required for Gcn1 binding are buried in the interface. This model strongly implies that Yih1 undergoes a large conformational rearrangement from a latent closed state to a primed open state to bind Gcn1. Our study provides structural insight into the interactions of Yih1 with partner molecules.

Tripartite transporters as mechanotransmitters in periplasmic alternating-access mechanisms

To remove xenobiotics from the periplasmic space, Gram-negative bacteria utilise unique tripartite efflux systems in which a molecular engine in the plasma membrane connects to periplasmic and outer membrane subunits. Substrates bind to periplasmic sections of the engine or sometimes to the periplasmic subunits. Then, the tripartite machines undergo conformational changes that allow the movement of the substrates down the substrate translocation pathway to the outside of the cell. The transmembrane (TM) domains of the tripartite resistance-nodulation-drug-resistance (RND) transporters drive these conformational changes by converting proton motive force into mechanical motion. Similarly, the TM domains of tripartite ATP-binding cassette (ABC) transporters transmit mechanical movement associated with nucleotide binding and hydrolysis at the nucleotide-binding domains to the relevant subunits in the periplasm. In this way, metabolic energy is coupled to periplasmic alternating-access mechanisms to achieve substrate transport across the outer membrane.

A Key Motif in the Cholesterol-Dependent Cytolysins Reveals a Large Family of Related Proteins

The cholesterol-dependent cytolysins (CDCs) are bacterial, β-barrel, pore-forming toxins. A central enigma of the pore-forming mechanism is how completion of the prepore is sensed to initiate its conversion to the pore. We identified a motif that is conserved between the CDCs and a diverse family of nearly 300 uncharacterized proteins present in over 220 species that span at least 10 bacterial and 2 eukaryotic phyla. Except for this motif, these proteins exhibit little similarity to the CDCs at the primary structure level. Studies herein show this motif is a critical component of the sensor that initiates the prepore-to-pore transition in the CDCs. We further show by crystallography, single particle analysis, and biochemical studies of one of these CDC-like (CDCL) proteins from Elizabethkingia anophelis, a commensal of the malarial mosquito midgut, that a high degree of structural similarity exists between the CDC and CDCL monomer structures and both form large oligomeric pore complexes. Furthermore, the conserved motif in the E. anophelis CDCL crystal structure occupies a nearly identical position and makes similar contacts to those observed in the structure of the archetype CDC, perfringolysin O (PFO). This suggests a common function in the CDCs and CDCLs and may explain why only this motif is conserved in the CDCLs. Hence, these studies identify a critical component of the sensor involved in initiating the prepore-to-pore transition in the CDCs, which is conserved in a large and diverse group of distant relatives of the CDCs.IMPORTANCE The cholesterol-dependent cytolysins' pore-forming mechanism relies on the ability to sense the completion of the oligomeric prepore structure and initiate the insertion of the β-barrel pore from the assembled prepore structure. These studies show that a conserved motif is an important component of the sensor that triggers the prepore-to-pore transition and that it is conserved in a large family of previously unidentified CDC-like proteins, the genes for which are present in a vast array of microbial species that span most terrestrial environments, as well as most animal and human microbiomes. These studies establish the foundation for future investigations that will probe the contribution of this large family of CDC-like proteins to microbial survival and human disease.

Methionine aminopeptidases with short sequence inserts within the catalytic domain are differentially inhibited: Structural and biochemical studies of three proteins from Vibrio spp

Methionine aminopeptidases (MetAPs) have been recognized as drug targets and have been extensively studied for discovery of selective inhibitors. MetAPs are essential enzymes in all living cells. While most prokaryotes contain a single gene, some prokaryotes and all eukaryotes including human have redundancy. Due to the similarity in the active sites of the MetAP enzyme between the pathogens and human limited the success of discovering selective inhibitors. We recently have discovered that MetAPs with small inserts within the catalytic domain to have different susceptibilities against some inhibitors compared to those that do not have. Using this clue we used bioinformatic tools to identify new variants of MetAPs with inserts in pathogenic species. Two new isoforms were identified in Vibrio species with two and three inserts in addition to an isoform without any insert. Multiple sequence alignment suggested that inserts are conserved in several of the Vibrio species. Two of the three inserts are common between two and three insert isoforms. One of the inserts is identified to have "NNKNN" motif that is similar to well-characterized quorum sensing peptide, "NNWNN". Another insert is predicted to have a posttranslational modification site. Three Vibrio proteins were cloned, expressed, purified, enzyme kinetics established and inhibitor screening has been performed. Several of the pyridinylpyrimidine derivatives selectively inhibited MetAPs with inserts compared to those that do not have, including the human enzyme. Crystal structure and molecular modeling studies provide the molecular basis for selective inhibition.

Mechanisms for target recognition and cleavage by the Cas12i RNA-guided endonuclease

Cas12i is a recently identified type V CRISPR-Cas endonuclease that predominantly cleaves the non-target strand of a double-stranded DNA substrate. This nicking activity of Cas12i could potentially be used for genome editing with high specificity. To elucidate its mechanisms for target recognition and cleavage, we determined cryo-EM structures of Cas12i in multiple functional states. Cas12i pre-orders a seven-nucleotide seed sequence of the crRNA for target recognition and undergoes a two-step activation through crRNA-DNA hybridization. Formation of 14 base pairs activates the nickase activity, and 28-bp hybridization promotes cleavage of the target strand. The atomic structures and mechanistic insights gained should facilitate the manipulation of Cas12i for genome editing applications.

Structure of the DOCK2-ELMO1 complex provides insights into regulation of the auto-inhibited state

DOCK (dedicator of cytokinesis) proteins are multidomain guanine nucleotide exchange factors (GEFs) for RHO GTPases that regulate intracellular actin dynamics. DOCK proteins share catalytic (DOCKDHR2) and membrane-associated (DOCKDHR1) domains. The structurally-related DOCK1 and DOCK2 GEFs are specific for RAC, and require ELMO (engulfment and cell motility) proteins for function. The N-terminal RAS-binding domain (RBD) of ELMO (ELMORBD) interacts with RHOG to modulate DOCK1/2 activity. Here, we determine the cryo-EM structures of DOCK2-ELMO1 alone, and as a ternary complex with RAC1, together with the crystal structure of a RHOG-ELMO2RBD complex. The binary DOCK2-ELMO1 complex adopts a closed, auto-inhibited conformation. Relief of auto-inhibition to an active, open state, due to a conformational change of the ELMO1 subunit, exposes binding sites for RAC1 on DOCK2DHR2, and RHOG and BAI GPCRs on ELMO1. Our structure explains how up-stream effectors, including DOCK2 and ELMO1 phosphorylation, destabilise the auto-inhibited state to promote an active GEF.

Real time structural search of the Protein Data Bank

Detection of protein structure similarity is a central challenge in structural bioinformatics. Comparisons are usually performed at the polypeptide chain level, however the functional form of a protein within the cell is often an oligomer. This fact, together with recent growth of oligomeric structures in the Protein Data Bank (PDB), demands more efficient approaches to oligomeric assembly alignment/retrieval. Traditional methods use atom level information, which can be complicated by the presence of topological permutations within a polypeptide chain and/or subunit rearrangements. These challenges can be overcome by comparing electron density volumes directly. But, brute force alignment of 3D data is a compute intensive search problem. We developed a 3D Zernike moment normalization procedure to orient electron density volumes and assess similarity with unprecedented speed. Similarity searching with this approach enables real-time retrieval of proteins/protein assemblies resembling a target, from PDB or user input, together with resulting alignments (http://shape.rcsb.org).

Structure, Dynamics and Cellular Insight Into Novel Substrates of the Legionella pneumophila Type II Secretion System

Legionella pneumophila is a Gram-negative bacterium that is able to replicate within a broad range of aquatic protozoan hosts. L. pneumophila is also an opportunistic human pathogen that can infect macrophages and epithelia in the lung and lead to Legionnaires’ disease. The type II secretion system is a key virulence factor of L. pneumophila and is used to promote bacterial growth at low temperatures, regulate biofilm formation, modulate host responses to infection, facilitate bacterial penetration of mucin gels and is necessary for intracellular growth during the initial stages of infection. The L. pneumophila type II secretion system exports at least 25 substrates out of the bacterium and several of these, including NttA to NttG, contain unique amino acid sequences that are generally not observed outside of the Legionella genus. NttA, NttC, and NttD are required for infection of several amoebal species but it is unclear what influence other novel substrates have within their host. In this study, we show that NttE is required for optimal infection of Acanthamoeba castellanii and Vermamoeba vermiformis amoeba and is essential for the typical colony morphology of L. pneumophila. In addition, we report the atomic structures of NttA, NttC, and NttE and through a combined biophysical and biochemical hypothesis driven approach we propose novel functions for these substrates during infection. This work lays the foundation for future studies into the mechanistic understanding of novel type II substrate functions and how these relate to L. pneumophila ecology and disease.

The Legionella kinase LegK7 exploits the Hippo pathway scaffold protein MOB1A for allostery and substrate phosphorylation

During infection, the bacterial pathogen Legionella pneumophila manipulates a variety of host cell signaling pathways, including the Hippo pathway which controls cell proliferation and differentiation in eukaryotes. Our previous studies revealed that L. pneumophila encodes the effector kinase LegK7 which phosphorylates MOB1A, a highly conserved scaffold protein of the Hippo pathway. Here, we show that MOB1A, in addition to being a substrate of LegK7, also functions as an allosteric activator of its kinase activity. A crystallographic analysis of the LegK7-MOB1A complex revealed that the N-terminal half of LegK7 is structurally similar to eukaryotic protein kinases, and that MOB1A directly binds to the LegK7 kinase domain. Substitution of interface residues critical for complex formation abrogated allosteric activation of LegK7 both in vitro and within cells and diminished MOB1A phosphorylation. Importantly, the N-terminal extension (NTE) of MOB1A not only regulated complex formation with LegK7 but also served as a docking site for downstream substrates such as the transcriptional coregulator YAP1. Deletion of the NTE from MOB1A or addition of NTE peptides as binding competitors attenuated YAP1 recruitment to and phosphorylation by LegK7. By providing mechanistic insight into the formation and regulation of the LegK7-MOB1A complex, our study unravels a sophisticated molecular mimicry strategy that is used by L. pneumophila to take control of the host cell Hippo pathway.

Structural and functional evidence of bacterial antiphage protection by Thoeris defense system via NAD+ degradation

The intense arms race between bacteria and phages has led to the development of diverse antiphage defense systems in bacteria. Unlike well-known restriction-modification and CRISPR-Cas systems, recently discovered systems are poorly characterized. One such system is the Thoeris defense system, which consists of two genes, thsA and thsB. Here, we report structural and functional analyses of ThsA and ThsB. ThsA exhibits robust NAD⁺ cleavage activity and a two-domain architecture containing sirtuin-like and SLOG-like domains. Mutation analysis suggests that NAD⁺ cleavage is linked to the antiphage function of Thoeris. ThsB exhibits a structural resemblance to TIR domain proteins such as nucleotide hydrolases and Toll-like receptors, but no enzymatic activity is detected in our in vitro assays. These results further our understanding of the molecular mechanism underlying the Thoeris defense system, highlighting a unique strategy for bacterial antiphage resistance via NAD⁺ degradation.

The Thoeris defense system is a recently discovered bacterial defense system that protects bacteria against phage infection and consists of the two genes thsA and thsB. Here, the authors present the crystal structures of Bacillus cereus ThsA and ThsB and show that ThsA is a NAD⁺ cleaving enzyme.

Atomic structures of the RNA end-healing 5'-OH kinase and 2',3'-cyclic phosphodiesterase domains of fungal tRNA ligase: conformational switches in the kinase upon binding of the GTP phosphate donor

Fungal tRNA ligase (Trl1) rectifies RNA breaks with 2′,3′-cyclic-PO₄ and 5′-OH termini. Trl1 consists of three catalytic modules: an N-terminal ligase (LIG) domain; a central polynucleotide kinase (KIN) domain; and a C-terminal cyclic phosphodiesterase (CPD) domain. Trl1 enzymes found in all human fungal pathogens are untapped targets for antifungal drug discovery. Here we report a 1.9 Å crystal structure of Trl1 KIN-CPD from the pathogenic fungus Candida albicans, which adopts an extended conformation in which separate KIN and CPD domains are connected by an unstructured linker. CPD belongs to the 2H phosphotransferase superfamily by dint of its conserved central concave β sheet and interactions of its dual HxT motif histidines and threonines with phosphate in the active site. Additional active site motifs conserved among the fungal CPD clade of 2H enzymes are identified. We present structures of the Candida Trl1 KIN domain at 1.5 to 2.0 Å resolution—as apoenzyme and in complexes with GTP•Mg²⁺, IDP•PO₄, and dGDP•PO₄—that highlight conformational switches in the G-loop (which recognizes the guanine base) and lid-loop (poised over the nucleotide phosphates) that accompany nucleotide binding.

The large family of PC4-like domains - similar folds and functions throughout all kingdoms of life

RNA- and DNA-binding domains are essential building blocks for specific regulation of gene expression. While a number of canonical nucleic acid binding domains share sequence and structural conservation, others are less obviously linked by evolutionary traits. In this review, we describe a protein fold of about 150 aa in length, bearing a conserved β-β-β-β-α-linker-β-β-β-β-α topology and similar nucleic acid binding properties but no apparent sequence conservation. The same overall fold can also be achieved by dimerization of two proteins, each bearing a β-β-β-β-α topology. These proteins include but are not limited to the transcription factors PC4 and P24 from humans and plants, respectively, the human RNA-transport factor Pur-α (also termed PURA), as well as the ssDNA-binding SP_0782 protein from Streptococcus pneumonia and the bacteriophage coat proteins PP7 and MS2. Besides their common overall topology, these proteins share common nucleic acids binding surfaces and thus functional similarity. We conclude that these PC4-like domains include proteins from all kingdoms of life and are much more abundant than previously known.

Novel Immunoglobulin Domain Proteins Provide Insights into Evolution and Pathogenesis of SARS-CoV-2-Related Viruses

A novel coronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was recently identified as the causative agent for the coronavirus disease 2019 (COVID-19) outbreak that has generated a global health crisis. We use a combination of genomic analysis and sensitive profile-based sequence and structure analysis to understand the potential pathogenesis determinants of this virus. As a result, we identify several fast-evolving genomic regions that might be at the interface of virus-host interactions, corresponding to the receptor binding domain of the Spike protein, the three tandem Macro fold domains in ORF1a, and the uncharacterized protein ORF8. Further, we show that ORF8 and several other proteins from alpha- and beta-CoVs belong to novel families of immunoglobulin (Ig) proteins. Among them, ORF8 is distinguished by being rapidly evolving, possessing a unique insert, and having a hypervariable position among SARS-CoV-2 genomes in its predicted ligand-binding groove. We also uncover numerous Ig domain proteins from several unrelated metazoan viruses, which are distinct in sequence and structure but share comparable architectures to those of the CoV Ig domain proteins. Hence, we propose that SARS-CoV-2 ORF8 and other previously unidentified CoV Ig domain proteins fall under the umbrella of a widespread strategy of deployment of Ig domain proteins in animal viruses as pathogenicity factors that modulate host immunity. The rapid evolution of the ORF8 Ig domain proteins points to a potential evolutionary arms race between viruses and hosts, likely arising from immune pressure, and suggests a role in transmission between distinct host species.IMPORTANCE The ongoing COVID-19 pandemic strongly emphasizes the need for a more complete understanding of the biology and pathogenesis of its causative agent SARS-CoV-2. Despite intense scrutiny, several proteins encoded by the genomes of SARS-CoV-2 and other SARS-like coronaviruses remain enigmatic. Moreover, the high infectivity and severity of SARS-CoV-2 in certain individuals make wet-lab studies currently challenging. In this study, we used a series of computational strategies to identify several fast-evolving regions of SARS-CoV-2 proteins which are potentially under host immune pressure. Most notably, the hitherto-uncharacterized protein encoded by ORF8 is one of them. Using sensitive sequence and structural analysis methods, we show that ORF8 and several other proteins from alpha- and beta-coronavirus comprise novel families of immunoglobulin domain proteins, which might function as potential immune modulators to delay or attenuate the host immune response against the viruses.

Structural and Functional Analyses of the FAM46C/Plk4 Complex

FAM46C, a non-canonical poly(A) polymerase, is frequently mutated in multiple myeloma. Loss of function of FAM46C promotes cell survival of multiple myeloma, suggesting a tumor-suppressive role. FAM46C is also essential for fastening sperm head and flagellum, indispensable for male fertility. The molecular mechanisms of these functions of FAM46C remain elusive. We report the crystal structure of FAM46C to provide the basis for its poly(A) polymerase activity and rationalize mutations associated with multiple myeloma. In addition, we found that FAM46C interacts directly with the serine/threonine kinase Plk4, the master regulator of centrosome duplication. We present the structure of FAM46C in complex with the Cryptic Polo-Box 1-2 domains of Plk4. Our structure-based mutational analyses show that the interaction with Plk4 recruits FAM46C to centrosomes. Our data suggest that Plk4-mediated localization of FAM46C enables its regulation of centrosome structure and functions, which may underlie the roles for FAM46C in cell proliferation and sperm development.

Mycobacterial DNA polymerase I: activities and crystal structures of the POL domain as apoenzyme and in complex with a DNA primer-template and of the full-length FEN/EXO-POL enzyme

Mycobacterial Pol1 is a bifunctional enzyme composed of an N-terminal DNA flap endonuclease/5' exonuclease domain (FEN/EXO) and a C-terminal DNA polymerase domain (POL). Here we document additional functions of Pol1: FEN activity on the flap RNA strand of an RNA:DNA hybrid and reverse transcriptase activity on a DNA-primed RNA template. We report crystal structures of the POL domain, as apoenzyme and as ternary complex with 3'-dideoxy-terminated DNA primer-template and dNTP. The thumb, palm, and fingers subdomains of POL form an extensive interface with the primer-template and the triphosphate of the incoming dNTP. Progression from an open conformation of the apoenzyme to a nearly closed conformation of the ternary complex entails a disordered-to-ordered transition of several segments of the thumb and fingers modules and an inward motion of the fingers subdomain-especially the O helix-to engage the primer-template and dNTP triphosphate. Distinctive structural features of mycobacterial Pol1 POL include a manganese binding site in the vestigial 3' exonuclease subdomain and a non-catalytic water-bridged magnesium complex at the protein-DNA interface. We report a crystal structure of the bifunctional FEN/EXO-POL apoenzyme that reveals the positions of two active site metals in the FEN/EXO domain.

Substrate Interaction with the EssC Coupling Protein of the Type VIIb Secretion System

Staphylococcus aureus employs the type VIIb secretion system (T7SSb) to secrete effector proteins that either have antibacterial activities or promote bacterial persistence in mouse infection models. Here, we present the crystal structure of the ATPase domain D3 of the EssC coupling protein from S. aureus USA300_FPR3757, an integral component of the T7SSb complex, resolved at a 1.7-Å resolution. EssC-D3 shares structural homology with FtsK/SpoIII-like ATPase domains of T7SSa and T7SSb and exhibits a conserved pocket on the surface with differential amino acid composition. In T7SSa, substrate EsxB interacts with the D3 domain through this pocket. Here, we identify amino acids in this pocket that are essential for effector protein secretion in the T7SSb. Our results reveal that the adjacent ATPase domain D2 is a substrate binding site on EssC and that substrates bound to D2 require domain D3 for further transport. Point mutations in the Walker B motif of domain D3 have diametric effects on secretion activity, either abolishing or boosting it, pointing to a critical role of domain D3 in the substrate transport. Finally, we identify ATPase domain D3 as a virulence determinant of S. aureus USA300_FPR3757 using an invertebrate in vivo infection model.IMPORTANCE The emergence of antibiotic-resistant bacteria poses a rising problem in antibiotic treatment (S. Boyle-Vavra and R. S. Daum, Lab Invest 87:3-9, 2007, https://doi.org/10.1038/labinvest.3700501). We have used the multidrug-resistant S. aureus USA300_FPR3757 as a model organism to study the T7SSb. Effector proteins of this system have been associated with abscess formation and bacterial persistence in mouse models (M. L. Burts, A. C. DeDent, and D. M. Missiakas, Mol Microbiol 69:736-746, 2008, https://doi.org/10.1111/j.1365-2958.2008.06324.x; M. L. Burts, W. A. Williams, K. DeBord, and D. M. Missiakas, Proc Natl Acad Sci U S A 102:1169-1174, 2005, https://doi.org/10.1073/pnas.0405620102). We determined the structure of the essential ATPase domain D3 of the T7SSb at atomic resolution and validated a surface-exposed pocket as a potential drug target to block secretion. Furthermore, our study provides new mechanistic insights into the T7SSb substrate transport.

Novel Immunoglobulin Domain Proteins Provide Insights into Evolution and Pathogenesis Mechanisms of SARS-Related Coronaviruses

A novel coronavirus (SARS-CoV-2) is the causative agent of an emergent severe respiratory disease (COVID-19) in humans that is threatening to result in a global health crisis. By using genomic, sequence, structural and evolutionary analysis, we show that Alpha- and Beta-CoVs possess several novel families of immunoglobulin (Ig) domain proteins, including ORF8 and ORF7a from SARS-related coronaviruses and two protein groups from certain Alpha-CoVs. Among them, ORF8 is distinguished in being rapidly evolving, possessing a unique insert and a hypervariable position among SARS-CoV-2 genomes in its predicted ligand-binding groove. We also uncover many Ig proteins from several metazoan viruses which are distinct in sequence and structure but share an architecture comparable to that of CoV Ig domain proteins. Hence, we propose that deployment of Ig domain proteins is a widely-used strategy by viruses, and SARS-CoV-2 ORF8 is a potential pathogenicity factor which evolves rapidly to counter the immune response and facilitate the transmission between hosts.

Highly regulated, diversifying NTP-dependent biological conflict systems with implications for the emergence of multicellularity

Social cellular aggregation or multicellular organization pose increased risk of transmission of infections through the system upon infection of a single cell. The generality of the evolutionary responses to this outside of Metazoa remains unclear. We report the discovery of several thematically unified, remarkable biological conflict systems preponderantly present in multicellular prokaryotes. These combine thresholding mechanisms utilizing NTPase chaperones (the MoxR-vWA couple), GTPases and proteolytic cascades with hypervariable effectors, which vary either by using a reverse transcriptase-dependent diversity-generating system or through a system of acquisition of diverse protein modules, typically in inactive form, from various cellular subsystems. Conciliant lines of evidence indicate their deployment against invasive entities, like viruses, to limit their spread in multicellular/social contexts via physical containment, dominant-negative interactions or apoptosis. These findings argue for both a similar operational 'grammar' and shared protein domains in the sensing and limiting of infections during the multiple emergences of multicellularity.

Clonorchis sinensis MF6p/HDM (CsMF6p/HDM) induces pro-inflammatory immune response in RAW 264.7 macrophage cells via NF-κB-dependent MAPK pathways

Background

MF6p/host defense molecules (HDMs) are a broad family of small proteins secreted by helminth parasites. Although the physiological role of MF6p/HDMs in trematode parasites is not fully understood, their potential biological function in maintaining heme homeostasis and modulating host immune response has been proposed.

Methods

A gene encoding the MF6p/HDM of Clonorchis sinensis (CsMF6p/HDM) was cloned. Recombinant CsMF6p/HDM (rCsMF6p/HDM) was expressed in Escherichia coli. The biochemical and immunological properties of rCsMF6/HDM were analyzed. CsMF6p/HDM induced pro-inflammatory response in RAW 264.7 cells was analyzed by cytokine array assay, reverse transcription polymerase chain reaction, and enzyme-linked immunosorbent assay. The structural feature of CsMF6p/HDM was analyzed by three-dimensional modeling and molecular docking simulations.

Results

The CsMF6p/HDM shares a high level of amino acid sequence similarity with orthologs from other trematodes and is expressed in diverse developmental stages of the parasite. The rCsMF6p/HDM bound to bacteria-derived lipopolysaccharide (LPS), without effectively neutralizing LPS-induced inflammatory response in RAW 264.7 macrophage cells. Rather, the rCsMF6p/HDM induced pro-inflammatory immune response, which is characterized by the expression of TNF-α and IL-6, in RAW 264.7 cells. The rCsMF6p/HDM-induced pro-inflammatory immune response was regulated by JNK and p38 MAPKs, and was effectively down-regulated via inhibition of NF-κB. The structural analysis of CsMF6p/HDM and the docking simulation with LPS suggested insufficient capture of LPS by CsMF6p/HDM, which suggested that rCsMF6p/HDM could not effectively neutralize LPS-induced inflammatory response in RAW 264.7 cells.

Conclusions

Although rCsMF6p/HDM binds to LPS, the binding affinity may not be sufficient to maintain a stable complex of rCsMF6p/HDM and LPS. Moreover, the rCsMF6p/HDM-induced pro-inflammatory response is characterized by the release of IL-6 and TNF-α in RAW 264.7 macrophage cells. The pro-inflammatory response induced by rCsMF6p/HDM is mediated via NF-κB-dependent MAPK signaling pathway. These results collectively suggest that CsMF6p/HDM mediates C. sinensis-induced inflammation cascades that eventually lead to hepatobiliary diseases.

Using Dali for Protein Structure Comparison

The exponential growth in the number of newly solved protein structures makes correlating and classifying the data an important task. Distance matrix alignment (Dali) is used routinely by crystallographers worldwide to screen the database of known structures for similarity to newly determined structures. Dali is easily accessible through the web server ( http://ekhidna.biocenter.helsinki.fi/dali ). Alternatively, the program may be downloaded and pairwise comparisons performed locally on Linux computers.

DALI and the persistence of protein shape

DALI is a popular resource for comparing protein structures. The software is based on distance-matrix alignment. The associated web server provides tools to navigate, integrate and organize some data pushed out by genomics and structural genomics. The server has been running continuously for the past 25 years. Structural biologists routinely use DALI to compare a new structure against previously known protein structures. If significant similarities are discovered, it may indicate a distant homology, that is, that the structures are of shared origin. This may be significant in determining the molecular mechanisms, as these may remain very similar from a distant predecessor to the present day, for example, from the last common ancestor of humans and bacteria. Meta-analysis of independent reference-based evaluations of alignment accuracy and fold discrimination shows DALI at top rank in six out of 12 studies. The web server and standalone software are available from http://ekhidna2.biocenter.helsinki.fi/dali.

Crystal structure of the translation recovery factor Trf from Sulfolobus solfataricus

During translation initiation, the heterotrimeric archaeal translation initiation factor 2 (aIF2) recruits the initiator tRNA_i to the small ribosomal subunit. In the stationary growth phase and/or during nutrient stress, Sulfolobus solfataricus aIF2 has a second function: It protects leaderless mRNAs against degradation by binding to their 5'-ends. The S. solfataricus protein Sso2509 is a translation recovery factor (Trf) that interacts with aIF2 and is responsible for the release of aIF2 from bound mRNAs, thereby enabling translation re-initiation. It is a member of the domain of unknown function 35 (DUF35) protein family and is conserved in Sulfolobales as well as in other archaea. Here, we present the X-ray structure of S. solfataricus Trf solved to a resolution of 1.65 Å. Trf is composed of an N-terminal rubredoxin-like domain containing a bound zinc ion and a C-terminal oligosaccharide/oligonucleotide binding fold domain. The Trf structure reveals putative mRNA binding sites in both domains. Surprisingly, the Trf protein is structurally but not sequentially very similar to proteins linked to acyl-CoA utilization-for example, the Sso2064 protein from S. solfataricus-as well as to scaffold proteins found in the acetoacetyl-CoA thiolase/high-mobility group-CoA synthase complex of the archaeon Methanothermococcus thermolithotrophicus and in a steroid side-chain-cleaving aldolase complex from the bacterium Thermomonospora curvata. This suggests that members of the DUF35 protein family are able to act as scaffolding and binding proteins in a wide variety of biological processes.

Binding site matching in rational drug design: algorithms and applications

Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.

The modular structure of α/β-hydrolases

The α/β-hydrolase fold family is highly diverse in sequence, structure and biochemical function. To investigate the sequence-structure-function relationships, the Lipase Engineering Database (https://led.biocatnet.de) was updated. Overall, 280 638 protein sequences and 1557 protein structures were analysed. All α/β-hydrolases consist of the catalytically active core domain, but they might also contain additional structural modules, resulting in 12 different architectures: core domain only, additional lids at three different positions, three different caps, additional N- or C-terminal domains and combinations of N- and C-terminal domains with caps and lids respectively. In addition, the α/β-hydrolases were distinguished by their oxyanion hole signature (GX-, GGGX- and Y-types). The N-terminal domains show two different folds, the Rossmann fold or the β-propeller fold. The C-terminal domains show a β-sandwich fold. The N-terminal β-propeller domain and the C-terminal β-sandwich domain are structurally similar to carbohydrate-binding proteins such as lectins. The classification was applied to the newly discovered polyethylene terephthalate (PET)-degrading PETases and MHETases, which are core domain α/β-hydrolases of the GX- and the GGGX-type respectively. To investigate evolutionary relationships, sequence networks were analysed. The degree distribution followed a power law with a scaling exponent γ = 1.4, indicating a highly inhomogeneous network which consists of a few hubs and a large number of less connected sequences. The hub sequences have many functional neighbours and therefore are expected to be robust toward possible deleterious effects of mutations. The cluster size distribution followed a power law with an extrapolated scaling exponent τ = 2.6, which strongly supports the connectedness of the sequence space of α/β-hydrolases. DATABASE: Supporting data about domains from other proteins with structural similarity to the N- or C-terminal domains of α/β-hydrolases are available in Data Repository of the University of Stuttgart (DaRUS) under doi: https://doi.org/10.18419/darus-458.

Structures of single-layer β-sheet proteins evolved from β-hairpin repeats

Free-standing single-layer β-sheets are extremely rare in naturally occurring proteins, even though β-sheet motifs are ubiquitous. Here we report the crystal structures of three homologous, single-layer, anti-parallel β-sheet proteins, comprised of three or four twisted β-hairpin repeats. The structures reveal that, in addition to the hydrogen bond network characteristic of β-sheets, additional hydrophobic interactions mediated by small clusters of residues adjacent to the turns likely play a significant role in the structural stability and compensate for the lack of a compact hydrophobic core. These structures enabled identification of a family of secreted proteins that are broadly distributed in bacteria from the human gut microbiome and are putatively involved in the metabolism of complex carbohydrates. A conserved surface patch, rich in solvent-exposed tyrosine residues, was identified on the concave surface of the β-sheet. These new modular single-layer β-sheet proteins may serve as a new model system for studying folding and design of β-rich proteins.

Structure and two-metal mechanism of fungal tRNA ligase

Fungal tRNA ligase (Trl1) is an essential enzyme that repairs RNA breaks with 2′,3′-cyclic-PO₄ and 5′-OH ends inflicted during tRNA splicing and non-canonical mRNA splicing in the fungal unfolded protein response. Trl1 is composed of C-terminal cyclic phosphodiesterase (CPD) and central GTP-dependent polynucleotide kinase (KIN) domains that heal the broken ends to generate the 3′-OH,2′-PO₄ and 5′-PO₄ termini required for sealing by an N-terminal ATP-dependent ligase domain (LIG). Here we report crystal structures of the Trl1-LIG domain from Chaetomium thermophilum at two discrete steps along the reaction pathway: the covalent LIG-(lysyl-Nζ)–AMP•Mn²⁺ intermediate and a LIG•ATP•(Mn²⁺)₂ Michaelis complex. The structures highlight a two-metal mechanism whereby a penta-hydrated metal complex stabilizes the transition state of the ATP α phosphate and a second metal bridges the β and γ phosphates to help orient the pyrophosphate leaving group. A LIG-bound sulfate anion is a plausible mimetic of the essential RNA terminal 2′-PO₄. Trl1-LIG has a distinctive C-terminal domain that instates fungal Trl1 as the founder of an Rnl6 clade of ATP-dependent RNA ligase. We discuss how the Trl1-LIG structure rationalizes the large body of in vivo structure–function data for Saccharomyces cerevisiae Trl1.

In silico evidence of de novo interactions between ribosomal and Epstein - Barr virus proteins

Background

Association of Epstein-Barr virus (EBV) encoded latent gene products with host ribosomal proteins (RPs) has not been fully explored, despite their involvement in the aetiology of several human cancers. To gain an insight into their plausible interactions, we employed a computational approach that encompasses structural alignment, gene ontology analysis, pathway analysis, and molecular docking.

Results

In this study, the alignment analysis based on structural similarity allows the prediction of 48 potential interactions between 27 human RPs and the EBV proteins EBNA1, LMP1, LMP2A, and LMP2B. Gene ontology analysis of the putative protein-protein interactions (PPIs) reveals their probable involvement in RNA binding, ribosome biogenesis, metabolic and biosynthetic processes, and gene regulation. Pathway analysis shows their possible participation in viral infection strategies (viral translation), as well as oncogenesis (Wnt and EGFR signalling pathways). Finally, our molecular docking assay predicts the functional interactions of EBNA1 with four RPs individually: EBNA1-eS10, EBNA1-eS25, EBNA1-uL10 and EBNA1-uL11.

Conclusion

These interactions have never been revealed previously via either experimental or in silico approach. We envisage that the calculated interactions between the ribosomal and EBV proteins herein would provide a hypothetical model for future experimental studies on the functional relationship between ribosomal proteins and EBV infection.

Electronic supplementary material

The online version of this article (10.1186/s12860-019-0219-y) contains supplementary material, which is available to authorized users.

Estimating statistical significance of local protein profile-profile alignments

BACKGROUND

Alignment of sequence families described by profiles provides a sensitive means for establishing homology between proteins and is important in protein evolutionary, structural, and functional studies. In the context of a steadily growing amount of sequence data, estimating the statistical significance of alignments, including profile-profile alignments, plays a key role in alignment-based homology search algorithms. Still, it is an open question as to what and whether one type of distribution governs profile-profile alignment score, especially when profile-profile substitution scores involve such terms as secondary structure predictions.

RESULTS

This study presents a methodology for estimating the statistical significance of this type of alignments. The methodology rests on a new algorithm developed for generating random profiles such that their alignment scores are distributed similarly to those obtained for real unrelated profiles. We show that improvements in statistical accuracy and sensitivity and high-quality alignment rate result from statistically characterizing alignments by establishing the dependence of statistical parameters on various measures associated with both individual and pairwise profile characteristics. Implemented in the COMER software, the proposed methodology yielded an increase of up to 34.2% in the number of true positives and up to 61.8% in the number of high-quality alignments with respect to the previous version of the COMER method.

CONCLUSIONS

The more accurate estimation of statistical significance is implemented in the COMER method, which is now more sensitive and provides an increased rate of high-quality profile-profile alignments. The results of the present study also suggest directions for future research.

mTM-align: a server for fast protein structure database search and multiple protein structure alignment

With the rapid increase of the number of protein structures in the Protein Data Bank, it becomes urgent to develop algorithms for efficient protein structure comparisons. In this article, we present the mTM-align server, which consists of two closely related modules: one for structure database search and the other for multiple structure alignment. The database search is speeded up based on a heuristic algorithm and a hierarchical organization of the structures in the database. The multiple structure alignment is performed using the recently developed algorithm mTM-align. Benchmark tests demonstrate that our algorithms outperform other peering methods for both modules, in terms of speed and accuracy. One of the unique features for the server is the interplay between database search and multiple structure alignment. The server provides service not only for performing fast database search, but also for making accurate multiple structure alignment with the structures found by the search. For the database search, it takes about 2-5 min for a structure of a medium size (∼300 residues). For the multiple structure alignment, it takes a few seconds for ∼10 structures of medium sizes. The server is freely available at: http://yanglab.nankai.edu.cn/mTM-align/.