Predicting functional residues of the Solanum lycopersicum aspartic protease inhibitor (SLAPI) by combining sequence and structural analysis with molecular docking

The Solanum lycopersicum aspartic protease inhibitor (SLAPI), which belongs to the STI-Kunitz family, is an effective inhibitor of the aspartic proteases human cathepsin D and Saccharomyces proteinase A. However, in contrast with the large number of studies on the inhibition mechanism of the serine proteases by the STI-Kunitz inhibitors, the structural aspects of the inhibition mechanism of aspartic proteases from this family of inhibitors are poorly understood. In the present study, we have combined sequence and structural analysis methods with protein-protein docking to gain a better understanding of the SLAPI inhibition mechanism of the proteinase A. The results suggest that: i) SLAPI loop L9 may be involved in the inhibitor interaction with the proteinase A′s active site, and ii) the residues I144, V148, L149, P151, F152 and R154 are implicated in the difference in the potency shown previously by SLAPI and another STI-Kunitz inhibitor isolated from Solanum tuberosum to inhibit proteinase A. These results will be useful in the design of site directed mutagenesis experiments to understand more thoroughly the aspartic protease inhibition mechanism of SLAPI and other related STI-Kunitz inhibitors.     

Molecular docking identifies the binding of 3-chloropyridine moieties specifically to the S1 pocket of SARS-CoV Mpro

The 3C-like main proteinase of the severe acute respiratory syndrome (SARS) coronavirus, SARS-CoV M(pro), is widely considered to be a major drug target for the development of anti-SARS treatment. Based on the chemical structure of a lead compound from a previous screening, we have designed and synthesized a number of non-peptidyl inhibitors, some of which have shown significantly improved inhibitory activity against SARS-CoV M(pro) with IC(50) values of approximately 60 nM. In the absence of SARS-CoV M(pro) crystal structures in complex with these synthetic inhibitors, molecular docking tools have been employed to study possible interactions between these inhibitors and SARS-CoV M(pro). The docking results suggest two major modes for the initial binding of these inhibitors to the active site of SARS-CoV M(pro). They also establish a structural basis for the 'core design' of these inhibitors by showing that the 3-chloropyridine functions common to all of the present inhibitors tend to cluster in the S1 specificity pocket. In addition, intrinsic flexibility in the S4 pocket allows for the accommodation of bulky groups such as benzene rings, suggesting that this structural plasticity can be further exploited for optimizing inhibitor-enzyme interactions that should promote a tighter binding mode. Most importantly, our results provide the structural basis for rational design of wide-spectrum antiviral drugs targeting the chymotrypsin-like cysteine proteinases from coronaviruses and picornaviruses.     

Understanding microscopic binding of macrophage migration inhibitory factor with phenolic hydrazones by molecular docking, molecular dynamics simulations and free energy calculations

Macrophage migration inhibitory factor (MIF), an immunoregulatory protein, is a potential target for a number of inflammatory diseases. In the current work, the interactions between MIF and a series of phenolic hydrazones were studied by molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decomposition analysis to determine the structural requirement for achieving favorable biological activity of phenolic hydrazones. First, molecular docking was used to predict the binding modes of inhibitors in the binding site of MIF. The good correlation between the predicted docking scores and the experimental activities shows that the binding conformations of the inhibitors in the active site of MIF are well predicted. Moreover, our results suggest that the flexibility of MIF is essential in ligand binding process. Then, MD simulations and MM/GBSA free energy calculations were employed to determine the dynamic binding process and compare the binding modes of the inhibitors with different activities. The predicted binding free energies given by MM/GBSA are not well correlated with the experimental activities for the two subsets of the inhibitors; however, for each subset, a good correlation between the predicted binding free energies and the experimental activities is achieved. The MM/GBSA free energy decomposition analysis highlights the importance of hydrophobic residues for the MIF binding of the studied inhibitors. Based on the essential factors for MIF-inhibitor interactions derived from the theoretical predictions, some derivatives were designed and the higher inhibitory activities of several candidates were confirmed by molecular docking studies. The structural insights obtained from our study are useful for designing potent inhibitors of MIF.     

Identifying inhibitors of the SARS coronavirus proteinase

