BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Tsvi Tlusty

Tsvi Tlusty contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Biological PhysicsBiomoleculescond-mat.softcond-mat.stat-mechMolecular NetworksInformation Theorymath.ITNeurons and Cognitioncond-mat.dis-nnphysics.chem-phphysics.flu-dynPopulations and EvolutionQuantitative MethodsGenomicsArtificial IntelligenceCell BehaviorComputation and Languagecond-mat.mtrl-scieess.ASmath-phmath.MPnlin.AOphysics.data-anphysics.soc-phSoundSubcellular Processes

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
44works
0followers
26topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Acceleration of enzymatic catalysis by active hydrodynamic fluctuations

The cellular milieu is teeming with biochemical nano-machines whose activity is a strong source of correlated non-thermal fluctuations termed active noise. Essential elements of this circuitry are enzymes, catalysts that speed up the rate of metabolic reactions by orders of magnitude, thereby making life possible. Here, we examine the possibility that active noise in the cell, or in vitro, affects enzymatic catalytic rate by accelerating or decelerating the crossing rate of energy barriers during the reaction. Considering hydrodynamic perturbations induced by biochemical activity as a source of active noise, we evaluate their impact on the enzymatic cycle using a combination of analytic and numerical methods. Our estimates show that the fast component of the active noise spectrum enhances the rate of enzymes, while reactions remain practically unaffected by the slow noise spectrum. Revisiting the physics of barrier crossing under the influence of active hydrodynamic fluctuations suggests that the biochemical activity of macromolecules such as enzymes is coupled to active noise. Thus, we propose that enzymatic catalysis is a collective, many-body process in which enzymes may affect each other's activity via long-range hydrodynamic interaction, with potential impact on biochemical networks in living and artificial systems alike.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spontaneous Chiral Symmetry Breaking in a Random Driven Chemical System

Living systems have evolved to efficiently consume available energy sources using an elaborate circuitry of chemical reactions which, puzzlingly, bear a strict restriction to asymmetric chiral configurations. While autocatalysis is known to promote such chiral symmetry breaking, whether a similar phenomenon may also be induced in a more general class of configurable chemical systems -- via energy exploitation -- is a sensible yet underappreciated possibility. This work examines this question within a model of randomly generated complex chemical networks. We show that chiral symmetry breaking may occur spontaneously and generically by harnessing energy sources from external environmental drives. Key to this transition are intrinsic fluctuations of achiral-to-chiral reactions and tight matching of system configurations to the environmental drives, which together amplify and sustain diverged enantiomer distributions. These asymmetric states emerge through steep energetic transitions from the corresponding symmetric states and sharply cluster as highly-dissipating states. The results thus demonstrate a generic mechanism in which energetic drives may give rise to homochirality in an otherwise totally symmetrical environment, and from an early-life perspective, might emerge as a competitive, energy-harvesting advantage.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inverse Design of Non-Equilibrium Steady-States: A Large Deviation Approach

The design of small scale non-equilibrium steady states (NESS) is a challenging, open ended question. While similar equilibrium problems are tractable using standard thermodynamics, a generalized description for non-equilibrium systems is lacking, making the design problem particularly difficult. Here we show we can exploit the large deviation behavior of a Brownian particle and design a variety of geometrically complex steady-state density distributions and flux field flows. We achieve this design target from direct knowledge of the joint large deviation functional for the empirical density and flow, and a "relaxation" algorithm on the desired target states via adjustable force field parameters. We validate the method by replicating analytical results, and demonstrate its capacity to yield complex prescribed targets, such as rose-curve or polygonal shapes on the plane. We consider this dynamical fluctuation approach a first step towards the design of more complex NESS where general frameworks are otherwise still lacking.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Positional information as a universal predictor of freezing

Variation of positional information, measured by the two-body excess entropy $\mathsf{S}_\mathrm{2}$, is studied across the liquid-solid equilibrium transition in a simple two-dimensional system. Analysis reveals a master relation between $\mathsf{S}_\mathrm{2}$ and the freezing temperature $T_{\mathrm{f}}$, from which a scaling law is extracted: $-\mathsf{S}_\mathrm{2}\sim |T_{\mathrm{f}} - T| ^{-1/3}$. Theoretical and practical implications of the observed universality are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

$αβ$ DCA method identifies unspecific binding but specific disruption of the group I intron by the StpA chaperone

Chaperone protein - the most disordered among all protein groups - help RNAs fold into their functional structure by destabilizing misfolded configurations or stabilizing the functional ones. But disentangling the mechanism underlying RNA chaperoning is challenging, mostly due to inherent disorder of the chaperones and the transient nature of their interactions with RNA. In particular, it is unclear how specific the interactions are and what role is played by amino acid charge and polarity patterns. Here, we address these questions in the RNA chaperone StpA. We adapted direct coupling analysis (DCA) into the $αβ$DCA method that can treat in tandem sequences written in two alphabets, nucleotides and amino acids. With $αβ$DCA, we could analyze StpA-RNA interactions and show consistency with a previously proposed two-pronged mechanism: StpA disrupts specific positions in the group I intron while globally and loosely binding to the entire structure. Moreover, the interactions are strongly associated with the charge pattern: negatively charged regions in the destabilizing StpA N-terminal affect a few specific positions in the RNA, located in stems and in the pseudoknot. In contrast, positive regions in the C-terminal contain strongly coupled amino acids that promote non-specific or weakly-specific binding to the RNA. The present study opens new avenues to examine the functions of disordered proteins and to design disruptive proteins based on their charge patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cross-cultural data shows musical scales evolved to maximise imperfect fifths

