Structure-based characterization of multiprotein complexes

What can AlphaFold do for antimicrobial amyloids?

Amyloids, protein, and peptide assemblies in various organisms are crucial in physiological and pathological processes. Their intricate structures, however, present significant challenges, limiting our understanding of their functions, regulatory mechanisms, and potential applications in biomedicine and technology. This study evaluated the AlphaFold2 ColabFold method's structure predictions for antimicrobial amyloids, using eight antimicrobial peptides (AMPs), including those with experimentally determined structures and AMPs known for their distinct amyloidogenic morphological features. Additionally, two well-known human amyloids, amyloid-β and islet amyloid polypeptide, were included in the analysis due to their disease relevance, short sequences, and antimicrobial properties. Amyloids typically exhibit tightly mated β-strand sheets forming a cross-β configuration. However, certain amphipathic α-helical subunits can also form amyloid fibrils adopting a cross-α structure. Some AMPs in the study exhibited a combination of cross-α and cross-β amyloid fibrils, adding complexity to structure prediction. The results showed that the AlphaFold2 ColabFold models favored α-helical structures in the tested amyloids, successfully predicting the presence of α-helical mated sheets and a hydrophobic core resembling the cross-α configuration. This implies that the AI-based algorithms prefer assemblies of the monomeric state, which was frequently predicted as helical, or capture an α-helical membrane-active form of toxic peptides, which is triggered upon interaction with lipid membranes.

Bdellovibrio bacteriovorus uses chimeric fibre proteins to recognize and invade a broad range of bacterial hosts

Predatory bacteria, like the model endoperiplasmic bacterium Bdellovibrio bacteriovorus, show several adaptations relevant to their requirements for locating, entering and killing other bacteria. The mechanisms underlying prey recognition and handling remain obscure. Here we use complementary genetic, microscopic and structural methods to address this deficit. During invasion, the B. bacteriovorus protein CpoB concentrates into a vesicular compartment that is deposited into the prey periplasm. Proteomic and structural analyses of vesicle contents reveal several fibre-like proteins, which we name the mosaic adhesive trimer (MAT) superfamily, and show localization on the predator surface before prey encounter. These dynamic proteins indicate a variety of binding capabilities, and we confirm that one MAT member shows specificity for surface glycans from a particular prey. Our study shows that the B. bacteriovorus MAT protein repertoire enables a broad means for the recognition and handling of diverse prey epitopes encountered during bacterial predation and invasion.

Crystal structure of domain of unknown function 507 (DUF507) reveals a new protein fold

The crystal structure of the domain of unknown function family 507 protein from Aquifex aeolicus is reported (AaDUF507, UniProt O67633, 183 residues). The structure was determined in two space groups (C2221 and P3221) at 1.9 Å resolution. The phase problem was solved by molecular replacement using an AlphaFold model as the search model. AaDUF507 is a Y-shaped α-helical protein consisting of an anti-parallel 4-helix bundle base and two helical arms that extend 30-Å from the base. The two crystal structures differ by a 25° rigid body rotation of the C-terminal arm. The tertiary structure exhibits pseudo-twofold symmetry. The structural symmetry mirrors internal sequence similarity: residues 11-57 and 102-148 are 30% identical and 53% similar with an E-value of 0.002. In one of the structures, electron density for an unknown ligand, consistent with nicotinamide or similar molecule, may indicate a functional site. Docking calculations suggest potential ligand binding hot spots in the region between the helical arms. Structure-based query of the Protein Data Bank revealed no other protein with a similar tertiary structure, leading us to propose that AaDUF507 represents a new protein fold.

Cryo-electron microscopy and image classification reveal the existence and structure of the coxsackievirus A6 virion

Coxsackievirus A6 (CV-A6) has recently overtaken enterovirus A71 and CV-A16 as the primary causative agent of hand, foot, and mouth disease worldwide. Virions of CV-A6 were not identified in previous structural studies, and it was speculated that the virus is unique among enteroviruses in using altered particles with expanded capsids to infect cells. In contrast, the virions of other enteroviruses are required for infection. Here we used cryo-electron microscopy (cryo-EM) to determine the structures of the CV-A6 virion, altered particle, and empty capsid. We show that the CV-A6 virion has features characteristic of virions of other enteroviruses, including a compact capsid, VP4 attached to the inner capsid surface, and fatty acid-like molecules occupying the hydrophobic pockets in VP1 subunits. Furthermore, we found that in a purified sample of CV-A6, the ratio of infectious units to virions is 1 to 500. Therefore, it is likely that virions of CV-A6 initiate infection, like those of other enteroviruses. Our results provide evidence that future vaccines against CV-A6 should target its virions instead of the antigenically distinct altered particles. Furthermore, the structure of the virion provides the basis for the rational development of capsid-binding inhibitors that block the genome release of CV-A6.

