Characterizing the potentially neuronal acetylcholinesterase reactivity toward chiral pyraclofos: Enantioselective insights from spectroscopy, in silico docking, molecular dynamics simulation and per-residue energy decomposition studies.

Chiral organophosphorus agents are distributed ubiquitously in the environment, but the neuroactivity of these asymmetric chemicals to humans remains uncertain. This scenario was to explore the stereoselective neurobiological response of human acetylcholinesterase (AChE) to chiral pyraclofos at the enantiomeric scale, and then decipher the microscopic basis of enantioselective neurotoxicity of pyraclofos enantiomers. The results indicated that (R)-/(S)-pyraclofos can form the bioconjugates with AChE with a stoichiometric ratio of 1:1, but the neuronal affinity of (R)-pyraclofos (K = 6.31 × 104 M-1) with AChE was larger than that of (S)-pyraclofos (K = 1.86 × 104 M-1), and significant enantioselectivity was existed in the biochemical reaction. The modes of neurobiological action revealed that pyraclofos enantiomers were situated at the substrate binding domain, and the strength of the overall noncovalent bonds between (S)-pyraclofos and the residues was weaker than that of (R)-pyraclofos, resulting in the high inhibitory effect of (R)-pyraclofos toward the activity of AChE. Dynamic enantioselective biointeractions illustrated that the intervention of inherent conformational flexibility in the AChE-(R)-pyraclofos was greater than that of the AChE-(S)-pyraclofos, which arises from the big spatial displacement and the conformational flip of the binding domain composed of the residues Thr-64～Asn-89, Gly-122～Asp-134, and Thr-436～Tyr-449. Energy decomposition exhibited that the Gibbs free energies of the AChE-(R)-/(S)-pyraclofos were ΔG° = ï¼37.4/－30.2 kJ mol-1, respectively, and the disparity comes from the electrostatic energy during the stereoselective neurochemical reactions. Quantitative conformational analysis further confirmed the atomic-scale computational chemistry conclusions, and the perturbation of (S)-pyraclofos on the AChE's ordered conformation was lower than that of (R)-pyraclofos, which is germane to the interaction energies of the crucial residues, e.g. Tyr-124, Tyr-337, Asp-74, Trp-86, and Tyr-119. Evidently, this attempt will contribute mechanistic information to uncovering the neurobiological effects of chiral organophosphates on the body.

Elucidating the potential neurotoxicity of chiral phenthoate: Molecular insight from experimental and computational studies.

Chiral organophosphorus pollutants are existed ubiquitously in the ecological environment, but the enantioselective toxicities of these nerve agents to humans and their molecular bases have not been fully elucidated. Using experimental and computational approaches, this story was to explore the neurotoxic response process of the target acetylcholinesterase (AChE) to chiral phenthoate and further decipher the microscopic mechanism of such toxicological effect at the enantiomeric level. The results showed that the toxic reaction of AChE with chiral phenthoate exhibited significant enantioselectivity, and (R)-phenthoate (K＝1.486 × 105 M-1) has a bioaffinity for the nerve enzyme nearly three times that of (S)-phenthoate (K＝4.503 × 104 M-1). Dynamic research outcomes interpreted the wet experiments, and the inherent conformational flexibility of the target enzyme has a great influence on the enantioselective neurotoxicological action processes, especially reflected in the conformational changes of the three key loop regions (i.e. residues His-447, Gly-448, and Tyr-449; residues Gly-122, Phe-123, and Tyr-124; and residues Thr-75, Leu-76, and Tyr-77) around the reaction patch. This was supported by the quantitative results of conformational studies derived from circular dichroism spectroscopy (α-helix: 34.7%→30.2%/31.6%; ß-sheet: 23.6%→19.5%/20.7%; turn: 19.2%→22.4%/21.9%; and random coil: 22.5%→27.9%/25.8%). Meanwhile, via analyzing the modes of toxic action and free energies, we can find that (R)-phenthoate has a strong inhibitory effect on the enzymatic activity of AChE, as compared with (S)-phenthoate, and electrostatic energy (-23.79/－17.77 kJ mol-1) played a critical role in toxicological reactions. These points were the underlying causes of chiral phenthoate displaying different degrees of enantioselective neurotoxicity.

Estimating the potential toxicity of chiral diclofop-methyl: Mechanistic insight into the enantioselective behavior.