Musical scales are used throughout the world, but the question of how they evolved remains open. Some suggest that scales based on the harmonic series are inherently pleasant, while others propose that scales are chosen that are easy to communicate. However, testing these theories has been hindered by the sparseness of empirical evidence. Here, we assimilate data from diverse ethnomusicological sources into a cross-cultural database of scales. We generate populations of scales based on multiple theories and assess their similarity to empirical distributions from the database. Most scales tend to include intervals which are close in size to perfect fifths (``imperfect fifths''), and packing arguments explain the salient features of the distributions. Scales are also preferred if their intervals are compressible, which may facilitate efficient communication and memory of melodies. While scales appear to evolve according to various selection pressures, the simplest h imperfect-fifths packing model best fits the empirical data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinetic Proofreading and the Limits of Thermodynamic Uncertainty

To mitigate errors induced by the cell's heterogeneous noisy environment, its main information channels and production networks utilize the kinetic proofreading (KPR) mechanism. Here, we examine two extensively-studied KPR circuits, DNA replication by the T7 DNA polymerase and translation by the E. coli ribosome. Using experimental data, we analyze the performance of these two vital systems in light of the fundamental bounds set by the recently-discovered thermodynamic uncertainty relation (TUR), which places an inherent trade-off between the precision of a desirable output and the amount of energy dissipation required. We show that the DNA polymerase operates close to the TUR lower bound, while the ribosome operates $\sim5$ times farther from this bound. This difference originates from the enhanced binding discrimination of the polymerase which allows it to operate effectively as a reduced reaction cycle prioritizing correct product formation. We show that approaching this limit also decouples the thermodynamic uncertainty factor from speed and error, thereby relaxing the accuracy-speed trade-off of the system. Altogether, our results show that operating near this reduced cycle limit not only minimizes thermodynamic uncertainty, but also results in global performance enhancement of KPR circuits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Master curve of boosted diffusion for ten catalytic enzymes

Molecular agitation more rapid than thermal Brownian motion is reported for cellular environments, motor proteins, synthetic molecular motors, enzymes, and common chemical reactions, yet that chemical activity couples to molecular motion contrasts with generations of accumulated knowledge about diffusion at equilibrium. To test the limits of this idea, a critical testbed is mobility of catalytically active enzymes. Sentiment is divided about reality of enhanced enzyme diffusion with evidence for and against. Here a master curve shows that enzyme diffusion coefficient increases in proportion to the energy release rate, the product of Michaelis-Menten reaction rate and Gibbs free energy change with the highly satisfactory correlation coefficient of 0.97. For ten catalytic enzymes (urease, acetylcholinesterase, seven enzymes from the glucose cascade cycle, and another), our measurements span from roughly 40% enhanced diffusion coefficient at high turnover rate and negative Gibbs free energy to no enhancement at slow turnover rate and positive Gibbs free energy. Moreover, two independent measures of mobility show consistency, provided that one avoids undesirable fluorescence photophysics. The master curve presented here quantifies the limits of both ideas, that enzymes display enhanced diffusion and that they do not within instrumental resolution, and has possible implications for understanding enzyme mobility in cellular environments. The striking linear dependence for the exergonic enzymes (negative Gibbs free energy) together with the vanishing effect for endergonic enzyme (positive Gibbs free energy) are consistent with a physical picture where the mechanism boosting the diffusion is an active one, utilizing the available work from the chemical reaction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reaching and violating thermodynamic uncertainty bounds in information engines

Thermodynamic uncertainty relations (TURs) set fundamental bounds on the fluctuation and dissipation of stochastic systems. Here, we examine these bounds, in experiment and theory, by exploring the entire phase space of a cyclic information engine operating in a non-equilibrium steady state. Close to its maximal efficiency, we find that the engine violates the original TUR. This first experimental demonstration of TUR violation agrees with recently proposed softer bounds: The engine satisfies two generalized TUR bounds derived from the detailed fluctuation theorem with feedback control and another bound linking fluctuation and dissipation to mutual information and Renyi divergence. We examine how the interplay of work fluctuation and dissipation shapes the information conversion efficiency of the engine, and find that dissipation is minimal at a finite noise level, where the original TUR is violated.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Efficiency fluctuations and noise induced refrigerator-to-heater transition in information engines

