Publications

Extracting consistent statistics between relevant free energy minima of a molecular system is essential for physics, chemistry, and biology. Molecular dynamics (MD) simulations can aid in this task but are computationally expensive, especially for systems that require quantum accuracy. To overcome this challenge, we developed an approach combining enhanced sampling with deep generative models and active learning of a machine learning potential (MLP). We introduce an adaptive Markov chain Monte Carlo framework that enables the training of one normalizing flow (NF) and one MLP per state, achieving rapid convergence toward the Boltzmann distribution. Leveraging the trained NF and MLP models, we compute thermodynamic observables such as free energy differences and optical spectra. We apply this method to study the isomerization of an ultrasmall silver nanocluster belonging to a set of systems with diverse applications in the fields of medicine and catalysis.

We consider the analysis of singular waveguides separating insulating phases in two-space dimensions. The insulating domains are modeled by a massive Schrödinger equation and the singular waveguide by appropriate jump conditions along the one-dimensional interface separating the insulators. We present an integral formulation of the problem and analyze its mathematical properties. We also implement a fast multipole and sweeping-accelerated iterative algorithm for solving the integral equations, and demonstrate numerically the fast convergence of this method. Several numerical examples of solutions and scattering effects illustrate our theory.

Pituitary adenomas cause hormonal dysregulation and severe morbidity. The absence of clear biomarkers for prognosis/treatment and the high risk of recurrence make management challenging. Non-functioning gonadotroph pituitary adenomas (NFPAs) represent proliferation of gonadotroph lineage cells without an increase in secretion of gonadotropins. Detection is often delayed until mass effects cause visual defects. Higher resolution molecular study may improve classification, diagnosis, treatment, and provide insight into adenoma biology. We analyzed four NFPAs using unbiased genome-wide same-cell single nucleus 10X Genomics multiome assay (transcriptome and chromatin accessibility) to define individual cell identity/states in comparison with normal gonadotrope cell references from twelve age/sex-matched human post-mortem pituitaries. In total, we generated high quality multiomics single-cell data for 33,832 adenoma cells and 5,451 normal gonadotropes. We identified tumor, immune, vascular, and proliferative cells confirming cellular heterogeneity. Multidimensional scaling analysis using Manhattan distances between cell types showed that NFPA gonadotropes were closer in gene expression to healthy gonadotropes among all pituitary cell types. In addition, we observed differences in the macrophage and endothelial cell populations between tumors and healthy pituitaries. Next, we studied gene regulatory circuitry (regulated gene, modulating transcription factor and its chromatin interaction site) using the new Control of Regulation Extracted from Multiomics Assays (CREMA, [1]) multiome analysis method which showed high variation between normal and tumor gonadotropes. Within each adenoma, the NFPA tumor cells encompassed two main subgroups of cells, one group (GT+ marker tumor cells) expressing established gonadotrope markers and another group not (GT- marker tumor cells). GT+ and GT- tumor cells formed a cluster distinct from healthy gonadotropes. GT+ and GT- tumor cells differed by their level of expression of mitochondrial genes, suggesting decreased aerobic capacity associated with the GT- tumor subgroup. Interestingly, differential analyses between the healthy and the NFPA gonadotropes from same-cell multiome data highlighted dissimilar gene expression and chromatin accessibility patterns.Our newly collected single cell data characterize the difference between NFPA tumor cells and normal gonadotrope cells. Our results suggest a Warburg-effect like difference in aerobic respiratory capacity among tumor cells within the same pituitary adenoma. Overall, this study brings new insights into the molecular characteristics of NFPAs, including tumor cell heterogeneity, and represents a new resource for the field.

Habituation—a phenomenon in which a dynamical system exhibits a diminishing response to repeated stimulations that eventually recovers when the stimulus is withheld—is universally observed in living systems from animals to unicellular organisms. Despite its prevalence, generic mechanisms for this fundamental form of learning remain poorly defined. Drawing inspiration from prior work on systems that respond adaptively to step inputs, we study habituation from a nonlinear dynamics perspective. This approach enables us to formalize classical hallmarks of habituation that have been experimentally identified in diverse organisms and stimulus scenarios. We use this framework to investigate distinct dynamical circuits capable of habituation. In particular, we show that driven linear dynamics of a memory variable with static nonlinearities acting at the input and output can implement numerous hallmarks in a mathematically interpretable manner. This work establishes a foundation for understanding the dynamical substrates of this primitive learning behavior and offers a blueprint for the identification of habituating circuits in biological systems.

Tensor networks are a compressed format for multi-dimensional data. One-dimensional tensor networks -- often referred to as tensor trains (TT) or matrix product states (MPS) -- are increasingly being used as a numerical ansatz for continuum functions by "quantizing" the inputs into discrete binary digits. Here we demonstrate the power of more general tree tensor networks for this purpose. We provide direct constructions of a number of elementary functions as generic tree tensor networks and interpolative constructions for more complicated functions via a generalization of the tensor cross interpolation algorithm. For a range of multi-dimensional functions we show how more structured tree tensor networks offer a significantly more efficient ansatz than the commonly used tensor train. We demonstrate an application of our methods to solving multi-dimensional, non-linear Fredholm equations, providing a rigorous bound on the rank of the solution which, in turn, guarantees exponentially scaling accuracy with the size of the tree tensor network for certain problems.

