GRID-Based Three-Dimensional Pharmacophores I: FLAPpharm, a Novel Approach for Pharmacophore Elucidation

Pharmacophore modeling using site-identification by ligand competitive saturation (SILCS) with multiple probe molecules

Receptor-based pharmacophore modeling is an efficient computer-aided drug design technique that uses the structure of the target protein to identify novel leads. However, most methods consider protein flexibility and desolvation effects in a very approximate way, which may limit their use in practice. The Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling protocol (SILCS-Pharm) was introduced recently to address these issues, as SILCS naturally takes both protein flexibility and desolvation effects into account by using full molecular dynamics simulations to determine 3D maps of the functional group-affinity patterns on a target receptor. In the present work, the SILCS-Pharm protocol is extended to use a wider range of probe molecules including benzene, propane, methanol, formamide, acetaldehyde, methylammonium, acetate and water. This approach removes the previous ambiguity brought by using water as both the hydrogen-bond donor and acceptor probe molecule. The new SILCS-Pharm protocol is shown to yield improved screening results, as compared to the previous approach based on three target proteins. Further validation of the new protocol using five additional protein targets showed improved screening compared to those using common docking methods, further indicating improvements brought by the explicit inclusion of additional feature types associated with the wider collection of probe molecules in the SILCS simulations. The advantage of using complementary features and volume constraints, based on exclusion maps of the protein defined from the SILCS simulations, is presented. In addition, reranking using SILCS-based ligand grid free energies is shown to enhance the diversity of identified ligands for the majority of targets. These results suggest that the SILCS-Pharm protocol will be of utility in rational drug design.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.

BioGPS descriptors for rational engineering of enzyme promiscuity and structure based bioinformatic analysis

A new bioinformatic methodology was developed founded on the Unsupervised Pattern Cognition Analysis of GRID-based BioGPS descriptors (Global Positioning System in Biological Space). The procedure relies entirely on three-dimensional structure analysis of enzymes and does not stem from sequence or structure alignment. The BioGPS descriptors account for chemical, geometrical and physical-chemical features of enzymes and are able to describe comprehensively the active site of enzymes in terms of “pre-organized environment” able to stabilize the transition state of a given reaction. The efficiency of this new bioinformatic strategy was demonstrated by the consistent clustering of four different Ser hydrolases classes, which are characterized by the same active site organization but able to catalyze different reactions. The method was validated by considering, as a case study, the engineering of amidase activity into the scaffold of a lipase. The BioGPS tool predicted correctly the properties of lipase variants, as demonstrated by the projection of mutants inside the BioGPS “roadmap”.

Playing with opening and closing of heterocycles: using the cusmano-ruccia reaction to develop a novel class of oxadiazolothiazinones, active as calcium channel modulators and P-glycoprotein inhibitors

As a result of the ring-into-ring conversion of nitrosoimidazole derivatives, we obtained a molecular scaffold that, when properly decorated, is able to decrease inotropy by blocking L-type calcium channels. Previously, we used this scaffold to develop a quantitative structure-activity relationship (QSAR) model, and we used the most potent oxadiazolothiazinone as a template for ligand-based virtual screening. Here, we enlarge the diversity of chemical decorations, present the synthesis and in vitro data for 11 new derivatives, and develop a new 3D-QSAR model with recent in silico techniques. We observed a key role played by the oxadiazolone moiety: given the presence of positively charged calcium ions in the transmembrane channel protein, we hypothesize the formation of a ternary complex between the oxadiazolothiazinone, the Ca2+ ion and the protein. We have supported this hypothesis by means of pharmacophore generation and through the docking of the pharmacophore into a homology model of the protein. We also studied with docking experiments the interaction with a homology model of P-glycoprotein, which is inhibited by this series of molecules, and provided further evidence toward the relevance of this scaffold in biological interactions.

Setting the record straight: the origin of the pharmacophore concept

