Understanding nucleic acid-ion interactions.

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid-ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid-ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid-ion interactions.

Ab initio generalized Langevin equation.

We introduce a machine learning-based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori-Zwanzig formalism, under the constraint of the fluctuation-dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

An encompassed representation of timescale hierarchies in first-order reaction network.

Complex networks are pervasive in various fields such as chemistry, biology, and sociology. In chemistry, first-order reaction networks are represented by a set of first-order differential equations, which can be constructed from the underlying energy landscape. However, as the number of nodes increases, it becomes more challenging to understand complex kinetics across different timescales. Hence, how to construct an interpretable, coarse-graining scheme that preserves the underlying timescales of overall reactions is of crucial importance. Here, we develop a scheme to capture the underlying hierarchical subsets of nodes, and a series of coarse-grained (reduced-dimensional) rate equations between the subsets as a function of time resolution from the original reaction network. Each of the coarse-grained representations guarantees to preserve the underlying slow characteristic timescales in the original network. The crux is the construction of a lumping scheme incorporating a similarity measure in deciphering the underlying timescale hierarchy, which does not rely on the assumption of equilibrium. As an illustrative example, we apply the scheme to four-state Markovian models and Claisen rearrangement of allyl vinyl ether (AVE), and demonstrate that the reduced-dimensional representation accurately reproduces not only the slowest but also the faster timescales of overall reactions although other reduction schemes based on equilibrium assumption well reproduce the slowest timescale but fail to reproduce the second-to-fourth slowest timescales with the same accuracy. Our scheme can be applied not only to the reaction networks but also to networks in other fields, which helps us encompass their hierarchical structures of the complex kinetics over timescales.

A theory of phonon-induced friction on molecular adsorbates.

In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface's atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).

Phonon screening and dissociation of excitons at finite temperatures from first principles.

The properties of excitons, or correlated electron-hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron-hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio [Formula: see text]-Bethe-Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and [Formula: see text]. We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron-hole pairs.

Topological wetting states of microdroplets on closed-loop structured surfaces: Breakdown of the Gibbs equation at the microscale.

Microdroplets are a class of soft matter that has been extensively employed for chemical, biochemical, and industrial applications. However, fabricating microdroplets with largely controllable contact-area shape and apparent contact angle, a key prerequisite for their applications, is still a challenge. Here, by engineering a type of surface with homocentric closed-loop microwalls/microchannels, we can achieve facile size, shape, and contact-angle tunability of microdroplets on the textured surfaces by design. More importantly, this class of surface topologies (with universal genus value = 1) allows us to reveal that the conventional Gibbs equation (widely used for assessing the edge effect on the apparent contact angle of macrodroplets) seems no longer applicable for water microdroplets or nanodroplets (evidenced by independent molecular dynamics simulations). Notably, for the flat surface with the intrinsic contact angle ~0°, we find that the critical contact angle on the microtextured counterparts (at edge angle 90°) can be as large as >130°, rather than 90° according to the Gibbs equation. Experiments show that the breakdown of the Gibbs equation occurs for microdroplets of different types of liquids including alcohol and hydrocarbon oils. Overall, the microtextured surface design and topological wetting states not only offer opportunities for diverse applications of microdroplets such as controllable chemical reactions and low-cost circuit fabrications but also provide testbeds for advancing the fundamental surface science of wetting beyond the Gibbs equation.

Capturing the mechanosensitivity of cell proliferation in models of epithelium.

Despite the primary role of cell proliferation in tissue development and homeostatic maintenance, the interplay between cell density, cell mechanoresponse, and cell growth and division is not yet understood. In this article, we address this issue by reporting on an experimental investigation of cell proliferation on all time- and length-scales of the development of a model tissue, grown on collagen-coated glass or deformable substrates. Through extensive data analysis, we demonstrate the relation between mechanoresponse and probability for cell division, as a function of the local cell density. Motivated by these results, we construct a minimal model of cell division in tissue environment that can recover the data. By parameterizing the growth and the dividing phases of the cell cycle, and introducing such a proliferation model in dissipative particle dynamics simulations, we recover the mechanoresponsive, time-dependent density profiles in 2D tissues growing to macroscopic scales. The importance of separating the cell population into growing and dividing cells, each characterized by a particular time scale, is further emphasized by calculations of density profiles based on adapted Fisher-Kolmogorov equations. Together, these results show that the mechanoresponse on the level of a constitutive cell and its proliferation results in a matrix-sensitive active pressure. The latter evokes massive cooperative displacement of cells in the invading tissue and is a key factor for developing large-scale structures in the steady state.

