Search

Clonal propagation of plants by induction of adventitious roots (ARs) from stem cuttings is a requisite step in breeding programs. Nevertheless, a major barrier exists for propagating valuable plants that naturally have low capacity to form ARs. Due to the central role of auxin in organogenesis, indole-3-butyric acid (IBA) is often utilized, yet many recalcitrant plants do not form ARs in response to such treatment. We describe the synthesis and screening of a focused library of synthetic auxin conjugates inEucalyptus grandiscuttings, highlighting 4-chlorophenoxyacetic acid-L-tryptophan-OMe as a competent enhancer of adventitious rooting in a number of recalcitrant woody plants. Comprehensive metabolic and functional analyses revealed that this activity is engendered by prolonged auxin signaling due to initial fast uptake and slow release and clearance of the free auxin 4-chlorophenoxyacetic acid. This work highlights the utility of a slow-release strategy for bioactive compounds and provides an exemplar for further rational development of more effective plant-growth regulators for agriculture.

The emergence of novel proteins, beyond these that can be readily made by duplication and recombination of preexisting domains, is elusive. De novo emergence from random sequences is unlikely because the vast majority of random chains would not even fold, let alone function. An alternative explanation is that novel proteins emerge by duplication and fusion of pre-existing polypeptide segments. In this case, traces of such ancient events may remain within contemporary proteins in the form of reused segments. Together with the late Dan Tawfik, we detected such similar segments, far shorter than intact protein domains, which are found in different environments. The detection of these, "bridging themes," was based on a unique search strategy, where in addition to searching for similarity of shared fragments, so-called "themes," we also explicitly searched for cases in which the sequence segments before and after the theme are dissimilar (both in sequence and structure). Here, using a similar strategy, we further expanded the search and discovered almost 500 additional "bridging themes," linking domains that are often from ancient folds. The themes, of 20 residues or more (average 53), do not retain their structure despite sharing 37% sequence identity on average. Indeed, conformation flexibility may confer an evolutionary advantage, in that it fits in multiple environments. We elaborate on two interesting themes, shared between Rossmann/Trefoil-Plexin-like domains and a β-propeller-like domain. FOR A BROAD AUDIENCE: A fundamental question in molecular evolution is how protein domains emerged. Similar segments shared between domains of seemingly distinct origins, may offer clues, as these may be remnants of the evolutionary process through which these domains emerged. However, finding such cases is difficult. Here, we expand the set of such cases which we curated previously, adding segments shared between domains that are considered ancient.

Accurate protein structure predictors use clusters of homologues, which disregard sequence specific effects. In this issue of Structure, Weißenow and colleagues report a deep learning-based tool, EMBER2, that efficiently predicts the distances in a protein structure from its amino acid sequence only. This approach should enable the analysis of mutation effects.

Currently known G protein-coupled receptors (GPCRs) have a single transmembrane domain. Many GPCRs form dimers that have two transmembrane domains (one per protein), and there are indications that this dimeric interaction is functionally meaningful. Here, based on sequence analysis and structure predictions, we report the existence of 57 proteins with two, three, or four GPCR domains within the same protein chain. We analyze the structures of these multi-GPCRs and show that almost all have DRY/NPxxY motifs, a strong indication of signaling activity. By homology, most of the multi-GPCRs that we identified are olfactory-related; a few are chemokine-related. Multi-GPCR candidates are found in various Chordata species including fish, camel, marmite, Chinese hamster, and new world monkeys. The discovery of receptors with multiple transmembrane domains suggests the possibility for signal regulation and amplification within an individual receptor, revealing another step in GPCR evolution and a new layer of complexity in signal transduction.

Influenza A virus is the most virulent influenza subtype and is associated with large-scale global pandemics characterized by high levels of morbidity and mortality. Developing simple and sensitive molecular methods for detecting influenza viruses is critical. Neuraminidase, an exo-glycosidase displayed on the surface of influenza virions, is responsible for the release of the virions and their spread in the infected host. Here, we present a new phenoxy-dioxetane chemiluminescent probe (CLNA) that can directly detect neuraminidase activity. The probe exhibits an effective turn-on response upon reaction with neuraminidase and produces a strong emission signal at 515 nm with an extremely high signal-to-noise ratio. Comparison measurements of our new probe with previously reported analogous neuraminidase optical probes showed superior detection capability in terms of response time and sensitivity. Thus, as far as we know, our probe is the most sensitive neuraminidase probe known to date. The chemiluminescence turn-on response produced by our neuraminidase probe enables rapid screening for small molecules that inhibit viral replication through different mechanisms as validated directly in influenza A-infected mammalian cells using the known inhibitors oseltamivir and amantadine. We expect that our new chemiluminescent neuraminidase probe will prove useful for various applications requiring neuraminidase detection including drug discovery assays against various influenza virus strains in mammalian cells.

