Search

Becoming mathematically literate means not only knowledge of the content of mathematics but also understanding of the nature of mathematics and the literate practices involved in the creation, communication, and consumption of its content. This case study examined one mathematician’s view of the nature of mathematics and his literate practices. Data collected include semi-structured interviews with the mathematician, observations of his daily work routines, and his think-alouds during the reading of a disciplinary text. These data were analyzed qualitatively through an iterative process involving multiple readings and identification and refinement of codes. The analysis revealed that the mathematician (a) viewed mathematics as rigorous and demanding, both theoretical and practical, relatively stable but highly rewarding, interconnected with other disciplines, and involving discipline-legitimated discursive practices; (b) engaged in extensive reading/viewing and writing, valued learning from repeated trials and errors, and collaborated with an international network of peers in research; (c) used a range of reading strategies (e.g. close reading, summarizing, questioning, storying, evaluating, annotating) to help him make sense of the disciplinary text; and (d) marshaled both verbal and visual resources to create specialized knowledge, engage in rigorous reasoning, develop logical argument, and construct professional identity. These findings provide important insights that can help teachers design activities that are authentic to mathematics practices and effective for promoting mathematics literacy.

Colchicine site tubulin inhibitors are currently developed as vascular disrupting agents (VDAs). However, they were found to have cardiotoxicity in clinical trials. To overcome the problem, we developed a stilbene derivative, cis-3, 4', 5-trimethoxy-3'-aminostilbene (stilbene 5c), which is highly potent and has no bone marrow and cardiac toxicity in mice. Here we attempt to optimize stilbene 5c using computer-based drug design and synthesize derivatives with benzimidazole or indole group. Biological evaluation showed that they are weaker than stilbene 5c without better water solubility. Alternative approach was thus adopted to make prodrugs of stilbene 5c. A water-soluble prodrug PD7 was synthesized by addition of a morpholino group with carbamate linkage to the amino group of stilbene 5c. In vitro studies show that PD7 induces mitotic arrest and disrupts microtubule similar to stilbene 5c. The cell signaling events in Cdc2, p53, Akt, and aurora kinase are similar in cells treated with stilbene 5c, CA4 or PD7, suggesting that they share the same mechanism. Although PD7 is less effective than stilbene 5c in vitro, the biological activity of PD7 as a single agent is similar to that of stilbene 5c. Combination of PD7 with VEGF inhibitor bevacizumab significantly enhances the therapeutic efficacy of PD7 in mouse xenograft model. These data suggest that PD7 could be a good candidate for further pre-clinical and clinical development as a new VDA for cancer therapy. [ABSTRACT FROM AUTHOR]

A method is described to deconstruct the network of hydropathic interactions within and between a protein's sidechain and its environment into residue-based three-dimensional maps. These maps encode favorable and unfavorable hydrophobic and polar interactions, in terms of spatial positions for optimal interactions, relative interaction strength, as well as character. In addition, these maps are backbone angle-dependent. After map calculation and clustering, a finite number of unique residue sidechain interaction maps exist for each backbone conformation, with the number related to the residue's size and interaction complexity. Structures for soluble proteins (~749,000 residues) and membrane proteins (~387,000 residues) were analyzed, with the latter group being subdivided into three subsets related to the residue's position in the membrane protein: soluble domain, core-facing transmembrane domain, and lipid-facing transmembrane domain. This work suggests that maps representing residue types and their backbone conformation can be reassembled to optimize the medium-to-high resolution details of a protein structure. In particular, the information encoded in maps constructed from the lipid-facing transmembrane residues appears to paint a clear picture of the protein-lipid interactions that are difficult to obtain experimentally.

Sphingolipids are a diverse class of lipids with essential functions as determinants of membrane physical properties and as intra- and intercellular signaling agents. Disruption of the normal biochemical processes that establish the levels of individual sphingolipids is associated with a variety of human diseases including cancer, cardiovascular disease, metabolic disease, skin diseases, and lysosomal storage diseases. A unique aspect of this metabolic network is that there is a single enzymatic step that initiates the biosynthetic pathway for all sphingolipids. This step is catalyzed by the enzyme serine palmitoyltranserase (SPT). Under most circumstances SPT condenses serine and the 16-carbon acyl-CoA, palmitoyl-CoA to produce the precursor of all sphingolipids. SPT, a four-subunit protein complex, is subject to classic feedback regulation: when cellular sphingolipids are elevated, SPT activity is inhibited. Ceramide is the sphingolipid sensed by this system and it regulates SPT by directly binding to the complex. The ceramide binding site in the SPT complex, and how ceramide binding results in SPT inhibition, has now been determined in vertebrates, plants, and yeast using molecular modeling and cryo-electron microscopy. Here we discuss the similarities and differences revealed by these resolved structures and the surprising result that ceramide binds at almost identical positions in the SPT complex of these divergent organisms, but accomplishes SPT regulation in very different ways., Competing Interests: Declaration of competing interest The authors have no conflict of interests to declare., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Human sirtuins play important roles in various cellular events including DNA repair, gene silencing, mitochondrial biogenesis, insulin secretion and apoptosis. They regulate a wide array of protein and enzyme targets through their NAD + -dependent deacetylase activities. Sirtuins are also thought to mediate the beneficial effects of low-calorie intake to extend longevity in diverse organisms from yeast to mammals. Small molecules mimicking calorie restriction to stimulate sirtuin activity are attractive therapeutics against age-related disorders such as cardiovascular diseases, diabetes and neurodegeneration. Little is known about one of the mitochondrial sirtuins, SIRT5. SIRT5 has emerged as a critical player in maintaining cardiac health and neuronal viability upon stress and functions as a tumour suppressor in a context-specific manner. Much has been debated about whether SIRT5 has evolved away from being a deacetylase because of its weak catalytic activity, especially in the in vitro testing. We have, for the first time, identified a SIRT5-selective allosteric activator, nicotinamide riboside (NR). It can increase SIRT5 catalytic efficiency with different synthetic peptide substrates. The mechanism of action was further explored using a combination of molecular biology and biochemical strategies. Based on the existing structural biology information, the NR binding site was also mapped out. These activators are powerful chemical probes for the elucidation of cellular regulations and biological functions of SIRT5. The knowledge gained in this study can be used to guide the design and synthesis of more potent, isotype-selective SIRT5 activators and to develop them into therapeutics for metabolic disorders and age-related diseases., (© 2023 The Authors. The FEBS Journal published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

The Department of Medicinal Chemistry, together with the Institute for Structural Biology, Drug Discovery and Development, at Virginia Commonwealth University (VCU) has evolved, organically with quite a bit of bootstrapping, into a unique drug discovery ecosystem in response to the environment and culture of the university and the wider research enterprise. Each faculty member that joined the department and/or institute added a layer of expertise, technology and most importantly, innovation, that fertilized numerous collaborations within the University and with outside partners. Despite moderate institutional support with respect to a typical drug discovery enterprise, the VCU drug discovery ecosystem has built and maintained an impressive array of facilities and instrumentation for drug synthesis, drug characterization, biomolecular structural analysis and biophysical analysis, and pharmacological studies. Altogether, this ecosystem has had major impacts on numerous therapeutic areas, such as neurology, psychiatry, drugs of abuse, cancer, sickle cell disease, coagulopathy, inflammation, aging disorders and others. Novel tools and strategies for drug discovery, design and development have been developed at VCU in the last five decades; e.g., fundamental rational structure-activity relationship (SAR)-based drug design, structure-based drug design, orthosteric and allosteric drug design, design of multi-functional agents towards polypharmacy outcomes, principles on designing glycosaminoglycans as drugs, and computational tools and algorithms for quantitative SAR (QSAR) and understanding the roles of water and the hydrophobic effect., Competing Interests: Declaration of Competing Interests The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Inc. All rights reserved.)

