A comparative study of the single crystal X-ray determination and molecular modelling of the binding of oligomycin to ATP Synthase

Recently published X-ray structures of three common forms, A, B and C, of oligomycin, including absolute configurations, are investigated to examine their binding to ATP Synthase. The X-ray studies reveal regions with differences in three-dimensional structure and hydrogen bonding propensity between the oligomycins, which may be associated with their potential to bind to sites on ATP Synthase. Computational docking studies carried out using MOE with the X-ray structures and an homology model of the F_O domain of ATP Synthase from Escherichia coli, are used to derive an induced fit pocket. Docking of all oligomycins to this pocket indicate that the B and C forms bind more tightly than the A form. Consideration of the single crystal X-ray data alone indicate the B form may be the best inhibitor and that O(24) is the most important ligating group for binding, this is supported by the docking data. The latter reveals Asn214 and other key proton translocating residues to be the main residues contacted by the inhibitor. These data allow the binding modes of different forms of oligomycin to be deduced from X-ray single crystal data supported by molecular modelling and computational docking studies.     

Dynamic mechanism for the autophosphorylation of CheA histidine kinase: molecular dynamics simulations

The two-component system (TCS) is an important signal transduction component for most bacteria. This signaling pathway is mediated by histidine kinases via autophosphorylation between P1 and P4 domains. Taking chemotaxis protein CheA as a model of TCS, the autophosphorylation mechanism of the TCS histidine kinases has been investigated in this study by using a computational approach integrated homology modeling, ligand-protein docking, protein-protein docking, and molecular dynamics (MD) simulations. Four nanosecond-scale MD simulations were performed on the free P4 domain, P4-ATP, P4-TNPATP, and P1-P4-ATP complexes, respectively. Upon its binding to the binding pocket of P4 with a folded conformation, ATP gradually extends to an open state with help from a water molecule. Meanwhile, ATP forms two hydrogen bonds with His413 and Lys494 at this state. Because of the lower energy of the folded conformations, ATP shrinks back to its folded conformations, leading to the rupture of the hydrogen bond between ATP and Lys494. Consequently, Lys494 moves away from the pocket entrance, resulting in an open of the ATP lid of P4. It is the open state of P4 that can bind tightly to P1, where the His45 of P1 occupies a favorable position for its autophosphorylation from ATP. This indicates that ATP is not only a phosphoryl group donor but also an activator for CheA phosphorylation. Accordingly, a mechanism of the autophosphorylation of CheA is proposed as that the ATP conformational switch triggers the opening of the ATP lid of P4, leading to P1 binding tightly, and subsequently autophosphorylation from ATP to P1.     

Computational Studies and Drug Design for HIV-1 Reverse Transcriptase Inhibitors of 3′,4′-di-O-(S)-camphanoyl-(+)-cis-Khellactone (DCK) Analogs

Molecular docking and molecular dynamics simulation were applied to study the binding mode of 3',4'-di-O-(S)-camphanoyl-(+)-cis-khellactone (DCK) analogs anti-HIV inhibitors with HIV-1 RT. The results suggest that there is a strong hydrogen bond between DCK O16 and NH of Lys101, and that DCK analogues might act similarly as other types of HIV-1 RT inhibitors. The investigation about drug resistance for DCK shows no remarkable influence on the most frequently observed mutation K103N of HIV-1 RT. Based on the proposed mechanism, some new structures were designed and predicted by a SVM model. All compounds exhibited potent inhibitory activities against HIV replication in H9 lymphocytes with EC50 values lower than 1.95 microM. The rationality of the method was validated by experimental results.     

Pharmacophore generation,atom-based 3D-QSAR and molecular dynamics simulation analyses of pyridine-3-carboxamide-6-yl-urea analogues as potential gyrase B inhibitors

