Search

A hallmark of Alzheimer's disease (AD) is the extracellular aggregation of toxic amyloid-beta (Aβ) peptides in form of plaques. Here, we identify netoglitazone, an antidiabetic compound previously tested in humans, as an Aβ aggregation antagonist. Netoglitazone improved cognition and reduced microglia activity in a mouse model of AD. Using quantitative whole-brain three-dimensional histology (Q3D), we precisely identified brain regions where netoglitazone reduced the number and size of Aβ plaques. We demonstrate the utility of Q3D in preclinical drug evaluation for AD by providing a high-resolution brain-wide view of drug efficacy. Applying Q3D has the potential to improve pre-clinical drug evaluation by providing information that can help identify mechanisms leading to brain region-specific drug efficacy.

Delivering macromolecules across biological barriers such as the blood-brain barrier limits their application in vivo. Previous work has demonstrated that Toxoplasma gondii, a parasite that naturally travels from the human gut to the central nervous system (CNS), can deliver proteins to host cells. Here we engineered T. gondii's endogenous secretion systems, the rhoptries and dense granules, to deliver multiple large (>100 kDa) therapeutic proteins into neurons via translational fusions to toxofilin and GRA16. We demonstrate delivery in cultured cells, brain organoids and in vivo, and probe protein activity using imaging, pull-down assays, scRNA-seq and fluorescent reporters. We demonstrate robust delivery after intraperitoneal administration in mice and characterize 3D distribution throughout the brain. As proof of concept, we demonstrate GRA16-mediated brain delivery of the MeCP2 protein, a putative therapeutic target for Rett syndrome. By characterizing the potential and current limitations of the system, we aim to guide future improvements that will be required for broader application., (© 2024. The Author(s).)

In recent years, advances in computing hardware and computational methods have prompted a wealth of activities for solving inverse problems in physics. These problems are often described by systems of partial differential equations (PDEs). The advent of machine learning has reinvigorated the interest in solving inverse problems using neural networks (NNs). In these efforts, the solution of the PDEs is expressed as NNs trained through the minimization of a loss function involving the PDE. Here, we show how to accelerate this approach by five orders of magnitude by deploying, instead of NNs, conventional PDE approximations. The framework of optimizing a discrete loss (ODIL) minimizes a cost function for discrete approximations of the PDEs using gradient-based and Newton's methods. The framework relies on grid-based discretizations of PDEs and inherits their accuracy, convergence, and conservation properties. The implementation of the method is facilitated by adopting machine-learning tools for automatic differentiation. We also propose a multigrid technique to accelerate the convergence of gradient-based optimizers. We present applications to PDE-constrained optimization, optical flow, system identification, and data assimilation. We compare ODIL with the popular method of physics-informed neural networks and show that it outperforms it by several orders of magnitude in computational speed while having better accuracy and convergence rates. We evaluate ODIL on inverse problems involving linear and nonlinear PDEs including the Navier-Stokes equations for flow reconstruction problems. ODIL bridges numerical methods and machine learning and presents a powerful tool for solving challenging, inverse problems across scientific domains., (© The Author(s) 2024. Published by Oxford University Press on behalf of National Academy of Sciences.)

We present a potent computational method for the solution of inverse problems in fluid mechanics. We consider inverse problems formulated in terms of a deterministic loss function that can accommodate data and regularization terms. We introduce a multigrid decomposition technique that accelerates the convergence of gradient-based methods for optimization problems with parameters on a grid. We incorporate this multigrid technique to the Optimizing a DIscrete Loss (ODIL) framework. The multiresolution ODIL (mODIL) accelerates by an order of magnitude the original formalism and improves the avoidance of local minima. Moreover, mODIL accommodates the use of automatic differentiation for calculating the gradients of the loss function, thus facilitating the implementation of the framework. We demonstrate the capabilities of mODIL on a variety of inverse and flow reconstruction problems: solution reconstruction for the Burgers equation, inferring conductivity from temperature measurements, and inferring the body shape from wake velocity measurements in three dimensions. We also provide a comparative study with the related, popular Physics-Informed Neural Networks (PINNs) method. We demonstrate that mODIL has three to five orders of magnitude lower computational cost than PINNs in benchmark problems including simple PDEs and lid-driven cavity problems. Our results suggest that mODIL is a very potent, fast and consistent method for solving inverse problems in fluid mechanics., (© 2023. The Author(s), under exclusive licence to EDP Sciences, SIF and Springer-Verlag GmbH Germany, part of Springer Nature.)

The stress-free state (SFS) of red blood cells (RBCs) is a fundamental reference configuration for the calibration of computational models, yet it remains unknown. Current experimental methods cannot measure the SFS of cells without affecting their mechanical properties, whereas computational postulates are the subject of controversial discussions. Here, we introduce data-driven estimates of the SFS shape and the visco-elastic properties of RBCs. We employ data from single-cell experiments that include measurements of the equilibrium shape of stretched cells and relaxation times of initially stretched RBCs. A hierarchical Bayesian model accounts for these experimental and data heterogeneities. We quantify, for the first time, the SFS of RBCs and use it to introduce a transferable RBC (t-RBC) model. The effectiveness of the proposed model is shown on predictions of unseen experimental conditions during the inference, including the critical stress of transitions between tumbling and tank-treading cells in shear flow. Our findings demonstrate that the proposed t-RBC model provides predictions of blood flows with unprecedented accuracy and quantified uncertainties., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The predictive capabilities of turbulent flow simulations, critical for aerodynamic design and weather prediction, hinge on the choice of turbulence models. The abundance of data from experiments and simulations and the advent of machine learning have provided a boost to turbulence modeling efforts. However, simulations of turbulent flows remain hindered by the inability of heuristics and supervised learning to model the near-wall dynamics. We address this challenge by introducing scientific multi-agent reinforcement learning (SciMARL) for the discovery of wall models for large-eddy simulations (LES). In SciMARL, discretization points act also as cooperating agents that learn to supply the LES closure model. The agents self-learn using limited data and generalize to extreme Reynolds numbers and previously unseen geometries. The present simulations reduce by several orders of magnitude the computational cost over fully-resolved simulations while reproducing key flow quantities. We believe that SciMARL creates unprecedented capabilities for the simulation of turbulent flows., (© 2022. The Author(s).)