Scientific machine learning for modeling and simulating complex fluids.

The formulation of rheological constitutive equations-models that relate internal stresses and deformations in complex fluids-is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these "digital fluid twins" to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation-a task that is not achievable using any previously developed data-driven rheological equation of state.

Fast protein folding is governed by memory-dependent friction.

When described by a low-dimensional reaction coordinate, the folding rates of most proteins are determined by a subtle interplay between free-energy barriers, which separate folded and unfolded states, and friction. While it is commonplace to extract free-energy profiles from molecular trajectories, a direct evaluation of friction is far more elusive and typically relies on fits of measured reaction rates to memoryless reaction-rate theories. Here, using memory-kernel extraction methods founded on a generalized Langevin equation (GLE) formalism, we directly calculate the time-dependent friction acting on the fraction of native contacts reaction coordinate Q, evaluated for eight fast-folding proteins, taken from a published set of large-scale molecular dynamics protein simulations. Our results reveal that, across the diverse range of proteins represented in this dataset, friction is more influential than free-energy barriers in determining protein folding rates. We also show that proteins fold in a regime where the finite decay time of friction significantly reduces the folding times, in some instances by as much as a factor of 10, compared to predictions based on memoryless friction.

Theoretical isotherm equation for adsorption-induced structural transition on flexible metal-organic frameworks.

Flexible metal-organic frameworks (MOFs) exhibit an adsorption-induced structural transition known as "gate opening" or "breathing," resulting in an S-shaped adsorption isotherm. This unique feature of flexible MOFs offers significant advantages, such as a large working capacity, high selectivity, and intrinsic thermal management capability, positioning them as crucial candidates for revolutionizing adsorption separation processes. Therefore, the interest in the industrial applications of flexible MOFs is increasing, and the adsorption engineering for flexible MOFs is becoming important. However, despite the establishment of the theoretical background for adsorption-induced structural transitions, no theoretical equation is available to describe S-shaped adsorption isotherms of flexible MOFs. Researchers rely on various empirical equations for process simulations that can lead to unreliable outcomes or may overlook insights into improving material performance owing to parameters without physical meaning. In this study, we derive a theoretical equation based on statistical mechanics that could be a standard for the structural transition type adsorption isotherms, as the Langmuir equation represents type I isotherms. The versatility of the derived equation is shown through four examples of flexible MOFs that exhibit gate opening and breathing. The consistency of the formula with existing theories, including the osmotic free energy analysis and intrinsic thermal management capabilities, is also discussed. The developed theoretical equation may lead to more reliable and insightful outcomes in adsorption separation processes, further advancing the direction of industrial applications of flexible MOFs.

In-context operator learning with data prompts for differential equation problems.

This paper introduces the paradigm of "in-context operator learning" and the corresponding model "In-Context Operator Networks" to simultaneously learn operators from the prompted data and apply it to new questions during the inference stage, without any weight update. Existing methods are limited to using a neural network to approximate a specific equation solution or a specific operator, requiring retraining when switching to a new problem with different equations. By training a single neural network as an operator learner, rather than a solution/operator approximator, we can not only get rid of retraining (even fine-tuning) the neural network for new problems but also leverage the commonalities shared across operators so that only a few examples in the prompt are needed when learning a new operator. Our numerical results show the capability of a single neural network as a few-shot operator learner for a diversified type of differential equation problems, including forward and inverse problems of ordinary differential equations, partial differential equations, and mean-field control problems, and also show that it can generalize its learning capability to operators beyond the training distribution.

Correlated surface-charging behaviors of two electrodes in an electrochemical cell.

