Multiple Alignment of Protein Structures in Three Dimensions

Structural and Functional Characterization of a Gene Cluster Responsible for Deglycosylation of C-glucosyl Flavonoids and Xanthonoids by Deinococcus aerius

Plant C-glycosylated aromatic polyketides are important for plant and animal health. These are specialized metabolites that perform functions both within the plant, and in interaction with soil or intestinal microbes. Despite the importance of these plant compounds, there is still limited knowledge of how they are metabolized. The Gram-positive aerobic soil bacterium Deinococcus aerius strain TR0125 and other Deinococcus species thrive in a wide range of harsh environments. In this work, we identified a C-glycoside deglycosylation gene cluster in the genome of D. aerius. The cluster includes three genes coding for a GMC-type oxidoreductase (DaCGO1) that oxidizes the glucosyl C3 position in aromatic C-glucosyl compounds, which in turn provides the substrate for the C-glycoside deglycosidase (DaCGD; composed of α+β subunits) that cleaves the glucosyl-aglycone C-C bond. Our results from size-exclusion chromatography, single particle cryo-electron microscopy and X-ray crystallography show that DaCGD is an α2β2 heterotetramer, which represents a novel oligomeric state among bacterial CGDs. Importantly, the high-resolution X-ray structure of DaCGD provides valuable insights into the activation of the catalytic hydroxide ion by Lys261. DaCGO1 is specific for the 6-C-glucosyl flavones isovitexin, isoorientin and the 2-C-glucosyl xanthonoid mangiferin, and the subsequent C-C-bond cleavage by DaCGD generated apigenin, luteolin and norathyriol, respectively. Of the substrates tested, isovitexin was the preferred substrate (DaCGO1, Km 0.047 mM, kcat 51 min-1; DaCGO1/DaCGD, Km 0.083 mM, kcat 0.42 min-1).

Structural and Biochemical Characterization of MppQ, an L-Enduracididine Biosynthetic Enzyme from Streptomyces hygroscopicus

MppQ is an enzyme of unknown function from Streptomyces hygroscopicus (ShMppQ) that operates in the biosynthesis of the nonproteinogenic amino acid L-enduracididine (L-End). Since L-End is a component of several peptides showing activity against antibiotic-resistant pathogens, understanding its biosynthetic pathway could facilitate the development of chemoenzymatic routes to novel antibiotics. Herein, we report on the crystal structures of ShMppQ complexed with pyridoxal-5'-phosphate (PLP) and pyridoxamine-5'-phosphate (PMP). ShMppQ is similar to fold-type I PLP-dependent aminotransferases like aspartate aminotransferase. The tertiary structure of ShMppQ is composed of an N-terminal extension, a large domain, and a small domain. The active site is placed at the junction of the large and small domains and includes residues from both protomers of the homodimer. We also report the first functional characterization of MppQ, which we incubated with the enzymatically produced 2-ketoenduracidine and observed the conversion to L-End, establishing ShMppQ as the final enzyme in L-End biosynthesis. Additionally, we have observed that MppQ has a relatively high affinity for 2-keto-5-guanidinovaleric acid (i.e., 2-ketoarginine), a shunt product of MppP, indicating the potential role of MppQ in increasing the efficiency of L-End biosynthesis by converting 2-ketoarginine back to the starting material, l-arginine. A panel of potential amino-donor substrates was tested for the transamination activity against a saturating concentration of 2-ketoarginine in end-point assays. Most l-Arg was produced with l-ornithine as the donor substrate. Steady-state kinetic analysis of the transamination reaction with l-Orn and 2-ketoarginine shows that the kinetic constants are in line with those for the amino donor substrate of other fold-type I aminotransferases.

Functional annotation of haloacid dehalogenase superfamily structural genomics proteins

