Improving the efficiency and convergence of geometry optimization with the polarizable continuum model: new energy gradients and molecular surface tessellation

New equations are derived and implemented for efficient and accurate computation of solvation energy derivatives for the conductor-like polarizable continuum model (C-PCM) and the isotropic integral equation formalism polarizable continuum model (IEF-PCM). Two new molecular surface tessellation procedures GEPOL-RT and GEPOL-AS that generate near continuous potential energy surfaces are proposed for PCM geometry optimization. The combined use of these new techniques leads to efficient and convergent geometry optimizations with the PCMs.     

Separation of isomeric disaccharides by traveling wave ion mobility mass spectrometry using CO2 as drift gas

The use of CO₂ as a massive and polarizable drift gas is shown to greatly improve peak‐to‐peak resolution (R_p‐p), as compared with N₂, for the separation of disaccharides in a Synapt G2 traveling wave ion mobility cell. Near or baseline R_p‐p was achieved for three pairs of sodiated molecules of disaccharide isomers, that is, cellobiose and sucrose (R_p‐p = 0.76), maltose and sucrose (R_p‐p = 1.04), and maltose and lactose (R_p‐p = 0.74). Ion mobility mass spectrometry using CO₂ as the drift gas offers therefore an attractive alternative for fast and efficient separation of isomeric disaccharides. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparative assessment of computational methods for the determination of solvation free energies in alcohol‐based molecules

The determination of differences in solvation free energies between related drug molecules remains an important challenge in computational drug optimization, when fast and accurate calculation of differences in binding free energy are required. In this study, we have evaluated the performance of five commonly used polarized continuum model (PCM) methodologies in the determination of solvation free energies for 53 typical alcohol and alkane small molecules. In addition, the performance of these PCM methods, of a thermodynamic integration (TI) protocol and of the Poisson–Boltzmann (PB) and generalized Born (GB) methods, were tested in the determination of solvation free energies changes for 28 common alkane‐alcohol transformations, by the substitution of an hydrogen atom for a hydroxyl substituent. The results show that the solvation model D (SMD) performs better among the PCM‐based approaches in estimating solvation free energies for alcohol molecules, and solvation free energy changes for alkane‐alcohol transformations, with an average error below 1 kcal/mol for both quantities. However, for the determination of solvation free energy changes on alkane‐alcohol transformation, PB and TI yielded better results. TI was particularly accurate in the treatment of hydroxyl groups additions to aromatic rings (0.53 kcal/mol), a common transformation when optimizing drug‐binding in computer‐aided drug design. © 2013 Wiley Periodicals, Inc.     

Valence state parameters of all transition metal atoms in metalloproteins—development of ABEEMσπ fluctuating charge force field

To promote accuracy of the atom‐bond electronegativity equalization method (ABEEMσπ) fluctuating charge polarizable force fields, and extend it to include all transition metal atoms, a new parameter, the reference charge is set up in the expression of the total energy potential function. We select over 700 model molecules most of which model metalloprotein molecules that come from Protein Data Bank. We set reference charges for different apparent valence states of transition metals and calibrate the parameters of reference charges, valence state electronegativities, and valence state hardnesses for ABEEMσπ through linear regression and least square method. These parameters can be used to calculate charge distributions of metalloproteins containing transition metal atoms (Sc‐Zn, Y‐Cd, and Lu‐Hg). Compared the results of ABEEMσπ charge distributions with those obtained by ab initio method, the quite good linear correlations of the two kinds of charge distributions are shown. The reason why the STO‐3G basis set in Mulliken population analysis for the parameter calibration is specially explained in detail. Furthermore, ABEEMσπ method can also quickly and quite accurately calculate dipole moments of molecules. Molecular dynamics optimizations of five metalloproteins as the examples show that their structures obtained by ABEEMσπ fluctuating charge polarizable force field are very close to the structures optimized by the ab initio MP2/6–311G method. This means that the ABEEMσπ/MM can now be applied to molecular dynamics simulations of systems that contain metalloproteins with good accuracy. © 2014 Wiley Periodicals, Inc.     

Modeling of Solute-Solvent Interactions Using an External Electric Field—From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute–solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.     

1,3,5-Triaza-7-Phosphaadamantane (PTA) as a 31P NMR Probe for Organometallic Transition Metal Complexes in Solution

Due to the rigid structure of 1,3,5-triaza-7-phosphaadamantane (PTA), its ³¹P chemical shift solely depends on non-covalent interactions in which the molecule is involved. The maximum range of change caused by the most common of these, hydrogen bonding, is only 6 ppm, because the active site is one of the PTA nitrogen atoms. In contrast, when the PTA phosphorus atom is coordinated to a metal, the range of change exceeds 100 ppm. This feature can be used to support or reject specific structural models of organometallic transition metal complexes in solution by comparing the experimental and Density Functional Theory (DFT) calculated values of this ³¹P chemical shift. This approach has been tested on a variety of the metals of groups 8–12 and molecular structures. General recommendations for appropriate basis sets are reported.     

Three‐body expansion of the fragment molecular orbital method combined with density‐functional tight‐binding

The three‐body fragment molecular orbital (FMO3) method is formulated for density‐functional tight‐binding (DFTB). The energy, analytic gradient, and Hessian are derived in the gas phase, and the energy and analytic gradient are also derived for polarizable continuum model. The accuracy of FMO3‐DFTB is evaluated for five proteins, sodium cation in explicit solvent, and three isomers of polyalanine. It is shown that FMO3‐DFTB is considerably more accurate than FMO2‐DFTB. Molecular dynamics simulations for sodium cation in water are performed for 100 ps, yielding radial distribution functions and coordination numbers. © 2017 Wiley Periodicals, Inc.     

Substituent and solvent effects on the structural bioactivity and anticancer characteristic of catechin as a bioactive constituent of green tea

In this study, quantum chemical calculations using density functional theory and 6‐311G (d) basis set have been applied to analyze the substituent effect on the electronic structural properties including thermochemical parameters as well as anticancer characteristics of catechin as a bioactive constituent of green tea. It has been found out that different substituents with two kinds of electron donating and electron withdrawing groups including ? NO₂, ? NH₂, ? Cl, ? OCH₃, and ? CH₃ leaded to different structural stabilities. The catechin with chloro substituents has been referred to the most stable catechin derivatives based on the minimum energy of investigated compounds. In this study, the calculated values have been compared in gas phase and different solvent media with a wide range of solvent dielectrics using the self‐consistent reaction field method with polarizable continuum model. Relative interaction energies between solute and solvent as well as dipole moment values have been analyzed. The results of solvent effect study revealed that the formation of hydrogen bonds accounts for the stabilization of catechin. These calculated data were found to be a logical way in predicting the relative stability of catechin derivatives and can be useful for drug design of anticancer drugs and would open a new door for the researchers who are interested in natural products. The results presented in this article will be helpful to improve existing model and will open up a whole new arena of study to understand the anticancer characteristic of these antioxidant drugs derived from natural products against the photosensitizing DNA damage induce by catechin extracts and enable a much clearer understanding of how anticancer drugs mediate their effects on the specific spot. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011