Surprising similarities in structure comparison

Crystal structure and activation mechanism of DR3 death domain

Death receptor 3 (DR3) (a.k.a. tumor necrosis factor receptor superfamily 25) plays a key role in the immune system by activating nuclear factor kappa-light-chain-enhancer of activated B cells signaling pathway. Here we present the crystal structures of human and mouse DR3 intracellular death domain (DD) at 2.7 and 2.5 Å resolutions, respectively. The mouse DR3 DD adopts a classical six-helix bundle structure while human DR3 DD displays an extended fold. Though there is one-amino-acid difference in the linker between maltose-binding protein (MBP) tag and DR3 DD, according to our self-interaction analysis, the hydrophobic interface discovered in MBP-hDR3 DD crystal structure is responsible for both hDR3 DD and mDR3 DD homotypic interaction. Furthermore, our biochemical analysis indicates that the sequence variation between human and mouse DR3 DD does not affect its structure and function. Small-angle X-ray scattering analysis shows the averaged solution structures of both human and mouse MBP-DR3 DD are the combination of different conformations with different proportion. Through switching to the open conformation, DR3 DD could improve the interaction with downstream element TNFR-associated death domain (TRADD). Here we propose an activation-dependent structural rearrangement model: the DD region is folded as the six-helix bundles in the resting state, while upon extracellular ligand engagement, it switches to the open conformation, which facilitates its self-association and the recruitment of TRADD. Our results provide detailed insights into the architecture of DR3 DD and the molecular mechanism of activation. DATABASES: All refined structure coordinates as well as the corresponding structure factors have been deposited in the PDB under the accession codes 5YGS, 5YEV, 5YGP, 5ZNY, 5ZNZ.

RUPEE: A fast and accurate purely geometric protein structure search

Given the close relationship between protein structure and function, protein structure searches have long played an established role in bioinformatics. Despite their maturity, existing protein structure searches either use simplifying assumptions or compromise between fast response times and quality of results. These limitations can prevent the easy and efficient exploration of relationships between protein structures, which is the norm in other areas of inquiry. To address these limitations we have developed RUPEE, a fast and accurate purely geometric structure search combining techniques from information retrieval and big data with a novel approach to encoding sequences of torsion angles. Comparing our results to the output of mTM, SSM, and the CATHEDRAL structural scan, it is clear that RUPEE has set a new bar for purely geometric big data approaches to protein structure searches. RUPEE in top-aligned mode produces equal or better results than the best available protein structure searches, and RUPEE in fast mode demonstrates the fastest response times coupled with high quality results. The RUPEE protein structure search is available at https://ayoubresearch.com. Code and data are available at https://github.com/rayoub/rupee.

Unsupervised Discovery of Geometrically Common Structural Motifs and Long-Range Contacts in Protein 3D Structures

The essential role of small evolutionarily conserved structural units in proteins has been extensively researched and validated. A popular example are serine proteases, where the peptide cleavage reaction is realized by a configuration of only three residues. Brought to spatial proximity during the protein folding process, such structural motifs are often long-range contacts and usually hard to detect at sequence level. Due to the constantly increasing resource of protein 3D structure data, the computational identification of structural motifs can contribute significantly to the understanding of protein fold and function. Thus, we propose a method to discover structural motifs of high geometrical similarity and desired sequence separation in protein 3D structure data. By utilizing methods originated from data mining, no a priori knowledge is required. The applicability of the method is demonstrated by the identification of the catalytic unit of serine proteases and the ion-coordination center of cupredoxins. Furthermore, large-scale analysis of the entire Protein Data Bank points towards the presence of ubiquitous structural motifs, independent of any specific fold or function. We envision that our method is suitable to uncover functional mechanisms and to derive fingerprint libraries of structural motifs, which could be used to assess protein family association.

Characterization of the housekeeping sortase from the human pathogen Propionibacterium acnes: first investigation of a class F sortase

Sortase enzymes play an important role in Gram-positive bacteria. They are responsible for the covalent attachment of proteins to the surface of the bacteria and perform this task via a highly sequence-specific transpeptidation reaction. Since these immobilized proteins are often involved in pathogenicity of Gram-positive bacteria, characterization of this type of enzyme is also of medical relevance. Different classes of sortases (A-F) have been found, which recognize characteristic recognition sequences present in substrate proteins. Up to date, sortase A from Staphylococcus aureus, a housekeeping class A sortase, is the most thoroughly studied representative of the sortase family of enzymes. Here we report the in-depth characterization of the class F sortase from Propionibacterium acnes, a class of sortases that has not been investigated before. As Sortase F is the only transpeptidase found in the P. acnes genome, it is the housekeeping sortase of this organism. Sortase F from P. acnes shows a behavior similar to sortases from class A in terms of pH dependence, recognition sequence and catalytic activity; furthermore, its activity is independent of bivalent ions, which contrasts to sortase A from S. aureus We demonstrate that sortase F is useful for protein engineering applications, by producing a site-specifically conjugated homogenous antibody-drug conjugate with a potency similar to that of a conjugate prepared with sortase A. Thus, the detailed characterization presented here will not only enable the development of anti-virulence agents targeting P. acnes but also provides a powerful alternative to sortase A for protein engineering applications.

Scalable Extraction of Big Macromolecular Data in Azure Data Lake Environment

Calculation of structural features of proteins, nucleic acids, and nucleic acid-protein complexes on the basis of their geometries and studying various interactions within these macromolecules, for which high-resolution structures are stored in Protein Data Bank (PDB), require parsing and extraction of suitable data stored in text files. To perform these operations on large scale in the face of the growing amount of macromolecular data in public repositories, we propose to perform them in the distributed environment of Azure Data Lake and scale the calculations on the Cloud. In this paper, we present dedicated data extractors for PDB files that can be used in various types of calculations performed over protein and nucleic acids structures in the Azure Data Lake. Results of our tests show that the Cloud storage space occupied by the macromolecular data can be successfully reduced by using compression of PDB files without significant loss of data processing efficiency. Moreover, our experiments show that the performed calculations can be significantly accelerated when using large sequential files for storing macromolecular data and by parallelizing the calculations and data extractions that precede them. Finally, the paper shows how all the calculations can be performed in a declarative way in U-SQL scripts for Data Lake Analytics.

