Kinetic Proofreading.

Biochemistry and molecular biology rely on the recognition of structural complementarity between molecules. Molecular interactions must be both quickly reversible, i.e., tenuous, and specific. How the cell reconciles these conflicting demands is the subject of this article. The problem and its theoretical solution are discussed within the wider theoretical context of the thermodynamics of stochastic processes (stochastic thermodynamics). The solution-an irreversible reaction cycle that decreases internal error at the expense of entropy export into the environment-is shown to be widely employed by biological processes that transmit genetic and regulatory information.

Fundamental behaviors emerge from simulations of a living minimal cell.

We present a whole-cell fully dynamical kinetic model (WCM) of JCVI-syn3A, a minimal cell with a reduced genome of 493 genes that has retained few regulatory proteins or small RNAs. Cryo-electron tomograms provide the cell geometry and ribosome distributions. Time-dependent behaviors of concentrations and reaction fluxes from stochastic-deterministic simulations over a cell cycle reveal how the cell balances demands of its metabolism, genetic information processes, and growth, and offer insight into the principles of life for this minimal cell. The energy economy of each process including active transport of amino acids, nucleosides, and ions is analyzed. WCM reveals how emergent imbalances lead to slowdowns in the rates of transcription and translation. Integration of experimental data is critical in building a kinetic model from which emerges a genome-wide distribution of mRNA half-lives, multiple DNA replication events that can be compared to qPCR results, and the experimentally observed doubling behavior.

Antibody agonists trigger immune receptor signaling through local exclusion of receptor-type protein tyrosine phosphatases.

Antibodies can block immune receptor engagement or trigger the receptor machinery to initiate signaling. We hypothesized that antibody agonists trigger signaling by sterically excluding large receptor-type protein tyrosine phosphatases (RPTPs) such as CD45 from sites of receptor engagement. An agonist targeting the costimulatory receptor CD28 produced signals that depended on antibody immobilization and were sensitive to the sizes of the receptor, the RPTPs, and the antibody itself. Although both the agonist and a non-agonistic anti-CD28 antibody locally excluded CD45, the agonistic antibody was more effective. An anti-PD-1 antibody that bound membrane proximally excluded CD45, triggered Src homology 2 domain-containing phosphatase 2 recruitment, and suppressed systemic lupus erythematosus and delayed-type hypersensitivity in experimental models. Paradoxically, nivolumab and pembrolizumab, anti-PD-1-blocking antibodies used clinically, also excluded CD45 and were agonistic in certain settings. Reducing these agonistic effects using antibody engineering improved PD-1 blockade. These findings establish a framework for developing new and improved therapies for autoimmunity and cancer.

Size-Dependent Segregation Controls Macrophage Phagocytosis of Antibody-Opsonized Targets.

Macrophages protect the body from damage and disease by targeting antibody-opsonized cells for phagocytosis. Though antibodies can be raised against antigens with diverse structures, shapes, and sizes, it is unclear why some are more effective at triggering immune responses than others. Here, we define an antigen height threshold that regulates phagocytosis of both engineered and cancer-specific antigens by macrophages. Using a reconstituted model of antibody-opsonized target cells, we find that phagocytosis is dramatically impaired for antigens that position antibodies >10 nm from the target surface. Decreasing antigen height drives segregation of antibody-bound Fc receptors from the inhibitory phosphatase CD45 in an integrin-independent manner, triggering Fc receptor phosphorylation and promoting phagocytosis. Our work shows that close contact between macrophage and target is a requirement for efficient phagocytosis, suggesting that therapeutic antibodies should target short antigens in order to trigger Fc receptor activation through size-dependent physical segregation.

EGFR Ligands Differentially Stabilize Receptor Dimers to Specify Signaling Kinetics.

Epidermal growth factor receptor (EGFR) regulates many crucial cellular programs, with seven different activating ligands shaping cell signaling in distinct ways. Using crystallography and other approaches, we show how the EGFR ligands epiregulin (EREG) and epigen (EPGN) stabilize different dimeric conformations of the EGFR extracellular region. As a consequence, EREG or EPGN induce less stable EGFR dimers than EGF-making them partial agonists of EGFR dimerization. Unexpectedly, this weakened dimerization elicits more sustained EGFR signaling than seen with EGF, provoking responses in breast cancer cells associated with differentiation rather than proliferation. Our results reveal how responses to different EGFR ligands are defined by receptor dimerization strength and signaling dynamics. These findings have broad implications for understanding receptor tyrosine kinase (RTK) signaling specificity. Our results also suggest parallels between partial and/or biased agonism in RTKs and G-protein-coupled receptors, as well as new therapeutic opportunities for correcting RTK signaling output.

