A tellurium-based cathepsin B inhibitor: Molecular structure, modelling, molecular docking and biological evaluation

The crystallographically determined structure of biologically active 4,4-dichloro-1,3-diphenyl-4-telluraoct-2-en-1-one, 3, shows the coordination geometry for Te to be distorted psi-pentagonal bipyramidal based on a C2OCl3(lone pair) donor set. Notable is the presence of an intramolecular axial Te center dot center dot center dot O (carbonyl) interaction, a design element included to reduce hydrolysis. Raman and molecular modelling studies indicate the persistence of the Te center dot center dot center dot O(carbonyl) interaction in the solution (CHCl3) and gasphases, respectively. Docking studies of 3' (i.e. original 3 less one chloride) with Cathepsin B reveals a change in the configuration about the vinyl C = C bond. i.e. to E from Z (crystal structure). This isomerism allows the optimisation of interactions in the complex which features a covalent Te-SGCys29 bond. Crucially, the E configuration observed for 3' allows for the formation of a hypervalent Te center dot center dot center dot O interaction as well as an O center dot center dot center dot H-O hydrogen bond with the Gly27 and Glu122 residues, respectively. Additional stabilisation is afforded by a combination of interactions spanning the S1, S2, S1' and S2' sub-sites of Cathepsin B. The greater experimental inhibitory activity of 3 compared with analogues is rationalised by the additional interactions formed between 3' and the His110 and His111 residues in the occluding loop, which serve to hinder the entrance to the active site. (C) 2012 Elsevier B.V. All rights reserved.

Computational methods for prediction of T-cell epitopes - a framework for modelling, testing, and applications

Computational models complement laboratory experimentation for efficient identification of MHC-binding peptides and T-cell epitopes. Methods for prediction of MHC-binding peptides include binding motifs, quantitative matrices, artificial neural networks, hidden Markov models, and molecular modelling. Models derived by these methods have been successfully used for prediction of T-cell epitopes in cancer, autoimmunity, infectious disease, and allergy. For maximum benefit, the use of computer models must be treated as experiments analogous to standard laboratory procedures and performed according to strict standards. This requires careful selection of data for model building, and adequate testing and validation. A range of web-based databases and MHC-binding prediction programs are available. Although some available prediction programs for particular MHC alleles have reasonable accuracy, there is no guarantee that all models produce good quality predictions. In this article, we present and discuss a framework for modelling, testing, and applications of computational methods used in predictions of T-cell epitopes. (C) 2004 Elsevier Inc. All rights reserved.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.

A facile and single pot strategy for the synthesis of novel naphthyridine derivatives under microwave irradiation conditions using ZnCl2 as catalyst, evaluation of AChE inhibitory activity, and molecular modeling studies

A series of novel naphthyridine derivatives 3 and 4 was prepared from substituted pyridine 2 and ketones using ZnCl2 as catalyst under microwave irradiation conditions. All the compounds were evaluated for AChE inhibitory activity and promising compounds 3d, 3e, 4b, and 4g was identified. Representative compounds 3d and 3e were found to show insignificant THLE-2 liver cell viability/toxicity. The binding mode between X-ray crystal structure of human AChE and compounds was studied using molecular docking method and fitness scores were found to be in good correlation with the activity data.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Delineation and modelling of a nucleolar retention signal in the coronavirus nucleocapsid protein

Unlike nuclear localization signals, there is no obvious consensus sequence for the targeting of proteins to the nucleolus. The nucleolus is a dynamic subnuclear structure which is crucial to the normal operation of the eukaryotic cell. Studying nucleolar trafficking signals is problematic as many nucleolar retention signals (NoRSs) are part of classical nuclear localization signals (NLSs). In addition, there is no known consensus signal with which to inform a study. The avian infectious bronchitis virus (IBV), coronavirus nucleocapsid (N) protein, localizes to the cytoplasm and the nucleolus. Mutagenesis was used to delineate a novel eight amino acid motif that was necessary and sufficient for nucleolar retention of N protein and colocalize with nucleolin and fibrillarin. Additionally, a classical nuclear export signal (NES) functioned to direct N protein to the cytoplasm. Comparison of the coronavirus NoRSs with known cellular and other viral NoRSs revealed that these motifs have conserved arginine residues. Molecular modelling, using the solution structure of severe acute respiratory (SARS) coronavirus N-protein, revealed that this motif is available for interaction with cellular factors which may mediate nucleolar localization. We hypothesise that the N-protein uses these signals to traffic to and from the nucleolus and the cytoplasm.

