Efficient determination of accurate atomic polarizabilities for polarizeable embedding calculations

We evaluate embedding potentials, obtained via various methods, used for polarizable embedding computations of excitation energies of para‐nitroaniline in water and organic solvents as well as of the green fluorescent protein. We found that isotropic polarizabilities derived from DFTD3 dispersion coefficients correlate well with those obtained via the LoProp method. We show that these polarizabilities in conjunction with appropriately derived point charges are in good agreement with calculations employing static multipole moments up to quadrupoles and anisotropic polarizabilities for both computed systems. The (partial) use of these easily‐accessible parameters drastically reduces the computational effort to obtain accurate embedding potentials especially for proteins. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

The large quadrupole of water molecules

Many quantum mechanical calculations indicate water molecules in the gas and liquid phase have much larger quadrupole moments than any of the common site models of water for computer simulations. Here, comparisons of multipoles from quantum mechanical∕molecular mechanical (QM∕MM) calculations at the MP2∕aug-cc-pVQZ level on a B3LYP∕aug-cc-pVQZ level geometry of a waterlike cluster and from various site models show that the increased square planar quadrupole can be attributed to the p-orbital character perpendicular to the molecular plane of the highest occupied molecular orbital as well as a slight shift of negative charge toward the hydrogens. The common site models do not account for the p-orbital type electron density and fitting partial charges of TIP4P- or TIP5P-type models to the QM∕MM dipole and quadrupole give unreasonable higher moments. Furthermore, six partial charge sites are necessary to account reasonably for the large quadrupole, and polarizable site models will not remedy the problem unless they account for the p-orbital in the gas phase since the QM calculations show it is present there too. On the other hand, multipole models by definition can use the correct multipoles and the electrostatic potential from the QM∕MM multipoles is much closer than that from the site models to the potential from the QM∕MM electron density. Finally, Monte Carlo simulations show that increasing the quadrupole in the soft-sticky dipole-quadrupole-octupole multipole model gives radial distribution functions that are in good agreement with experiment.     

Accuracy of distributed multipoles and polarizabilities: comparison between the LoProp and MpProp models

Localized multipole moments up to the fifth moment as well as localized dipole polarizabilities are calculated with the MpProp and the newly developed LoProp methods for a total of 20 molecules, predominantly derived from amino acids. A comparison of electrostatic potentials calculated from the multipole expansion obtained by the two methods with ab initio results shows that both methods reproduce the electrostatic interaction with an elementary charge with a mean absolute error of approximately 1.5 kJ/mol at contact distance and less than 0.1 kJ/mol at distances 2 A further out when terms up to the octupole moments are included. The polarizabilities are tested with homogenous electric fields and are found to have similar accuracy. The MpProp method gives better multipole moments unless diffuse basis sets are used, whereas LoProp gives better polarizabilities.     

The coupling constant polarizability and hyperpolarizabilty of 1J(NH) in N-methylacetamide, and its application for the multipole spin-spin coupling constant polarizability/reaction field approach to solvation

We present a benchmark study of a combined multipole spin-spin coupling constant (SSCC) polarizability/reaction field (MJP/RF) approach to the calculation of both specific and bulk solvation effects on SSCCs of solvated molecules. The MJP/RF scheme is defined by an expansion of the SSCCs of the solvated molecule in terms of coupling constant dipole and quadrupole polarizabilities and hyperpolarizabilities derived from single molecule ab initio calculations. The solvent electric field and electric field gradient are calculated based on data derived from molecular dynamics (MD) simulations thereby accounting for solute-solvent dynamical effects. The MJP/RF method is benchmarked against polarizable QM/MM calculations for the one-bond N-H coupling constant in N-methylacetamide. The best agreement between the MJP/RF and QM/MM approaches is found by truncating the electric field expansion in the MJP/RF approach at the linear electric field level. In addition, we investigate the sensitivity of the results due to the choice of one-electron basis set in the ab initio calculations of the coupling constant (hyper-)polarizabilities and find that they are affected by the basis set in a way similar to the coupling constants themselves.     

Fragmentation-based QM/MM simulations: length dependence of chain dynamics and hydrogen bonding of polyethylene oxide and polyethylene in aqueous solutions

Molecular fragmentation quantum mechanics (QM) calculations have been combined with molecular mechanics (MM) to construct the fragmentation QM/MM method for simulations of dilute solutions of macromolecules. We adopt the electrostatics embedding QM/MM model, where the low-cost generalized energy-based fragmentation calculations are employed for the QM part. Conformation energy calculations, geometry optimizations, and Born-Oppenheimer molecular dynamics simulations of poly(ethylene oxide), PEO(n) (n = 6-20), and polyethylene, PE(n) ( n = 9-30), in aqueous solution have been performed within the framework of both fragmentation and conventional QM/MM methods. The intermolecular hydrogen bonding and chain configurations obtained from the fragmentation QM/MM simulations are consistent with the conventional QM/MM method. The length dependence of chain conformations and dynamics of PEO and PE oligomers in aqueous solutions is also investigated through the fragmentation QM/MM molecular dynamics simulations.     

Conformational simulations of aqueous solvated alpha-conotoxin GI and its single disulfide analogues using a polarizable force field model

