Randy J. Zauhar

Evidence for a strong sulfur–aromatic interaction derived from crystallographic data

We have uncovered new evidence for a significant interaction between divalent sulfur atoms and aromatic rings. Our study involves a statistical analysis of interatomic distances and other geometric descriptors derived from entries in the Cambridge Crystallographic Database (F. H. Allen and O. Kennard, Chem. Design Auto. News, 1993, Vol. 8, pp. 1 and 31-37). A set of descriptors was defined sufficient in number and type so as to elucidate completely the preferred geometry of interaction between six-membered aromatic carbon rings and divalent sulfurs for all crystal structures of nonmetal-bearing organic compounds present in the database. In order to test statistical significance, analogous probability distributions for the interaction of the moiety X-CH(2)-X with aromatic rings were computed, and taken a priori to correspond to the null hypothesis of no significant interaction. Tests of significance were carried our pairwise between probability distributions of sulfur-aromatic interaction descriptors and their CH(2)-aromatic analogues using the Smirnov-Kolmogorov nonparametric test (W. W. Daniel, Applied Nonparametric Statistics, Houghton-Mifflin: Boston, New York, 1978, pp. 276-286), and in all cases significance at the 99% confidence level or better was observed. Local maxima of the probability distributions were used to define a preferred geometry of interaction between the divalent sulfur moiety and the aromatic ring. Molecular mechanics studies were performed in an effort to better understand the physical basis of the interaction. This study confirms observations based on statistics of interaction of amino acids in protein crystal structures (R. S. Morgan, C. E. Tatsch, R. H. Gushard, J. M. McAdon, and P. K. Warme, International Journal of Peptide Protein Research, 1978, Vol. 11, pp. 209-217; R. S. Morgan and J. M. McAdon, International Journal of Peptide Protein Research, 1980, Vol. 15, pp. 177-180; K. S. C. Reid, P. F. Lindley, and J. M. Thornton, FEBS Letters, 1985, Vol. 190, pp. 209-213), as well as studies involving molecular mechanics (G. Nemethy and H. A. Scheraga, Biochemistry and Biophysics Research Communications, 1981, Vol. 98, pp. 482-487) and quantum chemical calculations (B. V. Cheney, M. W. Schulz, and J. Cheney, Biochimica Biophysica Acta, 1989, Vol. 996, pp.116-124; J. Pranata, Bioorganic Chemistry, 1997, Vol. 25, pp. 213-219)-all of which point to the possible importance of the sulfur-aromatic interaction. However, the preferred geometry of the interaction, as determined from our analysis of the small-molecule crystal data, differs significantly from that found by other approaches.

University bioinformatics programs on the rise

Fueled by strong demand from students and industrys need for trained bioinformaticists, universities are increasing their offerings in this fast-growing field.

Fragment-based Shape Signatures: a new tool for virtual screening and drug discovery

Since its introduction in 2003, the Shape Signatures method has been successfully applied in a number of drug design projects. Because it uses a ray-tracing approach to directly measure molecular shape and properties (as opposed to relying on chemical structure), it excels at scaffold hopping, and is extraordinarily easy to use. Despite its advantages, a significant drawback of the method has hampered its application to certain classes of problems; namely, when the chemical structures considered are large and contain heterogeneous ring-systems, the method produces descriptors that tend to merely measure the overall size of the molecule, and begin to lose selective power. To remedy this, the approach has been reformulated to automatically decompose compounds into fragments using ring systems as anchors, and to likewise partition the ray-trace in accordance with the fragment assignments. Subsequently, descriptors are generated that are fragment-based, and query and target molecules are compared by mapping query fragments onto target fragments in all ways consistent with the underlying chemical connectivity. This has proven to greatly extend the selective power of the method, while maintaining the ease of use and scaffold-hopping capabilities that characterized the original implementation. In this work, we provide a full conceptual description of the next generation Shape Signatures, and we underline the advantages of the method by discussing its practical applications to ligand-based virtual screening. The new approach can also be applied in receptor-based mode, where protein-binding sites (partitioned into subsites) can be matched against the new fragment-based Shape Signatures descriptors of library compounds.

Application of Screening Methods, Shape Signatures and Engineered Biosensors in Early Drug Discovery Process

