Molecular dynamics simulation study and hybrid pharmacophore model development in human LTA4H inhibitor design

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

In silico investigations of some Cyperus rotundus compounds as potential anti-inflammatory inhibitors of 5-LO and LTA4H enzymes

The present study aimed to experimentally identify the essential oil of Algerian Cyperus rotundus L. and to model the interaction of some known anti-inflammatory molecules with two key enzymes involved in inflammation, 5-Lypoxygenase (5-LO) and leukotriene A4 hydrolase (LTA4H). Gas chromatography/gas chromatography-mass spectrometry (GC/GC-MS) revealed that 92.7% of the essential oil contains 35 compounds, including oxygenated sesquiterpenes (44.2%), oxygenated monoterpenes (30.2%), monoterpene hydrocarbons (11.8%) and sesquiterpene hydrocarbons (6.5%). The major identified oxygenated terpenes are humulene oxide II, caryophyllene oxide, khusinol, agarospirol, spathulinol and trans-pinocarveol Myrtenol and α-terpineol are known to exhibit anti-inflammatory activities. Several complexes obtained after docking the natural terpenes with 5-LO and LTA4H have shown strong hydrogen bonding interactions. The best docking energies were found with α-terpineol, Myrtenol and khusinol. The interaction between the natural products and amino-acid residues HIS367, ILE673 and GLN363 appears to be critical for 5-LO inhibition, while the interaction with residues GLU271, HIS295, TYR383, TYR378, GLU318, GLU296 and ASP375 is critical for LTA4H inhibition. Molecular dynamics (MD) trajectories of the selected docked complexes showed stable backbone root mean square deviation (RMSD), supporting the stability of the natural product-enzyme interaction.Communicated by Ramaswamy H. Sarma.

The Pharmacophore Concept and Its Applications in Computer-Aided Drug Design

Pharmacophore-based techniques currently are an integral part of many computer-aided drug design workflows and have been successfully and extensively applied for tasks such as virtual screening, de novo design, and lead optimization. Pharmacophore models can be derived both in a receptor-based and in a ligand-based manner, and provide an abstract description of essential non-bonded interactions that typically occur between small-molecule ligands and macromolecular targets. Due to their simplistic and abstract nature, pharmacophores are both perfectly suited for efficient computer processing and easy to comprehend by life and physical scientists. As a consequence, they have also proven to be a valuable tool for communicating between computational and medicinal chemists.This chapter aims to provide a short overview of the pharmacophore concept and its applications in modern computer-aided drug design. The chapter is divided into three distinct parts. The first section contains a brief introduction to the pharmacophore concept. The second section provides a description of the most common nonbonded interaction types and their representation as pharmacophoric features. Furthermore, it gives an overview of the various methods for pharmacophore generation and important pharmacophore-based techniques in drug design. This part concludes with examples for recent pharmacophore concept-related research and development. The last section is dedicated to a review of research in the field of natural product chemistry as carried out by employing pharmacophore-based drug design methods.

QM/MM Molecular Dynamics Investigations of the Substrate Binding of Leucotriene A4 Hydrolase: Implication for the Catalytic Mechanism

LTA4H is a monozinc bifunctional enzyme which exhibits both aminopeptidase and epoxide hydrolase activities. Its dual functions in anti- and pro-inflammatory roles have attracted wide attention to the inhibitor design. In this work, we tried to construct Michaelis complexes of LTA4H with both a native peptide substrate and LTA4 molecule using combined quantum mechanics and molecular mechanics molecular dynamics simulations. First of all, the zinc ion is coordinated by H295, H299, and E318. For its aminopeptidase activity, similar to conventional peptidases, the fourth ligand to the zinc ion is suggested to be an active site water, which is further hydrogen bonded with a downstream glutamic acid, E296. For the epoxide hydrolase activity, the fourth ligand to the zinc ion is found to be an epoxy oxygen atom. The potential of mean force calculation indicates about an 8.5 kcal/mol activation barrier height for the ring-opening reaction, which will generate a metastable carbenium intermediate. Subsequent frontier molecular orbital analyses suggest that the next step would be the nucleophilic attacking reaction at the C12 atom by a water molecule activated by D375. Our simulations also analyzed functions of several important residues like R563, K565, E271, Y383, and Y378 in the binding of peptide and LTA4.

Drug Repurposing of Histone Deacetylase Inhibitors That Alleviate Neutrophilic Inflammation in Acute Lung Injury and Idiopathic Pulmonary Fibrosis via Inhibiting Leukotriene A4 Hydrolase and Blocking LTB4 Biosynthesis