For over a century since the early 1900s, Paul Ehrlich was credited with originating the concept of pharmacophores. This was challenged by John Van Drie in 2007 due to the fact that Ehrlich did not use the word "pharmacophore" in his writings. Van Drie claimed that the attribution of the pharmacophore concept to Ehrlich was due to an erroneous citation made by Ariëns in a 1966 paper, and instead he claimed, Lemont B. Kier developed the pharmacophore concept (in the modern sense, as defined by the IUPAC) during 1967-1971. There are two separate issues that may have triggered this conflict. The first one is the shift in the meaning of pharmacophore from "chemical groups" to patterns of "abstract features" of a molecule that are responsible for a biological effect. Indeed, the original use of the term is different than the current definition proposed by the IUPAC. The term was redefined in 1960 by Schueler, and this modification formed the basis of IUPAC's modern definition. The second issue is the origin of the "concept" of pharmacophore. While Ehrlich's contemporaries have consistently attributed the origin of the concept to him, the issue is further complicated by the fact that Ehrlich did not use the term pharmacophore in his papers. He, instead, referred to the features of a molecule that are responsible for biological effects as toxophores, while his contemporaries were using the term pharmacophore for the same features. In this paper, we resolve any doubts about the origins of the pharmacophore concept. Our research points to Paul Ehrlich's 1898 paper for originating the concept, which identifies peripheral chemical groups in molecules responsible for binding that leads to the subsequent biological effect, and to Schueler's 1960 book that extends the concept to the modern definition where spatial patterns of abstract features of a molecule define the pharmacophore and are ultimately responsible for the biological effect.

Site-Identification by Ligand Competitive Saturation (SILCS) assisted pharmacophore modeling

Database screening using receptor-based pharmacophores is a computer-aided drug design technique that uses the structure of the target molecule (i.e. protein) to identify novel ligands that may bind to the target. Typically receptor-based pharmacophore modeling methods only consider a single or limited number of receptor conformations and map out the favorable binding patterns in vacuum or with a limited representation of the aqueous solvent environment, such that they may suffer from neglect of protein flexibility and desolvation effects. Site-Identification by Ligand Competitive Saturation (SILCS) is an approach that takes into account these, as well as other, properties to determine 3-dimensional maps of the functional group-binding patterns on a target receptor (i.e. FragMaps). In this study, a method to use the FragMaps to automatically generate receptor-based pharmacophore models is presented. It converts the FragMaps into SILCS pharmacophore features including aromatic, aliphatic, hydrogen-bond donor and acceptor chemical functionalities. The method generates multiple pharmacophore hypotheses that are then quantitatively ranked using SILCS grid free energies. The pharmacophore model generation protocol is validated using three different protein targets, including using the resulting models in virtual screening. Improved performance and efficiency of the SILCS derived pharmacophore models as compared to published docking studies, as well as a recently developed receptor-based pharmacophore modeling method is shown, indicating the potential utility of the approach in rational drug design.

Isozyme-specific ligands for O-acetylserine sulfhydrylase, a novel antibiotic target

The last step of cysteine biosynthesis in bacteria and plants is catalyzed by O-acetylserine sulfhydrylase. In bacteria, two isozymes, O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B, have been identified that share similar binding sites, although the respective specific functions are still debated. O-acetylserine sulfhydrylase plays a key role in the adaptation of bacteria to the host environment, in the defense mechanisms to oxidative stress and in antibiotic resistance. Because mammals synthesize cysteine from methionine and lack O-acetylserine sulfhydrylase, the enzyme is a potential target for antimicrobials. With this aim, we first identified potential inhibitors of the two isozymes via a ligand- and structure-based in silico screening of a subset of the ZINC library using FLAP. The binding affinities of the most promising candidates were measured in vitro on purified O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B from Salmonella typhimurium by a direct method that exploits the change in the cofactor fluorescence. Two molecules were identified with dissociation constants of 3.7 and 33 µM for O-acetylserine sulfhydrylase-A and O-acetylserine sulfhydrylase-B, respectively. Because GRID analysis of the two isoenzymes indicates the presence of a few common pharmacophoric features, cross binding titrations were carried out. It was found that the best binder for O-acetylserine sulfhydrylase-B exhibits a dissociation constant of 29 µM for O-acetylserine sulfhydrylase-A, thus displaying a limited selectivity, whereas the best binder for O-acetylserine sulfhydrylase-A exhibits a dissociation constant of 50 µM for O-acetylserine sulfhydrylase-B and is thus 8-fold selective towards the former isozyme. Therefore, isoform-specific and isoform-independent ligands allow to either selectively target the isozyme that predominantly supports bacteria during infection and long-term survival or to completely block bacterial cysteine biosynthesis.