Search

The fragment molecular orbital (FMO) scheme is one of the popular fragmentation‐based methods and has the potential advantage of making the circuit shallow for quantum chemical calculations on quantum computers. In this study, we used a GPU‐accelerated quantum simulator (cuQuantum) to perform the electron correlation part of the FMO calculation as unitary coupled‐cluster singles and doubles (UCCSD) with the variational quantum eigensolver (VQE) for hydrogen‐bonded (FH) 3 and (FH) 2‐H 2O systems with the STO‐3G basis set. VQE‐UCCSD calculations were performed using both canonical and localized MO sets, and the results were examined from the point of view of size‐consistency and orbital‐invariance affected by the Trotter error. It was found that the use of localized MO leads to better results, especially for (FH) 2‐H 2O. The GPU acceleration was substantial for the simulations with larger numbers of qubits, and was about a factor of 6.7–7.7 for 18 qubit systems. [ABSTRACT FROM AUTHOR]

Abstract Molecular fragmentation is an effective suite of approaches to reduce the formal computational complexity of quantum chemistry calculations while enhancing their algorithmic parallelisability. However, the practical applicability of fragmentation techniques remains hindered by a dearth of automation and effective metrics to assess the quality of a fragmentation scheme. In this article, we present the Quick Fragmentation via Automated Genetic Search (QFRAGS), a novel automated fragmentation algorithm that uses a genetic optimisation procedure to generate molecular fragments that yield low energy errors when adopted in Many Body Expansions (MBEs). Benchmark testing of QFRAGS on protein systems with less than 500 atoms, using two-body (MBE2) and three-body (MBE3) MBE calculations at the HF/6-31G* level, reveals mean absolute energy errors (MAEE) of 20.6 and 2.2 kJ $$\hbox {mol}^{-1}$$ mol - 1 , respectively. For larger protein systems exceeding 500 atoms, MAEEs are 181.5 kJ $$\hbox {mol}^{-1}$$ mol - 1 for MBE2 and 24.3 kJ $$\hbox {mol}^{-1}$$ mol - 1 for MBE3. Furthermore, when compared to three manual fragmentation schemes on a 40-protein dataset, using both MBE and Fragment Molecular Orbital techniques, QFRAGS achieves comparable or often lower MAEEs. When applied to a 10-lipoglycan/glycolipid dataset, MAEs of 7.9 and 0.3 kJ $$\hbox {mol}^{-1}$$ mol - 1 were observed at the MBE2 and MBE3 levels, respectively. Scientific Contribution This Article presents the Quick Fragmentation via Automated Genetic Search (QFRAGS), an innovative molecular fragmentation algorithm that significantly improves upon existing molecular fragmentation approaches by specifically addressing their lack of automation and effective fragmentation quality metrics. With an evolutionary optimisation strategy, QFRAGS actively pursues high quality fragments, generating fragmentation schemes that exhibit minimal energy errors on systems with hundreds to thousands of atoms. The advent of QFRAGS represents a significant advancement in molecular fragmentation, greatly improving the accessibility and computational feasibility of accurate quantum chemistry calculations.

Molecular fragmentation is an effective suite of approaches to reduce the formal computational complexity of quantum chemistry calculations while enhancing their algorithmic parallelisability. However, the practical applicability of fragmentation techniques remains hindered by a dearth of automation and effective metrics to assess the quality of a fragmentation scheme. In this article, we present the Quick Fragmentation via Automated Genetic Search (QFRAGS), a novel automated fragmentation algorithm that uses a genetic optimisation procedure to generate molecular fragments that yield low energy errors when adopted in Many Body Expansions (MBEs). Benchmark testing of QFRAGS on protein systems with less than 500 atoms, using two-body (MBE2) and three-body (MBE3) MBE calculations at the HF/6-31G* level, reveals mean absolute energy errors (MAEE) of 20.6 and 2.2 kJ mol - 1 , respectively. For larger protein systems exceeding 500 atoms, MAEEs are 181.5 kJ mol - 1 for MBE2 and 24.3 kJ mol - 1 for MBE3. Furthermore, when compared to three manual fragmentation schemes on a 40-protein dataset, using both MBE and Fragment Molecular Orbital techniques, QFRAGS achieves comparable or often lower MAEEs. When applied to a 10-lipoglycan/glycolipid dataset, MAEs of 7.9 and 0.3 kJ mol - 1 were observed at the MBE2 and MBE3 levels, respectively. Scientific Contribution This Article presents the Quick Fragmentation via Automated Genetic Search (QFRAGS), an innovative molecular fragmentation algorithm that significantly improves upon existing molecular fragmentation approaches by specifically addressing their lack of automation and effective fragmentation quality metrics. With an evolutionary optimisation strategy, QFRAGS actively pursues high quality fragments, generating fragmentation schemes that exhibit minimal energy errors on systems with hundreds to thousands of atoms. The advent of QFRAGS represents a significant advancement in molecular fragmentation, greatly improving the accessibility and computational feasibility of accurate quantum chemistry calculations. [ABSTRACT FROM AUTHOR]

