Ligand recognition by chloroperoxidase using molecular interaction fields and quantum chemistry calculations

We recently reported that chloroperoxidase (CPO) from Caldariomyces fumago showed atypical kinetic behavior for the oxidation of 4,6 dimethyl dibenzothiophene (DMDBT). In this work, we undertake the theoretical study of DMDBT docking into CPO's active site, in order to clarify its binding capacity and affinity using molecular interaction fields and quantum mechanical procedure. The results revealed that CPO has two substrate binding sites with similar affinities for DMDBT. This finding suggests that the atypical kinetic behavior of CPO for the oxidation of DMDBT might be due to the simultaneous binding of two DMDBT molecules during its reaction cycle. Finally, we extend these results to carbazole, a nitrogen-containing PAH, through experimental and theoretical studies.     

Control of adenosylmethionine-dependent radical generation in biotin synthase: a kinetic and thermodynamic analysis of substrate binding to active and inactive forms of BioB

Biotin synthase (BS) is an AdoMet-dependent radical enzyme that catalyzes the insertion of sulfur into saturated C6 and C9 atoms in the substrate dethiobiotin. To facilitate sulfur insertion, BS catalyzes the reductive cleavage of AdoMet to methionine and 5'-deoxyadenosyl radicals, which then abstract hydrogen atoms from the C6 and C9 positions of dethiobiotin. The enzyme from Escherichia coli is purified as a dimer that contains one [2Fe-2S]2+ cluster per monomer and can be reconstituted in vitro to contain an additional [4Fe-4S]2+ cluster per monomer. Since each monomer contains each type of cluster, the dimeric enzyme could contain one active site per monomer, or could contain a single active site at the dimer interface. To address these possibilities, and to better understand the manner in which biotin synthase controls radical generation and reactivity, we have examined the binding of AdoMet and DTB to reconstituted biotin synthase. We find that both the [2Fe-2S]2+ cluster and the [4Fe-4S]2+ cluster must be present for tight substrate binding. Further, substrate binding is highly cooperative, with the affinity for AdoMet increasing >20-fold in the presence of DTB, while DTB binds only in the presence of AdoMet. The stoichiometry of binding is ca. 2:1:1 AdoMet:DTB:BS dimer, suggesting that biotin synthase has a single functional active site per dimer. AdoMet binding, either in the presence or in the absence of DTB, leads to a decrease in the magnitude of the UV-visible absorption band at approximately 400 nm that we attribute to changes in the coordination environment of the [4Fe-4S]2+ cluster. Using these spectral changes as a probe, we have examined the kinetics of AdoMet and DTB binding, and propose an ordered binding mechanism that is followed by a conformational change in the enzyme-substrate complex. This kinetic analysis suggests that biotin synthase is evolved to bind AdoMet both weakly and slowly in the absence of DTB, while both the rate of binding and the affinity for AdoMet are increased in the presence of DTB. Cooperative binding of AdoMet and DTB may be an important mechanism for limiting the production of 5'-deoxyadenosyl radicals in the absence of the correct substrate.     

Substrate and Inhibitor-induced Dimerization and Cooperativity in Caspase-1 but Not Caspase-3

Caspases are intracellular cysteine-class proteases with aspartate specificity that is critical for driving processes as diverse as the innate immune response and apoptosis, exemplified by caspase-1 and caspase-3, respectively. Interestingly, caspase-1 cleaves far fewer cellular substrates than caspase-3 and also shows strong positive cooperativity between the two active sites of the homodimer, unlike caspase-3. Biophysical and kinetic studies here present a molecular basis for this difference. Analytical ultracentrifugation experiments show that mature caspase-1 exists predominantly as a monomer under physiological concentrations that undergoes dimerization in the presence of substrate; specifically, substrate binding shifts the K_D for dimerization by 20-fold. We have created a hemi-active site-labeled dimer of caspase-1, where one site is blocked with the covalent active site inhibitor, benzyloxycarbonyl-Val-Ala-Asp-fluoromethylketone. This hemi-labeled enzyme is about 9-fold more active than the apo-dimer of caspase-1. These studies suggest that substrate not only drives dimerization but also, once bound to one site in the dimer, promotes an active conformation in the other monomer. Steady-state kinetic analysis and modeling independently support this model, where binding of one substrate molecule not only increases substrate binding in preformed dimers but also drives the formation of heterodimers. Thus, the cooperativity in caspase-1 is driven both by substrate-induced dimerization as well as substrate-induced activation. Substrate-induced dimerization and activation seen in caspase-1 and not in caspase-3 may reflect their biological roles. Whereas caspase-1 cleaves a dramatically smaller number of cellular substrates that need to be concentrated near inflammasomes, caspase-3 is a constitutively active dimer that cleaves many more substrates located diffusely throughout the cell.     

