Bioenergetic costs and the evolution of noise regulation by microRNAs.

Noise control, together with other regulatory functions facilitated by microRNAs (miRNAs), is believed to have played important roles in the evolution of multicellular eukaryotic organisms. miRNAs can dampen protein fluctuations via enhanced degradation of messenger RNA (mRNA), but this requires compensation by increased mRNA transcription to maintain the same expression levels. The overall mechanism is metabolically expensive, leading to questions about how it might have evolved in the first place. We develop a stochastic model of miRNA noise regulation, coupled with a detailed analysis of the associated metabolic costs. Additionally, we calculate binding free energies for a range of miRNA seeds, the short sequences which govern target recognition. We argue that natural selection may have fine-tuned the Michaelis-Menten constant [Formula: see text] describing miRNA-mRNA affinity and show supporting evidence from analysis of experimental data. [Formula: see text] is constrained by seed length, and optimal noise control (minimum protein variance at a given energy cost) is achievable for seeds of 6 to 7 nucleotides in length, the most commonly observed types. Moreover, at optimality, the degree of noise reduction approaches the theoretical bound set by the Wiener-Kolmogorov linear filter. The results illustrate how selective pressure toward energy efficiency has potentially shaped a crucial regulatory pathway in eukaryotes.

Diverse mutant selection windows shape spatial heterogeneity in evolving populations.

Mutant selection windows (MSWs), the range of drug concentrations that select for drug-resistant mutants, have long been used as a model for predicting drug resistance and designing optimal dosing strategies in infectious disease. The canonical MSW model offers comparisons between two subtypes at a time: drug-sensitive and drug-resistant. In contrast, the fitness landscape model with N alleles, which maps genotype to fitness, allows comparisons between N genotypes simultaneously, but does not encode continuous drug response data. In clinical settings, there may be a wide range of drug concentrations selecting for a variety of genotypes in both cancer and infectious diseases. Therefore, there is a need for a more robust model of the pathogen response to therapy to predict resistance and design new therapeutic approaches. Fitness seascapes, which model genotype-by-environment interactions, permit multiple MSW comparisons simultaneously by encoding genotype-specific dose-response data. By comparing dose-response curves, one can visualize the range of drug concentrations where one genotype is selected over another. In this work, we show how N-allele fitness seascapes allow for N * 2N-1 unique MSW comparisons. In spatial drug diffusion models, we demonstrate how fitness seascapes reveal spatially heterogeneous MSWs, extending the MSW model to more fully reflect the selection of drug resistant genotypes. Furthermore, using synthetic data and empirical dose-response data in cancer, we find that the spatial structure of MSWs shapes the evolution of drug resistance in an agent-based model. By simulating a tumor treated with cyclic drug therapy, we find that mutant selection windows introduced by drug diffusion promote the proliferation of drug resistant cells. Our work highlights the importance and utility of considering dose-dependent fitness seascapes in evolutionary medicine.

Catch bonds in sickle cell disease: Shear-enhanced adhesion of red blood cells to laminin.

Could the phenomenon of catch bonding-force-strengthened cellular adhesion-play a role in sickle cell disease, where abnormal red blood cell (RBC) adhesion obstructs blood flow? Here, we investigate the dynamics of sickle RBCs adhering to a surface functionalized with the protein laminin (a component of the extracellular matrix around blood vessels) under physiologically relevant microscale flow. First, using total internal reflectance microscopy we characterize the spatial fluctuations of the RBC membrane above the laminin surface before detachment. The complex dynamics we observe suggest the possibility of catch bonding, where the mean detachment time of the cell from the surface initially increases to a maximum and then decreases as a function of shear force. We next conduct a series of shear-induced detachment experiments on blood samples from 25 sickle cell disease patients, quantifying the number and duration of adhered cells under both sudden force jumps and linear force ramps. The experiments reveal that a subset of patients does indeed exhibit catch bonding. By fitting the data to a theoretical model of the bond dynamics, we can extract the mean bond lifetime versus force for each patient. The results show a striking heterogeneity among patients, both in terms of the qualitative behavior (whether or not there is catch bonding) and in the magnitudes of the lifetimes. Patients with large bond lifetimes at physiological forces are more likely to have certain adverse clinical features, like a diagnosis of pulmonary arterial hypertension and intracardiac shunts. By introducing an in vitro platform for fully characterizing RBC-laminin adhesion dynamics, our approach could contribute to the development of patient-specific antiadhesive therapies for sickle cell disease. The experimental setup is also easily generalizable to studying adhesion dynamics in other cell types, for example, leukocytes or cancer cells, and can incorporate disease-relevant environmental conditions like oxygen deprivation.

