Jens Antony

A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu

The method of dispersion correction as an add-on to standard Kohn-Sham density functional theory (DFT-D) has been refined regarding higher accuracy, broader range of applicability, and less empiricism. The main new ingredients are atom-pairwise specific dispersion coefficients and cutoff radii that are both computed from first principles. The coefficients for new eighth-order dispersion terms are computed using established recursion relations. System (geometry) dependent information is used for the first time in a DFT-D type approach by employing the new concept of fractional coordination numbers (CN). They are used to interpolate between dispersion coefficients of atoms in different chemical environments. The method only requires adjustment of two global parameters for each density functional, is asymptotically exact for a gas of weakly interacting neutral atoms, and easily allows the computation of atomic forces. Three-body nonadditivity terms are considered. The method has been assessed on standard benchmark sets for inter- and intramolecular noncovalent interactions with a particular emphasis on a consistent description of light and heavy element systems. The mean absolute deviations for the S22 benchmark set of noncovalent interactions for 11 standard density functionals decrease by 15%-40% compared to the previous (already accurate) DFT-D version. Spectacular improvements are found for a tripeptide-folding model and all tested metallic systems. The rectification of the long-range behavior and the use of more accurate C(6) coefficients also lead to a much better description of large (infinite) systems as shown for graphene sheets and the adsorption of benzene on an Ag(111) surface. For graphene it is found that the inclusion of three-body terms substantially (by about 10%) weakens the interlayer binding. We propose the revised DFT-D method as a general tool for the computation of the dispersion energy in molecules and solids of any kind with DFT and related (low-cost) electronic structure methods for large systems.

Density functional theory including dispersion corrections for intermolecular interactions in a large benchmark set of biologically relevant molecules.

Density functional theory including dispersion corrections (DFT-D) is applied to calculate intermolecular interaction energies in an extensive benchmark set consisting mainly of DNA base pairs and amino acid pairs, for which CCSD(T) complete basis set limit estimates are available (JSCH-2005 database). The three generalized gradient approximation (GGA) density functionals B-LYP, PBE and the new B97-D are tested together with the popular hybrid functional B3-LYP. The DFT-D interaction energies deviate on average by less than 1 kcal mol(-1) or 10% from the reference values. In only six out of 161 cases, the deviation exceeds 2 kcal mol(-1). With one exception, the few larger deviations occur for non-equilibrium structures extracted from experimental geometries. The largest absolute deviations are observed for pairs of oppositely charged amino acids which are, however, not significant on a relative basis due to the huge interaction energies > 100 kcal mol(-1) involved. The counterpoise (CP) correction for the basis set superposition error with the applied triple-zeta AO basis sets varies between 0 and -1 kcal mol(-1) (<5% of the interaction energy in most cases) except for four complexes, where it is up to -1.4 kcal mol(-1). It is thus suggested to skip the laborious calculation of the CP correction in DFT-D treatments with reasonable basis sets. The three dispersion corrected GGAs considered differ mainly for the interactions of the hydrogen-bonded DNA base pairs, which are systematically too small by 0.6 kcal mol(-1) in case of B97-D, while for PBE-D they are too high by 1.5 kcal mol(-1), and for B-LYP-D by 0.5 kcal mol(-1). The all in all excellent results that have been obtained at affordable computational costs suggest the DFT-D method to be a routine tool for many applications in organic chemistry or biochemistry.

Anharmonic midinfrared vibrational spectra of benzoic acid monomer and dimer

Anharmonic vibrational calculations for the benzoic acid monomer and dimer in the mid-IR regime (500-1800 cm(-1)) are reported. Harmonic frequencies and intensities are obtained at the DFT/B3LYP level of theory employing D95(d,p) and cc-pVTZ basis sets. Anharmonic corrections obtained from standard perturbation theory lead to redshifts of 1%-3%. In almost all cases, the resulting frequencies deviate by less than 1% from previous measurements [Bakker et al., J. Chem. Phys. 119, 11180 (2003)]. Calculated intensities are in qualitative agreement with the absorption experiment, with the cc-pVTZ values being superior to the D95(d,p) ones for a few modes of the dimer. The antisymmetric out-of-plane bending mode of the dimer, which is strongly blueshifted with respect to the monomer frequency, represents a remarkable exception: The harmonic frequencies obtained for the two basis sets differ notably from each other, and the anharmonically corrected frequencies deviate from the experimental value by 8% [D95(d,p)] or 3% (cc-pVTZ). Nonperturbative calculations in reduced dimensionality reveal that the relatively small total anharmonic shift (few tens of cm(-1)) comprises of partly much larger contributions (few hundreds of cm(-1)) which are mostly canceling each other. Many of the individual anharmonic couplings are beyond the validity of second-order perturbation theory based on cubic and semidiagonal quartic force constants only. This emphasizes the need for high-dimensional, nonperturbative anharmonic calculations at high quantum-chemical level when accurate frequencies of H-atom vibrations in double hydrogen bonds are sought for.

