Ab Initio Molecular Orbital Methods Research Articles

Combining the Fragment Molecular Orbital and GRID Approaches for the Prediction of Ligand-Metalloenzyme Binding Affinity: The Case Study of hCA II Inhibitors.

Polarization and charge-transfer interactions play an important role in ligand-receptor complexes containing metals, and only quantum mechanics methods can adequately describe their contribution to the binding energy. In this work, we selected a set of benzenesulfonamide ligands of human Carbonic Anhydrase II (hCA II)-an important druggable target containing a Zn2+ ion in the active site-as a case study to predict the binding free energy in metalloprotein-ligand complexes and designed specialized computational methods that combine the ab initio fragment molecular orbital (FMO) method and GRID approach. To reproduce the experimental binding free energy in these systems, we adopted a machine-learning approach, here named formula generator (FG), considering different FMO energy terms, the hydrophobic interaction energy (computed by GRID) and logP. The main advantage of the FG approach is that it can find nonlinear relations between the energy terms used to predict the binding free energy, explicitly showing their mathematical relation. This work showed the effectiveness of the FG approach, and therefore, it might represent an important tool for the development of new scoring functions. Indeed, our scoring function showed a high correlation with the experimental binding free energy (R2 = 0.76-0.95, RMSE = 0.34-0.18), revealing a nonlinear relation between energy terms and highlighting the relevant role played by hydrophobic contacts. These results, along with the FMO characterization of ligand-receptor interactions, represent important information to support the design of new and potent hCA II inhibitors.

Performance and Mechanism of Chlorine Dioxide on BTEX Removal in Liquid and Indoor Air.

With the development of the chemical industry, benzene, toluene, ethylbenzene, and xylene (BTEX) have gradually become the major indoor air pollutants. Various gas treatment techniques are widely used to prevent the physical and mental health hazards of BTEX in semi-enclosed spaces. Chlorine dioxide (ClO2) is an alternative to chlorine as a secondary disinfectant with a strong oxidation ability, a wide range of action, and no carcinogenic effects. In addition, ClO2 has a unique permeability which allows it to eliminate volatile contaminants from the source. However, little attention has been paid to the removal of BTEX by ClO2, due to the difficulty of removing BTEX in semi-enclosed areas and the lack of testing methods for the reaction intermediates. Therefore, this study explored the performance of ClO2 advanced oxidation technology on both liquid and gaseous benzene, toluene, o-xylene, and m-xylene. The results showed that ClO2 was efficient in the removal of BTEX. The byproducts were detected by gas chromatography-mass spectrometry (GC-MS) and the reaction mechanism was speculated using the ab initio molecular orbital calculations method. The results demonstrated that ClO2 could remove the BTEX from the water and the air without causing secondary pollution.

Theoretical Study of Si/C Equally Mixed Dodecahedrane Analogues

To gain insight into the effect of Si/C arrangement on the molecular framework of strained polyhedral compounds, dodecahedrane analogues containing equal numbers of carbon and silicon (Si/C equally mixed dodecahedrane analogues) were investigated using the ab initio molecular orbital method. There are 1648 isomers for which the Si/C arrangement on the molecular framework is different. Based on the geometry optimization of all the isomers as well as the carbon and silicon analogues, the characteristics of the structure, relative energy, and strain energy of the Si/C equally mixed compounds are presented. Then, important factors controlling the relative energy, such as strain energy, are proposed through regression analysis. Also discussed is the correlation between the relative energy and the indices of Si/C dispersion, such as the number of skeletal C-Si single bonds and condensed five-membered rings constituting the polyhedral structure.

Numerical Analysis for Light Absorption Spectra of the Base of DNA-Wrapped Single-Walled Carbon Nanotubes

