Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

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.

Relieving the first bottleneck in the drug discovery pipeline:using array technologies to rationalize membrane protein production

The slow down in the drug discovery pipeline is, in part, owing to a lack of structural and functional information available for new drug targets. Membrane proteins, the targets of well over 50% of marketed pharmaceuticals, present a particular challenge. As they are not naturally abundant, they must be produced recombinantly for the structural biology that is a prerequisite to structure-based drug design. Unfortunately, however, obtaining high yields of functional, recombinant membrane proteins remains a major bottleneck in contemporary bioscience. While repeated rounds of trial-and-error optimization have not (and cannot) reveal mechanistic details of the biology of recombinant protein production, examination of the host response has provided new insights. To this end, we published an early transcriptome analysis that identified genes implicated in high-yielding yeast cell factories, which has enabled the engineering of improved production strains. These advances offer hope that the bottleneck of membrane protein production can be relieved rationally.

Catalytically active Dengue virus NS3 protease forms aggregates that are separable by size exclusion chromatography

An active form of the Dengue virus protease NS3 (CF40.Gly.NS3pro) was expressed in Escherichia coli. This construct consists of a critical 40 amino acid cofactor domain from NS2B fused to the N-terminal 184 amino acid protease domain of NS3 via a flexible, covalent linker (Gly(4)SerGly(4)). The recombinantly produced protein is soluble and has a hexa-histidine tag engineered at the N-terminus for ease of purification using metal affinity chromatography. However, the presence of lower molecular weight impurities after affinity chromatography indicated the need for additional purification steps. The consistent appearance of these impurities suggested that they may be the products of proteolysis and/or auto-proteolysis. The latter possibility was subsequently excluded by the observation of the same impurities in a purified, catalytically inactive form of the recombinant protease (CF40.Gly.NS3pro.SA). Further analysis indicated that these impurities may represent premature translation termination products. Regardless of their origin, they were shown to form various sized aggregates with full-length CF40.Gly.NS3pro that can be separated by size exclusion chromatography, yielding fractions of active protease of sufficient purity for crystallisation trials. The ultimate goal of these studies is to obtain a crystal structure of a catalytically active form of the Dengue virus NS3 protease for structure-based drug design. (C) 2002 Elsevier Science (USA). All rights reserved.

X-ray structures of Sap1 and Sap5: structural comparison of the secreted aspartic proteinases from Candida albicans.

Proteolytic activity is an important virulence factor for Candida albicans (C. albicans). It is attributed to the family of the secreted aspartic proteinases (Saps) from C. albicans with a minimum of 10 members. Saps show controlled expression and regulation for the individual stages of the infection process. Distinct isoenzymes can be responsible for adherence and tissue damage of local infections, while others cause systemic diseases. Earlier, only the structures of Sap2 and Sap3 were known. In our research, we have now succeeded in solving the X-ray crystal structures of the apoenzyme of Sap1 and Sap5 in complex with pepstatin A at 2.05 and 2.5 A resolution, respectively. With the structure of Sap1, we have completed the set of structures of isoenzyme subgroup Sap1-3. Of subgroup Sap4-6, the structure of the enzyme Sap5 is the first structure that has been described up to now. This facilitates comparison of structural details as well as inhibitor binding modes among the different subgroup members. Structural analysis reveals a highly conserved overall secondary structure of Sap1-3 and Sap5. However, Sap5 clearly differs from Sap1-3 by its electrostatic overall charge as well as through structural conformation of its entrance to the active site cleft. Design of inhibitors specific for Sap5 should concentrate on the S4 and S3 pockets, which significantly differ from Sap1-3 in size and electrostatic charge. Both Sap1 and Sap5 seem to play a major part in superficial Candida infections. Determination of the isoenzymes' structures can contribute to the development of new Sap-specific inhibitors for the treatment of superficial infections with a structure-based drug design program.

Immersive Projection for biochemical pattern matching

In all biological processes, protein molecules and other small molecules interact to function and form transient macromolecular complexes. This interaction of two or more molecules can be described by a docking event. Docking is an important phase for structure-based drug design strategies, as it can be used as a method to simulate protein-ligand interactions. Various docking programs exist that allow automated docking, but most of them have limited visualization and user interaction. It would be advantageous if scientists could visualize the molecules participating in the docking process, manipulate their structures and manually dock them before submitting the new conformations to an automated docking process in an immersive environment, which can help stimulate the design/docking process. This also could greatly reduce docking time and resources. To achieve this, we propose a new virtual modelling/docking program, whereby the advantages of virtual modelling programs and the efficiency of the algorithms in existing docking programs will be merged.