A long-lasting goal of computational biochemists, medicinal chemists, and structural biologists has been the development of tools capable of deciphering the molecule-molecule interaction code that produces a rich variety of complex biomolecular assemblies comprised of the many different simple and biological molecules of life: water, small metabolites, cofactors, substrates, proteins, DNAs, and RNAs. Software applications that can mimic the interactions amongst all of these species, taking account of the laws of thermodynamics, would help gain information for understanding qualitatively and quantitatively key determinants contributing to the energetics of the bimolecular recognition process. This, in turn, would allow the design of novel compounds that might bind at the intermolecular interface by either preventing or reinforcing the recognition. HINT, hydropathic interaction, was a model and software code developed from a deceptively simple idea of Donald Abraham with the close collaboration with Glen Kellogg at Virginia Commonwealth University. HINT is based on a function that scores atom-atom interaction using LogP, the partition coefficient of any molecule between two phases; here, the solvents are water that mimics the cytoplasm milieu and octanol that mimics the protein internal hydropathic environment. This review summarizes the results of the extensive and successful collaboration between Abraham and Kellogg at VCU and the group at the University of Parma for testing HINT in a variety of different biomolecular interactions, from proteins with ligands to proteins with DNA., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Kellogg, Marabotti, Spyrakis and Mozzarelli.)

The aliphatic hydrophobic amino acid residues-alanine, isoleucine, leucine, proline and valine-are among the most common found in proteins. Their structural role in proteins is seemingly obvious: engage in hydrophobic interactions to stabilize secondary, and to a lesser extent, tertiary and quaternary structure. However, favorable hydrophobic interactions involving the sidechains of these residue types are generally less significant than the unfavorable set arising from interactions with polar atoms. Importantly, the constellation of interactions between residue sidechains and their environments can be recorded as three-dimensional maps that, in turn, can be clustered. The clustered average map sets compose a library of interaction profiles encoding interaction strengths, interaction types and the optimal 3D position for the interacting partners. This library is backbone angle-dependent and suggests solvent and lipid accessibility for each unique interaction profile. In this work, in addition to analysis of soluble proteins, a large set of membrane proteins that contained optimized artificial lipids were evaluated by parsing the structures into three distinct components: soluble extramembrane domain, lipid facing transmembrane domain, core transmembrane domain. The aliphatic residues were extracted from each of these sets and passed through our calculation protocol. Notable observations include: the roles of aliphatic residues in soluble proteins and in the membrane protein's soluble domains are nearly identical, although the latter are slightly more solvent accessible; by comparing maps calculated with sidechain-lipid interactions to maps ignoring those interactions, the potential extent of residue-lipid and residue-interactions can be assessed and likely exploited in structure prediction and modeling; amongst these residue types, the levels of lipid engagement show isoleucine as the most engaged, while the other residues are largely interacting with neighboring helical residues., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 AL Mughram, Catalano, Herrington, Safo and Kellogg.)

Pyridoxal 5'-phosphate (PLP), the active form of vitamin B6, serves as a cofactor for scores of B6-dependent (PLP-dependent) enzymes involved in many cellular processes. One such B6 enzyme is dopa decarboxylase (DDC), which is required for the biosynthesis of key neurotransmitters, e.g., dopamine and serotonin. PLP-dependent enzymes are biosynthesized as apo-B6 enzymes and then converted to the catalytically active holo-B6 enzymes by Schiff base formation between the aldehyde of PLP and an active site lysine of the protein. In eukaryotes, PLP is made available to the B6 enzymes through the activity of the B6-salvage enzymes, pyridoxine 5'-phosphate oxidase (PNPO) and pyridoxal kinase (PLK). To minimize toxicity, the cell keeps the content of free PLP (unbound) very low through dephosphorylation and PLP feedback inhibition of PNPO and PLK. This has led to a proposed mechanism of complex formation between the B6-salvage enzymes and apo-B6 enzymes prior to the transfer of PLP, although such complexes are yet to be characterized at the atomic level, presumably due to their transient nature. A computational study, for the first time, was used to predict a likely PNPO and DDC complex, which suggested contact between the allosteric PLP tight-binding site on PNPO and the active site of DDC. Using isothermal calorimetry and/or surface plasmon resonance, we also show that PNPO binds both apoDDC and holoDDC with dissociation constants of 0.93 ± 0.07 μM and 2.59 ± 0.11 μM, respectively. Finally, in the presence of apoDDC, the tightly bound PLP on PNPO is transferred to apoDDC, resulting in the formation of about 35% holoDDC.

