The ultra-high affinity transport proteins of ubiquitous marine bacteria.

SAR11 bacteria are the most abundant microorganisms in the surface ocean1 and have global biogeochemical importance2-4. To thrive in their competitive oligotrophic environment, these bacteria rely heavily on solute-binding proteins that facilitate uptake of specific substrates via membrane transporters5,6. The functions and properties of these transport proteins are key factors in the assimilation of dissolved organic matter and biogeochemical cycling of nutrients in the ocean, but they have remained largely inaccessible to experimental investigation. Here we performed genome-wide experimental characterization of all solute-binding proteins in a prototypical SAR11 bacterium, revealing specific functions and general trends in their properties that contribute to the success of SAR11 bacteria in oligotrophic environments. We found that the solute-binding proteins of SAR11 bacteria have extremely high binding affinity (dissociation constant >20 pM) and high binding specificity, revealing molecular mechanisms of oligotrophic adaptation. Our functional data have uncovered new carbon sources for SAR11 bacteria and enable accurate biogeographical analysis of SAR11 substrate uptake capabilities throughout the ocean. This study provides a comprehensive view of the substrate uptake capabilities of ubiquitous marine bacteria, providing a necessary foundation for understanding their contribution to assimilation of dissolved organic matter in marine ecosystems.

Insertions and deletions mediated functional divergence of Rossmann fold enzymes.

Nucleobase-containing coenzymes are hypothesized to be relics of an early RNA-based world that preceded the emergence of proteins. Despite the importance of coenzyme-protein synergisms, their emergence and evolution remain understudied. An excellent target to address this issue is the Rossmann fold, the most catalytically diverse and abundant protein architecture in nature. We investigated two main Rossmann lineages: the nicotinamide adenine dinucleotide phosphate (NAD(P)) and the S-adenosyl methionine (SAM)- binding superfamilies. To identify the evolutionary changes that lead to a coenzyme specificity switch on these superfamilies, we performed structural and sequence-based Hidden Markov model analysis to systematically search for key motifs in their coenzyme-binding pockets. Our analyses revealed that through insertions and deletions (InDels) and a residue substitution, the ancient ß1-loop-α1 coenzyme-binding structure of NAD(P) could be reshaped into the SAM-binding ß1-loop-α1 structure. To experimentally prove this obsevation, we removed three amino acids from the NAD(P)-binding pocket and solved the structure of the resulting mutant, revealing the characteristic loop features of the SAM-binding pocket. To confirm the binding to SAM, we performed isothermal titration calorimetry measurements. Molecular dynamics simulations also corroborated the role of InDels in abolishing NAD binding and acquiring SAM binding. Our results uncovered how nature may have utilized insertions and deletions to optimize the different coenzyme-binding pockets and the distinct functionalities observed for Rossmann superfamilies. This work also proposes a general mechanism by which protein templates could have been recycled through the course of evolution to adopt different coenzymes and confer distinct chemistries.

Chemotactic Interactions Drive Migration of Membraneless Active Droplets.

In nature, chemotactic interactions are ubiquitous and play a critical role in driving the collective behavior of living organisms. Reproducing these interactions in vitro is still a paramount challenge due to the complexity of mimicking and controlling cellular features, such as tangled metabolic networks, cytosolic macromolecular crowding, and cellular migration, on a microorganism size scale. Here, we generate enzymatically active cell-sized droplets able to move freely, and by following a chemical gradient, able to interact with the surrounding droplets in a collective manner. The enzyme within the droplets generates a pH gradient that extends outside the edge of the droplets. We discovered that the external pH gradient triggers droplet migration and controls its directionality, which is selectively toward the neighboring droplets. Hence, by changing the enzyme activity inside the droplet, we tuned the droplet migration speed. Furthermore, we showed that these cellular-like features can facilitate the reconstitution of a simple and linear protometabolic pathway and increase the final reaction product generation. Our work suggests that simple and stable membraneless droplets can reproduce complex biological phenomena, opening new perspectives as bioinspired materials and synthetic biology tools.

Automated Structure- and Sequence-Based Design of Proteins for High Bacterial Expression and Stability.

Upon heterologous overexpression, many proteins misfold or aggregate, thus resulting in low functional yields. Human acetylcholinesterase (hAChE), an enzyme mediating synaptic transmission, is a typical case of a human protein that necessitates mammalian systems to obtain functional expression. We developed a computational strategy and designed an AChE variant bearing 51 mutations that improved core packing, surface polarity, and backbone rigidity. This variant expressed at ∼2,000-fold higher levels in E. coli compared to wild-type hAChE and exhibited 20°C higher thermostability with no change in enzymatic properties or in the active-site configuration as determined by crystallography. To demonstrate broad utility, we similarly designed four other human and bacterial proteins. Testing at most three designs per protein, we obtained enhanced stability and/or higher yields of soluble and active protein in E. coli. Our algorithm requires only a 3D structure and several dozen sequences of naturally occurring homologs, and is available at http://pross.weizmann.ac.il.

