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

Natural product drupacine acting on a novel herbicidal target shikimate dehydrogenase

Natural products are one of the important sources for the creation of new pesticides. Drupacine ((1R,11S,12S,13R,15S)-13-methoxy-5,7,21-trioxa-19-azahexacyclo[11.7.1.02,10.04,8.011,15.015,19]henicosa-2,4(8),9-trien-12-ol), isolated from Cephalotaxus sinensis (Chinese plum-yew), is a potent herbicidal compound containing an oxo-bridged oxygen bond structure. However, its molecular target still remains unknown. In this study, the targets of drupacine in Amaranthus retroflexus were identified by combining drug affinity responsive target stability (DARTS), cellular thermal shift assay coupled with mass spectrometry (CETSA MS), RNA-seq transcriptomic, and TMT proteomic analyses. Fifty-one and sixty-eight main binding proteins were identified by DARTS and CETSA MS, respectively, including nine co-existing binding proteins. In drupacine-treated A. retroflexus seedlings we identified 1389 up-regulated genes and 442 down-regulated genes, 34 up-regulated proteins, and 194 down-regulated proteins, respectively. Combining the symptoms and the biochemical profiles, Profilin, Shikimate dehydrogenase (SkDH), and Zeta-carotene desaturase were predicted to be the drupacine potential target proteins. At the same time, drupacine was found to bind SkDH stronger by molecular docking, and its inhibition on ArSkDH increased with the treatment concentration increase. Our results suggest that the molecular target of drupacine is SkDH, a new herbicide target, which lay a foundation for the rational design of herbicides based on new targets from natural products and enrich the target resources for developing green herbicides.

Shikimate Kinase Inhibitors: An Update on Promising Strategy against Mycobacterium tuberculosis

Humanity has been battling with tuberculosis (TB) for a long period, and despite the availability of drugs well-known to act against the deadly microbe, the menace is still very far from reaching its end. Moreover, problems related to TB chemotherapy, such as lengthy treatment periods leading to poor patient compliance, increasing drug resistance, and association with another deadlier disease HIV-AIDS, make the situation alarming, thereby pressing the need for the discovery of new potent drugs urgently. Therefore, a drug target that is essential for survival and exclusive to M. tuberculosis presents a promising platform to explore novel molecules against the microorganism for better pathogen clearance with minimal toxicity. The shikimate pathway that leads to the synthesis of essential aromatic amino acids is one such attractive target. Shikimate kinase, the fifth enzyme of this pathway, converts shikimate to shikimate-3-phosphate by using ATP as a cosubstrate. Targeting shikimate kinase could be an effective strategy in light of its essentiality and absence of any homologue in mammals. This review discusses different strategies adopted for discovering novel compounds or scaffolds targeting M. tuberculosis shikimate kinase (MtSK) in vitro. The application of substrate analogues, their structure, and ligand-based approach for screening a library of anti-mycobacterial compounds, marine-derived molecules, and commercially available libraries have yielded promising MtSK inhibitors exhibiting micro-molar activities. To develop these leads into future drugs with minimum off-target effects on the host microenvironment, the molecules need to be structurally optimized for improved activities against enzymes and whole-cell organisms.

A dual-target herbicidal inhibitor of lysine biosynthesis

Herbicides with novel modes of action are urgently needed to safeguard global agricultural industries against the damaging effects of herbicide-resistant weeds. We recently developed the first herbicidal inhibitors of lysine biosynthesis, which provided proof-of-concept for a promising novel herbicide target. In this study, we expanded upon our understanding of the mode of action of herbicidal lysine biosynthesis inhibitors. We previously postulated that these inhibitors may act as proherbicides. Here, we show this is not the case. We report an additional mode of action of these inhibitors, through their inhibition of a second lysine biosynthesis enzyme, and investigate the molecular determinants of inhibition. Furthermore, we extend our herbicidal activity analyses to include a weed species of global significance.

Translatable pathways classification (TransPath-C) for inferring processes germane to human biology from animal studies data: example application in neurobiology