Evolutionary dynamics of origin and loss in the deep history of phospholipase D toxin genes

Background

Venom-expressed sphingomyelinase D/phospholipase D (SMase D/PLD) enzymes evolved from the ubiquitous glycerophosphoryl diester phosphodiesterases (GDPD). Expression of GDPD-like SMaseD/PLD toxins in both arachnids and bacteria has inspired consideration of the relative contributions of lateral gene transfer and convergent recruitment in the evolutionary history of this lineage. Previous work recognized two distinct lineages, a SicTox-like (ST-like) clade including the arachnid toxins, and an Actinobacterial-toxin like (AT-like) clade including the bacterial toxins and numerous fungal homologs.

Results

Here we expand taxon sampling by homology detection to discover new GDPD-like SMase D/PLD homologs. The ST-like clade now includes homologs in a wider variety of arthropods along with a sister group in Cnidaria; the AT-like clade now includes additional fungal phyla and proteobacterial homologs; and we report a third clade expressed in diverse aquatic metazoan taxa, a few single-celled eukaryotes, and a few aquatic proteobacteria. GDPD-like SMaseD/PLDs have an ancient presence in chelicerates within the ST-like family and ctenophores within the Aquatic family. A rooted phylogenetic tree shows that the three clades derived from a basal paraphyletic group of proteobacterial GDPD-like SMase D/PLDs, some of which are on mobile genetic elements. GDPD-like SMase D/PLDs share a signature C-terminal motif and a shortened βα1 loop, features that distinguish them from GDPDs. The three major clades also have active site loop signatures that distinguish them from GDPDs and from each other. Analysis of molecular phylogenies with respect to organismal relationships reveals a dynamic evolutionary history including both lateral gene transfer and gene duplication/loss.

Conclusions

The GDPD-like SMaseD/PLD enzymes derive from a single ancient ancestor, likely proteobacterial, and radiated into diverse organismal lineages at least in part through lateral gene transfer.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1302-2) contains supplementary material, which is available to authorized users.

The application of Hadoop in structural bioinformatics

The paper reviews the use of the Hadoop platform in structural bioinformatics applications. For structural bioinformatics, Hadoop provides a new framework to analyse large fractions of the Protein Data Bank that is key for high-throughput studies of, for example, protein-ligand docking, clustering of protein-ligand complexes and structural alignment. Specifically we review in the literature a number of implementations using Hadoop of high-throughput analyses and their scalability. We find that these deployments for the most part use known executables called from MapReduce rather than rewriting the algorithms. The scalability exhibits a variable behaviour in comparison with other batch schedulers, particularly as direct comparisons on the same platform are generally not available. Direct comparisons of Hadoop with batch schedulers are absent in the literature but we note there is some evidence that Message Passing Interface implementations scale better than Hadoop. A significant barrier to the use of the Hadoop ecosystem is the difficulty of the interface and configuration of a resource to use Hadoop. This will improve over time as interfaces to Hadoop, e.g. Spark improve, usage of cloud platforms (e.g. Azure and Amazon Web Services (AWS)) increases and standardised approaches such as Workflow Languages (i.e. Workflow Definition Language, Common Workflow Language and Nextflow) are taken up.

Fast design of arbitrary length loops in proteins using InteractiveRosetta

Background

With increasing interest in ab initio protein design, there is a desire to be able to fully explore the design space of insertions and deletions. Nature inserts and deletes residues to optimize energy and function, but allowing variable length indels in the context of an interactive protein design session presents challenges with regard to speed and accuracy.

Results

Here we present a new module (INDEL) for InteractiveRosetta which allows the user to specify a range of lengths for a desired indel, and which returns a set of low energy backbones in a matter of seconds. To make the loop search fast, loop anchor points are geometrically hashed using C α-C α and C β-C β distances, and the hash is mapped to start and end points in a pre-compiled random access file of non-redundant, protein backbone coordinates. Loops with superposable anchors are filtered for collisions and returned to InteractiveRosetta as poly-alanine for display and selective incorporation into the design template. Sidechains can then be added using RosettaDesign tools.

Conclusions

INDEL was able to find viable loops in 100% of 500 attempts for all lengths from 3 to 20 residues. INDEL has been applied to the task of designing a domain-swapping loop for T7-endonuclease I, changing its specificity from Holliday junctions to paranemic crossover (PX) DNA.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2345-5) contains supplementary material, which is available to authorized users.

Carcinogenic Pesticide Control via Hijacking Endosymbiosis; The Paradigm of DSB-A from Wolbachia pipientis for the Management of Otiorhynchus singularis

BACKGROUND/AIM

Pesticides have little, if any specificity, to the pathogen they target in most cases. Wide spectrum toxic chemicals are being used to remove pestcides and salvage crops and economies linked to agriculture. The burden on the environment, public health and economy is huge. Traditional pestcide control is based on administering heavy loads of highly toxic compounds and elements that essentially strip all life from the field. Those chemicals are a leading cause of increased cancer related deaths in countryside. Herein, the Trojan horse of endosymbiosis was used, in an effort to control pests using high specificity compounds in reduced quantities.

MATERIALS AND METHODS

Our pipeline has been applied on the case of Otiorhynchus singularis, which is a very widespread pest, whose impact is devastating on a repertoire of crops. To date, there is no specific pesticide nor agent to control it. The deployed strategy involves the inhibition of the key DSB-A enzyme of its endosymbiotic Wolbachia pipientis bacterial strain.

RESULTS

Our methodology, provides the means to design, test and identify highly specific pestcide control substances that minimize the impact of toxic chemicals on health, economy and the environment.

CONCLUSION

All in all, in this study a radical computer-based pipeline is proposed that could be adopted under many other similar scenarios and pave the way for precision agriculture via optimized pest control.

