Insight into catalytically relevant correlated motions in human purine nucleoside phosphorylase

The catalytic site of the homotrimeric enzyme human purine nucleoside phosphorylase enzyme (hPNP) features residue F200 and the 241-265 loop directly skirting the purine base and a residue belonging to the adjacent monomer, F159, immediately conterminous to the ribosyl moiety. Crystallographic B-factors of apo human purine nucleoside phosphorylase, and hPNP complexed with substrate/transition state (TS) analogues, show that residue E250 is the centroid of a highly mobile loop region. Furthermore, superimposition of apo hPNP and hPNP complexed with TS analogue Immucillin-H shows a tightening of the active site, caused by the ligand-dependent 241-265 loop rearrangement taking place upon substrate/inhibitor binding, suggesting a putative dynamic role of the loop in binding/catalysis. However, crystallographic structures reveal only average atomic positions, and more detailed information is needed to discern the dynamic behavior of hPNP. The Essential Dynamics (ED) method is used here to investigate the existence of correlated motions in hPNP and consequently proposes mutagenesis assays to estimate the relative importance of these motions in the phosphorolytic efficiency of the reaction catalyzed by hPNP. We compare the concerted motions obtained from multiple molecular dynamics simulations of apo and Michaelis complex of hPNP both in vacuo and in solution. The results of the principal component analysis for the apo hPNP indicate the existence of strong correlations predominantly in the vicinity of residue F159. However, for the Michaelis complex, concerted motions are seen mostly around both active site residue F200 and loop residue E250. Additionally, for a simulation depicting the relaxation of tight complexed hPNP with a TS analogue, toward its relaxed apo form (after removal of the TS analog), a combination of the apo hPNP and Michaelis complex motions is found, with prominent concerted modes centered around neighboring residues F159, F200, and E250. Finally, we probed the extent to which these concerted motions bear an intrinsic catalytic role by performing experimental site-directed mutagenesis on some residues, followed by kinetic analysis. The F159G and F200G mutants displayed a strong increase in K(M) and modest decrease in k(cat), suggesting that these concerted motions may provide dynamical roles in substrate binding and/or catalysis. However, further structural data for the hPNP mutants are needed to confirm our hypothesis.     

NMR and EPR structural delineation of copper(II) complexes formed by kanamycin A in water

The complexes formed by kanamycin A at three different pH values (5.5, 7.4 and 12.0) were investigated by NMR and EPR spectroscopy. Paramagnetic relaxation contributions to proton relaxation rates were measured using a combination of the TOCSY sequence with the inversion recovery experiment in order to gain signal resolution in the bulk region. Measured contributions were converted into distances and used for structural determination by restrained simulated annealing where all possible chair and boat conformations of the rings were taken into account. The interaction of the Cu(II) ion with the nitrogen of the C ring is apparent at all pH values. At higher pH also the amino group of ring A starts to be involved in the metal coordination sphere. This is accompanied by a switch in conformation of ring C. Structures are consistent with the involvement in the coordination sphere either of the 2' or 4' hydroxyl oxygens at pH 5.5 and the 5 and the 6' hydroxyl oxygens (or the ring oxygen) at pH 12.0.     

Proton-coupled electron transfer versus hydrogen atom transfer in benzyl/toluene,methoxyl/methanol,and phenoxyl/phenol self-exchange reactions

