Formation of neutral molecules of potential stellar interest by neutralisation of negative ions in a mass spectrometer - The application of experiment and molecular modelling in concert

This paper is a modified version of a lecture which describes the synthesis, structure and reactivity of some neutral molecules of stellar significance. The neutrals are formed in the collision cell of a mass spectrometer following vertical Franck-Condon one electron oxidation of anions of known bond connectivity. Neutrals are characterised by conversion to positive ions and by extensive theoretical studies at the CCSD(T)/aug-cc-pVDZ//B3LYP/6-31G(d) level of theory. Four systems are considered in detail, viz (i) the formation of linear C-4 and its conversion to the rhombus C-4, (ii) linear C-5 and the atom scrambling of this system when energised, (iii) the stable cumulene oxide CCCCCO, and (iv) the elusive species O2C-CO. This paper is not intended to be a review of interstellar chemistry: examples are selected from our own work in this area. (C) 2002 Elsevier Science Inc. All rights reserved.

Molecular Modelling of Estrogen-Producing Enzymes CYP450 Aromatase and 17beta-Hydroxysteroid Dehydrogenase Type 1

Breast cancer is the most common cancer in women in the western countries. Approximately two-thirds of breast cancer tumours are hormone dependent, requiring estrogens to grow. Estrogens are formed in the human body via a multistep route starting from cholesterol. The final steps in the biosynthesis include the CYP450 aromatase enzyme, converting the male hormones androgens (preferred substrate androstenedione ASD) into estrogens(estrone E1), and the 17beta-HSD1 enzyme, converting the biologically less active E1 into the active hormone 17beta-hydroxyestradiol E2. E2 is bound to the nuclear estrogen receptors causing a cascade of biochemical reactions leading to cell proliferation in normal tissue, and to tumour growth in cancer tissue. Aromatase and 17beta-HSD1 are expressed in or near the breast tumour, locally providing the tissue with estrogens. One approach in treating hormone dependent breast tumours is to block the local estrogen production by inhibiting these two enzymes. Aromatase inhibitors are already on the market in treating breast cancer, despite the lack of an experimentally solved structure. The structure of 17beta-HSD1, on the other hand, has been solved, but no commercial drugs have emerged from the drug discovery projects reported in the literature. Computer-assisted molecular modelling is an invaluable tool in modern drug design projects. Modelling techniques can be used to generate a model of the target protein and to design novel inhibitors for them even if the target protein structure is unknown. Molecular modelling has applications in predicting the activities of theoretical inhibitors and in finding possible active inhibitors from a compound database. Inhibitor binding at atomic level can also be studied with molecular modelling. To clarify the interactions between the aromatase enzyme and its substrate and inhibitors, we generated a homology model based on a mammalian CYP450 enzyme, rabbit progesterone 21-hydroxylase CYP2C5. The model was carefully validated using molecular dynamics simulations (MDS) with and without the natural substrate ASD. Binding orientation of the inhibitors was based on the hypothesis that the inhibitors coordinate to the heme iron, and were studied using MDS. The inhibitors were dietary phytoestrogens, which have been shown to reduce the risk for breast cancer. To further validate the model, the interactions of a commercial breast cancer drug were studied with MDS and ligand–protein docking. In the case of 17beta-HSD1, a 3D QSAR model was generated on the basis of MDS of an enzyme complex with active inhibitor and ligand–protein docking, employing a compound library synthesised in our laboratory. Furthermore, four pharmacophore hypotheses with and without a bound substrate or an inhibitor were developed and used in screening a commercial database of drug-like compounds. The homology model of aromatase showed stable behaviour in MDS and was capable of explaining most of the results from mutagenesis studies. We were able to identify the active site residues contributing to the inhibitor binding, and explain differences in coordination geometry corresponding to the inhibitory activity. Interactions between the inhibitors and aromatase were in agreement with the mutagenesis studies reported for aromatase. Simulations of 17beta-HSD1 with inhibitors revealed an inhibitor binding mode with hydrogen bond interactions previously not reported, and a hydrophobic pocket capable of accommodating a bulky side chain. Pharmacophore hypothesis generation, followed by virtual screening, was able to identify several compounds that can be used in lead compound generation. The visualisation of the interaction fields from the QSAR model and the pharmacophores provided us with novel ideas for inhibitor development in our drug discovery project.

New structural insights and molecular-modelling studies of 4-methyl-5-beta-hydroxyethylthiazole kinase from Pyrococcus horikoshii OT3 (PhThiK)

4-Methyl-5-beta-hydroxyethylthiazole kinase (ThiK) catalyses the phosphorylation of the hydroxyl group of 4-methyl-5-beta-hydroxyethylthiazole. This work reports the first crystal structure of an archaeal ThiK: that from Pyrococcus horikoshii OT3 (PhThiK) at 1.85 angstrom resolution with a phosphate ion occupying the position of the beta-phosphate of the nucleotide. The topology of this enzyme shows the typical ribokinase fold of an alpha/beta protein. The overall structure of PhThiK is similar to those of Bacillus subtilis ThiK (BsThiK) and Enterococcus faecalis V583 ThiK (EfThiK). Sequence analysis of ThiK enzymes from various sources indicated that three-quarters of the residues involved in interfacial regions are conserved. It also revealed that the amino-acid residues in the nucleotide-binding, magnesium ion-binding and substrate-binding sites are conserved. Binding of the nucleotide and substrate to the ThiK enzyme do not influence the quaternary association (trimer) as revealed by the crystal structure of PhThiK.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Molecular modelling of saccharide-lipid interactions