In silico screening of HIV-1 non-nucleoside reverse transcriptase and protease inhibitors

Two targets, reverse transcriptase (RT) and protease from HIV-1, were used during the past two decades to the discovery of non-nucleoside reverse transcriptase inhibitors (NNRTI) and protease inhibitors (PI) that belong to the arsenal of the antiretroviral therapy. Herein these enzymes were chosen as templates for conducting a computer-aided ligand design. Ligand and structure-based drug designs were the starting points to select compounds from a database bearing more than five million compounds by means of cheminformatic tools. New promising lead structures are retrieved from the database, which are open to acquisition and test. Classes of molecules already described as NNRTI or PI in the literature also came out and were useful to prove the reliability of the workflow, and thus validating the work carried out so far. (c) 2007 Elsevier Masson SAS. All rights reserved.

MPO Inhibitors Selected by Virtual Screening

The hemeprotein myeloperoxidase (MPO) participates in innate immune defense through its ability to generate potent microbicidal oxidants. However, these oxidants are also key mediators of the tissue damage associated with many inflammatory diseases. Thus, there is considerable interest in developing therapeutically useful MPO inhibitors. Here, we used structure-based drug design (SBDD) and ligand-based drug design (LBDD) to select for potentially new and selective MPO inhibitors. A pharmacophore model was developed based on the crystal structure of human MPO in complex with salicylhydroxamic acid (SHA), a known inhibitor of the enzyme. The pharmacophore model was used to screen the ZINC database for potential ligands, which were further filtered on the basis of their physical-chemical properties and docking score. The filtered compounds were visually inspected, and nine were purchased for experimental studies. Surprisingly, almost all of the selected compounds belonged to the aromatic hydrazide class, which had been previously described as MPO inhibitors. The compounds selected by virtual screening were shown to inhibit the chlorinating activity of MPO; the top four compounds displayed IC(50) values ranging from 1.0 to 2.8 mM. MPO inactivation by the most effective compound was shown to be irreversible. Overall, our results show that SBDD and LBDD may be useful for the rational development of new MPO inhibitors.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Tuberculosis (TB) is one of the most common infectious diseases known to man and responsible for millions of human deaths in the world. The increasing incidence of TB in developing countries, the proliferation of multidrug resistant strains, and the absence of resources for treatment have highlighted the need of developing new drugs against TB. The shikimate pathway leads to the biosynthesis of chorismate, a precursor of aromatic amino acids. This pathway is absent from mammals and shown to be essential for the survival of Mycobacterium tuberculosis, the causative agent of TB. Accordingly, enzymes of aromatic amino acid biosynthesis pathway represent promising targets for structure-based drug design. The first reaction in phenylalanine biosynthesis involves the conversion of chorismate to prephenate, catalyzed by chorismate mutase. The second reaction is catalyzed by prephenate dehydratase (PDT) and involves decarboxylation and dehydratation of prephenate to form phenylpyruvate, the precursor of phenylalanine. Here, we describe utilization of different techniques to infer the structure of M. tuberculosis PDT (MtbPDT) in solution. Small angle X-ray scattering and ultracentrifugation analysis showed that the protein oligomeric state is a tetramer and MtbPDT is a flat disk protein. Bioinformatics tools were used to infer the structure of MtbPDT A molecular model for MtbPDT is presented and molecular dynamics simulations indicate that MtbPDT i.s stable. Experimental and molecular modeling results were in agreement and provide evidence for a tetrameric state of MtbPDT in solution.

Molecular models of protein targets from Mycobacterium tuberculosis

Structural characterization of enzymes that belong to microbial metabolic pathways is very important for structure-based drug design since some of these proteins may be present in the bacterial genome, but absent in humans. Thus, metabolic pathways became potential targets for drug design. The motivation of this work is the fact that Mycobacterium tuberculosis is the cause of the deaths of millions of people in the world, so that the structural characterization of protein targets to propose new drugs has become essential. DBMODELING is a relational database, created to highlight the importance of methods of molecular modeling applied to the Mycobacterium tuberculosis genome with the aim of proposing protein-ligand docking analysis. There are currently more than 300 models for proteins from Mycobacterium tuberculosis genome in the database. The database contains a detailed description of the reaction catalyzed by each enzyme and their atomic coordinates. Information about structures, a tool for animated gif image, a table with a specification of the metabolic pathway, modeled protein, inputs used in modeling, and analysis methods used in this project are available in the database for download. The search tool can be used for reseachers to find specific pathways or enzymes.

