ABSTRACT
Diabetic neuropathy (DN) is a long-term complication of diabetes mellitus, affecting different periphery nerve systems including sensory and motor neurons. Hyperglycemia is the major cause of DN with symptoms such as weakness of balance or coordination, insensitivity to sensation, weakness of the muscles as well as numbness and pain in limbs Analgesic drug such as opioids can be effective to relief neuropathy pain but there is no effective treatment. Adiponectin is an anti-diabetic adipokine, which possesses insulin-sensitizing and neuroprotective effects. In this project, we aim to identify an agent which is dual acting to opioid and adiponectin receptors. Within a virtual screening repositioning campaign, a large collection of compounds with different structures comprehensive of adipoRon-like piperidine derivatives was screened by docking. Recently developed opioid receptor benzomorphanic agonists finally emerged as good ligands to adiponectin receptors showing some 2D and 3D structural similarities with AdipoRon. Particularly, we have identified (+)-MML1017, which has high affinity to the same binding domain of AdipoR1 and AdipoR2 as AdipoRon. Our western blot results indicate (+)-MML1017 activates AMPK phosphorylation through both adipoR1 and adipoR2 in neuronal cell line. Moreover, pretreatment of (+)-MML1017 can improve the cell viability with motor neurons under hyperglycermic conditions. The (+)-MML1017 also activates µ-opioid receptor cells in a concentration-dependent manner. Our study identified a novel compound having dual activity on opioid receptors and adiponectin receptors that may have analgesic effects and neuroprotective effects to treat diabetic neuropathy.
Subject(s)
Diabetes Mellitus , Diabetic Neuropathies , Neuroprotective Agents , Humans , Receptors, Adiponectin/metabolism , Analgesics, Opioid , Diabetic Neuropathies/drug therapy , Receptors, Opioid , Adiponectin/metabolismABSTRACT
Endogenous peptides as part of physiological processes are targets of interest when it comes to finding desirable therapeutics which are able to modulate molecular interactions. The major limits presented by peptides when they are used as drugs have motivated the research of the synthesis of peptidomimetics obtained through chemical modification and the use of in silico approaches. Here recent works on the discovery of peptidomimetics by computational methods are reported. Together with molecular dynamic simulations, the use of pharmacophore research simulations helps to gain insight into and understand the molecular determinants underlying the physiological processes.
Subject(s)
Molecular Dynamics Simulation , Peptides/chemistry , Peptidomimetics , Peptides/chemical synthesis , SoftwareABSTRACT
The irreversible inhibitors of monoamine oxidases (MAO) slow neurotransmitter metabolism in depression and neurodegenerative diseases. After oxidation by MAO, hydrazines, cyclopropylamines and propargylamines form a covalent adduct with the flavin cofactor. To assist the design of new compounds to combat neurodegeneration, we have updated the kinetic parameters defining the interaction of these established drugs with human MAO-A and MAO-B and analyzed the required features. The Ki values for binding to MAO-A and molecular models show that selectivity is determined by the initial reversible binding. Common to all the irreversible inhibitor classes, the non-covalent 3D-chemical interactions depend on a H-bond donor and hydrophobic-aromatic features within 5.7 angstroms apart and an ionizable amine. Increasing hydrophobic interactions with the aromatic cage through aryl halogenation is important for stabilizing ligands in the binding site for transformation. Good and poor inactivators were investigated using visible spectroscopy and molecular dynamics. The initial binding, close and correctly oriented to the FAD, is important for the oxidation, specifically at the carbon adjacent to the propargyl group. The molecular dynamics study also provides evidence that retention of the allenyl imine product oriented towards FADH- influences the formation of the covalent adduct essential for effective inactivation of MAO.
Subject(s)
Monoamine Oxidase Inhibitors/chemistry , Monoamine Oxidase Inhibitors/pharmacology , Monoamine Oxidase/chemistry , Binding Sites , Dose-Response Relationship, Drug , Enzyme Activation/drug effects , Humans , Inhibitory Concentration 50 , Kinetics , Molecular Conformation , Molecular Docking Simulation , Molecular Dynamics Simulation , Molecular Structure , Oxidation-Reduction , Protein Binding , Structure-Activity Relationship , Substrate Specificity , Time FactorsABSTRACT
Background: Central and peripheral analgesia without adverse effects relies on the identification of µ-opioid agonists that are able to activate 'basal' antinociceptive pathways. Recently developed µ-selective benzomorphan agonists that are not antagonized by naloxone do not activate G-proteins and ß-arrestins. Which pathways do µ receptors activate? How can each of them be selectively activated? What role is played by allosteric binding sites? Methodology & results: Molecular modeling studies characterize the amino acid residues involved in the interaction with various classes of endogenous and exogenous ligands and with agonists and antagonists. Conclusions: Critical binding differences between various classes of agonists with different pharmacological profiles have been identified. MML series binding poses may be relevant in the search for an antinociception agent without side effects.
Subject(s)
Analgesics, Opioid/pharmacology , Molecular Dynamics Simulation , Receptors, Opioid, mu/agonists , Receptors, Opioid, mu/antagonists & inhibitors , Analgesics, Opioid/chemistry , Binding Sites/drug effects , Humans , Ligands , Molecular StructureABSTRACT
Functional antitumor vaccine constructs are the basis for active tumor immunotherapy, which is useful in the treatment of many types of cancers. MUC1 is one key glycoprotein for targeting and designing new strategies for multicomponent vaccines. Two self-adjuvant tetravalent vaccine candidates were prepared by clustering four or eight PDTRP MUC1 core epitope sequences on calixarene scaffolds. In this work, the different activities of two molecules with calix[4]arene and calix[8]arene skeleton are rationalized. Quantum mechanics, docking, and molecular dynamics structural optimization were first carried out followed by metadynamics to calculate the energy profiles. Further insights were obtained by complementarity studies of molecular fields. The molecular modeling results are in strong agreement with the experimental in vivo immunogenicity data. In conclusion, the overall data shows that, in the designing of anticancer vaccines, scaffold flexibility has a pivotal role in obtaining a suitable electrostatic, hydrophobic, and steric complementarity with the biological target.
Subject(s)
Calixarenes , Neoplasms , Vaccines , Humans , Molecular Dynamics Simulation , Mucin-1 , Static ElectricityABSTRACT
Modeling G-Protein Coupled Receptors (GPCRs) is an emergent field of research, since utility of high-quality models in receptor structure-based strategies might facilitate the discovery of interesting drug candidates. The findings from a quantitative analysis of eighteen resolved structures of rhodopsin family "A" receptors crystallized with antagonists and 153 pairs of structures are described. A strategy termed endeca-amino acids fragmentation was used to analyze the structures models aiming to detect the relationship between sequence identity and Root Mean Square Deviation (RMSD) at each trans-membrane-domain. Moreover, we have applied the leave-one-out strategy to study the shiftiness likelihood of the helices. The type of correlation between sequence identity and RMSD was studied using the aforementioned set receptors as representatives of membrane proteins and 98 serine proteases with 4753 pairs of structures as representatives of globular proteins. Data analysis using fragmentation strategy revealed that there is some extent of correlation between sequence identity and global RMSD of 11AA width windows. However, spatial conservation is not always close to the endoplasmic side as was reported before. A comparative study with globular proteins shows that GPCRs have higher standard deviation and higher slope in the graph with correlation between sequence identity and RMSD. The extracted information disclosed in this paper could be incorporated in the modeling protocols while using technique for model optimization and refinement.
Subject(s)
Models, Molecular , Protein Conformation , Rhodopsin/chemistry , Amino Acid Sequence , Humans , Molecular Dynamics Simulation , Protein Binding , Protein Interaction Domains and Motifs , Receptors, G-Protein-Coupled/chemistry , Structure-Activity RelationshipABSTRACT
G protein-coupled receptors (GPCRs) are a super-family of membrane proteins that attract great pharmaceutical interest due to their involvement in almost every physiological activity, including extracellular stimuli, neurotransmission, and hormone regulation. Currently, structural information on many GPCRs is mainly obtained by the techniques of computer modelling in general and by homology modelling in particular. Based on a quantitative analysis of eighteen antagonist-bound, resolved structures of rhodopsin family "A" receptors - also used as templates to build 153 homology models - it was concluded that a higher sequence identity between two receptors does not guarantee a lower RMSD between their structures, especially when their pair-wise sequence identity (within trans-membrane domain and/or in binding pocket) lies between 25 % and 40 %. This study suggests that we should consider all template receptors having a sequence identity ≤50 % with the query receptor. In fact, most of the GPCRs, compared to the currently available resolved structures of GPCRs, fall within this range and lack a correlation between structure and sequence. When testing suitability for structure-based drug design, it was found that choosing as a template the most similar resolved protein, based on sequence resemblance only, led to unsound results in many cases. Molecular docking analyses were carried out, and enrichment factors as well as attrition rates were utilized as criteria for assessing suitability for structure-based drug design.
Subject(s)
Receptors, G-Protein-Coupled/chemistry , Receptors, G-Protein-Coupled/metabolism , Receptors, Serotonin/chemistry , Receptors, Serotonin/metabolism , Amino Acid Sequence , Animals , Drug Design , Humans , Models, Molecular , Molecular Docking Simulation/methods , Rhodopsin/chemistry , Sequence Homology, Amino AcidABSTRACT
Agricultural practices are usually supported by several chemical substances, such as herbicides. Linuron and chlorbromuron are phenylurea herbicides largely used to protect crops from weeds, blocking photosynthesis by inhibition of the photosystem II complex. The former, also commercially known as lorox or afalon, is selectively used to protect bean and French bean plants, fennels, and celeriacs; the second, commercially known as maloran, is selectively used for carrots, peas, potatoes, soy sprouts, and sunflowers. Considering the widespread use of herbicides and, more generally, pesticides, it is important to clarify their involvement on human health, one of them concerning the possible direct or indirect effect on the genome of exposed populations. Here, we show that these herbicides are endowed by mutagenic properties, as demonstrated by an increased number of chromosomal aberrations (CAs) in two exposed Chinese hamster cell lines derived from ovary and epithelial liver, respectively. This was also confirmed by sister chromatid exchange (SCE) and micronucleus (MN) assays. Our present and previously obtained data clearly indicate that phenylurea herbicides must be used with great caution, especially for agricultural workers who use large amounts of herbicides during their work, and particular attention should be given to residues of these herbicides and their involvement in environmental pollution.
Subject(s)
Biomarkers/analysis , Chromosome Aberrations/drug effects , Herbicides/toxicity , Linuron/toxicity , Methylurea Compounds/toxicity , Mutagens/toxicity , Phenylurea Compounds/toxicity , Animals , CHO Cells , Cell Line , Cricetinae , Cricetulus , Female , Micronucleus Tests , Sister Chromatid Exchange/drug effectsABSTRACT
The human 5-HT7 receptor is expressed in both the central nervous system and peripheral tissues and is a potential drug target in behavioral and psychiatric disorders. We examined molecular determinants of ligand binding and G protein activation by the human 5-HT7(a) receptor. The role of several key residues in the 7th transmembrane domain (TMD) and helix 8 were elucidated combining in silico and experimental mutagenesis. Several single and two double point mutations of the 5-HT7(a) wild type receptor were made (W7.33V, E7.35T, E7.35R, E7.35D, E7.35A, R7.36V, Y7.43A, Y7.43F, Y7.43T, R8.52D, D8.53K; E7.35T-R7.36V, R8.52D-D8.53K), and their effects upon ligand binding were assessed by radioligand binding using a potent agonist (5-CT) and a potent antagonist (SB269970). In addition, the ability of the mutated 5-HT7(a) receptors to activate G protein after 5-HT-stimulation was determined through activation of adenylyl cyclase. In silico investigation on mutated receptors substantiated the predicted importance of TM7 and showed critical roles of residues E7.35, W7.33, R7.36 and Y7.43 in agonist and antagonist binding and conformational changes of receptor structure affecting adenylyl cyclase activation. Experimental data showed that mutants E7.35T and E7.35R were incapable of ligand binding and adenylyl cyclase activation, consistent with a requirement for a negatively charged residue at this position. The mutant R8.52D was unable to activate adenylyl cyclase, despite unaffected ligand binding, consistent with the R8.52 residue playing an important role in the receptor-G protein interface. The mutants Y7.43A and Y7.43T displayed reduced agonist binding and AC agonist potency, not seen in Y7.43F, consistent with a requirement for an aromatic residue at this position. Knowledge of the molecular interactions important in h5-HT7 receptor ligand binding and G protein activation will aid the design of selective h5-HT7 receptor ligands for potential pharmacological use.
ABSTRACT
The presence of polybrominated flame retardants in the environment seems to be increasing in the past decade. Considering the toxic effects of these pollutants, it is important evaluating the potential interaction with biological membranes for a risk assessment. In this study low and high brominated biphenyls and biphenyl ethers were used to investigate their interaction with biological membrane models constituted by liposomes, using differential scanning calorimetry (DSC) technique. The medium influence on membrane absorption was also assessed. The findings indicate that membrane interaction is controlled by compound structural characteristics. The membrane absorption is allowed by lipophilic medium; instead hydrophilic medium prevents membrane permeation.
Subject(s)
Calorimetry, Differential Scanning/methods , Flame Retardants/metabolism , Liposomes/metabolism , Polybrominated Biphenyls/metabolism , Absorption, Physicochemical , Environmental Pollutants/analysis , Hydrophobic and Hydrophilic Interactions , Membrane Lipids/metabolism , Membranes, Artificial , Phenyl Ethers/chemistry , Phenyl Ethers/metabolism , Polybrominated Biphenyls/chemistryABSTRACT
The human histamine H4 receptor (hH4R), a member of the G-protein coupled receptors (GPCR) family, is an increasingly attractive drug target. It plays a key role in many cell pathways and many hH4R ligands are studied for the treatment of several inflammatory, allergic and autoimmune disorders, as well as for analgesic activity. Due to the challenging difficulties in the experimental elucidation of hH4R structure, virtual screening campaigns are normally run on homology based models. However, a wealth of information about the chemical properties of GPCR ligands has also accumulated over the last few years and an appropriate combination of these ligand-based knowledge with structure-based molecular modeling studies emerges as a promising strategy for computer-assisted drug design. Here, two chemoinformatics techniques, the Intelligent Learning Engine (ILE) and Iterative Stochastic Elimination (ISE) approach, were used to index chemicals for their hH4R bioactivity. An application of the prediction model on external test set composed of more than 160 hH4R antagonists picked from the chEMBL database gave enrichment factor of 16.4. A virtual high throughput screening on ZINC database was carried out, picking â¼ 4000 chemicals highly indexed as H4R antagonists' candidates. Next, a series of 3D models of hH4R were generated by molecular modeling and molecular dynamics simulations performed in fully atomistic lipid membranes. The efficacy of the hH4R 3D models in discrimination between actives and non-actives were checked and the 3D model with the best performance was chosen for further docking studies performed on the focused library. The output of these docking studies was a consensus library of 11 highly active scored drug candidates. Our findings suggest that a sequential combination of ligand-based chemoinformatics approaches with structure-based ones has the potential to improve the success rate in discovering new biologically active GPCR drugs and increase the enrichment factors in a synergistic manner.
Subject(s)
Drug Discovery/methods , Informatics/methods , Molecular Docking Simulation , Molecular Dynamics Simulation , Receptors, G-Protein-Coupled/antagonists & inhibitors , Receptors, G-Protein-Coupled/metabolism , Receptors, Histamine/metabolism , Databases, Pharmaceutical , Humans , Inhibitory Concentration 50 , Ligands , Protein Binding , Protein Conformation , Receptors, G-Protein-Coupled/chemistry , Receptors, Histamine/chemistry , Receptors, Histamine H4 , Sequence Homology, Amino Acid , Thermodynamics , Time FactorsABSTRACT
Amylin is a 37-residue peptide hormone produced by the islet ß-cells of pancreas and the formation of amylin aggregates is strongly associated with ß-cell degeneration in type 2 diabetes, as demonstrated by more than 95% of patients exhibiting amylin amyloid upon autopsy. It is widely recognized that metal ions such as copper(II) have been implicated in the aggregation process of amyloidogenic peptides such as Aß and α-synuclein and there is evidence that amylin self-assembly is also largely affected by copper(II). For this reason, in this work, the role of copper(II) in the aggregation of amylin has been investigated by several different experimental approaches. Mass spectrometric investigations show that copper(II) induces significant changes in the amylin structure, which decrease the protein fibrillogenesis as observed by ThT measurements. Accordingly, solid-state NMR experiments together with computational analysis carried out on a model amylin fragment confirmed the non-fibrillogenic nature of the copper(II) induced aggregated structure. Finally, the presence of copper(II) is also shown to have a major influence on amylin proneness to be degraded by proteases and cytotoxicity studies on different cell cultures are reported.
Subject(s)
Copper/metabolism , Islet Amyloid Polypeptide/chemistry , Islet Amyloid Polypeptide/metabolism , Protein Aggregates , Amino Acid Sequence , Cell Line, Tumor , Cell Survival , Humans , Models, Molecular , Molecular Sequence Data , Protein Structure, Secondary , ProteolysisABSTRACT
Polycyclic aromatic hydrocarbons are a family of ubiquitous pollutants whose environmental behavior has been widely studied. Different bacterial species are able to decompose hydrocarbons by using them as a food source. One of the best-studied enzymes is naphthalene 1,2-dioxygenase (NDO). A practical way to optimize the degradation process is by mutating the protein involved, increasing both the degradation capacity of the enzyme and its ability to work under extreme environmental conditions of high temperature and low pH. Herein, we describe the study of NDO using molecular dynamics and docking calculations to discover new mutants with high degrading capabilities. We modeled eleven new mutants of NDO. The results indicate that increasing the size of the active site cavity in the mutants allowed for the insertion of high molecular weight PAHs. Additionally, the physicochemical properties of the NDO active sites make the sites well suited to interactions with PAHs, so most amino-acid modifications should not result in significantly altered behavior of NDO.
Subject(s)
Dioxygenases/chemistry , Dioxygenases/genetics , Molecular Docking Simulation , Molecular Dynamics Simulation , Binding Sites , Catalytic Domain , Mutant Proteins/chemistry , Mutation , Polycyclic Aromatic Hydrocarbons/chemistry , Protein Structure, Secondary , ThermodynamicsABSTRACT
Imidazoline ligands in I2-type binding sites in the brain alter monoamine turnover and release. One example of an I2 binding site characterized by binding studies, kinetics, and crystal structure has been described in monoamine oxidase B (MAO B). MAO A also binds imidazolines but has a different active site structure. Docking and molecular dynamics were used to explore how 2-(2-benzofuranyl)-2-imidazoline hydrochloride (2-BFI) binds to MAO A and to explain why tranylcypromine increases tight binding to MAO B. The energy for 2-BFI binding to MAO A was comparable to that for tranylcypromine-modified MAO B, but the location of 2-BFI in the MAO A could be anywhere in the monopartite substrate cavity. Binding to the tranylcypromine-modified MAO B was with high affinity and in the entrance cavity as in the crystal structure, but the energies of interaction with the native MAO B were less favorable. Molecular dynamics revealed that the entrance cavity of MAO B after tranylcypromine modification is both smaller and less flexible. This change in the presence of tranylcypromine may be responsible for the greater affinity of tranylcypromine-modified MAO B for imidazoline ligands.
Subject(s)
Imidazolines/metabolism , Monoamine Oxidase/metabolism , Binding Sites , Humans , Imidazolines/chemistry , Molecular Dynamics SimulationABSTRACT
In recent years, the number of studies in the field of bioremediation has been growing steadily. Although a large number of studies provide information that is highly detailed and offer great amounts of knowledge on a given subject, the downside is that the hunt for more information requires the combined efforts of researchers from many areas, which are becoming increasingly difficult to attain. In this review, we present an overview of recent work investigating enzyme degradation of polycyclic aromatic hydrocarbons. In the first part, this review examines several of the new enzymes able to degrade pollutants, with special attention being given to those with a well-resolved structure. The second part explores some of the most recent work in which computational approaches, such as molecular dynamics, docking, density functional theory and database retrieval, have been employed to study enzymes with specific bioremediation activities.
Subject(s)
Computer Simulation , Environmental Restoration and Remediation , Models, Theoretical , Polycyclic Aromatic Hydrocarbons/chemistry , Biocatalysis , Environmental Pollutants/chemistry , Environmental Pollutants/metabolism , Enzymes/chemistry , Enzymes/metabolism , Polycyclic Aromatic Hydrocarbons/metabolismABSTRACT
In this paper, the techniques of modelling, docking and molecular dynamics were used to study eight single amino acid mutations of the enzyme PhnI to optimise its enzymatic degradation capability. The eight mutants were first equilibrated to avoid deformations of the secondary and tertiary structure and to minimise alterations in the functionality of the chimera enzymes that were obtained. For this purpose, we monitored the potential energy of the systems and the fluctuations of the backbone of the enzymes. The structures of mutant enzymes, at equilibrium, were subjected to docking calculations with selected PAHs. The results indicated a significant increase in the PAH-enzyme interaction with respect to the wild-type protein. The considerable computing resources offered by the GRID computing system made it possible to perform calculations on the entire enzyme system, consisting of six protein subunits, as highlighted in the recent literature.
Subject(s)
Enzymes/chemistry , Molecular Dynamics Simulation , Polycyclic Aromatic Hydrocarbons/chemistry , Enzymes/metabolism , Molecular Conformation , Polycyclic Aromatic Hydrocarbons/metabolism , Protein BindingABSTRACT
The hydrosoluble resveratrol derivative 3-O-phosphorylresveratrol was shown to be more cytotoxic against DU 145 prostate cancer cells than its analog 4'-O-phosphorylresveratrol. In an attempt to unveil the molecular determinants that lye at the root of their different biological effects, here we investigate the interactions of the two resveratrol derivatives with DMPC model membranes by using DSC, membrane permeation/poration assays and molecular dynamics. The results show that the 3-O-derivative interacts with DMPC membranes and diffuses across them. The 4'-O-derivative lies preferentially onto the surface of membrane. The MD simulations provide a molecular interpretation of the experiments and highlight that, in order to maximize the apolar interactions, the 3-O-derivative is embedded in the lipid hydrophobic region. This topographical position of the 3-O resveratrol analog perturbs the liquid-crystalline order of the lipid bilayer promoting membrane curvature and partial lipid loss from the vesicle. This finding reconciles with the lowering of the enthalpy of the lipid phase transition and the ability of the molecule to diffuse across membranes. The present data contribute to explain the different biological activity of the two molecules and evidence that membrane permeability is a key requirement for effective design of resveratrol derivatives to be used for therapeutic purposes.
Subject(s)
Stilbenes/chemistry , Stilbenes/pharmacology , Antineoplastic Agents, Phytogenic/chemistry , Antineoplastic Agents, Phytogenic/pharmacology , Calorimetry, Differential Scanning , Cell Line, Tumor , Cell Membrane Permeability , Dimyristoylphosphatidylcholine/chemistry , Drug Design , Humans , Lipid Bilayers/chemistry , Male , Membranes, Artificial , Models, Molecular , Molecular Dynamics Simulation , Phosphorylation , ResveratrolABSTRACT
The massive computational resources available in the framework of a grid paradigm approach represent an emerging tool in the bioinformatics field. In this paper, we used the above approach in the rapid determination of the interactions between the ring-hydroxylating dioxygenase, comprised six enzymatic subunits, and polycyclic aromatic hydrocarbons (PAHs) in their optimal positions. The results were obtained by simulating enzyme dynamics at 300 K through molecular dynamics calculations. For the first time, the equilibrated structure of the dioxygenase revealed a network of channels throughout the enzyme that were sufficiently large to allow a flow of small ions or molecules from the inner core of the complex to its exterior surface. The ring-hydroxylating dioxygenase was able to interact with some of the studied PAHs. Additionally, not only the number of aromatic rings but also the PAH shape were critical in predicting the ability of the dioxygenase to interact with these types of molecules. Docking calculations shed light on a new possible binding site that is far from the enzymatic one, which is potentially interesting in considering the stability of the enzyme itself.
Subject(s)
Bacterial Proteins/chemistry , Dioxygenases/chemistry , Polycyclic Aromatic Hydrocarbons/chemistry , Sphingomonas/enzymology , Bacterial Proteins/metabolism , Catalytic Domain , Dioxygenases/metabolism , Enzyme Stability , Models, Molecular , Molecular Dynamics Simulation , Molecular Structure , Polycyclic Aromatic Hydrocarbons/metabolism , Protein ConformationABSTRACT
The intensive use of herbicides over the last few decades has caused a general increase of environmental pollution. It is thus very important to evaluate the possible genotoxic properties of these chemical compounds as well as identifying their mode of action. Phenylurea herbicides are selective agents widely used for the control of infestant plants. Of these herbicides, which are widely used in agriculture, we analysed four of the less intensively studied molecules. More precisely, we investigated the genotoxic effects of fenuron, chlorotoluron, diuron, and difenoxuron by analyses of chromosomal aberrations (CAs) and sister chromatid exchange (SCE) in exposed mammalian cells. We used the Chinese hamster ovary (CHO) and epithelial liver (CHEL) cell lines, endowed with the absence or the presence, respectively, of an enzymatic system to activate pro-mutagenic compounds. Our results show that all herbicides tested induce, at high concentrations, an increasing number of CAs in non-metabolising CHO cells. Instead, in the exposed CHEL cell line, the four herbicides induced CAs also at the lowest dose-level. In the CHEL cells, a statistically significant increase of SCE was also observed. The phenylurea herbicides showed direct genotoxic activity, but the cytogenetic effects were greatly enhanced after metabolic conversion. These data, together with other information on phenylurea herbicides, are of great interest from the environmental point of view, and for human health. In fact, intensive use of herbicides contaminates soil, surface water, groundwater and agricultural products, and thus should be taken in particular consideration not only for those initiatives to specifically protect exposed workers, but also to safeguard the health of consumers of agricultural products.