A cryo-EM structure for the three conformers of coxsackievirus A6 provides insight into the infection process of this enterovirus, which is responsible for numerous cases of hand, foot, and mouth disease.

Structure of the ciliogenesis-associated CPLANE complex

Dysfunctional cilia cause pleiotropic human diseases termed ciliopathies. These hereditary maladies are often caused by defects in cilia assembly, a complex event that is regulated by the ciliogenesis and planar polarity effector (CPLANE) proteins Wdpcp, Inturned, and Fuzzy. CPLANE proteins are essential for building the cilium and are mutated in multiple ciliopathies, yet their structure and molecular functions remain elusive. Here, we show that mammalian CPLANE proteins comprise a bona fide complex and report the near-atomic resolution structures of the human Wdpcp-Inturned-Fuzzy complex and of the mouse Wdpcp-Inturned-Fuzzy complex bound to the small guanosine triphosphatase Rsg1. Notably, the crescent-shaped CPLANE complex binds phospholipids such as phosphatidylinositol 3-phosphate via multiple modules and a CPLANE ciliopathy mutant exhibits aberrant lipid binding. Our study provides critical structural and functional insights into an enigmatic ciliogenesis-associated complex as well as unexpected molecular rationales for ciliopathies.

Cellular homologs of the double jelly-roll major capsid proteins clarify the origins of an ancient virus kingdom

Viruses are the most abundant biological entities on Earth and ubiquitous parasites of cellular life forms. The general scenario for the origin of viruses involves evolution from nonviral replicators, such as plasmids and transposons, via recruitment of host proteins for virion formation. One of the most common virion core components, the double jelly-roll major capsid protein of a broad variety of viruses with double-stranded DNA genomes, so far has been thought to represent a virus innovation. However, we present evidence, obtained by protein structure comparison, that this type of virus capsid protein also evolved from a cellular ancestor, a distinct family of carbohydrate-active enzymes. These findings reinforce the chimeric scenario of virus origin.

Viruses are a distinct type of replicators that encode structural proteins encasing virus genomes in virions. For some of the widespread virus capsid proteins and other major components of virions, likely ancestors encoded by cellular life forms are identifiable. In particular, one of the most common capsid proteins, with the single jelly-roll (SJR) fold, appears to have evolved from a particular family of cellular carbohydrate-binding proteins. However, the double jelly-roll major capsid protein (DJR-MCP), the hallmark of the enormously diverse viruses of the kingdom Bamfordvirae within the realm Varidnaviria, which includes bacterial and archaeal icosahedral viruses as well as eukaryotic giant viruses, has been perceived as a virus innovation that evolved by duplication and fusion of the SJR capsid proteins. Here we employ protein structure comparison to show that the DJR fold is represented in several widespread families of cellular proteins, including several groups of carbohydrate-active enzymes. We show that DJR-MCPs share a common ancestry with a distinct family of bacterial DJR proteins (DUF2961) involved in carbohydrate metabolism. Based on this finding, we propose a scenario in which bamfordviruses evolved from nonviral replicators, in particular plasmids, by recruiting a host protein for capsid formation. This sequence of events appears to be the general route of virus origin. The results of this work indicate that virus kingdoms Bamfordvirae, with the DJR-MCPs, and Helvetiavirae that possess two SJR-MCPs, have distinct origins, suggesting a reappraisal of the realm Varidnaviria.

Unlimited Cooperativity of Betatectivirus SSB, a Novel DNA Binding Protein Related to an Atypical Group of SSBs From Protein-Primed Replicating Bacterial Viruses

Bam35 and related betatectiviruses are tail-less bacteriophages that prey on members of the Bacillus cereus group. These temperate viruses replicate their linear genome by a protein-primed mechanism. In this work, we have identified and characterized the product of the viral ORF2 as a single-stranded DNA binding protein (hereafter B35SSB). B35SSB binds ssDNA with great preference over dsDNA or RNA in a sequence-independent, highly cooperative manner that results in a non-specific stimulation of DNA replication. We have also identified several aromatic and basic residues, involved in base-stacking and electrostatic interactions, respectively, that are required for effective protein-ssDNA interaction. Although SSBs are essential for DNA replication in all domains of life as well as many viruses, they are very diverse proteins. However, most SSBs share a common structural domain, named OB-fold. Protein-primed viruses could constitute an exception, as no OB-fold DNA binding protein has been reported. Based on databases searches as well as phylogenetic and structural analyses, we showed that B35SSB belongs to a novel and independent group of SSBs. This group contains proteins encoded by protein-primed viral genomes from unrelated viruses, spanning betatectiviruses and Φ29 and close podoviruses, and they share a conserved pattern of secondary structure. Sensitive searches and structural predictions indicate that B35SSB contains a conserved domain resembling a divergent OB-fold, which would constitute the first occurrence of an OB-fold-like domain in a protein-primed genome.