Malaria parasite translocon structure and mechanism of effector export

The putative Plasmodium translocon of exported proteins (PTEX) is essential for transport of malarial effector proteins across a parasite-encasing vacuolar membrane into host erythrocytes, but the mechanism of this process remains unknown. Here we show that PTEX is a bona fide translocon by determining structures of the PTEX core complex at near-atomic resolution using cryo-electron microscopy. We isolated the endogenous PTEX core complex containing EXP2, PTEX150 and HSP101 from Plasmodium falciparum in the 'engaged' and 'resetting' states of endogenous cargo translocation using epitope tags inserted using the CRISPR-Cas9 system. In the structures, EXP2 and PTEX150 interdigitate to form a static, funnel-shaped pseudo-seven-fold-symmetric protein-conducting channel spanning the vacuolar membrane. The spiral-shaped AAA+ HSP101 hexamer is tethered above this funnel, and undergoes pronounced compaction that allows three of six tyrosine-bearing pore loops lining the HSP101 channel to dissociate from the cargo, resetting the translocon for the next threading cycle. Our work reveals the mechanism of P. falciparum effector export, and will inform structure-based design of drugs targeting this unique translocon.

Crystal structure of the Streptococcus agalactiae CAMP factor provides insights into its membrane-permeabilizing activity

Streptococcus agalactiae is an important human opportunistic pathogen that can cause serious health problems, particularly among newborns and older individuals. S. agalactiae contains the CAMP factor, a pore-forming toxin first identified in this bacterium. The CAMP reaction is based on the co-hemolytic activity of the CAMP factor and is commonly used to identify S. agalactiae in the clinic. Closely related proteins are present also in other Gram-positive pathogens. Although the CAMP toxin was discovered more than a half century ago, no structure from this toxin family has been reported, and the mechanism of action of this toxin remains unclear. Here, we report the first structure of this toxin family, revealing a structural fold composed of 5 + 3-helix bundles. Further analysis by protein truncation and site-directed mutagenesis indicated that the N-terminal 5-helix bundle is responsible for membrane permeabilization, whereas the C-terminal 3-helix bundle is likely responsible for host receptor binding. Interestingly, the C-terminal domain inhibited the activity of both full-length toxin and its N-terminal domain. Moreover, we observed that the linker region is highly conserved and has a conserved DLXXXDXAT sequence motif. Structurally, this linker region extensively interacted with both terminal CAMP factor domains, and mutagenesis disclosed that the conserved sequence motif is required for CAMP factor's co-hemolytic activity. In conclusion, our results reveal a unique structure of this bacterial toxin and help clarify the molecular mechanism of its co-hemolytic activity.

Functional and Structural Diversity of Acyl-coA Binding Proteins in Oil Crops

Diversities in structure and function of ACBP were discussed in this review. ACBP are important proteins that could transport newly synthesized fatty acid, activated into -coA, from plastid to endoplasmic reticulum, where oil in the form of triacylglycerol occurs. ACBP were detected in various animal and plants species, which indicated their importance in biological function. In fact, involvement of ACBP in important process such as lipid metabolism, regulation of enzyme and gene expression, and in response to plant stresses has been proven in several studies. In this review, findings on ACBP of 11 well-known oil crops were reviewed to comprehend diversity, comparative analyses on ACBP structure were made, and link between structure and function, tissue expression and subcellular location of ACBP were also observed. Incomplete reports in some species were mentioned, which might be encouraging to start or to perform deeper studies. Similar characteristics were found in paralogs ACBP, and orthologs ACBP had different functions, despite the high identity in amino acid sequence. At the end, it is confirmed that ortholog proteins could not necessarily display the same function, even from closely related species.

GPU-Based Point Cloud Superpositioning for Structural Comparisons of Protein Binding Sites

In this paper, we present a novel approach to solve the labeled point cloud superpositioning problem for performing structural comparisons of protein binding sites. The solution is based on a parallel evolution strategy that operates on large populations and runs on GPU hardware. The proposed evolution strategy reduces the likelihood of getting stuck in a local optimum of the multimodal real-valued optimization problem represented by labeled point cloud superpositioning. The performance of the GPU-based parallel evolution strategy is compared to a previously proposed CPU-based sequential approach for labeled point cloud superpositioning, indicating that the GPU-based parallel evolution strategy leads to qualitatively better results and significantly shorter runtimes, with speed improvements of up to a factor of 1,500 for large populations. Binary classification tests based on the ATP, NADH, and FAD protein subsets of CavBase, a database containing putative binding sites, show average classification rate improvements from about 92 percent (CPU) to 96 percent (GPU). Further experiments indicate that the proposed GPU-based labeled point cloud superpositioning approach can be superior to traditional protein comparison approaches based on sequence alignments.

Architecture of the human GATOR1 and GATOR1-Rag GTPases complexes

Nutrients, such as amino acids and glucose, signal through the Rag GTPases to activate mTORC1. The GATOR1 protein complex-comprising DEPDC5, NPRL2 and NPRL3-regulates the Rag GTPases as a GTPase-activating protein (GAP) for RAGA; loss of GATOR1 desensitizes mTORC1 signalling to nutrient starvation. GATOR1 components have no sequence homology to other proteins, so the function of GATOR1 at the molecular level is currently unknown. Here we used cryo-electron microscopy to solve structures of GATOR1 and GATOR1-Rag GTPases complexes. GATOR1 adopts an extended architecture with a cavity in the middle; NPRL2 links DEPDC5 and NPRL3, and DEPDC5 contacts the Rag GTPase heterodimer. Biochemical analyses reveal that our GATOR1-Rag GTPases structure is inhibitory, and that at least two binding modes must exist between the Rag GTPases and GATOR1. Direct interaction of DEPDC5 with RAGA inhibits GATOR1-mediated stimulation of GTP hydrolysis by RAGA, whereas weaker interactions between the NPRL2-NPRL3 heterodimer and RAGA execute GAP activity. These data reveal the structure of a component of the nutrient-sensing mTORC1 pathway and a non-canonical interaction between a GAP and its substrate GTPase.