The dopamine D4 receptor (D4R) is a promising therapeutic target in widespread diseases, and the search for novel agonists and antagonists appears to be clinically relevant. The mechanism of binding to the receptor (R) for antagonists and agonists varies. In the present study, we conducted an in-depth computational study, teasing out key similarities and differences in binding modes, complex dynamics, and binding energies for D4R agonists and antagonists. The dynamic network method was applied to investigate the communication paths between the ligand (L) and G-protein binding site (GBS) of human D4R. Finally, the fragment molecular orbitals with pair interaction energy decomposition analysis (FMO/PIEDA) scheme was used to estimate the binding energies of L–R complexes. We found that a strong salt bridge with D3.32 initiates the inhibition of the dopamine D4 receptor. This interaction also occurs in the binding of agonists, but the change in the receptor conformation to the active state starts with interaction with cysteine C3.36. Such a mechanism may arise in the case of agonists unable to form a hydrogen bond with the serine S5.46, considered, so far, to be crucial in the activation of GPCRs. The energy calculations using the FMO/PIEDA method indicate that antagonists show higher residue occupancy of the receptor binding site than agonists, suggesting they could form relatively more stable complexes. Additionally, antagonists were characterized by repulsive interactions with S5.46 distinguishing them from agonists. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD) is a neurological disease, and its signs and symptoms appear slowly over time. Although current Alzheimer's disease treatments can alleviate symptoms, they cannot prevent the disease from progressing. To accurately diagnose and treat Alzheimer's disease, it is therefore necessary to establish effective methods for diagnosis. Apolipoprotein E4 (ApoE4), the most frequent genetic risk factor for AD, is expressed in more than half of patients with AD, making it an attractive target for AD therapy. We used molecular docking simulations, classical molecular mechanics optimizations, and ab initio fragment molecular orbital (FMO) calculations to investigate the specific interactions between ApoE4 and the naturally occurring compounds found in the plant Moringa Oleifera. According to the FMO calculations, quercetin had the highest binding affinity to ApoE4 among the sixteen compounds because its hydroxyl groups generated strong hydrogen bonds with the ApoE4 residues Trp11, Asp12, Arg15, and Asp130. As a result, we proposed various quercetin derivatives by introducing a hydroxyl group into quercetin and studied their ApoE4 binding properties. The FMO data clearly showed that adding a hydroxyl group to quercetin improved its binding capacity to ApoE4. Furthermore, ApoE4 Trp11, Asp12, Arg15, and Asp130 residues were discovered to be required for significant interactions between ApoE4 and quercetin derivatives. They had a higher ApoE4 binding affinity than our previously proposed epicatechin derivatives. Accordingly, the current results evaluated using the ab initio FMO method will be useful for designing potent ApoE4 inhibitors that can be used as a candidate agent for AD treatment. [ABSTRACT FROM AUTHOR]

The dopamine D4 receptor (D4R) is a promising therapeutic target in widespread diseases, and the search for novel agonists and antagonists appears to be clinically relevant. The mechanism of binding to the receptor (R) for antagonists and agonists varies. In the present study, we conducted an in-depth computational study, teasing out key similarities and differences in binding modes, complex dynamics, and binding energies for D4R agonists and antagonists. The dynamic network method was applied to investigate the communication paths between the ligand (L) and G-protein binding site (GBS) of human D4R. Finally, the fragment molecular orbitals with pair interaction energy decomposition analysis (FMO/PIEDA) scheme was used to estimate the binding energies of L–R complexes. We found that a strong salt bridge with D3.32 initiates the inhibition of the dopamine D4 receptor. This interaction also occurs in the binding of agonists, but the change in the receptor conformation to the active state starts with interaction with cysteine C3.36. Such a mechanism may arise in the case of agonists unable to form a hydrogen bond with the serine S5.46, considered, so far, to be crucial in the activation of GPCRs. The energy calculations using the FMO/PIEDA method indicate that antagonists show higher residue occupancy of the receptor binding site than agonists, suggesting they could form relatively more stable complexes. Additionally, antagonists were characterized by repulsive interactions with S5.46 distinguishing them from agonists.

Alzheimer’s disease (AD) is a neurological disease, and its signs and symptoms appear slowly over time. Although current Alzheimer’s disease treatments can alleviate symptoms, they cannot prevent the disease from progressing. To accurately diagnose and treat Alzheimer’s disease, it is therefore necessary to establish effective methods for diagnosis. Apolipoprotein E4 (ApoE4), the most frequent genetic risk factor for AD, is expressed in more than half of patients with AD, making it an attractive target for AD therapy. We used molecular docking simulations, classical molecular mechanics optimizations, and ab initio fragment molecular orbital (FMO) calculations to investigate the specific interactions between ApoE4 and the naturally occurring compounds found in the plant Moringa Oleifera. According to the FMO calculations, quercetin had the highest binding affinity to ApoE4 among the sixteen compounds because its hydroxyl groups generated strong hydrogen bonds with the ApoE4 residues Trp11, Asp12, Arg15, and Asp130. As a result, we proposed various quercetin derivatives by introducing a hydroxyl group into quercetin and studied their ApoE4 binding properties. The FMO data clearly showed that adding a hydroxyl group to quercetin improved its binding capacity to ApoE4. Furthermore, ApoE4 Trp11, Asp12, Arg15, and Asp130 residues were discovered to be required for significant interactions between ApoE4 and quercetin derivatives. They had a higher ApoE4 binding affinity than our previously proposed epicatechin derivatives. Accordingly, the current results evaluated using the ab initio FMO method will be useful for designing potent ApoE4 inhibitors that can be used as a candidate agent for AD treatment.

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated. [ABSTRACT FROM AUTHOR]

