Energy decomposition analysis in solution based on the fragment molecular orbital method

We develop the pair interaction energy decomposition analysis (PIEDA) in solution by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The solvent screening of the electrostatic interaction and the desolvation penalty in complex formation are described by this approach from ab initio calculations of fragments and their pairs. The applications to the complex of solvated sodium and chlorine ions, as well as to lysine and aspartic acid, show how the analysis helps reveal the physical picture. The PIEDA/PCM method is also applied to a small protein chignolin (PDB: 1UAO), and the solvent screening of the pair interactions is discussed.     

Pair interaction energy decomposition analysis

The energy decomposition analysis (EDA) by Kitaura and Morokuma was redeveloped in the framework of the fragment molecular orbital method (FMO). The proposed pair interaction energy decomposition analysis (PIEDA) can treat large molecular clusters and the systems in which fragments are connected by covalent bonds, such as proteins. The interaction energy in PIEDA is divided into the same contributions as in EDA: the electrostatic, exchange-repulsion, and charge transfer energies, to which the correlation (dispersion) term was added. The careful comparison to the ab initio EDA interaction energies for water clusters with 2-16 molecules revealed that PIEDA has the error of at most 1.2 kcal/mol (or about 1%). The analysis was applied to (H2O)1024, the alpha helix, beta turn, and beta strand of polyalanine (ALA)10, as well as to the synthetic protein (PDB code 1L2Y) with 20 residues. The comparative aspects of the polypeptide isomer stability are discussed in detail.     

Energy gradients in combined fragment molecular orbital and polarizable continuum model (FMO/PCM) calculation

The analytic energy gradients for the combined fragment molecular orbital and polarizable continuum model (FMO/PCM) method are derived and implemented. Applications of FMO/PCM geometry optimization to polyalanine show that the structures obtained with the FMO/PCM method are very close to those obtained with the corresponding full ab initio PCM methods. FMO/PCM (RHF/6‐31G* level) is used to optimize the solution structure of the 304‐atom Trp‐cage miniprotein and the result is in agreement with NMR experiments. The key factors determining the relative stability of the α‐helix, β‐turn and the extended form in solution are elucidated for polyalanine. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Geometry optimizations of open-shell systems with the fragment molecular orbital method

The ability to perform geometry optimizations on large molecular systems is desirable for both closed- and open-shell species. In this work, the restricted open-shell Hartree-Fock (ROHF) gradients for the fragment molecular orbital (FMO) method are presented. The accuracy of the gradients is tested, and the ability of the method to reproduce adiabatic excitation energies is also investigated. Timing comparisons between the FMO method and full ab initio calculations are also performed, demonstrating the efficiency of the FMO method in modeling large open-shell systems.     

Analytic gradient for second order Møller-Plesset perturbation theory with the polarizable continuum model based on the fragment molecular orbital method

A new energy expression is proposed for the fragment molecular orbital method interfaced with the polarizable continuum model (FMO/PCM). The solvation free energy is shown to be more accurate on a set of representative polypeptides with neutral and charged residues, in comparison to the original formulation at the same level of the many-body expansion of the electrostatic potential determining the apparent surface charges. The analytic first derivative of the energy with respect to nuclear coordinates is formulated at the second-order M?ller-Plesset (MP2) perturbation theory level combined with PCM, for which we derived coupled perturbed Hartree-Fock equations. The accuracy of the analytic gradient is demonstrated on test calculations in comparison to numeric gradient. Geometry optimization of the small Trp-cage protein (PDB: 1L2Y) is performed with FMO/PCM/6-31(+)G(d) at the MP2 and restricted Hartree-Fock with empirical dispersion (RHF/D). The root mean square deviations between the FMO optimized and NMR experimental structure are found to be 0.414 and 0.426 A? for RHF/D and MP2, respectively. The details of the hydrogen bond network in the Trp-cage protein are revealed.     

