Mukund Thattai

The Origins of Specificity in Polyketide Synthase Protein Interactions

Polyketides, a diverse group of heteropolymers with antibiotic and antitumor properties, are assembled in bacteria by multiprotein chains of modular polyketide synthase (PKS) proteins. Specific protein–protein interactions determine the order of proteins within a multiprotein chain, and thereby the order in which chemically distinct monomers are added to the growing polyketide product. Here we investigate the evolutionary and molecular origins of protein interaction specificity. We focus on the short, conserved N- and C-terminal docking domains that mediate interactions between modular PKS proteins. Our computational analysis, which combines protein sequence data with experimental protein interaction data, reveals a hierarchical interaction specificity code. PKS docking domains are descended from a single ancestral interacting pair, but have split into three phylogenetic classes that are mutually noninteracting. Specificity within one such compatibility class is determined by a few key residues, which can be used to define compatibility subclasses. We identify these residues using a novel, highly sensitive co-evolution detection algorithm called CRoSS (correlated residues of statistical significance). The residue pairs selected by CRoSS are involved in direct physical interactions in a docked-domain NMR structure. A single PKS system can use docking domain pairs from multiple classes, as well as domain pairs from multiple subclasses of any given class. The termini of individual proteins are frequently shuffled, but docking domain pairs straddling two interacting proteins are linked as an evolutionary module. The hierarchical and modular organization of the specificity code is intimately related to the processes by which bacteria generate new PKS pathways.

Odor representations in the rat olfactory bulb change smoothly with morphing stimuli.

Summary Many species of mammals are very good at categorizing odors. One model for how this is achieved involves the formation of “attractor” states in the olfactory processing pathway, which converge to stable representations for the odor. We analyzed the responses of rat olfactory bulb mitral/tufted (M/T) cells using stimuli “morphing” from one odor to another through intermediate mixtures. We then developed a phenomenological model for the representation of odors and mixtures by M/T cells and show that >80% of odorant responses to different concentrations and mixtures can be expressed in terms of smoothly summing responses to air and the two pure odorants. Furthermore, the model successfully predicts M/T cell responses to odor mixtures when respiration dependence is eliminated. Thus, odor mixtures are represented in the bulb through summation of components, rather than distinct attractor states. We suggest that our olfactory coding model captures many aspects of single and mixed odor representation in M/T cells.

Ancient dynamin segments capture early stages of host-mitochondrial integration.

Significance Mitochondria were originally free-living bacteria with their own division machinery, which took up residence within another cell 2 billion y ago. The host cell subsequently tamed mitochondrial division using dynamin, a membrane-pinching protein. We have found that a single ancient dynamin at the root of extant eukaryotes, capable of pinching both mitochondria and vesicles, duplicated independently in plants and animals into specialized mitochondrial and vesicle variants. A “living fossil” of this ancient bifunctional dynamin still survives in scattered eukaryotic species, along with the original bacterial FtsZ division protein. The mitochondria of these organisms, preserved as if in amber, might teach us how the fateful partnership between host and endosymbiont was first established. Eukaryotic cells use dynamins—mechano-chemical GTPases—to drive the division of endosymbiotic organelles. Here we probe early steps of mitochondrial and chloroplast endosymbiosis by tracing the evolution of dynamins. We develop a parsimony-based phylogenetic method for protein sequence reconstruction, with deep time resolution. Using this, we demonstrate that dynamins diversify through the punctuated transformation of sequence segments on the scale of secondary-structural elements. We find examples of segments that have remained essentially unchanged from the 1.8-billion-y-old last eukaryotic common ancestor to the present day. Stitching these together, we reconstruct three ancestral dynamins: The first is nearly identical to the ubiquitous mitochondrial division dynamins of extant eukaryotes, the second is partially preserved in the myxovirus-resistance-like dynamins of metazoans, and the third gives rise to the cytokinetic dynamins of amoebozoans and plants and to chloroplast division dynamins. The reconstructed sequences, combined with evolutionary models and published functional data, suggest that the ancestral mitochondrial division dynamin also mediated vesicle scission. This bifunctional protein duplicated into specialized mitochondrial and vesicle variants at least three independent times—in alveolates, green algae, and the ancestor of fungi and metazoans—accompanied by the loss of the ancient prokaryotic mitochondrial division protein FtsZ. Remarkably, many extant species that retain FtsZ also retain the predicted ancestral bifunctional dynamin. The mitochondrial division apparatus of such organisms, including amoebozoans, red algae, and stramenopiles, seems preserved in a near-primordial form.