Structural and functional modeling of viral protein 5 of Infectious Bursal Disease Virus

Infectious Bursal Disease (IBD) is an acute, highly contagious and immunosuppressive disease of young chicken. The causative virus (IBDV) is a bi-segmented, double-stranded RNA virus. The virus encodes five major proteins, viral protein (VP) 1-5. VPs 1-3 have been characterized crystallographically. Albeit a rise in the number of studies reporting successful heterologous expression of VP5 in recent times, challenging the notion that rapid death of host cells overexpressing VP5 disallows obtaining sufficiently pure preparations of the protein for crystallographic studies, the structure of VP5 remains unknown and its function controversial. Our study describes the first 3D model of IBD VP5 obtained through an elaborate computational workflow. Based on the results of the study, IBD VP5 can be predicted to be a structural analog of the leucine-rich repeat (LRR) family of proteins. Functional implications arising from structural similarity of VP5 with host Toll-like receptor (Tlr) 3 also satisfy the previously reported opposing roles of the protein in first abolishing and later inducing host-cell apoptosis.

The Classification of Protein Domains

The significant expansion in protein sequence and structure data that we are now witnessing brings with it a pressing need to bring order to the protein world. Such order enables us to gain insights into the evolution of proteins, their function and the extent to which the functional repertoire can vary across the three kingdoms of life. This has lead to the creation of a wide range of protein family classifications that aim to group proteins based upon their evolutionary relationships.In this chapter we discuss the approaches and methods that are frequently used in the classification of proteins, with a specific emphasis on the classification of protein domains. The construction of both domain sequence and domain structure databases is considered and we show how the use of domain family annotations to assign structural and functional information is enhancing our understanding of genomes.

The crystal structure of pseudokinase PEAK1 (Sugen kinase 269) reveals an unusual catalytic cleft and a novel mode of kinase fold dimerization

The pseudokinase group encompasses some 10% of protein kinases, but pseudokinases diverge from canonical kinases in key motifs. The two members of the small new kinase family 3 (NKF3) group are considered pseudokinases. These proteins, pseudopodium-enriched atypical kinase 1 (PEAK1, Sugen kinase 269, or SgK269) and pragmin (Sugen kinase 223 or SgK223), act as scaffolds in growth factor signaling pathways, and both contain a kinase fold with degraded kinase motifs at their C termini. These kinases may harbor regions that mediate oligomerization or control other aspects of signal transduction, but a lack of structural information has precluded detailed investigations into their functional roles. In this study, we determined the X-ray crystal structure of the PEAK1 pseudokinase domain to 2.3 Å resolution. The structure revealed that the PEAK1 kinase-like domain contains a closed nucleotide-binding cleft that in this conformation may deleteriously affect nucleotide binding. Moreover, we found that N- and C-terminal extensions create a highly unusual all α-helical split-dimerization region, termed here the split helical dimerization (SHED) region. Sequence conservation analysis suggested that this region facilitates a dimerization mode that is conserved between PEAK1 and pragmin. Finally, we observed structural similarities between the PEAK1 SHED region and the C-terminal extension of the Parkinson's disease-associated kinase PINK1. In summary, PEAK1's kinase cleft is occluded, and its newly identified SHED region may promote an unexpected dimerization mode. Similarities of PEAK1 with the active kinase PINK1 may reclassify the latter as a member of the new kinase family 3 group.

Bayesian comparison of protein structures using partial Procrustes distance

An important topic in bioinformatics is the protein structure alignment. Some statistical methods have been proposed for this problem, but most of them align two protein structures based on the global geometric information without considering the effect of neighbourhood in the structures. In this paper, we provide a Bayesian model to align protein structures, by considering the effect of both local and global geometric information of protein structures. Local geometric information is incorporated to the model through the partial Procrustes distance of small substructures. These substructures are composed of β-carbon atoms from the side chains. Parameters are estimated using a Markov chain Monte Carlo (MCMC) approach. We evaluate the performance of our model through some simulation studies. Furthermore, we apply our model to a real dataset and assess the accuracy and convergence rate. Results show that our model is much more efficient than previous approaches.

Modeling of RAS complexes supports roles in cancer for less studied partners

Background

RAS protein interactions have predominantly been studied in the context of the RAF and PI3kinase oncogenic pathways. Structural modeling and X-ray crystallography have demonstrated that RAS isoforms bind to canonical downstream effector proteins in these pathways using the highly conserved switch I and II regions. Other non-canonical RAS protein interactions have been experimentally identified, however it is not clear whether these proteins also interact with RAS via the switch regions.

Results

To address this question we constructed a RAS isoform-specific protein-protein interaction network and predicted 3D complexes involving RAS isoforms and interaction partners to identify the most probable interaction interfaces. The resulting models correctly captured the binding interfaces for well-studied effectors, and additionally implicated residues in the allosteric and hyper-variable regions of RAS proteins as the predominant binding site for non-canonical effectors. Several partners binding to this new interface (SRC, LGALS1, RABGEF1, CALM and RARRES3) have been implicated as important regulators of oncogenic RAS signaling. We further used these models to investigate competitive binding and multi-protein complexes compatible with RAS surface occupancy and the putative effects of somatic mutations on RAS protein interactions.

Conclusions

We discuss our findings in the context of RAS localization to the plasma membrane versus within the cytoplasm and provide a list of RAS protein interactions with possible cancer-related consequences, which could help guide future therapeutic strategies to target RAS proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s13628-017-0037-6) contains supplementary material, which is available to authorized users.

Development of a Novel Lead that Targets M. tuberculosis Polyketide Synthase 13

