Tatsusada Yoshida

Correlation Analyses on Binding Affinity of Sialic Acid Analogues with Influenza Virus Neuraminidase-1 Using ab Initio MO Calculations on Their Complex Structures

We carried out full ab initio molecular orbital calculations on complexes between neuraminidase-1 (N1-NA) in the influenza A virus and a series of eight sialic acid analogues including oseltamivir (Tamiflu) in order to quantitatively examine the binding mechanism and variation in the inhibitory potency at the atomic and electronic levels. FMO-MP2-IFIE (interfragment interaction energy at the MP2 level of ab initio fragment molecular orbital calculations) analyses quantitatively revealed (1) that the complex formation is driven by strong electrostatic interactions of charged functional groups in the analogues with ionized amino acid residues and water molecules in the active site of N1-NA, and (2) that the variation in the inhibitory potency among the eight analogues is determined by the dispersion and/or hydrophobic interaction energies of the 3-pentyl ether and charged amino moieties in oseltamivir with certain residues and water molecules in the active site of N1-NA. The current results will be useful for the development of new antiinfluenza drugs with high potency against various subtypes of wild-type and drug-resistant NAs.

Novel Quantitative Structure-Activity Studies of HIV-1 Protease Inhibitors of the Cyclic Urea Type Using Descriptors Derived from Molecular Dynamics and Molecular Orbital Calculations

Human immunodeficiency virus type 1 protease (HIV-1 PR) is an essential enzyme for the replication cycle of HIV-1. HIV-1 PR inhibitors have been extensively investigated as anti-AIDS drugs. For developments of HIV-1 PR inhibitors more promising than those utilized at the moment, the construction of reliable QSAR models that can elucidate the inhibitory mechanism as consistently as possible should be one of the most significant issues. Garg, Kurup, and their groups published comprehensive QSAR studies using past structure-activity data for HIV-1 PR inhibitors, and summarized some physicochemical structural factors of inhibitors that govern variations in the inhibitory activity for various structural types. There seem to exist much to be clarified further, especially for effects of electronic structure of inhibitors. It is also expected to incorporate structural and physicochemical information about the enzyme protein into the QSAR model. In this article, we reviewed our own QSAR study on a series of cyclic urea inhibitors with newly proposed QSAR descriptors. We performed molecular dynamics simulations of HIV-1 PR-inhibitor complexes to provide the accurate geometry to the fragment molecular orbital (FMO) calculations as well as to the estimation of an accessible surface area descriptor for inhibitors and amino acid residues. With the FMO procedure to cover full electronic feature of three-dimensional structure of protease-inhibitor complexes, we derived electronic descriptors for inhibitors and amino acid residues. The successful results are believed to provide a new insight into QSAR and understanding of binding mechanism of inhibitors with HIV-1 PR at atomic and electronic levels.

Correlation Analyses on Binding Affinity of Sialic Acid Analogues and Anti-Influenza Drugs with Human Neuraminidase Using ab Initio MO Calculations on Their Complex Structures – LERE-QSAR Analysis (IV)

We carried out full ab initio fragment molecular orbital (FMO) calculations for complexes comprising human neuraminidase-2 (hNEU2) and sialic acid analogues including anti-influenza drugs zanamivir (Relenza) and oseltamivir (Tamiflu) in order to examine the variation in the observed inhibitory activity toward hNEU2 at the atomic and electronic levels. We recently proposed the LERE (linear expression by representative energy terms)-QSAR (quantitative structure-activity relationship) procedure. LERE-QSAR analysis quantitatively revealed that the complex formation is driven by hydrogen-bonding and electrostatic interaction of hNEU2 with sialic acid analogues. The most potent inhibitory activity, that of zanamivir, is attributable to the strong electrostatic interaction of a positively charged guanidino group in zanamivir with negatively charged amino acid residues in hNEU2. After we confirmed that the variation in the observed inhibitory activity among sialic acid analogues is excellently reproducible with the LERE-QSAR equation, we examined the reason for the remarkable difference between the inhibitory potencies of oseltamivir as to hNEU2 and influenza A virus neuraminidase-1 (N1-NA). Several amino acid residues in close contact with a positively charged amino group in oseltamivir are different between hNEU2 and N1-NA. FMO-IFIE (interfragment interaction energy) analysis showed that the difference in amino acid residues causes a remarkably large difference between the overall interaction energies of oseltamivir with hNEU2 and N1-NA. The current results will be useful for the development of new anti-influenza drugs with high selectivity and without the risk of adverse side effects.

