Protein-protein docking on molecular models of Aspergillus niger RNase and human actin: novel target for anticancer therapeutics

The 3D models of human actin protein and A.niger RNase were designed using the templates ACTBIND (PDB ID: 3D3Z) and crystalline profilin-beta-actin (PDB ID: 2BTF), respectively in Modeller9v5. These models are testified using several validation methods including PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The stereo-chemical quality of the models was judged by Ramachandran plot with PROCHECK. The total quality G-factor −0.2, shows a good quality model. The ERRAT score for the human actin and A.niger RNase models are 86.104 and 84.615, respectively, fit well within the range of a high quality model. The ERRAT score for the templates 2BTF and 3D3Z are 91.111 and 97.391, respectively. The WHAT-IF evaluation justifies a reasonable homology model structure as none of the scores for each residue in the homology model is lower than −5.0. The energy-minimized model of human actin with PROSA reveals the Z-score value −10.52 between native conformations of the crystal structures. The VERIFY 3D average score is 0.36. All evidence suggests that the geometric quality of the backbone conformation, the residue interaction, the residue contact and the energy profile of the structures were well within the limits of reliable structures. The interaction energy of docking was calculated using the HEX server. The Etotal, lowest docked energy, and calculated RMSD values were −1.608 kcal mol⁻¹, -8.369 kcal mol⁻¹ and 0.617 Å, respectively. The study presented in the current project may be useful to design molecules that may have anticancer activity.

     

Homology modelling,virtual screening and molecular dynamics of the MARK3 KA1 domain for cancer drug design

Structural homology modelling was done with the software AMPS, MODELLER, PROCHECK, WHATIF AND VERIFY-3D to generate a quality model of human MARK3. Macromolecular docking simulations seem to confirm recent data in the literature and in MARK3 there does not occur intramolecular interactions between the associated kinase domain KA1 and the catalytic domain. Using virtual screening, we were able to identify and suggest the principal residues of MARK3 which interact with the ligands in addition to those reported in the literature. The pharmacophoric model obtained from Discovery Studio coincides with those obtained by molecular interaction fields, indicating the principal ligand residues of the MARK3 KA1 domain. Using virtual screening with pharmacophoric constraints as well as molecular dynamics, the most stable compounds in the ligand site as well as their potential toxicities were used to select potential inhibitors for further in vitro and in vivo investigations of human MARK3 KA1 domain, which could eventually pass to the market to be used for the treatment of head and neck cancer.     

In silico analysis of Myoglobin in Channa striata

Myoglobin is a cytoplasmic hemoprotein, expressed solely in cardiac myocytes and oxidative skeletal muscle fibers, that reversibly binds O₂ by its heme residue. Myoglobin is an essential oxygen-storage hemoprotein capable of facilitating oxygen transport and modulating nitric oxide homeostasis within cardiac and skeletal myocytes. Functionally, myoglobin is well accepted as an O₂- storage protein in muscle, capable of releasing O₂ during periods of hypoxia or anoxia. There is no evidence available regarding active sites, ligand binding sites, antigenic determinants and the ASA value of myoglobin in Channa striata. We further document the predicted active sites in the structural model with solvent exposed ASA residues. During this study, the model was built by CPH program and validated through PROCHECK, Verify 3D, ERRAT and ProSA for reliability. The active sites were predicted in the model with further ASA analysis of active site residues. The discussed information thus provides the predicted active sites, ligand binding sites, antigenic determinants and ASA values of myoglobin model in Channa striata.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.     

Computational Analysis of N-acetyl transferase in Tribolium castaneum

