1
|
Fedorov DG. Use of caps in the auxiliary basis set formulation of the fragment molecular orbital method. J Comput Chem 2024; 45:1540-1551. [PMID: 38490813 DOI: 10.1002/jcc.27345] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 02/19/2024] [Accepted: 02/29/2024] [Indexed: 03/17/2024]
Abstract
An auxiliary polarization formulation of the fragment molecular orbital (FMO) method is developed, combining a basis set correction computed for capped isolated fragments with a polarization obtained from uncapped fragments. For a set of organic and inorganic test systems, it is shown that the total energy and atomic charges are accurately reproduced with respect to full unfragmented calculations. It is demonstrated that the method is accurate for computing electronic excited states. The developed approach is applied to rank the isomers of chignolin from experimental NMR data (PDB: 1UAO) according to their relative energy. Contributions of polarization and basis set effects to pair interactions between fragments are elucidated.
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), Tsukuba, Japan
| |
Collapse
|
2
|
Matsuoka S, Sakakura K, Akinaga Y, Akisawa K, Okuwaki K, Doi H, Mochizuki Y. Enhancement of energy decomposition analysis in fragment molecular orbital calculations. J Comput Chem 2024; 45:898-902. [PMID: 38158621 DOI: 10.1002/jcc.27297] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2023] [Revised: 12/10/2023] [Accepted: 12/13/2023] [Indexed: 01/03/2024]
Abstract
Energy decomposition analysis is one of the most attractive features of fragment molecular orbital (FMO) calculations from the point of view of practical applications. Here we report some enhancements for PIEDA in the ABINIT-MP program. One is a separation of the dispersion-type stabilization from the electron correlation energy, traditionally referred to as the "dispersion interaction" (DI). Another is an alternative evaluation of the electrostatic (ES) interaction using the restrained electrostatic potential (RESP) charges. The GA:CT stacked base pair and the Trp-Cage miniprotein were used as illustrative examples.
Collapse
Affiliation(s)
- Sota Matsuoka
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, Tokyo, Japan
| | | | | | - Kazuki Akisawa
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, Tokyo, Japan
| | - Koji Okuwaki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, Tokyo, Japan
- JSOL Corp., Kudan-Kaikan Terrace, Tokyo, Japan
| | - Hideo Doi
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, Tokyo, Japan
| | - Yuji Mochizuki
- Department of Chemistry and Research Center for Smart Molecules, Faculty of Science, Rikkyo University, Tokyo, Japan
- Institute of Industrial Science, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
3
|
Khire SS, Gadre SR. Development and testing of an algorithm for efficient MP2/CCSD(T) energy estimation of molecular clusters with the 2-body approach. J Comput Chem 2023; 44:261-267. [PMID: 35514315 DOI: 10.1002/jcc.26881] [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: 03/03/2022] [Revised: 04/08/2022] [Accepted: 04/13/2022] [Indexed: 12/31/2022]
Abstract
This work reports the development and testing of an automated algorithm for estimating the energies of weakly bound molecular clusters employing correlated theory. Firstly, the monomers and dimers of (homo/hetero) clusters are identified, and the sum of one-body and two-body contributions to correlation energy is calculated. The addition of this contribution to the Hartree-Fock full calculation (FC) energies provides a good estimate of the total energies at Møller-Plesset second-order perturbation theory (MP2)/coupled-cluster method with singles and doubles (CCSD) (T)-level theory using augmented Dunning basis sets. The estimated energies for several test clusters show an excellent agreement with their FC counterparts, with a substantial wall-clock time saving employing off-the-shelf hardware. Furthermore, the complete basis set (CBS) limit for MP2 energy computed using the two-body approach also agrees with its CBS energy with its FC counterpart.
Collapse
Affiliation(s)
- Subodh S Khire
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune, India
| |
Collapse
|
4
|
Fedorov DG. Parametrized quantum-mechanical approaches combined with the fragment molecular orbital method. J Chem Phys 2022; 157:231001. [PMID: 36550057 DOI: 10.1063/5.0131256] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
Fast parameterized methods such as density-functional tight-binding (DFTB) facilitate realistic calculations of large molecular systems, which can be accelerated by the fragment molecular orbital (FMO) method. Fragmentation facilitates interaction analyses between functional parts of molecular systems. In addition to DFTB, other parameterized methods combined with FMO are also described. Applications of FMO methods to biochemical and inorganic systems are reviewed.
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
|
5
|
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
|
6
|
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
|
7
|
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
|
8
|
Chen J, Chan B, Shao Y, Ho J. How accurate are approximate quantum chemical methods at modelling solute-solvent interactions in solvated clusters? Phys Chem Chem Phys 2020; 22:3855-3866. [PMID: 32022044 PMCID: PMC7394230 DOI: 10.1039/c9cp06792b] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
In this paper, the performance of a wide range of DFT methods is assessed for the calculation of interaction energies of thermal clusters of a solute in water. Three different charge states (neutral, proton transfer transition state and zwitterion) of glycine were solvated by 1 to 40 water molecules as sampled from molecular dynamics simulations. While some ab initio composite methods that employ insufficiently large basis sets incurred significant errors even for a cluster containing only 5 water molecules relative to the W1X-2 benchmark, the DLPNO-CCSD(T)/CBS and DSD-PBEP86 (triple zeta basis set) levels of theory predicted very accurate interaction energies. These levels of theory were used to benchmark the performance of 16 density functionals from different rungs of Jacob's Ladder. Of the Rung 4 functionals examined, the ωB97M-V and ωB97X-V functionals stood out for predicting absolute interaction energies in 40-water clusters with mean absolute deviations (MAD) ∼4 kJ mol-1. The B3LYP-D3(BJ) functional performed exceptionally well with a MAD ∼1.7 kJ mol-1 and is the overall best performing method. Calculations of relative interaction energies allow for cancellation of systematic errors, including basis set truncation and superposition errors, and the ωB97M-V and B3LYP-D3(BJ) double zeta basis set calculations yielded relative interaction energies that are within ∼3 kJ mol-1 of the benchmark. The ONIOM approximation provides another strategy for accelerating the calculation of accurate absolute interaction energies provided that the calculations have converged with respect to the size of the "high-level-layer".
Collapse
Affiliation(s)
- Junbo Chen
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-Machi 1-14, Nagasaki 852-8521, Japan.
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Junming Ho
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia.
| |
Collapse
|
9
|
Geometry Optimization, Transition State Search, and Reaction Path Mapping Accomplished with the Fragment Molecular Orbital Method. Methods Mol Biol 2020. [PMID: 32016888 DOI: 10.1007/978-1-0716-0282-9_6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register]
Abstract
Recent development of the fragment molecular orbital (FMO) method related to energy gradients, geometry optimization, transition state search, and chemical reaction mapping is summarized. The frozen domain formulation of FMO is introduced in detail, and the structure of related GAMESS input files for FMO is described.
Collapse
|
10
|
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
|
11
|
Townsend-Nicholson A, Altwaijry N, Potterton A, Morao I, Heifetz A. Computational prediction of GPCR oligomerization. Curr Opin Struct Biol 2019; 55:178-184. [PMID: 31170578 DOI: 10.1016/j.sbi.2019.04.005] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2019] [Revised: 04/16/2019] [Accepted: 04/17/2019] [Indexed: 01/08/2023]
Abstract
There has been a recent and prolific expansion in the number of GPCR crystal structures being solved: in both active and inactive forms and in complex with ligand, with G protein and with each other. Despite this, there is relatively little experimental information about the precise configuration of GPCR oligomers during these different biologically relevant states. While it may be possible to identify the experimental conditions necessary to crystallize a GPCR preferentially in a specific structural conformation, computational approaches afford a potentially more tractable means of describing the probability of formation of receptor dimers and higher order oligomers. Ensemble-based computational methods based on structurally determined dimers, coupled with a computational workflow that uses quantum mechanical methods to analyze the chemical nature of the molecular interactions at a GPCR dimer interface, will generate the reproducible and accurate predictions needed to predict previously unidentified GPCR dimers and to inform future advances in our ability to understand and begin to precisely manipulate GPCR oligomers in biological systems. It may also provide information needed to achieve an increase in the number of experimentally determined oligomeric GPCR structures.
Collapse
Affiliation(s)
- Andrea Townsend-Nicholson
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom.
| | - Nojood Altwaijry
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom
| | - Andrew Potterton
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London WC1E 6BT, United Kingdom; Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| |
Collapse
|
12
|
Heifetz A, James T, Southey M, Morao I, Aldeghi M, Sarrat L, Fedorov DG, Bodkin MJ, Townsend-Nicholson A. Characterising GPCR-ligand interactions using a fragment molecular orbital-based approach. Curr Opin Struct Biol 2019; 55:85-92. [PMID: 31022570 DOI: 10.1016/j.sbi.2019.03.021] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2019] [Revised: 02/19/2019] [Accepted: 03/14/2019] [Indexed: 10/27/2022]
Abstract
There has been fantastic progress in solving GPCR crystal structures. However, the ability of X-ray crystallography to guide the drug discovery process for GPCR targets is limited by the availability of accurate tools to explore receptor-ligand interactions. Visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum mechanical approaches (QM) are often too computationally expensive, but the fragment molecular orbital (FMO) method offers an excellent solution that combines accuracy, speed and the ability to reveal key interactions that would otherwise be hard to detect. Integration of GPCR crystallography or homology modelling with FMO reveals atomistic details of the individual contributions of each residue and water molecule towards ligand binding, including an analysis of their chemical nature.
Collapse
Affiliation(s)
- Alexander Heifetz
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom.
| | - Tim James
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Michelle Southey
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Inaki Morao
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Matteo Aldeghi
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
| | - Laurie Sarrat
- Evotec (France) SAS, 195 Route d' Espagne, 31036 Toulouse, France
| | - 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
| | - Mike J Bodkin
- Evotec (UK) Ltd., 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4RZ, United Kingdom
| | - Andrea Townsend-Nicholson
- Institute of Structural & Molecular Biology, Research Department of Structural & Molecular Biology, Division of Biosciences, University College London, London,WC1E 6BT, United Kingdom
| |
Collapse
|
13
|
Northey T, Kirrander A. Ab Initio Fragment Method for Calculating Molecular X-ray Diffraction. J Phys Chem A 2019; 123:3395-3406. [PMID: 30892904 DOI: 10.1021/acs.jpca.9b00621] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A fragment-based approach for the prediction of elastic X-ray scattering is presented. The total diffraction pattern is assembled from anisotropic form factors calculated for individual molecular fragments, optionally including corrections for pairwise interactions between fragments. The approach is evaluated against full ab initio scattering calculations in the peptide diphenylalanine, and the optimal selection of fragments is examined in the ethanol molecule. The approach is found to improve significantly on the independent atom model while remaining conceptually simple and computationally efficient. It is expected to be particularly useful for macromolecules with repeated subunits, such as peptides, proteins, DNA, or RNA and other polymers, where it is straightforward to define appropriate fragments.
Collapse
Affiliation(s)
- Thomas Northey
- EaStCHEM, School of Chemistry , University of Edinburgh , David Brewster Road , Edinburgh EH9 3FJ , United Kingdom
| | - Adam Kirrander
- EaStCHEM, School of Chemistry , University of Edinburgh , David Brewster Road , Edinburgh EH9 3FJ , United Kingdom
| |
Collapse
|