Ataxia telangiectasia mutated (ATM) is a serine-threonine protein kinase and important regulator of the DNA damage response (DDR). One critical ATM target is the structural subunit A (PR65-S401) of protein phosphatase 2A (PP2A), known to regulate diverse cellular processes such as mitosis and cell growth as well as dephosphorylating many proteins during the recovery from the DDR. We generated mouse embryonic fibroblasts expressing PR65-WT, -S401A (cannot be phosphorylated), and -S401D (phospho-mimetic) transgenes. Significantly, S401 mutants exhibited extensive chromosomal aberrations, impaired DNA double-strand break (DSB) repair and underwent increased mitotic catastrophe after radiation. Both S401A and the S401D cells showed impaired DSB repair (nonhomologous end joining and homologous recombination repair) and exhibited delayed DNA damage recovery, which was reflected in reduced radiation survival. Furthermore, S401D cells displayed increased ERK and AKT signaling resulting in enhanced growth rate further underscoring the multiple roles ATM-PP2A signaling plays in regulating prosurvival responses. Time-lapse video and cellular localization experiments showed that PR65 was exported to the cytoplasm after radiation by CRM1, a nuclear export protein, in line with the very rapid pleiotropic effects observed. A putative nuclear export sequence (NES) close to S401 was identified and when mutated resulted in aberrant PR65 shuttling. Our study demonstrates that the phosphorylation of a single, critical PR65 amino acid (S401) by ATM fundamentally controls the DDR, and balances DSB repair quality, cell survival and growth by spatiotemporal PR65 nuclear-cytoplasmic shuttling mediated by the nuclear export receptor CRM1., (© 2022. The Author(s).)

Evaluation of the intramolecular stability of proteins plays a key role in the comprehension of their biological behavior and mechanism of action. Small structural alterations such as mutations induced by single nucleotide polymorphism can impact biological activity and pharmacological modulation. Covid-19 mutations, that affect viral replication and the susceptibility to antibody neutralization, and the action of antiviral drugs, are just one example. In this work, the intramolecular stability of mutated proteins, like Spike glycoprotein and its complexes with the human target, is evaluated through hydropathic intramolecular energy scoring originally conceived by Abraham and Kellogg based on the "Extension of the fragment method to calculate amino acid zwitterion and side-chain partition coefficients" by Abraham and Leo in Proteins: Struct. Funct. Genet. 1987, 2:130 - 52. HINT is proposed as a fast and reliable tool for the stability evaluation of any mutated system. This work has been written in honor of Prof. Donald J. Abraham (1936-2021)., (© 2022. The Author(s).)

Background: Deregulated translation initiation is implicated extensively in cancer initiation and progression. It is actively pursued as a viable target that circumvents the dependency on oncogenic signaling, a significant factor in current strategies. Eukaryotic translation initiation factor (eIF) 4A plays an essential role in translation initiation by unwinding the secondary structure of messenger RNA (mRNA) upstream of the start codon, enabling active ribosomal recruitment on the downstream genes. Several natural product molecules with similar scaffolds, such as Rocaglamide A (RocA), targeting eIF4A have been reported in the last decade. However, their clinical utilization is still elusive due to several pharmacological limitations. In this study we identified new eIF4A1 inhibitors and their possible mechanisms., Methods: In this report, we conducted a pharmacophore-based virtual screen of RocA complexed with eIF4A and a polypurine RNA strand for novel eIF4A inhibitors from commercially available compounds in the MolPort Database. We performed target-based screening and optimization of active pharmacophores. We assessed the effects of novel compounds on biochemical and cell-based assays for efficacy and mechanistic evaluation., Results: We validated three new potent eIF4A inhibitors, RBF197, RBF 203, and RBF 208, which decreased diffuse large B-cell lymphoma (DLBCL) cell viability. Biochemical and cellular studies, molecular docking, and functional assays revealed that thosenovel compounds clamp eIF4A into mRNA in an ATP-independent manner. Moreover, we found that RBF197 and RBF208 significantly depressed eIF4A-dependent oncogene expression as well as the colony formation capacity of DLBCL. Interestingly, exposure of these compounds to non-malignant cells had only minimal impact on their growth and viability., Conclusions: Identified compounds suggest a new strategy for designing novel eIF4A inhibitors., (© 2022. The Author(s).)

Knowledge of three-dimensional protein structure is integral to most modern drug discovery efforts. Recent advancements have highlighted new techniques for 3D protein structure determination and, where structural data cannot be collected experimentally, prediction of protein structure. We have undertaken a major effort to use existing protein structures to collect, characterize, and catalogue the inter-atomic interactions that define and compose 3D structure by mapping hydropathic interaction environments as maps in 3D space. This work has been performed on a residue-by-residue basis, where we have seen evidence for relationships between environment character, residue solvent-accessible surface areas and their secondary structures. In this graphical review, we apply principles from our earlier studies and expand the scope to all common amino acid residue types in both soluble and membrane proteins. Key to this analysis is parsing the Ramachandran plot to an 8-by-8 chessboard to define secondary structure bins. Our analysis yielded a number of quantitative discoveries: 1) increased fraction of hydrophobic residues (alanine, isoleucine, leucine, phenylalanine and valine) in membrane proteins compared to their fractions in soluble proteins; 2) less burial coupled with significant increases in favorable hydrophobic interactions for hydrophobic residues in membrane proteins compared to soluble proteins; and 3) higher burial and more favorable polar interactions for polar residues now preferring the interior of membrane proteins. These observations and the supporting data should provide benchmarks for current studies of protein residues in different environments and may be able to guide future protein structure prediction efforts., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2021 The Author(s).)

Aspartic acid, glutamic acid and histidine are ionizable residues occupying various protein environments and perform many different functions in structures. Their roles are tied to their acid/base equilibria, solvent exposure, and backbone conformations. We propose that the number of unique environments for ASP, GLU and HIS is quite limited. We generated maps of these residue's environments using a hydropathic scoring function to record the type and magnitude of interactions for each residue in a 2703-protein structural dataset. These maps are backbone-dependent and suggest the existence of new structural motifs for each residue type. Additionally, we developed an algorithm for tuning these maps to any pH, a potentially useful element for protein design and structure building. Here, we elucidate the complex interplay between secondary structure, relative solvent accessibility, and residue ionization states: the degree of protonation for ionizable residues increases with solvent accessibility, which in turn is notably dependent on backbone structure., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Herrington and Kellogg.)