The Severe Acute Respiratory Syndrome (SARS) is a serious respiratory illness that has recently been reported in parts of Asia and Canada. In this study, we use molecular dynamics (MD) simulations and docking techniques to screen 29 approved and experimental drugs against the theoretical model of the SARS CoV proteinase as well as the experimental structure of the transmissible gastroenteritis virus (TGEV) proteinase. Our predictions indicate that existing HIV-1 protease inhibitors, L-700,417 for instance, have high binding affinities and may provide good starting points for designing SARS CoV proteinase inhibitors.     

Exploring the conformational and binding properties of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 through docking and molecular dynamics simulations

B‐cell lymphoma (Bcl‐2) is commonly associated with the progression and preservation of cancer and certain lymphomas; therefore, it is considered as a biological target against cancer. Nevertheless, evidence of all its structural binding sites has been hidden because of the lack of a complete Bcl‐2 model, given the presence of a flexible loop domain (FLD), which is responsible for its complex behavior. FLD region has been implicated in phosphorylation, homotrimerization, and heterodimerization associated with Bcl‐2 antiapoptotic function. In this contribution, homology modeling, molecular dynamics (MD) simulations in the microsecond (µs) time‐scale and docking calculations were combined to explore the conformational complexity of unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 systems. Conformational ensembles generated through MD simulations allowed for identifying the most populated unphosphorylated/phosphorylated monomeric conformations, which were used as starting models to obtain trimeric complexes through protein–protein docking calculations, also submitted to µs MD simulations. Principal component analysis showed that FLD represents the main contributor to total Bcl‐2 mobility, and is affected by phosphorylation and oligomerization. Subsequently, based on the most representative unphosphorylated/phosphorylated monomeric and trimeric Bcl‐2 conformations, docking studies were initiated to identify the ligand binding site of several known Bcl‐2 inhibitors to explain their influence in homo‐complex formation and phosphorylation. Docking studies showed that the different conformational states experienced by FLD, such as phosphorylation and oligomerization, play an essential role in the ability to make homo and hetero‐complexes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 393–413, 2016.     

Structural insights into the binding of MMP9 inhibitors

Matrix Metalloproteinase are family of enzymes responsible for degradation of extracellular matrix. MMP9 (gelatinase B) is one of the common matrix metalloproteinase that is associated with tissue destruction in a number of disease states such as rheumatoid arthiritis, fibrotic lung disease, dilated cardiomyopathy, as well as cancer invasion and metastasis. Recent study demonstrates that increased expression of MMP9 results in augmentation of myopathy with increased inflammation and fibernecrosis. Previous studies do not provide any conclusive information related to structural specificity of MMP9 inhibitors towards its active site, but with the availability of experimental structures it is now possible to study the structural specificity of MMP9 inhibitors. In light of availability of this information, we have applied docking and molecular dynamics approach to study the binding of inhibitors to the active site of MMP9. Three categories of inhibitor consisting of sulfonamide hydroxamate, thioester, and carboxylic moieties as zinc binding groups (ZBG) were chosen in the present study. Our docking results demonstrate that thioester based zinc binding group gives favourable docking scores as compared to other two groups. Molecular Dynamics simulations further reveal that tight binding conformation for thioester group has high specificity for MMP9 active site. Our study provides valuable insights on inhibitor specificity of MMP9 which provides valuable hints for future design of potent inhibitors and drugs.     

Direct interaction of ONO-5046 with human neutrophil elastase through ¹H NMR and molecular docking

Human neutrophil elastase (HNE) has been implicated as a major contributor in the pathogenesis of diseases, such as pulmonary emphysema, acute lung injury (ALI), acute respiratory distress syndrome (ARDS), and other inflammatory diseases. Therefore, searching for appropriate and potential human neutrophil elastase inhibitors (HNEI) that would restore the balance between the free enzyme and the endogenous inhibitors would be of therapeutic interest. ONO-5046 is the first specific HNEI to improve respiratory function and protect lung tissues against various lung injuries. However, the mechanism of ONO-5046 to HNE is still unclear. In this study, the binding properties of ONO-5046 were investigated through (1)H NMR, molecular docking, and bioassay methods to understand the effect of ONO-5046 to HNE. The proton spin-lattice relaxation rate and molecular rotational correlation time results indicated that ONO-5046 has higher affinity with HNE. The molecular docking study showed that ONO-5046 is perfectly matched for the primary enzyme specificity pocket (S1 pocket), and is tightly bound to this pocket of HNE through hydrophobic and hydrogen bonding interactions. The results of both methods were validated through analysis of the HNE inhibitory activity bioassay of ONO-5046 with an IC(50) value of 87.05 nM. Our data suggested that ONO-5046 could bind to HNE through direct interaction, and that molecular docking and NMR methods are valid approaches to survey new HNEI.     