Chiral pollutants are widely distributed in the environment; however, the enantioselective toxic effects of these chemicals have still not fully been clarified. Using wet experiments and computational toxicology, this story was to explore the static and dynamic toxic reactions between chiral diclofop-methyl and target protein at the enantiomeric level, and further unveil the microscopic mechanism of enantioselective toxicity of chiral pesticide. Steady-state and time-resolved results indicated that both (R)-/(S)-enantiomers can form the stable toxic conjugates with target protein and the bioaffinities were 1.156 × 104 M-1/1.734 × 104 M-1, respectively, and significant enantioselectivity was occurred in the reaction. Results of the modes of toxic action revealed that diclofop-methyl enantiomers located in the subdomain IIA, and the strength of important noncovalent interactions between (S)-diclofop-methyl and the residues was greater than that of (R)-diclofop-methyl. The Gibbs free energies of the chiral reactions were －26.89/－29.40 kJ mol-1 and －25.79/－30.08 kJ mol-1, respectively, which was consistent with the outcomes of photochemistry and site-specific competitive assay. Dynamic enantioselective processes explained that the impact of intrinsic protein conformational flexibility on the toxic reaction of (R)-diclofop-methyl was lower than that of (S)-diclofop-methyl, which originates from the conformational changes and spatial displacement of the four loop regions (i.e. h6↔h7, h1↔h2, h5↔h6, and h8↔h9). The quantitative data of circular dichroism spectra confirmed such results. Energy decomposition displayed that the electrostatic energy of the target protein-(S)-diclofop-methyl system (－25.86 kJ mol-1) was higher than that of the target protein-(R)-diclofop-methyl complex (－18.21 kJ mol-1). Some crucial residues such as Lys-195, Lys-199, Ser-202, and Trp-214 have been shown to be of different importance for the enantioselective toxicity of chiral diclofop-methyl. Obviously this scenario will contribute mechanistic clues to assessing the potential hazards of chiral environmental pollutants to the body.

Molecular basis for the resistance of American sloughgrass to aryloxyphenoxypropionic acid pesticides and its environmental relevance: A combined experimental and computational study.

Organic pesticides are one of the main environmental pollutants, and how to reduce their environmental risks is an important issue. In this contribution, we disclose the molecular basis for the resistance of American sloughgrass to aryloxyphenoxypropionic acid pesticides using site-directed mutagenesis and molecular modeling and then construct an effective screening model. The results indicated that the target-site mutation (Trp-1999-Leu) in acetyl-coenzyme A carboxylase (ACCase) can affect the effectiveness of the pesticides (clodinafop, fenoxaprop, cyhalofop, and metamifop), and the plant resistance to fenoxaprop, clodinafop, cyhalofop, and metamifop was found to be 564, 19.5, 10, and 0.19 times, respectively. The established computational models (i.e. wild-type/mutant ACCase models) could be used for rational screening and evaluation of the resistance to pesticides. The resistance induced by target gene mutation can markedly reduce the bioreactivity of the ACCase-clodinafop/fenoxaprop adducts, and the magnitudes are 10 and 102, respectively. Such event will seriously aggravate environmental pollution. However, the biological issue has no distinct effect on cyhalofop (RI＝10), and meanwhile it may markedly increase the bioefficacy of metamifop (RI＝0.19). We could selectively adopt the two chemicals so as to decrease the residual pesticides in the environment. Significantly, research findings from the computational screening models were found to be negatively correlated with the resistance level derived from the bioassay testing, suggesting that the screening models can be used to guide the usage of pesticides. Obviously, this story may shed novel insight on the reduction of environmental risks of pesticides and other organic pollutants.

Understanding the bioconjugation reaction of phenthoate with human serum albumin: New insights from experimental and computational approaches.

Organophosphates are chemical pollutants that are existed widely in the environment, but the reactions of these agents with blood proteins are still not fully clarified. The current story was to analyze the static and dynamic interactions between human serum albumin (HSA) and phenthoate and then uncover the impact of the conjugations on the acetylcholinesterase (AChE) activity at the microscopic scale. Experimental results revealed clearly that the bioconjugate of the HSA-phenthoate was yielded and the conformation of HSA can produce autoregulation during the reaction. Dynamic reaction processes suggested that the conformational flexibility of the specific protein domain was changed significantly in equilibrium, and the electrostatic interaction energy played a major role in total energy of the biosystems, which matches the results of wet experiment and molecular docking. We also found that the modes of homologous proteins-phenthoate have obvious distinctions, and this point is related closely to the local dynamic flexibility of biomolecular structures. Additionally, the degree of bioconjugation of the HSA-phenthoate is positively associated with the enzymatic activity of target AChE, which may be attributed to the competitive reactions between HSA and AChE. Evidently, this scenario could provide useful molecular information for the systematic exploration of the toxicokinetics of organophosphorus compounds.