Rational Design of Resveratrol O-methyltransferase for the Production of Pinostilbene

Pinostilbene is a monomethyl ether analog of the well-known nutraceutical resveratrol. Both compounds have health-promoting properties, but the latter undergoes rapid metabolization and has low bioavailability. O-methylation improves the stability and bioavailability of resveratrol. In plants, these reactions are performed by O-methyltransferases (OMTs). Few efficient OMTs that monomethylate resveratrol to yield pinostilbene have been described so far. Here, we report the engineering of a resveratrol OMT from Vitis vinifera (VvROMT), which has the highest catalytic efficiency in di-methylating resveratrol to yield pterostilbene. In the absence of a crystal structure, we constructed a three-dimensional protein model of VvROMT and identified four critical binding site residues by applying different in silico approaches. We performed point mutations in these positions generating W20A, F24A, F311A, and F318A variants, which greatly reduced resveratrol's enzymatic conversion. Then, we rationally designed eight variants through comparison of the binding site residues with other stilbene OMTs. We successfully modified the native substrate selectivity of VvROMT. Variant L117F/F311W showed the highest conversion to pinostilbene, and variant L117F presented an overall increase in enzymatic activity. Our results suggest that VvROMT has potential for the tailor-made production of stilbenes.

How Coiled-Coil Assemblies Accommodate Multiple Aromatic Residues

Rational protein design requires understanding the contribution of each amino acid to a targeted protein fold. For a subset of protein structures, namely, α-helical coiled coils (CCs), knowledge is sufficiently advanced to allow the rational de novo design of many structures, including entirely new protein folds. Current CC design rules center on using aliphatic hydrophobic residues predominantly to drive the folding and assembly of amphipathic α helices. The consequences of using aromatic residues-which would be useful for introducing structural probes, and binding and catalytic functionalities-into these interfaces are not understood. There are specific examples of designed CCs containing such aromatic residues, e.g., phenylalanine-rich sequences, and the use of polar aromatic residues to make buried hydrogen-bond networks. However, it is not known generally if sequences rich in tyrosine can form CCs, or what CC assemblies these would lead to. Here, we explore tyrosine-rich sequences in a general CC-forming background and resolve new CC structures. In one of these, an antiparallel tetramer, the tyrosine residues are solvent accessible and pack at the interface between the core and the surface. In another more complex structure, the residues are buried and form an extended hydrogen-bond network.

A rotary mechanism for allostery in bacterial hybrid malic enzymes

Bacterial hybrid malic enzymes (MaeB grouping, multidomain) catalyse the transformation of malate to pyruvate, and are a major contributor to cellular reducing power and carbon flux. Distinct from other malic enzyme subtypes, the hybrid enzymes are regulated by acetyl-CoA, a molecular indicator of the metabolic state of the cell. Here we solve the structure of a MaeB protein, which reveals hybrid enzymes use the appended phosphotransacetylase (PTA) domain to form a hexameric sensor that communicates acetyl-CoA occupancy to the malic enzyme active site, 60 Å away. We demonstrate that allostery is governed by a large-scale rearrangement that rotates the catalytic subunits 70° between the two states, identifying MaeB as a new model enzyme for the study of ligand-induced conformational change. Our work provides the mechanistic basis for metabolic control of hybrid malic enzymes, and identifies inhibition-insensitive variants that may find utility in synthetic biology.

Discovering and efficiently promoting the extracellular secretory expression of Thermobacillus sp. ZCTH02-B1 sucrose phosphorylase in Escherichia coli

Sucrose phosphorylase (SPase, EC2.4.1.7) is a promising transglycosylation biocatalyst used for producing glycosylated compounds that are widely used in the food, cosmetics, and pharmaceutical industries. In this study, a recombinant SPase from the Thermobacillus sp. ZCTH02-B1 (rTSPase), which was previously reported to have high thermostability and the catalytic ability to synthesize ascorbic acid 2-glucoside, was attempted to be extracellularly expressed in Escherichia coli BL21(DE3) by fusion of endogenous osmotically-inducible protein Y. Unexpectedly, the rTSPase itself was produced outside the cells with an underestimated performance, although no typical signal peptide was predicted. Further N- and C-terminal truncation experiments revealed that both termini of rTSPase have an important role in protein folding and enzymatic activity, while its secretion was N-terminus associated. Extracellular protein concentration and rTSPase activity achieved 1.8 mg/mL and 6.2 U/mL after induction of 36 h in a 5-L fermenter. High-level extracellular rTSPase production could also be obtained from E. coli within 24 h by inducing overexpression of D, D-carboxypeptidase for cell lysis.

TopMatch-web: pairwise matching of large assemblies of protein and nucleic acid chains in 3D