Cytochrome P450 (CYP) enzymes play essential roles in the synthesis and metabolic activation of physiologically active substances. CYP has a prosthetic heme (iron protoporphyrin IX) in its active center, where Fe ion (heme-Fe) is deeply involved in enzymatic reactions of CYP. To precisely describe the structure and electronic states around heme-Fe, we modified the force fields (FFs) around heme-Fe in molecular mechanics (MM) simulations and conducted ab initio fragment molecular orbital (FMO) calculations for the CYP–ligand complex. To describe the coordination bond between heme-Fe and its coordinated ligand (ketoconazole), we added FF between heme-Fe and the N atom of ketoconazole, and then the structure of the complex was optimized using the modified FF. Its adequacy was confirmed by comparing the MM-optimized structure with the X-ray crystal one of the CYP–ketoconazole complex. We also performed 100 ns molecular dynamics simulations and revealed that the coordination bonds around heme-Fe were maintained even at 310 K and that the CYP–ketoconazole structure was kept similar to the X-ray structure. Furthermore, we investigated the electronic states of the complex using the ab initio FMO method to identify the CYP residues and parts of ketoconazole that contribute to strong binding between CYP and ketoconazole. The present procedure of constructing FF between heme-Fe and ketoconazole can be applicable to other CYP–ligand complexes, and the modified FF can provide their accurate structures useful for predicting the specific interactions between CYP and its ligands. Optimized conformation of ketoconazole coordinated to heme-Fe using MM method with bonded FF. The distance between heme-Fe and the N atom of ketoconazole is indicated. [Display omitted] • We studied structure and electronic states of cytochrome P450 (CYP) and ketoconazole complex. • Force fields around heme-Fe of CYP were modified to describe precisely coordination bonds around Fe. • MM optimized structure of CYP−ketoconazole complex was comparable to the X-ray structure. • Stability of the MM optimized structure was confirmed by MD simulations. • Specific interactions between ketoconazole and CYP residues and heme were investigated using ab initio fragment MO method. [ABSTRACT FROM AUTHOR]

Matrix metalloproteinases (MMPs) are calcium-dependent zinc-containing endopeptidases involved in multiple cellular processes. Among the MMP isoforms, MMP-9 regulates cancer invasion, rheumatoid arthritis, and osteoarthritis by degrading extracellular matrix proteins present in the tumor microenvironment and cartilage and promoting angiogenesis. Here, we identified two potent natural product inhibitors of the non-catalytic hemopexin domain of MMP-9 using a novel quantum mechanical fragment molecular orbital (FMO)-based virtual screening workflow. The workflow integrates qualitative pharmacophore modeling, quantitative binding affinity prediction, and a raw material search of natural product inhibitors with the BMDMS-NP library. In binding affinity prediction, we made a scoring function with the FMO method and applied the function to two protein targets (acetylcholinesterase and fibroblast growth factor 1 receptor) from DUD-E benchmark sets. In the two targets, the FMO method outperformed the Glide docking score and MM/PBSA methods. By applying this workflow to MMP-9, we proposed two potent natural product inhibitors (laetanine 9 and genkwanin 10) that interact with hotspot residues of the hemopexin domain of MMP-9. Laetanine 9 and genkwanin 10 bind to MMP-9 with a dissociation constant (KD) of 21.6 and 0.614 μM, respectively. Overall, we present laetanine 9 and genkwanin 10 for MMP-9 and demonstrate that the novel FMO-based workflow with a quantum mechanical approach is promising to discover potent natural product inhibitors of MMP-9, satisfying the pharmacophore model and good binding affinity. [ABSTRACT FROM AUTHOR]

The intrinsic inductive properties of atoms or functional groups depend on the chemical properties of either electron-withdrawing groups (EWGs) or electron-donating groups (EDGs). This study aimed to evaluate in silico methods to determine whether changes in chemical properties of the compound by single atomic substitution affect the biological activity of target proteins and whether the results depend on the properties of the functional groups. We found an imidazo[4,5-b]pyridine-based PAK4 inhibitor, compound 1, as an initial hit compound with the well-defined binding mode for PAK4. In this study, we used both experimental and in silico methods to investigate the effect of atomic substitution on biological activity to optimize the initial hit compound. In biological assays, in the case of EWG, as the size of the halogen atom became smaller and the electronegativity increased, the biological activity IC50 value ranged from 5150 nM to inactive; in the case of EDG, biological activity was inactive. Furthermore, we analyzed the interactions of PAK4 with compounds, focusing on the hinge region residues, L398 and E399, and gatekeeper residues, M395 and K350, of the PAK4 protein using molecular docking studies and fragment molecular orbital (FMO) methods to determine the differences between the effect of EWG and EDG on the activity of target proteins. These results of the docking score and binding energy did not explain the differences in biological activity. However, the pair-interaction energy obtained from the results of the FMO method indicated that there was a difference in the interaction energy between the EWG and EDG in the hinge region residues, L398 and E399, as well as in M395 and K350. The two groups with different properties exhibited opposite electrostatic energy and charge transfer energy between L398 and E399. Additionally, we investigated the electron distribution of the parts interacting with the hinge region by visualizing the molecular electrostatic potential (MEP) surface of the compounds. In conclusion, we described the properties of functional groups that affect biological activity using an in silico method, FMO. [ABSTRACT FROM AUTHOR]

