Structure-based discovery of inhibitors of the YycG histidine kinase: New chemical leads to combat Staphylococcus epidermidis infections

Background

Coagulase-negative Staphylococcus epidermidis has become a major frequent cause of infections in relation to the use of implanted medical devices. The pathogenicity of S. epidermidis has been attributed to its capacity to form biofilms on surfaces of medical devices, which greatly increases its resistance to many conventional antibiotics and often results in chronic infection. It has an urgent need to design novel antibiotics against staphylococci infections, especially those can kill cells embedded in biofilm.

Results

In this report, a series of novel inhibitors of the histidine kinase (HK) YycG protein of S. epidermidis were discovered first using structure-based virtual screening (SBVS) from a small molecular lead-compound library, followed by experimental validation. Of the 76 candidates derived by SBVS targeting of the homolog model of the YycG HATPase_c domain of S. epidermidis, seven compounds displayed significant activity in inhibiting S. epidermidis growth. Furthermore, five of them displayed bactericidal effects on both planktonic and biofilm cells of S. epidermidis. Except for one, the compounds were found to bind to the YycG protein and to inhibit its auto-phosphorylation in vitro, indicating that they are potential inhibitors of the YycG/YycF two-component system (TCS), which is essential in S. epidermidis. Importantly, all these compounds did not affect the stability of mammalian cells nor hemolytic activities at the concentrations used in our study.

Conclusion

These novel inhibitors of YycG histidine kinase thus are of potential value as leads for developing new antibiotics against infecting staphylococci. The structure-based virtual screening (SBVS) technology can be widely used in screening potential inhibitors of other bacterial TCSs, since it is more rapid and efficacious than traditional screening technology.     

Fast identification of folded human protein domains expressed in <Emphasis Type="Italic">E. coli</Emphasis> suitable for structural analysis

Background  

High-throughput protein structure analysis of individual protein domains requires analysis of large numbers of expression clones to identify suitable constructs for structure determination. For this purpose, methods need to be implemented for fast and reliable screening of the expressed proteins as early as possible in the overall process from cloning to structure determination.     

Protein kinase substrate identification on functional protein arrays

Background  

Over the last decade, kinases have emerged as attractive therapeutic targets for a number of different diseases, and numerous high throughput screening efforts in the pharmaceutical community are directed towards discovery of compounds that regulate kinase function. The emerging utility of systems biology approaches has necessitated the development of multiplex tools suitable for proteomic-scale experiments to replace lower throughput technologies such as mass spectroscopy for the study of protein phosphorylation. Recently, a new approach for identifying substrates of protein kinases has applied the miniaturized format of functional protein arrays to characterize phosphorylation for thousands of candidate protein substrates in a single experiment. This method involves the addition of protein kinases in solution to arrays of immobilized proteins to identify substrates using highly sensitive radioactive detection and hit identification algorithms.     

Scoring functions and enrichment: a case study on Hsp90

Background  

The need for fast and accurate scoring functions has been driven by the increased use of in silico virtual screening twinned with high-throughput screening as a method to rapidly identify potential candidates in the early stages of drug development. We examine the ability of some the most common scoring functions (GOLD, ChemScore, DOCK, PMF, BLEEP and Consensus) to discriminate correctly and efficiently between active and non-active compounds among a library of ~3,600 diverse decoy compounds in a virtual screening experiment against heat shock protein 90 (Hsp90).     

Ligand scaffold hopping combining 3D maximal substructure search and molecular similarity

Background  

Virtual screening methods are now well established as effective to identify hit and lead candidates and are fully integrated in most drug discovery programs. Ligand-based approaches make use of physico-chemical, structural and energetics properties of known active compounds to search large chemical libraries for related and novel chemotypes. While 2D-similarity search tools are known to be fast and efficient, the use of 3D-similarity search methods can be very valuable to many research projects as integration of "3D knowledge" can facilitate the identification of not only related molecules but also of chemicals possessing distant scaffolds as compared to the query and therefore be more inclined to scaffolds hopping. To date, very few methods performing this task are easily available to the scientific community.     