Influence of pH on the Mn2+ activation of and binding to yeast enolase: a functional study.

The influence of pH on the activation of yeast enolase by Mn2+ was measured by steady-state kinetics. The pH influence on the binding of Mn2+ to apoenolase and the enolase-substrate complex was measured by EPR spectroscopy. At pH values above 6.6, activation by Mn2+ is fit by Michaelis-Menten kinetics, but at higher concentrations of Mn2+, inhibition is observed. Under conditions analogous to the kinetic studies, the enzyme binds two Mn2+ per dimer with a Kd in the micromolar range. In the presence of the substrate 2-phosphoglycerate, three thermodynamically distinct cation binding sites per monomer are detected and the binding constants are determined by a fit to the data. As the pH decreases, the reaction velocity decreases and the cation inhibition becomes minimal. Under these conditions, only two Mn2+ binding sites per monomer are observed; the third site must be the inhibitory site. The velocity and kinetic constants are minimally affected by buffer except at pH 5.8 with PIPES. Under these conditions, the velocity is only about 40% that observed with other buffers and only a single binding site for Mn2+ per monomer is detected in the presence or absence of substrate. A direct role in the catalytic mechanism by the second cation is called to question. The binding constant for Mn2+ at site I is independent of pH over the range from 7.5 to 5.2, and the binding at site II increases only slightly over this same pH range. These results indicate that the cation sites at positions I and II contain ligands that are pH independent over this range.(ABSTRACT TRUNCATED AT 250 WORDS)     

Cooperative effects of CTP on calf liver CTP synthetase.

In all previous kinetics studies of calf liver CTP synthetase, simple Michaelis-Menten hyperbolic plots were obtained. In this study it was shown that calf liver CTP synthetase could generate sigmoidal kinetic plots as a function of the substrate UTP when in the presence of the product of the reaction, CTP. The Hill number was estimated to be 2.8. The enzyme did not generate sigmoidal plots as a function of the other substrates (L-glutamine and ATP) either in the presence or absence of CTP. Thus, CTP apparently induced changes in the liver enzyme which altered the binding of UTP to the enzyme by acting at a site distinct from the UTP binding site (allosteric site). This concept was further strengthened by the fact that 3-deazaUTP, a known competitive inhibitor of the liver enzyme, did not induce sigmoidal kinetic plots. It was also shown that CTP had no effect upon the dimerization of the enzyme, thus ruling out monomer to dimer transitions as a potential mechanism for the observed sigmoidal kinetics.     

Preparation of monomeric cytochrome C oxidase: its kinetics differ from those of the dimeric enzyme

Bovine cytochrome c oxidase in 0.1% dodecylmaltoside, 50 mM KCl and 10 mM Tris-HCl, pH 7.4 is monodisperse with an apparent Mr 360,000 (dimer) as estimated by filtration on Ultrogel AcA 34. In the absence of added KCl the apparent Mr is 160,000 (monomer). The dimeric enzyme has a high and a low affinity site for cytochrome c; the monomeric, only the high affinity site. The results are consistent with the existence of one active site per monomer, having high affinity for cytochrome c. Since in a dimer the two sites are in close proximity, the binding of the first molecule of cytochrome c to the first site hinders the binding of the second molecule to the second site. The kinetic data fit with a model of homotropic negative cooperativity. The effect of salts on the cytochrome c oxidase kinetics is also present in isolated bovine heart mitochondria.     

Kinetics and thermodynamics of ligand binding by cytochrome P450 3A4