Degenerate hydrogen atom exchange reactions have been studied using calculations, based on density functional theory (DFT), for (i) benzyl radical plus toluene, (ii) phenoxyl radical plus phenol, and (iii) methoxyl radical plus methanol. The first and third reactions occur via hydrogen atom transfer (HAT) mechanisms. The transition structure (TS) for benzyl/toluene hydrogen exchange has C(2)(h)() symmetry and corresponds to the approach of the 2p-pi orbital on the benzylic carbon of the radical to a benzylic hydrogen of toluene. In this TS, and in the similar C(2) TS for methoxyl/methanol hydrogen exchange, the SOMO has significant density in atomic orbitals that lie along the C-H vectors in the former reaction and nearly along the O-H vectors in the latter. In contrast, the SOMO at the phenoxyl/phenol TS is a pi symmetry orbital within each of the C(6)H(5)O units, involving 2p atomic orbitals on the oxygen atoms that are essentially orthogonal to the O.H.O vector. The transferring hydrogen in this reaction is a proton that is part of a typical hydrogen bond, involving a sigma lone pair on the oxygen of the phenoxyl radical and the O-H bond of phenol. Because the proton is transferred between oxygen sigma orbitals, and the electron is transferred between oxygen pi orbitals, this reaction should be described as a proton-coupled electron transfer (PCET). The PCET mechanism requires the formation of a hydrogen bond, and so is not available for benzyl/toluene exchange. The preference for phenoxyl/phenol to occur by PCET while methoxyl/methanol exchange occurs by HAT is traced to the greater pi donating ability of phenyl over methyl. This results in greater electron density on the oxygens in the PCET transition structure for phenoxyl/phenol, as compared to the PCET hilltop for methoxyl/methanol, and the greater electron density on the oxygens selectively stabilizes the phenoxyl/phenol TS by providing a larger binding energy of the transferring proton.     

Exploring nucleoside hydrolase catalysis in silico: molecular dynamics study of enzyme-bound substrate and transition state

The mechanism of action of inosine-uridine nucleoside hydrolase has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Five MD studies have been performed with enzyme substrate complex (E.S), enzyme substrate complex with protonated His241 (EH.S), enzyme transition state complex (E.TS), enzyme transition state complex with protonated His241 (EH.TS), and His241Ala transition state complex E(H241A).TS. Special attention has been given to the role of His241, which has been considered as the general acid catalyst to assist departure of the leaving nucleobase on the basis of its location in the active site in the X-ray crystal structure (). Yet on the basis of the location in the active site, Tyr229 is closer to the aniline ring of pAPIR as compared to His241. On initiation of MD simulations, His241 does not approach the nucleobase in the structures of EH.S, E.S, EH.TS, and E.TS. In the solvated enzyme, Tyr229, which is a member of the hydrogen bonding network inosine O2'.Asp14.His241.Tyr229.inosine N7, serves as a proton source to the leaving nucleobase. The loss of significant activity of His241Ala mutant is shown to be related to the disruption of the above hydrogen bonded network and the distancing of Tyr229 from inosine N7. The structures of the enzyme complexes with substrate or TS are not visibly altered on protonation of His241, a most unusual outcome. The bell-shaped pH dependence upon pK(app)'s of 7.1 and 9.1 may be attributed to the necessity of the dissociation of Asp10 or Asp15 and the acid form of Tyr229, respectively. In TS, the residue Ile81 migrated closer, whereas Arg233 moved away from the nucleobase. The probability of ribooxocarbenium ion stabilization by Asn168 and Asp14 is discussed. The Asp14-CO(2)(-) is hydrogen bonded to the ribose 2'-OH for 96% of the MD simulation time. Nucleophilic addition of water138 to ribooxocarbenium ion is suggested to be assisted by the proton shuttle from water138 --> Asp10 --> Asp15 --> water pool. An anticorrelation motion between Tyr229-OH and Asn168-OD1 in EH.S and E.S is observed. The relationship of this anticorrelated motion to mechanism, if any, deserves further exploration, perhaps the formation of a near attack conformation.     

Theoretical Analysis of Structural Characteristics of Morin Rearrangement for Cephalosporins Prepared from Penicillin Sulfoxide

