Counterpoise-corrected interaction energy analysis based on the fragment molecular orbital scheme

Use of caps in the auxiliary basis set formulation of the fragment molecular orbital method

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.

Hydrophobic core formation and secondary structure elements in uranyl(VI)-binding peptides

Cyclic peptides as well as a modified EF-hand motif of calmodulin have been newly designed to achieve high affinity towards uranyl(VI). Cyclic peptides may be engineered to bind uranyl(VI) to its backbone under acidic conditions, which may enhance its selectivity. For the modified EF-hand motif of calmodulin, strong electrostatic interactions between uranyl(VI) and negatively charged side chains play an important role in achieving high affinity; however, it is also essential to have a secondary structure element and formation of hydrophobic cores in the metal-bound state of the peptide.

Interaction analyses of SARS-CoV-2 spike protein based on fragment molecular orbital calculations

At the stage of SARS-CoV-2 infection in human cells, the spike protein consisting of three chains, A, B, and C, with a total of 3300 residues plays a key role, and thus its structural properties and the binding nature of receptor proteins to host human cells or neutralizing antibodies has attracted considerable interest. Here, we report on interaction analyses of the spike protein in both closed (PDB-ID: 6VXX) and open (6VYB) structures, based on large-scale fragment molecular orbital (FMO) calculations at the level of up to the fourth-order Møller–Plesset perturbation with singles, doubles, and quadruples (MP4(SDQ)). Inter-chain interaction energies were evaluated for both structures, and a mutual comparison indicated considerable losses of stabilization energies in the open structure, especially in the receptor binding domain (RBD) of chain-B. The role of charged residues in inter-chain interactions was illuminated as well. By two separate calculations for the RBD complexes with angiotensin-converting enzyme 2 (ACE2) (6M0J) and B38 Fab antibody (7BZ5), it was found that the binding with ACE2 or antibody partially compensated for this stabilization loss of RBD.

Visualized IFIE results seen from chain-B of spike protein.

Electron-correlated fragment-molecular-orbital calculations for biomolecular and nano systems

Recent developments in the fragment molecular orbital (FMO) method for theoretical formulation, implementation, and application to nano and biomolecular systems are reviewed. The FMO method has enabled ab initio quantum-mechanical calculations for large molecular systems such as protein-ligand complexes at a reasonable computational cost in a parallelized way. There have been a wealth of application outcomes from the FMO method in the fields of biochemistry, medicinal chemistry and nanotechnology, in which the electron correlation effects play vital roles. With the aid of the advances in high-performance computing, the FMO method promises larger, faster, and more accurate simulations of biomolecular and related systems, including the descriptions of dynamical behaviors in solvent environments. The current status and future prospects of the FMO scheme are addressed in these contexts.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.

Open-shell pair interaction energy decomposition analysis (PIEDA): formulation and application to the hydrogen abstraction in tripeptides

An open-shell extension of the pair interaction energy decomposition analysis (PIEDA) within the framework of the fragment molecular orbital (FMO) method is developed. The open-shell PIEDA method allows the analysis of inter- and intramolecular interactions in terms of electrostatic, exchange-repulsion, charge-transfer, dispersion, and optional polarization energies for molecular systems with a radical or high-spin fragment. Taking into account the low computational cost and scalability of the FMO and PIEDA methods, the new scheme provides a means to characterize the stabilization of radical and open-shell sites in biologically relevant species. The open-shell PIEDA is applied to the characterization of intramolecular interactions in capped trialanine upon hydrogen abstraction (HA) at various sites on the peptide. Hydrogen abstraction reaction is the first step in the oxidative pathway initiated by reactive oxygen or nitrogen species, associated with oxidative stress. It is found that HA results in significant geometrical reorganization of the trialanine peptide. Depending on the HA site, terminal interactions in the radical fold conformers may become weaker or stronger compared to the parent molecule, and often change the character of the non-covalent bonding from amide stacking to hydrogen bonding.

An Efficient Method to Evaluate Intermolecular Interaction Energies in Large Systems Using Overlapping Multicenter ONIOM and the Fragment Molecular Orbital Method

We propose an approach based on the overlapping multicenter ONIOM to evaluate intermolecular interaction energies in large systems and demonstrate its accuracy on several representative systems in the complete basis set limit at the MP2 and CCSD(T) level of theory. In the application to the intermolecular interaction energy between insulin dimer and 4'-hydroxyacetanilide at the MP2/CBS level, we use the fragment molecular orbital method for the calculation of the entire complex assigned to the lowest layer in three-layer ONIOM. The developed method is shown to be efficient and accurate in the evaluation of the protein-ligand interaction energies.

Reassessment of Hammett σ as an effective parameter representing intermolecular interaction energy-links between traditional and modern QSAR approaches

The Hammett σ constant has for a long time been known to be one of most important linear free-energy related parameters that correlate with biological activity. It is a conventionally used electronic parameter in studies of enzymatic quantitative structure-activity relationships (QSAR). However, it is not necessarily obvious why σ represents variations in the free-energy change associated with the complex formation between a congeneric series of ligands with their target protein. So far, several powerful molecular calculations, such as the ab initio fragment molecular orbital (FMO) one, that are directly applicable to ligand-protein complexes have emerged. In this study, we comprehensively reevaluate experimentally derived parameter σ confirming it represents intermolecular interaction energy terms, by applying molecular orbital (MO) calculations to a simple ligand-protein complex model. The current results provide a rational and quantitative basis for bridging the gap between the traditional QSAR approach and 'the modern QSAR one', which involves the molecular calculations to evaluate the overall free-energy change for complex formation.