N-acetyl transferase (NAT) is responsible to catalyze the transfer of acetyl groups to arylamines from acetyl-CoA. Aralkylamine Nacetyl transferase (AANAT), which belongs to GCN5-related N-acetyl transferase member, is a globular 23-kDa cytosolic protein that forms a reversible regulatory complex with 14-3-3 proteins, AANAT regulates the daily cycle of melatonin biosynthesis in mammals, making it an attractive target for therapeutic control of abnormal melatonin production in mood and sleep disorders. There is no evidence available regarding α and β subunits, active site and their ASA value in Dopamine N-acetyl transferase. Therefore, we describe the development of Dopamine N-acetyl transferase model in Tribolium castaneum. We further document the predicted active sites in the structural model with solvent exposed ASA residues. During this study, the model was built by CPH program and validated through PROCHECK, Verify 3D, ERRAT and ProSA for reliability. The active sites were predicted in the model with further ASA analysis of active site residues. The discussed information thus provides insight to the predicted active site and ASA values of Dopamine N-acetyl transferase model in Tribolium castaneum.     

Molecular docking studies of chloroquine and its derivatives against P23pro-zbd domain of chikungunya virus: Implication in designing of novel therapeutic strategies

The arthropod-transmitted chikungunya virus has emerged as an epidemic menace that causes debilitating polyarthritis. With this life-threatening impact on humans, the possible treatment requires to cure the viral infectivity. But, devoid of any vaccine against the chikungunya virus (CHIKV), there is a need to develop a novel chemotherapeutic strategy to treat this noxious infection. CHIKV carries highly compact P23^pro-zbd structure that possesses potential RNA-binding surface domains which extremely influences the use of RNA template during genome replication at the time of infection and pathogenesis. Therefore, computational approaches were used to explore the novel small molecule inhibitors targeting P23^pro-zbd domain. The tertiary structure was modeled and optimized using in silico approaches. The results obtained from PROCHECK (93.1% residues in favored regions), ERRAT (87.480 overall model quality) and ProSA (Z-score: −11.72) revealed the reliability of the proposed model. Interestingly, a previously reported inhibitor, chloroquine possesses good binding affinities with the target domain. In-depth analysis revealed that chloroquine derivatives such as didesethyl chloroquine hydroxyacetamide, cletoquine, hydroxychloroquine exhibited a better binding affinity. Notably, MD simulation analysis exhibited that Thr1312, Ala1355, Ala1356, Asn1357, Asp1364, Val1366, Cys1367, Ala1401, Gly1403, Ser1443, Tyr1444, Gly1445, Asn1459, and Thr1463 residues are the key amino acid responsible for stable ligand-protein interaction. The results obtained from this study provide new insights and advances the understanding to develop a new approach to consider effective and novel drug against chikungunya. However, a detailed in vivo study is required to explore its drug likeliness against this life-threatening disease.     

Exploring the molecular basis for selective binding of Mycobacterium tuberculosis Asp kinase toward its natural substrates and feedback inhibitors: A docking and molecular dynamics study

Tuberculosis (TB) is still a major public health problem, compounded by the human immunodeficiency virus (HIV)-TB co-infection and recent emergence of multidrug-resistant (MDR) and extensively drug resistant (XDR)-TB. In this context, aspartokinase of mycobacterium tuberculosis has drawn attention for designing novel anti-TB drugs. Asp kinase is an enzyme responsible for the synthesis of 4-phospho-L-aspartate from L-aspartate and involved in the branched biosynthetic pathway leading to the synthesis of amino acids lysine, threonine, methionine and isoleucine. An intermediate of lysine biosynthetic branch, mesodiaminopimelate is also a component of the peptidoglycan which is a component of bacterial cell wall. To interfere with the production of all these amino acids and cell wall, it is possible to inhibit Asp kinase activity. This can be achieved using Asp kinase inhibitors. In order to design novel Asp kinase inhibitors as effective anti-TB drugs, it is necessary to have an understanding of the binding sites of Asp kinase. As no crystal structure of the enzyme has yet been published, we built a homology model of Asp kinase using the crystallized Asp kinase from M. Jannaschii, as template structures (2HMF and 3C1M). After the molecular dynamics refinement, the optimized homology model was assessed as a reliable structure by PROCHECK, ERRAT, WHAT-IF, PROSA2003 and VERIFY-3D. The results of molecular docking studies with natural substrates, products and feedback inhibitors are in agreement with the published data and showed that ACT domain plays an important role in binding to ligands. Based on the docking conformations, pharmacophore model can be developed by probing the common features of ligands. By analyzing the results, ACT domain architecture, certain key residues that are responsible for binding to feedback inhibitors and natural substrates were identified. This would be very helpful in understanding the blockade mechanism of Asp kinase and providing insights into rational design of novel Asp kinase inhibitors for M.tuberculosis.     