The structure of penicillin sulfoxide is rearranged to cephalosporins by the Morin rearrangement. It is a unit reaction for the preparation of various types of cephalosporins. In order to make better use of the reaction and in view of the shortage of the reaction theory, this study used m062x/6-311++G(d, p) to explore the possible ring-opening reaction of the penicillin sulfoxide. It is found that the isomer of(S)-sulfoxide is a necessary structure. At the same time, the intramolecular hydrogen bonding effect between the side-chain amide proton(-CONH-) and the sulfinyl oxygen(-SO) is the decisive structure factor for the formation of alkenyl in ring-opening reaction, and the best reaction path is S0- TS2- IN1 channel. The main effect of acid catalysis is to catalyze the dehydration reaction of sulfenic acid to form sulfur cations for subsequently ring closing reaction.     

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.     

Catalytic mechanism of glycosyltransferases: hybrid quantum mechanical/molecular mechanical study of the inverting N-acetylglucosaminyltransferase I

The Golgi glycosyltransferase, N-acetylglucosaminyltransferase I (GnT-I), catalyzes the transfer of a GlcNAc residue from the donor UDP-GlcNAc to the C2-hydroxyl group of a mannose residue in the trimannosyl core of the Man5GlcNAc2-Asn-X oligosaccharide. The catalytic mechanism of GnT-I was investigated using a hybrid quantum mechanical/molecular mechanical (QM/MM) method with a QM part containing 88 atoms treated with density functional theory (DFT) at the BP/TZP level. The remaining parts of a GnT-I complex, altogether 5633 atoms, were modeled using the AMBER molecular force field. A theoretical model of a Michaelis complex was built using the X-ray structure of GnT-I in complex with the donor having geometrical features consistent with kinetic studies. The QM(DFT)/MM model identified a concerted SN2-type of transition state with D291 as the catalytic base for the reaction in the enzyme active site. The TS model features nearly simultaneous nucleophilic addition and dissociation steps accompanied by the transfer of the nucleophile proton Hb2 to the catalytic base D291. The structure of the TS model is characterized by the Ob2-C1 and C1-O1 bond distances of 1.912 and 2.542 A, respectively. The activation energy for the proposed reaction mechanism was estimated to be approximately 19 kcal mol-1. The calculated alpha-deuterium kinetic isotope effect of 1.060 is consistent with the proposed reaction mechanism. Theoretical results also identified interactions between the Hb6 and beta-phosphate oxygen of the UDP and a low-barrier hydrogen bond between the nucleophile and the catalytic base D291. It is proposed that these interactions contribute to a stabilization of TS. This modeling study provided detailed insight into the mechanism of the GlcNAc transfer catalyzed by GnT-I, which is the first step in the conversion of high mannose oligosaccharides to complex and hybrid N-glycan structures.     

Coupling of electron and proton transport in photosynthetic membranes: molecular mechanism

Using the method of Modified Neglect of Diatomic Overlap (MNDO), the electronic structure of plastoquinol (PQH(2)) and plastoquinone (PQ) in neutral, singly (PQ(-)) and doubly (PQ(2-)) reduced states is studied. The conformational analysis performed on these molecules shows that in the lowest energy conformation, the angle between the first link of the tail backbone and the ring plane of neutral and singly reduced PQ and plastoquinol is nearly the same and differs by 15 degrees from that of doubly reduced PQ. Nevertheless, for all states of plastoquinone and for plastoquinol, the total energy changes by less than 0.2 eV when the studied angle is varied from 0 degrees to 180 degrees. As in Rhodobacter sphaeroides, the oxygen of the PQ ring is reported to form a hydrogen bond with a nitrogen in the ring of Histidine (His) L 190. The energy of the PQ-His complex was calculated for different redox states of PQ and for several values of the distance between the molecules (N-O distance from 0.2 to 0.5 nm). For every considered complexes, the total energy dependence on the proton position on the line connecting the N and O atoms was determined, to see if the hydrogen bond is formed. It is shown that for only singly reduced PQ this dependence has a symmetric two-well form, i.e. the hydrogen bond is formed. For neutral and doubly reduced PQ, the curve is also two-well but asymmetric, so that the proton is bound to His or to PQ, correspondingly.On the basis of these results, we propose the following scheme of electron-proton coupling. Negatively charged oxygens of PQ form H-bonds with proton donor groups of the surrounding protein and fix PQ in its pocket. While the negative charges of oxygens increase after quinone reduction, protons shift to PQ oxygens and form strong hydrogen bonds with them. Upon second PQ reduction, protons are torn away from surrounding amino acids and form covalent bonds with the quinol. Resulting PQH(2) detaches from its binding place and is replaced by a neutral PQ. The lacking protons on amino acids in the Q(B) pocket are replaced by a step-by-step transfer from the stroma bulk through the proton channels.     

Novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: free energy calculation for the Beckmann rearrangement promoted by proton transfers in the supercritical water

The Beckmann rearrangement of acetone oxime promoted by proton transfers in the supercritical water has been investigated by means of the hybrid quantum mechanical/molecular mechanical approach combined with the theory of energy representation (QM/MM-ER) recently developed. The transition state (TS) structures have been explored by ab initio calculations for the reaction of hydrated acetone oxime on the assumption that the reaction is catalyzed by proton transfers along the hydrogen bonds connecting the solute and the solvent water molecules. Up to two water molecules have been considered as reactants that take part in the proton transfers. As a result of the density functional theory calculations with B3LYP functional and aug-cc-pVDZ basis set, it has been found that participation of two water molecules in the reaction reduces the activation free energy by -12.3 kcal/mol. Furthermore, the QM/MM-ER simulations have revealed that the TS is more stabilized than the reactant state in the supercritical water by 2.7 kcal/mol when two water molecules are involved in the reaction. Solvation free energies of the reactant and the TS have been decomposed into terms due to the electronic polarization of the solute, electron density fluctuation, and others to elucidate the origin of the stabilization of the TS as compared with the reactant. It has been revealed that the promotion of the chemical reaction due to the hydration mainly originates from the interaction between the nonpolarized solute and the solvent water molecules at the supercritical state.     

Exploring the optimal site for modifications of pyranmycins with the extended arm approach

[reaction: see text] Continuing from the syntheses and the antibacterial studies of a library of pyranmycins, we further probed the proximity around ring III of pyranmycin by introducing an "extended arm" that has hydroxyethyl or aminoethyl groups at the O-2' ', O-3' ', or O-4' ' positions. The results from the antibacterial studies reveal the optimal structural motif is the attachment of an extended arm with a terminal hydroxyl group at the O-3' ' position.     

Binding thermodynamics and interaction patterns of human purine nucleoside phosphorylase-inhibitor complexes from extensive free energy calculations

Journal of Computer-Aided Molecular Design - Human purine nucleoside phosphorylase (hPNP) plays a significant role in the catabolism of deoxyguanosine. The trimeric protein is an important target...     

Microporous crystal 12CaO x 7Al(2)O(3) encaging abundant O(-) radicals

Extremely high concentrations (>1020 cm-3) of active oxygenic radicals, O- and O2-, have been created in the zeolitic crystal, 12CaO.7Al2O3 (C12A7), which can accommodate anions in its cavities. An increase in oxygen pressure and a decrease in water vapor pressure at high temperature enhance the formation of the radicals. The oxidation of Pt is observed on the surface of the material as a result of reaction with the active oxygens.     

乙基乙酰胺热解反应机理的MO研究

