Simulations of Raman Spectra Using the Fragment Molecular Orbital Method

Combining fragmentation method and high-performance computing: Geometry optimization and vibrational spectra of proteins

Exploring the structures and spectral features of proteins with advanced quantum chemical methods is an uphill task. In this work, a fragment-based molecular tailoring approach (MTA) is appraised for the CAM-B3LYP/aug-cc-pVDZ-level geometry optimization and vibrational infrared (IR) spectra calculation of ten real proteins containing up to 407 atoms and 6617 basis functions. The use of MTA and the inherently parallel nature of the fragment calculations enables a rapid and accurate calculation of the IR spectrum. The applicability of MTA to optimize the protein geometry and evaluate its IR spectrum employing a polarizable continuum model with water as a solvent is also showcased. The typical errors in the total energy and IR frequencies computed by MTA vis-à-vis their full calculation (FC) counterparts for the studied protein are 5-10 millihartrees and 5 cm-1, respectively. Moreover, due to the independent execution of the fragments, large-scale parallelization can also be achieved. With increasing size and level of theory, MTA shows an appreciable advantage in computer time as well as memory and disk space requirement over the corresponding FCs. The present study suggests that the geometry optimization and IR computations on the biomolecules containing ∼1000 atoms and/or ∼15 000 basis functions using MTA and HPC facility can be clearly envisioned in the near future.

Analytic Gradient for Time-Dependent Density Functional Theory Combined with the Fragment Molecular Orbital Method

The analytic energy gradient of energy with respect to nuclear coordinates is derived for the fragment molecular orbital (FMO) method combined with time-dependent density functional theory (TDDFT). The response terms arising from the use of a polarizable embedding are derived. The obtained analytic FMO-TDDFT gradient is shown to be accurate in comparison to both numerical FMO-TDDFT and unfragmented TDDFT gradients, at the level of two- and three-body expansions. The gradients are used for geometry optimizations, molecular dynamics, vibrational calculations, and simulations of IR and Raman spectra of excited states. The developed method is used to optimize the geometry of the ground and excited electronic states of the photoactive yellow protein (PDB: 2PHY).

Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method

Fast parameterized methods such as density-functional tight-binding (DFTB) facilitate realistic calculations of large molecular systems, which can be accelerated by the fragment molecular orbital (FMO) method. Fragmentation facilitates interaction analyses between functional parts of molecular systems. In addition to DFTB, other parameterized methods combined with FMO are also described. Applications of FMO methods to biochemical and inorganic systems are reviewed.

Partitioning of the Vibrational Free Energy

Vibrational energies are partitioned into the contributions of molecular parts called segments, for instance, residues in proteins. The fragment molecular orbital method is used to facilitate vibrational calculations of large systems at the DFTB and HF-3c levels. The vibrational analysis is combined with the partitioning of the electronic energy, yielding free-energy contributions of segments to the binding energy, pinpointing hot spots for drug discovery and other studies. The analysis is illustrated on two protein-ligand complexes in solution.

Computational Methods for Biochemical Simulations Implemented in GAMESS

Computational methods for modeling biochemical processes implemented in GAMESS package are reviewed; in particular, quantum mechanics combined with molecular mechanics (QM/MM), semi-empirical, and fragmentation approaches. A detailed summary of capabilities is provided for the QM/MM implementation in QuanPol program and the fragment molecular orbital (FMO) method. Molecular modeling and visualization packages useful for biochemical simulations with GAMESS are described. GAMESS capabilities with corresponding references are tabulated for reader's convenience.

Geometry Optimization, Transition State Search, and Reaction Path Mapping Accomplished with the Fragment Molecular Orbital Method

Recent development of the fragment molecular orbital (FMO) method related to energy gradients, geometry optimization, transition state search, and chemical reaction mapping is summarized. The frozen domain formulation of FMO is introduced in detail, and the structure of related GAMESS input files for FMO is described.

Fantasy versus reality in fragment-based quantum chemistry

