Determination of fuzzy logic membership functions using genetic algorithms: application to structure–odor modeling

Fuzzy logic has been used as a tool in structure–camphoraceous odor relationships. The data base studied included 99 molecules. The rules used to discriminate between camphor and non camphor molecules lead to 77% correct discrimination. Such rules account for the shape and the size of the molecule. Their adjustment by means of genetic algorithms led to 84% correct discrimination between camphor and non-camphor molecules.Figure Membership function for the chosen variables     

A computational approach to the synthesis of dirithromycin

Dirithromycin is a macrolide antibiotic derived from erythromycin A. Dirithromycin is synthesized by the condensation of 9(S)-erythromycylamine with 2-(2-methoxyethoxy)-acetaldehyde. To gain insight into the synthesis, the condensation mechanism has been analyzed computationally by the AM1 method in the gas phase. First, the formation of the Schiff bases of dirithromycin and epidirithromycin from 9(S)-erythromycylamine and 2-(2-methoxyethoxy)-acetaldehyde were modeled. Then, the tautomerization of the Schiff bases to dirithromycin and epidirithromycin were considered. Finally, the epimerization of the Schiff base of epidirithromycin to the Schiff base of dirithromycin was investigated. Our results show that, even though carbinolamine forms faster for epidirithromycin than the corresponding structure for dirithromycin, dirithromycin is the major product of the synthesis. Figure Synthesis of dirithromycin     

Quantitative structure–activity relationship study on some benzodiazepine derivatives as anti-Alzheimer agents

A QSAR study was performed in an attempt to explore the pharmacophore of some benzodiazepine derivatives as anti-Alzheimer agents for the inhibition of -secretase. The study, which used the electrotopological state atom (ETSA) index, which encodes electronic and topological information, reveals the importance of two phenyl rings—one substituted and another unsubstituted, for the inhibition of the enzyme. Fluorine substitution on the substituted phenyl ring has an important contribution to the activity. R-configurations of the aliphatic chain substituents provide the exact conformation of the molecules to enter into the binding pockets of the receptor(s).Figure General structure of benzodiazepine containing -secretase inhibitors     

An ab initio study of di- and trifluorobenzene–benzene complexes as relevant to carbonic anhydrase II–drug interactions

Ab initio calculations at the MP2/6-31G* level have shown that variously substituted di- and trifluorobenzenes form non-covalent complexes with benzene that adopt either aromatic–aromatic or H---F binding, the choice being determined by the pattern of fluorination. The binding energies of these structures are from 3.4 to 4.5 kcal mol^–1. This range is large enough to account for observed variations in the binding affinity of a library of fluoroaromatic inhibitors of carbonic anhydrase. This enzyme has an aromatic amino acid at a central position in the active site. The diverse modes of binding of the dimers also suggest that aggregates of fluorobenzenes might adopt specified 3-dimensional shapes in the solid state. Figure Observed structure for 1,2-difluorobenzene     

Characterization of molecular structures and properties of polyurethanes using molecular dynamics simulations

Thermotropic polyurethanes with mesogenic groups in side chains were prepared from two diisocyanates and four diols with stoichiometric ratios of reactive isocyanate (NCO) and hydroxy (OH) groups. Their thermal behavior was determined by differential scanning calorimetry. The effect of structure modifications of the diisocyanates and diols, in particular changes in the mesogen, were investigated. Introduction of mesogenic segments into the polymers suppresses the ordering. Stiff end substituents (phenyl and alkoxy groups) of the mesogens stabilize the mesophases to such an extent that the negative influence of long polymer chains is compensated and the liquid-crystalline properties are recovered. All-atom molecular dynamics simulations in the Cerius² modeling environment were carried out to characterize the structures of the polymers. Analysis of the dynamic trajectories at 20, 100, 120 and 170 °C revealed changes in conformation of macromolecules, which correlate with DSC measurements.Figure Example of structure relaxation of D4/TDI molecule at indicated simulation times (temperature 20 °C): a complete structure; b backbone structure; c top view of molecule     

Theoretical studies of the conformational behavior of chain molecules containing polar groups: simulations of a poly(vinylidene fluoride) model

