Prediction of Binding Pose and Affinity of Nelfinavir, a SARS-CoV-2 Main Protease Repositioned Drug, by Combining Docking, Molecular Dynamics, and Fragment Molecular Orbital Calculations.

A novel drug design procedure is described targeting the Main protease (Mpro) of the SARS-CoV-2 virus. The procedure combines molecular docking, molecular dynamics (MD), and fragment molecular orbital (FMO) calculations. The binding structure and properties of Mpro were predicted for Nelfinavir (NFV), which had been identified as a candidate compound through drug repositioning, targeting Mpro.

Dynamic Cooperativity of Ligand-Residue Interactions Evaluated with the Fragment Molecular Orbital Method.

By the splendid advance in computation power realized with the Fugaku supercomputer, it has become possible to perform ab initio fragment molecular orbital (FMO) calculations for thousands of dynamic structures of protein-ligand complexes in a parallel way. We thus carried out electron-correlated FMO calculations for a complex of the 3C-like (3CL) main protease (M) of the new coronavirus (SARS-CoV-2) and its inhibitor N3 incorporating the structural fluctuations sampled by classical molecular dynamics (MD) simulation in hydrated conditions. Along with a statistical evaluation of the interfragment interaction energies (IFIEs) between the N3 ligand and the surrounding amino-acid residues for 1000 dynamic structure samples, in this study we applied a novel approach based on principal component analysis (PCA) and singular value decomposition (SVD) to the analysis of IFIE data in order to extract the dynamically cooperative interactions between the ligand and the residues.

Ca-ATPase Molecules as a Calcium-Sensitive Membrane-Endoskeleton of Sarcoplasmic Reticulum.

The Ca-transport ATPase of sarcoplasmic reticulum (SR) is an integral, transmembrane protein. It sequesters cytoplasmic calcium ions released from SR during muscle contraction, and causes muscle relaxation. Based on negative staining and transmission electron microscopy of SR vesicles isolated from rabbit skeletal muscle, we propose that the ATPase molecules might also be a calcium-sensitive membrane-endoskeleton.

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.

Fragment Molecular Orbital Based Interaction Analyses on COVID-19 Main Protease - Inhibitor N3 Complex (PDB ID: 6LU7).

The worldwide spread of COVID-19 (new coronavirus found in 2019) is an emergent issue to be tackled. In fact, a great amount of works in various fields have been made in a rather short period. Here, we report a fragment molecular orbital (FMO) based interaction analysis on a complex between the SARS-CoV-2 main protease (Mpro) and its peptide-like inhibitor N3 (PDB ID: 6LU7).

Taking Water into Account with the Fragment Molecular Orbital Method.

This chapter describes the current status of development of the fragment molecular orbital (FMO) method for analyzing the electronic state and intermolecular interactions of biomolecular systems in solvent. The orbital energies and the inter-fragment interaction energies (IFIEs) for a specific molecular structure can be obtained directly by performing FMO calculations by exposing water molecules and counterions around biomolecular systems. Then, it is necessary to pay attention to the thickness of the water shell surrounding the biomolecules.

Cm/Eu induced structural, mechanistic and functional implications for calmodulin.

Trivalent actinides and their lanthanide homologues are being scrutinized for their potential health risk when ingested as a result of a range of industrial activities such as mining. Importantly, these ions are known to exhibit high affinity towards calmodulin (CaM). In case of their inadvertent uptake, the holoproteins that are occupied by these cations may block signal transduction pathways or increase the concentration of these ions in intact cells, which could lead to accumulation in human organs.

Potentiating Antigen-Specific Antibody Production with Peptides Obtained from In Silico Screening for High-Affinity against MHC-II.

Monoclonal antibodies with high affinity and specificity are essential for research and clinical purposes, yet remain difficult to produce. Agretope peptides that can potentiate antigen-specific antibody production have been reported recently. Here, we screened in silico for peptides with higher affinity against the agretope binding pocket in the MHC-II.

Hydration of ligands of influenza virus neuraminidase studied by the fragment molecular orbital method.

The fragment molecular orbital (FMO) method was applied to quantum chemical calculations of neuramic acid, the natural substrate of the influenza virus neuraminidase, and two of its competitive inhibitors, Oseltamivir (Tamiful(®)) and Zanamivir (Relenza(®)), to investigate their hydrated structures and energetics. Each of the three ligands was immersed in an explicit water solvent, geometry-optimized by classical MM and QM/MM methods, and subjected to FMO calculations with 2-, 3-, and 4-body corrections under several fragmentation options. The important findings were that QM/MM optimization was preferable to obtain reliable hydrated structures of the ligands, that the 3-body correction was important for quantitative evaluation of the solvation energy, and that the dehydration effect was most remarkable near the hydrophobic sections of the ligands.

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.