Effects of an Anticarcinogenic Bowman-Birk Protease Inhibitor on Purified 20S Proteasome and MCF-7 Breast Cancer Cells

Proteasome inhibitors have been described as an important target for cancer therapy due to their potential to regulate the ubiquitin-proteasome system in the degradation pathway of cellular proteins. Here, we reported the effects of a Bowman-Birk-type protease inhibitor, the Black-eyed pea Trypsin/Chymotrypsin Inhibitor (BTCI), on proteasome 20S in MCF-7 breast cancer cells and on catalytic activity of the purified 20S proteasome from horse erythrocytes, as well as the structural analysis of the BTCI-20S proteasome complex. In vitro experiments and confocal microscopy showed that BTCI readily crosses the membrane of the breast cancer cells and co-localizes with the proteasome in cytoplasm and mainly in nucleus. Indeed, as indicated by dynamic light scattering, BTCI and 20S proteasome form a stable complex at temperatures up to 55°C and at neutral and alkaline pHs. In complexed form, BTCI strongly inhibits the proteolytic chymotrypsin-, trypsin- and caspase-like activities of 20S proteasome, indicated by inhibition constants of 10⁻⁷ M magnitude order. Besides other mechanisms, this feature can be associated with previously reported cytostatic and cytotoxic effects of BTCI in MCF-7 breast cancer cells by means of apoptosis.     

Structural insights into IKKβ inhibition by natural products staurosporine and quercetin

This Letter describes the results of two combined approaches: homology modeling and molecular docking studies, in order to propose the molecular basis of IKKβ inhibition by staurosporine and quercetin as ATP-competitive inhibitors. The results provides a rationale and structural frameworks for designing potent ATP binding-site inhibitors of IKKβ, which is an attractive drug target for inflammatory diseases and has been found to be responsible for some of the already observed pharmacological effects for marketed drugs.     

BTCI enhances guanylin-induced natriuresis and promotes renal glomerular and tubular effects.

Guanylin and uroguanylin are small cysteine-rich peptides involved in the regulation of fluid and electrolyte homeostasis through binding and activation of guanylyl cyclases signaling molecules expressed in intestine and kidney. Guanylin is less potent than uroguanylin as a natriuretic agent and is degraded in vitro by chymotrypsin due to unique structural features in the bioactive moiety of the peptide. Thus, the aim of this study was to verify whether or not guanylin is degraded by chymotrypsin-like proteases present in the kidney brush-border membranes. The isolated perfused rat kidney assay was used in this regard. Guanylin (0.2 microM) induced no changes in kidney function. However, when pretreated by the black-eyed pea trypsin and chymotrypsin inhibitor (BTCI - 1.0 microM; guanylin - 0.2 microM) it promoted increases in urine flow (DeltaUF of 0.25 +/- 0.09 mL.g(-1)/min, P < 0.05) and Na+ excretion (% Delta ENa+ of 18.20 +/- 2.17, P < 0.05). BTCI (1.0 microM) also increased %ENa+ (from 22.8 +/- 1.30 to 34.4 +/- 3.48, P < 0.05, 90 minutes). Furthermore, BTCI (3.0 microM) induced increases in glomerular filtration rate (GFR; from 0.96 +/- 0.02 to 1.28 0.02 mL.g(-1)/min, P < 0.05, 60 minutes). The present paper strongly suggests that chymotrypsin-like proteases play a role in renal metabolism of guanylin and describes for the first time renal effects induced by a member of the Bowman-Birk family of protease inhibitors.     

基于作用机制的氨基恶唑啉呫吨类BACE1抑制剂的分子设计

