1
|
Gómez S, Cappelli C. The Role of Hydrogen Bonding in the Raman Spectral Signals of Caffeine in Aqueous Solution. Molecules 2024; 29:3035. [PMID: 38998986 PMCID: PMC11243038 DOI: 10.3390/molecules29133035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2024] [Revised: 06/22/2024] [Accepted: 06/23/2024] [Indexed: 07/14/2024] Open
Abstract
The identification and quantification of caffeine is a common need in the food and pharmaceutical industries and lately also in the field of environmental science. For that purpose, Raman spectroscopy has been used as an analytical technique, but the interpretation of the spectra requires reliable and accurate computational protocols, especially as regards the Resonance Raman (RR) variant. Herein, caffeine solutions are sampled using Molecular Dynamics simulations. Upon quantification of the strength of the non-covalent intermolecular interactions such as hydrogen bonding between caffeine and water, UV-Vis, Raman, and RR spectra are computed. The results provide general insights into the hydrogen bonding role in mediating the Raman spectral signals of caffeine in aqueous solution. Also, by analyzing the dependence of RR enhancement on the absorption spectrum of caffeine, it is proposed that the sensitivity of the RR technique could be exploited at excitation wavelengths moderately far from 266 nm, yet achieving very low detection limits in the quantification caffeine content.
Collapse
Affiliation(s)
- Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
2
|
Liu C, Lv C, Yao YY, Du X, Zhao DX, Yang ZZ. Water-Mediated Oxidation of Guanine by a Repair Enzyme: Simulation Using the ABEEM Polarizable Force Field. J Chem Theory Comput 2021; 17:3525-3538. [PMID: 34018392 DOI: 10.1021/acs.jctc.1c00107] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The recognition mechanism of oxidative damage in organisms has long been a research hotspot. Water is an important medium in the recognition process, but its specific role remains unknown. There is a need to develop a suitable force field that can adequately describe the electrostatic, hydrogen bond, and other interactions among the molecules in the complex system of the repair enzyme and oxidized base. The developing ABEEM polarizable force field (PFF) has been used to simulate the repaired enzyme hOGG1 and oxidized DNA (PDB ID: 1EBM) in a biological environment, and the corresponding results are better than those of the fixed-charge force fields OPLS/AA and AMBER OL15. 8-Oxo-G is recognized by Gln315 of hOGG1 mainly through hydrogen bonds mediated by continuous exchange of 2 water molecules. Phe319 and Cys253 are stacked on both sides of the π planes of bases to form sandwich structures. The charge polarization effect gives an important signal to drive the exchange of water molecules and maintains the recognition of oxidation bases by enzymes. The mediated main water molecule A and mediated auxiliary water molecule B together pull Gln315 to recognize 8-oxo-G by hydrogen bond interactions. Then, the charge polarization signal of solvent water molecule C with a large absolute charge causes the absolute charge of O atoms in water molecule A or B to increase by approximately 0.2 e, and water molecule A or B leaves Gln315 and 8-oxo-G. The other water molecule and water molecule C synergistically recognize 8-oxo-G with Gln315. Even though the water molecules between Gln315 and 8-oxo-G are removed, the MD simulation results show that water molecules appear between Gln315 and 8-oxo-G in a very short time (<2 ps). The dwell time of each water molecule is approximately 60 ps. The radial distribution function and dwell time support the correctness of the above mechanism. These polarization effects and hydrogen bonding interactions cannot be simulated by a fixed-charge force field.
Collapse
Affiliation(s)
- Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Change Lv
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Yu-Ying Yao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Xue Du
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, China
| |
Collapse
|
3
|
Raynal M, Li Y, Troufflard C, Przybylski C, Gontard G, Maistriaux T, Idé J, Lazzaroni R, Bouteiller L, Brocorens P. Experimental and computational diagnosis of the fluxional nature of a benzene-1,3,5-tricarboxamide-based hydrogen-bonded dimer. Phys Chem Chem Phys 2021; 23:5207-5221. [PMID: 33625418 DOI: 10.1039/d0cp06128j] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Precise characterization of the hydrogen bond network present in discrete self-assemblies of benzene-1,3,5-tricarboxamide monomers derived from amino-esters (ester BTAs) is crucial for the construction of elaborated functional co-assemblies. For all ester BTA dimeric structures previously reported, ester carbonyls in the side chain acted as hydrogen bond acceptors, yielding well-defined dimers stabilized by six hydrogen bonds. The ester BTA monomer derived from glycine (BTA Gly) shows a markedly different self-assembly behaviour. We report herein a combined experimental and computational investigation aimed at determining the nature of the dimeric species formed by BTA Gly. Two distinct dimeric structures were characterized by single-crystal X-ray diffraction measurements. Likewise, a range of spectroscopic and scattering techniques as well as molecular modelling were employed to diagnose the nature of dynamic dimeric structures in toluene. Our results unambiguously establish that both ester and amide carbonyls are involved in the hydrogen bond network of the discrete dimeric species formed by BTA Gly. The participation of roughly 4.5 ester carbonyls and 1.5 amide carbonyls per dimer as determined by FT-IR spectroscopy implies that several conformations coexist in solution. Moreover, NMR analysis and modelling data reveal rapid interconversion between these different conformers leading to a symmetric structure on the NMR timescale. Rapid hydrogen bond shuffling between conformers having three (three), two (four), one (five) and zero (six) amide carbonyl groups (ester carbonyl groups, respectively) as hydrogen bond acceptors is proposed to explain the magnetic equivalence of the amide N-H on the NMR timescale. When compared to other ester BTA derivatives in which only ester carbonyls act as hydrogen bond acceptors, the fluxional behaviour of the hydrogen-bonded dimers of BTA Gly likely originates from a larger range of energetically favorable conformations accessible through rotation of the BTA side chains.
Collapse
Affiliation(s)
- M Raynal
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - Y Li
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - C Troufflard
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - C Przybylski
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - G Gontard
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - T Maistriaux
- Service de Chimie des Matériaux Nouveaux, Institut de Recherche sur les Matériaux, Université de Mons, Place du Parc, 20, B-7000, Mons, Belgium.
| | - J Idé
- Service de Chimie des Matériaux Nouveaux, Institut de Recherche sur les Matériaux, Université de Mons, Place du Parc, 20, B-7000, Mons, Belgium.
| | - R Lazzaroni
- Service de Chimie des Matériaux Nouveaux, Institut de Recherche sur les Matériaux, Université de Mons, Place du Parc, 20, B-7000, Mons, Belgium.
| | - L Bouteiller
- Sorbonne Université, CNRS, IPCM, UMR 8232, 4 Place Jussieu, 75252 Paris Cedex 05, France.
| | - P Brocorens
- Service de Chimie des Matériaux Nouveaux, Institut de Recherche sur les Matériaux, Université de Mons, Place du Parc, 20, B-7000, Mons, Belgium.
| |
Collapse
|
4
|
Shi H, Gong LD, Liu C, Lu LN, Yang ZZ. ABEEM/MM OH - Models for OH -(H 2O) n Clusters and Aqueous OH -: Structures, Charge Distributions, and Binding Energies. J Phys Chem A 2020; 124:5963-5978. [PMID: 32520555 DOI: 10.1021/acs.jpca.0c03941] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Based on the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM), two fluctuating charge models of OH--water system were proposed. The difference between these two models is whether there is charge transfer between OH- and its first-shell water molecules. The structures, charge distributions, charge transfer, and binding energies of the OH-(H2O)n (n = 1-8, 10, 15, 23) clusters were studied by these two ABEEM/MM models, the OPLS/AA force field, the OPLS-SMOOTH/AA force field, and the QM methods. The results demonstrate that two ABEEM/MM models can search out all stable structures just as the QM methods, and the structures and charge distributions agree well with those from the QM calculations. The structures, the charge transfer, and the strength of hydrogen bonds in the first hydration shell are closely related to the coordination number of OH-. Molecular dynamics simulations on the aqueous OH- solution are performed at 298 and 278 K using ABEEM/MM-I model. The MD results show that the populations of three-, four-, and five-coordinated OH- are 29.6%, 67.1%, and 3.4% at 298 K, respectively, and those of two-, three-, four-, and five-coordinated OH- are 10.8%, 44.9%, 39.2%, and 4.9% at 278 K, respectively; the average hydrogen bond lengths and the hydrogen bond angle in the first shell increase with the temperature decreasing.
Collapse
Affiliation(s)
- Hua Shi
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China.,School of Marine Science and Environment, Dalian Ocean University, Dalian 116023, People's Republic of China
| | - Li-Dong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian, 116029, People's Republic of China
| |
Collapse
|
5
|
Li H, Wang D, Zhao X, Lu LN, Liu C, Gong LD, Zhao DX, Yang ZZ. Reaction mechanism of NO with hydrolysates of NAMI-A: an MD simulation by combining the QM/MM(ABEEM) with the MD-FEP method. J Comput Chem 2019; 40:1141-1150. [PMID: 30375671 DOI: 10.1002/jcc.25734] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Revised: 09/27/2018] [Accepted: 09/28/2018] [Indexed: 11/06/2022]
Abstract
Nitrosylation reaction mechanisms of the hydrolysates of NAMI-A and hydrolysis reactions of ruthenium nitrosyl complexes were investigated in the triplet state and the singlet state. Activation free energies were calculated by combining the QM/MM(ABEEM) method with free energy perturbation theory, and the explicit solvent environment was simulated by an ABEEMσπ polarizable force field. Our results demonstrate that nitrosylation reactions of the hydrolysates of NAMI-A occur in both the triplet and the singlet states. The Ru-N-O angle of the triplet ruthenium nitrosyl complexes is in the range of 132.0°-138.2°. However, all the ruthenium nitrosyl complexes at the singlet state show an almost linear Ru-N-O angle. The nitrosylation reaction happens prior to the hydrolysis reaction for the first-step hydrolysates. The activation free energies of the nitrosylation reactions show that the H2 O-NO exchange reaction of [RuCl4 (Im)(H2 O)] in the singlet spin sate is the most likely one. Comparing with the activation free energies of the hydrolysis reactions of the ruthenium nitrosyl complexes, the results indicate that the rate of the DMSO-H2 O exchange reaction of [RuCl3 (NO)(Im)(DMSO)] is faster than that of [RuCl3 (H2 O)(Im)(DMSO)] in both the triplet spin state and the singlet spin state. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hui Li
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China.,Department of Chemistry, Bohai University, Jinzhou 121013, China
| | - Di Wang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Xin Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Li-Nan Lu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Cui Liu
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Li-Dong Gong
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| | - Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering, Liaoning Normal University, Dalian 116029, People's Republic of China
| |
Collapse
|
6
|
Liu C, Li Y, Han BY, Gong LD, Lu LN, Yang ZZ, Zhao DX. Development of the ABEEMσπ Polarization Force Field for Base Pairs with Amino Acid Residue Complexes. J Chem Theory Comput 2017; 13:2098-2111. [DOI: 10.1021/acs.jctc.6b01206] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Cui Liu
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Yue Li
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Bing-Yu Han
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Li-Dong Gong
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Li-Nan Lu
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Zhong-Zhi Yang
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| | - Dong-Xia Zhao
- School of Chemistry
and Chemical
Engineering, Liaoning Normal University, Dalian 116029, China
| |
Collapse
|
7
|
Tuttle T. Computational Approaches to Understanding the Self-assembly of Peptide-based Nanostructures. Isr J Chem 2015. [DOI: 10.1002/ijch.201400188] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
|
8
|
Wu Y, Hu N, Yue L, Wei L, Guan W. Effects of polarizability on the structural and thermodynamics properties of [Cnmim][Gly] ionic liquids (n = 1–4) using EEM/MM molecular dynamic simulations. J Chem Phys 2015; 142:064503. [DOI: 10.1063/1.4907281] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Yang Wu
- Key Laboratory of Green Synthesis and Preparative Chemistry of Advanced Materials, College of Chemistry, Liaoning University, Shenyang 110036, China
| | - Na Hu
- Key Laboratory of Green Synthesis and Preparative Chemistry of Advanced Materials, College of Chemistry, Liaoning University, Shenyang 110036, China
| | - Lili Yue
- Key Laboratory of Green Synthesis and Preparative Chemistry of Advanced Materials, College of Chemistry, Liaoning University, Shenyang 110036, China
| | - Lihong Wei
- Key Laboratory of Green Synthesis and Preparative Chemistry of Advanced Materials, College of Chemistry, Liaoning University, Shenyang 110036, China
| | - Wei Guan
- Key Laboratory of Green Synthesis and Preparative Chemistry of Advanced Materials, College of Chemistry, Liaoning University, Shenyang 110036, China
| |
Collapse
|
9
|
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
|
10
|
Yang ZZ, Wang JJ, Zhao DX. Valence state parameters of all transition metal atoms in metalloproteins-development of ABEEMσπ fluctuating charge force field. J Comput Chem 2014; 35:1690-706. [DOI: 10.1002/jcc.23676] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Revised: 06/01/2014] [Accepted: 06/18/2014] [Indexed: 01/08/2023]
Affiliation(s)
- Zhong-Zhi Yang
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| | - Jian-Jiang Wang
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| | - Dong-Xia Zhao
- School of Chemistry and Chemical Engineering; Department of Chemistry, Liaoning Normal University; Dalian China 116029
| |
Collapse
|
11
|
Development of ABEEMσπ polarizable force field for oxidized adenine base pairs: investigation of the interaction and mutagenic mechanism. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1469-6] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
12
|
Development and applications of the ABEEM fluctuating charge molecular force field in the ion-containing systems. Sci China Chem 2012. [DOI: 10.1007/s11426-012-4787-3] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
|