The human Na+/H+ antiporter NHA2 (SLC9B2) transports Na+ or Li+ across the plasma membrane in exchange for protons, and is implicated in various pathologies. It is a 537 amino acids protein with an 82 residues long hydrophilic cytoplasmic N-terminus followed by a transmembrane part comprising 14 transmembrane helices. We optimized the functional expression of HsNHA2 in the plasma membrane of a salt-sensitive Saccharomyces cerevisiae strain and characterized a set of mutated or truncated versions of HsNHA2 in terms of their substrate specificity, transport activity, localization and protein stability. We identified a highly conserved proline 246, located in the core of the protein, as being crucial for ion selectivity. The replacement of P246 with serine or threonine resulted in antiporters with altered substrate specificity and increased resistance to the HsNHA2-specific inhibitor phloretin that were not only highly active at an acidic pH of 4.0 (like the native antiporter), but also at neutral pH. We also experimentally confirmed the importance of a putative salt bridge between E215 and R432 for antiporter function and structural integrity. Truncations of the first 50 - 70 residues of the N-terminus doubled the transport activity of HsNHA2, whilst changes in the charge at positions E47, E56, K57, or K58 decreased the antiporter’s transport activity. Thus, the hydrophilic N-terminal part of the protein appears to allosterically autoinhibit its cation transport. Our data also show this in vivo approach to be useful for a rapid screening of SNP’s effect on HsNHA2 activity.

Significance Membrane-embedded cation/proton antiporters (CPAs) exchange monovalent cations with protons by alternating between two extreme states: inward- and outward-facing. Although atomic-resolution structures are available, there is still an ongoing debate regarding these antiporters’ transport mechanism. Here, we use computer simulations to explore the dynamics of two bacterial antiporters along two axes of motion: a tilting movement and a vertical translation of the antiporters’ two functional domains. By analyzing these dynamic changes, rather than comparing the two extreme states that CPAs adopt, we determine that it is the tilting movement of the antiporters’ mobile domain that drives the relevant conformational changes. Finally, by applying the knowledge gained from these simulations, we predict an unknown outward-facing conformation and verify it experimentally., There is ongoing debate regarding the mechanism through which cation/proton antiporters (CPAs), like Thermus thermophilus NapA (TtNapA) and Escherichia coli NapA (EcNhaA), alternate between their outward- and inward-facing conformations in the membrane. CPAs comprise two domains, and it is unclear whether the transition is driven by their rocking-bundle or elevator motion with respect to each other. Here we address this question using metadynamics simulations of TtNapA, where we bias conformational sampling along two axes characterizing the two proposed mechanisms: angular and translational motions, respectively. By applying the bias potential for the two axes simultaneously, as well as to the angular, but not the translational, axis alone, we manage to reproduce each of the two known states of TtNapA when starting from the opposite state, in support of the rocking-bundle mechanism as the driver of conformational change. Next, starting from the inward-facing conformation of EcNhaA, we sample what could be its long-sought-after outward-facing conformation and verify it using cross-linking experiments.