Widespread resistance to first-line TB drugs is a major problem that will likely only be resolved through the development of new drugs with novel mechanisms of action. We have used structure-guided methods to develop a lead molecule that targets the thioesterase activity of polyketide synthase Pks13, an essential enzyme that forms mycolic acids, required for the cell wall of Mycobacterium tuberculosis. Our lead, TAM16, is a benzofuran class inhibitor of Pks13 with highly potent in vitro bactericidal activity against drug-susceptible and drug-resistant clinical isolates of M. tuberculosis. In multiple mouse models of TB infection, TAM16 showed in vivo efficacy equal to the first-line TB drug isoniazid, both as a monotherapy and in combination therapy with rifampicin. TAM16 has excellent pharmacological and safety profiles, and the frequency of resistance for TAM16 is ∼100-fold lower than INH, suggesting that it can be developed as a new antitubercular aimed at the acute infection. PAPERCLIP.

Structure Reveals Mechanisms of Viral Suppressors that Intercept a CRISPR RNA-Guided Surveillance Complex

Genetic conflict between viruses and their hosts drives evolution and genetic innovation. Prokaryotes evolved CRISPR-mediated adaptive immune systems for protection from viral infection, and viruses have evolved diverse anti-CRISPR (Acr) proteins that subvert these immune systems. The adaptive immune system in Pseudomonas aeruginosa (type I-F) relies on a 350 kDa CRISPR RNA (crRNA)-guided surveillance complex (Csy complex) to bind foreign DNA and recruit a trans-acting nuclease for target degradation. Here, we report the cryo-electron microscopy (cryo-EM) structure of the Csy complex bound to two different Acr proteins, AcrF1 and AcrF2, at an average resolution of 3.4 Å. The structure explains the molecular mechanism for immune system suppression, and structure-guided mutations show that the Acr proteins bind to residues essential for crRNA-mediated detection of DNA. Collectively, these data provide a snapshot of an ongoing molecular arms race between viral suppressors and the immune system they target.

Structure and Dimerization of IreB, a Negative Regulator of Cephalosporin Resistance in Enterococcus faecalis

Enterococcus faecalis, a leading cause of hospital-acquired infections, exhibits intrinsic resistance to most cephalosporins, which are antibiotics in the beta-lactam family that target cell-wall biosynthesis. A comprehensive understanding of the underlying genetic and biochemical mechanisms of cephalosporin resistance in E. faecalis is lacking. We previously determined that a transmembrane serine/threonine kinase (IreK) and its cognate phosphatase (IreP) reciprocally regulate cephalosporin resistance in E. faecalis, dependent on the kinase activity of IreK. Other than IreK itself, thus far the only known substrate for reversible phosphorylation by IreK and IreP is IreB, a small protein of unknown function that is well conserved in low-GC Gram-positive bacteria. We previously showed that IreB acts as a negative regulator of cephalosporin resistance in E. faecalis. However, the biochemical mechanism by which IreB modulates cephalosporin resistance remains unknown. As a next step toward an understanding of the mechanism by which IreB regulates resistance, we initiated a structure-function study on IreB. The NMR solution structure of IreB was determined, revealing that IreB adopts a unique fold and forms a dimer in vitro. Dimerization of IreB was confirmed in vivo. Substitutions at the dimer interface impaired IreB function and stability in vivo, indicating that dimerization is functionally important for the biological activity of IreB. Hence, these studies provide new insights into the structure and function of a widely conserved protein of unknown function that is an important regulator of antimicrobial resistance in E. faecalis.

Structure of the Toll/Interleukin-1 Receptor (TIR) Domain of the B-cell Adaptor That Links Phosphoinositide Metabolism with the Negative Regulation of the Toll-like Receptor (TLR) Signalosome

Ligand binding to Toll-like receptors (TLRs) results in dimerization of their cytosolic Toll/interleukin-1 receptor (TIR) domains and recruitment of post-receptor signal transducers into a complex signalosome. TLR activation leads to the production of transcription factors and pro-inflammatory molecules and the activation of phosphoinositide 3-kinases (PI3K) in a process that requires the multimodular B-cell adaptor for phosphoinositide 3-kinase (BCAP). BCAP has a sequence previously proposed as a "cryptic" TIR domain. Here, we present the structure of the N-terminal region of human BCAP and show that it possesses a canonical TIR fold. Dimeric BCAP associates with the TIR domains of TLR2/4 and MAL/TIRAP, suggesting that it is recruited to the TLR signalosome by multitypic TIR-TIR interactions. BCAP also interacts with the p85 subunit of PI3K and phospholipase Cγ, enzymes that deplete plasma membrane phosphatidylinositol 4,5-bisphosphate (PIP2), and these interactions provide a molecular explanation for BCAP-mediated down-regulation of inflammatory signaling.

Protein Function Prediction

Protein function is a concept that can have different interpretations in different biological contexts, and the number and diversity of novel proteins identified by large-scale "omics" technologies poses increasingly new challenges. In this review we explore current strategies used to predict protein function focused on high-throughput sequence analysis, as for example, inference based on sequence similarity, sequence composition, structure, and protein-protein interaction. Various prediction strategies are discussed together with illustrative workflows highlighting the use of some benchmark tools and knowledge bases in the field.

Mechanism of arginine sensing by CASTOR1 upstream of mTORC1

The mechanistic Target of Rapamycin Complex 1 (mTORC1) is a major regulator of eukaryotic growth that coordinates anabolic and catabolic cellular processes with inputs such as growth factors and nutrients, including amino acids. In mammals arginine is particularly important, promoting diverse physiological effects such as immune cell activation, insulin secretion, and muscle growth, largely mediated through activation of mTORC1 (refs 4, 5, 6, 7). Arginine activates mTORC1 upstream of the Rag family of GTPases, through either the lysosomal amino acid transporter SLC38A9 or the GATOR2-interacting Cellular Arginine Sensor for mTORC1 (CASTOR1). However, the mechanism by which the mTORC1 pathway detects and transmits this arginine signal has been elusive. Here, we present the 1.8 Å crystal structure of arginine-bound CASTOR1. Homodimeric CASTOR1 binds arginine at the interface of two Aspartate kinase, Chorismate mutase, TyrA (ACT) domains, enabling allosteric control of the adjacent GATOR2-binding site to trigger dissociation from GATOR2 and downstream activation of mTORC1. Our data reveal that CASTOR1 shares substantial structural homology with the lysine-binding regulatory domain of prokaryotic aspartate kinases, suggesting that the mTORC1 pathway exploited an ancient, amino-acid-dependent allosteric mechanism to acquire arginine sensitivity. Together, these results establish a structural basis for arginine sensing by the mTORC1 pathway and provide insights into the evolution of a mammalian nutrient sensor.