Alkaline protease aeruginolysin (APR) is an important virulence factor in the evasion of the immune system by Pseudomonas aeruginosa (P. aeruginosa). The P. aeruginosa genome also encodes the highly potent and specific APR peptide inhibitor (APRin). However, the structural reason for the significant inhibition has not been revealed. Using ab initio molecular simulations, we here investigated the specific interactions between APR and APRin to elucidate which amino acid residues of APRin and APR contribute strongest to the inhibition. Since APR has a Zn2+ ion at the ligand-binding site and histidine and glutamic acid residues are coordinated with Zn2+, it is essential to precisely describe these coordination bonds to elucidate the specific interactions between APR and APRin. Therefore, we employed the ab initio fragment molecular orbital method to investigate the specific interactions at an electronic level. The results revealed that Ser1 and Ser2 at the N-terminus of APRin significantly contribute to the binding between APRin and APR. In particular, Ser1 binds strongly to Zn2+ as well as to the sidechains of His176(Hid), His180(Hid), and His186(Hid) in APR. This is the main reason for the strong interaction between APR and APRin. The results also elucidated significant contributions of the positively charged Arg83 and Arg90 residues of APRin to the binding with APR. These findings may provide information useful for the design of novel small agents as potent APR inhibitors. [ABSTRACT FROM AUTHOR]

Analyzing the electronic states and inter‐/intra‐molecular interactions of amyloid oligomers expand our understanding of the molecular basis of Alzheimer's disease and other amyloid diseases. In the current study, several high‐resolution crystal structures of oligomeric assemblies of Aβ‐derived peptides have been studied by the ab initio fragment molecular orbital (FMO) method. The FMO method provides comprehensive details of the molecular interactions between the residues of the amyloid oligomers at the quantum mechanical level. Based on the calculations, two sequential aromatic residues (F19 and F20) and negatively charged E22 on the central region of Aβ have been identified as key residues in oligomer stabilization and potential interesting pharmacophores for preventing oligomer formation. [ABSTRACT FROM AUTHOR]

To elucidate structural changes in the retinoic acid receptor-related orphan receptor gamma (RORγt) induced by the binding of an agonist or an inverse agonist, we conducted molecular dynamics (MD) simulations in explicit water. In addition, ab initio fragment molecular orbital calculations were carried out for certain characteristic structures obtained from the MD simulations to reveal important interactions between the amino acid residues of RORγt, and to distinguish the different effects in the binding of an agonist and an inverse agonist on the structure of RORγt. The results elucidate that the hydrogen bond between His479 of helix11 (H11) and Tyr502 of helix12 (H12) is important to keep the H12 conformation in the agonist-bound RORγt. In contrast, in the inverse-agonist-bound RORγt, the side chain of His479 rotates, significantly weakening the interaction between His479 and Tyr502, leading to a conformational change in H12. Therefore, the present molecular simulations clearly indicate that the conformational change in the side chain of His479 in the inverse-agonist-bound RORγt is the main reason for the H12 destabilization induced by the binding of the inverse agonist. Such a conformational change does not occur on the binding of the agonist in RORγt, owing to the strong hydrogen bond between His479 and Tyr502.

Water molecules play a crucial role in protein–ligand binding, and many tools exist that aim to predict the position and relative energies of these important, but challenging participants of biomolecular recognition. The available tools are, in general, capable of predicting the location of water molecules. However, predicting the effects of their displacement is still very challenging. In this work, a linear-scaling quantum mechanics-based approach was used to assess water network energetics and the changes in network stability upon ligand structural modifications. This approach offers a valuable way to improve understanding of SAR data and help guide compound design. [ABSTRACT FROM AUTHOR]

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.

Matrix metalloproteinases (MMPs) are calcium-dependent zinc-containing endopeptidases involved in multiple cellular processes. Among the MMP isoforms, MMP-9 regulates cancer invasion, rheumatoid arthritis, and osteoarthritis by degrading extracellular matrix proteins present in the tumor microenvironment and cartilage and promoting angiogenesis. Here, we identified two potent natural product inhibitors of the non-catalytic hemopexin domain of MMP-9 using a novel quantum mechanical fragment molecular orbital (FMO)-based virtual screening workflow. The workflow integrates qualitative pharmacophore modeling, quantitative binding affinity prediction, and a raw material search of natural product inhibitors with the BMDMS-NP library. In binding affinity prediction, we made a scoring function with the FMO method and applied the function to two protein targets (acetylcholinesterase and fibroblast growth factor 1 receptor) from DUD-E benchmark sets. In the two targets, the FMO method outperformed the Glide docking score and MM/PBSA methods. By applying this workflow to MMP-9, we proposed two potent natural product inhibitors (laetanine 9 and genkwanin 10) that interact with hotspot residues of the hemopexin domain of MMP-9. Laetanine 9 and genkwanin 10 bind to MMP-9 with a dissociation constant (KD) of 21.6 and 0.614 μM, respectively. Overall, we present laetanine 9 and genkwanin 10 for MMP-9 and demonstrate that the novel FMO-based workflow with a quantum mechanical approach is promising to discover potent natural product inhibitors of MMP-9, satisfying the pharmacophore model and good binding affinity.

The intrinsic inductive properties of atoms or functional groups depend on the chemical properties of either electron-withdrawing groups (EWGs) or electron-donating groups (EDGs). This study aimed to evaluate in silico methods to determine whether changes in chemical properties of the compound by single atomic substitution affect the biological activity of target proteins and whether the results depend on the properties of the functional groups. We found an imidazo[4,5-b]pyridine-based PAK4 inhibitor, compound 1, as an initial hit compound with the well-defined binding mode for PAK4. In this study, we used both experimental and in silico methods to investigate the effect of atomic substitution on biological activity to optimize the initial hit compound. In biological assays, in the case of EWG, as the size of the halogen atom became smaller and the electronegativity increased, the biological activity IC50 value ranged from 5150 nM to inactive; in the case of EDG, biological activity was inactive. Furthermore, we analyzed the interactions of PAK4 with compounds, focusing on the hinge region residues, L398 and E399, and gatekeeper residues, M395 and K350, of the PAK4 protein using molecular docking studies and fragment molecular orbital (FMO) methods to determine the differences between the effect of EWG and EDG on the activity of target proteins. These results of the docking score and binding energy did not explain the differences in biological activity. However, the pair-interaction energy obtained from the results of the FMO method indicated that there was a difference in the interaction energy between the EWG and EDG in the hinge region residues, L398 and E399, as well as in M395 and K350. The two groups with different properties exhibited opposite electrostatic energy and charge transfer energy between L398 and E399. Additionally, we investigated the electron distribution of the parts interacting with the hinge region by visualizing the molecular electrostatic potential (MEP) surface of the compounds. In conclusion, we described the properties of functional groups that affect biological activity using an in silico method, FMO.

