Fast 3D shape screening of large chemical databases through alignment-recycling

<p>Abstract</p> <p>Background</p> <p>Large chemical databases require fast, efficient, and simple ways of looking for similar structures. Although such tasks are now fairly well resolved for graph-based similarity queries, they remain an issue for 3D approaches, particularly for those based on 3D shape overlays. Inspired by a recent technique developed to compare molecular shapes, we designed a hybrid methodology, alignment-recycling, that enables efficient retrieval and alignment of structures with similar 3D shapes.</p> <p>Results</p> <p>Using a dataset of more than one million PubChem compounds of limited size (< 28 heavy atoms) and flexibility (< 6 rotatable bonds), we obtained a set of a few thousand diverse structures covering entirely the 3D shape space of the conformers of the dataset. Transformation matrices gathered from the overlays between these diverse structures and the 3D conformer dataset allowed us to drastically (100-fold) reduce the CPU time required for shape overlay. The alignment-recycling heuristic produces results consistent with <it>de novo </it>alignment calculation, with better than 80% hit list overlap on average.</p> <p>Conclusion</p> <p>Overlay-based 3D methods are computationally demanding when searching large databases. Alignment-recycling reduces the CPU time to perform shape similarity searches by breaking the alignment problem into three steps: selection of diverse shapes to describe the database shape-space; overlay of the database conformers to the diverse shapes; and non-optimized overlay of query and database conformers using common reference shapes. The precomputation, required by the first two steps, is a significant cost of the method; however, once performed, querying is two orders of magnitude faster. Extensions and variations of this methodology, for example, to handle more flexible and larger small-molecules are discussed.</p

PubChem3D: Similar conformers

<p>Abstract</p> <p>Background</p> <p>PubChem is a free and open public resource for the biological activities of small molecules. With many tens of millions of both chemical structures and biological test results, PubChem is a sizeable system with an uneven degree of available information. Some chemical structures in PubChem include a great deal of biological annotation, while others have little to none. To help users, PubChem pre-computes "neighboring" relationships to relate similar chemical structures, which may have similar biological function. In this work, we introduce a "Similar Conformers" neighboring relationship to identify compounds with similar 3-D shape and similar 3-D orientation of functional groups typically used to define pharmacophore features.</p> <p>Results</p> <p>The first two diverse 3-D conformers of 26.1 million PubChem Compound records were compared to each other, using a shape Tanimoto (ST) of 0.8 or greater and a color Tanimoto (CT) of 0.5 or greater, yielding 8.16 billion conformer neighbor pairs and 6.62 billion compound neighbor pairs, with an average of 253 "Similar Conformers" compound neighbors per compound. Comparing the 3-D neighboring relationship to the corresponding 2-D neighboring relationship ("Similar Compounds") for molecules such as caffeine, aspirin, and morphine, one finds unique sets of related chemical structures, providing additional significant biological annotation. The PubChem 3-D neighboring relationship is also shown to be able to group a set of non-steroidal anti-inflammatory drugs (NSAIDs), despite limited PubChem 2-D similarity.</p> <p>In a study of 4,218 chemical structures of biomedical interest, consisting of many known drugs, using more diverse conformers per compound results in more 3-D compound neighbors per compound; however, the overlap of the compound neighbor lists per conformer also increasingly resemble each other, being 38% identical at three conformers and 68% at ten conformers. Perhaps surprising is that the average count of conformer neighbors per conformer increases rather slowly as a function of diverse conformers considered, with only a 70% increase for a ten times growth in conformers per compound (a 68-fold increase in the conformer pairs considered).</p> <p>Neighboring 3-D conformers on the scale performed, if implemented naively, is an intractable problem using a modest sized compute cluster. Methodology developed in this work relies on a series of filters to prevent performing 3-D superposition optimization, when it can be determined that two conformers cannot possibly be a neighbor. Most filters are based on Tanimoto equation volume constraints, avoiding incompatible conformers; however, others consider preliminary superposition between conformers using reference shapes.</p> <p>Conclusion</p> <p>The "Similar Conformers" 3-D neighboring relationship locates similar small molecules of biological interest that may go unnoticed when using traditional 2-D chemical structure graph-based methods, making it complementary to such methodologies. The computational cost of 3-D similarity methodology on a wide scale, such as PubChem contents, is a considerable issue to overcome. Using a series of efficient filters, an effective throughput rate of more than 150,000 conformers per second per processor core was achieved, more than two orders of magnitude faster than without filtering.</p

Bcl-2 protein family: Implications in vascular apoptosis and atherosclerosis

