Solid-body trajectoids shaped to roll along desired pathways.

In everyday life, rolling motion is typically associated with cylindrical (for example, car wheels) or spherical (for example, billiard balls) bodies tracing linear paths. However, mathematicians have, for decades, been interested in more exotically shaped solids such as the famous oloids1, sphericons2, polycons3, platonicons4 and two-circle rollers5 that roll downhill in curvilinear paths (in contrast to cylinders or spheres) yet indefinitely (in contrast to cones, Supplementary Video 1). The trajectories traced by such bodies have been studied in detail6-9, and can be useful in the context of efficient mixing10,11 and robotics, for example, in magnetically actuated, millimetre-sized sphericon-shaped robots12,13, or larger sphericon- and oloid-shaped robots translocating by shifting their centre of mass14,15. However, the rolling paths of these shapes are all sinusoid-like and their diversity ends there. Accordingly, we were intrigued whether a more general problem is solvable: given an infinite periodic trajectory, find the shape that would trace this trajectory when rolling down a slope. Here, we develop an algorithm to design such bodies-which we call 'trajectoids'-and then validate these designs experimentally by three-dimensionally printing the computed shapes and tracking their rolling paths, including those that close onto themselves such that the body's centre of mass moves intermittently uphill (Supplementary Video 2). Our study is motivated largely by fundamental curiosity, but the existence of trajectoids for most paths has unexpected implications for quantum and classical optics, as the dynamics of qubits, spins and light polarization can be exactly mapped to trajectoids and their paths16.

The ribosome as an optimal decoder: a lesson in molecular recognition.

The ribosome is a complex molecular machine that, in order to synthesize proteins, has to decode mRNAs by pairing their codons with matching tRNAs. Decoding is a major determinant of fitness and requires accurate and fast selection of correct tRNAs among many similar competitors. However, it is unclear whether the modern ribosome, and in particular its large conformational changes during decoding, are the outcome of adaptation to its task as a decoder or the result of other constraints. Here, we derive the energy landscape that provides optimal discrimination between competing substrates and thereby optimal tRNA decoding. We show that the measured landscape of the prokaryotic ribosome is sculpted in this way. This model suggests that conformational changes of the ribosome and tRNA during decoding are means to obtain an optimal decoder. Our analysis puts forward a generic mechanism that may be utilized broadly by molecular recognition systems.

Programmable Potentials Choreograph Defects in a Colloidal Crystal Shell.

Crystallization on spherical surfaces is obliged by topology to induce lattice defects. But controlling the organization of such defects remains a great challenge due to the long-range constraints of the curved geometry. Here, we report on DNA-coated colloids whose programmable interaction potentials can be used to regulate the arrangement of defects and even achieve perfect icosahedral order on a sphere. Combined simulations and theoretical analysis show how the potential can be tuned by changing the temperature, thereby controlling the number of defects. An explicit expression for the effective potential is derived, allowing us to distinguish the effects of entropic repulsion and enthalpic attraction. Altogether, the present findings provide insights into the physics of crystallization on curved spaces and may be used for designing desired crystal geometries.

Gene birth in a model of non-genic adaptation.

BACKGROUND: Over evolutionary timescales, genomic loci can switch between functional and non-functional states through processes such as pseudogenization and de novo gene birth. Particularly, de novo gene birth is a widespread process, and many examples continue to be discovered across diverse evolutionary lineages. However, the general mechanisms that lead to functionalization are poorly understood, and estimated rates of de novo gene birth remain contentious. Here, we address this problem within a model that takes into account mutations and structural variation, allowing us to estimate the likelihood of emergence of new functions at non-functional loci. RESULTS: Assuming biologically reasonable mutation rates and mutational effects, we find that functionalization of non-genic loci requires the realization of strict conditions. This is in line with the observation that most de novo genes are localized to the vicinity of established genes. Our model also provides an explanation for the empirical observation that emerging proto-genes are often lost despite showing signs of adaptation. CONCLUSIONS: Our work elucidates the properties of non-genic loci that make them fertile for adaptation, and our results offer mechanistic insights into the process of de novo gene birth.

General Theory of Specific Binding: Insights from a Genetic-Mechano-Chemical Protein Model.