Since the introduction of the fragment molecular orbital method 20 years ago, fragment-based approaches have occupied a small but growing niche in quantum chemistry. These methods decompose a large molecular system into subsystems small enough to be amenable to electronic structure calculations, following which the subsystem information is reassembled in order to approximate an otherwise intractable supersystem calculation. Fragmentation sidesteps the steep rise (with respect to system size) in the cost of ab initio calculations, replacing it with a distributed cost across numerous computer processors. Such methods are attractive, in part, because they are easily parallelizable and therefore readily amenable to exascale computing. As such, there has been hope that distributed computing might offer the proverbial "free lunch" in quantum chemistry, with the entrée being high-level calculations on very large systems. While fragment-based quantum chemistry can count many success stories, there also exists a seedy underbelly of rarely acknowledged problems. As these methods begin to mature, it is time to have a serious conversation about what they can and cannot be expected to accomplish in the near future. Both successes and challenges are highlighted in this Perspective.

Simulations of infrared and Raman spectra in solution using the fragment molecular orbital method

Calculation of IR and Raman spectra in solution for large molecular systems made possible with analytic FMO/PCM Hessians.

Quantum mechanical force fields for condensed phase molecular simulations

Molecular simulations are powerful tools for providing atomic-level details into complex chemical and physical processes that occur in the condensed phase. For strongly interacting systems where quantum many-body effects are known to play an important role, density-functional methods are often used to provide the model with the potential energy used to drive dynamics. These methods, however, suffer from two major drawbacks. First, they are often too computationally intensive to practically apply to large systems over long time scales, limiting their scope of application. Second, there remain challenges for these models to obtain the necessary level of accuracy for weak non-bonded interactions to obtain quantitative accuracy for a wide range of condensed phase properties. Quantum mechanical force fields (QMFFs) provide a potential solution to both of these limitations. In this review, we address recent advances in the development of QMFFs for condensed phase simulations. In particular, we examine the development of QMFF models using both approximate and ab initio density-functional models, the treatment of short-ranged non-bonded and long-ranged electrostatic interactions, and stability issues in molecular dynamics calculations. Example calculations are provided for crystalline systems, liquid water, and ionic liquids. We conclude with a perspective for emerging challenges and future research directions.

Vibrational infrared and Raman spectra of polypeptides: Fragments-in-fragments within molecular tailoring approach

The present work reports the calculation of vibrational infrared (IR) and Raman spectra of large molecular systems employing molecular tailoring approach (MTA). Further, it extends the grafting procedure for the accurate evaluation of IR and Raman spectra of large molecular systems, employing a new methodology termed as Fragments-in-Fragments (FIF), within MTA. Unlike the previous MTA-based studies, the accurate estimation of the requisite molecular properties is achieved without performing any full calculations (FC). The basic idea of the grafting procedure is implemented by invoking the nearly basis-set-independent nature of the MTA-based error vis-à-vis the respective FCs. FIF has been tested out for the estimation of the above molecular properties for three isomers, viz., β-strand, 310- and α-helix of acetyl(alanine)nNH2 (n = 10, 15) polypeptides, three conformers of doubly protonated gramicidin S decapeptide and trpzip2 protein (PDB id: 1LE1), respectively, employing BP86/TZVP, M06/6-311G**, and M05-2X/6-31G** levels of theory. For most of the cases, a maximum difference of 3 cm(-1) is achieved between the grafted-MTA frequencies and the corresponding FC values. Further, a comparison of the BP86/TZVP level IR and Raman spectra of α-helical (alanine)20 and its N-deuterated derivative shows an excellent agreement with the existing experimental spectra. In view of the requirement of only MTA-based calculations and the ability of FIF to work at any level of theory, the current methodology provides a cost-effective solution for obtaining accurate spectra of large molecular systems.

Accurate vibrational spectra via molecular tailoring approach: a case study of water clusters at MP2 level

In spite of the recent advents in parallel algorithms and computer hardware, high-level calculation of vibrational spectra of large molecules is still an uphill task. To overcome this, significant effort has been devoted to the development of new algorithms based on fragmentation methods. The present work provides the details of an efficient and accurate procedure for computing the vibrational spectra of large clusters employing molecular tailoring approach (MTA). The errors in the Hessian matrix elements and dipole derivatives arising due to the approximation nature of MTA are reduced by grafting the corrections from a smaller basis set. The algorithm has been tested out for obtaining vibrational spectra of neutral and charged water clusters at Møller-Plesset second order level of theory, and benchmarking them against the respective full calculation (FC) and/or experimental results. For (H2O)16 clusters, the estimated vibrational frequencies are found to differ by a maximum of 2 cm(-1) with reference to the corresponding FC values. Unlike the FC, the MTA-based calculations including grafting procedure can be performed on a limited hardware, yet take a fraction of the FC time. The present methodology, thus, opens a possibility of the accurate estimation of the vibrational spectra of large molecular systems, which is otherwise impossible or formidable.

