Publications

The quantum Fourier transform (QFT) is a key component of many important quantum algorithms, most famously being the essential ingredient in Shor’s algorithm for factoring products of primes. Given its remarkable capability, one would think it can introduce large entanglement to qubit systems and would be difficult to simulate classically. While early results showed the QFT indeed has maximal operator entanglement, we show that this is entirely due to the bit reversal in the QFT. The core part of the QFT has Schmidt coefficients decaying exponentially quickly, and thus it can generate only a constant amount of entanglement regardless of the number of qubits. In addition, we show the entangling power of the QFT is the same as the time evolution of a Hamiltonian with exponentially decaying interactions, and thus a variant of the area law for dynamics can be used to understand the low entanglement intuitively. Using the low entanglement property of the QFT, we show that classical simulations of the QFT on a matrix product state with low bond dimension take time linear in the number of qubits, providing a potential speedup over the classical fast Fourier transform on many classes of functions. We demonstrate this speedup in test calculations on some simple functions. For data vectors of length 10^6 – 10^8, the speedup can be a few orders of magnitude.

To check the accuracy of Bayesian computations, it is common to use rank-based simulation-based calibration (SBC). However, SBC has drawbacks: The test statistic is somewhat ad-hoc, interactions are difficult to examine, multiple testing is a challenge, and the resulting p-value is not a divergence metric. We propose to replace the marginal rank test with a flexible classification approach that learns test statistics from data. This measure typically has a higher statistical power than the SBC rank test and returns an interpretable divergence measure of miscalibration, computed from classification accuracy. This approach can be used with different data generating processes to address likelihood-free inference or traditional inference methods like Markov chain Monte Carlo or variational inference. We illustrate an automated implementation using neural networks and statistically-inspired features, and validate the method with numerical and real data experiments.

Recent studies have highlighted the significance of the spindle midzone – the region positioned between chromosomes – in ensuring proper chromosome segregation. By combining advanced 3D electron tomography and cutting-edge light microscopy we have discovered a previously unknown role of the regulation of microtubule dynamics within the spindle midzone of C. elegans. Using Fluorescence recovery after photobleaching and a combination of second harmonic generation and two-photon fluorescence microscopy, we found that the length of the antiparallel microtubule overlap zone in the spindle midzone is constant throughout anaphase, and independent of cortical pulling forces as well as the presence of the microtubule bundling protein SPD-1. Further investigations of SPD-1 and the chromokinesin KLP-19 in C. elegans suggest that KLP-19 regulates the overlap length and functions independently of SPD-1. Our data shows that KLP-19 plays an active role in regulating the length and turn-over of microtubules within the midzone as well as the size of the antiparallel overlap region throughout mitosis. Depletion of KLP-19 in mitosis leads to an increase in microtubule length in the spindle midzone, which also leads to increased microtubule – microtubule interaction, thus building up a more robust microtubule network. The spindle is globally stiffer and more stable, which has implications for the transmission of forces within the spindle affecting chromosome segregation dynamics. Our data shows that by localizing KLP-19 to the spindle midzone in anaphase microtubule dynamics can be locally controlled allowing the formation of a functional midzone.

The design of protein interaction inhibitors is a promising approach to address aberrant protein interactions that cause disease. One strategy in designing inhibitors is to use peptidomimetic scaffolds that mimic the natural interaction interface. A central challenge in using peptidomimetics as protein interaction inhibitors, however, is determining how best the molecular scaffold aligns to the residues of the interface it is attempting to mimic. Here we present the Scaffold Matcher algorithm that aligns a given molecular scaffold onto hotspot residues from a protein interaction interface. To optimize the degrees of freedom of the molecular scaffold we implement the covariance matrix adaptation evolution strategy (CMA-ES), a state-of-the-art derivative-free optimization algorithm in Rosetta. To evaluate the performance of the CMA-ES, we used 26 peptides from the FlexPepDock Benchmark and compared with three other algorithms in Rosetta, specifically, Rosetta's default minimizer, a Monte Carlo protocol of small backbone perturbations, and a Genetic algorithm. We test the algorithms' performance on their ability to align a molecular scaffold to a series of hotspot residues (i.e., constraints) along native peptides. Of the 4 methods, CMA-ES was able to find the lowest energy conformation for all 26 benchmark peptides. Additionally, as a proof of concept, we apply the Scaffold Match algorithm with CMA-ES to align a peptidomimetic oligooxopiperazine scaffold to the hotspot residues of the substrate of the main protease of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our implementation of CMA-ES into Rosetta allows for an alternative optimization method to be used on macromolecular modeling problems with rough energy landscapes. Finally, our Scaffold Matcher algorithm allows for the identification of initial conformations of interaction inhibitors that can be further designed and optimized as high-affinity reagents.

Vertebrate groups have evolved strikingly diverse color patterns. However, it remains unknown to what extent the diversification of such patterns has been shaped by the proximate, developmental mechanisms that regulate their formation. While these developmental mechanisms have long been inaccessible empirically, here we take advantage of recent insights into rodent pattern formation to investigate the role of development in shaping pattern diversification across rodents. Based on a broad survey of museum specimens, we first establish that various rodents have independently evolved diverse patterns consisting of longitudinal stripes, varying across species in number, color, and relative positioning. We then interrogate this diversity using a simple model that incorporates recent molecular and developmental insights into stripe formation in African striped mice. Our results suggest that, on the one hand, development has facilitated pattern diversification: The diversity of patterns seen across species can be generated by a single developmental process, and small changes in this process suffice to recapitulate observed evolutionary changes in pattern organization. On the other hand, development has constrained diversification: Constraints on stripe positioning limit the scope of evolvable patterns, and although pattern organization appears at first glance phylogenetically unconstrained, development turns out to impose a cryptic constraint. Altogether, this work reveals that pattern diversification in rodents can in part be explained by the underlying development and illustrates how pattern formation models can be leveraged to interpret pattern evolution.

