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.

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.

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.

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.

Mechanisms of protein evolution

How do proteins evolve? How do changes in sequence mediate changes in protein structure, and in turn in function? This question has multiple angles, ranging from biochemistry and biophysics to evolutionary biology. This review provides a brief integrated view of some key mechanistic aspects of protein evolution. First, we explain how protein evolution is primarily driven by randomly acquired genetic mutations and selection for function, and how these mutations can even give rise to completely new folds. Then, we also comment on how phenotypic protein variability, including promiscuity, transcriptional and translational errors, may also accelerate this process, possibly via "plasticity-first" mechanisms. Finally, we highlight open questions in the field of protein evolution, with respect to the emergence of more sophisticated protein systems such as protein complexes, pathways, and the emergence of pre-LUCA enzymes.

Structural dynamics shape the fitness window of alanine:glyoxylate aminotransferase

The conformational landscape of a protein is constantly expanded by genetic variations that have a minimal impact on the function(s) while causing subtle effects on protein structure. The wider the conformational space sampled by these variants, the higher the probabilities to adapt to changes in environmental conditions. However, the probability that a single mutation may result in a pathogenic phenotype also increases. Here we present a paradigmatic example of how protein evolution balances structural stability and dynamics to maximize protein adaptability and preserve protein fitness. We took advantage of known genetic variations of human alanine:glyoxylate aminotransferase (AGT1), which is present as a common major allelic form (AGT‐Ma) and a minor polymorphic form (AGT‐Mi) expressed in 20% of Caucasian population. By integrating crystallographic studies and molecular dynamics simulations, we show that AGT‐Ma is endowed with structurally unstable (frustrated) regions, which become disordered in AGT‐Mi. An in‐depth biochemical characterization of variants from an anticonsensus library, encompassing the frustrated regions, correlates this plasticity to a fitness window defined by AGT‐Ma and AGT‐Mi. Finally, co‐immunoprecipitation analysis suggests that structural frustration in AGT1 could favor additional functions related to protein–protein interactions. These results expand our understanding of protein structural evolution by establishing that naturally occurring genetic variations tip the balance between stability and frustration to maximize the ensemble of conformations falling within a well‐defined fitness window, thus expanding the adaptability potential of the protein.

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.

Proteins are critical components of living systems because they facilitate most biological processes like protein synthesis, DNA replication, chemical catalysis, etc. Proteins are encoded in their genes. During evolution, genes accumulate mutations that get translated at the protein level. These mutations can be “neutral” if they do not affect the protein function immediately and directly; otherwise, mutations can be functional if they directly modify protein function. An event that provides an opportunity to study protein function is gene duplication namely, when two copies of a gene encoding the same protein appear. One copy of the protein often retains the same function while the other is free to diverge and specialize to a different function. This work sheds light on an alternative outcome of gene duplication that might be critical to discern between neutral and functional mutations. By looking at 88 fish genomes, we found proteins in which the evolution of their sequences does not follow the expected pattern of divergence after gene duplication. In this case, the protein sequence of a subgroup of species diverges in the copy expected to retain its function, while the sequence is retained in the expectedly divergent one. We called this event “inter-paralog amino acid inversion”. Our data shows that this “inversion” event is correlated to function, and its detection has to be considered for assigning protein functions correctly.

Correction to "Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases"

[This corrects the article DOI: 10.1021/jacsau.1c00464.].

Implications of divergence of methionine adenosyltransferase in archaea

Methionine adenosyltransferase (MAT) catalyzes the biosynthesis of S-adenosyl methionine from l-methionine and ATP. MAT enzymes are ancient, believed to share a common ancestor, and are highly conserved in all three domains of life. However, the sequences of archaeal MATs show considerable divergence compared with their bacterial and eukaryotic counterparts. Furthermore, the structural significance and functional significance of this sequence divergence are not well understood. In the present study, we employed structural analysis and ancestral sequence reconstruction to investigate archaeal MAT divergence. We observed that the dimer interface containing the active site (which is usually well conserved) diverged considerably between the bacterial/eukaryotic MATs and archaeal MAT. A detailed investigation of the available structures supports the sequence analysis outcome: The protein domains and subdomains of bacterial and eukaryotic MAT are more similar than those of archaea. Finally, we resurrected archaeal MAT ancestors. Interestingly, archaeal MAT ancestors show substrate specificity, which is lost during evolution. This observation supports the hypothesis of a common MAT ancestor for the three domains of life. In conclusion, we have demonstrated that archaeal MAT is an ideal system for studying an enzyme family that evolved differently in one domain compared with others while maintaining the same catalytic activity.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

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

The tRNA modification m¹G37, 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 m¹G37 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 m¹G37 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 tRNA^Pro is aminoacylated less efficiently than m¹G37-modified tRNA^Pro, and that growth of trmD mutant strains can be largely restored by single mutations in proS that restore aminoacylation of G37-unmodified tRNA^Pro. These results show that inefficient aminoacylation of tRNA^Pro 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 tRNA^Pro 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.

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.

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.

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.

•

A new computational method is used to stabilize five recalcitrant proteins

