Compass: A shape-based machine learning tool for drug design

Summary Building predictive models for iterative drug design in the absence of a known target protein structure is an important challenge. We present a novel technique, Compass, that removes a major obstacle to accurate prediction by automatically selecting conformations and alignments of molecules without the benefit of a characterized active site. The technique combines explicit representation of molecular shape with neural network learning methods to produce highly predictive models, even across chemically distinct classes of molecules. We apply the method to predicting human perception of musk odor and show how the resulting models can provide graphical guidance for chemical modifications.     

The multi-copy simultaneous search methodology: a fundamental tool for structure-based drug design

Fragment-based ligand design approaches, such as the multi-copy simultaneous search (MCSS) methodology, have proven to be useful tools in the search for novel therapeutic compounds that bind pre-specified targets of known structure. MCSS offers a variety of advantages over more traditional high-throughput screening methods, and has been applied successfully to challenging targets. The methodology is quite general and can be used to construct functionality maps for proteins, DNA, and RNA. In this review, we describe the main aspects of the MCSS method and outline the general use of the methodology as a fundamental tool to guide the design of de novo lead compounds. We focus our discussion on the evaluation of MCSS results and the incorporation of protein flexibility into the methodology. In addition, we demonstrate on several specific examples how the information arising from the MCSS functionality maps has been successfully used to predict ligand binding to protein targets and RNA.     

In silico fragment-mapping method: a new tool for fragment-based/structure-based drug discovery

Here, we propose an in silico fragment-mapping method as a potential tool for fragment-based/structure-based drug discovery (FBDD/SBDD). For this method, we created a database named Canonical Subsite–Fragment DataBase (CSFDB) and developed a knowledge-based fragment-mapping program, Fsubsite. CSFDB consists of various pairs of subsite–fragments derived from X-ray crystal structures of known protein–ligand complexes. Using three-dimensional similarity-matching between subsites on one protein and another, Fsubsite compares the surface of a target protein with all subsites in CSFDB. When a local topography similar to the subsite is found on the surface, Fsubsite places a fragment combined with the subsite in CSFDB on the target protein. For validation purposes, we applied the method to the apo-structure of cyclin-dependent kinase 2 (CDK2) and identified four compounds containing three mapped fragments that existed in the list of known inhibitors of CDK2. Next, the utility of our fragment-mapping method for fragment-growing was examined on the complex structure of tRNA-guanine transglycosylase with a small ligand. Fsubsite mapped appropriate fragments on the same position as the binding ligand or in the vicinity of the ligand. Finally, a 3D-pharmacophore model was constructed from the fragments mapped on the apo-structure of heat shock protein 90-α (HSP90α). Then, 3D pharmacophore-based virtual screening was carried out using a commercially available compound database. The resultant hit compounds were very similar to a known ligand of HSP90α. As a result of these findings, this in silico fragment-mapping method seems to be a useful tool for computational FBDD and SBDD.     

In search of new lead compounds for trypanosomiasis drug design: A protein structure-based linked-fragment approach

Summary A modular method for pursuing structure-based inhibitor design in the framework of a design cycle is presented. The approach entails four stages: (1) a design pathway is defined in the three-dimensional structure of a target protein; (2) this pathway is divided into subregions; (3) complementary building blocks, also called fragments, are designed in each subregion; complementarity is defined in terms of shape, hydrophobicity, hydrogen bond properties and electrostatics; and (4) fragments from different subregions are linked into potential lead compounds. Stages (3) and (4) are qualitatively guided by force-field calculations. In addition, the designed fragments serve as entries for retrieving existing compounds from chemical databases. This linked-fragment approach has been applied in the design of potentially selective inhibitors of triosephosphate isomerase from Trypanosoma brucei, the causative agent of sleeping sickness.     

Discovery of a monomeric green fluorescent protein sensor for chloride by structure-guided bioinformatics

