The ProFunc Function Prediction Server

BeEM: fast and faithful conversion of mmCIF format structure files to PDB format

BACKGROUND

Although mmCIF is the current official format for deposition of protein and nucleic acid structures to the protein data bank (PDB) database, the legacy PDB format is still the primary supported format for many structural bioinformatics tools. Therefore, reliable software to convert mmCIF structure files to PDB files is needed. Unfortunately, existing conversion programs fail to correctly convert many mmCIF files, especially those with many atoms and/or long chain identifies.

RESULTS

This study proposed BeEM, which converts any mmCIF format structure files to PDB format. BeEM conversion faithfully retains all atomic and chain information, including chain IDs with more than 2 characters, which are not supported by any existing mmCIF to PDB converters. The conversion speed of BeEM is at least ten times faster than existing converters such as MAXIT and Phenix. Part of the reason for the speed improvement is the avoidance of conversion between numerical values and text strings.

CONCLUSION

BeEM is a fast and accurate tool for mmCIF-to-PDB format conversion, which is a common procedure in structural biology. The source code is available under the BSD licence at https://github.com/kad-ecoli/BeEM/ .

BESFA: bioinformatics based evolutionary, structural & functional analysis of prostrate, Placenta, Ovary, Testis, and Embryo (POTE) paralogs

The POTE family comprises 14 paralogues and is primarily expressed in Prostrate, Placenta, Ovary, Testis, Embryo (POTE), and cancerous cells. The prospective function of the POTE protein family under physiological conditions is less understood. We systematically analyzed their cellular localization and molecular docking analysis to elucidate POTE proteins' structure, function, and Adaptive Divergence. Our results suggest that group three POTE paralogs (POTEE, POTEF, POTEI, POTEJ, and POTEKP (a pseudogene)) exhibits significant variation among other members could be because of their Adaptive Divergence. Furthermore, our molecular docking studies on POTE protein revealed the highest binding affinity with NCI-approved anticancer compounds. Additionally, POTEE, POTEF, POTEI, and POTEJ were subject to an explicit molecular dynamic simulation for 50ns. MM-GBSA and other essential electrostatics were calculated that showcased that only POTEE and POTEF have absolute binding affinities with minimum energy exploitation. Thus, this study’s outcomes are expected to drive cancer research to successful utilization of POTE genes family as a new biomarker, which could pave the way for the discovery of new therapies.

Substrate size-dependent conformational changes of bacterial pectin-binding protein crucial for chemotaxis and assimilation

Gram-negative Sphingomonas sp. strain A1 exhibits positive chemotaxis toward acidic polysaccharide pectin. SPH1118 has been identified as a pectin-binding protein involved in both pectin chemotaxis and assimilation. Here we show tertiary structures of SPH1118 with six different conformations as determined by X-ray crystallography. SPH1118 consisted of two domains with a large cleft between the domains and substrates bound to positively charged and aromatic residues in the cleft through hydrogen bond and stacking interactions. Substrate-free SPH1118 adopted three different conformations in the open form. On the other hand, the two domains were closed in substrate-bound form and the domain closure ratio was changed in response to the substrate size, suggesting that the conformational change upon binding to the substrate triggered the expression of pectin chemotaxis and assimilation. This study first clarified that the solute-binding protein with dual functions recognized the substrate through flexible conformational changes in response to the substrate size.

Structural and functional characterization of β-cyanoalanine synthase from Tetranychus urticae

Tetranychus urticae is a polyphagous spider mite that can feed on more than 1100 plant species including cyanogenic plants. The herbivore genome contains a horizontally acquired gene tetur10g01570 (TuCAS) that was previously shown to participate in cyanide detoxification. To understand the structure and determine the function of TuCAS in T. urticae, crystal structures of the protein with lysine conjugated pyridoxal phosphate (PLP) were determined. These structures reveal extensive TuCAS homology with the β-substituted alanine synthase family, and they show that this enzyme utilizes a similar chemical mechanism involving a stable α-aminoacrylate intermediate in β-cyanoalanine and cysteine synthesis. We demonstrate that TuCAS is more efficient in the synthesis of β-cyanoalanine, which is a product of the detoxification reaction between cysteine and cyanide, than in the biosynthesis of cysteine. Also, the enzyme carries additional enzymatic activities that were not previously described. We show that TuCAS can detoxify cyanide using O-acetyl-L-serine as a substrate, leading to the direct formation of β-cyanoalanine. Moreover, it catalyzes the reaction between the TuCAS-bound α-aminoacrylate intermediate and aromatic compounds with a thiol group. In addition, we have tested several compounds as TuCAS inhibitors. Overall, this study identifies additional functions for TuCAS and provides new molecular insight into the xenobiotic metabolism of T. urticae.