Haloacid dehalogenases (HAD) are members of a large superfamily that includes many Structural Genomics proteins with poorly characterized functionality. This superfamily consists of multiple types of enzymes that can act as sugar phosphatases, haloacid dehalogenases, phosphonoacetaldehyde hydrolases, ATPases, or phosphate monoesterases. Here, we report on predicted functional annotations and experimental testing by direct biochemical assay for Structural Genomics proteins from the HAD superfamily. To characterize the functions of HAD superfamily members, nine representative HAD proteins and 21 structural genomics proteins are analyzed. Using techniques based on computed chemical and electrostatic properties of individual amino acids, the functions of five structural genomics proteins from the HAD superfamily are predicted and validated by biochemical assays. A dehalogenase-like hydrolase, RSc1362 (Uniprot Q8XZN3, PDB 3UMB) is predicted to be a dehalogenase and dehalogenase activity is confirmed experimentally. Four proteins predicted to be sugar phosphatases are characterized as follows: a sugar phosphatase from Thermophilus volcanium (Uniprot Q978Y6) with trehalose-6-phosphate phosphatase and fructose-6-phosphate phosphatase activity; haloacid dehalogenase-like hydrolase from Bacteroides thetaiotaomicron (Uniprot Q8A2F3; PDB 3NIW) with fructose-6-phosphate phosphatase and sucrose-6-phosphate phosphatase activity; putative phosphatase from Eubacterium rectale (Uniprot D0VWU2; PDB 3DAO) as a sucrose-6-phosphate phosphatase; and hypothetical protein from Geobacillus kaustophilus (Uniprot Q5L139; PDB 2PQ0) as a fructose-6-phosphate phosphatase. Most of these sugar phosphatases showed some substrate promiscuity.

A previously unrecognized superfamily of macro-conotoxins includes an inhibitor of the sensory neuron calcium channel Cav2.3

Animal venom peptides represent valuable compounds for biomedical exploration. The venoms of marine cone snails constitute a particularly rich source of peptide toxins, known as conotoxins. Here, we identify the sequence of an unusually large conotoxin, Mu8.1, which defines a new class of conotoxins evolutionarily related to the well-known con-ikot-ikots and 2 additional conotoxin classes not previously described. The crystal structure of recombinant Mu8.1 displays a saposin-like fold and shows structural similarity with con-ikot-ikot. Functional studies demonstrate that Mu8.1 curtails calcium influx in defined classes of murine somatosensory dorsal root ganglion (DRG) neurons. When tested on a variety of recombinantly expressed voltage-gated ion channels, Mu8.1 displayed the highest potency against the R-type (Cav2.3) calcium channel. Ca2+ signals from Mu8.1-sensitive DRG neurons were also inhibited by SNX-482, a known spider peptide modulator of Cav2.3 and voltage-gated K+ (Kv4) channels. Our findings highlight the potential of Mu8.1 as a molecular tool to identify and study neuronal subclasses expressing Cav2.3. Importantly, this multidisciplinary study showcases the potential of uncovering novel structures and bioactivities within the largely unexplored group of macro-conotoxins.

Staphylococcal Periscope proteins Aap, SasG, and Pls project noncanonical legume-like lectin adhesin domains from the bacterial surface

Staphylococcus aureus and Staphylococcus epidermidis are frequently associated with medical device infections that involve establishment of a bacterial biofilm on the device surface. Staphylococcal surface proteins Aap, SasG, and Pls are members of the Periscope Protein class and have been implicated in biofilm formation and host colonization; they comprise a repetitive region ("B region") and an N-terminal host colonization domain within the "A region," predicted to be a lectin domain. Repetitive E-G5 domains (as found in Aap, SasG, and Pls) form elongated "stalks" that would vary in length with repeat number, resulting in projection of the N-terminal A domain variable distances from the bacterial cell surface. Here, we present the structures of the lectin domains within A regions of SasG, Aap, and Pls and a structure of the Aap lectin domain attached to contiguous E-G5 repeats, suggesting the lectin domains will sit at the tip of the variable length rod. We demonstrate that these isolated domains (Aap, SasG) are sufficient to bind to human host desquamated nasal epithelial cells. Previously, proteolytic cleavage or a deletion within the A domain had been reported to induce biofilm formation; the structures suggest a potential link between these observations. Intriguingly, while the Aap, SasG, and Pls lectin domains bind a metal ion, they lack the nonproline cis peptide bond thought to be key for carbohydrate binding by the lectin fold. This suggestion of noncanonical ligand binding should be a key consideration when investigating the host cell interactions of these bacterial surface proteins.