Chloride is an essential anion for all forms of life. Beyond electrolyte balance, an increasing body of evidence points to new roles for chloride in normal physiology and disease. Over the last two decades, this understanding has been advanced by chloride-sensitive fluorescent proteins for imaging applications in living cells. To our surprise, these sensors have primarily been engineered from the green fluorescent protein (GFP) found in the jellyfish Aequorea victoria. However, the GFP family has a rich sequence space that could already encode for new sensors with desired properties, thereby minimizing protein engineering efforts and accelerating biological applications. To efficiently sample this space, we present and validate a stepwise bioinformatics strategy focused first on the chloride binding pocket and second on a monomeric oligomerization state. Using this, we identified GFPxm163 from GFPxm found in the jellyfish Aequorea macrodactyla. In vitro characterization shows that the binding of chloride as well as bromide, iodide, and nitrate rapidly tunes the ground state chromophore equilibrium from the phenolate to the phenol state generating a pH-dependent, turn-off fluorescence response. Furthermore, live-cell fluorescence microscopy reveals that GFPxm163 provides a reversible, yet indirect readout of chloride transport via iodide exchange. With this demonstration, we anticipate that the pairing of bioinformatics with protein engineering methods will provide an efficient methodology to discover and design new chloride-sensitive fluorescent proteins for cellular applications.

We developed a workflow to identify and apply GFPxm163 as a new green fluorescent protein-based sensor for chloride.     

Coordinated design of cofactor and active site structures in development of new protein catalysts

New methods for the synthesis of artificial metalloenzymes are important for the construction of novel biocatalysts and biomaterials. Recently, we reported new methodology for the synthesis of artificial metalloenzymes by reconstituting apo-myoglobin with metal complexes (Ohashi, M. et al., Angew Chem., Int. Ed. 2003, 42, 1005-1008). However, it has been difficult to improve their reactivity, since their crystal structures were not available. In this article, we report the crystal structures of M(III)(Schiff base).apo-A71GMbs (M = Cr and Mn). The structures suggest that the position of the metal complex in apo-Mb is regulated by (i) noncovalent interaction between the ligand and surrounding peptides and (ii) the ligation of the metal ion to proximal histidine (His93). In addition, it is proposed that specific interactions of Ile107 with 3- and 3'-substituent groups on the salen ligand control the location of the Schiff base ligand in the active site. On the basis of these results, we have successfully controlled the enantioselectivity in the sulfoxidation of thioanisole by changing the size of substituents at the 3 and 3' positions. This is the first example of an enantioselective enzymatic reaction regulated by the design of metal complex in the protein active site.     

Force mode atomic force microscopy as a tool for protein folding studies

The advent of a new class of force microscopes designed specifically to “pull” biomolecules has allowed non-specialists to use force microscopy as a tool to study single-molecule protein unfolding. This powerful new technique has the potential to explore regions of the protein energy landscape that are not accessible in conventional bulk studies. It has the added advantage of allowing direct comparison with single-molecule simulation experiments. However, as with any new technique, there is currently no well described consensus for carrying out these experiments. Adoption of standard schemes of data selection and analysis will facilitate comparison of data from different laboratories and on different proteins. In this review, some guidelines and principles, which have been adopted by our laboratories, are suggested. The issues associated with collecting sufficient high quality data and the analysis of those data are discussed. In single-molecule studies, there is an added complication since an element of judgement has to be applied in selecting data to analyse; we propose criteria to make this process more objective. The principal sources of error are identified and standardised methods of selecting and analysing the data are proposed. The errors associated with the kinetic parameters obtained from such experiments are evaluated. The information that can be obtained from dynamic force experiments is compared, both quantitatively and qualitatively to that derived from conventional protein folding studies.     

pH-gradient ion-exchange chromatography: an analytical tool for design and optimization of protein separations

This work demonstrates that a highly linear, controllable and wide-ranged pH-gradient can be generated through an ion-exchange chromatography (IEC) column. Such a pH-gradient anion-exchange chromatography was evaluated with 17 model proteins and found that acidic (pI<6) and basic (pI>8) proteins elute roughly at their pI, whereas neutral proteins (pI 6-8) elute at pH 8-9 regardless their pI values. Because of the flat nature of protein titration curves from pH approximately 6 to approximately 9, neutral proteins indeed exhibit nearly zero net charge at pH approximately 9. The elution-pH in pH-gradient IEC or the titration curve, but not the pI, was identified as the key parameter for pH optimization of preparative IEC in a fast and rational way. The pH-gradient IEC was also applied and found to be an excellent analytical tool for the fractionation of crude protein mixtures.     

