1
|
Abstract
We review different models for introducing electric polarization in force fields, with special focus on methods where polarization is modelled at the atomic charge level. While electric polarization has been included in several force fields, the common approach has been to focus on atomic dipole polarizability. Several approaches allow modelling electric polarization by using charge-flow between charge sites instead, but this has been less exploited, despite that atomic charges and charge-flow is expected to be more important than atomic dipoles and dipole polarizability. A number of challenges are required to be solved for charge-flow models to be incorporated into polarizable force fields, for example how to parameterize the models and how to make them computational efficient.
Collapse
Affiliation(s)
- Frank Jensen
- Department of Chemistry, Aarhus University, Denmark.
| |
Collapse
|
2
|
Munárriz J, Laplaza R, Martín Pendás A, Contreras-García J. A first step towards quantum energy potentials of electron pairs. Phys Chem Chem Phys 2019; 21:4215-4223. [PMID: 30747171 DOI: 10.1039/c8cp07509c] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
A first step towards the construction of a quantum force field for electron pairs in direct space is taken. Making use of topological tools (Interacting Quantum Atoms and the Electron Localisation Function), we have analysed the dependency of electron pairs electrostatic, kinetic and exchange-correlation energies upon bond stretching. Simple correlations were found, and can be explained with elementary models such as the homogeneous electron gas. The resulting energy model is applicable to various bonding regimes: from homopolar to highly polarized and even to non-conventional bonds. Overall, this is a fresh approach for developing real space-based force fields including an exchange-correlation term. It provides the relative weight of each of the contributions, showing that, in common Lewis structures, the exchange correlation contribution between electron pairs is negligible. However, our results reveal that classical approximations progressively fail for delocalised electrons, including lone pairs. This theoretical framework justifies the success of the classic Bond Charge Model (BCM) approach in solid state systems and sets the basis of its limits. Finally, this approach opens the door towards the development of quantitative rigorous energy models based on the ELF topology.
Collapse
Affiliation(s)
- Julen Munárriz
- Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, 50009, Zaragoza, Spain.
| | | | | | | |
Collapse
|
3
|
El Khoury L, Naseem-Khan S, Kwapien K, Hobaika Z, Maroun RG, Piquemal JP, Gresh N. Importance of explicit smeared lone-pairs in anisotropic polarizable molecular mechanics. Torture track angular tests for exchange-repulsion and charge transfer contributions. J Comput Chem 2017; 38:1897-1920. [DOI: 10.1002/jcc.24830] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Revised: 03/18/2017] [Accepted: 04/03/2017] [Indexed: 12/15/2022]
Affiliation(s)
- 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
| | - Sehr Naseem-Khan
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC; UMR7616 CNRS Paris France
| | - Karolina Kwapien
- Chemistry and Biology, Nucleo(s)tides and Immunology for Therapy (CBNIT); UMR 8601 CNRS, UFR Biomédicale Paris France
- Institut Charles-Gerhardt, UMR 5253, CNRS-UM2-UM1-ENSM; Montpellier France
| | - 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; Texas 78712
| | - Nohad Gresh
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC; UMR7616 CNRS Paris France
| |
Collapse
|
4
|
Varfolomeeva VV, Terentev AV. Weak hydrogen bonds in adsorption of nonrigid molecules on graphitized thermal carbon black. J STRUCT CHEM+ 2017. [DOI: 10.1134/s0022476617030180] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
5
|
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
|
6
|
Reinscheid F, Reinscheid U. Stereochemical analysis of menthol and menthylamine isomers using calculated and experimental optical rotation data. J Mol Struct 2016. [DOI: 10.1016/j.molstruc.2015.09.013] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
7
|
O’Meara MJ, Leaver-Fay A, Tyka M, Stein A, Houlihan K, DiMaio F, Bradley P, Kortemme T, Baker D, Snoeyink J, Kuhlman B. Combined covalent-electrostatic model of hydrogen bonding improves structure prediction with Rosetta. J Chem Theory Comput 2015; 11:609-22. [PMID: 25866491 PMCID: PMC4390092 DOI: 10.1021/ct500864r] [Citation(s) in RCA: 169] [Impact Index Per Article: 18.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Interactions between polar atoms are challenging to model because at very short ranges they form hydrogen bonds (H-bonds) that are partially covalent in character and exhibit strong orientation preferences; at longer ranges the orientation preferences are lost, but significant electrostatic interactions between charged and partially charged atoms remain. To simultaneously model these two types of behavior, we refined an orientation dependent model of hydrogen bonds [Kortemme et al. J. Mol. Biol. 2003, 326, 1239] used by the molecular modeling program Rosetta and then combined it with a distance-dependent Coulomb model of electrostatics. The functional form of the H-bond potential is physically motivated and parameters are fit so that H-bond geometries that Rosetta generates closely resemble H-bond geometries in high-resolution crystal structures. The combined potentials improve performance in a variety of scientific benchmarks including decoy discrimination, side chain prediction, and native sequence recovery in protein design simulations and establishes a new standard energy function for Rosetta.
Collapse
Affiliation(s)
- Matthew J. O’Meara
- Department of Computer Science, University of North Carolina, 201 S Columbia St. Chapel Hill, North Carolina 27599, United States
| | - Andrew Leaver-Fay
- Department of Biochemistry and Biophysics, University of North Carolina, 120 Mason Farm Rd Chapel Hill, North Carolina 27599, United States
| | - Mike Tyka
- Google Inc., 1600 Amphitheatre Parkway Mountain View, California 94043, United States
| | - Amelie Stein
- Department of Bioengineering and Therapeutic Science, University of California San Francisco, 513 Parnassus Avenue San Francisco, California 94143, United States
| | - Kevin Houlihan
- Department of Biochemistry and Biophysics, University of North Carolina, 120 Mason Farm Rd Chapel Hill, North Carolina 27599, United States
| | - Frank DiMaio
- Department of Biochemistry, University of Washington, 1705 North East Pacific Street Seattle Washington 98195, United States
| | - Philip Bradley
- Fred Hutchinson Cancer Research Center, 1100 Fairview Avenue North, Seattle Washington 98109, United States
| | - Tanja Kortemme
- Department of Bioengineering and Therapeutic Science, University of California San Francisco, 513 Parnassus Avenue San Francisco, California 94143, United States
| | - David Baker
- Department of Biochemistry, University of Washington, 1705 North East Pacific Street Seattle Washington 98195, United States
| | - Jack Snoeyink
- Department of Computer Science, University of North Carolina, 201 S Columbia St. Chapel Hill, North Carolina 27599, United States
| | - Brian Kuhlman
- Department of Biochemistry and Biophysics, University of North Carolina, 120 Mason Farm Rd Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
8
|
Taha HA, Richards MR, Lowary TL. Conformational Analysis of Furanoside-Containing Mono- and Oligosaccharides. Chem Rev 2012; 113:1851-76. [DOI: 10.1021/cr300249c] [Citation(s) in RCA: 99] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Hashem A. Taha
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| | - Michele R. Richards
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| | - Todd L. Lowary
- Alberta Glycomics Centre and Department of Chemistry, Gunning−Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2
| |
Collapse
|
9
|
Carter M, Rappé AK, Ho PS. Scalable Anisotropic Shape and Electrostatic Models for Biological Bromine Halogen Bonds. J Chem Theory Comput 2012; 8:2461-73. [DOI: 10.1021/ct3001969] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Megan Carter
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| | - Anthony K. Rappé
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| | - P. Shing Ho
- Department
of Biochemistry and Molecular Biology and ‡Department of Chemistry, Colorado State University, Fort Collins,
Colorado 80523, United States
| |
Collapse
|
10
|
Lu Z, Zhou N, Wu Q, Zhang Y. Directional Dependence of Hydrogen Bonds: a Density-based Energy Decomposition Analysis and Its Implications on Force Field Development. J Chem Theory Comput 2011; 7:4038-4049. [PMID: 22267958 PMCID: PMC3259744 DOI: 10.1021/ct2003226] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
One well-known shortcoming of widely-used biomolecular force fields is the description of the directional dependence of hydrogen bonding (HB). Here we aim to better understand the origin of this difficulty and thus provide some guidance for further force field development. Our theoretical approaches center on a novel density-based energy decomposition analysis (DEDA) method [J. Chem. Phys., 131, 164112 (2009)], in which the frozen density energy is variationally determined through constrained search. This unique and most significant feature of DEDA enables us to find that the frozen density interaction term is the key factor in determining the HB orientation, while the sum of polarization and charge-transfer components shows very little HB directional dependence. This new insight suggests that the difficulty for current non-polarizable force fields to describe the HB directional dependence is not due to the lack of explicit polarization or charge-transfer terms. Using the DEDA results as reference, we further demonstrate that the main failure coming from the atomic point charge model can be overcome largely by introducing extra charge sites or higher order multipole moments. Among all the electrostatic models explored, the smeared charge distributed multipole model (up to quadrupole), which also takes account of charge penetration effects, gives the best agreement with the corresponding DEDA results. Meanwhile, our results indicate that the van der Waals interaction term needs to be further improved to better model directional hydrogen bonding.
Collapse
Affiliation(s)
- Zhenyu Lu
- Department of Chemistry, New York University, New York 10003
| | - Nengjie Zhou
- Department of Chemistry, New York University, New York 10003
| | - Qin Wu
- Center for Functional Nanomaterials, Brookhaven National Laboratory, Upton, New York 11973
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York 10003
| |
Collapse
|
11
|
Taha HA, Roy PN, Lowary TL. Theoretical Investigations on the Conformation of the β-d-Arabinofuranoside Ring. J Chem Theory Comput 2010; 7:420-32. [DOI: 10.1021/ct100450s] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Hashem A. Taha
- Department of Chemistry and Alberta Ingenuity Centre for Carbohydrate Science, Gunning-Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2 and Department of Chemistry, University of Waterloo, Waterloo, ON, Canada N2L 3G1
| | - Pierre-Nicholas Roy
- Department of Chemistry and Alberta Ingenuity Centre for Carbohydrate Science, Gunning-Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2 and Department of Chemistry, University of Waterloo, Waterloo, ON, Canada N2L 3G1
| | - Todd L. Lowary
- Department of Chemistry and Alberta Ingenuity Centre for Carbohydrate Science, Gunning-Lemieux Chemistry Centre, University of Alberta, Edmonton, AB, Canada T6G 2G2 and Department of Chemistry, University of Waterloo, Waterloo, ON, Canada N2L 3G1
| |
Collapse
|
12
|
Selenoglycosides in silico: ab initio-derived reparameterization of MM4, conformational analysis using histo-blood group ABH antigens and lectin docking as indication for potential of bioactivity. J Comput Aided Mol Des 2010; 24:1009-21. [PMID: 20976527 DOI: 10.1007/s10822-010-9392-y] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2010] [Accepted: 10/06/2010] [Indexed: 10/18/2022]
Abstract
The identification of glycan epitopes such as the histo-blood group ABH determinants as docking sites for bacterial/viral infections and signals in growth regulation fuels the interest to develop non-hydrolysable mimetics for therapeutic applications. Inevitably, the required substitution of the linkage oxygen atom will alter the derivative's topology. Our study addresses the question of the impact of substitution of oxygen by selenium. In order to characterize spatial parameters and flexibility of selenoglycosides, we first performed ab initio calculations on model compounds to refine the MM4 force field. The following application of the resulting MM4R version appears to reduce the difference to ab initio data when compared to using the MM4 estimator. Systematic conformational searches on the derivatives of histo-blood group ABH antigens revealed increased flexibility with acquisition of additional low-energy conformer(s), akin to the behavior of S-glycosides. Docking analysis using the Glide program for eight test cases indicated potential for bioactivity, giving further experimental investigation a clear direction to testing Se-glycosides as lectin ligands.
Collapse
|
13
|
Stortz CA, Johnson GP, French AD, Csonka GI. Comparison of different force fields for the study of disaccharides. Carbohydr Res 2009; 344:2217-28. [PMID: 19758584 DOI: 10.1016/j.carres.2009.08.019] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2009] [Revised: 08/13/2009] [Accepted: 08/18/2009] [Indexed: 11/30/2022]
Abstract
Eighteen empirical force fields and the semi-empirical quantum method PM3CARB-1 were compared for studying beta-cellobiose, alpha-maltose, and alpha-galabiose [alpha-D-Galp-(1-->4)-alpha-D-Galp]. For each disaccharide, the energies of 54 conformers with differing hydroxymethyl, hydroxyl, and glycosidic linkage orientations were minimized by the different methods, some at two dielectric constants. By comparing these results and the available crystal structure data and/or higher level density functional theory results, it was concluded that the newer parameterizations for force fields (GROMOS, GLYCAM06, OPLS-2005 and CSFF) give results that are reasonably similar to each other, whereas the older parameterizations for Amber, CHARMM or OPLS were more divergent. However, MM3, an older force field, gave energy and geometry values comparable to those of the newer parameterizations, but with less sensitivity to dielectric constant values. These systems worked better than MM2 variants, which were still acceptable. PM3CARB-1 also gave adequate results in terms of linkage and exocyclic torsion angles. GROMOS, GLYCAM06, and MM3 appear to be the best choices, closely followed by MM4, CSFF, and OPLS-2005. With GLYCAM06 and to a lesser extent, CSFF, and OPLS-2005, a number of the conformers that were stable with MM3 changed to other forms.
Collapse
Affiliation(s)
- Carlos A Stortz
- Departamento de Química Orgánica-CIHIDECAR, FCEyN-Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina.
| | | | | | | |
Collapse
|