Cytochrome P450 (P450) 3A4, the major catalyst involved in human drug oxidation, displays substrate- and reaction-dependent homotropic and heterotropic cooperative behavior. Although several models have been proposed, these mainly rely on steady-state kinetics and do not provide information on the contribution of the individual steps of P450 catalytic cycle to the observed cooperativity. In this work, we focused on the kinetics of substrate binding, and the fluorescent properties of bromocriptine and alpha-naphthoflavone allowed analysis of an initial ligand-P450 3A4 interaction that does not cause a perturbation of the heme spectrum. The binding stoichiometry for bromocriptine was determined to be unity using isothermal titration calorimetry and equilibrium dialysis methods, suggesting that the ligand bound to the peripheral site during the initial encounter dissociates subsequently. A three-step substrate binding model is proposed, based on absorbance and fluorescence stopped-flow kinetic data and equilibrium binding data obtained with bromocriptine, and evaluated using kinetic modeling. The results are consistent with the substrate molecule binding at a site peripheral to the active site and subsequently moving toward the active site to bind to the heme and resulting in a low to high spin iron shift. The last step is attributed to a conformational change in the enzyme active site. The later steps of binding were shown to have rate constants comparable with the subsequent steps of the catalytic cycle. The P450 3A4 binding process is more complex than a two-state system, and the overlap of rates of some of the events with subsequent steps is proposed to underlie the observed cooperativity.     

Tyrosyl-tRNA synthetase acts as an asymmetric dimer in charging tRNA. A rationale for half-of-the-sites activity

Tyrosyl-tRNA synthetase from Bacillus stearothermophilus is a classical example of an enzyme with half-of-the-sites activity. The enzyme crystallizes as a symmetrical dimer that is composed of identical subunits, each having a complete active site. In solution, however, tyrosyl-tRNA synthetase binds tightly, and activates rapidly, only 1 mol of Tyr/mol of dimer. It has recently been shown that the half-of-the-sites activity results from an inherent asymmetry of the enzyme. Only one subunit catalyzes formation of Tyr-AMP, and interchange of activity between subunits is not detectable over a long time scale. Paradoxically, however, the kinetics of tRNA charging are biphasic with respect to [Tyr], suggesting that both subunits of the dimer are catalytically active. This paradox has now been resolved by kinetic analysis of heterodimeric enzymes containing different mutations in each subunit. Biphasic kinetics with unchanged values of KM for Tyr are maintained when one of the two tRNA-binding domains is removed and also when the affinity of the "inactive" site for Try is reduced by 2-58-fold. The biphasic kinetics do not result from catalysis at both active sites, but instead appear to result from two molecules of Tyr binding sequentially to the same site. A second molecule of Tyr perhaps aids the dissociation of Tyr-tRNA by displacing the tyrosyl moiety from its binding site. A monomer of the enzyme is probably too small to allow both recognition and aminoacylation of a tRNA molecule. This could explain the requirement for the enzyme to function as an asymmetric dimer.     

Kinetic and physical properties of Co2+ enolase

The activation of yeast enolase by cobaltous ion in 0.1 M KCl is characterized by an activation constant of 1 microM and an inhibition constant of 18 microM. Measurements of binding of Co2+ to the apoenzyme show that a maximum of four Co2+ ions are bound per dimer in the presence or absence of substrate although binding is far tighter in the presence of substrate. Ultraviolet spectral titrations show evidence for a conformational change due exclusively to the binding of the first two ions of Co2+. Both visible and EPR spectra confirm that the environment of the first pair of cobalt ions ("conformational sites") is markedly different from that of the second pair in the "catalytic" sites. Cobalt at the conformational site appears to be a tetragonally distorted octahedral complex while the second pair of metal ions appears to be in a more regular tetrahedral symmetry. Addition of either Mg2+ or substrate to the enzyme with only one pair of cobalt ions per dimer causes striking changes in the metal ion environment. The conformational metal sites appear sufficiently shielded from solvent to be inaccessible to oxidation by H2O2, in contrast to the second pair of cobaltous ions whose ready oxidation by H2O2 inactivates the enzyme. Comparison of kinetic and binding data suggests that only one site of the dimeric enzyme can be active, since activity requires more than two metals bound per dimer and inactivation results from the binding of the fourth ion per dimer.     

Catalytic acid-base groups in yeast pyruvate decarboxylase. 3. A steady-state kinetic model consistent with the behavior of both wild-type and variant enzymes at all relevant pH values.

