Structure-based virtual screening to identify inhibitors against Staphylococcus aureus MurD enzyme

Structural Chemistry - The MurD enzyme of Staphylococcus aureus is an attractive drug target as it is essential and ubiquitous in bacteria but absent in mammalian cells. In the present study, we...     

Improved electron-conformational method of pharmacophore identification and bioactivity prediction. Application to angiotensin converting enzyme inhibitors

The electron-conformational (EC) method of pharmacophore (Pha) identification and bioactivity prediction, suggested earlier, is given here two major improvements. First, an atomic index of orbital and charge controlled interaction is introduced to better represent the ligand (substrate) in its interaction with the bioreceptor. Second, the multiconformational problem is considered in view of ligand-receptor binding states, resulting in essential simplification of the expression of bioactivity. The details of the improved EC method are demonstrated in application to the problem of angiotensin converting enzyme (ACE) inhibitors. The Pha of the latter is identified by separation of the heavily populated conformations of the chosen 51 compounds (the training set), calculation of the electronic structure, construction of their EC matrixes of congruity, and processing of the latter in comparison with the activities to reveal a common submatrix of all the active only compounds that describes the Pha. The latter contains three oxygen atoms plus a fourth atom X = S, N, O at certain interatomic distances and with restricted electronic parameters (within assumed tolerances), the position of the atom X being more changeable from one active compound to another. For quantitative prediction of the bioactivity, an expression is deduced which takes into account the duly parametrized influence of auxiliary groups (AG) which, being positioned outside the Pha, either diminish the activity (antipharmacophore shielding) or enhance it. It is shown that in case of many conformations of the same compound only one of them, that of the lowest energy which has the Pha, should be parametrized. The 15 parameters chosen to represent the AG in case of ACE inhibitors are weighted by variational (adjustable) coefficients which are determined from a regression treatment of the calculated versus known activities in the training set. Then the formulas with known coefficients are used to validate the method by calculating the bioactivity of other compounds not used in the training set. The prediction of the activity proved to be more than 90% (within experimental error and available compounds) qualitatively (yes, no) and about 60%-70% quantitatively.     

Structure-based virtual screening for novel potential selective inhibitors of class IIa histone deacetylases for cancer treatment

The fundamental cause of human cancer is strongly influenced by down- or up-regulations of epigenetic factors. Upregulated histone deacetylases (HDAC) have been shown to be effectively neutralized by the action of HDACs inhibitors (HDACi). However, cytotoxicity has been reported in normal cells because of non-specificity of several available HDACis that are in clinical use or at different phases of clinical trials. Because of the high amino acid sequence and structural similarity among HDAC enzymes, it is believed to be a challenging task to obtain isoform-selectivity. The essential aim of the present research work was to identify isoform-selective inhibitors against class IIa HDACs via structure-based drug design. Based on the highest binding affinity and isoform-selectivity, the top-ranked inhibitors were in silico tested for their absorption, distribution, metabolism, elimination, and toxicity (ADMET) properties, which were classified as drug-like compounds. Later, molecular dynamics simulation (MD) was carried out for all compound-protein complexes to evaluate the structural stability and the biding mode of the inhibitors, which showed high stability throughout the 100 ns simulation. Free binding energy predictions by MM-PBSA method showed the high binding affinity of the identified compounds toward their respective targets. Hence, these inhibitors could be used as drug candidates or as lead compounds for more in silico or in vitro optimization to design safe isoform-selective HDACs inhibitors.     

Antioxidant activity of peptide-based angiotensin converting enzyme inhibitors

Angiotensin converting enzyme (ACE) inhibitors are important for the treatment of hypertension as they can decrease the formation of vasopressor hormone angiotensin II (Ang II) and elevate the levels of vasodilating hormone bradykinin. It is observed that bradykinin contains a Ser-Pro-Phe motif near the site of hydrolysis. The selenium analogues of captopril represent a novel class of ACE inhibitors as they also exhibit significant antioxidant activity. In this study, several di- and tripeptides containing selenocysteine and cysteine residues at the N-terminal were synthesized. Hydrolysis of angiotensin I (Ang I) to Ang II by ACE was studied in the presence of these peptides. It is observed that the introduction of L-Phe to Sec-Pro and Cys-Pro peptides significantly increases the ACE inhibitory activity. On the other hand, the introduction of L-Val or L-Ala decreases the inhibitory potency of the parent compounds. The presence of an L-Pro moiety in captopril analogues appears to be important for ACE inhibition as the replacement of L-Pro by L-piperidine 2-carboxylic acid decreases the ACE inhibition. The synthetic peptides were also tested for their ability to scavenge peroxynitrite (PN) and to exhibit glutathione peroxidase (GPx)-like activity. All the selenium-containing peptides exhibited good PN-scavenging and GPx activities.     