Crashing ocean waves, cappuccino froths, and microfluidic bubble crystals are examples of foamy flows. Foamy flows are critical in numerous natural and industrial processes and remain notoriously difficult to compute as they involve coupled, multiscale physical processes. Computations need to resolve the interactions of the bubbles separated by stable thin liquid films. We present the multilayer volume-of-fluid method (Multi-VOF) that advances the state of the art in simulation capabilities of foamy flows. The method introduces a scheme to handle multiple bubbles that do not coalesce. Multi-VOF is verified and validated with experimental results and complemented with open-source software. We demonstrate capturing of crystalline structures of bubbles in realistic microfluidics devices and foamy flows involving tens of thousands of bubbles in a waterfall. The present technique extends the classical volume-of-fluid methodology and allows for large-scale predictive simulations of flows with multiple interfaces.

The extreme liquid transport properties of carbon nanotubes present new opportunities for surpassing conventional technologies in water filtration and purification. We demonstrate that carbon nanotubes with wettability surface patterns act as nanopumps for the ultrafast transport of picoliter water droplets without requiring externally imposed pressure gradients. Large-scale molecular dynamics simulations evidence unprecedented speeds and accelerations on the order of 10 10 g of droplet propulsion caused by interfacial energy gradients. This phenomenon is persistent for nanotubes of varying sizes, stepwise pattern configurations, and initial conditions. We present a scaling law for water transport as a function of wettability gradients through simple models for the droplet dynamic contact angle and friction coefficient. Our results show that patterned nanotubes are energy-efficient nanopumps offering a realistic path toward ultrafast water nanofiltration and precision drug delivery.

Simulations are vital for understanding and predicting the evolution of complex molecular systems. However, despite advances in algorithms and special purpose hardware, accessing the time scales necessary to capture the structural evolution of biomolecules remains a daunting task. In this work, we present a novel framework to advance simulation time scales by up to 3 orders of magnitude by learning the effective dynamics (LED) of molecular systems. LED augments the equation-free methodology by employing a probabilistic mapping between coarse and fine scales using mixture density network (MDN) autoencoders and evolves the non-Markovian latent dynamics using long short-term memory MDNs. We demonstrate the effectiveness of LED in the Müller-Brown potential, the Trp cage protein, and the alanine dipeptide. LED identifies explainable reduced-order representations, i.e., collective variables, and can generate, at any instant, all-atom molecular trajectories consistent with the collective variables. We believe that the proposed framework provides a dramatic increase to simulation capabilities and opens new horizons for the effective modeling of complex molecular systems.

The tunable polarity of water can be exploited in emerging technologies including catalysis, gas storage, and green chemistry. Recent experimental and theoretical studies have shown that water can be rendered into an effectively apolar solvent under nanoconfinement. We furthermore demonstrate, through molecular simulations, that the static dielectric constant of water can be modified by changing the wettability of the confining material. We find the out-of-plane dielectric response to be highly sensitive to the level of confinement and can be reduced up to 40× , in accordance with experimental data. By altering the surface wettability from superhydrophilic to superhydrophobic, we observe a 36% increase for the out-of-plane and a 31% decrease for the in-plane dielectric constants. Our findings demonstrate the feasibility of tunable water polarity, a phenomenon with great potential for scientific and technological impact.

Efficient point-to-point navigation in the presence of a background flow field is important for robotic applications such as ocean surveying. In such applications, robots may only have knowledge of their immediate surroundings or be faced with time-varying currents, which limits the use of optimal control techniques. Here, we apply a recently introduced Reinforcement Learning algorithm to discover time-efficient navigation policies to steer a fixed-speed swimmer through unsteady two-dimensional flow fields. The algorithm entails inputting environmental cues into a deep neural network that determines the swimmer's actions, and deploying Remember and Forget Experience Replay. We find that the resulting swimmers successfully exploit the background flow to reach the target, but that this success depends on the sensed environmental cue. Surprisingly, a velocity sensing approach significantly outperformed a bio-mimetic vorticity sensing approach, and achieved a near 100% success rate in reaching the target locations while approaching the time-efficiency of optimal navigation trajectories., (© 2021. The Author(s).)

We present high-fidelity numerical simulations of expiratory biosol transport during normal breathing under indoor, stagnant air conditions with and without a facile mask. We investigate mask efficacy to suppress the spread of saliva particles that is underpinnings existing social distancing recommendations. The present simulations incorporate the effect of human anatomy and consider a spectrum of saliva particulate sizes that range from 0.1 to 10  μ m while also accounting for their evaporation. The simulations elucidate the vorticity dynamics of human breathing and show that without a facile mask, saliva particulates could travel over 2.2 m away from the person. However, a non-medical grade face mask can drastically reduce saliva particulate propagation to 0.72 m away from the person. This study provides new quantitative evidence that facile masks can successfully suppress the spreading of saliva particulates due to normal breathing in indoor environments., (© 2021 Author(s).)