How to translate insights gained from studies in one organismal species for what is most likely to be germane in another species, such as from mice to humans, is a ubiquitous challenge in basic biology as well as biomedicine. This is an especially difficult problem when there are few molecular features that are obviously important in both species for a given phenotype of interest. Neuropathologies are a prominent realm of this complication. Schizophrenia is complex psychiatric disorder that affects 1% of the population. Many genetic factors have been proposed to drive the development of schizophrenia, and the 22q11 microdeletion (MD) syndrome has been shown to dramatically increase this risk. Due to heterogeneity of presentation of symptoms, diagnosis and formulation of treatment options for patients can often be delayed, and there is an urgent need for novel therapeutics directed toward the treatment of schizophrenia. Here, we present a novel computational approach, Translational Pathways Classification (TransPath-C), that can be used to identify shared pathway dysregulation between mouse models and human schizophrenia cohorts. This method uses variation of pathway activation in the mouse model to predict both mouse and human disease phenotype. Analysis of shared dysregulated pathways called out by both the mouse and human classifiers of TransPath-C can identify pathways that can be targeted in both preclinical and human cohorts of schizophrenia. In application to the 22q11 MD mouse model, our findings suggest that PAR1 pathway activation found upregulated in this mouse phenotype is germane for the corresponding human schizophrenia cohort such that inhibition of PAR1 may offer a novel therapeutic target.

Uncovering Flexible Active Site Conformations of SARS-CoV-2 3CL Proteases through Protease Pharmacophore Clusters and COVID-19 Drug Repurposing

The infectious SARS-CoV-2 causes COVID-19, which is now a global pandemic. Aiming for effective treatments, we focused on the key drug target, the viral 3C-like (3CL) protease. We modeled a big dataset with 42 SARS-CoV-2 3CL protease-ligand complex structures from ∼98.7% similar SARS-CoV 3CL protease with abundant complex structures. The diverse flexible active site conformations identified in the dataset were clustered into six protease pharmacophore clusters (PPCs). For the PPCs with distinct flexible protease active sites and diverse interaction environments, we identified pharmacophore anchor hotspots. A total of 11 "PPC consensus anchors" (a distinct set observed in each PPC) were observed, of which three "PPC core anchors" EHV2, HV1, and V3 are strongly conserved across PPCs. The six PPC cavities were then applied in virtual screening of 2122 FDA drugs for repurposing, using core anchor-derived "PPC scoring S" to yield seven drug candidates. Experimental testing by SARS-CoV-2 3CL protease inhibition assay and antiviral cytopathic effect assays discovered active hits, Boceprevir and Telaprevir (HCV drugs) and Nelfinavir (HIV drug). Specifically, Boceprevir showed strong protease inhibition with micromolar IC₅₀ of 1.42 μM and an antiviral activity with EC₅₀ of 49.89 μM, whereas Telaprevir showed moderate protease inhibition only with an IC₅₀ of 11.47 μM. Nelfinavir solely showed antiviral activity with a micromolar EC₅₀ value of 3.28 μM. Analysis of binding mechanisms of protease inhibitors revealed the role of PPC core anchors. Our PPCs revealed the flexible protease active site conformations, which successfully enabled drug repurposing.

Targeting shikimate pathway: In silico analysis of phosphoenolpyruvate derivatives as inhibitors of EPSP synthase and DAHP synthase

The shikimate pathway consists of seven enzymatic steps involved in the conversion of erythrose-4-phosphate and phosphoenolpyruvate to chorismate and also responsible to the production of aromatic amino acids, such as phenylalanine, tyrosine, and tryptophan which are essential to the bacterial metabolism. The 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase (DAHPS) and 5-enolpyruvylshikimate 3-phosphate synthase (EPSPS) catalyze important steps in the shikimate pathway using as substrate the phosphoenolpyruvate (PEP). Due to the importance of PEP in shikimate pathway, its structure has been investigated to develop new bioinspired competitive inhibitors against DAHPS and EPSPS. In the present study, we perform a literature survey of 28 PEP derivatives, then we analyzed the selectivity and affinity of these compounds against the EPSPS and DAHPS structures using consensual molecular docking, pharmacophore prediction, molecular dynamics (MD) simulations, and binding free energy calculations. Here, we propose consistent binding modes of the selected ligands and indicate that their structures show interesting pharmacophoric properties related to multi-targets inhibitors for both enzymes. Our computational results are supported by previous experimental findings related to the interactions of PEP derivatives with DAHPS and EPSPS structures.

A site-moiety map and virtual screening approach for discovery of novel 5-LOX inhibitors