Biological effects of α-adrenergic phentolamine on erythrocyte hemeprotein: Molecular insights from biorecognition behavior, protein dynamics and flexibility.

Phentolamine is one of the most representative nonselective α-adrenoreceptor blocking agents, which have been proved to be owned various pharmacological actions. Unfortunately, whether erythrocytes in the veins intervene in biological behaviors of such drug are largely obscured. With the aid of multiple biophysical techniques, this scenario was to detailed explore the potential biorecognition between phentolamine and the hemeprotein in the cytosol of erythrocytes, and the influences of dynamic characters of protein during the bioreaction. Steady-state and time-resolved fluorescence data manifested that the biomolecular recognition of phentolamine by hemeprotein was processed through the biopolymer-drug adduct with a moderate strength of 104M-1. Such procedure causes a reduction in fluorescence intensity of the aromatic tryptophan (Trp) residues, and the R-T transition of the globular protein occurred concurrently. Circular dichroism demonstrated the conclusions of fluorescence essays, viz. biorecognition can induce fairly structural transformation (self-regulation) of protein conformation. Furthermore, one could find that a specific domain for phentolamine is located at the polypeptide chains α1ß2 interface, and hydrogen bonds, π-conjugated and hydrophobic effects are discovered to be held the lowest energy state of the biomacromolecule-drug biosystem, which overtly matches the outcomes of wet experiments. Meanwhile, several crucial residues such as Trp-37 and Arg-40 were confirmed to have directly noncovalent interactions with phentolamine, and the effect of the heme group on the biomolecule-drug recognition is minimal. Further analyses of molecular dynamics simulation supported that the inherent protein flexibility may notably elicit alterations in some key noncovalent bonds between biomacromolecule and drug during the dynamic biointeraction, which might primarily be attributed to the torsion of drug structure and the conformational changes of essential residues. Undoubtedly, this research will not only help to thoroughly unearth the pharmacological profiles of phentolamine, but to elaborate the impacts of the intrinsic features (i.e. dynamics and flexibility) of critically cellular proteins on the biological conducts of active α-adrenergic blockers.

Biophysical exploration of protein-flavonol recognition: effects of molecular properties and conformational flexibility.

The current work explores the biomolecular recognition of a series of flavonols by a protein and then uncovers the influences of the structural features of flavonols and the protein's own characteristics, e.g. the dynamics and flexibility, on the bioavailability of flavonols by using the pivotal biomacromolecule hemoglobin as a model. The experimental results revealed that flavonol may lead to a notable decrease in the steady-state fluorescence intensity of the ß-37 Trp residue, and in the meantime the R-T transition of the protein transpired. Such noncovalent recognition forms the ground-state adduct, with an association intensity of 3.991 × 10(4) M(-1) in the reaction process, which has already been authenticated by the detailed analysis of time-resolved fluorescence and UV/vis absorption spectra. Furthermore, flavonol can form hydrogen bonds and π-conjugation effects with several amino acid residues on the polypeptide chain, for example, Trp-37, Arg-40, Asp-99 and Asn-102, and this event would induce self-regulation of the compact, regular conformation of the protein to a certain extent, which explicitly corroborates the results of circular dichroism. According to the study of molecular docking and structure-activity relationships, we could see that the recognition capacities of the protein-flavonols are inversely interrelated with the C log P values of the flavonol molecules. Moreover, the properties of the substituents in the structural B-ring unit of flavonols, i.e. polarity, position and number, will also prominently affect the degree of affinity and bioavailability of the protein-flavonol complexes. The analytical results of molecular dynamics (MD) simulation testified that the discussions of the structure-activity relationships are entirely logical, and the conformations of the amino acid residues forming noncovalent interactions tend to be stable in the MD simulation, as further elucidated from the dynamics data. Plainly, molecular recognition of the protein-flavonols might noticeably cause relatively large changes in protein flexibility, and then manifest different recognition strengths and corresponding biological activities. This issue will be carefully validated by the interpretation of root-mean-square fluctuation.