Understanding noisy information engines is a fundamental problem of non-equilibrium physics, particularly in biomolecular systems agitated by thermal and active fluctuations in the cell. By the generalized second law of thermodynamics, the efficiency of these engines is bounded by the mutual information passing through their noisy feedback loop. Yet, direct measurement of the interplay between mutual information and energy has so far been elusive. To allow such examination, we explore here the entire phase-space of a noisy colloidal information engine, and study efficiency fluctuations due to the stochasticity of the mutual information and extracted work. We find that the average efficiency is maximum for non-zero noise level, at which the distribution of efficiency switches from bimodal to unimodal, and the stochastic efficiency often exceeds unity. We identify a line of anomalous, noise-driven equilibrium states that defines a refrigerator-to-heater transition, and test the generalized integral fluctuation theorem for continuous engines.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Enhanced diffusion and enzyme dissociation

The concept that catalytic enzymes can act as molecular machines transducing chemical activity into motion has conceptual and experimental support, but much of the claimed support comes from experimental conditions where the substrate concentration is higher than biologically relevant and accordingly exceeds kM, the Michaelis-Menten constant. Moreover, many of the enzymes studied experimentally to date are oligomeric. Urease, a hexamer of subunits, has been considered to be the gold standard demonstrating enhanced diffusion. Here we show that urease and certain other oligomeric enzymes of high catalytic activity above kM dissociate into their smaller subunit fragments that diffuse more rapidly, thus providing a simple physical mechanism of enhanced diffusion in this regime of concentrations. Mindful that this conclusion may be controversial, our findings are sup-ported by four independent analytical techniques, static light scattering, dynamic light scattering (DLS), size-exclusion chroma-tography (SEC), and fluorescence correlation spectroscopy (FCS). Data for urease are presented in the main text and the con-clusion is validated for hexokinase and acetylcholinesterase with data presented in supplementary information. For substrate concentration regimes below kM at which these enzymes do not dissociate, our findings from both FCS and DLS validate that enzymatic catalysis does lead to the enhanced diffusion phenomenon. INTRODUCT

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Large-Scale Survey of Cell-Differentiation Programs in a Generative Model Reveals Regeneration as an Epiphenomenon of Development

