An official website of the United States government
Monte Carlo (MC) methods are important computational tools for molecular structure optimizations and predictions. When solvent effects are explicitly considered, MC methods become very expensive due to the large degree of freedom associated with the water molecules and mobile ions. Alternatively implicit-solvent MC can largely reduce the computational cost by applying a mean field approximation to solvent effects and meanwhile maintains the atomic detail of the target molecule. The two most popular implicit-solvent models are the Poisson-Boltzmann (PB) model and the Generalized Born (GB) model in a way such that the GB model is an approximation to the PB model but is much faster in simulation time. In this work, we develop a machine learning-based implicit-solvent Monte Carlo (MLIMC) method by combining the advantages of both implicit solvent models in accuracy and efficiency. Specifically, the MLIMC method uses a fast and accurate PB-based machine learning (PBML) scheme to compute the electrostatic solvation free energy at each step. We validate our MLIMC method by using a benzene-water system and a protein-water system. We show that the proposed MLIMC method has great advantages in speed and accuracy for molecular structure optimization and prediction.
more »
« less
A generalized Kirkwood implicit solvent for the polarizable AMOEBA protein model
Computational simulation of biomolecules can provide important insights into protein design, protein-ligand binding interactions, and ab initio biomolecular folding, among other applications. Accurate treatment of the solvent environment is essential in such applications, but the use of explicit solvents can add considerable cost. Implicit treatment of solvent effects using a dielectric continuum model is an attractive alternative to explicit solvation since it is able to describe solvation effects without the inclusion of solvent degrees of freedom. Previously, we described the development and parameterization of implicit solvent models for small molecules. Here, we extend the parameterization of the generalized Kirkwood (GK) implicit solvent model for use with biomolecules described by the AMOEBA force field via the addition of corrections to the calculation of effective radii that account for interstitial spaces that arise within biomolecules. These include element-specific pairwise descreening scale factors, a short-range neck contribution to describe the solvent-excluded space between pairs of nearby atoms, and finally tanh-based rescaling of the overall descreening integral. We then apply the AMOEBA/GK implicit solvent to a set of ten proteins and achieve an average coordinate root mean square deviation for the experimental structures of 2.0 Å across 500 ns simulations. Overall, the continued development of implicit solvent models will help facilitate the simulation of biomolecules on mechanistically relevant timescales.
more »
« less
- Award ID(s):
- 1751688
- PAR ID:
- 10437302
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 159
- Issue:
- 5
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Solvation effects profoundly influence the characteristics and behavior of chemical systems in liquid solutions. The interaction between solute and solvent molecules intricately impacts solubility, reactivity, stability, and various chemical processes. Continuum solvation models gained prominence in quantum chemistry by implicitly capturing these interactions and enabling efficient investigations of diverse chemical systems in solution. In comparison, continuum solvation models in condensed matter simulation are very recent. Among these, the self‐consistent continuum solvation (SCCS) and the soft‐sphere continuum solvation models (SSCS) have been among the first to be successfully parameterized and extended to model periodic systems in aqueous solutions and electrolytes. As most continuum approaches, these models depend on a number of parameters that are linked to experimental or theoretical properties of the solvent, or that can be tuned based on reference data. Here, we present a systematic parameterization of the SSCS model for over 100 nonaqueous solvents. We validate the model's efficacy across diverse solvent environments by leveraging experimental solvation‐free energies and partition coefficients from comprehensive databases. The average root means square error over all the solvents was calculated as 0.85 kcal/mol which is below the chemical accuracy (1 kcal/mol). Similarly to what has been reported by Hille et al. (J. Chem. Phys.2019,150, 041710.) for the SCCS model, a single‐parameter model accurately reproduces experimental solvation energies, showcasing the transferability and predictive power of these continuum approaches. Our findings underscore the potential for a unified approach to predict solvation properties, paving the way for enhanced computational studies across various chemical environments.more » « less
-
Nanoporous carbons play an important role in different electrochemical applications such as being utilized as electrodes in supercapacitors. Application of electric potential to a porous electrode in electrolyte solution stimulates adsorption or desorption of ions on the electrode surface. Electrosorption causes appearance of solvation pressure in the pores and results in electrode deformation. In this work, using molecular dynamics simulations and the continuum theory based on the modified Poisson-Boltzmann equation, we studied the structure of the electrical double layer in slit graphitic micropores filled with a NaCl aqueous solution, and solvation pressure in these pores. We focused on the behavior of the solvation pressure as a function of pore width and surface charge density. Within molecular dynamics simulations, two different water models were used -- an explicit model based on SPC/E water molecules and an implicit model, i.e., structureless background with fixed dielectric permittivity. The latter allows us to relate molecular dynamics simulations to the continuum theory. Simulations with explicit water show a qualitatively different behavior of the solvation pressure in the 1 and 2 nm pores as a function of the surface charge density. We demonstrated that the value of the solvation pressure is defined by a delicate balance between Van der Waals and electrostatic contributions. We demonstrated that the theory predicts the dependence of the solvation pressure on the pore width, which matches the results of simulations using the implicit water model. Finally, we adapted the continuum theory, developed for adsorption-induced deformation to estimate the deformation of a carbon electrode due to electrosorption. Our results can be used in the further development of nanoporous actuators working based on electrosorption-induced deformation.more » « less
-
n→π* interactions between consecutive carbonyls stabilize the α-helix and polyproline II helix (PPII) conformations in proteins. n→π* interactions have been suggested to provide significant conformational biases to the disordered states of proteins. To understand the roles of solvation on the strength of n→π* interactions, computational investigations were conducted on a model n→π* interaction, the twisted-parallel-offset formaldehyde dimer, as a function of explicit solvation of the donor and acceptor carbonyls, using water and HF. In addition, the effects of urea, thiourea, guanidinium, and monovalent cations on n→π* interaction strength were examined. Solvation of the acceptor carbonyl significantly strengthens the n→π* interaction, while solvation of the donor carbonyl only modestly weakens the n→π* interaction. The n→π* interaction strength was maximized with two solvent molecules on the acceptor carbonyl. Urea stabilized the n→π* interaction via simultaneous engagement of both oxygen lone pairs on the acceptor carbonyl. Solvent effects were further investigated in the model peptides Ac-Pro-NMe 2 , Ac-Ala-NMe 2 , and Ac-Pro 2 -NMe 2 . Solvent effects in peptides were similar to those in the formaldehyde dimer, with solvation of the acceptor carbonyl increasing n→π* interaction strength and resulting in more compact conformations, in both the proline endo and exo ring puckers, as well as a reduction in the energy difference between these ring puckers. Carbonyl solvation leads to an energetic preference for PPII over both the α-helix and β/extended conformations, consistent with experimental data that protic solvents and protein denaturants both promote PPII. Solvation of the acceptor carbonyl weakens the intraresidue C5 hydrogen bond that stabilizes the β conformation.more » « less
-
Abstract Although cyclic voltammetry (CV) measurements in solution have been widely used to determine the highest occupied molecular orbital energy (EHOMO) of semiconducting organic molecules, an understanding of the experimentally observed discrepancies due to the solvent used is lacking. To explain these differences, we investigate the solvent effects onEHOMOby combining density functional theory and molecular dynamics calculations for four donor molecules with a common backbone moiety. We compare the experimentalEHOMOvalues to the calculated values obtained from either implicit or first solvation shell theories. We find that the first solvation shell method can capture theEHOMOvariation arising from the functional groups in solution, unlike the implicit method. We further applied the first solvation shell method to other semiconducting organic molecules measured in solutions for different solvents. We find that theEHOMOobtained using an implicit method is insensitive to solvent choice. The first solvation shell, however, producesEHOMOvalues that are sensitive to solvent choices and agrees with published experimental results. The solvent sensitivity arises from a hierarchy of three effects: (1) the solute electronic state within a surrounding dielectric continuum, (2) ambient temperature or solvent atoms changing the solute geometry, and (3) electronic interactions between the solute and solvents. The implicit method, on the other hand, only captures the effect of a dielectric environment. Our findings suggest thatEHOMOobtained by CV measurements should account for the influence of solvent when the results are reported, interpreted, or compared to other molecules.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1063/5.0158914