Measuring the bioactivity and molecular conformation of typically globular proteins with phenothiazine-derived methylene blue in solid and in solution: A comparative study using photochemistry and computational chemistry.

Methylene blue is a phenothiazine agent, that possesses a diversity of biomedical and biological therapeutic purpose, and it has also become the lead compound for the exploitation of other pharmaceuticals such as chlorpromazine and the tricyclic antidepressants. However, the U.S. Food and Drug Administration has acquired cases of detrimental effects of methylene blue toxicities such as hemolytic anemia, methemoglobinemia and phototoxicity. In this work, the molecular recognition of methylene blue by two globular proteins, hemoglobin and lysozyme was characterized by employing fluorescence, circular dichroism (CD) along with molecular modeling at the molecular scale. The recognition of methylene blue with proteins appears fluorescence quenching via static type, this phenomenon does cohere with time-resolved fluorescence lifetime decay that nonfluorescent protein-drug conjugate formation has a strength of 10(4)M(-1), and the primary noncovalent bonds, that is hydrogen bonds, π-conjugated effects and hydrophobic interactions were operated and remained adduct stable. Meantime, the results of far-UV CD and synchronous fluorescence suggest that the α-helix of hemoglobin/lysozyme decreases from 78.2%/34.7% (free) to 58.7%/23.8% (complex), this elucidation agrees well with the elaborate description of three-dimensional fluorescence showing the polypeptide chain of proteins partially destabilized upon conjugation with methylene blue. Furthermore, both extrinsic fluorescent indicator and molecular modeling clearly exhibit methylene blue is situated within the cavity constituted by α1, ß2 and α2 subunits of hemoglobin, while it was located at the deep fissure on the lysozyme surface and Trp-62 and Trp-63 residues are nearby. With the aid of computational analyses and combining the wet experiments, it can evidently be found that the recognition ability of proteins for methylene blue is patterned upon the following sequence: lysozyme

Realizing the recognition features of model antipsychotic compounds by important protein: Photochemical and computational studies.

Phenothiazine and its derivatives are the most effective antipsychotic drugs. They have been used in the treatment of serious mental and emotional symptoms including bipolar disorder, organic psychoses, psychotic depression and schizophrenia. However, these drugs cause serious side effects such as akathisia, hyperprolactinaemia and neuroleptic malignant syndrome. In this work we investigated the molecular recognition of two typical phenothiazine compounds, phenosafranin and safranin O by the most pivotal heme protein hemoglobin using steady state and time-resolved fluorescence, extrinsic 8-anilino-1-naphthalenesulfonic acid (ANS) fluorescent probe, circular dichroism (CD) along with computational modeling. Results show phenothiazines complex with protein via formation of adducts at 298 K with moderate strengths of 3.555×10(4) M(-1) and 2.567×10(4) M(-1) for the hemoglobin-phenosafranin and hemoglobin-safranin O, respectively. We also found phenothiazines were effectors at the protein allosteric site, which affects the allosteric equilibrium. Further, time-resolved fluorescence and hydrophobic ANS experiments showed the static mechanism is dominated for the shrinkage in the fluorescence intensity of ß-37 Trp residue at the α1ß2 interface. The stoichiometric proportion of the protein-drug adduct is 1:1, as derived from Job's plot. Several crucial noncovalent bonds, including hydrogen bonds, π-π stacking and hydrophobic interactions played a major role in stabilizing the noncovalent conjugates. Based on three-dimensional fluorescence, we concluded that the conformation of hemoglobin is partially destabilized after recognition with phenothiazines. These alterations were confirmed by far-UV CD spectra that showed the α-helix of protein decreased from 78.3% in free hemoglobin to 62.8% and 64.8% in hemoglobin-phenosafranin and hemoglobin-safranin O, respectively. Computer-aided molecular docking was consistent, indicating that both phenothiazines are situated within the pocket composed of α1 and ß2 subunits. Affinity of hemoglobin to phenosafranin is superior compared with safranin O. This difference may be explained by the methyl group substituent on A- and C-rings, and by the different molecular volume between phenosafranin and safranin O. Our data provides further explanation of the overall pharmacokinetics of phenothiazines and sheds light on the allosteric regulation of heme proteins.

Renal protein reactivity and stability of antibiotic amphenicols: structure and affinity.