Atomic-resolution protein structural models are prerequisites for many downstream activities like structure-function studies or structure-based drug discovery. Unfortunately, this data is often unavailable for some of the most interesting and therapeutically important proteins. Thus, computational tools for building native-like structural models from less-than-ideal experimental data are needed. To this end, interaction homology exploits the character, strength and loci of the sets of interactions that define a structure. Each residue type has its own limited set of backbone angle-dependent interaction motifs, as defined by their environments. In this work, we characterize the interactions of serine, cysteine and S-bridged cysteine in terms of 3D hydropathic environment maps. As a result, we explore several intriguing questions. Are the environments different between the isosteric serine and cysteine residues? Do some environments promote the formation of cystine S-S bonds? With the increasing availability of structural data for water-insoluble membrane proteins, are there environmental differences for these residues between soluble and membrane proteins? The environments surrounding serine and cysteine residues are dramatically different: serine residues are about 50% solvent exposed, while cysteines are only 10% exposed; the latter are more involved in hydrophobic interactions although there are backbone angle-dependent differences. Our analysis suggests that one driving force for -S-S- bond formation is a rather substantial increase in burial and hydrophobic interactions in cystines. Serine and cysteine become less and more, respectively, solvent-exposed in membrane proteins. 3D hydropathic environment maps are an evolving structure analysis tool showing promise as elements in a new protein structure prediction paradigm., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: 1910 Genetics, Inc. provided support for preliminary studies but had no role in the reporting of data or derived conclusions., (© 2021 The Author(s).)

Previously, we discovered that FOSL1 facilitates the metastasis of head and neck squamous cell carcinoma (HNSCC) cancer stem cells in a spontaneous mouse model. However, the molecular mechanisms remained unclear. Here, we demonstrated that FOSL1 serves as the dominant activating protein 1 (AP1) family member and is significantly upregulated in HNSCC tumor tissues and correlated with metastasis of HNSCC. Mechanistically, FOSL1 exerts its function in promoting tumorigenicity and metastasis predominantly via selective association with Mediators to establish super-enhancers (SEs) at a cohort of cancer stemness and pro-metastatic genes, such as SNAI2 and FOSL1 itself. Depletion of FOSL1 led to disruption of SEs and expression inhibition of these key oncogenes, which resulted in the suppression of tumor initiation and metastasis. We also revealed that the abundance of FOSL1 is positively associated with the abundance of SNAI2 in HNSCC and the high expression levels of FOSL1 and SNAI2 are associated with short overall disease-free survival. Finally, the administration of the FOSL1 inhibitor SR11302 significantly suppressed tumor growth and lymph node metastasis of HNSCC in a patient-derived xenograft model. These findings indicate that FOSL1 is a master regulator that promotes the metastasis of HNSCC through a SE-driven transcription program that may represent an attractive target for therapeutic interventions., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The American Society of Gene and Cell Therapy. Published by Elsevier Inc. All rights reserved.)

In our continuing efforts to develop novel neuroprotectants for Alzheimer's disease (AD), a series of analogs based on a lead compound that was recently shown to target the mitochondrial complex I were designed, synthesized and biologically characterized to understand the structure features that are important for neuroprotective activities. The results from a cellular AD model highlighted the important roles of the 4-OH on the phenyl ring and the 5-OCH 3 on the indole ring of the lead compound. The results also demonstrated that the β-keto moiety can be modified to retain or improve the neuroprotective activity. Docking studies of selected analogs to the FMN site of mitochondrial complex I also supported the observed neuroprotective activities. Collectively, the results provide further information to guide optimization and development of analogs based on this chemical scaffold as neuroprotectants with a novel mechanism of action for AD., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Three-dimensional (3D) maps of the hydropathic environments of protein amino acid residues are information-rich descriptors of preferred conformations, interaction types and energetics, and solvent accessibility. The interactions made by each residue are the primary factor for rotamer selection and the secondary, tertiary, and even quaternary protein structure. Our evolving basis set of environmental data for each residue type can be used to understand the protein structure. This work focuses on the aromatic residues phenylalanine, tyrosine, and tryptophan and their structural roles. We calculated and analyzed side chain-to-environment 3D maps for over 70,000 residues of these three types that reveal, with respect to hydrophobic and polar interactions, the environment around each. After binning with backbone ϕ/ψ and side chain χ 1 , we clustered each bin by 3D similarities between map-map pairs. For each of the three residue types, four bins were examined in detail: one in the β-pleat, two in the right-hand α-helix, and one in the left-hand α-helix regions of the Ramachandran plot. For high degrees of side chain burial, encapsulation of the side chain by hydrophobic interactions is ubiquitous. The more solvent-exposed side chains are more likely to be involved in polar interactions with their environments. Evidence for π-π interactions was observed in about half of the residues surveyed [phenylalanine (PHE): 53.3%, tyrosine (TYR): 34.1%, and tryptophan (TRP): 55.7%], but on an energy basis, this contributed to only ∼4% of the total. Evidence for π-cation interactions was observed in 14.1% of PHE, 8.3% of TYR, and 26.8% of TRP residues, but on an energy basis, this contributed to only ∼1%. This recognition of even these subtle interactions in the 3D hydropathic environment maps is key support for our interaction homology paradigm of protein structure elucidation and possibly prediction.

Mitochondrial dysfunction has been recognized as an essential contributor to many human diseases including neurodegenerative disorders. However, the exact pathological role of mitochondrial dysfunction, especially in mitochondrial reactive oxygen species-associated oxidative stress, remains elusive, partially due to the lack of chemical probes with well-defined mechanisms of action. Herein, we describe the characterization and discovery of a rationally designed small molecule ZCM-I-1 as a selective modulator of the production of reactive oxygen species from mitochondrial complex I that does not alter mitochondrial membrane potential and bioenergetics. Chemical biology studies employing photoaffinity probes derived from ZCM-I-1 demonstrated its novel mechanism of action of modulating complex I via interactions with the flavin mononucleotide site, proximal in the reaction pathway within complex I.