Mycobacterium tuberculosis puromycin hydrolase displays a prolyl oligopeptidase fold and an acyl aminopeptidase activity

Puromycin-hydrolizing peptidases have been described as members of the prolyl oligopeptidase peptidase family. These enzymes are present across all domains of life but still little is known of the homologs found in the pathogenic bacterium Mycobacterium tuberculosis. The crystal structure of a M. tuberculosis puromycin hydrolase peptidase has been determined at 3 Angstrom resolution, revealing a conserved prolyl oligopeptidase fold, defined by α/β-hydrolase and β-propeller domains with two distinctive loops that occlude access of large substrates to the active site. The enzyme displayed amino peptidase activity with a substrate specificity preference for hydrophobic residues in the decreasing order of phenylalanine, leucine, alanine and proline. The enzyme's active site is lined by residues Glu564 for the coordination of the substrates amino terminal moiety and His561, Val608, Tyr78, Trp306, Phe563 and Ty567 for the accommodation of hydrophobic substrates. The availability of a crystal structure for puromycin hydrolase of M. tuberculosis shall facilitate the development of inhibitors with therapeutic applications.

Cryo-EM structures reveal intricate Fe-S cluster arrangement and charging in Rhodobacter capsulatus formate dehydrogenase

Metal-containing formate dehydrogenases (FDH) catalyse the reversible oxidation of formate to carbon dioxide at their molybdenum or tungsten active site. They display a diverse subunit and cofactor composition, but structural information on these enzymes is limited. Here we report the cryo-electron microscopic structures of the soluble Rhodobacter capsulatus FDH (RcFDH) as isolated and in the presence of reduced nicotinamide adenine dinucleotide (NADH). RcFDH assembles into a 360 kDa dimer of heterotetramers revealing a putative interconnection of electron pathway chains. In the presence of NADH, the RcFDH structure shows charging of cofactors, indicative of an increased electron load.

Rhodobacter capsulatus NAD⁺ dependent formate dehydrogenase (RcFDH) is a molybdoenzyme that catalyses the reversible oxidation of formate to carbon dioxide, and is of interest for biotechnological applications. Here the authors present the cryo-EM structures of RcFDH as isolated from R. capsulatus and in the reduced state with bound NADH, and discuss the enzyme mechanism.

Binding Options for the Small Subunit-Like Domain of Cyanobacteria to Rubisco

Two proteins found in cyanobacteria contain a C-terminal domain with homology to the small subunit of rubisco (RbcS). These small subunit-like domains (SSLDs) are important features of CcmM, a protein involved in the biogenesis of carboxysomes found in all β-cyanobacteria, and a rubisco activase homolog [activase-like protein of cyanobacteria (ALC)] found in over a third of sequenced cyanobacterial genomes. Interaction with rubisco is crucial to the function of CcmM and is believed to be important to ALC as well. In both cases, the SSLD aggregates rubisco, and this nucleation event may be important in regulating rubisco assembly and activity. Recently, two independent studies supported the conclusion that the SSLD of CcmM binds equatorially to L₈S₈ holoenzymes of rubisco rather than by displacing an RbcS, as its structural homology would suggest. We use sequence analysis and homology modeling to examine whether the SSLD from the ALC could bind the large subunit of rubisco either via an equatorial interaction or in an RbcS site, if available. We suggest that the SSLD from the ALC of Fremyella diplosiphon could bind either in a vacant RbcS site or equatorially. Our homology modeling takes into account N-terminal residues not represented in available cryo-electron microscopy structures that potentially contribute to the interface between the large subunit of rubisco (RbcL) and RbcS. Here, we suggest the perspective that binding site variability as a means of regulation is plausible and that the dynamic interaction between the RbcL, RbcS, and SSLDs may be important for carboxysome assembly and function.

A bipartite structural organization defines the SERINC family of HIV-1 restriction factors

The human integral membrane protein SERINC5 potently restricts HIV-1 infectivity and sensitizes the virus to antibody-mediated neutralization. Here, using cryo-EM, we determine the structures of human SERINC5 and its orthologue from Drosophila melanogaster at subnanometer and near-atomic resolution, respectively. The structures reveal a novel fold comprised of ten transmembrane helices organized into two subdomains and bisected by a long diagonal helix. A lipid binding groove and clusters of conserved residues highlight potential functional sites. A structure-based mutagenesis scan identified surface-exposed regions and the interface between the subdomains of SERINC5 as critical for HIV-1-restriction activity. The same regions are also important for viral sensitization to neutralizing antibodies, directly linking the antiviral activity of SERINC5 with remodeling of the HIV-1 envelope glycoprotein.