DNA gyrase subunit B (GyrB) is an attractive drug target for the development of antibacterial agents with therapeutic potential. In the present study, computational studies based on pharmacophore modelling, atom-based QSAR, molecular docking, free binding energy calculation and dynamics simulation were performed on a series of pyridine-3-carboxamide-6-yl-urea derivatives. A pharmacophore model using 49 molecules revealed structural and chemical features necessary for these molecules to inhibit GyrB. The best fitted model AADDR.13 was generated with a coefficient of determination (r²) of 0.918. This model was validated using test set molecules and had a good r² of 0.78. 3D contour maps generated by the 3D atom-based QSAR revealed the key structural features responsible for the GyrB inhibitory activity. Extra precision molecular docking showed hydrogen bond interactions with key amino acid residues of ATP-binding pocket, important for inhibitor binding. Further, binding free energy was calculated by the MM-GBSA rescoring approach to validate the binding affinity. A 10 ns MD simulation of inhibitor #47 showed the stability of the predicted binding conformations. We identified 10 virtual hits by in silico high-throughput screening. A few new molecules were also designed as potent GyrB inhibitors. The information obtained from these methodologies may be helpful to design novel inhibitors of GyrB.     

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.     

Investigating the hydrogen-bond acceptor site of the nicotinic pharmacophore model: a computational and experimental study using epibatidine-related molecular probes

The binding mode of nicotinic agonists has been thoroughly investigated in the last decades. It is now accepted that the charged amino group is bound by a cation-π interaction to a conserved tryptophan residue, and that the aromatic moiety is projected into a hydrophobic pocket deeply located inside the binding cleft. A hydrogen bond donor/acceptor, maybe a water molecule solvating this receptor subsite, contributes to further stabilize the nicotinic ligands. The position of this water molecule has been established by several X-ray structures of the acetylcholine-binding protein. In this study, we computationally analyzed the role of this water molecule as a putative hydrogen bond donor/acceptor moiety in the agonist binding site of the three most relevant heteromeric (α4β2, α3β4) and homomeric (α7) neuronal nicotinic acetylcholine receptor (nAChR) subtypes. Our theoretical investigation made use of epibatidine 1 and deschloroepibatidine 2 as molecular probes, and was then extended to their analogues 3 and 4, which were subsequently synthesized and tested at the three target receptor subtypes. Although the pharmacological data for the new ligands 3 and 4 indicated a reduction of the affinity at the studied nAChRs with respect to reference agonists, a variation of the selectivity profile was clearly evidenced.     

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25相似文献   

Use of MM-PBSA in reproducing the binding free energies to HIV-1 RT of TIBO derivatives and predicting the binding mode to HIV-1 RT of efavirenz by docking and MM-PBSA.

In this work, a new ansatz is presented that combines molecular dynamics simulations with MM-PBSA (Molecular Mechanics Poisson-Boltzmann/surface area) to rank the binding affinities of 12 TIBO-like HIV-1 RT inhibitors. Encouraging results have been obtained not only for the relative binding free energies, but also for the absolute ones, which have a root-mean-square deviation of 1.0 kcal/mol (the maximum error is 1.89 kcal/mol). Since the root-mean-square error is rather small, this approach can be reliably applied in ranking the ligands from the databases for this important target. Encouraged by the results, we decided to apply MM-PBSA combined with molecular docking to determine the binding mode of efavirenz SUSTIVA(TM) another promising HIV-1 RT inhibitor for which no ligand-protein crystal structure had been published at the time of this work. To proceed, we define the following ansatz: Five hundred picosecond molecular dynamics simulations were first performed for the five binding modes suggested by DOCK 4.0, and then MM-PBSA was carried out for the collected snapshots. MM-PBSA successfully identified the correct binding mode, which has a binding free energy about 7 kcal/mol more favorable than the second best mode. Moreover, the calculated binding free energy (-13.2 kcal/mol) is in reasonable agreement with experiment (-11.6 kcal/mol). In addition, this procedure was also quite successful in modeling the complex and the structure of the last snapshot was quite close to that of the measured 2,3 A resolution crystal (structure the root-mean-square deviation of the 54 C(alpha) around the binding site and the inhibitor is 1.1 A). We want to point out that this result was achieved without prior knowledge of the structure of the efavirenz/RT complex. Therefore, molecular docking combined with MD simulations followed by MM-PBSA analysis is an attractive approach for modeling protein complexes a priori.     

