1
|
Sarangi R, Maity S, Acharya A. Machine Learning Approach to Vertical Energy Gap in Redox Processes. J Chem Theory Comput 2024; 20:6747-6755. [PMID: 39044422 PMCID: PMC11325558 DOI: 10.1021/acs.jctc.4c00715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
A straightforward approach to calculating the free energy change (ΔG) and reorganization energy of a redox process is linear response approximation (LRA). However, accurate prediction of redox properties is still challenging due to difficulties in conformational sampling and vertical energy-gap sampling. Expensive hybrid quantum mechanical/molecular mechanical (QM/MM) calculations are typically employed in sampling energy gaps using conformations from simulations. To alleviate the computational cost associated with the expensive QM method in the QM/MM calculation, we propose machine learning (ML) methods to predict the vertical energy gaps (VEGs). We tested several ML models to predict the VEGs and observed that simple models like linear regression show excellent performance (mean absolute error ∼0.1 eV) in predicting VEGs in all test systems, even when using features extracted from cheaper semiempirical methods. Our best ML model (extra trees regressor) shows a mean absolute error of around 0.1 eV while using features from the cheapest QM method. We anticipate our approach can be generalized to larger macromolecular systems with more complex redox centers.
Collapse
Affiliation(s)
- Ronit Sarangi
- Department of Chemistry, Syracuse University, Syracuse, New York 13244, United States
| | - Suman Maity
- Department of Chemistry, Syracuse University, Syracuse, New York 13244, United States
| | - Atanu Acharya
- Department of Chemistry, Syracuse University, Syracuse, New York 13244, United States
- BioInspired Syracuse, Syracuse University, Syracuse, New York 13244, United States
| |
Collapse
|
2
|
Choudhury A, Santra S, Ghosh D. Understanding the Photoprocesses in Biological Systems: Need for Accurate Multireference Treatment. J Chem Theory Comput 2024; 20:4951-4964. [PMID: 38864715 DOI: 10.1021/acs.jctc.4c00027] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Light-matter interaction is crucial to life itself and revolves around many of the central processes in biology. The need for understanding these photochemical and photophysical processes cannot be overemphasized. Interaction of light with biological systems starts with the absorption of light and subsequent phenomena that occur in the excited states of the system. However, excited states are typically difficult to understand within the mean field approximation of quantum chemical methods. Therefore, suitable multireference methods and methodologies have been developed to understand these phenomena. In this Perspective, we will describe a few methods and methodologies suitable for these descriptions and discuss some persisting difficulties.
Collapse
Affiliation(s)
- Arpan Choudhury
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Supriyo Santra
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Debashree Ghosh
- School of Chemical Sciences, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| |
Collapse
|
3
|
Spencer RJ, Zhanserkeev AA, Yang EL, Steele RP. The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings. J Am Chem Soc 2024; 146:15376-15392. [PMID: 38771156 DOI: 10.1021/jacs.4c03198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.
Collapse
Affiliation(s)
- Ryan J Spencer
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Asylbek A Zhanserkeev
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Emily L Yang
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Ryan P Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| |
Collapse
|
4
|
Chung MKJ, Ponder JW. Beyond isotropic repulsion: Classical anisotropic repulsion by inclusion of p orbitals. J Chem Phys 2024; 160:174118. [PMID: 38748037 PMCID: PMC11078554 DOI: 10.1063/5.0203678] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 04/15/2024] [Indexed: 05/19/2024] Open
Abstract
Accurate modeling of intermolecular repulsion is an integral component in force field development. Although repulsion can be explicitly calculated by applying the Pauli exclusion principle, this approach is computationally viable only for systems of limited sizes. Instead, it has previously been shown that repulsion can be reformulated in a "classical" picture: the Pauli exclusion principle prohibits electrons from occupying the same state, leading to a depletion of electronic charge between atoms, giving rise to an enhanced nuclear-nuclear electrostatic repulsion. This classical picture is called the isotropic S2/R approximation, where S is the overlap and R is the interatomic distance. This approximation accurately captures the repulsion of isotropic atoms such as noble gas dimers; however, a key deficiency is that it fails to capture the angular dependence of the repulsion of anisotropic molecules. To include directionality, the wave function must at least be a linear combination of s and p orbitals. We derive a new anisotropic S2/R repulsion model through the inclusion of the anisotropic p orbital term in the total wave function. Because repulsion is pairwise and decays rapidly, it can be truncated at a short range, making it amenable for efficient calculation of energy and forces in complex biomolecular systems. We present a parameterization of the S101 dimer database against the ab initio benchmark symmetry-adapted perturbation theory, which yields an rms error of only 0.9 kcal/mol. The importance of the anisotropic term is demonstrated through angular scans of water-water dimers and dimers involving halobenzene. Simulation of liquid water shows that the model can be computed efficiently for realistic system sizes.
Collapse
Affiliation(s)
| | - Jay W. Ponder
- Author to whom correspondence should be addressed: . Tel.: 314-935-4275
| |
Collapse
|
5
|
Slipchenko LV. Detangling Solvatochromic Effects by the Effective Fragment Potential Method. J Phys Chem A 2024; 128:656-669. [PMID: 38193780 DOI: 10.1021/acs.jpca.3c06194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.
Collapse
Affiliation(s)
- Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| |
Collapse
|
6
|
Kim S, Conrad JA, Tow GM, Maginn EJ, Boatz JA, Gordon MS. Intermolecular interactions in clusters of ethylammonium nitrate and 1-amino-1,2,3-triazole. Phys Chem Chem Phys 2023; 25:30428-30457. [PMID: 37917371 DOI: 10.1039/d3cp02407e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
The intermolecular interaction energies, including hydrogen bonds (H-bonds), of clusters of the ionic liquid ethylammonium nitrate (EAN) and 1-amino-1,2,3-triazole (1-AT) based deep eutectic propellants (DeEP) are examined. 1-AT is introduced as a neutral hydrogen bond donor (HBD) to EAN in order to form a eutectic mixture. The effective fragment potential (EFP) is used to examine the bonding interactions in the DeEP clusters. The resolution of the Identity (RI) approximated second order Møller-Plesset perturbation theory (RI-MP2) and coupled cluster theory (RI-CCSD(T)) are used to validate the EFP results. The EFP method predicts that there are significant polarization and charge transfer effects in the EAN:1-AT complexes, along with Coulombic, dispersion and exchange repulsion interactions. The EFP interaction energies are in good agreement with the RI-MP2 and RI-CCSD(T) results. The quasi-atomic orbital (QUAO) bonding and kinetic bond order (KBO) analyses are additionally used to develop a conceptual and semi-quantitative understanding of the H-bonding interactions as a function of the size of the system. The QUAO and KBO analyses suggest that the H-bonds in the examined clusters follow the characteristic hydrogen bonding three-center four electron interactions. The strongest H-bonding interactions between the (EAN)1:(1-AT)n and (EAN)2:(1-AT)n (n = 1-5) complexes are observed internally within EAN; that is, between the ethylammonium cation [EA]+ and the nitrate anion ([NO3]-). The weakest H-bonding interactions occur between [NO3]- and 1-AT. Consequently, the average strengths of the H-bonds within a given (EAN)x:(1-AT)n complex decrease as more 1-AT molecules are introduced into the EAN monomer and EAN dimer. The QUAO bonding analysis suggests that 1-AT in (EAN)x:(1-AT)n can act as both a HBD and a hydrogen bond acceptor simultaneously. It is observed that two 1-AT molecules can form H-bonds to each other. Although the KBOs that correspond to H-bonding interactions in [EA]+:1-AT, [NO3]-:1-AT and between two 1-AT molecules are weaker than the H-bonds in EAN, those weak H-bond networks with 1-AT could be important to form a stable DeEP.
Collapse
Affiliation(s)
- Shinae Kim
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, USA
| | - Justin A Conrad
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
| | - Garrett M Tow
- Department of Chemical & Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Edward J Maginn
- Department of Chemical & Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Jerry A Boatz
- Aerospace Systems Directorate, Air Force Research Laboratory, Edwards Air Force Base, California 93524, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, IA 50011, USA.
| |
Collapse
|
7
|
Mashkovtsev D, Orimoto Y, Aoki Y. Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation. J Chem Theory Comput 2023; 19:5548-5562. [PMID: 37471461 DOI: 10.1021/acs.jctc.3c00266] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.
Collapse
Affiliation(s)
- Denis Mashkovtsev
- Department of Molecular and Material Sciences, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| |
Collapse
|
8
|
Kim Y, Mitchell Z, Lawrence J, Morozov D, Savikhin S, Slipchenko LV. Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example. J Phys Chem Lett 2023; 14:7038-7044. [PMID: 37524046 DOI: 10.1021/acs.jpclett.3c01461] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Zach Mitchell
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Jack Lawrence
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
9
|
Church JR, Olsen JMH, Schapiro I. Induction effects on the absorption maxima of photoreceptor proteins. Biophys Physicobiol 2023; 20:e201007. [PMID: 38362325 PMCID: PMC10865876 DOI: 10.2142/biophysico.bppb-v20.s007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2022] [Accepted: 01/20/2023] [Indexed: 01/24/2023] Open
Abstract
Multiscale simulations have been established as a powerful tool to calculate and predict excitation energies in complex systems such as photoreceptor proteins. In these simulations the chromophore is typically treated using quantum mechanical (QM) methods while the protein and surrounding environment are described by a classical molecular mechanics (MM) force field. The electrostatic interactions between these regions are often treated using electrostatic embedding where the point charges in the MM region polarize the QM region. A more sophisticated treatment accounts also for the polarization of the MM region. In this work, the effect of such a polarizable embedding on excitation energies was benchmarked and compared to electrostatic embedding. This was done for two different proteins, the lipid membrane-embedded jumping spider rhodopsin and the soluble cyanobacteriochrome Slr1393g3. It was found that the polarizable embedding scheme produces absorption maxima closer to experimental values. The polarizable embedding scheme was also benchmarked against expanded QM regions and found to be in qualitative agreement. Treating individual residues as polarizable recovered between 50% and 71% of the QM improvement in the excitation energies, depending on the system. A detailed analysis of each amino acid residue in the chromophore binding pocket revealed that aromatic residues result in the largest change in excitation energy compared to the electrostatic embedding. Furthermore, the computational efficiency of polarizable embedding allowed it to go beyond the binding pocket and describe a larger portion of the environment, further improving the results.
Collapse
Affiliation(s)
- Jonathan R. Church
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | | | - Igor Schapiro
- Fritz Haber Center for Molecular Dynamics Research, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| |
Collapse
|
10
|
Makoś MZ, Gurunathan PK, Raugei S, Kowalski K, Glezakou VA, Rousseau R. Modeling Absolute Redox Potentials of Ferrocene in the Condensed Phase. J Phys Chem Lett 2022; 13:10005-10010. [PMID: 36264148 DOI: 10.1021/acs.jpclett.2c02447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
Absolute thermodynamic quantities for critical chemical reactions are needed to determine the role of solvents and reactive environments in catalysis and electrocatalysis. Theoretical methods can provide such quantification but are often hindered by the innate complexity of electron correlation and dynamic relaxation of solvent environments. We present and validate a protocol for calculating the redox potentials of the ferrocene/ferrocenium redox pair in acetonitrile. Equation-of-motion and effective fragment potential (EFP) methods are used to characterize the adiabatic and vertical ionization potentials as well as the electron affinity processes. We benchmark molecular mechanics against the EFP model to show the differences in the ferrocene electronic polarizability in two redox states. Our best estimate of the redox potential (4.94 eV) agrees well with the experimental value (4.93 eV). This demonstrates the ability of modern computational methods to predict absolute redox potentials quantitatively and to quantify the correlation of dynamic effects, which underlie their origin.
Collapse
Affiliation(s)
- Małgorzata Zofia Makoś
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Pradeep Kumar Gurunathan
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Simone Raugei
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | | | - Roger Rousseau
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
11
|
Nakamura T, Yokaichiya T, Fedorov DG. Analysis of Guest Adsorption on Crystal Surfaces Based on the Fragment Molecular Orbital Method. J Phys Chem A 2022; 126:957-969. [PMID: 35080391 DOI: 10.1021/acs.jpca.1c10229] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
For gaining insights into interactions in periodic systems, an analysis is developed based on the fragment molecular orbital method combined with periodic boundary conditions. The adsorption energy is decomposed into guest and surface polarization and deformation energy, guest-surface and guest-guest interactions, and the vibrational free energy. The analysis is applied to the adsorption of guest molecules to Ih (001) ice surface. The cooperativity effects result in a non-linear change in the adsorption energy with coverage due to many-body effects. The role of dispersion is found to be dominant for guests with long hydrophobic tails. A rule is proposed relating the length of the alkyl tail with the formation of the guest layer. The computed binding enthalpies are in good agreement with experimental values. For high coverage, adsorbed molecules can form an ordered layer known as self-assembled monolayer.
Collapse
Affiliation(s)
- Taiji Nakamura
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Tomoko Yokaichiya
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
12
|
Chen CG, Nardi AN, Amadei A, D’Abramo M. Theoretical Modeling of Redox Potentials of Biomolecules. Molecules 2022; 27:1077. [PMID: 35164342 PMCID: PMC8838479 DOI: 10.3390/molecules27031077] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 01/21/2022] [Accepted: 01/25/2022] [Indexed: 11/28/2022] Open
Abstract
The estimation of the redox potentials of biologically relevant systems by means of theoretical-computational approaches still represents a challenge. In fact, the size of these systems typically does not allow a full quantum-mechanical treatment needed to describe electron loss/gain in such a complex environment, where the redox process takes place. Therefore, a number of different theoretical strategies have been developed so far to make the calculation of the redox free energy feasible with current computational resources. In this review, we provide a survey of such theoretical-computational approaches used in this context, highlighting their physical principles and discussing their advantages and limitations. Several examples of these approaches applied to the estimation of the redox potentials of both proteins and nucleic acids are described and critically discussed. Finally, general considerations on the most promising strategies are reported.
Collapse
Affiliation(s)
- Cheng Giuseppe Chen
- Department of Chemistry, Sapienza University of Rome, 00185 Rome, Italy; (C.G.C.); (A.N.N.)
| | | | - Andrea Amadei
- Department of Chemical and Technological Sciences, Tor Vergata University, 00133 Rome, Italy;
| | - Marco D’Abramo
- Department of Chemistry, Sapienza University of Rome, 00185 Rome, Italy; (C.G.C.); (A.N.N.)
| |
Collapse
|
13
|
Wang X, Li X, He X, Zhang JZH. A fixed multi-site interaction charge model for an accurate prediction of the QM/MM interactions. Phys Chem Chem Phys 2021; 23:21001-21012. [PMID: 34522933 DOI: 10.1039/d1cp02776j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A fixed multi-site interaction charge (FMIC) model was proposed for the accurate prediction of intermolecular electrostatic interactions based on the quantum mechanical linear response of a molecule to an external electric field. In such a model, some additional off-center interaction sites were added for capturing multipole interactions for a given molecule. By multivariate least-square fitting analysis of the calculated QM/MM interactions of a given molecule with the electrostatic environment and the electrostatic potentials of the environment at the pre-defined distributed interaction sites, the FMIC of the molecule was obtained. The model system of CO in myoglobin (Mb) was utilized to demonstrate the derivation of the FMIC. The accuracy of FMIC in predicting the electrostatic interactions between CO and the Mb environment was investigated using 10 000 different Mb-CO configurations generated from the 400 ps QM/MM MD simulation. In comparison to the QM/MM calculations at the B3LYP/aug-cc-pVTZ/ff99SB level, the mean unsigned error (MUE) of the results based on the FMIC model was merely 0.10 kcal mol-1, and the root mean square error (RMSE) was only 0.13 kcal mol-1, which are significantly lower than the results predicted by the ESP charge model (MUE = 1.45 kcal mol-1, and RMSE = 1.7 kcal mol-1, respectively). The transferability of FMIC was tested by applying the obtained FMIC in the wild type Mb-CO system to the mutants of V68F and H64L Mb-CO systems. The MUEs of the obtained results for 10 000 different configurations are both smaller than 0.2 kcal mol-1 for the V68F and H64L Mb-CO systems in comparison to the B3LYP/aug-cc-pVTZ/ff99SB calculations, and the RMSEs are also lower than 0.2 kcal mol-1 for both mutants. The applications of FMIC were extended to model the electrostatic interactions between a hydrogen fluoride molecule and 492 waters in a truncated octahedron box; our study showed that the FMIC could give satisfactory results with a MUE of 0.12 kcal mol-1 and a RMSE of 0.16 kcal mol-1 in comparison to the B3LYP/aug-cc-pVDZ/TIP3P calculations for 10 000 different configurations generated using the 10 ns classical MD simulation. Therefore, the FMIC method provides an accurate and efficient tool for predicting intermolecular electrostatic interactions, which can be utilized in the future development of molecular force fields.
Collapse
Affiliation(s)
- Xianwei Wang
- College of Science, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China. .,Shanghai Engineering Research Center for Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry and Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.
| | - Xilong Li
- College of Computer Science and Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310023, China
| | - Xiao He
- Shanghai Engineering Research Center for Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry and Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China. .,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.
| | - John Z H Zhang
- Shanghai Engineering Research Center for Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry and Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China. .,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China. .,Department of Chemistry, New York University, New York, New York 10003, USA
| |
Collapse
|
14
|
Nakamura T, Yokaichiya T, Fedorov DG. Quantum-Mechanical Structure Optimization of Protein Crystals and Analysis of Interactions in Periodic Systems. J Phys Chem Lett 2021; 12:8757-8762. [PMID: 34478310 DOI: 10.1021/acs.jpclett.1c02510] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A fast quantum-mechanical approach, density-functional tight-binding combined with the fragment molecular orbital method and periodic boundary conditions, is used to optimize atomic coordinates and cell parameters for a set of protein crystals: 1ETL, 5OQZ, 3Q8J, 1CBN, and 2VB1. Good agreement between experimental and calculated structures is obtained for both atomic coordinates and cell parameters. Sterical clashes present in the experimental structures are corrected by simulations. The partition analysis is extended to treat periodic boundary conditions and applied to analyze protein-solvent interactions in crystals.
Collapse
Affiliation(s)
- Taiji Nakamura
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Tomoko Yokaichiya
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
15
|
Karnaukh EA, Bravaya KB. The redox potential of a heme cofactor in Nitrosomonas europaea cytochrome c peroxidase: a polarizable QM/MM study. Phys Chem Chem Phys 2021; 23:16506-16515. [PMID: 34017969 PMCID: PMC11178132 DOI: 10.1039/d0cp06632j] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
Redox reactions are crucial to biological processes that protect organisms against oxidative stress. Metalloenzymes, such as peroxidases which reduce excess reactive oxygen species into water, play a key role in detoxification mechanisms. Here we present the results of a polarizable QM/MM study of the reduction potential of the electron transfer heme in the cytochrome c peroxidase of Nitrosomonas europaea. We have found that environment polarization does not substantially affect the computed value of the redox potential. Particular attention has been given to analyzing the role of electrostatic interactions within the protein environment and the solvent on tuning the redox potential of the heme co-factor. We have found that the electrostatic interactions predominantly explain the fluctuations of the vertical ionization/attachment energies of the heme for the sampled configurations, and that the long range electrostatic interactions (up to 40 Å) contribute substantially to the absolute values of the vertical energy gaps.
Collapse
|
16
|
Wang Z, Liu W. iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem. J Chem Theory Comput 2021; 17:4831-4845. [PMID: 34240856 DOI: 10.1021/acs.jctc.1c00445] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.
Collapse
Affiliation(s)
- Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| |
Collapse
|
17
|
Abstract
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.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
18
|
Nochebuena J, Naseem-Khan S, Cisneros GA. Development and application of quantum mechanics/molecular mechanics methods with advanced polarizable potentials. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1515. [PMID: 34367343 PMCID: PMC8341087 DOI: 10.1002/wcms.1515] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Accepted: 12/19/2020] [Indexed: 01/02/2023]
Abstract
Quantum mechanics/molecular mechanics (QM/MM) simulations are a popular approach to study various features of large systems. A common application of QM/MM calculations is in the investigation of reaction mechanisms in condensed-phase and biological systems. The combination of QM and MM methods to represent a system gives rise to several challenges that need to be addressed. The increase in computational speed has allowed the expanded use of more complicated and accurate methods for both QM and MM simulations. Here, we review some approaches that address several common challenges encountered in QM/MM simulations with advanced polarizable potentials, from methods to account for boundary across covalent bonds and long-range effects, to polarization and advanced embedding potentials.
Collapse
Affiliation(s)
- Jorge Nochebuena
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - Sehr Naseem-Khan
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| | - G Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas, USA
| |
Collapse
|
19
|
Glick ZL, Koutsoukas A, Cheney DL, Sherrill CD. Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles. J Chem Phys 2021; 154:224103. [PMID: 34241239 DOI: 10.1063/5.0050444] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole-multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule's electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.
Collapse
Affiliation(s)
- Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Alexios Koutsoukas
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - Daniel L Cheney
- Molecular Structure and Design, Bristol Myers Squibb Company, P.O. Box 5400, Princeton, New Jersey 08543, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
20
|
Nottoli M, Cupellini L, Lipparini F, Granucci G, Mennucci B. Multiscale Models for Light-Driven Processes. Annu Rev Phys Chem 2021; 72:489-513. [PMID: 33561359 DOI: 10.1146/annurev-physchem-090419-104031] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.
Collapse
Affiliation(s)
- Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Giovanni Granucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy;
| |
Collapse
|
21
|
Błasiak B, Bednarska JD, Chołuj M, Góra RW, Bartkowiak W. Ab initio effective one-electron potential operators: Applications for charge-transfer energy in effective fragment potentials. J Comput Chem 2021; 42:398-411. [PMID: 33349929 DOI: 10.1002/jcc.26462] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 11/18/2020] [Indexed: 12/19/2022]
Abstract
The concept of effective one-electron potentials (EOPs) has proven to be extremely useful in efficient description of electronic structure of chemical systems, especially extended molecular aggregates such as interacting molecules in condensed phases. Here, a general method for EOP-based elimination of electron repulsion integrals is presented, that is tuned toward the fragment-based calculation methodologies such as the second generation of the effective fragment potentials (EFP2) method. Two general types of the EOP operator matrix elements are distinguished and treated either via the distributed multipole expansion or the extended density fitting (DF) schemes developed in this work. The EOP technique is then applied to reduce the high computational costs of the effective fragment charge-transfer (CT) terms being the bottleneck of EFP2 potentials. The alternative EOP-based CT energy model is proposed, derived within the framework of intermolecular perturbation theory with Hartree-Fock noninteracting reference wavefunctions, compatible with the original EFP2 formulation. It is found that the computational cost of the EFP2 total interaction energy calculation can be reduced by up to 38 times when using the EOP-based formulation of CT energy, as compared to the original EFP2 scheme, without compromising the accuracy for a wide range of weakly interacting neutral and ionic molecular fragments. The proposed model can thus be used routinely within the EFP2 framework.
Collapse
Affiliation(s)
- Bartosz Błasiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Joanna D Bednarska
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Marta Chołuj
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Robert W Góra
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| | - Wojciech Bartkowiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland
| |
Collapse
|
22
|
Mroginski MA, Adam S, Amoyal GS, Barnoy A, Bondar AN, Borin VA, Church JR, Domratcheva T, Ensing B, Fanelli F, Ferré N, Filiba O, Pedraza-González L, González R, González-Espinoza CE, Kar RK, Kemmler L, Kim SS, Kongsted J, Krylov AI, Lahav Y, Lazaratos M, NasserEddin Q, Navizet I, Nemukhin A, Olivucci M, Olsen JMH, Pérez de Alba Ortíz A, Pieri E, Rao AG, Rhee YM, Ricardi N, Sen S, Solov'yov IA, De Vico L, Wesolowski TA, Wiebeler C, Yang X, Schapiro I. Frontiers in Multiscale Modeling of Photoreceptor Proteins. Photochem Photobiol 2021; 97:243-269. [PMID: 33369749 DOI: 10.1111/php.13372] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 11/01/2020] [Indexed: 02/06/2023]
Abstract
This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.
Collapse
Affiliation(s)
| | - Suliman Adam
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Gil S Amoyal
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Avishai Barnoy
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ana-Nicoleta Bondar
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Veniamin A Borin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Jonathan R Church
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Tatiana Domratcheva
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Department Biomolecular Mechanisms, Max Planck Institute for Medical Research, Heidelberg, Germany
| | - Bernd Ensing
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Francesca Fanelli
- Department of Life Sciences, Center for Neuroscience and Neurotechnology, Università degli Studi di Modena e Reggio Emilia, Modena, Italy
| | | | - Ofer Filiba
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Laura Pedraza-González
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | - Ronald González
- Institut für Chemie, Technische Universität Berlin, Berlin, Germany
| | | | - Rajiv K Kar
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Lukas Kemmler
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Seung Soo Kim
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Odense, Denmark
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| | - Yigal Lahav
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel.,MIGAL - Galilee Research Institute, S. Industrial Zone, Kiryat Shmona, Israel
| | - Michalis Lazaratos
- Freie Universität Berlin, Department of Physics, Theoretical Molecular Biophysics Group, Berlin, Germany
| | - Qays NasserEddin
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Isabelle Navizet
- MSME, Univ Gustave Eiffel, CNRS UMR 8208, Univ Paris Est Creteil, Marne-la-Vallée, France
| | - Alexander Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia.,Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow, Russia
| | - Massimo Olivucci
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy.,Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark.,Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, Tromsø, Norway
| | - Alberto Pérez de Alba Ortíz
- Van 't Hoff Institute for Molecular Science and Amsterdam Center for Multiscale Modeling, University of Amsterdam, Amsterdam, The Netherlands
| | - Elisa Pieri
- Aix-Marseille Univ, CNRS, ICR, Marseille, France
| | - Aditya G Rao
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology, Daejeon, Korea
| | - Niccolò Ricardi
- Département de Chimie Physique, Université de Genève, Genève, Switzerland
| | - Saumik Sen
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Ilia A Solov'yov
- Department of Physics, Carl von Ossietzky University Oldenburg, Oldenburg, Germany
| | - Luca De Vico
- Department of Biotechnology, Chemistry and Pharmacy, Università degli Studi di Siena, Siena, Italy
| | | | - Christian Wiebeler
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| | - Xuchun Yang
- Chemistry Department, Bowling Green State University, Bowling Green, OH, USA
| | - Igor Schapiro
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel
| |
Collapse
|
23
|
Abstract
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.
Collapse
|
24
|
Osella S. Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead? NANOMATERIALS 2021; 11:nano11020299. [PMID: 33498961 PMCID: PMC7911014 DOI: 10.3390/nano11020299] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 01/19/2021] [Accepted: 01/20/2021] [Indexed: 12/13/2022]
Abstract
A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.
Collapse
Affiliation(s)
- Silvio Osella
- Chemical and Biological Systems Simulation Lab, Center of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
| |
Collapse
|
25
|
Kim Y, Bui Y, Tazhigulov RN, Bravaya KB, Slipchenko LV. Effective Fragment Potentials for Flexible Molecules: Transferability of Parameters and Amino Acid Database. J Chem Theory Comput 2020; 16:7735-7747. [PMID: 33236635 DOI: 10.1021/acs.jctc.0c00758] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
An accurate but efficient description of noncovalent interactions is a key to predictive modeling of biological and materials systems. The effective fragment potential (EFP) is an ab initio-based force field that provides a physically meaningful decomposition of noncovalent interactions of a molecular system into Coulomb, polarization, dispersion, and exchange-repulsion components. An EFP simulation protocol consists of two steps, preparing parameters for molecular fragments by a series of ab initio calculations on each individual fragment, and calculation of interaction energy and properties of a total molecular system based on the prepared parameters. As the fragment parameters (distributed multipoles, polarizabilities, localized wave function, etc.) depend on a fragment geometry, straightforward application of the EFP method requires recomputing parameters of each fragment if its geometry changes, for example, during thermal fluctuations of a molecular system. Thus, recomputing fragment parameters can easily become both computational and human bottlenecks and lead to a loss of efficiency of a simulation protocol. An alternative approach, in which fragment parameters are adjusted to different fragment geometries, referred to as "flexible EFP", is explored here. The parameter adjustment is based on translations and rotations of local coordinate frames associated with fragment atoms. The protocol is validated on extensive benchmark of amino acid dimers extracted from molecular dynamics snapshots of a cryptochrome protein. A parameter database for standard amino acids is developed to automate flexible EFP simulations in proteins. To demonstrate applicability of flexible EFP in large-scale protein simulations, binding energies and vertical electron ionization and electron attachment energies of a lumiflavin chromophore of the cryptochrome protein are computed. The results obtained with flexible EFP are in a close agreement with the standard EFP procedure but provide a significant reduction in computational cost.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Yen Bui
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| | - Ruslan N Tazhigulov
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Ksenia B Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
26
|
Abstract
High-order charge transfer is incorporated into the fragment molecular orbital (FMO) method using a charge transfer state with fractional charges. This state is used for a partition analysis of properties based on segments that may be different from fragments in FMO. The partition analysis is also formulated for calculations without fragmentation. All development in this work is limited to density-functional tight-binding. The analysis is applied to a water cluster, crambin (PDB: 1CBN), and two complexes of Trp-cage (1L2Y) with ligands. The contributions of functional groups in ligands are obtained, providing useful information for drug discovery.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
27
|
Viquez Rojas CI, Slipchenko LV. Exchange Repulsion in Quantum Mechanical/Effective Fragment Potential Excitation Energies: Beyond Polarizable Embedding. J Chem Theory Comput 2020; 16:6408-6417. [PMID: 32786899 DOI: 10.1021/acs.jctc.9b01156] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Hybrid quantum mechanical and molecular mechanical (QM/MM) approaches facilitate computational modeling of large biological and materials systems. Typically, in QM/MM, a small region of the system is modeled with an accurate quantum mechanical method and its surroundings with a more efficient alternative, such as a classical force field or the effective fragment potential (EFP). The reliability of QM/MM calculations depends largely on the treatment of interactions between the two subregions, also known as embedding. The polarizable embedding, which allows mutual polarization between solvent and solute, is considered to be essential for describing electronic excitations in polar solvents. In this work, we employ the QM/EFP model and extend the polarizable embedding by incorporating two short-range terms-a charge penetration correction to the electrostatic term and the exchange-repulsion term-both of which are modeled with one-electron contributions to the quantum Hamiltonian. We evaluate the accuracy of these terms by computing excitation energies across 37 molecular clusters consisting of biologically relevant chromophores surrounded by polar solvent molecules. QM/EFP excitation energies are compared to the fully quantum mechanical calculations with the configuration interaction singles (CIS) method. We find that the charge penetration correction diminishes the accuracy of the QM/EFP calculations. On the other hand, while the effect of exchange-repulsion is negligible for most ππ* transitions, the exchange-repulsion significantly improves description of nπ* transitions with blue solvatochromic shifts. As a result, addition of the exchange-repulsion term improves the overall accuracy of QM/EFP. Performances of QM/EFP models remain similar when excitation energies are modeled with cc-pVDZ and aug-cc-pVDZ basis sets.
Collapse
Affiliation(s)
- Claudia I Viquez Rojas
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47906, United States
| |
Collapse
|
28
|
Fedorov DG. Three-Body Energy Decomposition Analysis Based on the Fragment Molecular Orbital Method. J Phys Chem A 2020; 124:4956-4971. [DOI: 10.1021/acs.jpca.0c03085] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
29
|
Kim Y, Morozov D, Stadnytskyi V, Savikhin S, Slipchenko LV. Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex. J Phys Chem Lett 2020; 11:1636-1643. [PMID: 32013435 DOI: 10.1021/acs.jpclett.9b03486] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
High efficiency of light harvesting in photosynthetic pigment-protein complexes is governed by evolutionary-perfected protein-assisted tuning of individual pigment properties and interpigment interactions. Due to the large number of spectrally overlapping pigments in a typical photosynthetic complex, experimental methods often fail to unambiguously identify individual chromophore properties. Here, we report a first-principles-based modeling protocol capable of predicting properties of pigments in protein environment to a high precision. The technique was applied to successfully uncover electronic properties of the Fenna-Matthews-Olson (FMO) pigment-protein complex. Each of the three subunits of the FMO complex contains eight strongly coupled bacteriochlorophyll a (BChl a) pigments. The excitonic structure of FMO can be described by an electronic Hamiltonian containing excitation (site) energies of BChl a pigments and electronic couplings between them. Several such Hamiltonians have been developed in the past based on the information from various spectroscopic measurements of FMO; however, fine details of the excitonic structure and energy transfer in FMO, especially assignments of short-lived high-energy sites, remain elusive. Utilizing polarizable embedding quantum mechanics/molecular mechanics with the effective fragment potentials, we computed the electronic Hamiltonian of FMO that is in general agreement with previously reported empirical Hamiltonians and quantitatively reproduces experimental absorption and circular dichroism spectra of the FMO protein. The developed computational protocol is sufficiently simple and can be utilized for predictive modeling of other wild-type and mutated photosynthetic pigment-protein complexes.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, Jyväskylä 40014, Finland
| | - Valentyn Stadnytskyi
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
- Laboratory of Chemical Physics, National Institute of Diabetes, Digestion and Kidney Diseases, National Institutes of Health, 5 Memorial Drive, Bethesda, Maryland 20892, United States
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
30
|
Herbers S, Fritz SM, Mishra P, Kim Y, Slipchenko L, Zwier TS. The unusual symmetry of hexafluoro-o-xylene-A microwave spectroscopy and computational study. J Chem Phys 2020; 152:064302. [PMID: 32061218 DOI: 10.1063/1.5142169] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The rotational constants and quartic centrifugal distortion coefficients of hexafluoro-o-xylene were precisely derived from the 8 GHz to 18 GHz gas phase microwave spectrum. In addition, the rotational constants of all singly substituted 13C isotopologues were determined. Instead of the intuitively expected symmetry of C2v, as in o-xylene, calculations with a variety of methods (B3LYP, CAM-B3LYP, ωB97XD, MP2, and coupled-cluster singles, doubles, and perturbative triples) predict a C2 symmetry structure in which the two CF3 groups rotate in opposite directions by about 16°. While the experimental results in this study are not capable of proving this unusual symmetry, they can support future microwave, gas phase electron diffraction or nuclear magnetic resonance studies. From the presented data, a preliminary r0 structure was determined, reproducing the experimental rotational constants with deviations of no more than 15 kHz. Analysis of the interactions between the two CF3 groups using an effective fragment potential approach identifies two major contributions to their interaction, due to exchange-repulsion and electrostatic repulsion, with electrostatic repulsion responsible for the barrier at the C2v geometry.
Collapse
Affiliation(s)
- Sven Herbers
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Sean M Fritz
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Piyush Mishra
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Lyudmila Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Timothy S Zwier
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| |
Collapse
|
31
|
Sattasathuchana T, Xu P, Gordon MS. An Accurate Quantum-Based Approach to Explicit Solvent Effects: Interfacing the General Effective Fragment Potential Method with Ab Initio Electronic Structure Theory. J Phys Chem A 2019; 123:8460-8475. [PMID: 31365250 DOI: 10.1021/acs.jpca.9b05801] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
An interface between ab initio quantum mechanics (QM) methods and the general effective fragment potential (EFP2) method, QM-EFP2, is implemented in which the intermolecular interactions between a QM molecule and EFP fragments consist of Coulomb, polarization, exchange repulsion (exrep), and dispersion components. In order to ensure accuracy in the QM-EFP2 exrep interaction energy, the EFP2-EFP2 spherical Gaussian overlap (SGO) approximation is abandoned and replaced with the exact electron repulsion integrals (ERI) that are evaluated with a direct method to reduce disk usage. A Gaussian damping function for the QM-EFP2 Coulomb component damps both the EFP nuclear and electronic charges. A new overlap damping function has been implemented for the QM-EFP2 dispersion component. The current QM-EFP2 implementation has been benchmarked with the S22 and S66 data sets and demonstrates excellent agreement with symmetry-adapted perturbation theory (SAPT) for component energies and with coupled cluster theory [CCSD(T)] for the total interaction energies. Water clusters of various sizes (up to 256 water molecules) have been tested; it is shown that the QM-EFP2 method has an accuracy that is comparable to that of EFP2-EFP2. It has been shown previously that the accuracy of EFP2-EFP2 intermolecular interactions is comparable to that of second-order perturbation theory (MP2) or better. The implementation of the distributed data interface (DDI) parallelization scheme significantly improves the efficiency of QM-EFP2 calculations. The time to form the QM-EFP2 Fock operator per SCF iteration for water clusters scales linearly with the number EFP basis functions.
Collapse
Affiliation(s)
- Tosaporn Sattasathuchana
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| | - Peng Xu
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| | - Mark S Gordon
- Department of Chemistry , Iowa State University and Ames Laboratory Ames , Iowa 50011 , United States
| |
Collapse
|
32
|
Izsák R. Single‐reference coupled cluster methods for computing excitation energies in large molecules: The efficiency and accuracy of approximations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1445] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Róbert Izsák
- Max‐Planck‐Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|
33
|
Fedorov DG. Solvent Screening in Zwitterions Analyzed with the Fragment Molecular Orbital Method. J Chem Theory Comput 2019; 15:5404-5416. [PMID: 31461277 DOI: 10.1021/acs.jctc.9b00715] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Based on induced solvent charges, a new model of solvent screening is developed in the framework of the fragment molecular orbital combined with the polarizable continuum model. The developed model is applied to analyze interactions in a prototypical zwitterionic system, sodium chloride in water, and it is shown that the large underestimation of the interaction in the original solvent screening based on local charges is successfully corrected. The model is also applied to a complex of the Trp-cage (PDB: 1L2Y ) miniprotein with an anionic ligand, and the physical factors determined protein-ligand binding in solution are unraveled.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| |
Collapse
|
34
|
Conrad JA, Kim S, Gordon MS. Ionic liquids from a fragmented perspective. Phys Chem Chem Phys 2019; 21:16878-16888. [PMID: 31359024 DOI: 10.1039/c9cp02836f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
The efficacy of using fragmentation methods, such as the effective fragment potential, the fragment molecular orbital and the effective fragment molecular orbital methods is discussed. The advantages and current limitations of these methods are considered, potential improvements are suggested, and a prognosis for the future is provided.
Collapse
Affiliation(s)
- Justin A Conrad
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| | - Shinae Kim
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| | - Mark S Gordon
- Department of Chemistry, Iowa State University, Ames, IA 50014, USA.
| |
Collapse
|
35
|
Fedorov DG, Brekhov A, Mironov V, Alexeev Y. Molecular Electrostatic Potential and Electron Density of Large Systems in Solution Computed with the Fragment Molecular Orbital Method. J Phys Chem A 2019; 123:6281-6290. [DOI: 10.1021/acs.jpca.9b04936] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan
| | - Anton Brekhov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, 119991, Russian Federation
| | - Yuri Alexeev
- Argonne Leadership Computing Facility and Computational Science Division, Argonne National Laboratory, Argonne, Illinois, 60439, United States
| |
Collapse
|
36
|
Bhattacharyya D, Zhang Y, Elles CG, Bradforth SE. Electronic Structure of Liquid Methanol and Ethanol from Polarization-Dependent Two-Photon Absorption Spectroscopy. J Phys Chem A 2019; 123:5789-5804. [DOI: 10.1021/acs.jpca.9b04040] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Dhritiman Bhattacharyya
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Yuyuan Zhang
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Christopher G. Elles
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Stephen E. Bradforth
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
37
|
Tazhigulov RN, Gurunathan PK, Kim Y, Slipchenko LV, Bravaya KB. Polarizable embedding for simulating redox potentials of biomolecules. Phys Chem Chem Phys 2019; 21:11642-11650. [PMID: 31116217 PMCID: PMC6611676 DOI: 10.1039/c9cp01533g] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana. The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (-0.15 V).
Collapse
Affiliation(s)
- Ruslan N Tazhigulov
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA.
| | | | | | | | | |
Collapse
|
38
|
Wolff FE, Höfener S, Elstner M, Wesołowski TA. Origin of the Solvatochromism in Organic Fluorophores with Flexible Side Chains: A Case Study of Flugi-2. J Phys Chem A 2019; 123:4581-4587. [DOI: 10.1021/acs.jpca.9b02474] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | | | | | - Tomasz A. Wesołowski
- Département de Chimie Physique, Université de Genève, 30 quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
39
|
Nemukhin AV, Grigorenko BL, Khrenova MG, Krylov AI. Computational Challenges in Modeling of Representative Bioimaging Proteins: GFP-Like Proteins, Flavoproteins, and Phytochromes. J Phys Chem B 2019; 123:6133-6149. [DOI: 10.1021/acs.jpcb.9b00591] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Affiliation(s)
- Alexander V. Nemukhin
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
- Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow 119334, Russia
| | - Bella L. Grigorenko
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
- Emanuel Institute of Biochemical Physics, Russian Academy of Sciences, Moscow 119334, Russia
| | - Maria G. Khrenova
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
- Federal Research Center of Biotechnology, Bach Institute of Biochemistry, Russian Academy of Sciences, Moscow 119071, Russian
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
40
|
Haldar S, Riplinger C, Demoulin B, Neese F, Izsak R, Dutta AK. Multilayer Approach to the IP-EOM-DLPNO-CCSD Method: Theory, Implementation, and Application. J Chem Theory Comput 2019; 15:2265-2277. [DOI: 10.1021/acs.jctc.8b01263] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| | | | - Baptiste Demoulin
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Robert Izsak
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| |
Collapse
|
41
|
Sterling CM, Bjornsson R. Multistep Explicit Solvation Protocol for Calculation of Redox Potentials. J Chem Theory Comput 2018; 15:52-67. [PMID: 30511855 DOI: 10.1021/acs.jctc.8b00982] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The calculation of molecular redox potentials in aqueous solution presents a challenge to quantum chemistry due to the need to calculate charged, open-shell species experiencing large solvent effects. Traditionally, redox potentials are calculated via the use of density functional theory and continuum solvation methods, but such protocols have been found to often suffer from large errors, particularly in the case of aqueous solution. While explicit solvation models hold promise of higher accuracy to describe solvent effects in general, their complicated use and lack of well-defined, reliable protocols has hindered their adoption. In this study, we present an explicit-solvation-based approach for the calculation of molecular redox potentials. We combine the use of affordable semiempirical QM/MM molecular dynamics (making use of the recently proposed GFN-xTB method by Grimme et al.) for both redox states and use the linear response approximation to relate vertical ionization energies to the adiabatic redox potential. Simulation length, averaging over snapshots, and accounting for bulk and polarization effects are systematically evaluated using phenol as a working example. We find that it is crucial to reliably account for bulk solvation effects in these calculations, as well as polarization effects which we divide up into short-range and long-range contributions. The short-range polarization contribution is accounted for via QM-region expansion, while the long-range contribution is accounted for via Drude-polarizable QM/MM. Our multistep protocol has been coded to be used in a fully automatic way in a local version of Chemshell. It has been evaluated on a test set of oxidation potentials of organic molecules and found to give gas-solution redox shifts with a mean absolute error of 0.13 eV with respect to experiment, compared to mean absolute errors of 0.26 and 0.21 eV with CPCM and SMD continuum models, respectively.
Collapse
Affiliation(s)
- Cody M Sterling
- Science Institute , University of Iceland , Dunhagi 3 , 107 Reykjavík , Iceland
| | - Ragnar Bjornsson
- Science Institute , University of Iceland , Dunhagi 3 , 107 Reykjavík , Iceland.,Department of Inorganic Spectroscopy , Max-Planck-Institut für Chemische Energiekonversion , Stiftstrasse 34-36 , 45470 Mülheim an der Ruhr , Germany
| |
Collapse
|
42
|
Nanda KD, Krylov AI. The effect of polarizable environment on two-photon absorption cross sections characterized by the equation-of-motion coupled-cluster singles and doubles method combined with the effective fragment potential approach. J Chem Phys 2018; 149:164109. [DOI: 10.1063/1.5048627] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
43
|
Nakata H, Fedorov DG. Analytic second derivatives for the efficient electrostatic embedding in the fragment molecular orbital method. J Comput Chem 2018; 39:2039-2050. [DOI: 10.1002/jcc.25360] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 04/27/2018] [Accepted: 04/29/2018] [Indexed: 01/09/2023]
Affiliation(s)
- Hiroya Nakata
- Department of Fundamental Technology Research; Research and Development Center Kagoshima, Kyocera, 1-4 Kokubu Yamashita-cho; Kirishima-shi Kagoshima, 899-4312 Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono; Tsukuba Ibaraki, 305-8568 Japan
| |
Collapse
|
44
|
Viquez Rojas CI, Fine J, Slipchenko LV. Exchange-repulsion energy in QM/EFP. J Chem Phys 2018; 149:094103. [PMID: 30195305 DOI: 10.1063/1.5043107] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
The effective fragment potential (EFP) is a quantum mechanics (QM)-based model designed to accurately describe intermolecular interactions. Hybrid QM/EFP calculations combine quantum mechanical methods with an EFP embedding to study complex systems in which many-body effects are relevant. As in EFP-only calculations, non-bonded interactions between the QM region and EFP fragments are computed as a sum of electrostatic, polarization, dispersion, and exchange-repulsion energies. The exchange-repulsion term is a computational bottleneck of the EFP calculations. Here, we present a general procedure for computing the QM/EFP exchange-repulsion interactions based on one-electron contributions to the QM Hamiltonian, by using Gaussian functions to represent localized molecular orbitals of the effective fragments. The accuracy of the exchange-repulsion and total QM/EFP interaction energies is evaluated on a diverse set of dimers, including complexes from the S22 dataset of non-covalent interactions. In most cases, the QM/EFP energies are at least as accurate as corresponding EFP energies. A simple and computationally efficient form of the introduced QM/EFP exchange-repulsion term will facilitate further developments and applications of QM/EFP methods.
Collapse
Affiliation(s)
| | - Jonathan Fine
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | | |
Collapse
|
45
|
Empirical corrections and pair interaction energies in the fragment molecular orbital method. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.06.025] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
|
46
|
Gökcan H, Kratz E, Darden TA, Piquemal JP, Cisneros GA. QM/MM Simulations with the Gaussian Electrostatic Model: A Density-based Polarizable Potential. J Phys Chem Lett 2018; 9:3062-3067. [PMID: 29775314 PMCID: PMC6069983 DOI: 10.1021/acs.jpclett.8b01412] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
The use of advanced polarizable potentials in quantum mechanical/molecular mechanical (QM/MM) simulations has been shown to improve the overall accuracy of the calculation. We have developed a density-based potential called the Gaussian electrostatic model (GEM), which has been shown to provide very accurate environments for QM wave functions in QM/MM. In this contribution we present a new implementation of QM/GEM that extends our implementation to include all components (Coulomb, exchange-repulsion, polarization, and dispersion) for the total intermolecular interaction energy in QM/MM calculations, except for the charge-transfer term. The accuracy of the method is tested using a subset of water dimers from the water dimer potential energy surface reported by Babin et al. ( J. Chem. Theory Comput. 2013 9, 5395-5403). Additionally, results of the new implementation are contrasted with results obtained with the classical AMOEBA potential. Our results indicate that GEM provides an accurate MM environment with average root-mean-square error <0.15 kcal/mol for every intermolecular interaction energy component compared with SAPT2+3/aug-cc-pVTZ reference calculations.
Collapse
Affiliation(s)
- Hatice Gökcan
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| | - Eric Kratz
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| | - Thomas A. Darden
- OpenEye Scientific Software, Santa Fe, New Mexico 87508, United States
| | - Jean-Philip Piquemal
- Department of Chemistry, Sorbonne Université, Paris 75005, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
- Insitute Universitaire de France, Paris 75231, France
| | - G. Andrés Cisneros
- Department of Chemistry, University of North Texas, Denton, Texas 76201, United States
| |
Collapse
|
47
|
Fedorov DG. Analysis of solute-solvent interactions using the solvation model density combined with the fragment molecular orbital method. Chem Phys Lett 2018. [DOI: 10.1016/j.cplett.2018.05.002] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
48
|
Odinokov AV, Dubinets NO, Bagaturyants AA. pyEFP: Automatic decomposition of the complex molecular systems into rigid polarizable fragments. J Comput Chem 2018; 39:807-814. [PMID: 29280158 DOI: 10.1002/jcc.25149] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 11/08/2017] [Accepted: 12/08/2017] [Indexed: 11/06/2022]
Abstract
We present an open source tool able to describe intermolecular electrostatic interactions within the framework of the effective fragment potential (EFP) method. Complex molecular structure is subdivided into compact rigid fragments and parameters of their interactions are obtained from ab initio calculations. Automatic procedure allows for searching of these parameters into the existing database and merge new fragments into it. A set of standard fragments useful for the studies of organic semiconductors is also provided. Input files both for purely EFP and hybrid QM/MM calculations can be generated. The program is written in python and freely available on GitHub: https://github.com/ale-odinokov/pyEFP © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Alexey V Odinokov
- Photochemistry Center RAS, Novatorov str., 7a bld. 1, Moscow, 119421, Russia
| | - Nikita O Dubinets
- Photochemistry Center RAS, Novatorov str., 7a bld. 1, Moscow, 119421, Russia.,National Research Nuclear University "MEPhI", Kashirskoe shosse, 31, Moscow, 115409, Russia
| | - Alexander A Bagaturyants
- Photochemistry Center RAS, Novatorov str., 7a bld. 1, Moscow, 119421, Russia.,National Research Nuclear University "MEPhI", Kashirskoe shosse, 31, Moscow, 115409, Russia
| |
Collapse
|
49
|
Morzan UN, Alonso de Armiño DJ, Foglia NO, Ramírez F, González Lebrero MC, Scherlis DA, Estrin DA. Spectroscopy in Complex Environments from QM–MM Simulations. Chem Rev 2018; 118:4071-4113. [DOI: 10.1021/acs.chemrev.8b00026] [Citation(s) in RCA: 97] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Uriel N. Morzan
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Diego J. Alonso de Armiño
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Nicolás O. Foglia
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Francisco Ramírez
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Mariano C. González Lebrero
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Damián A. Scherlis
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| | - Darío A. Estrin
- Departamento de Química Inorgánica, Analítica y Química Física/INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, Pab. II, C1428EHA Buenos Aires, Argentina
| |
Collapse
|
50
|
Nishimoto Y, Fedorov DG. Adaptive frozen orbital treatment for the fragment molecular orbital method combined with density-functional tight-binding. J Chem Phys 2018; 148:064115. [DOI: 10.1063/1.5012935] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Yoshio Nishimoto
- Fukui Institute for Fundamental Chemistry, Kyoto University, 34-4 Takano Nishihiraki-cho, Sakyo-ku, Kyoto 606-8103, Japan
| | - Dmitri G. Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| |
Collapse
|