Frozen density functional free energy simulations of redox proteins: computational studies of the reduction potential of plastocyanin and rusticyanin

The evaluation of reduction potentials of proteins by ab initio approaches presents a major challenge for computational chemistry. This is addressed in the present investigation by reporting detailed calculations of the reduction potentials of the blue copper proteins plastocyanin and rusticyanin using the QM/MM all-atom frozen density functional theory, FDFT, method. The relevant ab initio free energies are evaluated by using a classical reference potential. This approach appears to provide a general consistent and effective way for reproducing the configurational ensemble needed for consistent ab initio free energy calculations. The FDFT formulation allows us to treat a large part of the protein quantum mechanically by a consistently coupled QM/QM/MM embedding method while still retaining a proper configurational sampling. To establish the importance of proper configurational sampling and the need for a complete representation of the protein+solvent environment, we also consider several classical approaches. These include the semi-macroscopic PDLD/S-LRA method and classical all-atom simulations with and without a polarizable force field. The difference between the reduction potentials of the two blue copper proteins is reproduced in a reasonable way, and its origin is deduced from the different calculations. It is found that the protein permanent dipole tunes down the reduction potential for plastocyanin compared to the active site in regular water solvent, whereas in rusticyanin it is instead tuned up. This electrostatic environment, which is the major effect determining the reduction potential, is a property of the entire protein and solvent system and cannot be ascribed to any particular single interaction.     

