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

Modification of MM force fields around heme-Fe in the CYP-ligand complex and ab initio FMO calculations for the complex

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.

The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH4)2 and Ne2

We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.

Protein-Protein Interaction Modelling with the Fragment Molecular Orbital Method

Fragment molecular orbital (FMO) method enables ab initio quantum-chemical calculations for biomolecular systems with high accuracy and moderate computational cost. Through this analysis we can evaluate the inter-fragment interaction energies (IFIEs) that provide useful measures for effective interactions between the fragments representing amino-acid residues and ligand molecules. Here I describe how to prepare the input structures and perform the FMO calculations for protein-protein complex system. In addition to the pre-processing, some useful tools for the post-processing analysis are also illustrated.

Structure-based drug design of novel M. tuberculosis InhA inhibitors based on fragment molecular orbital calculations

2-trans enoyl-acyl carrier protein reductase (InhA) is a promising target for developing novel chemotherapy agents for tuberculosis, and their inhibitory effects on InhA activity were widely investigated by the physicochemical experiments. However, the reason for the wide range of their inhibitory effects induced by similar agents was not explained by only the difference in their chemical structures. In our previous molecular simulations, a series of heteroaryl benzamide derivatives were selected as candidate inhibitors against InhA, and their binding properties with InhA were investigated to propose novel derivatives with higher binding affinity to InhA. In the present study, we extended the simulations for a series of 4-hydroxy-2-pyridone derivatives to search widely for more potent inhibitors against InhA. Using ab initio fragment molecular orbital (FMO) calculations, we elucidated the specific interactions between InhA residues and the derivatives at an electronic level and highlighted key interactions between InhA and the derivatives. The FMO results clearly indicated that the most potent inhibitor has strong hydrogen bonds with the backbones of Tyr158, Thr196, and NADH of InhA. This finding may provide informative structural concepts for designing novel 4-hydroxy-2-pyridone derivatives with higher binding affinity to InhA. Our previous and present molecular simulations could provide important guidelines for the rational design of more potent InhA inhibitors.

Natural inhibitors for severe acute respiratory syndrome coronavirus 2 main protease from Moringa oleifera, Aloe vera, and Nyctanthes arbor-tristis: molecular docking and ab initio fragment molecular orbital calculations

The novel coronavirus 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread rapidly worldwide, and new drug treatments for COVID-19 are urgently required. To find the potential inhibitors against the main protease (Mpro) of SARS-CoV-2, we investigated the inhibitory potential of naturally occurring compounds from the plants Moringa oleifera, Aloe vera, and Nyctanthes arbor-tristis, using molecular docking, classical molecular mechanics optimizations, and ab initio fragment molecular orbital (FMO) calculations. Of the 35 compounds that we simulated, feralolide from Aloe vera exhibited the highest binding affinity against Mpro. Therefore, we proposed novel compounds based on the feralolide and investigated their binding properties to Mpro. The FMO results indicated that the introduction of a hydroxyl group into feralolide significantly enhances its binding affinity to Mpro. These results provide useful information for developing potent Mpro inhibitors.

Specific interactions between the alkaline protease of P. aeruginosa and its natural peptide inhibitor: ab initio molecular simulations

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 Zn²⁺ ion at the ligand-binding site and histidine and glutamic acid residues are coordinated with Zn²⁺, 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 Zn²⁺ 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.

Proposal of selective inhibitor for bacterial zinc metalloprotease: Molecular mechanics and ab initio molecular orbital calculations