Cheater suppression and stochastic clearance through quorum sensing.

The evolutionary consequences of quorum sensing in regulating bacterial cooperation are not fully understood. In this study, we reveal unexpected effects of regulating public good production through quorum sensing on bacterial population dynamics, showing that quorum sensing can be a collectively harmful alternative to unregulated production. We analyze a birth-death model of bacterial population dynamics accounting for public good production and the presence of non-producing cheaters. Our model demonstrates that when demographic noise is a factor, the consequences of controlling public good production according to quorum sensing depend on the cost of public good production and the growth rate of populations in the absence of public goods. When public good production is inexpensive, quorum sensing is a destructive alternative to unconditional production, in terms of the mean population extinction time. When costs are higher, quorum sensing becomes a constructive strategy for the producing strain, both stabilizing cooperation and decreasing the risk of population extinction.

Integrating deep learning with microfluidics for biophysical classification of sickle red blood cells adhered to laminin.

Sickle cell disease, a genetic disorder affecting a sizeable global demographic, manifests in sickle red blood cells (sRBCs) with altered shape and biomechanics. sRBCs show heightened adhesive interactions with inflamed endothelium, triggering painful vascular occlusion events. Numerous studies employ microfluidic-assay-based monitoring tools to quantify characteristics of adhered sRBCs from high resolution channel images. The current image analysis workflow relies on detailed morphological characterization and cell counting by a specially trained worker. This is time and labor intensive, and prone to user bias artifacts. Here we establish a morphology based classification scheme to identify two naturally arising sRBC subpopulations-deformable and non-deformable sRBCs-utilizing novel visual markers that link to underlying cell biomechanical properties and hold promise for clinically relevant insights. We then set up a standardized, reproducible, and fully automated image analysis workflow designed to carry out this classification. This relies on a two part deep neural network architecture that works in tandem for segmentation of channel images and classification of adhered cells into subtypes. Network training utilized an extensive data set of images generated by the SCD BioChip, a microfluidic assay which injects clinical whole blood samples into protein-functionalized microchannels, mimicking physiological conditions in the microvasculature. Here we carried out the assay with the sub-endothelial protein laminin. The machine learning approach segmented the resulting channel images with 99.1±0.3% mean IoU on the validation set across 5 k-folds, classified detected sRBCs with 96.0±0.3% mean accuracy on the validation set across 5 k-folds, and matched trained personnel in overall characterization of whole channel images with R2 = 0.992, 0.987 and 0.834 for total, deformable and non-deformable sRBC counts respectively. Average analysis time per channel image was also improved by two orders of magnitude (∼ 2 minutes vs ∼ 2-3 hours) over manual characterization. Finally, the network results show an order of magnitude less variance in counts on repeat trials than humans. This kind of standardization is a prerequisite for the viability of any diagnostic technology, making our system suitable for affordable and high throughput disease monitoring.

The 2019 mathematical oncology roadmap.

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

Modeling the Growth of Organisms Validates a General Relation between Metabolic Costs and Natural Selection.

Metabolism and evolution are closely connected: if a mutation incurs extra energetic costs for an organism, there is a baseline selective disadvantage that may or may not be compensated for by other adaptive effects. A long-standing, but to date unproven, hypothesis is that this disadvantage is equal to the fractional cost relative to the total resting metabolic expenditure. We validate this result from physical principles through a general growth model and show it holds to excellent approximation for experimental parameters drawn from a wide range of species.

Unraveling the mechanism of the cadherin-catenin-actin catch bond.

The adherens junctions between epithelial cells involve a protein complex formed by E-cadherin, ß-catenin, α-catenin and F-actin. The stability of this complex was a puzzle for many years, since in vitro studies could reconstitute various stable subsets of the individual proteins, but never the entirety. The missing ingredient turned out to be mechanical tension: a recent experiment that applied physiological forces to the complex with an optical tweezer dramatically increased its lifetime, a phenomenon known as catch bonding. However, in the absence of a crystal structure for the full complex, the microscopic details of the catch bond mechanism remain mysterious. Building on structural clues that point to α-catenin as the force transducer, we present a quantitative theoretical model for how the catch bond arises, fully accounting for the experimental lifetime distributions. The underlying hypothesis is that force induces a rotational transition between two conformations of α-catenin, overcoming a significant energy barrier due to a network of salt bridges. This transition allosterically regulates the energies at the interface between α-catenin and F-actin. The model allows us to predict these energetic changes, as well as highlighting the importance of the salt bridge rotational barrier. By stabilizing one of the α-catenin states, this barrier could play a role in how the complex responds to additional in vivo binding partners like vinculin. Since significant conformational energy barriers are a common feature of other adhesion systems that exhibit catch bonds, our model can be adapted into a general theoretical framework for integrating structure and function in a variety of force-regulated protein complexes.