The accumulation of abnormal prion protein (PrPSc) produced by the structure conversion of PrP (PrPC) in the brain induces prion disease. Although the conversion process of the protein is still not fully elucidated, it has been known that the intramolecular chemical bridging in the most fragile pocket of PrP, known as the "hot spot," stabilizes the structure of PrPC and inhibits the conversion process. Using our original structure-based drug discovery algorithm, we identified the low molecular weight compounds that predicted binding to the hot spot. NPR-130 and NPR-162 strongly bound to recombinant PrP in vitro, and fragment molecular orbital (FMO) analysis indicated that the high affinity of those candidates to the PrP is largely dependent on nonpolar interactions, such as van der Waals interactions. Those NPRs showed not only significant reduction of the PrPSc levels but also remarkable decrease of the number of aggresomes in persistently prion-infected cells. Intriguingly, treatment with those candidate compounds significantly prolonged the survival period of prion-infected mice and suppressed prion disease-specific pathological damage, such as vacuole degeneration, PrPSc accumulation, microgliosis, and astrogliosis in the brain, suggesting their possible clinical use. Our results indicate that in silico drug discovery using NUDE/DEGIMA may be widely useful to identify candidate compounds that effectively stabilize the protein. [ABSTRACT FROM AUTHOR]

Dengue infection is caused by a mosquito-borne virus, particularly in children, which may even cause death. No effective prevention or therapeutic agents to cure this disease are available up to now. The dengue viral envelope (E) protein was discovered to be a promising target for inhibition in several steps of viral infection. Structure-based virtual screening has become an important technique to identify first hits in a drug screening process, as it is possible to reduce the number of compounds to be assayed, allowing to save resources. In the present study, pharmacophore models were generated using the common hits approach (CHA), starting from trajectories obtained from molecular dynamics (MD) simulations of the E protein complexed with the active inhibitor, flavanone (FN5Y). Subsequently, compounds presented in various drug databases were screened using the LigandScout 4.2 program. The obtained hits were analyzed in more detail by molecular docking, followed by extensive MD simulations of the complexes. The highest-ranked compound from this procedure was then synthesized and tested on its inhibitory efficiency by experimental assays. [ABSTRACT FROM AUTHOR]

Tuberculosis (TB), the second leading infectious killer, causes serious public health problems worldwide. To develop novel anti-TB agents, many biochemical studies have targeted the subunit B of DNA gyrase (GyrB), which captures a second DNA segment and responses for ATP hydrolysis. Here, we investigated specific interactions between GyrB residues and existing pyrrolamide derivatives at an electronic level using ab initio fragment molecular orbital (FMO) calculations and designed potent inhibitors against GyrB. The evaluated binding affinities between GyrB and pyrrolamides were confirmed to be consistent with the IC 50 values obtained from previous experiments. Thus, we employed the most potent pyrrolamide (compound 1 ) as a lead compound and proposed novel pyrrolamide derivatives. The specific interactions between GyrB and these derivatives were investigated using molecular mechanic optimizations and FMO calculations. The results revealed that our proposed derivatives had strong hydrogen bonds with Asp79 and Arg141 and exhibited electrostatic interactions with Glu56 and Ile84 of GyrB. In addition, the binding affinity between GyrB and compound 1 was enhanced significantly by the replacement at the R3 site of compound 1 . The present results may provide structural concepts for the rational design of potent GyrB inhibitors as anti-TB agents.Communicated by Ramaswamy H. Sarma.

Aggregates of amyloid-beta proteins (Aβ) have been recognised to be intimately related to pathogenesis of Alzheimer's disease (AD). Indeed, Aβ aggregates of various sizes from dimers to fibrils were found in the brains of AD patients, and these aggregates can be self-organised. Since abnormal accumulation of metal ions such as Zn, Cu and Fe was also observed in the brains, the association between Aβ aggregations and these metal ions has been studied widely. In the present study, to elucidate the influence of Zn ions on the stability of Aβ aggregates, we performed molecular dynamics (MD) simulations and ab initio fragment molecular orbital (FMO) calculations on the Aβ nonamers with and without Zn ions and investigated the change in its structure and electronic states induced by Zn ions at atomic and electronic levels. The MD simulations revealed that Aβ nonamer cannot keep its symmetry structure, whereas Aβ nonamer with Zn ions keeps the structure. The FMO results indicated that electrostatic interactions among the charged amino-acid residues of Aβ nonamer are significantly changed by the influence of Zn ions to stabilise Aβ nonamer. These results provide useful information for proposing novel compounds, which binds specifically to Aβ and inhibits the Aβ aggregation. [ABSTRACT FROM AUTHOR]