Between Protein Fold and Nucleophile Identity: Multiscale Modeling of the TEV Protease Enzyme-Substrate Complex

The cysteine protease from the tobacco etch virus (TEVp) is a well-known and widely utilized enzyme. TEVp's chymotrypsin-like fold is generally associated with serine catalytic triads that differ in terms of a reaction mechanism from the most well-studied papain-like cysteine proteases. The question of what dominates the TEVp mechanism, nucleophile identity, or structural composition has never been previously addressed. Here, we use enhanced sampling multiscale modeling to uncover that TEVp combines the features of two worlds in such a way that potentially hampers its activity. We show that TEVp cysteine is strictly in the anionic form in a free enzyme similar to papain. Peptide binding shifts the equilibrium toward the nucleophile's protonated form, characteristic of chymotrypsin-like proteases, although the cysteinyl anion form is still present and interconversion is rapid. This way cysteine protonation generates enzyme states that are a diversion from the most effective course of action, with only 13.2% of Michaelis complex sub-states able to initiate the reaction. As a result, we propose an updated view on the reaction mechanism catalyzed by TEVp. We also demonstrate that AlphaFold is able to construct protease-substrate complexes with high accuracy. We propose that our findings open a way for its industrious use in enzymological tasks. Unique features of TEVp discovered in this work open a discussion on the evolutionary history and trade-offs of optimizing serine triad-associated folds to cysteine as a nucleophile.

AlphaFold predicts the most complex protein knot and composite protein knots

The computer artificial intelligence system AlphaFold has recently predicted previously unknown three‐dimensional structures of thousands of proteins. Focusing on the subset with high‐confidence scores, we algorithmically analyze these predictions for cases where the protein backbone exhibits rare topological complexity, that is, knotting. Amongst others, we discovered a 7₁‐knot, the most topologically complex knot ever found in a protein, as well several six‐crossing composite knots comprised of two methyltransferase or carbonic anhydrase domains, each containing a simple trefoil knot. These deeply embedded composite knots occur evidently by gene duplication and interconnection of knotted dimers. Finally, we report two new five‐crossing knots including the first 5₁‐knot. Our list of analyzed structures forms the basis for future experimental studies to confirm these novel‐knotted topologies and to explore their complex folding mechanisms.

Functional Characterization of Structural Genomics Proteins in the Crotonase Superfamily

Members of the Crotonase superfamily, a mechanistically diverse family of proteins that share a conserved quaternary structure, can often catalyze more than one reaction. However, the spectrum of activity for its members has not been well studied. We report on measured crotonase and hydrolase activity for eight structural genomics (SG) proteins from the Crotonase superfamily plus two previously characterized proteins, intended as controls: human enoyl CoA hydratase (ECH) and Anabaena β-diketone hydrolase. Like most of the 15,000+ SG protein structures deposited in the Protein Data Bank (PDB), the eight SG proteins are of unknown or uncertain biochemical function. The functional characterization of the eight SG proteins is guided by the Structurally Aligned Local Sites of Activity (SALSA), a local-structure-based computational approach to functional annotation. For human ECH, the turnover number for hydrolase activity is threefold higher than that for ECH activity, although the catalytic efficiency is 160-fold higher for ECH. Three SG proteins originally annotated as ECHs were predicted by SALSA to be hydrolases and are observed to have higher catalytic efficiencies for hydrolase activity than for ECH activity, on par with the previously characterized hydrolase. Among the five SG proteins predicted by SALSA to be ECHs, all but one also show some hydrolase activity; all five exhibit lower ECH activity than the human ECH with respect to the crotonyl-CoA substrate. Here, we show examples demonstrating that SALSA can correct functional misannotations even within enzyme families that display promiscuous activity.

