TIP4P-Ew Water Research Articles

Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP’s conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP’s coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

Influence of Ionic Strength on Hydrophobic Interactions in Water: Dependence on Solute Size and Shape.

Hydrophobicity is a phenomenon of great importance in biology, chemistry, and biochemistry. It is defined as the interaction between nonpolar molecules or groups in water and their low solubility. Hydrophobic interactions affect many processes in water, for example, complexation, surfactant aggregation, and coagulation. These interactions play a pivotal role in the formation and stability of proteins or biological membranes. In the present study, we assessed the effect of ionic strength, solute size, and shape on hydrophobic interactions between pairs of nonpolar particles. Pairs of methane, neopentane, adamantane, fullerene, ethane, propane, butane, hexane, octane, and decane were simulated by molecular dynamics in AMBER 16.0 force field. As a solvent, TIP3P and TIP4PEW water models were used. Potential of mean force (PMF) plots of these dimers were determined at four values of ionic strength, 0, 0.04, 0.08, and 0.40 mol/dm3, to observe its impact on hydrophobic interactions. The characteristic shape of PMFs with three extrema (contact minimum, solvent-separated minimum, and desolvation maximum) was observed for most of the compounds for hydrophobic interactions. Ionic strength affected hydrophobic interactions. We observed a tendency to deepen contact minima with an increase in ionic strength value in the case of spherical and spheroidal molecules. Additionally, two-dimensional distribution functions describing water density and average number of hydrogen bonds between water molecules were calculated in both water models for adamantane and hexane. It was observed that the density of water did not significantly change with the increase in ionic strength, but the average number of hydrogen bonds changed. The latter tendency strongly depends on the water model used for simulations.

Electron Paramagnetic Resonance Measurements of Four Nitroxide Probes in Supercooled Water Explained by Molecular Dynamics Simulations.

Electron paramagnetic resonance (EPR) measurements of the rotational diffusion of small nitroxide probes have been demonstrated to be a powerful technique for experimentally investigating the properties of supercooled liquids, such as water. However, since only the rotational diffusion of the probe molecules is measured and EPR measurements are indirect, it is not clear what the relationship between the behavior of water and the probe molecule is. To address this, we have performed molecular dynamics simulations of four nitroxide probes in TIP4P-Ew and OPC water models to directly compare with EPR experiments and to determine the behavior of the water and the underlying microscopic coupling between the water and the probes. In all, 200 ns simulations were run for 23 temperatures between 253 and 283 K for all four probes with each water model for an aggregate of 36.8 μs of simulation time. Simulations for both water models systematically underestimated the rotational diffusion coefficients for both water and probes, though OPC simulations were generally in better agreement with the experiments than TIP4P-Ew simulations. Despite this, when the temperature dependence of the data was fit to a power law, fit parameters for TIP4P-Ew were generally in better agreement with the experiments than OPC. For probe molecules, the singular temperature was found to be T0 = 226.5 ± 0.4 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 215 ± 2 K for TIP4P-Ew water. While for water molecules, the singular temperature was found to be T0 = 220.3 ± 0.2 K from experiments, T0 = 208 ± 2 K for OPC water, and T0 = 220 ± 1 K for TIP4P-Ew water. Systematic underestimation of the rotational diffusion coefficients was most pronounced at lower temperatures and was clearly observed in changes to the Arrhenius activation energy. Above the maximum density temperature of Tρmax = 277 K, an activation energy of EA ≈ 16.7 kJ/mol was observed for the probes from experiments, while OPC had EA ≈ 15.2 kJ/mol and TIP4P-Ew had EA ≈ 14.6 kJ/mol. Below the maximum density temperature, the activation energy jumped to EA ≈ 32.5 kJ/mol for experiments but only EA ≈ 23 kJ/mol for OPC and EA ≈ 22 kJ/mol for TIP4P-Ew. In all cases, we saw good agreement between the behavior of the probe molecules and water. To understand why, we calculated the average number of hydrogen bonds between the probe molecules and water. From this, we were able to explain the rotational diffusion times for all of the probes. These results show that current molecular models are sufficient to capture physical phenomena observed with EPR and to help elucidate why the probes provide accurate insights into the behavior of supercooled water.

2,4-SUBSTITUTED QUINAZOLINE AS JAK2 INHIBITOR: DOCKING AND MOLECULAR DYNAMICS STUDY