Frequently, the complete functional units of biological molecules are assemblies of protein and nucleic acid chains. Stunning examples are the complex structures of ribosomes. Here, we present TopMatch-web, a computational tool for the study of the three-dimensional structure, function and evolution of such molecules. The unique feature of TopMatch is its ability to match the protein as well as nucleic acid chains of complete molecular assemblies simultaneously. The resulting structural alignments are visualized instantly using the high-performance molecular viewer NGL. We use the mitochondrial ribosomes of human and yeast as an example to demonstrate the capabilities of TopMatch-web. The service responds immediately, enabling the interactive study of many pairwise alignments of large molecular assemblies in a single session. TopMatch-web is freely accessible at https://topmatch.services.came.sbg.ac.at.

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

Mechanism of conditional partner selectivity in MITF/TFE family transcription factors with a conserved coiled coil stammer motif

Interrupted dimeric coiled coil segments are found in a broad range of proteins and generally confer selective functional properties such as binding to specific ligands. However, there is only one documented case of a basic-helix–loop–helix leucine zipper transcription factor—microphthalmia-associated transcription factor (MITF)—in which an insertion of a three-residue stammer serves as a determinant of conditional partner selectivity. To unravel the molecular principles of this selectivity, we have analyzed the high-resolution structures of stammer-containing MITF and an engineered stammer-less MITF variant, which comprises an uninterrupted symmetric coiled coil. Despite this fundamental difference, both MITF structures reveal identical flanking in-phase coiled coil arrangements, gained by helical over-winding and local asymmetry in wild-type MITF across the stammer region. These conserved structural properties allow the maintenance of a proper functional readout in terms of nuclear localization and binding to specific DNA-response motifs regardless of the presence of the stammer. By contrast, MITF heterodimer formation with other bHLH-Zip transcription factors is only permissive when both factors contain either the same type of inserted stammer or no insert. Our data illustrate a unique principle of conditional partner selectivity within the wide arsenal of transcription factors with specific partner-dependent functional readouts.

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

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

Structural modeling of protein complexes: Current capabilities and challenges

Proteins frequently interact with each other, and the knowledge of structures of the corresponding protein complexes is necessary to understand how they function. Computational methods are increasingly used to provide structural models of protein complexes. Not surprisingly, community-wide Critical Assessment of protein Structure Prediction (CASP) experiments have recently started monitoring the progress in this research area. We participated in CASP13 with the aim to evaluate our current capabilities in modeling of protein complexes and to gain a better understanding of factors that exert the largest impact on these capabilities. To model protein complexes in CASP13, we applied template-based modeling, free docking and hybrid techniques that enabled us to generate models of the topmost quality for 27 of 42 multimers. If templates for protein complexes could be identified, we modeled the structures with reasonable accuracy by straightforward homology modeling. If only partial templates were available, it was nevertheless possible to predict the interaction interfaces correctly or to generate acceptable models for protein complexes by combining template-based modeling with docking. If no templates were available, we used rigid-body docking with limited success. However, in some free docking models, despite the incorrect subunit orientation and missed interface contacts, the approximate location of protein binding sites was identified correctly. Apparently, our overall performance in docking was limited by the quality of monomer models and by the imperfection of scoring methods. The impact of human intervention on our results in modeling of protein complexes was significant indicating the need for improvements of automatic methods.

Evolutionary relationship between the cysteine and histidine rich domains (CHORDs) and Btk-type zinc fingers

Summary

Cysteine and histidine rich domains (CHORDs), implicated in immunity and disease resistance signaling in plants, and in development and signal transduction in muscles and tumorigenesis in animals, are seen to have a cylindrical three-dimensional structure stabilized by the tetrahedral chelation of two zinc ions. CHORDs are regarded as novel zinc-binding domains and classified independently in Pfam and ECOD. Our sequence and structure analysis reveals that both the zinc-binding sites in CHORD possess a zinc ribbon fold and are likely related to each other by duplication and circular permutation. Interestingly, we also detect an evolutionary relationship between each of the CHORD zinc fingers (ZFs) and the Bruton's tyrosine kinase (Btk)-type ZF of the zinc ribbon fold group. Btk_ZF is found in eukaryotic Tec kinase family proteins that are also implicated in signaling pathways in several lineages of hematopoietic cells involved in mammalian immunity. Our analysis suggests that the unique zinc-stabilized fold seen only in the CHORD and Btk_ZFs likely emerged specifically in eukaryotes to mediate diverse signaling pathways.

Supplementary information

Supplementary data are available at Bioinformatics online.

Navigating Among Known Structures in Protein Space