The widely quoted kinetic model for the mechanism of yeast pyruvate decarboxylase (YPDC, EC 4.1.1.1), an enzyme subject to substrate activation, is based on data for the wild-type enzyme under optimal experimental conditions. The major feature of the model is the obligatory binding of substrate in the regulatory site prior to substrate binding at the catalytic site. The activated monomer would complete the cycle by irreversible decarboxylation of the substrate and product (acetaldehyde) release. Our recent kinetic studies of YPDC variants substituted at positions D28 and E477 at the active center necessitate some modification of the mechanism. It was found that enzyme without substrate activation apparently is still catalytically competent. Further, substrate-dependent inhibition of D28-substituted variants leads to an enzyme form with nonzero activity at full saturation, requiring a second major branch point in the mechanism. Kinetic data for the E477Q variant suggest that three consecutive substrate binding steps may be needed to release product acetaldehyde, unlikely if YPDC monomer is the minimal catalytic unit with only two binding sites for substrate. A model to account for all kinetic observations involves a functional dimer operating through alternation of active sites. In the context of this mechanism, roles are suggested for the active center acid-base groups D28, E477, H114, and H115. The results underline once more the enormous importance that both aromatic rings of the thiamin diphosphate, rather than only the thiazolium ring, have in catalysis, a fact little appreciated prior to the availability of the 3-dimensional structure of these enzymes.     

Evidence of dimerisation among class D beta-lactamases: kinetics of OXA-14 beta-lactamase

OXA-14 enzyme, a class D beta-lactamase, gave biphasic kinetics with all penicillin and cephalosporin substrates tested, such that the catalytic rate declined more swiftly than was explicable by substrate depletion. This biphasic behaviour was independent of temperature or extraneous protein but was lost if the enzyme was diluted to occupy almost the total assay volume before addition of a small amount of concentrated substrate. The presence of substrate could partially protect the enzyme against conversion to the less active form, with protection greatest at substrate concentration above the K(m). These observations are compatible with the hypothesis that the biphasic kinetics depended on the enzyme existing as a highly active dimer at high concentration and as a less active monomer at low concentration. Direct evidence supporting this hypothesis came from the observation that gel exclusion chromatography indicated a higher molecular weight for concentrated enzyme than for dilute. Biphasic kinetics are not so universal for different substrates amongst beta-lactamases (OXA-10, -11, -13, -16 and -17) that differ from OXA-14 by only one to two amino acid substitutions. It may be that the monomer:dimer equilibrium is more rapidly achieved with these enzymes than with OXA-14, or that the kinetic properties of the dimers and monomers of these enzymes are similar, masking any biphasic trait.     

Peroxide oxidation of indole to oxindole by chloroperoxidase catalysis.

In the presence of chloroperoxidase, indole was oxidized by H2O2 to give oxindole as the major product. Under most conditions oxindole was the only product formed, and under optimal conditions the conversion was quantitative. This reaction displayed maximal activity at pH 4.6, although appreciable activity was observed throughout the entire pH range investigated, namely pH 2.5-6.0. Enzyme saturation by indole could not be demonstrated, up to the limit of indole solubility in the buffer. The oxidation kinetics were first-order with respect to indole up to 8 mM, which was the highest concentration of indole that could be investigated. On the other hand, 2-methylindole was not affected by H2O2 and chloroperoxidase, but was a strong inhibitor of indole oxidation. The isomer 1-methylindole was a poor substrate for chloroperoxidase oxidation, and a weak inhibitor of indole oxidation. These results suggest the possibility that chloroperoxidase oxidation of the carbon atom adjacent to the nitrogen atom in part results from hydrogen-bonding of the substrate N-H group to the enzyme active site.     

Modulation of dimer stability in yeast pyrophosphatase by mutations at the subunit interface and ligand binding to the active site

