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:Auxetic metamaterials possess unnatural properties, such as a negative Poisson's ratio, which offers interesting features when combined with traditional materials. This paper describes the deformation behavior of a gel consisting of spherical auxetic inclusions when embedded in a conventional matrix. The auxetic inclusions and conventional matrix were modeled as spherical objects with a controlled pore shape. The auxetic particle had a reentrant honeycomb, and the conventional phase contained honeycomb-shaped pores. The deformation behavior was simulated using various existing models based on continuum mechanics. For the continuum mechanics models-the simplest of which are the Mori-Tanaka theory and self-consistent field mechanics models-the auxetic particle was homogenized as a solid element with Young's modulus and Poisson's ratio and compared with the common composite gel filled with rigid spheres. The finite element analysis simulations using these models were performed for two cases: (1) a detailed model of one particle and its surroundings in which the structure included the design of both the reentrant and conventional honeycombs; and (2) a multiparticle face-centered cubic lattice where both the classic matrix and auxetic particle were homogenized. Our results suggest that auxetic inclusion-filled gels provide an unsurpassed balance of low density and enhanced stiffness.
| S-EPMC9689618 | biostudies-literature
Project description:This article shows that the gas-phase entropy of molecules is proportional to the area of the molecules, with corrections for the different curvatures of the molecular surface. The ability to estimate gas-phase entropy by the area law also allows us to calculate molecular entropy faster and more accurately than currently popular methods of estimating molecular entropy with harmonic oscillator approximation. The speed and accuracy of our method will open up new possibilities for the explicit inclusion of entropy in various computational biology methods.
| S-EPMC11353761 | 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:Parvalbumin (PV) is a globular calcium (Ca(2+))-selective protein expressed in a variety of biological tissues. Our computational studies of the rat β-parvalbumin (β-PV) isoform seek to elucidate the molecular thermodynamics of Ca(2+) versus magnesium (Mg(2+)) binding at the protein's two EF-hand motifs. Specifically, we have utilized molecular dynamics (MD) simulations and a mean-field electrolyte model (mean spherical approximation (MSA) theory) to delineate how the EF-hand scaffold controls the "local" thermodynamics of Ca(2+) binding selectivity over Mg(2+). Our MD simulations provide the probability density of metal-chelating oxygens within the EF-hand scaffolds for both Ca(2+) and Mg(2+), as well the conformational strain induced by Mg(2+) relative to Ca(2+) binding. MSA theory utilizes the binding domain oxygen and charge distributions to predict the chemical potential of ion binding, as well as their corresponding concentrations within the binding domain. We find that the electrostatic and steric contributions toward ion binding were similar for Mg(2+) and Ca(2+), yet the latter was 5.5 kcal/mol lower in enthalpy when internal strain within the EF hand was considered. We therefore speculate that beyond differences in dehydration energies for the Ca(2+) versus Mg(2+), strain induced in the β-PV EF hand by cation binding significantly contributes to the nearly 10,000-fold difference in binding affinity reported in the literature. We further complemented our analyses of local factors governing cation binding selectivity with whole-protein (global) contributions, such as interhelical residue-residue contacts and solvent exposure of hydrophobic surface. These contributions were found to be comparable for both Ca(2+)- and Mg(2+)-bound β-PV, which may implicate local factors, EF-hand strain, and dehydration, in providing the primary means of selectivity. We anticipate these methods could be used to estimate metal binding thermodynamics across a broad range of PV sequence homologues and EF-hand-containing, Ca(2+) binding proteins.
| S-EPMC8867401 | 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:Numerous biological systems contain vesicle-like biomolecular compartments without membranes, which contribute to diverse functions including gene regulation, stress response, signaling, and skin barrier formation. Coacervation, as a form of liquid-liquid phase separation (LLPS), is recognized as a representative precursor to the formation and assembly of membrane-less vesicle-like structures, although their formation mechanism remains unclear. In this study, a coacervation-driven membrane-less vesicle-like structure is constructed using two proteins, GG1234 (an anionic intrinsically disordered protein) and bhBMP-2 (a bioengineered human bone morphogenetic protein 2). GG1234 formed both simple coacervates by itself and complex coacervates with the relatively cationic bhBMP-2 under acidic conditions. Upon addition of dissolved bhBMP-2 to the simple coacervates of GG1234, a phase transition from spherical simple coacervates to vesicular condensates occurred via the interactions between GG1234 and bhBMP-2 on the surface of the highly viscoelastic GG1234 simple coacervates. Furthermore, the shell structure in the outer region of the GG1234/bhBMP-2 vesicular condensates exhibited gel-like properties, leading to the formation of multiphasic vesicle-like compartments. A potential mechanism is proposed for the formation of the membrane-less GG1234/bhBMP-2 vesicle-like compartments. This study provides a dynamic process underlying the formation of biomolecular multiphasic condensates, thereby enhancing the understanding of these biomolecular structures.
| S-EPMC10870063 | 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