Abstract 2-trans enoyl-acyl carrier protein reductase (InhA) has been identified as a promising target for the development of novel chemotherapy for tuberculosis. In the present study, a series of heteroaryl benzamide derivatives were selected as potent inhibitors against InhA, and their binding properties with InhA were investigated at atomic and electronic levels by ab initio molecular simulations based on protein-ligand docking, classical molecular mechanics optimizations and ab initio fragment molecular orbital (FMO) calculations. The results evaluated by FMO highlight some key interactions between InhA and the derivatives, indicating that the most potent derivative has strong hydrogen bonds with the Met98 side chain of InhA and strong electrostatic interactions with the nicotinamide adenine dinucleotide cofactor. These findings provide informative structural concepts for designing novel heteroaryl benzamide derivatives with higher binding affinity to InhA. Graphical abstract Interacting structures between the InhA residues and the derivative A4/B4. Red, blue and black lines indicate hydrogen bonding, electrostatic and π-π interactions, respectively. Image 1 Highlights • Interactions between InhA and its inhibitors were studied by molecular simulations. • Ab initio fragment molecular orbital (FMO) method was used for calculations. • The most potent inhibitor has strong hydrogen bonds with Met98 residue of InhA. • The NH group of the inhibitor is important for forming the hydrogen bonds with Met98. • Novel InhA inhibitors with large binding energy were proposed by FMO calculations. [ABSTRACT FROM AUTHOR]

Graphical abstract Interacting structures between hVDR residues and ligand 2 depending on the protonation states of His305 and His397 Highlights • We investigated specific interactions between vitamin D receptor (VDR) and its ligand. • Ab initio fragment molecular orbital (FMO) method was used for calculations. • Stable protonation states of His residues in VDR were precisely determined by FMO. • The difference in chirality of ligands causes the difference in His protonations. Abstract Vitamin D is recognized to play important roles in the onset of immunological diseases as well as the regulation of the amount of Ca in the blood. Since these physiological actions caused by active vitamin D are triggered by the specific interaction between the vitamin D receptor (VDR) and active vitamin D, many types of compounds have been developed as potent ligands against VDR. It was found that the binding affinity between VDR and its ligand depends significantly on the chirality of the ligand. However, the reason for the dependence has, thus far, not been elucidated. In the present study, we investigated the specific interactions between VDR and some ligands with different chirality, using ab initio fragment molecular orbital (FMO) calculations. The FMO results reveal that two histidine residues of VDR contribute significantly to the binding between VDR and ligand and that their protonation states can affect the specific interactions between VDR and ligand. We therefore considered other possible protonation states of these histidine residues and determined their most stable states, using the ab initio FMO calculations. The results illustrate the possibility that the difference in the chirality of a ligand can induce the change in protonation states of the histidine residues of VDR existing near the ligand. This finding provides an important warning that the protonation states of histidine residues existing near the ligand should be considered more precisely in the molecular simulations for investigating the specific interactions between protein and ligand. [ABSTRACT FROM AUTHOR]

Unveiling the events leading to the formation of prion particles is a nowadays challenge in the field of neurochemistry. Pathogenic mutants of prion protein (PrP) are characterized by both an intrinsic tendency to aggregation and scrapie conversion propensity. However, the question about a possible correlation between these two events lasts still unanswered. Here, a multilayered computational workflow was employed to investigate structure, stability, and molecular interaction properties of a dimer of PrPC‐E200K, a well‐known mutant of the PrP that represents a reduced model of early aggregates of this protein. Based on the combination of molecular dynamics and quantum mechanical approaches, this study provided for an in depth insight of PrPC‐E200K dimer in terms of residue‐residue interactions. Assembly hypotheses for the early aggregation of PrPC‐E200K are paved and compared with PrPSc models reported in the literature to find a structural link between early and late (scrapie) aggregates of this protein. [ABSTRACT FROM AUTHOR]

The inhibition of a bacterial cell division protein, filamentous temperature-sensitive Z (FtsZ), prevents the reproduction of Mycobacteria. To propose potent inhibitors of FtsZ, the binding properties of FtsZ with various derivatives of Zantrin ZZ3 were investigated at an electronic level, using molecular simulations. We here employed protein–ligand docking, classical molecular mechanics (MM) optimizations, and ab initio fragment molecular orbital (FMO) calculations. Based on the specific interactions between FtsZ and the derivatives, as determined by FMO calculations, we proposed novel ligands, which can strongly bind to FtsZ and inhibit its aggregations. The introduction of a hydroxyl group into ZZ3 was found to enhance its binding affinity to FtsZ.

The analytic second derivatives of the energy with respect to nuclear coordinates are developed for restricted Hartree–Fock and density functional theory, based on the two‐body fragment molecular orbital method (FMO) and combined with the electrostatic embedding potential, self‐consistently determined by point charges for far separated fragments and electron densities for near fragments. The accuracy of the method is established with respect to FMO using the exact embedding potential based on electron densities and to full calculations without fragmentation. The computational efficiency of parallelization is measured on the K supercomputer and the method is applied to simulate infrared spectra of two proteins, Trp‐cage (PDB: 1L2Y) and crambin (1CRN). The nature of the vibrations in the Amide I peak of crambin and the Tyr symmetric stretch peak in Trp‐cage are analyzed in terms of localized vibrations. © 2018 Wiley Periodicals, Inc. Vibrational frequencies and infrared spectra of large molecular systems can be evaluated using the developed quantum‐mechanical method. In the application of two proteins, Trp‐cage and crambin, broad bands corresponding to the coupling of multiple quasi‐degenerate normal modes are decoupled using the vibration mode delocalization, yielding detailed information on the vibrations in large molecular systems normally hidden by the broad vibrational bands. A map of the couplings between localized modes is shown. [ABSTRACT FROM AUTHOR]