对于烷基乙酰胺的初始热解反应机理, 通常认为与酯类的热解反应相类似。Maccoll和Nagra通过对该热解反应的动力学研究, 认为两者存在不同。差异之一, 烷基乙酰胺存在两种可能的热解途径【参见本文(129页)前述反应方程(1),(2)】。而在酯类热解反应中(2)的活化能过高, 且四元环过渡态极不稳定。差异之二, 极性溶剂(比如乙酸)对酰胺热解反应的催化作用不明显, 而对酯类等气相热解反应的催化作用是十分显著的。为此, 我们用MINDO/3分子轨道法对乙基乙酰胺的初始热解反应进行了较全面的研究。用能量梯度法对此反应的反应物、中间体和生成物的平衡几何构型进行了全优化。(如图1所示)用极小能量途径法分别寻找反应(1)和反应(2)的初始过渡态, 继而用Powell法全优化过渡态的几何构型, 计算所得的过渡态TS1、TS2和TS3分别见图2a, 图3a和图4a。为了确证这些过渡态, 进行了振动分析研究, 结果表明这些过渡态的力常数矩阵的诸本征值中均只有一个负值, 且虚振动模式展示了走向各自的反应物和生成物的趋势, (如图2b,图3b和图4b所示)。它们的总能量及反应(1)和反应(2)的活化能列于表1. 对整个热解反应(1)作了内禀反应坐标(IRC)理论分析, 反应历程见图5所示. 与IRC相应的总能偶极矩以及部分关键的键长和原子净电荷变化一并列于表2.本文研究结果表明, 在乙基乙酰胺的初始反应中主反应即反应(1)与酯类反应相类似, ...     

Identification of the transition states in the inversion of 1,4-benzodiazepines with the ab initio replica path method

The inversion of four 1,4-benzodiazepines was investigated with the ab initio "replica path method" with density functional theory at the B3LYP/6-31G* level. The reaction path constructed with this method for the inversion provides an approximate transition state (TS) geometry, which, upon further TS optimization, leads to the TS geometry characterized by a single vibrational frequency. 1,4-Benzodiazepines lacking a 5-phenyl ring have a single reaction path for the inversion with Cs symmetry at the TS. In contrast, the inversion of benzodiazepines with a 5-phenyl ring, such as the peripheral benzodiazepine receptor ligand 4'-chlorodiazepam (Ro5-4864) and its N1-desmethyl analog (Ro5-2752), can proceed through multiple reaction paths having a TS with or without Cs symmetry. Notably, the replica path method found a path connecting two asymmetric TSs of 4'-chlorodiazepam via a symmetrical TS, suggesting that these inversion paths can be readily crossed from one to another. The stabilization energies gained by 4'-chlorodiazepam and its N1-desmethyl analog from the breaking of Cs symmetry at the TS were calculated to be 0.10 and 0.07 kcal/mol, respectively. The origin of the broken symmetry of Cs was traced to the coupling of the puckering of the diazepine ring with the rotation of the chlorophenyl ring. These results show the advantages of the replica path method for locating the TSs as well as for constructing the reaction paths for the inversion of 1,4-benzodiazepines.     

Brønsted Acid Catalyzed Morita–Baylis–Hillman Reaction: A New Mechanistic View for Thioureas Revealed by ESI‐MS(/MS) Monitoring and DFT Calculations

A Morita–Baylis–Hillman (MBH) reaction catalyzed by thiourea was monitored by ESI‐MS(/MS) and key intermediates were intercepted and characterized. These intermediates suggest that thiourea acts as an organocatalyst in all steps of the MBH reaction cycle, including the rate‐limiting proton‐transfer step. DFT calculations, performed for a model MBH reaction between formaldehyde and acrolein with trimethylamine as base and in the presence or the absence of thiourea, suggest that thiourea accelerates MBH reactions by decreasing the transition‐state (TS) energies through bidentate hydrogen bonding throughout the whole catalytic cycle. In the rate‐limiting proton‐transfer step, the thiourea acts not as a proton shuttle, but as a Brønsted acid stabilizing the basic oxygen center that is formed in the TS.     

Computational study of phosphatase activity in soluble epoxide hydrolase: high efficiency through a water bridge mediated proton shuttle