Proteins need to selectively interact with specific targets among a multitude of similar molecules in the cell. However, despite a firm physical understanding of binding interactions, we lack a general theory of how proteins evolve high specificity. Here, we present such a model that combines chemistry, mechanics, and genetics and explains how their interplay governs the evolution of specific protein-ligand interactions. The model shows that there are many routes to achieving molecular discrimination-by varying degrees of flexibility and shape/chemistry complementarity-but the key ingredient is precision. Harder discrimination tasks require more collective and precise coaction of structure, forces, and movements. Proteins can achieve this through correlated mutations extending far from a binding site, which fine-tune the localized interaction with the ligand. Thus, the solution of more complicated tasks is enabled by increasing the protein size, and proteins become more evolvable and robust when they are larger than the bare minimum required for discrimination. The model makes testable, specific predictions about the role of flexibility and shape mismatch in discrimination, and how evolution can independently tune affinity and specificity. Thus, the proposed theory of specific binding addresses the natural question of "why are proteins so big?". A possible answer is that molecular discrimination is often a hard task best performed by adding more layers to the protein.

AlphaFold2 Can Predict Single-Mutation Effects.

AlphaFold2 (AF) is a promising tool, but is it accurate enough to predict single mutation effects? Here, we report that the localized structural deformation between protein pairs differing by only 1-3 mutations-as measured by the effective strain-is correlated across 3901 experimental and AF-predicted structures. Furthermore, analysis of ∼11 000 proteins shows that the local structural change correlates with various phenotypic changes. These findings suggest that AF can predict the range and magnitude of single-mutation effects on average, and we propose a method to improve precision of AF predictions and to indicate when predictions are unreliable.

Bona fide stochastic resonance under nonGaussian active fluctuations.

We report on the experimental observation of stochastic resonance (SR) in a nonGaussian active bath without any periodic modulation. A Brownian particle hopping in a nanoscale double-well potential under the influence of nonGaussian correlated noise, with mean interval τP and correlation time τc, shows a series of equally-spaced peaks in the residence time distribution at integral multiples of τP. The strength of the first peak is found to be maximum when the mean residence time d matches the double condition, 4τc ≈ τP ≈ d/2, demonstrating a new type of bona fide SR. The experimental findings agree with a simple model that explains the emergence of SR without periodic modulation of the double-well potential. Additionally, we show that generic SR under periodic modulation, known to degrade in strongly correlated continuous noise, is recovered by the discrete nonGaussian kicks.

Dimensional reduction in complex living systems: Where, why, and how.

The unprecedented prowess of measurement techniques provides a detailed, multi-scale look into the depths of living systems. Understanding these avalanches of high-dimensional data-by distilling underlying principles and mechanisms-necessitates dimensional reduction. We propose that living systems achieve exquisite dimensional reduction, originating from their capacity to learn, through evolution and phenotypic plasticity, the relevant aspects of a non-random, smooth physical reality. We explain how geometric insights by mathematicians allow one to identify these genuine hallmarks of life and distinguish them from universal properties of generic data sets. We illustrate these principles in a concrete example of protein evolution, suggesting a simple general recipe that can be applied to understand other biological systems.

Master curve of boosted diffusion for 10 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 coupled to molecular motion contrasts with generations of accumulated knowledge about diffusion at equilibrium. To test the limits of this idea, a critical testbed is the mobility of catalytically active enzymes. Sentiment is divided about the reality of enhanced enzyme diffusion, with evidence for and against. Here a master curve shows that the enzyme diffusion coefficient increases in proportion to the energy release rate-the product of Michaelis-Menten reaction rate and Gibbs free energy change (ΔG)-with a highly satisfactory correlation coefficient of 0.97. For 10 catalytic enzymes (urease, acetylcholinesterase, seven enzymes from the glucose cascade cycle, and one other), our measurements span from a roughly 40% enhanced diffusion coefficient at a high turnover rate and negative ΔG to no enhancement at a slow turnover rate and positive ΔG 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 of ΔG for the exergonic enzymes (ΔG <0), together with the vanishing effect for endergonic enzyme (ΔG >0), are consistent with a physical picture in which the mechanism boosting the diffusion is an active one, utilizing the available work from the chemical reaction.

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

Chaperone proteins-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 because of 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 amino-terminal affect a few specific positions in the RNA, located in stems and in the pseudoknot. In contrast, positive regions in the carboxy-terminal contain strongly coupled amino acids that promote nonspecific 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.

Gauging Nanoswimmer Dynamics via the Motion of Large Bodies.

Nanoswimmers are ubiquitous in biotechnology and nanotechnology but are extremely challenging to measure due to their minute size and driving forces. A simple method is proposed for detecting the elusive physical features of nanoswimmers by observing how they affect the motion of much larger, easily traceable particles. Modeling the swimmers as hydrodynamic force dipoles, we find direct, easy-to-calibrate relations between the observable power spectrum and diffusivity of the tracers and the dynamic characteristics of the swimmers-their force dipole moment and correlation times.

A comprehensive survey of developmental programs reveals a dearth of tree-like lineage graphs and ubiquitous regeneration.