Effective intervention strategies for epidemics rely on the identification of their origin and on the robustness of the predictions made by network disease models. We introduce a Bayesian uncertainty quantification framework to infer model parameters for a disease spreading on a network of communities from limited, noisy observations; the state-of-the-art computational framework compensates for the model complexity by exploiting massively parallel computing architectures. Using noisy, synthetic data, we show the potential of the approach to perform robust model fitting and additionally demonstrate that we can effectively identify the disease origin via Bayesian model selection. As disease-related data are increasingly available, the proposed framework has broad practical relevance for the prediction and management of epidemics., Competing Interests: We declare we have no competing interests., (© 2021 The Authors.)

The systematic identification of infected individuals is critical for the containment of the COVID-19 pandemic. Currently, the spread of the disease is mostly quantified by the reported numbers of infections, hospitalisations, recoveries and deaths; these quantities inform epidemiology models that provide forecasts for the spread of the epidemic and guide policy making. The veracity of these forecasts depends on the discrepancy between the numbers of reported, and unreported yet infectious, individuals. We combine Bayesian experimental design with an epidemiology model and propose a methodology for the optimal allocation of limited testing resources in space and time, which maximises the information gain for such unreported infections. The proposed approach is applicable at the onset and spread of the epidemic and can forewarn of a possible recurrence of the disease after relaxation of interventions. We examine its application in Switzerland; the open source software is, however, readily adaptable to countries around the world. We find that following the proposed methodology can lead to vastly less uncertain predictions for the spread of the disease, thus improving estimates of the effective reproduction number and the future number of unreported infections. This information can provide timely and systematic guidance for the effective identification of infectious individuals and for decision-making regarding lockdown measures and the distribution of vaccines.

