Fast Simulation of Lipid Vesicle Deformation Using Spherical Harmonic Approximation.
Ontology highlight
ABSTRACT:
SUBMITTER: Mikucki M
PROVIDER: S-EPMC5552105 | biostudies-other | 2017 Jan
REPOSITORIES: biostudies-other
Similar Datasets
Project description:In cognitive neuroscience, computational modeling can formally adjudicate between theories and affords quantitative fits to behavioral/brain data. Pragmatically, however, the space of plausible generative models considered is dramatically limited by the set of models with known likelihood functions. For many models, the lack of a closed-form likelihood typically impedes Bayesian inference methods. As a result, standard models are evaluated for convenience, even when other models might be superior. Likelihood-free methods exist but are limited by their computational cost or their restriction to particular inference scenarios. Here, we propose neural networks that learn approximate likelihoods for arbitrary generative models, allowing fast posterior sampling with only a one-off cost for model simulations that is amortized for future inference. We show that these methods can accurately recover posterior parameter distributions for a variety of neurocognitive process models. We provide code allowing users to deploy these methods for arbitrary hierarchical model instantiations without further training.
| S-EPMC8102064 | biostudies-literature
Project description:Stochastic channel gating is the major source of intrinsic neuronal noise whose functional consequences at the microcircuit- and network-levels have been only partly explored. A systematic study of this channel noise in large ensembles of biophysically detailed model neurons calls for the availability of fast numerical methods. In fact, exact techniques employ the microscopic simulation of the random opening and closing of individual ion channels, usually based on Markov models, whose computational loads are prohibitive for next generation massive computer models of the brain. In this work, we operatively define a procedure for translating any Markov model describing voltage- or ligand-gated membrane ion-conductances into an effective stochastic version, whose computer simulation is efficient, without compromising accuracy. Our approximation is based on an improved Langevin-like approach, which employs stochastic differential equations and no Montecarlo methods. As opposed to an earlier proposal recently debated in the literature, our approximation reproduces accurately the statistical properties of the exact microscopic simulations, under a variety of conditions, from spontaneous to evoked response features. In addition, our method is not restricted to the Hodgkin-Huxley sodium and potassium currents and is general for a variety of voltage- and ligand-gated ion currents. As a by-product, the analysis of the properties emerging in exact Markov schemes by standard probability calculus enables us for the first time to analytically identify the sources of inaccuracy of the previous proposal, while providing solid ground for its modification and improvement we present here.
| S-EPMC3053314 | biostudies-literature
Project description:We present hierarchical spherical deformation for a group-wise shape correspondence to address template selection bias and to minimize registration distortion. In this work, we aim at a continuous and smooth deformation field to guide accurate cortical surface registration. In conventional spherical registration methods, a global rigid alignment and local deformation are independently performed. Motivated by the composition of precession and intrinsic rotation, we simultaneously optimize global rigid rotation and non-rigid local deformation by utilizing spherical harmonics interpolation of local composite rotations in a single framework. To this end, we indirectly encode local displacements by such local composite rotations as functions of spherical locations. Furthermore, we introduce an additional regularization term to the spherical deformation, which maximizes its rigidity while reducing registration distortion. To improve surface registration performance, we employ the second order approximation of the energy function that enables fast convergence of the optimization. In the experiments, we validate our method on healthy normal subjects with manual cortical surface parcellation in registration accuracy and distortion. We show an improved shape correspondence with high accuracy in cortical surface parcellation and significantly low registration distortion in surface area and edge length. In addition to validation, we discuss parameter tuning, optimization, and implementation design with potential acceleration.
| S-EPMC6733638 | biostudies-literature
Project description:PurposeTo describe and implement a strategy for dynamic slice-by-slice and multiband B0 shimming using spherical harmonic shims in the human brain at 7T.TheoryFor thin axial slices, spherical harmonic shims can be divided into pairs of shims (z-degenerate and non-z-degenerate) that are spatially degenerate, such that only ½ of the shims (non-z-degenerate) are required for single slice optimizations. However, when combined, the pairs of shims can be used to simultaneously generate the same in-plane symmetries but with different amplitudes as a function of their z location. This enables multiband shimming equivalent to that achievable by single slice-by-slice optimization.MethodsAll data were acquired at 7T using a spherical harmonic shim insert enabling shimming up through 4th order with two additional 5th order shims (1st-4th+). Dynamic shim updating was achieved using a 10A shim power supply with 2 ms ramps and constrained optimizations to minimize eddy currents.ResultsIn groups of eight subjects, we demonstrated that: 1) dynamic updating using 1st-4th+ order shims reduced the SD of the B0 field over the whole brain from 32.4 ± 2.6 and 24.9 ± 2 Hz with 1st-2nd and 1st-4th+ static global shimming to 15.1 ± 1.7 Hz; 2) near equivalent performance was achieved when dynamically updating only the non-z-degenerate shims (14.3 ± 1.5 Hz), or when a using multiband shim factor of 2, MBs = 2, and all shims (14.4 ± 2.0 Hz).ConclusionHigh order spherical harmonics provide substantial improvements over static global shimming and enable dynamic multiband shimming with near equivalent performance to that of dynamic slice-by-slice shimming. This reduces distortion in echo planar imaging.
| S-EPMC8098050 | biostudies-literature
Project description:This study addresses the nonlinear structure-extended cavity interaction simulation using a new version of the multilevel residue harmonic balance method. This method has only been adopted once to solve a nonlinear beam problem. This is the first study to use this method to solve a nonlinear structural acoustic problem. This study has two focuses: 1) the new version of the multilevel residue harmonic balance method can generate the higher-level nonlinear solutions ignored in the previous version and 2) the effect of the extended cavity, which has not been considered in previous studies, is examined. The cavity length of a panel-cavity system is sometimes longer than the panel length. However, many studies have adopted a model in which the cavity length is equal to the panel length. The effects of excitation magnitude, cavity depth, damping and number of structural modes on sound and vibration responses are investigated for various panel cases. In the simulations, the present harmonic balance solutions agree reasonably well with those obtained from the classical harmonic balance method. There are two important findings. First, the nonlinearity of a structural acoustic system highly depends on the cavity size. If the cavity size is smaller, the nonlinearity is higher. A large cavity volume implies a low stiffness or small acoustic pressure transmitted from the source panel to the nonlinear panel. In other words, the additional volume in an extended cavity affects the nonlinearity, sound and vibration responses of a structural acoustic system. Second, if an acoustic resonance couples with a structural resonance, nonlinearity is amplified and thus the insertion loss is adversely affected.
| S-EPMC6029766 | biostudies-literature
Project description:Styryl dyes and genetically encoded pH-sensitive fluorescent proteins like pHluorin are well-established tools for the optical analysis of synaptic vesicle (SV) recycling at presynaptic boutons. Here, we describe the development of a new class of fluorescent probes based on pH-sensitive organic dyes covalently bound to lipids, providing a promising complementary assay to genetically encoded fluorescent probes. These new optical tracers allow a pure read out of membrane turnover during synaptic activity and visualization of multiple rounds of stimulation-dependent SV recycling without genetic perturbation. Measuring the incorporation efficacy of different dye-labeled lipids into budding SVs, we did not observe an enrichment of lipids with affinity for liquid ordered membrane domains. But most importantly, we found no evidence for a static segregation of SVs into recycling and resting pools. A small but significant fraction of SVs that is reluctant to release during a first round of evoked activity can be exocytosed during a second bout of stimulation, showing fast intermixing of SV pools within seconds. Furthermore, we found that SVs recycling spontaneously have a higher chance to re-occupy release sites than SVs recycling during high-frequency evoked activity. In summary, our data provide strong evidence for a highly dynamic and use-dependent control of the fractions of releasable or resting SVs.
| S-EPMC5801418 | biostudies-literature
Project description:Model-based analyses of natural selection often categorize sites into a relatively small number of site classes. Forcing each site to belong to one of these classes places unrealistic constraints on the distribution of selection parameters, which can result in misleading inference due to model misspecification. We present an approximate hierarchical Bayesian method using a Markov chain Monte Carlo (MCMC) routine that ensures robustness against model misspecification by averaging over a large number of predefined site classes. This leaves the distribution of selection parameters essentially unconstrained, and also allows sites experiencing positive and purifying selection to be identified orders of magnitude faster than by existing methods. We demonstrate that popular random effects likelihood methods can produce misleading results when sites assigned to the same site class experience different levels of positive or purifying selection--an unavoidable scenario when using a small number of site classes. Our Fast Unconstrained Bayesian AppRoximation (FUBAR) is unaffected by this problem, while achieving higher power than existing unconstrained (fixed effects likelihood) methods. The speed advantage of FUBAR allows us to analyze larger data sets than other methods: We illustrate this on a large influenza hemagglutinin data set (3,142 sequences). FUBAR is available as a batch file within the latest HyPhy distribution (http://www.hyphy.org), as well as on the Datamonkey web server (http://www.datamonkey.org/).
| S-EPMC3670733 | biostudies-literature
Project description:We performed molecular dynamics simulation of nanoindentation on Cu/Ni nanotwinned multilayer films using a spherical indenter, aimed to investigate the effects of hetero-twin interface and twin thickness on hardness. We found that both twinning partial slip (TPS) and partial slip parallel with twin boundary (PSPTB) can reduce hardness and therefore should not be ignored when evaluating mechanical properties at nanoscale. There is a critical range of twin thickness ? (~25?Å?<???<?~31?Å), in which hardness of the multilayer films is maximized. At a smaller ?, TPSs appear due to the reaction between partial dislocations and twin boundary accounts for the softening-dominated mechanism. We also found that the combination of the lowered strengthening due to confined layer slips and the softening due to TPSs and PSPTBs results in lower hardness at a larger ?.
| S-EPMC5073369 | biostudies-literature
Project description:We develop an analytical theory to explain the experimentally observed morphological transitions of quasispherical giant vesicles induced by alternating electric fields. The model treats the inner and suspending media as lossy dielectrics, and the membrane as an impermeable flexible incompressible-fluid sheet. The vesicle shape is obtained by balancing electric, hydrodynamic, bending, and tension stresses exerted on the membrane. Our approach, which is based on force balance, also allows us to describe the time evolution of the vesicle deformation, in contrast to earlier works based on energy minimization, which are able to predict only stationary shapes. Our theoretical predictions for vesicle deformation are consistent with experiment. If the inner fluid is more conducting than the suspending medium, the vesicle always adopts a prolate shape. In the opposite case, the vesicle undergoes a transition from a prolate to oblate ellipsoid at a critical frequency, which the theory identifies with the inverse membrane charging time. At frequencies higher than the inverse Maxwell-Wagner polarization time, the electrohydrodynamic stresses become too small to alter the vesicle's quasispherical rest shape. The model can be used to rationalize the transient and steady deformation of biological cells in electric fields.
| S-EPMC2712034 | biostudies-literature
Project description:In this manuscript we expand significantly on our earlier communication by investigating the bilayer self-assembly of eight different types of phospholipids in unbiased molecular dynamics (MD) simulations using three widely used all-atom lipid force fields. Irrespective of the underlying force field, the lipids are shown to spontaneously form stable lamellar bilayer structures within 1 microsecond, the majority of which display properties in satisfactory agreement with the experimental data. The lipids self-assemble via the same general mechanism, though at formation rates that differ both between lipid types, force fields and even repeats on the same lipid/force field combination. In addition to zwitterionic phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipids, anionic phosphatidylserine (PS) and phosphatidylglycerol (PG) lipids are represented. To our knowledge this is the first time bilayer self-assembly of phospholipids with negatively charged head groups is demonstrated in all-atom MD simulations.
| S-EPMC5091658 | biostudies-literature