Phylogenetic analysis, homology modelling, molecular dynamics and docking studies of caffeoyl–CoA-O- methyl transferase (CCoAOMT 1 and 2) isoforms isolated from subabul (Leucaena leucocephala)

Caffeoyl coenzyme A O-methyltransferase (CCoAOMT) is an important enzyme that participates in lignin biosynthesis especially in the formation of cell wall ferulic esters of plants. It plays a pivotal role in the methylation of the 3-hydroxyl group of caffeoyl CoA. Two cDNA clones that code CCoAOMT were isolated earlier from subabul and in the present study; 3D models of CCoAOMT1 and CCoAOMT2 enzymes were built using the MODELLER7v7 software to find out the substrate binding sites. These two proteins differed only in two amino acids and may have little or no functional redundancy. Refined models of the proteins were obtained after energy minimization and molecular dynamics in a solvated water layer. The models were further assessed by PROCHECK, WHATCHECK, Verify_3D and ERRAT programs and the results indicated that these models are reliable for further active site and docking analysis. The refined models showed that the two proteins have 9 and 10 α-helices, 6 and 7 β-sheets respectively. The models were used for docking the substrates CoA, SAM, SAH, caffeoyl CoA, feruloyl CoA, 5-hydroxy feruloyl CoA and sinapyl CoA which showed that CoA and caffeoyl CoA are binding with high affinity with the enzymes in the presence and absence of SAM. It appears therefore that caffeoyl CoA is the substrate for both the isoenzymes. The results also indicated that CoA and caffeoyl CoA are binding with higher affinity to CCoAOMT2 than CCoAOMT1. Therefore, CCoAOMT2 conformation is thought to be the active form that exists in subabul. Docking studies indicated that conserved active site residues Met58, Thr60, Val63, Glu82, Gly84, Ser90, Asp160, Asp162, Thr169, Asn191 and Arg203 in CCoAOMT1 and CCoAOMT2 enzymes create the positive charge to balance the negatively charged caffeoyl CoA and play an important role in maintaining a functional conformation and are directly involved in donor-substrate binding.     

Insights from Streptococcus pneumoniae glucose kinase structural model

Streptococcus pneumonia is the common cause of sepsis and meningitis. Emergence of multiple antibiotic resistant strains in the community‐acquired bacterium is catastrophic. Glucose kinase (GLK) is a regulatory enzyme capable of adding phosphate group to glucose in the first step of streptomycin biosynthesis. The activity of glucose kinase was regulated by the Carbon Catabolite Repression (CCR) system. Therefore, it is important to establish the structure‐function relation of GLK in S. pneumoniae. However, a solved structure for S. pneumoniae GLK is not available at the protein data bank (PDB). Therefore, we created a model of GLK from S. pnemoniae using the X‐ray structure of Glk from E. faecalis as template with MODELLER (a comparative modeling program). The model was validated using protein structure checking tools such as PROCHECK, WHAT IF and ProSA for reliability. The active site amino acid Asp114 in the template is retained in S. pneumoniae GLK model (Asp115). Solvent accessible surface area (ASA) analysis of the GLK model showed that known key residues playing important role in active site for ligand binding and metal ion binding are buried and hence not accessible to solvent. The information thus discussed provides insight to the molecular understanding of glucose kinase in S. pneumoniae.     