The guide-RNA sequence dictates the slicing kinetics and conformational dynamics of the Argonaute silencing complex.

The RNA-induced silencing complex (RISC), which powers RNA interference (RNAi), consists of a guide RNA and an Argonaute protein that slices target RNAs complementary to the guide. We find that, for different guide-RNA sequences, slicing rates of perfectly complementary bound targets can be surprisingly different (>250-fold range), and that faster slicing confers better knockdown in cells. Nucleotide sequence identities at guide-RNA positions 7, 10, and 17 underlie much of this variation in slicing rates. Analysis of one of these determinants implicates a structural distortion at guide nucleotides 6-7 in promoting slicing. Moreover, slicing directed by different guide sequences has an unanticipated, 600-fold range in 3'-mismatch tolerance, attributable to guides with weak (AU-rich) central pairing requiring extensive 3' complementarity (pairing beyond position 16) to more fully populate the slicing-competent conformation. Together, our analyses identify sequence determinants of RISC activity and provide biochemical and conformational rationale for their action.

Genome-wide quantification of RNA flow across subcellular compartments reveals determinants of the mammalian transcript life cycle.

Dissecting the regulatory mechanisms controlling mammalian transcripts from production to degradation requires quantitative measurements of mRNA flow across the cell. We developed subcellular TimeLapse-seq to measure the rates at which RNAs are released from chromatin, exported from the nucleus, loaded onto polysomes, and degraded within the nucleus and cytoplasm in human and mouse cells. These rates varied substantially, yet transcripts from genes with related functions or targeted by the same transcription factors and RNA-binding proteins flowed across subcellular compartments with similar kinetics. Verifying these associations uncovered a link between DDX3X and nuclear export. For hundreds of RNA metabolism genes, most transcripts with retained introns were degraded by the nuclear exosome, while the remaining molecules were exported with stable cytoplasmic lifespans. Transcripts residing on chromatin for longer had extended poly(A) tails, whereas the reverse was observed for cytoplasmic mRNAs. Finally, machine learning identified molecular features that predicted the diverse life cycles of mRNAs.

Systematic conformation-to-phenotype mapping via limited deep sequencing of proteins.

Non-native conformations drive protein-misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

E3 ligase autoinhibition by C-degron mimicry maintains C-degron substrate fidelity.

E3 ligase recruitment of proteins containing terminal destabilizing motifs (degrons) is emerging as a major form of regulation. How those E3s discriminate bona fide substrates from other proteins with terminal degron-like sequences remains unclear. Here, we report that human KLHDC2, a CRL2 substrate receptor targeting C-terminal Gly-Gly degrons, is regulated through interconversion between two assemblies. In the self-inactivated homotetramer, KLHDC2's C-terminal Gly-Ser motif mimics a degron and engages the substrate-binding domain of another protomer. True substrates capture the monomeric CRL2KLHDC2, driving E3 activation by neddylation and subsequent substrate ubiquitylation. Non-substrates such as NEDD8 bind KLHDC2 with high affinity, but its slow on rate prevents productive association with CRL2KLHDC2. Without substrate, neddylated CRL2KLHDC2 assemblies are deactivated via distinct mechanisms: the monomer by deneddylation and the tetramer by auto-ubiquitylation. Thus, substrate specificity is amplified by KLHDC2 self-assembly acting like a molecular timer, where only bona fide substrates may bind before E3 ligase inactivation.

Transcription activation is enhanced by multivalent interactions independent of phase separation.

Transcription factors (TFs) consist of a DNA-binding domain and an activation domain (AD) that are frequently considered to be independent and exchangeable modules. However, recent studies report that the physicochemical properties of the AD can control TF assembly at chromatin by driving phase separation into transcriptional condensates. Here, we dissected transcription activation by comparing different synthetic TFs at a reporter gene array with real-time single-cell fluorescence microscopy. In these experiments, binding site occupancy, residence time, and coactivator recruitment in relation to multivalent TF interactions were compared. While phase separation propensity and activation strength of the AD were linked, the actual formation of liquid-like TF droplets had a neutral or inhibitory effect on transcription activation. We conclude that multivalent AD-mediated interactions enhance the transcription activation capacity of a TF by increasing its residence time in the chromatin-bound state and facilitating the recruitment of coactivators independent of phase separation.