天冬氨酰蛋白酶(β-site amyloid precursor protein cleaving enzyme 1, BACE1)作为治疗阿尔兹海默症的潜在靶点,其抑制剂的开发已成为医学领域的重要研究方向。本文以59个氨基恶唑啉呫吨类BACE1抑制剂为研究对象,运用比较分子相似性指数(comparative molecular similarity index, CoMSIA)和分子对接方法,深入挖掘影响抑制剂活性的特征结构,以及抑制剂与BACE1间的结合模式和作用力类型,并以此为基础设计新型抑制剂并预测其活性。CoMSIA模拟结果表明,由立体场、静电场、疏水场和氢键供体场4个场组合建立的构效关系模型具有较强的预测能力,交叉验证相关系数Q2=0.48, 非交叉验证相关系数Rncv2=0.94, 外部预测相关系数Rpre2=0.85;通过分子对接,发现抑制剂占据了靶标的S3、S1和S2'位点,与BACE1之间的结合主要是通过氢键作用力和π-π堆积作用实现的;占据S2'位点的R取代基是立体场、静电场和疏水场影响的敏感区域,氨基恶唑啉核心官能团是氢键供体场的敏感区域。基于以上分析获得的抑制剂特征结构信息及其与蛋白质受体的作用机制,成功设计出了新的分子并预测了抑制活性。实验所得模型和信息,为后续新型BACE1抑制剂的结构优化和改造提供了重要理论依据     

Searching the conformational complexity and binding properties of HDAC6 through docking and molecular dynamic simulations

Histone deacetylases (HDACs) are a family of proteins involved in the deacetylation of histones and other non-histones substrates. HDAC6 belongs to class II and shares similar biological functions with others of its class. Nevertheless, its three-dimensional structure that involves the catalytic site remains unknown for exploring the ligand recognition properties. Therefore, in this contribution, homology modeling, 100-ns-long Molecular Dynamics (MD) simulation and docking calculations were combined to explore the conformational complexity and binding properties of the catalytic domain 2 from HDAC6 (DD2-HDAC6), for which activity and affinity toward five different ligands have been reported. Clustering analysis allowed identifying the most populated conformers present during the MD simulation, which were used as starting models to perform docking calculations with five DD2-HDAC6 inhibitors: Cay10603 (CAY), Rocilinostat (RCT), Tubastatin A (TBA), Tubacin (TBC), and Nexturastat (NXT), and then were also submitted to 100-ns-long MD simulations. Docking calculations revealed that the five inhibitors bind at the DD2-HDAC6 binding site with the lowest binding free energy, the same binding mode is maintained along the 100-ns-long MD simulations. Overall, our results provide structural information about the molecular flexibility of apo and holo DD2-HDAC6 states as well as insight of the map of interactions between DD2-HDAC6 and five well-known DD2-HDAC6 inhibitors allowing structural details to guide the drug design. Finally, we highlight the importance of combining different theoretical approaches to provide suitable structural models for structure-based drug design.     

Ligand-based pharmacophore modeling and docking studies on vitamin D receptor inhibitors

In recent years, pharmacophore modeling and molecular docking approaches have been extensively used to characterize the structural requirements and explore the conformational space of a ligand in the binding pocket of the selected target protein. Herein, we report a pharmacophore modeling and molecular docking of 45 compounds comprising of the indole scaffold as vitamin D receptor (VDR) inhibitors. Based on the selected best hypothesis (DRRRR.61), an atom-based three-dimensional quantitative structure-activity relationships model was developed to rationalize the structural requirement of biological activity modulating components. The developed model predicted the binding affinity for the training set and test set with R²_(training) = 0.8869 and R²_(test) = 0.8139, respectively. Furthermore, molecular docking and dynamics simulation were performed to understand the underpinning of binding interaction and stability of selected VDR inhibitors in the binding pocket. In conclusion, the results presented here, in the form of functional and structural data, agreed well with the proposed pharmacophores and provide further insights into the development of novel VDR inhibitors with better activity.     

Gene ontology improves template selection in comparative protein docking

Structural characterization of protein-protein interactions is essential for our ability to study life processes at the molecular level. Computational modeling of protein complexes (protein docking) is important as the source of their structure and as a way to understand the principles of protein interaction. Rapidly evolving comparative docking approaches utilize target/template similarity metrics, which are often based on the protein structure. Although the structural similarity, generally, yields good performance, other characteristics of the interacting proteins (eg, function, biological process, and localization) may improve the prediction quality, especially in the case of weak target/template structural similarity. For the ranking of a pool of models for each target, we tested scoring functions that quantify similarity of Gene Ontology (GO) terms assigned to target and template proteins in three ontology domains—biological process, molecular function, and cellular component (GO-score). The scoring functions were tested in docking of bound, unbound, and modeled proteins. The results indicate that the combined structural and GO-terms functions improve the scoring, especially in the twilight zone of structural similarity, typical for protein models of limited accuracy.     