Abstract: The crystal structure of metabotropic glutamate receptor 1 (mGluR1) complexed with 4‐fluoro‐N‐(4‐(6‐(isopropylamino)pyrimidin‐4‐yl)thiazol‐2‐yl)‐N‐methylbenzamide (FITM, a negative allosteric modulator) and its twelve close structural analogs with a broad spectrum of affinities (2.4 nM < IC50 > 10 000 nM) were investigated using quantum mechanical methods. The our own N‐layered integrated molecular orbital and molecular mechanics (ONIOM) was used to optimize the molecular geometries of the receptor with complexed ligands, which were then used to perform the ab initio calculations using the fragment molecular orbitals method with energy decomposition analysis (FMO‐EDA). The results clearly showed that residues Q6603.28 and/or Y8056.55 were the anchoring points for all the studied analogs of FITM, while the H‐bond with T8157.38 determined only the orientation of very active molecules containing an amino substituent in the pyrimidine moiety (e.g., FITM). The orientation of the other parts of ligands resulted from hydrophobic interactions mainly with L7575.44, F8016.51, or W7986.48. The applied ONIOM/FMO–EDA approach facilitated the study of effects related to very small changes in the ligand structure and led to conclusions regarding the significance of individual interactions in the allosteric binding pocket of mGluR1. [ABSTRACT FROM AUTHOR]

Vitamin D plays an important role in the regulation of the calcium and phosphorus metabolism as well as in bone formation. These physiological actions caused by vitamin D are triggered by the specific binding of vitamin D to its receptor (VDR). Here we investigated the specific interactions and binding affinities between VDR and vitamin D derivatives, using ab initio fragment molecular orbital (FMO) calculations. The FMO results elucidate that relative position of the two hydroxyl groups of the derivatives is essential for the strong binding affinity between the derivative and Arg274 residue of VDR. It is therefore expected that novel potent ligands, which have a great binding affinity for VDR, are developed by adjusting the positions of the hydroxyl groups in the derivatives in such a way as these groups form strong hydrogen bonds with VDR residues. We proposed these novel derivatives and investigated their specific interactions with VDR at atomic and electronic levels to obtain a more potent ligand for VDR. [ABSTRACT FROM AUTHOR]

Spectrophotometric titration of buffered solutions of gamma cyclodextrin ( γ CD) and 1,1 ′ -diethyl,2,2 ′ -dicarbocyanine (DDI) demonstrates extension of the known 1:2 host:guest complex to form a previously unreported 2:2 complex near the γ CD solubility limit. Though DDI is predominantly hosted as a non-fluorescent H-aggregate, both complexes exist in respective equilibria with two secondary complexes hosting unaggregated DDI as 1:1 and 2:1 complexes. The 2:1 complex exhibits significant fluorescence emission, with a quantum yield six times that of DDI in organic solvents, but ten times lower than that of an analogous indodicarbocyanine. Fragment Molecular Orbital calculations suggest that the 2:1 complex has the tail-to-tail conformation, and that solvent access to the dye strongly favors photoisomerization. In the host-guest complex, γ CD limits solvent access to the dye and hinders rotation of the quinolyl terminal groups, but nevertheless pairwise rotation of methine carbons within the γ CD cavity likely remains as a significant nonradiative relaxation pathway for the excited state. [ABSTRACT FROM AUTHOR]

Filamentous temperature-sensitive Z (FtsZ) protein plays essential role in bacteria cell division, and its inhibition prevents Mycobacteria reproduction. Here we adopted curcumin derivatives as candidates of novel inhibitors and investigated their specific interactions with FtsZ, using ab initio molecular simulations based on protein–ligand docking, classical molecular mechanics and ab initio fragment molecular orbital (FMO) calculations. Based on FMO calculations, we specified the most preferable site of curcumin binding to FtsZ and highlighted the key amino acid residues for curcumin binding at an electronic level. The result will be useful for proposing novel inhibitors against FtsZ based on curcumin derivatives. [ABSTRACT FROM AUTHOR]

Accumulation of amyloid-β (Aβ) peptides in a brain is closely related with the pathogenesis of Alzheimer's disease. To suppress the production of Aβ peptides, we propose novel curcumin derivatives and investigate their binding properties with the amyloid precursor protein (APP), using protein-ligand docking as well as ab initio molecular simulations. Our proposed derivative (curcumin XIV) is found to have a large binding energy with APP and interacts strongly with the cleavage site Ala19 by secretase. It is thus expected that curcumin XIV can protect APP from the secretase attack and be a potent inhibitor against the production of Aβ peptides. [ABSTRACT FROM AUTHOR]

Bacteria secrete the enzyme pseudolysin (PLN) to degrade exocellular proteins, and the produced peptides are used as a nutrient for the bacteria. It is thus expected that inhibition of PLN can suppress bacterial growth. Since such inhibitors do not attack to bacteria directly, the risk of producing drug-resistance bacteria is less. However, endogenous proteinases such as the matrix metalloproteinases (MMPs) have active site similar to that of PLN, and there is a possibility that PLN inhibitors also inhibit the activity of MMPs, resulting in a loss of substrate degradation by these proteinases. Therefore, agents that inhibit the activity of only PLN and not MMPs are required. In the present study, we select two compounds (ARP101 and LM2) and investigate their specific interactions with PLN and MMPs by use of ab initio molecular simulations. Based on the results, we propose several novel compounds as candidates for potent PLN inhibition and investigate their binding properties with PLN, elucidating that the compound, in which a toluene group is introduced into LM2, has larger binding energy with PLN compared with the pristine LM2. Therefore, this compound is suggested to be a potent PLN inhibitor. [ABSTRACT FROM AUTHOR]

