Publications

Molecular electronics break-junction experiments are widely used to investigate fundamental physics and chemistry at the nanoscale. Reproducibility in these experiments relies on measuring conductance on thousands of freshly formed molecular junctions, yielding a broad histogram of conductance events. Experiments typically focus on the most probable conductance, while the information content of the conductance histogram has remained unclear. Here we develop a microscopic theory for the conductance histogram by merging the theory of force-spectroscopy with molecular conductance. The procedure yields analytical equations that accurately fit the conductance histogram of a wide range of molecular junctions and augments the information content that can be extracted from them. Our formulation captures contributions to the conductance dispersion due to conductance changes during the mechanical elongation inherent to the experiments. In turn, the histogram shape is determined by the non-equilibrium stochastic features of junction rupture and formation. The microscopic parameters in the theory capture the junction’s electromechanical properties and can be isolated from separate conductance and rupture force (or junction-lifetime) measurements. The predicted behavior can be used to test the range of validity of the theory, understand the conductance histograms, design molecular junction experiments with enhanced resolution and molecular devices with more reproducible conductance properties.

We study bacterial diffusion in disordered porous media. Interactions with obstacles, at unknown locations, make this problem challenging. We approach it by abstracting the environment to cell states with memoryless transitions. With this, we derive an effective diffusivity that agrees well with simulations in explicit geometries. The diffusivity is non-monotonic, and we solve the optimal run length. We also find a rescaling that causes all of the theory and simulations to collapse. Our results indicate that a small set of microscopic features captures bacterial diffusion in disordered media.

The propulsion of mammalian spermatozoa relies on the spontaneous periodic oscillation of their flagella. These oscillations are driven internally by the coordinated action of ATP-powered dynein motors that exert sliding forces between microtubule doublets, resulting in bending waves that propagate along the flagellum and enable locomotion. We present an integrated chemomechanical model of a freely swimming spermatozoon that uses a sliding-control model of the axoneme capturing the two-way feedback between motor kinetics and elastic deformations while accounting for detailed fluid mechanics around the moving cell. We develop a robust computational framework that solves a boundary integral equation for the passive sperm head alongside the slender-body equation for the deforming flagellum described as a geometrically nonlinear internally actuated Euler-Bernoulli beam, and captures full hydrodynamic interactions. Nonlinear simulations are shown to produce spontaneous oscillations with realistic beating patterns and trajectories, which we analyze as a function of sperm number and motor activity. Our results indicate that the swimming velocity does not vary monotonically with dynein activity, but instead displays two maxima corresponding to distinct modes of swimming, each characterized by qualitatively different wave forms and trajectories. Our model also provides an estimate for the efficiency of swimming, which peaks at low sperm number.

A growing list of diverse biological systems and their equally diverse functionalities provides realizations of a paradigm of emergent behavior. In each of these biological systems, pervasive ensembles of weak, short-lived, spatially local interactions act autonomously to convey functionalities at larger spatial and temporal scales. In this article, a range of diverse systems and functionalities are presented in a cursory manner with literature citations for further details. Then two systems and their properties are discussed in more detail: yeast chromosome biology and human respiratory mucus

Understanding the transport properties of microorganisms and self-propelled particles in porous media has important implications for human health as well as microbial ecology. In free space, most microswimmers perform diffusive random walks as a result of the interplay of self-propulsion and orientation decorrelation mechanisms such as run-and-tumble dynamics or rotational diffusion. In an unstructured porous medium, collisions with the microstructure result in a decrease in the effective spatial diffusivity of the particles from its free-space value. Here, we analyze this problem for a simple model system consisting of noninteracting point particles performing run-and-tumble dynamics through a two-dimensional disordered medium composed of a random distribution of circular obstacles, in the absence of Brownian diffusion or hydrodynamic interactions. The particles are assumed to collide with the obstacles as hard spheres and subsequently slide on the obstacle surface with no frictional resistance while maintaining their orientation, until they either escape or tumble. We show that the variations in the long-time diffusivity can be described by a universal dimensionless hindrance function f (φ,Pe) of the obstacle area fraction φ and Péclet number Pe, or ratio of the swimmer run length to the obstacle size. We analytically derive an asymptotic expression for the hindrance function valid for dilute media (Peφ 1), and its extension to denser media is obtained using stochastic simulations. As we explain, the model is also easily generalized to describe dispersion in three dimensions.