The zinc metalloprotease pseudolysin (PLN) secreted from Pseudomonas aeruginosa degrades extracellular proteins to produce bacterial nutrition, and various types of PLN inhibitors have been developed to suppress the bacterial growth. However, as the structure of the ligand-binding pocket of PLN has large similarities to those of human matrix metalloproteinases (MMPs) and other human zinc metalloprotease, there is a risk that PLN inhibitors also inhibit human zinc proteases. In this study, we propose a novel agent that may bind stronger to PLN than to MMPs. The compound is proposed based on the specific molecular interactions between existing agents and PLN/MMP metalloproteases evaluated by the present molecular simulations. First, we confirmed that the binding energies of PLN agents evaluated using the ab initio fragment molecular orbital method were comparable to the IC₅₀ values obtained through previous experiments. In addition, the specific molecular interactions between these agents and MMP-9 were investigated to elucidate the fact that some of the agents bind weaker to MMP than PLN. Based on the results, we proposed a novel agent having a succinimide group introduce by a hydroxamic acid group and investigated its binding properties with PLN and MMP. The results may provide useful information for the development of potent inhibitors for PLN with few potential side effects in human bodies.

Proposal of novel natural inhibitors of severe acute respiratory syndrome coronavirus 2 main protease: Molecular docking and ab initio fragment molecular orbital calculations

This paper proposes natural drug candidate compounds for the treatment of coronavirus disease 2019 (COVID-19). We investigated the binding properties between the compounds in the Moringa oleifera plant and the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 using molecular docking and ab initio fragment molecular orbital calculations. Among the 12 compounds, niaziminin was found to bind the strongest to Mpro. We furthermore proposed novel compounds based on niaziminin and investigated their binding properties to Mpro. The results reveal that the introduction of a hydroxyl group into niaziminin enhances its binding affinity to Mpro. These niaziminin derivatives can be promising candidate drugs for the treatment of COVID-19.

Proposal of novel inhibitors for vitamin-D receptor: Molecular docking, molecular mechanics and ab initio molecular orbital simulations

The specific binding of active vitamin-D to the vitamin-D receptor (VDR) is closely related to the onset of immunological diseases. To inhibit the binding, various compounds have been developed as potent inhibitors against VDR. Among them, a compound NS-54c, which was developed based on the first VDR antagonist TEI-9647 (25-dehydro-1α-hydroxyvitamin D₃-26,23-lactone), was revealed to posse almost 1000-fold improved antagonistic activity over the original TEI-9647. However, the reason for this significant improvement has not been elucidated. In the present study, we investigated the specific interactions between VDR and these inhibitors, using molecular simulations based on molecular docking, molecular mechanics and ab initio fragment molecular orbital calculations. Based on the results simulated, we furthermore proposed novel inhibitors and investigated their binding properties to VDR. The results elucidate that the replacement of propyl group at the 24th site of NS-54c by a phenethyl group can enhance the binding affinity of the inhibitor to VDR. This finding provides useful information for developing novel potent inhibitors against VDR.

Proposal of Potent Inhibitors for a Bacterial Cell Division Protein FtsZ: Molecular Simulations Based on Molecular Docking and ab Initio Molecular Orbital Calculations

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.

Identification of correlated inter-residue interactions in protein complex based on the fragment molecular orbital method

A theoretical scheme to systematically describe correlated (network-like) interactions between molecular fragments is proposed within the framework of the fragment molecular orbital (FMO) method. The method is mathematically based on the singular value decomposition (SVD) of the inter-fragment interaction energy (IFIE) matrix obtained by the FMO calculation, and can be applied to a comprehensive description of protein-protein interactions in the context of molecular recognition. In the present study we apply the proposed method to a complex of measles virus hemagglutinin and human SLAM receptor, thus finding a usefulness for efficiently eliciting the correlated interactions among the amino-acid residues involved in the two proteins. Additionally, collective interaction networks by amino-acid residues important for mutation experiments can be clearly visualized.

Structural change of retinoic-acid receptor-related orphan receptor induced by binding of inverse-agonist: Molecular dynamics and ab initio molecular orbital simulations

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.

Protonation states of central amino acids in a zinc metalloprotease complexed with inhibitor: Molecular mechanics optimizations and ab initio molecular orbital calculations