Objective: The involvement of Janus kinase2/signal transducer and activator of transcription (JAK2/STAT3) pathway reported in various solid tumors made authors study the conformational changes of JAK2-3e complex which was previously reported with a moderate percentage of In-vitro JAK2 inhibition.
 Methods: In this present study Compound 3e was reported with a moderate percentage of inhibition of JAK2 protein selected for performing molecular docking and molecular dynamics studies to elucidate the conformational changes with JAK2-3e complex. Docking studies were performed using ChemSketch to draw the structure of the compound and optimized/energy minimized using the Ligprep module of Schrodinger suite, employing optimized potentials for liquid simulations (OPLS-2005) force field. Molecular dynamics simulations were performed for 10 ns for complex using TIP4PEW water solvent model and neutralized by adding sodium ions.
 Results: Docking studies of Compound 3e which has been reported as one of the effective cytotoxic agents and a moderate percentage of In-vitro JAK2 inhibition among the series, showed H-bond interaction with leucine 855, serine936, aspartine994. Dock score and Ligand binding energy with protein suggested compound 3e has shown-4.049,-66.003 kcal/mol respectively. Molecular dynamics simulations elucidated the mechanistic insight of JAK-2 inhibition. The Root means square deviation (RMSD) pattern of both protein and ligands in the JAK2-3e complex observed to be different over 10 ns simulation. In the JAK2-3e complex, an exponential increase in RMSD of Cα and side-chain amino acids is observed during the first 1-3 ns simulation and is stabilized till 10 ns. During the 10 ns simulation, ligand 3e seems to be stable in the complex with an overall deviation<1 Å, despite a drastic increase between 1-3 ns. The ligand RMSD plot suggests that the ligand 3e remained intact within the binding site of the protein and longer time period simulation may elucidate the binding pattern and fate of ligand 3e.
 Conclusion: Results from molecular dynamics simulations elucidated the mechanistic insight of JAK-2 inhibition by 2, 4 disubstituted quinazoline compound that is N’(2-(4-nitrophenyl)quinazoline-4-yl) isonicotinohydrazide) and their binding phenomenon. Molecular docking studies further supported the elucidation of binding patterns of the molecules in the JAK-2 protein environment. Further simulations with a longer time period may provide deeper insights into ligand interactions in the protein environment. It is noteworthy to use compound 3e as a new scaffold for further development of multifunctional compounds.

On the composition dependence of the microscopic structure, thermodynamic, dynamic and dielectric properties of water-dimethyl formamide model mixtures. Molecular dynamics simulation results

Isothermal-isobaric molecular dynamics simulations have been performed to examine an ample set of properties of the model water-N,N-dimethylformamide (DMF) mixture as a function of composition. The SPC-E and TIP4P-Ew water models together with two united atom models for DMF [Chalaris M., Samios J., J. Chem. Phys., 2000, 112, 8581; Cordeiro J., Int. J. Quantum Chem., 1997, 65, 709] were used. Our principal analyses concern the behaviour of structural properties in terms of radial distribution functions, and the number of hydrogen bonds between molecules of different species as well as thermodynamic properties. Namely, we explore the density, excess mixing molar volume and enthalpy, the heat capacity and excess mixing heat capacity. Finally, the self-diffusion coefficients of species and the dielectric constant of the system are discussed. In addition, surface tension of water-DMF mixtures has been calculated and analyzed.

Reduced Point Charge Models of Proteins: Effect of Protein–Water Interactions in Molecular Dynamics Simulations of Ubiquitin Systems

We investigate the influence of various solvent models on the structural stability and protein-water interface of three ubiquitin complexes (PDB access codes: 1Q0W , 2MBB , 2G3Q ) modeled using the Amber99sb force field (FF) and two different point charge distributions. A previously developed reduced point charge model (RPCM), wherein each amino acid residue is described by a limited number of point charges, is tested and compared to its all-atom (AA) version. The complexes are solvated in TIP4P-Ew or TIP3P type water molecules, involving either the scaling of the Lennard-Jones protein-Owater interaction parameters, or the coarse-grain (CG) SIRAH water description. The best agreements between the RPCM and AA models were obtained for structural, protein-water, and ligand-ubiquitin properties when using the TIP4P-Ew water FF with a scaling factor γ of 0.7. At the RPCM level, a decrease in γ, or the inclusion of SIRAH particles, allows weakening of the protein-water interactions. It results in a slight collapse of the protein structure and a less compact hydration shell and, thus, in a decrease in the number of protein-water and water-water H-bonds. The dynamics of the surface protein atoms and of the water shell molecules are also slightly refrained, which allow the generation of stable RPCM trajectories.