It is revealed herein that surface-charging behaviors of the two electrodes constituting an electrochemical cell cannot be described independently by their respective electric double-layer (EDL) properties. Instead, they are correlated in such a way that the surface-charging behavior of each electrode is determined by the EDL and the reaction kinetics at both electrodes. Two fundamental equations describing the correlated surface-charging behaviors are derived, and approximate analytical solutions are obtained at low and high current densities, respectively, to facilitate transparent understanding. Important implications of the presented conceptual analysis for theoretical and computational electrochemistry are discussed. A strategy of modulating the activity of one electrode by tuning EDL parameters of the other in a two-electrode electrochemical cell is demonstrated.

Gray and White Matter Metrics Demonstrate Distinct and Complementary Prediction of Differences in Cognitive Performance in Children: Findings from ABCD (N = 11,876).

Individual differences in cognitive performance in childhood are a key predictor of significant life outcomes such as educational attainment and mental health. Differences in cognitive ability are governed in part by variations in brain structure. However, studies commonly focus on either gray or white matter metrics in humans, leaving open the key question as to whether gray or white matter microstructure plays distinct or complementary roles supporting cognitive performance. To compare the role of gray and white matter in supporting cognitive performance, we used regularized structural equation models to predict cognitive performance with gray and white matter measures. Specifically, we compared how gray matter (volume, cortical thickness, and surface area) and white matter measures (volume, fractional anisotropy, and mean diffusivity) predicted individual differences in cognitive performance. The models were tested in 11,876 children (ABCD Study; 5,680 female, 6,196 male) at 10 years old. We found that gray and white matter metrics bring partly nonoverlapping information to predict cognitive performance. The models with only gray or white matter explained respectively 15.4 and 12.4% of the variance in cognitive performance, while the combined model explained 19.0%. Zooming in, we additionally found that different metrics within gray and white matter had different predictive power and that the tracts/regions that were most predictive of cognitive performance differed across metrics. These results show that studies focusing on a single metric in either gray or white matter to study the link between brain structure and cognitive performance are missing a key part of the equation.

Linear and Nonlinear Behaviors of the Photoreceptor Coupled Network.

Photoreceptors are electrically coupled to one another, and the spatiotemporal properties of electrical synapses in a two-dimensional retinal network are still not well studied, because of the limitation of the single electrode or pair recording techniques which do not allow simultaneously measuring responses of multiple photoreceptors at various locations in the retina. A multiple electrode recording system is needed. In this study, we investigate the network properties of the two-dimensional rod coupled array of the salamander retina (both sexes were used) by using the newly available multiple patch electrode system that allows simultaneous recordings from up to eight cells and to determine the electrical connectivity among multiple rods. We found direct evidence that voltage signal spread in the rod-rod coupling network in the absence of I h (mediated by HCN channels) is passive and follows the linear cable equation. Under physiological conditions, I h shapes the network signal by progressively shortening the response time-to-peak of distant rods, compensating the time loss of signal traveling from distant rods to bipolar cell somas and facilitating synchronization of rod output signals. Under voltage-clamp conditions, current flow within the coupled rods follows Ohm's law, supporting the idea that nonlinear behaviors of the rod network are dependent on membrane voltage. Rod-rod coupling is largely symmetrical in the 2D array, and voltage-clamp blocking the next neighboring rod largely suppresses rod signal spread into the second neighboring rod, suggesting that indirect coupling pathways play a minor role in rod-rod coupling.

Mendelian randomization analysis using multiple biomarkers of an underlying common exposure.

Mendelian randomization (MR) analysis is increasingly popular for testing the causal effect of exposures on disease outcomes using data from genome-wide association studies. In some settings, the underlying exposure, such as systematic inflammation, may not be directly observable, but measurements can be available on multiple biomarkers or other types of traits that are co-regulated by the exposure. We propose a method for MR analysis on latent exposures (MRLE), which tests the significance for, and the direction of, the effect of a latent exposure by leveraging information from multiple related traits. The method is developed by constructing a set of estimating functions based on the second-order moments of GWAS summary association statistics for the observable traits, under a structural equation model where genetic variants are assumed to have indirect effects through the latent exposure and potentially direct effects on the traits. Simulation studies show that MRLE has well-controlled type I error rates and enhanced power compared to single-trait MR tests under various types of pleiotropy. Applications of MRLE using genetic association statistics across five inflammatory biomarkers (CRP, IL-6, IL-8, TNF-α, and MCP-1) provide evidence for potential causal effects of inflammation on increasing the risk of coronary artery disease, colorectal cancer, and rheumatoid arthritis, while standard MR analysis for individual biomarkers fails to detect consistent evidence for such effects.

