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

A linear reciprocal relationship between robustness and plasticity in homeostatic biological networks

In physics of living systems, a search for relationships of a few macroscopic variables that emerge from many microscopic elements is a central issue. We evolved gene regulatory networks so that the expression of core genes (partial system) is insensitive to environmental changes. Then, we found the expression levels of the remaining genes autonomously increase to provide a plastic (sensitive) response. A feedforward structure from the non-core to core genes evolved autonomously. Negative proportionality was observed between the average changes in core and non-core genes, reflecting reciprocity between the macroscopic robustness of homeostatic genes and plasticity of regulator genes. The proportion coefficient between those genes is represented by their number ratio, as in the "lever principle", whereas the decrease in the ratio results in a transition from perfect to partial adaptation, in which only a portion of the core genes exhibits robustness against environmental changes. This reciprocity between robustness and plasticity was satisfied throughout the evolutionary course, imposing an evolutionary constraint. This result suggests a simple macroscopic law for the adaptation characteristic in evolved complex biological networks.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

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.

Physical Constraints on Epistasis

Living systems evolve one mutation at a time, but a single mutation can alter the effect of subsequent mutations. The underlying mechanistic determinants of such epistasis are unclear. Here, we demonstrate that the physical dynamics of a biological system can generically constrain epistasis. We analyze models and experimental data on proteins and regulatory networks. In each, we find that if the long-time physical dynamics is dominated by a slow, collective mode, then the dimensionality of mutational effects is reduced. Consequently, epistatic coefficients for different combinations of mutations are no longer independent, even if individually strong. Such epistasis can be summarized as resulting from a global nonlinearity applied to an underlying linear trait, that is, as global epistasis. This constraint, in turn, reduces the ruggedness of the sequence-to-function map. By providing a generic mechanistic origin for experimentally observed global epistasis, our work suggests that slow collective physical modes can make biological systems evolvable.

Allostery and Epistasis: Emergent Properties of Anisotropic Networks

Understanding the underlying mechanisms behind protein allostery and non-additivity of substitution outcomes (i.e., epistasis) is critical when attempting to predict the functional impact of mutations, particularly at non-conserved sites. In an effort to model these two biological properties, we extend the framework of our metric to calculate dynamic coupling between residues, the Dynamic Coupling Index (DCI) to two new metrics: (i) EpiScore, which quantifies the difference between the residue fluctuation response of a functional site when two other positions are perturbed with random Brownian kicks simultaneously versus individually to capture the degree of cooperativity of these two other positions in modulating the dynamics of the functional site and (ii) DCI_asym, which measures the degree of asymmetry between the residue fluctuation response of two sites when one or the other is perturbed with a random force. Applied to four independent systems, we successfully show that EpiScore and DCI_asym can capture important biophysical properties in dual mutant substitution outcomes. We propose that allosteric regulation and the mechanisms underlying non-additive amino acid substitution outcomes (i.e., epistasis) can be understood as emergent properties of an anisotropic network of interactions where the inclusion of the full network of interactions is critical for accurate modeling. Consequently, mutations which drive towards a new function may require a fine balance between functional site asymmetry and strength of dynamic coupling with the functional sites. These two tools will provide mechanistic insight into both understanding and predicting the outcome of dual mutations.

Allosteric communication in molecular machines via information exchange: what can be learned from dynamical modeling

Allosteric regulation is crucial for the operation of protein machines and molecular motors. A major challenge is to characterize and quantify the information exchange underlying allosteric communication between remote functional sites in a protein, and to identify the involved relevant pathways. We review applications of two topical approaches of dynamical protein modeling, a kinetic-based single-molecule stochastic model, which employs information thermodynamics to quantify allosteric interactions, and structure-based coarse-grained modeling to characterize intra-molecular couplings in terms of conformational motions and propagating mechanical strain. Both descriptions resolve the directionality of allosteric responses within a protein, emphasizing the concept of causality as the principal hallmark of protein allostery. We discuss the application of techniques from information thermodynamics to dynamic protein elastic networks and evolutionary designed model structures, and the ramifications for protein allostery.

Direct coupling analysis of epistasis in allosteric materials

In allosteric proteins, the binding of a ligand modifies function at a distant active site. Such allosteric pathways can be used as target for drug design, generating considerable interest in inferring them from sequence alignment data. Currently, different methods lead to conflicting results, in particular on the existence of long-range evolutionary couplings between distant amino-acids mediating allostery. Here we propose a resolution of this conundrum, by studying epistasis and its inference in models where an allosteric material is evolved in silico to perform a mechanical task. We find in our model the four types of epistasis (Synergistic, Sign, Antagonistic, Saturation), which can be both short or long-range and have a simple mechanical interpretation. We perform a Direct Coupling Analysis (DCA) and find that DCA predicts well the cost of point mutations but is a rather poor generative model. Strikingly, it can predict short-range epistasis but fails to capture long-range epistasis, in consistence with empirical findings. We propose that such failure is generic when function requires subparts to work in concert. We illustrate this idea with a simple model, which suggests that other methods may be better suited to capture long-range effects.