The study of G-protein coupled receptor oligomerization with computational modeling and bioinformatics

To achieve a structural context for the analysis of G-protein coupled receptor (GPCR) oligomers, molecular modeling must be used to predict the corresponding interaction interfaces. The task is complicated by the paucity of detailed structural data at atomic resolution, and the large number of possible modes in which the bundles of seven transmembrane (TM) segments of the interacting GPCR monomers can be packed together into dimers and/or higher-order oligomers. Approaches and tools offered by bioinformatics can be used to reduce the complexity of this task and, combined with computational modeling, can serve to yield testable predictions for the structural properties of oligomers. Most of the bioinformatics methods take advantage of the evolutionary relation that exists among GPCRs, as expressed in their sequences and measurable in the common elements of their structural and functional features. These common elements are responsible for the presence of detectable patterns of motifs and correlated mutations evident from the alignment of the sequences of these complex biological systems. The decoding of these patterns in terms of structural and functional determinants can provide indications about the most likely interfaces of dimerization/oligomerization of GPCRs. We review here the main approaches from bioinformatics, enhanced by computational molecular modeling, that have been used to predict likely interfaces of dimerization/oligomerization of GPCRs, and compare results from their application to rhodopsin-like GPCRs. A compilation of the most frequently predicted GPCR oligomerization interfaces points to specific regions of TMs 4-6.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Serine proteinase inhibitors from nematodes and the arms race between host and pathogen

Serine proteinase inhibitors are encoded by a large gene family of long evolutionary standing. Recent discoveries of parasite proteins that inhibit human serine proteinases, together with the complete genomic sequence from Caenorhabditis elegans, have provided a set of new serine proteinase inhibitors from more primitive metazoan animals such as nematodes. The structural features (e.g. reactive centre residues), gene organization (including intron arrangements) and inhibitory function and targets (e.g. inflammatory and coagulation pathway proteinase) all contribute important new insights into proteinase inhibitor evolution. Some parasite products have evolved that block enzymes in the mammalian host, but the human host responds with a significant immune response to the parasite inhibitors. Thus, infection produces a finely balanced conflict between host and pathogen at the molecular level, and this might have accelerated the evolution of these proteins in parasitic species as well as their hosts.     

Homology modeling,molecular dynamic simulation and docking studies of cyclin dependent kinase 1

In order to develop promising cyclin dependent kinase 1 inhibitors, homology modeling, docking and molecular dynamic simulation techniques were applied to get insight into the functional and structural properties of cyclin dependent kinase 1 (CDK1). Since there is no reported CDK1 crystal structural data, the three dimensional structure of CDK1 was constructed based on homology modeling. An extensive dynamic simulation was also performed on a Flavopiridol-CDK1 complex for probing the binding pattern of Flavopiridol in the active site of CDK1. The binding modes of other inhibitors to CDK1 were also proposed by molecular docking. The structural requirement for developing more potent CDK1 inhibitors was obtained by the above-mentioned molecular simulations and pharmacophore modeling.     

Molecular design of new aggrecanases-2 inhibitors

Aggrecanases-2 is a very important potential drug target for the treatment of osteoarthritis. In this study, a series of known aggrecanases-2 inhibitors was analyzed by the technologies of three-dimensional quantitative structure–activity relationships (3D-QSAR) and molecular docking. Two 3D-QSAR models, which based on comparative molecular field analysis (CoMFA) and comparative molecular similarity analysis (CoMSIA) methods, were established. Molecular docking was employed to explore the details of the interaction between inhibitors and aggrecanases-2 protein. According to the analyses for these models, several new potential inhibitors with higher activity predicted were designed, and were supported by the simulation of molecular docking. This work propose the fast and effective approach to design and prediction for new potential inhibitors, and the study of the interaction mechanism provide a better understanding for the inhibitors binding into the target protein, which will be useful for the structure-based drug design and modifications.