Water interactions with hydrophobic groups: assessment and recalibration of semiempirical molecular orbital methods

Hydration effects on the vibrational properties of carboxylates: From continuum models to QM/MM simulations

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.

A specific MNDO parameterization for water

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.

Semi-Empirical Born-Oppenheimer Molecular Dynamics (SEBOMD) within the Amber Biomolecular Package

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.

Thermally induced hopping model for long-range triplet excitation energy transfer in DNA

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.

Amber-Compatible Parametrization Procedure for Peptide-like Compounds: Application to 1,4- and 1,5-Substituted Triazole-Based Peptidomimetics

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.

Why many semiempirical molecular orbital theories fail for liquid water and how to fix them

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.

Enhanced semiempirical QM methods for biomolecular interactions

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.