Development combines three basic processes asymmetric --- cell division, signaling and gene regulation --- in a multitude of ways to create an overwhelming diversity of multicellular life-forms. Here, we attempt to chart this diversity using a generative model. We sample millions of biologically feasible developmental schemes, allowing us to comment on the statistical properties of cell-differentiation trajectories they produce. Our results indicate that, in contrast to common views, cell-type lineage graphs are unlikely to be tree-like. Instead, they are more likely to be directed acyclic graphs, with multiple lineages converging on the same terminal cell-type. Additionally, in line with the hypothesis that whole body regeneration is an epiphenomenon of development, a majority of the `organism' generated by our model can regenerate using pluripotent cells. The generative framework is modular and flexible, and can be adapted to test additional hypotheses about general features of development.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Universality, scaling and collapse in supercritical fluids

The Supercritical Fluid (SCF) is known to exhibit salient dynamic and thermodynamic crossovers and inhomogeneous molecular distribution. But the question as to what basic physics underlies these microscopic and macroscopic anomalies remains open. Here, using an order parameter extracted by machine learning, the fraction of gas-like (or liquid-like) molecules, we find simplicity and universality in SCF: First, all isotherms of a given fluid collapse onto a single master curve described by a scaling relation. The observed power law holds from the high-temperature and pressure regime down to the critical point where it diverges. Second, phase diagrams of different compounds collapse onto their master curves by the same scaling exponent, thereby demonstrating a putative law of corresponding supercritical states in simple fluids. The reported results support a model of the SCF as a mixture of two interchangeable microstates, whose spatiotemporal dynamics gives rise to unique macroscopic properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Self reference in word definitions

Dictionaries are inherently circular in nature. A given word is linked to a set of alternative words (the definition) which in turn point to further descendants. Iterating through definitions in this way, one typically finds that definitions loop back upon themselves. The graph formed by such definitional relations is our object of study. By eliminating those links which are not in loops, we arrive at a core subgraph of highly connected nodes. We observe that definitional loops are conveniently classified by length, with longer loops usually emerging from semantic misinterpretation. By breaking the long loops in the graph of the dictionary, we arrive at a set of disconnected clusters. We find that the words in these clusters constitute semantic units, and moreover tend to have been introduced into the English language at similar times, suggesting a possible mechanism for language evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A Coarse-Grained Biophysical Model of E. coli and Its Application to Perturbation of the rRNA Operon Copy Number

We propose a biophysical model of Escherichia coli that predicts growth rate and an effective cellular composition from an effective, coarse-grained representation of its genome. We assume that E. coli is in a state of balanced exponential steadystate growth, growing in a temporally and spatially constant environment, rich in resources. We apply this model to a series of past measurements, where the growth rate and rRNA-to-protein ratio have been measured for seven E. coli strains with an rRNA operon copy number ranging from one to seven (the wild-type copy number). These experiments show that growth rate markedly decreases for strains with fewer than six copies. Using the model, we were able to reproduce these measurements. We show that the model that best fits these data suggests that the volume fraction of macromolecules inside E. coli is not fixed when the rRNA operon copy number is varied. Moreover, the model predicts that increasing the copy number beyond seven results in a cytoplasm densely packed with ribosomes and proteins. Assuming that under such overcrowded conditions prolonged diffusion times tend to weaken binding affinities, the model predicts that growth rate will not increase substantially beyond the wild-type growth rate, as indicated by other experiments. Our model therefore suggests that changing the rRNA operon copy number of wild-type E. coli cells growing in a constant rich environment does not substantially increase their growth rate. Other observations regarding strains with an altered rRNA operon copy number, such as nucleoid compaction and the rRNA operon feedback response, appear to be qualitatively consistent with this model. In addition, we discuss possible design principles suggested by the model and propose further experiments to test its validity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A colorful origin for the genetic code: Information theory, statistical mechanics and the emergence of molecular codes

The genetic code maps the sixty-four nucleotide triplets (codons) to twenty amino-acids. While the biochemical details of this code were unraveled long ago, its origin is still obscure. We review information-theoretic approaches to the problem of the code's origin and discuss the results of a recent work that treats the code in terms of an evolving, error-prone information channel. Our model - which utilizes the rate-distortion theory of noisy communication channels - suggests that the genetic code originated as a result of the interplay of the three conflicting evolutionary forces: the needs for diverse amino-acids, for error-tolerance and for minimal cost of resources. The description of the code as an information channel allows us to mathematically identify the fitness of the code and locate its emergence at a second-order phase transition when the mapping of codons to amino-acids becomes nonrandom. The noise in the channel brings about an error-graph, in which edges connect codons that are likely to be confused. The emergence of the code is governed by the topology of the error-graph, which determines the lowest modes of the graph-Laplacian and is related to the map coloring problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A model for the emergence of the genetic code as a transition in a noisy information channel

The genetic code maps the sixty-four nucleotide triplets (codons) to twenty amino-acids. Some argue that the specific form of the code with its twenty amino-acids might be a 'frozen accident' because of the overwhelming effects of any further change. Others see it as a consequence of primordial biochemical pathways and their evolution. Here we examine a scenario in which evolution drives the emergence of a genetic code by selecting for an amino-acid map that minimizes the impact of errors. We treat the stochastic mapping of codons to amino-acids as a noisy information channel with a natural fitness measure. Organisms compete by the fitness of their codes and, as a result, a genetic code emerges at a supercritical transition in the noisy channel, when the mapping of codons to amino-acids becomes nonrandom. At the phase transition, a small expansion is valid and the emergent code is governed by smooth modes of the Laplacian of errors. These modes are in turn governed by the topology of the error-graph, in which codons are connected if they are likely to be confused. This topology sets an upper bound - which is related to the classical map-coloring problem - on the number of possible amino-acids. The suggested scenario is generic and may describe a mechanism for the formation of other error-prone biological codes, such as the recognition of DNA sites by proteins in the transcription regulatory network.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A rate-distortion scenario for the emergence and evolution of noisy molecular codes

We discuss, in terms of rate-distortion theory, the fitness of molecular codes as the problem of designing an optimal information channel. The fitness is governed by an interplay between the cost and quality of the channel, which induces smoothness in the code. By incorporating this code fitness into population dynamics models, we suggest that the emergence and evolution of molecular codes may be explained by simple channel design considerations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A relation between the multiplicity of the second eigenvalue of a graph Laplacian, Courant's nodal line theorem and the substantial dimension of tight polyhedral surfaces

This note discusses a relation between the multiplicity m of the second eigenvalue λ2 of a Laplacian on a graph G, tight mappings of G and a discrete analogue of Courant's nodal line theorem. For a certain class of graphs, we show that the m-dimensional eigenspace of λ2 is tight and thus defines a tight mapping of G into an m-dimensional Euclidean space. The tightness of the mapping is shown to set Colin de Verdière's upper bound on the maximal λ2-multiplicity, where chr(γ(G)) is the chromatic number and γ(G) is the genus of G.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A simple model for the evolution of molecular codes driven by the interplay of accuracy, diversity and cost

Molecular codes translate information written in one type of molecules into another molecular language. We introduce a simple model that treats molecular codes as noisy information channels. An optimal code is a channel that conveys information accurately and efficiently while keeping down the impact of errors. The equipoise of the three conflicting needs, for minimal error-load, minimal cost of resources and maximal diversity of vocabulary, defines the fitness of the code. The model suggests a mechanism for the emergence of a code when evolution varies the parameters that control this equipoise and the mapping between the two molecular languages becomes non-random. This mechanism is demonstrated by a simple toy model that is formally equivalent to a mean-field Ising magnet.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anomalous Microfluidic Phonons Induced by the Interplay of Hydrodynamic Screening and Incompressibility

We investigate the acoustic normal modes ("phonons") of a 1D microfluidic droplet crystal at the crossover between 2D flow and confined 1D plug flow. The unusual phonon spectra of the crystal, which arise from long-range hydrodynamic interactions, change anomalously under confinement. The boundaries induce weakening and screening of the interactions, but when approaching the 1D limit we measure a marked increase in the crystal sound velocity, a sign of interaction strengthening. This non-monotonous behavior of the phonon spectra is explained theoretically by the interplay of screening and plug flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Casting Polymer Nets to Optimize Noisy Molecular Codes

Life relies on the efficient performance of molecular codes, which relate symbols and meanings via error-prone molecular recognition. We describe how optimizing a code to withstand the impact of molecular recognition noise may be approximated by the statistics of a two-dimensional network made of polymers. The noisy code is defined by partitioning the space of symbols into regions according to their meanings. The "polymers" are the boundaries between these regions and their statistics defines the cost and the quality of the noisy code. When the parameters that control the cost-quality balance are varied, the polymer network undergoes a first-order transition, where the number of encoded meanings rises discontinuously. Effects of population dynamics on the evolution of molecular codes are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coding limits on the number of transcription factors

Transcription factor proteins bind specific DNA sequences to control the expression of genes. They contain DNA binding domains which belong to several super-families, each with a specific mechanism of DNA binding. The total number of transcription factors encoded in a genome increases with the number of genes in the genome. Here, we examined the number of transcription factors from each super-family in diverse organisms. We find that the number of transcription factors from most super-families appears to be bounded. For example, the number of winged helix factors does not generally exceed 300, even in very large genomes. The magnitude of the maximal number of transcription factors from each super-family seems to correlate with the number of DNA bases effectively recognized by the binding mechanism of that super-family. Coding theory predicts that such upper bounds on the number of transcription factors should exist, in order to minimize cross-binding errors between transcription factors. This theory further predicts that factors with similar binding sequences should tend to have similar biological effect, so that errors based on mis-recognition are minimal. We present evidence that transcription factors with similar binding sequences tend to regulate genes with similar biological functions, supporting this prediction. The present study suggests limits on the transcription factor repertoire of cells, and suggests coding constraints that might apply more generally to the mapping between binding sites and biological function.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Conformational Proofreading: The Impact of Conformational Changes on the Specificity of Molecular Recognition

To perform recognition, molecules must locate and specifically bind their targets within a noisy biochemical environment with many look-alikes. Molecular recognition processes, especially the induced-fit mechanism, are known to involve conformational changes. This arises a basic question: does molecular recognition gain any advantage by such conformational changes? By introducing a simple statistical-mechanics approach, we study the effect of conformation and flexibility on the quality of recognition processes. Our model relates specificity to the conformation of the participant molecules and thus suggests a possible answer: Optimal specificity is achieved when the ligand is slightly off target, that is a conformational mismatch between the ligand and its main target improves the selectivity of the process. This indicates that deformations upon binding serve as a conformational proofreading mechanism, which may be selected for via evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cross-species analysis traces adaptation of Rubisco towards optimality in a low dimensional landscape

Rubisco, probably the most abundant protein in the biosphere, performs an essential part in the process of carbon fixation through photosynthesis thus facilitating life on earth. Despite the significant effect that Rubisco has on the fitness of plants and other photosynthetic organisms, this enzyme is known to have a remarkably low catalytic rate and a tendency to confuse its substrate, carbon dioxide, with oxygen. This apparent inefficiency is puzzling and raises questions regarding the roles of evolution versus biochemical constraints in shaping Rubisco. Here we examine these questions by analyzing the measured kinetic parameters of Rubisco from various organisms in various environments. The analysis presented here suggests that the evolution of Rubisco is confined to an effectively one-dimensional landscape, which is manifested in simple power law correlations between its kinetic parameters. Within this one dimensional landscape, which may represent biochemical and structural constraints, Rubisco appears to be tuned to the intracellular environment in which it resides such that the net photosynthesis rate is nearly optimal. Our analysis indicates that the specificity of Rubisco is not the main determinant of its efficiency but rather the tradeoff between the carboxylation velocity and CO2 affinity. As a result, the presence of oxygen has only moderate effect on the optimal performance of Rubisco, which is determined mostly by the local CO2 concentration. Rubisco appears as an experimentally testable example for the evolution of proteins subject both to strong selection pressure and to biochemical constraints which strongly confine the evolutionary plasticity to a low dimensional landscape.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Development of input connections in neural cultures

We introduce a novel approach for the quantitative assessment of the connectivity in neuronal cultures, based on the statistical mechanics of percolation on a graph. This allows us to follow the development of the culture and see the emergence of connectivity in the network. The culture becomes fully connected at a time equivalent to full term. The spontaneous bursting activity that characterizes cultures develops in parallel with the connectivity. The average number of inputs per neuron can be quantitatively determined in units of $m_0$, the number of activated inputs needed to excite the neuron. For $m_0\sim 10$ we find that hippocampal neurons have on average $\sim 40-80$ inputs while cortical neurons have $\sim 50-100$, depending on neuronal density. The ratio of excitatory to inhibitory neurons is determined using the GABA$_\text{A}$ antagonist bicuculine. This ratio changes during development and reaches the final value at day $7-8$, coinciding with the expected time of the GABA switch. For hippocampal cultures the inhibitory cells comprise about $30\%$ of the neurons in the culture while for cortical cultures they are about $20\%$. Such detailed global information on the connectivity of networks in neuronal cultures is at present inaccessible by any electrophysiological or other technique.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High fidelity of RecA-catalyzed recombination: a watchdog of genetic diversity

Homologous recombination plays a key role in generating genetic diversity, while maintaining protein functionality. The mechanisms by which RecA enables a single-stranded segment of DNA to recognize a homologous tract within a whole genome are poorly understood. The scale by which homology recognition takes place is of a few tens of base pairs, after which the quest for homology is over. To study the mechanism of homology recognition, RecA-promoted homologous recombination between short DNA oligomers with different degrees of heterology was studied in vitro, using fluorescence resonant energy transfer. RecA can detect single mismatches at the initial stages of recombination, and the efficiency of recombination is strongly dependent on the location and distribution of mismatches. Mismatches near the 5' end of the incoming strand have a minute effect, whereas mismatches near the 3' end hinder strand exchange dramatically. There is a characteristic DNA length above which the sensitivity to heterology decreases sharply. Experiments with competitor sequences with varying degrees of homology yield information about the process of homology search and synapse lifetime. The exquisite sensitivity to mismatches and the directionality in the exchange process support a mechanism for homology recognition that can be modeled as a kinetic proofreading cascade.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-Fidelity DNA Sensing by Protein Binding Fluctuations

One of the major functions of RecA protein in the cell is to bind single-stranded DNA exposed upon damage, thereby triggering the SOS repair response.We present fluorescence anisotropy measurements at the binding onset, showing enhanced DNA length discrimination induced by adenosine triphosphate consumption. Our model explains the observed DNA length sensing as an outcome of out-of equilibrium binding fluctuations, reminiscent of microtubule dynamic instability. The cascade architecture of the binding fluctuations is a generalization of the kinetic proofreading mechanism. Enhancement of precision by an irreversible multistage pathway is a possible design principle in the noisy biological environment.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Leaders of neuronal cultures in a quorum percolation model

We present a theoretical framework using quorum-percolation for describing the initiation of activity in a neural culture. The cultures are modeled as random graphs, whose nodes are excitatory neurons with kin inputs and kout outputs, and whose input degrees kin = k obey given distribution functions pk. We examine the firing activity of the population of neurons according to their input degree (k) classes and calculate for each class its firing probability Φ_k(t) as a function of t. The probability of a node to fire is found to be determined by its in-degree k, and the first-to-fire neurons are those that have a high k. A small minority of high-k classes may be called "Leaders", as they form an inter-connected subnetwork that consistently fires much before the rest of the culture. Once initiated, the activity spreads from the Leaders to the less connected majority of the culture. We then use the distribution of in-degree of the Leaders to study the growth rate of the number of neurons active in a burst, which was experimentally measured to be initially exponential. We find that this kind of growth rate is best described by a population that has an in-degree distribution that is a Gaussian centered around k = 75 with width σ = 31 for the majority of the neurons, but also has a power law tail with exponent -2 for ten percent of the population. Neurons in the tail may have as many as k = 4, 700 inputs. We explore and discuss the correspondence between the degree distribution and a dynamic neuronal threshold, showing that from the functional point of view, structure and elementary dynamics are interchangeable. We discuss possible geometric origins of this distribution, and comment on the importance of size, or of having a large number of neurons, in the culture.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Messenger RNA Fluctuations and Regulatory RNAs Shape the Dynamics of Negative Feedback Loop

Single cell experiments of simple regulatory networks can markedly differ from cell population experiments. Such differences arise from stochastic events in individual cells that are averaged out in cell populations. For instance, while individual cells may show sustained oscillations in the concentrations of some proteins, such oscillations may appear damped in the population average. In this paper we investigate the role of RNA stochastic fluctuations as a leading force to produce a sustained excitatory behavior at the single cell level. Opposed to some previous models, we build a fully stochastic model of a negative feedback loop that explicitly takes into account the RNA stochastic dynamics. We find that messenger RNA random fluctuations can be amplified during translation and produce sustained pulses of protein expression. Motivated by the recent appreciation of the importance of non--coding regulatory RNAs in post--transcription regulation, we also consider the possibility that a regulatory RNA transcript could bind to the messenger RNA and repress translation. Our findings show that the regulatory transcript helps reduce gene expression variability both at the single cell level and at the cell population level.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Molecular Model of the Contractile Ring

We present a model for the actin contractile ring of adherent animal cells. The model suggests that the actin concentration within the ring and consequently the power that the ring exerts both increase during contraction. We demonstrate the crucial role of actin polymerization and depolymerization throughout cytokinesis, and the dominance of viscous dissipation in the dynamics. The physical origin of two phases in cytokinesis dynamics ("biphasic cytokinesis") follows from a limitation on the actin density. The model is consistent with a wide range of measurements of the midzone of dividing animal cells.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Molecular Recognition as an Information Channel: The Role of Conformational Changes

Molecular recognition, which is essential in processing information in biological systems, takes place in a crowded noisy biochemical environment and requires the recognition of a specific target within a background of various similar competing molecules. We consider molecular recognition as a transmission of information via a noisy channel and use this analogy to gain insights on the optimal, or fittest, molecular recognizer. We focus on the optimal structural properties of the molecules such as flexibility and conformation. We show that conformational changes upon binding, which often occur during molecular recognition, may optimize the detection performance of the recognizer. We thus suggest a generic design principle termed 'conformational proofreading' in which deformation enhances detection. We evaluate the optimal flexibility of the molecular recognizer, which is analogous to the stochasticity in a decision unit. In some scenarios, a flexible recognizer, i.e., a stochastic decision unit, performs better than a rigid, deterministic one. As a biological example, we discuss conformational changes during homologous recombination, the process of genetic exchange between two DNA strands.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimal Design of a Molecular Recognizer: Molecular Recognition as a Bayesian Signal Detection Problem

Numerous biological functions-such as enzymatic catalysis, the immune response system, and the DNA-protein regulatory network-rely on the ability of molecules to specifically recognize target molecules within a large pool of similar competitors in a noisy biochemical environment. Using the basic framework of signal detection theory, we treat the molecular recognition process as a signal detection problem and examine its overall performance. Thus, we evaluate the optimal properties of a molecular recognizer in the presence of competition and noise. Our analysis reveals that the optimal design undergoes a "phase transition" as the structural properties of the molecules and interaction energies between them vary. In one phase, the recognizer should be complementary in structure to its target (like a lock and a key), while in the other, conformational changes upon binding, which often accompany molecular recognition, enhance recognition quality. Using this framework, the abundance of conformational changes may be explained as a result of increasing the fitness of the recognizer. Furthermore, this analysis may be used in future design of artificial signal processing devices based on biomolecules.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pearling in cells: A clue to understanding cell shape

Gradual disruption of the actin cytoskeleton induces a series of structural shape changes in cells leading to a transformation of cylindrical cell extensions into a periodic chain of "pearls". Quantitative measurements of the pearling instability give a square-root behavior for the wavelength as a function of drug concentration. We present a theory that explains these observations in terms of the interplay between rigidity of the submembranous actin shell and tension that is induced by boundary conditions set by adhesion points. The theory allows estimation of the rigidity and thickness of this supporting shell. The same theoretical considerations explain the shape of nonadherent edges in the general case of untreated cells.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Percolation Approach to Study Connectivity in Living Neural Networks

We study neural connectivity in cultures of rat hippocampal neurons. We measure the neurons' response to an electric stimulation for gradual lower connectivity, and characterize the size of the giant cluster in the network. The connectivity undergoes a percolation transition described by the critical exponent $β\simeq 0.65$. We use a theoretic approach based on bond.percolation on a graph to describe the process of disintegration of the network and extract its statistical properties. Together with numerical simulations we show that the connectivity in the neural culture is local, characterized by a gaussian degree distribution and not a power law on

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Percolation in living neural networks

We study living neural networks by measuring the neurons' response to a global electrical stimulation. Neural connectivity is lowered by reducing the synaptic strength, chemically blocking neurotransmitter receptors. We use a graph-theoretic approach to show that the connectivity undergoes a percolation transition. This occurs as the giant component disintegrates, characterized by a power law with critical exponent $β\simeq 0.65$ is independent of the balance between excitatory and inhibitory neurons and indicates that the degree distribution is gaussian rather than scale free

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phonons in a one-dimensional microfluidic crystal

The development of a general theoretical framework for describing the behaviour of a crystal driven far from equilibrium has proved difficult1. Microfluidic crystals, formed by the introduction of droplets of immiscible fluid into a liquid-filled channel, provide a convenient means to explore and develop models to describe non-equilibrium dynamics2, 3, 4, 5, 6, 7, 8, 9, 10, 11. Owing to the fact that these systems operate at low Reynolds number (Re), in which viscous dissipation of energy dominates inertial effects, vibrations are expected to be over-damped and contribute little to their dynamics12, 13, 14. Against such expectations, we report the emergence of collective normal vibrational modes (equivalent to acoustic 'phonons') in a one-dimensional microfluidic crystal of water-in-oil droplets at Reapprox10-4. These phonons propagate at an ultra-low sound velocity of approx100 mum s-1 and frequencies of a few hertz, exhibit unusual dispersion relations markedly different to those of harmonic crystals, and give rise to a variety of crystal instabilities that could have implications for the design of commercial microfluidic systems. First-principles theory shows that these phonons are an outcome of the symmetry-breaking flow field that induces long-range inter-droplet interactions, similar in nature to those observed in many other systems including dusty plasma crystals15, 16, vortices in superconductors17, 18, active membranes19 and nucleoprotein filaments20.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protein-DNA computation by stochastic assembly cascade

The assembly of RecA on single-stranded DNA is measured and interpreted as a stochastic finite-state machine that is able to discriminate fine differences between sequences, a basic computational operation. RecA filaments efficiently scan DNA sequence through a cascade of random nucleation and disassembly events that is mechanistically similar to the dynamic instability of microtubules. This iterative cascade is a multistage kinetic proofreading process that amplifies minute differences, even a single base change. Our measurements suggest that this stochastic Turing-like machine can compute certain integral transforms.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quorum Percolation in Living Neural Networks

Cooperative effects in neural networks appear because a neuron fires only if a minimal number $m$ of its inputs are excited. The multiple inputs requirement leads to a percolation model termed {\it quorum percolation}. The connectivity undergoes a phase transition as $m$ grows, from a network--spanning cluster at low $m$ to a set of disconnected clusters above a critical $m$. Both numerical simulations and the model reproduce the experimental results well. This allows a robust quantification of biologically relevant quantities such as the average connectivity $\kbar$ and the distribution of connections $p_k$

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

RecA-mediated homology search as a nearly optimal signal detection system

Homologous recombination facilitates the exchange of genetic material between homologous DNA molecules. This crucial process requires detecting a specific homologous DNA sequence within a huge variety of heterologous sequences. The detection is mediated by RecA in E. coli, or members of its superfamily in other organisms. Here we examine how well is the RecA-DNA interaction adjusted to its task. By formulating the DNA recognition process as a signal detection problem, we find the optimal value of binding energy that maximizes the ability to detect homologous sequences. We show that the experimentally observed binding energy is nearly optimal. This implies that the RecA-induced deformation and the binding energetics are fine-tuned to ensure optimal sequence detection. Our analysis suggests a possible role for DNA extension by RecA, in which deformation enhances detection. The present signal detection approach provides a general recipe for testing the optimality of other molecular recognition systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rules for biological regulation based on error minimization

The control of gene expression involves complex mechanisms that show large variation in design. For example, genes can be turned on either by the binding of an activator (positive control) or the unbinding of a repressor (negative control). What determines the choice of mode of control for each gene? This study proposes rules for gene regulation based on the assumption that free regulatory sites are exposed to nonspecific binding errors, whereas sites bound to their cognate regulators are protected from errors. Hence, the selected mechanisms keep the sites bound to their designated regulators for most of the time, thus minimizing fitness-reducing errors. This offers an explanation of the empirically demonstrated Savageau demand rule: Genes that are needed often in the natural environment tend to be regulated by activators, and rarely needed genes tend to be regulated by repressors; in both cases, sites are bound for most of the time, and errors are minimized. The fitness advantage of error minimization appears to be readily selectable. The present approach can also generate rules for multi-regulator systems. The error-minimization framework raises several experimentally testable hypotheses. It may also apply to other biological regulation systems, such as those involving protein-protein interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Screening by symmetry of long-range hydrodynamic interactions of polymers confined in sheets

Hydrodynamic forces may significantly affect the motion of polymers. In sheet-like cavities, such as the cell's cytoplasm and microfluidic channels, the hydrodynamic forces are long-range. It is therefore expected that that hydrodynamic interactions will dominate the motion of polymers in sheets and will be manifested by Zimm-like scaling. Quite the opposite, we note here that although the hydrodynamic forces are long-range their overall effect on the motion of polymers vanishes due to the symmetry of the two-dimensional flow. As a result, the predicted scaling of experimental observables such as the diffusion coefficient or the rotational diffusion time is Rouse-like, in accord with recent experiments. The effective screening validates the use of the non-interacting blobs picture for polymers confined in a sheet.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The physical language of molecular codes: A rate-distortion approach to the evolution and emergence of biological codes

The function of the organism hinges on the performance of its information-processing networks, which convey information via molecular recognition. Many paths within these networks utilize molecular codebooks, such as the genetic code, to translate information written in one class of molecules into another molecular "language" . The present paper examines the emergence and evolution of molecular codes in terms of rate-distortion theory and reviews recent results of this approach. We discuss how the biological problem of maximizing the fitness of an organism by optimizing its molecular coding machinery is equivalent to the communication engineering problem of designing an optimal information channel. The fitness of a molecular code takes into account the interplay between the quality of the channel and the cost of resources which the organism needs to invest in its construction and maintenance. We analyze the dynamics of a population of organisms that compete according to the fitness of their codes. The model suggests a generic mechanism for the emergence of molecular codes as a phase transition in an information channel. This mechanism is put into biological context and demonstrated in a simple example.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Physics of Living Neural Networks

Improvements in technique in conjunction with an evolution of the theoretical and conceptual approach to neuronal networks provide a new perspective on living neurons in culture. Organization and connectivity are being measured quantitatively along with other physical quantities such as information, and are being related to function. In this review we first discuss some of these advances, which enable elucidation of structural aspects. We then discuss two recent experimental models that yield some conceptual simplicity. A one-dimensional network enables precise quantitative comparison to analytic models, for example of propagation and information transport. A two-dimensional percolating network gives quantitative information on connectivity of cultured neurons. The physical quantities that emerge as essential characteristics of the network in vitro are propagation speeds, synaptic transmission, information creation and capacity. Potential application to neuronal devices is discussed.

0 citations0 reviewsNo code linkedOpen access