Emergent gene order in a model of modular polyketide synthases

Polyketides are a class of biologically active heteropolymers produced by assembly line-like multiprotein complexes of modular polyketide synthases (PKS). The polyketide product is encoded in the order of the PKS proteins in the assembly line, suggesting that polyketide diversity derives from combinatorial rearrangement of these PKS complexes. Remarkably, the order of PKS genes on the chromosome follows the order of PKS proteins in the assembly line: This fact is commonly referred to as “collinearity”. Here we propose an evolutionary origin for collinearity and demonstrate the mechanism by using a computational model of PKS evolution in a population. Assuming continuous evolutionary pressure for novel polyketides, and that new polyketide pathways are formed by horizontal transfer/recombination of PKS-encoding DNA, we demonstrate the existence of a broad range of parameters for which collinearity emerges spontaneously. Collinearity confers no fitness advantage in our model; it is established and maintained through a “secondary selection” mechanism, as a trait which increases the probability of forming long, novel PKS complexes through recombination. Consequently, collinearity hitchhikes on the successful genotypes which periodically sweep through the evolving population. In addition to computer simulation of a simplified model of PKS evolution, we provide a mathematical framework describing the secondary selection mechanism, which generalizes beyond the context of the present model.

Prediction by Promoter Logic in Bacterial Quorum Sensing

Quorum-sensing systems mediate chemical communication between bacterial cells, coordinating cell-density-dependent processes like biofilm formation and virulence-factor expression. In the proteobacterial LuxI/LuxR quorum sensing paradigm, a signaling molecule generated by an enzyme (LuxI) diffuses between cells and allosterically stimulates a transcriptional regulator (LuxR) to activate its cognate promoter (pR). By expressing either LuxI or LuxR in positive feedback from pR, these versatile systems can generate smooth (monostable) or abrupt (bistable) density-dependent responses to suit the ecological context. Here we combine theory and experiment to demonstrate that the promoter logic of pR – its measured activity as a function of LuxI and LuxR levels – contains all the biochemical information required to quantitatively predict the responses of such feedback loops. The interplay of promoter logic with feedback topology underlies the versatility of the LuxI/LuxR paradigm: LuxR and LuxI positive-feedback systems show dramatically different responses, while a dual positive/negative-feedback system displays synchronized oscillations. These results highlight the dual utility of promoter logic: to probe microscopic parameters and predict macroscopic phenotype.

On the Archaeal Origins of Eukaryotes and the Challenges of Inferring Phenotype from Genotype.

If eukaryotes arose through a merger between archaea and bacteria, what did the first true eukaryotic cell look like? A major step toward an answer came with the discovery of Lokiarchaeum, an archaeon whose genome encodes small GTPases related to those used by eukaryotes to regulate membrane traffic. Although ‘Loki’ cells have yet to be seen, their existence has prompted the suggestion that the archaeal ancestor of eukaryotes engulfed the future mitochondrion by phagocytosis. We propose instead that the archaeal ancestor was a relatively simple cell, and that eukaryotic cellular organization arose as the result of a gradual transfer of bacterial genes and membranes driven by an ever-closer symbiotic partnership between a bacterium and an archaeon.

New Organelles by Gene Duplication in a Biophysical Model of Eukaryote Endomembrane Evolution