Present-day protein space is the result of 3.7 billion years of evolution, constrained by the underlying physicochemical qualities of the proteins. It is difficult to differentiate between evolutionary traces and effects of physicochemical constraints. Nonetheless, as a rule of thumb, instances of structural reuse, or focusing on structural similarity, are likely attributable to physicochemical constraints, whereas sequence reuse, or focusing on sequence similarity, may be more indicative of evolutionary relationships. Both types of relationships have been studied and can provide meaningful insights to protein biophysics and evolution, which in turn can lead to better algorithms for protein search, annotation, and maybe even design.In broad strokes, studies of protein space vary in the entities they represent, the similarity measure comparing these entities, and the representation used. The entities can be, for example, protein chains, domains, supra-domains, or smaller protein sub-parts denoted themes. The measures of similarity between the entities can be based on sequence, structure, function, or any combination of these. The representation can be global, encompassing the whole space, or local, focusing on a particular region surrounding protein(s) of interest. Global representations include lists of grouped proteins, protein networks, and maps. Networks are the abstraction that is derived most directly from the similarity data: each node is the protein entity (e.g., a domain), and edges connect similar domains. Selecting the entities, the similarity measure, and the abstraction are three intertwined decisions: the similarity measures allow us to identify the entities, and the selection of entities influences what is a meaningful similarity measure. Similarly, we seek entities that are related to each other in a way, for which a simple representation describes their relationships succinctly and accurately. This chapter will cover studies that rely on different entities, similarity measures, and a range of representations to better understand protein structure space. Scholars may use publicly available navigators offering a global representation, and in particular the hierarchical classifications SCOP, CATH, and ECOD, or a local representation, which encompass structural alignment algorithms. Alternatively, scholars can configure their own navigator using existing tools. To demonstrate this DIY (do it yourself) approach for navigating in protein space, we investigate substrate-binding proteins. By presenting sequence similarities among this large and diverse protein family as a network, we can infer that one member (pdb ID 4ntl; of yet unknown function) may bind methionine and suggest a putative binding mechanism.

Maintaining and breaking symmetry in homomeric coiled-coil assemblies

In coiled-coil (CC) protein structures α-helices wrap around one another to form rope-like assemblies. Most natural and designed CCs have two–four helices and cyclic (C_n) or dihedral (D_n) symmetry. Increasingly, CCs with five or more helices are being reported. A subset of these higher-order CCs is of interest as they have accessible central channels that can be functionalised; they are α-helical barrels. These extended cavities are surprising given the drive to maximise buried hydrophobic surfaces during protein folding and assembly in water. Here, we show that α-helical barrels can be maintained by the strategic placement of β-branched aliphatic residues lining the lumen. Otherwise, the structures collapse or adjust to give more-complex multi-helix assemblies without C_n or D_n symmetry. Nonetheless, the structural hallmark of CCs—namely, knobs-into-holes packing of side chains between helices—is maintained leading to classes of CCs hitherto unobserved in nature or accessed by design.

Higher order coiled coils with five or more helices can form α-helical barrels. Here the authors show that placing β-branched aliphatic residues along the lumen yields stable and open α-helical barrels, which is of interest for the rational design of functional proteins; whereas, the absence of β-branched side chains leads to unusual low-symmetry α-helical bundles.

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.

Structural changes of homodimers in the PDB

Protein complexes are involved in various biological phenomena. These complexes are intrinsically flexible, and structural changes are essential to their functions. To perform a large-scale automated analysis of the structural changes of complexes, we combined two original methods. An application, SCPC, compares two structures of protein complexes and decides the match of binding mode. Another application, Motion Tree, identifies rigid-body motions in various sizes and magnitude from the two structural complexes with the same binding mode. This approach was applied to all available homodimers in the Protein Data Bank (PDB). We defined two complex-specific motions: interface motion and subunit-spanning motion. In the former, each subunit of a complex constitutes a rigid body, and the relative movement between subunits occurs at the interface. In the latter, structural parts from distinct subunits constitute a rigid body, providing the relative movement spanning subunits. All structural changes were classified and examined. It was revealed that the complex-specific motions were common in the homodimers, detected in around 40% of families. The dimeric interfaces were likely to be small and flat for interface motion, while large and rugged for subunit-spanning motion. Interface motion was accompanied by a drastic change in contacts at the interface, while the change in the subunit-spanning motion was moderate. These results indicate that the interface properties of homodimers correlated with the type of complex-specific motion. The study demonstrates that the pipeline of SCPC and Motion Tree is useful for the massive analysis of structural change of protein complexes.

Evolutionary analysis of a novel zinc ribbon in the N-terminal region of threonine synthase