This study numerically demonstrates the light absorption spectra of each base of DNA-wrapped single-walled carbon nanotubes (SWCNTs). Previous experimental and theoretical studies show that the optical properties of these composites are different from the bare SWCNTs. In this work, we investigated the bases of DNA that influence optical properties. To obtain stable molecular states for studying optical properties, molecular dynamics calculations were performed. Additionally, light absorption spectra in the ultraviolet-to-near-infrared region of one type of base-wrapped (e.g., adenine-, thymine-, cytosine-, or guanine-wrapped) SWCNTs were investigated by utilizing the semi-empirical molecular orbital theory using SCIGRESS commercial software. This method can significantly reduce the calculation time compared to the ab initio molecular orbital method, making the handling of composites of bases and SWCNTs possible. We found that the largest peaks appear at a wavelength of around 300 nm for all the composites. Furthermore, we found that the light absorption spectra above 570 nm are strongly influenced by adenine and cytosine. Thus, our computational results provide insight into the optical properties and the effects of base-SWCNTs that are difficult to investigate experimentally under the influence of solvents and various molecules.

Predicting Kinetics and Dynamics of Spin-Dependent Processes

ConspectusPredicting mechanisms and rates of nonadiabatic spin-dependent processes including photoinduced intersystem crossings, thermally activated spin-forbidden reactions, and spin crossovers in metal centers is a very active field of research. These processes play critical roles in transition-metal-based and metalloenzymatic catalysis, molecular magnets, light-harvesting materials, organic light-emitting diodes, photosensitizers for photodynamic therapy, and many other applications. Therefore, accurate modeling of spin-dependent processes in complex systems and on different time scales is important for many problems in chemistry, biochemistry, and materials sciences.Nonadiabatic statistical theory (NAST) and nonadiabatic molecular dynamics (NAMD) are two complementary approaches to modeling the kinetics and dynamics of spin-dependent processes. NAST predicts the probabilities and rate constants of nonradiative transitions between electronic states with different spin multiplicities using molecular properties at only few critical points on the potential energy surfaces (PESs), including the reactant minimum and the minimum energy crossing point (MECP) between two spin states. This makes it possible to obtain molecular properties for NAST calculations using accurate but often computationally expensive electronic structure methods, which is critical for predicting the rate constants of spin-dependent processes. Alternatively, NAST can be used to study spin-dependent processes in very large complex molecular systems using less computationally expensive electronic structure methods. The nuclear quantum effects, such as zero-point vibrational energy, tunneling, and interference between reaction paths can be easily incorporated. However, the statistical and local nature of NAST makes it more suitable for large systems and slow kinetics. In contrast, NAMD explores entire PESs of interacting electronic states, making it ideal for modeling fast barrierless spin-dependent processes. Because the knowledge of large portions of PESs is often needed, the simulations require a very large number of electronic structure calculations, which limits the NAMD applicability to relatively small molecular systems and ultrafast kinetics.In this Account, we discuss our contribution to the development of the NAST and NAMD approaches for predicting the rates and mechanism of spin-dependent processes. First, we briefly describe our NAST and NAMD implementations. The NAST implementation is an extension of the transition state theory to the processes involving two crossing potential energy surfaces of different spin multiplicities. The NAMD approach includes the trajectory surface hopping (TSH) and ab initio multiple spawning (AIMS) methods. Second, we discuss several applications of NAST and NAMD to model spin-dependent processes in different systems. The NAST applicability to large complex systems is demonstrated by the studies of the spin-forbidden isomerization of the active sites of metal-sulfur proteins. Our implementation of the MECP search algorithm within the fully ab initio fragment molecular orbital method allows applying NAST to systems with thousands of atoms, such as the solvated protein rubredoxin. Applications of NAMD to ultrafast spin-dependent processes are represented by the generalized AIMS simulations utilizing the fast GPU-based TeraChem electronic structure program to gain insight into the complex photoexcited state relaxation in 2-cyclopentenone.

The FMO2 analysis of the ligand-receptor binding energy: the Biscarbene-Gold(I)/DNA G-Quadruplex case study.

