Baroreflex Curve Fitting Using a WYSIWYG Boltzmann Sigmoidal Equation
Ontology highlight
ABSTRACT:
SUBMITTER: Heusser K
PROVIDER: S-EPMC8519000 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Sigmoidal and double-sigmoidal dynamics are commonly observed in many areas of biology. Here we present sicegar, an R package for the automated fitting and classification of sigmoidal and double-sigmoidal data. The package categorizes data into one of three categories, "no signal," "sigmoidal," or "double-sigmoidal," by rigorously fitting a series of mathematical models to the data. The data is labeled as "ambiguous" if neither the sigmoidal nor double-sigmoidal model fit the data well. In addition to performing the classification, the package also reports a wealth of metrics as well as biologically meaningful parameters describing the sigmoidal or double-sigmoidal curves. In extensive simulations, we find that the package performs well, can recover the original dynamics even under fairly high noise levels, and will typically classify curves as "ambiguous" rather than misclassifying them. The package is available on CRAN and comes with extensive documentation and usage examples.
| S-EPMC5774301 | biostudies-literature
Project description:Quantitative real-time PCR has revolutionized many aspects of genetic research, biomedical diagnostics and pathogen detection. Nevertheless, the full potential of this technology has yet to be realized, primarily due to the limitations of the threshold-based methodologies that are currently used for quantitative analysis. Prone to errors caused by variations in reaction preparation and amplification conditions, these approaches necessitate construction of standard curves for each target sequence, significantly limiting the development of high-throughput applications that demand substantive levels of reliability and automation. In this study, an alternative approach based upon fitting of fluorescence data to a four-parametric sigmoid function is shown to dramatically increase both the utility and reliability of quantitative real-time PCR. By mathematically modeling individual amplification reactions, quantification can be achieved without the use of standard curves and without prior knowledge of amplification efficiency. Combined with provision of quantitative scale via optical calibration, sigmoidal curve-fitting could confer the capability for fully automated quantification of nucleic acids with unparalleled accuracy and reliability.
| S-EPMC545475 | biostudies-literature
Project description:Microtubules are highly negatively charged proteins which have been shown to behave as bio-nanowires capable of conducting ionic currents. The electrical characteristics of microtubules are highly complicated and have been the subject of previous work; however, the impact of the ionic concentration of the buffer solution on microtubule electrical properties has often been overlooked. In this work we use the non-linear Poisson Boltzmann equation, modified to account for a variable permittivity and a Stern Layer, to calculate counterion concentration profiles as a function of the ionic concentration of the buffer. We find that for low-concentration buffers ([KCl] from 10 μM to 10 mM) the counterion concentration is largely independent of the buffer's ionic concentration, but for physiological-concentration buffers ([KCl] from 100 to 500 mM) the counterion concentration varies dramatically with changes in the buffer's ionic concentration. We then calculate the conductivity of microtubule-counterion complexes, which are found to be more conductive than the buffer when the buffer's ionic concentrations is less than ≈100 mM and less conductive otherwise. These results demonstrate the importance of accounting for the ionic concentration of the buffer when analyzing microtubule electrical properties both under laboratory and physiological conditions. We conclude by calculating the basic electrical parameters of microtubules over a range of ionic buffer concentrations applicable to nanodevice and medical applications.
| S-EPMC8027483 | biostudies-literature
Project description:Temperature measurements with thin thermocouples embedded in ultrasound fields are strongly subjected to a viscous heating artifact (VHA). The artifact contribution decays over time; therefore, it can be minimized at late temperature readings. However, previous studies have failed to demonstrate a rigorous method for determining the optimal time point at which the artifact contribution is negligible. In this study, we present an iterative processing method based on successive curve fittings using an artifact-independent model. The fitting starting point moves at each iteration until the maximum R2 indicates where the viscous heating is minimum. A solution of the bioheat transfer equation is used to account for blood perfusion, thus enabling in vivo measurements. Three T-type thermocouples with different diameters and sensitivities were assessed in an excised canine liver and in the mouse brain in vivo. We found that the artifact constitutes up to 81% ± 5% of wire thermocouple readings. The best-fit time varied in the liver samples ( n = 3 ) from 0 to 3.335 ± 0.979 s and in the mouse brain ( n = 5 ) from 0 to 0.498 ± 0.457 s at variable experimental conditions, which clearly demonstrates the need of the method for finding the appropriate starting time point of the fit. This study introduces a statistical method to determine the best time to fit a curve that can back-estimate temperature in tissues under ultrasound exposure using thermocouples. This method allows temperature evaluation in vivo and in vitro during a validation and safety assessment of a wide range of therapeutic and diagnostic ultrasound modalities.
| S-EPMC6944748 | biostudies-literature
Project description:The research of adsorption theory has recently gained renewed attention due to its critical relevance to a number of trending industrial applications, hydrogen storage and shale gas exploration for instance. The existing theoretical foundation, laid mostly in the early twentieth century, was largely based on simple heuristic molecular interaction models and static interaction potential which, although being insightful in illuminating the fundamental mechanisms, are insufficient for computations with realistic adsorbent structure and adsorbate hydrodynamics, both critical for real-life applications. Here we present and validate a novel lattice Boltzmann model incorporating both adsorbate-adsorbate and adsorbate-adsorbent interactions with hydrodynamics which, for the first time, allows adsorption to be computed with real-life details. Connection with the classic Ono-Kondo lattice theory is established and various adsorption isotherms, both within and beyond the IUPAC classification are observed as a pseudo-potential is varied. This new approach not only enables an important physical to be simulated for real-life applications, but also provides an enabling theoretical framework within which the fundamentals of adsorption can be studied.
| S-EPMC4891696 | biostudies-other
Project description:Experimental results have demonstrated that the numbers of counterions surrounding nucleic acids differ from those predicted by the nonlinear Poisson-Boltzmann equation, NLPBE. Some studies have fit these data against the ion size in the size-modified Poisson-Boltzmann equation, SMPBE, but the present study demonstrates that other parameters, such as the Stern layer thickness and the molecular surface definition, can change the number of bound ions by amounts comparable to varying the ion size. These parameters will therefore have to be fit simultaneously against experimental data. In addition, the data presented here demonstrate that the derivative, SK, of the electrostatic binding free energy, ΔGel, with respect to the logarithm of the salt concentration is sensitive to these parameters, and experimental measurements of SK could be used to parameterize the model. However, although better values for the Stern layer thickness and ion size and better molecular surface definitions could improve the model's predictions of the numbers of ions around biomolecules and SK, ΔGel itself is more sensitive to parameters, such as the interior dielectric constant, which in turn do not significantly affect the distributions of ions around biomolecules. Therefore, improved estimates of the ion size and Stern layer thickness to use in the SMPBE will not necessarily improve the model's predictions of ΔGel.
| S-EPMC3977794 | biostudies-literature
Project description:We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L(∞) estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme are demonstrated with FETK through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.
| S-EPMC3178275 | biostudies-literature
Project description:BackgroundFitting four-parameter sigmoidal models is one of the methods established in the analysis of quantitative real-time PCR (qPCR) data. We had observed that these models are not optimal in the fitting outcome due to the inherent constraint of symmetry around the point of inflection. Thus, we found it necessary to employ a mathematical algorithm that circumvents this problem and which utilizes an additional parameter for accommodating asymmetrical structures in sigmoidal qPCR data.ResultsThe four-parameter models were compared to their five-parameter counterparts by means of nested F-tests based on the residual variance, thus acquiring a statistical measure for higher performance. For nearly all qPCR data we examined, five-parameter models resulted in a significantly better fit. Furthermore, accuracy and precision for the estimation of efficiencies and calculation of quantitative ratios were assessed with four independent dilution datasets and compared to the most commonly used quantification methods. It could be shown that the five-parameter model exhibits an accuracy and precision more similar to the non-sigmoidal quantification methods.ConclusionThe five-parameter sigmoidal models outperform the established four-parameter model with high statistical significance. The estimation of essential PCR parameters such as PCR efficiency, threshold cycles and initial template fluorescence is more robust and has smaller variance. The model is implemented in the qpcR package for the freely available statistical R environment. The package can be downloaded from the author's homepage.
| S-EPMC2386824 | biostudies-literature
Project description:To investigate the dosimetric impact of the heterogeneity dose calculation Acuros (Transpire Inc., Gig Harbor, WA), a grid-based Boltzmann equation solver (GBBS), for brachytherapy in a cohort of cervical cancer patients.The impact of heterogeneities was retrospectively assessed in treatment plans for 26 patients who had previously received (192)Ir intracavitary brachytherapy for cervical cancer with computed tomography (CT)/magnetic resonance-compatible tandems and unshielded colpostats. The GBBS models sources, patient boundaries, applicators, and tissue heterogeneities. Multiple GBBS calculations were performed with and without solid model applicator, with and without overriding the patient contour to 1 g/cm(3) muscle, and with and without overriding contrast materials to muscle or 2.25 g/cm(3) bone. Impact of source and boundary modeling, applicator, tissue heterogeneities, and sensitivity of CT-to-material mapping of contrast were derived from the multiple calculations. American Association of Physicists in Medicine Task Group 43 (TG-43) guidelines and the GBBS were compared for the following clinical dosimetric parameters: Manchester points A and B, International Commission on Radiation Units and Measurements (ICRU) report 38 rectal and bladder points, three and nine o'clock, and (D2cm3) to the bladder, rectum, and sigmoid.Points A and B, D(2) cm(3) bladder, ICRU bladder, and three and nine o'clock were within 5% of TG-43 for all GBBS calculations. The source and boundary and applicator account for most of the differences between the GBBS and TG-43 guidelines. The D(2cm3) rectum (n = 3), D(2cm3) sigmoid (n = 1), and ICRU rectum (n = 6) had differences of >5% from TG-43 for the worst case incorrect mapping of contrast to bone. Clinical dosimetric parameters were within 5% of TG-43 when rectal and balloon contrast were mapped to bone and radiopaque packing was not overridden.The GBBS has minimal impact on clinical parameters for this cohort of patients with unshielded applicators. The incorrect mapping of rectal and balloon contrast does not have a significant impact on clinical parameters. Rectal parameters may be sensitive to the mapping of radiopaque packing.
| S-EPMC3955059 | biostudies-literature
Project description:Grid-based solvers of the Poisson-Boltzmann, PB, equation are routinely used to estimate electrostatic binding, ??Gel, and solvation, ?Gel, free energies. The accuracies of such estimates are subject to grid discretization errors from the finite difference approximation to the PB equation. Here, we show that the grid discretization errors in ??Gel are more significant than those in ?Gel, and can be divided into two parts: (i) errors associated with the relative positioning of the grid and (ii) systematic errors associated with grid spacing. The systematic error in particular is significant for methods, such as the molecular mechanics PB surface area, MM-PBSA, approach that predict electrostatic binding free energies by averaging over an ensemble of molecular conformations. Although averaging over multiple conformations can control for the error associated with grid placement, it will not eliminate the systematic error, which can only be controlled by reducing grid spacing. The present study indicates that the widely-used grid spacing of 0.5 Å produces unacceptable errors in ??Gel, even though its predictions of ?Gel are adequate for the cases considered here. Although both grid discretization errors generally increase with grid spacing, the relative sizes of these errors differ according to the solute-solvent dielectric boundary definition. The grid discretization errors are generally smaller on the Gaussian surface used in the present study than on either the solvent-excluded or van der Waals surfaces, which both contain more surface discontinuities (e.g., sharp edges and cusps). Additionally, all three molecular surfaces converge to very different estimates of ??Gel.
| S-EPMC3756143 | biostudies-literature