Crystal structure of a hypothetical T2SS effector Lpg0189 from Legionella pneumophila reveals a novel protein fold

Lpg0189 is a type II secretion system-dependent extracellular protein with unknown function from Legionella pneumophila. Herein, we determined the crystal structure of Lpg0189 at 1.98 Å resolution by using single-wavelength anomalous diffraction (SAD). Lpg0189 folds into a novel chair-shaped architecture, with two sheets roughly perpendicular to each other. Bioinformatics analysis suggests Lpg0189 and its homologues are unique to Legionellales and evolved divergently. The interlinking structural and bioinformatics studies provide a better understanding of this hypothetical protein.

Structural and functional characterization of an intradiol ring-cleavage dioxygenase from the polyphagous spider mite herbivore Tetranychus urticae Koch

Genome analyses of the polyphagous spider mite herbivore Tetranychus urticae (two-spotted spider mite) revealed the presence of a set of 17 genes that code for secreted proteins belonging to the "intradiol dioxygenase-like" subgroup. Phylogenetic analyses indicate that this novel enzyme family has been acquired by horizontal gene transfer. In order to better understand the role of these proteins in T. urticae, we have structurally and functionally characterized one paralog (tetur07g02040). It was demonstrated that this protein is indeed an intradiol ring-cleavage dioxygenase, as the enzyme is able to cleave catechol between two hydroxyl-groups using atmospheric dioxygen. The enzyme was characterized functionally and structurally. The active site of the T. urticae enzyme contains an Fe³⁺ cofactor that is coordinated by two histidine and two tyrosine residues, an arrangement that is similar to those observed in bacterial homologs. However, the active site is significantly more solvent exposed than in bacterial proteins. Moreover, the mite enzyme is monomeric, while almost all structurally characterized bacterial homologs form oligomeric assemblies. Tetur07g02040 is not only the first spider mite dioxygenase that has been characterized at the molecular level, but is also the first structurally characterized intradiol ring-cleavage dioxygenase originating from a eukaryote.

Function Prediction for G Protein-Coupled Receptors through Text Mining and Induction Matrix Completion

G protein-coupled receptors (GPCRs) constitute the key component of cellular signal transduction. Accurately annotating the biological functions of GPCR proteins is vital to the understanding of the physiological processes they involve in. With the rapid development of text mining technologies and the exponential growth of biomedical literature, it becomes urgent to explore biological functional information from various literature for systematically and reliably annotating these known GPCRs. We design a novel three-stage approach, TM-IMC, using text mining and inductive matrix completion, for automated prediction of the gene ontology (GO) terms of the GPCR proteins. Large-scale benchmark tests show that inductive matrix completion models contribute to GPCR-GO association prediction for both molecular function and biological process aspects. Moreover, our detailed data analysis shows that information extracted from GPCR-associated literature indeed contributes to the prediction of GPCR-GO associations. The study demonstrated a new avenue to enhance the accuracy of GPCR function annotation through the combination of text mining and induction matrix completion over baseline methods in critical assessment of protein function annotation algorithms and literature-based GO annotation methods. Source codes of TM-IMC and the involved datasets can be freely downloaded from https://zhanglab.ccmb.med.umich.edu/TM-IMC for academic purposes.

Genome-wide transcriptional analyses in Anopheles mosquitoes reveal an unexpected association between salivary gland gene expression and insecticide resistance

BACKGROUND

To combat malaria transmission, the Ugandan government has embarked upon an ambitious programme of indoor residual spraying (IRS) with a carbamate class insecticide, bendiocarb. In preparation for this campaign, we characterized bendiocarb resistance and associated transcriptional variation among Anopheles gambiae s.s. mosquitoes from two sites in Uganda.

RESULTS

Gene expression in two mosquito populations displaying some resistance to bendiocarb (95% and 79% An. gambiae s.l. WHO tube bioassay mortality in Nagongera and Kihihi, respectively) was investigated using whole-genome microarrays. Significant overexpression of several genes encoding salivary gland proteins, including D7r2 and D7r4, was detected in mosquitoes from Nagongera. In Kihihi, D7r4, two detoxification-associated genes (Cyp6m2 and Gstd3) and an epithelial serine protease were among the genes most highly overexpressed in resistant mosquitoes. Following the first round of IRS in Nagongera, bendiocarb-resistant mosquitoes were collected, and real-time quantitative PCR analyses detected significant overexpression of D7r2 and D7r4 in resistant mosquitoes. A single nucleotide polymorphism located in a non-coding transcript downstream of the D7 genes was found at a significantly higher frequency in resistant individuals. In silico modelling of the interaction between D7r4 and bendiocarb demonstrated similarity between the insecticide and serotonin, a known ligand of D7 proteins. A meta-analysis of published microarray studies revealed a recurring association between D7 expression and insecticide resistance across Anopheles species and locations.

