skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A Transfer Free Energy Based Implicit Solvent Model for Protein Simulations in Solvent Mixtures: Urea-Induced Denaturation as a Case Study
Award ID(s):
1716956
PAR ID:
10356899
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
126
Issue:
24
ISSN:
1520-6106
Page Range / eLocation ID:
4472 to 4482
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We report a Monte Carlo simulation study of length-scale dependent density fluctuations in cavities in the coarse-grained mW representation of water at ambient conditions. Specifically, we use a combination of test particle insertion and umbrella sampling techniques to examine the full range of water occupation states in spherical cavities up to 6.3 Å in radius in water. As has previously been observed, water density fluctuations are found to be effectively Gaussian in nature for atomic-scale cavities, but as the cavities get larger they exhibit a non-Gaussian “fat-tail” distribution for lower occupancy states. We introduce a new statistical thermodynamic approach to analyze non-Gaussian fluctuations based on the radial distribution of waters about cavities with varying numbers of waters within its boundaries. It is shown that the onset on these non-Gaussian fluctuations is a result of the formation of a bubble within the cavity as it is emptied that is accompanied by the adsorption of waters onto its interior surface. We revisit a theoretical framework we previously introduced to describe Gaussian fluctuations within cavities to now incorporate bubble formation by including surface tension contributions. This modified theory accurately describes density fluctuations within both atomic and meso-scale cavities. Moreover, the theory predicts the transition from Gaussian to non-Gaussian fluctuations at a specific cavity occupancy in excellent agreement with simulation observations. 
    more » « less
  2. Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters. 
    more » « less
  3. This computational study characterises charge-transfer-to-solvent (CTTS) states of aqueous thiocyanate anion using equation-of-motion coupled-cluster methods combined with electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) scheme. Equilibrium sampling was carried out using classical molecular dynamics (MD) with standard force-fields and QM/MM ab initio molecular dynamics (AIMD) using density functional theory. The two calculations yield significantly different local structure around solvated SCN− . Because of the diffuse character of CTTS states, they are very sensitive to the local structure of solvent around the solute and its dynamic fluctuations. Owing to this sensitivity, the spectra computed using MD and AIMD based snapshots differ considerably. This sensitivity suggests that the spectroscopy exploiting CTTS transitions can provide an experimental handle for assessing the quality of force-fields and density functionals. By combining CTTS-based spectroscopies with reliable theoretical modeling, detailed microscopic information of the solvent structure can be obtained. We present a robust computational protocol for modeling spectra of solvated anions and emphasise the use of an ab initio characterization of individual electronic transitions as CTTS or local excitations. 
    more » « less
  4. Wetting experiments show pure graphene to be weakly hydrophilic, but its contact angle (CA) also reflects the character of the supporting material. Measurements and Molecular Dynamics simulations on suspended and supported graphene often reveal a CA reduction due to the presence of the supporting substrate. A similar reduction is consistently observed when graphene is wetted from both sides. The effect has been attributed to transparency to molecular interactions across the graphene sheet, however, the possibility of substrate-induced graphene polarization has also been considered. Computer simulations of CA on graphene have so far been determined by ignoring the material’s conducting properties. We improve the graphene model by incorporating its conductivity according to the Constant Applied Potential Molecular Dynamics. Using this method, we compare the wettabilities of suspended graphene and graphene supported by water by measuring the CA of cylindrical water drops on the sheets. The inclusion of graphene conductivity and concomitant polarization effects lead to a lower CA on suspended graphene but the CA reduction is significantly bigger when the sheets are also wetted from the opposite side. The stronger adhesion is accompanied by a profound change in the correlations among water molecules across the sheet. While partial charges on water molecules interacting across an insulator sheet attract charges of the opposite sign, apparent attraction among like charges is manifested across the conducting graphene. The change is associated with graphene polarization, as the image charges inside the conductor attract equally signed partial charges of water molecules on both sides of the sheet. Additionally, by using a non-polar liquid (diiodomethane), we affirm a detectable wetting translucency when liquid-liquid forces are dominated by dispersive interactions. Our findings are important for predictive modeling toward a variety of applications including sensors, fuel cell membranes, water filtration, and graphene-based electrode materials in high-performance supercapacitors. 
    more » « less