Extant eukaryotic cells have a dynamic traffic network that consists of diverse membrane-bound organelles exchanging matter via vesicles. This endomembrane system arose and diversified during a period characterized by massive expansions of gene families involved in trafficking after the acquisition of a mitochondrial endosymbiont by a prokaryotic host cell >1.8 billion years ago. Here we investigate the mechanistic link between gene duplication and the emergence of new nonendosymbiotic organelles, using a minimal biophysical model of traffic. Our model incorporates membrane-bound compartments, coat proteins and adaptors that drive vesicles to bud and segregate cargo from source compartments, and SNARE proteins and associated factors that cause vesicles to fuse into specific destination compartments. In simulations, arbitrary numbers of compartments with heterogeneous initial compositions segregate into a few compositionally distinct subsets that we term organelles. The global structure of the traffic system (i.e., the number, composition, and connectivity of organelles) is determined completely by local molecular interactions. On evolutionary timescales, duplication of the budding and fusion machinery followed by loss of cross-interactions leads to the emergence of new organelles, with increased molecular specificity being necessary to maintain larger organellar repertoires. These results clarify potential modes of early eukaryotic evolution as well as more recent eukaryotic diversification.

Evolutionary cell biology: Lessons from diversity

Novel perspectives emerge from a recent conference on the origins of eukaryotic cells, which covered phylogenetics, population genetics and evolutionary consequences of energy requirements and host–pathogen interactions.

Interactions among Quorum Sensing Inhibitors

Many pathogenic bacteria use quorum sensing (QS) systems to regulate the expression of virulence genes in a density-dependent manner. In one widespread QS paradigm the enzyme LuxI generates a small diffusible molecule of the acyl-homoserine lactone (AHL) family; high cell densities lead to high AHL levels; AHL binds the transcription factor LuxR, triggering it to activate gene expression at a virulence promoter. The emergence of antibiotic resistance has generated interest in alternative anti-microbial therapies that target QS. Inhibitors of LuxI and LuxR have been developed and tested in vivo, and can act at various levels: inhibiting the synthesis of AHL by LuxI, competitively or non-competitively inhibiting LuxR, or increasing the turnover of LuxI, LuxR, or AHL. Here use an experimentally validated computational model of LuxI/LuxR QS to study the effects of using inhibitors individually and in combination. The model includes the effect of transcriptional feedback, which generates highly non-linear responses as inhibitor levels are increased. For the ubiquitous LuxI-feedback virulence systems, inhibitors of LuxI are more effective than those of LuxR when used individually. Paradoxically, we find that LuxR competitive inhibitors, either individually or in combination with other inhibitors, can sometimes increase virulence by weakly activating LuxR. For both LuxI-feedback and LuxR-feedback systems, a combination of LuxR non-competitive inhibitors and LuxI inhibitors act multiplicatively over a broad parameter range. In our analysis, this final strategy emerges as the only robust therapeutic option.

Universal Poisson Statistics of mRNAs with Complex Decay Pathways

Messenger RNA (mRNA) dynamics in single cells are often modeled as a memoryless birth-death process with a constant probability per unit time that an mRNA molecule is synthesized or degraded. This predicts a Poisson steady-state distribution of mRNA number, in close agreement with experiments. This is surprising, since mRNA decay is known to be a complex process. The paradox is resolved by realizing that the Poisson steady state generalizes to arbitrary mRNA lifetime distributions. A mapping between mRNA dynamics and queueing theory highlights an identifiability problem: a measured Poisson steady state is consistent with a large variety of microscopic models. Here, I provide a rigorous and intuitive explanation for the universality of the Poisson steady state. I show that the mRNA birth-death process and its complex decay variants all take the form of the familiar Poisson law of rare events, under a nonlinear rescaling of time. As a corollary, not only steady-states but also transients are Poisson distributed. Deviations from the Poisson form occur only under two conditions, promoter fluctuations leading to transcriptional bursts or nonindependent degradation of mRNA molecules. These results place severe limits on the power of single-cell experiments to probe microscopic mechanisms, and they highlight the need for single-molecule measurements.