Docking and small angle X-ray scattering studies of purine nucleoside phosphorylase

Docking simulations have been used to assess protein complexes with some success. Small angle X-ray scattering (SAXS) is a well-established technique to investigate protein spatial configuration. This work describes the integration of geometric docking with SAXS to investigate the quaternary structure of recombinant human purine nucleoside phosphorylase (PNP). This enzyme catalyzes the reversible phosphorolysis of N-ribosidic bonds of purine nucleosides and deoxynucleosides. A genetic deficiency due to mutations in the gene encoding for PNP causes gradual decrease in T-cell immunity. Inappropriate activation of T-cells has been implicated in several clinically relevant human conditions such as transplant rejection, rheumatoid arthritis, lupus, and T-cell lymphomas. PNP is therefore a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. The present analysis confirms the trimeric structure observed in the crystal. The potential application of the present procedure to other systems is discussed. (C) 2003 Elsevier B.V. All rights reserved.

Determining the structural basis for specificity of ligands using crystallographic screening

Crystallographic screening has been used to identify new inhibitors for potential target for drug development. Here, we describe the application of the crystallographic screening to assess the structural basis of specificity of ligands against a protein target. The method is efficient and results in detailed crystallographic information. The utility of the method is demonstrated in the study of the structural basis for specificity of ligands for human purine nucleoside phosphorylase (PNP). Purine nucleoside phosphorylase catalyzes the phosphorolysis of the N-ribosidic bonds of purine nucleosides and deoxynucleosides. This enzyme is a target for inhibitor development aiming at T-cell immune response modulation and has been submitted to extensive structure-based drug design. This methodology may help in the future development of a new generation of PNP inhibitors.

Molecular models of NS3 protease variants of the Hepatitis C virus

Background: Hepatitis C virus (HCV) currently infects approximately three percent of the world population. In view of the lack of vaccines against HCV, there is an urgent need for an efficient treatment of the disease by an effective antiviral drug. Rational drug design has not been the primary way for discovering major therapeutics. Nevertheless, there are reports of success in the development of inhibitor using a structure-based approach. One of the possible targets for drug development against HCV is the NS3 protease variants. Based on the three-dimensional structure of these variants we expect to identify new NS3 protease inhibitors. In order to speed up the modeling process all NS3 protease variant models were generated in a Beowulf cluster. The potential of the structural bioinformatics for development of new antiviral drugs is discussed.Results: the atomic coordinates of crystallographic structure 1CU1 and 1DY9 were used as starting model for modeling of the NS3 protease variant structures. The NS3 protease variant structures are composed of six subdomains, which occur in sequence along the polypeptide chain. The protease domain exhibits the dual beta-barrel fold that is common among members of the chymotrypsin serine protease family. The helicase domain contains two structurally related beta-alpha-beta subdomains and a third subdomain of seven helices and three short beta strands. The latter domain is usually referred to as the helicase alpha-helical subdomain. The rmsd value of bond lengths and bond angles, the average G-factor and Verify 3D values are presented for NS3 protease variant structures.Conclusions: This project increases the certainty that homology modeling is an useful tool in structural biology and that it can be very valuable in annotating genome sequence information and contributing to structural and functional genomics from virus. The structural models will be used to guide future efforts in the structure-based drug design of a new generation of NS3 protease variants inhibitors. All models in the database are publicly accessible via our interactive website, providing us with large amount of structural models for use in protein-ligand docking analysis.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Crystallization and preliminary X-ray diffraction analysis of prephenate dehydratase from Mycobacterium tuberculosis H37Rv

Tuberculosis remains the leading cause of mortality arising from a bacterial pathogen ( Mycobacterium tuberculosis). There is an urgent need for the development of new antimycobacterial agents. The aromatic amino-acid pathway is essential for the survival of this pathogen and represents a target for structure-based drug design. Accordingly, the M. tuberculosis prephenate dehydratase has been cloned, expressed, purified and crystallized by the hanging-drop vapour-diffusion method using PEG 400 as a precipitant. The crystal belongs to the orthorhombic space group I222 or I2(1)2(1)2(1), with unit-cell parameters a = 98.26, b = 133.22, c = 225.01 angstrom, and contains four molecules in the asymmetric unit. A complete data set was collected to 3.2 angstrom resolution using a synchrotron-radiation source.