Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCoupling nonpolar and polar solvation free energies in implicit solvent models
76
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Recent studies on the solvation of atomistic and nanoscale solutes indicate that a strong coupling exists between the hydrophobic, dispersion, and electrostatic contributions to the solvation free energy, a facet not considered in current implicit solvent models. We suggest a theoretical formalism which accounts for coupling by minimizing the Gibbs free energy of the solvent with respect to a solvent volume exclusion function. The resulting differential equation is similar to the Laplace-Young equation for the geometrical description of capillary interfaces, but is extended to microscopic scales by explicitly considering curvature corrections as well as dispersion and electrostatic contributions. Unlike existing implicit solvent approaches, the solvent accessible surface is an output of our model. The presented formalism is illustrated on spherically or cylindrically symmetrical systems of neutral or charged solutes on different length scales. The results are in agreement with computer simulations and, most importantly, demonstrate that our method captures the strong sensitivity of solvent expulsion and dewetting to the particular form of the solvent-solute interactions.
rate research
Read More
In this work, a systematic protocol is proposed to automatically parametrize implicit solvent models with polar and nonpolar components. The proposed protocol utilizes the classical Poisson model or the Kohn-Sham density functional theory (KSDFT) based polarizable Poisson model for modeling polar solvation free energies. For the nonpolar component, either the standard model of surface area, molecular volume, and van der Waals interactions, or a model with atomic surface areas and molecular volume is employed. Based on the assumption that similar molecules have similar parametrizations, we develop scoring and ranking algorithms to classify solute molecules. Four sets of radius parameters are combined with four sets of charge force fields to arrive at a total of 16 different parametrizations for the Poisson model. A large database with 668 experimental data is utilized to validate the proposed protocol. The lowest leave-one-out root mean square (RMS) error for the database is 1.33k cal/mol. Additionally, five subsets of the database, i.e., SAMPL0-SAMPL4, are employed to further demonstrate that the proposed protocol offers some of the best solvation predictions. The optimal RMS errors are 0.93, 2.82, 1.90, 0.78, and 1.03 kcal/mol, respectively for SAMPL0, SAMPL1, SAMPL2, SAMPL3, and SAMPL4 test sets. These results are some of the best, to our best knowledge.
Continuum models to handle solvent and electrolyte effects in an effective way have a long tradition in quantum-chemistry simulations and are nowadays also being introduced in computational condensed-matter and materials simulations. A key ingredient of continuum models is the choice of the solute cavity, i.e. the definition of the sharp or smooth boundary between the regions of space occupied by the quantum-mechanical (QM) system and the continuum embedding environment. Although most of the solute-based approaches developed lead to models with comparable and high accuracy when applied to small organic molecules, they can introduce significant artifacts when complex systems are considered. As an example, condensed-matter simulations often deal with supports that present open structures. Similarly, unphysical pockets of continuum solvent may appear in systems featuring multiple molecular components. Here, we introduce a solvent-aware approach to eliminate the unphysical effects where regions of space smaller than the size of a single solvent molecule could still be filled with a continuum environment. We do this by defining a smoothly varying solute cavity that overcomes several of the limitations of straightforward solute-based definitions. This new approach applies to any smooth local definition of the continuum interface, being it based on the electronic density or the atomic positions of the QM system. It produces boundaries that are continuously differentiable with respect to the QM degrees of freedom, leading to accurate forces and/or Kohn-Sham potentials. Benchmarks on semiconductor substrates and on explicit water substrates confirm the flexibility and the accuracy of the approach and provide a general set of parameters for condensed-matter systems featuring open structures and/or explicit liquid components.
We demonstrate that with two small modifications, the popular dielectric continuum model is capable of predicting, with high accuracy, ion solvation thermodynamics in numerous polar solvents, and ion solvation free energies in water--co-solvent mixtures. The first modification involves perturbing the macroscopic dielectric-flux interface condition at the solute--solvent interface with a nonlinear function of the local electric field, giving what we have called a solvation-layer interface condition (SLIC). The second modification is a simple treatment of the microscopic interface potential (static potential). We show that the resulting model exhibits high accuracy without the need for fitting solute atom radii in a state-dependent fashion. Compared to experimental results in nine water--co-solvent mixtures, SLIC predicts transfer free energies to within 2.5 kJ/mol. The co-solvents include both protic and aprotic species, as well as biologically relevant denaturants such as urea and dimethylformamide. Furthermore, our results indicate that the interface potential is essential to reproduce entropies and heat capacities. The present work, together with previous studies of SLIC illustrating its accuracy for biomolecules in water, indicates it as a promising dielectric continuum model for accurate predictions of molecular solvation in a wide range of conditions.
We consider the phase diagram of self-avoiding walks (SAW) on the simple cubic lattice subject to surface and bulk interactions, modeling an adsorbing surface and variable solvent quality for a polymer in dilute solution, respectively. We simulate SAWs at specific interaction strengths to focus on locating certain transitions and their critical behavior. By collating these new results with previous results we sketch the complete phase diagram and show how the adsorption transition is affected by changing the bulk interaction strength. This expands on recent work considering how adsorption is affected by solvent quality. We demonstrate that changes in the adsorption crossover exponent coincide with phase boundaries.
The Jarzynski equality and the fluctuation theorem relate equilibrium free energy differences to non-equilibrium measurements of the work. These relations extend to single-molecule experiments that have probed the finite-time thermodynamics of proteins and nucleic acids. The effects of experimental error and instrument noise have not previously been considered. Here, we present a Bayesian formalism for estimating free-energy changes from non-equilibrium work measurements that compensates for instrument noise and combines data from multiple driving protocols. We reanalyze a recent set of experiments in which a single RNA hairpin is unfolded and refolded using optical tweezers at three different rates. Interestingly, the fastest and farthest-from-equilibrium measurements contain the least instrumental noise, and therefore provide a more accurate estimate of the free energies than a few slow, more noisy, near-equilibrium measurements. The methods we propose here will extend the scope of single-molecule experiments; they can be used in the analysis of data from measurements with AFM, optical, and magnetic tweezers.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