Efficient Exploration of Sequence Space by Sequence-Guided Protein Engineering and Design.

The rapid growth of sequence databases over the past two decades means that protein engineers faced with optimizing a protein for any given task will often have immediate access to a vast number of related protein sequences. These sequences encode information about the evolutionary history of the protein and the underlying sequence requirements to produce folded, stable, and functional protein variants. Methods that can take advantage of this information are an increasingly important part of the protein engineering tool kit. In this Perspective, we discuss the utility of sequence data in protein engineering and design, focusing on recent advances in three main areas: the use of ancestral sequence reconstruction as an engineering tool to generate thermostable and multifunctional proteins, the use of sequence data to guide engineering of multipoint mutants by structure-based computational protein design, and the use of unlabeled sequence data for unsupervised and semisupervised machine learning, allowing the generation of diverse and functional protein sequences in unexplored regions of sequence space. Altogether, these methods enable the rapid exploration of sequence space within regions enriched with functional proteins and therefore have great potential for accelerating the engineering of stable, functional, and diverse proteins for industrial and biomedical applications.

Nucleoside-Driven Specificity of DNA Methyltransferase.

We have studied the adenosine binding specificities of two bacterial DNA methyltransferases, Taq methyltransferase (M.TaqI), and HhaI methyltransferase (M.HhaI). While they have similar cofactor binding pocket interactions, experimental data showed different specificity for novel S-nucleobase-l-methionine cofactors (SNMs; N=guanosyl, cytidyl, uridyl). Protein dynamics corroborate the experimental data on the cofactor specificities. For M.TaqI the specificity for S-adenosyl-l-methionine (SAM) is governed by the tight binding on the nucleoside part of the cofactor, while for M.HhaI the degree of freedom of the nucleoside chain allows the acceptance of other bases. The experimental data prove catalytically productive methylation by the M.HhaI binding pocket for all the SNMs. Our results suggest a new route for successful design of unnatural SNM analogues for methyltransferases as a tool for cofactor engineering.

Inter-paralog amino acid inversion events in large phylogenies of duplicated proteins.

Connecting protein sequence to function is becoming increasingly relevant since high-throughput sequencing studies accumulate large amounts of genomic data. In order to go beyond the existing database annotation, it is fundamental to understand the mechanisms underlying functional inheritance and divergence. If the homology relationship between proteins is known, can we determine whether the function diverged? In this work, we analyze different possibilities of protein sequence evolution after gene duplication and identify "inter-paralog inversions", i.e., sites where the relationship between the ancestry and the functional signal is decoupled. The amino acids in these sites are masked from being recognized by other prediction tools. Still, they play a role in functional divergence and could indicate a shift in protein function. We develop a method to specifically recognize inter-paralog amino acid inversions in a phylogeny and test it on real and simulated datasets. In a dataset built from the Epidermal Growth Factor Receptor (EGFR) sequences found in 88 fish species, we identify 19 amino acid sites that went through inversion after gene duplication, mostly located at the ligand-binding extracellular domain. Our work uncovers an outcome of protein duplications with direct implications in protein functional annotation and sequence evolution. The developed method is optimized to work with large protein datasets and can be readily included in a targeted protein analysis pipeline.

Evolutionary repair reveals an unexpected role of the tRNA modification m1G37 in aminoacylation.

The tRNA modification m1G37, introduced by the tRNA methyltransferase TrmD, is thought to be essential for growth in bacteria because it suppresses translational frameshift errors at proline codons. However, because bacteria can tolerate high levels of mistranslation, it is unclear why loss of m1G37 is not tolerated. Here, we addressed this question through experimental evolution of trmD mutant strains of Escherichia coli. Surprisingly, trmD mutant strains were viable even if the m1G37 modification was completely abolished, and showed rapid recovery of growth rate, mainly via duplication or mutation of the proline-tRNA ligase gene proS. Growth assays and in vitro aminoacylation assays showed that G37-unmodified tRNAPro is aminoacylated less efficiently than m1G37-modified tRNAPro, and that growth of trmD mutant strains can be largely restored by single mutations in proS that restore aminoacylation of G37-unmodified tRNAPro. These results show that inefficient aminoacylation of tRNAPro is the main reason for growth defects observed in trmD mutant strains and that proS may act as a gatekeeper of translational accuracy, preventing the use of error-prone unmodified tRNAPro in translation. Our work shows the utility of experimental evolution for uncovering the hidden functions of essential genes and has implications for the development of antibiotics targeting TrmD.

Binding of single-mutant epidermal growth factor (EGF) ligands alters the stability of the EGF receptor dimer and promotes growth signaling.