Fragment quantum chemical approach to geometry optimization and vibrational spectrum calculation of proteins

Geometry optimization and vibrational spectra (infrared and Raman spectra) calculations of proteins are carried out by a quantum chemical approach using the EE-GMFCC (electrostatically embedded generalized molecular fractionation with conjugate caps) method (J. Phys. Chem. A, 2013, 117, 7149).

Elucidating the fundamental forces in protein crystal formation: the case of crambin

This study demonstrates the feasibility of periodic all-electron hybrid density functional theory calculations in the description of protein crystals, using crambin as a test case.

Molecular simulations of proteins have been usually accomplished through empirical or semi-empirical potentials, due to the large size and inherent complexity of these biological systems. On the other hand, a theoretical description of proteins based on quantum-mechanical methods would however provide an unbiased characterization of their electronic properties, possibly offering a link between these and the ultimate biological activity. Yet, such approaches have been historically hindered by the large amount of requested computational power. Here we demonstrate the feasibility of periodic all-electron density functional theory calculations in the description of the crystal of the protein crambin (46 aminoacids), which is determined with exceptional structural accuracy. We have employed the hybrid B3LYP functional, coupled to an empirical description of London interactions (D*) to simulate the crambin crystal with an increasing amount of lattice water molecules in the cell (up to 172H₂O per cell). The agreement with the experiment is good for both protein geometry and protein–water interactions. The energetics was computed to predict crystal formation energies, protein–water and protein–protein interaction energies. We studied the role of dispersion interactions which are crucial for holding the crambin crystal in place. B3LYP-D* electrostatic potential and dipole moment of crambin as well as the electronic charge flow from crambin to the solvating water molecules (0.0015e per H₂O) have also been predicted. These results proved that quantum-mechanical simulations of small proteins, both free and in their crystalline state, are now feasible in a reasonable amount of time, by programs capable of exploiting high performance computing architectures, allowing the study of protein properties not easily amenable through classical force fields.

The use of many-body expansions and geometry optimizations in fragment-based methods

Conspectus Chemists routinely work with complex molecular systems: solutions, biochemical molecules, and amorphous and composite materials provide some typical examples. The questions one often asks are what are the driving forces for a chemical phenomenon? How reasonable are our views of chemical systems in terms of subunits, such as functional groups and individual molecules? How can one quantify the difference in physicochemical properties of functional units found in a different chemical environment? Are various effects on functional units in molecular systems additive? Can they be represented by pairwise potentials? Are there effects that cannot be represented in a simple picture of pairwise interactions? How can we obtain quantitative values for these effects? Many of these questions can be formulated in the language of many-body effects. They quantify the properties of subunits (fragments), referred to as one-body properties, pairwise interactions (two-body properties), couplings of two-body interactions described by three-body properties, and so on. By introducing the notion of fragments in the framework of quantum chemistry, one obtains two immense benefits: (a) chemists can finally relate to quantum chemistry, which now speaks their language, by discussing chemically interesting subunits and their interactions and (b) calculations become much faster due to a reduced computational scaling. For instance, the somewhat academic sounding question of the importance of three-body effects in water clusters is actually another way of asking how two hydrogen bonds affect each other, when they involve three water molecules. One aspect of this is the many-body charge transfer (CT), because the charge transfers in the two hydrogen bonds are coupled to each other (not independent). In this work, we provide a generalized view on the use of many-body expansions in fragment-based methods, focusing on the general aspects of the property expansion and a contraction of a many-body expansion in a formally two-body series, as exemplified in the development of the fragment molecular orbital (FMO) method. Fragment-based methods have been very successful in delivering the properties of fragments, as well as the fragment interactions, providing insights into complex chemical processes in large molecular systems. We briefly review geometry optimizations performed with fragment-based methods and present an efficient geometry optimization method based on the combination of FMO with molecular mechanics (MM), applied to the complex of a subunit of protein kinase 2 (CK2) with a ligand. FMO results are discussed in comparison with experimental and MM-optimized structures.