Isobaric-isothermal molecular dynamics computer simulations of the properties of water-1,2-dimethoxyethane model mixtures

Isothermal-isobaric molecular dynamics simulations have been performed to examine a broad set of properties of the model water-1,2-dimethoxyethane (DME) mixture as a function of composition. The SPC-E and TIP4P-Ew water models and the modified TraPPE model for DME were applied. Our principal focus was to explore the trends of behaviour of the structural properties in terms of the radial distribution functions, coordination numbers and number of hydrogen bonds between molecules of different species, and of conformations of DME molecules. Thermodynamic properties, such as density, molar volume, enthalpy of mixing and heat capacity at constant pressure have been examined. Finally, the self-diffusion coefficients of species and the dielectric constant of the system were calculated and analyzed.

Comparison of simulation and experimental results for a model aqueous tert-butanol solution.

Molecular dynamics simulations are used to investigate the behavior of aqueous tert-butanol (TBA) solutions for a range of temperatures, using the CHARMM generalized force field (CGenFF) to model TBA and the TIP4P/2005 or TIP4P-Ew water model. Simulation results for the density, isothermal compressibility, constant pressure heat capacity, and self-diffusion coefficients are in good accord with experimental measurements. Agreement with the experiment is particularly good at low TBA concentration, where experiments have revealed anomalies in a number of thermodynamic properties. Importantly, the CGenFF model does not exhibit liquid-liquid demixing at temperatures between 290 and 320 K (for systems of 32 000 molecules), in contrast with the situation for several other common TBA models [R. Gupta and G. N. Patey, J. Chem. Phys. 137, 034509 (2012)]. However, whereas real water and TBA are miscible at all temperatures where the liquid is stable, we observe some evidence of demixing at 340 K and above. To evaluate the structural properties at low concentrations, we compare with both neutron scattering and recent spectroscopic measurements. This reveals that while the CGenFF model is a definite improvement over other models that have been considered, the TBA molecules still exhibit a tendency to associate at low concentrations that is somewhat stronger than that indicated by experiments. Finally, we discuss the range and decay times of the long-range correlations, providing an indication of the system size and simulation times that are necessary in order to obtain reliable results for certain properties.

Reparametrization of Protein Force Field Nonbonded Interactions Guided by Osmotic Coefficient Measurements from Molecular Dynamics Simulations.

There is a small, but growing, body of literature describing the use of osmotic coefficient measurements to validate and reparametrize simulation force fields. Here we have investigated the ability of five very commonly used force field and water model combinations to reproduce the osmotic coefficients of seven neutral amino acids and five small molecules. The force fields tested include AMBER ff99SB-ILDN, CHARMM36, GROMOS54a7, and OPLS-AA, with the first of these tested in conjunction with the TIP3P and TIP4P-Ew water models. In general, for both the amino acids and the small molecules, the tested force fields produce computed osmotic coefficients that are lower than experiment; this is indicative of excessively favorable solute-solute interactions. The sole exception to this general trend is provided by GROMOS54a7 when applied to amino acids: in this case, the computed osmotic coefficients are consistently too high. Importantly, we show that all of the force fields tested can be made to accurately reproduce the experimental osmotic coefficients of the amino acids when minor modifications-some previously reported by others and some that are new to this study-are made to the van der Waals interactions of the charged terminal groups. Special care is required, however, when simulating Proline with a number of the force fields, and a hydroxyl-group specific modification is required in order to correct Serine and Threonine when simulated with AMBER ff99SB-ILDN. Interestingly, an alternative parametrization of the van der Waals interactions in the latter force field, proposed by the Nerenberg and Head-Gordon groups, is shown to immediately produce osmotic coefficients that are in excellent agreement with experiment. Overall, this study reinforces the idea that osmotic coefficient measurements can be used to identify general shortcomings in commonly used force fields' descriptions of solute-solute interactions and further demonstrates that modifications to van der Waals parameters provide a simple route to optimizing agreement with experiment.

Sensitivity of Protein Glass Transition to the Choice of Water Model.