The epidermal growth factor receptor (EGFR) is a membrane-anchored tyrosine kinase that is able to selectively respond to multiple extracellular stimuli. Previous studies have indicated that the modularity of this system may be caused by ligand-induced differences in the stability of the receptor dimer. However, this hypothesis has not been explored using single-mutant ligands thus far. Herein, we developed a new approach to identify residues responsible for functional divergence by selecting residues in the epidermal growth factor (EGF) ligand that are conserved among orthologs yet divergent between paralogs. Then, we mutated these residues and assessed the mutants' effects on the receptor using a combination of molecular dynamics (MD) and biochemical techniques. Although the EGF mutants had binding affinities for the EGFR comparable with the WT ligand, the EGF mutants showed differential patterns of receptor phosphorylation and cell growth in multiple cell lines. The MD simulations of the EGF mutants indicated that mutations had long-range effects on the receptor dimer interface. This study shows for the first time that a single mutation in the EGF is sufficient to alter the activation of the EGFR signaling pathway at the cellular level. These results also support that biased ligand-receptor signaling in the tyrosine kinase receptor system can lead to differential downstream outcomes and demonstrate a promising new method to study ligand-receptor interactions.

Cycloserine enantiomers are reversible inhibitors of human alanine:glyoxylate aminotransferase: implications for Primary Hyperoxaluria type 1.

Peroxisomal alanine:glyoxylate aminotransferase (AGT) is responsible for glyoxylate detoxification in human liver and utilizes pyridoxal 5'-phosphate (PLP) as coenzyme. The deficit of AGT leads to Primary Hyperoxaluria Type I (PH1), a rare disease characterized by calcium oxalate stones deposition in the urinary tract as a consequence of glyoxylate accumulation. Most missense mutations cause AGT misfolding, as in the case of the G41R, which induces aggregation and proteolytic degradation. We have investigated the interaction of wild-type AGT and the pathogenic G41R variant with d-cycloserine (DCS, commercialized as Seromycin), a natural product used as a second-line treatment of multidrug-resistant tuberculosis, and its synthetic enantiomer l-cycloserine (LCS). In contrast with evidences previously reported on other PLP-enzymes, both ligands are AGT reversible inhibitors showing inhibition constants in the micromolar range. While LCS undergoes half-transamination generating a ketimine intermediate and behaves as a classical competitive inhibitor, DCS displays a time-dependent binding mainly generating an oxime intermediate. Using a mammalian cellular model, we found that DCS, but not LCS, is able to promote the correct folding of the G41R variant, as revealed by its increased specific activity and expression as a soluble protein. This effect also translates into an increased glyoxylate detoxification ability of cells expressing the variant upon treatment with DCS. Overall, our findings establish that DCS could play a role as pharmacological chaperone, thus suggesting a new line of intervention against PH1 based on a drug repositioning approach. To a widest extent, this strategy could be applied to other disease-causing mutations leading to AGT misfolding.

Rossmann-Fold Methyltransferases: Taking a "ß-Turn" around Their Cofactor, S-Adenosylmethionine.

Methyltransferases (MTases) are superfamilies of enzymes that catalyze the transfer of a methyl group from S-adenosylmethionine (SAM), a nucleoside-based cofactor, to a wide variety of substrates such as DNA, RNA, proteins, small molecules, and lipids. Depending upon their structural features, the MTases can be further classified into different classes; we consider exclusively the largest class of MTases, the Rossmann-fold MTases. It has been shown that the nucleoside cofactor-binding Rossmann enzymes, particularly the nicotinamide adenine dinucleotide (NAD)-, flavin adenine dinucleotide (FAD)-, and SAM-binding MTases enzymes, share common binding motifs that include a Gly-rich loop region that interacts with the cofactor and a highly conserved acidic residue (Asp/Glu) that interacts with the ribose moiety of the cofactor. Here, we observe that the Gly-rich loop region of the Rossmann MTases adapts a specific type II' ß-turn in the proximity of the cofactor (<4 Å), and it appears to be a key feature of these superfamilies. Additionally, we demonstrate that the conservation of this ß-turn could play a critical role in the enzyme-cofactor interaction, thereby shedding new light on the structural conformation of the Gly-rich loop region from Rossmann MTases.

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors.

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second ß-strand (ß2-Asp/Glu). We identified a canonical motif, defined by the ß2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose's ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint--geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the ß2-Asp/Glu motif.

Spontaneous Emergence of S-Adenosylmethionine and the Evolution of Methylation.

S-Adenosylmethionine (SAM) is an essential methylation cofactor. The origins of SAM methylation are complex, seemingly demanding the simultaneous emergence of an enzyme that makes SAM and enzyme(s) that utilize it. We report that both ATP and adenosine spontaneously react with methionine to yield SAM, thus suggesting that SAM could have emerged by chance. SAM methylation thus exemplifies how metabolites and pathways can co-emerge through the gradual recruitment of individual enzymes in reverse order.