The immune system works in conjunction with inflammation. Excessive inflammation underlies various human diseases, such as asthma, diabetes and heart disease. Previous studies found that 5-lipoxygenase (5-LOX) plays a crucial role in metabolizing arachidonic acid into inflammatory mediators and is a potential therapeutic target. In this study, we performed an in silico approach to establish a site-moiety map (SiMMap) to screen for new 5-LOX inhibitors. The map is composed of several anchors that contain key residues, moiety preferences, and their interaction types (i.e., electrostatic (E), hydrogen-bonding (H), and van der Waals (V) interactions) within the catalytic site. In total, we identified one EH, one H, and five V anchors, within the 5-LOX catalytic site. Based on the SiMMap, three 5-LOX inhibitors (YS1, YS2, and YS3) were identified. An enzyme-based assay validated inhibitory activity of YS1, YS2, and YS3 against 5-LOX with an IC₅₀ value of 2.7, 4.2, and 5.3 μM, respectively. All three inhibitors significantly decrease LPS-induced TNF-α and IL-6 production, which suggests its potential use an anti-inflammatory agent. In addition, the identified 5-LOX inhibitors contain a novel scaffold. The discovery of these inhibitors presents an opportunity for designing specific anti-inflammatory drugs.

Zika Virus NS3 Protease Pharmacophore Anchor Model and Drug Discovery

Zika virus (ZIKV) of the flaviviridae family, is the cause of emerging infections characterized by fever, Guillain-Barré syndrome (GBS) in adults and microcephaly in newborns. There exists an urgent unmet clinical need for anti-ZIKV drugs for the treatment of infected individuals. In the current work, we aimed at the promising virus drug target, ZIKV NS3 protease and constructed a Pharmacophore Anchor (PA) model for the active site. The PA model reveals a total of 12 anchors (E, H, V) mapped across the active site subpockets. We further identified five of these anchors to be critical core anchors (CEH1, CH3, CH7, CV1, CV3) conserved across flaviviral proteases. The ZIKV protease PA model was then applied in anchor-enhanced virtual screening yielding 14 potential antiviral candidates, which were tested by in vitro assays. We discovered FDA drugs Asunaprevir and Simeprevir to have potent anti-ZIKV activities with EC₅₀ values 4.7 µM and 0.4 µM, inhibiting the viral protease with IC₅₀ values 6.0 µM and 2.6 µM respectively. Additionally, the PA model anchors aided in the exploration of inhibitor binding mechanisms. In conclusion, our PA model serves as a promising guide map for ZIKV protease targeted drug discovery and the identified ‘previr’ FDA drugs are promising for anti-ZIKV treatments.

Alternative splicing in human cancer cells is modulated by the amiloride derivative 3,5-diamino-6-chloro-N-(N-(2,6-dichlorobenzoyl)carbamimidoyl)pyrazine-2-carboxide

Alternative splicing (AS) is a process that enables the generation of multiple protein isoforms with different biological properties from a single mRNA. Cancer cells often use the maneuverability conferred by AS to produce proteins that contribute to growth and survival. In our previous studies, we identified that amiloride modulates AS in cancer cells. However, the effective concentration of amiloride required to modulate AS is too high for use in cancer treatment. In this study, we used computational algorithms to screen potential amiloride derivatives for their ability to regulate AS in cancer cells. We found that 3,5-diamino-6-chloro-N-(N-(2,6-dichlorobenzoyl)carbamimidoyl)pyrazine-2-carboxamide (BS008) can regulate AS of apoptotic gene transcripts, including HIPK3, SMAC, and BCL-X, at a lower concentration than amiloride. This splicing regulation involved various splicing factors, and it was accompanied by a change in the phosphorylation state of serine/arginine-rich proteins (SR proteins). RNA sequencing was performed to reveal that AS of many other apoptotic gene transcripts, such as AATF, ATM, AIFM1, NFKB1, and API5, was also modulated by BS008. In vivo experiments further indicated that treatment of tumor-bearing mice with BS008 resulted in a marked decrease in tumor size. BS008 also had inhibitory effects in vitro, either alone or in a synergistic combination with the cytotoxic chemotherapeutic agents sorafenib and nilotinib. BS008 enabled sorafenib dose reduction without compromising antitumor activity. These findings suggest that BS008 may possess therapeutic potential for cancer treatment.

Synthesis, biological activity and molecular modelling studies of shikimic acid derivatives as inhibitors of the shikimate dehydrogenase enzyme of Escherichia coli