PurposeIn this study, two unreported estrogen antagonists were identified using a combination of computational screening and a simple bacterial estrogen sensor.MethodsMolecules here presented were initially part of a group obtained from a library of over a half million chemical compounds, using the Shape Signatures method. The structures within this group were then clustered and compared to known antagonists based on their physico-chemical parameters, and possible binding modes of the compounds to the Estrogen Receptor α (ERα) were analyzed. Finally, thirteen candidate compounds were purchased, and two of them were shown to behave as potential subtype-selective estrogen antagonists using a set of bacterial estrogen biosensors, which included sensors for ERα, ERβ, and a negative control thyroid hormone β biosensor. These activities were then analyzed using an ELISA assay against activated ERα in human MCF-7 cell extract.ResultsTwo new estrogen receptor antagonists were detected using in silico Shape Signatures method with an engineered subtype-selective bacterial estrogen biosensor and commercially available ELISA assay. Additional thyroid biosensor control experiments confirmed no compounds interacted with human thyroid receptor β.ConclusionsThis work demonstrates an effective combination of computational analysis and simple bacterial screens for rapid identification of potential hormone-like therapeutics.

Computational analysis of EBNA1 “druggability” suggests novel insights for Epstein-Barr virus inhibitor design

The Epstein-Barr Nuclear Antigen 1 (EBNA1) is a critical protein encoded by the Epstein-Barr Virus (EBV). During latent infection, EBNA1 is essential for DNA replication and transcription initiation of viral and cellular genes and is necessary to immortalize primary B-lymphocytes. Nonetheless, the concept of EBNA1 as drug target is novel. Two EBNA1 crystal structures are publicly available and the first small-molecule EBNA1 inhibitors were recently discovered. However, no systematic studies have been reported on the structural details of EBNA1 “druggable” binding sites. We conducted computational identification and structural characterization of EBNA1 binding pockets, likely to accommodate ligand molecules (i.e. “druggable” binding sites). Then, we validated our predictions by docking against a set of compounds previously tested in vitro for EBNA1 inhibition (PubChem AID-2381). Finally, we supported assessments of pocket druggability by performing induced fit docking and molecular dynamics simulations paired with binding affinity predictions by Molecular Mechanics Generalized Born Surface Area calculations for a number of hits belonging to druggable binding sites. Our results establish EBNA1 as a target for drug discovery, and provide the computational evidence that active AID-2381 hits disrupt EBNA1:DNA binding upon interacting at individual sites. Lastly, structural properties of top scoring hits are proposed to support the rational design of the next generation of EBNA1 inhibitors.

Computer-aided identification of novel 3,5-substituted rhodanine derivatives with activity against Staphylococcus aureus DNA gyrase.

Methicillin resistant Staphylococcus aureus (MRSA) is among the major drug resistant bacteria that persist in both the community and clinical settings due to resistance to commonly used antimicrobials. This continues to fuel the need for novel compounds that are active against this organism. For this purpose we have targeted the type IIA bacterial topoisomerase, DNA gyrase, an essential enzyme involved in bacterial replication, through the ATP-dependent supercoiling of DNA. The virtual screening tool Shape Signatures was applied to screen a large database for agents with shape similar to Novobiocin, a known gyrase B inhibitor. The binding energetics of the top hits from this initial screen were further validated by molecular docking. Compounds with the highest score on available crystal structure of homologous DNA gyrase from Thermus thermophilus were selected. From this initial set of compounds, several rhodanine-substituted derivatives had the highest antimicrobial activity against S. aureus, as determined by minimal inhibitory concentration assays, with Novobiocin as the positive control. Further activity validation of the rhodanine compounds through biochemical assays confirmed their inhibition of both the supercoiling and the ATPase activity of DNA gyrase. Subsequent docking and molecular dynamics on the crystal structure of DNA gyrase from S. aureus when it became available, provides further rationalization of the observed biochemical activity and understanding of the receptor-ligand interactions. A regression model for MIC prediction against S. aureus is generated based on the current molecules studied as well as other rhodanines derivatives found in the literature.

Modeling androgen receptor flexibility: a binding mode hypothesis of CYP17 inhibitors/antiandrogens for prostate cancer therapy.

