Motivated by growing evidence for pathway heterogeneity and alternative functions of molecular machines, we demonstrate a computational approach for investigating two questions: (1) Are there multiple mechanisms (state-space pathways) by which a machine can perform a given function, such as cotransport across a membrane? (2) How can additional functionality, such as proofreading/error-correction, be built into machine function using standard biochemical processes? Answers to these questions will aid both the understanding of molecular-scale cell biology and the design of synthetic machines. Focusing on transport in this initial study, we sample a variety of mechanisms by employing Metropolis Markov chain Monte Carlo. Trial moves adjust transition rates among an automatically generated set of conformational and binding states while maintaining fidelity to thermodynamic principles and a user-supplied fitness/functionality goal. Each accepted move generates a new model. The simulations yield both single and mixed reaction pathways for cotransport in a simple environment with a single substrate along with a driving ion. In a "competitive" environment including an additional decoy substrate, several qualitatively distinct reaction pathways are found which are capable of extremely high discrimination coupled to a leak of the driving ion, akin to proofreading. The array of functional models would be difficult to find by intuition alone in the complex state-spaces of interest.

Sorting and degradation of receptors and associated signaling molecules maintain homeostasis of conserved signaling pathways during cell specification and tissue development. Yet, whether machineries that sort signaling proteins act preferentially on different receptors and ligands in different contexts remains mysterious. Here, we show that Vacuolar protein sorting 25, Vps25, a component of ESCRT-II (Endosomal Sorting Complex Required for Transport II), directs preferential endosome-mediated modulation of FGF signaling in limbs. By ENU-induced mutagenesis, we isolated a polydactylous mouse line carrying a hypomorphic mutation of Vps25 (Vps25(ENU)). Unlike Vps25-null embryos we generated, Vps25(ENU/ENU) mutants survive until late gestation. Their limbs display FGF signaling enhancement and consequent hyperactivation of the FGF-SHH feedback loop causing polydactyly, whereas WNT and BMP signaling remain unperturbed. Notably, Vps25(ENU/ENU) Mouse Embryonic Fibroblasts exhibit aberrant FGFR trafficking and degradation; however, SHH signaling is unperturbed. These studies establish that the ESCRT-II machinery selectively limits FGF signaling in vertebrate skeletal patterning.

Multiple coronaviruses have emerged independently in the past 20 years that cause lethal human diseases. Although vaccine development targeting these viruses has been accelerated substantially, there remain patients requiring treatment who cannot be vaccinated or who experience breakthrough infections. Understanding the common host factors necessary for the life cycles of coronaviruses may reveal conserved therapeutic targets. Here, we used the known substrate specificities of mammalian protein kinases to deconvolute the sequence of phosphorylation events mediated by three host protein kinase families (SRPK, GSK-3, and CK1) that coordinately phosphorylate a cluster of serine and threonine residues in the viral N protein, which is required for viral replication. We also showed that loss or inhibition of SRPK1/2, which we propose initiates the N protein phosphorylation cascade, compromised the viral replication cycle. Because these phosphorylation sites are highly conserved across coronaviruses, inhibitors of these protein kinases not only may have therapeutic potential against COVID-19 but also may be broadly useful against coronavirus-mediated diseases.