In this work, the ab initio fragment molecular orbital (FMO) method was applied to calculate and analyze the binding energy of two biscarbene-Au(I) derivatives, [Au(9-methylcaffein-8-ylidene)2]+ and [Au(1,3-dimethylbenzimidazol-2-ylidene)2]+, to the DNA G-Quadruplex structure. The FMO2 binding energy considers the ligand-receptor complex as well as the isolated forms of energy-minimum state of ligand and receptor, providing a better description of ligand-receptor affinity compared with simple pair interaction energies (PIE). Our results highlight important features of the binding process of biscarbene-Au(I) derivatives to DNA G-Quadruplex, indicating that the total deformation-polarization energy and desolvation penalty of the ligands are the main terms destabilizing the binding. The pair interaction energy decomposition analysis (PIEDA) between ligand and nucleobases suggest that the main interaction terms are electrostatic and charge-transfer energies supporting the hypothesis that Au(I) ion can be involved in π-cation interactions further stabilizing the ligand-receptor complex. Moreover, the presence of polar groups on the carbene ring, as C = O, can improve the charge-transfer interaction with K+ ion. These findings can be employed to design new powerful biscarbene-Au(I) DNA-G quadruplex binders as promising anticancer drugs. The procedure described in this work can be applied to investigate any ligand-receptor system and is particularly useful when the binding process is strongly characterized by polarization, charge-transfer and dispersion interactions, properly evaluated by ab initio methods.

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 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.

Interaction Analysis on the SARS-CoV-2 Spike Protein Receptor Binding Domain Using Visualization of the Interfacial Electrostatic Complementarity.

Visualization of the interfacial electrostatic complementarity (VIINEC) is a recently developed method for analyzing protein–protein interactions using electrostatic potential (ESP) calculated via the ab initio fragment molecular orbital method. In this Letter, the molecular interactions of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with human angiotensin-converting enzyme 2 (ACE2) and B38 neutralizing antibody were examined as an illustrative application of VIINEC. The results of VIINEC revealed that the E484 of RBD has a role in making a local electrostatic complementary with ACE2 at the protein–protein interface, while it causes a considerable repulsive electrostatic interaction. Furthermore, the calculated ESP map at the interface of the RBD/B38 complex was significantly different from that of the RBD/ACE2 complex, which is discussed herein in association with the mechanism of the specificity of the antibody binding to the target protein.

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 IC50 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.

Hot Spot Analysis of YAP-TEAD Protein-Protein Interaction Using the Fragment Molecular Orbital Method and Its Application for Inhibitor Discovery.

Simple SummaryYAP-TEAD protein-protein interaction is a crucial feature in the Hippo pathway for cancer therapeutics. In this study, to investigate novel YAP-TEAD protein-protein interaction inhibitors, the fragment molecular orbital method for hot spot analysis of YAP-TEAD PPI and FMO-based pharmacophore construction were introduced. We performed pharmacophore based virtual screening through an in-house library. Virtual hit compounds were derived, and the biological activity of the hit compounds was confirmed, through in vitro experiments. Finally, we proposed a compound BY03 as a novel promising PPI inhibitor for the YAP-TEAD interaction.The Hippo pathway is an important signaling pathway modulating growth control and cancer cell proliferation. Dysregulation of the Hippo pathway is a common feature of several types of cancer cells. The modulation of the interaction between yes-associated protein (YAP) and transcriptional enhancer associated domain (TEAD) in the Hippo pathway is considered an attractive target for cancer therapeutic development, although the inhibition of PPI is a challenging task. In order to investigate the hot spots of the YAP and TEAD1 interacting complex, an ab initio Fragment Molecular Orbital (FMO) method was introduced. With the hot spots, pharmacophores for the inhibitor design were constructed, then virtual screening was performed to an in-house library. Next, we performed molecular docking simulations and FMO calculations for screening results to study the binding modes and affinities between PPI inhibitors and TEAD1. As a result of the virtual screening, three compounds were selected as virtual hit compounds. In order to confirm their biological activities, cellular (luciferase activity, proximity ligation assay and wound healing assay in A375 cells, qRT-PCR in HEK 293T cells) and biophysical assays (surface plasmon resonance assays) were performed. Based on the findings of the study, we propose a novel PPI inhibitor BY03 and demonstrate a profitable strategy to analyze YAP–TEAD PPI and discover novel PPI inhibitors.

Identification of key stabilizing interactions of amyloid-β oligomers based on fragment molecular orbital calculations on macrocyclic β-hairpin peptides.

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.