Apoptosis has been recognized as a central component in the pathogenesis of atherosclerosis, in addition to the other human pathologies such as cancer and diabetes. The pathophysiology of atherosclerosis is complex, involving both apoptosis and proliferation at different phases of its progression. Oxidative modification of lipids and inflammation differentially regulate the apoptotic and proliferative responses of vascular cells during progression of the atherosclerotic lesion. Bcl-2 proteins act as the major regulators of extrinsic and intrinsic apoptosis signalling pathways and more recently it has become evident that they mediate the apoptotic response of vascular cells in response to oxidation and inflammation either in a provocative or an inhibitory mode of action. Here we address Bcl-2 proteins as major therapeutic targets for the treatment of atherosclerosis and underscore the need for the novel preventive and therapeutic interventions against atherosclerosis, which should be designed in the light of molecular mechanisms regulating apoptosis of vascular cells in atherosclerotic lesions

Mutual Regulation of Bcl-2 Proteins Independent of the BH3 Domain as Shown by the BH3-Lacking Protein Bcl-xAK

The BH3 domain of Bcl-2 proteins was regarded as indispensable for apoptosis induction and for mutual regulation of family members. We recently described Bcl-xAK, a proapoptotic splice product of the bcl-x gene, which lacks BH3 but encloses BH2, BH4 and a transmembrane domain. It remained however unclear, how Bcl-xAK may trigger apoptosis

The Anti-Apoptotic Bcl-xL Protein, a New Piece in the Puzzle of Cytochrome C Interactome

A structural model of the adduct between human cytochrome c and the human anti-apoptotic protein Bcl-xL, which defines the protein-protein interaction surface, was obtained from solution NMR chemical shift perturbation data. The atomic level information reveals key intermolecular contacts identifying new potentially druggable areas on cytochrome c and Bcl-xL. Involvement of residues on cytochrome c other than those in its complexes with electron transfer partners is apparent. Key differences in the contact area also exist between the Bcl-xL adduct with the Bak peptide and that with cytochrome c. The present model provides insights to the mechanism by which cytochrome c translocated to cytosol can be intercepted, so that the apoptosome is not assembled

The MCL-1 BH3 helix is an exclusive MCL-1 inhibitor and apoptosis sensitizer