Background Coronavirus disease 2019 (COVID-19) and dengue fever are difficult to distinguish given shared clinical and laboratory features. Failing to consider COVID-19 due to false-positive dengue serology can have serious implications. We aimed to assess this possible cross-reactivity. Methods We analyzed clinical data and serum samples from 55 individuals with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. To assess dengue serology status, we used dengue-specific antibodies by means of lateral-flow rapid test, as well as enzyme-linked immunosorbent assay (ELISA). Additionally, we tested SARS-CoV-2 serology status in patients with dengue and performed in-silico protein structural analysis to identify epitope similarities. Results Using the dengue lateral-flow rapid test we detected 12 positive cases out of the 55 (21.8%) COVID-19 patients versus zero positive cases in a control group of 70 healthy individuals (P = 2.5E−5). This includes 9 cases of positive immunoglobulin M (IgM), 2 cases of positive immunoglobulin G (IgG), and 1 case of positive IgM as well as IgG antibodies. ELISA testing for dengue was positive in 2 additional subjects using envelope protein directed antibodies. Out of 95 samples obtained from patients diagnosed with dengue before September 2019, SARS-CoV-2 serology targeting the S protein was positive/equivocal in 21 (22%) (16 IgA, 5 IgG) versus 4 positives/equivocal in 102 controls (4%) (P = 1.6E−4). Subsequent in-silico analysis revealed possible similarities between SARS-CoV-2 epitopes in the HR2 domain of the spike protein and the dengue envelope protein. Conclusions Our findings support possible cross-reactivity between dengue virus and SARS-CoV-2, which can lead to false-positive dengue serology among COVID-19 patients and vice versa. This can have serious consequences for both patient care and public health.

Sulfur is essential for biological processes such as amino acid biogenesis, iron-sulfur cluster formation, and redox homeostasis. To acquire sulfur-containing compounds from the environment, bacteria have evolved high-affinity uptake systems, predominant among which is the ABC transporter family. Theses membrane-embedded enzymes use the energy of ATP hydrolysis for transmembrane transport of a wide range of biomolecules against concentration gradients. Three distinct bacterial ABC import systems of sulfur-containing compounds have been identified, but the molecular details of their transport mechanism remain poorly characterized. Here we provide results from a biochemical analysis of the purified Escherichia coli YecSC-FliY cysteine/cystine import system. We found that the substrate-binding protein FliY binds l-cystine, l-cysteine, and d-cysteine with micromolar affinities. However, binding of the l- and d-enantiomers induced different conformational changes of FliY, where the l- enantiomer-substrate-binding protein complex interacted more efficiently with the YecSC transporter. YecSC had low basal ATPase activity that was moderately stimulated by apo FliY, more strongly by d-cysteine-bound FliY, and maximally by l-cysteine- or l-cystine-bound FliY. However, at high FliY concentrations, YecSC reached maximal ATPase rates independent of the presence or nature of the substrate. These results suggest that FliY exists in a conformational equilibrium between an open, unliganded form that does not bind to the YecSC transporter and closed, unliganded and closed, liganded forms that bind this transporter with variable affinities but equally stimulate its ATPase activity. These findings differ from previous observations for similar ABC transporters, highlighting the extent of mechanistic diversity in this large protein family.

Significance How proteins evolved to recognize and bind their ligands is a key mystery in protein function evolution. To explore this mystery, we study how proteins bind adenine, an ancient fragment. We characterize physicochemical patterns of protein–adenine interactions and link these to proteins’ evolutionary origins. In conflict with previous findings, we see that all of adenine’s hydrogen donors and acceptors have been used to bind proteins, and that adenine binding is likely to have emerged multiple times in evolution. To identify adenine-binding sites of shared origin, we use “themes”: short amino acid segments suggested to constitute evolutionary building blocks. We detect specific themes that are engaged in adenine binding; the detection of these in a protein’s sequence might reveal its function., Proteins’ interactions with ancient ligands may reveal how molecular recognition emerged and evolved. We explore how proteins recognize adenine: a planar rigid fragment found in the most common and ancient ligands. We have developed a computational pipeline that extracts protein–adenine complexes from the Protein Data Bank, structurally superimposes their adenine fragments, and detects the hydrogen bonds mediating the interaction. Our analysis extends the known motifs of protein–adenine interactions in the Watson–Crick edge of adenine and shows that all of adenine’s edges may contribute to molecular recognition. We further show that, on the proteins' side, binding is often mediated by specific amino acid segments (“themes”) that recur across different proteins, such that different proteins use the same themes when binding the same adenine-containing ligands. We identify numerous proteins that feature these themes and are thus likely to bind adenine-containing ligands. Our analysis suggests that adenine binding has emerged multiple times in evolution.