Selective sweeps in SARS-CoV-2 variant competition.

The main mathematical result in this paper is that change of variables in the ordinary differential equation (ODE) for the competition of two infections in a Susceptible-Infected-Removed (SIR) model shows that the fraction of cases due to the new variant satisfies the logistic differential equation, which models selective sweeps. Fitting the logistic to data from the Global Initiative on Sharing All Influenza Data (GISAID) shows that this correctly predicts the rapid turnover from one dominant variant to another. In addition, our fitting gives sensible estimates of the increase in infectivity. These arguments are applicable to any epidemic modeled by SIR equations.

Anomalous heating in a colloidal system.

We report anomalous heating in a colloidal system, an experimental observation of the inverse Mpemba effect, where for two initial temperatures lower than the temperature of the thermal bath, the colder of the two systems heats up faster when coupled to the same thermal bath. For an overdamped, Brownian colloidal particle moving in a tilted double-well potential, we find a nonmonotonic dependence of the heating times on the initial temperature of the system. Entropic effects make the inverse Mpemba effect generically weaker-harder to observe-than the usual Mpemba effect (anomalous cooling). We also observe a strong version of anomalous heating, where a cold system heats up exponentially faster than systems prepared under slightly different conditions.

A method to predict the response to directional selection using a Kalman filter.

Predicting evolution remains challenging. The field of quantitative genetics provides predictions for the response to directional selection through the breeder's equation, but these predictions can have errors. The sources of these errors include omission of traits under selection, inaccurate estimates of genetic variance, and nonlinearities in the relationship between genetic and phenotypic variation. Previous research showed that the expected value of these prediction errors is often not zero, so predictions are systematically biased. Here, we propose that this bias, rather than being a nuisance, can be used to improve the predictions. We use this to develop a method to predict evolution, which is built on three key innovations. First, the method predicts change as the breeder's equation plus a bias term. Second, the method combines information from the breeder's equation and from the record of past changes in the mean to predict change using a Kalman filter. Third, the parameters of the filter are fitted in each generation using a learning algorithm on the record of past changes. We compare the method to the breeder's equation in two artificial selection experiments, one using the wing of the fruit fly and another using simulations that include a complex mapping of genotypes to phenotypes. The proposed method outperforms the breeder's equation, particularly when traits under selection are omitted from the analysis, when data are noisy, and when additive genetic variance is estimated inaccurately or not estimated at all. The proposed method is easy to apply, requiring only the trait means over past generations.

Likelihood-based non-Markovian models from molecular dynamics.

SignificanceThe analysis of complex systems with many degrees of freedom generally involves the definition of low-dimensional collective variables more amenable to physical understanding. Their dynamics can be modeled by generalized Langevin equations, whose coefficients have to be estimated from simulations of the initial high-dimensional system. These equations feature a memory kernel describing the mutual influence of the low-dimensional variables and their environment. We introduce and implement an approach where the generalized Langevin equation is designed to maximize the statistical likelihood of the observed data. This provides an efficient way to generate reduced models to study dynamical properties of complex processes such as chemical reactions in solution, conformational changes in biomolecules, or phase transitions in condensed matter systems.

Using single-cell models to predict the functionality of synthetic circuits at the population scale.

SignificanceAt the single-cell level, biochemical processes are inherently stochastic. For many natural systems, the resulting cell-to-cell variability is exploited by microbial populations. In synthetic biology, however, the interplay of cell-to-cell variability and population processes such as selection or growth often leads to circuits not functioning as predicted by simple models. Here we show how multiscale stochastic kinetic models that simultaneously track single-cell and population processes can be obtained based on an augmentation of the chemical master equation. These models enable us to quantitatively predict complex population dynamics of a yeast optogenetic differentiation system from a specification of the circuit's components and to demonstrate how cell-to-cell variability can be exploited to purposefully create unintuitive circuit functionality.