Blind Prediction of Binding Affinities for Charged Supramolecular Host–Guest Systems: Achievements and Shortcomings of DFT-D3

Association free energies ΔGa are calculated for two different types of host-guest systems, the rigid cucurbit[7]uril (CB7) and the basket shaped octa-acid (OA), and a number of charged guest molecules each by quantum chemical methods from first principles in the context of a recent blind test challenge (SAMPL4). For CB7, the overall agreement between theory and experiment is excellent. In comparison with all other submitted calculated relative ΔGa,rel values for this part of the blind test, our results ranked on top. Modeling the binding free energy in the case of the OA host mainly suffers from the problem that the binding situation is undefined with respect to the charge state and due to its intrinsic flexibility the host-guest complex is not represented well by a single configuration, but qualitative features of the binding process such as the proper binding orientation and the order of magnitude of ΔGa are represented in accord with the experimental expectations even though an accurate ranking is not possible.

Fully ab initio protein‐ligand interaction energies with dispersion corrected density functional theory

Dispersion corrected density functional theory (DFT‐D3) is used for fully ab initio protein‐ligand (PL) interaction energy calculation via molecular fractionation with conjugated caps (MFCC) and applied to PL complexes from the PDB comprising 3680, 1798, and 1060 atoms. Molecular fragments with n amino acids instead of one in the original MFCC approach are considered, thereby allowing for estimating the three‐body and higher many‐body terms. n > 1 is recommended both in terms of accuracy and efficiency of MFCC. For neutral protein side‐chains, the computed PL interaction energy is visibly independent of the fragment length n. The MFCC fractionation error is determined by comparison to a full‐system calculation for the 1060 atoms containing PL complex. For charged amino acid side‐chains, the variation of the MFCC result with n is increased. For these systems, using a continuum solvation model with a dielectricity constant typical for protein environments (ϵ = 4) reduces both the variation with n and improves the stability of the DFT calculations considerably. The PL interaction energies for two typical complexes obtained ab initio for the first time are found to be rather large (−30 and −54 kcal/mol).

Complexes of thiomandelate and captopril mercaptocarboxylate inhibitors to metallo‐β‐lactamase by polarizable molecular mechanics. Validation on model binding sites by quantum chemistry

Using the polarizable molecular mechanics method SIBFA, we have performed a search for the most stable binding modes of D‐ and L‐thiomandelate to a 104‐residue model of the metallo‐β‐lactamase from B. fragilis, an enzyme involved in the acquired resistance of bacteria to antibiotics. Energy balances taking into account solvation effects computed with a continuum reaction field procedure indicated the D‐isomer to be more stably bound than the L‐one, conform to the experimental result. The most stably bound complex has the S− ligand bridging monodentately the two Zn(II) cations and one carboxylate O− H‐bonded to the Asn193 side chain. We have validated the SIBFA energy results by performing additional SIBFA as well as quantum chemical (QC) calculations on small (88 atoms) model complexes extracted from the 104‐residue complexes, which include the residues involved in inhibitor binding. Computations were done in parallel using uncorrelated (HF) as well as correlated (DFT, LMP2, MP2) computations, and the comparisons extended to corresponding captopril complexes (Antony et al., J Comput Chem 2002, 23, 1281). The magnitudes of the SIBFA intermolecular interaction energies were found to correctly reproduce their QC counterparts and their trends for a total of twenty complexes.

Binding of D‐ and L‐captopril inhibitors to metallo‐β‐lactamase studied by polarizable molecular mechanics and quantum mechanics