FMO-MD simulations on the hydration of formaldehyde in water solution with constraint dynamics.

Full-quantum mechanical fragment molecular orbital-based molecular dynamics (FMO-MD) simulations were applied to the hydration reaction of formaldehyde in water solution under neutral conditions. Two mechanisms, a concerted and a stepwise one, were considered with respect to the nucleophilic addition and the proton transfer. Preliminary molecular orbital calculations by means of polarized continuum model reaction field predicted that the hydration prefers a concerted mechanism.

Fragment Molecular Orbital method-based Molecular Dynamics (FMO-MD) as a simulator for chemical reactions in explicit solvation.

Fragment Molecular Orbital based-Molecular Dynamics (FMO-MD, Komeiji et al., Chem Phys Lett 2003, 372, 342) is an ab initio MD method suitable for large molecular systems. Here, FMO-MD was implemented to conduct full quantum simulations of chemical reactions in explicit solvation.

Water-mediated interactions in the CRP-cAMP-DNA complex: does water mediate sequence-specific binding at the DNA primary-kink site?

The cyclic AMP receptor protein (CRP) of Escherichia coli binds preferentially to DNA sequences possessing a T:A base pair at position 6 (at which the DNA becomes kinked), but with which it does not form any direct interactions. It has been proposed that indirect readout is involved in CRP-DNA binding, in which specificity for this base pair is primarily related to sequence effects on the energetic susceptibility of the DNA to kink formation. In the current study, the possibility of contributions to indirect readout by water-mediated hydrogen bonding of CRP with the T:A base pair was investigated.

Visualization analysis of inter-fragment interaction energies of CRP-cAMP-DNA complex based on the fragment molecular orbital method.

A visualization method for inter-fragment interaction energies (IFIEs) of biopolymers is presented on the basis of the fragment molecular orbital (FMO) method. The IFIEs appropriately illustrate the information about the interaction energies between the fragments consisting of amino acids, nucleotides and other molecules. The IFIEs are usually analyzed in a matrix form called an IFIE matrix.

Change in a protein's electronic structure induced by an explicit solvent: an ab initio fragment molecular orbital study of ubiquitin.

The effect of solvation on the electronic structure of the ubiquitin protein was analyzed using the ab initio fragment molecular orbital (FMO) method. FMO calculations were performed for the protein in vacuo, and the protein was immersed in an explicit solvent shell as thick as 12 A at the HF or MP2 level by using the 6-31G* basis set. The protein's physical properties examined were the net charge, the dipole moment, the internal energy, and the solvent interaction energy.

Intra- and intermolecular interactions between cyclic-AMP receptor protein and DNA: ab initio fragment molecular orbital study.

The ab initio fragment molecular orbital (FMO) calculations were performed for the cAMP receptor protein (CRP) complexed with a cAMP and DNA duplex to elucidate their sequence-specific binding and the stability of the DNA duplex, as determined by analysis of their inter- and intramolecular interactions. Calculations were performed with the AMBER94 force field and at the HF and MP2 levels with several basis sets. The interfragment interaction energies (IFIEs) were analyzed for interactions of CRP-cAMP with each base pair, DNA duplex with each amino acid residue, and each base pair with each residue.

Ab initio fragment molecular orbital (FMO) method applied to analysis of the ligand-protein interaction in a pheromone-binding protein.

Full quantum computation of the electronic state of proteins has recently become possible by the advent of the ab initio fragment molecular orbital (FMO) method. We applied this method to the analysis of the interaction between the Bombyx mori pheromone-binding protein and its ligand, bombykol. The protein-ligand interaction of this molecular complex was minutely analyzed by the FMO method, and the analysis revealed several important interactions between the ligand and amino acid residues.

PEACH 4 with ABINIT-MP: a general platform for classical and quantum simulations of biological molecules.

A program package for molecular simulations of biological molecules was developed. The package, "PEACH version 4 with ABINIT-MP version 20021029," was constructed by incorporating ABINIT-MP, a program for the fragment molecular orbital (FMO) method [Chem.Phys.

Molecular dynamics simulations revealed Ca(2+)-dependent conformational change of Calmodulin.

Molecular dynamics simulations were performed to simulate Ca(2+)-dependent conformational change of calmodulin (CaM). Simulations of the fully Ca(2+)-bound form of CaM (Holo-CaM) and the Ca(2+)-free form (Apo-CaM) were performed in solution for 4 ns starting from the X-ray crystal structure of Holo-CaM. A striking difference was observed between the trajectories of Holo-CaM and Apo-CaM: the central helix remained straight in the former but became largely bent in the latter.