The dynamic mechanism of 4E-BP1 recognition and phosphorylation by mTORC1.

The activation of cap-dependent translation in eukaryotes requires multisite, hierarchical phosphorylation of 4E-BP by the 1 MDa kinase mammalian target of rapamycin complex 1 (mTORC1). To resolve the mechanism of this hierarchical phosphorylation at the atomic level, we monitored by NMR spectroscopy the interaction of intrinsically disordered 4E binding protein isoform 1 (4E-BP1) with the mTORC1 subunit regulatory-associated protein of mTOR (Raptor). The N-terminal RAIP motif and the C-terminal TOR signaling (TOS) motif of 4E-BP1 bind separate sites in Raptor, resulting in avidity-based tethering of 4E-BP1. This tethering orients the flexible central region of 4E-BP1 toward the mTORC1 kinase site for phosphorylation. The structural constraints imposed by the two tethering interactions, combined with phosphorylation-induced conformational switching of 4E-BP1, explain the hierarchy of 4E-BP1 phosphorylation by mTORC1. Furthermore, we demonstrate that mTORC1 recognizes both free and eIF4E-bound 4E-BP1, allowing rapid phosphorylation of the entire 4E-BP1 pool and efficient activation of translation. Finally, our findings provide a mechanistic explanation for the differential rapamycin sensitivity of the 4E-BP1 phosphorylation sites.

Human mitochondrial carriers of the SLC25 family function as monomers exchanging substrates with a ping-pong kinetic mechanism.

Members of the SLC25 mitochondrial carrier family link cytosolic and mitochondrial metabolism and support cellular maintenance and growth by transporting compounds across the mitochondrial inner membrane. Their monomeric or dimeric state and kinetic mechanism have been a matter of long-standing debate. It is believed by some that they exist as homodimers and transport substrates with a sequential kinetic mechanism, forming a ternary complex where both exchanged substrates are bound simultaneously. Some studies, in contrast, have provided evidence indicating that the mitochondrial ADP/ATP carrier (SLC25A4) functions as a monomer, has a single substrate binding site, and operates with a ping-pong kinetic mechanism, whereby ADP is imported before ATP is exported. Here we reanalyze the oligomeric state and kinetic properties of the human mitochondrial citrate carrier (SLC25A1), dicarboxylate carrier (SLC25A10), oxoglutarate carrier (SLC25A11), and aspartate/glutamate carrier (SLC25A13), all previously reported to be dimers with a sequential kinetic mechanism. We demonstrate that they are monomers, except for dimeric SLC25A13, and operate with a ping-pong kinetic mechanism in which the substrate import and export steps occur consecutively. These observations are consistent with a common transport mechanism, based on a functional monomer, in which a single central substrate-binding site is alternately accessible.

Phylogenic Determinants of Cardiovascular Frailty, Focus on Hemodynamics and Arterial Smooth Muscle Cells.

The evolution of the circulatory system from invertebrates to mammals has involved the passage from an open system to a closed in-parallel system via a closed in-series system, accompanying the increasing complexity and efficiency of life's biological functions. The archaic heart enables pulsatile motion waves of hemolymph in invertebrates, and the in-series circulation in fish occurs with only an endothelium, whereas mural smooth muscle cells appear later. The present review focuses on evolution of the circulatory system. In particular, we address how and why this evolution took place from a closed, flowing, longitudinal conductance at low pressure to a flowing, highly pressurized and bifurcating arterial compartment. However, although arterial pressure was the latest acquired hemodynamic variable, the general teleonomy of the evolution of species is the differentiation of individual organ function, supported by specific fueling allowing and favoring partial metabolic autonomy. This was achieved via the establishment of an active contractile tone in resistance arteries, which permitted the regulation of blood supply to specific organ activities via its localized function-dependent inhibition (active vasodilation). The global resistance to viscous blood flow is the peripheral increase in frictional forces caused by the tonic change in arterial and arteriolar radius, which backscatter as systemic arterial blood pressure. Consequently, the arterial pressure gradient from circulating blood to the adventitial interstitium generates the unidirectional outward radial advective conductance of plasma solutes across the wall of conductance arteries. This hemodynamic evolution was accompanied by important changes in arterial wall structure, supported by smooth muscle cell functional plasticity, including contractility, matrix synthesis and proliferation, endocytosis and phagocytosis, etc. These adaptive phenotypic shifts are due to epigenetic regulation, mainly related to mechanotransduction. These paradigms actively participate in cardio-arterial pathologies such as atheroma, valve disease, heart failure, aneurysms, hypertension, and physiological aging.