The serotonin 5-HT 7 G protein-coupled receptor (GPCR) is a proposed pharmacotherapeutic target for a variety of central and peripheral indications, albeit, there are no approved drugs selective for binding 5-HT 7 . We previously reported that a lead analog based on the 5-substituted-N,N-disubstituted-1,2,3,4-tetrahydronaphthalen-2-amine (5-substituted-2-aminotetralin, 5-SAT) scaffold binds with high affinity at the 5-HT 7 GPCR, and can treat symptoms of autism in mouse models; subsequently, the lead was found to have high affinity at the 5-HT 1A GPCR. Herein, we report the synthesis of novel 5-SAT analogs to develop a 3-dimensional quantitative structure-affinity relationship (3D-QSAR) at the human 5-HT 7 receptor for comparison with similar studies at the highly homologous 5-HT 1A receptor. We report 35 new 5-SAT ligands, some with very high affinity (K i  ≤ 1 nM) and stereoselectivity at 5-HT 7 + or 5-HT 1A receptors, several with modest selectivity (up to 12-fold) for binding at 5-HT 7 , and, several ligands with high selectivity (up to 40-fold) at the 5-HT 1A receptor. 3D-QSAR results indicate that steric extensions at the C(5)-position improve selectivity for the 5-HT 7 over 5-HT 1A receptor, while steric and hydrophobic extensions at the chiral C(2)-amino position impart 5-HT 1A selectivity. In silico receptor homology modeling studies, supplemented with molecular dynamics simulations and binding free energy calculations, were used to rationalize experimentally-determined receptor selectivity and stereoselective affinity results. The data from these studies indicate that the 5-SAT chemotype, previously shown to be safe and efficacious in rodent paradigms of neurodevelopmental and neuropsychiatric disorders, is amenable to structural modification to optimize affinity at serotonin 5-HT 7 vs. 5-HT 1A GPCRs, as may be required for successful clinical translation., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Phosphatidylinositol-transfer proteins (PITPs) are key regulators of lipid signaling in eukaryotic cells. These proteins both potentiate the activities of phosphatidylinositol (PtdIns) 4-OH kinases and help channel production of specific pools of phosphatidylinositol 4-phosphate (PtdIns(4)P) dedicated to specific biological outcomes. In this manner, PITPs represent a major contributor to the mechanisms by which the biological outcomes of phosphoinositide are diversified. The two-ligand priming model proposes that the engine by which Sec14-like PITPs potentiate PtdIns kinase activities is a heterotypic lipid-exchange cycle where PtdIns is a common exchange substrate among the Sec14-like PITP family, but the second exchange ligand varies with the PITP. A major prediction of this model is that second-exchangeable ligand identity will vary from PITP to PITP. To address the heterogeneity in the second exchange ligand for Sec14-like PITPs, we used structural, computational, and biochemical approaches to probe the diversities of the lipid-binding cavity microenvironments of the yeast Sec14-like PITPs. The collective data report that yeast Sec14-like PITP lipid-binding pockets indeed define diverse chemical microenvironments that translate into differential ligand-binding specificities across this protein family., (© 2019 Tripathi et al.)

Analyses of the hydropathic environments of protein amino acid residues reveal structural information on multiple levels. The interactions made by each residue are the basis for sidechain (rotamer) conformation and ultimately for secondary, tertiary and even quaternary protein structure. By identifying and characterizing the interactions for each residue type, we are developing a basis set of environmental data that can be used to understand protein structure. This work focuses alanine and its roles. We calculated and analyzed separately backbone-to-environment and sidechain-to-environment 3D maps for over 57,000 alanines that, with respect to hydrophobic and polar interactions, show the environment around each. After binning by backbone ϕ and ψ angles, we clustered each bin with k-means based on calculated map similarities between map-map pairs. Four bins were examined in detail: one in the β-pleat region, two in the right-hand α-helix (RHα) region and one in the left-hand α-helix region of the Ramachandran plot. All regions indicated a common map motif of hydrophobic-hydrophobic interactions along the CA-CB axis, accounting for 62% in the β-pleat bin, about one-third in the two RHα bins and 42% in the LHα bin. Another shared motif shows no interactions along the CA-CB axis; this was uncommon (8%) in β-pleat, but >30% elsewhere. The maps calculated for the two RHα bins are extremely similar (pairwise >0.9787), which suggests that the hydropathic interaction sets or motifs found around each residue are conserved. Altogether, these results are integral to a new paradigm for understanding protein structure and function., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The pneumococcal surface antigen A (PsaA) metal transporter protein provides manganese to bacterial cells. The X-ray crystal structures of PsaA, in both closed (Mn bound) and open (metal free) conformations, were explored with virtual screening to identify potential inhibitors of manganese transport. We pursued three strategies for inhibition: i) targeting a cavity close to the bound Mn to keep the metal in place; ii) targeting the metal-free Mn site to prevent metal uptake; and iii) targeting a potentially druggable allosteric site involving loops that translate between the conformations. Tiered assays were used to test the resulting 170 acquired hits: i) assay 1 tested the compounds' growth inhibition of the TIGR4 S. pneumoniae strain (ΔPsaA mutant control), yielding 80 compounds (MIC≤250 μm); ii) assay 2 tested if the addition of 20 μm Mn to inhibited cell cultures restored growth, yielding 21 compounds; and iii) assay 3 confirmed that the restored bacterial growth was Mn concentration dependent, as was the restoration of ΔPsaA growth, yielding 12 compounds with MICs of 125 μm or greater. It may be possible for a small molecule to inhibit PsaA, but we have not yet identified a compound with exemplary properties., (© 2018 Wiley-VHCA AG, Zurich, Switzerland.)

Meso-Diaminopimelic acid (meso-2,6-diamino-heptanedioic acid, DAP) is an important component of the cell wall of many bacteria. Meso-diaminopimelate dehydrogenase (m-Ddh) is a critical enzyme in the process of converting tetrahydrodipicolinate to DAP. Here, we are proposing that DAP analogs targeting m-Ddh may be considered as potential antibiotics. Four DAP analogs without significant structural change from DAP have been obtained and their inhibitory potencies against m-Ddh from the P. gingivalis strain W83 show significant differences from that of DAP. However, their inhibitory mechanisms as for how simple structural change influences the inhibitory potency remain unknown. Therefore, we employed molecular modeling methods to obtain insight into the inhibitory mechanisms of DAP and analogs with m-Ddh. The predicted binding mode of DAP was highly consistent with the experimental structural data and disclosed the important roles played by the binding pocket residues. According to our predictions, the isoxazoline ring of compounds 1 and 2 and the double bonds in compounds 3 and 4 had distinct influences on these compounds' binding to m-Ddh. This enriched understanding of the inhibitory mechanisms of DAP and these four analogs to m-Ddh has provided new and relevant information for future rational development of potent inhibitors targeting m-Ddh., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Two diastereomeric analogs (1 and 2) of diaminopimelic acid (DAP) bearing an isoxazoline moiety were synthesized and evaluated for their inhibitory activities against meso-diaminopimelate dehydrogenase (m-Ddh) from the periodontal pathogen, Porphyromonas gingivalis. Compound 2 showed promising inhibitory activity against m-Ddh with an IC 50 value of 14.9µM at pH 7.8. The two compounds were further tested for their antibacterial activities against a panel of periodontal pathogens, and compound 2 was shown to be selectively potent to P. gingivalis strains W83 and ATCC 33277 with minimum inhibitory concentration (MIC) values of 773µM and 1.875mM, respectively. Molecular modeling studies revealed that the inversion of chirality at the C-5 position of these compounds was the primary reason for their different biological profiles. Based on these preliminary results, we believe that compound 2 has properties consistent with it being a lead compound for developing novel pathogen selective antibiotics to treat periodontal diseases., (Published by Elsevier Ltd.)