Analysis of the molecular association of rifampicin with hydroxypropyl-²-cyclodextrin

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Protease-inhibiting, molecular modeling and antimicrobial activities of extracts and constituents from Helichrysum foetidum and Helichrysum mechowianum (compositae)

Background Helichrysum species are used extensively for stress-related ailments and as dressings for wounds normally encountered in circumcision rites, bruises, cuts and sores. It has been reported that Helichysum species are used to relief abdominal pain, heart burn, cough, cold, wounds, female sterility, menstrual pain. Results From the extracts of Helichrysum foetidum (L.) Moench, six known compounds were isolated and identified. They were 7, 4′-dihydroxy-5-methoxy-flavanone (1), 6′-methoxy-2′,4, 4′-trihydroxychalcone (2), 6′-methoxy-2′,4-dihydroxychalcone -4′-O-β-D-glucoside (3), apigenin (4), apigenin-7-O-β-D-glucoside (5), kaur-16-en-18-oic acid (6) while two known compounds 3,5,7-trihydroxy-8-methoxyflavone (12), 4,5-dicaffeoyl quinic acid (13) together with a mixture of phytosterol were isolated from the methanol extract of Helichrysum mechowianum Klatt. All the compounds were characterized by spectroscopic and mass spectrometric methods, and by comparison with literature data. Both extracts and all the isolates were screened for the protease inhibition, antibacterial and antifungal activities. In addition, the phytochemical profiles of both species were investigated by ESI-MS experiments. Conclusions These results showed that the protease inhibition assay of H. foetidum could be mainly attributed to the constituents of flavonoids glycosides (3, 5) while the compound (13) from H. mechowianum contributes to the stomach protecting effects. In addition, among the antibacterial and antifungal activities of all the isolates, compound (6) was found to possess a potent inhibitor effect against the tested microorganisms. The heterogeneity of the genus is also reflected in its phytochemical diversity. The differential bioactivities and determined constituents support the traditional use of the species. Molecular modelling was carried out by computing selected descriptors related to drug absorption, distribution, metabolism, excretion and toxicity (ADMET).

Untersuchungen zu COX/LOX-Inhibitoren mit Hydroxylradikalfängereigenschaften – Synthese, Testung und Modelling