Fully analytic energy gradient in the fragment molecular orbital method

The Z-vector equations are derived and implemented for solving the response term due to the external electrostatic potentials, and the corresponding contribution is added to the energy gradients in the framework of the fragment molecular orbital (FMO) method. To practically solve the equations for large molecules like proteins, the equations are decoupled by taking advantage of the local nature of fragments in the FMO method and establishing the self-consistent Z-vector method. The resulting gradients are compared with numerical gradients for the test molecular systems: (H(2)O)(64), alanine decamer, hydrated chignolin with the protein data bank (PDB) ID of 1UAO, and a Trp-cage miniprotein construct (PDB ID: 1L2Y). The computation time for calculating the response contribution is comparable to or less than that of the FMO self-consistent charge calculation. It is also shown that the energy gradients for the electrostatic dimer approximation are fully analytic, which significantly reduces the computational costs. The fully analytic FMO gradient is parallelized with an efficiency of about 98% on 32 nodes.     

Exploring chemistry with the fragment molecular orbital method

The fragment molecular orbital (FMO) method makes possible nearly linear scaling calculations of large molecular systems, such as water clusters, proteins and DNA. In particular, FMO has been widely used in biochemical applications involving protein-ligand binding and drug design. The method has been efficiently parallelized suitable for petascale computing. Many commonly used wave functions and solvent models have been interfaced with FMO. We review the historical background of FMO, and summarize its method development and applications.     

Rapid and accurate assessment of GPCR–ligand interactions Using the fragment molecular orbital‐based density‐functional tight‐binding method

The reliable and precise evaluation of receptor–ligand interactions and pair‐interaction energy is an essential element of rational drug design. While quantum mechanical (QM) methods have been a promising means by which to achieve this, traditional QM is not applicable for large biological systems due to its high computational cost. Here, the fragment molecular orbital (FMO) method has been used to accelerate QM calculations, and by combining FMO with the density‐functional tight‐binding (DFTB) method we are able to decrease computational cost 1000 times, achieving results in seconds, instead of hours. We have applied FMO‐DFTB to three different GPCR–ligand systems. Our results correlate well with site directed mutagenesis data and findings presented in the published literature, demonstrating that FMO‐DFTB is a rapid and accurate means of GPCR–ligand interactions. © 2017 Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Application of PCM thermal energy storage system to reduce building energy consumption

The building sector is known to make a large contribution to total energy consumption and CO₂ emissions. Phase change materials (PCMs) have been considered for thermal energy storage (TES) in buildings. They can balance out the discrepancies between energy demand and energy supply, which are temporally out of phase. However, traditional PCMs need special latent storage devices or containers to encapsulate the PCM, in order to store and release the latent heat of the PCM. The proper design of TES systems using a PCM requires quantitative information and knowledge about the heat transfer and phase change processes in the PCM. In Korea, radiant floor heating systems, which have traditionally been used in residential buildings, consume approximately 55% of the total residential building energy consumption in heating. This article reviews the development of available latent heat thermal energy storage technologies and discusses PCM application methods for residential building using radiant floor heating systems with the goal of reducing energy consumption.     

Copper(II) and cobalt(II) complexes of 2,6-diacetylpyridine bis(<Emphasis Type="Italic">O</Emphasis>-methyloxime): A theoretical investigation

The molecular geometries and vibrational frequencies of the title compounds in the ground state are calculated using the Hartree-Fock (HF) and density functional theory (DFT/B3LYP) methods with the LANL2DZ basis set and compared with the experimental data. The calculated results show that the optimized geometries can well reproduce the crystal structural parameters, and the theoretical vibrational frequencies show good agreement with the experimental values. The energetic behavior of the title compounds in solvent media is examined using the B3LYP method with the LANL2DZ basis set by applying the Onsager and polarizable continuum model (PCM). In addition, molecular electrostatic potential (MEP) and frontier molecular orbital (FMO) analyses of the title compounds are investigated by theoretical calculations.     