Structure based drug designing of a novel antiflaviviral inhibitor for nonstructural 3 protein

Literature shows that Flaviviruses cause a variety of diseases, including fevers, encephalitis, and hemorrhagic fevers. NS3 is a multifunctional protein with an Nterminal protease domain (NS3pro) that is responsible for proteolytic processing of the viral polyprotein, and a C-terminal region that contains an RNA triphosphatase, RNA helicase and RNA-stimulated NTPase domain that are essential for RNA replication. Therefore, NS3 protein is the preferential choice for inhibition to stop the proteolytic processing. Hence, the 3D structure of NS3 protein was modeled using homology modeling by MODELLER 9v7. Evaluation of the constructed NS3 protein models were done by PROCHECK, VERYFY3D and through ProSA calculations. Ligands for the catalytic triad were designed using LIGBUILDER. The NS3 protein's catalytic triad was explored to find out the critical interactions pattern for inhibitor binding using molecular docking methodology using AUTODOCK Vina. It should be noted that these predicted data should be validated using suitable assays for further consideration. ABBREVIATIONS: DOPE - Discrete optimized protein energy, WHO - World Health Organization, ADME/T - Absorption, Distribution, Metabolism, Excretion and Toxicity.     

Structure modeling and antidiabetic activity of a seed protein of Momordica charantia in non-obese diabetic (NOD) mice

Momordica charantia is a well known medicinal plant used in the traditional medicinal system for the treatment of various diseases including diabetes mellitus. Recently, a novel protein termed as ADMc1 from the seed extract of M. charantia has been identified and isolated showing significant antihyperglycemic activity in type 1 diabetic rats in which diabetes was induced. However, the structure of this protein has not yet been analyzed. Homology modeling approach was used to generate a high quality protein 3D structure for the amino acid sequence of the ADMc1 protein in this study. The comparative assessment of secondary structures revealed ADMc1 as an all-alpha helix protein with random coils. Tertiary structure predicted on the template structure of Napin of B. Napus (PDB ID: 1SM7) with which the ADMc1 showed significant sequence similarity, was validated using protein structure validation tools like PROCHECK, WHAT_CHECK, VERIFY3D and ProSA. Arrangement of disulfide bridges formed by cysteine residues were predicted by the Dianna 1.1 server. The presence of multiple disulfide bond confers the stable nature of the ADMc1 protein. Further, the biological activity of the ADMc1 was assessed in non-obese diabetic (NOD) mice which are spontaneous model of type 1 diabetes. Significant reduction in the blood glucose levels of NOD mice was observed up to 8 h post administration of the rADMc1 protein. Overall, the structural characterizations with antihyperglycemic activity of this seed protein of Momordica charantia demonstrate its potential as an antidiabetic agent.     

Homology modelling of a sensor histidine kinase from Aeromonas hydrophila

Aeromonas hydrophila has been implicated in extra-intestinal infection and diarrhoea in humans. Targetting unique effectors of bacterial pathogens is considered a powerful strategy for drug design against bacterial variations to drug resistance. The two-component bacterial system involving sensor histidine kinase (SHK) and its response regulators is considered a lucrative target for drug design. This is the first report describing a three-dimensional (3D) structure for SHK of A. hydrophila. The model was constructed by homology modelling using the X-ray structure of PleD—a response regulator—in conjunction with cdiGMP (PDB code 1W25) and HemAT sensor domain (PDB code 1OR4)—a globin coupled sensor. A combination of homology modelling methodology and molecular dynamics (MD) simulations was applied to obtain a reasonable structure to understand the dynamic behaviour of SHK. Homology modelling was performed using MODELLER9v2 software. The structure was relaxed to eliminate bad atomic contacts. The final model obtained by molecular mechanics and dynamics methods was assessed using PROCHECK and VERIFY 3D graph, which confirmed that the final refined model is reliable. Until complete biochemical and structural data of SHK are determined by experimental means, this model can serve as a valuable reference for characterising the protein and could be explored for drug targetting by design of suitable inhibitors.     

