Variational approach for nonpolar solvation analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Chen Z
PROVIDER: S-EPMC3436914 | biostudies-other | 2012 Aug
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Using precomputed near neighbor or proximal distribution functions (pDFs) that approximate solvent density about atoms in a chemically bonded context one can estimate the solvation structures around complex solutes and the corresponding solute-solvent energetics. In this contribution, we extend this technique to calculate the solvation free energies (ΔG) of a variety of solutes. In particular we use pDFs computed for small peptide molecules to estimate ΔG for larger peptide systems. We separately compute the non polar (ΔGvdW) and electrostatic (ΔGelec) components of the underlying potential model. Here we show how the former can be estimated by thermodynamic integration using pDF-reconstructed solute-solvent interaction energy. The electrostatic component can be approximated with Linear Response theory as half of the electrostatic solute-solvent interaction energy. We test the method by calculating the solvation free energies of butane, propanol, polyalanine, and polyglycine and by comparing with traditional free energy simulations. Results indicate that the pDF-reconstruction algorithm approximately reproduces ΔGvdW calculated by benchmark free energy simulations to within ∼ kcal/mol accuracy. The use of transferable pDFs for each solute atom allows for a rapid estimation of ΔG for arbitrary molecular systems.
| S-EPMC5398929 | biostudies-literature
Project description:We incorporate the Poisson-Boltzmann (PB) theory of electrostatics into our variational implicit-solvent model (VISM) for the solvation of charged molecules in an aqueous solvent. In order to numerically relax the VISM free-energy functional by our level-set method, we develop highly accurate methods for solving the dielectric PB equation and for computing the dielectric boundary force. We also apply our VISM-PB theory to analyze the solvent potentials of mean force and the effect of charges on the hydrophobic hydration for some selected molecular systems. These include some single ions, two charged particles, two charged plates, and the host-guest system Cucurbit[7]uril and Bicyclo[2.2.2]octane. Our computational results show that VISM with PB theory can capture well the sensitive response of capillary evaporation to the charge in hydrophobic confinement and the polymodal hydration behavior and can provide accurate estimates of binding affinity of the host-guest system. We finally discuss several issues for further improvement of VISM.
| S-EPMC3985794 | biostudies-other
Project description:A goal in computational chemistry is computing hydration free energies of nonpolar and charged solutes accurately, but with much greater computational speeds than in today's explicit-water simulations. Here, we take one step in that direction: a simple model of solvating waters that is analytical and thus essentially instantaneous to compute. Each water molecule is a 2-dimensional dipolar hydrogen-bonding disk that interacts around small circular solutes with different nonpolar and charge interactions. The model gives good qualitative agreement with experiments. As a function of the solute radius, it gives the solvation free energy, enthalpy and entropy as a function of temperature for the inert gas series Ne, Ar, Kr, and Xe. For anions and cations, it captures relatively well the trends versus ion radius. This approach should be readily generalizable to three dimensions.
| S-EPMC7958497 | biostudies-literature
Project description:Continuum solvation models provide appealing alternatives to explicit solvent methods because of their ability to reproduce solvation effects while alleviating the need for expensive sampling. Our previous work has demonstrated that Poisson-Boltzmann methods are capable of faithfully reproducing polar explicit solvent forces for dilute protein systems; however, the popular solvent-accessible surface area model was shown to be incapable of accurately describing nonpolar solvation forces at atomic-length scales. Therefore, alternate continuum methods are needed to reproduce nonpolar interactions at the atomic scale. In the present work, we address this issue by supplementing the solvent-accessible surface area model with additional volume and dispersion integral terms suggested by scaled particle models and Weeks-Chandler-Andersen theory, respectively. This more complete nonpolar implicit solvent model shows very good agreement with explicit solvent results and suggests that, although often overlooked, the inclusion of appropriate dispersion and volume terms are essential for an accurate implicit solvent description of atomic-scale nonpolar forces.
| S-EPMC1482494 | biostudies-literature
Project description:Detection-nondetection data are often used to investigate species range dynamics using Bayesian occupancy models which rely on the use of Markov chain Monte Carlo (MCMC) methods to sample from the posterior distribution of the parameters of the model. In this article we develop two Variational Bayes (VB) approximations to the posterior distribution of the parameters of a single-season site occupancy model which uses logistic link functions to model the probability of species occurrence at sites and of species detection probabilities. This task is accomplished through the development of iterative algorithms that do not use MCMC methods. Simulations and small practical examples demonstrate the effectiveness of the proposed technique. We specifically show that (under certain circumstances) the variational distributions can provide accurate approximations to the true posterior distributions of the parameters of the model when the number of visits per site (K) are as low as three and that the accuracy of the approximations improves as K increases. We also show that the methodology can be used to obtain the posterior distribution of the predictive distribution of the proportion of sites occupied (PAO).
| S-EPMC4771718 | biostudies-literature
Project description:Solvation is a fundamental driving force in many biological processes including biomolecular recognition and self-assembly, not to mention protein folding, dynamics, and function. The variational implicit solvent method (VISM) is a theoretical tool currently developed and optimized to estimate solvation free energies for systems of very complex topology, such as biomolecules. VISM's theoretical framework makes it unique because it couples hydrophobic, van der Waals, and electrostatic interactions as a functional of the solvation interface. By minimizing this functional, VISM produces the solvation interface as an output of the theory. In this work, we push VISM to larger scale applications by combining it with coarse-grained solute Hamiltonians adapted from the MARTINI framework, a well-established mesoscale force field for modeling large-scale biomolecule assemblies. We show how MARTINI-VISM (MVISM) compares with atomistic VISM (AVISM) for a small set of proteins differing in size, shape, and charge distribution. We also demonstrate MVISM's suitability to study the solvation properties of an interesting encounter complex, barnase-barstar. The promising results suggest that coarse-graining the protein with the MARTINI force field is indeed a valuable step to broaden VISM's and MARTINI's applications in the near future.
| S-EPMC5740479 | biostudies-literature
Project description:In evolutionary biology, niche construction is sometimes described as a genuine evolutionary process whereby organisms, through their activities and regulatory mechanisms, modify their environment such as to steer their own evolutionary trajectory, and that of other species. There is ongoing debate, however, on the extent to which niche construction ought to be considered a bona fide evolutionary force, on a par with natural selection. Recent formulations of the variational free-energy principle as applied to the life sciences describe the properties of living systems, and their selection in evolution, in terms of variational inference. We argue that niche construction can be described using a variational approach. We propose new arguments to support the niche construction perspective, and to extend the variational approach to niche construction to current perspectives in various scientific fields.
| S-EPMC5938575 | biostudies-literature
Project description:Previous work describes a computational solvation model called semi-explicit assembly (SEA). The SEA water model computes the free energies of solvation of nonpolar and polar solutes in water with good efficiency and accuracy. However, SEA gives systematic errors in the solvation free energies of ions and charged solutes. Here, we describe field-SEA, an improved treatment that gives accurate solvation free energies of charged solutes, including monatomic and polyatomic ions and model dipeptides, as well as nonpolar and polar molecules. Field-SEA is computationally inexpensive for a given solute because explicit-solvent model simulations are relegated to a precomputation step and because it represents solvating waters in terms of a solute's free-energy field. In essence, field-SEA approximates the physics of explicit-model simulations within a computationally efficient framework. A key finding is that an atom's solvation shell inherits characteristics of a neighboring atom, especially strongly charged neighbors. Field-SEA may be useful where there is a need for solvation free-energy computations that are faster than explicit-solvent simulations and more accurate than traditional implicit-solvent simulations for a wide range of solutes.
| S-EPMC4065164 | biostudies-literature
Project description:We propose a novel Riemannian geometric framework for variational inference in Bayesian models based on the nonparametric Fisher-Rao metric on the manifold of probability density functions. Under the square-root density representation, the manifold can be identified with the positive orthant of the unit hypersphere S ∞ in L2 , and the Fisher-Rao metric reduces to the standard L2 metric. Exploiting such a Riemannian structure, we formulate the task of approximating the posterior distribution as a variational problem on the hypersphere based on the α-divergence. This provides a tighter lower bound on the marginal distribution when compared to, and a corresponding upper bound unavailable with, approaches based on the Kullback-Leibler divergence. We propose a novel gradient-based algorithm for the variational problem based on Fréchet derivative operators motivated by the geometry of S ∞, and examine its properties. Through simulations and real data applications, we demonstrate the utility of the proposed geometric framework and algorithm on several Bayesian models.
| S-EPMC7540671 | biostudies-literature
Project description:Testing the urine nonpolar sulfateome can enable discovery of xenobiotics that are most likely to be bioactive. This is based on the fact that nonpolar xenobiotics are more likely to enter cells where they tend to undergo metabolism, in part, to sulfates that are then largely excreted into the urine.The following sequence of steps, with conditions that achieve high reproducibility, was applied to large human urine samples: (1) competitive nonpolar extraction with a porous extraction paddle; (2) weak anion-exchange extraction with strong organic washing; and (3) ultrahigh-performance liquid chromatography (UHPLC)/negative ion matrix-assisted laser desorption ionization tandem time-of-flight mass spectrometery (MALDI-TOF/TOF-MS) with recording of ions with signal-to-noise (S/N) ??20 that yielded M-1-80 (loss of SO3 ) or m/z 97 (HSO4- ) upon fragmentation.From a collection of urine samples from six pregnant women, the masses of 1129 putative sulfates were measured. Three lists of candidate compounds (preliminary hits) from these masses were formed by searching METLIN, especially via MATLAB, yielding putative xenobiotic contaminants (35 compounds), steroids (122), and flavonoids (1582).A new way to reveal some of the nonpolar xenobiotic exposome has been developed that applies to urine samples. The value of the method is to suggest xenobiotics for subsequent targeted analysis in the population of people under study, in order to relate the environment to health and disease. Copyright © 2016 John Wiley & Sons, Ltd.
| S-EPMC5070653 | biostudies-literature