Shikimic acid (SA) pathway is the common route used by bacteria, plants, fungi, algae, and certain Apicomplexa parasites for the biosynthesis of aromatic amino acids and other secondary metabolites. As this essential pathway is absent in mammals designing inhibitors against implied enzymes may lead to the development of antimicrobial and herbicidal agents harmless to humans. Shikimate dehydrogenase (SDH) is the fourth enzyme of the SA pathway. In this contribution, a series of SA amide derivatives were synthesised and evaluated for in vitro SDH inhibition and antibacterial activity against Escherichia coli. All tested compounds showed to be mixed type inhibitors; diamide derivatives displayed more inhibitory activity than synthesised monoamides. Among the evaluated compounds, molecules called 4a and 4b were the most active derivatives with IC50 588 and 589 µM, respectively. Molecular modelling studies suggested two different binding modes of monoamide and diamide derivatives to the SDH enzyme of E. coli.

Pharmacophore anchor models of flaviviral NS3 proteases lead to drug repurposing for DENV infection

Background

Viruses of the flaviviridae family are responsible for some of the major infectious viral diseases around the world and there is an urgent need for drug development for these diseases. Most of the virtual screening methods in flaviviral drug discovery suffer from a low hit rate, strain-specific efficacy differences, and susceptibility to resistance. It is because they often fail to capture the key pharmacological features of the target active site critical for protein function inhibition. So in our current work, for the flaviviral NS3 protease, we summarized the pharmacophore features at the protease active site as anchors (subsite-moiety interactions).

Results

For each of the four flaviviral NS3 proteases (i.e., HCV, DENV, WNV, and JEV), the anchors were obtained and summarized into ‘Pharmacophore anchor (PA) models’. To capture the conserved pharmacophore anchors across these proteases, were merged the four PA models. We identified five consensus core anchors (CEH1, CH3, CH7, CV1, CV3) in all PA models, represented as the “Core pharmacophore anchor (CPA) model” and also identified specific anchors unique to the PA models. Our PA/CPA models complied with 89 known NS3 protease inhibitors. Furthermore, we proposed an integrated anchor-based screening method using the anchors from our models for discovering inhibitors. This method was applied on the DENV NS3 protease to screen FDA drugs discovering boceprevir, telaprevir and asunaprevir as promising anti-DENV candidates. Experimental testing against DV2-NGC virus by in-vitro plaque assays showed that asunaprevir and telaprevir inhibited viral replication with EC₅₀ values of 10.4 μM & 24.5 μM respectively. The structure-anchor-activity relationships (SAAR) showed that our PA/CPA model anchors explained the observed in-vitro activities of the candidates. Also, we observed that the CEH1 anchor engagement was critical for the activities of telaprevir and asunaprevir while the extent of inhibitor anchor occupation guided their efficacies.

Conclusion

These results validate our NS3 protease PA/CPA models, anchors and the integrated anchor-based screening method to be useful in inhibitor discovery and lead optimization, thus accelerating flaviviral drug discovery.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1957-5) contains supplementary material, which is available to authorized users.

Crystal structure of 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase from the ESKAPE pathogen Acinetobacter baumannii

The enzyme 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase catalyzes the sixth step of the seven-step shikimate pathway. Chorismate, the product of the pathway, is a precursor for the biosynthesis of aromatic amino acids, siderophores and metabolites such as folate, ubiquinone and vitamin K. The shikimate pathway is present in bacteria, fungi, algae, plants and apicomplexan parasites, but is absent in humans. The EPSP synthase enzyme produces 5-enolpyruvylshikimate 3-phosphate and phosphate from phosphoenolpyruvate and shikimate 3-phosphate via a transferase reaction, and is the target of the herbicide glyphosate. The Acinetobacter baumannii gene encoding EPSP synthase, aroA, has previously been demonstrated to be essential during host infection for the growth and survival of this clinically important drug-resistant ESKAPE pathogen. Prephenate dehydrogenase is also encoded by the bifunctional A. baumannii aroA gene, but its activity is dependent upon EPSP synthase since it operates downstream of the shikimate pathway. As part of an effort to evaluate new antimicrobial targets, recombinant A. baumannii EPSP (AbEPSP) synthase, comprising residues Ala301-Gln756 of the aroA gene product, was overexpressed in Escherichia coli, purified and crystallized. The crystal structure, determined to 2.37 Å resolution, is described in the context of a potential antimicrobial target and in comparison to EPSP synthases that are resistant or sensitive to the herbicide glyphosate.