The Androgen Receptor (AR) is a family of nuclear receptor proteins and a ligand-activated transcription factor. Since its abnormal activation can cause the progression of prostate cancer, numerous types of antagonists against AR have been developed as promising agents for treating prostate cancers. We here investigated the specific interactions between AR and several types of non-steroid agents at an electronic level, using ab initio molecular simulations based on molecular mechanics and ab initio fragment molecular orbital (FMO) methods From the results obtained by FMO, we proposed novel agents as potent ligands against AR and investigated the binding properties between AR and these agents to confirm that some of them can bind more strongly with AR than the existing non-steroid agents and can be strongly effective ligands against AR. [ABSTRACT FROM AUTHOR]

Vitamin D is recognized to play important roles not only in the bone metabolism and the regulation of Ca amount in the blood but also in the onset of immunological diseases. These physiological actions caused by vitamin D are triggered by the specific interaction between vitamin D receptor (VDR) and vitamin D. In the present study, we investigated the interactions between VDR and vitamin D derivatives using ab initio molecular simulation, in order to elucidate the reason for the significant difference in their effects on VDR activity. Based on the results simulated, we elucidated which parts of the derivatives and which residues of VDR mainly contribute to the specific binding between VDR and the derivatives at an electronic level. This finding will be helpful for proposing new vitamin D derivatives as a potent modulator or inhibitor against VDR. [ABSTRACT FROM AUTHOR]

The accumulation of amyloid-beta (Aβ) aggregates in brain contributes to the onset of Alzheimer's disease (AD). Recent structural analysis for the tissue obtained from AD patients revealed that Aβ aggregates have a single structure with three-fold symmetry. To explain why this structure possesses significant stability, we here investigated the specific interactions between Aβ peptides in the aggregate, using ab initio fragment molecular orbital calculations. The results indicate that the interactions between the Aβ peptides of the stacked Aβ pair are stronger than those between the Aβ peptides of the trimer with three-fold symmetry and that the charged amino-acids are important. [ABSTRACT FROM AUTHOR]

A series of fragment molecular orbital (FMO) calculations on DNA models was performed to assess the reliability of the scaled third-order Møller-Plesset perturbation method (MP2.5) (Pitoňák et al., Chem. Phys. Chem., 10 (2009) 282) for the incorporation of higher-order electron correlation. The Watson-Crick H-bonding (horizontal interaction) and the 1,2-stacking (vertical interaction) of four combinations of two base pairs (AT:AT, AT:TA, AT:GC, and CG:GC) were calculated at the levels of MP2, MP3 (nonscaled), MP2.5, and coupled cluster singles and doubles with perturbative triples (CCSD(T)). The comparison showed sufficient agreement between MP2.5 and the reference CCSD(T), and indicated that MP2.5 is an inexpensive but precise alternative to costly CCSD(T). The influence of basis sets including the model core potential (MCP) was also examined for the two base pairs. The Watson-Crick H-bonding interaction energies were fairly uninfluenced by these methods. In contrast, the values of the stacking interactions among bases quantitatively varied depending on the calculation options. As a cost-effective procedure, the FMO-MP2.5/MCP calculations were applied to a couple of realistic DNA models explicitly solvated by water: a 12-base-pair model and a 10-base-pair model cross-linked with cis- platin (Pt(NH 3 ) 2 ). The H-bonding and stacking interactions in them were discussed from a biochemical viewpoint. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD), one of the most common neurodegenerative diseases, is a major factor contributing to cognitive impairment in older adults. Current therapeutic treatments can only relieve the symptoms of AD, but they cannot stop the progression of the disease because it takes a long time for clinical symptoms to manifest. Therefore, it is essential to develop effective diagnostic strategies for early detection and treatment of AD. As the most common genetic risk factor for AD, apolipoprotein E4 (ApoE4) is present in more than half of patients with AD, and it can be a target protein for AD therapy. We used molecular docking, classical molecular mechanics optimizations, and ab initio fragment molecular orbital (FMO) calculations to investigate the specific interactions between ApoE4 and Cinnamon-derived compounds. Of the 10 compounds, epicatechin was found to have the highest binding affinity to ApoE4 because the hydroxyl groups of epicatechin form strong hydrogen bonds with the Asp130 and Asp12 residues of ApoE4. Therefore, we proposed some epicatechin derivatives by adding a hydroxyl group to epicatechin and studied their interactions with ApoE4. The FMO results indicate that the addition of a hydroxyl group to epicatechin increases its binding affinity to ApoE4. It is also revealed that the Asp130 and Asp12 residues of ApoE4 are important for the binding between ApoE4 and the epicatechin derivatives. These findings will help propose potent inhibitors against ApoE4, leading to a proposal for effective therapeutic candidates for AD. Outline for the design of novel natural ApoE4 inhibitors from the spice Cinnamon based on ab initio molecular simulations. [Display omitted] • The interactions between ApoE4 and the natural compounds in Cinnamon were investigated. • We performed molecular docking and ab initio fragment molecular orbital calculations. • The most potent compound against ApoE4 was epicatechin. • We designed novel epicatechin derivatives and examined their interactions with ApoE4. • Asp12 and Asp130 of ApoE4 were crucial for its strong interaction with ligand. [ABSTRACT FROM AUTHOR]