The zinc-metalloprotease pseudolysin (PLN) secreted from bacteria degrades extracellular proteins to produce bacterial nutrition. Since PLN has a Zn ion at the inhibitor-binding site, the interactions between Zn and PLN residues as well as inhibitor can be significantly changed depending on the protonation states of PLN residues at the inhibitor-binding site. To determine stable protonation states of these residues, we here considered different protonation states for Glu and His residues located around Zn and investigated the electronic states of the PLN + inhibitor complex, using ab initio molecular simulations. The protonation state of His223 was found to significantly affect the specific interactions between PLN and the inhibitor.

Specific interactions between tau protein and curcumin derivatives: Molecular docking and ab initio molecular orbital simulations

Alzheimer's disease (AD) is the most common neurodegenerative disorder in the world, and there is currently no potent medicine for the treatment of ADs. Curcumin, a primary chemical contained in the ancient Indian herb known as turmeric, has been extensively studied and shown to be effective in inhibiting the aggregations of amyloid-β and tau proteins, both of which are observed in the brains of AD patients. In the present study, we focused on the tau protein and investigated its specific interactions with curcumin derivatives, using molecular simulations based on molecular docking, molecular mechanics and ab initio fragment molecular orbital calculations. Based on the results, we attempted to propose novel potent inhibitors against the tau protein aggregation. Our molecular simulations provide useful information for developing novel medicines for the treatment of ADs.

Characterization of PD-L1 binding sites by a combined FMO/GRID-DRY approach

The programmed cell death protein 1 (PD-1) and its ligand, PD-L1, constitute an important co-inhibitory immune checkpoint leading to downregulation of immune system. Tumor cells developed a strategy to trigger PD-1/PD-L1 pathway reducing the T cell anticancer activity. Anti-PD-L1 small drugs, generally with improved pharmacokinetic and technological profiles than monoclonal antibodies, became an attractive research topic. Nevertheless, still few works have been published on the chemical features of possible binding sites. In this work, we applied a novel computational protocol based on the combination of the ab initio Fragment Molecular Orbital (FMO) method and a newly developed GRID-DRY approach in order to characterize the PD-L1 binding sites, starting from PD-1/PD-L1 and PD-L1/BMS-ligands (Bristol-Mayers Squibb ligands) complexes. The FMO method allows the calculation of the pair-residues as well as the ligand-residues interactions with ab initio accuracy, whereas the GRID-DRY approach is an effective tool to investigate hydrophobic interactions, not easily detectable by ab initio methods. The present GRID-DRY protocol is able to determine the energy contributions of each ligand atoms to each hydrophobic interaction, both qualitatively and quantitatively. We were also able to identify the three specific hot regions involved in PD-1/PD-L1 protein-protein interaction and in PD-L1/BMS-ligand interactions, in agreement with preceding theoretical/experimental results, and to suggest a specific pharmacophore for PD-L1 inhibitors.

Design of galardine analogs as putative psudolysin inhibitors based on ab initio fragment molecular orbital calculations

Pseudolysin (PLN) is a metalloproteinase secreted from bacteria that degrades extracellular proteins to produce bacterial nutrition. It is thus expected that inhibitors against PLN can suppress the growth of bacteria and their pandemic spread. In addition, since these inhibitors do not attack to bacteria directly, there is a reduced risk for producing drug-resistant bacteria. On the other hand, as PLN has large structural similarity in the active sites with human matrix-metalloproteinases (MMPs), there is a possibility that the inhibitors for PLN also inhibit MMP activity, resulting in a loss of necessary nutrients to be produced by MMPs. Therefore, it is required the agents inhibiting the activity of only PLN not MMPs. In the present study, we employed a hydroxamate compound galardin, which has a significant inhibition effect against PLN and MMP, and investigated its specific interactions with PLN/MMP at atomic and electronic levels, by use of ab initio molecular simulations. Based on the results, we proposed several derivatives of galardin and elucidated which derivatives that can bind more strongly to PLN and be putative antimicrobial agents capable of inhibiting the PLN activity.Communicated by Ramaswamy H. Sarma.

