1
|
Yamasaki K, Tsuzuki S, Tateno H. Stabilization of the Protein Structure by the Many-Body Cooperative Effect in the NH/π Hydrogen-bonding Tryptophan Triad. J Phys Chem B 2024; 128:7401-7406. [PMID: 39018377 DOI: 10.1021/acs.jpcb.4c02391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/19/2024]
Abstract
The indole ring of tryptophan can form NH/π hydrogen bonds, acting both as a hydrogen donor at the NH group in the pyrrole subring and as a hydrogen acceptor at the benzene subring. In the structural core of the trimeric stable protein Pholiota squarrosa lectin (PhoSL), three indoles are symmetrically arranged and form NH/π hydrogen bonds among each other. Here, we conducted quantum chemical calculations on this indole triad by using various methods and basis sets. The analyses revealed cooperativity in triad formation, with the many-body effect contributing approximately -2 kcal mol-1, which significantly stabilizes this protein. Symmetry-adapted perturbation theory ascribed this effect to the induced polarization. The electrostatic potential and atomic charges indeed revealed a charge redistribution through the NH/π hydrogen bond, which was favorable for triad formation.
Collapse
Affiliation(s)
- Kazuhiko Yamasaki
- Biomedical Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Higashi, Tsukuba 305-8566, Japan
| | - Seiji Tsuzuki
- Department of Applied Physics, The University of Tokyo, 7-3-1 Hongo, Tokyo 113-8656, Japan
| | - Hiroaki Tateno
- Cellular and Molecular Biotechnology Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Higashi, Tsukuba 305-8566, Japan
| |
Collapse
|
2
|
Santis GD, Herman KM, Heindel JP, Xantheas SS. Descriptors of water aggregation. J Chem Phys 2024; 160:054306. [PMID: 38341703 DOI: 10.1063/5.0179815] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2023] [Accepted: 01/05/2024] [Indexed: 02/13/2024] Open
Abstract
We rely on a total of 23 (cluster size, 8 structural, and 14 connectivity) descriptors to investigate structural patterns and connectivity motifs associated with water cluster aggregation. In addition to the cluster size n (number of molecules), the 8 structural descriptors can be further categorized into (i) one-body (intramolecular): covalent OH bond length (rOH) and HOH bond angle (θHOH), (ii) two-body: OO distance (rOO), OHO angle (θOHO), and HOOX dihedral angle (ϕHOOX), where X lies on the bisector of the HOH angle, (iii) three-body: OOO angle (θOOO), and (iv) many-body: modified tetrahedral order parameter (q) to account for two-, three-, four-, five-coordinated molecules (qm, m = 2, 3, 4, 5) and radius of gyration (Rg). The 14 connectivity descriptors are all many-body in nature and consist of the AD, AAD, ADD, AADD, AAAD, AAADD adjacencies [number of hydrogen bonds accepted (A) and donated (D) by each water molecule], Wiener index, Average Shortest Path Length, hydrogen bond saturation (% HB), and number of non-short-circuited three-membered cycles, four-membered cycles, five-membered cycles, six-membered cycles, and seven-membered cycles. We mined a previously reported database of 4 948 959 water cluster minima for (H2O)n, n = 3-25 to analyze the evolution and correlation of these descriptors for the clusters within 5 kcal/mol of the putative minima. It was found that rOH and % HB correlated strongly with cluster size n, which was identified as the strongest predictor of energetic stability. Marked changes in the adjacencies and cycle count were observed, lending insight into changes in the hydrogen bond network upon aggregation. A Principal Component Analysis (PCA) was employed to identify descriptor dependencies and group clusters into specific structural patterns across different cluster sizes. The results of this study inform our understanding of how water clusters evolve in size and what appropriate descriptors of their structural and connectivity patterns are with respect to system size, stability, and similarity. The approach described in this study is general and can be easily extended to other hydrogen-bonded systems.
Collapse
Affiliation(s)
- Garrett D Santis
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Joseph P Heindel
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
- Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MSIN J7-10, Richland, Washington 99352, USA
- Computational and Theoretical Chemistry Institute (CTCI), Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| |
Collapse
|
3
|
Li XL, Li CM, Zhu JY, Zhou Z, Hao Q, Wang CS. A scheme for rapid evaluation of the intermolecular three-body polarization effect in water clusters. J Comput Chem 2023; 44:677-686. [PMID: 36408852 DOI: 10.1002/jcc.27032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Revised: 09/22/2022] [Accepted: 10/24/2022] [Indexed: 11/22/2022]
Abstract
The ability to accurately and rapidly evaluate the intermolecular many-body polarization effect of the water system is very important for computer simulations of biomolecule in aqueous. In this paper, a scheme is proposed based on the polarizable dipole-dipole interaction model and used to rapidly estimate the intermolecular many-body polarization effect in water clusters. We use a bond-dipole-based polarization function to evaluate the polarization energy. We regard two OH bonds of a water molecule as two bond-dipoles and set the permanent OH bond-dipole moment of a water molecule to be 1.51 Debye. We estimate the induced OH bond-dipole moment via a simple formula in which only one correction factor is needed. This scheme is then applied to tens of water clusters to calculate the three- and four-body interaction energies. The three-body interaction energies of 93 water clusters produced by our scheme are compared with those produced by the counterpoise-corrected CCSD(T)/aug-cc-pVDZ, MP2/aug-cc-pVDZ, M06-2X/jul-cc-pVTZ methods, by the AMOEBApro13, iAMOEBA, AMOEBA+, AMOEBA+(CF) methods, and by the MB-pol method. The four-body interaction energies of 47 water clusters yielded by our scheme are compared with those yielded by the counterpoise-corrected MP2/aug-cc-pVDZ and M06-2X/ jul-cc-pVTZ methods, by the AMOEBApro13, AMOEBA+, AMOEBA+(CF) methods, and by the MB-pol method. The comparison results show that the scheme proposed in this paper can reproduce the counterpoise-corrected CCSD(T)/aug-cc-pVDZ three-body interaction energies and reproduce the counterpoise-corrected MP2/aug-cc-pVDZ four-body interaction energies both accurately and efficiently. We anticipate the scheme proposed here can be useful for computer simulations of liquid water and aqueous solutions.
Collapse
Affiliation(s)
- Xiao-Lei Li
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| | - Chao-Ming Li
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| | - Jia-Yi Zhu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| | - Zhan Zhou
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| | - Qiang Hao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| | - Chang-Sheng Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, People's Republic of China
| |
Collapse
|
4
|
Perez FP, Bandeira JP, Perez Chumbiauca CN, Lahiri DK, Morisaki J, Rizkalla M. Multidimensional insights into the repeated electromagnetic field stimulation and biosystems interaction in aging and age-related diseases. J Biomed Sci 2022; 29:39. [PMID: 35698225 PMCID: PMC9190166 DOI: 10.1186/s12929-022-00825-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Accepted: 06/07/2022] [Indexed: 11/25/2022] Open
Abstract
We provide a multidimensional sequence of events that describe the electromagnetic field (EMF) stimulation and biological system interaction. We describe this process from the quantum to the molecular, cellular, and organismal levels. We hypothesized that the sequence of events of these interactions starts with the oscillatory effect of the repeated electromagnetic stimulation (REMFS). These oscillations affect the interfacial water of an RNA causing changes at the quantum and molecular levels that release protons by quantum tunneling. Then protonation of RNA produces conformational changes that allow it to bind and activate Heat Shock Transcription Factor 1 (HSF1). Activated HSF1 binds to the DNA expressing chaperones that help regulate autophagy and degradation of abnormal proteins. This action helps to prevent and treat diseases such as Alzheimer's and Parkinson's disease (PD) by increasing clearance of pathologic proteins. This framework is based on multiple mathematical models, computer simulations, biophysical experiments, and cellular and animal studies. Results of the literature review and our research point towards the capacity of REMFS to manipulate various networks altered in aging (Reale et al. PloS one 9, e104973, 2014), including delay of cellular senescence (Perez et al. 2008, Exp Gerontol 43, 307-316) and reduction in levels of amyloid-β peptides (Aβ) (Perez et al. 2021, Sci Rep 11, 621). Results of these experiments using REMFS at low frequencies can be applied to the treatment of patients with age-related diseases. The use of EMF as a non-invasive therapeutic modality for Alzheimer's disease, specifically, holds promise. It is also necessary to consider the complicated and interconnected genetic and epigenetic effects of the REMFS-biological system's interaction while avoiding any possible adverse effects.
Collapse
Affiliation(s)
- Felipe P Perez
- Indiana University School of Medicine, Indianapolis, IN, USA.
- Division of General Internal Medicine and Geriatrics, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, USA.
| | - Joseph P Bandeira
- Indiana University School of Medicine, Indianapolis, IN, USA
- Division of General Internal Medicine and Geriatrics, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, USA
| | - Cristina N Perez Chumbiauca
- Indiana University School of Medicine, Indianapolis, IN, USA
- Division of Rheumatology, Department of Medicine, Indiana University School of Medicine, Indianapolis, IN, USA
| | - Debomoy K Lahiri
- Indiana University School of Medicine, Indianapolis, IN, USA
- Department of Psychiatry, Institute of Psychiatric Research, Neuroscience Research Center, Indiana University School of Medicine, Indianapolis, IN, USA
- Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, USA
| | - Jorge Morisaki
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Maher Rizkalla
- Department of Electrical and Computer Engineering, Indiana University-Purdue University, Indianapolis, IN, USA
| |
Collapse
|
5
|
Fedorov DG. Three-Body Energy Decomposition Analysis Based on the Fragment Molecular Orbital Method. J Phys Chem A 2020; 124:4956-4971. [DOI: 10.1021/acs.jpca.0c03085] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
6
|
Combining Structural Data with Computational Methodologies to Investigate Structure-Function Relationships in TRP Channels. Methods Mol Biol 2019; 1987:65-82. [PMID: 31028674 DOI: 10.1007/978-1-4939-9446-5_5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/08/2022]
Abstract
Since the emergence of high-resolution three-dimensional structures of membrane proteins, and the increasing availability of state-of-the-art algorithms and high-performance-computing facilities, classical molecular dynamics (MD) simulations have become a routine device to explore the molecular behavior of these proteins. The rise of cryo-electron microscopy (cryo-EM) as a credible experimental tool to resolve structures at an atomic level has revolutionized structural biology in recent years, culminating in the disclosure of the first high-resolution three-dimensional structure of a transient receptor potential (TRP) channel, the vanilloid receptor 1 (TRPV1). As a result, the number of research articles investigating the molecular behavior of TRP channels using macromolecular simulation techniques has proliferated. This review provides an overview of the current state of this field, including our understanding of TRP channel structure, the framework of classical MD simulations, and how to perform such simulations to investigate structure-function relationships in TRP channels.
Collapse
|
7
|
Heßelmann A. Correlation effects and many-body interactions in water clusters. Beilstein J Org Chem 2018; 14:979-991. [PMID: 29977369 PMCID: PMC6009095 DOI: 10.3762/bjoc.14.83] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2018] [Accepted: 04/11/2018] [Indexed: 11/23/2022] Open
Abstract
Background: The quantum-chemical description of the interactions in water clusters is an essential basis for deriving accurate and physically sound models of the interaction potential for water to be used in molecular simulations. In particular, the role of many-body interactions beyond the two-body interactions, which are often not explicitly taken into account by empirical force fields, can be accurately described by quantum chemistry methods on an adequate level, e.g., random-phase approximation electron correlation methods. The relative magnitudes of the different interaction energy contributions obtained by accurate ab initio calculations can therefore provide useful insights that can be exploited to develop enhanced force field methods. Results: In line with earlier theoretical studies of the interactions in water clusters, it has been found that the main contribution to the many-body interactions in clusters with a size of up to N = 13 molecules are higher-order polarisation interaction terms. Compared to this, many-body dispersion interactions are practically negligible for all studied sytems. The two-body dispersion interaction, however, plays a significant role in the formation of the structures of the water clusters and their stability, since it leads to a distinct compression of the cluster sizes compared to the structures optimized on an uncorrelated level. Overall, the many-body interactions amount to about 13% of the total interaction energy, irrespective of the cluster size. The electron correlation contribution to these, however, amounts to only about 30% to the total many-body interactions for the largest clusters studied and is repulsive for all structures considered in this work. Conclusion: While this shows that three- and higher-body interactions can not be neglected in the description of water complexes, the electron correlation contributions to these are much smaller in comparison to the two-body electron correlation effects. Efficient quantum chemistry approaches for describing intermolecular interactions between water molecules may therefore describe higher-body interactions on an uncorrelated Hartree-Fock level without a serious loss in accuracy.
Collapse
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, 91058 Erlangen, Germany
| |
Collapse
|
8
|
Abstract
Hydrogen bond directionality in the water dimer is explained on the basis of symmetry-adapted intermolecular perturbation theory which directly separates the intermolecular interaction energy into four physically interpretable components: electrostatics, exchange-repulsion, dispersion, and induction. Analysis of these four main contributions to the binding energy allows a deeper understanding of the dominant factors ruling the mutual arrangement of the two monomers. A preference for the linear configuration is shown to be due to a subtle interplay of all four energy components. While the first-order terms, electrostatic and exchange-repulsion, almost perfectly cancel each other near the equilibrium geometry of the dimer, the importance of the second- and higher-order terms, induction and dispersion, becomes evident.
Collapse
Affiliation(s)
- Maxim Tafipolsky
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg, Campus Hubland Nord , Emil-Fischer-Strasse 42, D-97074 Würzburg, Germany
| |
Collapse
|
9
|
Tan SYS, Izgorodina EI. Comparison of the Effective Fragment Potential Method with Symmetry-Adapted Perturbation Theory in the Calculation of Intermolecular Energies for Ionic Liquids. J Chem Theory Comput 2016; 12:2553-68. [DOI: 10.1021/acs.jctc.6b00141] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Samuel Y. S. Tan
- School of Chemistry, Monash University, 17 Rainforest
Walk, Clayton, Victoria 3800, Australia
| | - Ekaterina I. Izgorodina
- School of Chemistry, Monash University, 17 Rainforest
Walk, Clayton, Victoria 3800, Australia
| |
Collapse
|
10
|
Misquitta AJ, Stone AJ. Accurate Induction Energies for Small Organic Molecules: 1. Theory. J Chem Theory Comput 2015; 4:7-18. [PMID: 26619975 DOI: 10.1021/ct700104t] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The induction energy often plays a very important role in determining the structure and properties of clusters of organic molecules, but only in recent years has an effort been made to include this energy in such calculations, notably in the field of organic crystal structure prediction. In this paper and the following one in this issue we provide ab initio methods suitable for the accurate inclusion of the induction energy for molecules containing as many as 30 atoms or so. These techniques are based on Symmetry-Adapted Perturbation Theory using Density Functional Theory [SAPT(DFT)] and use distributed polarizabilities computed using the recently developed density-fitting algorithm with constrained refinement. With this approach we are able to obtain induction models of varying complexity and study the effects of overlap and related numerical issues. Basis set effects on the exact and asymptotic induction energies are investigated, and the roles of higher-order induction energies and many-body effects are explored.
Collapse
Affiliation(s)
- Alston J Misquitta
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| | - Anthony J Stone
- University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, U.K., and University College London, 20 Gordon Street, London WC1H 0AJ, U.K
| |
Collapse
|
11
|
Saha S, Sastry GN. Quantifying cooperativity in water clusters: an attempt towards obtaining a generalised equation. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1072648] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
12
|
Raston PL, Jäger W. Microwave spectroscopy of the seeded binary and ternary clusters CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)(N=1,2). J Chem Phys 2015; 142:144308. [PMID: 25877579 DOI: 10.1063/1.4917420] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We report the Fourier transform microwave spectra of the a-type J = 1-0 transitions of the binary and ternary CO-(pH2)2, CO-pH2-He, CO-HD, and CO-(oD2)N=1,2 clusters. In addition to the normal isotopologue of CO for all clusters, we observed the transitions of the minor isotopologues, (13)C(16)O, (12)C(18)O, and (13)C(18)O, for CO-(pH2)2 and CO-pH2-He. All transitions lie within 335 MHz of the experimentally or theoretically predicted values. In comparison to previously reported infrared spectra [Moroni et al., J. Chem. Phys. 122, 094314 (2005)], we are able to tentatively determine the vibrational shift for CO-pH2-He, in addition to its b-type J = 1-0 transition frequency. The a-type frequency of CO-pH2-He is similar to that of CO-He2 [Surin et al., Phys. Rev. Lett. 101, 233401 (2008)], suggesting that the pH2 molecule has a strong localizing effect on the He density. Perturbation theory analysis of CO-oD2 reveals that it is approximately T-shaped, with an anisotropy of the intermolecular potential amounting to ∼9 cm(-1).
Collapse
Affiliation(s)
- Paul L Raston
- Department of Chemistry, University of Adelaide, SA 5005, Australia
| | - Wolfgang Jäger
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| |
Collapse
|
13
|
Temelso B, Renner CR, Shields GC. Importance and Reliability of Small Basis Set CCSD(T) Corrections to MP2 Binding and Relative Energies of Water Clusters. J Chem Theory Comput 2015; 11:1439-48. [DOI: 10.1021/ct500944v] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Berhane Temelso
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - Carla R. Renner
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| | - George C. Shields
- Dean’s
Office, College of Arts and Sciences, and Department of Chemistry, Bucknell University, Lewisburg, Pennsylvania 17837, United States
| |
Collapse
|
14
|
Theoretical investigation of hydrogen bonding interaction in H3O+(H2O)9 complex. J Mol Model 2014; 20:2480. [DOI: 10.1007/s00894-014-2480-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2014] [Accepted: 09/28/2014] [Indexed: 11/24/2022]
|
15
|
Ji C, Mei Y. Some practical approaches to treating electrostatic polarization of proteins. Acc Chem Res 2014; 47:2795-803. [PMID: 24883956 DOI: 10.1021/ar500094n] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Conspectus Electrostatic interaction plays a significant role in determining many properties of biomolecules, which exist and function in aqueous solution, a highly polar environment. For example, proteins are composed of amino acids with charged, polar, and nonpolar side chains and their specific electrostatic properties are fundamental to the structure and function of proteins. An important issue that arises in computational study of biomolecular interaction and dynamics based on classical force field is lack of polarization. Polarization is a phenomenon in which the charge distribution of an isolated molecule will be distorted when interacting with another molecule or presented in an external electric field. The distortion of charge distribution is intended to lower the overall energy of the molecular system, which is counter balanced by the increased internal energy of individual molecules due to the distorted charge distributions. The amount of the charge redistribution, which characterizes the polarizability of a molecule, is determined by the level of the charge distortion. Polarization is inherently quantum mechanical, and therefore classical force fields with fixed atomic charges are incapable of capturing this important effect. As a result, simulation studies based on popular force fields, AMBER, CHARMM, etc., lack the polarization effect, which is a widely known deficiency in most computational studies of biomolecules today. Many efforts have been devoted to remedy this deficiency, such as adding additional movable charge on the atom, allowing atomic charges to fluctuate, or including induced multipoles. Although various successes have been achieved and progress at various levels has been reported over the past decades, the issue of lacking polarization in force field based simulations is far from over. For example, some of these methods do not always give converged results, and other methods require huge computational cost. This Account reviews recent work on developing polarized and polarizable force fields based on fragment quantum mechanical calculations for proteins. The methods described here are based on quantum mechanical calculations of proteins in solution, but with a different level of rigor and different computational efficiency for the molecular dynamics applications. In the general approach, a fragment quantum mechanical calculation for protein with implicit solvation is carried out to derive a polarized protein-specific charge (PPC) for any given protein structure. The PPC correctly reflects the polarization state of the protein in a given conformation, and it can also be dynamically changed as the protein changes conformation in dynamics simulations. Another approach that is computationally more efficient is the effective polarizable bond method in which only polar bonds or groups can be polarized and their polarizabilities are predetermined from quantum mechanical calculations of these groups in external electric fields. Both methods can be employed for applications in various situations by taking advantage of their unique features.
Collapse
Affiliation(s)
- Changge Ji
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
- Institute for Advanced Interdisciplinary Research, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
16
|
Fedorov DG, Asada N, Nakanishi I, Kitaura K. The use of many-body expansions and geometry optimizations in fragment-based methods. Acc Chem Res 2014; 47:2846-56. [PMID: 25144610 DOI: 10.1021/ar500224r] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Conspectus Chemists routinely work with complex molecular systems: solutions, biochemical molecules, and amorphous and composite materials provide some typical examples. The questions one often asks are what are the driving forces for a chemical phenomenon? How reasonable are our views of chemical systems in terms of subunits, such as functional groups and individual molecules? How can one quantify the difference in physicochemical properties of functional units found in a different chemical environment? Are various effects on functional units in molecular systems additive? Can they be represented by pairwise potentials? Are there effects that cannot be represented in a simple picture of pairwise interactions? How can we obtain quantitative values for these effects? Many of these questions can be formulated in the language of many-body effects. They quantify the properties of subunits (fragments), referred to as one-body properties, pairwise interactions (two-body properties), couplings of two-body interactions described by three-body properties, and so on. By introducing the notion of fragments in the framework of quantum chemistry, one obtains two immense benefits: (a) chemists can finally relate to quantum chemistry, which now speaks their language, by discussing chemically interesting subunits and their interactions and (b) calculations become much faster due to a reduced computational scaling. For instance, the somewhat academic sounding question of the importance of three-body effects in water clusters is actually another way of asking how two hydrogen bonds affect each other, when they involve three water molecules. One aspect of this is the many-body charge transfer (CT), because the charge transfers in the two hydrogen bonds are coupled to each other (not independent). In this work, we provide a generalized view on the use of many-body expansions in fragment-based methods, focusing on the general aspects of the property expansion and a contraction of a many-body expansion in a formally two-body series, as exemplified in the development of the fragment molecular orbital (FMO) method. Fragment-based methods have been very successful in delivering the properties of fragments, as well as the fragment interactions, providing insights into complex chemical processes in large molecular systems. We briefly review geometry optimizations performed with fragment-based methods and present an efficient geometry optimization method based on the combination of FMO with molecular mechanics (MM), applied to the complex of a subunit of protein kinase 2 (CK2) with a ligand. FMO results are discussed in comparison with experimental and MM-optimized structures.
Collapse
Affiliation(s)
- Dmitri G. Fedorov
- NRI, National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Naoya Asada
- Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Isao Nakanishi
- Department of Pharmaceutical Sciences, Kinki University, 3-4-1,
Kowakae, Higashi-Osaka, Osaka 577-8502, Japan
| | - Kazuo Kitaura
- Graduate
School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan
| |
Collapse
|
17
|
Ouyang JF, Cvitkovic MW, Bettens RPA. Trouble with the Many-Body Expansion. J Chem Theory Comput 2014; 10:3699-707. [DOI: 10.1021/ct500396b] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- John F. Ouyang
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| | - Milan W. Cvitkovic
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| | - Ryan P. A. Bettens
- Department of Chemistry, National University of Singapore, 3 Science
Drive 3, Singapore 117543
| |
Collapse
|
18
|
Jackson NE, Chen LX, Ratner MA. Solubility of Nonelectrolytes: A First-Principles Computational Approach. J Phys Chem B 2014; 118:5194-202. [DOI: 10.1021/jp5024197] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Nicholas E. Jackson
- Department
of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Lin X. Chen
- Department
of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Mark A. Ratner
- Department
of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| |
Collapse
|
19
|
|
20
|
Semerci F, Yeşilel OZ, Ölmez H, Büyükgüngör O. Supramolecular assemblies of copper(II)–pyridine-2,3-dicarboxylate complexes with N-donor ligands and clustered water molecules. Inorganica Chim Acta 2014. [DOI: 10.1016/j.ica.2013.10.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
21
|
Guevara-Vela JM, Chávez-Calvillo R, García-Revilla M, Hernández-Trujillo J, Christiansen O, Francisco E, Martín Pendás Á, Rocha-Rinza T. Hydrogen-Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach. Chemistry 2013; 19:14304-15. [DOI: 10.1002/chem.201300656] [Citation(s) in RCA: 73] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2013] [Revised: 05/11/2013] [Indexed: 11/08/2022]
|
22
|
|
23
|
Salmi T, Sälli E, Halonen L. A Nine-Dimensional Calculation of the Vibrational OH Stretching and HOH Bending Spectrum of the Water Trimer. J Phys Chem A 2012; 116:5368-74. [DOI: 10.1021/jp3017584] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Teemu Salmi
- Laboratory
of Physical Chemistry, Department of Chemistry,
A.I. Virtasen aukio 1 (P.O. BOX 55), FI-00014, University of Helsinki, Finland
| | - Elina Sälli
- Laboratory
of Physical Chemistry, Department of Chemistry,
A.I. Virtasen aukio 1 (P.O. BOX 55), FI-00014, University of Helsinki, Finland
| | - Lauri Halonen
- Laboratory
of Physical Chemistry, Department of Chemistry,
A.I. Virtasen aukio 1 (P.O. BOX 55), FI-00014, University of Helsinki, Finland
| |
Collapse
|
24
|
Li QZ, Li R, Guo P, Li H, Li WZ, Cheng JB. Competition of chalcogen bond, halogen bond, and hydrogen bond in SCSHOX and SeCSeHOX (X=Cl and Br) complexes. COMPUT THEOR CHEM 2012. [DOI: 10.1016/j.comptc.2011.11.019] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
25
|
Kharat B, Deshmukh V, Chaudhari A. Hydrogen-bonding interactions in acetonitrile oligomers using density functional theory method. Struct Chem 2011. [DOI: 10.1007/s11224-011-9909-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
|
26
|
Kharat B, Deshmukh V, Chaudhari A. Cyclic and ladder hydrogen bonded cyanamide oligomers: a density functional theory and many-body analysis approach. Struct Chem 2011. [DOI: 10.1007/s11224-011-9841-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
27
|
Abstract
The origin of the intermolecular interaction, especially the many-body interaction, in eight low-lying water hexamer structures (prism, cage, book-1, book-2, cyclic-chair, bag, cyclic-boat-1, and cyclic-boat-2) is unraveled using the localized molecular orbital energy decomposition analysis (LMO-EDA) method at the second-order Møller-Plesset perturbation (MP2) level of theory with a large basis set. It is found that the relative stabilities of these hexamer structures are determined by delicate balances between different types of interaction. According to LMO-EDA, electrostatic and exchange interactions are strictly pairwise additive. Dispersion interaction in these water hexamer structures is almost pairwise additive, with many-body effects varying from -0.13 to +0.05 kcal/mol. Repulsion interaction is roughly pairwise additive, with many-body effects varying from -0.84 to -0.62 kcal/mol. Polarization interaction is not pairwise additive, with many-body effects varying from -13.10 to -8.85 kcal/mol.
Collapse
Affiliation(s)
- Yiming Chen
- Department of Chemistry, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | | |
Collapse
|
28
|
Huang Z, Dai Y, Wang H, Yu L. Microsolvation of aminoethanol: a study using DFT combined with QTAIM. J Mol Model 2011; 17:2781-96. [DOI: 10.1007/s00894-011-0973-z] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2010] [Accepted: 01/12/2011] [Indexed: 11/30/2022]
|
29
|
Chaudhari A, Meraj GA, Lee SL. Many-body energies during proton transfer in an aqueous system. J Mol Model 2010; 16:1559-66. [DOI: 10.1007/s00894-010-0675-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2009] [Accepted: 02/02/2010] [Indexed: 11/29/2022]
Affiliation(s)
- Ajay Chaudhari
- School of Physical Sciences, Swami Ramanand Teerth Marathwada University, Nanded, 431 606, India.
| | | | | |
Collapse
|
30
|
Nolan SJ, Bygrave PJ, Allan NL, Manby FR. Comparison of the incremental and hierarchical methods for crystalline neon. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2010; 22:074201. [PMID: 21386379 DOI: 10.1088/0953-8984/22/7/074201] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
We present a critical comparison of the incremental and hierarchical methods for the evaluation of the static cohesive energy of crystalline neon. Both of these schemes make it possible to apply the methods of molecular electronic structure theory to crystalline solids, offering a systematically improvable alternative to density functional theory. Results from both methods are compared with previous theoretical and experimental studies of solid neon and potential sources of error are discussed. We explore the similarities of the two methods and demonstrate how they may be used in tandem to study crystalline solids.
Collapse
Affiliation(s)
- S J Nolan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, BS8 1TS, UK
| | | | | | | |
Collapse
|
31
|
Sebetci A, Beran GJO. Spatially Homogeneous QM/MM for Systems of Interacting Molecules with on-the-Fly ab Initio Force-Field Parametrization. J Chem Theory Comput 2009; 6:155-67. [DOI: 10.1021/ct900545v] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ali Sebetci
- Department of Chemistry, University of California, Riverside, California 92521
| | - Gregory J. O. Beran
- Department of Chemistry, University of California, Riverside, California 92521
| |
Collapse
|
32
|
Wu J, Yan H, Cao W. Study on the nature of unusual hydrogen-bonded complexes FClHX and halogen-bonded complexes FClClF. ACTA ACUST UNITED AC 2009. [DOI: 10.1016/j.theochem.2008.09.039] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
33
|
Szatyłowicz H, Krygowski TM, Panek JJ, Jezierska A. H-Bonded Complexes of Aniline with HF/F− and Anilide with HF in Terms of Symmetry-Adapted Perturbation, Atoms in Molecules, and Natural Bond Orbitals Theories. J Phys Chem A 2008; 112:9895-905. [DOI: 10.1021/jp803592v] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Halina Szatyłowicz
- Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3, 00-664 Warsaw, Poland
| | | | - Jarosław J. Panek
- National Institute of Chemistry, Hajdrihova 19, 1001 Ljubljana, Slovenia, and Faculty of Chemistry, University of Wrocław, F. Joliot-Curie 14, 50-383 Wrocław, Poland
| | - Aneta Jezierska
- National Institute of Chemistry, Hajdrihova 19, 1001 Ljubljana, Slovenia, and Faculty of Chemistry, University of Wrocław, F. Joliot-Curie 14, 50-383 Wrocław, Poland
| |
Collapse
|
34
|
Kumar A, Kołaski M, Lee HM, Kim KS. Photoexcitation and Photoionization Dynamics of Water Photolysis. J Phys Chem A 2008; 112:5502-8. [DOI: 10.1021/jp711485b] [Citation(s) in RCA: 47] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Anupriya Kumar
- Center for Superfunctional Materials and Basic Science Research Institute, Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Maciej Kołaski
- Center for Superfunctional Materials and Basic Science Research Institute, Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Han Myoung Lee
- Center for Superfunctional Materials and Basic Science Research Institute, Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| | - Kwang S. Kim
- Center for Superfunctional Materials and Basic Science Research Institute, Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Namgu, Pohang 790-784, Korea
| |
Collapse
|
35
|
Du S, Francisco JS, Schenter GK, Garrett BC. Many-body decomposition of the binding energies for OH⋅(H2O)2 and OH⋅(H2O)3 complexes. J Chem Phys 2008; 128:084307. [DOI: 10.1063/1.2828522] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Sediki A, Lebsir F, Martiny L, Dauchez M, Krallafa A. Ab initio investigation of the topology and properties of three-dimensional clusters of water (H2O)n. Food Chem 2008. [DOI: 10.1016/j.foodchem.2007.03.080] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
37
|
van der Avoird A, Szalewicz K. Water trimer torsional spectrum from accurate ab initio and semiempirical potentials. J Chem Phys 2008; 128:014302. [DOI: 10.1063/1.2812556] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
38
|
Kim D, Lee EC, Kim KS, Tarakeshwar P. Cation−π−Anion Interaction: A Theoretical Investigation of the Role of Induction Energies. J Phys Chem A 2007; 111:7980-6. [PMID: 17655210 DOI: 10.1021/jp073337x] [Citation(s) in RCA: 97] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Cation-pi and the corresponding anion-pi interactions have in general been investigated as binary complexes despite their association with counterions. However, a recent study of the ammonia channel highlights the important but overlooked role of anions in cation-pi interactions. In an effort to examine the structural and energetic consequences of the presence of counterions, we have carried out detailed ab initio calculations on some model cation-pi-anion ternary complexes and evaluated the nonpair potential terms, three-body contributions, and attractive and repulsive energy components of the interaction energy. The presence of the anion in the vicinity of the pi system leads to a large redistribution of electron density and hence leads to an inductive stabilization. The resulting electronic and geometrical changes have important consequences in both chemical and biological systems. Compared to cation-pi-anion ternary complexes, the magnitude of the cation-pi interaction in pi-cation-anion ternary complexes is markedly lower because of charge transfer from the anion to the cation.
Collapse
Affiliation(s)
- Dongwook Kim
- Center for Superfunctional Materials, Department of Chemistry, Pohang University of Science and Technology, San 31, Hyojadong, Pohang 790-784, Korea
| | | | | | | |
Collapse
|
39
|
Dahlke EE, Truhlar DG. Assessment of the pairwise additive approximation and evaluation of many-body terms for water clusters. J Phys Chem B 2007; 110:10595-601. [PMID: 16771303 DOI: 10.1021/jp061039e] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We have tested the ability of four commonly used density functionals (three of which are semilocal and one of which is nonlocal) to outperform accurate pairwise additive approximations in the prediction of binding energies for a series of water clusters ranging in size from dimer to pentamer. Comparison to results obtained with the Weizmann-1 (W1) level of wave function theory shows that while all density functionals are capable of outperforming the accurate pairwise data, the choice of basis set used is crucial to the performance of the method, and if a poor choice of basis set is made the errors obtained with density functional theory (DFT) can be larger than those obtained with the simple pairwise approximation. We have also compared the binding energies and many-body terms determined with DFT to those obtained with W1, and have found that all semilocal functionals have significant errors in the many-body components of the full interactions energy. Despite this limitation, however, we have found that, of the four functionals tested, PBE1W/MG3S is the most accurate for predicting the binding energies of the clusters.
Collapse
|
40
|
Ruckenstein E, Shulgin IL, Shulgin LI. Cooperativity in Ordinary Ice and Breaking of Hydrogen Bonds. J Phys Chem B 2007; 111:7114-21. [PMID: 17550284 DOI: 10.1021/jp071582a] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
The total interaction energy between two H-bonded water molecules in a condensed phase is composed of a binding energy between them and an energy due to a cooperative effect. An approximate simple expression is suggested for the dependence of the interaction energy between two H-bonded water molecules on the number of neighboring water molecules with which they are H-bonded. Using this expression, the probabilities of breaking a H bond with various numbers of H-bonded neighbors are estimated. These probabilities are used in computer simulations of the breaking of specified fractions of H bonds in an ordinary (hexagonal) ice. A large "piece" of hexagonal ice (up to 8 millions molecules) is built up, and various percentages of H bonds are considered broken. It is shown that 62-63% of H bonds must be broken in order to disintegrate the "piece" of ice into disconnected clusters. This value is only a little larger than the percolation threshold (61%) predicted both by the percolation theory for tetrahedral ice and by simulations in which all H bonds were considered equally probable to be broken. When the percentage of broken bonds is smaller than 62-63%, there is a network of H-bonded molecules which contains the overwhelming majority of water molecules. This result contradicts some models of water which consider that water consists of a mixture of water clusters of various sizes. The distribution of water molecules with unequal probabilities for breaking is compared with the simulation involving equal probabilities for breaking. It was found that in the former case, there is an enhanced number of water monomers without H bonds, that the numbers of 2- and 3-bonded molecules are smaller, and the number of 4-bonded molecules is larger than in the latter case.
Collapse
Affiliation(s)
- Eli Ruckenstein
- Department of Chemical and Biological Engineering, State University of New York at Buffalo, Amherst, New York 14260, USA. feaeliru@acsu. buffalo.edu
| | | | | |
Collapse
|
41
|
Theoretical investigation on intermolecular interactions between HCN and HNC: The nature and thermodynamic properties. ACTA ACUST UNITED AC 2007. [DOI: 10.1016/j.theochem.2006.11.030] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
42
|
|
43
|
Wu J, Zhang J, Wang Z, Cao W. Theoretical Investigations of the Nature of Interaction of ClF with HF, H2O, and NH3. J Chem Theory Comput 2006; 3:95-102. [DOI: 10.1021/ct600229n] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Junyong Wu
- Institute of Modern Catalysis, Beijing University of Chemical Technology, State Key Laboratory of Chemical Resource Engineering, Beijing 100029, China
| | - Jingchang Zhang
- Institute of Modern Catalysis, Beijing University of Chemical Technology, State Key Laboratory of Chemical Resource Engineering, Beijing 100029, China
| | - Zhaoxu Wang
- Institute of Modern Catalysis, Beijing University of Chemical Technology, State Key Laboratory of Chemical Resource Engineering, Beijing 100029, China
| | - Weiliang Cao
- Institute of Modern Catalysis, Beijing University of Chemical Technology, State Key Laboratory of Chemical Resource Engineering, Beijing 100029, China
| |
Collapse
|
44
|
Mrázek J, Burda JV. Can the pH value of water solutions be estimated by quantum chemical calculations of small water clusters? J Chem Phys 2006; 125:194518. [PMID: 17129134 DOI: 10.1063/1.2363383] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In the study, various water clusters were explored from the point of view of the proton transfer between H-bonded neighbors. A relatively modest approach--the MP26-31++G(d,p) level--was chosen as acceptable considering the fact that also larger systems have to be included. The tight ion-pair model (with usually three fixed O-O distances) was adopted for the autodissociation process. First, cluster-estimated pH values rapidly decrease as cluster size increases from 2 to 6. For larger clusters in gas phase, the topology of H bonds plays an important role, varying pH from 7 to 13 in hexamers and from 5 to 15 in octamer clusters. The relationship energy/distance was quantified, too. Enhancing our model with the conductorlike screening model (COSMO) approach brought significant improvement in description of the autodissociation reaction with a stable zwitterionic structure. However, when the O-O restrictions were released, the small barrier for backward reaction disappeared, reforming neutral cluster spontaneously. Also Boltzmann weighting procedure was applied, and for the explored clusters in vacuo, the series of pH 25-18-14-13-10 was obtained for cluster sizes n=2, 3, 4, 6, 8. Using the COSMO approach, the analogous series is 15-14-12-10-9. The limit of the series is still about two to three units above the experimentally known pH. In order to reach the size-independent (bulk) value, larger clusters are needed. However, the situation is far from hopeless since (as it was proven in the study) four-coordinated molecules are not involved in the proton transfer process directly; they can only be a part of the surrounding environment.
Collapse
Affiliation(s)
- Jirí Mrázek
- Department of Chemical Physics and Optics, Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 121 16 Prague 2, Czech Republic
| | | |
Collapse
|
45
|
Wang ZX, Zhang JC, Cao WL. Study on the Cation-π Interactions between Ammonium Ion and Aromatic π Systems. CHINESE J CHEM 2006. [DOI: 10.1002/cjoc.200690287] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
46
|
Vilela AF, Barreto PR, Gargano R, Cunha CR. Ab initio studies of hydrogen-bonded complexes: The H2O dimer, trimer and H2OCO. Chem Phys Lett 2006. [DOI: 10.1016/j.cplett.2006.06.040] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
47
|
Torheyden M, Jansen G. A new potential energy surface for the water dimer obtained from separate fits ofab initioelectrostatic, induction, dispersion and exchange energy contributions. Mol Phys 2006. [DOI: 10.1080/00268970600679188] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
48
|
Vaupel S, Brutschy B, Tarakeshwar P, Kim KS. Characterization of Weak NH−π Intermolecular Interactions of Ammonia with Various Substituted π-Systems. J Am Chem Soc 2006; 128:5416-26. [PMID: 16620113 DOI: 10.1021/ja056454j] [Citation(s) in RCA: 100] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Among the several weak intermolecular interactions pervading chemistry and biology, the NH-pi interaction is one of the most widely known. Nevertheless its weak nature makes it one of the most poorly understood and characterized interactions. The present study details the results obtained on gas-phase complexes of ammonia with various substituted pi systems using both laser vibrational spectroscopy and ab initio calculations. The spectroscopic measurements carried out by applying one-color resonant two-photon ionization (R2PI) and IR-vibrational predissociation spectroscopy in the region of the NH stretches yield the first experimental NH stretching shifts of ammonia upon its interaction with various kinds of pi-systems. The experiments were complemented by ab initio calculations and energy decompositions, carried out at the second-order Møller-Plesset (MP2) level of theory. The observed complexes show characteristic vibrational spectra which are very similar to the calculated ones, thereby allowing an in-depth analysis of the interaction forces and energies. The interaction energy of the conformers responsible for the observed vibrational spectra has the maximum contribution from dispersion energies. This implies that polarizabilities of the pi-electron systems play a very important role in governing the nature and geometry of the NH-pi interaction. The larger polarizability of ammonia as compared to water and the tendency to maximize the dispersion energy implies that the characteristics of the NH-pi interactions are markedly different from that of the corresponding OH-pi interactions.
Collapse
Affiliation(s)
- Sascha Vaupel
- Institut für Physikalische und Theoretische Chemie, J. W. Goethe-Universität Frankfurt, Marie-Curie-Str. 11, D-60439 Frankfurt/Main, Germany
| | | | | | | |
Collapse
|
49
|
Piquemal JP, Marquez A, Parisel O, Giessner-Prettre C. A CSOV study of the difference between HF and DFT intermolecular interaction energy values: The importance of the charge transfer contribution. J Comput Chem 2005; 26:1052-62. [PMID: 15898112 DOI: 10.1002/jcc.20242] [Citation(s) in RCA: 95] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Intermolecular interaction energy decompositions using the Constrained Space Orbital Variation (CSOV) method are carried out at the Hartree-Fock level on the one hand and using DFT with usual GGA functionals on the other for a number of model complexes to analyze the role of electron correlation in the intermolecular stabilization energy. In addition to the overall stabilization, the results provide information on the variation, with respect to the computational level, of the different contributions to the interaction energy. The complexes studied are the water linear dimer, the N-methylformamide dimer, the nucleic acid base pairs, the benzene-methane and benzene-N2 van der Waals complexes, [Cu+ -(ImH)3]2, where "ImH" stands for the Imidazole ligand, and ImH-Zn++. The variation of the frozen core energy (the sum of the intermolecular electrostatic energy and the Pauli repulsion energy) calculated from the unperturbed orbitals of the interacting entities indicates that the intramolecular correlation contributions can be stabilizing as well as destabilizing, and that general trends can be derived from the results obtained using usual density functionals. The most important difference between the values obtained from HF and DFT computations concerns the charge transfer contribution, which, in most cases, undergoes the largest increase. The physical meaning of these results is discussed. The present work gives reference calculations that might be used to parametrize new correlated molecular mechanics potentials.
Collapse
Affiliation(s)
- Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, UMR 7616-CNRS/UPMC, Université P. & M. Curie, Case courrier 137, 4 place Jussieu, F. 75252 Paris Cedex 05, France
| | | | | | | |
Collapse
|
50
|
Ejtehadi MR, Avall SP, Plotkin SS. Three-body interactions improve the prediction of rate and mechanism in protein folding models. Proc Natl Acad Sci U S A 2004; 101:15088-93. [PMID: 15469920 PMCID: PMC524050 DOI: 10.1073/pnas.0403486101] [Citation(s) in RCA: 81] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2004] [Indexed: 11/18/2022] Open
Abstract
Here we study the effects of many-body interactions on rate and mechanism in protein folding by using the results of molecular dynamics simulations on numerous coarse-grained Calpha-model single-domain proteins. After adding three-body interactions explicitly as a perturbation to a Gō-like Hamiltonian with native pairwise interactions only, we have found (i) a significantly increased correlation with experimental phi values and folding rates, (ii) a stronger correlation of folding rate with contact order, matching the experimental range in rates when the fraction of three-body energy in the native state is approximately 20%, and (iii) a considerably larger amount of three-body energy present in chymotripsin inhibitor than in the other proteins studied.
Collapse
Affiliation(s)
- M R Ejtehadi
- Department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada V6T 1Z1.
| | | | | |
Collapse
|