Achieving linear-scaling computational cost for the polarizable continuum model of solvation

This work describes a new and low-scaling implementation of the polarizable continuum model (PCM) for computing the self-consistent solvent reaction field. The PCM approach is both general and accurate. It is applicable in the framework of both quantum and classical calculations, and also to hybrid quantum/classical methods. In order to further extend the range of applicability of PCM we addressed the problem of its computational cost. The generation of the finite-elements molecular cavity has been reviewed and reimplemented, achieving linear scaling for systems containing up to 500 atoms. Linear scaling behavior has been achieved also for the iterative solution of the PCM equations, by exploiting the fast multipole method (FMM) for computing electrostatic interactions. Numerical results for large (both linear and globular) chemical systems are discussed.     

The extension of the fragment molecular orbital method with the many-particle Green's function

By using the many-particle Green's function (GF) the extension of the fragment molecular orbital (FMO) method by Kitaura et al. [Chem. Phys. Lett. 313, 701 (1999)] is proposed. It is shown that the partial summation of the cluster expansion of GF reproduces the same extrapolation formula as that of FMO. Therefore we can determine the excitation energy, the transition moment, and the linear response of a molecule from GF approximated with the FMO procedure. It is also shown that no wave function exists which is consistent to the FMO results. The perturbation expansion in which the self-consistent charge approximation defines the unperturbed state is reported. By using it the three-body effects missing in the pair approximation of FMO are analyzed and the corrections to the energy and the reduced density matrices are proposed. In contrast to the previous works these new corrections are not expressed as the addition or the subtraction of the energies of fragments. They are size extensive and require only the quantities available by the FMO calculation. The accuracy of these corrections is validated with the extended Hubbard model and the several test molecules.     

Covalent bond fragmentation suitable to describe solids in the fragment molecular orbital method

To improve the accuracy of the fragment molecular orbital method (FMO), we introduce a new fragmentation scheme based on using frozen orbitals to describe fractioned bonds. By applying this scheme to a set of polyalanine systems of up to 40 residues for the alpha-helix and beta-strand isomers, we established its accuracy, which is considerably improved compared to the original hybrid orbital projection method used for detaching bonds in FMO. For instance, at the two-body FMO expansion with the 6-311G* basis set, the error was typically reduced 2-4 times, and for 6-31G* the accuracy increase was even larger (10 times in terms of the maximum error). For the Trp-cage protein (PDB file 1L2Y) with many charged residues, a fairly large error was observed, which was shown to become small with a larger fragment size or at the three-body level. Consequently, we applied the new scheme to the adsorption of toluene and phenol on a faujasite zeolite, and we demonstrated that good accuracy can be achieved in reproducing ab initio results.     

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.     

Analytic gradient and molecular dynamics simulations using the fragment molecular orbital method combined with effective potentials

The completely analytic energy gradients are derived and implemented for the two-body fragment molecular orbital (FMO2) method combined with the model core potentials (MCP) and effective fragment potentials (EFP). The many-body terms in EFP require solving coupled-perturbed Hartree-Fock equations, which are derived and implemented. The molecular dynamics (MD) simulations are performed using the FMO2/MCP method for the capped alanine decamer and with the FMO2/EFP method for the zwitterionic conformer of glycine tetramer immersed in the water layer of 6.0 Å (135 water molecules). The results of the MD simulations using the FMO2/EFP and FMO2/MCP gradients show that the total energy is conserved at the time steps less than 1 fs.     

相变蓄热材料研究进展

相变蓄热材料(恒温潜热热能储存材料)是目前最热门的功能材料之一。它在发生相变时储存、放出的热量能够帮助所在系统进行能量的储存,同时可以一定程度上缓解双方在时间、强度及地点上的不匹配程度。相变蓄热材料优点突出,其中包括在使用过程中自身温度变化较小、有很好的稳定性、储热能力较强等。此类材料对环境友好,响应了国家近年来节能环保的政策,同时也可以极大地优化所在系统的运行效率。本文综述了近年来几类相变蓄热材料的种类、特点及国内外学者应对于不同缺陷做出的改进及其应用于行业的研究现状,并对未来的发展进行了探讨与展望。     

