Partial geometry optimization with FMO-MP2 gradient: Application to TrpCage

Halogenated Baicalein as a Promising Antiviral Agent toward SARS-CoV-2 Main Protease

The coronavirus disease pandemic is a constant reminder that global citizens are in imminent danger of exposure to emerging infectious diseases. Therefore, developing a technique for inhibitor discovery is essential for effective drug design. Herein, we proposed fragment molecular orbital (FMO)-based virtual screening to predict the molecular binding energy of potential severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease inhibitors. The integration of quantum mechanical approaches and trajectory analysis from a microsecond molecular dynamics simulation was used to identify potential inhibitors. We identified brominated baicalein as a potent inhibitor of the SARS-CoV-2 main protease and confirmed its inhibitory activity in an in vitro assay. Brominated baicalein did not demonstrate significant toxicity in either in vitro or in vivo studies. The pair interaction energy from FMO-RIMP2/PCM and inhibitory constants based on the protease enzyme assay suggested that the brominated baicalein could be further developed into novel SARS-CoV-2 protease inhibitors.

Interaction of 8-anilinonaphthalene-1-sulfonate with SARS-CoV-2 main protease and its application as a fluorescent probe for inhibitor identification

The 3C-like main protease of SARS-CoV-2 (3CLPro) is responsible for the cleavage of the viral polyprotein. This process is essential for the viral life cycle. Therefore, 3CLPro is a promising target to develop antiviral drugs for COVID-19 prevention and treatment. Traditional enzymatic assays for the identification of 3CLPro inhibitors rely on peptide-based colorimetric or fluorogenic substrates. However, the COVID-19 pandemic has limit or delay access to these substrates, especially for researchers in developing countries attempting to screen natural product libraries. We explored the use of the fluorescent probe 8-anilinonaphthalene-1-sulfonate (ANS) as an alternative assay for inhibitor identification. Fluorescence enhancement upon binding of ANS to 3CLPro was observed, and this interaction was competitive with a peptide substrate. The utility of ANS-based competitive binding assay to identify 3CLPro inhibitors was demonstrated with the flavonoid natural products baicalein and rutin. The molecular nature of ANS and rutin interaction with 3CLPro was explored with molecular modeling. Our results suggested that ANS could be employed in a competitive binding assay to facilitate the identification of novel SARS-CoV-2 antiviral compounds.

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

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

Theoretical Analysis of Activity Cliffs among Benzofuranone-Class Pim1 Inhibitors Using the Fragment Molecular Orbital Method with Molecular Mechanics Poisson-Boltzmann Surface Area (FMO+MM-PBSA) Approach

Significant activity changes due to small structural changes (i.e., activity cliffs) of serine/threonine kinase Pim1 inhibitors were studied theoretically using the fragment molecular orbital method with molecular mechanics Poisson-Boltzmann surface area (FMO+MM-PBSA) approach. This methodology enables quantum-chemical calculations for large biomolecules with solvation. In the course of drug discovery targeting Pim1, six benzofuranone-class inhibitors were found to differ only in the position of the indole-ring nitrogen atom. By comparing the various qualities of complex structures based on X-ray, classical molecular mechanics (MM)-optimized, and quantum/molecular mechanics (QM/MM)-optimized structures, we found that the QM/MM-optimized structures provided the best correlation (R² = 0.85) between pIC₅₀ and the calculated FMO+MM-PBSA binding energy. Combining the classical solvation energy with the QM binding energy was important to increase the correlation. In addition, decomposition of the interaction energy into various physicochemical components by pair interaction energy decomposition analysis suggested that CH-π and electrostatic interactions mainly caused the activity differences.

Subsystem Analysis for the Fragment Molecular Orbital Method and Its Application to Protein-Ligand Binding in Solution

A subsystem analysis is derived incorporating interfragment interactions into the fragment properties, such as energies or charges. The relative stabilities of three alanine isomers, the α-helix, the β-turn, and the extended form are studied and the differences in fragment properties are elucidated. The analysis is further elaborated for studies of binding energies. The binding of the Trp-cage protein (PDB: 1L2Y ) to two ligands is studied in detail. Binding energies defined for each fragment can be used as a convenient descriptor for analyzing contributions to binding in solution.

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.

CH-π hydrogen bonds in biological macromolecules

This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.

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

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