The solution conformation of alpha-conotoxin GI and its two single disulfide analogues are simulated using a polarizable force field in combination with the molecular fragmentation quantum chemical calculation. The polarizability is explicitly described by allowing the partial charges and fragment dipole moments to be variables, with values coming from the linear-scaling energy-based molecular fragmentation calculations at the B3LYP/6-31G(d) level. In comparison with the full quantum chemical calculations, the fragmentation approaches can yield precise ground-state energies, dipole moments, and static polarizabilities for peptides. The B3LYP/6-31G(d) charges and fragment-centered dipole moments are introduced in calculations of electrostatic terms in both AmberFF03 and OPLS force fields. Our test calculations on the gas-phase glucagon (PDB code: 1gcn) and solvated alpha-conotoxin GI (PDB code: 1not) demonstrate that the present polarization model is capable of describing the structural properties (such as the relative conformational energies, intramolecular hydrogen bonds, and disulfide bonds) with accuracy comparable to some other polarizable force fields (ABEEM/MM and OPLS-PFF) and the quantum mechanics/molecular mechanics (QM/MM) hybrid model. The employment of fragment-centered dipole moments in calculations of dipole-dipole interactions can save computational time in comparison with those polarization models using atom-centered dipole moments without much loss of accuracy. The molecular dynamics simulations using the polarizable force field demonstrate that two single disulfide GI analogues are more flexible and less structured than the native alpha-conotoxin GI, in agreement with NMR experiments. The polarization effect is important in simulations of the folding/unfolding process of solvated proteins.     

Density functional self-consistent quantum mechanics/molecular mechanics theory for linear and nonlinear molecular properties: Applications to solvated water and formaldehyde

A combined quantum mechanics/molecular mechanics (QM/MM) method is described, where the polarization between the solvent and solute is accounted for using a self-consistent scheme linear in the solvent polarization. The QM/MM method is implemented for calculation of energies and molecular response properties including the calculation of linear and quadratic response functions using the density-functional theory (DFT) and the Hartree-Fock (HF) theory. Sample calculations presented for ground-state energies, first-order ground-state properties, excitation energies, first-order excited state properties, polarizabilities, first-hyperpolarizabilities, and two-photon absorptions strengths of formaldehyde suggests that DFT may in some cases be a sufficiently reliable alternative to high-level theory, such as coupled-cluster (CC) theory, in modeling solvent shifts, whereas results obtained with the HF wave function deviate significantly from the CC results. Calculations carried out on water gives results that also are comparable with CC calculations in accuracy for ground-state and first-order properties. However, to obtain such accuracy an exchange-correlation functional capable of describing the diffuse Rydberg states must be chosen.     

Contributions of pauli repulsions to the energetics and physical properties computed in QM/MM methods

Conventional combined quantum mechanical/molecular mechanical (QM/MM) methods lack explicit treatment of Pauli repulsions between the quantum‐mechanical and molecular‐mechanical subsystems. Instead, classical Lennard‐Jones (LJ) potentials between QM and MM nuclei are used to model electronic Pauli repulsion and long‐range London dispersion, despite the fact that the latter two are inherently of quantum nature. Use of the simple LJ potential in QM/MM methods can reproduce minimal geometries and energies of many molecular clusters reasonably well, as compared to full QM calculations. However, we show here that the LJ potential cannot correctly describe subtle details of the electron density of the QM subsystem because of the neglect of Pauli repulsions between the QM and MM subsystems. The inaccurate electron density subsequently affects the calculation of electronic and magnetic properties of the QM subsystem. To explicitly consider Pauli interactions with QM/MM methods, we propose a method to use empirical effective potentials on the MM atoms. The test case of the binding energy and magnetic properties of a water dimer shows promising results for the general application of effective potentials to mimic Pauli repulsions in QM/MM calculations. © 2013 Wiley Periodicals, Inc.     

Converging ligand‐binding free energies obtained with free‐energy perturbations at the quantum mechanical level

In this article, the convergence of quantum mechanical (QM) free‐energy simulations based on molecular dynamics simulations at the molecular mechanics (MM) level has been investigated. We have estimated relative free energies for the binding of nine cyclic carboxylate ligands to the octa‐acid deep‐cavity host, including the host, the ligand, and all water molecules within 4.5 Å of the ligand in the QM calculations (158–224 atoms). We use single‐step exponential averaging (ssEA) and the non‐Boltzmann Bennett acceptance ratio (NBB) methods to estimate QM/MM free energy with the semi‐empirical PM6‐DH2X method, both based on interaction energies. We show that ssEA with cumulant expansion gives a better convergence and uses half as many QM calculations as NBB, although the two methods give consistent results. With 720,000 QM calculations per transformation, QM/MM free‐energy estimates with a precision of 1 kJ/mol can be obtained for all eight relative energies with ssEA, showing that this approach can be used to calculate converged QM/MM binding free energies for realistic systems and large QM partitions. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Computational linear dependence in molecular electronic structure calculations using universal basis sets

Distributed universal even‐tempered basis sets have been developed over recent years that are capable of supporting Hartree–Fock energies to an accuracy approaching the sub‐μHartree level. These basis sets have also been exploited in correlation studies, in applications to polyatomic molecules, and in the calculation of electric properties, such as multipole moments, polarizabilities, and hyperpolarizabilities. Jorge and coworkers have also developed universal basis sets and have recently reported applications to diatomic molecular systems. In this article, we compare the molecular calculations reported by Jorge and coworkers with our previous studies. Particular attention is given to the degree of computational linear dependence associated with the various basis sets employed and the consequential effects of the accuracy of the calculated energies. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005