Predicted structural mimicry of spike receptor-binding motifs from highly pathogenic human coronaviruses

•

Potential coronavirus spike protein mimicry revealed by structural comparison.

•

Human and non-human protein potential interactions with virus identified.

•

Predicted structural mimicry corroborated by protein–protein docking.

•

Epitope-based alignments may help guide vaccine efforts.

Viruses often encode proteins that mimic host proteins in order to facilitate infection. Little work has been done to understand the potential mimicry of the SARS-CoV-2, SARS-CoV, and MERS-CoV spike proteins, particularly the receptor-binding motifs, which could be important in determining tropism and druggability of the virus. Peptide and epitope motifs have been detected on coronavirus spike proteins using sequence homology approaches; however, comparing the three-dimensional shape of the protein has been shown as more informative in predicting mimicry than sequence-based comparisons. Here, we use structural bioinformatics software to characterize potential mimicry of the three coronavirus spike protein receptor-binding motifs. We utilize sequence-independent alignment tools to compare structurally known protein models with the receptor-binding motifs and verify potential mimicked interactions with protein docking simulations. Both human and non-human proteins were returned for all three receptor-binding motifs. For example, all three were similar to several proteins containing EGF-like domains: some of which are endogenous to humans, such as thrombomodulin, and others exogenous, such as Plasmodium falciparum MSP-1. Similarity to human proteins may reveal which pathways the spike protein is co-opting, while analogous non-human proteins may indicate shared host interaction partners and overlapping antibody cross-reactivity. These findings can help guide experimental efforts to further understand potential interactions between human and coronavirus proteins.

Comparison of Hyaluronic Acid Biosynthetic Genes From Different Strains of Pasteurella multocida

Pasteurella multocida produces a capsule composed of different polysaccharides according to the capsular serotype (A, B, D, E, and F). Hyaluronic acid (HA) is a component of certain capsular types of this bacterium, especially capsular type A. Previously, 2 HA biosynthetic genes from a capsular type A strain were studied for the industrial-scale improvement of HA production. Molecular comparison of these genes across different capsular serotypes of P multocida has not been reported. This study aimed to compare 8 HA biosynthetic genes (pgi, pgm, galU, hyaC, glmS, glmM, glmU, and hyaD) of 22 P multocida strains (A:B:D:F = 6:6:6:4) with those of other organisms using sequence and structural bioinformatics analyses. These 8 genes showed a high level of within-species similarity (98%-99%) compared with other organisms. Only the last gene of 4 strains with capsular type F (HN07, PM70, HNF01, and HNF02) significantly differed from those of other strains (82%). Analysis of amino acid patterns together with phylogenetic results showed that the HA biosynthetic genes of the type A were closely related within the group. The genes in the capsular type F strain were notably similar to those of the capsular type A strain. Protein structural analysis supported structural similarities of the encoded enzymes between the strains of capsular types A, B, D, and F, except for the Pgm, GlmS, GlmU, and HyaD proteins. Our bioinformatics analytic workflow proposed that variations observed within these genes could be useful for genetic engineering–based improvement of hyaluronic acid–producing enzymes.

A Bacteriophage DNA Mimic Protein Employs a Non-specific Strategy to Inhibit the Bacterial RNA Polymerase

DNA mimicry by proteins is a strategy that employed by some proteins to occupy the binding sites of the DNA-binding proteins and deny further access to these sites by DNA. Such proteins have been found in bacteriophage, eukaryotic virus, prokaryotic, and eukaryotic cells to imitate non-coding functions of DNA. Here, we report another phage protein Gp44 from bacteriophage SPO1 of Bacillus subtilis, employing mimicry as part of unusual strategy to inhibit host RNA polymerase. Consisting of three simple domains, Gp44 contains a DNA binding motif, a flexible DNA mimic domain and a random-coiled domain. Gp44 is able to anchor to host genome and interact bacterial RNA polymerase via the β and β' subunit, resulting in bacterial growth inhibition. Our findings represent a non-specific strategy that SPO1 phage uses to target different bacterial transcription machinery regardless of the structural variations of RNA polymerases. This feature may have potential applications like generation of genetic engineered phages with Gp44 gene incorporated used in phage therapy to target a range of bacterial hosts.