Atomistic Monte Carlo simulations have been conducted to elucidate the conformational behavior of a single chain molecule containing polar functional groups. Here, we resort to an atomistic poly(vinylidene fluoride) (PVDF) chain model as a representative example. The model is modified in such a way that bond lengths and bond angles are fixed, aiming to manifest the role of dipolar interactions. For a given chain length, chain conformation is sensitive to two environmental parameters, temperature and dielectric constant. The mean chain size increases when temperature and/or dielectric constant are increased. The conformational behavior is further characterized by chain size distribution function, and our findings show that temperature induced conformational transition for a chain molecule can be discrete or continuous, depending on its chain length. Also, the dipolar interactions in PVDF are effectively attractive, and enhance chain contraction. As a result, when the strength of dipolar interactions is increased, the discrete conformational transition shifts toward longer chains; and for a given chain length, such a transition occurs at higher temperatures.Figure Variation of R² with temperature for different dielectric constants =1 and 8, denoted by dotted and solid lines, respectively, and for different chain lengths M=8 and 12, as marked. Lines are meant for eye guidance     

Predicting anti-HIV activity of 2,3-diaryl-1,3-thiazolidin-4-ones: computational approach using reformed eccentric connectivity index

The eccentric connectivity index, which has recently been employed successfully for the development of numerous mathematical models for the prediction of biological activities of diverse nature, has been reformed to overcome its limitations caused by degeneracy and insensitivity towards heteroatoms. The reformed eccentric connectivity index, termed the eccentric connectivity topochemical index, overcomes the limitations of the eccentric connectivity index by exhibiting very low degeneracy and displaying sensitivity to both the presence and relative position of heteroatoms without compromizing the discriminating power of the eccentric connectivity index. The relationship of the eccentric connectivity topochemical index, eccentric connectivity index and Wieners index with regard to the anti-HIV activity of 2, 3-diaryl-1, 3-thiazolidin-4-one derivatives was subsequently investigated. The values of the eccentric connectivity topochemical index, the eccentric connectivity index and Wieners index of each of 31 analogues comprizing the data set were computed using in-house computer program. Resultant data was analyzed and suitable models developed after identification of active ranges. Subsequently, each derivative was assigned a biological activity using these models, which was then compared with the reported anti-HIV activity. The accuracy of prediction using these models was found to vary from 81 to 90%. The proposed index offers a vast potential for virtual screening of combinatorial libraries, structure property/activity studies and drug design.Figure Basic structure of 2,3-diaryl-1, 3-thiazoidin-4-ones.     

Reinvestigation of molecular structure and barrier to internal rotation of pyridinium N-phenolate betaine dye

In the present paper, the results of a systematic theoretical study of the molecular structure of 4-(1-pyridinium-1-yl)phenolate betaine are reported. The ground-state molecular structure and the barrier to internal rotation of the betaine dye molecule were calculated ab inito (with Hartree–Fock theory and the second-order of Möller–Plesset method) and with density functional theory (DFT). In order to estimate the complete basis set limit, the calculations of barriers to internal rotations were performed using correlation–consistent basis sets with a maximal cardinal number of four. It was determined that electron correlation is crucial in order to obtain reliable geometries and rotational barriers of the molecule investigated. For the sake of comparison, the results of calculations using the AM1 Hamiltonian are also presented.Figure Investigated betaine dye.     

AIPAR: ab initio parametrization of intermolecular potentials for computer simulations

An unambiguous, fully ab initio and automated technique denoted AIPAR (ab initio parametrization) implemented in the SJBR program has been proposed to yield intermolecular interaction potentials between polar molecules and water. The AIPAR procedure has been applied to several organic molecules covering a wide range of structure and functional groups, namely methanol, acetone (propanone), methanethiol (methyl mercaptan), imidazole (1,3-diazole), oxazole and furan. The AIPAR-derived sets of parameters compare well with the empirical OPLS ones, mainly when the all-atoms model is employed in the OPLS procedure. Monte Carlo simulations were performed for an aqueous solution of methanol and for an equimolar binary mixture methanol–water using the AIPAR and OPLS parameters. The thermodynamic and geometric results obtained with the parameters obtained with the AIPAR procedure compare favorably with the OPLS simulations, even for the binary mixture, demonstrating the precision, robustness and transferability of the parameters obtained with the AIPAR procedure.Figure Superimposed configurations of water (without the hydrogen atoms) around the methanol molecule obtained with the AGOA procedure.     