Directly measuring single-molecule heterogeneity using force spectroscopy.

One of the most intriguing results of single-molecule experiments on proteins and nucleic acids is the discovery of functional heterogeneity: the observation that complex cellular machines exhibit multiple, biologically active conformations. The structural differences between these conformations may be subtle, but each distinct state can be remarkably long-lived, with interconversions between states occurring only at macroscopic timescales, fractions of a second or longer. Although we now have proof of functional heterogeneity in a handful of systems-enzymes, motors, adhesion complexes-identifying and measuring it remains a formidable challenge. Here, we show that evidence of this phenomenon is more widespread than previously known, encoded in data collected from some of the most well-established single-molecule techniques: atomic force microscopy or optical tweezer pulling experiments. We present a theoretical procedure for analyzing distributions of rupture/unfolding forces recorded at different pulling speeds. This results in a single parameter, quantifying the degree of heterogeneity, and also leads to bounds on the equilibration and conformational interconversion timescales. Surveying 10 published datasets, we find heterogeneity in 5 of them, all with interconversion rates slower than 10 s(-1) Moreover, we identify two systems where additional data at realizable pulling velocities is likely to find a theoretically predicted, but so far unobserved crossover regime between heterogeneous and nonheterogeneous behavior. The significance of this regime is that it will allow far more precise estimates of the slow conformational switching times, one of the least understood aspects of functional heterogeneity.

Force-dependent switch in protein unfolding pathways and transition-state movements.

Although it is known that single-domain proteins fold and unfold by parallel pathways, demonstration of this expectation has been difficult to establish in experiments. Unfolding rate, [Formula: see text], as a function of force f, obtained in single-molecule pulling experiments on src SH3 domain, exhibits upward curvature on a [Formula: see text] plot. Similar observations were reported for other proteins for the unfolding rate [Formula: see text]. These findings imply unfolding in these single-domain proteins involves a switch in the pathway as f or [Formula: see text] is increased from a low to a high value. We provide a unified theory demonstrating that if [Formula: see text] as a function of a perturbation (f or [Formula: see text]) exhibits upward curvature then the underlying energy landscape must be strongly multidimensional. Using molecular simulations we provide a structural basis for the switch in the pathways and dramatic shifts in the transition-state ensemble (TSE) in src SH3 domain as f is increased. We show that a single-point mutation shifts the upward curvature in [Formula: see text] to a lower force, thus establishing the malleability of the underlying folding landscape. Our theory, applicable to any perturbation that affects the free energy of the protein linearly, readily explains movement in the TSE in a ß-sandwich (I27) protein and single-chain monellin as the denaturant concentration is varied. We predict that in the force range accessible in laser optical tweezer experiments there should be a switch in the unfolding pathways in I27 or its mutants.

Phenomenological and microscopic theories for catch bonds.

Lifetimes of bound states of protein complexes or biomolecule folded states typically decrease when subject to mechanical force. However, a plethora of biological systems exhibit the counter-intuitive phenomenon of catch bonding, where non-covalent bonds become stronger under externally applied forces. The quest to understand the origin of catch-bond behavior has led to the development of phenomenological and microscopic theories that can quantitatively recapitulate experimental data. Here, we assess the successes and limitations of such theories in explaining experimental data. The most widely applied approach is a phenomenological two-state model, which fits all of the available data on a variety of complexes: actomyosin, kinetochore-microtubule, selectin-ligand, and cadherin-catenin binding to filamentous actin. With a primary focus on the selectin family of cell-adhesion complexes, we discuss the positives and negatives of phenomenological models and the importance of evaluating the physical relevance of fitting parameters. We describe a microscopic theory for selectins, which provides a structural basis for catch bonds and predicts a crucial allosteric role for residues Asn82-Glu88. We emphasize the need for new theories and simulations that can mimic experimental conditions, given the complex response of cell adhesion complexes to force and their potential role in a variety of biological contexts.

Extreme sensitivity biosensing platform based on hyperbolic metamaterials.