Specific interactions between 2-trans enoyl-acyl carrier protein reductase and its ligand: Protein-ligand docking and ab initio fragment molecular orbital calculations

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.

Specific interactions between zinc metalloproteinase and its inhibitors: Ab initio fragment molecular orbital calculations

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.

Investigation of the molecular similarity in closely related protein systems: The PrP case study

The amyloid conversion is a massive detrimental modification affecting several proteins upon specific physical or chemical stimuli characterizing a plethora of diseases. In many cases, the amyloidogenic stimuli induce specific structural features to the protein conferring the propensity to misfold and form amyloid deposits. The investigation of mutants, structurally similar to their native isoform but inherently prone to amyloid conversion, may be a viable strategy to elucidate the structural features connected with amyloidogenesis. In this article, we present a computational protocol based on the combination of molecular dynamics (MD) and grid-based approaches suited for the pairwise comparison of closely related protein structures. This method was applied on the cellular prion protein (PrP(C)) as a case study and, in particular, addressed to the quali/quantification of the structural features conferred by either E200K mutations and treatment with CaCl(2), both able to induce the scrapie conversion of PrP. Several schemes of comparison were developed and applied to this case study, and made up suitable of application to other protein systems. At this purpose an in-house python codes has been implemented that, together with the parallelization of the GRID force fields program, will spread the applicability of the proposed computational procedure.

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.

Comparison of the local structural stabilities of mammalian prion protein (PrP) by fragment molecular orbital calculations

Bovine spongiform encephalopathy (BSE), a member of the prion diseases, is a fatal neurodegenerative disorder suspected to be caused by a malfunction of prion protein (PrP). Although BSE prions have been reported to be transmitted to a wide range of animal species, dogs and hamsters are known to be BSE-resistant animals. Analysis of canine and hamster PrP could elucidate the molecular mechanisms supporting the species barriers to BSE prion transmission. The structural stability of 6 mammalian PrPs, including human, cattle, mouse, hamster, dog and cat, was analyzed. We then evaluated intramolecular interactions in PrP by fragment molecular orbital (FMO) calculations. Despite similar backbone structures, the PrP side-chain orientations differed among the animal species examined. The pair interaction energies between secondary structural elements in the PrPs varied considerably, indicating that the local structural stabilities of PrP varied among the different animal species. Principal component analysis (PCA) demonstrated that different local structural stability exists in bovine PrP compared with the PrP of other animal species examined. The results of the present study suggest that differences in local structural stabilities between canine and bovine PrP link diversity in susceptibility to BSE prion infection.

Molecular dynamics simulation of the induced-fit binding process of DNA aptamer and L-argininamide

Aptamers are rare functional nucleic acids with binding affinity to and specificity for target ligands. Recent experiments have lead to the proposal of an induced-fit binding mechanism for L-argininamide (Arm) and its binding aptamer. However, at the molecular level, this mechanism between the aptamer and its coupled ligand is still poorly understood. The present study used explicit solvent molecular dynamics (MD) simulations to examine the critical bases involved in aptamer-Arm binding and the induced-fit binding process at atomic resolution. The simulation results revealed that the Watson-Crick pair (G10-C16), C9, A12, and C17 bases play important roles in aptamer-Arm binding, and that binding of Arm results in an aptamer conformation optimized through a general induced-fit process. In an aqueous solution, the mechanism has the following characteristic stages: (a) adsorption stage, the Arm anchors to the binding site of aptamer with strong electrostatic interaction; (b) binding stage, the Arm fits into the binding site of aptamer by hydrogen-bond formation; and (c) complex stabilization stage, the hydrogen bonding and electrostatic interactions cooperatively stabilize the complex structure. This study provides dynamics information on the aptamer-ligand induced-fit binding mechanism. The critical bases in aptamer-ligand binding may provide a guideline in aptamer design for molecular recognition engineering.