skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods

Abstract

In this work, we present a study of the ability of different semiempirical methods to describe intermolecular interactions in water solution. In particular, we focus on methods based on the Neglect of Diatomic Differential Overlap approximation. Significant improvements of these methods have been reported in the literature in the past years regarding the description of non-covalent interactions. In particular, a broad range of methodologies has been developed to deal with the properties of hydrogen-bonded systems, with varying degrees of success. In contrast, the interactions between water and a molecule containing hydrophobic groups have been little analyzed. Indeed, by considering the potential energy surfaces obtained using different semiempirical Hamiltonians for the intermolecular interactions of model systems, we found that none of the available methods provides an entirely satisfactory description of both hydrophobic and hydrophilic interactions in water. In addition, a vibrational analysis carried out in a model system for these interactions, a methane clathrate cluster, showed that some recent methods cannot be used to carry out studies of vibrational properties. Following a procedure established in our group [M. I. Bernal-Uruchurtu, M. T. C. Martins-Costa, C. Millot, and M. F. Ruiz-López, J. Comput. Chem. 21, 572 (2000); W. Harb, M. I. Bernal-Uruchurtu,more » and M. F. Ruiz-López, Theor. Chem. Acc. 112, 204 (2004)], we developed new parameters for the core-core interaction terms based on fitting potential energy curves obtained at the MP2 level for our model system. We investigated the transferability of the new parameters to describe a system, having both hydrophilic and hydrophobic groups, interacting with water. We found that only by introducing two different sets of parameters for hydrophilic and hydrophobic hydrogen atom types we are able to match the features of the ab initio calculated properties. Once this assumption is made, a good agreement with the MP2 reference is achieved. The results reported in this work provide therefore a direction for future developments of semiempirical approaches that are still required to investigate chemical processes in biomolecules and in large disordered systems.« less

Authors:
; ; ;  [1];  [2]
  1. Université de Lorraine, SRSMC UMR 7565, Vandoeuvre-lès-Nancy F-54506 (France)
  2. (France)
Publication Date:
OSTI Identifier:
22419878
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Chemical Physics; Journal Volume: 141; Journal Issue: 3; Other Information: (c) 2014 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
74 ATOMIC AND MOLECULAR PHYSICS; APPROXIMATIONS; ATOMS; CLATHRATES; HAMILTONIANS; INTERACTIONS; METHANE; MOLECULAR ORBITAL METHOD; MOLECULES; POTENTIAL ENERGY; WATER

Citation Formats