The sensitivity of protein glass transition temperature and dynamics of Trp-cage miniprotein to the choice of water models have been investigated by performing molecular dynamics simulations at different temperatures between 130 and 300 K using mTIP3P, TIP4P, and TIP4P-Ew water models with the CHARMM22 force field. Analysis of the structure of the protein and the ordering of the hydration shell water molecules revealed correlated dynamic crossovers of the protein and its hydration water. It is demonstrated that the parametrization of different water models is robust enough to probe the onset temperatures for the two dynamic crossovers of the peptide: the first one around 150 K which is attributed to active methyl group rotations and the second one around 200-230 K specifically pertaining to glass transition of the protein. The behavior of the protein has been found to be qualitatively similar to the three water models employed, while the structural ordering and number of hydration shell waters vary significantly among the 3-site and 4-site water models. The changes in the probability density map of methyl hydrogens of the protein with temperature also show sensitivity to the choice of water model.

Molecular Dynamics Simulation Study of Parallel Telomeric DNA Quadruplexes at Different Ionic Strengths: Evaluation of Water and Ion Models.

Most molecular dynamics (MD) simulations of DNA quadruplexes have been performed under minimal salt conditions using the Åqvist potential parameters for the cation with the TIP3P water model. Recently, this combination of parameters has been reported to be problematic for the stability of quadruplex DNA, especially caused by the ion interactions inside or near the quadruplex channel. Here, we verify how the choice of ion parameters and water model can affect the quadruplex structural stability and the interactions with the ions outside the channel. We have performed a series of MD simulations of the human full-parallel telomeric quadruplex by neutralizing its negative charge with K(+) ions. Three combinations of different cation potential parameters and water models have been used: (a) Åqvist ion parameters, TIP3P water model; (b) Joung and Cheatham ion parameters, TIP3P water model; and (c) Joung and Cheatham ion parameters, TIP4Pew water model. For the combinations (b) and (c), the effect of the ionic strength has been evaluated by adding increasing amounts of KCl salt (50, 100, and 200 mM). Two independent simulations using the Åqvist parameters with the TIP3P model show that this combination is clearly less suited for the studied quadruplex with K(+) as counterions. In both simulations, one ion escapes from the channel, followed by significant deformation of the structure, leading to deviating conformation compared to that in the reference crystallographic data. For the other combinations of ion and water potentials, no tendency is observed for the channel ions to escape from the quadruplex channel. In addition, the internal mobility of the three loops, torsion angles, and counterion affinity have been investigated at varied salt concentrations. In summary, the selection of ion and water models is crucial as it can affect both the structure and dynamics as well as the interactions of the quadruplex with its counterions. The results obtained with the TIP4Pew model are found to be closest to the experimental data at all of the studied ion concentrations.

Osmotic Pressure Simulations of Amino Acids and Peptides Highlight Potential Routes to Protein Force Field Parameterization.

Recent molecular dynamics (MD) simulations of proteins have suggested that common force fields overestimate the strength of amino acid interactions in aqueous solution. In an attempt to determine the causes of these effects, we have measured the osmotic coefficients of a number of amino acids using the AMBER ff99SB-ILDN force field with two popular water models, and compared the results with available experimental data. With TIP4P-Ew water, interactions between aliphatic residues agree well with experiment, but interactions of the polar residues serine and threonine are found to be excessively attractive. For all tested amino acids, the osmotic coefficients are lower when the TIP3P water model is used. Additional simulations performed on charged amino acids indicate that the osmotic coefficients are strongly dependent on the parameters assigned to the salt ions, with a reparameterization of the sodium/carboxylate interaction reported by the Aksimentiev group significantly improving description of the osmotic coefficient for glutamate. For five neutral amino acids, we also demonstrate a decrease in solute-solute attractions using the recently reported TIP4P-D water model and using the KBFF force field. Finally, we show that for four two-residue peptides improved agreement with experiment can be achieved by rederiving the partial charges for each peptide.

Comparison of hydration behavior and conformational preferences of the Trp-cage mini-protein in different rigid-body water models.