When a founder cell and its progeny divide with incomplete cytokinesis, a network forms in which each intercellular bridge corresponds to a past mitotic event. Such networks are required for gamete production in many animals, and different species have evolved diverse final network topologies. Although mechanisms regulating network assembly have been identified in particular organisms, we lack a quantitative framework to understand network assembly and inter-species variability. Motivated by cell networks responsible for oocyte production in invertebrates, where the final topology is typically invariant within each species, we devised a mathematical model for generating cell networks, in which each node is an oscillator and, after a full cycle, the node produces a daughter to which it remains connected. These cell cycle oscillations are transient and coupled via diffusion over the edges of the network. By variation of three biologically motivated parameters, our model generates nearly all such networks currently reported across invertebrates. Furthermore, small parameter variations can rationalize cases of intra-species variation. Because cell networks outside of the ovary often form less deterministically, we propose model generalizations to account for sources of stochasticity.

High wall shear stress (WSS) is associated with risk of atherosclerotic plaque rupture, but there are numerous gaps in validating ultrasound-derived measurements of the parameter. Two major challenges are using simple models of stenosis and only evaluating WSS along a single 2D plane. To overcome these limitations, a novel simulation framework is herein demonstrated. The framework first models volumetric blood flow in actual stenosed human carotid artery geometries (using an immersed interface method (IIM) fluid structure interaction solver) and calculates the associated WSS. Then, the framework projects the modeled blood flow onto scatterers in Field II simulations of its ultrasonic interrogation. Volumetric ultrasound vector Doppler (VD) imaging using an elevationally swept L7-4 linear array was simulated in Field II, with variations in transmit sequences and flow conditions. In a ~55% stenosed human carotid artery under 600 mL/min flow, Bland-Altman analysis showed that a 3-angle plane wave (PW) transmit scheme estimated WSS with 0.04±0.64 Pa error (bias ± 95% CI) relative to the IIM ground truth, whereas transmitting with 5 angles increased accuracy, but decreased precision to -0.01±1.07 Pa, due to aliasing. These findings illustrate that the simulation framework enables direct comparison of data acquisition and processing methods for efficient development, validation, and refinement of WSS estimation methods in realistic clinical environments.

Theranostic materials research is experiencing rapid growth driven by the interest in integrating both therapeutic and diagnostic modalities. These materials offer the unique capability to not only provide treatment but also track the progression of a disease. However, to create an ideal theranostic biomaterial without compromising drug encapsulation, diagnostic imaging must be optimized for improved sensitivity and spatial localization. Herein, we create a protein-engineered fluorinated coiled-coil fiber, Q2TFL, capable of improved sensitivity to 19F magnetic resonance spectroscopy (MRS) detection. Leveraging residue-specific noncanonical amino acid incorporation of trifluoroleucine (TFL) into the coiled-coil, Q2, which self-assembles into nanofibers, we generate Q2TFL. We demonstrate that fluorination results in a greater increase in thermostability and 19F magnetic resonance detection compared to the nonfluorinated parent, Q2. Q2TFL also exhibits linear ratiometric 19F MRS thermoresponsiveness, allowing it to act as a temperature probe. Furthermore, we explore the ability of Q2TFL to encapsulate the anti-inflammatory small molecule, curcumin (CCM), and its impact on the coiled-coil structure. Q2TFL also provides hyposignal contrast in 1H MRI, echogenic signal with high-frequency ultrasound and sensitive detection by 19F MRS in vivo illustrating fluorination of coiled-coils for supramolecular assembly and their use with 1H MRI, 19F MRS and high frequency ultrasound as multimodal theranostic agents.

Transient receptor potential (TRP) ion channels are among the most well-studied classes of temperature-sensing molecules. Yet, the molecular mechanism and thermodynamic basis for the temperature sensitivity of TRP channels remains to this day poorly understood. One hypothesis is that the temperature-sensing mechanism can simply be described by a difference in heat capacity between the closed and open channel states. While such a two-state model may be simplistic it nonetheless has descriptive value, in the sense that it can be used to compare overall temperature sensitivity between different channels and mutants. Here, we introduce a mathematical framework based on the two-state model to reliably extract temperature-dependent thermodynamic potentials and heat capacities from measurements of equilibrium constants at different temperatures. Our framework is implemented in an open-source data analysis package that provides a straightforward way to fit both linear and nonlinear van ’t Hoff plots, thus avoiding some of the previous, potentially erroneous, assumptions when extracting thermodynamic variables from TRP channel electrophysiology data.