The bacterial Zn2+ metallo‐β‐lactamase from B. fragilis is a zinc‐enzyme with two potential metal ion binding sites. It cleaves the lactam ring of antibiotics, thus contributing to the acquired resistance of bacteria against antibiotics. The present study bears on the binuclear form of the enzyme. We compare several possible binding modes of captopril, a mercaptocarboxamide inhibitor of several zinc‐metalloenzymes. Two diastereoisomers of captopril were considered, with either a D‐ or an L‐proline residue. We have used the polarizable molecular mechanics procedure SIBFA (Sum of Interactions Between Fragments ab initio computed). Two β‐lactamase models were considered, encompassing 104 and 188 residues, respectively. The energy balances included the inter and intramolecular interaction energies as well as the contribution from solvation computed using a continuum reaction field procedure. The thiolate ion of the inhibitor is binding to both metal ions, expelling the bridging solvent molecule from the uncomplexed enzyme. Different competing binding modes of captopril were considered, either where the inhibitor binds in a monodentate mode to the zinc cations only with its thiolate ion, or in bidentate modes involving additional zinc binding by its carboxylate or ketone carbonyl groups. The additional coordination by the inhibitors carboxylate or carbonyl group always occurs at the zinc ion, which is bound by a histidine, a cysteine, and an aspartate side chain. For both diastereomers, the energy balances favor monodentate binding of captopril via S−. The preference over bidentate binding is small. The interaction energies were recomputed in model sites restricted to captopril, the Zn2+ cations, and their coordinating end side chains from β‐lactamase (98 atoms). The interaction energies and their ranking among competing arrangements were consistent with those computed by ab initio HF and DFT procedures.

On the role of the axial ligand in heme-based catalysis of the peroxidase and P450 type

Abstract The present commentary focusses on the role of the axial ligand in peroxidase- and P450-type catalysis. Based on molecular orbital calculations and the experimental evidence available, it is argued that the ligand of a heme-containing enzyme may be a factor in setting the relative chance, although not the intrinsic capability, of the enzyme to catalyse a specific type of heme-based reaction chemistry. The ligand can do so by influencing the electrophilicity, i.e. the redox potential of the high-valency iron-oxo complex, and also by influencing the energy barrier for a reaction pathway through delocalization of valence electrons along the axial ligands, thereby, in the case of a cysteinate but not a histidine axial ligand, stabilizing oxygen transfer pathways.

First principle calculations of 113 Cd chemical shifts for proteins and model systems

Abstract113Cd isotropic NMR shieldings are calculated for a number of metal ion binding sites in proteins, using the GIAO-B3LYP and GIAO-HF methods with the uncontracted (19s15p9d4f) polarized basis set of Kellö and Sadlej on cadmium and 6-31G(d) on the ligands. The results compare favorably with experimental data, indicating that first principle calculations are a useful tool for structural interpretation of 113Cd chemical shift data from metal ion containing proteins. The effect of different ligand types (thiolate, imidazole, water, and monodentate carboxylate), coordination number, and deviations of the coordination geometry from ideal structures is evaluated. In particular, the ligand type and coordination number are important factors, but also changes in cadmium–ligand bond lengths may cause significant changes of the chemical shift.

Nonadiabatic effects on peptide vibrational dynamics induced by conformational changes

Quantum dynamical simulations of vibrational spectroscopy have been carried out for glycine dipeptide (CH(3)-CO-NH-CH(2)-CO-NH-CH(3)). Conformational structure and dynamics are modeled in terms of the two Ramachandran dihedral angles of the molecular backbone. Potential energy surfaces and harmonic frequencies are obtained from electronic structure calculations at the density functional theory (DFT) [B3LYP/6-31+G(d)] level. The ordering of the energetically most stable isomers (C(7) and C(5)) is reversed upon inclusion of the quantum mechanical zero point vibrational energy. Vibrational spectra of various isomers show distinct differences, mainly in the region of the amide modes, thereby relating conformational structures and vibrational spectra. Conformational dynamics is modeled by propagation of quantum mechanical wave packets. Assuming a directed energy transfer to the torsional degrees of freedom, transitions between the C(7) and C(5) minimum energy structures occur on a sub-picosecond time scale (700...800 fs). Vibrationally nonadiabatic effects are investigated for the case of the coupled, fundamentally excited amide I states. Using a two state-two mode model, the resulting wave packet dynamics is found to be strongly nonadiabatic due to the presence of a seam of the two potential energy surfaces. Initially prepared adiabatic vibrational states decay upon conformational change on a time scale of 200...500 fs with population transfer of more than 50% between the coupled amide I states. Also the vibrational energy transport between localized (excitonic) amide I vibrational states is strongly influenced by torsional dynamics of the molecular backbone where both enhanced and reduced decay rates are found. All these observations should allow the detection of conformational changes by means of time-dependent vibrational spectroscopy.