BackgroundMolecular docking methods are commonly used for predicting binding modes and energies of ligands to proteins. For accurate complex geometry and binding energy estimation, an appropriate method for calculating partial charges is essential. AutoDockTools software, the interface for preparing input files for one of the most widely used docking programs AutoDock 4, utilizes the Gasteiger partial charge calculation method for both protein and ligand charge calculation. However, it has already been shown that more accurate partial charge calculation - and as a consequence, more accurate docking- can be achieved by using quantum chemical methods. For docking calculations quantum chemical partial charge calculation as a routine was only used for ligands so far. The newly developed Mozyme function of MOPAC2009 allows fast partial charge calculation of proteins by quantum mechanical semi-empirical methods. Thus, in the current study, the effect of semi-empirical quantum-mechanical partial charge calculation on docking accuracy could be investigated.ResultsThe docking accuracy of AutoDock 4 using the original AutoDock scoring function was investigated on a set of 53 protein ligand complexes using Gasteiger and PM6 partial charge calculation methods. This has enabled us to compare the effect of the partial charge calculation method on docking accuracy utilizing AutoDock 4 software. Our results showed that the docking accuracy in regard to complex geometry (docking result defined as accurate when the RMSD of the first rank docking result complex is within 2 Å of the experimentally determined X-ray structure) significantly increased when partial charges of the ligands and proteins were calculated with the semi-empirical PM6 method.Out of the 53 complexes analyzed in the course of our study, the geometry of 42 complexes were accurately calculated using PM6 partial charges, while the use of Gasteiger charges resulted in only 28 accurate geometries. The binding affinity estimation was not influenced by the partial charge calculation method - for more accurate binding affinity prediction development of a new scoring function for AutoDock is needed.ConclusionOur results demonstrate that the accuracy of determination of complex geometry using AutoDock 4 for docking calculation greatly increases with the use of quantum chemical partial charge calculation on both the ligands and proteins.

The application of natural product inspired synthesis has identified novel antibacterial tetramic acids which exhibit wide ranging antibacterial activity, and which provide potential lead structures for antibacterial drug discovery. Their phenotypic activity appears to correlate with action at two enzymes, UPPS and RNAP, which operate in independent metabolic pathways. SAR maps and identification of their relevant binding sites by molecular modelling has been achieved, and characterisation of the most active compounds suggests that these systems offer potential for topical antibiotics but that for oral and injectable use further optimisation is required.

Human P-glycoprotein (P-gp) is an ATP-binding cassette multidrug transporter that confers resistance to a wide range of chemotherapeutic agents in cancer cells by active efflux of the drugs from cells. P-gp also plays a key role in limiting oral absorption and brain penetration and in facilitating biliary and renal elimination of structurally diverse drugs. Thus, identification of drugs or new molecular entities to be P-gp substrates is of vital importance for predicting the pharmacokinetics, efficacy, safety, or tissue levels of drugs or drug candidates. At present, publicly available, reliable in silico models predicting P-gp substrates are scarce. In this study, a support vector machine (SVM) method was developed to predict P-gp substrates and P-gp-substrate interactions, based on a training data set of 197 known P-gp substrates and non-substrates collected from the literature. We showed that the SVM method had a prediction accuracy of approximately 80% on an independent external validation data set of 32 compounds. A homology model of human P-gp based on the X-ray structure of mouse P-gp as a template has been constructed. We showed that molecular docking to the P-gp structures successfully predicted the geometry of P-gp-ligand complexes. Our SVM prediction and the molecular docking methods have been integrated into a free web server (http://pgp.althotas.com), which allows the users to predict whether a given compound is a P-gp substrate and how it binds to and interacts with P-gp. Utilization of such a web server may prove valuable for both rational drug design and screening.

No abstract

Human serum albumin (HSA), the most abundant plasma protein is well known for its extraordinary binding capacity for both endogenous and exogenous substances, including a wide range of drugs. Interaction with the two principal binding sites of HSA in subdomain IIA (site 1) and in subdomain IIIA (site 2) controls the free, active concentration of a drug, provides a reservoir for a long duration of action and ultimately affects the ADME (absorption, distribution, metabolism, and excretion) profile. Due to the continuous demand to investigate HSA binding properties of novel drugs, drug candidates and drug-like compounds, a support vector machine (SVM) model was developed that efficiently predicts albumin binding. Our SVM model was integrated to a free, web-based prediction platform (http://albumin.althotas.com). Automated molecular docking calculations for prediction of complex geometry are also integrated into the web service. The platform enables the users (i) to predict if albumin binds the query ligand, (ii) to determine the probable ligand binding site (site 1 or site 2), (iii) to select the albumin X-ray structure which is complexed with the most similar ligand and (iv) to calculate complex geometry using molecular docking calculations. Our SVM model and the potential offered by the combined use of in silico calculation methods and experimental binding data is illustrated.

No abstract

No abstract

ABSTRACT:Previous studies in our laboratory showed that among cDNAexpressed human cytochrome P450 (P450) supersomes, CYP2C19 was the most active in methoxychlor-O-demethylation (Hu et al., 2004). However, based on the lack of inhibition of methoxychlor-O-demethylation by monoclonal anti-CYP2C19 antibodies in human liver microsomes (HLM), CYP2C19 did not seem to catalyze that reaction in HLM. By contrast, CYP2C9, much less active than CYP2C19 in supersomes, was the most active in HLM. The current study examines whether the lack of methoxychlor-O-demethylation by CYP2C19 in HLM was due to CYP2C19 exhibiting inferior competition for the NADPH-cytochrome P450 reductase (CPR) versus CYP2C9 and explores the interactions between CYP2C9 and CYP2C19 in a singular and binary complex of a reconstituted system. When reconstituted with CPR, cytochrome b 5 , and lipid, purified CYP2C19 and CYP2C9 catalyzed methoxychlor-O-demethylation.However, whereas equimolar CPR to CYP2C9 supported maximal rates of methoxychlor demethylation and diclofenac hydroxylation, the rate of methoxychlor demethylation by CYP2C19 was not fully saturated, even with a 9-fold molar excess of CPR over CYP2C19. This behavior of CYP2C19 was also observed with S-mephenytoin as the substrate. When a binary reconstitution system was prepared by mixing CYP2C9 and CYP2C19 enzymes, methoxychlor-O-demethylation and S-mephenytoin hydroxylation by CYP2C19 were dramatically inhibited. Inhibition depended on the amount of CPR and substrate used. By contrast, in the incubation containing CYP2C9, diclofenac hydroxylation was activated by the presence of CYP2C19. These results show that interactions among P450 enzymes can modulate their catalytic rates, which depend on the substrate undergoing metabolism.