Synthesis, characterization and antioxidant activity of angiotensin converting enzyme inhibitors

Angiotensin converting enzyme (ACE) catalyzes the conversion of angiotensin I (Ang I) to angiotensin II (Ang II). ACE also cleaves the terminal dipeptide of vasodilating hormone bradykinin (a nonapeptide) to inactivate this hormone. Therefore, inhibition of ACE is generally used as one of the methods for the treatment of hypertension. 'Oxidative stress' is another disease state caused by an imbalance in the production of oxidants and antioxidants. A number of studies suggest that hypertension and oxidative stress are interdependent. Therefore, ACE inhibitors having antioxidant property are considered beneficial for the treatment of hypertension. As selenium compounds are known to exhibit better antioxidant behavior than their sulfur analogues, we have synthesized a number of selenium analogues of captopril, an ACE inhibitor used as an antihypertensive drug. The selenium analogues of captopril not only inhibit ACE activity but also effectively scavenge peroxynitrite, a strong oxidant found in vivo.     

Structure-based discovery of novel non-nucleosidic DNA alkyltransferase inhibitors: virtual screening and in vitro and in vivo activities

The human DNA-repair O (6)-alkylguanine DNA alkyltransferase (MGMT or hAGT) protein protects DNA from environmental alkylating agents and also plays an important role in tumor resistance to chemotherapy treatment. Available inhibitors, based on pseudosubstrate analogs, have been shown to induce substantial bone marrow toxicity in vivo. These deficiencies and the important role of MGMT as a resistance mechanism in the treatment of some tumors with dismal prognosis like glioblastoma multiforme, the most common and lethal primary malignant brain tumor, are increasing the attention toward the development of improved MGMT inhibitors. Here, we report the identification for the first time of novel non-nucleosidic MGMT inhibitors by using docking and virtual screening techniques. The discovered compounds are shown to be active in both in vitro and in vivo cellular assays, with activities in the low to medium micromolar range. The chemical structures of these new compounds can be classified into two families according to their chemical architecture. The first family corresponds to quinolinone derivatives, while the second is formed by alkylphenyl-triazolo-pyrimidine derivatives. The predicted inhibitor protein interactions suggest that the inhibitor binding mode mimics the complex between the excised, flipped out damaged base and MGMT. This study opens the door to the development of a new generation of MGMT inhibitors.     

基于结构设计合成新型喹啉类分泌型PLA2抑制剂

磷脂酶A2(PLA2)在很多人类疾病的病理研究中起重要作用,是药物化学研究所的热点之一。因此,发展新型的PLA2抑制剂对生物有机研究和临床应用均有重要意义。我们设计合成分泌型PLA2的非底物类似物以寻找新型PLA2抑制剂。本文基于结构设计并合成了喹啉－4－乙酰胺作为PLA2抑制剂,目标化合物结构经1^HNMR,IR,MS和元素分析确认,初步活性检测显示该类化合物具有较好体外活性及动物活性。     

Structure-based virtual screening of dipeptidyl peptidase 4 inhibitors and their in vitro analysis

Type 2 diabetes mellitus (T2DM) is one of the most widely prevalent metabolic disorders with no cure to date thus remains the most challenging task in the current drug discovery. Therefore, the only strategy to control diabetes prevalence is to develop novel efficacious therapeutics. Dipeptidyl Peptidase 4 (DPP-4) inhibitors are currently used as anti-diabetic drugs for the inhibition of incretins. This study aims to construct the chemical feature based on pharmacophore models for dipeptidyl peptidase IV. The structure-based pharmacophore modeling has been employed to evaluate new inhibitors of DPP-4. A four-featured pharmacophore model was developed from crystal structure of DPP-4 enzyme with 4-(2-aminoethyl) benzenesulfonyl fluoride in its active site via pharmacophore constructing tool of Molecular Operating Environment (MOE) consisting F1 Hyd (hydrophobic region), F2 Hyd|Cat|Don (hydrophobic cationic and donor region), F3 Acc (acceptor region) and F4 Hyd (hydrophobic region). The generated pharmacophore model was used for virtual screening of in-house compound library (the available compounds which were used for initial screening to get the few compounds for the current studies). The resultant selected compounds, after virtual screening were further validated using in vitro assay. Furthermore, structure-activity relationship was carried out for the compounds possessing significant inhibition potential after docking studies. The binding free energy of analogs was evaluated via molecular mechanics generalized Born surface area (MM-GBSA) and Poisson-Boltzmann surface area (MM-PBSA) methods using AMBER 16 as a molecular dynamics (MD) simulation package. Based on potential findings, we report that selected candidates are more likely to be used as DPP-4 inhibitors or as starting leads for the development of novel and potent DPP-4 inhibitors.     