Certain saccharides, including trehalose, sucrose and glucose, stabilize lipid bilayers against dehydration. It has been suggested that these saccharides replace waters of hydration as the system is dried, thereby maintaining the headgroups at their hydrated spacing. The lipid acyl chains consequently have sufficient free volume to remain in the liquid crystallines state, and the processes that disrupt membrane integrity are inhibited. Initial molecular graphic investigations of a model trehalose/DMPC system supported this idea (Chandrasekhar, I. and Gaber, B.P. (1988) J. Biomol. Stereodyn, 5, 1163–1171). We have extended these studies to glucose and sucrose. A set of AMBER potential parameters has been established that reproduce simple saccharide conformations, including the anomeric effect. Extensive energy minimizations have been conducted on all three systems. The saccharide-lipid interaction energies become less stable in the order trehalose

Molecular modelling of ABC transporters: from ATP hydrolysis to substrate transport

Dissertação para a obtenção de grau de doutor em Bioquímica pelo Instituto de Tecnologia Química e Biológica da Universidade Nova de Lisboa

Metal complexes of a tetraazacyclophane: solution and molecular modelling studies

An alternative synthetic approach to yield the compound 2,3,5,6,8,9,11,14-octahydrobenzo[1][ 1,4,7,10]tetraazacyclotetradecine (bz[14]N-4) is presented. The protonation constants of bz[14]N-4 and the stability constants of its complexes with Ni2+, Cu2+, Zn2+, Cd2+, Pb2+ were determined in H2O at 25degreesC with ionic strength 0.10 mol dm(-3) in KNO3 and they were compared with structurally related macrocycles cyclam (1,4,8,11-tetraazacyclotetradecane) and cyclen (1,4,7,10-tetraazacyclododecane). These studies indicate that only 1 : 1 ( M : L) species are formed in solution, and the ligand exhibits a high affinity for larger ions such as Cd2+ and Pb2+. The X-ray study of [bz[14]N4H3](3+) shows that an inclusion compound with a chloride counter-anion is formed through NH...Cl hydrogen bonds. Spectroscopic data in solution ( electronic and NMR spectra) showed that the macrocycle adopts a planar arrangement upon metal complexation. Molecular mechanics calculations reveal that in spite of the presence of the benzene ring in the macrocyclic framework this ligand can encapsulate metal ions with different stereo-electronic sizes in square planar arrangements. Our results indicate that the presence of the benzene ring in the backbone of the bz[14]N-4 confers a coordination behaviour intermediate between that of cyclam and cyclen.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Application of Molecular Modelling studies to the search of novel approaches for neuroprotection and stem cells applications

Il lavoro di ricerca presentato in questa tesi di dottorato riguarda l'applicazione di studi di modellistica molecolare per l'individuazione di nuovi approcci farmacologici nel campo della neuroprotezione e del controllo della proliferazione di cellule staminali. Durante il mio dottorato di ricerca, mi sono concentrata sullo studio del sistema degli endocannabinoidi come target per lo sviluppo di nuovi trattamenti neuroprotettivi. In particolare, la mia ricerca ha avuto come obiettivo la modulazione dei livelli di 2-arachidonilglicerolo e arachidonil-etanolamide tramite l'inibizione degli enzimi MGL (monoglyceride lipase) e FAAH (fatty acid amide hydrolase). Il mio progetto di ricerca comprende anche studi di modellistica molecolare per l'individuazione di piccole molecole in grado di inibire il complesso proteina-proteina YAP-TEAD. Tale complesso, coinvolto nei sistemi di regolazione della proliferazione cellulare, rappresenta un target di cruciale importanza nel controllo della proliferazione e differenziazione di cellule staminali e, al tempo stesso, nel controllo dell'espansione tumorale

Molecular modelling of human CYP1B1 substrate interactions and investigation of allelic variant effects on metabolism

Molecular modelling of human CYP1B1 based on homology with the mammalian P450, CYP2C5, of known three-dimensional structure is reported. The enzyme model has been used to investigate the likely mode of binding for selected CYP1B1 substrates, particularly with regard to the possible effects of allelic variants of CYP1B1 on metabolism. In general, it appears that the CYP1B1 model is consistent with known substrate selectivity for the enzyme, and the sites of metabolism can be rationalized in terms of specific contacts with key amino acid residues within the CYP1B1 heme locus. Further-more, a mode of binding interaction for the inhibitor, a-naphthoflavone, is presented which accords with currently available information. The current paper shows that a combination of molecular modelling and experimental determinations on the substrate metabolism for CYP1B1 allelic variants can aid in the understanding of structure-function relationships within P450 enzymes. (C) 2003 Elsevier Science Ireland Ltd. All rights reserved.

Molecular modelling assisted design and synthesis of cyclothialidine derivatives as DNA gyrase inhibitors

Since cyclothialidine was discovered as the most active DNA gyrase inhibitor in 1994, enormous efforts have been devoted to make it into a commercial medicine by a number of pharmaceutical companies and research groups worldwide. However, no serious breakthrough has been made up to now. An essential problem involved with cyclothialidine is that though it demonstrated the potent inhibition of DNA gyrase, it showed little activity against bacteria. This probably is attributable to its inability to penetrate bacterial cell walls and membranes. We applied the TSAR programme to generate a QSAR equation to the gram-negative organisms. In that equation, LogP is profoundly indicated as the key factor influencing the cyclothialidine activity against bacteria. However, the synthesized new analogues have failed to prove that. In the structure based drug design stage, we designed a group of open chain cyclothialidine derivatives by applying the SPROUT programme and completed the syntheses. Improved activity is found in a few analogues and a 3D pharmacophore of the DNA gyrase B is proposed to lead to synthesis of the new derivatives for development of potent antibiotics.