Structural and Functional Investigation of the Periplasmic Arsenate-Binding Protein ArrX from Chrysiogenes arsenatis

The anaerobic bacterium Chrysiogenes arsenatis respires using the oxyanion arsenate (AsO₄^3-) as the terminal electron acceptor, where it is reduced to arsenite (AsO₃^3-) while concomitantly oxidizing various organic (e.g., acetate) electron donors. This respiratory activity is catalyzed in the periplasm of the bacterium by the enzyme arsenate reductase (Arr), with expression of the enzyme controlled by a sensor histidine kinase (ArrS) and a periplasmic-binding protein (PBP), ArrX. Here, we report for the first time, the molecular structure of ArrX in the absence and presence of bound ligand arsenate. Comparison of the ligand-bound structure of ArrX with other PBPs shows a high level of conservation of critical residues for ligand binding by these proteins; however, this suite of PBPs shows different structural alterations upon ligand binding. For ArrX and its homologue AioX (from Rhizobium sp. str. NT-26), which specifically binds arsenite, the structures of the substrate-binding sites in the vicinity of a conserved and critical cysteine residue contribute to the discrimination of binding for these chemically similar ligands.

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).

A pair of esterases from a commensal gut bacterium remove acetylations from all positions on complex β-mannans

β-mannans and xylans are important components of the plant cell wall and they are acetylated to be protected from degradation by glycoside hydrolases. β-mannans are widely present in human and animal diets as fiber from leguminous plants and as thickeners and stabilizers in processed foods. There are many fully characterized acetylxylan esterases (AcXEs); however, the enzymes deacetylating mannans are less understood. Here we present two carbohydrate esterases, RiCE2 and RiCE17, from the Firmicute Roseburia intestinalis, which together deacetylate complex galactoglucomannan (GGM). The three-dimensional (3D) structure of RiCE17 with a mannopentaose in the active site shows that the CBM35 domain of RiCE17 forms a confined complex, where the axially oriented C2-hydroxyl of a mannose residue points toward the Ser41 of the catalytic triad. Cavities on the RiCE17 surface may accept galactosylations at the C6 positions of mannose adjacent to the mannose residue being deacetylated (subsite -1 and +1). In-depth characterization of the two enzymes using time-resolved NMR, high-performance liquid chromatography (HPLC), and mass spectrometry demonstrates that they work in a complementary manner. RiCE17 exclusively removes the axially oriented 2-O-acetylations on any mannose residue in an oligosaccharide, including double acetylated mannoses, while the RiCE2 is active on 3-O-, 4-O-, and 6-O-acetylations. Activity of RiCE2 is dependent on RiCE17 removing 2-O-acetylations from double acetylated mannose. Furthermore, transacetylation of oligosaccharides with the 2-O-specific RiCE17 provided insight into how temperature and pH affects acetyl migration on manno-oligosaccharides.

Crystal structures of plant inorganic pyrophosphatase, an enzyme with a moonlighting autoproteolytic activity