Representation of target-bound drugs by computed conformers: implications for conformational libraries

Background  

The increasing number of known protein structures provides valuable information about pharmaceutical targets. Drug binding sites are identifiable and suitable lead compounds can be proposed. The flexibility of ligands is a critical point for the selection of potential drugs. Since computed 3D structures of millions of compounds are available, the knowledge of their binding conformations would be a great benefit for the development of efficient screening methods.     

Systematic comparison of SCOP and CATH: a new gold standard for protein structure analysis

Background  

SCOP and CATH are widely used as gold standards to benchmark novel protein structure comparison methods as well as to train machine learning approaches for protein structure classification and prediction. The two hierarchies result from different protocols which may result in differing classifications of the same protein. Ignoring such differences leads to problems when being used to train or benchmark automatic structure classification methods. Here, we propose a method to compare SCOP and CATH in detail and discuss possible applications of this analysis.     

MS-DOCK: Accurate multiple conformation generator and rigid docking protocol for multi-step virtual ligand screening

Background  

The number of protein targets with a known or predicted tri-dimensional structure and of drug-like chemical compounds is growing rapidly and so is the need for new therapeutic compounds or chemical probes. Performing flexible structure-based virtual screening computations on thousands of targets with millions of molecules is intractable to most laboratories nor indeed desirable. Since shape complementarity is of primary importance for most protein-ligand interactions, we have developed a tool/protocol based on rigid-body docking to select compounds that fit well into binding sites.     

Identification of type I and type II inhibitors of c-Yes kinase using in silico and experimental techniques

c-Yes kinase is considered as one of the attractive targets for anti-cancer drug design. The DFG (Asp-Phe-Gly) motif present in most of the kinases will adopt active and inactive conformations, known as DFG-in and DFG-out and their inhibitors are classified into type I and type II, respectively. In the present study, two screening protocols were followed for identification of c-Yes kinase inhibitors. (i) Structure-based virtual screening (SBVS) and (ii) Structure-based (SB) and Pharmacophore-based (PB) tandem screening. In SBVS, the c-Yes kinase structure was obtained from homology modeling and seven ensembles with different active site scaffolds through molecular dynamics (MD) simulations. For SB-PB tandem screening, we modeled ligand bound active and inactive conformations. Physicochemical properties of inhibitors of Src kinase family and c-Yes kinase were used to prepare target focused libraries for screenings. Our screening procedure along with docking showed 520 probable hits in SBVS and tandem screening (120 and 400, respectively). Out of 5000 compounds identified from different computational methods, 2410 were examined using kinase inhibition assays. It includes 266 compounds (5.32%) identified from our method. We observed that 14 compounds (12%) are identified by the present method out of 168 that showed > 30% inhibition. Among them, three compounds are novel, unique, and showed good inhibition. Further, we have studied the binding of these compounds at the DFG-in and DFG-out conformations and reported the probable class (type I or type II). Hence, we suggest that these compounds could be novel drug leads for regulation of colorectal cancer.     

Fast and accurate protein substructure searching with simulated annealing and GPUs

Background  

Searching a database of protein structures for matches to a query structure, or occurrences of a structural motif, is an important task in structural biology and bioinformatics. While there are many existing methods for structural similarity searching, faster and more accurate approaches are still required, and few current methods are capable of substructure (motif) searching.     

Multiplexed expression and screening for recombinant protein production in mammalian cells

Background  

A variety of approaches to understanding protein structure and function require production of recombinant protein. Mammalian based expression systems have advantages over bacterial systems for certain classes of protein but can be slower and more laborious. Thus the availability of a simple system for production and rapid screening of constructs or conditions for mammalian expression would be of great benefit. To this end we have coupled an efficient recombinant protein production system based on transient transfection in HEK-293 EBNA1 (HEK-293E) suspension cells with a dot blot method allowing pre-screening of proteins expressed in cells in a high throughput manner.     