Engineering and Directed Evolution of DNA Methyltransferases.

DNA methyltransferases (MTases) constitute an attractive target for protein engineering, thus opening the road to new ways of manipulating DNA in a unique and selective manner. Here, we review various aspects of MTase engineering, both methodological and conceptual, and also discuss future directions and challenges. Bacterial MTases that are part of restriction/modification (R/M) systems offer a convenient way for the selection of large gene libraries, both in vivo and in vitro. We review these selection methods, their strengths and weaknesses, and also the prospects for new selection approaches that will enable the directed evolution of mammalian DNA methyltransferases (Dnmts). We explore various properties of MTases that may be subject to engineering. These include engineering for higher stability and soluble expression (MTases, including bacterial ones, are prone to misfolding), engineering of the DNA target specificity, and engineering for the usage of S-adenosyl-L-methionine (AdoMet) analogs. Directed evolution of bacterial MTases also offers insights into how these enzymes readily evolve in nature, thus yielding MTases with a huge spectrum of DNA target specificities. Engineering for alternative cofactors, on the other hand, enables modification of DNA with various groups other than methyl and thus can be employed to map and redirect DNA epigenetic modifications.

DddD is a CoA-transferase/lyase producing dimethyl sulfide in the marine environment.

Dimethyl sulfide (DMS) is produced in oceans in vast amounts (>10(7) tons/year) and mediates a wide range of processes from regulating marine life forms to cloud formation. Nonetheless, none of the enzymes that produce DMS from dimethylsulfoniopropionate (DMSP) has been adequately characterized. We describe the expression and purification of DddD from the marine bacterium Marinomonas sp. MWYL1 and its biochemical characterization. We identified DMSP and acetyl-coenzyme A to be DddD's native substrates and Asp602 as the active site residue mediating the CoA-transferase prior to lyase activity. These findings shed light on the biochemical utilization of DMSP in the marine environment.

Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site.

Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.

Non-consecutive enzyme interactions within TCA cycle supramolecular assembly regulate carbon-nitrogen metabolism.

Enzymes of the central metabolism tend to assemble into transient supramolecular complexes. However, the functional significance of the interactions, particularly between enzymes catalyzing non-consecutive reactions, remains unclear. Here, by co-localizing two non-consecutive enzymes of the TCA cycle from Bacillus subtilis, malate dehydrogenase (MDH) and isocitrate dehydrogenase (ICD), in phase separated droplets we show that MDH-ICD interaction leads to enzyme agglomeration with a concomitant enhancement of ICD catalytic rate and an apparent sequestration of its reaction product, 2-oxoglutarate. Theory demonstrates that MDH-mediated clustering of ICD molecules explains the observed phenomena. In vivo analyses reveal that MDH overexpression leads to accumulation of 2-oxoglutarate and reduction of fluxes flowing through both the catabolic and anabolic branches of the carbon-nitrogen intersection occupied by 2-oxoglutarate, resulting in impeded ammonium assimilation and reduced biomass production. Our findings suggest that the MDH-ICD interaction is an important coordinator of carbon-nitrogen metabolism.

Toxoplasma gondii secretory proteins bind to sulfated heparin structures.

Toxoplasma gondii is the causative agent of toxoplasmosis, one of the most widespread infections in humans and animals, and is a major opportunistic pathogen in immunocompromised patients. Toxoplasma gondii is unique as it can invade virtually any nucleated cell, although the mechanisms are not completely understood. Parasite attachment to the host cell is a prerequisite for reorientation and penetration and likely requires the recognition of molecules at the host cell surface. It has been reported that the affinity of tachyzoites, the invasive form of T. gondii, for host cells can be inhibited by a variety of soluble-sulfated glycosaminoglycans (GAGs), such as heparan sulfate. Using heparin-functionalized zeolites in the absence of host cells, we visualized heparin-binding sites on the surface of tachyzoites by confocal and atomic force microscopy. Furthermore, we report that protein components of the parasite rhoptry, dense granule and surface bind GAGs. In particular, the proteins ROP2 and ROP4 from the rhoptry, GRA2 from the dense granules and the surface protein SAG1 were found to bind heparin. The binding specificities and affinities of individual parasite proteins for natural heparin and heparin oligosaccharides were determined by a combination of heparin oligosaccharide microarrays and surface plasmon resonance. Our results suggest that interactions between sulfated GAGs and parasite surface antigens contribute to T. gondii attachment to host cell surfaces as well as initiating the invasion process, while rhoptries and dense granule organelles may play an important role during the establishment of the infection and during the life of the parasite inside the parasitophorous vacuole.