The secondary structure conformational properties and hydration shell metrics of the Trp-cage mini-protein are examined in the folded and unfolded ensembles in mTIP3P, TIP4P, and TIP4P-Ew water models with the CHARMM22 force-field using molecular dynamics simulations at 250 K. Upon changing the water model, the conformational order metrics of the peptide show significant differences in the unfolded rather than in the folded ensemble. The unfolding temperatures for Trp-cage are observed to be around 460, 470, and 430 K in mTIP3P, TIP4P, and TIP4P-Ew, respectively. Upon comparing the results with a previous study on a 16-residue β-hairpin fragment of the 2GB1 protein, the same set of conformational order metrics are found to be insufficient in describing the free energy landscape of peptides having a distinct native secondary structure. However, the hydration shell properties of the peptide have been found to be independent of the sequence of the peptide and it changes in conformation upon unfolding. Our calculations reveal that for a particular water model, the secondary structure preferences in the unfolded ensembles of the two peptides are qualitatively different. The unfolded structures of Trp-cage prefer extended and compact structures in TIP4P-Ew and mTIP3P water, respectively, whereas the β-hairpin peptide prefers extended unfolded structures in mTIP3P. The conformational preferences of the unfolded peptide in a given water model have been found to depend on the peptide sequence, where the binding energies of the water molecules around the polar residues in the unfolded conformations show sensitivity to the multipole moments of the water models. The significance of an accurate description of peptide-solvent interactions in the parametrization of biomolecular force-fields, to obtain an accurate description of conformational preferences, in particular in the unfolded ensembles of proteins, is highlighted.

Cellulose Nanofibrils and Mechanism of their Mineralization in Biomimetic Synthesis of Hydroxyapatite/Native Bacterial Cellulose Nanocomposites: Molecular Dynamics Simulations.

Molecular dynamics (MD) simulation of a nanofibril of native bacterial cellulose (BC) in solutions of mineral ions is presented. The supersaturated calcium-phosphate (CP) solution with the ionic composition of hydroxyapatite and CaCl2 solutions with the concentrations below, equal to, and above the solubility limits are simulated. The influence of solvation models (TIP3P and TIP4P-ew water models) on structural characteristics of the simulated nanofibril and on the crystal nucleation process is assessed. The structural characteristics of cellulose nanofibrils (in particular, of the surface layer) are found to be nearly independent of the solvation models used in the simulation and on the presence of ions in the solutions. It is shown that ionic clusters are formed in the solution rather than on the fibril surface. The cluster sizes are slightly different for the two water models. The effect of the ion-ion interaction parameters on the results is discussed. The main conclusion is that the activity of hydroxyl groups on the BC fibril surface is not high enough to cause adsorption of Ca(2+) ions from the solution. Therefore, the nucleation of CP crystals takes place initially in solution, and then the crystallites formed can be adsorbed on BC nanofibril surfaces.

Cage occupancy of methane hydrates from Gibbs ensemble Monte Carlo simulations

Isobaric-isothermal Gibbs ensemble Monte Carlo simulations of sI methane hydrates in equilibrium with bulk methane are performed to calculate large and small cage occupancies. The OPLS united-atom Lennard-Jones potential, and a variation, are used to model methane, while the TIP4P/Ice model represents water. This model system produces systematically higher total cage occupancies than those of simulations by Henley and Lucia [1] using TIP4P-Ew water and TraPPE-UA methane. We also see higher total occupancies than those observed in prior Raman spectrometry studies by Uchida et al. [2], but find good agreement for large cage occupancies. The simulations provide comprehensive predictions of large and small cage occupancies for temperatures and pressures within the ranges 270–290 K and 20–400 bar, respectively. These predictions demonstrate the advantages of the Gibbs ensemble for the simulation of clathrate hydrates; the ability to control the temperature and pressure, while not constraining the chemical composition of the hydrate phase is the key advantage of this ensemble.

Composition dependence of thermodynamic, dynamic and dielectric properties of water–methanol model mixtures. Molecular dynamics simulation results with the OPLS-AA model for methanol

We have studied thermodynamic and dynamic properties, as well as the dielectric constant, of water–methanol model mixtures in the entire range of composition, by using constant pressure molecular dynamics simulations at room temperature and ambient pressure. The SPC-E and TIP4P-Ew water models are explored in combination with the OPLS-AA (all atom) model of methanol. Changes of the average number of hydrogen bonds between particles of different species and of the fractions of differently bonded molecules are put in correspondence with the evolution of excess mixing properties, of self-diffusion coefficients, reorientational relaxation times and power spectra. We were looking for the existence of correlations between properties of a different kind, in order to find reasons of satisfactory or dissatisfactory performance of the models applied. It seems that any improvement of the description of water–methanol mixtures requires a more sophisticated, carefully parametrized model for methanol and possibly, certain deviations in the combination rules.