Asialo-α1-acid glycoprotein resialylated with 9-amino-5-N-acetyl-d-neuraminic acid is resistant towards bacterial,viral and mammalian sialidases

We demonstrate that 9-amino-NeuAc transferred to asialo-₁-acid glycoprotein resists cleavage by bacterial, viral and mammalian sialidases. This is the first synthetic sialic acid analogue, which can be activated and transferred to glycoprotein, but is not a sialidase (EC 3.2.1.18) substrate.Abbreviations HPLC high performance liquid chromatography - BSA bovine serum albumin - NeuAc N-acetyl-d-neuraminic acid, 5-acetamido-3,5-dideoxy-d-glycero-d-galacto-non-2-ulosonic acid - 9-Amino-NeuAc 9-amino-5-N-acetyl-d-neuraminic acid, 5-acetamido-9-trideoxy-d-glycero-d-galacto-non-2-ulosonic acid - CMP-NeuAc cytidine-5-monophospho-N-acetyl-d-neuraminic acid - CMP-9-amino-NeuAc cytidine-5-monophospho-9-amino-5-N-acetyl-d-neuraminic acid - 9-azido-NeuAc 5-acetamido-9-azido-3,5,9-trideoxy-d-glycero-d-galacto-non-2-ulosonic acid. Enzymes EC 3.2.1.18 sialidase, acylneuraminylhydrolase - EC 2.4.99.1 Galß1-4GlcNAc a(2-6)-sialytransferase     

Observations on phase-locking within the response of primary muscle spindle afferents to pseudo-random stretch

In order to uncover encoder properties of primary muscle spindle afferent fibers, time coupling (phase-locking) of action potentials on cyclic muscle stretch was studied by means of pseudo-random noise. In cats Ia action potentials were recorded from dorsal root filaments and the gastrocnemius muscles of one hind leg were stretched. The stimulus time course was a determined sequence of randomly varying muscle length which could be applied repeatedly (sequence duration 0.6 or 20 s). The noise amplitude (standard deviation of displacements) was varied between 5 and 300 m, the upper cut-off frequency of noise f _c was varied between 20 and 100 Hz. The responses to the consecutive pseudo-random noise cycles were displayed as raster diagrams and cycle histograms. Phaselocking characterized the responses at all noise amplitudes outside the near threshold range (>10 m). The higher and f _c , the stronger was the phase-locking of impulses on the stretch. When and f _c were selected to achieve high mean stretch velocities of about 500 mm/s, phase-locking was as precise as 0.15 ms, measured as the variability of spike occurrences with respect to stretch. The rasters obtained with low noise amplitudes (<40 m) showed a loose phase-locking and this gave insight into underlying mechanisms: The elicitation of action potentials caused by dynamic stretch can be prevented by a post-spike depression of excitability. This disfacilitation was very effective in counteracting weak stretch components within the random sequence and less effective or even missing when relatively strong stretch components could force the spike elicitation. This led to the reestablishment of phase-locked patterns. The results were discussed in relation to the known encoder models.     

Lethal effect of lambda DNA terminase in recombination deficient Escherichia coli

Summary DNA terminase is the enzyme that catalyses the cleavage of DNA concatemers into genome-size molecules and packages them into the capsid. The cleavage (DNA maturation) takes place in a specific site in the phage DNA called cos. Either one of two Escherichia coli proteins, integration host factor (IHF) and terminase host factor (THF), is required, in addition to terminase, for maturation of wild-type DNA in vitro. In vivo, at least some cos cleavage is known to occur in mutants that are unable to synthesize active IHF. No THF-defective mutants have yet been isolated. In order to determine if IHF, THF or any other host protein is involved in DNA maturation in vivo, I devised a selection for host mutants that are unable to support cos cleavage. The selection is based on the assumption that DNA terminase will kill cells by cleaving chromosomally located cos sites. I found that DNA terminase will indeed kill cells provided that they contain a chromosomal cos site and provided also that they are defective in the host recA or recB genes. These two genes are required for certain pathways of genetic recombination and repair of damaged DNA, and I suggest that they prevent terminase-induced killing by repairing broken chromosomes. Interstingly, mutation in a related host gene, recD, did not render cells susceptible to terminase killing. recD and recB both encode subunits of exonuclease V, but recD mutants, unlike recB, remain proficient in genetic recombination and repair. I found mutants that survived the lethal effect of terminase in cos-containing E. coli recA at a frequency of about 5×10^-5. About 90% of these survivors were defective in terminase synthesis, and the rest were defective in IHF function. This result suggests that in the absence of IHF in vivo cos cleavage decreases to a level that permits repair of the damage, and therefore survival, even in recombination deficient cells. The absence of mutations in any other host gene suggests that IHF is the major accessory factor in DNA maturation in vivo. Alternatively, or in addition, mutations in other accessory factors are lethal.Abbreviations gp gene product: e.g. gpA, product of gene A - () prophage state - [] plasmid-carrier state     