Homology modeling, substrate docking, and molecular simulation studies of mycobacteriophage Che12 lysin A

Mycobacteriophages produce lysins that break down the host cell wall at the end of lytic cycle to release their progenies. The ability to lyse mycobacterial cells makes the lysins significant. Mycobacteriophage Che12 is the first reported temperate phage capable of infecting and lysogenising Mycobacterium tuberculosis. Gp11 of Che12 was found to have Chitinase domain that serves as endolysin (lysin A) for Che12. Structure of gp11 was modeled and evaluated using Ramachandran plot in which 98 % of the residues are in the favored and allowed regions. Che12 lysin A was predicted to act on NAG-NAM-NAG molecules in the peptidoglycan of cell wall. The tautomers of NAG-NAM-NAG molecule were generated and docked with lysin A. The stability and binding affinity of lysin A - NAG-NAM-NAG tautomers were studied using molecular dynamics simulations.

Loop-to-helix transition in the structure of multidrug regulator AcrR at the entrance of the drug-binding cavity

Multidrug transcription regulator AcrR from Salmonella enterica subsp. enterica serovar Typhimurium str. LT2 belongs to the tetracycline repressor family, one of the largest groups of bacterial transcription factors. The crystal structure of dimeric AcrR was determined and refined to 1.56Å resolution. The tertiary and quaternary structures of AcrR are similar to those of its homologs. The multidrug binding site was identified based on structural alignment with homologous proteins and has a di(hydroxyethyl)ether molecule bound. Residues from helices α4 and α7 shape the entry into this binding site. The structure of AcrR reveals that the extended helical conformation of helix α4 is stabilized by the hydrogen bond between Glu67 (helix α4) and Gln130 (helix α7). Based on the structural comparison with the closest homolog structure, the Escherichia coli AcrR, we propose that this hydrogen bond is responsible for control of the loop-to-helix transition within helix α4. This local conformational switch of helix α4 may be a key step in accessing the multidrug binding site and securing ligands at the binding site. Solution small-molecule binding studies suggest that AcrR binds ligands with their core chemical structure resembling the tetracyclic ring of cholesterol.

MS3ALIGN: an efficient molecular surface aligner using the topology of surface curvature

Background

Aligning similar molecular structures is an important step in the process of bio-molecular structure and function analysis. Molecular surfaces are simple representations of molecular structure that are easily constructed from various forms of molecular data such as 3D atomic coordinates (PDB) and Electron Microscopy (EM) data.

Methods

We present a Multi-Scale Morse-Smale Molecular-Surface Alignment tool, MS3ALIGN, which aligns molecular surfaces based on significant protrusions on the molecular surface. The input is a pair of molecular surfaces represented as triangle meshes. A key advantage of MS3ALIGN is computational efficiency that is achieved because it processes only a few carefully chosen protrusions on the molecular surface. Furthermore, the alignments are partial in nature and therefore allows for inexact surfaces to be aligned.

Results

The method is evaluated in four settings. First, we establish performance using known alignments with varying overlap and noise values. Second, we compare the method with SurfComp, an existing surface alignment method. We show that we are able to determine alignments reported by SurfComp, as well as report relevant alignments not found by SurfComp. Third, we validate the ability of MS3ALIGN to determine alignments in the case of structurally dissimilar binding sites. Fourth, we demonstrate the ability of MS3ALIGN to align iso-surfaces derived from cryo-electron microscopy scans.

Conclusions

We have presented an algorithm that aligns Molecular Surfaces based on the topology of surface curvature. A webserver and standalone software implementation of the algorithm available at http://vgl.serc.iisc.ernet.in/ms3align.

Electronic supplementary material

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

Efficient and automated large-scale detection of structural relationships in proteins with a flexible aligner

Background

The total number of known three-dimensional protein structures is rapidly increasing. Consequently, the need for fast structural search against complete databases without a significant loss of accuracy is increasingly demanding. Recently, TopSearch, an ultra-fast method for finding rigid structural relationships between a query structure and the complete Protein Data Bank (PDB), at the multi-chain level, has been released. However, comparable accurate flexible structural aligners to perform efficient whole database searches of multi-domain proteins are not yet available. The availability of such a tool is critical for a sustainable boosting of biological discovery.

Results

Here we report on the development of a new method for the fast and flexible comparison of protein structure chains. The method relies on the calculation of 2D matrices containing a description of the three-dimensional arrangement of secondary structure elements (angles and distances). The comparison involves the matching of an ensemble of substructures through a nested-two-steps dynamic programming algorithm. The unique features of this new approach are the integration and trade-off balancing of the following: 1) speed, 2) accuracy and 3) global and semiglobal flexible structure alignment by integration of local substructure matching. The comparison, and matching with competitive accuracy, of one medium sized (250-aa) query structure against the complete PDB database (216,322 protein chains) takes about 8 min using an average desktop computer. The method is at least 2–3 orders of magnitude faster than other tested tools with similar accuracy. We validate the performance of the method for fold and superfamily assignment in a large benchmark set of protein structures. We finally provide a series of examples to illustrate the usefulness of this method and its application in biological discovery.

Conclusions

The method is able to detect partial structure matching, rigid body shifts, conformational changes and tolerates substantial structural variation arising from insertions, deletions and sequence divergence, as well as structural convergence of unrelated proteins.

Electronic supplementary material

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

PASS2 database for the structure-based sequence alignment of distantly related SCOP domain superfamilies: update to version 5 and added features