Structure of glutathione S-transferase of the filarial parasite Wuchereria bancrofti: a target for drug development against adult worm

A three dimensional structural model of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. The three dimensional X-ray crystal structure of porcine -class GST with PDB ID: 2gsr-A chain protein with 42% sequential and functional homology was used as the template. The model of wbGST built by MODELLER6v2 was analyzed by the PROCHECK programs. Ramachandran plot analysis showed that 93.5% of the residues are in the core region followed by 5.4 and 1.1% residues in the allowed and generously allowed regions, respectively. None of the non-glycine residues is in disallowed regions. The PROSA II z-score and the energy graph for the final model further confirmed the quality of the modeled structure. The computationally modeled three-dimensional (3D) structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668). 1SFM was used for docking with GST inhibitors by Hex4.2 macromolecular docking using spherical polar Fourier correlations.Figure: A three-dimensional (3D) structure of Glutathione-S-transferase (GST) of the lymphatic filarial parasite Wuchereria bancrofti (wb) was constructed by homology modeling. This modeled 3D structure of wbGST has been submitted to the Protein Data Bank (PDB) (PDB ID: 1SFM and RCSB ID: RCSB021668).     

Modeling of three dimensional structure of human alpha-fetoprotein complexed with diethylstilbestrol: docking and molecular dynamics simulation study

It has been long experimentally demonstrated that human alpha-fetoprotein (HAFP) has an ability to bind immobilized estrogens with the most efficiency for synthetic estrogen analog - diethylstilbestrol (DES). However, the question remains why the human AFP (HAFP), unlike rodent AFP, cannot bind free estrogens. Moreover, despite the fact that AFP was first discovered more than 50 years ago and is presently recognized as a "golden standard" among onco-biomarkers, its three-dimensional (3D) structure has not been experimentally solved yet. In this work using MODELLER program, we generated 3D model of HAFP on the basis of homology with human serum albumin (HSA) and Vitamin D-binding protein (VTDB) with subsequent molecular docking of DES to the model structure and molecular dynamics (MD) simulation study of the complex obtained. The model constructed has U-shaped structure in which a cavity may be distinguished. In this cavity the putative estrogen-binding site is localized. Validation by RMSD calculation and with the use of PROCHECK program showed good quality of the model and stability of extended region of four alpha-helical structures that contains putative hormone-binding residues. Data extracted from MD simulation trajectory allow proposing two types of interactions between amino acid residues of HAFP and DES molecule: (1) hydrogen bonding with involvement of residues S445, R452, and E551; (2) hydrophobic interactions with participation of L138, M448, and M548 residues. A suggestion is made that immobilization of the hormone using a long spacer provides delivery of the estrogen molecule to the binding site and, thereby, facilitates interaction between HAFP and the hormone.     

Streptomyces-derived actinomycin D inhibits biofilm formation by Staphylococcus aureus and its hemolytic activity

Staphylococcus aureus is a versatile human pathogen that produces diverse virulence factors, and its biofilm cells are difficult to eradicate due to their inherent ability to tolerate antibiotics. The anti-biofilm activities of the spent media of 252 diverse endophytic microorganisms were investigated using three S. aureus strains. An attempt was made to identify anti-biofilm compounds in active spent media and to assess their anti-hemolytic activities and hydrophobicities in order to investigate action mechanisms. Unlike other antibiotics, actinomycin D (0.5 μg ml^?1) from Streptomyces parvulus significantly inhibited biofilm formation by all three S. aureus strains. Actinomycin D inhibited slime production in S. aureus and it inhibited hemolysis by S. aureus and caused S. aureus cells to become less hydrophobic, thus supporting its anti-biofilm effect. In addition, surface coatings containing actinomycin D prevented S. aureus biofilm formation on glass surfaces. Given these results, FDA-approved actinomycin D warrants further attention as a potential antivirulence agent against S. aureus infections.     