Automated high throughput analysis of antiretroviral drugs in dried blood spots

For therapeutic drug monitoring in remote settings, dried blood spots (DBS) are particularly advantageous, as blood sample collection and handling is uncomplicated. The aim of this study was to develop and validate an automated extraction method for the analysis of nevirapine, efavirenz and lopinavir in DBS samples. Automated extraction was performed with methanol : water (70 : 30 v /v ), using a DBS‐MS 500 autosampler coupled to a liquid chromatography tandem mass spectrometry system. The autosampler used digital images of each DBS to position the extraction head, sprayed 10 μl of internal standard onto each DBS and extracted a 4‐mm disc (Ø) from the centre of each spot by unilateral flow using 25‐μl extraction solvent. The analytes were baseline separated on a pentafluorophenyl column and analysed by using electrospray ionization with multiple reaction monitoring in positive polarity mode for nevirapine and lopinavir and in negative mode for efavirenz. The method was linear between 10 and 10 000 ng/ml for all analytes. Automated sample extraction resulted in consistent recoveries (nevirapine: 70 ± 6%, efavirenz: 63 ± 11% and lopinavir: 60 ± 10%) and matrix effects between different donors and concentration levels. Intra‐day and inter‐day accuracy and precision deviations were ≤15%. Manual and automated extractions of DBS samples collected within the framework of an adherence assessment study in rural Tanzania showed good agreements with deviations of less than 10%. Our study highlights that therapeutic drug monitoring samples obtained in the resource‐constrained setting of rural Africa can be reliably determined by automated extraction of DBS. Overall, automatization improved method sensitivity and facilitates analysis of large sample numbers. Copyright © 2017 John Wiley & Sons, Ltd.     

RK-805, an endothelial-cell-growth inhibitor produced by Neosartorya sp., and a docking model with methionine aminopeptidase-2

We found that a fungus Neosartorya sp. produced an angiogenesis inhibitor, RK-805. By spectroscopic analyses and semi-synthetic methods from fumagillin, the structure of RK-805 was identified as 6-oxo-6-deoxyfumagillol, which has not been reported as a natural product. RK-805 preferentially inhibited the growth of human umbilical vein endothelial cells (HUVECs) rather than that of human normal fibroblast in cell proliferation assays and blocked endothelial cell migration induced by vascular endothelial growth factor (VEGF). Moreover, RK-805 selectively inhibited methionine aminopeptidase-2 (MetAP2), but not methionine aminopeptidase-1 (MetAP1). The docked structure of RK-805 complexed with human MetAP2 indicated that not only a covalent bond between a nucleophilic imidazole nitrogen atom of His231 and the carbon of the reactive spirocyclic epoxide of RK-805, but also a hydrogen bond between NH (Asn329) and the carbonyl group of RK-805 at C-6 promote close contact in the binding pocket of the enzyme. Taken together, these results suggest that structure activity relationships of RK-805 derivatives at both C-4 and C-6, in comparison with ovalicin and TNP-470, would be useful for development of new angiogenesis inhibitors.     

白色念珠菌N-肉豆蔻酰基转移酶中药物结合位点的MCSS分析

采用多拷贝同时搜寻方法(MCSS)分析得到了CaNMT活性位点的疏水区域、氢键结合位点和负电性区域. MCSS计算结果显示, CaNMT活性位点有两个疏水性比较强的区域: 一个由Tyr107, Tyr109, Val108, Phe117, Phe123, Ala127, Phe176和Leu337等残基组成; 另一个由Phe115, Phe240和Phe339组成. CaNMT活性位点发现有两个氢键作用区域, 其中Tyr119, His227, Asn392和Leu451是与已有抑制剂的氢键结合位点, Tyr107, Asn175, Thr211和Asp412是新发现的氢键结合位点, 而且在NMT家族中高度稳定, 它们对设计新结构类型的CaNMT抑制剂具有重要作用. Leu451是负电性兼氢键作用位点, 是抑制剂设计时所必需考虑的位点.