New microtubule depolymerizing agents with potent cytotoxic activities have been prepared with a 5-cyano or 5-oximino group attached to a pyrrole core. The utilization of ortho activation of a bromopyrrole ester to facilitate successful Suzuki-Miyaura cross-coupling reactions was a key aspect of the synthetic methodology. This strategy allows for control of regiochemistry with the attachment of four completely different groups at the 2, 3, 4 and 5 positions of the pyrrole scaffold. Biological evaluations and molecular modeling studies are reported for these examples., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

C-terminal Binding Protein (CtBP) is a transcriptional co-regulator that downregulates the expression of many tumor-suppressor genes. Utilizing a crystal structure of CtBP with its substrate 4-methylthio-2-oxobutyric acid (MTOB) and NAD(+) as a guide, we have designed, synthesized, and tested a series of small molecule inhibitors of CtBP. From our first round of compounds, we identified 2-(hydroxyimino)-3-phenylpropanoic acid as a potent CtBP inhibitor (IC50=0.24μM). A structure-activity relationship study of this compound further identified the 4-chloro- (IC50=0.18μM) and 3-chloro- (IC50=0.17μM) analogues as additional potent CtBP inhibitors. Evaluation of the hydroxyimine analogues in a short-term cell growth/viability assay showed that the 4-chloro- and 3-chloro-analogues are 2-fold and 4-fold more potent, respectively, than the MTOB control. A functional cellular assay using a CtBP-specific transcriptional readout revealed that the 4-chloro- and 3-chloro-hydroxyimine analogues were able to block CtBP transcriptional repression activity. This data suggests that substrate-competitive inhibition of CtBP dehydrogenase activity is a potential mechanism to reactivate tumor-suppressor gene expression as a therapeutic strategy for cancer., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The fundamental pathophysiology of sickle cell disease is predicated by the polymerization of deoxygenated (T-state) sickle hemoglobin (Hb S) into fibers that distort red blood cells into the characteristic sickle shape. The crystal structure of deoxygenated Hb S (DeoxyHb S) and other studies suggest that the polymer is initiated by a primary interaction between the mutation βVal6 from one Hb S molecule, and a hydrophobic acceptor pocket formed by the residues βAla70, βPhe85 and βLeu88 of an adjacent located Hb S molecule. On the contrary, oxygenated or liganded Hb S does not polymerize or incorporate in the polymer. In this paper we present the crystal structure of carbonmonoxy-ligated sickle Hb (COHb S) in the quaternary classical R-state at 1.76Å. The overall structure and the pathological donor and acceptor environments of COHb S are similar to those of the isomorphous CO-ligated R-state normal Hb (COHb A), but differ significantly from DeoxyHb S as expected. More importantly, the packing of COHb S molecules does not show the typical pathological interaction between βVal6 and the βAla70, βPhe85 and βLeu88 hydrophobic acceptor pocket observed in DeoxyHb S crystal. The structural analysis of COHb S, COHb A and DeoxyHb S provides atomic level insight into why liganded hemoglobin does not form a polymer., (Copyright © 2016 Elsevier Inc. All rights reserved.)

A refined model of the colchicine site on tubulin was used to design an improved analog of the pyrrole parent compound, JG-03-14. The optimized compound, NT-7-16, was evaluated in biological assays that confirm that it has potent activities as a new colchicine site microtubule depolymerizer. NT-7-16 exhibits antiproliferative and cytotoxic activities against multiple cancer cell lines, with IC(50) values of 10-16 nM, and it is able to overcome drug resistance mediated by the expression of P-glycoprotein and the βIII isotype of tubulin. NT-7-16 initiated the concentration-dependent loss of cellular microtubules and caused the formation of abnormal mitotic spindles, leading to mitotic accumulation. The direct interaction of NT-7-16 with purified tubulin was confirmed, and it was more potent than combretastatin A-4 in these assays. Binding studies verified that NT-7-16 binds to tubulin within the colchicine site. The antitumor effects of NT-7-16 were evaluated in an MDA-MB-435 xenograft model and it had excellent activity at concentrations that were not toxic. A second compound, NT-9-21, which contains dichloro moieties in place of the 3,5-dibromo substituents of NT-7-16, had a poorer fit within the colchicine site as predicted by modeling and the Hydropathic INTeractions score. Biological evaluations showed that NT-9-21 has 10-fold lower potency than NT-7-16, confirming the modeling predictions. These studies highlight the value of the refined colchicine-site model and identify a new pyrrole-based colchicine-site agent with potent in vitro activities and promising in vivo antitumor actions., (Copyright © 2016 by The American Society for Pharmacology and Experimental Therapeutics.)

Species-specific antimicrobial therapy has the potential to combat the increasing threat of antibiotic resistance and alteration of the human microbiome. We therefore set out to demonstrate the beginning of a pathogen-selective drug discovery method using the periodontal pathogen Porphyromonas gingivalis as a model. Through our knowledge of metabolic networks and essential genes we identified a "druggable" essential target, meso-diaminopimelate dehydrogenase, which is found in a limited number of species. We adopted a high-throughput virtual screen method on the ZINC chemical library to select a group of potential small-molecule inhibitors. Meso-diaminopimelate dehydrogenase from P. gingivalis was first expressed and purified in Escherichia coli then characterized for enzymatic inhibitor screening studies. Several inhibitors with similar structural scaffolds containing a sulfonamide core and aromatic substituents showed dose-dependent inhibition. These compounds were further assayed showing reasonable whole-cell activity and the inhibition mechanism was determined. We conclude that the establishment of this target and screening strategy provides a model for the future development of new antimicrobials.