Correlation analyses on binding affinity of substituted benzenesulfonamides with carbonic anhydrase using ab initio MO calculations on their complex structures (II)

We proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as ab initio fragment molecular orbital and generalized Born/surface area ones. We applied LERE-QSAR to two datasets for the free-energy changes during complex formation between carbonic anhydrase and a series of substituted benzenesulfonamides. The first compound set (Set I) and the second one (Set II) include relatively small substituents and alkyl chains of different lengths in the benzene ring, respectively. Variation of the inhibitory activity in Set I is expressed as the combination of Hammett σ and the hydrophobic substituent constant π in classical QSAR, and variation in Set II only by π. LERE-QSAR analyses clearly revealed that effects of σ and π on the activity variations in Sets I and II are consistently explainable with the energy terms in the LERE formulation, and provide more detailed and direct information as to the binding mechanism. The proposed procedure was demonstrated to provide a quantitative basis for understanding ligand-protein interactions at the electronic and atomic levels.

Combined QM/MM (ONIOM) and QSAR Approach to the Study of Complex Formation of Matrix Metalloproteinase-9 with a Series of Biphenylsulfonamides–LERE-QSAR Analysis (V)

We previously proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as an ab initio fragment molecular orbital ones. In the present work, we applied LERE-QSAR to complex formation of matrix metalloproteinase-9 (MMP-9) with a series of substituted biphenylsulfonamides. The results shows that the overall free-energy change accompanying complex formation is due to predominantly the contribution from the electrostatic interaction with the zinc atom in the active site of MMP-9. Carbonic anhydrase (CA) belongs to the zinc-containing protease family. In contrast to the current case of MMP-9, the overall free-energy change during complex formation of CA with a series of benzenesulfonamides is due to the contributions from the solvation and dissociation free-energy changes, as previously reported. Comparison of the two sets of results indicates quantitative differences in the relative contributions of free-energy components to the overall free-energy change between the two data sets, corresponding with those in the respective classical QSAR equations. The LERE-QSAR procedure was demonstrated to quantitatively reveal differences in the binding mechanisms between the two cases involving similar but different zinc-containing proteins at the electronic and atomic levels.

Molecular dynamics and density functional studies on the metabolic selectivity of antipsychotic thioridazine by cytochrome P450 2D6: Connection with crystallographic and metabolic results

CYP2D6, a cytochrome P450 isoform, significantly contributes to the metabolism of many clinically important drugs. Thioridazine (THD) is one of the phenothiazine-type antipsychotics, which exhibit dopamine D2 antagonistic activity. THD shows characteristic metabolic profiles compared to other phenothiazine-type antipsychotics such as chlorpromazine. The sulfur atom attached to the phenothiazine ring is preferentially oxidized mainly by CYP2D6, that is, the 2-sulfoxide is a major metabolite, and interestingly this metabolite shows more potent activity against dopamine D2 receptors than THD. On the other hand, the formation of this metabolite causes many serious problems for its clinical use. Wójcikowski et al. (Drug Metab. Dispos. 2006, 34, 471) reported a kinetic study of THD formed by CYP2D6. Recently, Wang et al. (J. Biol. Chem. 2012, 287, 10834 and J. Biol. Chem. 2015, 290, 5092) revealed the crystallographic structure of THD with CYP2D6. In the current study, the binding and reaction mechanisms at the atomic and electronic levels were computationally examined based on the assumption as to whether or not the different crystallographic binding poses correspond to the different metabolites. The binding and oxidative reaction steps in the whole metabolic process were investigated using molecular dynamics and density functional theory calculations, respectively. The current study demonstrated the essential importance of the orientation of the substrate in the reaction center of CYP2D6 for the metabolic reaction.

A simple and efficient dispersion correction to the Hartree–Fock theory

One of the most challenging problems in computational chemistry and in drug discovery is the accurate prediction of the binding energy between a ligand and a protein receptor. It is well known that the binding energy calculated with the Hartree-Fock molecular orbital theory (HF) lacks the dispersion interaction energy that significantly affects the accuracy of the total binding energy of a large molecular system. We propose a simple and efficient dispersion energy correction to the HF theory (HF-Dtq). The performance of HF-Dtq was compared with those of several recently proposed dispersion corrected density functional theory methods (DFT-Ds) as to the binding energies of 68 small non-covalent complexes. The overall performance of HF-Dtq was found to be nearly equivalent to that of more sophisticated B3LYP-D3. HF-Dtq will thus be a useful and powerful method for accurately predicting the binding energy between a ligand and a protein, albeit it is a simple correction procedure based on HF.