We present a high-order boundary integral equation (BIE) method for the frequency-domain acoustic scattering of a point source by a singly-periodic, infinite, corrugated boundary. We apply it to the accurate numerical study of acoustic radiation in the neighborhood of a sound-hard two-dimensional staircase modeled after the El Castillo pyramid. Such staircases support trapped waves which travel along the surface and decay exponentially away from it. We use the array scanning method (Floquet--Bloch transform) to recover the scattered field as an integral over the family of quasiperiodic solutions parameterized by their on-surface wavenumber. Each such BIE solution requires the quasiperiodic Green's function, which we evaluate using an efficient integral representation of lattice sum coefficients. We avoid the singularities and branch cuts present in the array scanning integral by complex contour deformation. For each frequency, this enables a solution accurate to around 10 digits in a couple of seconds. We propose a residue method to extract the limiting powers carried by trapped modes far from the source. Finally, by computing the trapped mode dispersion relation, we use a simple ray model to explain an observed acoustic "raindrop" effect (chirp-like time-domain response).

During morphogenesis, epithelial monolayers actively alter their shape to create future body parts of the animal; this makes the epithelium one of the most active and critical structures in early animal development. While epithelia are often modeled as two-dimensional systems, real epithelia are not necessarily thin relative to cell cross section, and advances in 3D imaging have shown the possibility of substantial cell deformations in the third dimension, as well as differences in dynamics of the apical and basal surfaces indicative of three-dimensional coupling. With the importance of the third dimension in mind, we have developed a self-sculpting, three-dimensional model of epithelia whose dynamics are driven by active forces on its surface. We present a first, fundamental study for a reduced version of epithelia that investigates how surface forces affect its internal dynamics. Our model captures the 3D slab-like geometry of epithelia, viscoelasticity of tissue response, fluid surroundings, and driving from active surface forces. We represent epithelial tissue as a thick slab, a 3D continuum comprised of a Stokes fluid with an extra viscoelastic stress. Employing this model, we present both analytical and numerical solutions of the system and make quantitative predictions about cell shapes, cell dynamics, and the tissue's response to surface force in a three-dimensional setting. In particular, we investigate the implications of our model on the initiation of ventral furrow invagination and T1 transitions in Drosophila embryogenesis. In the former, we demonstrate the importance of fluid and geometric surroundings to drive invagination. In the latter, we show the limitations of surface forces alone to drive T1 transitions.

In this paper, we present a novel methodology for automatic adaptive weighting of Bayesian Physics-Informed Neural Networks (BPINNs), and we demonstrate that this makes it possible to robustly address multi-objective and multiscale problems. BPINNs are a popular framework for data assimilation, combining the constraints of Uncertainty Quantification (UQ) and Partial Differential Equation (PDE). The relative weights of the BPINN target distribution terms are directly related to the inherent uncertainty in the respective learning tasks. Yet, they are usually manually set a-priori, that can lead to pathological behavior, stability concerns, and to conflicts between tasks which are obstacles that have deterred the use of BPINNs for inverse problems with multiscale dynamics.

The present weighting strategy automatically tunes the weights by considering the multitask nature of target posterior distribution. We show that this remedies the failure modes of BPINNs and provides efficient exploration of the optimal Pareto front. This leads to better convergence and stability of BPINN training while reducing sampling bias. The determined weights moreover carry information about task uncertainties, reflecting noise levels in the data and adequacy of the PDE model.

We demonstrate this in numerical experiments in Sobolev training, and compare them to analytically ε-optimal baseline, and in a multiscale Lotka-Volterra inverse problem. We eventually apply this framework to an inpainting task and an inverse problem, involving latent field recovery for incompressible flow in complex geometries.

Interactions between commuting individuals can lead to large-scale spreading of rumors, ideas, or disease, even though the commuters have no net displacement. The emergent dynamics depend crucially on the commuting distribution of a population, that is how the probability to travel to a destination decays with distance from home. Applying this idea to epidemics, we will demonstrate the qualitatively different infection dynamics emerging from populations with different commuting distributions. If the commuting distribution is exponentially localized, then we recover a reaction-diffusion system and observe Fisher waves traveling at a speed proportional to the characteristic commuting distance. If the commuting distribution has a long tail, then no finite-velocity waves can form, but we show that, in some regimes, there is nontrivial spatial dependence that the well-mixed approximation neglects. We discuss how, in all cases, an initial dispersal-dominated regime can allow the disease to go undetected for a finite amount of time before exponential growth takes over. This “offset time” is a quantity of huge importance for epidemic surveillance and yet largely ignored in the literature.

Endurance exercise is an important health modifier. We studied cell-type specific adaptations of human skeletal muscle to acute endurance exercise using single-nucleus (sn) multiome sequencing in human vastus lateralis samples collected before and 3.5 hours after 40 min exercise at 70% VO2max in four subjects, as well as in matched time of day samples from two supine resting circadian controls. High quality same-cell RNA-seq and ATAC-seq data were obtained from 37,154 nuclei comprising 14 cell types. Among muscle fiber types, both shared and fiber-type specific regulatory programs were identified. Single-cell circuit analysis identified distinct adaptations in fast, slow and intermediate fibers as well as LUM-expressing FAP cells, involving a total of 328 transcription factors (TFs) acting at altered accessibility sites regulating 2,025 genes. These data and circuit mapping provide single-cell insight into the processes underlying tissue and metabolic remodeling responses to exercise.