Inorganic pyrophosphatases (PPases, EC 3.6.1.1), which hydrolyze inorganic pyrophosphate to phosphate in the presence of divalent metal cations, play a key role in maintaining phosphorus homeostasis in cells. DNA coding inorganic pyrophosphatases from Arabidopsis thaliana (AtPPA1) and Medicago truncatula (MtPPA1) were cloned into a bacterial expression vector and the proteins were produced in Escherichia coli cells and crystallized. In terms of their subunit fold, AtPPA1 and MtPPA1 are reminiscent of other members of Family I soluble pyrophosphatases from bacteria and yeast. Like their bacterial orthologs, both plant PPases form hexamers, as confirmed in solution by multi-angle light scattering and size-exclusion chromatography. This is in contrast with the fungal counterparts, which are dimeric. Unexpectedly, the crystallized AtPPA1 and MtPPA1 proteins lack ∼30 amino acid residues at their N-termini, as independently confirmed by chemical sequencing. In vitro, self-cleavage of the recombinant proteins is observed after prolonged storage or during crystallization. The cleaved fragment corresponds to a putative signal peptide of mitochondrial targeting, with a predicted cleavage site at Val31-Ala32. Site-directed mutagenesis shows that mutations of the key active site Asp residues dramatically reduce the cleavage rate, which suggests a moonlighting proteolytic activity. Moreover, the discovery of autoproteolytic cleavage of a mitochondrial targeting peptide would change our perception of this signaling process.

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the "open" hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine-rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.

Structure and Misfolding of the Flexible Tripartite Coiled-Coil Domain of Glaucoma-Associated Myocilin

Glaucoma-associated myocilin is a member of the olfactomedins, a protein family involved in neuronal development and human diseases. Molecular studies of the myocilin N-terminal coiled coil demonstrate a unique tripartite architecture: a Y-shaped parallel dimer-of-dimers with distinct tetramer and dimer regions. The structure of the dimeric C-terminal 7-heptad repeats elucidates an unexpected repeat pattern involving inter-strand stabilization by oppositely charged residues. Molecular dynamics simulations reveal an alternate accessible conformation in which the terminal inter-strand disulfide limits the extent of unfolding and results in a kinked configuration. By inference, full-length myocilin is also branched, with two pairs of C-terminal olfactomedin domains. Selected variants within the N-terminal region alter the apparent quaternary structure of myocilin but do so without compromising stability or causing aggregation. In addition to increasing our structural knowledge of naturally occurring extracellular coiled coils and biomedically important olfactomedins, this work broadens the scope of protein misfolding in the pathogenesis of myocilin-associated glaucoma.

The structure of Plasmodium falciparum 3D7_0606800 reveals a bi-lobed architecture that supports re-annotation as a Venus Flytrap protein

Plasmodium falciparum, the causative agent of malaria, employs a diverse array of surface displayed proteins to promote dissemination and establish infection in the human host. Of these, Pf3D7_0606800 is highly immunogenic and has been designated a potential top 10 candidate for inclusion in a multicomponent malarial vaccine. The role of Pf3D7_0606800 in parasite biology, however, is unknown and its characterization has been complicated by a lack of sequence identity with proteins of known structure or function. Towards elucidating Pf3D7_0606800 function, we determined its structure to a resolution of 2.35 Å using selenium single wavelength anomalous dispersion. A bi-lobed architecture displays the core structural hallmarks of Venus Flytrap (VFT) proteins prompting us to re-annotate Pf3D7_0606800 as PfVFT1. Structural analysis further revealed an extended inter-lobe groove that, when interrogated by molecular docking, appears well suited to bind peptide-based ligands. Collectively, our structural characterization of the highly antigenic P. falciparum surface protein PfVFT1 provides intriguing functional insight and establishes a structural template that could prove valuable for malaria vaccine engineering studies.

Nmd3 is a structural mimic of eIF5A, and activates the cpGTPase Lsg1 during 60S ribosome biogenesis

During ribosome biogenesis in eukaryotes, nascent subunits are exported to the cytoplasm in a functionally inactive state. 60S subunits are activated through a series of cytoplasmic maturation events. The last known events in the cytoplasm are the release of Tif6 by Efl1 and Sdo1 and the release of the export adapter, Nmd3, by the GTPase Lsg1. Here, we have used cryo-electron microscopy to determine the structure of the 60S subunit bound by Nmd3, Lsg1, and Tif6. We find that a central domain of Nmd3 mimics the translation elongation factor eIF5A, inserting into the E site of the ribosome and pulling the L1 stalk into a closed position. Additional domains occupy the P site and extend toward the sarcin-ricin loop to interact with Tif6. Nmd3 and Lsg1 together embrace helix 69 of the B2a intersubunit bridge, inducing base flipping that we suggest may activate the GTPase activity of Lsg1.