Comparison of group 14 elements in sp3 and sp2 environment by fragment structure energy analysis.

To explore the characteristics of group 14 elements in sp3 and sp2 hybrid environments systematically, the concept of fragment structure energy (FSE) was proposed through ab initio molecular orbital method. A fragment in this study means a building block of molecular structure and the energy of each kind of fragments (FSE) is obtained from the total electronic energy of different size of model compounds by just a simple algebraic method. The model compounds are the hydrogen-terminated diamond and graphene analogs as the typical compounds in the ideal sp3 - and sp2 -environment, respectively. As a result, the sp3 environment was found to be more stable than the sp2 for all the group 14 elements (C-Pb) in the ideal hybrid environment. Furthermore, a mixing effect of carbon and heavier group 14 elements (E=Si, Ge, Sn, and Pb) was investigated by comparing the E/C-alternately mixed- and the pure FSEs.

Proposal of novel potent inhibitors against androgen receptor based on ab initio molecular orbital calculations

The androgen receptor (AR), a family of nuclear receptor proteins, stimulates the transcription of androgen-responsive genes. As its abnormal activation can cause the progression of prostate cancer, numerous types of ligands for AR have been developed as promising antagonists for the treatment of prostate cancer. We previously investigated the specific interactions between AR and nine types of existing non-steroidal ligands, using molecular simulations based on molecular mechanics and ab initio fragment molecular orbital methods. The results were confirmed to be comparable to the binding affinities of these ligands observed in experiments. We here propose novel ligands as potent inhibitors against AR and investigate their binding properties to AR, using the same molecular simulations. The results indicate that the most promising ligand binds stronger to AR than the existing non-steroidal ligands, and that our proposed ligand binds strongly to a mutant-type AR, which has drug resistance to the existing non-steroidal ligands.

Acceleration of Environmental Electrostatic Potential Using Cholesky Decomposition with Adaptive Metric (CDAM) for Fragment Molecular Orbital (FMO) Method

Abstract The calculation speed of the ab initio fragment molecular orbital (FMO) method can and must be increased by applying approximations to the environmental electrostatic potential (ESP) and the dimer electrostatic potential (dimer-es). These approximations were previously implemented by introducing the Cholesky decomposition with adaptive metric (CDAM) method to the FMO (Okiyama et al., Chem. Phys. Lett.2010, 490, 84). In this study, a distributed memory algorithm of CDAM was introduced to reduce the necessary memory size. The improved version of CDAM was applied to the ESP approximation and was shown to give sufficiently precise energy values while halving the calculation time.

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.

Application of Fragment Molecular Orbital Method to investigate dopamine receptors

GPCRs are a vast family of seven-domain transmembrane proteins. This family includes dopamine receptors (D1, D2, D3, D4, and D5), which mediate the variety of dopamine-controlled physiological functions in the brain and periphery. Ligands of dopamine receptors are used for managing several neuropsychiatric disorders, including bipolar disorder, schizophrenia, anxiety, and Parkinson’s disease. Recent studies have revealed that dopamine receptors could be part of multiple signaling cascades, rather than of a single signaling pathway. For these targets, a variety of experimental and computational drug design techniques are utilized. In this work, dopamine receptors D2, D3, and D4 were investigated using molecular dynamic method as well as computational ab initio Fragment Molecular Orbital method (FMO), which can reveal atomistic details about ligand binding. The results provided useful insights into the significances of amino acid residues in ligand binding sites. Moreover, similarities and differences between active-sites of three studied types of receptors were examined.

Proposal of therapeutic curcumin derivatives for Alzheimer’s disease based on ab initio molecular simulations

Aggregation of amyloid-β peptides (Aβs) in a brain is deeply involved in the onset of Alzheimer’s disease, and various agents inhibiting the aggregation have been developed based on natural substances to suppress the onset. We here proposed novel curcumin derivatives as potent inhibitors against the aggregation and investigated their binding properties to Aβ, using protein-ligand docking and ab initio fragment molecular orbital methods. The results reveal that a curcumin derivative, in which COH3 group of the aromatic ring is replaced by OH, binds strongly to Aβ and can be a potent inhibitor against the Aβ aggregation.