Numerical study of the effects of nanofluids and phase-change materials in photovoltaic thermal (PVT) systems

In this paper, the effects of pure water, SiO₂/water nanofluid, and a phase-change material (PCM) as coolants on the performance of a photovoltaic thermal (PVT) system are numerically investigated. The simulations are performed on two modules of PVT with PCM (PVT/PCM module) and without (PVT module). Parameters including PV surface temperature, thermal, and electrical efficiencies of the systems are studied and compared with each other. Moreover, the results of nanofluid as a working fluid is compared with those obtained using pure water. The results show that in the water-based PVT/PCM, the average PV cell temperature is decreased by 16 °C compared to that of the PVT system. This results in an increase of 8% in the electrical efficiency and 25% in the thermal efficiency. In addition, using nanofluid (SiO₂ with 1 and 3 mass% mass fraction) as a coolant in the PVT/PCM system increases the thermal efficiency by 3.51% and 10.40%, for 1 and 3 mass%, respectively, compared to that of the PVT/PCM with pure water as a coolant. This study shows that increasing the melting temperature of the phase-change material leads to an increase in the thermal efficiency of the PVT/PCM system.

     

Assessment and acceleration of binding energy calculations for protein–ligand complexes by the fragment molecular orbital method

In the field of drug discovery, it is important to accurately predict the binding affinities between target proteins and drug applicant molecules. Many of the computational methods available for evaluating binding affinities have adopted molecular mechanics‐based force fields, although they cannot fully describe protein–ligand interactions. A noteworthy computational method in development involves large‐scale electronic structure calculations. Fragment molecular orbital (FMO) method, which is one of such large‐scale calculation techniques, is applied in this study for calculating the binding energies between proteins and ligands. By testing the effects of specific FMO calculation conditions (including fragmentation size, basis sets, electron correlation, exchange‐correlation functionals, and solvation effects) on the binding energies of the FK506‐binding protein and 10 ligand complex molecule, we have found that the standard FMO calculation condition, FMO2‐MP2/6‐31G(d), is suitable for evaluating the protein–ligand interactions. The correlation coefficient between the binding energies calculated with this FMO calculation condition and experimental values is determined to be R = 0.77. Based on these results, we also propose a practical scheme for predicting binding affinities by combining the FMO method with the quantitative structure–activity relationship (QSAR) model. The results of this combined method can be directly compared with experimental binding affinities. The FMO and QSAR combined scheme shows a higher correlation with experimental data (R = 0.91). Furthermore, we propose an acceleration scheme for the binding energy calculations using a multilayer FMO method focusing on the protein–ligand interaction distance. Our acceleration scheme, which uses FMO2‐HF/STO‐3G:MP2/6‐31G(d) at R_int = 7.0 Å, reduces computational costs, while maintaining accuracy in the evaluation of binding energy. © 2015 Wiley Periodicals, Inc.     

Unrestricted Hartree-Fock based on the fragment molecular orbital method: Energy and its analytic gradient

A consideration of the surrounding environment is necessary for a meaningful analysis of the reaction activity in large molecular systems. We propose an approach to perform unrestricted Hartree-Fock (UHF) calculations within the framework of the fragment molecular orbital (FMO) method (FMO-UHF) to study large systems with unpaired electrons. Prior to an energy analysis one has to optimize geometry, which requires an accurate analytic energy gradient. We derive the FMO-UHF energy and its analytic gradient and implement them into GAMESS. The performance of FMO-UHF is evaluated for a solvated organic molecule and a solvated metal complex, as well as for the active part of a protein, in terms of energy, gradient, and geometry optimization.