Structural Insights into the Tetrameric State of Aspartate-β-semialdehyde Dehydrogenases from Fungal Species

Aspartate-β-semialdehyde dehydrogenase (ASADH) catalyzes the second reaction in the aspartate pathway, a pathway required for the biosynthesis of one fifth of the essential amino acids in plants and microorganisms. Microarray analysis of a fungal pathogen T. rubrum responsible for most human dermatophytoses identified the upregulation of ASADH (trASADH) expression when the fungus is exposed to human skin, underscoring its potential as a drug target. Here we report the crystal structure of trASADH, revealing a tetrameric ASADH with a GAPDH-like fold. The tetramerization of trASADH was confirmed by sedimentation and SAXS experiments. Native PAGE demonstrated that this ASADH tetramerization is apparently universal in fungal species, unlike the functional dimer that is observed in all bacterial ASADHs. The helical subdomain in dimeric bacteria ASADH is replaced by the cover loop in archaeal/fungal ASADHs, presenting the determinant for this altered oligomerization. Mutations that disrupt the tetramerization of trASADH also abolish the catalytic activity, suggesting that the tetrameric state is required to produce the active fungal enzyme form. Our findings provide a basis to categorize ASADHs into dimeric and tetrameric enzymes, adopting a different orientation for NADP binding and offer a structural framework for designing drugs that can specifically target the fungal pathogens.

CombAlign: a code for generating a one-to-many sequence alignment from a set of pairwise structure-based sequence alignments

Background

In order to better define regions of similarity among related protein structures, it is useful to identify the residue-residue correspondences among proteins. Few codes exist for constructing a one-to-many multiple sequence alignment derived from a set of structure or sequence alignments, and a need was evident for creating such a tool for combining pairwise structure alignments that would allow for insertion of gaps in the reference structure.

Results

This report describes a new Python code, CombAlign, which takes as input a set of pairwise sequence alignments (which may be structure based) and generates a one-to-many, gapped, multiple structure- or sequence-based sequence alignment (MSSA). The use and utility of CombAlign was demonstrated by generating gapped MSSAs using sets of pairwise structure-based sequence alignments between structure models of the matrix protein (VP40) and pre-small/secreted glycoprotein (sGP) of Reston Ebolavirus and the corresponding proteins of several other filoviruses. The gapped MSSAs revealed structure-based residue-residue correspondences, which enabled identification of structurally similar versus differing regions in the Reston proteins compared to each of the other corresponding proteins.

Conclusions

CombAlign is a new Python code that generates a one-to-many, gapped, multiple structure- or sequence-based sequence alignment (MSSA) given a set of pairwise sequence alignments (which may be structure based). CombAlign has utility in assisting the user in distinguishing structurally conserved versus divergent regions on a reference protein structure relative to other closely related proteins. CombAlign was developed in Python 2.6, and the source code is available for download from the GitHub code repository.

Electronic supplementary material

The online version of this article (doi:10.1186/s13029-015-0039-1) contains supplementary material, which is available to authorized users.

Structural Characterization of the Complex between Hen Egg-White Lysozyme and Zr(IV) -Substituted Keggin Polyoxometalate as Artificial Protease