In der vorliegenden Arbeit werden 52 Verbindungen beschrieben, welche auf COX/LOX-Inhibition mit zusätzlichen Hydroxylradikalfängereigenschaften getestet worden sind. rnEs war möglich eine neue Synthesestrategie für noch nicht beschriebene 4,5-Diarylisoselenazole zu entwickeln und eine vorhandene Synthese für Isothiazoliumchloride von zwei Stufen, mit mäßigen Ausbeuten, auf eine Stufe, mit hoher Ausbeute, zu verkürzen.rnEs wurden mehrere COX-Inhibitoren identifiziert. MSD4a, MSD4h, MSD5a und MSD5h konnten als COX-1-, COX-2- und 5-LOX-Hemmer identifiziert werden. Besonders hervorzuheben ist die Verbindung MSD5h, die zusätzlich zur COX-1-, COX-2- und 5-LOX-Inhibition eine leichte Hemmung im Hydroxylradikalfänger-Assay zeigt, für die ein clog P-Wert von 2,65 berechnet wurde und die im XTT-Zytotoxizitätstestsystem, selbst bei einer Konzentration von 100 µM, kaum toxische Eigenschaften besitzt.rnWeiterhin war es möglich zu zeigen, dass Carbonsäuren gute Hydroxylradikalfängereigenschaften in unserem, auf der Fenton-Reaktion basierenden, Testsystem haben. Die Potenz der Carbonsäuren MSD8b und MSD11j im Vergleich zu den unwirksamen korrespondierenden Ester MSD8a und MSD11i führte zu Untersuchungen mit weiteren Carbonsäuren und deren Ester. Um den Wirkungsmechanismus zu erforschen wurde das Testsystem modifiziert, um eine Komplexierung der Eisenionen durch die Carbonsäuren auszuschließen. An Hand der Substanzen MSD8b und MSD11j wurde nachgewiesen, dass diese mit dem Hydroxylradikal reagieren, ohne zu decarboxylieren oder andere Zerfallsreaktionen einzugehen.rnZusätzlich zu den Untersuchungen der Enzym-Inhibition sowie des Hydroxylradikal-Scavenings wurden Molecular Modelling Studien durchgeführt. Die Ergebnisse der Dockingstudien in COX-1- (1eqg), COX-2- (1cx2) und in COX-1 mutierte COX-2-Kristallstrukturen (1cx2) führen zu einer kritischen Bewertung des folgenden Ansatzes: Es ist nicht unbedingt sinnvoll zuerst Strukturen mit dem Computer zu entwerfen und zu modeln und sie erst dann zu synthetisieren und in Enzym- oder Zellassays zu testen. Die Begründung dafür liegt in der Schwierigkeit einschätzen zu können, wie nah das gewählte Modell der Wirklichkeit ist. In den durchgeführten Dockingstudien konnte der sehr große Einfluss des kokristallisierten Liganden in der als Grundlage dienenden Kristallstruktur auf die Dockingergebnisse gezeigt werden. Durch einen zu kleinen kokristallisierten Liganden in der COX-1-Bindungstasche wurden als Ergebnis der Dockingstudie alle Verbindungen als nicht potent eingestuft, obwohl diese zum Teil im Enzymtestsystem wirksam waren. Dies konnte mit den Mutationsversuchen ausgeglichen werden. rnDeshalb kann man aus diesen Ergebnissen als Fazit ziehen, dass eine Strategie, Strukturen zu synthetisieren, in vitro zu testen und dabei die Strukturentwicklung mit Molecular Modelling Studien zu unterstützen, die Methode der Wahl darstellt.rn

Receptor activity-modifying proteins 2 and 3 generate adrenomedullin receptor subtypes with distinct molecular properties

Adrenomedullin (AM) is a peptide hormone with numerous effects in the vascular systems. AM signals through the AM1 and AM2 receptors formed by the obligate heterodimerization of a G protein-coupled receptor, the calcitonin receptor-like receptor (CLR), and receptor activity-modifying proteins (RAMP) 2 and 3, respectively. These different CLR-RAMP interactions yield discrete receptor pharmacology and physiological effects. The effective design of therapeutics that target the individual AM receptors is dependent on understanding the molecular details of the effects of RAMPs on CLR. To understand the role of RAMPs 2 and 3 on the activation and conformation of the CLR subunit of AM receptors we mutated 68 individual amino acids in the juxtamembrane region of CLR, a key region for activation of AM receptors and determined the effects on cAMP signalling. Sixteen CLR mutations had differential effects between the AM1 and AM2 receptors. Accompanying this, independent molecular modelling of the full-length AM-bound AM1 and AM2 receptors predicted differences in the binding pocket, and differences in the electrostatic potential of the two AM receptors. Druggability analysis indicated unique features that could be used to develop selective small molecule ligands for each receptor. The interaction of RAMP2 or RAMP3 with CLR induces conformational variation in the juxtamembrane region, yielding distinct binding pockets, probably via an allosteric mechanism. These subtype-specific differences have implications for the design of therapeutics aimed at specific AM receptors and for understanding the mechanisms by which accessory proteins affect G protein-coupled receptor function.

An investigation of the molecular pharmacology of G protein-coupled receptor 35