Structure-based 3D QSAR and design of novel acetylcholinesterase inhibitors

The paper describes the construction, validation and application of a structure-based 3D QSAR model of novel acetylcholinesterase (AChE) inhibitors. Initial use was made of four X-ray structures of AChE complexed with small, non-specific inhibitors to create a model of the binding of recently developed aminopyridazine derivatives. Combined automated and manual docking methods were applied to dock the co-crystallized inhibitors into the binding pocket. Validation of the modelling process was achieved by comparing the predicted enzyme-bound conformation with the known conformation in the X-ray structure. The successful prediction of the binding conformation of the known inhibitors gave confidence that we could use our model to evaluate the binding conformation of the aminopyridazine compounds. The alignment of 42 aminopyridazine compounds derived by the docking procedure was taken as the basis for a 3D QSAR analysis applying the GRID/GOLPE method. A model of high quality was obtained using the GRID water probe, as confirmed by the cross-validation method (q² _LOO=0.937, q² _{L50% O}=0.910). The validated model, together with the information obtained from the calculated AChE-inhibitor complexes, were considered for the design of novel compounds. Seven designed inhibitors which were synthesized and tested were shown to be highly active. After performing our modelling study the X-ray structure of AChE complexed with donepezil, an inhibitor structurally related to the developed aminopyirdazines, has been made available. The good agreement found between the predicted binding conformation of the aminopyridazines and the one observed for donepezil in the crystal structure further supports our developed model.     

Structure-based rational quest for potential novel inhibitors of human HMG-CoA reductase by combining CoMFA 3D QSAR modeling and virtual screening

3-Hydroxy-3-methylglutaryl-coenzyme A reductase (HMGR) catalyzes the formation of mevalonate. In many classes of organisms, this is the committed step leading to the synthesis of essential compounds, such as cholesterol. However, a high level of cholesterol is an important risk factor for coronary heart disease, for which an effective clinical treatment is to block HMGR using inhibitors like statins. Recently the structures of catalytic portion of human HMGR complexed with six different statins have been determined by a delicate crystallography study (Istvan and Deisenhofer Science 2001, 292, 1160-1164), which established a solid basis of structure and mechanism for the rational design, optimization, and development of even better HMGR inhibitors. In this study, three-dimensional quantitative structure-activity relationship (3D QSAR) with comparative molecular field analysis (CoMFA) was performed on a training set of up to 35 statins and statin-like compounds. Predictive models were established by using two different ways: (1) Models-fit, obtained by SYBYL conventional fit-atom molecular alignment rule, has cross-validated coefficients (q2) up to 0.652 and regression coefficients (r2) up to 0.977. (2) Models-dock, obtained by FlexE by docking compounds into the HMGR active site, has cross-validated coefficients (q2) up to 0.731 and regression coefficients (r2) up to 0.947. These models were further validated by an external testing set of 12 statins and statin-like compounds. Integrated with CoMFA 3D QSAR predictive models, molecular surface property (electrostatic and steric) mapping and structure-based (both ligand and receptor) virtual screening have been employed to explore potential novel hits for the HMGR inhibitors. A representative set of eight new compounds of non-statin-like structures but with high pIC(50) values were sorted out in the present study.     

新型血管紧张素转化酶抑制剂: 氮杂三肽类似物的合成及结构活性关系

设计与合成一系列氮杂三肽类似物, 测定了氮杂三肽类似物对血管紧张素转化酶的体外抑制活性, 探讨了它们的结构与抑制活性之间的关系, 结果表明氮杂丙氨酸是丙氨酸很好的替代物。     

Structure-based design of potent non-peptide MDM2 inhibitors

A successful structure-based design of a class of non-peptide small-molecule MDM2 inhibitors targeting the p53-MDM2 protein-protein interaction is reported. The most potent compound 1d binds to MDM2 protein with a Ki value of 86 nM and is 18 times more potent than a natural p53 peptide (residues 16-27). Compound 1d is potent in inhibition of cell growth in LNCaP prostate cancer cells with wild-type p53 and shows only a weak activity in PC-3 prostate cancer cells with a deleted p53. Importantly, 1d has a minimal toxicity to normal prostate epithelial cells. Our studies provide a convincing example that structure-based strategy can be employed to design highly potent, non-peptide, cell-permeable, small-molecule inhibitors to target protein-protein interaction, which remains a very challenging area in chemical biology and drug design.     