•

Designed variants show higher expression and stability with unmodified function

•

A designed human acetylcholinesterase variant expresses solubly in bacteria

•

The method is fully automated and implemented on a webserver

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.

6-methoxy-7-benzofuranoxy and 6-methoxy-7-indolyloxy analogues of 2-[2-(2,6-Dimethoxyphenoxy)ethyl]aminomethyl-1,4-benzodioxane (WB4101):1 discovery of a potent and selective α1D-adrenoceptor antagonist

Previous results have shown that replacement of one of the two o-methoxy groups at the phenoxy residue of the potent, but not subtype-selective, α1-AR antagonist (S)-WB4101 [(S)-1] by phenyl, or by ortho,meta-fused cyclohexane, or especially by ortho,meta-fused benzene preferentially elicits α1D-AR antagonist affinity. Such observations inspired the design of four new analogues of 1 bearing, in lieu of the 2,6-dimethoxyphenoxy residue, a 6-methoxy-substituted 7-benzofuranoxy or 7-indolyloxy group or, alternatively, their corresponding 2,3-dihydro form. Of these new compounds, which maintain, rigidified, the characteristic ortho heterodisubstituted phenoxy substructure of 1, the S enantiomer of the dihydrobenzofuranoxy derivative exhibited the highest α1D-AR antagonist affinity (pA2 9.58) with significant α1D/α1A and α1D/α1B selectivity. In addition, compared both to α1D-AR antagonists structurally related to 1 and to the well-known α1D-AR antagonist BMY7378, this derivative had modest 5-HT1A affinity and neutral α1-AR antagonist behavior.

Continuous-flow reactor-based synthesis of carbohydrate and dihydrolipoic acid-capped quantum dots

A detailed protocol for the large-scale synthesis of carbohydrate and dihydrolipoic acid (DHLA)-coated CdSe/ZnS and CdTe/ZnS nanoparticles using continuous flow reactors is described here. Three continuous flow microreaction systems, operating at three different temperatures, are used for the synthesis of mannose-, galactose- or DHLA-functionalized quantum dots (QDs). In the first step of synthesis, the CdSe and CdTe nanoparticles are prepared. The size and spectral properties of the CdSe core of the nanoparticles are controlled by adjustment of the residence time and the temperature. As a second step, the zinc sulfide capping under homogenous conditions is carried out at a substantially lower temperature than is required for nanoparticle growth in batch processes. Finally, the trioctylphosphine/oleic acid ligand is effectively replaced with either carbohydrate PEG-thiol moieties or DHLA at 60 °C. This new protocol allows the synthesis of biologically active fluorescent QDs in 4 d.

Detection of bacteria using glyco-dendronized polylysine prepared by continuous flow photofunctionalization

Biocompatible glyco-dendronized poly-l-lysine (PLL) polymers carry either three or nine mannose- or galactose-bearing dendrons that selectively bind, and thus can be used to detect, bacteria. Central to the synthesis of glyco-dendronized polymers was the development of a continuous flow [2 + 2] photocycloaddition reaction to connect the dendrons and PLL. Glycodendronized polymers cluster bacteria by binding to cell-surface carbohydrate receptors and thereby result in an easy read-out using microscopic analyses.

Design, synthesis and biological evaluation of carbohydrate-functionalized cyclodextrins and liposomes for hepatocyte-specific targeting

Targeting glycan-binding receptors is an attractive strategy for cell-specific drug and gene delivery. The C-type lectin asialoglycoprotein receptor (ASGPR) is particularly suitable for liver-specific delivery due to its exclusive expression by parenchymal hepatocytes. In this study, we designed and developed an efficient synthesis of carbohydrate-functionalized β-cyclodextrins (βCDs) and liposomes for hepatocyte-specific delivery. For targeting of ASGPR, rhodamine B-loaded βCDs were functionalized with glycodendrimers. Liposomes were equipped with synthetic glycolipids containing a terminal D-GalNAc residue to mediate binding to ASGPR. Uptake studies in the human hepatocellular carcinoma cell line HepG2 demonstrated that βCDs and liposomes displaying terminal D-Gal/D-GalNAc residues were preferentially endocytosed. In contrast, uptake of βCDs and liposomes with terminal d-Man or D-GlcNAc residues was markedly reduced. The d-Gal/d-GalNAc-functionalized βCDs and liposomes presented here enable hepatocyte-specific targeting. Gal-functionalized βCDs are efficient molecular carriers to deliver doxorubicin in vitro into hepatocytes and induce apoptosis.

Asymmetric reactions in continuous flow

An overview of asymmetric synthesis in continuous flow and microreactors is presented in this review. Applications of homogeneous and heterogeneous asymmetric catalysis as well as biocatalysis in flow are discussed.

In vitro imaging and in vivo liver targeting with carbohydrate capped quantum dots

PEGylated quantum dots (QDs) capped with d-mannose, d-galactose, and d-galactosamine have been synthesized. The stable, high quantum yield fluorescence of QDs was exploited to study specific carbohydrate-protein interactions in vitro and in vivo.