Successful co-crystallization of a noncovalent complex between hen egg-white lysozyme (HEWL) and the monomeric Zr(IV) -substituted Keggin polyoxometalate (POM) (Zr1 K1), (Et2 NH2)3 [Zr(PW11 O39)] (1), has been achieved, and its single-crystal X-ray structure has been determined. The dimeric Zr(IV) -substituted Keggin-type polyoxometalate (Zr1 K2), (Et2 NH2)10 [Zr(PW11 O39 )2] (2), has been previously shown to exhibit remarkable selectivity towards HEWL hydrolysis. The reported X-ray structure shows that the hydrolytically active Zr(IV) -substituted Keggin POM exists as a monomeric species. Prior to hydrolysis, this monomer interacts with HEWL in the vicinity of the previously identified cleavage sites found at Trp28-Val29 and Asn44-Arg45, through water-mediated H-bonding and electrostatic interactions. Three binding sites are observed at the interface of the negatively charged Keggin POM and the positively charged regions of HEWL at: 1) Gly16, Tyr20, and Arg21; 2) Asn44, Arg45, and Asn46; and 3) Arg128.

Thermostability of in vitro evolved Bacillus subtilis lipase A: a network and dynamics perspective

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure.

Common and distant structural characteristics of feruloyl esterase families from Aspergillus oryzae

Background

Feruloyl esterases (FAEs) are important biomass degrading accessory enzymes due to their capability of cleaving the ester links between hemicellulose and pectin to aromatic compounds of lignin, thus enhancing the accessibility of plant tissues to cellulolytic and hemicellulolytic enzymes. FAEs have gained increased attention in the area of biocatalytic transformations for the synthesis of value added compounds with medicinal and nutritional applications. Following the increasing attention on these enzymes, a novel descriptor based classification system has been proposed for FAEs resulting into 12 distinct families and pharmacophore models for three FAE sub-families have been developed.

Methodology/Principal Findings

The feruloylome of Aspergillus oryzae contains 13 predicted FAEs belonging to six sub-families based on our recently developed descriptor-based classification system. The three-dimensional structures of the 13 FAEs were modeled for structural analysis of the feruloylome. The three genes coding for three enzymes, viz., A.O.2, A.O.8 and A.O.10 from the feruloylome of A. oryzae, representing sub-families with unknown functional features, were heterologously expressed in Pichia pastoris, characterized for substrate specificity and structural characterization through CD spectroscopy. Common feature-based pharamacophore models were developed according to substrate specificity characteristics of the three enzymes. The active site residues were identified for the three expressed FAEs by determining the titration curves of amino acid residues as a function of the pH by applying molecular simulations.

Conclusions/Significance

Our findings on the structure-function relationships and substrate specificity of the FAEs of A. oryzae will be instrumental for further understanding of the FAE families in the novel classification system. The developed pharmacophore models could be applied for virtual screening of compound databases for short listing the putative substrates prior to docking studies or for post-processing docking results to remove false positives. Our study exemplifies how computational predictions can complement to the information obtained through experimental methods.

Implications for proteasome nuclear localization revealed by the structure of the nuclear proteasome tether protein Cut8

Degradation of nuclear proteins by the 26S proteasome is essential for cell viability. In yeast, the nuclear envelope protein Cut8 mediates nuclear proteasomal sequestration by an uncharacterized mechanism. Here we describe structures of Schizosaccharomyces pombe Cut8, which shows that it contains a unique, modular fold composed of an extended N-terminal, lysine-rich segment that when ubiquitinated binds the proteasome, a dimer domain followed by a six-helix bundle connected to a flexible C tail. The Cut8 six-helix bundle shows structural similarity to 14-3-3 phosphoprotein-binding domains, and binding assays show that this domain is necessary and sufficient for liposome and cholesterol binding. Moreover, specific mutations in the 14-3-3 regions corresponding to putative cholesterol recognition/interaction amino acid consensus motifs abrogate cholesterol binding. In vivo studies confirmed that the 14-3-3 region is necessary for Cut8 membrane localization and that dimerization is critical for its function. Thus, the data reveal the Cut8 organization at the nuclear envelope. Reconstruction of Cut8 evolution suggests that it was present in the last common ancestor of extant eukaryotes and accordingly that nuclear proteasomal sequestration is an ancestral eukaryotic feature. The importance of Cut8 for cell viability and its absence in humans suggests it as a possible target for the development of specific chemotherapeutics against invasive fungal infections.