available in PMC 2011 February 3.MCL-1 has emerged as a major oncogenic and chemoresistance factor. A screen of stapled peptide helices identified the MCL-1 BH3 domain as selectively inhibiting MCL-1 among the related anti-apoptotic Bcl-2 family members, providing insights into the molecular determinants of binding specificity and a new approach for sensitizing cancer cells to apoptosis.National Institutes of Health (U.S.) (NIH award 5RO1GM084181)National Institutes of Health (U.S.) (NIH grant 5P01CA92625)National Institutes of Health (U.S.) (Ruth L. Kirschstein National Research Service Award 1F31CA144566)Burroughs Wellcome Fund (Career Award

Crystal Structures of T. b. rhodesiense Adenosine Kinase Complexed with Inhibitor and Activator: Implications for Catalysis and Hyperactivation

Recently, we discovered that 4-[5-(4-phenoxyphenyl)-2H-pyrazol-3-yl]morpholine (compound 1) and its derivatives exhibit specific antitrypanosomal activity toward T. b. rhodesiense, the causative agent of the acute form of HAT. We found that compound 1 would target the parasite adenosine kinase (TbrAK), an important enzyme of the purine salvage pathway, by acting via hyperactivation of the enzyme. This represents a novel and hitherto unexplored strategy for the development of trypanocides. These findings prompted us to investigate the mechanism of action at the molecular level. The present study reports the first three-dimensional crystal structures of TbrAK in complex with the bisubstrate inhibitor AP5A, and in complex with the activator (compound 1). The subsequent structural analysis sheds light on substrate and activator binding, and gives insight into the possible mechanism leading to hyperactivation. Further structure-activity relationships in terms of TbrAK activation properties support the observed binding mode of compound 1 in the crystal structure and may open the field for subsequent optimization of this compound series

Selective regulation of IP3-receptor-mediated Ca2+ signaling and apoptosis by the BH4 domain of Bcl-2 versus Bcl-Xl

Antiapoptotic B-cell lymphoma 2 (Bcl-2) targets the inositol 1,4,5-trisphosphate receptor (IP3R) via its BH4 domain, thereby suppressing IP3R Ca2+-flux properties and protecting against Ca2+-dependent apoptosis. Here, we directly compared IP3R inhibition by BH4-Bcl-2 and BH4-Bcl-Xl. In contrast to BH4-Bcl-2, BH4-Bcl-Xl neither bound the modulatory domain of IP3R nor inhibited IP3-induced Ca2+ release (IICR) in permeabilized and intact cells. We identified a critical residue in BH4-Bcl-2 (Lys17) not conserved in BH4-Bcl-Xl (Asp11). Changing Lys17 into Asp in BH4-Bcl-2 completely abolished its IP3R-binding and -inhibitory properties, whereas changing Asp11 into Lys in BH4-Bcl-Xl induced IP3R binding and inhibition. This difference in IP3R regulation between BH4-Bcl-2 and BH4-Bcl-Xl controls their antiapoptotic action. Although both BH4-Bcl-2 and BH4-Bcl-Xl had antiapoptotic activity, BH4-Bcl-2 was more potent than BH4-Bcl-Xl. The effect of BH4-Bcl-2, but not of BH4-Bcl-Xl, depended on its binding to IP(3)Rs. In agreement with the IP3R-binding properties, the antiapoptotic activity of BH4-Bcl-2 and BH4-Bcl-Xl was modulated by the Lys/Asp substitutions. Changing Lys17 into Asp in full-length Bcl-2 significantly decreased its binding to the IP3R, its ability to inhibit IICR and its protection against apoptotic stimuli. A single amino-acid difference between BH4-Bcl-2 and BH4-Bcl-Xl therefore underlies differential regulation of IP(3)Rs and Ca2+-driven apoptosis by these functional domains. Mutating this residue affects the function of Bcl-2 in Ca2+ signaling and apoptosis

Dynamic Interaction of cBid with Detergents, Liposomes and Mitochondria

The BH3-only protein Bid plays a key role in the induction of mitochondrial apoptosis, but its mechanism of action is still not completely understood. Here we studied the two main activation events of Bid: Caspase-8 cleavage and interaction with the membrane bilayer. We found a striking reversible behaviour of the dissociation-association events between the Bid fragments p15 and p7. Caspase-8 cleavage does not induce per se separation of the two Bid fragments, which remain in a stable complex resembling the full length Bid. Detergents trigger a complete dissociation, which can be fully reversed by detergent removal in a range of protein concentrations from 100 µM down to 500 nM. Incubation of cBid with cardiolipin-containing liposomes leads to partial dissociation of the complex. Only p15 (tBid) fragments are found at the membrane, while p7 shows no tendency to interact with the bilayer, but complete removal of p7 strongly increases the propensity of tBid to become membrane-associated. Despite the striking structural similarities of inactive Bid and Bax, Bid does not form oligomers and reacts differently in the presence of detergents and membranes, highlighting clear differences in the modes of action of the two proteins. The partial dissociation of cBid triggered by the membrane is suggested to depend on the strong and specific interaction between p15 and p7. The reversible disassembly and re-assembly of the cBid molecules at the membrane was as well proven by EPR using spin labeled cBid in the presence of isolated mitochondria. The observed dynamic dissociation of the two Bid fragments could allow the assistance to the pore-forming Bax to occur repeatedly and may explain the proposed “hit-and-run" mode of action of Bid at the bilayer

Let’s not forget tautomers

A compound exhibits tautomerism if it can be represented by two structures that are related by an intramolecular movement of hydrogen from one atom to another. The different tautomers of a molecule usually have different molecular fingerprints, hydrophobicities and pKa’s as well as different 3D shape and electrostatic properties; additionally, proteins frequently preferentially bind a tautomer that is present in low abundance in water. As a result, the proper treatment of molecules that can tautomerize, ~25% of a database, is a challenge for every aspect of computer-aided molecular design. Library design that focuses on molecular similarity or diversity might inadvertently include similar molecules that happen to be encoded as different tautomers. Physical property measurements might not establish the properties of individual tautomers with the result that algorithms based on these measurements may be less accurate for molecules that can tautomerize—this problem influences the accuracy of filtering for library design and also traditional QSAR. Any 2D or 3D QSAR analysis must involve the decision of if or how to adjust the observed Ki or IC50 for the tautomerization equilibria. QSARs and recursive partitioning methods also involve the decision as to which tautomer(s) to use to calculate the molecular descriptors. Docking virtual screening must involve the decision as to which tautomers to include in the docking and how to account for tautomerization in the scoring. All of these decisions are more difficult because there is no extensive database of measured tautomeric ratios in both water and non-aqueous solvents and there is no consensus as to the best computational method to calculate tautomeric ratios in different environments