Binding interaction of SGLT with sugar and thiosugar by the molecular dynamics simulation

The human sodium-glucose co-transporter 2 (hSGLT2) is a transporter responsible for reabsorption of glucose in the proximal convoluted tubule of the kidney. hSGLT2 inhibitors, including luseogliflozin, have been developed as drugs for type 2 diabetes mellitus. Only luseogliflozin contains a thiosugar ring in its chemical structure, while other hSGLT2 inhibitors contain glucose rings. Consequently, we focused on the binding interactions of hSGLT2 with sugars and thiosugars. We first revealed that the binding affinities of thiosugars are stronger than those of sugars through molecular dynamics simulations of Vibrio parahaemolyticus, sodium-galactose co-transporter, and human hSGLT2. We then demonstrated that Na(+) dissociates from the protein to the cytoplasmic solution more slowly in the thiosugar system than in the sugar system. These differences between sugars and thiosugars are discussed on the basis of the different binding modes due to the atom at the 5-position of the sugar and thiosugar rings. Finally, as a result of Na(+) dissociation, we suggest that the dissociation of thiosugars is slower than that of sugars.

Connecting Classical QSAR and LERE Analyses Using Modern Molecular Calculations, LERE-QSAR (VI): Hydrolysis of Substituted Hippuric Acid Phenyl Esters by Trypsin

The reaction mechanism of trypsin was studied by applying DFT and ab initio molecular orbital (MO) calculations to complexes of trypsin with a congeneric series of eight para‐substituted hippuric acid phenyl esters, for which a previous quantitative structureactivity relationship (QSAR) study revealed nice linearity of Hammett substitution constant σ− with logarithmic values of the MichaelisMenten and catalytic rate constants. Based on the LERE procedure, we performed QSAR analyses on each elementary reaction step during the acylation process. The present calculations showed that the rate‐determining step during the acylation process is the transition state (TS) between the enzymesubstrate complex (ES) and tetrahedral intermediate (TET), and that the proton transfer occurs from Ser195 to His57, not between His57 and Asp102. The LERE‐QSAR analysis statistically suggested that the variation of overall free‐energy changes leading to formation of TS is governed mostly by that of activation energies required to form TS from ES. In spite of a very limited number of congeneric ligands in the current work, it is critically essential to clarify and verify physicochemical meanings of a typical QSAR/Chemoinformatics parameter, Hammett σ− based on quantum chemical calculations on the proteinligand kinetics; how Hammett σ− behaves in terms of proteinligand interaction energies.

Expression and molecular dynamics studies on effect of amino acid substitutions at Arg344 in human cathepsin A on the protein local conformation

Human lysosomal protective protein/cathepsin A (CathA) is a multifunctional protein that exhibits not only protective functions as to lysosomal glycosidases, i.e., neuraminidase 1 (NEU1) and beta-galactosidase (GLB), but also its own serine carboxypeptidase activity, and exhibits conserved structural similarity to yeast and wheat homologs (CPY and CPW). Our previous study revealed that the R344 (Arg344) residue in CathA could contribute to the binding and recognition of the serine peptidase inhibitor chymostatin. We examined here the effects of substitution of R344 with other amino acids, including A, D, E, G, I, K, M, N, P, Q, S, and V, denoted as R344X, including the wild-type CathA, on expression of CathA activity and intracellular processing. Among the mutant gene products, the 54-kDa precursor/zymogen with the R344D substitution was not processed to the 32/20-kDa mature form with CathA activity in a fibroblastic cell line derived from a galactosialidosis patient. Molecular dynamics (MD) simulations on the total twelve R344X mutants and the wild-type revealed that only R344D takes on a significantly different conformation of S293-D295 in the excision peptide (M285-R298) compared to the other R344X mutants; the side chains of S293 and D295 in R344D are exposed on the molecular surface, although those in the other twelve R344X mutants are buried inside the protein. The results of the current work strongly suggest that the distinct conformational change of the S293-D295 region in the R344D protein causes the processing defect of the 54-kDa precursor of the R344D mutant gene product in cultured cells.