BACKGROUND: Multicellular organisms are characterized by a wide diversity of forms and complexity despite a restricted set of key molecules and mechanisms at the base of organismal development. Development combines three basic processes-asymmetric cell division, signaling, and gene regulation-in a multitude of ways to create this overwhelming diversity of multicellular life forms. Here, we use a generative model to test the limits to which such processes can be combined to generate multiple differentiation paths during development, and attempt to chart the diversity of multicellular organisms generated. RESULTS: We sample millions of biologically feasible developmental schemes, allowing us to comment on the statistical properties of cell differentiation trajectories they produce. We characterize model-generated "organisms" using the graph topology of their cell type lineage maps. Remarkably, tree-type lineage differentiation maps are the rarest in our data. Additionally, a majority of the "organisms" generated by our model appear to be endowed with the ability to regenerate using pluripotent cells. CONCLUSIONS: 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. Furthermore, the high incidence of pluripotent cells in model-generated organisms stands in line with the long-standing hypothesis that whole body regeneration is an epiphenomenon of development. We discuss experimentally testable predictions of our model and some ways to adapt the generative framework to test additional hypotheses about general features of development.

A topological look into the evolution of developmental programs.

Rapid advance of experimental techniques provides an unprecedented in-depth view into complex developmental processes. Still, little is known on how the complexity of multicellular organisms evolved by elaborating developmental programs and inventing new cell types. A hurdle to understanding developmental evolution is the difficulty of even describing the intertwined network of spatiotemporal processes underlying the development of complex multicellular organisms. Nonetheless, an overview of developmental trajectories can be obtained from cell type lineage maps. Here, we propose that these lineage maps can also reveal how developmental programs evolve: the modes of evolving new cell types in an organism should be visible in its developmental trajectories and therefore in the geometry of its cell type lineage map. This idea is demonstrated using a parsimonious generative model of developmental programs, which allows us to reliably survey the universe of all possible programs and examine their topological features. We find that, contrary to belief, tree-like lineage maps are rare, and lineage maps of complex multicellular organisms are likely to be directed acyclic graphs in which multiple developmental routes can converge on the same cell type. Although cell type evolution prescribes what developmental programs come into existence, natural selection prunes those programs that produce low-functioning organisms. Our model indicates that additionally, lineage map topologies are correlated with such a functional property: the ability of organisms to regenerate.

Slowest-first protein translation scheme: Structural asymmetry and co-translational folding.

Proteins are translated from the N to the C terminus, raising the basic question of how this innate directionality affects their evolution. To explore this question, we analyze 16,200 structures from the Protein Data Bank (PDB). We find remarkable enrichment of α helices at the C terminus and ß strands at the N terminus. Furthermore, this α-ß asymmetry correlates with sequence length and contact order, both determinants of folding rate, hinting at possible links to co-translational folding (CTF). Hence, we propose the "slowest-first" scheme, whereby protein sequences evolved structural asymmetry to accelerate CTF: the slowest of the cooperatively folding segments are positioned near the N terminus so they have more time to fold during translation. A phenomenological model predicts that CTF can be accelerated by asymmetry in folding rate, up to double the rate, when folding time is commensurate with translation time; analysis of the PDB predicts that structural asymmetry is indeed maximal in this regime. This correspondence is greater in prokaryotes, which generally require faster protein production. Altogether, this indicates that accelerating CTF is a substantial evolutionary force whose interplay with stability and functionality is encoded in secondary structure asymmetry.

Catalytic enzymes are active matter.

Using a microscopic theory to analyze experiments, we demonstrate that enzymes are active matter. Superresolution fluorescence measurements-performed across four orders of magnitude of substrate concentration, with emphasis on the biologically relevant regime around or below the Michaelis-Menten constant-show that catalysis boosts the motion of enzymes to be superdiffusive for a few microseconds, enhancing their effective diffusivity over longer timescales. Occurring at the catalytic turnover rate, these fast ballistic leaps maintain direction over a duration limited by rotational diffusion, driving enzymes to execute wormlike trajectories by piconewton forces performing work of a few kBT against viscosity. The boosts are more frequent at high substrate concentrations, biasing the trajectories toward substrate-poor regions, thus exhibiting antichemotaxis, demonstrated here experimentally over a wide range of aqueous concentrations. Alternative noncatalytic, passive mechanisms that predict chemotaxis, cross-diffusion, and phoresis, are critically analyzed. We examine the physical interpretation of our findings, speculate on the underlying mechanism, and discuss the avenues they open with biological and technological implications. These findings violate the classical paradigm that chemical reaction and motility are distinct processes, and suggest reaction-motion coupling as a general principle of catalysis.