Allostery in proteins is the property of highly specific responses to ligand binding at a distant site. To inform protocols of de novo drug design, it is fundamental to understand the impact of mutations on allosteric regulation and whether it can be predicted from evolutionary correlations. In this work we consider allosteric architectures artificially evolved to optimize the cooperativity of binding at allosteric and active site. We first characterize the emergent pattern of epistasis as well as the underlying mechanical phenomena, finding the four types of epistasis (Synergistic, Sign, Antagonistic, Saturation), which can be both short or long-range. The numerical evolution of these allosteric architectures allows us to benchmark Direct Coupling Analysis, a method which relies on co-evolution in sequence data to infer direct evolutionary couplings, in connection to allostery. We show that Direct Coupling Analysis predicts quantitatively point mutation costs but underestimates strong long-range epistasis. We provide an argument, based on a simplified model, illustrating the reasons for this discrepancy. Our analysis suggests neural networks as more promising tool to measure epistasis.

Mechanics of Allostery: Contrasting the Induced Fit and Population Shift Scenarios

In allosteric proteins, binding a ligand can affect function at a distant location, for example, by changing the binding affinity of a substrate at the active site. The induced fit and population shift models, which differ by the assumed number of stable configurations, explain such cooperative binding from a thermodynamic viewpoint. Yet, understanding what mechanical principles constrain these models remains a challenge. Here, we provide an empirical study on 34 proteins supporting the idea that allosteric conformational change generally occurs along a soft elastic mode presenting extended regions of high shear. We argue, based on a detailed analysis of how the energy profile along such a mode depends on binding, that in the induced fit scenario, there is an optimal stiffness ka∗ ∼ 1/N for cooperative binding, where N is the number of residues. We find that the population shift scenario is more robust to mutations affecting stiffness because binding becomes more and more cooperative with stiffness up to the same characteristic value ka∗, beyond which cooperativity saturates instead of decaying. We numerically confirm these findings in a nonlinear mechanical model. Dynamical considerations suggest that a stiffness of order ka∗ is favorable in that scenario as well, supporting that for proper function, proteins must evolve a functional elastic mode that is softer as their size increases. In consistency with this view, we find a fair anticorrelation between the stiffness of the allosteric response and protein size in our data set.

Simple mechanics of protein machines

While belonging to the nanoscale, protein machines are so complex that tracing even a small fraction of their cycle requires weeks of calculations on supercomputers. Surprisingly, many aspects of their operation can be however already reproduced by using very simple mechanical models of elastic networks. The analysis suggests that, similar to other self-organized complex systems, functional collective dynamics in such proteins is effectively reduced to a low-dimensional attractive manifold.

Revealing evolutionary constraints on proteins through sequence analysis

Statistical analysis of alignments of large numbers of protein sequences has revealed "sectors" of collectively coevolving amino acids in several protein families. Here, we show that selection acting on any functional property of a protein, represented by an additive trait, can give rise to such a sector. As an illustration of a selected trait, we consider the elastic energy of an important conformational change within an elastic network model, and we show that selection acting on this energy leads to correlations among residues. For this concrete example and more generally, we demonstrate that the main signature of functional sectors lies in the small-eigenvalue modes of the covariance matrix of the selected sequences. However, secondary signatures of these functional sectors also exist in the extensively-studied large-eigenvalue modes. Our simple, general model leads us to propose a principled method to identify functional sectors, along with the magnitudes of mutational effects, from sequence data. We further demonstrate the robustness of these functional sectors to various forms of selection, and the robustness of our approach to the identification of multiple selected traits.

Catalytic enzymes are active matter

Catalysis and mobility of reactants in fluid are normally thought to be decoupled. Violating this classical paradigm, this paper presents the catalyst laws of motion. Comparing experimental data to the theory presented here, we conclude that part of the free energy released by chemical reaction is channeled into driving catalysts to execute wormlike trajectories by piconewton forces performing work of a few k_BT against fluid viscosity, where the rotational diffusion rate dictates the trajectory persistence length. This active motion agitates the fluid medium and produces antichemotaxis, the migration of catalyst down the gradient of the reactant concentration. Alternative explanations of enhanced catalyst mobility are examined critically.

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 k_BT 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.

Designed Elastic Networks: Models of Complex Protein Machinery

Recently, the design of mechanical networks with protein-inspired responses has become increasingly popular. Here, we review contributions which were motivated by studies of protein dynamics employing coarse-grained elastic network models. First, the concept of evolutionary optimization that we developed to design network structures which execute prescribed tasks is explained. We then review what presumably marks the origin of the idea to design complex functional networks which encode protein-inspired behavior, namely the design of an elastic network structure which emulates the cycles of ATP-powered conformational motion in protein machines. Two recent applications are reviewed. First, the construction of a model molecular motor, whose operation incorporates both the tight coupling power stroke as well as the loose coupling Brownian ratchet mechanism, is discussed. Second, the evolutionary design of network structures which encode optimal long-range communication between remote sites and represent mechanical models of allosteric proteins is presented. We discuss the prospects of designed protein-mimicking elastic networks as model systems to elucidate the design principles and functional signatures underlying the operation of complex protein machinery.