Acetate is a beneficial nutrient for E. coli at low glycolytic flux.

Acetate, a major by-product of glycolytic metabolism in Escherichia coli and many other microorganisms, has long been considered a toxic waste compound that inhibits microbial growth. This counterproductive auto-inhibition represents a major problem in biotechnology and has puzzled the scientific community for decades. Recent studies have however revealed that acetate is also a co-substrate of glycolytic nutrients and a global regulator of E. coli metabolism and physiology. Here, we used a systems biology strategy to investigate the mutual regulation of glycolytic and acetate metabolism in E. coli. Computational and experimental analyses demonstrate that decreasing the glycolytic flux enhances co-utilization of acetate with glucose. Acetate metabolism thus compensates for the reduction in glycolytic flux and eventually buffers carbon uptake so that acetate, rather than being toxic, actually enhances E. coli growth under these conditions. We validated this mechanism using three orthogonal strategies: chemical inhibition of glucose uptake, glycolytic mutant strains, and alternative substrates with a natively low glycolytic flux. In summary, acetate makes E. coli more robust to glycolytic perturbations and is a valuable nutrient, with a beneficial effect on microbial growth.

Prebiotic access to enantioenriched amino acids via peptide-mediated transamination reactions.

The kinetic resolution of racemic amino acids mediated by dipeptides and pyridoxal provides a prebiotically plausible route to enantioenriched proteinogenic amino acids. The enzymatic transamination cycles that are key to modern biochemical formation of enantiopure amino acids may have evolved from this half of the reversible reaction couple. Kinetic resolution of racemic precursors emerges as a general route to enantioenrichment under prebiotic conditions.

Temperature compensation through kinetic regulation in biochemical oscillators.

Nearly all circadian clocks maintain a period that is insensitive to temperature changes, a phenomenon known as temperature compensation (TC). Yet, it is unclear whether there is any common feature among different systems that exhibit TC. From a general timescale invariance, we show that TC relies on the existence of certain period-lengthening reactions wherein the period of the system increases strongly with the rates in these reactions. By studying several generic oscillator models, we show that this counterintuitive dependence is nonetheless a common feature of oscillators in the nonlinear (far-from-onset) regime where the oscillation can be separated into fast and slow phases. The increase of the period with the period-lengthening reaction rates occurs when the amplitude of the slow phase in the oscillation increases with these rates while the progression speed in the slow phase is controlled by other rates of the system. The positive dependence of the period on the period-lengthening rates balances its inverse dependence on other kinetic rates in the system, which gives rise to robust TC in a wide range of parameters. We demonstrate the existence of such period-lengthening reactions and their relevance for TC in all four model systems we considered. Theoretical results for a model of the Kai system are supported by experimental data. A study of the energy dissipation also shows that better TC performance requires higher energy consumption. Our study unveils a general mechanism by which a biochemical oscillator achieves TC by operating in parameter regimes far from the onset where period-lengthening reactions exist.

Discovering optimal kinetic pathways for self-assembly using automatic differentiation.