Prostate Cancer (PCa), a leading cause of cancer death worldwide (www.cancer.gov), is a complex malignancy where a spectrum of targets leads to a diversity of PCa forms. A widely pursued therapeutic target is the Androgen Receptor (AR). As a Steroid Hormone Receptor, AR serves as activator of transcription upon binding to androgens and plays a central role in the development of PCa. AR is a structurally flexible protein, and conformational plasticity of residues in the binding-pocket is a key to its ability to accommodate ligands from various chemical classes. Besides direct modulation of AR activity by antagonists, inhibition of cytochrome CYP17 (17α-hydroxylase/17,20-lyase), essential in androgen biosynthesis, has widely been considered an effective strategy against PCa. Interestingly, Handratta et al. (2005) discovered new, potent inhibitors of CYP17 (C-17 steroid derivatives) with pure AR antagonistic properties. Although the antiandrogenic activity of their lead compound (VN/124-1) has been experimentally proven both in vitro and in vivo, no structural data are currently available to elucidate the molecular determinants responsible for these desirable dual inhibitory properties. We implemented a Structure-based Drug Design (SBDD) approach to generate a valuable hypothesis as to the binding modes of steroidal CYP17 inhibitors/antiandrogens against the AR. To deal with the plasticity of residues buried in the Ligand Binding Domain (LBD), we developed a flexible-receptor Docking protocol based on Induced-Fit (IFD) methodology (www.schrodinger.com/). Our results constitute an ideal starting point for the rational design of next-generation analogues of CYP17 inhibitors/antiandrogens as well as an attractive tool to suggest novel chemical classes of AR antagonists.

An SH2 domain model of STAT5 in complex with phospho-peptides define "STAT5 Binding Signatures".

The signal transducer and activator of transcription 5 (STAT5) is a member of the STAT family of proteins, implicated in cell growth and differentiation. STAT activation is regulated by phosphorylation of protein monomers at conserved tyrosine residues, followed by binding to phospho-peptide pockets and subsequent dimerization. STAT5 is implicated in the development of severe pathological conditions, including many cancer forms. However, nowadays a few STAT5 inhibitors are known, and only one crystal structure of the inactive STAT5 dimer is publicly available. With a view to enabling structure-based drug design, we have: (1) analyzed phospho-peptide binding pockets on SH2 domains of STAT5, STAT1 and STAT3; (2) generated a model of STAT5 bound to phospho-peptides; (3) assessed our model by docking against a class of known STAT5 inhibitors (Müller et al. in ChemBioChem 9:723–727, 2008); (4) used molecular dynamics simulations to optimize the molecular determinants responsible for binding and (5) proposed unique “Binding Signatures” of STAT5. Our results put in place the foundations to address STAT5 as a target for rational drug design, from sequence, structural and functional perspectives.

Dual inhibition of Staphylococcus aureus DNA gyrase and topoisomerase IV activity by phenylalanine-derived (Z)-5-arylmethylidene rhodanines.

Methicillin resistant Staphylococcus aureus (MRSA) is a major drug resistant bacteria that persists in both community and clinical settings due to growing resistance to current drug regimens. Thus, there is a continued need for novel compounds that are active against this organism. Previously, we reported that various rhodanine derivatives inhibited the supercoiling activity of DNA gyrase. In this study, we determined the effect of new phenylalanine-derived (Z)-5-arylmethylidene rhodanines (which are efficacious against MRSA) on the activity of the two type II bacterial topoisomerases, DNA gyrase and topoisomerase IV (Topo IV). Compounds 1 and 9 showed the greatest efficacy against DNA gyrase with a minimal inhibitory concentration (MIC) of 5 μM while compounds 2 and 3 were the most efficacious against Topo IV with MIC values of 0.75 μM and 0.5 μM, respectively. Induced fit docking, using the crystallographic structures of the target enzymes, indicated that these rhodanine derivatives bind to the ATPase domain of gyrB and ParE subunits on DNA gyrase and Topo IV, respectively. These new compounds were efficacious against both DNA gyrase and Topo IV. The increased efficacy of these new rhodanine compounds, as compared to other rhodanine derivatives, results from their dual inhibition of DNA gyrase and Topo IV, thereby making them good candidates for further drug design and development.

RNA expression in human retina

Recent Genome-wide Association Studies (GWASs) for eye diseases/traits have delivered a number of novel findings across a diverse range of diseases, including age-related macular degeneration (AMD), glaucoma and refractive error. However, despite this astonishing rate of success, the major challenge still remains to not only confirm that the genes implicated in these studies are truly the genes conferring protection from or risk of disease but also to define the functional roles these genes play in disease. Ongoing evidence is accumulating that the single nucleotide polymorphisms (SNPs) used in GWAS and fine mapping studies have causal effects through their influence on gene expression rather than affecting protein function. The biological interpretation of SNP regulatory effects for a tissue requires knowledge of the transcriptome for that tissue. We summarize the reasons to characterize the complete retinal transcriptome as well as the evidence to include an assessment of differences in regional retinal expression.