The Escherichia coli yobA–yebZ–yebY (AZY) operon encodes the proteins YobA, YebZ, and YebY. YobA and YebZ are homologs of the CopC periplasmic copper-binding protein and the CopD putative copper importer, respectively, whereas YebY belongs to the uncharacterized Domain of Unknown Function 2511 family. Despite numerous studies of E. coli copper homeostasis and the existence of the AZY operon in a range of bacteria, the operon's proteins and their functional roles have not been explored. In this study, we present the first biochemical and functional studies of the AZY proteins. Biochemical characterization and structural modeling indicate that YobA binds a single Cu2+ ion with high affinity. Bioinformatics analysis shows that YebY is widespread and encoded either in AZY operons or in other genetic contexts unrelated to copper homeostasis. We also determined the 1.8 Å resolution crystal structure of E. coli YebY, which closely resembles that of the lantibiotic self-resistance protein MlbQ. Two strictly conserved cysteine residues form a disulfide bond, consistent with the observed periplasmic localization of YebY. Upon treatment with reductants, YebY binds Cu+ and Cu2+ with low affinity, as demonstrated by metal-binding analysis and tryptophan fluorescence. Finally, genetic manipulations show that the AZY operon is not involved in copper tolerance or antioxidant defense. Instead, YebY and YobA are required for the activity of the copper-related NADH dehydrogenase II. These results are consistent with a potential role of the AZY operon in copper delivery to membrane proteins.