The CHARMM–TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics,free energies,and excited state properties

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM–TURBOMOLE is well‐suited for extended QM/MM MD simulations using first principles methods with large (triple‐ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg²⁺(aq), where the MM outer sphere water molecules are represented using the SWM4‐NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg²⁺ and Zn²⁺ were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time‐averaged solution absorption spectrum of indole, the emission spectrum of the indole excited state, and the electronic circular dichroism spectrum of an oxacepham. © 2014 Wiley Periodicals, Inc.     

A polarizable continuum approach for the study of heterogeneous dielectric environments

We present a computational method, exploiting some features of the polarizable continuum model (PCM) to describe heterogeneous media; it belongs to the family of electrostatic embedding mixed methods, such as the more common quantum-mechanical (QM)/molecular mechanics approaches, with the electrostatic long range effects accounted for by a polarized continuum instead of atomic point charges. Provided effective dielectric constants are determined for the various parts of the system, the method is much faster than its atomistic counterpart, and allows for high-level QM calculations on the fragment of interest, using all the highly efficient computational tools developed for homogeneous PCM. Two case studies (the calculation of the pKa of solvent exposed acidic residues in a model protein, and the calculation of the electron spin resonance spectrum of a typical spin probe partially embedded in a membrane) are analyzed in some detail, to illustrate the application of the method to complex systems.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

Variational calculation of quantum mechanical/molecular mechanical free energy with electronic polarization of solvent

Quantum mechanical/molecular mechanical (QM/MM) free energy calculation presents a significant challenge due to an excessive number of QM calculations. A useful approach for reducing the computational cost is that based on the mean field approximation to the QM subsystem. Here, we describe such a mean-field QM/MM theory for electronically polarizable systems by starting from the Hartree product ansatz for the total system and invoking a variational principle of free energy. The MM part is then recast to a classical polarizable model by introducing the charge response kernel. Numerical test shows that the potential of mean force (PMF) thus obtained agrees quantitatively with that obtained from a direct QM/MM calculation, indicating the utility of self-consistent mean-field approximation. Next, we apply the obtained method to prototypical reactions in several qualitatively different solvents and make a systematic comparison of polarization effects. The results show that in aqueous solution the PMF does not depend very much on the water models employed, while in nonaqueous solutions the PMF is significantly affected by explicit polarization. For example, the free energy barrier for a phosphoryl dissociation reaction in acetone and cyclohexane is found to increase by more than 10 kcal/mol when switching the solvent model from an empirical to explicitly polarizable one. The reason for this is discussed based on the parametrization of empirical nonpolarizable models.     

Molecular dynamics simulations using the drude polarizable force field on GPUs with OpenMM: Implementation,validation, and benchmarks

Presented is the implementation of the Drude force field in the open‐source OpenMM simulation package allowing for access to graphical processing unit (GPU) hardware. In the Drude model, electronic degrees of freedom are represented by negatively charged particles attached to their parent atoms via harmonic springs, such that extra computational overhead comes from these additional particles and virtual sites representing lone pairs on electronegative atoms, as well as the associated thermostat and integration algorithms. This leads to an approximately fourfold increase in computational demand over additive force fields. However, by making the Drude model accessible to consumer‐grade desktop GPU hardware it will be possible to perform simulations of one microsecond or more in less than a month, indicating that the barrier to employ polarizable models has largely been removed such that polarizable simulations with the classical Drude model are readily accessible and practical.     

Polarization energies in oligoacene semiconductor crystals

Characterization of the electronically polarized environment and the nuclear relaxation that accompanies charge carriers is fundamental to charge transport in crystalline, polycrystalline, and amorphous organic solids. To study the polarization effects of localized charged carriers, we use quantum/classical QM/MM approaches with charge redistribution and polarizable force field schemes and apply them to crystals of naphthalene through pentacene. We describe the results of a comprehensive investigation of the electronic polarization energies in molecular crystal structures of these oligoacenes and discuss as well the evolution of the nuclear relaxation energies calculated for model oligoacene systems.     

Computational photochemistry of retinal proteins

High spectral tunability and quantum yield are the striking features of rhodopsin photochemistry. They rely on a strong and complex interaction of their chromophore, the protonated Schiff base of retinal, with its protein environment. In this article, we review the progress in the computational modeling of these systems, focusing on the optical properties and the excited state dynamics. While the earlier success of atomistic theoretical models was based on the breakthrough in X-ray crystallography and combined quantum mechanical molecular mechanical (QM/MM) methodology, recent advances point out the importance of high-level QM methods and the incorporation of effects that are neglected in conventional QM/MM or ONIOM schemes, like polarization and charge transfer.     

Theory and algorithms for the excited states of large molecules and molecular aggregates

This project aims to attack the frontiers of electronic structure calculations on the excited states of large molecules and molecular aggregates by developing novel theoretical and computational methods. The methodology development is especially based on the time-dependent density functional theory (TDDFT) and valence bond (VB) theory, and is expected to be computationally effective and accurate as well. Research works on the following related subjects will be performed: (1) The analytical energy-derivative approaches for electronically excited state within TDDFT will be developed to reduce bypass the computational costs in the calculation of molecular excited-state properties. (2) The ab initio methods for electronically excited state based on VB theory and hybrid TDDFT-VB method will be developed to overcome the limitations of current TDDFT in simulating photophysics and photochemistry. (3) For larger aggregates, neither ab initio methods nor TDDFT is applicable. We intend to build the effective model Hamiltonian by developing novel theoretical and computational methods to calculate the involved microscopic physical parameters from the first-principles methods. The constructed effective Hamiltonian is then used to describe the excitonic states and excitonic dynamics of the natural or artificial photosynthesized systems, organic or inorganic photovoltaic cell. (4) The condensed phase environment is taken into account by combining the developed theories and algorithms based on TDDFT and VB with the polarizable continuum solvent models (PCM), molecular mechanism (MM), classical electrodynamics (ED) or molecular dynamics (MD) theory. (5) Highly efficient software packages will be designed and developed.     

态定(线性响应)-极化连续模型/含时密度泛函方法研究一种有机发光材料的吸收和发射光谱

3-(二氰亚甲基)-5,5-二甲基-1-(3-[9-(2-乙基-己基)-咔唑基]-乙烯基)环己烷(DCDHCC)是一种用于光电器件中的有机染料,它具有良好的发光特性.我们使用含时密度泛函方法(TD-PBE0,TD-BMK和TD-M06)以及极化连续模型(PCM)计算了该材料在溶剂中的吸收和发射特性.计算中使用了线性响应(LR)、态定(SS)两种溶剂模型和6-31G(d)、6-31+G(d,p)两种基组.计算了DCDHCC在苯、四氢呋喃和丙酮溶剂中的吸收和发射光谱,并与实验观测进行了比较.结果表明:对于吸收光谱的计算,杂化函数的影响大于基组和溶剂模型,在三种函数中BMK更适于研究DCDHCC的吸收光谱;而对于发射光谱,基组的影响最大,基组通过影响激发态构型从而影响发射光谱,对于激发态构型的优化需要使用6-31+G(d,p)基组.我们希望这些研究能对今后设计类似的发光分子有帮助.     

Polarizable molecular mechanics studies of Cu(I)/Zn(II) superoxide dismutase: Bimetallic binding site and structured waters

The existence of a network of structured waters in the vicinity of the bimetallic site of Cu/Zn‐superoxide dismutase (SOD) has been inferred from high‐resolution X‐ray crystallography. Long‐duration molecular dynamics (MD) simulations could enable to quantify the lifetimes and possible interchanges of these waters between themselves as well as with a ligand diffusing toward the bimetallic site. The presence of several charged or polar ligands makes it necessary to resort to second‐generation polarizable potentials. As a first step toward such simulations, we benchmark in this article the accuracy of one such potential, sum of interactions between fragments Ab initio computed (SIBFA), by comparisons with quantum mechanics (QM) computations. We first consider the bimetallic binding site of a Cu/Zn‐SOD, in which three histidines and a water molecule are bound to Cu(I) and three histidines and one aspartate are bound to Zn(II). The comparisons are made for different His6 complexes with either one or both cations, and either with or without Asp and water. The total net charges vary from zero to three. We subsequently perform preliminary short‐duration MD simulations of 296 waters solvating Cu/Zn‐SOD. Six representative geometries are selected and energy‐minimized. Single‐point SIBFA and QM computations are then performed in parallel on model binding sites extracted from these six structures, each of which totals 301 atoms including the closest 28 waters from the Cu metal site. The ranking of their relative stabilities as given by SIBFA is identical to the QM one, and the relative energy differences by both approaches are fully consistent. In addition, the lowest‐energy structure, from SIBFA and QM, has a close overlap with the crystallographic one. The SIBFA calculations enable to quantify the impact of polarization and charge transfer in the ranking of the six structures. Five structural waters, which connect Arg141 and Glu131, are endowed with very high dipole moments (2.7–3.0 Debye), equal and larger than the one computed by SIBFA in ice‐like arrangements (2.7 D).     

Computational and data driven molecular material design assisted by low scaling quantum mechanics calculations and machine learning

Electronic structure methods based on quantum mechanics (QM) are widely employed in the computational predictions of the molecular properties and optoelectronic properties of molecular materials. The computational costs of these QM methods, ranging from density functional theory (DFT) or time-dependent DFT (TDDFT) to wave-function theory (WFT), usually increase sharply with the system size, causing the curse of dimensionality and hindering the QM calculations for large sized systems such as long polymer oligomers and complex molecular aggregates. In such cases, in recent years low scaling QM methods and machine learning (ML) techniques have been adopted to reduce the computational costs and thus assist computational and data driven molecular material design. In this review, we illustrated low scaling ground-state and excited-state QM approaches and their applications to long oligomers, self-assembled supramolecular complexes, stimuli-responsive materials, mechanically interlocked molecules, and excited state processes in molecular aggregates. Variable electrostatic parameters were also introduced in the modified force fields with the polarization model. On the basis of QM computational or experimental datasets, several ML algorithms, including explainable models, deep learning, and on-line learning methods, have been employed to predict the molecular energies, forces, electronic structure properties, and optical or electrical properties of materials. It can be conceived that low scaling algorithms with periodic boundary conditions are expected to be further applicable to functional materials, perhaps in combination with machine learning to fast predict the lattice energy, crystal structures, and spectroscopic properties of periodic functional materials.

Low scaling quantum mechanics calculations and machine learning can be employed to efficiently predict the molecular energies, forces, and optical and electrical properties of molecular materials and their aggregates.     

Coupled cluster theory with the polarizable continuum model of solvation

The environment may significantly affect molecular properties. Thus, it is desirable to account explicitly for these effects on the wave function and its derivatives, especially when the latter are evaluated with accurate methods, such as those belonging to coupled cluster (CC) theory. In this tutorial review, we discuss how to combine CC methods with the polarizable continuum model of solvation (PCM). We describe useful approximations that include the solvent response to the correlation and excited state equations while maintaining the computational cost comparable to in vacuo calculations. Although applied to PCM, the theoretical framework presented in this review is general and can be used with any polarizable embedding model. Representative applications of the CC-PCM method to ground and excited state properties of solvated molecules are presented, and comparisons with experiment, and between the full and approximate schemes are discussed.     

LICHEM: A QM/MM program for simulations with multipolar and polarizable force fields

We introduce an initial implementation of the LICHEM software package. LICHEM can interface with Gaussian, PSI4, NWChem, TINKER, and TINKER–HP to enable QM/MM calculations using multipolar/polarizable force fields. LICHEM extracts forces and energies from unmodified QM and MM software packages to perform geometry optimizations, single‐point energy calculations, or Monte Carlo simulations. When the QM and MM regions are connected by covalent bonds, the pseudo‐bond approach is employed to smoothly transition between the QM region and the polarizable force field. A series of water clusters and small peptides have been employed to test our initial implementation. The results obtained from these test systems show the capabilities of the new software and highlight the importance of including explicit polarization. © 2016 Wiley Periodicals, Inc.