Direct Coexistence Methods to Determine the Solubility of Salts in Water from Numerical Simulations. Test Case NaCl.

The solubility of NaCl, an equilibrium between a saturated solution of ions and a solid with a crystalline structure, was obtained from molecular dynamics simulations using the SPC/E and TIP4P-Ew water models. Four initial setups on supersaturated systems were tested on sodium chloride (NaCl) solutions to determine the equilibrium conditions and computational performance: (1) an ionic solution confined between two crystal plates of periodic NaCl, (2) a solution with all the ions initially distributed randomly, (3) a nanocrystal immersed in pure water, and (4) a nanocrystal immersed in an ionic solution. In some cases, the equilibration of the system can take several microseconds. The results from this work showed that the solubility of NaCl was the same, within simulation error, for the four setups, and in agreement with previously reported values from simulations with the setup (1). The system of a nanocrystal immersed in supersaturated solution was found to equilibrate faster than others. In agreement with laser-Doppler droplet measurements, at equilibrium with the solution the crystals in all the setups had a slight positive charge.

Temperature dependence of the microscopic structure and density anomaly of the SPC/E and TIP4P-Ew water models. Molecular dynamics simulation results

We have investigated temperature trends of the microscopic structure of the SPC/E and TIP4P-Ew water models in terms of the pair distribution functions, coordination numbers, the average number of hydrogen bonds, the distribution of bonding states of a single molecule as well as the angular distribution of molecules by using the constant pressure molecular dynamics simulations. The evolution of the structure is put in correspondence with the dependence of water density on high temperatures down to the region of temperatures where the system becomes supercooled. It is shown that the fraction of molecules with three and four bonds determine the maximum density for both models. Moreover, the temperature dependence of the dielectric constant is obtained and analyzed.

Bridging Calorimetry and Simulation through Precise Calculations of Cucurbituril-Guest Binding Enthalpies.

We used microsecond time scale molecular dynamics simulations to compute, at high precision, binding enthalpies for cucurbit[7]uril (CB7) with eight guests in aqueous solution. The results correlate well with experimental data from previously published isothermal titration calorimetry studies, and decomposition of the computed binding enthalpies by interaction type provides plausible mechanistic insights. Thus, dispersion interactions appear to play a key role in stabilizing these complexes, due at least in part to the fact that their packing density is greater than that of water. On the other hand, strongly favorable Coulombic interactions between the host and guests are compensated by unfavorable solvent contributions, leaving relatively modest electrostatic contributions to the binding enthalpies. The better steric fit of the aliphatic guests into the circular host appears to explain why their binding enthalpies tend to be more favorable than those of the more planar aromatic guests. The present calculations also bear on the validity of the simulation force field. Somewhat unexpectedly, the TIP3P water yields better agreement with experiment than the TIP4P-Ew water model, although the latter is known to replicate the properties of pure water more accurately. More broadly, the present results demonstrate the potential for computational calorimetry to provide atomistic explanations for thermodynamic observations.

Quantum mechanical study of solvent effects in a prototype SN2 reaction in solution: Cl- attack on CH3Cl.

The nucleophilic attack of a chloride ion on methyl chloride is an important prototype SN2 reaction in organic chemistry that is known to be sensitive to the effects of the surrounding solvent. Herein, we develop a highly accurate Specific Reaction Parameter (SRP) model based on the Austin Model 1 Hamiltonian for chlorine to study the effects of solvation into an aqueous environment on the reaction mechanism. To accomplish this task, we apply high-level quantum mechanical calculations to study the reaction in the gas phase and combined quantum mechanical/molecular mechanical simulations with TIP3P and TIP4P-ew water models and the resulting free energy profiles are compared with those determined from simulations using other fast semi-empirical quantum models. Both gas phase and solution results with the SRP model agree very well with experiment and provide insight into the specific role of solvent on the reaction coordinate. Overall, the newly parameterized SRP Hamiltonian is able to reproduce both the gas phase and solution phase barriers, suggesting it is an accurate and robust model for simulations in the aqueous phase at greatly reduced computational cost relative to comparably accurate ab initio and density functional models.