The ConSurf web server (https://consurf.tau.ac.il/) for using evolutionary data to detect functional regions is useful for analyzing proteins. The analysis is based on the premise that functional regions, which may for example facilitate ligand binding and catalysis, often evolve slowly. The analysis requires finding enough effective, i.e., non-redundant, sufficiently remote homologs. Indeed, the ConSurf pipeline, which is based on state-of-the-art protein sequence databases and analysis tools, is highly valuable for protein analysis. ConSurf also allows evolutionary analysis of RNA, but the analysis often fails due to insufficient data, particularly the inability of the current pipeline to detect enough effective RNA homologs. This is because the RNA search tools and databases offered are not as good as those used for protein analysis. Fortunately, ConSurf also allows importing external collections of homologs in the form of a multiple sequence alignment (MSA). Leveraging this, here we describe various protocols for constructing MSAs for successful ConSurf analysis of RNA queries. We report the level of success of these protocols on an exemplary set comprising a dozen RNA molecules of diverse structure and function. © 2021 Wiley Periodicals LLC. Basic Protocol 1: Standard ConSurf evolutionary conservation analysis of an RNA query. Basic Protocol 2: ConSurf evolutionary conservation analysis of an RNA query with external MSA. Support Protocol 1: Construction of an MSA for an RNA query using other online servers. Support Protocol 2: Construction of an MSA for an RNA query using nHMMER locally.

The Sigma-1 receptor (S1R) is a transmembrane protein with important roles in cellular homeostasis in normal physiology and in disease. Especially in neurodegenerative diseases, S1R activation has been shown to provide neuroprotection by modulating calcium signaling, mitochondrial function and reducing endoplasmic reticulum (ER) stress. S1R missense mutations are one of the causes of the neurodegenerative Amyotrophic Lateral Sclerosis and distal hereditary motor neuronopathies. Although the S1R has been studied intensively, basic aspects remain controversial, such as S1R topology and whether it reaches the plasma membrane. To address these questions, we have undertaken several approaches. C-terminal tagging with a small biotin-acceptor peptide and BirA biotinylation in cells suggested a type II membrane orientation (cytosolic N-terminus). However, N-terminal tagging gave an equal probability for both possible orientations. This might explain conflicting reports in the literature, as tags may affect the protein topology. Therefore, we studied untagged S1R using a protease protection assay and a glycosylation mapping approach, introducing N-glycosylation sites. Both methods provided unambiguous results showing that the S1R is a type II membrane protein with a short cytosolic N-terminal tail. Assessments of glycan processing, surface fluorescence-activated cell sorting, and cell surface biotinylation indicated ER retention, with insignificant exit to the plasma membrane, in the absence or presence of S1R agonists or of ER stress. These findings may have important implications for S1R-based therapeutic approaches.

Significance All currently known architectures of outer-membrane beta barrels (OMBBs) have only one barrel. While the vast majority function as oligomers, with barrels from different chains packing against each other in the membrane, it was assumed that these multiple chains are needed to form multibarrel structures. And yet, here we show that multibarrel chains exist. Using state-of-the-art sequence and structure analysis tools, we report the discovery of more than 30 multibarrel architectures from gram-negative bacteria. The discovery of these architectures reveals another interesting chapter in OMBB evolution and has implications for protein engineering. The evolutionary advantages of multibarrels are yet to be discovered., Outer-membrane beta barrels (OMBBs) are found in the outer membrane of gram-negative bacteria and eukaryotic organelles. OMBBs fold as antiparallel β-sheets that close onto themselves, forming pores that traverse the membrane. Currently known structures include only one barrel, of 8 to 36 strands, per chain. The lack of multi-OMBB chains is surprising, as most OMBBs form oligomers, and some function only in this state. Using a combination of sensitive sequence comparison methods and coevolutionary analysis tools, we identify many proteins combining multiple beta barrels within a single chain; combinations that include eight-stranded barrels prevail. These multibarrels seem to be the result of independent, lineage-specific fusion and amplification events. The absence of multibarrels that are universally conserved in bacteria with an outer membrane, coupled with their frequent de novo genesis, suggests that their functions are not essential but rather beneficial in specific environments. Adjacent barrels of complementary function within the same chain may allow for functions beyond those of the individual barrels.

Patterns observed by examining the evolutionary relationships among proteins of common origin can reveal the structural and functional importance of specific residue positions. In particular, amino acids that are highly conserved (i.e., their positions evolve at a slower rate than other positions) are particularly likely to be of biological importance, for example, for ligand binding. ConSurf is a bioinformatics tool for accurately estimating the evolutionary rate of each position in a protein family. Here we introduce a new release of ConSurf‐DB, a database of precalculated ConSurf evolutionary conservation profiles for proteins of known structure. ConSurf‐DB provides high‐accuracy estimates of the evolutionary rates of the amino acids in each protein. A reliable estimate of a query protein's evolutionary rates depends on having a sufficiently large number of effective homologues (i.e., nonredundant yet sufficiently similar). With current sequence data, ConSurf‐DB covers 82% of the PDB proteins. It will be updated on a regular basis to ensure that coverage remains high—and that it might even increase. Much effort was dedicated to improving the user experience. The repository is available at https://consurfdb.tau.ac.il/. Broader audience By comparing a protein to other proteins of similar origin, it is possible to determine the extent to which each amino acid position in the protein evolved slowly or rapidly. A protein's evolutionary profile can provide valuable insights: For example, amino acid positions that are highly conserved (i.e., evolved slowly) are particularly likely to be of structural and/or functional importance, for example, for ligand binding and catalysis. We introduce here a new and improved version of ConSurf‐DB, a continually updated database that provides precalculated evolutionary profiles of proteins with known structure.

Measuring evolutionary rates at the residue level is indispensable for gaining structural and functional insights into proteins. State-of-the-art tools for estimating rates take as input a large set of homologous proteins, a probabilistic model of evolution and a phylogenetic tree. However, a gap exists when only few or no homologous proteins can be found, e.g., orphan proteins. In addition, such tools do not take the three-dimensional (3D) structure of the protein into account. The association between the 3D structure and site-specific rates can be learned using machine-learning regression tools from a cohort of proteins for which both the structure and a large set of homologs exist. Here we present EvoRator, a user-friendly web server that implements a machine-learning regression algorithm to predict site-specific evolutionary rates from protein structures. We show that EvoRator outperforms predictions obtained using traditional physicochemical features, such as relative solvent accessibility and weighted contact number. We also demonstrate the application of EvoRator in three common scenarios that arise in protein evolution research: (1) orphan proteins for which no (or few) homologs exist; (2) When homologous sequences exist, our algorithm contrasts structure-based estimates of the evolutionary rates and the phylogeny-based estimates. This allows detecting sites that are likely conserved due to functional rather than structural constraints; (3) Algorithms that only rely on homologous sequence often fail to accurately measure the evolutionary rates of positions in gapped sequence alignments, which frequently occurs as a result of a clade-specific insertion. Our algorithm makes use of training data and known 3D structure of such gapped positions to predict their evolutionary rates. EvoRator is freely available for all users at: https://evorator.tau.ac.il/.