Image classification and denoising suffer from complementary issues of lack of robustness or partially ignoring conditioning information. We argue that they can be alleviated by unifying both tasks through a model of the joint probability of (noisy) images and class labels. Classification is performed with a forward pass followed by conditioning. Using the Tweedie-Miyasawa formula, we evaluate the denoising function with the score, which can be computed by marginalization and back-propagation. The training objective is then a combination of cross-entropy loss and denoising score matching loss integrated over noise levels. Numerical experiments on CIFAR-10 and ImageNet show competitive classification and denoising performance compared to reference deep convolutional classifiers/denoisers, and significantly improves efficiency compared to previous joint approaches. Our model shows an increased robustness to adversarial perturbations compared to a standard discriminative classifier, and allows for a novel interpretation of adversarial gradients as a difference of denoisers.

2D template matching (2DTM) can be used to detect molecules and their assemblies in cellular cryo-EM images with high positional and orientational accuracy. While 2DTM successfully detects spherical targets such as large ribosomal subunits, challenges remain in detecting smaller and more aspherical targets in various environments. In this work, a novel 2DTM metric, referred to as the 2DTM p-value, is developed to extend the 2DTM framework to more complex applications. The 2DTM p-value combines information from two previously used 2DTM metrics, namely the 2DTM signal-to-noise ratio (SNR) and z-score, which are derived from the cross-correlation coefficient between the target and the template. The 2DTM p-value demonstrates robust detection accuracies under various imaging and sample conditions and outperforms the 2DTM SNR and z-score alone. Specifically, the 2DTM p-value improves the detection of aspherical targets such as a modified artificial tubulin patch particle (500 kDa) and a much smaller clathrin monomer (193 kDa) in simulated data. It also accurately recovers mature 60S ribosomes in yeast lamellae samples, even under conditions of increased Gaussian noise. The new metric will enable the detection of a wider variety of targets in both purified and cellular samples through 2DTM.

2D template matching (2DTM) can be used to detect molecules and their assemblies in cellular cryo-EM images with high positional and orientational accuracy. While 2DTM successfully detects spherical targets such as large ribosomal subunits, challenges remain in detecting smaller and more aspherical targets in various environments. In this work, a novel 2DTM metric, referred to as the 2DTM p-value, is developed to extend the 2DTM framework to more complex applications. The 2DTM p-value combines information from two previously used 2DTM metrics, namely the 2DTM signal-to-noise ratio (SNR) and z-score, which are derived from the cross-correlation coefficient between the target and the template. The 2DTM p-value demonstrates robust detection accuracies under various imaging and sample conditions and outperforms the 2DTM SNR and z-score alone. Specifically, the 2DTM p-value improves the detection of aspherical targets such as a modified artificial tubulin patch particle (500 kDa) and a much smaller clathrin monomer (193 kDa) in simulated data. It also accurately recovers mature 60S ribosomes in yeast lamellae samples, even under conditions of increased Gaussian noise. The new metric will enable the detection of a wider variety of targets in both purified and cellular samples through 2DTM.

The terminal cells of the Drosophila larval tracheal system are perhaps the simplest delivery networks, providing an analogue for mammalian vascular growth and function in a system with many fewer components. These cells are a prime example of single-cell morphogenesis, branching significantly over time to adapt to the needs of the growing tissue they supply. While the genetic mechanisms governing local branching decisions have been studied extensively, an understanding of the emergence of a global network architecture is still lacking. Mapping out the full network architecture of populations of terminal cells at different developmental times of Drosophila larvae, we find that cell growth follows scaling laws relating the total edge length, supply area, and branch density. Using time-lapse imaging of individual terminal cells, we identify that the cells grow in three ways: by extending branches, by the side budding of new branches, and by internally growing existing branches. A generative model based on these modes of growth recapitulates statistical properties of the terminal cell network data. These results suggest that the scaling laws arise from the coupled contributions of branching and internal growth. This study establishes the terminal cell as a uniquely tractable model system for further studies of transportation and distribution networks.

In vivo molecular imaging tools hold immense potential to drive transformative breakthroughs by enabling researchers to visualize cellular and molecular interactions in real-time and/or at high resolution. These advancements will facilitate a deeper understanding of fundamental biological processes and their dysregulation in disease states. Here, we develop and characterize a self-assembling protein nanomicelle called collagen type I binding – thermoresponsive assembled protein (Col1-TRAP) that binds tightly to type I collagen in vitro with nanomolar affinity. For ex vivo visualization, Col1-TRAP is labeled with a near-infrared fluorescent dye (NIR-Col1-TRAP). Both Col1-TRAP and NIR-Col1-TRAP display approximately a 3.8-fold greater binding to type I collagen compared to TRAP when measured by surface plasmon resonance (SPR). We present a proof-of-concept study using NIR-Col1-TRAP to detect fibrotic type I collagen deposition ex vivo in the livers of mice with non-alcoholic steatohepatitis (NASH). We show that NIR-Col1-TRAP demonstrates significantly decreased plasma recirculation time as well as increased liver accumulation in the NASH mice compared to mice without disease over 4 hours. As a result, NIR-Col1-TRAP shows potential as an imaging probe for NASH with in vivo targeting performance after injection in mice.