Optical sensor technology offers significant opportunities in the field of medical research and clinical diagnostics, particularly for the detection of small numbers of molecules in highly diluted solutions. Several methods have been developed for this purpose, including label-free plasmonic biosensors based on metamaterials. However, the detection of lower-molecular-weight (<500 Da) biomolecules in highly diluted solutions is still a challenging issue owing to their lower polarizability. In this context, we have developed a miniaturized plasmonic biosensor platform based on a hyperbolic metamaterial that can support highly confined bulk plasmon guided modes over a broad wavelength range from visible to near infrared. By exciting these modes using a grating-coupling technique, we achieved different extreme sensitivity modes with a maximum of 30,000 nm per refractive index unit (RIU) and a record figure of merit (FOM) of 590. We report the ability of the metamaterial platform to detect ultralow-molecular-weight (244 Da) biomolecules at picomolar concentrations using a standard affinity model streptavidin-biotin.

Protein collapse is encoded in the folded state architecture.

Folded states of single domain globular proteins are compact with high packing density. The radius of gyration, Rg, of both the folded and unfolded states increase as Nν where N is the number of amino acids in the protein. The values of the Flory exponent ν are, respectively, ≈⅓ and ≈0.6 in the folded and unfolded states, coinciding with those for homopolymers. However, the extent of compaction of the unfolded state of a protein under low denaturant concentration (collapsibility), conditions favoring the formation of the folded state, is unknown. We develop a theory that uses the contact map of proteins as input to quantitatively assess collapsibility of proteins. Although collapsibility is universal, the propensity to be compact depends on the protein architecture. Application of the theory to over two thousand proteins shows that collapsibility depends not only on N but also on the contact map reflecting the native structure. A major prediction of the theory is that ß-sheet proteins are far more collapsible than structures dominated by α-helices. The theory and the accompanying simulations, validating the theoretical predictions, provide insights into the differing conclusions reached using different experimental probes assessing the extent of compaction of proteins. By calculating the criterion for collapsibility as a function of protein length we provide quantitative insights into the reasons why single domain proteins are small and the physical reasons for the origin of multi-domain proteins. Collapsibility of non-coding RNA molecules is similar ß-sheet proteins structures adding support to "Compactness Selection Hypothesis".

Plasticity of hydrogen bond networks regulates mechanochemistry of cell adhesion complexes.

Mechanical forces acting on cell adhesion receptor proteins regulate a range of cellular functions by formation and rupture of noncovalent interactions with ligands. Typically, force decreases the lifetimes of intact complexes ("slip bonds"), making the discovery that these lifetimes can also be prolonged ("catch bonds") a surprise. We created a microscopic analytic theory by incorporating the structures of selectin and integrin receptors into a conceptual framework based on the theory of stochastic equations, which quantitatively explains a wide range of experimental data (including catch bonds at low forces and slip bonds at high forces). Catch bonds arise due to force-induced remodeling of hydrogen bond networks, a finding that also accounts for unbinding in structurally unrelated integrin-fibronectin and actomyosin complexes. For the selectin family, remodeling of hydrogen bond networks drives an allosteric transition resulting in the formation of the maximum number of hydrogen bonds determined only by the structure of the receptor and independent of the ligand. A similar transition allows us to predict the increase in the number of hydrogen bonds in a particular allosteric state of α5ß1 integrin-fibronectin complex, a conformation which is yet to be crystallized. We also make a testable prediction that a single point mutation (Tyr51Phe) in the ligand associated with selectin should dramatically alter the nature of the catch bond compared with the wild type. Our work suggests that nature uses a ductile network of hydrogen bonds to engineer function over a broad range of forces.

Sequence-resolved free energy profiles of stress-bearing vimentin intermediate filaments.

Intermediate filaments (IFs) are key to the mechanical strength of metazoan cells. Their basic building blocks are dimeric coiled coils mediating hierarchical assembly of the full-length filaments. Here we use single-molecule force spectroscopy by optical tweezers to assess the folding and stability of coil 2B of the model IF protein vimentin. The coiled coil was unzipped from its N and C termini. When pulling from the C terminus, we observed that the coiled coil was resistant to force owing to the high stability of the C-terminal region. Pulling from the N terminus revealed that the N-terminal half is considerably less stable. The mechanical pulling assay is a unique tool to study and control seed formation and structure propagation of the coiled coil. We then used rigorous theory-based deconvolution for a model-free extraction of the energy landscape and local stability profiles. The data obtained from the two distinct pulling directions complement each other and reveal a tripartite stability of the coiled coil: a labile N-terminal half, followed by a medium stability section and a highly stable region at the far C-terminal end. The different stability regions provide important insight into the mechanics of IF assembly.