Yeast (Saccharomyces cerevisiae) pyrophosphatase (Y-PPase) is a tight homodimer with two active sites separated in space from the subunit interface. The present study addresses the effects of mutation of four amino acid residues at the subunit interface on dimer stability and catalytic activity. The W52S variant of Y-PPase is monomeric up to an enzyme concentration of 300 microm, whereas R51S, H87T, and W279S variants produce monomer only in dilute solutions at pH > or = 8.5, as revealed by sedimentation, gel electrophoresis, and activity measurements. Monomeric Y-PPase is considerably more sensitive to the SH reagents N-ethylmaleimide and p-hydroxymercurobenzosulfonate than the dimeric protein. Additionally, replacement of a single cysteine residue (Cys(83)), which is not part of the subunit interface or active site, with Ser resulted in insensitivity of the monomer to SH reagents and stabilization against spontaneous inactivation during storage. Active site ligands (Mg(2+) cofactor, P(i) product, and the PP(i) analog imidodiphosphate) stabilized the W279S dimer versus monomer predominantly by decreasing the rate of dimer to monomer conversion. The monomeric protein exhibited a markedly increased (5-9-fold) Michaelis constant, whereas k(cat) remained virtually unchanged, compared with dimer. These results indicate that dimerization of Y-PPase improves its substrate binding performance and, conversely, that active site adjustment through cofactor, product, or substrate binding strengthens intersubunit interactions. Both effects appear to be mediated by a conformational change involving the C-terminal segment that generally shields the Cys(83) residue in the dimer.     

Regulation and state of aggregation of Bacillus subtilis prephenate dehydratase in the presence of allosteric effectors.

Prephenate dehydratase from Bacillus subtilis was found to exist in three states of aggregation. A high molecular weight (210,000) species was fully active and the catalytic activity was unaffected by the effectors methionine or phenylalanine. Low concentrations of phenylalanine caused dissociation to a Mr = 55,000 dimer. Heating to 32 degrees C also caused dissociation, but cooling and adding substrate or methionine favored association. When no effectors were present the enzyme eluted from Sephadex columns as a monomer. Both methionine and phenylalanine shifted the equilibrium from the inactive monomer to the active dimeric enzyme. In the presence of a saturating methionine concentration, the dimer possessed the same high activity as did the 210,000-dalton form. Phenylalanine inhibited the dimer, but not the higher molecular weight form. A model involving only three types of sites (catalytic, association-activation, and inhibition) is consistent with the data. It is proposed that phenylalanine is the preferred metabolite for binding both effector sites on the dimer; it binds the association-activation site with higher affinity than the inhibition site, but binding at the latter site has a greater effect on the catalytic rate. Methionine, like phenylalanine, has a hydrophobic side chain but is accommodated only at the association-activation site.     

Reaction pathway of bovine aortic lysyl oxidase

The catalysis of amine oxidation by lysyl oxidase has been probed to assess for the likely order of substrate binding and product release and to discriminate between mechanistic alternatives previously proposed for other copper-dependent amine oxidases using molecular oxygen as a substrate. Lineweaver-Burk plots revealed a pattern of parallel lines when the oxidation of n-butylamine was followed at different fixed concentrations of oxygen consistent with a "ping-pong" kinetic mechanism in which the aldehyde is produced and released before the binding of oxygen, the second substrate. Initial burst experiments revealed the ability of lysyl oxidase to form and release n-butyraldehyde in amounts stoichiometric with functional active site content in the absence of oxygen, consistent with the ping-pong kinetics obtained. Reciprocal plots of n-butylamine oxidation in the presence of fixed concentrations of the reaction products were consistent with a Uni Uni Uni Bi ping-pong kinetic mechanism with the aldehyde being the first, H2O2 the second, and ammonia the last departing product. Moreover, spectral studies of the oxidation of p-hydroxybenzylamine by lysyl oxidase indicated that the enzyme does not process the amine substrate to a noncovalently bound p-hydroxybenzaldimine intermediate subsequently to be hydrolyzed to p-hydroxybenzaldehyde. The kinetic mechanism of lysyl oxidase thus appears to be similar to those described for diamine oxidase and pig plasma monoamine oxidase.     

Allosteric interactions within subsites of a monomeric enzyme: kinetics of fluorogenic substrates of PI-specific phospholipase C

Two novel water-soluble fluorescein myo-inositol phosphate (FLIP) substrates, butyl-FLIP and methyl-FLIP, were used to examine the kinetics and subsite interactions of Bacillus cereus phosphatidylinositol-specific phospholipase C. Butyl-FLIP exhibited sigmoidal kinetics when initial rates are plotted versus substrate concentration. The data fit a Hill coefficient of 1.2-1.5, suggesting an allosteric interaction between two sites. Two substrate molecules bind to this enzyme, one at the active site and one at a subsite, causing an increase in activity. The kinetic behavior is mathematically similar to that of well-known cooperative multimeric enzymes even though this phosphatidylinositol-specific phospholipase C is a small, monomeric enzyme. The less hydrophobic substrate, methyl-FLIP, binds only to the active site and not the activator site, and thus exhibits standard hyperbolic kinetics. An analytical expression is presented that accounts for the kinetics of both substrates in the absence and presence of a nonsubstrate short-chain phospholipid, dihexanoylphosphatidylcholine. The fluorogenic substrates detect activation at much lower concentrations of dihexanoylphosphatidylcholine than previously reported.     