Change in binding states between catabolite activating protein and DNA induced by ligand-binding: molecular dynamics and ab initio fragment molecular orbital calculations.

The transcription mechanism of genetic information from DNA to RNA is efficiently controlled by regulatory proteins, such as catabolite activator protein (CAP), and their ligands. When cyclic AMP (cAMP) binds to CAP, the complex forms a dimer and binds specifically to DNA to activate the transcription mechanism. On the other hand, when cyclic GMP (cGMP) binds to CAP, the complex has no marked effect on the mechanism. In our previous study, based on molecular dynamics (MD) and ab initio fragment molecular orbital (FMO) methods, we elucidated which residues of CAP are important for the specific interactions between CAP and DNA in the CAP-monomer+DNA + cAMP complex. However, this monomer model for CAP cannot describe real interactions between the CAP-dimer and DNA because CAPs form a dimer before binding to DNA. Accordingly, here, we investigated stable structures and their electronic states for the CAP-dimer+DNA complex with cAMP or cGMP ligand, to clarify the influence of ligand-binding on the interactions between CAP-dimer and DNA. The MD simulations elucidated that the DNA-binding domains of CAP-dimer behave differently depending on the ligand bound to the CAP-dimer. In addition, FMO calculations revealed that the binding energy between CAP-dimer and DNA for the CAP-dimer+DNA + cAMP complex is larger than that for the CAP-dimer+DNA + cGMP complex, being consistent with experiments. It was also highlighted that the Arg185 and Lys188 residues of CAP-dimer are important for the binding between CAP-dimer and DNA. These results provide useful information for proposing new compounds that efficiently control the transcription mechanism.

Ab Initio Molecular Orbital Study of the First Four Si/C Alternately Substituted Annulenes.

The ground and low-lying excited states of four alternating Si/C annulenes, H nSi n/2C n/2 with n = 4, 6, 8, and 10, have been investigated by ab initio molecular orbital methods and compared to those of their all-carbon and all-silicon analogues. In the ground state, all of the Si/C-mixed annulenes, except for the largest 10-membered annulene (H10Si5C5), assume equal-bond-length structures by adopting a closed-shell electronic structure in the possible highest symmetry. For the largest H10Si5C5, the trend of the bond delocalization still remains but the circular structure is considerably distorted and nonplanar due to severe angle strain. In the low-lying singlet (S1) and triplet (T1) states, the geometry of the compounds tends to be nonplanar as the excitations produce silyl radical character. Relative energies of the T1 and S1 states of the 6-membered ring, compared to those of the respective ground states (S0), are higher than those of the 4- and 8-membered rings, suggesting a special stability for H6Si3C3. The planar rhombus shape of the formally antiaromatic H4Si2C2 suggests that a synthetic effort is merited. Bonding analyses are given to support the conclusions reached on the basis of geometric structures and excited-state energetics.

Complexes between H2 and neutral oxyacid beryllium derivatives. The role of angular strain

ABSTRACTThe complexes between 14 different Be-salts of oxyacids from groups 13 to 16 with one and two molecules of H2 have been investigated by means of the MP2/aug-cc-pVTZ ab initio molecular orbital theory method, which was found to be reliable for the treatment of these weakly bound species. The main conclusion is that these Be-salts yield rather stable complexes with dihydrogen, with binding energies one order of magnitude larger than other typical H2 complexes reported in the literature. This strong binding is shown to be due to an enhancement of the electron-deficient nature of Be when attached to an oxyacid moiety, which depends more on the type of coordination of the central atom of the oxyacid moiety. The formation of these complexes is followed by a significant lengthening of the H2 internuclear distance and a concomitant red-shift of the H–H stretching frequency, which becomes a good indicator of the strength of the interaction. The charge shifting from the bonding region of the H2 molecule to the interboundary Be···H2 region is the physical phenomenon behind the stability of these complexes. Accordingly, the most important contributor to this stability is the inductive term, followed by the electrostatic interactions. The ability of Be to bind H2 is enhanced by the angular arrangement of the O–Be–O electron-acceptor group.