Structure-based sequence alignment is an essential step in assessing and analysing the relationship of distantly related proteins. PASS2 is a database that records such alignments for protein domain superfamilies and has been constantly updated periodically. This update of the PASS2 version, named as PASS2.5, directly corresponds to the SCOPe 2.04 release. All SCOPe structural domains that share less than 40% sequence identity, as defined by the ASTRAL compendium of protein structures, are included. The current version includes 1977 superfamilies and has been assembled utilizing the structure-based sequence alignment protocol. Such an alignment is obtained initially through MATT, followed by a refinement through the COMPARER program. The JOY program has been used for structural annotations of such alignments. In this update, we have automated the protocol and focused on inclusion of new features such as mapping of GO terms, absolutely conserved residues among the domains in a superfamily and inclusion of PDBs, that are absent in SCOPe 2.04, using the HMM profiles from the alignments of the superfamily members and are provided as a separate list. We have also implemented a more user-friendly manner of data presentation and options for downloading more features. PASS2.5 version is available at http://caps.ncbs.res.in/pass2/.

Database resources of the National Center for Biotechnology Information

The National Center for Biotechnology Information (NCBI) provides a large suite of online resources for biological information and data, including the GenBank(®) nucleic acid sequence database and the PubMed database of citations and abstracts for published life science journals. Additional NCBI resources focus on literature (PubMed Central (PMC), Bookshelf and PubReader), health (ClinVar, dbGaP, dbMHC, the Genetic Testing Registry, HIV-1/Human Protein Interaction Database and MedGen), genomes (BioProject, Assembly, Genome, BioSample, dbSNP, dbVar, Epigenomics, the Map Viewer, Nucleotide, Probe, RefSeq, Sequence Read Archive, the Taxonomy Browser and the Trace Archive), genes (Gene, Gene Expression Omnibus (GEO), HomoloGene, PopSet and UniGene), proteins (Protein, the Conserved Domain Database (CDD), COBALT, Conserved Domain Architecture Retrieval Tool (CDART), the Molecular Modeling Database (MMDB) and Protein Clusters) and chemicals (Biosystems and the PubChem suite of small molecule databases). The Entrez system provides search and retrieval operations for most of these databases. Augmenting many of the web applications are custom implementations of the BLAST program optimized to search specialized datasets. All of these resources can be accessed through the NCBI home page at www.ncbi.nlm.nih.gov.

Structural basis for leucine sensing by the Sestrin2-mTORC1 pathway

Eukaryotic cells coordinate growth with the availability of nutrients through the mechanistic target of rapamycin complex 1 (mTORC1), a master growth regulator. Leucine is of particular importance and activates mTORC1 via the Rag guanosine triphosphatases and their regulators GATOR1 and GATOR2. Sestrin2 interacts with GATOR2 and is a leucine sensor. Here we present the 2.7 angstrom crystal structure of Sestrin2 in complex with leucine. Leucine binds through a single pocket that coordinates its charged functional groups and confers specificity for the hydrophobic side chain. A loop encloses leucine and forms a lid-latch mechanism required for binding. A structure-guided mutation in Sestrin2 that decreases its affinity for leucine leads to a concomitant increase in the leucine concentration required for mTORC1 activation in cells. These results provide a structural mechanism of amino acid sensing by the mTORC1 pathway.

PhyreStorm: A Web Server for Fast Structural Searches Against the PDB

The identification of structurally similar proteins can provide a range of biological insights, and accordingly, the alignment of a query protein to a database of experimentally determined protein structures is a technique commonly used in the fields of structural and evolutionary biology. The PhyreStorm Web server has been designed to provide comprehensive, up-to-date and rapid structural comparisons against the Protein Data Bank (PDB) combined with a rich and intuitive user interface. It is intended that this facility will enable biologists inexpert in bioinformatics access to a powerful tool for exploring protein structure relationships beyond what can be achieved by sequence analysis alone. By partitioning the PDB into similar structures, PhyreStorm is able to quickly discard the majority of structures that cannot possibly align well to a query protein, reducing the number of alignments required by an order of magnitude. PhyreStorm is capable of finding 93 ± 2% of all highly similar (TM-score > 0.7) structures in the PDB for each query structure, usually in less than 60 s. PhyreStorm is available at http://www.sbg.bio.ic.ac.uk/phyrestorm/.

CAMPR3: a database on sequences, structures and signatures of antimicrobial peptides

Antimicrobial peptides (AMPs) are known to have family-specific sequence composition, which can be mined for discovery and design of AMPs. Here, we present CAMP_R3; an update to the existing CAMP database available online at www.camp3.bicnirrh.res.in. It is a database of sequences, structures and family-specific signatures of prokaryotic and eukaryotic AMPs. Family-specific sequence signatures comprising of patterns and Hidden Markov Models were generated for 45 AMP families by analysing 1386 experimentally studied AMPs. These were further used to retrieve AMPs from online sequence databases. More than 4000 AMPs could be identified using these signatures. AMP family signatures provided in CAMP_R3 can thus be used to accelerate and expand the discovery of AMPs. CAMP_R3 presently holds 10247 sequences, 757 structures and 114 family-specific signatures of AMPs. Users can avail the sequence optimization algorithm for rational design of AMPs. The database integrated with tools for AMP sequence and structure analysis will be a valuable resource for family-based studies on AMPs.

Characterization of two cation diffusion facilitators NpunF0707 and NpunF1794 in Nostoc punctiforme

AIMS

To characterize genes involved in maintaining homeostatic levels of zinc in the cyanobacterium Nostoc punctiforme.

METHODS AND RESULTS