Modeling of adsorption in nanopores

Adsorption in nonporous materials has been studied using Grand Canonical Monte Carlo simulations. We discuss three types of materials: (a) a model of cylindrical pores with smooth walls, representing MCM-41 like materials, (b) a model of cylindrical pores with regular structured walls (model of carbon nanotubes) and (c) a material with crystalline wall structure (zeolites). Typical problems related to the stability of adsorbed layers have been analyzed. We have shown that the mechanism of adsorption is strongly dependent on the structure of the pore walls. In the case of amorphous walls it may lead to metastable configurations. In nanotubes, the ordered corrugation structure of walls determines the low temperature structure of the adsorbed system. In 3D ordered porous system, such as zeolites, the mechanism of adsorption is mostly determined by characteristic sites of adsorption.Figure Adsorbed atoms and energy fluctuations at the pressure of the first layer formation of krypton atoms: (a) instantaneous numbers of adsorbed atoms (per nm² of the pore wall) as a function of the time of simulation (Monte Carlo steps) observed in a relatively long run, (b) the bimodal distribution of the energy fluctuations is a consequence of the behavior of the systems as shown in (a).     

An extended aqueous solvation model based on atom-weighted solvent accessible surface areas: SAWSA v2.0 model

A new method is proposed for calculating aqueous solvation free energy based on atom-weighted solvent accessible surface areas. The method, SAWSA v2.0, gives the aqueous solvation free energy by summing the contributions of component atoms and a correction factor. We applied two different sets of atom typing rules and fitting processes for small organic molecules and proteins, respectively. For small organic molecules, the model classified the atoms in organic molecules into 65 basic types and additionally. For small organic molecules we proposed a correction factor of hydrophobic carbon to account for the aggregation of hydrocarbons and compounds with long hydrophobic aliphatic chains. The contributions for each atom type and correction factor were derived by multivariate regression analysis of 379 neutral molecules and 39 ions with known experimental aqueous solvation free energies. Based on the new atom typing rules, the correlation coefficient (r) for fitting the whole neutral organic molecules is 0.984, and the absolute mean error is 0.40 kcal mol^–1, which is much better than those of the model proposed by Wang et al. and the SAWSA model previously proposed by us. Furthermore, the SAWSA v2.0 model was compared with the simple atom-additive model based on the number of atom types (NA). The calculated results show that for small organic molecules, the predictions from the SAWSA v2.0 model are slightly better than those from the atom-additive model based on NA. However, for macromolecules such as proteins, due to the connection between their molecular conformation and their molecular surface area, the atom-additive model based on the number of atom types has little predictive power. In order to investigate the predictive power of our model, a systematic comparison was performed on seven solvation models including SAWSA v2.0, GB/SA_1, GB/SA_2, PB/SA_1, PB/SA_2, AM1/SM5.2R and SM5.0R. The results showed that for organic molecules the SAWSA v2.0 model is better than the other six solvation models. For proteins, the model classified the atoms into 20 basic types and the predicted aqueous free energies of solvation by PB/SA were used for fitting. The solvation model based on the new parameters was employed to predict the solvation free energies of 38 proteins. The predicted values from our model were in good agreement with those from the PB/SA model and were much better than those given by the other four models developed for proteins.Figure The definition of hydrophobic carbons. Here CA, CB and CD are three carbon atoms; X represents a heteroatom. According to our definition, CB is a hydrophobic carbon, CA is not a hydrophobic carbon because a heteroatom is within four atoms and CD is not a hydrophobic carbon because CD is sp²- hydridized and in a six-member ring.Electronic Supplementary Material Supplementary material is available for this article at     