The reproduction number is broadly considered as a key indicator for the spreading of the COVID-19 pandemic. Its estimated value is a measure of the necessity and, eventually, effectiveness of interventions imposed in various countries. Here we present an online tool for the data-driven inference and quantification of uncertainties for the reproduction number, as well as the time points of interventions for 51 European countries. The study relied on the Bayesian calibration of the SIR model with data from reported daily infections from these countries. The model fitted the data, for most countries, without individual tuning of parameters. We also compared the results of SIR and SEIR models, which give different estimates of the reproduction number, and provided an analytical relationship between the respective numbers. We deployed a Bayesian inference framework with efficient sampling algorithms, to present a publicly available graphical user interface (https://cse-lab.ethz.ch/coronavirus) that allows the user to assess and compare predictions for pairs of European countries. The results quantified the rate of the disease’s spread before and after interventions, and provided a metric for the effectiveness of non-pharmaceutical interventions in different countries. They also indicated how geographic proximity and the times of interventions affected the progression of the epidemic.

Fish schooling implies an awareness of the swimmers for their companions. In flow mediated environments, in addition to visual cues, pressure and shear sensors on the fish body are critical for providing quantitative information that assists the quantification of proximity to other fish. Here we examine the distribution of sensors on the surface of an artificial swimmer so that it can optimally identify a leading group of swimmers. We employ Bayesian experimental design coupled with numerical simulations of the two-dimensional Navier Stokes equations for multiple self-propelled swimmers. The follower tracks the school using information from its own surface pressure and shear stress. We demonstrate that the optimal sensor distribution of the follower is qualitatively similar to the distribution of neuromasts on fish. Our results show that it is possible to identify accurately the center of mass and the number of the leading swimmers using surface only information.

Patient-specific modelling of haemodynamics in arterial networks has so far relied on parameter estimation for inexpensive or small-scale models. We describe here a Bayesian uncertainty quantification framework which makes two major advances: an efficient parallel implementation, allowing parameter estimation for more complex forward models, and a system for practical model selection, allowing evidence-based comparison between distinct physical models. We demonstrate the proposed methodology by generating simulated noisy flow velocity data from a branching arterial tree model in which a structural defect is introduced at an unknown location; our approach is shown to accurately locate the abnormality and estimate its physical properties even in the presence of significant observational and systemic error. As the method readily admits real data, it shows great potential in patient-specific parameter fitting for haemodynamical flow models., Competing Interests: The authors declare no competing interests., (© 2019 The Authors.)

Glioblastoma (GBM) is a highly invasive brain tumor, whose cells infiltrate surrounding normal brain tissue beyond the lesion outlines visible in the current medical scans. These infiltrative cells are treated mainly by radiotherapy. Existing radiotherapy plans for brain tumors derive from population studies and scarcely account for patient-specific conditions. Here, we provide a Bayesian machine learning framework for the rational design of improved, personalized radiotherapy plans using mathematical modeling and patient multimodal medical scans. Our method, for the first time, integrates complementary information from high-resolution MRI scans and highly specific FET-PET metabolic maps to infer tumor cell density in GBM patients. The Bayesian framework quantifies imaging and modeling uncertainties and predicts patient-specific tumor cell density with credible intervals. The proposed methodology relies only on data acquired at a single time point and, thus, is applicable to standard clinical settings. An initial clinical population study shows that the radiotherapy plans generated from the inferred tumor cell infiltration maps spare more healthy tissue thereby reducing radiation toxicity while yielding comparable accuracy with standard radiotherapy protocols. Moreover, the inferred regions of high tumor cell densities coincide with the tumor radioresistant areas, providing guidance for personalized dose-escalation. The proposed integration of multimodal scans and mathematical modeling provides a robust, non-invasive tool to assist personalized radiotherapy design.

We propose the S-leaping algorithm for the acceleration of Gillespie's stochastic simulation algorithm that combines the advantages of the two main accelerated methods; the [Formula: see text]-leaping and R-leaping algorithms. These algorithms are known to be efficient under different conditions; the [Formula: see text]-leaping is efficient for non-stiff systems or systems with partial equilibrium, while the R-leaping performs better in stiff system thanks to an efficient sampling procedure. However, even a small change in a system's set up can critically affect the nature of the simulated system and thus reduce the efficiency of an accelerated algorithm. The proposed algorithm combines the efficient time step selection from the [Formula: see text]-leaping with the effective sampling procedure from the R-leaping algorithm. The S-leaping is shown to maintain its efficiency under different conditions and in the case of large and stiff systems or systems with fast dynamics, the S-leaping outperforms both methods. We demonstrate the performance and the accuracy of the S-leaping in comparison with the [Formula: see text]-leaping and R-leaping on a number of benchmark systems involving biological reaction networks.

The necessity for accurate and computationally efficient representations of water in atomistic simulations that can span biologically relevant timescales has born the necessity of coarse-grained (CG) modeling. Despite numerous advances, CG water models rely mostly on a-priori specified assumptions. How these assumptions affect the model accuracy, efficiency, and in particular transferability, has not been systematically investigated. Here we propose a data driven comparison and selection for CG water models through a Hierarchical Bayesian framework. We examine CG water models that differ in their level of coarse-graining, structure, and number of interaction sites. We find that the importance of electrostatic interactions for the physical system under consideration is a dominant criterion for the model selection. Multi-site models are favored, unless the effects of water in electrostatic screening are not relevant, in which case the single site model is preferred due to its computational savings. The charge distribution is found to play an important role in the multi-site model's accuracy while the flexibility of the bonds/angles may only slightly improve the models. Furthermore, we find significant variations in the computational cost of these models. We present a data informed rationale for the selection of CG water models and provide guidance for future water model designs.

Fish in schooling formations navigate complex flow fields replete with mechanical energy in the vortex wakes of their companions. Their schooling behavior has been associated with evolutionary advantages including energy savings, yet the underlying physical mechanisms remain unknown. We show that fish can improve their sustained propulsive efficiency by placing themselves in appropriate locations in the wake of other swimmers and intercepting judiciously their shed vortices. This swimming strategy leads to collective energy savings and is revealed through a combination of high-fidelity flow simulations with a deep reinforcement learning (RL) algorithm. The RL algorithm relies on a policy defined by deep, recurrent neural nets, with long-short-term memory cells, that are essential for capturing the unsteadiness of the two-way interactions between the fish and the vortical flow field. Surprisingly, we find that swimming in-line with a leader is not associated with energetic benefits for the follower. Instead, "smart swimmer(s)" place themselves at off-center positions, with respect to the axis of the leader(s) and deform their body to synchronize with the momentum of the oncoming vortices, thus enhancing their swimming efficiency at no cost to the leader(s). The results confirm that fish may harvest energy deposited in vortices and support the conjecture that swimming in formation is energetically advantageous. Moreover, this study demonstrates that deep RL can produce navigation algorithms for complex unsteady and vortical flow fields, with promising implications for energy savings in autonomous robotic swarms., Competing Interests: The authors declare no conflict of interest.

The prediction of extreme events, from avalanches and droughts to tsunamis and epidemics, depends on the formulation and analysis of relevant, complex dynamical systems. Such dynamical systems are characterized by high intrinsic dimensionality with extreme events having the form of rare transitions that are several standard deviations away from the mean. Such systems are not amenable to classical order-reduction methods through projection of the governing equations due to the large intrinsic dimensionality of the underlying attractor as well as the complexity of the transient events. Alternatively, data-driven techniques aim to quantify the dynamics of specific, critical modes by utilizing data-streams and by expanding the dimensionality of the reduced-order model using delayed coordinates. In turn, these methods have major limitations in regions of the phase space with sparse data, which is the case for extreme events. In this work, we develop a novel hybrid framework that complements an imperfect reduced order model, with data-streams that are integrated though a recurrent neural network (RNN) architecture. The reduced order model has the form of projected equations into a low-dimensional subspace that still contains important dynamical information about the system and it is expanded by a long short-term memory (LSTM) regularization. The LSTM-RNN is trained by analyzing the mismatch between the imperfect model and the data-streams, projected to the reduced-order space. The data-driven model assists the imperfect model in regions where data is available, while for locations where data is sparse the imperfect model still provides a baseline for the prediction of the system state. We assess the developed framework on two challenging prototype systems exhibiting extreme events. We show that the blended approach has improved performance compared with methods that use either data streams or the imperfect model alone. Notably the improvement is more significant in regions associated with extreme events, where data is sparse., Competing Interests: The authors have declared that no competing interests exist.

We introduce a data-driven forecasting method for high-dimensional chaotic systems using long short-term memory (LSTM) recurrent neural networks. The proposed LSTM neural networks perform inference of high-dimensional dynamical systems in their reduced order space and are shown to be an effective set of nonlinear approximators of their attractor. We demonstrate the forecasting performance of the LSTM and compare it with Gaussian processes (GPs) in time series obtained from the Lorenz 96 system, the Kuramoto-Sivashinsky equation and a prototype climate model. The LSTM networks outperform the GPs in short-term forecasting accuracy in all applications considered. A hybrid architecture, extending the LSTM with a mean stochastic model (MSM-LSTM), is proposed to ensure convergence to the invariant measure. This novel hybrid method is fully data-driven and extends the forecasting capabilities of LSTM networks., Competing Interests: We declare we have no competing interests.

Background: The diagnostic assessment of fetal arrhythmias relies on the measurements of atrioventricular (AV) and ventriculoatrial (VA) time intervals. Pulsed Doppler over in- and outflow of the left ventricle and tissue Doppler imaging are well-described methods, while Doppler measurements between the left brachiocephalic vein and the aortic arch are less investigated. The aim of this study was to compare these methods of measurement, to find influencing factors on AV and VA times and their ratio, and to create reference ranges., Methods: Echocardiography was performed between 16 and 40 weeks of gestation in normal singleton pregnancies. Nomograms for the individual measurements were created using quantile regression with Matlab Data Analytics. Statistical analyses were performed with GraphPad version 5.0 for Windows., Results: A total of 329 pregnant women were enrolled. A significant correlation exists between AV and VA times and gestational age (GA) (p = 0.0104 to <0.0001, σ = 0.1412 to 0.3632). No correlation was found between the AV:VA ratio and GA (p = 0.08 to 0.60). All measurements differed significantly amongst the studied methods (p < 0.0001)., Conclusions: AV and VA intervals increase proportionally with GA; no other independent influencing factors could be identified. As significant differences exist between the three methods of assessment, it is crucial to use appropriate reference ranges to diagnose pathologies., (© 2017 S. Karger AG, Basel.)

The Lennard-Jones (LJ) potential is a cornerstone of Molecular Dynamics (MD) simulations and among the most widely used computational kernels in science. The LJ potential models atomistic attraction and repulsion with century old prescribed parameters (q = 6, p = 12, respectively), originally related by a factor of two for simplicity of calculations. We propose the inference of the repulsion exponent through Hierarchical Bayesian uncertainty quantification We use experimental data of the radial distribution function and dimer interaction energies from quantum mechanics simulations. We find that the repulsion exponent p ≈ 6.5 provides an excellent fit for the experimental data of liquid argon, for a range of thermodynamic conditions, as well as for saturated argon vapour. Calibration using the quantum simulation data did not provide a good fit in these cases. However, values p ≈ 12.7 obtained by dimer quantum simulations are preferred for the argon gas while lower values are promoted by experimental data. These results show that the proposed LJ 6-p potential applies to a wider range of thermodynamic conditions, than the classical LJ 6-12 potential. We suggest that calibration of the repulsive exponent in the LJ potential widens the range of applicability and accuracy of MD simulations.

The coordinated motion by multiple swimmers is a fundamental component in fish schooling. The flow field induced by the motion of each self-propelled swimmer implies non-linear hydrodynamic interactions among the members of a group. How do swimmers compensate for such hydrodynamic interactions in coordinated patterns? We provide an answer to this riddle though simulations of two, self-propelled, fish-like bodies that employ a learning algorithm to synchronise their swimming patterns. We distinguish between learned motion patterns and the commonly used a-priori specified movements, that are imposed on the swimmers without feedback from their hydrodynamic interactions. First, we demonstrate that two rigid bodies executing pre-specified motions, with an alternating leader and follower, can result in substantial drag-reduction and intermittent thrust generation. In turn, we study two self-propelled swimmers arranged in a leader-follower configuration, with a-priori specified body-deformations. These two self-propelled swimmers do not sustain their tandem configuration. The follower experiences either an increase or decrease in swimming speed, depending on the initial conditions, while the swimming of the leader remains largely unaffected. This indicates that a-priori specified patterns are not sufficient to sustain synchronised swimming. We then examine a tandem of swimmers where the leader has a steady gait and the follower learns to synchronize its motion, to overcome the forces induced by the leader's vortex wake. The follower employs reinforcement learning to adapt its swimming-kinematics so as to minimize its lateral deviations from the leader's path. Swimming in such a sustained synchronised tandem yields up to [Formula: see text] reduction in energy expenditure for the follower, in addition to a [Formula: see text] increase in its swimming-efficiency. The present results show that two self-propelled swimmers can be synchronised by adapting their motion patterns to compensate for flow-structure interactions. Moreover, swimmers can exploit the vortical structures of their flow field so that synchronised swimming is energetically beneficial.

We propose a hierarchical Bayesian framework to systematically integrate heterogeneous data for the calibration of force fields in Molecular Dynamics (MD) simulations. Our approach enables the fusion of diverse experimental data sets of the physico-chemical properties of a system at different thermodynamic conditions. We demonstrate the value of this framework for the robust calibration of MD force-fields for water using experimental data of its diffusivity, radial distribution function, and density. In order to address the high computational cost associated with the hierarchical Bayesian models, we develop a novel surrogate model based on the empirical interpolation method. Further computational savings are achieved by implementing a highly parallel transitional Markov chain Monte Carlo technique. The present method bypasses possible subjective weightings of the experimental data in identifying MD force-field parameters.

The increasing power density and the decreasing dimensions of transistors present severe thermal challenges to the design of modern microprocessors. Furthermore, new technologies such as three-dimensional chip-stack architectures require novel cooling solutions for their thermal management. Here, we demonstrate, through transient heat-dissipation simulations, that a covalently bonded graphene-carbon nanotube (G-CNT) hybrid immersed in water is a promising solution for the ultrafast cooling of such high-temperature and high heat-flux surfaces. The G-CNT hybrid offers a unique platform to integrate the superior axial heat transfer capability of individual CNTs via their parallel arrangement. The immersion of the G-CNT in water enables an additional heat dissipation path via the solid-liquid interaction, allowing for the sustainable cooling of the hot surface under a constant power input of up to 10 000 W cm -2 .

Growth of tissues is highly reproducible; yet, growth of individual cells in a tissue is highly variable, and neighboring cells can grow at different rates. We analyzed the growth of epidermal cell lineages in the Arabidopsis (Arabidopsis thaliana) sepal to determine how the growth curves of individual cell lineages relate to one another in a developing tissue. To identify underlying growth trends, we developed a continuous displacement field to predict spatially averaged growth rates. We showed that this displacement field accurately describes the growth of sepal cell lineages and reveals underlying trends within the variability of in vivo cellular growth. We found that the tissue, individual cell lineages, and cell walls all exhibit growth rates that are initially low, accelerate to a maximum, and decrease again. Accordingly, these growth curves can be represented by sigmoid functions. We examined the relationships among the cell lineage growth curves and surprisingly found that all lineages reach the same maximum growth rate relative to their size. However, the cell lineages are not synchronized; each cell lineage reaches this same maximum relative growth rate but at different times. The heterogeneity in observed growth results from shifting the same underlying sigmoid curve in time and scaling by size. Thus, despite the variability in growth observed in our study and others, individual cell lineages in the developing sepal follow similarly shaped growth curves., (© 2015 American Society of Plant Biologists. All Rights Reserved.)

The Kapitza resistance (RK) between few-layer graphene (FLG) and water was studied using molecular dynamics simulations. The RK was found to depend on the number of the layers in the FLG though, surprisingly, not on the water block thickness. This distinct size dependence is attributed to the large difference in the phonon mean free path between the FLG and water. Remarkably, RK is strongly dependent on the layering of water adjacent to the FLG, exhibiting an inverse proportionality relationship to the peak density of the first water layer, which is consistent with better acoustic phonon matching between FLG and water. These findings suggest novel ways to engineer the thermal transport properties of solid-liquid interfaces by controlling and regulating the liquid layering at the interface.

Rough surfaces immersed under water remain practically dry if the liquid-solid contact is on roughness peaks, while the roughness valleys are filled with gas. Mechanisms that prevent water from invading the valleys are well studied. However, to remain practically dry under water, additional mechanisms need consideration. This is because trapped gas (e.g. air) in the roughness valleys can dissolve into the water pool, leading to invasion. Additionally, water vapor can also occupy the roughness valleys of immersed surfaces. If water vapor condenses, that too leads to invasion. These effects have not been investigated, and are critically important to maintain surfaces dry under water. In this work, we identify the critical roughness scale, below which it is possible to sustain the vapor phase of water and/or trapped gases in roughness valleys - thus keeping the immersed surface dry. Theoretical predictions are consistent with molecular dynamics simulations and experiments.

Morphogenesis emerges from complex multiscale interactions between genetic and mechanical processes. To understand these processes, the evolution of cell shape, proliferation and gene expression must be quantified. This quantification is usually performed either in full 3D, which is computationally expensive and technically challenging, or on 2D planar projections, which introduces geometrical artifacts on highly curved organs. Here we present MorphoGraphX ( www.MorphoGraphX.org), a software that bridges this gap by working directly with curved surface images extracted from 3D data. In addition to traditional 3D image analysis, we have developed algorithms to operate on curved surfaces, such as cell segmentation, lineage tracking and fluorescence signal quantification. The software's modular design makes it easy to include existing libraries, or to implement new algorithms. Cell geometries extracted with MorphoGraphX can be exported and used as templates for simulation models, providing a powerful platform to investigate the interactions between shape, genes and growth.

The mouse retina has become a prominent model for studying angiogenesis. The easy access and well-known developmental progression have significantly propelled our ability to examine and manipulate blood vessels in vivo. Nonetheless, most studies have restricted their evaluations to the superficial plexus (an upper vascular layer in contact with the vitreous). Here we present experimental data and quantification for the developmental progression of the full retina including the intermediate and deeper plexus that sprouts from the superficial layer. We analyze the origin and advancement of vertical sprouting and present the progression of vascular perfusion within the tissue. Furthermore, we introduce the use of Minkowsky functionals to quantify remodeling in the superficial and deeper plexus. The work expands information on the retina towards a 3D structure. This is of particular interest, as recent data have demonstrated differential effects of gene deletion on the upper and deeper plexus, highlighting the concept of distinct operational pathways during sprouting angiogenesis.

For over five decades, molecular dynamics (MD) simulations have helped to elucidate critical mechanisms in a broad range of physiological systems and technological innovations. MD simulations are synergetic with experiments, relying on measurements to calibrate their parameters and probing "what if scenarios" for systems that are difficult to investigate experimentally. However, in certain systems, such as nanofluidics, the results of experiments and MD simulations differ by several orders of magnitude. This discrepancy may be attributed to the spatiotemporal scales and structural information accessible by experiments and simulations. Furthermore, MD simulations rely on parameters that are often calibrated semiempirically, while the effects of their computational implementation on their predictive capabilities have only been sporadically probed. In this work, we show that experimental and MD investigations can be consolidated through a rigorous uncertainty quantification framework. We employ a Bayesian probabilistic framework for large scale MD simulations of graphitic nanostructures in aqueous environments. We assess the uncertainties in the MD predictions for quantities of interest regarding wetting behavior and hydrophobicity. We focus on three representative systems: water wetting of graphene, the aggregation of fullerenes in aqueous solution, and the water transport across carbon nanotubes. We demonstrate that the dominant mode of calibrating MD potentials in nanoscale fluid mechanics, through single values of water contact angle on graphene, leads to large uncertainties and fallible quantitative predictions. We demonstrate that the use of additional experimental data reduces uncertainty, improves the predictive accuracy of MD models, and consolidates the results of experiments and simulations.

Understanding and accelerating the mechanisms of endothelial wound healing is of fundamental interest for biotechnology and of significant medical utility in repairing pathologic changes to the vasculature induced by invasive medical interventions. We report the fundamental mechanisms that determine the influence of substrate topography and flow on the efficiency of endothelial regeneration. We exposed endothelial monolayers, grown on topographically engineered substrates (gratings), to controlled levels of flow-induced shear stress. The wound healing dynamics were recorded and analyzed in various configurations, defined by the relative orientation of an inflicted wound, the topography and the flow direction. Under flow perpendicular to the wound, the speed of endothelial regeneration was significantly increased on substrates with gratings oriented in the direction of the flow when compared to flat substrates. This behavior is linked to the dynamic state of cell-to-cell adhesions in the monolayer. In particular, interactions with the substrate topography counteract Vascular Endothelial Cadherin phosphorylation induced by the flow and the wounding. This effect contributes to modulating the mechanical connection between migrating cells to an optimal level, increasing their coordination and resulting in coherent cell motility and preservation of the monolayer integrity, thus accelerating wound healing. We further demonstrate that the reduction of vascular endothelial cadherin phosphorylation, through specific inhibition of Src activity, enhances endothelial wound healing in flows over flat substrates., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Cell migration is commonly quantified by tracking the speed of the cell layer interface in wound healing assays. This quantification is often hampered by low signal to noise ratio, in particular when complex substrates are employed to emulate in vivo cell migration in geometrically complex environments. Moreover, information about the cell motion, readily available inside the migrating cell layers, is not usually harvested. We introduce Cell Image Velocimetry (CIV), a combination of cell layer segmentation and image velocimetry algorithms, to drastically enhance the quantification of cell migration by wound healing assays. The resulting software analyses the speed of the interface as well as the detailed velocity field inside the cell layers in an automated fashion. CIV is shown to be highly robust for images with low signal to noise ratio, low contrast and frame shifting and it is portable across various experimental settings. The modular design and parametrization of CIV is not restricted to wound healing assays and allows for the exploration and quantification of flow phenomena in any optical microscopy dataset. Here, we demonstrate the capabilities of CIV in wound healing assays over topographically engineered surfaces and quantify the relative merits of differently aligned gratings on cell migration.

We present a Bayesian probabilistic framework for quantifying and propagating the uncertainties in the parameters of force fields employed in molecular dynamics (MD) simulations. We propose a highly parallel implementation of the transitional Markov chain Monte Carlo for populating the posterior probability distribution of the MD force-field parameters. Efficient scheduling algorithms are proposed to handle the MD model runs and to distribute the computations in clusters with heterogeneous architectures. Furthermore, adaptive surrogate models are proposed in order to reduce the computational cost associated with the large number of MD model runs. The effectiveness and computational efficiency of the proposed Bayesian framework is demonstrated in MD simulations of liquid and gaseous argon.

Particle-mesh interpolations are fundamental operations for particle-in-cell codes, as implemented in vortex methods, plasma dynamics and electrostatics simulations. In these simulations, the mesh is used to solve the field equations and the gradients of the fields are used in order to advance the particles. The time integration of particle trajectories is performed through an extensive resampling of the flow field at the particle locations. The computational performance of this resampling turns out to be limited by the memory bandwidth of the underlying computer architecture. We investigate how mesh-particle interpolation can be efficiently performed on graphics processing units (GPUs) and multicore central processing units (CPUs), and we present two implementation techniques. The single-precision results for the multicore CPU implementation show an acceleration of 45-70×, depending on system size, and an acceleration of 85-155× for the GPU implementation over an efficient single-threaded C++ implementation. In double precision, we observe a performance improvement of 30-40× for the multicore CPU implementation and 20-45× for the GPU implementation. With respect to the 16-threaded standard C++ implementation, the present CPU technique leads to a performance increase of roughly 2.8-3.7× in single precision and 1.7-2.4× in double precision, whereas the GPU technique leads to an improvement of 9× in single precision and 2.2-2.8× in double precision.

We use the Dpp morphogen gradient in the Drosophila wing disc as a model to address the fundamental question of how a gradient of a growth factor can produce uniform growth. We first show that proper expression and subcellular localization of components in the Fat tumor-suppressor pathway, which have been argued to depend on Dpp activity differences, are not reliant on the Dpp gradient. We next analyzed cell proliferation in discs with uniformly high Dpp or uniformly low Fat signaling activity and found that these pathways regulate growth in a complementary manner. While the Dpp mediator Brinker inhibits growth in the primordium primarily in the lateral regions, Fat represses growth mostly in the medial region. Together, our results indicate that the activities of both signaling pathways are regulated in a parallel rather than sequential manner and that uniform proliferation is achieved by their complementary action on growth., (Copyright © 2011 Elsevier Inc. All rights reserved.)

This study investigated specific gene expression profiles in patients with donor-specific cytotoxic-hyporesponsiveness, reflected by cytotoxic T-lymphocyte precursor frequency (CTLpf). The effect of calcineurin inhibitor (CNI) withdrawal was studied on markers for cytotoxicity (perforin, granzyme B), apoptosis (Fas,FasL), Th1 and Th2 cytokines (IL-2, IL-10), Th1 and Th2 transcription factors (T-bet, GATA 3), Th17 transcription factor and cytokine (RORγt, IL-17), and for immune regulation/activation (CD25, FOXP3). Peripheral blood samples from renal allograft recipients (n = 18) more than two yr after transplantation with stable renal function were analyzed before and four months after CNI withdrawal. Additionally, systolic and diastolic blood pressure, cholesterol, serum creatinine and proteinuria were evaluated, and no significant differences were measured before and after CNI withdrawal. However, CNIs' discontinuation influenced peripheral gene expression profiles. After CNI withdrawal, the mRNA expression of Granzyme B, Perforin, Fas, FasL, T-bet, GATA3 and CD25 were significantly lower than during CNI treatment. After CNI discontinuation, donor-specific CTLpf decreased, while FOXP3 expression discriminated between detectable and non-detectable donor-specific cytolysis reactivity; FOXP3 transcript values were highest in absence of donor-specific cytotoxicity (p < 0.01). Our study shows CNI withdrawal in stable kidney transplant recipients two yr after transplantation is safe. Moreover, discontinuation of CNIs' treatment allows FOXP3+ regulatory T-cells development, resulting in a significant decrease of anti-donor immune reactivity., (© 2010 John Wiley & Sons A/S.)

We present a simple algorithm for the simulation of stiff, discrete-space, continuous-time Markov processes. The algorithm is based on the concept of flow averaging for the integration of stiff ordinary and stochastic differential equations and ultimately leads to a straightforward variation of the the well-known stochastic simulation algorithm (SSA). The speedup that can be achieved by the present algorithm [flow averaging integrator SSA (FLAVOR-SSA)] over the classical SSA comes naturally at the expense of its accuracy. The error of the proposed method exhibits a cutoff phenomenon as a function of its speed-up, allowing for optimal tuning. Two numerical examples from chemical kinetics are provided to illustrate the efficiency of the method.

Reaction-diffusion processes on complex deforming surfaces are fundamental to a number of biological processes ranging from embryonic development to cancer tumor growth and angiogenesis. The simulation of these processes using continuum reaction-diffusion models requires computational methods capable of accurately tracking the geometric deformations and discretizing on them the governing equations. We employ a Lagrangian level-set formulation to capture the deformation of the geometry and use an embedding formulation and an adaptive particle method to discretize both the level-set equations and the corresponding reaction-diffusion. We validate the proposed method and discuss its advantages and drawbacks through simulations of reaction-diffusion equations on complex and deforming geometries.

We study shear-induced instabilities of lipid bilayers immersed in water using coarse-grained molecular dynamics simulations. The shear imposed by the flow of the water induces initially microscopic structural changes of the membrane, starting with tilting of the molecules in the direction of the shear. The tilting propagates in the spanwise direction when the shear rate exceeds a critical value and the membrane undergoes a bucklinglike deformation in the direction perpendicular to the shear. The bucklinglike undulation continues until a localized Kelvin-Helmholtz-like instability leads to membrane rupture. We study the different modes of membrane undulation using membranes of different geometries and quantify the relative importance of the bucklinglike bending and the Kelvin-Helmholtz-like instability of the membrane.

Cytoplasmic transport of organelles, nucleic acids and proteins on microtubules is usually bidirectional with dynein and kinesin motors mediating the delivery of cargoes in the cytoplasm. Here we combine live cell microscopy, single virus tracking and trajectory segmentation to systematically identify the parameters of a stochastic computational model of cargo transport by molecular motors on microtubules. The model parameters are identified using an evolutionary optimization algorithm to minimize the Kullback-Leibler divergence between the in silico and the in vivo run length and velocity distributions of the viruses on microtubules. The present stochastic model suggests that bidirectional transport of human adenoviruses can be explained without explicit motor coordination. The model enables the prediction of the number of motors active on the viral cargo during microtubule-dependent motions as well as the number of motor binding sites, with the protein hexon as the binding site for the motors.

An exact method for stochastic simulation of chemical reaction networks, which accelerates the stochastic simulation algorithm (SSA), is proposed. The present "ER-leap" algorithm is derived from analytic upper and lower bounds on the multireaction probabilities sampled by SSA, together with rejection sampling and an adaptive multiplicity for reactions. The algorithm is tested on a number of well-quantified reaction networks and is found experimentally to be very accurate on test problems including a chaotic reaction network. At the same time ER-leap offers a substantial speedup over SSA with a simulation time proportional to the 23 power of the number of reaction events in a Galton-Watson process.

We present a multiscale algorithm to couple atomistic water models with continuum incompressible flow simulations via a Schwarz domain decomposition approach. The coupling introduces an inhomogeneity in the description of the atomistic domain and prevents the use of periodic boundary conditions. The use of a mass conserving specular wall results in turn to spurious oscillations in the density profile of the atomistic description of water. These oscillations can be eliminated by using an external boundary force that effectively accounts for the virial component of the pressure. In this Rapid Communication, we extend a control algorithm, previously introduced for monatomic molecules, to the case of atomistic water and demonstrate the effectiveness of this approach. The proposed computational method is validated for the cases of equilibrium and Couette flow of water.