Partially-supervised protein subclass discovery with simultaneous annotation of functional residues

Background  

The study of functional subfamilies of protein domain families and the identification of the residues which determine substrate specificity is an important question in the analysis of protein domains. One way to address this question is the use of clustering methods for protein sequence data and approaches to predict functional residues based on such clusterings. The locations of putative functional residues in known protein structures provide insights into how different substrate specificities are reflected on the protein structure level.     

Secondary structure spatial conformation footprint: a novel method for fast protein structure comparison and classification

Background  

Recently a new class of methods for fast protein structure comparison has emerged. We call the methods in this class projection methods as they rely on a mapping of protein structure into a high-dimensional vector space. Once the mapping is done, the structure comparison is reduced to distance computation between corresponding vectors. As structural similarity is approximated by distance between projections, the success of any projection method depends on how well its mapping function is able to capture the salient features of protein structure. There is no agreement on what constitutes a good projection technique and the three currently known projection methods utilize very different approaches to the mapping construction, both in terms of what structural elements are included and how this information is integrated to produce a vector representation.     

Improved prediction of critical residues for protein function based on network and phylogenetic analyses

Background  

Phylogenetic approaches are commonly used to predict which amino acid residues are critical to the function of a given protein. However, such approaches display inherent limitations, such as the requirement for identification of multiple homologues of the protein under consideration. Therefore, complementary or alternative approaches for the prediction of critical residues would be desirable. Network analyses have been used in the modelling of many complex biological systems, but only very recently have they been used to predict critical residues from a protein's three-dimensional structure. Here we compare a couple of phylogenetic approaches to several different network-based methods for the prediction of critical residues, and show that a combination of one phylogenetic method and one network-based method is superior to other methods previously employed.     

Biomedical word sense disambiguation with ontologies and metadata: automation meets accuracy

Background  

Ontology term labels can be ambiguous and have multiple senses. While this is no problem for human annotators, it is a challenge to automated methods, which identify ontology terms in text. Classical approaches to word sense disambiguation use co-occurring words or terms. However, most treat ontologies as simple terminologies, without making use of the ontology structure or the semantic similarity between terms. Another useful source of information for disambiguation are metadata. Here, we systematically compare three approaches to word sense disambiguation, which use ontologies and metadata, respectively.     

A protein domain interaction interface database: InterPare

Background  

Most proteins function by interacting with other molecules. Their interaction interfaces are highly conserved throughout evolution to avoid undesirable interactions that lead to fatal disorders in cells. Rational drug discovery includes computational methods to identify the interaction sites of lead compounds to the target molecules. Identifying and classifying protein interaction interfaces on a large scale can help researchers discover drug targets more efficiently.     

A topological algorithm for identification of structural domains of proteins

Background  

Identification of the structural domains of proteins is important for our understanding of the organizational principles and mechanisms of protein folding, and for insights into protein function and evolution. Algorithmic methods of dissecting protein of known structure into domains developed so far are based on an examination of multiple geometrical, physical and topological features. Successful as many of these approaches are, they employ a lot of heuristics, and it is not clear whether they illuminate any deep underlying principles of protein domain organization. Other well-performing domain dissection methods rely on comparative sequence analysis. These methods are applicable to sequences with known and unknown structure alike, and their success highlights a fundamental principle of protein modularity, but this does not directly improve our understanding of protein spatial structure.     

Improved estimation of structure predictor quality

Background  

Methods that can automatically assess the quality of computationally predicted protein structures are important, as they enable the selection of the most accurate structure from an ensemble of predictions. Assessment methods that determine the quality of a predicted structure by comparing it against the various structures predicted by different servers have been shown to outperform approaches that rely on the intrinsic characteristics of the structure itself.