Metal efflux transporters play a central role in maintaining homeostatic levels of trace elements such as zinc. Sequence analyses of the N. punctiforme genome identified two potential cation diffusion facilitator (CDF) metal efflux transporters, Npun_F0707 (Cdf31) and Npun_F1794 (Cdf33). Deletion of either Cdf31or Cdf33 resulted in increased zinc retention over 3 h. Interestingly, Cdf31(-) and Cdf33(-) mutants showed no change in sensitivity to zinc exposure in comparison with the wild type, suggesting some compensatory capacity for the loss of each other. Using qRT-PCR, a possible interaction was observed between the two cdf's, where the Cdf31(-) mutant had a more profound effect on cdf33 expression than Cdf33(-) did on cdf31. Over-expression of Cdf31 and Cdf33 in ZntA(-) - and ZitB(-) -deficient Escherichia coli revealed function similarities between the ZntA and ZitB of E. coli and the cyanobacterial transporters.

CONCLUSIONS

The data presented shed light on the function of two important transporters that regulate zinc homeostasis in N. punctiforme.

SIGNIFICANCE AND IMPACT OF THE STUDY

This study shows for the first time the functional characterization of two cyanobacterial zinc efflux proteins belonging to the CDF family.

Structural Perspectives on the Evolutionary Expansion of Unique Protein-Protein Binding Sites

Structures of protein complexes provide atomistic insights into protein interactions. Human proteins represent a quarter of all structures in the Protein Data Bank; however, available protein complexes cover less than 10% of the human proteome. Although it is theoretically possible to infer interactions in human proteins based on structures of homologous protein complexes, it is still unclear to what extent protein interactions and binding sites are conserved, and whether protein complexes from remotely related species can be used to infer interactions and binding sites. We considered biological units of protein complexes and clustered protein-protein binding sites into similarity groups based on their structure and sequence, which allowed us to identify unique binding sites. We showed that the growth rate of the number of unique binding sites in the Protein Data Bank was much slower than the growth rate of the number of structural complexes. Next, we investigated the evolutionary roots of unique binding sites and identified the major phyletic branches with the largest expansion in the number of novel binding sites. We found that many binding sites could be traced to the universal common ancestor of all cellular organisms, whereas relatively few binding sites emerged at the major evolutionary branching points. We analyzed the physicochemical properties of unique binding sites and found that the most ancient sites were the largest in size, involved many salt bridges, and were the most compact and least planar. In contrast, binding sites that appeared more recently in the evolution of eukaryotes were characterized by a larger fraction of polar and aromatic residues, and were less compact and more planar, possibly due to their more transient nature and roles in signaling processes.

Structure of a herpesvirus nuclear egress complex subunit reveals an interaction groove that is essential for viral replication

Herpesviruses require a nuclear egress complex (NEC) for efficient transit of nucleocapsids from the nucleus to the cytoplasm. The NEC orchestrates multiple steps during herpesvirus nuclear egress, including disruption of nuclear lamina and particle budding through the inner nuclear membrane. In the important human pathogen human cytomegalovirus (HCMV), this complex consists of nuclear membrane protein UL50, and nucleoplasmic protein UL53, which is recruited to the nuclear membrane through its interaction with UL50. Here, we present an NMR-determined solution-state structure of the murine CMV homolog of UL50 (M50; residues 1-168) with a strikingly intricate protein fold that is matched by no other known protein folds in its entirety. Using NMR methods, we mapped the interaction of M50 with a highly conserved UL53-derived peptide, corresponding to a segment that is required for heterodimerization. The UL53 peptide binding site mapped onto an M50 surface groove, which harbors a large cavity. Point mutations of UL50 residues corresponding to surface residues in the characterized M50 heterodimerization interface substantially decreased UL50-UL53 binding in vitro, eliminated UL50-UL53 colocalization, prevented disruption of nuclear lamina, and halted productive virus replication in HCMV-infected cells. Our results provide detailed structural information on a key protein-protein interaction involved in nuclear egress and suggest that NEC subunit interactions can be an attractive drug target.

Known structure, unknown function: An inquiry-based undergraduate biochemistry laboratory course

Undergraduate biochemistry laboratory courses often do not provide students with an authentic research experience, particularly when the express purpose of the laboratory is purely instructional. However, an instructional laboratory course that is inquiry- and research-based could simultaneously impart scientific knowledge and foster a student's research expertise and confidence. We have developed a year-long undergraduate biochemistry laboratory curriculum wherein students determine, via experiment and computation, the function of a protein of known three-dimensional structure. The first half of the course is inquiry-based and modular in design; students learn general biochemical techniques while gaining preparation for research experiments in the second semester. Having learned standard biochemical methods in the first semester, students independently pursue their own (original) research projects in the second semester. This new curriculum has yielded an improvement in student performance and confidence as assessed by various metrics. To disseminate teaching resources to students and instructors alike, a freely accessible Biochemistry Laboratory Education resource is available at http://biochemlab.org.

Computational Prediction of acyl-coA Binding Proteins Structure in Brassica napus

Acyl-coA binding proteins could transport acyl-coA esters from plastid to endoplasmic reticulum, prior to fatty acid biosynthesis, leading to the formation of triacylglycerol. The structure and the subcellular localization of acyl-coA binding proteins (ACBP) in Brassica napus were computationally predicted in this study. Earlier, the structure analysis of ACBPs was limited to the small ACBPs, the current study focused on all four classes of ACBPs. Physicochemical parameters including the size and the length, the intron-exon structure, the isoelectric point, the hydrophobicity, and the amino acid composition were studied. Furthermore, identification of conserved residues and conserved domains were carried out. Secondary structure and tertiary structure of ACBPs were also studied. Finally, subcellular localization of ACBPs was predicted. The findings indicated that the physicochemical parameters and subcellular localizations of ACBPs in Brassica napus were identical to Arabidopsis thaliana. Conserved domain analysis indicated that ACBPs contain two or three kelch domains that belong to different families. Identical residues in acyl-coA binding domains corresponded to eight amino acid residues in all ACBPs of B. napus. However, conserved residues of common ACBPs in all species of animal, plant, bacteria and fungi were only inclusive in small ACBPs. Alpha-helixes were displayed and conserved in all the acyl-coA binding domains, representing almost the half of the protein structure. The findings confirm high similarities in ACBPs between A. thaliana and B. napus, they might share the same functions but loss or gain might be possible.