Sidechain rotamer libraries are obtained through exhaustive statistical analysis of existing crystallographic structures of proteins and have been applied in multiple aspects of structural biology, for example, crystallography of relatively low-resolution structures, in homology model building and in biomolecular NMR. Little is known, however, about the driving forces that lead to the preference or suitability of one rotamer over another. Construction of 3D hydropathic interaction maps for nearly 30,000 tyrosines reveals the environment around each, in terms of hydrophobic (π-π stacking, etc.) and polar (hydrogen bonding, etc.) interactions. After partitioning the tyrosines into backbone-dependent (ϕ, ψ) bins, a map similarity metric based on the correlation coefficient was applied to each map-map pair to build matrices suitable for clustering with k-means. The first bin (-200° ≤ ϕ < -155°; -205° ≤ ψ < -160°), representing 631 tyrosines, reduced to 14 unique hydropathic environments, with most diversity arising from favorable hydrophobic interactions with many different residue partner types. Polar interactions for tyrosine include surprisingly ubiquitous hydrogen bonding with the phenolic OH and a handful of unique environments surrounding the tyrosine backbone. The memberships of all but one of the 14 environments are dominated (>50%) by a single χ(1)/χ(2) rotamer. The last environment has weak or no interactions with the tyrosine ring and its χ(1)/χ(2) rotamer is indeterminate, which is consistent with it being composed of mostly surface residues. Each tyrosine residue attempts to fulfill its hydropathic valence and thus, structural water molecules are seen in a variety of roles throughout protein structure., (© 2015 Wiley Periodicals, Inc.)