Discovery of potential M2 channel inhibitors based on the amantadine scaffold via virtual screening and pharmacophore modeling

The M2 channel protein on the influenza A virus membrane has become the main target of the anti-flu drugs amantadine and rimantadine. The structure of the M2 channel proteins of the H3N2 (PDB code 2RLF) and 2009-H1N1 (Genbank accession number GQ385383) viruses may help researchers to solve the drug-resistant problem of these two adamantane-based drugs and develop more powerful new drugs against influenza A virus. In the present study, we searched for new M2 channel inhibitors through a combination of different computational methodologies, including virtual screening with docking and pharmacophore modeling. Virtual screening was performed to calculate the free energies of binding between receptor M2 channel proteins and 200 new designed ligands. After that, pharmacophore analysis was used to identify the important M2 protein-inhibitor interactions and common features of top binding compounds with M2 channel proteins. Finally, the two most potential compounds were determined as novel leads to inhibit M2 channel proteins in both H3N2 and 2009-H1N1 influenza A virus.     

Studies on angiotensin converting enzyme inhibitors. VI. Synthesis and angiotensin converting enzyme inhibitory activities of the dicarboxylic acid derivative of imidapril and its diastereoisomers.

All possible diastereoisomers of the dicarboxylic acid (10a), the biologically active form of imidapril (1), were synthesized, and their inhibitory activity against angiotensin converting enzyme (ACE) was examined. The in vitro ACE inhibitory activity of these compounds greatly depended on the configurations of the three asymmetric carbons in each molecule. The (S,S,S) isomer (10a) showed much more potent activity than the others.     

亲和色谱法筛选中药中血管紧张素转化酶抑制剂

以壳聚糖微球为载体、戊二醛（glutaraldehyde,GA）为交联剂对血管紧张素转化酶（Angiotensin converting enzyme,ACE）进行固定化．用固化的ACE作为亲和介质,利用血管紧张素转化酶抑制剂（Angiotensin convertingenzym einhibitor,ACEI）与ACE之间的亲和作用,结合高效液相色谱对亲和前后的体系进行检测,比较两者各组分色谱峰的差异,以此实现快速筛选复杂体系中的ACE抑制剂．应用赖诺普利（Lisinopril）、九肽抑制剂、依那普利（Enalapril）、培哚普利（Perindopril）、卡托普利（Captopril）等已上市的ACEI对方法进行验证,反映方法具有高度选择性．将方法应用于中药地龙及山楂筛选,发现共有5个组分与ACE有亲和作用,并且都能抑制ACE酶活性,它们对酶活性抑制的IC50值在0．45～4．62μg／mL范围．通过对亲和方法重现性考察,6次测定的相对标准偏差小于1％,说明方法可靠．提出的亲和色谱．色谱指纹差异法非常适合于从中药及天然产物等复杂混合物库中快速筛选靶点活性物质．     

Refined pharmacophore features for virtual screening of human thromboxane A2 receptor antagonists

For a long time, the structural basis of TXA2 receptor is limited due to the lack of crystal structure information, till the release of the crystal structure of TXA2 receptor, which deepens our understanding about ligand recognition and selectivity mechanisms of this physiologically important receptor. In this research, we report the successful implementation in the discovery of an optimal pharmacophore model of human TXA2 receptor antagonists through virtual screening. Structure-based pharmacophore models were generated based on two crystal structures of human TXA2 receptor (PDB entry 6IIU and 6IIV). Docking simulation revealed interaction modes of the virtual screening hits against TXA2 receptor, which was validated through molecular dynamics simulation and binding free energy calculation. ADMET properties were also analyzed to evaluate the toxicity and physio-chemical characteristics of the hits. The research would provide valuable insight into the binding mechanisms of TXA2 receptor antagonists and thus be helpful for designing novel antagonists.     

Discovery of novel human acrosin inhibitors by virtual screening

Human acrosin is an attractive target for the discovery of male contraceptive drugs. For the first time, structure-based drug design was applied to discover structurally diverse human acrosin inhibitors. A parallel virtual screening strategy in combination with pharmacophore-based and docking-based techniques was used to screen the SPECS database. From 16 compounds selected by virtual screening, a total of 10 compounds were found to be human acrosin inhibitors. Compound 2 was found to be the most potent hit (IC₅₀ = 14 μM) and its binding mode was investigated by molecular dynamics simulations. The hit interacted with human acrosin mainly through hydrophobic and hydrogen-bonding interactions, which provided a good starting structure for further optimization studies.