In the present work, the molecular recognition of the oldest active amphenicols by the most popular renal carrier, lysozyme, was deciphered by using fluorescence, circular dichroism (CD) and molecular modeling at the molecular scale. Steady state fluorescence data showed that the recognition of amphenicol by lysozyme yields a static type of fluorescence quenching. This corroborates time-resolved fluorescence results that lysozyme-amphenicol adduct formation has a moderate affinity of 10(4) M(-1), and the driving forces were found to be chiefly hydrogen bonds, hydrophobic interactions and π stacking. Far-UV CD spectra confirmed that the spatial structure of lysozyme was slightly changed with a distinct reduction of α-helices in the presence of amphenicol, suggesting partial destabilization of the protein. Furthermore, via the extrinsic 8-anilino-1-naphthalenesulfonic acid fluorescence spectral properties and molecular modeling, one could see that the amphenicol binding site was situated at the deep crevice on the protein surface, and the ligand was also near to several crucial amino acid residues, such as Trp-62, Trp-63 and Arg-73. Simultaneously, contrastive studies of protein-amphenicols revealed clearly that some substituting groups, e.g. nitryl in the molecular structure of ligands, may be vitally important for the recognition activity of amphenicols with lysozyme. Due to the connection of amphenicols with fatal detrimental effects and because lysozyme has been applied as a drug carrier for proximal tubular targeting, the discussion herein is necessary for rational antibiotic use, development of safe antibiotics and particularly a better appraisal of the risks associated with human exposure to toxic agrochemicals.

Evaluation of the biointeraction of colorant flavazin with human serum albumin: insights from multiple spectroscopic studies, in silico docking and molecular dynamics simulation.