Modelling of compound combination effects and applications to efficacy and toxicity: state-of-the-art, challenges and perspectives

The development of treatments involving combinations of drugs is a promising approach towards combating complex or multifactorial disorders. However, the large number of compound combinations that can be generated, even from small compound collections, means that exhaustive experimental testing is infeasible. The ability to predict the behaviour of compound combinations in biological systems, whittling down the number of combinations to be tested, is therefore crucial. Here, we review the current state-of-the-art in the field of compound combination modelling, with the aim to support the development of approaches that, as we hope, will finally lead to an integration of chemical with systems-level biological information for predicting the effect of chemical mixtures.

Structure of shikimate kinase, an in vivo essential metabolic enzyme in the nosocomial pathogen Acinetobacter baumannii, in complex with shikimate

Acinetobacter baumannii is an opportunistic Gram-negative pathogen that is an important cause of healthcare-associated infections exhibiting high mortality rates. Clinical isolates of multidrug-resistant (MDR) and extremely drug-resistant (XDR) A. baumannii strains are increasingly being observed. Compounding this concern is the dearth of new antibacterial agents in late-stage development that are effective against MDR and XDR A. baumannii. As part of an effort to address these concerns, two genes (aroA and aroC) of the shikimate pathway have previously been determined to be essential for the growth and survival of A. baumannii during host infection (i.e. to be essential in vivo). This study expands upon these results by demonstrating that the A. baumannii aroK gene, encoding shikimate kinase (SK), is also essential in vivo in a rat soft-tissue infection model. The crystal structure of A. baumannii SK in complex with the substrate shikimate and a sulfate ion that mimics the binding interactions expected for the β-phosphate of ATP was then determined to 1.91 Å resolution and the enzyme kinetics were characterized. The flexible shikimate-binding domain and LID region are compared with the analogous regions in other SK crystal structures. The impact of structural differences and sequence divergence between SKs from pathogenic bacteria that may influence antibiotic-development efforts is discussed.

Computational and Biochemical Discovery of RSK2 as a Novel Target for Epigallocatechin Gallate (EGCG)

The most active anticancer component in green tea is epigallocatechin-3-gallate (EGCG). Protein interaction with EGCG is a critical step for mediating the effects of EGCG on the regulation of various key molecules involved in signal transduction. By using computational docking screening methods for protein identification, we identified a serine/threonine kinase, 90-kDa ribosomal S6 kinase (RSK2), as a novel molecular target of EGCG. RSK2 includes two kinase catalytic domains in the N-terminal (NTD) and the C-terminal (CTD) and RSK2 full activation requires phosphorylation of both terminals. The computer prediction was confirmed by an in vitro kinase assay in which EGCG inhibited RSK2 activity in a dose-dependent manner. Pull-down assay results showed that EGCG could bind with RSK2 at both kinase catalytic domains in vitro and ex vivo. Furthermore, results of an ATP competition assay and a computer-docking model showed that EGCG binds with RSK2 in an ATP-dependent manner. In RSK2^+/+ and RSK2^-/- murine embryonic fibroblasts, EGCG decreased viability only in the presence of RSK2. EGCG also suppressed epidermal growth factor-induced neoplastic cell transformation by inhibiting phosphorylation of histone H3 at Ser10. Overall, these results indicate that RSK2 is a novel molecular target of EGCG.

Anchor-based classification and type-C inhibitors for tyrosine kinases

Tyrosine kinases regulate various biological processes and are drug targets for cancers. At present, the design of selective and anti-resistant inhibitors of kinases is an emergent task. Here, we inferred specific site-moiety maps containing two specific anchors to uncover a new binding pocket in the C-terminal hinge region by docking 4,680 kinase inhibitors into 51 protein kinases, and this finding provides an opportunity for the development of kinase inhibitors with high selectivity and anti-drug resistance. We present an anchor-based classification for tyrosine kinases and discover two type-C inhibitors, namely rosmarinic acid (RA) and EGCG, which occupy two and one specific anchors, respectively, by screening 118,759 natural compounds. Our profiling reveals that RA and EGCG selectively inhibit 3% (EGFR and SYK) and 14% of 64 kinases, respectively. According to the guide of our anchor model, we synthesized three RA derivatives with better potency. These type-C inhibitors are able to maintain activities for drug-resistant EGFR and decrease the invasion ability of breast cancer cells. Our results show that the type-C inhibitors occupying a new pocket are promising for cancer treatments due to their kinase selectivity and anti-drug resistance.