Threonine synthase (TS) catalyzes the terminal reaction in the biosynthetic pathway of threonine and requires pyridoxal phosphate as a cofactor. TSs share a common catalytic domain with other fold type II PALP dependent enzymes. TSs are broadly grouped into two classes based on their sequence, quaternary structure, and enzyme regulation. We report the presence of a novel zinc ribbon domain in the N-terminal region preceding the catalytic core in TS. The zinc ribbon domain is present in TSs belonging to both classes. Our sequence analysis reveals that archaeal TSs possess all zinc chelating residues to bind a metal ion that are lacking in the structurally characterized homologs. Phylogenetic analysis suggests that TSs with an N-terminal zinc ribbon likely represents the ancestral state of the enzyme while TSs without a zinc ribbon must have diverged later in specific lineages. The zinc ribbon and its N- and C-terminal extensions are important for enzyme stability, activity and regulation. It is likely that the zinc ribbon domain is involved in higher order oligomerization or mediating interactions with other biomolecules leading to formation of larger metabolic complexes.

Structural insights into the catalytic mechanism of a sacrificial sulfur insertase of the N-type ATP pyrophosphatase family, LarE

The lar operon in Lactobacillus plantarum encodes five Lar proteins (LarA/B/C/D/E) that collaboratively synthesize and incorporate a niacin-derived Ni-containing cofactor into LarA, an Ni-dependent lactate racemase. Previous studies have established that two molecules of LarE catalyze successive thiolation reactions by donating the sulfur atom of their exclusive cysteine residues to the substrate. However, the catalytic mechanism of this very unusual sulfur-sacrificing reaction remains elusive. In this work, we present the crystal structures of LarE in ligand-free and several ligand-bound forms, demonstrating that LarE is a member of the N-type ATP pyrophosphatase (PPase) family with a conserved N-terminal ATP PPase domain and a unique C-terminal domain harboring the putative catalytic site. Structural analysis, combined with structure-guided mutagenesis, leads us to propose a catalytic mechanism that establishes LarE as a paradigm for sulfur transfer through sacrificing its catalytic cysteine residue.

Structural and functional mapping of Rtg2p determinants involved in retrograde signaling and aging of Saccharomyces cerevisiae

In Saccharomyces cerevisiae mitochondrial dysfunction induces retrograde signaling, a pathway of communication from mitochondria to the nucleus that promotes a metabolic remodeling to ensure sufficient biosynthetic precursors for replication. Rtg2p is a positive modulator of this pathway that is also required for cellular longevity. This protein belongs to the ASKHA superfamily, and contains a putative N-terminal ATP-binding domain, but there is no detailed structural and functional map of the residues in this domain that accounts for their contribution to retrograde signaling and aging. Here we use Decomposition of Residue Correlation Networks and site-directed mutagenesis to identify Rtg2p structural determinants of retrograde signaling and longevity. We found that most of the residues involved in retrograde signaling surround the ATP-binding loops, and that Rtg2p N-terminus is divided in three regions whose mutants have different aging phenotypes. We also identified E137, D158 and S163 as possible residues involved in stabilization of ATP at the active site. The mutants shown here may be used to map other Rtg2p activities that crosstalk to other pathways of the cell related to genomic stability and aging.

Classification of the treble clef zinc finger: noteworthy lessons for structure and function evolution

Treble clef (TC) zinc fingers constitute a large fold-group of structural zinc-binding protein domains that mediate numerous cellular functions. We have analysed the sequence, structure, and function relationships among all TCs in the Protein Data Bank. This led to the identification of novel TCs, such as lsr2, YggX and TFIIIC τ 60 kDa subunit, and prediction of a nuclease-like function for the DUF1364 family. The structural malleability of TCs is evident from the many examples with variations to the core structural elements of the fold. We observe domains wherein the structural core of the TC fold is circularly permuted, and also some examples where the overall fold resembles both the TC motif and another unrelated fold. All extant TC families do not share a monophyletic origin, as several TC proteins are known to have been present in the last universal common ancestor and the last eukaryotic common ancestor. We identify several TCs where the zinc-chelating site and residues are not merely responsible for structure stabilization but also perform other functions, such as being redox active in C1B domain of protein kinase C, a nucleophilic acceptor in Ada and catalytic in organomercurial lyase, MerB.

An atypical segment swap in the DN and DC domains of the Acr_tran family resistance-nodulation-cell division pump

Domain/segment swapping is an exchange of equivalent secondary structure element(s) among two or more protein domains resulting in the reconstitution of the original fold while simultaneously causing oligomerization. Here we report an example of the outer membrane factor docking region of the Acr_tran family (PF00873) resistance-nodulation-cell division pump, in which a swapped, misfolded state, of the ferredoxin-like fold of the DN and DC domains, effectuates oligomerization. The atypical segment swap and the associated displacement of a region of the ferredoxin-like fold leads to a topology that is distinct from the original fold. To our knowledge, such segment swaps and associated fold change are rare. This exemplifies the role of functional constraints including oligomerization that determine the interplay between sequence and the three-dimensional structure of proteins.

De Novo modeling in cryo-EM density maps with Pathwalking