Immobilization of chloroperoxidase on silica-based materials for 4,6-dimethyl dibenzothiophene oxidation

Silica-based materials have been used as effective supports for the immobilization of enzymes. Moreover, the understanding on the oxidation of sulfur compounds by immobilized chloroperoxidase represents a step further in the development of a biocatalytic desulfurization process of fossil fuels. Here, chloroperoxidase from Caldariomyces fumago was immobilized on amorphous and structured silica-based materials either physically or covalently using an organosilane derivative for the oxidation of a recalcitrant organosulfur compound currently found in gas oil and diesel, such as 4,6-dimethyldibenzothiophene (4,6-DMDBT). Such materials were characterized by FTIR, N₂-adsorption, XRD, SEM and TEM. We have found that the chemical differences on the silanol/siloxane groups of SG/67 and SBA15 mesoporous materials deeply modify the enzymatic load, activity, thermal stability and reusability. The physical immobilization of CPO was characterized by a high adsorption capacity (q_m) and affinity constants (K_m) when compared to the covalent approach, but it resulted more sensitive to temperature than free, the silanized and covalently immobilized enzyme. The thermal residual activity as well as reusability of CPO were first improved by silanization, then by covalent immobilization in a support with a large pore size and high silanol/siloxane ratio.     

Mutation of Glu-166 Blocks the Substrate-Induced Dimerization of SARS Coronavirus Main Protease

The maturation of SARS coronavirus involves the autocleavage of polyproteins 1a and 1ab by the main protease (Mpro) and a papain-like protease; these represent attractive targets for the development of anti-SARS drugs. The functional unit of Mpro is a homodimer, and each subunit has a His-41?Cys-145 catalytic dyad. Current thinking in this area is that Mpro dimerization is essential for catalysis, although the influence of the substrate binding on the dimer formation has never been explored. Here, we delineate the contributions of the peptide substrate to Mpro dimerization. Enzyme kinetic assays indicate that the monomeric mutant R298A/L exhibits lower activity but in a cooperative manner. Analytical ultracentrifugation analyses indicate that in the presence of substrates, the major species of R298A/L shows a significant size shift toward the dimeric form and the monomer-dimer dissociation constant of R298A/L decreases by 12- to 17-fold, approaching that for wild-type. Furthermore, this substrate-induced dimerization was found to be reversible after substrates were removed. Based on the crystal structures, a key residue, Glu-166, which is responsible for recognizing the Gln-P1 of the substrate and binding to Ser-1 of another protomer, will interact with Asn-142 and block the S1 subsite entrance in the monomer. Our studies indicate that mutation of Glu-166 in the R298A mutant indeed blocks the substrate-induced dimerization. This demonstrates that Glu-166 plays a pivotal role in connecting the substrate binding site with the dimer interface. We conclude that protein-ligand and protein-protein interactions are closely correlated in Mpro.     

Atypical kinetic behavior of chloroperoxidase-mediated oxidative halogenation of polycyclic aromatic hydrocarbons

We have identified an atypical kinetic behavior for the oxidative halogenation of several polycyclic aromatic hydrocarbons (PAHs) by chloroperoxidase (CPO) from Caldariomyces fumago. This behavior resembles the capacity of some members of the P450 family to simultaneously recognize several substrate molecules at their active sites. Indeed, fluorometric studies showed that PAHs exist in solution as monomers and π-π dimers that interact to different extents with CPO. The dissociation constants of dimerization were evaluated for every single PAH by spectrofluorometry. Furthermore, docking studies also suggest that CPO might recognize either one or two substrate molecules in its active site. The atypical sigmoidal kinetic behavior of CPO in the oxidative halogenation of PAHs is explained in terms of different kinetic models for non-heteroatomic PAHs (naphthalene, anthracene and pyrene). The results suggest that the actual substrate for CPO in this study was the π-π dimer for all evaluated PAHs.