Selective Mycobacterium tuberculosis Shikimate Kinase Inhibitors as Potential Antibacterials

Owing to the persistence of tuberculosis (TB) as well as the emergence of multidrug-resistant and extensively drug-resistant (XDR) forms of the disease, the development of new antitubercular drugs is crucial. Developing inhibitors of shikimate kinase (SK) in the shikimate pathway will provide a selective target for antitubercular agents. Many studies have used in silico technology to identify compounds that are anticipated to interact with and inhibit SK. To a much more limited extent, SK inhibition has been evaluated by in vitro methods with purified enzyme. Currently, there are no data on in vivo activity of Mycobacterium tuberculosis shikimate kinase (MtSK) inhibitors available in the literature. In this review, we present a summary of the progress of SK inhibitor discovery and evaluation with particular attention toward development of new antitubercular agents.

The shikimate dehydrogenase family: functional diversity within a conserved structural and mechanistic framework

Shikimate dehydrogenase (SDH) catalyzes the NADPH-dependent reduction of 3-deydroshikimate to shikimate, an essential reaction in the biosynthesis of the aromatic amino acids and a large number of other secondary metabolites in plants and microbes. The indispensible nature of this enzyme makes it a potential target for herbicides and antimicrobials. SDH is the archetypal member of a large protein family, which contains at least four additional functional classes with diverse metabolic roles. The different members of the SDH family share a highly similar three-dimensional structure and utilize a conserved catalytic mechanism, but exhibit distinct substrate preferences, making the family a particularly interesting system for studying modes of substrate recognition used by enzymes. Here, we review our current understanding of the biochemical and structural properties of each of the five previously identified SDH family functional classes.

Identification of Novel Polyphenolic Inhibitors of Shikimate Dehydrogenase (AroE)

Shikimate dehydrogenase (AroE) is an attractive target for herbicides and antimicrobial agents due to its conserved and essential nature in plants, fungi, and bacteria. Here, we have performed an in vitro screen using a collection of more than 5500 compounds and identified 24 novel inhibitors of AroE from Pseudomonas putida. The IC50 values for the two most potent inhibitors we identified, epigallocatechin gallate (EGCG) and epicatechin gallate (ECG), were 3.0 ± 0.2 µM and 3.7 ± 0.5 µM, respectively. Based on the high level of structural conservation between AroE orthologs, we predicted that the identified compounds would also inhibit AroE enzymes from other organisms. Consistent with this hypothesis, we found that EGCG and ECG inhibit the AroE domain of the bifunctional dehydroquinate dehydratase-shikimate dehydrogenase (DHQ-SDH) from Arabidopsis thaliana with IC50 values of 2.1 ± 0.3 µM and 2.0 ± 0.2 µM, respectively.

Synergistic Synthetic Biology: Units in Concert

Synthetic biology aims at translating the methods and strategies from engineering into biology in order to streamline the design and construction of biological devices through standardized parts. Modular synthetic biology devices are designed by means of an adequate elimination of cross-talk that makes circuits orthogonal and specific. To that end, synthetic constructs need to be adequately optimized through in silico modeling by choosing the right complement of genetic parts and by experimental tuning through directed evolution and craftsmanship. In this review, we consider an additional and complementary tool available to the synthetic biologist for innovative design and successful construction of desired circuit functionalities: biological synergies. Synergy is a prevalent emergent property in biological systems that arises from the concerted action of multiple factors producing an amplification or cancelation effect compared with individual actions alone. Synergies appear in domains as diverse as those involved in chemical and protein activity, polypharmacology, and metabolic pathway complementarity. In conventional synthetic biology designs, synergistic cross-talk between parts and modules is generally attenuated in order to verify their orthogonality. Synergistic interactions, however, can induce emergent behavior that might prove useful for synthetic biology applications, like in functional circuit design, multi-drug treatment, or in sensing and delivery devices. Synergistic design principles are therefore complementary to those coming from orthogonal design and may provide added value to synthetic biology applications. The appropriate modeling, characterization, and design of synergies between biological parts and units will allow the discovery of yet unforeseeable, novel synthetic biology applications.