Yen-Fu Chen

iGEMDOCK: a graphical environment of enhancing GEMDOCK using pharmacological interactions and post-screening analysis

BackgroundPharmacological interactions are useful for understanding ligand binding mechanisms of a therapeutic target. These interactions are often inferred from a set of active compounds that were acquired experimentally. Moreover, most docking programs loosely coupled the stages (binding-site and ligand preparations, virtual screening, and post-screening analysis) of structure-based virtual screening (VS). An integrated VS environment, which provides the friendly interface to seamlessly combine these VS stages and to identify the pharmacological interactions directly from screening compounds, is valuable for drug discovery.ResultsWe developed an easy-to-use graphic environment, i GEMDOCK, integrating VS stages (from preparations to post-screening analysis). For post-screening analysis, i GEMDOCK provides biological insights by deriving the pharmacological interactions from screening compounds without relying on the experimental data of active compounds. The pharmacological interactions represent conserved interacting residues, which often form binding pockets with specific physico-chemical properties, to play the essential functions of a target protein. Our experimental results show that the pharmacological interactions derived by i GEMDOCK are often hot spots involving in the biological functions. In addition, i GEMDOCK provides the visualizations of the protein-compound interaction profiles and the hierarchical clustering dendrogram of the compounds for post-screening analysis.ConclusionsWe have developed i GEMDOCK to facilitate steps from preparations of target proteins and ligand libraries toward post-screening analysis. i GEMDOCK is especially useful for post-screening analysis and inferring pharmacological interactions from screening compounds. We believe that i GEMDOCK is useful for understanding the ligand binding mechanisms and discovering lead compounds. i GEMDOCK is available at http://gemdock.life.nctu.edu.tw/dock/igemdock.php.

Aurintricarboxylic acid inhibits influenza virus neuraminidase

Abstract There is a continuing threat that the highly pathogenic avian influenza virus will cause future influenza pandemics. In this study, we screened a library of compounds that are biologically active and structurally diverse for inhibitory activity against influenza neuraminidase (NA). We found that aurintricarboxylic acid (ATA) is a potent inhibitor of NA activity of both group-1 and group-2 influenza viruses with IC50s (effective concentration to inhibit NA activity by 50%) values at low micromolar concentrations. ATA was equally potent in inhibiting the NA activity derived from wild-type NA and its H274Y mutant which renders NA resistance to inhibition by oseltamivir. Although ATA is structurally distinct from sialic acid, molecular modeling experiments suggested that ATA binds to NA at the enzyme’s substrate binding site. These results indicate that ATA may be a good starting material for the design of a novel class of NA inhibitors for the treatment influenza viruses.

Time‐Evolving Self‐Organization and Autonomous Structural Adaptation of Cobalt(II)–Organic Framework Materials with scu and pts Nets

Self-organization is a process, in which an internal system spontaneously opens a new route to increase system complexity without being guided by an external source. The concept of self-organization is central to the understanding of living organisms, biominerals, and new supramolecular materials. For chemistry, self-organizing equilibrium conditions can be controlled by changing a few critical factors (concentration, template, pH, temperature, solvent system, etc.) to generate desirable compounds. However, these explorations seem not to be completely applied in a few particular supramolecular systems. Inspired by biology, to construct a high-order architecture from individual building components, various driving forces may competitively predominate at certain stages of the self-assembly process. A subtle thermodynamic/kinetic balance may control and tune the materials growth delicately. Namely, self-organization processes can be operative if the building components are sufficient and in close proximity, under suitable conditions. If the supply of building units is depleted or reduced, the original equilibrium conditions will change, and a new self-organization process will take place. These intriguing phenomena of self-organization are triggered by an internal stimulus and seem to be easily understood in biology, but the phenomena has not been addressed in the synthesis system of metal–organic framework (MOF) materials. As part of our ongoing efforts in the design and synthesis of functional crystalline materials, we report herein on an intriguing supramolecular system that involves a distinct self-organization process, in which the product structures adapt to autonomous dynamic changes in the ratio of build-

SiMMap: a web server for inferring site-moiety map to recognize interaction preferences between protein pockets and compound moieties

The protein–ligand interacting mechanism is essential to biological processes and drug discovery. The SiMMap server statistically derives site-moiety map with several anchors, which describe the relationship between the moiety preferences and physico-chemical properties of the binding site, from the interaction profiles between query target protein and its docked (or co-crystallized) compounds. Each anchor includes three basic elements: a binding pocket with conserved interacting residues, the moiety composition of query compounds and pocket–moiety interaction type (electrostatic, hydrogen bonding or van der Waals). We provide initial validation of the site-moiety map on three targets, thymidine kinase, and estrogen receptors of antagonists and agonists. Experimental results show that an anchor is often a hot spot and the site-moiety map can help to assemble potential leads by optimal steric, hydrogen bonding and electronic moieties. When a compound highly agrees with anchors of site-moiety map, this compound often activates or inhibits the target protein. We believe that the site-moiety map is useful for drug discovery and understanding biological mechanisms. The SiMMap web server is available at http://simfam.life.nctu.edu.tw/.

Core Site-Moiety Maps Reveal Inhibitors and Binding Mechanisms of Orthologous Proteins by Screening Compound Libraries

Members of protein families often share conserved structural subsites for interaction with chemically similar moieties despite low sequence identity. We propose a core site-moiety map of multiple proteins (called CoreSiMMap) to discover inhibitors and mechanisms by profiling subsite-moiety interactions of immense screening compounds. The consensus anchor, the subsite-moiety interactions with statistical significance, of a CoreSiMMap can be regarded as a “hot spot” that represents the conserved binding environments involved in biological functions. Here, we derive the CoreSiMMap with six consensus anchors and identify six inhibitors (IC50<8.0 µM) of shikimate kinases (SKs) of Mycobacterium tuberculosis and Helicobacter pylori from the NCI database (236,962 compounds). Studies of site-directed mutagenesis and analogues reveal that these conserved interacting residues and moieties contribute to pocket-moiety interaction spots and biological functions. These results reveal that our multi-target screening strategy and the CoreSiMMap can increase the accuracy of screening in the identification of novel inhibitors and subsite-moiety environments for elucidating the binding mechanisms of targets.