Acute lung injury (ALI) and idiopathic pulmonary fibrosis (IPF) are both serious public health problems with high incidence and mortality rate in adults, and with few drugs available for the efficient treatment in clinic. In this study, we identified that two known histone deacetylase (HDAC) inhibitors, suberanilohydroxamic acid (SAHA, 1) and its analogue 4-(dimethylamino)-N-[7-(hydroxyamino)-7-oxoheptyl]benzamide (2), are effective inhibitors of Leukotriene A4 hydrolase (LTA4H), a key enzyme in the biosynthesis of leukotriene B4 (LTB4), across a panel of 18 HDAC inhibitors, using enzymatic assay, thermofluor assay, and X-ray crystallographic investigation. Importantly, both 1 and 2 markedly diminish early neutrophilic inflammation in mouse models of ALI and IPF under a clinical safety dose. Detailed mechanisms of down-regulation of proinflammatory cytokines by 1 or 2 were determined in vivo. Collectively, 1 and 2 would provide promising agents with well-known clinical safety for potential treatment in patients with ALI and IPF via pharmacologically inhibiting LAT4H and blocking LTB4 biosynthesis.

Common Hits Approach: Combining Pharmacophore Modeling and Molecular Dynamics Simulations

We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hit-list. The score for a particular molecule is calculated based on the number of representative pharmacophore models which classified it as active. Hence, the method is called common hits approach (CHA). The steps between the MD simulation and the final hit-list are performed automatically and without user interaction. We test the performance of CHA for virtual screening using screening databases with active and inactive compounds for 40 protein-ligand systems. The results of the CHA are compared to the (i) median screening performance of all representative pharmacophore models of protein-ligand systems, as well as to the virtual screening performance of (ii) a random classifier, (iii) the pharmacophore model derived from the experimental structure in the PDB, and (iv) the representative pharmacophore model appearing most frequently during the MD simulation. For the 34 (out of 40) protein-ligand complexes, for which at least one of the approaches was able to perform better than a random classifier, the highest enrichment was achieved using CHA in 68% of the cases, compared to 12% for the PDB pharmacophore model and 20% for the representative pharmacophore model appearing most frequently. The availabilithy of diverse sets of different pharmacophore models is utilized to analyze some additional questions of interest in 3D pharmacophore-based virtual screening.

Relationship between human LTA4H polymorphisms and extra-pulmonary tuberculosis in an ethnic Han Chinese population in Eastern China

Two single nucleotide polymorphisms in Leukotriene A4 hydrolase (LTA4H) gene were reported to be associated with protection from pulmonary tuberculosis in Vietnamese population. But these associations were not found in the Russians. To investigate the association of LTA4H polymorphisms with tuberculosis in a Han Chinese population in Eastern China, we genotyped 5 SNPs of LTA4H gene in 743 of pulmonary tuberculosis patients, 372 of extra-pulmonary tuberculosis patients and 888 of healthy controls individuals. The CC and TT homozygotes of rs1978331 and rs2540474 were identified to have higher rates (P < 0.01) and be risk factors in the patients with extra-pulmonary tuberculosis (OR = 1.412; 95% CI = 1.104-1.804 and(OR = 1.380; 95% CI = 1.080-1.764). However, no significant association was found between any of the SNPs and pulmonary tuberculosis. In the extra-pulmonary tuberculosis subgroups. LTA4H gene were significantly associated with tuberculous meningitis, lymph node tuberculosis, bone tuberculosis and other extra-pulmonary tuberculosis except for pleural tuberculosis. The present findings suggest that polymorphisms in the LTA4H gene may affect susceptibility to extra-pulmonary tuberculosis and change the risk of developing the disease in the Han nationality in the East China.

Structural origins for the loss of catalytic activities of bifunctional human LTA4H revealed through molecular dynamics simulations

Human leukotriene A4 hydrolase (hLTA4H), which is the final and rate-limiting enzyme of arachidonic acid pathway, converts the unstable epoxide LTA4 to a proinflammatory lipid mediator LTB4 through its hydrolase function. The LTA4H is a bi-functional enzyme that also exhibits aminopeptidase activity with a preference over arginyl tripeptides. Various mutations including E271Q, R563A, and K565A have completely or partially abolished both the functions of this enzyme. The crystal structures with these mutations have not shown any structural changes to address the loss of functions. Molecular dynamics simulations of LTA4 and tripeptide complex structures with functional mutations were performed to investigate the structural and conformation changes that scripts the observed differences in catalytic functions. The observed protein-ligand hydrogen bonds and distances between the important catalytic components have correlated well with the experimental results. This study also confirms based on the structural observation that E271 is very important for both the functions as it holds the catalytic metal ion at its location for the catalysis and it also acts as N-terminal recognition residue during peptide binding. The comparison of binding modes of substrates revealed the structural changes explaining the importance of R563 and K565 residues and the required alignment of substrate at the active site. The results of this study provide valuable information to be utilized in designing potent hLTA4H inhibitors as anti-inflammatory agents.