G protein-coupled receptors (GPCRs) are seven-pass integral membrane proteins that act as transducers of extracellular signals across the lipid bilayer. Their location and involvement in basic and pathological physiological processes has secured their role as key targets for pharmaceutical intervention. GPCRs are targeted by many of the best-selling drugs on the market and there are a substantial number of GPCRs that are yet to be characterised; these could offer interest for therapeutic targeting. GPR35 is one such receptor that, as a result of gene knockout and genome wide association studies, has attracted interest through its association with cardiovascular and gastrointestinal disease. Elucidation of the basic physiological function of GPR35 has, however, been difficult due a paucity of potent and selective ligands in addition to a lack of consensus on the endogenous ligand. Herein, a focussed drug discovery effort was carried out to identify agonists of GPR35. Various in vitro cellular assays were employed in conjunction with N- or C-terminally manipulated forms of the receptor to investigate GPR35’s signalling profile and to provide an assay format suitable for the characterisation of newly identified ligands. Although GPR35 associates with both Gαi/o and Gα13 families of small heterotrimeric G proteins, the G protein-independent β-arrestin-2 recruitment format was found to be the most suited to drug screening efforts. Small molecule compound screening, carried out in conjunction with the Medical Research Council Technology, identified compound 1 as the most potent ligand of human GPR35 reported at that time. However, the lower efficacy and potency of compound 1 at the rodent species orthologues of GPR35 prevented its use in in vivo studies. A subsequent effort, carried out with Novartis, focused on mast cell stabilisers as putative agonists of GPR35, revealed lodoxamide and bufrolin as highly potent agonists that activated human and rat GPR35 with equal potency. This finding offered–for the first time–the opportunity to employ the same GPR35 ligand between species at a similar concentration, an important factor to consider when translating rodent in vivo functional studies to those in man. Additionally, using molecular modelling and site directed mutagenesis studies, these newly identified compounds were used to aid characterisation of the ligand binding pockets of human and rat GPR35 to reveal the molecular basis of species selectivity at this receptor. In summary, this research effort presents GPR35 tool compounds that can now be used to dissect the basic biology of GPR35 and investigate its contribution to disease.

Homozygous Inactivating Mutation In Nanos3 In Two Sisters With Primary Ovarian Insufficiency.

Despite the increasing understanding of female reproduction, the molecular diagnosis of primary ovarian insufficiency (POI) is seldom obtained. The RNA-binding protein NANOS3 poses as an interesting candidate gene for POI since members of the Nanos family have an evolutionarily conserved function in germ cell development and maintenance by repressing apoptosis. We performed mutational analysis of NANOS3 in a cohort of 85 Brazilian women with familial or isolated POI, presenting with primary or secondary amenorrhea, and in ethnically-matched control women. A homozygous p.Glu120Lys mutation in NANOS3 was identified in two sisters with primary amenorrhea. The substituted amino acid is located within the second C2HC motif in the conserved zinc finger domain of NANOS3 and in silico molecular modelling suggests destabilization of protein-RNA interaction. In vitro analyses of apoptosis through flow cytometry and confocal microscopy show that NANOS3 capacity to prevent apoptosis was impaired by this mutation. The identification of an inactivating missense mutation in NANOS3 suggests a mechanism for POI involving increased primordial germ cells (PGCs) apoptosis during embryonic cell migration and highlights the importance of NANOS proteins in human ovarian biology.

A computational analysis of substrate binding strength by phosphorylase kinase and protein kinase A

Protein kinases exhibit various degrees of substrate specificity. The large number of different protein kinases in the eukaryotic proteomes makes it impractical to determine the specificity of each enzyme experimentally. To test if it were possible to discriminate potential substrates from non-substrates by simple computational techniques, we analysed the binding enthalpies of modelled enzyme-substrate complexes and attempted to correlate it with experimental enzyme kinetics measurements. The crystal structures of phosphorylase kinase and cAMP-dependent protein kinase were used to generate models of the enzyme with a series of known peptide substrates and non-substrates, and the approximate enthalpy of binding assessed following energy minimization. We show that the computed enthalpies do not correlate closely with kinetic measurements, but the method can distinguish good substrates from weak substrates and non-substrates. Copyright (C) 2002 John Wiley Sons, Ltd.

Computational tools for the study of allergens

Allergy is a major cause of morbidity worldwide. The number of characterized allergens and related information is increasing rapidly creating demands for advanced information storage, retrieval and analysis. Bioinformatics provides useful tools for analysing allergens and these are complementary to traditional laboratory techniques for the study of allergens. Specific applications include structural analysis of allergens, identification of B- and T-cell epitopes, assessment of allergenicity and cross-reactivity, and genome analysis. In this paper, the most important bioinformatic tools and methods with relevance to the study of allergy have been reviewed.