Discrete Step Sizes of Molecular Motors Lead to Bimodal Non-Gaussian Velocity Distributions under Force.

Fluctuations in the physical properties of biological machines are inextricably linked to their functions. Distributions of run lengths and velocities of processive molecular motors, like kinesin-1, are accessible through single-molecule techniques, but rigorous theoretical models for these probabilities are lacking. Here, we derive exact analytic results for a kinetic model to predict the resistive force (F)-dependent velocity [P(v)] and run length [P(n)] distribution functions of generic finitely processive molecular motors. Our theory quantitatively explains the zero force kinesin-1 data for both P(n) and P(v) using the detachment rate as the only parameter. In addition, we predict the F dependence of these quantities. At nonzero F, P(v) is non-Gaussian and is bimodal with peaks at positive and negative values of v, which is due to the discrete step size of kinesin-1. Although the predictions are based on analyses of kinesin-1 data, our results are general and should hold for any processive motor, which walks on a track by taking discrete steps.

Design principles governing the motility of myosin V.

The molecular motor myosin V (MyoV) exhibits a wide repertoire of pathways during the stepping process, which is intimately connected to its biological function. The best understood of these is the hand-over-hand stepping by a swinging lever arm movement toward the plus end of actin filaments. Single-molecule experiments have also shown that the motor "foot stomps," with one hand detaching and rebinding to the same site, and back-steps under sufficient load. The complete taxonomy of MyoV's load-dependent stepping pathways, and the extent to which these are constrained by motor structure and mechanochemistry, are not understood. Using a polymer model, we develop an analytical theory to describe the minimal physical properties that govern motor dynamics. We solve the first-passage problem of the head reaching the target-binding site, investigating the competing effects of backward load, strain in the leading head biasing the diffusion in the direction of the target, and the possibility of preferential binding to the forward site due to the recovery stroke. The theory reproduces a variety of experimental data, including the power stroke and slow diffusive search regimes in the mean trajectory of the detached head, and the force dependence of the forward-to-backward step ratio, run length, and velocity. We derive a stall force formula, determined by lever arm compliance and chemical cycle rates. By exploring the MyoV design space, we predict that it is a robust motor whose dynamical behavior is not compromised by reasonable perturbations to the reaction cycle and changes in the architecture of the lever arm.

From mechanical folding trajectories to intrinsic energy landscapes of biopolymers.

In single-molecule laser optical tweezer (LOT) pulling experiments, a protein or RNA is juxtaposed between DNA handles that are attached to beads in optical traps. The LOT generates folding trajectories under force in terms of time-dependent changes in the distance between the beads. How to construct the full intrinsic folding landscape (without the handles and beads) from the measured time series is a major unsolved problem. By using rigorous theoretical methods--which account for fluctuations of the DNA handles, rotation of the optical beads, variations in applied tension due to finite trap stiffness, as well as environmental noise and limited bandwidth of the apparatus--we provide a tractable method to derive intrinsic free-energy profiles. We validate the method by showing that the exactly calculable intrinsic free-energy profile for a generalized Rouse model, which mimics the two-state behavior in nucleic acid hairpins, can be accurately extracted from simulated time series in a LOT setup regardless of the stiffness of the handles. We next apply the approach to trajectories from coarse-grained LOT molecular simulations of a coiled-coil protein based on the GCN4 leucine zipper and obtain a free-energy landscape that is in quantitative agreement with simulations performed without the beads and handles. Finally, we extract the intrinsic free-energy landscape from experimental LOT measurements for the leucine zipper.

Evidence of disorder in biological molecules from single molecule pulling experiments.

Heterogeneity in biological molecules, resulting in molecule-to-molecule variations in their dynamics and function, is an emerging theme. To elucidate the consequences of heterogeneous behavior at the single molecule level, we propose an exactly solvable model in which the unfolding rate due to mechanical force depends parametrically on an auxiliary variable representing an entropy barrier arising from fluctuations in internal dynamics. When the rate of fluctuations--a measure of dynamical disorder--is comparable to or smaller than the rate of force-induced unbinding, we show that there are two experimentally observable consequences: nonexponential survival probability at constant force, and a heavy-tailed rupture force distribution at constant loading rate. By fitting our analytical expressions to data from single molecule pulling experiments on proteins and DNA, we quantify the extent of disorder. We show that only by analyzing data over a wide range of forces and loading rates can the role of disorder due to internal dynamics be quantitatively assessed.