Incorporating protein flexibility in structure-based drug discovery: using HIV-1 protease as a test case

We have developed a receptor-based pharmacophore method which utilizes a collection of protein structures to account for inherent protein flexibility in structure-based drug design. Several procedures were systematically evaluated to derive the most general protocol for using multiple protein structures. Most notably, incorporating more protein flexibility improved the performance of the method. The pharmacophore models successfully discriminate known inhibitors from drug-like non-inhibitors. Furthermore, the models correctly identify the bound conformations of some ligands. We used unliganded HIV-1 protease to develop and validate this method. Drug design is always initiated with a protein-ligand structure, and such success with unbound protein structures is remarkable - particularly in the case of HIV-1 protease, which has a large conformational change upon binding. This technique holds the promise of successful computer-based drug design before bound crystal structures are even discovered, which can mean a jump-start of 1-3 years in tackling some medically relevant systems with computational methods.     

A new tool for the rational design of methylbiotin hosts

To study the binding mode of biotin related compounds with artificial hosts, we have developed a new tool to be used as a guide to test their behaviour prior to their synthesis. In that way, we have considered a set of 23 complexes comprising biotin and urea derivatives with synthetic hosts to develop a Partial Least Squares Cross-Validated (PLS-CV) model. The data, for such a model, are the binding constants (K_b) of each complex and the interaction energies (−E_min) calculated by molecular mechanics with AMBER and OPLS force fields. The predictive power of the model has been verified.     

Immunoaffinity capillary electrophoresis: a new versatile tool for determining protein biomarkers in inflammatory processes

Many diseases caused by inflammatory processes can progress to a chronic state causing deterioration in the quality of life and a poor prognosis for long-term survival. To address inflammatory diseases effectively, early detection and novel therapeutics are required. However, this can be challenging, in part because of the lack of early predictive biomarkers and the limited availability of adequate technologies capable of the identification/characterization of key predictive biomarkers present in biological materials, especially those found at picomolar concentrations and below. This review highlights the need for state-of-the art methodologies, with high-sensitivity and high-throughput capabilities, for determination of multiple biomarkers. Although many new biomarkers have been discovered recently, existing technology has failed to successfully bring this advancement to the patient's bedside. We present an overview of the various advances available today to extend the discovery of predictive biomarkers of inflammatory diseases; in particular, we review the technology of immunoaffinity capillary electrophoresis (IACE), which combines the use of antibodies as highly selective capture agents with the high resolving power of capillary electrophoresis. This two-dimensional hybrid technology permits the quantification and characterization of several protein biomarkers simultaneously, including subtle structural changes such as variants, isoforms, peptide fragments, and post-translational modifications. Furthermore, the results are rapid, sensitive, can be performed at a relatively low cost, without the introduction of false positive or false negative data. The IACE instrumentation can have relevance to medical, pharmaceutical, environmental, military, cultural heritage (authenticity of art work), forensic science, industrial and research fields, and in particular as a point-of-care biomarker analyzer in translational medicine.     

GREEN: A program package for docking studies in rational drug design

Summary A program package, GREEN, has been developed that enables docking studies between ligand molecules and a protein molecule. Based on the structure of the protein molecule, the physical and chemical environment of the ligand-binding site is expressed as three-dimensional grid-point data. The grid-point data are used for the real-time evaluation of the protein-ligand interaction energy, as well as for the graphical representation of the binding-site environment. The interactive docking operation is facilitated by various built-in functions, such as energy minimization, energy contribution analysis and logging of the manipulation trajectory. Interactive modeling functions are incorporated for designing new ligand molecules while considering the binding-site environment and the protein-ligand interaction. As an example of the application of GREEN, a docking study is presented on the complex between trypsin and a synthetic trypsin inhibitor. The program package will be useful for rational drug design, based on the 3D structure of the target protein.     

Aroma WaterLOGSY: a fast and sensitive screening tool for drug discovery

