101
|
Jing Z, Qi R, Thibonnier M, Ren P. Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides. J Chem Theory Comput 2019; 15:6422-6432. [PMID: 31553600 PMCID: PMC6889957 DOI: 10.1021/acs.jctc.9b00519] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712
| | - Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712
| | | | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712
| |
Collapse
|
102
|
Bagwill CS, Wireduaah S, Cusworth B, Korba J, Kirkpatrick CC, Lewis M, Stalcup AM. Use of HPLC retention to investigate new P descriptors designed to represent ion-π interactions. J LIQ CHROMATOGR R T 2019. [DOI: 10.1080/10826076.2019.1667822] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
| | | | - Brian Cusworth
- School of Medicine, Washington University, St. Louis, MO, USA
| | | | | | - Michael Lewis
- Department of Chemistry, Saint Louis University, St. Louis, MO, USA
| | - Apryll M. Stalcup
- School of Chemical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland
| |
Collapse
|
103
|
Robalo JR, Streacker LM, Mendes de Oliveira D, Imhof P, Ben-Amotz D, Verde AV. Hydrophobic but Water-Friendly: Favorable Water–Perfluoromethyl Interactions Promote Hydration Shell Defects. J Am Chem Soc 2019; 141:15856-15868. [DOI: 10.1021/jacs.9b06862] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- João R. Robalo
- Department of Theory & Bio-systems, Max Planck Institute for Colloids and Interfaces, Science Park, Potsdam 14476, Germany
| | - Louis M. Streacker
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | | | - Petra Imhof
- Institute for Theoretical Physics, Free University of Berlin, Arnimallee 14, 14195 Berlin, Germany
| | - Dor Ben-Amotz
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ana Vila Verde
- Department of Theory & Bio-systems, Max Planck Institute for Colloids and Interfaces, Science Park, Potsdam 14476, Germany
| |
Collapse
|
104
|
Derricotte WD. Symmetry-Adapted Perturbation Theory Decomposition of the Reaction Force: Insights into Substituent Effects Involved in Hemiacetal Formation Mechanisms. J Phys Chem A 2019; 123:7881-7891. [PMID: 31429558 DOI: 10.1021/acs.jpca.9b06865] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The decomposition of the reaction force based on symmetry-adapted perturbation theory (SAPT) has been proposed. This approach was used to investigate the substituent effects along the reaction coordinate pathway for the hemiacetal formation mechanism between methanol and substituted aldehydes of the form CX3CHO (X = H, F, Cl, and Br), providing a quantitative evaluation of the reaction-driving and reaction-retarding force components. Our results highlight the importance of more favorable electrostatic and induction effects in the reactions involving halogenated aldehydes that leads to lower activation energy barriers. These substituent effects are further elucidated by applying the functional-group partition of symmetry-adapted perturbation theory (F-SAPT). The results show that the reaction is largely driven by favorable direct noncovalent interactions between the CX3 group on the aldehyde and the OH group on methanol.
Collapse
Affiliation(s)
- Wallace D Derricotte
- Department of Chemistry , Morehouse College , Atlanta , Georgia 30314 , United States
| |
Collapse
|
105
|
Qi R, Walker B, Jing Z, Yu M, Stancu G, Edupuganti R, Dalby KN, Ren P. Computational and Experimental Studies of Inhibitor Design for Aldolase A. J Phys Chem B 2019; 123:6034-6041. [PMID: 31268712 DOI: 10.1021/acs.jpcb.9b04551] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Glycolytic enzyme fructose-bisphosphate aldolase A is an emerging therapeutic target in cancer. Recently, we have solved the crystal structure of murine aldolase in complex with naphthalene-2,6-diyl bisphosphate (ND1) that served as a template of the design of bisphosphate-based inhibitors. In this work, a series of ND1 analogues containing difluoromethylene (-CF2), methylene (-CH2), or aldehyde substitutions were designed. All designed compounds were studied using molecular dynamics (MD) simulations with the AMOEBA force field. Both energetics and structural analyses have been done to understand the calculated binding free energies. The average distances between ligand and protein atoms for ND1 were very similar to those for the ND1 crystal structure, which indicates that our MD simulation is sampling the correct conformation well. CF2 insertion lowers the binding free energy by 10-15 kcal/mol, while CF2 substitution slightly increases the binding free energy, which matches the experimental measurement. In addition, we found that NDB with two CF2 insertions, the strongest binder, is entropically driven, while others including NDA with one CF2 insertion are all enthalpically driven. This work provides insights into the mechanisms underlying protein-phosphate binding and enhances the capability of applying computational and theoretical frameworks to model, predict, and design diagnostic strategies targeting cancer.
Collapse
Affiliation(s)
| | | | | | - Maiya Yu
- Department of Biochemistry and Mathematics , University of Michigan , Ann Arbor , Michigan 48109 , United States
| | | | | | | | | |
Collapse
|
106
|
Bootsma AN, Doney AC, Wheeler SE. Predicting the Strength of Stacking Interactions between Heterocycles and Aromatic Amino Acid Side Chains. J Am Chem Soc 2019; 141:11027-11035. [DOI: 10.1021/jacs.9b00936] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Andrea N. Bootsma
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Analise C. Doney
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Steven E. Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
107
|
Ang SJ, Ser CT, Wong MW. Modeling halogen bonding with planewave density functional theory: Accuracy and challenges. J Comput Chem 2019; 40:1829-1835. [DOI: 10.1002/jcc.25835] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2019] [Revised: 02/20/2019] [Accepted: 03/18/2019] [Indexed: 12/12/2022]
Affiliation(s)
- Shi Jun Ang
- NUS Graduate School for Integrative Sciences and Engineering University Hall, Tan Chin Tuan Wing, #04‐02, 21 Lower Kent Ridge, Singapore 119077 Singapore
| | - Cher Tian Ser
- Department of ChemistryNational University of Singapore 3 Science Drive 3, Singapore 117543 Singapore
| | - Ming Wah Wong
- NUS Graduate School for Integrative Sciences and Engineering University Hall, Tan Chin Tuan Wing, #04‐02, 21 Lower Kent Ridge, Singapore 119077 Singapore
- Department of ChemistryNational University of Singapore 3 Science Drive 3, Singapore 117543 Singapore
| |
Collapse
|
108
|
Maji R, Ugale H, Wheeler SE. Understanding the Reactivity and Selectivity of Fluxional Chiral DMAP-Catalyzed Kinetic Resolutions of Axially Chiral Biaryls. Chemistry 2019; 25:4452-4459. [PMID: 30657217 DOI: 10.1002/chem.201806068] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 01/17/2019] [Indexed: 12/17/2022]
Abstract
Fluxional chiral DMAP-catalyzed kinetic resolutions of axially chiral biaryls were examined using density functional theory. Computational analyses lead to a revised understanding of this reaction in which the interplay of numerous non-covalent interactions control the conformation and flexibility of the active catalyst, the preferred mechanism, and the stereoselectivity. Notably, while the DMAP catalyst itself is confirmed to be highly fluxional, electrostatically driven π⋅⋅⋅π+ interactions render the active, acylated form of the catalyst highly rigid, explaining its pronounced stereoselectivity.
Collapse
Affiliation(s)
- Rajat Maji
- Department of Chemistry, Texas A&M University, College Station, TX, 77842, USA
| | - Heena Ugale
- Department of Chemistry, Texas A&M University, College Station, TX, 77842, USA
| | - Steven E Wheeler
- Department of Chemistry, Texas A&M University, College Station, TX, 77842, USA.,Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, GA, 30602, USA
| |
Collapse
|
109
|
Study of the Wilcox torsion balance in solution for a Tröger's base derivative with hexyl-and heptyl substituents using a combined molecular mechanics and quantum chemistry approach. J Mol Model 2019; 25:69. [PMID: 30778689 DOI: 10.1007/s00894-019-3935-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] [Received: 10/31/2018] [Accepted: 01/21/2019] [Indexed: 02/06/2023]
Abstract
The folding equilibrium of the Wilcox torsion balance in solution has been studied using a molecular mechanics method for sampling the conformational space and semi-empirical and density-functional quantum chemistry methods for characterizing the relative stabilities of various solute-solvent clusters extracted with the aid of the MD-quench technique from the different simulations that were performed. The role of the solvent environment has been analyzed by choosing four solvents of different polarities, namely water, acetone, tetrachloromethane, and n-hexane. In all cases, it is found that the attractive intramolecular interactions in folded conformations are strongly compensated by the increase of the solute-solvent interaction energies when the molecule unfolds. The latter can be well explained by the larger number of solvent molecules that can bind to the Wilcox molecule when in an unfolded conformation. The results of this work therefore support the experimental results of Yang et al. (Nature Chem 5:1006, 2013) that the folding free energy of the Wilcox balance is strongly reduced in solution as compared to the gas phase.
Collapse
|
110
|
Hapka M, Przybytek M, Pernal K. Second-Order Dispersion Energy Based on Multireference Description of Monomers. J Chem Theory Comput 2019; 15:1016-1027. [PMID: 30525591 DOI: 10.1021/acs.jctc.8b01058] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
We propose a method for calculating a second-order dispersion energy for weakly interacting multireference systems in arbitrary electronic states. It is based on response properties obtained from extended random phase approximation equations. The introduced formalism is general and requires only one- and two-particle reduced density matrices of monomers. We combine the new method with either generalized valence bond perfect pairing (GVB) or complete active space (CAS) self-consistent field description of the interacting systems. In addition to a general scheme, three approximations, leading to significant reduction of the computational cost, are developed by exploiting Dyall partitioning of the monomer Hamiltonians. For model multireference systems (H2···H2 and Be···Be) the method is accurate, unlike its single-reference-based counterpart. Neither GVB nor CAS description of single-reference monomers improves the dispersion energy with respect to the Hartree-Fock-based results.
Collapse
Affiliation(s)
- Michał Hapka
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland.,Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Michał Przybytek
- Faculty of Chemistry , University of Warsaw , ul. L. Pasteura 1 , 02-093 Warsaw , Poland
| | - Katarzyna Pernal
- Institute of Physics , Lodz University of Technology , ul. Wolczanska 219 , 90-924 Lodz , Poland
| |
Collapse
|
111
|
Varadwaj A, Marques HM, Varadwaj PR. Is the Fluorine in Molecules Dispersive? Is Molecular Electrostatic Potential a Valid Property to Explore Fluorine-Centered Non-Covalent Interactions? Molecules 2019; 24:E379. [PMID: 30678158 PMCID: PMC6384640 DOI: 10.3390/molecules24030379] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2018] [Revised: 01/15/2019] [Accepted: 01/18/2019] [Indexed: 11/23/2022] Open
Abstract
Can two sites of positive electrostatic potential localized on the outer surfaces of two halogen atoms (and especially fluorine) in different molecular domains attract each other to form a non-covalent engagement? The answer, perhaps counterintuitive, is yes as shown here using the electronic structures and binding energies of the interactions for a series of 22 binary complexes formed between identical or different atomic domains in similar or related halogen-substituted molecules containing fluorine. These were obtained using various computational approaches, including density functional and ab initio first-principles theories with M06-2X, RHF, MP2 and CCSD(T). The physical chemistry of non-covalent bonding interactions in these complexes was explored using both Quantum Theory of Atoms in Molecules and Symmetry Adapted Perturbation Theories. The surface reactivity of the 17 monomers was examined using the Molecular Electrostatic Surface Potential approach. We have demonstrated inter alia that the dispersion term, the significance of which is not always appreciated, which emerges either from an energy decomposition analysis, or from a correlated calculation, plays a structure-determining role, although other contributions arising from electrostatic, exchange-repulsion and polarization effects are also important. The 0.0010 a.u. isodensity envelope, often used for mapping the electrostatic potential is found to provide incorrect information about the complete nature of the surface reactive sites on some of the isolated monomers, and can lead to a misinterpretation of the results obtained.
Collapse
Affiliation(s)
- Arpita Varadwaj
- Department of Chemical System Engineering, School of Engineering, The University of Tokyo 7-3-1, Hongo, Bunkyo-ku 113-8656, Japan.
- National Institute of Advanced Industrial Science and Technology, 1 Chome-1-1 Umezono, Tsukuba, Ibaraki Prefecture, Ibaraki 305-8560, Japan.
| | - Helder M Marques
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, Johannesburg 2050, South Africa.
| | - Pradeep R Varadwaj
- Department of Chemical System Engineering, School of Engineering, The University of Tokyo 7-3-1, Hongo, Bunkyo-ku 113-8656, Japan.
- National Institute of Advanced Industrial Science and Technology, 1 Chome-1-1 Umezono, Tsukuba, Ibaraki Prefecture, Ibaraki 305-8560, Japan.
| |
Collapse
|
112
|
Antonijević IS, Malenov DP, Hall MB, Zarić SD. Study of stacking interactions between two neutral tetrathiafulvalene molecules in Cambridge Structural Database crystal structures and by quantum chemical calculations. ACTA CRYSTALLOGRAPHICA SECTION B-STRUCTURAL SCIENCE CRYSTAL ENGINEERING AND MATERIALS 2019; 75:1-7. [DOI: 10.1107/s2052520618015494] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/08/2018] [Accepted: 11/02/2018] [Indexed: 11/10/2022]
Abstract
Tetrathiafulvalene (TTF) and its derivatives are very well known as electron donors with widespread use in the field of organic conductors and superconductors. Stacking interactions between two neutral TTF fragments were studied by analysing data from Cambridge Structural Database crystal structures and by quantum chemical calculations. Analysis of the contacts found in crystal structures shows high occurrence of parallel displaced orientations of TTF molecules. In the majority of the contacts, two TTF molecules are displaced along their longer C
2 axis. The most frequent geometry has the strongest TTF–TTF stacking interaction, with CCSD(T)/CBS energy of −9.96 kcal mol−1. All the other frequent geometries in crystal structures are similar to geometries of the minima on the calculated potential energy surface.
Collapse
|
113
|
Wolf ME, Zhang B, Turney JM, Schaefer HF. A comparison between hydrogen and halogen bonding: the hypohalous acid–water dimers, HOX⋯H2O (X = F, Cl, Br). Phys Chem Chem Phys 2019; 21:6160-6170. [DOI: 10.1039/c9cp00422j] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Hypohalous acids (HOX) are a class of molecules that play a key role in the atmospheric seasonal depletion of ozone and have the ability to form both hydrogen and halogen bonds.
Collapse
Affiliation(s)
- Mark E. Wolf
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- USA
| | - Boyi Zhang
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- USA
| | - Justin M. Turney
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry
- University of Georgia
- Athens
- USA
| |
Collapse
|
114
|
Harding DP, Bootsma AN, Wheeler SE. Better Sensing through Stacking: The Role of Non-Covalent Interactions in Guanine-Binding Sensors. J Phys Chem B 2018; 123:487-495. [DOI: 10.1021/acs.jpcb.8b12158] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Drew P. Harding
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Andrea N. Bootsma
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Steven E. Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
- Center for Computational Quantum Chemistry, Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
115
|
Bootsma AN, Wheeler SE. Tuning Stacking Interactions between Asp-Arg Salt Bridges and Heterocyclic Drug Fragments. J Chem Inf Model 2018; 59:149-158. [PMID: 30507185 DOI: 10.1021/acs.jcim.8b00563] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Stacking interactions can play an integral role in the strength and selectivity of protein-drug binding and are of particular interest given the ubiquity and variety of heterocyclic fragments in drugs. In addition to traditional stacking interactions between aromatic rings, stacking interactions involving heterocyclic drug fragments and protein salt bridges have also been observed. These "salt-bridge stacking interactions" can be quite strong but are not well understood. We studied stacked dimers of the acetate···guanidinium ion pair with a diverse set of 63 heterocycles using robust ab initio methods. The computed interaction energies span more than 10 kcal mol-1, indicating the sensitivity of these salt-bridge stacking interactions to heterocycle features. Trends in both the strength and preferred geometry of these interactions can be understood through analyses of the electrostatic potentials and electric fields above the heterocycles. We have developed new heterocycle descriptors that quantify these effects and used them to create robust predictors of the strength of salt-bridge stacking interactions both in the gas phase and a protein-like dielectric environment. These predictive tools, combined with a set of qualitative guidelines, should facilitate the choice of heterocycles that maximize salt-bridge stacking interactions in drug binding sites.
Collapse
Affiliation(s)
- Andrea N Bootsma
- Department of Chemistry , Texas A&M University , College Station , Texas 77842 , United States.,Center for Computational Quantum Chemistry, Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| | - Steven E Wheeler
- Center for Computational Quantum Chemistry, Department of Chemistry , University of Georgia , Athens , Georgia 30602 , United States
| |
Collapse
|
116
|
Pastorczak E, Pernal K. Molecular interactions in electron-groups embedding generalized valence bond picture. Theor Chem Acc 2018. [DOI: 10.1007/s00214-018-2378-x] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
117
|
Viquez Rojas CI, Fine J, Slipchenko LV. Exchange-repulsion energy in QM/EFP. J Chem Phys 2018; 149:094103. [PMID: 30195305 DOI: 10.1063/1.5043107] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The effective fragment potential (EFP) is a quantum mechanics (QM)-based model designed to accurately describe intermolecular interactions. Hybrid QM/EFP calculations combine quantum mechanical methods with an EFP embedding to study complex systems in which many-body effects are relevant. As in EFP-only calculations, non-bonded interactions between the QM region and EFP fragments are computed as a sum of electrostatic, polarization, dispersion, and exchange-repulsion energies. The exchange-repulsion term is a computational bottleneck of the EFP calculations. Here, we present a general procedure for computing the QM/EFP exchange-repulsion interactions based on one-electron contributions to the QM Hamiltonian, by using Gaussian functions to represent localized molecular orbitals of the effective fragments. The accuracy of the exchange-repulsion and total QM/EFP interaction energies is evaluated on a diverse set of dimers, including complexes from the S22 dataset of non-covalent interactions. In most cases, the QM/EFP energies are at least as accurate as corresponding EFP energies. A simple and computationally efficient form of the introduced QM/EFP exchange-repulsion term will facilitate further developments and applications of QM/EFP methods.
Collapse
Affiliation(s)
| | - Jonathan Fine
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | | |
Collapse
|
118
|
Rackers JA, Liu C, Ren P, Ponder JW. A physically grounded damped dispersion model with particle mesh Ewald summation. J Chem Phys 2018; 149:084115. [PMID: 30193468 DOI: 10.1063/1.5030434] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Accurate modeling of dispersion is critical to the goal of predictive biomolecular simulations. To achieve this accuracy, a model must be able to correctly capture both the short-range and asymptotic behavior of dispersion interactions. We present here a damped dispersion model based on the overlap of charge densities that correctly captures both regimes. The overlap damped dispersion model represents a classical physical interpretation of dispersion: the interaction between the instantaneous induced dipoles of two distinct charge distributions. This model is shown to be an excellent fit with symmetry adapted perturbation theory dispersion energy calculations, yielding an RMS error on the S101x7 database of 0.5 kcal/mol. Moreover, the damping function used in this model is wholly derived and parameterized from the electrostatic dipole-dipole interaction, making it not only physically grounded but transferable as well.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Jay W Ponder
- Program in Computational and Molecular Biophysics, Washington University School of Medicine, Saint Louis, Missouri 63110, USA
| |
Collapse
|
119
|
Awoonor-Williams E, Rowley CN. The hydration structure of methylthiolate from QM/MM molecular dynamics. J Chem Phys 2018; 149:045103. [PMID: 30068187 DOI: 10.1063/1.5038010] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Thiols are widely present in biological systems, most notably as the side chain of cysteine amino acids in proteins. Thiols can be deprotonated to form a thiolate which affords a diverse range of enzymatic activity and modes for chemical modification of proteins. Parameters for modeling thiolates using molecular mechanical force fields have not yet been validated, in part due to the lack of structural data on thiolate solvation. Here, the CHARMM36 and Amber models for thiolates in aqueous solutions are assessed using free energy perturbation and hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations. The hydration structure of methylthiolate was calculated from 1 ns of QM/MM MD (PBE0-D3/def2-TZVP//TIP3P), which shows that the water-S- distances are approximately 2 Å with a coordination number near 6. The CHARMM thiolate parameters predict a thiolate S radius close to the QM/MM value and predict a hydration Gibbs energy of -329.2 kJ/mol, close to the experimental value of -318 kJ/mol. The cysteine thiolate model in the Amber force field underestimates the thiolate radius by 0.2 Å and overestimates the thiolate hydration energy by 119 kJ/mol because it uses the same Lennard-Jones parameters for thiolates as for thiols. A recent Drude polarizable model for methylthiolate with optimized thiolate parameters also performs well. SAPT2+ [Symmetry Adapted Perturbation Theory (SAPT)] analysis indicates that exchange repulsion is larger for the methylthiolate, consistent with it having a more diffuse electron density distribution in comparison with the parent thiol. These data demonstrate that it is important to define distinct non-bonded parameters for the protonated/deprotonated states of amino acid side chains in molecular mechanical force fields.
Collapse
Affiliation(s)
- Ernest Awoonor-Williams
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador A1C 5S7, Canada
| | - Christopher N Rowley
- Department of Chemistry, Memorial University of Newfoundland, St. John's, Newfoundland and Labrador A1C 5S7, Canada
| |
Collapse
|
120
|
Heindel JP, Knodel ES, Schofield DP. Origin of Many-Body Vibrational Frequency Shifts in Water Clusters. J Phys Chem A 2018; 122:6724-6735. [DOI: 10.1021/acs.jpca.8b04380] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Joseph P. Heindel
- Department of Chemistry and Biochemistry, Seattle Pacific University, Seattle, Washington 98119, United States
| | - Elizabeth S. Knodel
- Department of Chemistry and Biochemistry, Seattle Pacific University, Seattle, Washington 98119, United States
| | - Daniel P. Schofield
- Department of Chemistry and Biochemistry, Seattle Pacific University, Seattle, Washington 98119, United States
| |
Collapse
|
121
|
Stasyuk OA, Sedlak R, Guerra CF, Hobza P. Comparison of the DFT-SAPT and Canonical EDA Schemes for the Energy Decomposition of Various Types of Noncovalent Interactions. J Chem Theory Comput 2018; 14:3440-3450. [PMID: 29926727 DOI: 10.1021/acs.jctc.8b00034] [Citation(s) in RCA: 58] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Interaction energies computed with density functional theory can be divided into physically meaningful components by symmetry-adapted perturbation theory (DFT-SAPT) or the canonical energy decomposition analysis (EDA). In this work, the decomposition results obtained by these schemes were compared for more than 200 hydrogen-, halogen-, and pnicogen-bonded, dispersion-bound, and mixed complexes to investigate their similarity in the evaluation of the nature of noncovalent interactions. BLYP functional with D3(BJ) correction was used for the EDA scheme, whereas asymptotically corrected PBE0 functional for DFT-SAPT provided some of the best combinations for description of noncovalent interactions. Both schemes provide similar results concerning total interaction energies and insight into the individual energy components. For most complexes, the dominant energetic term was identified equally by both decomposition schemes. Because the canonical EDA is computationally less demanding than the DFT-SAPT, the former can be especially used in cases where the systems investigated are very large.
Collapse
Affiliation(s)
- Olga A Stasyuk
- Institute of Organic Chemistry and Biochemistry , Academy of Sciences of the Czech Republic , 166 10 Prague 6, Czech Republic
| | - Robert Sedlak
- Institute of Organic Chemistry and Biochemistry , Academy of Sciences of the Czech Republic , 166 10 Prague 6, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry , Palacký University , 771 46 Olomouc , Czech Republic
| | - Célia Fonseca Guerra
- Department of Theoretical Chemistry and Amsterdam Center for Multiscale Modeling , VU Amsterdam , De Boelelaan 1083 , 1081 HV Amsterdam , The Netherlands.,Leiden Institute of Chemistry, Gorlaeus Laboratories , Leiden University , P.O. Box 9502, 2300 RA Leiden , The Netherlands
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry , Academy of Sciences of the Czech Republic , 166 10 Prague 6, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry , Palacký University , 771 46 Olomouc , Czech Republic
| |
Collapse
|
122
|
Tu YJ, Lin Z, Allen MJ, Cisneros GA. Molecular dynamics investigation of water-exchange reactions on lanthanide ions in water/1-ethyl-3-methylimidazolium trifluoromethylsufate ([EMIm][OTf]). J Chem Phys 2018; 148:024503. [PMID: 29331119 DOI: 10.1063/1.4997008] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
We report a kinetic study of the water exchange on lanthanide ions in water/[1-ethyl-3-methylimidazolium][trifluoromethylsufate] (water/[EMIm][OTf]). The results from 17O-NMR measurements show that the water-exchange rates in water/[EMIm][OTf] increase with decreasing size of the lanthanide ions. This trend for water-exchange is similar to the previously reported trend in water/1-ethyl-3-methylimidazolium ethyl sulfate (water/[EMIm][EtSO4]) but opposite to that in water. To gain atomic-level insight into these water-exchange reactions, molecular dynamics simulations for lanthanide ions in water/[EMIm][OTf] have been performed using the atomic-multipole-optimized-energetics-for-biomolecular-application polarizable force field. Our molecular dynamics simulations reproduce the experimental water-exchange rates in terms of the trend and provide possible explanations for the observed experimental behavior. The smaller lanthanide ions in water/[EMIm][OTf] undergo faster water exchange because the smaller lanthanide ions coordinate to the first shell [OTf]- anions more tightly, resulting in a stronger screening effect for the second-shell water. The screening effect weakens the interaction of the lanthanide ions with the second-shell water molecules, facilitating the dissociation of water from the second-shell and subsequent association of water molecules from the outer solvation shells.
Collapse
Affiliation(s)
- Yi-Jung Tu
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - Zhijin Lin
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - Matthew J Allen
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, USA
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, USA
| |
Collapse
|
123
|
Abstract
Even though (H2O)2 and (HF)2 are arguably the most thoroughly characterized prototypes for hydrogen bonding, their heterogeneous analogue H2O···HF has received relatively little attention. Here we report that the experimental dissociation energy ( D0) of this important paradigm for heterogeneous hydrogen bonding is too large by 2 kcal mol-1 or 30% relative to our computed value of 6.3 kcal mol-1. For reference, computational procedures similar to those employed here to compute D0 (large basis set CCSD(T) computations with anharmonic corrections from second-order vibrational perturbation theory) provide results within 0.1 kcal mol-1 of the experimental values for (H2O)2 and (HF)2. Near the CCSD(T) complete basis set limit, the electronic dissociation energy for H2O···HF is ∼4 kcal mol-1 larger than those for (H2O)2 and (HF)2 (∼9 kcal mol-1 for the heterogeneous dimer vs ∼5 kcal mol-1 for the homogeneous dimers). Results reported here from symmetry-adapted perturbation theory computations suggest that this large difference is primarily due to the induction contribution to the interaction energy.
Collapse
Affiliation(s)
- Thomas More Sexton
- Department of Chemistry and Biochemistry , University of Mississippi , University , Mississippi 38677-1848 , United States
| | - J Coleman Howard
- Department of Chemistry , Virginia Tech , Blacksburg , Virginia 24061 , United States
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry , University of Mississippi , University , Mississippi 38677-1848 , United States
| |
Collapse
|
124
|
Lao KU, Herbert JM. Atomic Orbital Implementation of Extended Symmetry-Adapted Perturbation Theory (XSAPT) and Benchmark Calculations for Large Supramolecular Complexes. J Chem Theory Comput 2018; 14:2955-2978. [DOI: 10.1021/acs.jctc.8b00058] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
125
|
Intermolecular dispersion energies from coupled exact-exchange Kohn-Sham excitation energies and vectors. COMPUT THEOR CHEM 2018. [DOI: 10.1016/j.comptc.2018.02.019] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
126
|
Zakharov BA, Michalchuk AAL, Morrison CA, Boldyreva EV. Anisotropic lattice softening near the structural phase transition in the thermosalient crystal 1,2,4,5-tetrabromobenzene. Phys Chem Chem Phys 2018. [PMID: 29537423 DOI: 10.1039/c7cp08609a] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The thermosalient effect (crystal jumping on heating) attracts much attention as both an intriguing academic phenomenon and in relation to its potential for the development of molecular actuators but its mechanism remains unclear. 1,2,4,5-Tetrabromobenzene (TBB) is one of the most extensively studied thermosalient compounds that has been shown previously to undergo a phase transition on heating, accompanied by crystal jumping and cracking. The difference in the crystal structures and intermolecular interaction energies of the low- and high-temperature phases is, however, too small to account for the large stress that arises over the course of the transformation. The energy is released spontaneously, and crystals jump across distances that exceed the crystal size by orders of magnitude. In the present work, the anisotropy of lattice strain is followed across the phase transition by single-crystal X-ray diffraction, focusing on the structural evolution from 273 to 343 K. A pronounced lattice softening is observed close to the transition point, with the structure becoming more rigid immediately after the phase transition. The diffraction studies are further supported by theoretical analysis of pairwise intermolecular energies and zone-centre lattice vibrations. Only three modes are found to monotonically soften up to the phase transition, with complex behaviour exhibited by the remaining lattice modes. The thermosalient effect is delayed with respect to the structural transformation itself. This can originate from the martensitic mechanism of the transformation, and the accumulation of stress associated with vibrational switching across the phase transition. The finding of this study sheds more light on the nature of the thermosalient effect in 1,2,4,5-tetrabromobenzene and can be applicable also to other thermosalient compounds.
Collapse
Affiliation(s)
- Boris A Zakharov
- Institute of Solid State Chemistry and Mechanochemistry, Siberian Branch of the Russian Academy of Sciences, Kutateladze Str. 18, Novosibirsk, 630128, Russian Federation.
| | | | | | | |
Collapse
|
127
|
Heßelmann A. DFT-SAPT Intermolecular Interaction Energies Employing Exact-Exchange Kohn–Sham Response Methods. J Chem Theory Comput 2018; 14:1943-1959. [DOI: 10.1021/acs.jctc.7b01233] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Andreas Heßelmann
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstrasse 3, D-91058 Erlangen, Germany
| |
Collapse
|
128
|
Rackers JA, Wang Q, Liu C, Piquemal JP, Ren P, Ponder JW. An optimized charge penetration model for use with the AMOEBA force field. Phys Chem Chem Phys 2018; 19:276-291. [PMID: 27901142 DOI: 10.1039/c6cp06017j] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The principal challenge of using classical physics to model biomolecular interactions is capturing the nature of short-range interactions that drive biological processes from nucleic acid base stacking to protein-ligand binding. In particular most classical force fields suffer from an error in their electrostatic models that arises from an ability to account for the overlap between charge distributions occurring when molecules get close to each other, known as charge penetration. In this work we present a simple, physically motivated model for including charge penetration in the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) force field. With a function derived from the charge distribution of a hydrogen-like atom and a limited number of parameters, our charge penetration model dramatically improves the description of electrostatics at short range. On a database of 101 biomolecular dimers, the charge penetration model brings the error in the electrostatic interaction energy relative to the ab initio SAPT electrostatic interaction energy from 13.4 kcal mol-1 to 1.3 kcal mol-1. The model is shown not only to be robust and transferable for the AMOEBA model, but also physically meaningful as it universally improves the description of the electrostatic potential around a given molecule.
Collapse
Affiliation(s)
- Joshua A Rackers
- Program in Computational & Molecular Biophysics, Washington University, School of Medicine, Saint Louis, Missouri 63110, USA
| | - Qiantao Wang
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Paris 06, UMR 7616, case courrier 137, 4 place Jussieu, F-75005, Paris, France
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Jay W Ponder
- Department of Chemistry, Washington University in Saint Louis, Saint Louis, Missouri 63130, USA.
| |
Collapse
|
129
|
Parrish RM, Thompson KC, Martínez TJ. Large-Scale Functional Group Symmetry-Adapted Perturbation Theory on Graphical Processing Units. J Chem Theory Comput 2018; 14:1737-1753. [DOI: 10.1021/acs.jctc.7b01053] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Robert M. Parrish
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Keiran C. Thompson
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
130
|
Zhang Y, Wang W, Wang YB. The nature of the noncovalent interactions between fullerene C60 and aromatic hydrocarbons. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.11.004] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
|
131
|
Shirkov L, Sladek V. Benchmark CCSD-SAPT study of rare gas dimers with comparison to MP-SAPT and DFT-SAPT. J Chem Phys 2017; 147:174103. [DOI: 10.1063/1.4997569] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Leonid Shirkov
- Faculty of Chemistry, Adam Mickiewicz University, Umultowska 89b, 61-614 Poznań, Poland
| | - Vladimir Sladek
- Institute of Chemistry–Centre for Glycomics, Slovak Academy of Sciences, 845 38 Bratislava, Slovakia
| |
Collapse
|
132
|
Pinto de Magalhães H, Togni A, Lüthi HP. Importance of Nonclassical σ-Hole Interactions for the Reactivity of λ3-Iodane Complexes. J Org Chem 2017; 82:11799-11805. [DOI: 10.1021/acs.joc.7b01716] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Affiliation(s)
- Halua Pinto de Magalhães
- Departement of Chemistry
and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Antonio Togni
- Departement of Chemistry
and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Hans Peter Lüthi
- Departement of Chemistry
and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
133
|
Begum S, Subramanian R. Bonding and spectroscopic analyses of N 2O-CS 2 and N 2O-OCS heterodimer complexes and their atmospheric consequences. Phys Chem Chem Phys 2017; 19:26412-26422. [PMID: 28944392 DOI: 10.1039/c7cp03936k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Different isomers of the valence isoelectronic pairs of the heterodimers N2O-SCS and N2O-OCS were investigated using MP2 and CCSD(T) methods with the aug-cc-pVXZ (X = D, T) basis set with anharmonic frequency calculations. Interaction energies of both the heterodimers were estimated at the CBS limit and for CCSD(T)/aug-cc-pVTZ. One isomer was located for the N2O-SCS heterodimer. In the case of the N2O-OCS heterodimer, three planar slipped-parallel isomers were characterized. There was no significant change in electron density at the bond critical points (BCPs) of various intra-bonds but the presence of inter-bonds with a weak electron density at the BCPs in the complexes was found as per the results from the quantum theory of atoms in molecules (QTAIM) analysis. The interorbital interactions of the monomer in the heterodimers are discussed herein in terms of the natural bond orbital (NBO) charge transfer. Symmetry-adapted perturbation analysis was performed for all the isomers. In the case of the stable complex of N2O-SCS, induction dominated over the dispersion energy. Among the N2O-OCS isomers, the dispersion effect predominates for the most stable isomer, while for the other two isomers, the induction effect is dominant. Frequency calculations for these complexes showed a number of interactions induced by the low frequency modes in the far IR region. The atmospheric implications are also discussed for these complexes.
Collapse
Affiliation(s)
- Samiyara Begum
- Department of Chemistry, Indian Institute of Technology, Patna, 800013, India.
| | | |
Collapse
|
134
|
Bull CL, Flowitt-Hill G, de Gironcoli S, Küçükbenli E, Parsons S, Pham CH, Playford HY, Tucker MG. ζ-Glycine: insight into the mechanism of a polymorphic phase transition. IUCRJ 2017; 4:569-574. [PMID: 28989714 PMCID: PMC5619850 DOI: 10.1107/s205225251701096x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/08/2017] [Accepted: 07/25/2017] [Indexed: 05/30/2023]
Abstract
Glycine is the simplest and most polymorphic amino acid, with five phases having been structurally characterized at atmospheric or high pressure. A sixth form, the elusive ζ phase, was discovered over a decade ago as a short-lived intermediate which formed as the high-pressure ∊ phase transformed to the γ form on decompression. However, its structure has remained unsolved. We now report the structure of the ζ phase, which was trapped at 100 K enabling neutron powder diffraction data to be obtained. The structure was solved using the results of a crystal structure prediction procedure based on fully ab initio energy calculations combined with a genetic algorithm for searching phase space. We show that the fate of ζ-glycine depends on its thermal history: although at room temperature it transforms back to the γ phase, warming the sample from 100 K to room temperature yielded β-glycine, the least stable of the known ambient-pressure polymorphs.
Collapse
Affiliation(s)
- Craig L. Bull
- ISIS Facility, Rutherford–Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0QX, UK
| | - Giles Flowitt-Hill
- ISIS Facility, Rutherford–Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0QX, UK
- School of Chemistry and Centre for Science at Extreme Conditions, The University of Edinburgh, King’s Buildings, W. Mains Road, Edinburgh EH9 3FJ, UK
| | - Stefano de Gironcoli
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| | - Emine Küçükbenli
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| | - Simon Parsons
- School of Chemistry and Centre for Science at Extreme Conditions, The University of Edinburgh, King’s Buildings, W. Mains Road, Edinburgh EH9 3FJ, UK
| | - Cong Huy Pham
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Helen Y. Playford
- ISIS Facility, Rutherford–Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0QX, UK
| | - Matthew G. Tucker
- ISIS Facility, Rutherford–Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, Oxfordshire OX11 0QX, UK
- Oak Ridge National Laboratory, PO Box 2008, Oak Ridge, TN 37831, USA
| |
Collapse
|
135
|
Deng S, Wang Q, Ren P. Estimating and modeling charge transfer from the SAPT induction energy. J Comput Chem 2017; 38:2222-2231. [PMID: 28766729 DOI: 10.1002/jcc.24864] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2016] [Revised: 06/04/2017] [Accepted: 06/06/2017] [Indexed: 02/05/2023]
Abstract
Recent studies using quantum mechanics energy decomposition methods, for example, SAPT and ALMO, have revealed that the charge transfer energy may play an important role in short ranged inter-molecular interactions, and have a different distance dependence comparing with the polarization energy. However, the charge transfer energy component has been ignored in most current polarizable or non-polarizable force fields. In this work, first, we proposed an empirical decomposition of SAPT induction energy into charge transfer and polarization energy that mimics the regularized SAPT method (ED-SAPT). This empirical decomposition is free of the divergence issue, hence providing a good reference for force field development. Then, we further extended this concept in the context of AMOEBA polarizable force field, proposed a consistent approach to treat the charge transfer phenomenon. Current results show a promising application of this charge transfer model in future force field development. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Shi Deng
- Department of Urology, Institute of Urology, West China Hospital, Sichuan University, Chengdu, 610041, China
| | - Qiantao Wang
- Key Laboratory of Drug Targeting and Drug Delivery System of Education Ministry, West China School of Pharmacy, Sichuan University, Chengdu, 610041, China
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas, 78712
| |
Collapse
|
136
|
Parrish RM, Burns LA, Smith DGA, Simmonett AC, DePrince AE, Hohenstein EG, Bozkaya U, Sokolov AY, Di Remigio R, Richard RM, Gonthier JF, James AM, McAlexander HR, Kumar A, Saitow M, Wang X, Pritchard BP, Verma P, Schaefer HF, Patkowski K, King RA, Valeev EF, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability. J Chem Theory Comput 2017; 13:3185-3197. [PMID: 28489372 PMCID: PMC7495355 DOI: 10.1021/acs.jctc.7b00174] [Citation(s) in RCA: 758] [Impact Index Per Article: 108.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.
Collapse
Affiliation(s)
- Robert M Parrish
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Daniel G A Smith
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Andrew C Simmonett
- National Institutes of Health , National Heart, Lung and Blood Institute, Laboratory of Computational Biology, 5635 Fishers Lane, T-900 Suite, Rockville, Maryland 20852, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University , Tallahassee, Florida 32306-4390, United States
| | - Edward G Hohenstein
- Department of Chemistry and Biochemistry, The City College of New York , New York, New York 10031, United States
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe University , Ankara 06800, Turkey
| | - Alexander Yu Sokolov
- Division of Chemistry and Chemical Engineering, California Institute of Technology , Pasadena, California 91125, United States
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, UiT, The Arctic University of Norway , N-9037 Tromsø, Norway
| | - Ryan M Richard
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Jérôme F Gonthier
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Andrew M James
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Harley R McAlexander
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Masaaki Saitow
- Department of Chemistry and Research Center for Smart Molecules, Rikkyo University , 3-34-1 Nishi-ikebukuro, Toshima-ku, Tokyo 171-8501, Japan
| | - Xiao Wang
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - Benjamin P Pritchard
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| | - Prakash Verma
- Department of Chemistry, Emory University , Atlanta, Georgia 30322, United States
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University , Auburn, Alabama 36849, United States
| | - Rollin A King
- Department of Chemistry, Bethel University , St. Paul, Minnesota 55112, United States
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | | | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia , Athens, Georgia 30602, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech , Blacksburg, Virginia 24061, United States
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology , Atlanta, Georgia 30332-0400, United States
| |
Collapse
|
137
|
Kwapien K, Damergi M, Nader S, El Khoury L, Hobaika Z, Maroun RG, Piquemal JP, Gavara L, Berthomieu D, Hernandez JF, Gresh N. Calibration of 1,2,4-Triazole-3-Thione, an Original Zn-Binding Group of Metallo-β-Lactamase Inhibitors. Validation of a Polarizable MM/MD Potential by Quantum Chemistry. J Phys Chem B 2017; 121:6295-6312. [DOI: 10.1021/acs.jpcb.7b01053] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Karolina Kwapien
- Chemistry
and Biology, Nucléo(s)tides and Immunology for Therapy (CBNIT),
UMR 8601, CNRS, UFR Biomédicale, Paris, France
- Institut Charles-Gerhardt, MACS, UMR 5253 CNRS-ENSCM-UM, 8 rue de l’Ecole Normale, 34296 Montpellier Cedex 5, France
| | - Mirna Damergi
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
- Centre
d’Analyses et de Recherche, UR EGFEM, LSIM, Faculté
des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut 1116-2050, Lebanon
| | - Serge Nader
- Chemistry
and Biology, Nucléo(s)tides and Immunology for Therapy (CBNIT),
UMR 8601, CNRS, UFR Biomédicale, Paris, France
- Centre
d’Analyses et de Recherche, UR EGFEM, LSIM, Faculté
des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut 1116-2050, Lebanon
| | - Léa El Khoury
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
- Centre
d’Analyses et de Recherche, UR EGFEM, LSIM, Faculté
des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut 1116-2050, Lebanon
| | - Zeina Hobaika
- Centre
d’Analyses et de Recherche, UR EGFEM, LSIM, Faculté
des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut 1116-2050, Lebanon
| | - Richard G. Maroun
- Centre
d’Analyses et de Recherche, UR EGFEM, LSIM, Faculté
des Sciences, Saint Joseph University of Beirut, BP 11-514, Riad El Solh, Beirut 1116-2050, Lebanon
| | - Jean-Philip Piquemal
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
- Institut Universitaire de France, Paris Cedex 05, 75231, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Laurent Gavara
- Institut des Biomolécules Max Mousseron,
UMR 5247 CNRS, Faculté de Pharmacie, Université de Montpellier, ENSCM, 15 avenue Charles Flahault, 34093 Montpellier, France
| | - Dorothée Berthomieu
- Institut Charles-Gerhardt, MACS, UMR 5253 CNRS-ENSCM-UM, 8 rue de l’Ecole Normale, 34296 Montpellier Cedex 5, France
| | - Jean-François Hernandez
- Institut des Biomolécules Max Mousseron,
UMR 5247 CNRS, Faculté de Pharmacie, Université de Montpellier, ENSCM, 15 avenue Charles Flahault, 34093 Montpellier, France
| | - Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC, UMR7616 CNRS, Paris, France
| |
Collapse
|
138
|
Reany O, Li A, Yefet M, Gilson MK, Keinan E. Attractive Interactions between Heteroallenes and the Cucurbituril Portal. J Am Chem Soc 2017; 139:8138-8145. [PMID: 28532152 PMCID: PMC5581494 DOI: 10.1021/jacs.6b13005] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In this paper, we report on the noteworthy attractive interaction between organic azides and the portal carbonyls of cucurbiturils. Five homologous bis-α,ω-azidoethylammonium alkanes were prepared, where the number of methylene groups between the ammonium groups ranges from 4 to 8. Their interactions with cucurbit[6]uril were studied by NMR spectroscopy, IR spectroscopy, X-ray crystallography, and computational methods. Remarkably, while the distance between the portal plane and most atoms at the guest end groups increases progressively with the molecular size, the β-nitrogen atoms maintain a constant distance from the portal plane in all homologues, pointing at a strong attractive interaction between the azide group and the portal. Both crystallography and NMR support a specific electrostatic interaction between the carbonyl and the azide β-nitrogen, which stabilizes the canonical resonance form with positive charge on the β-nitrogen and negative charge on the γ-nitrogen. Quantum computational analyses strongly support electrostatics, in the form of orthogonal dipole-dipole interaction, as the main driver for this attraction. The alternative mechanism of n → π* orbital delocalization does not seem to play a significant role in this interaction. The computational studies also indicate that the interaction is not limited to azides, but generalizes to other isoelectronic heteroallene functions, such as isocyanate and isothiocyanate. This essentially unexploited attractive interaction could be more broadly utilized as a tool not only in relation to cucurbituril chemistry, but also for the design of novel supramolecular architectures.
Collapse
Affiliation(s)
- Ofer Reany
- The Avinoam Adam Department of Natural Sciences, The Open
University of Israel, 1 University Road, Ra’anana 43537, Israel
| | - Amanda Li
- Skaggs School of Pharmaceutic and Pharmaceutical Sciences,
University of California San Diego, 9500 Gilman Drive, La Jolla, California,
USA
| | - Maayan Yefet
- The Schulich faculty of Chemistry, Technion-Israel Institute
of Technology, Technion city, 32000 Haifa, Israel
| | - Michael K. Gilson
- Skaggs School of Pharmaceutic and Pharmaceutical Sciences,
University of California San Diego, 9500 Gilman Drive, La Jolla, California,
USA
| | - Ehud Keinan
- The Schulich faculty of Chemistry, Technion-Israel Institute
of Technology, Technion city, 32000 Haifa, Israel
| |
Collapse
|
139
|
MONDAL ANIRBAN, BALASUBRAMANIAN SUNDARAM. Molecular Dynamics Investigation of Efficient SO2 Absorption by Anion-Functionalized Ionic Liquids. J CHEM SCI 2017. [DOI: 10.1007/s12039-017-1236-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
140
|
Unravelling the influence of quantum proton delocalization on electronic charge transfer through the hydrogen bond. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.04.034] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
141
|
Kumar S, Singh SK, Vaishnav JK, Hill JG, Das A. Interplay among Electrostatic, Dispersion, and Steric Interactions: Spectroscopy and Quantum Chemical Calculations of π-Hydrogen Bonded Complexes. Chemphyschem 2017; 18:828-838. [PMID: 28124829 DOI: 10.1002/cphc.201601405] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2016] [Revised: 01/24/2017] [Indexed: 11/07/2022]
Abstract
π-Hydrogen bonding interactions are ubiquitous in both materials and biology. Despite their relatively weak nature, great progress has been made in their investigation by experimental and theoretical methods, but this becomes significantly more complicated when secondary intermolecular interactions are present. In this study, the effect of successive methyl substitution on the supramolecular structure and interaction energy of indole⋅⋅⋅methylated benzene (ind⋅⋅⋅n-mb, n=1-6) complexes is probed through a combination of supersonic jet experiments and benchmark-quality quantum chemical calculations. It is demonstrated that additional secondary interactions introduce a subtle interplay among electrostatic and dispersion forces, as well as steric repulsion, which fine-tunes the overall structural motif. Resonant two-photon ionization and IR-UV double-resonance spectroscopy techniques are used to probe jet-cooled ind⋅⋅⋅n-mb (n=2, 3, 6) complexes, with redshifting of the N-H IR stretching frequency showing that increasing the degree of methyl substitution increases the strength of the primary N-H⋅⋅⋅π interaction. Ab initio harmonic frequency and binding energy calculations confirm this trend for all six complexes. Electronic spectra of the three dimers are broad and structureless, with quantum chemical calculations revealing that this is likely to be due to multiple tilted conformations of each dimer possessing similar stabilization energies.
Collapse
Affiliation(s)
- Sumit Kumar
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Dr. Homi Bhabha Road, Pashan, Pune, 411008, Maharashtra, India.,Present address: Department of Dynamics at Surfaces, Max Planck Institute of Biophysical Chemistry, Am Faßberg 11, 37077, Gottingen, Germany
| | - Santosh K Singh
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Dr. Homi Bhabha Road, Pashan, Pune, 411008, Maharashtra, India
| | - Jamuna K Vaishnav
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Dr. Homi Bhabha Road, Pashan, Pune, 411008, Maharashtra, India.,Present address: Indian Institute of Technology (IIT) Indore, Khandwa Rd, Simrol, Madhya, Pradesh, 452020, India
| | - J Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield, S3 7HF, UK
| | - Aloke Das
- Department of Chemistry, Indian Institute of Science Education and Research (IISER), Dr. Homi Bhabha Road, Pashan, Pune, 411008, Maharashtra, India
| |
Collapse
|
142
|
Beker W, van der Kamp MW, Mulholland AJ, Sokalski WA. Rapid Estimation of Catalytic Efficiency by Cumulative Atomic Multipole Moments: Application to Ketosteroid Isomerase Mutants. J Chem Theory Comput 2017; 13:945-955. [PMID: 28103023 DOI: 10.1021/acs.jctc.6b01131] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We propose a simple atomic multipole electrostatic model to rapidly evaluate the effects of mutation on enzyme activity and test its performance on wild-type and mutant ketosteroid isomerase. The predictions of our atomic multipole model are similar to those obtained with symmetry-adapted perturbation theory at a fraction of the computational cost. We further show that this approach is relatively insensitive to the precise amino acid side chain conformation in mutants and may thus be useful in computational enzyme (re)design.
Collapse
Affiliation(s)
- Wiktor Beker
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wrocław University of Science and Technology , Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Marc W van der Kamp
- School of Biochemistry, Biomedical Sciences Building, University Walk , Bristol BS8 1TD, United Kingdom.,BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, Tyndall Avenue, University of Bristol , Bristol BS8 1TQ, United Kingdom.,Centre of Computational Chemistry, School of Chemistry, Cantock's Close, University of Bristol , Bristol BS8 1TS, United Kingdom
| | - Adrian J Mulholland
- BrisSynBio Synthetic Biology Research Centre, Life Sciences Building, Tyndall Avenue, University of Bristol , Bristol BS8 1TQ, United Kingdom.,Centre of Computational Chemistry, School of Chemistry, Cantock's Close, University of Bristol , Bristol BS8 1TS, United Kingdom
| | - W Andrzej Sokalski
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wrocław University of Science and Technology , Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|
143
|
Stasyuk OA, Jakubec D, Vondrášek J, Hobza P. Noncovalent Interactions in Specific Recognition Motifs of Protein-DNA Complexes. J Chem Theory Comput 2017; 13:877-885. [PMID: 27992205 DOI: 10.1021/acs.jctc.6b00775] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
In view of the importance of protein-DNA interactions in biological processes, we extracted from the Protein Data Bank several one-to-one complexes of amino acids with nucleotides that matched certain geometric and energetic specificity criteria and investigated them using quantum chemistry methods. The CCSD(T)/CBS interaction energies were used as a benchmark to compare the performance of the MP2.5, MP2-F12, DFT-D3, and PM6-D3H4 methods. All methods yielded good agreement with the reference values, with declining accuracy from MP2.5 to PM6-D3H4. Regardless of the site of interaction, the minima found after full optimization in implicit solvent with high dielectric constant were close to the structures experimentally detected in protein-DNA complexes. According to DFT-SAPT analysis, the nature of noncovalent interactions strongly depends on the type of amino acid. The negatively charged sugar-phosphate backbone of DNA heavily influences the strength of interactions and must be included in the computational model, especially in the case of interactions with charged amino acids.
Collapse
Affiliation(s)
- Olga A Stasyuk
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - David Jakubec
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University in Prague , Albertov 6, 128 43 Prague, Czech Republic
| | - Jiří Vondrášek
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic , Flemingovo nám. 2, 166 10 Prague, Czech Republic.,Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Palacký University , 771 46 Olomouc, Czech Republic
| |
Collapse
|
144
|
Blagojević JP, Veljković DŽ, Zarić SD. Stacking interactions between hydrogen-bridged and aromatic rings: study of crystal structures and quantum chemical calculations. CrystEngComm 2017. [DOI: 10.1039/c6ce02045c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
|
145
|
Harczuk I, Nagy B, Jensen F, Vahtras O, Ågren H. Local decomposition of imaginary polarizabilities and dispersion coefficients. Phys Chem Chem Phys 2017; 19:20241-20250. [DOI: 10.1039/c7cp02399e] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present a new way to compute the two-body contribution to the dispersion energy using ab initio theory.
Collapse
Affiliation(s)
- Ignat Harczuk
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| | - Balazs Nagy
- Department of Chemistry
- Aarhus University
- Langelandsgade 140
- 8000 Aarhus C
- Denmark
| | - Frank Jensen
- Department of Chemistry
- Aarhus University
- Langelandsgade 140
- 8000 Aarhus C
- Denmark
| | - Olav Vahtras
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| | - Hans Ågren
- KTH Royal Institute of Technology
- School of Biotechnology
- Division of Theoretical Chemistry and Biology
- SE-106 91 Stockholm
- Sweden
| |
Collapse
|
146
|
Vandenbrande S, Waroquier M, Speybroeck VV, Verstraelen T. The Monomer Electron Density Force Field (MEDFF): A Physically Inspired Model for Noncovalent Interactions. J Chem Theory Comput 2016; 13:161-179. [DOI: 10.1021/acs.jctc.6b00969] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Steven Vandenbrande
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Michel Waroquier
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| | - Toon Verstraelen
- Center for Molecular Modeling
(CMM), QCMM Ghent−Brussels Alliance, Ghent University, Technologiepark
903, B9000 Ghent, Belgium
| |
Collapse
|
147
|
Sherman MC, Ams MR, Jordan KD. Symmetry-Adapted Perturbation Theory Energy Analysis of Alkyl Fluorine-Aromatic Interactions in Torsion Balance Systems. J Phys Chem A 2016; 120:9292-9298. [PMID: 27797186 DOI: 10.1021/acs.jpca.6b09193] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Symmetry-adapted perturbation theory (SAPT) calculations are carried out to elucidate the intermolecular interactions present between fluorinated and nonfluorinated alkyl chain groups and aromatic π systems in the folded and unfolded conformers of Wilcox torsion balance systems. The calculations predict the folded conformers to be 2.0-2.3 kcal/mol more stable than the unfolded conformers, with the preference for the folded conformer being greater in the fluorinated alkyl chain case. We also establish that a simple electrostatic analysis, based on atomic charges, is inadequate for understanding the conformational preferences of these systems. In the folded conformers, there are sizable charge penetration contributions that are not recovered by point charge models. Additionally, the SAPT analysis reveals that exchange-repulsion interactions make a significant contribution to the relative stability of the folded and unfolded conformer.
Collapse
Affiliation(s)
- Mary C Sherman
- Department of Chemistry, University of Pittsburgh , Pittsburgh, Pennsylvania 15260, United States
| | - Mark R Ams
- Department of Chemistry, Allegheny College , Meadville, Pennsylvania 16335, United States
| | - Kenneth D Jordan
- Department of Chemistry, University of Pittsburgh , Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
148
|
Bagwill C, Anderson C, Sullivan E, Manohara V, Murthy P, Kirkpatrick CC, Stalcup A, Lewis M. Predicting the Strength of Anion−π Interactions of Substituted Benzenes: the Development of Anion−π Binding Substituent Constants. J Phys Chem A 2016; 120:9235-9243. [DOI: 10.1021/acs.jpca.6b06276] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Christina Bagwill
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Christa Anderson
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Elizabeth Sullivan
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Varun Manohara
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Prithvi Murthy
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Charles C. Kirkpatrick
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| | - Apryll Stalcup
- Irish
Separation Science Cluster, National Centre for Sensor Research, Dublin City University, Glasnevin, Dublin 9, Ireland
| | - Michael Lewis
- Department
of Chemistry, Saint Louis University, 3501 Laclede Avenue, St. Louis, Missouri 63103, United States
| |
Collapse
|
149
|
Qi R, Wang Q, Ren P. General van der Waals potential for common organic molecules. Bioorg Med Chem 2016; 24:4911-4919. [PMID: 27519463 DOI: 10.1016/j.bmc.2016.07.062] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2016] [Revised: 07/27/2016] [Accepted: 07/28/2016] [Indexed: 11/18/2022]
Abstract
This work presents a systematic development of a new van der Waals potential (vdW2016) for common organic molecules based on symmetry-adapted perturbation theory (SAPT) energy decomposition. The Buf-14-7 function, as well as Cubic-mean and Waldman-Hagler mixing rules were chosen given their best performance among other popular potentials. A database containing 39 organic molecules and 108 dimers was utilized to derive a general set of vdW parameters, which were further validated on nucleobase stacking systems and testing organic dimers. The vdW2016 potential is anticipated to significantly improve the accuracy and transferability of new generations of force fields for organic molecules.
Collapse
Affiliation(s)
- Rui Qi
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, United States
| | - Qiantao Wang
- West China School of Pharmacy, Sichuan University, Chengdu 610041, China
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX 78712, United States.
| |
Collapse
|
150
|
Li MM, Wang YB, Zhang Y, Wang W. The Nature of the Noncovalent Interactions between Benzene and C60 Fullerene. J Phys Chem A 2016; 120:5766-72. [DOI: 10.1021/acs.jpca.6b06492] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ming-Ming Li
- Department
of Chemistry and Key Laboratory of Guizhou High Performance Computational
Chemistry, Guizhou University, Guiyang 550025, P.R. China
| | - Yi-Bo Wang
- Department
of Chemistry and Key Laboratory of Guizhou High Performance Computational
Chemistry, Guizhou University, Guiyang 550025, P.R. China
| | - Yu Zhang
- College
of Chemistry and Chemical Engineering and Henan Key Laboratory of
Function-Oriented Porous Materials, Luoyang Normal University, Luoyang 471934, P.R. China
| | - Weizhou Wang
- College
of Chemistry and Chemical Engineering and Henan Key Laboratory of
Function-Oriented Porous Materials, Luoyang Normal University, Luoyang 471934, P.R. China
| |
Collapse
|