Recently, a new branch of fatty acid metabolism has been opened by the novel phosphatase activity found in the N-terminal domain of the, hence bifunctional, soluble epoxide hydrolase (sEH). Importantly, this finding has also provided a new site for drug targeting in sEH's activity regulation. Classical MD and hybrid Car-Parrinello QM/MM calculations have been performed to investigate the reaction mechanism of the phosphoenzyme intermediate formation in the first step of the catalysis. The results support a concerted multi-event reaction mechanism: (1) a dissociative in-line nucleophilic substitution for the phosphoryl transfer reaction; (2) a double proton transfer involved in the formation of a good leaving group in the transition state. The presence of a water bridge in the substrate/enzyme complex allowed an efficient proton shuttle, showing its key role in speeding up the catalysis. The calculated free energy of the favored catalytic pathway is approximately 19 kcal/mol, in excellent agreement with experimental data.     

19F NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA‐Catalyzed GTP Hydrolysis

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use ¹⁹F NMR chemical shifts in the MgF₃^? transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ‐PO₃^? oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ‐PO₃^? group based on the structure of a RhoA/GAP‐GDP‐MgF₃^? TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85′ side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln‐63 carboxamide by Arg85′ to preclude direct hydrogen bonding from water to the target γ‐PO₃^? group.     

Theoretical studies for excited-state tautomerization in the 7-azaindole-(CH3OH)n (n = 1 and 2) complexes in the gas phase

The excited-state tautomerization of 7-azaindole (7AI) complexes bonded with either one or two methanol molecule(s) was studied by systematic quantum mechanical calculations in the gas phases. Electronic structures and energies for the reactant, transition state (TS), and product were computed at the complete active space self-consistent field (CASSCF) levels with the second-order multireference perturbation theory (MRPT2) to consider the dynamic electron correlation. The time-dependent density functional theory (TDDFT) was also used for comparison. The excited-state double proton transfer (ESDPT) in 7AI-CH(3)OH occurs in a concerted but asynchronous mechanism. Similarly, such paths are also found in the two transition states during the excited-state triple proton transfer (ESTPT) of the 7AI-(CH(3)OH)(2) complex. In the first TS, the pyrrole ring proton first migrated to methanol, while in the second the methanol proton moved first to the pyridine ring. The CASSCF level with the MRPT2 correction showed that the former path was much preferable to the latter, and the ESDPT is much slower than the ESTPT. Additionally, the vibrational-mode enhanced tautomerization in the 7AI-(CH(3)OH)(2) complex was also studied. We found that the excitation of the low-frequency mode shortens the reaction path to increase the tautomerization rate. Overall, most TDDFT methods used in this study predicted different TS structures and barriers from the CASSCF methods with MRPT2 corrections.     

Cα-Cβ and Cα-N bond cleavage in the dissociation of protonated N-benzyllactams: dissociative proton transfer and intramolecular proton-transport catalysis

In mass spectrometry of protonated N-benzylbutyrolactams, the added proton is initially localized on the carbonyl oxygen, which is the thermodynamically preferred protonation site. Upon collisional activation, dissociative proton transfer takes place leading to the occurrence of fragmentation reactions. The major fragmentations observed are the cleavages of C(α)-C(β) and C(α)-N bonds on the two sides of the methylene linker, which is different to the cleavage of the amide bond itself seen in most amide cases. Theoretical calculations and isotopic labeling experiments demonstrate that the phenyl ring regulates the proton transfer reactions. The proton directly migrates to the C(β) position via a 1,5-H shift leading to the efficient loss of benzene, while it stepwise migrates to the amide nitrogen resulting in the formation of a benzyl cation. The stepwise proton transfer is achieved via intramolecular proton-transport catalysis. The C(γ) position accepts the proton from the carbonyl oxygen via a 1,6-H shift, and then donates it to the amide nitrogen via a 1,4-H shift. The general 1,3-H shift from the carbonyl oxygen to the amide nitrogen can be excluded in this case due to its significant energy barrier. The substituent effects are also applied to explore the reaction mechanism, and it proves that both C(β) and C(γ) are involved in the dissociative proton transfer processes. For monosubstituted N-benzylbutyrolactams, the abundance ratios of the two competing product ions are well correlated with the nature of the substituents.