One-dimensional NMR spectroscopy has proven to be a powerful technique for screening compound libraries in drug discovery. We report a novel water ligand-observed gradient spectroscopy (WaterLOGSY) pulse sequence, named Aroma WaterLOGSY, that selectively detects aromatic WaterLOGSY signals from compounds or ligands. In the Aroma WaterLOGSY, water magnetization is untouched after water excitation and utilizes the whole period of the remaining pulse sequence to relax back to the +z direction. Due to the phase cycling design, the water magnetization is allowed to relax for the period of two full scans before it gets inverted again. Therefore, the recycle delay can be significantly shortened. Within similar experimental time, Aroma WaterLOGSY shows approximately two times higher sensitivity than the standard scheme. This method also allows the use of non-deuterated reagents, thereby accelerating experimental set-up time for ligand-binding studies.     

Preparation and in vivo studies of a new drug delivery system

Polyhexylcyanoacrylate nanoparticles have been prepared with vincamine as the model drug. These particles had an average size of 200 nm and adsorbed approximately 435 of vincamine. The adsorption of vincamine to nanoparticles modified the distribution of vincamine in tissues. After iv injection the distribution volumes were increased in comparison with an aqueous solution of drug. In comparison with an aqueous solution of drug, the absolute bioavailability of vincamine was also increased after an oral administration of nanoparticles.     

Rhodotorula aurantiaca penicillin V acylase: Active site characterization and fluorometric studies

Penicillin V acylase (PVA), a member of newly evolved Ntn-hydrolase superfamily, is a pharmaceutically important enzyme to produce 6-aminopenicillanic acid. Active site characterization of recently purified monomeric PVA from Rhodotorula aurantiaca (Ra-PVA), the yeast source, showed the involvement of serine and tryptophan in the enzyme activity. Modification of the protein with serine and tryptophan specific reagents such as PMSF and NBS showed partial loss of PVA activity and substrate protection. Ra-PVA was found to be a multi-tryptophan protein exhibiting one tryptophan, in native and, four in its denatured condition. Various solute quenchers and substrate were used to probe the microenvironment of the putative reactive tryptophan through fluorescence quenching. The results obtained indicate that the tryptophan residues of Ra-PVA were largely buried in hydrophobic core of the protein matrix. Quenching of the fluorescence by acrylamide was collisional. Acrylamide was the most effective quencher amongst all the used quenchers, which quenched 71.6% of the total intrinsic fluorescence of the protein, at a very less final concentration of 0.1 M. Surface tryptophan residues were found to have predominantly more electropositively charged amino acids around them, however differentially accessible for ionic quenchers. Denaturation led to shift in λ_max from 336, in native state, to 357 nm and more exposed to the solvent, consequently increase in fluorescence quenching with all quenchers. This is an attempt towards the conformational studies of Ra-PVA.     

Conditional protein splicing: a new tool to control protein structure and function in vitro and in vivo

Protein splicing is a naturally occurring process in which an intervening intein domain excises itself out of a precursor polypeptide in an autocatalytic fashion with concomitant linkage of the two flanking extein sequences by a native peptide bond. We have recently reported an engineered split VMA intein whose splicing activity in trans between two polypeptides can be triggered by the small molecule rapamycin. In this report, we show that this conditional protein splicing (CPS) system can be used in mammalian cells. Two model constructs harboring maltose-binding protein (MBP) and a His-tag as exteins were expressed from a constitutive promoter after transient transfection. The splicing product MBP-His was detected by Western blotting and immunoprecipitation in cells treated with rapamycin or a nontoxic analogue thereof. No background splicing in the absence of the small-molecule inducer was observed over a 24-h time course. Product formation could be detected within 10 min of addition of rapamycin, indicating the advantage of the posttranslational nature of CPS for quick responses. The level of protein splicing was dose dependent and could be competitively attenuated with the small molecule ascomycin. In related studies, the geometric flexibility of the CPS components was investigated with a series of purified proteins. The FKBP and FRB domains, which are dimerized by rapamycin and thereby induce the reconstitution of the split intein, were fused to the extein sequences of the split intein halves. CPS was still triggered by rapamycin when FKBP and FRB occupied one or both of the extein positions. This finding suggests yet further applications of CPS in the area of proteomics. In summary, CPS holds great promise to become a powerful new tool to control protein structure and function in vitro and in living cells.