As electron cryo-microscopy (cryo-EM) can now frequently achieve near atomic resolution, accurate interpretation of these density maps in terms of atomistic detail has become paramount in deciphering macromolecular structure and function. However, there are few software tools for modeling protein structure from cryo-EM density maps in this resolution range. Here, we present an extension of our original Pathwalking protocol, which can automatically trace a protein backbone directly from a near-atomic resolution (3-6Å) density map. The original Pathwalking approach utilized a Traveling Salesman Problem solver for backbone tracing, but manual adjustment was still required during modeling. In the new version, human intervention is minimized and we provide a more robust approach for backbone modeling. This includes iterative secondary structure identification, termini detection and the ability to model multiple subunits without prior segmentation. Overall, the new Pathwalking procedure provides a more complete and robust tool for annotating protein structure function in near-atomic resolution density maps.

Evolutionary relationships between heme-binding ferredoxin α + β barrels

BACKGROUND

The α + β barrel superfamily of the ferredoxin-like fold consists of a functionally diverse group of evolutionarily related proteins. The barrel architecture of these proteins is formed by either homo-/hetero-dimerization or duplication and fusion of ferredoxin-like domains. Several members of this superfamily bind heme in order to carry out their functions.

RESULTS

We analyze the heme-binding sites in these proteins as well as their barrel topologies. Our comparative structural analysis of these heme-binding barrels reveals two distinct modes of packing of the ferredoxin-like domains to constitute the α + β barrel, which is typified by the Type-1/IsdG-like and Type-2/OxdA-like proteins, respectively. We examine the heme-binding pockets and explore the versatility of the α + β barrels ability to accommodate heme or heme-related moieties, such as siroheme, in at least three different sites, namely, the mode seen in IsdG/OxdA, Cld/DyP/EfeB/HemQ and siroheme decarboxylase barrels.

CONCLUSIONS

Our study offers insights into the plausible evolutionary relationships between the two distinct barrel packing topologies and relate the observed heme-binding sites to these topologies.

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.

Structural and Energetic Characterization of the Ankyrin Repeat Protein Family

Ankyrin repeat containing proteins are one of the most abundant solenoid folds. Usually implicated in specific protein-protein interactions, these proteins are readily amenable for design, with promising biotechnological and biomedical applications. Studying repeat protein families presents technical challenges due to the high sequence divergence among the repeating units. We developed and applied a systematic method to consistently identify and annotate the structural repetitions over the members of the complete Ankyrin Repeat Protein Family, with increased sensitivity over previous studies. We statistically characterized the number of repeats, the folding of the repeat-arrays, their structural variations, insertions and deletions. An energetic analysis of the local frustration patterns reveal the basic features underlying fold stability and its relation to the functional binding regions. We found a strong linear correlation between the conservation of the energetic features in the repeat arrays and their sequence variations, and discuss new insights into the organization and function of these ubiquitous proteins.

Some natural proteins are formed with repetitions of similar amino acid stretches. Ankyrin-repeat proteins constitute one of the most abundant families of this class of proteins that serve as model systems to analyze how variations in sequences exert effects in structures and biological functions. We present an in-depth analysis of the ankyrin repeat protein family, characterizing the variations in the repeating arrays both at the structural and energetic level. We introduce a consistent annotation for the repeat characteristics and describe how the structural differences are related to the sequences by their underlying energetic signatures.

Structural principles for computational and de novo design of 4Fe-4S metalloproteins

Iron-sulfur centers in metalloproteins can access multiple oxidation states over a broad range of potentials, allowing them to participate in a variety of electron transfer reactions and serving as catalysts for high-energy redox processes. The nitrogenase FeMoCO cluster converts di-nitrogen to ammonia in an eight-electron transfer step. The 2(Fe4S4) containing bacterial ferredoxin is an evolutionarily ancient metalloprotein fold and is thought to be a primordial progenitor of extant oxidoreductases. Controlling chemical transformations mediated by iron-sulfur centers such as nitrogen fixation, hydrogen production as well as electron transfer reactions involved in photosynthesis are of tremendous importance for sustainable chemistry and energy production initiatives. As such, there is significant interest in the design of iron-sulfur proteins as minimal models to gain fundamental understanding of complex natural systems and as lead-molecules for industrial and energy applications. Herein, we discuss salient structural characteristics of natural iron-sulfur proteins and how they guide principles for design. Model structures of past designs are analyzed in the context of these principles and potential directions for enhanced designs are presented, and new areas of iron-sulfur protein design are proposed. This article is part of a Special issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, protein networks, edited by Ronald L. Koder and J.L Ross Anderson.

The UBR-box and its relationship to binuclear RING-like treble clef zinc fingers

Background

The N-end rule pathway is a part of the ubiquitin–dependent proteolytic system wherein N-recognin proteins recognize the amino terminal degradation signals (N-degrons) of the substrate. The type 1 N-degron recognizing UBR-box domain of the eukaryotic Arg/N-end rule pathway is known to possess a novel three-zinc-stabilized heart-shaped fold.

