1
|
Fast calculation of hydrogen-bond strengths and free energy of hydration of small molecules. Sci Rep 2023; 13:4143. [PMID: 36914670 PMCID: PMC10011384 DOI: 10.1038/s41598-023-30089-x] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2022] [Accepted: 02/15/2023] [Indexed: 03/16/2023] Open
Abstract
Hydrogen bonding is an interaction of great importance in drug discovery and development as it may significantly affect chemical and biological processes including the interaction of small molecules with other molecules, proteins, and membranes. In particular, hydrogen bonding can impact drug-like properties such as target affinity and oral availability which are critical to developing effective pharmaceuticals, and therefore, numerous methods for the calculation of properties such as hydrogen-bond strengths, free energy of hydration, or water solubility have been proposed over time. However, the accessibility to efficient methods for the predictions of such properties is still limited. Here, we present the development of Jazzy, an open-source tool for the prediction of hydrogen-bond strengths and free energies of hydration of small molecules. Jazzy also allows the visualisation of hydrogen-bond strengths with atomistic resolution to support the design of compounds with desired properties and the interpretation of existing data. The tool is described in its implementation, parameter fitting, and validation against two data sets of experimental hydration free energies. Jazzy is also applied against two chemical series of bioactive compounds to show that hydrogen-bond strengths can be used to understand their structure-activity relationships. Results from the validations highlight the strengths and limitations of Jazzy, and suggest its suitability for interactive design, screening, and machine-learning featurisation.
Collapse
|
2
|
Rahimi AM, Jamali S, Bardhan JP, Lustig SR. Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model. J Chem Theory Comput 2022; 18:5539-5558. [PMID: 36001344 DOI: 10.1021/acs.jctc.2c00248] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. In this work, we develop and assess two multiscale solvation models for computing solvation thermodynamic properties. The new SLIC/CDC model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with new statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). Given the structures of 500 solutes, the SLIC/CDC model predicts Gibbs energies of solvation in water with an average accuracy better than 1 kcal/mol, when compared to experimental measurements, and better than 0.8 kcal/mol, when compared to explicit-solvent molecular dynamics simulations. The individual SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The previously published SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent-accessible surface area (SASA) nonpolar energy mode. With our new, improved parametrization method, the SLIC/SASA model now predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems, compared to experimental and explicit-solvent molecular dynamics, and better than 1.6 kcal/mol average accuracy in ionic liquids, compared to explicit-solvent molecular dynamics. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.
Collapse
Affiliation(s)
- Ali Mehdizadeh Rahimi
- Department of Mechanical and Industrial Engineering, Northeastern University, 360 Huntington Ave., Boston Massachusetts 02115, United States
| | - Safa Jamali
- Department of Mechanical and Industrial Engineering, Northeastern University, 360 Huntington Ave., Boston Massachusetts 02115, United States
| | - Jaydeep P Bardhan
- Pacific Northwest National Laboratory, 902 Battelle Blvd., Richland, Washington 99354, United States
| | - Steven R Lustig
- Department of Chemical Engineering, Northeastern University, 360 Huntington Ave., Boston, Massachusetts 02115, United States
| |
Collapse
|
3
|
Molavi Tabrizi A, Goossens S, Mehdizadeh Rahimi A, Knepley M, Bardhan JP. Predicting solvation free energies and thermodynamics in polar solvents and mixtures using a solvation-layer interface condition. J Chem Phys 2017. [DOI: 10.1063/1.4977037] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Amirhossein Molavi Tabrizi
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, USA
| | - Spencer Goossens
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, USA
| | - Ali Mehdizadeh Rahimi
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, USA
| | - Matthew Knepley
- Department of Computational and Applied Mathematics, Rice University, Houston, Texas 77005, USA
| | - Jaydeep P. Bardhan
- Department of Mechanical and Industrial Engineering, Northeastern University, Boston, Massachusetts 02115, USA
| |
Collapse
|
4
|
Sandberg L. Predicting hydration free energies with chemical accuracy: the SAMPL4 challenge. J Comput Aided Mol Des 2014; 28:211-9. [PMID: 24550133 DOI: 10.1007/s10822-014-9725-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2013] [Accepted: 01/30/2014] [Indexed: 12/14/2022]
Abstract
An implicit solvent model described by a non-simple dielectric medium is used for the prediction of hydration free energies on the dataset of 47 molecules in the SAMPL4 challenge. The solute is represented by a minimal parameter set model based on a new all atom force-field, named the liquid simulation force-field. The importance of a first solvation shell correction to the hydration free energy prediction is discussed and two different approaches are introduced to address it: either with an empirical correction to a few functional groups (alcohol, ether, ester, amines and aromatic nitrogen), or an ab initio correction based on the formation of a solute/explicit water complex. Both approaches give equally good predictions with an average unsigned error <1 kcal/mol. Chemical accuracy is obtained.
Collapse
Affiliation(s)
- Lars Sandberg
- Division of Biological Chemistry and Drug Discovery College of Life Sciences, University of Dundee, Dundee, UK,
| |
Collapse
|
5
|
Parameterization of the hydration free energy computations for organic solutes in the framework of the implicit solvent model with the nonuniform dielectric function. COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2012.12.027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
6
|
|
7
|
Xiong R, Sandler SI, Vlachos DG. Molecular screening of alcohol and polyol adsorption onto MFI-type zeolites. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2012; 28:4491-4499. [PMID: 22320250 DOI: 10.1021/la204710j] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Configurational-bias grand canonical Monte Carlo (CB-GCMC) simulations and expanded ensemble (EE)-CB-GCMC simulations were performed to obtain adsorption isotherms of alcohols and polyols onto MFI-type zeolites from the gas phase and aqueous solution. In adsorption from both phases, Henry's constants and heats of adsorption at infinite dilution for straight-chain alcohols, diols, and triols in silicalite-1 are found to increase, and the saturation loadings decrease with increasing carbon number. Adsorption of straight-chain alcohols is more favorable than that of branched-chain alcohols. Henry's constants increase with increasing number of hydroxyl groups for gas-phase adsorption but decrease for adsorption from aqueous solution due to the strong hydrophilic solvent effect of water. The location of the hydroxyls does not affect significantly the adsorption from aqueous solution but does so in gas-phase adsorption. The saturation pressures for gas-phase adsorption decrease by orders of magnitude from the alcohols to the triols. Nonframework cations increase the adsorption of the small alcohols by an order magnitude at low concentrations (<1 mg/mL), but result in only a factor of 2 increase for larger alcohols like butanol at low concentrations (<0.03 mg/mL), and then decrease the adsorption at higher concentrations. Overall, the simulated results are in reasonable agreement with available experimental data.
Collapse
Affiliation(s)
- Ruichang Xiong
- Department of Chemical Engineering, University of Delaware, Newark, Delaware 19716, United States
| | | | | |
Collapse
|
8
|
Kongsted J, Söderhjelm P, Ryde U. How accurate are continuum solvation models for drug-like molecules? J Comput Aided Mol Des 2009; 23:395-409. [DOI: 10.1007/s10822-009-9271-6] [Citation(s) in RCA: 63] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2009] [Accepted: 04/15/2009] [Indexed: 12/01/2022]
|
9
|
Gong H, Freed KF. Langevin-Debye model for nonlinear electrostatic screening of solvated ions. PHYSICAL REVIEW LETTERS 2009; 102:057603. [PMID: 19257555 DOI: 10.1103/physrevlett.102.057603] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2008] [Revised: 01/06/2009] [Indexed: 05/27/2023]
Abstract
Ion-ion interactions are central to numerous phenomena in geology, biology, and material science, but current understanding is based on linear theories of limited physical applicability. Using the Langevin-Debye model, we investigate how nonlinear dielectric saturation alters screening of electrostatic interactions between solvated charges and show that charge screening in liquids strongly depends on the charge magnitudes but negligibly on the ion radii. Qualitatively different universal screening functions are predicted for ions with charges of the same or opposite signs.
Collapse
Affiliation(s)
- Haipeng Gong
- The James Franck Institute, University of Chicago, Chicago, Illinois 60637, USA
| | | |
Collapse
|
10
|
Bondesson L, Rudberg E, Luo Y, Sałek P. Basis set dependence of solute-solvent interaction energy of benzene in water: A HF/DFT study. J Comput Chem 2008; 29:1725-32. [DOI: 10.1002/jcc.20930] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
11
|
Jha AK, Freed KF. Solvation effect on conformations of 1,2:dimethoxyethane: charge-dependent nonlinear response in implicit solvent models. J Chem Phys 2008; 128:034501. [PMID: 18205504 PMCID: PMC2717614 DOI: 10.1063/1.2815764] [Citation(s) in RCA: 45] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The physical content of and, in particular, the nonlinear contributions from the Langevin-Debye model are illustrated using two applications. First, we provide an improvement in the Langevin-Debye model currently used in some implicit solvent models for computer simulations of solvation free energies of small organic molecules, as well as of biomolecular folding and binding. The analysis is based on the implementation of a charge-dependent Langevin-Debye (qLD) model that is modified by subsequent corrections due to Onsager and Kirkwood. Second, the physical content of the model is elucidated by discussing the general treatment within the LD model of the self-energy of a charge submerged in a dielectric medium for three different limiting conditions and by considering the nonlinear response of the medium. The modified qLD model is used to refine an implicit solvent model (previously applied to protein dynamics). The predictions of the modified implicit solvent model are compared with those from explicit solvent molecular dynamics simulations for the equilibrium conformational populations of 1,2-dimethoxyethane (DME), which is the shortest ether molecule to reproduce the local conformational properties of polyethylene oxide, a polymer with tremendous technological importance and a wide variety of applications. Because the conformational population preferences of DME change dramatically upon solvation, DME is a good test case to validate our modified qLD model. The present analysis of the modified qLD model provides the motivation and tools for studying a wide variety of other interesting systems with heterogeneous dielectric properties and spatial anisotropy.
Collapse
Affiliation(s)
- Abhishek K Jha
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA.
| | | |
Collapse
|
12
|
Bondesson L, Rudberg E, Luo Y, Sałek P. A Linear Scaling Study of Solvent−Solute Interaction Energy of Drug Molecules in Aqua Solution. J Phys Chem B 2007; 111:10320-8. [PMID: 17676891 DOI: 10.1021/jp072621l] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Solvent-solute interaction energies for three well-known drug molecules in water solution are computed at the Hartree-Fock and B3LYP density functional theory levels using a linear scaling technique, which allows one to explicitly include in the model water molecules up to 14 A away from the solute molecule. The dependence of calculated interaction energies on the amount of included solvent has been examined. It is found that it is necessary to account for water molecules within an 8 A radius around the drug molecule to reach the saturated solvent interaction level. Effects of electron correlation and basis set on solvent-solute interaction energies are discussed.
Collapse
Affiliation(s)
- Laban Bondesson
- Department of Theoretical Chemistry, Royal Institute of Technology, SE-10691 Stockholm, Sweden
| | | | | | | |
Collapse
|
13
|
Abstract
A force field for the computer simulation of aqueous solutions of amides is presented. The force field is designed to reproduce the experimentally observed density and Kirkwood-Buff integrals for N-methylacetamide (NMA), allowing for an accurate description of the NMA activity. Other properties such as the translational diffusion constant and heat of mixing are also well reproduced. The force field is then extended to include N,N'-dimethylacetamide and acetamide with good success. Analysis of the simulations of low concentrations of NMA in water indicates a high degree of solvation with only 15% of the NMA molecules involved in solute-solute hydrogen bonding. There is only a weak angular dependence of the solute-solute hydrogen bonding interaction with a minimum at an angle of 65 degrees for the N-H and C=O dipole vectors. The models presented here provide a basis for an accurate force field for peptides and proteins.
Collapse
Affiliation(s)
- Myungshim Kang
- Department of Chemistry, Kansas State University, 111 Willard Hall, Manhattan, Kansas 66506-3701, USA
| | | |
Collapse
|
14
|
Minezawa N, Kato S. Efficient implementation of three-dimensional reference interaction site model self-consistent-field method: Application to solvatochromic shift calculations. J Chem Phys 2007; 126:054511. [PMID: 17302489 DOI: 10.1063/1.2431809] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The authors present an implementation of the three-dimensional reference interaction site model self-consistent-field (3D-RISM-SCF) method. First, they introduce a robust and efficient algorithm for solving the 3D-RISM equation. The algorithm is a hybrid of the Newton-Raphson and Picard methods. The Jacobian matrix is analytically expressed in a computationally useful form. Second, they discuss the solute-solvent electrostatic interaction. For the solute to solvent route, the electrostatic potential (ESP) map on a 3D grid is constructed directly from the electron density. The charge fitting procedure is not required to determine the ESP. For the solvent to solute route, the ESP acting on the solute molecule is derived from the solvent charge distribution obtained by solving the 3D-RISM equation. Matrix elements of the solute-solvent interaction are evaluated by the direct numerical integration. A remarkable reduction in the computational time is observed in both routes. Finally, the authors implement the first derivatives of the free energy with respect to the solute nuclear coordinates. They apply the present method to "solute" water and formaldehyde in aqueous solvent using the simple point charge model, and the results are compared with those from other methods: the six-dimensional molecular Ornstein-Zernike SCF, the one-dimensional site-site RISM-SCF, and the polarizable continuum model. The authors also calculate the solvatochromic shifts of acetone, benzonitrile, and nitrobenzene using the present method and compare them with the experimental and other theoretical results.
Collapse
Affiliation(s)
- Noriyuki Minezawa
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa, Sakyo-ku, Kyoto 606-8502, Japan
| | | |
Collapse
|
15
|
Westergren J, Lindfors L, Höglund T, Lüder K, Nordholm S, Kjellander R. In Silico Prediction of Drug Solubility: 1. Free Energy of Hydration. J Phys Chem B 2007; 111:1872-82. [PMID: 17266351 DOI: 10.1021/jp064220w] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
As a first step in the computational prediction of drug solubility the free energy of hydration, DeltaG*(vw) in TIP4P water has been computed for a data set of 48 drug molecules using the free energy of perturbation method and the optimized potential for liquid simulations all-atom force field. The simulations were performed in two steps, where first the Coulomb and then the Lennard-Jones interactions between the solute and the water molecules were scaled down from full to zero strength to provide physical understanding and simpler predictive models. The results have been interpreted using a theory assuming DeltaG*(vw) = A(MS)gamma + E(LJ) + E(C)/2 where A(MS) is the molecular surface area, gamma is the water-vapor surface tension, and E(LJ) and E(C) are the solute-water Lennard-Jones and Coulomb interaction energies, respectively. It was found that by a proper definition of the molecular surface area our results as well as several results from the literature were found to be in quantitative agreement using the macroscopic surface tension of TIP4P water. This is in contrast to the surface tension for water around a spherical cavity that previously has been shown to be dependent on the size of the cavity up to a radius of approximately 1 nm. The step of scaling down the electrostatic interaction can be represented by linear response theory.
Collapse
Affiliation(s)
- Jan Westergren
- Pharmaceutical and Analytical R&D, AstraZeneca R&D, Mölndal, SE-431 83 Mölndal, Sweden
| | | | | | | | | | | |
Collapse
|
16
|
Muthukrishnan A, Sangaranarayanan M. Hydration energies of C60 and C70 fullerenes – A novel Monte Carlo simulation study. Chem Phys 2007. [DOI: 10.1016/j.chemphys.2006.10.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
17
|
Amovilli C, Filippi C, Floris FM. Coupling Quantum Monte Carlo to a Nonlinear Polarizable Continuum Model for Spherical Solutes. J Phys Chem B 2006; 110:26225-31. [PMID: 17181280 DOI: 10.1021/jp066006d] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Starting from the nonlinear dielectric response model of Sandberg and Edholm, we derive an analytical expression of the polarization contribution to the solvation free energy in terms of the electronic density of the solute and the dielectric properties of the solvent. The solvent inhomogeneity is taken into account with the use of a smooth switching function whose spacial variation is established on the basis of how the solvent is arranged around the solute. An explicit form of a local potential representing the solvent effect on the solute is thus obtained by functional analysis. This effective potential can be combined with density functional or quantum chemical methods for the quantum mechanical treatment of the solute. Here, we use quantum Monte Carlo techniques for the solute and apply the method to the hydration of atomic ions finding very good agreement with experimental data.
Collapse
Affiliation(s)
- Claudio Amovilli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, I-56126 Pisa, Italy.
| | | | | |
Collapse
|
18
|
|
19
|
Basilevsky MV, Leontyev IV, Luschekina SV, Kondakova OA, Sulimov VB. Computation of hydration free energies of organic solutes with an implicit water model. J Comput Chem 2006; 27:552-70. [PMID: 16463371 DOI: 10.1002/jcc.20332] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
A new approach for computing hydration free energies DeltaG(solv) of organic solutes is formulated and parameterized. The method combines a conventional PCM (polarizable continuum model) computation for the electrostatic component DeltaG(el) of DeltaG(solv) and a specially detailed algorithm for treating the complementary nonelectrostatic contributions (DeltaG(nel)). The novel features include the following: (a) two different cavities are used for treating DeltaG(el) and DeltaG(nel). For the latter case the cavity is larger and based on thermal atomic radii (i.e., slightly reduced van der Waals radii). (b) The cavitation component of DeltaG(nel) is taken to be proportional to the volume of the large cavity. (c) In the treatment of van der Waals interactions, all solute atoms are counted explicitly. The corresponding interaction energies are computed as integrals over the surface of the larger cavity; they are based on Lennard Jones (LJ) type potentials for individual solute atoms. The weighting coefficients of these LJ terms are considered as fitting parameters. Testing this method on a collection of 278 uncharged organic solutes gave satisfactory results. The average error (RMSD) between calculated and experimental free energy values varies between 0.15 and 0.5 kcal/mol for different classes of solutes. The larger deviations found for the case of oxygen compounds are probably due to a poor approximation of H-bonding in terms of LJ potentials. For the seven compounds with poorest fit to experiment, the error exceeds 1.5 kcal/mol; these outlier points were not included in the parameterization procedure. Several possible origins of these errors are discussed.
Collapse
Affiliation(s)
- Mikhail V Basilevsky
- Karpov Institute of Physical Chemistry, ul. Vorontsovo Pole 10, 105064 Moscow, Russia
| | | | | | | | | |
Collapse
|
20
|
|
21
|
Affiliation(s)
- Jacopo Tomasi
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento 35, 56126 Pisa, Italy.
| | | | | |
Collapse
|
22
|
Sudha V, Harinipriya S, Sangaranarayanan M. Dehydration energies of alkaline earth metal halides – a novel simulation methodology. Chem Phys 2005. [DOI: 10.1016/j.chemphys.2004.10.030] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
23
|
Sudha V, Harinipriya S, Sangaranarayanan MV. A simple simulation methodology for estimation of dehydration energies and surface potentials of concentrated NaCl solutions. J Colloid Interface Sci 2004; 280:139-48. [PMID: 15476784 DOI: 10.1016/j.jcis.2004.07.036] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2004] [Accepted: 07/21/2004] [Indexed: 10/26/2022]
Abstract
A novel simulation procedure for estimating the dehydration energies of NaCl solutions spanning a wide range of concentrations, which incorporates ionic and molecular sizes, ion-pair formation, etc., is proposed on geometric and phenomenological considerations. The extent of dehydration during each movement of the hydrated molecule is evaluated using the expected and actual displacement of the species and mean nearest-neighbor distances. The interdependence between the size of the simulation box, number of molecules, and electrolyte concentration is pointed out.
Collapse
Affiliation(s)
- V Sudha
- Department of Chemistry, Indian Institute of Technology, Madras 600-036, India
| | | | | |
Collapse
|
24
|
Sudha V, Harinipriya S, Sangaranarayanan M. Hydration energies of trihalides of lanthanide and actinide series—a novel simulation methodology. ACTA ACUST UNITED AC 2004. [DOI: 10.1016/j.theochem.2004.07.008] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
25
|
Nagy PI. Effects of the Solute Model and Concentration on the Calculated Free Energy of Hydration in Explicit Solvent Solution. J Phys Chem B 2004. [DOI: 10.1021/jp049249j] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- P. I. Nagy
- Department of Medicinal and Biological Chemistry, The University of Toledo, Toledo, Ohio 43606-3390
| |
Collapse
|
26
|
Nagy PI, Takács-Novák K. Tautomeric and conformational equilibria of biologically important (hydroxyphenyl)alkylamines in the gas phase and in aqueous solution. Phys Chem Chem Phys 2004. [DOI: 10.1039/b314924b] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|