1
|
Zhang S, Giese TJ, Lee TS, York DM. Alchemical Enhanced Sampling with Optimized Phase Space Overlap. J Chem Theory Comput 2024; 20:3935-3953. [PMID: 38666430 PMCID: PMC11157682 DOI: 10.1021/acs.jctc.4c00251] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2024]
Abstract
An alchemical enhanced sampling (ACES) method has recently been introduced to facilitate importance sampling in free energy simulations. The method achieves enhanced sampling from Hamiltonian replica exchange within a dual topology framework while utilizing new smoothstep softcore potentials. A common sampling problem encountered in lead optimization is the functionalization of aromatic rings that exhibit distinct conformational preferences when interacting with the protein. It is difficult to converge the distribution of ring conformations due to the long time scale of ring flipping events; however, the ACES method addresses this issue by modeling the syn and anti ring conformations within a dual topology. ACES thereby samples the conformer distributions by alchemically tunneling between states, as opposed to traversing a physical pathway with a high rotational barrier. We demonstrate the use of ACES to overcome conformational sampling issues involving ring flipping in ML300-derived noncovalent inhibitors of SARS-CoV-2 Main Protease (Mpro). The demonstrations explore how the use of replica exchange and the choice of softcore selection affects the convergence of the ring conformation distributions. Furthermore, we examine how the accuracy of the calculated free energies is affected by the degree of phase space overlap (PSO) between adjacent states (i.e., between neighboring λ-windows) and the Hamiltonian replica exchange acceptance ratios. Both of these factors are sensitive to the spacing between the intermediate states. We introduce a new method for choosing a schedule of λ values. The method analyzes short "burn-in" simulations to construct a 2D map of the nonlocal PSO. The schedule is obtained by optimizing an alchemical pathway on the 2D map that equalizes the PSO between the λ intervals. The optimized phase space overlap λ-spacing method (Opt-PSO) leads to more numerous end-to-end single passes and round trips due to the correlation between PSO and Hamiltonian replica exchange acceptance ratios. The improved exchange statistics enhance the efficiency of ACES method. The method has been implemented into the FE-ToolKit software package, which is freely available.
Collapse
Affiliation(s)
- Shi Zhang
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Timothy J. Giese
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Tai-Sung Lee
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, NJ 08854, USA
| |
Collapse
|
2
|
York DM. Modern Alchemical Free Energy Methods for Drug Discovery Explained. ACS PHYSICAL CHEMISTRY AU 2023; 3:478-491. [PMID: 38034038 PMCID: PMC10683484 DOI: 10.1021/acsphyschemau.3c00033] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 09/12/2023] [Accepted: 09/13/2023] [Indexed: 12/02/2023]
Abstract
This Perspective provides a contextual explanation of the current state-of-the-art alchemical free energy methods and their role in drug discovery as well as highlights select emerging technologies. The narrative attempts to answer basic questions about what goes on "under the hood" in free energy simulations and provide general guidelines for how to run simulations and analyze the results. It is the hope that this work will provide a valuable introduction to students and scientists in the field.
Collapse
Affiliation(s)
- Darrin M. York
- Laboratory for Biomolecular
Simulation Research, Institute for Quantitative Biomedicine, and Department
of Chemistry and Chemical Biology, Rutgers
University, Piscataway, New Jersey 08854, United States
| |
Collapse
|
3
|
Stones AE, Aarts DGAL. Measuring many-body distribution functions in fluids using test-particle insertion. J Chem Phys 2023; 159:194502. [PMID: 37975484 DOI: 10.1063/5.0172664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Accepted: 10/19/2023] [Indexed: 11/19/2023] Open
Abstract
We derive a hierarchy of equations, which allow a general n-body distribution function to be measured by test-particle insertion of between 1 and n particles. We apply it to measure the pair and three-body distribution functions in a simple fluid using snapshots from Monte Carlo simulations in the grand canonical ensemble. The resulting distribution functions obtained from insertion methods are compared with the conventional distance-histogram method: the insertion approach is shown to overcome the drawbacks of the histogram method, offering enhanced structural resolution and a more straightforward normalization. At high particle densities, the insertion method starts breaking down, which can be delayed by utilizing the underlying hierarchical structure of the insertion method. Our method will be especially useful in characterizing the structure of inhomogeneous fluids and investigating closure approximations in liquid state theory.
Collapse
Affiliation(s)
- Adam Edward Stones
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Dirk G A L Aarts
- Department of Chemistry, Physical and Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
4
|
Boulougouris GC. Accessible Molecular System Creator: Building Molecular Configurations Based on the Inaccessible Molecular Volume and Accessible Molecular Surface via Static Monte Carlo Sampling. J Phys Chem B 2023; 127:9520-9531. [PMID: 37883744 DOI: 10.1021/acs.jpcb.3c03670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
Monte Carlo (MC) stochastic sampling is a powerful tool in classical molecular simulations that directly connects the observable macroscopic properties of matter and the underlying atomistic interactions. This connection operates within the framework of the statistical mechanics proposed by Gibbs. Most MC simulations are "dynamic," creating statistical ensembles of microstates via a Markovian chain, where each microstate in the ensemble depends only on its previous microstate. Herein, we re-examine an alternative form of MC that generates ensemble members through a "static" approach, building molecular systems stepwise. The basic theory for such an approach traces back to Rosenbluth and Rosenbluth, who proposed "static" stepwise sampling of a polymeric chain. It is almost as old as the Metropolis importance sampling approach used in dynamic MC, although the latter has been considerably more popular than the former. Herein, we address the main obstacle in static MC that has hindered the widespread adoption of Rosenbluth-based approaches in atomistic simulations. The obstacle lies in mapping the molecular accessible volume for adding a molecule in a Rosenbluth-like static sampling of atomistic configurations. We demonstrate a breakthrough by leveraging the ability to analytically map the inaccessible molecular volume and the accessible molecular surface owing to interatomically excluded volume interactions. This advance substantially enhances the ability to create molecular samples using a Rosenbluth-like static building process. The proposed approach can be used as a tool for creating initial configurations in MC or molecular dynamics simulations─a field where Rosenbluth-like static building has been applied. Additionally, this approach can be used as the first step in a perturbation scheme that accurately estimates free energy differences by estimating the chemical work related to molecule addition, removal, or reinsertion within the context of free energy perturbation schemes employed in molecular simulations.
Collapse
Affiliation(s)
- Georgios C Boulougouris
- Laboratory of Computational Physical Chemistry, Department of Molecular Biology and Genetics, University of Thrace, GR 681 00 Alexandroupoulis, Greece
| |
Collapse
|
5
|
Kournopoulos S, Santos MS, Ravipati S, Haslam AJ, Jackson G, Economou IG, Galindo A. The Contribution of the Ion-Ion and Ion-Solvent Interactions in a Molecular Thermodynamic Treatment of Electrolyte Solutions. J Phys Chem B 2022; 126:9821-9839. [PMID: 36395498 PMCID: PMC9720728 DOI: 10.1021/acs.jpcb.2c03915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
Developing molecular equations of state to treat electrolyte solutions is challenging due to the long-range nature of the Coulombic interactions. Seminal approaches commonly used are the mean spherical approximation (MSA) and the Debye-Hückel (DH) theory to account for ion-ion interactions and, often, the Born theory of solvation for ion-solvent interactions. We investigate the accuracy of the MSA and DH approaches using each to calculate the contribution of the ion-ion interactions to the chemical potential of NaCl in water, comparing these with newly computer-generated simulation data; the ion-ion contribution is isolated by selecting an appropriate primitive model with a Lennard-Jones force field to describe the solvent. A study of mixtures with different concentrations and ionic strengths reveals that the calculations from both MSA and DH theories are of similar accuracy, with the MSA approach resulting in marginally better agreement with the simulation data. We also demonstrate that the Born theory provides a good qualitative description of the contribution of the ion-solvent interactions; we employ an explicitly polar water model in these simulations. Quantitative agreement up to moderate salt concentrations and across the relevant range of temperature is achieved by adjusting the Born radius using simulation data of the free energy of solvation. We compute the radial and orientational distribution functions of the systems, thereby providing further insight on the differences observed between the theory and simulation. We thus provide rigorous benchmarks for use of the MSA, DH, and Born theories as perturbation approaches, which will be of value for improving existing models of electrolyte solutions, especially in the context of equations of state.
Collapse
Affiliation(s)
- Spiros Kournopoulos
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering,
and Institute for Molecular Science and Engineering, Imperial College, London, London SW7 2AZ, United Kingdom
| | - Mirella Simões Santos
- Laboratoire
de Chimie, École Normale Supérieure
de Lyon, 46 Allée d’Italie, 69364 Lyon, France,Australian
Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Srikanth Ravipati
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering,
and Institute for Molecular Science and Engineering, Imperial College, London, London SW7 2AZ, United Kingdom
| | - Andrew J. Haslam
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering,
and Institute for Molecular Science and Engineering, Imperial College, London, London SW7 2AZ, United Kingdom
| | - George Jackson
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering,
and Institute for Molecular Science and Engineering, Imperial College, London, London SW7 2AZ, United Kingdom
| | - Ioannis G. Economou
- Chemical
Engineering Program, Texas A&M University
at Qatar, Doha 23874, Qatar
| | - Amparo Galindo
- Department
of Chemical Engineering, Sargent Centre for Process Systems Engineering,
and Institute for Molecular Science and Engineering, Imperial College, London, London SW7 2AZ, United Kingdom,
| |
Collapse
|
6
|
Sabatino SJ, Paluch AS. Predicting octanol/water partition coefficients using molecular simulation for the SAMPL7 challenge: comparing the use of neat and water saturated 1-octanol. J Comput Aided Mol Des 2021; 35:1009-1024. [PMID: 34495430 DOI: 10.1007/s10822-021-00415-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Accepted: 08/09/2021] [Indexed: 11/29/2022]
Abstract
Blind predictions of octanol/water partition coefficients at 298.15 K for 22 drug-like compounds were made for the SAMPL7 challenge. The octanol/water partition coefficients were predicted using solvation free energies computed using molecular dynamics simulations, wherein we considered the use of both pure and water-saturated 1-octanol to model the octanol-rich phase. Water and 1-octanol were modeled using TIP4P and TrAPPE-UA, respectively, which have been shown to well reproduce the experimental mutual solubility, and the solutes were modeled using GAFF. After the close of the SAMPL7 challenge, we additionally made predictions using TIP4P/2005 water. We found that the predictions were sensitive to the choice of water force field. However, the effect of water in the octanol-rich phase was found to be even more significant and non-negligible. The effect of inclusion of water was additionally sensitive to the chemical structure of the solute.
Collapse
Affiliation(s)
- Spencer J Sabatino
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, OH, 45056, USA
| | - Andrew S Paluch
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, OH, 45056, USA.
| |
Collapse
|
7
|
Schultz AJ, Kofke DA. Identifying and estimating bias in overlap-sampling free-energy calculations. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2020.1758695] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Andrew J. Schultz
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, NY, USA
| | - David A. Kofke
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, NY, USA
| |
Collapse
|
8
|
Rahbari A, Hens R, Ramdin M, Moultos OA, Dubbeldam D, Vlugt TJH. Recent advances in the continuous fractional component Monte Carlo methodology. MOLECULAR SIMULATION 2020. [DOI: 10.1080/08927022.2020.1828585] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| | - D. Dubbeldam
- Van 't Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, Netherlands
| |
Collapse
|
9
|
Stephan S, Thol M, Vrabec J, Hasse H. Thermophysical Properties of the Lennard-Jones Fluid: Database and Data Assessment. J Chem Inf Model 2019; 59:4248-4265. [DOI: 10.1021/acs.jcim.9b00620] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Simon Stephan
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| | - Monika Thol
- Thermodynamics, Ruhr University Bochum, 44801 Bochum, Germany
| | - Jadran Vrabec
- Thermodynamics and Process Engineering, TU Berlin, 10587 Berlin, Germany
| | - Hans Hasse
- Laboratory of Engineering Thermodynamics (LTD), TU Kaiserslautern, 67663 Kaiserslautern, Germany
| |
Collapse
|
10
|
Giese TJ, York DM. Development of a Robust Indirect Approach for MM → QM Free Energy Calculations That Combines Force-Matched Reference Potential and Bennett's Acceptance Ratio Methods. J Chem Theory Comput 2019; 15:5543-5562. [PMID: 31507179 DOI: 10.1021/acs.jctc.9b00401] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We use the PBE0/6-31G* density functional method to perform ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations under periodic boundary conditions with rigorous electrostatics using the ambient potential composite Ewald method in order to test the convergence of MM → QM/MM free energy corrections for the prediction of 17 small-molecule solvation free energies and eight ligand binding free energies to T4 lysozyme. The "indirect" thermodynamic cycle for calculating free energies is used to explore whether a series of reference potentials improve the statistical quality of the predictions. Specifically, we construct a series of reference potentials that optimize a molecular mechanical (MM) force field's parameters to reproduce the ab initio QM/MM forces from a QM/MM simulation. The optimizations form a systematic progression of successively expanded parameters that include bond, angle, dihedral, and charge parameters. For each reference potential, we calculate benchmark quality reference values for the MM → QM/MM correction by performing the mixed MM and QM/MM Hamiltonians at 11 intermediate states, each for 200 ps. We then compare forward and reverse application of Zwanzig's relation, thermodynamic integration (TI), and Bennett's acceptance ratio (BAR) methods as a function of reference potential, simulation time, and the number of simulated intermediate states. We find that Zwanzig's equation is inadequate unless a large number of intermediate states are explicitly simulated. The TI and BAR mean signed errors are very small even when only the end-state simulations are considered, and the standard deviations of the TI and BAR errors are decreased by choosing a reference potential that optimizes the bond and angle parameters. We find a robust approach for the data sets of fairly rigid molecules considered here is to use bond + angle reference potential together with the end-state-only BAR analysis. This requires QM/MM simulations to be performed in order to generate reference data to parametrize the bond + angle reference potential, and then this same simulation serves a dual purpose as the full QM/MM end state. The convergence of the results with respect to time suggests that computational resources may be used more efficiently by running multiple simulations for no more than 50 ps, rather than running one long simulation.
Collapse
Affiliation(s)
- Timothy J Giese
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| | - Darrin M York
- Laboratory for Biomolecular Simulation Research, Center for Integrative Proteomics Research and Department of Chemistry and Chemical Biology , Rutgers University , Piscataway , New Jersey 08854-8087 , United States
| |
Collapse
|
11
|
Roese SN, Margulis GV, Schmidt AJ, Uzat CB, Heintz JD, Paluch AS. A Simple Method to Predict and Interpret the Formation of Azeotropes in Binary Systems Using Conventional Solvation Free Energy Calculations. Ind Eng Chem Res 2019. [DOI: 10.1021/acs.iecr.9b03694] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Sydnee N. Roese
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Griffin V. Margulis
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Alexa J. Schmidt
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06520, United States
| | - Cole B. Uzat
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Justin D. Heintz
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Andrew S. Paluch
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| |
Collapse
|
12
|
Stones AE, Dullens RPA, Aarts DGAL. Model-Free Measurement of the Pair Potential in Colloidal Fluids Using Optical Microscopy. PHYSICAL REVIEW LETTERS 2019; 123:098002. [PMID: 31524476 DOI: 10.1103/physrevlett.123.098002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Indexed: 06/10/2023]
Abstract
We report a straightforward, model-free approach for measuring pair potentials from particle-coordinate data, based on enforcing consistency between the pair distribution function measured separately by the distance-histogram and test-particle insertion routes. We demonstrate the method's accuracy and versatility in simulations of simple fluids, before applying it to an experimental system composed of superparamagnetic colloidal particles. The method will enable experimental investigations into many-body interactions and allow for effective coarse graining of interactions from simulations.
Collapse
Affiliation(s)
- Adam Edward Stones
- Department of Chemistry, Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Roel P A Dullens
- Department of Chemistry, Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| | - Dirk G A L Aarts
- Department of Chemistry, Physical & Theoretical Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
13
|
Josephson TR, Singh R, Minkara MS, Fetisov EO, Siepmann JI. Partial molar properties from molecular simulation using multiple linear regression. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1648898] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Tyler R. Josephson
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| | - Ramanish Singh
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| | - Mona S. Minkara
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
| | - Evgenii O. Fetisov
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Chemical Physics and Analysis, Pacific Northwest National Laboratory, Richland, Washington, USA
| | - J. Ilja Siepmann
- Department of Chemistry and Chemical Theory Center, Minneapolis, MN, USA
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, MN, USA
| |
Collapse
|
14
|
Long GE, Dhakal P, Redeker BN, Paluch AS. Using limiting activity coefficients to efficiently evaluate the ability of fixed-charge force fields to model miscible water plus cosolvent mixtures. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1531399] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Garrett E. Long
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Pratik Dhakal
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Bryce N. Redeker
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Andrew S. Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| |
Collapse
|
15
|
Dasari S, Mallik BS. Association of Nucleobases in Hydrated Ionic Liquid from Biased Molecular Dynamics Simulations. J Phys Chem B 2018; 122:9635-9645. [PMID: 30260229 DOI: 10.1021/acs.jpcb.8b05778] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
We employed metadynamics-based classical molecular dynamics simulations to methylated adenine-thymine (mA-mT) and guanine-cytosine (mG-mC) base pairs to see favorable conformations in various concentrations of hydrated 1-ethyl, 3-methyl imidazolium acetate. We investigated various stacked and hydrogen-bonded conformations of association of base pairs through appropriately chosen collective variables. Stacked conformations more favored in water for both base pairs, whereas Watson-Crick (WC) hydrogen-bonding conformations are favored in pure and hydrated ionic liquids (ILs) except for 0.75 mol fraction IL. We observe that EMIm cations surround the base pairs in WC conformations creating a kind of hydrophobic cavity and protect the hydrogen bonds between base pairs. However, the five-membered heteroaromatic rings of cations stack with the nucleobases in the cation-base-cation (π-π-π) model, which resembles the base-base-base stacking in a DNA duplex. Interestingly, from additional simulations of 0.5 mol fraction hydrated choline dihydrogen phosphate IL, we observe that the stacked conformations become more favored than the WC conformation due to the absence of π-bonds in cations. The calculated values of relative solubility of base pairs in pure and hydrated ionic liquids compared to those in pure water correlate well with the free energy values of WC and stacked conformations.
Collapse
Affiliation(s)
- Sathish Dasari
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi , Sangareddy 502285 , Telangana , India
| | - Bhabani S Mallik
- Department of Chemistry , Indian Institute of Technology Hyderabad , Kandi , Sangareddy 502285 , Telangana , India
| |
Collapse
|
16
|
Perego C, Valsson O, Parrinello M. Chemical potential calculations in non-homogeneous liquids. J Chem Phys 2018; 149:072305. [PMID: 30134712 DOI: 10.1063/1.5024631] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
The numerical computation of chemical potential in dense non-homogeneous fluids is a key problem in the study of confined fluid thermodynamics. To this day, several methods have been proposed; however, there is still need for a robust technique, capable of obtaining accurate estimates at large average densities. A widely established technique is the Widom insertion method, which computes the chemical potential by sampling the energy of insertion of a test particle. Non-homogeneity is accounted for by assigning a density dependent weight to the insertion points. However, in dense systems, the poor sampling of the insertion energy is a source of inefficiency, hampering a reliable convergence. We have recently presented a new technique for the chemical potential calculation in homogeneous fluids. This novel method enhances the sampling of the insertion energy via well-tempered metadynamics, reaching accurate estimates at very large densities. In this paper, we extend the technique to the case of non-homogeneous fluids. The method is successfully tested on a confined Lennard-Jones fluid. In particular, we show that, thanks to the improved sampling, our technique does not suffer from a systematic error that affects the classic Widom method for non-homogeneous fluids, providing a precise and accurate result.
Collapse
Affiliation(s)
- Claudio Perego
- Department of Polymer Theory, Max-Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Omar Valsson
- Department of Polymer Theory, Max-Planck Institute for Polymer Research, Ackermannweg 10, D-55128 Mainz, Germany
| | - Michele Parrinello
- Department of Chemistry and Applied Biosciences, ETH Zurich, c/o USI Campus, Via Giuseppe Buffi 13, CH-6900 Lugano, Switzerland
| |
Collapse
|
17
|
Li L, Totton T, Frenkel D. Computational methodology for solubility prediction: Application to sparingly soluble organic/inorganic materials. J Chem Phys 2018; 149:054102. [PMID: 30089373 DOI: 10.1063/1.5040366] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The solubility of a crystalline material can be estimated from the absolute free energy of the solid and the excess solvation free energy. In the earlier work, we presented a general-purpose molecular-dynamics-based methodology enabling solubility predictions of crystalline compounds, yielding accurate estimates of the aqueous solubilities of naphthalene at various pressures and temperatures. In the present work, we investigate a number of prototypical complex materials, including phenanthrene, calcite, and aragonite over a range of temperatures and pressures. Our results provide stronger evidence for the power of the methodology for universal solubility predictions.
Collapse
Affiliation(s)
- Lunna Li
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Tim Totton
- BP Exploration Operating Co. Ltd., Sunbury-on-Thames TW16 7LN, United Kingdom
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| |
Collapse
|
18
|
A new interpretation of chemical potential in adsorption systems and the vapour–liquid interface. ADSORPTION 2018. [DOI: 10.1007/s10450-018-9957-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
19
|
Li L, Totton T, Frenkel D. Computational methodology for solubility prediction: Application to the sparingly soluble solutes. J Chem Phys 2018; 146:214110. [PMID: 28595415 DOI: 10.1063/1.4983754] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The solubility of a crystalline substance in the solution can be estimated from its absolute solid free energy and excess solvation free energy. Here, we present a numerical method, which enables convenient solubility estimation of general molecular crystals at arbitrary thermodynamic conditions where solid and solution can coexist. The methodology is based on standard alchemical free energy methods, such as thermodynamic integration and free energy perturbation, and consists of two parts: (1) systematic extension of the Einstein crystal method to calculate the absolute solid free energies of molecular crystals at arbitrary temperatures and pressures and (2) a flexible cavity method that can yield accurate estimates of the excess solvation free energies. As an illustration, via classical Molecular Dynamic simulations, we show that our approach can predict the solubility of OPLS-AA-based (Optimized Potentials for Liquid Simulations All Atomic) naphthalene in SPC (Simple Point Charge) water in good agreement with experimental data at various temperatures and pressures. Because the procedure is simple and general and only makes use of readily available open-source software, the methodology should provide a powerful tool for universal solubility prediction.
Collapse
Affiliation(s)
- Lunna Li
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Tim Totton
- BP Exploration Operating Co. Ltd., Sunbury-on-Thames TW16 7LN, United Kingdom
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| |
Collapse
|
20
|
Heidari M, Kremer K, Cortes-Huerto R, Potestio R. Spatially Resolved Thermodynamic Integration: An Efficient Method To Compute Chemical Potentials of Dense Fluids. J Chem Theory Comput 2018; 14:3409-3417. [DOI: 10.1021/acs.jctc.8b00002] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Maziar Heidari
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
| | | | - Raffaello Potestio
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz, Germany
- Physics Department, University of Trento, via Sommarive 14 Povo, Trento 38123, Italy
- INFN-TIFPA, Trento Institute for Fundamental Physics and Applications, I-38123 Trento, Italy
| |
Collapse
|
21
|
Rahbari A, Hens R, Nikolaidis IK, Poursaeidesfahani A, Ramdin M, Economou IG, Moultos OA, Dubbeldam D, Vlugt TJH. Computation of partial molar properties using continuous fractional component Monte Carlo. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1451663] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- A. Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - R. Hens
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - I. K. Nikolaidis
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research Demokritos, Institute of Nanoscience and Nanotechnology , Attikis, Greece
| | - A. Poursaeidesfahani
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - M. Ramdin
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - I. G. Economou
- Molecular Thermodynamics and Modeling of Materials Laboratory, National Center for Scientific Research Demokritos, Institute of Nanoscience and Nanotechnology , Attikis, Greece
- Chemical Engineering Program, Texas A&M University at Qatar , Doha, Qatar
| | - O. A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| | - D. Dubbeldam
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam , Amsterdam, The Netherlands
| | - T. J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering , Delft, The Netherlands
| |
Collapse
|
22
|
Rahbari A, Poursaeidesfahani A, Torres-Knoop A, Dubbeldam D, Vlugt TJH. Chemical potentials of water, methanol, carbon dioxide and hydrogen sulphide at low temperatures using continuous fractional component Gibbs ensemble Monte Carlo. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1391385] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Ahmadreza Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ali Poursaeidesfahani
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| | - Ariana Torres-Knoop
- Van’t Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - David Dubbeldam
- Van’t Hoff Institute for Molecular Sciences, University of Amsterdam, Amsterdam, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
23
|
Noroozi J, Paluch AS. Microscopic Structure and Solubility Predictions of Multifunctional Solids in Supercritical Carbon Dioxide: A Molecular Simulation Study. J Phys Chem B 2017; 121:1660-1674. [DOI: 10.1021/acs.jpcb.6b12390] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Javad Noroozi
- Department
of Chemical and Petroleum Engineering, Sharif University of Technology, Tehran, Iran
| | - Andrew S. Paluch
- Department
of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| |
Collapse
|
24
|
Phifer JR, Cox CE, da Silva LF, Nogueira GG, Barbosa AKP, Ley RT, Bozada SM, O’Loughlin EJ, Paluch AS. Predicting the equilibrium solubility of solid polycyclic aromatic hydrocarbons and dibenzothiophene using a combination of MOSCED plus molecular simulation or electronic structure calculations. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1284356] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Jeremy R. Phifer
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Courtney E. Cox
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Larissa Ferreira da Silva
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
- Departamento de Engenharia Química, Universidade Federal de São João del-Rei, Ouro Branco, Brazil
| | - Gabriel Gonçalves Nogueira
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
- Departamento de Engenharia Química, Universidade Federal de São João del-Rei, Ouro Branco, Brazil
| | - Ana Karolyne Pereira Barbosa
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
- Instituto de Ciência e Tecnologia, Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, Brazil
| | - Ryan T. Ley
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Samantha M. Bozada
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | | | - Andrew S. Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| |
Collapse
|
25
|
Cox CE, Phifer JR, Ferreira da Silva L, Gonçalves Nogueira G, Ley RT, O’Loughlin EJ, Pereira Barbosa AK, Rygelski BT, Paluch AS. Combining MOSCED with molecular simulation free energy calculations or electronic structure calculations to develop an efficient tool for solvent formulation and selection. J Comput Aided Mol Des 2017; 31:183-199. [DOI: 10.1007/s10822-016-0001-6] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Accepted: 12/16/2016] [Indexed: 11/28/2022]
|
26
|
Caudle MR, Cox CE, Ley RT, Paluch AS. A molecular study of the wastewater contaminants atenolol and atrazine in 1-n-butyl-3-methylimidazolium based ionic liquids for potential treatment applications. Mol Phys 2017. [DOI: 10.1080/00268976.2016.1278478] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Miranda R. Caudle
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Courtney E. Cox
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Ran T. Ley
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| | - Andrew S. Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, OH, USA
| |
Collapse
|
27
|
Ustinov EA. Thermodynamics and simulation of hard-sphere fluid and solid: Kinetic Monte Carlo method versus standard Metropolis scheme. J Chem Phys 2017; 146:034110. [DOI: 10.1063/1.4974141] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
28
|
Affiliation(s)
- Stefan Boresch
- University of Vienna, Faculty of Chemistry, Department of Computational Biological Chemistry, Vienna, Austria
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, FL, USA
| |
Collapse
|
29
|
Ley RT, Paluch AS. Understanding the large solubility of lidocaine in 1-n-butyl-3-methylimidazolium based ionic liquids using molecular simulation. J Chem Phys 2016; 144:084501. [PMID: 26931706 DOI: 10.1063/1.4942025] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Room temperature ionic liquids have been proposed as replacement solvents in a wide range of industrial separation processes. Here, we focus on the use of ionic liquids as solvents for the pharmaceutical compound lidocaine. We show that the solubility of lidocaine in seven common 1-n-butyl-3-methylimidazolium based ionic liquids is greatly enhanced relative to water. The predicted solubility is greatest in [BMIM](+)[CH3CO2](-), which we find results from favorable hydrogen bonding between the lidocaine amine hydrogen and the [CH3CO2](-) oxygen, favorable electrostatic interactions between the lidocaine amide oxygen with the [BMIM](+) aromatic ring hydrogens, while lidocaine does not interfere with the association of [BMIM](+) with [CH3CO2](-). Additionally, by removing functional groups from the lidocaine scaffold while maintaining the important amide group, we found that as the van der Waals volume increases, solubility in [BMIM](+)[CH3CO2](-) relative to water increases.
Collapse
Affiliation(s)
- Ryan T Ley
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, Ohio 45056, USA
| | - Andrew S Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University, Oxford, Ohio 45056, USA
| |
Collapse
|
30
|
McDonnell MT, Greeley DA, Kit KM, Keffer DJ. Molecular Dynamics Simulations of Hydration Effects on Solvation, Diffusivity, and Permeability in Chitosan/Chitin Films. J Phys Chem B 2016; 120:8997-9010. [DOI: 10.1021/acs.jpcb.6b05999] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Marshall T. McDonnell
- Department of Chemical and Biomolecular
Engineering,
and ‡Department of Materials
Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Duncan A. Greeley
- Department of Chemical and Biomolecular
Engineering,
and ‡Department of Materials
Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Kevin M. Kit
- Department of Chemical and Biomolecular
Engineering,
and ‡Department of Materials
Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - David J. Keffer
- Department of Chemical and Biomolecular
Engineering,
and ‡Department of Materials
Science and Engineering, University of Tennessee, Knoxville, Tennessee 37996, United States
| |
Collapse
|
31
|
Padmanabhan P, Martinez-Veracoechea F, Escobedo FA. Computation of Free Energies of Cubic Bicontinuous Phases for Blends of Diblock Copolymer and Selective Homopolymer. Macromolecules 2016. [DOI: 10.1021/acs.macromol.6b00123] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
- Poornima Padmanabhan
- School of Chemical & Biomolecular Engineering, Cornell University, Ithaca, New York 14853, United States
| | | | - Fernando A. Escobedo
- School of Chemical & Biomolecular Engineering, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
32
|
Ley RT, Fuerst GB, Redeker BN, Paluch AS. Developing a Predictive Form of MOSCED for Nonelectrolyte Solids Using Molecular Simulation: Application to Acetanilide, Acetaminophen, and Phenacetin. Ind Eng Chem Res 2016. [DOI: 10.1021/acs.iecr.5b04807] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ryan T. Ley
- Department
of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Georgia B. Fuerst
- Department
of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Bryce N. Redeker
- Department
of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Andrew S. Paluch
- Department
of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| |
Collapse
|
33
|
Paluch AS, Lourenço TC, Han F, Costa LT. Understanding the Solubility of Acetaminophen in 1-n-Alkyl-3-methylimidazolium-Based Ionic Liquids Using Molecular Simulation. J Phys Chem B 2016; 120:3360-9. [PMID: 26974037 DOI: 10.1021/acs.jpcb.5b11648] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
During the manufacturing of pharmaceutical compounds, solvent mixtures are commonly used, where the addition of a cosolvent allows for the tuning of the intermolecular interactions present in the system. Here we demonstrate how a similar effect can be accomplished using a room temperature ionic liquid. The pharmaceutical compound acetaminophen is studied in 21 common ionic liquids composed of a 1-n-alkyl-3-methylimidazolium cation with 1 of 7 anions. Using the acetate anion, we predict a large enhancement in solubility of acetaminophen relative to water. We show how this is caused by a synergistic effect of favorable interactions between the ionic liquid and the phenyl, hydroxyl and amide groups of acetaminophen, demonstrating how the ionic liquid cation and anion may be chosen to preferentially solvate different functional groups of complex pharmaceutical compounds. Additionally, while the use of charge scaling in ionic liquid force fields has previously been found to have a minute effect on ionic liquid structural properties, we find it appreciably affects the computed solvation free energy of acetaminophen, which in turn affects the predicted solubility.
Collapse
Affiliation(s)
- Andrew S Paluch
- Department of Chemical, Paper and Biomedical Engineering, Miami University , Oxford, Ohio 45056, United States
| | - Tuanan C Lourenço
- Instituto de Química, Universidade Federal Fluminense , Outeiro de São João Batista, s/n CEP: 24020-141, Niterói, RJ, Brazil
| | - Fenglin Han
- Department of Chemical, Paper and Biomedical Engineering, Miami University , Oxford, Ohio 45056, United States
| | - Luciano T Costa
- Instituto de Química, Universidade Federal Fluminense , Outeiro de São João Batista, s/n CEP: 24020-141, Niterói, RJ, Brazil
| |
Collapse
|
34
|
Mehandzhiyski AY, Grimes BA. Calculation of the probability for ionic association and dissociation reactions by molecular dynamics and umbrella sampling. Mol Phys 2016. [DOI: 10.1080/00268976.2016.1155776] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
| | - Brian A. Grimes
- Department of Chemical Engineering, Norwegian University of Science and Technology, Trondheim, Norway
| |
Collapse
|
35
|
Fuerst GB, Ley RT, Paluch AS. Calculating the Fugacity of Pure, Low Volatile Liquids via Molecular Simulation with Application to Acetanilide, Acetaminophen, and Phenacetin. Ind Eng Chem Res 2015. [DOI: 10.1021/acs.iecr.5b01827] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Georgia B. Fuerst
- Department of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Ryan T. Ley
- Department of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| | - Andrew S. Paluch
- Department of Chemical, Paper
and Biomedical Engineering, Miami University, Oxford, Ohio 45056, United States
| |
Collapse
|
36
|
Paluch AS, Parameswaran S, Liu S, Kolavennu A, Mobley DL. Predicting the excess solubility of acetanilide, acetaminophen, phenacetin, benzocaine, and caffeine in binary water/ethanol mixtures via molecular simulation. J Chem Phys 2015; 142:044508. [PMID: 25637996 PMCID: PMC4312346 DOI: 10.1063/1.4906491] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2014] [Accepted: 01/12/2015] [Indexed: 01/13/2023] Open
Abstract
We present a general framework to predict the excess solubility of small molecular solids (such as pharmaceutical solids) in binary solvents via molecular simulation free energy calculations at infinite dilution with conventional molecular models. The present study used molecular dynamics with the General AMBER Force Field to predict the excess solubility of acetanilide, acetaminophen, phenacetin, benzocaine, and caffeine in binary water/ethanol solvents. The simulations are able to predict the existence of solubility enhancement and the results are in good agreement with available experimental data. The accuracy of the predictions in addition to the generality of the method suggests that molecular simulations may be a valuable design tool for solvent selection in drug development processes.
Collapse
Affiliation(s)
- Andrew S Paluch
- Department of Chemical, Paper, and Biomedical Engineering, Miami University, Oxford, Ohio 45056, USA
| | - Sreeja Parameswaran
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - Shuai Liu
- Department of Pharmaceutical Sciences, University of California Irvine, Irvine, California 92697, USA
| | - Anasuya Kolavennu
- Department of Chemistry, University of New Orleans, New Orleans, Louisiana 70148, USA
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California Irvine, Irvine, California 92697, USA
| |
Collapse
|
37
|
Vasileiadis M, Pantelides CC, Adjiman CS. Prediction of the crystal structures of axitinib, a polymorphic pharmaceutical molecule. Chem Eng Sci 2015. [DOI: 10.1016/j.ces.2014.08.058] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
38
|
Schultz AJ, Kofke DA. Fifth to eleventh virial coefficients of hard spheres. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:023301. [PMID: 25215845 DOI: 10.1103/physreve.90.023301] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2014] [Indexed: 06/03/2023]
Abstract
Virial coefficients B(n) of three-dimensional hard spheres are reported for n=5 to 11, with precision exceeding that presently available in the literature. Calculations are performed using the recursive method due to Wheatley, and a binning approach is proposed to allow more flexibility in where computational effort is directed in the calculations. We highlight the difficulty as a general measure that quantifies performance of an algorithm that computes a stochastic average and show how it can be used as the basis for optimizing such calculations.
Collapse
Affiliation(s)
- Andrew J Schultz
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| | - David A Kofke
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York 14260-4200, USA
| |
Collapse
|
39
|
Boulougouris GC. Free energy calculations, enhanced by a Gaussian ansatz, for the "chemical work" distribution. J Comput Chem 2014; 35:1024-35. [PMID: 24664967 DOI: 10.1002/jcc.23590] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 01/18/2014] [Accepted: 02/26/2014] [Indexed: 01/09/2023]
Abstract
The evaluation of the free energy is essential in molecular simulation because it is intimately related with the existence of multiphase equilibrium. Recently, it was demonstrated that it is possible to evaluate the Helmholtz free energy using a single statistical ensemble along an entire isotherm by accounting for the "chemical work" of transforming each molecule, from an interacting one, to an ideal gas. In this work, we show that it is possible to perform such a free energy perturbation over a liquid vapor phase transition. Furthermore, we investigate the link between a general free energy perturbation scheme and the novel nonequilibrium theories of Crook's and Jarzinsky. We find that for finite systems away from the thermodynamic limit the second law of thermodynamics will always be an inequality for isothermal free energy perturbations, resulting always to a dissipated work that may tend to zero only in the thermodynamic limit. The work, the heat, and the entropy produced during a thermodynamic free energy perturbation can be viewed in the context of the Crooks and Jarzinsky formalism, revealing that for a given value of the ensemble average of the "irreversible" work, the minimum entropy production corresponded to a Gaussian distribution for the histogram of the work. We propose the evaluation of the free energy difference in any free energy perturbation based scheme on the average irreversible "chemical work" minus the dissipated work that can be calculated from the variance of the distribution of the logarithm of the work histogram, within the Gaussian approximation. As a consequence, using the Gaussian ansatz for the distribution of the "chemical work," accurate estimates for the chemical potential and the free energy of the system can be performed using much shorter simulations and avoiding the necessity of sampling the computational costly tails of the "chemical work." For a more general free energy perturbation scheme that the Gaussian ansatz may not be valid, the free energy calculation can be expressed in terms of the moment generating function of the "chemical work" distribution.
Collapse
Affiliation(s)
- Georgios C Boulougouris
- Department of Molecular Biology and Genetics, Democritius University, Alexandroupolis, 68100, Greece
| |
Collapse
|
40
|
A hybrid stochastic-deterministic algorithm for lattice-gas models of catalytic reactions and the computation of TPD spectra. Comput Chem Eng 2014. [DOI: 10.1016/j.compchemeng.2013.08.010] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
41
|
Liu H, Zhang Z, Bara JE, Turner CH. Electrostatic Potential within the Free Volume Space of Imidazole-Based Solvents: Insights into Gas Absorption Selectivity. J Phys Chem B 2013; 118:255-64. [DOI: 10.1021/jp410143j] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Haining Liu
- Department of Chemical & Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| | - Zhongtao Zhang
- Department of Chemical & Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| | - Jason E. Bara
- Department of Chemical & Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| | - C. Heath Turner
- Department of Chemical & Biological Engineering, The University of Alabama, Tuscaloosa, Alabama 35487-0203, United States
| |
Collapse
|
42
|
Ashton DJ, Sánchez-Gil V, Wilding NB. Monte Carlo methods for estimating depletion potentials in highly size-asymmetrical hard sphere mixtures. J Chem Phys 2013; 139:144102. [DOI: 10.1063/1.4824137] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- D J Ashton
- Department of Physics, University of Bath, Bath BA2 7AY, United Kingdom
| | | | | |
Collapse
|
43
|
Moustafa SG, Schultz AJ, Kofke DA. A comparative study of methods to compute the free energy of an ordered assembly by molecular simulation. J Chem Phys 2013; 139:084105. [DOI: 10.1063/1.4818990] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Gai L, Maerzke K, Cummings PT, McCabe C. A Wang-Landau study of a lattice model for lipid bilayer self-assembly. J Chem Phys 2013; 137:144901. [PMID: 23061859 DOI: 10.1063/1.4754536] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The Wang-Landau (WL) Monte Carlo method has been applied to simulate the self-assembly of a lipid bilayer on a 3D lattice. The WL method differs from conventional Monte Carlo methods in that a complete density of states is obtained directly for the system, from which properties, such as the free energy, can be derived. Furthermore, from a single WL simulation, continuous curves of the average energy and heat capacity can be determined, which provide a complete picture of the phase behavior. The lipid model studied consists of 3 or 5 coarse-grained segments on lattices of varying sizes, with the empty lattice sites representing water. A bilayer structure is found to form at low temperatures, with phase transitions to clusters as temperature increases. For 3-segment chains, varying lattice sizes were studied, with the observation that the ratio of chain number to lattice area (i.e., area per lipid) affects the phase transition temperature. At small ratios, only one phase transition occurs between the bilayer and cluster phases, while at high lipid ratios the phase transition occurs in a two-step process with a stable intermediate phase. This second phase transition was not observed in conventional Metropolis Monte Carlo simulations on the same model, demonstrating the advantage of being able to perform a complete scan of the whole temperature range with the WL method. For longer 5-segment chains similar phase transitions are also observed with changes in temperature. In the WL method, due to the extensive nature of the energy, the number of energy bins required to represent the density of states increases as the system size increases and so limits its practical application to larger systems. To improve this, an extension of the WL algorithm, the statistical-temperature Monte Carlo method that allows simulations with larger energy bin sizes, has recently been proposed and is implemented in this work for the 3-segment lattice model. The results obtained are in good agreement with the original WL method and appear to be independent of the energy bin size used.
Collapse
Affiliation(s)
- Lili Gai
- Department of Chemical and Biomolecular Engineering, Vanderbilt University, Nashville, Tennessee 37235, USA
| | | | | | | |
Collapse
|
45
|
Dowdle JR, Buldyrev SV, Stanley HE, Debenedetti PG, Rossky PJ. Temperature and length scale dependence of solvophobic solvation in a single-site water-like liquid. J Chem Phys 2013; 138:064506. [DOI: 10.1063/1.4789981] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
46
|
Hong QJ, van de Walle A. Direct first-principles chemical potential calculations of liquids. J Chem Phys 2013; 137:094114. [PMID: 22957562 DOI: 10.1063/1.4749287] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.
Collapse
Affiliation(s)
- Qi-Jun Hong
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA.
| | | |
Collapse
|
47
|
Schultz AJ, Barlow NS, Chaudhary V, Kofke DA. Mayer Sampling Monte Carlo calculation of virial coefficients on graphics processors. Mol Phys 2012. [DOI: 10.1080/00268976.2012.730642] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
48
|
|
49
|
Luzhkov VB. Free-energy molecular simulations of the inclusion complex of Ne with fullerene C 60in water. MOLECULAR SIMULATION 2012. [DOI: 10.1080/08927022.2011.632417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
50
|
Eshet H, Khaliullin RZ, Kühne TD, Behler J, Parrinello M. Microscopic origins of the anomalous melting behavior of sodium under high pressure. PHYSICAL REVIEW LETTERS 2012; 108:115701. [PMID: 22540486 DOI: 10.1103/physrevlett.108.115701] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/31/2011] [Indexed: 05/31/2023]
Abstract
X-ray diffraction experiments have shown that sodium exhibits a dramatic pressure-induced drop in melting temperature, which extends from 1000 K at ~30 GPa to as low as room temperature at ~120 GPa. Despite significant theoretical effort to understand the anomalous melting, its origins are still debated. In this work, we reconstruct the sodium phase diagram by using an ab initio quality neural-network potential. Furthermore, we demonstrate that the reentrant behavior results from the screening of interionic interactions by conduction electrons, which at high pressure induces a softening in the short-range repulsion.
Collapse
Affiliation(s)
- Hagai Eshet
- Department of Chemistry and Applied Biosciences, ETH Zürich, Zürich, Switzerland.
| | | | | | | | | |
Collapse
|