1
|
Porwal VK, Carof A, Ingrosso F. Hydration effects on the vibrational properties of carboxylates: From continuum models to QM/MM simulations. J Comput Chem 2023. [PMID: 37300426 DOI: 10.1002/jcc.27171] [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: 09/12/2022] [Revised: 05/15/2023] [Accepted: 05/21/2023] [Indexed: 06/12/2023]
Abstract
The presence of carboxyl groups in a molecule delivers an affinity to metal cations and a sensitivity to the chemical environment, especially for an environment that can give rise to intermolecular hydrogen bonds. Carboxylate groups can also induce intramolecular interactions, such as the formation of hydrogen bonds with donor groups, leading to an impact on the conformational space of biomolecules. In the latter case, the protonation state of the amino groups plays an important role. In order to provide an accurate description of the modifications induced in a carboxylated molecule by the formation of hydrogen bonds, one needs a compromise between a quantum chemical description of the system and the necessity to take into account explicit solvent molecules. In this work, we propose a bottom-up approach to study the conformational space and the carboxylate stretching band of (bio)organic anions. Starting from the anions in a continuum solvent, we then move to calculations using a microsolvation approach including one explicit water molecule per polar group, immersed in a continuum. Finally, we run QM/MM molecular dynamics simulations to analyze the solvation properties and to explore the anions conformational space. The results thus obtained are in good agreement with the description given by the microsolvation approach and they bring a more detailed description of the solvation shell and of the intermolecular hydrogen bonds.
Collapse
Affiliation(s)
- Vishal Kumar Porwal
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques UMR 7019, Nancy, France
| | - Antoine Carof
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques UMR 7019, Nancy, France
| | - Francesca Ingrosso
- Université de Lorraine and CNRS, Laboratoire de Physique et Chimie Théoriques UMR 7019, Nancy, France
| |
Collapse
|
2
|
Abstract
A modified neglect of differential overlap has been parameterized specifically for water and its oligomers with the addition of polarization functions on both hydrogen and oxygen, Feynman dispersion, and a slight modification of the treatment of the hydrogen nucleus. The results show that it is possible to easily obtain good geometries and energies for hydrogen-bonded water aggregates. Data from the Benchmark Energy and Geometry Database water-cluster database were used to parameterize the new Hamiltonian for water clusters from the dimer to the decamer using MP2/aug-cc-pVDZ optimized geometries and CCSD(T)/CBS oligomerization energies. Seventy five oligomerization and rearrangement energies derived from the parameterization data are reproduced with a root mean-square error (RMSE) of 0.79 kcal mol-1 and the geometries of 38 oligomers with an RMSE of 0.17 Å. Interestingly, the Feynman dispersion term adopts a role different from that intended and tunes the atomic polarizability. The implications of these results in terms of future dedicated neglect of diatomic differential overlap Hamiltonians and those that use force-field-like atom types are discussed.
Collapse
Affiliation(s)
- Matthias Hennemann
- Computer-Chemistry-Center, Department of Chemistry and Pharmacy, Friedrich-Alexander-University Erlangen-Nuernberg, Naegelsbachstr. 25, 91052 Erlangen, Germany
| | - Timothy Clark
- Computer-Chemistry-Center, Department of Chemistry and Pharmacy, Friedrich-Alexander-University Erlangen-Nuernberg, Naegelsbachstr. 25, 91052 Erlangen, Germany
| |
Collapse
|
3
|
Examination of the performance of semiempirical methods in QM/MM studies of the SN2-like reaction of an adenylyl group transfer catalysed by ANT4′. Theor Chem Acc 2019. [DOI: 10.1007/s00214-019-2507-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
4
|
Marion A, Gokcan H, Monard G. Semi-Empirical Born-Oppenheimer Molecular Dynamics (SEBOMD) within the Amber Biomolecular Package. J Chem Inf Model 2019; 59:206-214. [PMID: 30433776 DOI: 10.1021/acs.jcim.8b00605] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Semi-empirical quantum methods from the neglect of differential diatomic overlap (NDDO) family such as MNDO, AM1, or PM3 are fast albeit approximate quantum methods. By combining them with linear scaling methods like the divide & conquer (D&C) method, it is possible to quickly evaluate the energy of systems containing hundreds to thousands of atoms. We here present our implementation in the Amber biomolecular package of a SEBOMD module that provides a way to run semi-empirical Born-Oppenheimer molecular dynamics. At each step of a SEBOMD, a fully converged self-consistent field (SCF) calculation is performed to obtain the semiempirical quantum potential energy of a molecular system encaged or not in periodic boundary conditions. We describe the implementation and the features of our SEBOMD implementation. We show the requirements to conserve the total energy in NVE simulations, and how to accelerate SCF convergence through density matrix extrapolation. Specific ways of handling periodic boundary conditions using mechanical embedding or electrostatic embedding through a tailored quantum Ewald summation is developed. The parallel performance of SEBOMD simulations using the D&C scheme are presented for liquid water systems of various sizes, and a comparison between the traditional full diagonalization scheme and the D&C approach for the reproduction of the structure of liquid water illustrates the potentiality of SEBOMD to simulate molecular systems containing several hundreds of atoms for hundreds of picoseconds with a quantum mechanical potential in a reasonable amount of CPU time.
Collapse
Affiliation(s)
- Antoine Marion
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France.,Department of Chemistry , Middle East Technical University , 06800 , Ankara , Turkey
| | - Hatice Gokcan
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France.,Department of Chemistry , University of North Texas , Denton , Texas 76201 , United States
| | - Gerald Monard
- Université de Lorraine, CNRS, LPCT , F-54000 Nancy , France
| |
Collapse
|
5
|
Kartashynska E, Vysotsky Y, Fainerman V, Vollhardt D, Miller R. Quantum-chemical analysis of condensed monolayer phases of N-alkanoyl-substituted alanine at the air/water interface. Colloids Surf A Physicochem Eng Asp 2018. [DOI: 10.1016/j.colsurfa.2018.03.027] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
6
|
Blancafort L, Voityuk AA. Thermally induced hopping model for long-range triplet excitation energy transfer in DNA. Phys Chem Chem Phys 2018; 20:4997-5000. [DOI: 10.1039/c7cp07811k] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Long-range triplet excitation energy transfer in DNA occurs through hops between T bases, providing it with a shallow distance dependence. Hops over two thymines on the same strand are preferred.
Collapse
Affiliation(s)
- Lluís Blancafort
- Institut de Química Computacional i Catàlisi and Departament de Química
- Universitat de Girona
- Facultat de Ciències
- Girona 17003
- Spain
| | - Alexander A. Voityuk
- Institut de Química Computacional i Catàlisi and Departament de Química
- Universitat de Girona
- Facultat de Ciències
- Girona 17003
- Spain
| |
Collapse
|
7
|
Marion A, Góra J, Kracker O, Fröhr T, Latajka R, Sewald N, Antes I. Amber-Compatible Parametrization Procedure for Peptide-like Compounds: Application to 1,4- and 1,5-Substituted Triazole-Based Peptidomimetics. J Chem Inf Model 2017; 58:90-110. [PMID: 29112399 DOI: 10.1021/acs.jcim.7b00305] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Peptidomimetics are molecules of particular interest in the context of drug design and development. They are proteolytically and metabolically more stable than their natural peptide counterparts but still offer high specificity toward their biological targets. In recent years, 1,4- and 1,5-disubstituted 1,2,3-triazole-based peptidomimetics have emerged as promising lead compounds for the design of various inhibitory and tumor-targeting molecules as well as for the synthesis of peptide analogues. The growing popularity of triazole-based peptidomimetics and a constantly broadening range of their application generated a demand for elaborate theoretical investigations by classical molecular dynamics simulations and molecular docking. Despite this rising interest, accurate and coherent force field parameters for triazole-based peptidomimetics are still lacking. Here, we report the first complete set of parameters dedicated to this group of compounds, named TZLff. This parametrization is compatible with the latest version of the AMBER force field (ff14SB) and can be readily applied for the modeling of pure triazole-based peptidomimetics as well as natural peptide sequences containing one or more triazole-based modifications in their backbone. The parameters were optimized to reproduce HF/6-31G* electrostatic potentials as well as MP2/cc-pVTZ equilibrium Hessian matrices and conformational potential energy surfaces through the use of a genetic algorithm-based search and least-squares fitting. Following the standards of AMBER, we introduce residue building units, thus allowing the user to define any given sequence of triazole-based peptidomimetics. Validation of the parameter set against ab initio- and NMR-based reference systems shows that we obtain fairly accurate results, which properly capture the conformational features of triazole-based peptidomimetics. The successful and efficient parametrization strategy developed in this work is general enough to be applied in a straightforward manner for parametrization of other peptidomimetics and, potentially, any polymeric assemblies.
Collapse
Affiliation(s)
- Antoine Marion
- Center for Integrated Protein Science Munich at the Department of Biosciences, Technische Universität München , Emil-Erlenmeyer-Forum 8, D-85354 Freising, Germany
| | - Jerzy Góra
- Organic and Bioorganic Chemistry, Department of Chemistry, Bielefeld University , Universitätsstraße 25, D-33615 Bielefeld, Germany.,Department of Organic and Pharmaceutical Technology, Faculty of Chemistry, Wroclaw University of Science and Technology , Wybrzeze Wyspianskiego 27, PL-50-370 Wroclaw, Poland
| | - Oliver Kracker
- Organic and Bioorganic Chemistry, Department of Chemistry, Bielefeld University , Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Tanja Fröhr
- Organic and Bioorganic Chemistry, Department of Chemistry, Bielefeld University , Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Rafał Latajka
- Department of Organic and Pharmaceutical Technology, Faculty of Chemistry, Wroclaw University of Science and Technology , Wybrzeze Wyspianskiego 27, PL-50-370 Wroclaw, Poland
| | - Norbert Sewald
- Organic and Bioorganic Chemistry, Department of Chemistry, Bielefeld University , Universitätsstraße 25, D-33615 Bielefeld, Germany
| | - Iris Antes
- Center for Integrated Protein Science Munich at the Department of Biosciences, Technische Universität München , Emil-Erlenmeyer-Forum 8, D-85354 Freising, Germany
| |
Collapse
|
8
|
Alcaraz-Torres A, Gamboa-Suárez A, Bernal-Uruchurtu MI. Is Br2 hydration hydrophobic? J Chem Phys 2017; 146:084501. [DOI: 10.1063/1.4975688] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- A. Alcaraz-Torres
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Morelos, Avenida Universidad 1001, Cuernavaca, Morelos 62209,
México
| | - A. Gamboa-Suárez
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Morelos, Avenida Universidad 1001, Cuernavaca, Morelos 62209,
México
| | - M. I. Bernal-Uruchurtu
- Centro de Investigaciones Químicas, Universidad Autónoma del Estado de Morelos, Avenida Universidad 1001, Cuernavaca, Morelos 62209,
México
| |
Collapse
|
9
|
Theoretical description of 2D-cluster formation of nonionic surfactants at the air/water interface. Colloid Polym Sci 2015. [DOI: 10.1007/s00396-015-3630-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
10
|
Effects of monohydration on an adenine–thymine base pair. Theor Chem Acc 2015. [DOI: 10.1007/s00214-015-1686-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
11
|
Welborn M, Chen J, Wang LP, Van Voorhis T. Why many semiempirical molecular orbital theories fail for liquid water and how to fix them. J Comput Chem 2015; 36:934-9. [PMID: 25766721 DOI: 10.1002/jcc.23887] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2014] [Revised: 12/16/2014] [Accepted: 02/08/2015] [Indexed: 11/06/2022]
Abstract
Water is an extremely important liquid for chemistry and the search for more accurate force fields for liquid water continues unabated. Neglect of diatomic differential overlap (NDDO) molecular orbital methods provide and intriguing generalization of classical force fields in this regard because they can account both for bond breaking and electronic polarization of molecules. However, we show that most standard NDDO methods fail for water because they give an incorrect description of hydrogen bonding, water's key structural feature. Using force matching, we design a reparameterized NDDO model and find that it qualitatively reproduces the experimental radial distribution function of water, as well as various monomer, dimer, and bulk properties that PM6 does not. This suggests that the apparent limitations of NDDO models are primarily due to poor parameterization and not to the NDDO approximations themselves. Finally, we identify the physical parameters that most influence the condensed phase properties. These results help to elucidate the chemistry that a semiempirical molecular orbital picture of water must capture. We conclude that properly parameterized NDDO models could be useful for simulations that require electronically detailed explicit solvent, including the calculation of redox potentials and simulation of charge transfer and photochemistry.
Collapse
Affiliation(s)
- Matthew Welborn
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts, 02139
| | | | | | | |
Collapse
|
12
|
Yilmazer ND, Korth M. Enhanced semiempirical QM methods for biomolecular interactions. Comput Struct Biotechnol J 2015; 13:169-75. [PMID: 25848495 PMCID: PMC4372622 DOI: 10.1016/j.csbj.2015.02.004] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/16/2015] [Revised: 02/17/2015] [Accepted: 02/19/2015] [Indexed: 12/21/2022] Open
Abstract
Recent successes and failures of the application of 'enhanced' semiempirical QM (SQM) methods are reviewed in the light of the benefits and backdraws of adding dispersion (D) and hydrogen-bond (H) correction terms. We find that the accuracy of SQM-DH methods for non-covalent interactions is very often reported to be comparable to dispersion-corrected density functional theory (DFT-D), while computation times are about three orders of magnitude lower. SQM-DH methods thus open up a possibility to simulate realistically large model systems for problems both in life and materials science with comparably high accuracy.
Collapse
Affiliation(s)
| | - Martin Korth
- Institute of Theoretical Chemistry, Ulm University, D-89069 Ulm, Germany
| |
Collapse
|
13
|
Jiang T, Boereboom JM, Michel C, Fleurat-Lessard P, Bulo RE. Proton Transfer in Aqueous Solution: Exploring the Boundaries of Adaptive QM/MM. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2015. [DOI: 10.1007/978-3-319-21626-3_2] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|