Sequence analysis and homology modeling of peroxidase from Medicago sativa

Plant peroxidases are one of the most extensively studied group of enzymes which find applications in the environment, health, pharmaceutical, chemical and biotechnological processes. Class III secretary peroxidase from alfalfa (Medicago sativa) has been characterized using bioinformatics approach Physiochemical properties and topology of alfalfa peroxidase were compared with that of soybean and horseradish peroxidase, two most popular commercially available peroxidase preparations. Lower value of instability index as predicted by ProtParam and presence of extra disulphide linkages as predicted by Cys_REC suggested alfalfa peroxidase to be more stable than either of the commercial preparations. Multiple Sequence Alignment (MSA) with other functionally similar proteins revealed the presence of highly conserved catalytic residues. Three dimensional model of alfalfa peroxidase was constructed based on the crystal structure of soybean peroxidase (PDB Id: 1FHF A) by homology modelling approach. The model was checked for stereo chemical quality by PROCHECH, VERIFY 3D, WHAT IF, ERRAT, 3D MATCH AND ProSA servers. The best model was selected, energy minimized and used to analyze structure function relationship with substrate hydrogen peroxide by Autodock 4.0. The enzyme substrate complex was viewed with Swiss PDB viewer and one residue ASP43 was found to stabilize the interaction by hydrogen bonds. The results of the study may be a guiding point for further investigations on alfalfa peroxidase.     

A homology model for <Emphasis Type="Italic">Clostridium difficile</Emphasis> methionyl tRNA synthetase: active site analysis and docking interactions

Treatment of C. difficile infection is one of the most difficult biomedical challenges. To develop novel antibacterials, researchers have been targeting bacterial molecular functions that are essential for its growth. The methionyl tRNA synthetase (MetRS) is strictly required for protein biosynthesis and success was reported in developing antibacterials to inhibit this enzyme. The present study was aimed at building and analyzing a homology model for C. difficile MetRS in the context of drug design. A homology model of C. difficile MetRS was constructed using Molecular Operating Environment (MOE) software. A. aeolicus MetRS was the main template while the query zinc binding domain was modeled using T. thermophilus MetRS. The model has been assessed and compared to its main template (Ramachandran, ERRAT and ProSA). The active site of the query protein has been predicted from its sequence using a detailed conservation analysis (ClustalW2). Using MOE software, suitable ligands were docked in the constructed model, including a C. difficile MetRS inhibitor REP3123 and the enzyme natural substrate, and the key active site residues and interactions were identified. These docking studies have validated the active site conformation in the constructed model and identified binding interactions.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Non-toxic plant metabolites regulate Staphylococcus viability and biofilm formation: a natural therapeutic strategy useful in the treatment and prevention of skin infections

In the present study, the efficacy of generally recognised as safe (GRAS) antimicrobial plant metabolites in regulating the growth of Staphylococcus aureus and S. epidermidis was investigated. Thymol, carvacrol and eugenol showed the strongest antibacterial action against these microorganisms, at a subinhibitory concentration (SIC) of ≤ 50 μg ml^?1. Genistein, hydroquinone and resveratrol showed antimicrobial effects but with a wide concentration range (SIC = 50–1,000 μg ml^?1), while catechin, gallic acid, protocatechuic acid, p-hydroxybenzoic acid and cranberry extract were the most biologically compatible molecules (SIC ≥ 1000 μg ml^?1). Genistein, protocatechuic acid, cranberry extract, p-hydroxybenzoic acid and resveratrol also showed anti-biofilm activity against S. aureus, but not against S. epidermidis in which, surprisingly, these metabolites stimulated biofilm formation (between 35% and 1,200%). Binary combinations of cranberry extract and resveratrol with genistein, protocatechuic or p-hydroxibenzoic acid enhanced the stimulatory effect on S. epidermidis biofilm formation and maintained or even increased S. aureus anti-biofilm activity.