CONCLUSIONS

A whole-genome microarray approach identified an association between novel insecticide resistance candidates and bendiocarb resistance in Uganda. In addition, a single nucleotide polymorphism associated with this resistance mechanism was discovered. The use of such impartial screening methods allows for discovery of resistance candidates that have no previously-ascribed function in insecticide binding or detoxification. Characterizing these novel candidates will broaden our understanding of resistance mechanisms and yield new strategies for combatting widespread insecticide resistance among malaria vectors.

Heligmosomoides polygyrus Venom Allergen-like Protein-4 (HpVAL-4) is a sterol binding protein

•

Heligmosomoides polygyrus Venom Allergen-like Protein-4 (HpVAL-4) was produced in plants as a glycosylated protein.

•

The crystal structure of HpVAL-4 was solved and reveals three distinct cavities.

•

These cavities are the central cavity; the sterol-binding caveolin-binding motif (CBM); and the palmitate-binding cavity.

•

The central cavity of Hp-VAL-4 lacks the characteristic histidines that coordinate divalent cations.

•

Hp-VAL-4 binds sterol in vivo and in vitro.

Heligmosomoides polygyrus bakeri is a model parasitic hookworm used to study animal and human helminth diseases. During infection, the parasite releases excretory/secretory products that modulate the immune system of the host. The most abundant protein family in excretory/secretory products comprises the venom allergen-like proteins (VALs), which are members of the SCP/TAPS (sperm-coating protein/Tpx/antigen 5/pathogenesis related-1/Sc7) superfamily. There are >30 secreted Heligmosomoides polygyrus VAL proteins (HpVALs) and these proteins are characterised by having either one or two 15 kDa CAP (cysteine-rich secretory protein (CRISP)/antigen 5/pathogenesis related-1) domains. The first known HpVAL structure, HpVAL-4, refined to 1.9 Å is reported. HpVAL-4 was produced as a homogeneously glycosylated protein in leaves of Nicotiana benthamiana infiltrated with recombinant plasmids, making this plant expression platform amenable for the production of biological products. The overall topology of HpVAL-4 is a three layered αβα sandwich between a short N-terminal loop and a C-terminal cysteine rich extension. The C-terminal cysteine rich extension has two strands stabilized by two disulfide bonds and superposes well with the previously reported extension from the human hookworm Necator americanus Ancylostoma secreted protein-2 (Na-ASP-2). The N-terminal loop is connected to alpha helix 2 via a disulfide bond previously observed in Na-ASP-2. HpVAL-4 has a central cavity that is more similar to the N-terminal CAP domain of the two CAP Na-ASP-1 from Necator americanus. Unlike Na-ASP-2, mammalian CRISP, and the C-terminal CAP domain of Na-ASP-1, the large central cavity of HpVAL-4 lacks the two histidines required to coordinate divalent cations. HpVAL-4 has both palmitate-binding and sterol-binding cavities and is able to complement the in vivo sterol export phenotype of yeast mutants lacking their endogenous CAP proteins. More studies are required to determine endogenous binding partners of HpVAL-4 and unravel the possible impact of sterol binding on immune-modulatory functions.

Crystal Structure of Borrelia turicatae protein, BTA121, a differentially regulated gene in the tick-mammalian transmission cycle of relapsing fever spirochetes

Tick-borne relapsing fever (RF) borreliosis is a neglected disease that is often misdiagnosed. RF species circulating in the United States include Borrelia turicatae, which is transmitted by argasid ticks. Environmental adaptation by RF Borrelia is poorly understood, however our previous studies indicated differential regulation of B. turicatae genes localized on the 150 kb linear megaplasmid during the tick-mammalian transmission cycle, including bta121. This gene is up-regulated by B. turicatae in the tick versus the mammal, and the encoded protein (BTA121) is predicted to be surface localized. The structure of BTA121 was solved by single-wavelength anomalous dispersion (SAD) using selenomethionine-derivative protein. The topology of BTA121 is unique with four helical domains organized into two helical bundles. Due to the sequence similarity of several genes on the megaplasmid, BTA121 can serve as a model for their tertiary  structures. BTA121 has large interconnected tunnels and cavities that can accommodate ligands, notably long parallel helices, which have a large hydrophobic central pocket. Preliminary in-vitro studies suggest that BTA121 binds lipids, notably palmitate with a similar order of binding affinity as tablysin-15, a known palmitate-binding protein. The reported data will guide mechanistic studies to determine the role of BTA121 in the tick-mammalian transmission cycle of B. turicatae.