Enzyme leaps fuel antichemotaxis.

There is mounting evidence that enzyme diffusivity is enhanced when the enzyme is catalytically active. Here, using superresolution microscopy [stimulated emission-depletion fluorescence correlation spectroscopy (STED-FCS)], we show that active enzymes migrate spontaneously in the direction of lower substrate concentration ("antichemotaxis") by a process analogous to the run-and-tumble foraging strategy of swimming microorganisms and our theory quantifies the mechanism. The two enzymes studied, urease and acetylcholinesterase, display two families of transit times through subdiffraction-sized focus spots, a diffusive mode and a ballistic mode, and the latter transit time is close to the inverse rate of catalytic turnover. This biochemical information-processing algorithm may be useful to design synthetic self-propelled swimmers and nanoparticles relevant to active materials. Executed by molecules lacking the decision-making circuitry of microorganisms, antichemotaxis by this run-and-tumble process offers the biological function to homogenize product concentration, which could be significant in situations when the reactant concentration varies from spot to spot.

Green function of correlated genes in a minimal mechanical model of protein evolution.

The function of proteins arises from cooperative interactions and rearrangements of their amino acids, which exhibit large-scale dynamical modes. Long-range correlations have also been revealed in protein sequences, and this has motivated the search for physical links between the observed genetic and dynamic cooperativity. We outline here a simplified theory of protein, which relates sequence correlations to physical interactions and to the emergence of mechanical function. Our protein is modeled as a strongly coupled amino acid network with interactions and motions that are captured by the mechanical propagator, the Green function. The propagator describes how the gene determines the connectivity of the amino acids and thereby, the transmission of forces. Mutations introduce localized perturbations to the propagator that scatter the force field. The emergence of function is manifested by a topological transition when a band of such perturbations divides the protein into subdomains. We find that epistasis-the interaction among mutations in the gene-is related to the nonlinearity of the Green function, which can be interpreted as a sum over multiple scattering paths. We apply this mechanical framework to simulations of protein evolution and observe long-range epistasis, which facilitates collective functional modes.

Enhanced Diffusion and Oligomeric Enzyme Dissociation.

The concept that catalytic enzymes can act as molecular machines transducing chemical activity into motion has conceptual and experimental support, but experimental support has involved oligomeric enzymes, often studied under conditions where the substrate concentration is higher than biologically relevant and accordingly exceeds kM, the Michaelis constant. 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 dissociate above kM into their subunits that diffuse more rapidly, thus providing a simple physical mechanism that contributes to enhanced diffusion in this regime of concentrations. Mindful that this conclusion may be controversial, our findings are supported by four independent analytical techniques: static light scattering, dynamic light scattering (DLS), size-exclusion chromatography (SEC), and fluorescence correlation spectroscopy (FCS). Data for urease are emphasized and the conclusion is validated for hexokinase, acetylcholinesterase, and aldolase. For hexokinase and aldolase no enhanced diffusion is observed except under conditions when these oligomeric enzymes dissociate. At substrate concentration regimes below kM at which acetylcholinesterase and urease do not dissociate, our finding showing up to 10% enhancement of the diffusion coefficient is consistent with various theoretical scenarios in the literature.

Strain analysis of protein structures and low dimensionality of mechanical allosteric couplings.

In many proteins, especially allosteric proteins that communicate regulatory states from allosteric to active sites, structural deformations are functionally important. To understand these deformations, dynamical experiments are ideal but challenging. Using static structural information, although more limited than dynamical analysis, is much more accessible. Underused for protein analysis, strain is the natural quantity for studying local deformations. We calculate strain tensor fields for proteins deformed by ligands or thermal fluctuations using crystal and NMR structure ensembles. Strains-primarily shears-show deformations around binding sites. These deformations can be induced solely by ligand binding at distant allosteric sites. Shears reveal quasi-2D paths of mechanical coupling between allosteric and active sites that may constitute a widespread mechanism of allostery. We argue that strain-particularly shear-is the most appropriate quantity for analysis of local protein deformations. This analysis can reveal mechanical and biological properties of many proteins.

Lossless Brownian Information Engine.

We report on a lossless information engine that converts nearly all available information from an error-free feedback protocol into mechanical work. Combining high-precision detection at a resolution of 1 nm with ultrafast feedback control, the engine is tuned to extract the maximum work from information on the position of a Brownian particle. We show that the work produced by the engine achieves a bound set by a generalized second law of thermodynamics, demonstrating for the first time the sharpness of this bound. We validate a generalized Jarzynski equality for error-free feedback-controlled information engines.