Marion, Antoine, Monard, Gérald, Ruiz-López, Manuel F., E-mail: Manuel.Ruiz@univ-lorraine.fr, Ingrosso, Francesca, E-mail: Francesca.Ingrosso@univ-lorraine.fr, and CNRS, SRSMC UMR 7565, Vandoeuvre-lès-Nancy F-54506. Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods. United States: N. p., 2014. Web. doi:10.1063/1.4886655.
Marion, Antoine, Monard, Gérald, Ruiz-López, Manuel F., E-mail: Manuel.Ruiz@univ-lorraine.fr, Ingrosso, Francesca, E-mail: Francesca.Ingrosso@univ-lorraine.fr, & CNRS, SRSMC UMR 7565, Vandoeuvre-lès-Nancy F-54506. Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods. United States. doi:10.1063/1.4886655.
Marion, Antoine, Monard, Gérald, Ruiz-López, Manuel F., E-mail: Manuel.Ruiz@univ-lorraine.fr, Ingrosso, Francesca, E-mail: Francesca.Ingrosso@univ-lorraine.fr, and CNRS, SRSMC UMR 7565, Vandoeuvre-lès-Nancy F-54506. Mon . "Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods". United States. doi:10.1063/1.4886655.
@article{osti_22419878,
title = {Water interactions with hydrophobic groups: Assessment and recalibration of semiempirical molecular orbital methods},
author = {Marion, Antoine and Monard, Gérald and Ruiz-López, Manuel F., E-mail: Manuel.Ruiz@univ-lorraine.fr and Ingrosso, Francesca, E-mail: Francesca.Ingrosso@univ-lorraine.fr and CNRS, SRSMC UMR 7565, Vandoeuvre-lès-Nancy F-54506},
abstractNote = {In this work, we present a study of the ability of different semiempirical methods to describe intermolecular interactions in water solution. In particular, we focus on methods based on the Neglect of Diatomic Differential Overlap approximation. Significant improvements of these methods have been reported in the literature in the past years regarding the description of non-covalent interactions. In particular, a broad range of methodologies has been developed to deal with the properties of hydrogen-bonded systems, with varying degrees of success. In contrast, the interactions between water and a molecule containing hydrophobic groups have been little analyzed. Indeed, by considering the potential energy surfaces obtained using different semiempirical Hamiltonians for the intermolecular interactions of model systems, we found that none of the available methods provides an entirely satisfactory description of both hydrophobic and hydrophilic interactions in water. In addition, a vibrational analysis carried out in a model system for these interactions, a methane clathrate cluster, showed that some recent methods cannot be used to carry out studies of vibrational properties. Following a procedure established in our group [M. I. Bernal-Uruchurtu, M. T. C. Martins-Costa, C. Millot, and M. F. Ruiz-López, J. Comput. Chem. 21, 572 (2000); W. Harb, M. I. Bernal-Uruchurtu, and M. F. Ruiz-López, Theor. Chem. Acc. 112, 204 (2004)], we developed new parameters for the core-core interaction terms based on fitting potential energy curves obtained at the MP2 level for our model system. We investigated the transferability of the new parameters to describe a system, having both hydrophilic and hydrophobic groups, interacting with water. We found that only by introducing two different sets of parameters for hydrophilic and hydrophobic hydrogen atom types we are able to match the features of the ab initio calculated properties. Once this assumption is made, a good agreement with the MP2 reference is achieved. The results reported in this work provide therefore a direction for future developments of semiempirical approaches that are still required to investigate chemical processes in biomolecules and in large disordered systems.},
doi = {10.1063/1.4886655},
journal = {Journal of Chemical Physics},
number = 3,
volume = 141,
place = {United States},
year = {Mon Jul 21 00:00:00 EDT 2014},
month = {Mon Jul 21 00:00:00 EDT 2014}
}
  • Molecular recognition, via non-covalent processes such as hydrogen bonding, {pi}-{pi}, and hydrophobic interactions, is an important biological phenomenon for guests, such as drugs, proteins, and other important biological molecules with, for example, host DNA/RNA. We have studied a novel molecular recognition process using guests that encompass aromatic and aliphatic amino acids [L-alanine, L-glutamine (L-Gln), L-histidine, L-isoleucine(L-Ile), L-leucine(L-Leu), L-phenylalanine(L-Phe), L-proline, L-tryptophan(L-Trp), L-valine(L-Val)], substituted aromatic carboxylic acids o-, m-, p-aminobenzoic acids (G1-3), benzoic acid (G4), phenylacetic acid (G5), p-methoxyphenylacetic acid (G6), o-methyoxybenozoic acid (G9), o-nitrobenzoic acid (G10), and aliphatic carboxylic acids [cyclohexylacetic acid (G7), 1-adamantanecarboxylic acid (G8)] with supramolecular, bioorganometallic hosts, ({eta}{supmore » 5}-pentamethylcyclopentadienyl)rhodium (Cp{sup *}Rh)-nucleobase, nucleoside, and nucleotide cyclic trimer complexes in aqueous solution at pH 7, utilizing {sup 1}H NMR, NOE, and molecular modeling techniques, and, as well, determining association constants (K{sub a}) and free energies of complexation ({Delta}{degree}G). The host-guest complexation occurs predominantly via non-covalent {pi}-{pi}, hydrophobic, and possible subtle H-bonding interactions, with steric, electronic, and molecular conformational parameters as important criteria. 8 refs., 6 figs., 3 tabs.« less
  • The authors present experimental and computational determination of the excited-state properties in several perylene dyestuffs that are potential candidates for use as laser dyes. Attention is focused on the following species derived from 3,4,9,10-perylenetetracarboxylic acid dianhydride: the bis ((2,6-dimethylphenyl)imide) (1b, DXP); the bis(methylimide) (1c, DMP); both 1,6,7,12-tetrachloro- (1d, Cl{sub 4}DMP) and 1,2,5,6,7,8,11,12-octachloro- (1e, Cl{sub 8}DMP) derivatives. Soluble derivatives of seven-ringed or larger aromatic systems are produced by the introduction of relatively rigid out-of-plane substituents that prevent intermolecular close packing. Chiral distortions in ring-chlorinated perylene derivatives significantly alter the shape of the absorption bands, a consequence of symmetry breaking. Triplet-triplet absorptionsmore » and oxygen-quenching rates are observed under the conditions found in a dye laser cavity. Semiempirical molecular orbital calculations (INDO/S) provide detailed mapping of the singlet and triplet excited state manifolds of DMP, Cl{sub 4}DMP, and Cl{sub 8}DMP. Computed transition energies and intensities are used in the interpretation of the spectral features, in particular the observed T{sub 1} {yields} T{sub n} and potential S{sub 1} {yields} S{sub n} absorptions. The authors conclude the perylene-3,4,9,10-tetracarboxylic acid diimide chromophores may be solubilized and utilized in laser materials exhibiting superior performance in terms of power output, tuning range, and light stability. 42 refs., 5 figs., 7 tabs.« less
  • The analytic first derivative with respect to nuclear coordinates is formulated and implemented in the framework of the three-body fragment molecular orbital (FMO) method. The gradient has been derived and implemented for restricted Hartree-Fock, second-order Møller-Plesset perturbation, and density functional theories. The importance of the three-body fully analytic gradient is illustrated through the failure of the two-body FMO method during molecular dynamics simulations of a small water cluster. The parallel implementation of the fragment molecular orbital method, its parallel efficiency, and its scalability on the Blue Gene/Q architecture up to 262,144 CPU cores, are also discussed.
  • Deuterium quadrupole coupling constants (DQCC) for deuterium nuclei in a variety of bonding situations have been obtained by means of semiempirical molecular orbital wavefunctions in the INDO approximation. All integrals of the operator entering the electronic contributions have been included with no approximations in their evaluation. Analytical expressions are tabulated for the relevant two-center integrals; three-center integrals were obtained by numerical integration techniques. The factors responsible for the very strong dependence of deuterium QCC on internuclear distance, and the relative constancy of the parameters for a given bonding situation are discussed with regard to the types of terms which entermore » the integral expressions and the molecular electronic distributions. Some interesting and possibly useful predictions are also included in the tabulated results; the DQCC for the bridging deuterium in diborane is predicted to have a value of -110 kHz, which may be the first negative value for a DQCC; the effects of hydrogen bonding can be significant, the result of lengthening of the D--X bond on hydrogen bond formation.« less
  • PM3 and AM1 semiempirical molecular orbital techniques are used to establish a model for internal rotation about the N{sub 2}-N{sub 3} axis of 1-aryl-3,3-dialkyltriazines. The PM3 method is satisfactory for obtaining agreement between the experimental and calculated results, but the AM1 method has an artifact in the potential energy curve of internal rotation about the N{sub 2}-N{sub 3} axis. 24 refs., 6 figs., 5 tabs.