A series of 17-cyclopropylmethyl-3,14β-dihydroxy-4,5α-epoxy-6α-(isoquinoline-3'-carboxamido)morphinan (NAQ) analogues were synthesized and pharmacologically characterized to study their structure-activity relationship at the mu opioid receptor (MOR). The competition binding assay showed two-atom spacer and aromatic side chain were optimal for MOR selectivity. Meanwhile, substitutions at the 1'- and/or 4'-position of the isoquinoline ring retained or improved MOR selectivity over the kappa opioid receptor while still possessing above 20-fold MOR selectivity over the delta opioid receptor. In contrast, substitutions at the 6'- and/or 7'-position of the isoquinoline ring reduced MOR selectivity as well as MOR efficacy. Among this series of ligands, compound 11 acted as an antagonist when challenged with morphine in warm-water tail immersion assay and produced less significant withdrawal symptoms compared to naltrexone in morphine-pelleted mice. Compound 11 also antagonized the intracellular Ca(2+) increase induced by DAMGO. Molecular dynamics simulation studies of 11 in three opioid receptors indicated orientation of the 6'-nitro group varied significantly in the different 'address' domains of the receptors and played a crucial role in the observed binding affinities and selectivity. Collectively, the current findings provide valuable insights for future development of NAQ-based MOR selective ligands., (Published by Elsevier Ltd.)

L-Threonine aldolases (TAs), a family of enzymes belonging to the fold-type I pyridoxal 5'-phosphate (PLP) dependent enzymes, play a role in catalyzing the reversible cleavage of l-3-hydroxy-α-amino acids to glycine and the corresponding aldehydes. Threonine aldolases have great biotechnological potential for the syntheses of pharmaceutically relevant drug molecules because of their stereospecificity. The pH-dependency of their catalytic activity, affecting reaction intermediates, led us to study the effect of low-pH on Escherichia coli TA (eTA) structure. We report here a low-pH crystal structure of eTA at 2.1 Å resolution, with a non-covalently bound uncleaved l-serine substrate, and a PLP cofactor bound as an internal aldimine. This structure contrasts with other eTA structures obtained at physiological pH that show products or substrates bound as PLP-external aldimines. The non-productive binding at low-pH is due to an unusual substrate serine binding orientation in which the α-amino group and carboxylate group are in the wrong positions (relative to the active site residues) as a result of protonation of the α-amino group of the serine, as well as the active site histidines, His83 and His126. Protonation of these residues prevents the characteristic nucleophilic attack of the α-amino group of substrate serine on C4' of PLP to form the external aldimine. Our study shows that at low pH the change in charge distribution at the active site can result in substrates binding in a non-productive orientation., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Human cytomegalovirus (CMV) causes significant disease in immunocompromised patients and serious birth defects if acquired in utero. Available CMV antivirals target the viral DNA polymerase, have significant toxicities, and suffer from resistance. New drugs targeting different pathways would be beneficial. The anthraquinone emodin is proposed to inhibit herpes simplex virus by blocking the viral nuclease. Emodin and related anthraquinones are also reported to inhibit CMV. In the present study, emodin reduced CMV infectious yield with an EC50 of 4.9μM but was cytotoxic at concentrations only twofold higher. Related anthraquinones acid blue 40 and alizarin violet R inhibited CMV at only high concentrations (238-265μM) that were also cytotoxic. However, atanyl blue PRL inhibited infectious yield of CMV with an EC50 of 6.3μM, significantly below its 50% cytotoxic concentration of 216μM. Atanyl blue PRL reduced CMV infectivity and inhibited spread. When added up to 1h after infection, it dramatically reduced CMV immediate early protein expression and blocked viral DNA synthesis. However, it had no antiviral activity when added 24h after infection. Interestingly, atanyl blue PRL inhibited nuclease activities of purified CMV UL98 protein with IC50 of 4.5 and 9.3μM. These results indicate that atanyl blue PRL targets very early post-entry events in CMV replication and suggest it may act through inhibition of UL98, making it a novel CMV inhibitor. This compound may provide valuable insights into molecular events that occur at the earliest times post-infection and serve as a lead structure for antiviral development., (Copyright © 2014 Elsevier B.V. All rights reserved.)

Cystalysin from Treponema denticola is a pyridoxal 5'-phosphate dependent lyase that catalyzes the formation of pyruvate, ammonia, and sulfide from cysteine. It is a virulence factor in adult periodontitis because its reaction contributes to hemolysis, which sustains the pathogen. Therefore, it was proposed as a potential antimicrobial target. To identify specific inhibitors by structure-based in silico methods, we first validated the crystal structure of cystalysin as a reliable starting point for the design of ligands. By using single-crystal absorption microspectrophotometry, we found that the enzyme in the crystalline state, with respect to that in solution, exhibits: 1) the same absorption spectra for the catalytic intermediates, 2) a close pKa value for the residue controlling the keto enamine ionization, and 3) similar reactivity with glycine, L-serine, L-methionine, and the nonspecific irreversible inhibitor aminoethoxyvinylglycine. Next, we screened in silico a library of 9357 compounds with the Fingerprints for Ligands and Proteins (FLAP) software, by using the three-dimensional structure of cystalysin as a template. From the library, 17 compounds were selected and experimentally evaluated by enzyme assays and spectroscopic methods. Two compounds were found to competitively inhibit recombinant T. denticola cystalysin, with inhibition constant (Ki ) values of 25 and 37 μM. One of them exhibited a minimum inhibitory concentration (MIC) value of 64 μg mL(-1) on Moraxella catarrhalis ATCC 23246, which proves its ability to cross bacterial membranes., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Characterizing the nature of interaction between proteins that have not been experimentally cocrystallized requires a computational docking approach that can successfully predict the spatial conformation adopted in the complex. In this work, the Hydropathic INTeractions (HINT) force field model was used for scoring docked models in a data set of 30 high-resolution crystallographically characterized "dry" protein-protein complexes and was shown to reliably identify native-like models. However, most current protein-protein docking algorithms fail to explicitly account for water molecules involved in bridging interactions that mediate and stabilize the association of the protein partners, so we used HINT to illuminate the physical and chemical properties of bridging waters and account for their energetic stabilizing contributions. The HINT water Relevance metric identified the "truly" bridging waters at the 30 protein-protein interfaces and we utilized them in "solvated" docking by manually inserting them into the input files for the rigid body ZDOCK program. By accounting for these interfacial waters, a statistically significant improvement of ∼24% in the average hit-count within the top-10 predictions the protein-protein dataset was seen, compared to standard "dry" docking. The results also show scoring improvement, with medium and high accuracy models ranking much better than incorrect ones. These improvements can be attributed to the physical presence of water molecules that alter surface properties and better represent native shape and hydropathic complementarity between interacting partners, with concomitantly more accurate native-like structure predictions., (© 2013 Wiley Periodicals, Inc.)

Cannabinoid receptors are a family of G-protein coupled receptors that are involved in a wide variety of physiological processes and diseases. One of the key regulators that are unique to cannabinoid receptors is the cannabinoid receptor interacting proteins (CRIPs). Among them CRIP1a was found to decrease the constitutive activity of the cannabinoid type-1 receptor (CB1R). The aim of this study is to gain an understanding of the interaction between CRIP1a and CB1R through using different computational techniques. The generated model demonstrated several key putative interactions between CRIP1a and CB1R, including the critical involvement of Lys130 in CRIP1a., (Published by Elsevier Ltd.)

Several biologically active compounds have been identified from Tacca species, including glycosides, diarylheptanoids, saponins, withanolides, and the taccalonolide class of microtubule stabilizers. More recently, two cytotoxic retro-dihydrochalcones named evelynin A (7) and taccabulin A (6) were isolated and their biological activities characterized, including the finding that taccabulin has microtubule destabilizing effects. Here we describe the identification and characterization of five new retro-chalcones, named taccabulins B-E (1-4) and evelynin B (5) from Tacca sp. extracts. Their structures were determined using 1D and 2D NMR as well as mass spectroscopic data and modeled into the colchicine binding site of tubulin. The antiproliferative and microtubule effects of each compound were determined experimentally and found to be well correlated with modeling studies. The isolation and biological characterization of several retro-dihydrochalcones facilitated preliminary structure-activity relationships for this compound class concerning its antiproliferative and microtubule depolymerizing activities.

The importance of protein-protein interactions (PPIs) is becoming increasingly appreciated, as these interactions lie at the core of virtually every biological process. Small molecule modulators that target PPIs are under exploration as new therapies. One of the greatest obstacles faced in crystallographically determining the 3D structures of proteins is coaxing the proteins to form "artificial" PPIs that lead to uniform crystals suitable for X-ray diffraction. This work compares interactions formed naturally, i.e., "biological", with those artificially formed under crystallization conditions or "non-biological". In particular, a detailed analysis of water molecules at the interfaces of high-resolution (≤2.30 Å) X-ray crystal structures of protein-protein complexes, where 140 are biological protein-protein complex structures and 112 include non-biological protein-protein interfaces, was carried out using modeling tools based on the HINT forcefield. Surprisingly few and relatively subtle differences were observed between the two types of interfaces: (i) non-biological interfaces are more polar than biological interfaces, yet there is better organized hydrogen bonding at the latter; (ii) biological associations rely more on water-mediated interactions with backbone atoms compared to non-biological associations; (iii) aromatic/planar residues play a larger role in biological associations with respect to water, and (iv) Lys has a particularly large role at non-biological interfaces. A support vector machines (SVMs) classifier using descriptors from this study was devised that was able to correctly classify 84% of the two interface types., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Highly selective opioid receptor antagonists are essential pharmacological probes in opioid receptor structural characterization and opioid agonist functional studies. Currently, there is no highly selective, nonpeptidyl and reversible mu opioid receptor antagonist available. Among a series of naltrexamine derivatives that have been designed and synthesized, two compounds, NAP and NAQ, were previously identified as novel leads for this purpose based on their in vitro and in vivo pharmacological profiles. Both compounds displayed high binding affinity and selectivity to the mu opioid receptor. To further study the interaction of these two ligands with the three opioid receptors, the recently released opioid receptor crystal structures were employed in docking studies to further test our original hypothesis that the ligands recognize a unique 'address' domain in the mu opioid receptor involving Trp318 that facilitates their selectivity. These modeling results were supported by site-directed mutagenesis studies on the mu opioid receptor, where the mutants Y210A and W318A confirmed the role of the latter in binding. Such work not only enriched the 'message-address' concept, also facilitated our next generation ligand design and development., (Published by Elsevier Ltd.)

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.

αβ-Tubulin colchicine site inhibitors (CSIs) from four scaffolds that we previously tested for antiproliferative activity were modeled to better understand their effect on microtubules. Docking models, constructed by exploiting the SAR of a pyrrole subset and HINT scoring, guided ensemble docking of all 59 compounds. This conformation set and two variants having progressively less structure knowledge were subjected to CoMFA, CoMFA+HINT, and CoMSIA 3D-QSAR analyses. The CoMFA+HINT model (docked alignment) showed the best statistics: leave-one-out q(2) of 0.616, r(2) of 0.949, and r(2)pred (internal test set) of 0.755. An external (tested in other laboratories) collection of 24 CSIs from eight scaffolds were evaluated with the 3D-QSAR models, which correctly ranked their activity trends in 7/8 scaffolds for CoMFA+HINT (8/8 for CoMFA). The combination of SAR, ensemble docking, hydropathic analysis, and 3D-QSAR provides an atomic-scale colchicine site model more consistent with a target structure resolution much higher than the ~3.6 Å available for αβ-tubulin.