Structures of Helicobacter pylori shikimate kinase reveal a selective inhibitor-induced-fit mechanism.

Shikimate kinase (SK), which catalyzes the specific phosphorylation of the 3-hydroxyl group of shikimic acid in the presence of ATP, is the enzyme in the fifth step of the shikimate pathway for biosynthesis of aromatic amino acids. This pathway is present in bacteria, fungi, and plants but absent in mammals and therefore represents an attractive target pathway for the development of new antimicrobial agents, herbicides, and antiparasitic agents. Here we investigated the detailed structure–activity relationship of SK from Helicobacter pylori (HpSK). Site-directed mutagenesis and isothermal titration calorimetry studies revealed critical conserved residues (D33, F48, R57, R116, and R132) that interact with shikimate and are therefore involved in catalysis. Crystal structures of HpSK·SO4, R57A, and HpSK•shikimate-3-phosphate•ADP show a characteristic three-layer architecture and a conformationally elastic region consisting of F48, R57, R116, and R132, occupied by shikimate. The structure of the inhibitor complex, E114A•162535, was also determined, which revealed a dramatic shift in the elastic LID region and resulted in conformational locking into a distinctive form. These results reveal considerable insight into the active-site chemistry of SKs and a selective inhibitor-induced-fit mechanism.

Orientation-dependent conductance study of pentacene nanocrystals by conductive atomic force microscopy

Oriented pentacene nanocrystals with long molecular axis either parallel or perpendicular to a Au substrate were prepared on a bare Au surface or a self-assembled monolayer (SAM)-modified Au surface, respectively. The conductance across the differently oriented pentacene crystals were measured by conductive atomic force microscopy in a similar device configuration of Au/SAM/pentacene/Au-tip and Au/pentacene/SAM-modified-Au-tip, respectively. Rectifying current was observed depending on the location of the SAM in the device. With an average thickness of 50nm, the conductance along the C–H⋯π stacking direction (a-b plane) was nearly five orders of magnitude larger than along the layer direction (c axis).

Pathway-based screening strategy for multitarget inhibitors of diverse proteins in metabolic pathways.

Many virtual screening methods have been developed for identifying single-target inhibitors based on the strategy of “one–disease, one–target, one–drug”. The hit rates of these methods are often low because they cannot capture the features that play key roles in the biological functions of the target protein. Furthermore, single-target inhibitors are often susceptible to drug resistance and are ineffective for complex diseases such as cancers. Therefore, a new strategy is required for enriching the hit rate and identifying multitarget inhibitors. To address these issues, we propose the pathway-based screening strategy (called PathSiMMap) to derive binding mechanisms for increasing the hit rate and discovering multitarget inhibitors using site-moiety maps. This strategy simultaneously screens multiple target proteins in the same pathway; these proteins bind intermediates with common substructures. These proteins possess similar conserved binding environments (pathway anchors) when the product of one protein is the substrate of the next protein in the pathway despite their low sequence identity and structure similarity. We successfully discovered two multitarget inhibitors with IC50 of <10 µM for shikimate dehydrogenase and shikimate kinase in the shikimate pathway of Helicobacter pylori. Furthermore, we found two selective inhibitors (IC50 of <10 µM) for shikimate dehydrogenase using the specific anchors derived by our method. Our experimental results reveal that this strategy can enhance the hit rates and the pathway anchors are highly conserved and important for biological functions. We believe that our strategy provides a great value for elucidating protein binding mechanisms and discovering multitarget inhibitors.

TSCC: Two-Stage Combinatorial Clustering for virtual screening using protein-ligand interactions and physicochemical features

BackgroundThe increasing numbers of 3D compounds and protein complexes stored in databases contribute greatly to current advances in biotechnology, being employed in several pharmaceutical and industrial applications. However, screening and retrieving appropriate candidates as well as handling false positives presents a challenge for all post-screening analysis methods employed in retrieving therapeutic and industrial targets.ResultsUsing the TSCC method, virtually screened compounds were clustered based on their protein-ligand interactions, followed by structure clustering employing physicochemical features, to retrieve the final compounds. Based on the protein-ligand interaction profile (first stage), docked compounds can be clustered into groups with distinct binding interactions. Structure clustering (second stage) grouped similar compounds obtained from the first stage into clusters of similar structures; the lowest energy compound from each cluster being selected as a final candidate.ConclusionBy representing interactions at the atomic-level and including measures of interaction strength, better descriptions of protein-ligand interactions and a more specific analysis of virtual screening was achieved. The two-stage clustering approach enhanced our post-screening analysis resulting in accurate performances in clustering, mining and visualizing compound candidates, thus, improving virtual screening enrichment.

Direct visualization of triplex DNA molecular dynamics by fluorescence resonance energy transfer and atomic force microscopy measurements

We have detected the dynamics of 17-mer DNA triplex dissociation mechanism at the molecular level. Fluorescence resonance energy transfer (FRET) was used as an indicator of intermolecular interaction in nanometer range, whereas atomic force microscopy (AFM) was employed to address single molecule with sub-angstrom precision. The maximum rupture force of DNA triplex was found at pH 4.65, consistent with macroscopic observations. These results indicated that the FRET together with an AFM detection system could be used to reveal the DNA triplex interaction in nanometer scale unambiguously.