Azo compounds are the largest chemical class of agents frequently used as colorants in a variety of consumer goods and farm produce; therefore, they may become a hazard to public health, because numerous azo compounds and their metabolites are proven to be carcinogens and mutagens. Herein several qualitative and quantitative analytical techniques, including steady state and time-resolved fluorescence, circular dichroism (CD), computer-aided molecular docking as well as molecular dynamics simulation, were employed to ascertain the molecular recognition between the principal vehicle of ligands in human plasma, albumin and a model azo compound, flavazin. The results show that the albumin spatial structure was changed in the presence of flavazin with a decrease of α-helix suggesting partial protein destabilization/self-regulation, as derived from steady state fluorescence, far-UV CD and detailed analyses of three-dimensional fluorescence spectra. Time-resolved fluorescence further evinced that the recognition mechanism is related to albumin-flavazin adduct formation with an association intensity of 10(4) M(-1), and the driving forces were found to be chiefly π-π interactions, hydrophobic interactions and hydrogen bonds. The specific binding domain of flavazin in protein was defined from molecular docking; subdomain IIA (Sudlow's site I) was found to retain high affinity for the ligand flavazin. This finding corroborates the results of competitive ligand displacement experiments, a hydrophobic 8-anilino-1-naphthalenesulfonic acid probe study and protein denaturation results, placing flavazin at the warfarin-azapropazone site. Based on molecular dynamics simulation, it can be said with certainty that the results of molecular docking are credible, and the key amino acid residues participating in the molecular recognition of flavazin by protein are clearly Trp-214, Arg-222 and Lys-436. The outcomes presented here will help to further comprehend the molecular recognition of azo compounds by protein and the possible toxicological profiles of other compounds that have configurations analogous to azo chemicals.

Bioavailability and activity of natural food additive triterpenoids as influenced by protein.

Triterpenoids were thought to be biologically ineffective for a very long time, but aggregating proof on their widely ranging pharmacological activities paired with a dubious toxicity portrait has motivated regenerated attraction for human health and disease. In the current contribution, our central goal was to integratively dissect the biointeraction of two typical triterpenoids, ursolic acid and oleanolic acid, by the most fundamental macromolecule bovine serum albumin (BSA) by employing molecular modeling, steady state and time-resolved fluorescence, and circular dichroism spectra at the molecular scale. Based on molecular modeling, subdomain IIA, which matches Sudlow's site I, was allocated to retain high affinity for triterpenoids, but the affinity of ursolic acid with subdomain IIA is somewhat inferior compared to that of oleanolic acid, probably because the affinity differentiation arises from the different positions of the methyl group on the E-ring in the two triterpenoids. This sustains the site-specific ligands, and hydrophobic 8-anilino-1-naphthalenesulfonic acid probe results in arranging the triterpenoids at the warfarin-azapropazone site. The data of steady state and time-resolved fluorescence indicated that the recognition of triterpenoids by BSA produced quenching by a static type, in other words, the ground state BSA-triterpenoid complex formation with the affinities of 1.507/1.734, 1.042/1.186, and 0.8395/0.9863 × 10(4) M(-1) at 298, 304, and 310 K for ursolic acid/oleanolic acid, respectively. Thermodynamic analyses show that the basic forces acting between BSA and triterpenoids are hydrogen bonds, van der Waals forces, and hydrophobic interactions; this occurrence provoked the alterations of the BSA spatial structure with a noticeable decline of α-helix evoking perturbation of the protein, as stemmed from circular dichroism, synchronous fluorescence, and three-dimensional fluorescence measurements. We anticipate that the complexation of plant triterpenoids with protein delineated here may be exploited as a biologically relevant model for evaluating the physiologically applicable noncovalent complexes in in vivo examination of triterpenoid properties such as accumulation, bioavailability, and distribution.

Molecular recognition of malachite green by hemoglobin and their specific interactions: insights from in silico docking and molecular spectroscopy.

Malachite green is an organic compound that can be widely used as a dyestuff for various materials; it has also emerged as a controversial agent in aquaculture. Since malachite green is proven to be carcinogenic and mutagenic, it may become a hazard to public health. For this reason, it is urgently required to analyze this controversial dye in more detail. In our current research, the interaction between malachite green and hemoglobin under physiological conditions was investigated by the methods of molecular modeling, fluorescence spectroscopy, circular dichroism (CD) as well as hydrophobic ANS displacement experiments. From the molecular docking, the central cavity of hemoglobin was assigned to possess high-affinity for malachite green, this result was corroborated by time-resolved fluorescence and hydrophobic ANS probe results. The recognition mechanism was found to be of static type, or rather the hemoglobin-malachite green complex formation occurred via noncovalent interactions such as π-π interactions, hydrogen bonds and hydrophobic interactions with an association constant of 10(4) M(-1). Moreover, the results also show that the spatial structure of the biopolymer was changed in the presence of malachite green with a decrease of the α-helix and increase of the ß-sheet, turn and random coil suggesting protein damage, as derived from far-UV CD and three-dimensional fluorescence. Results of this work will help to further comprehend the molecular recognition of malachite green by the receptor protein and the possible toxicological profiles of other compounds, which are the metabolites and ramifications of malachite green.

Binding patterns and structure-affinity relationships of food azo dyes with lysozyme: a multitechnique approach.

Food dyes serve to beguile consumers: they are often used to imitate the presence of healthful, colorful food produce such as fruits and vegetables. But considering the hurtful impact of these chemicals on the human body, it is time to thoroughly uncover the toxicity of these food dyes at the molecular level. In the present contribution, we have examined the molecular reactions of protein lysozyme with model food azo compound Color Index (C.I.) Acid Red 2 and its analogues C.I. Acid Orange 52, Solvent Yellow 2, and the core structure of azobenzene using a combination of biophysical methods at physiological conditions. Fluorescence, circular dichroism (CD), time-resolved fluorescence, UV-vis absorption as well as computer-aided molecular modeling were used to analyze food dye affinity, binding mode, energy transfer, and the effects of food dye complexation on lysozyme stability and conformation. Fluorescence emission spectra indicate complex formation at 10(-5) M dye concentration, and this corroborates time-resolved fluorescence results showing the diminution in the tryptophan (Trp) fluorescence mainly via a static type (KSV = 1.505 × 10(4) M(-1)) and Förster energy transfer. Structural analysis displayed the participation of several amino acid residues in food dye protein adducts, with hydrogen bonds, π-π and cation-π interactions, but the conformation of lysozyme was unchanged in the process, as derived from fluorescence emission, far-UV CD, and synchronous fluorescence spectra. The overall affinity of food dye is 10(4) M(-1) and there exists only one kind of binding domain in protein for food dye. These data are consistent with hydrophobic probe 8-anilino-1-naphthalenesulfonic acid (ANS) displacement, and molecular modeling manifesting the food dye binding patch was near to Trp-62 and Trp-63 residues of lysozyme. On the basis of the computational analyses, we determine that the type of substituent on the azobenzene structure has a powerful influence on the toxicity of food dyes. Results from this work testify that model protein, though an indirect method, provides a more comprehensive profile of the essence of toxicity evaluation of food dyes.