GTP-binding proteins regulate high conductance anion channels in rat bile duct epithelial cells

Summary Epithelial cells from the intrahepatic bile duct contribute to bile formation, but little is known of the cellular mechanisms responsible. In these studies, we have characterized the endogenous GTP-binding proteins (G proteins) present in these cells and evaluated their role in regulation of high conductance anion channels. G proteins were identified in purified plasma membranes of isolated bile duct epithelial cells using specific antisera on Western blots, and ion channel activity was measured in excised inside-out membrane patches using patch-clamp recording techniques. In patches without spontaneous channel activity, addition of cholera toxin to the cytoplasmic surface had no effect (n=10). Addition of pertussis toxin caused an activation of channels in 13/34 (38%) attempts, as detected by an increase in channel open probability. Activated channels were anion selective (gluconate/Cl^– permeability ratio of 0.17±0.04) and had a unitary conductance of 380 pS. Channel open probability was also increased by the nonhydrolyzable GDP analogue guanosine 5-0-(2-thiodiphosphate) in 8/14 (57%) attempts. In contrast, channel open probability was rapidly and reversibly decreased by the nonhydrolyzable analogue of GTP 5guanylylimidodiphosphate in 7/9 (78%) attempts. Western blotting with specific antisera revealed that both G_i –2 and G_i –3 were present in significant amounts, whereas G_i –1 and G_o were not detected. These studies indicate that in bile duct epithelial cells, high conductance anion channels are inhibited, in a membrane-delimited manner, by PTX-sensitive G proteins.We gratefully acknowledge the assistance of Marwan Farouk, M.D. in the preparation of bile duct epithelial cells, Lucy Seger in the identification of the G proteins, C.F. Starmer in channel analysis, and P.J. Casey for the gift of bacteria expressing the different G-protein -subunits. This work was supported in part by grants from the National Institutes of Health DK43278 (to J.G.F.), DK42486 (to T.W.G.), and DK07568 (to J.M.M.); American Gastroenterological Association/G.D. Searle Research Scholar Award (to J.G.F.) and an American Gastroenterological Association Advanced Research Training Award (J.M.M.).     

Direct ab initio dynamics studies of the hydrogen abstraction reactions of hydrogen atom with n-propyl radical and isopropyl radical

The kinetics of the hydrogen abstraction reactions of hydrogen atom with n-propyl radical and isopropyl radical were studied using the direct ab initio dynamics approach. BHandHLYP/cc-pVDZ method was employed to optimize the geometries of stationary points as well as the points on the minimum energy path (MEP). The energies of all the points for the two reactions were further refined at the QCISD(T)/cc-pVTZ level of theory. No barrier was found at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory for both reactions. The forward and reverse rate constants were evaluated with both canonical variational transition state theory (CVT) and microcanonical variational transition state theory ( VT) in the temperature range of 300–2,500 K. The fitted three-parameter Arrhenius expression of the calculated CVT rate constants at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory are k^CVT (n – C₃H₇)=1.68×10^–14 T^0.84 e^(319.5/T) cm³ molecule^–1 s^–1 and k^CVT (iso-C₃H₇)=4.99×10^–14 T^0.90 e^(159.5/T) cm³ molecule^–1 s^–1 for reactions of n-C₃H₇ + H and iso-C₃H₇ + H, respectively, which are in good agreement with available literature data. The variational effects were analysed.Figure Comparison of the calculated forward rate constants at the QCISD(T)/cc-pVTZ//BHandHLYP/cc-pVDZ level of theory and the available experimental and theoretical data of the reaction vs 1,000/T for the two reactions.     

Binding affinity of hydroxamate inhibitors of matrix metalloproteinase-2

Here we report molecular dynamics (MD) and free energy perturbation (FEP) simulations applied to hydroxamate-matrix metalloproteinase-2 (MMP-2) complex systems. We have developed some new force field parameters for the hydroxamate functional group that were not included in the AMBER94 force field but were necessary in our simulations. For the representation of the active zinc center, a bonded model was adopted in which restrained electrostatic potential fitting (RESP) charges were used as the electrostatic representation of this model. Using the resulted bonded model, FEP simulations predict the relative binding free energy in good agreement with the experimental value. By analyzing the molecular dynamics (MD) trajectories of the two complex systems, we can provide an explanation of why one of the two inhibitors is favored over the other. The results provide a chemical insight into the interactions between inhibitor and enzyme, and can indicate changes in the inhibitor that would enhance inhibitor–enzyme interactions.Figure The scheme of the binding site     