Macromolecular complexes are often composed of diverse subunits. The self-assembly of these subunits is inherently nonequilibrium and must avoid kinetic traps to achieve high yield over feasible timescales. We show how the kinetics of self-assembly benefits from diversity in subunits because it generates an expansive parameter space that naturally improves the "expressivity" of self-assembly, much like a deeper neural network. By using automatic differentiation algorithms commonly used in deep learning, we searched the parameter spaces of mass-action kinetic models to identify classes of kinetic protocols that mimic biological solutions for productive self-assembly. Our results reveal how high-yield complexes that easily become kinetically trapped in incomplete intermediates can instead be steered by internal design of rate-constants or external and active control of subunits to efficiently assemble. Internal design of a hierarchy of subunit binding rates generates self-assembly that can robustly avoid kinetic traps for all concentrations and energetics, but it places strict constraints on selection of relative rates. External control via subunit titration is more versatile, avoiding kinetic traps for any system without requiring molecular engineering of binding rates, albeit less efficiently and robustly. We derive theoretical expressions for the timescales of kinetic traps, and we demonstrate our optimization method applies not just for design but inference, extracting intersubunit binding rates from observations of yield-vs.-time for a heterotetramer. Overall, we identify optimal kinetic protocols for self-assembly as a powerful mechanism to achieve efficient and high-yield assembly in synthetic systems whether robustness or ease of "designability" is preferred.

Kinetic coevolutionary models predict the temporal emergence of HIV-1 resistance mutations under drug selection pressure.

Drug resistance in HIV type 1 (HIV-1) is a pervasive problem that affects the lives of millions of people worldwide. Although records of drug-resistant mutations (DRMs) have been extensively tabulated within public repositories, our understanding of the evolutionary kinetics of DRMs and how they evolve together remains limited. Epistasis, the interaction between a DRM and other residues in HIV-1 protein sequences, is key to the temporal evolution of drug resistance. We use a Potts sequence-covariation statistical-energy model of HIV-1 protein fitness under drug selection pressure, which captures epistatic interactions between all positions, combined with kinetic Monte-Carlo simulations of sequence evolutionary trajectories, to explore the acquisition of DRMs as they arise in an ensemble of drug-naive patient protein sequences. We follow the time course of 52 DRMs in the enzymes protease, RT, and integrase, the primary targets of antiretroviral therapy. The rates at which DRMs emerge are highly correlated with their observed acquisition rates reported in the literature when drug pressure is applied. This result highlights the central role of epistasis in determining the kinetics governing DRM emergence. Whereas rapidly acquired DRMs begin to accumulate as soon as drug pressure is applied, slowly acquired DRMs are contingent on accessory mutations that appear only after prolonged drug pressure. We provide a foundation for using computational methods to determine the temporal evolution of drug resistance using Potts statistical potentials, which can be used to gain mechanistic insights into drug resistance pathways in HIV-1 and other infectious agents.

Protein structure and aggregation: a marriage of necessity ruled by aggregation gatekeepers.

Protein aggregation propensity is a pervasive and seemingly inescapable property of proteomes. Strikingly, a significant fraction of the proteome is supersaturated, meaning that, for these proteins, their native conformation is less stable than the aggregated state. Maintaining the integrity of a proteome under such conditions is precarious and requires energy-consuming proteostatic regulation. Why then is aggregation propensity maintained at such high levels over long evolutionary timescales? Here, we argue that the conformational stability of the native and aggregated states are correlated thermodynamically and that codon usage strengthens this correlation. As a result, the folding of stable proteins requires kinetic control to avoid aggregation, provided by aggregation gatekeepers. These unique residues are evolutionarily selected to kinetically favor native folding, either on their own or by coopting chaperones.

RIG-I Uses an ATPase-Powered Translocation-Throttling Mechanism for Kinetic Proofreading of RNAs and Oligomerization.

RIG-I has a remarkable ability to specifically select viral 5'ppp dsRNAs for activation from a pool of cytosolic self-RNAs. The ATPase activity of RIG-I plays a role in RNA discrimination and activation, but the underlying mechanism was unclear. Using transient-state kinetics, we elucidated the ATPase-driven "kinetic proofreading" mechanism of RIG-I activation and RNA discrimination, akin to DNA polymerases, ribosomes, and T cell receptors. Even in the autoinhibited state of RIG-I, the C-terminal domain kinetically discriminates against self-RNAs by fast off rates. ATP binding facilitates dsRNA engagement but, interestingly, makes RIG-I promiscuous, explaining the constitutive signaling by Singleton-Merten syndrome-linked mutants that bind ATP without hydrolysis. ATP hydrolysis dissociates self-RNAs faster than 5'ppp dsRNA but, more importantly, drives RIG-I oligomerization through translocation, which we show to be regulated by helicase motif IVa. RIG-I translocates directionally from the dsRNA end into the stem region, and the 5'ppp end "throttles" translocation to provide a mechanism for threading and building a signaling-active oligomeric complex.