Results

Using sequence and structure analysis we argue that the UBR-box fold emerged from a binuclear RING-like treble clef zinc finger. The RING-like core is preserved in the UBR-box and the metal-chelating motifs display significant sequence and structural similarity to B-box and ZZ domains. UBR-box domains retrieved in our analysis co-occur with a variety of other protein domains, suggestive of its involvement in diverse biological roles. The UBR-box is a unique family of RING-like treble clefs as it displays a distinct circular permutation at the zinc-knuckle of the first zinc-binding site unlike other documented permutations of the RING-like domains which occur at the second zinc-binding site. The circular permutation of the RING-like treble clef scaffold has possibly aided the gain of a novel and relatively deep cleft suited for binding N-degrons. The N- and C-terminal extensions to the circularly permuted RING-like region bind a third zinc ion, which likely provides additional stability to the domain by keeping the two halves of the permuted zinc-knuckle together.

Conclusions

Structural modifications and extensions to the RING-like core have resulted in a novel UBR-box fold, which can recognize and target the type 1 N-degron containing proteins for ubiquitin-mediated proteolysis. The UBR-box appears to have emerged during the expansion of ubiquitin system pathway-related functions in eukaryotes, but is also likely to have other non-N-recognin functions as well.

Reviewers

This article was reviewed by Eugene Koonin, Balaji Santhanam, Kira S. Makarova.

Electronic supplementary material

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

Linking structural features of protein complexes and biological function

Protein-protein interaction (PPI) establishes the central basis for complex cellular networks in a biological cell. Association of proteins with other proteins occurs at varying affinities, yet with a high degree of specificity. PPIs lead to diverse functionality such as catalysis, regulation, signaling, immunity, and inhibition, playing a crucial role in functional genomics. The molecular principle of such interactions is often elusive in nature. Therefore, a comprehensive analysis of known protein complexes from the Protein Data Bank (PDB) is essential for the characterization of structural interface features to determine structure-function relationship. Thus, we analyzed a nonredundant dataset of 278 heterodimer protein complexes, categorized into major functional classes, for distinguishing features. Interestingly, our analysis has identified five key features (interface area, interface polar residue abundance, hydrogen bonds, solvation free energy gain from interface formation, and binding energy) that are discriminatory among the functional classes using Kruskal-Wallis rank sum test. Significant correlations between these PPI interface features amongst functional categories are also documented. Salt bridges correlate with interface area in regulator-inhibitors (r = 0.75). These representative features have implications for the prediction of potential function of novel protein complexes. The results provide molecular insights for better understanding of PPIs and their relation to biological functions.

BCSearch: fast structural fragment mining over large collections of protein structures

Resources to mine the large amount of protein structures available today are necessary to better understand how amino acid variations are compatible with conformation preservation, to assist protein design, engineering and, further, the development of biologic therapeutic compounds. BCSearch is a versatile service to efficiently mine large collections of protein structures. It relies on a new approach based on a Binet-Cauchy kernel that is more discriminative than the widely used root mean square deviation criterion. It has statistics independent of size even for short fragments, and is fast. The systematic mining of large collections of structures such as the complete SCOPe protein structural classification or comprehensive subsets of the Protein Data Bank can be performed in few minutes. Based on this new score, we propose four innovative applications: BCFragSearch and BCMirrorSearch, respectively, search for fragments similar and anti-similar to a query and return information on the diversity of the sequences of the hits. BCLoopSearch identifies candidate fragments of fixed size matching the flanks of a gaped structure. BCSpecificitySearch analyzes a complete protein structure and returns information about sites having few similar fragments. BCSearch is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/BCSearch.

Taking structure searches to the next dimension

Structure comparisons are now the first step when a new experimental high-resolution protein structure has been determined. In this issue of Structure, Wiederstein and colleagues describe their latest tool for comparing structures, which gives us the unprecedented power to discover crucial structural connections between whole complexes of proteins in the full structural database in real time.

The structure of the cyanobactin domain of unknown function from PatG in the patellamide gene cluster

Patellamides are members of the cyanobactin family of ribosomally synthesized and post-translationally modified cyclic peptide natural products, many of which, including some patellamides, are biologically active. A detailed mechanistic understanding of the biosynthetic pathway would enable the construction of a biotechnological `toolkit' to make novel analogues of patellamides that are not found in nature. All but two of the protein domains involved in patellamide biosynthesis have been characterized. The two domains of unknown function (DUFs) are homologous to each other and are found at the C-termini of the multi-domain proteins PatA and PatG. The domain sequence is found in all cyanobactin-biosynthetic pathways characterized to date, implying a functional role in cyanobactin biosynthesis. Here, the crystal structure of the PatG DUF domain is reported and its binding interactions with plausible substrates are investigated.