Normal coordinate analyses of 3,5-dichlorophenylcyanamide

The structure of 3,5-dichlorophenylcyanamide c-C₆H₃Cl₂–NHCN was investigated by DFT-B3LYP and ab initio MP2 calculations with the 6-311+G** basis set. The planar to perpendicular rotational barrier was calculated to be of about 5 kcal mol^–1 at both levels of calculation. The stability of the planar structure of the molecule was explained on the basis of conjugation effects between the cyanamide–NHCN moiety and the phenyl c-C₆H₅ ring in agreement with earlier NMR results. The CNC and the HNC bond angles were calculated to be about 120° especially by MP2 calculation, which is consistent with sp² (planar –NH–CN group) and not sp³ (pyramidal –NH–CN group) structure. The vibrational frequencies of the d₀, d₁ and d₃ species of 3,5-dichlorophenylcyanamide and the potential energy distributions among symmetry coordinates of the normal modes of the parent molecule were computed at the DFT-B3LYP level. The calculated infrared and Raman spectra of the molecule were plotted. Complete vibrational assignments were made on the basis of isotopic substitution and normal coordinate calculations.Figure Potential curves for the internal rotation in 3,5-dichlorophenylcyanamide as determined by DFT-B3LYP/6-311+G** (solid) and MP2/6-311+G** (dotted) calculations     

The Mutation Spectrum of the CFTR Gene in Cystic Fibrosis Patients from Bashkortostan

Mutations of CFTR were studied in patients with cystic fibrosis (CF) from Bashkortostan. In total, 15 mutations were observed and 51% of all mutant alleles identified. The most diagnostically significant mutations were delF508 (33.8%), 394delTT (3.52%), CFTRdele2,3(21kb) (1.41%), R334W (1.41%), 3849 + 10kbC T (1.41%), and N1303K (1.41%). Mutations G542X, 2184insA, S1196X, and W1282X were each found in less than 1% patients. Five new mutations and two neutral substitutions were revealed. These were I488M (exon 10), 1811 + 12A C (intron 11), T663S (exon 13), I1226R (exon 19), 4005 + 9A C (intron 20), 2097A C (A655A, exon 13), and 3996G C (V1288V, exon 20). Bashkortostan was shown to differ in the CFTR mutation spectrum from other regions of Russia. The results will allow direct DNA diagnostics of CF in far more families. Molecular screening of probands" relatives will contribute to identification and medical genetic counseling of heterozygous carriers, which is essential for CF prevention.     

A quantitative structure–antifungal activity relationship study of oxygenated aromatic essential oil compounds using data structuring and PLS regression analysis

Twenty two oxygenated aromatic essential oil compounds were chosen for the study of the antifungal activity against two wood-decaying fungi, the white-rot Trametes versicolor, which mainly metabolizes lignin, and the brown-rot Coniophoha puteana, which digests cellulose in plant cell walls. Minimal inhibitory concentrations (MICs) were determined by the agar dilution method, using dimethyl sulfoxide (DMSO) as the solvent for the selected compounds and potato-dextrose agar (PDA) as the growth medium for both fungi. The MICs were then used to generate a tree structure, which represents the structuring of the essential oil compounds by the nature and position of the substituents in their aromatic rings, and as dependent variables (log(1/MIC)) in the QSAR analysis. Data structuring proved that a relationship between the molecular structures of the essential oil compounds and their antifungal activity exists, and the hypotheses derived therefrom were complemented by performing a QSAR analysis using the partial least squares (PLS) method. Statistically significant PLS models were obtained with the 1-octanol–water partition coefficient (C log P), the energy of the highest occupied molecular orbital (E _HOMO), and the number of hydrogen-bond donor atoms in the molecules of the compounds studied (Donor) for T. versicolor and with C log P and the fractional negative surface area (FNSA1) for C. puteana.Figure Tree structure representing the structuring of the oxygenated aromatic essential-oil compounds by the position and nature of their substituents in the aromatic ring