Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies.

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Structural determinants of specificity and regulation of activity in the allosteric loop network of human KLK8/neuropsin.

Human KLK8/neuropsin, a kallikrein-related serine peptidase, is mostly expressed in skin and the hippocampus regions of the brain, where it regulates memory formation by synaptic remodeling. Substrate profiles of recombinant KLK8 were analyzed with positional scanning using fluorogenic tetrapeptides and the proteomic PICS approach, which revealed the prime side specificity. Enzyme kinetics with optimized substrates showed stimulation by Ca2+ and inhibition by Zn2+, which are physiological regulators. Crystal structures of KLK8 with a ligand-free active site and with the inhibitor leupeptin explain the subsite specificity and display Ca2+ bound to the 75-loop. The variants D70K and H99A confirmed the antagonistic role of the cation binding sites. Molecular docking and dynamics calculations provided insights in substrate binding and the dual regulation of activity by Ca2+ and Zn2+, which are important in neuron and skin physiology. Both cations participate in the allosteric surface loop network present in related serine proteases. A comparison of the positional scanning data with substrates from brain suggests an adaptive recognition by KLK8, based on the tertiary structures of its targets. These combined findings provide a comprehensive picture of the molecular mechanisms underlying the enzyme activity of KLK8.

Kallikrein-related peptidase 5 and seasonal influenza viruses, limitations of the experimental models for activating proteases.

Every year, influenza A virus (IAV) affects and kills many people worldwide. The viral hemagglutinin (HA) is a critical actor in influenza virus infectivity which needs to be cleaved by host serine proteases to exert its activity. KLK5 has been identified as an activating protease in humans with a preference for the H3N2 IAV subtype. We investigated the origin of this preference using influenza A/Puerto Rico/8/34 (PR8, H1N1) and A/Scotland/20/74 (Scotland, H3N2) viruses. Pretreatment of noninfectious virions with human KLK5 increased infectivity of Scotland IAV in MDCK cells and triggered influenza pneumonia in mice. These effects were not observed with the PR8 IAV. Molecular modeling and in vitro enzymatic studies of peptide substrates and recombinant HAs revealed that the sequences around the cleavage site do not represent the sole determinant of the KLK5 preference for the H3N2 subtype. Using mouse Klk5 and Klk5-deficient mice, we demonstrated in vitro and in vivo that the mouse ortholog protease is not an IAV activating enzyme. This may be explained by unfavorable interactions between H3 HA and mKlk5. Our data highlight the limitations of some approaches used to identify IAV-activating proteases.

Specificity profiling of human trypsin-isoenzymes.

In humans, three different trypsin-isoenzymes have been described. Of these, trypsin-3 appears to be functionally different from the others. In order to systematically study the specificity of the trypsin-isoenzymes, we utilized proteome-derived peptide libraries and quantitative proteomics. We found similar specificity profiles dominated by the well-characterized preference for cleavage after lysine and arginine. Especially, trypsin-1 slightly favored lysine over arginine in this position, while trypsin-3 did not discriminate between them. In the P1' position, which is the residue C-terminal to the cleavage site, we noticed a subtle enrichment of alanine and glycine for all three trypsins and for trypsin-3 there were additional minor P1' and P2' preferences for threonine and aspartic acid, respectively. These findings were confirmed by FRET peptide substrates showing different susceptibility to cleavage by different trypsins. The preference of trypsin-3 for aspartic acid in P2' is explained by salt bridge formation with the unique Arg193. This salt bridge enables and stabilizes a canonical oxyanion conformation by the amides of Ser195 and Arg193, thus manifesting a selective substrate-assisted catalysis. As trypsin-3 has been proposed to be a therapeutic target and marker for cancers, our results may aid the development of specific inhibitors for cancer therapy and diagnostic probes.

Structural analyses of Arabidopsis thaliana legumain γ reveal differential recognition and processing of proteolysis and ligation substrates.

Legumain is a dual-function protease-peptide ligase whose activities are of great interest to researchers studying plant physiology and to biotechnological applications. However, the molecular mechanisms determining the specificities for proteolysis and ligation are unclear because structural information on the substrate recognition by a fully activated plant legumain is unavailable. Here, we present the X-ray structure of Arabidopsis thaliana legumain isoform γ (AtLEGγ) in complex with the covalent peptidic Ac-YVAD chloromethyl ketone (CMK) inhibitor targeting the catalytic cysteine. Mapping of the specificity pockets preceding the substrate-cleavage site explained the known substrate preference. The comparison of inhibited and free AtLEGγ structures disclosed a substrate-induced disorder-order transition with synergistic rearrangements in the substrate-recognition sites. Docking and in vitro studies with an AtLEGγ ligase substrate, sunflower trypsin inhibitor (SFTI), revealed a canonical, protease substrate-like binding to the active site-binding pockets preceding and following the cleavage site. We found the interaction of the second residue after the scissile bond, P2'-S2', to be critical for deciding on proteolysis versus cyclization. cis-trans-Isomerization of the cyclic peptide product triggered its release from the AtLEGγ active site and prevented inadvertent cleavage. The presented integrative mechanisms of proteolysis and ligation (transpeptidation) explain the interdependence of legumain and its preferred substrates and provide a rational framework for engineering optimized proteases, ligases, and substrates.

N-Terminomics identifies HtrA1 cleavage of thrombospondin-1 with generation of a proangiogenic fragment in the polarized retinal pigment epithelial cell model of age-related macular degeneration.

Age-related macular degeneration (AMD) is the leading cause of irreversible blindness in the elderly population. Variants in the HTRA1-ARMS2 locus have been linked to increased AMD risk. In the present study we investigated the impact of elevated HtrA1 levels on the retina pigment epithelial (RPE) secretome using a polarized culture system. Upregulation of HtrA1 alters the abundance of key proteins involved in angiogenesis and extracellular matrix remodeling. Thrombospondin-1, an angiogenesis modulator, was identified as a substrate for HtrA1 using terminal amine isotope labeling of substrates in conjunction with HtrA1 specificity profiling. HtrA1 cleavage of thrombospondin-1 was further corroborated by in vitro cleavage assays and targeted proteomics together with small molecule inhibition of HtrA1. While thrombospondin-1 is anti-angiogenic, the proteolytically released N-terminal fragment promotes the formation of tube-like structure by endothelial cells. Taken together, our findings suggest a mechanism by which increased levels of HtrA1 may contribute to AMD pathogenesis. The proteomic data has been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the data set identifier. For quantitative secretome analysis, project accession: PXD007691, username: reviewer45093@ebi.ac.uk, password: 1FUpS6Yq. For TAILS analysis, project accession: PXD007139, username: reviewer76731@ebi.ac.uk, password: sNbMp7xK.

Distinct Roles of Catalytic Cysteine and Histidine in the Protease and Ligase Mechanisms of Human Legumain As Revealed by DFT-Based QM/MM Simulations.

The cysteine protease enzyme legumain hydrolyzes peptide bonds with high specificity after asparagine and under more acidic conditions after aspartic acid [Baker E. N.J. Mol. Biol.1980, 141, 441-484; Baker E. N.; J. Mol. Biol.1977, 111, 207-210; Drenth J.; Biochemistry1976, 15, 3731-3738; Menard R.; J. Cell. Biochem.1994, 137; Polgar L.Eur. J. Biochem.1978, 88, 513-521; Storer A. C.; Methods Enzymol.1994, 244, 486-500. Remarkably, legumain additionally exhibits ligase activity that prevails at pH > 5.5. The atomic reaction mechanisms including their pH dependence are only partly understood. Here we present a density functional theory (DFT)-based quantum mechanics/molecular mechanics (QM/MM) study of the detailed reaction mechanism of both activities for human legumain in solution. Contrasting the situation in other papain-like proteases, our calculations reveal that the active site Cys189 must be present in the protonated state for a productive nucleophilic attack and simultaneous rupture of the scissile peptide bond, consistent with the experimental pH profile of legumain-catalyzed cleavages. The resulting thioester intermediate (INT1) is converted by water attack on the thioester into a second intermediate, a diol (INT2), which is released by proton abstraction by Cys189. Surprisingly, we found that ligation is not the exact reverse of the proteolysis but can proceed via two distinct routes. Whereas the transpeptidation route involves aminolysis of the thioester (INT1), at pH 6 a cysteine-independent, histidine-assisted ligation route was found. Given legumain's important roles in immunity, cancer, and neurodegenerative diseases, our findings open up possibilities for targeted drug design in these fields.

Discovery of a Potent Inhibitor Class with High Selectivity toward Clostridial Collagenases.

Secreted virulence factors like bacterial collagenases are conceptually attractive targets for fighting microbial infections. However, previous attempts to develop potent compounds against these metalloproteases failed to achieve selectivity against human matrix metalloproteinases (MMPs). Using a surface plasmon resonance-based screening complemented with enzyme inhibition assays, we discovered an N-aryl mercaptoacetamide-based inhibitor scaffold that showed sub-micromolar affinities toward collagenase H (ColH) from the human pathogen Clostridium histolyticum. Moreover, these inhibitors also efficiently blocked the homologous bacterial collagenases, ColG from C. histolyticum, ColT from C. tetani, and ColQ1 from the Bacillus cereus strain Q1, while showing negligible activity toward human MMPs-1, -2, -3, -7, -8, and -14. The most active compound displayed a more than 1000-fold selectivity over human MMPs. This selectivity can be rationalized by the crystal structure of ColH with this compound, revealing a distinct non-primed binding mode to the active site. The non-primed binding mode presented here paves the way for the development of selective broad-spectrum bacterial collagenase inhibitors with potential therapeutic application in humans.

Endolysosomal Degradation of Allergenic Ole e 1-Like Proteins: Analysis of Proteolytic Cleavage Sites Revealing T Cell Epitope-Containing Peptides.

Knowledge of the susceptibility of proteins to endolysosomal proteases provides valuable information on immunogenicity. Though Ole e 1-like proteins are considered relevant allergens, little is known about their immunogenic properties and T cell epitopes. Thus, six representative molecules, i.e., Ole e 1, Fra e 1, Sal k 5, Che a 1, Phl p 11 and Pla l 1, were investigated. Endolysosomal degradation and peptide generation were simulated using microsomal fractions of JAWS II dendritic cells. Kinetics and peptide patterns were evaluated by gel electrophoresis and mass spectrometry. In silico MHC (major histocompatibility complex) class II binding prediction was performed with ProPred. Cleavage sites were assigned to the primary and secondary structure, and in silico docking experiments between the protease cathepsin S and Ole e 1 were performed. Different kinetics during endolysosomal degradation were observed while similar peptide profiles especially at the C-termini were detected. Typically, the identified peptide clusters comprised the previously-reported T cell epitopes of Ole e 1, consistent with an in silico analysis of the T cell epitopes. The results emphasize the importance of the fold on allergen processing, as also reflected by conserved cleavage sites located within the large flexible loop. In silico docking and mass spectrometry results suggest that one of the first Ole e 1 cleavages might occur at positions 107-108. Our results provided kinetic and structural information on endolysosomal processing of Ole e 1-like proteins.

Targeting of a Helix-Loop-Helix Transcriptional Regulator by a Short Helical Peptide.

The Id proteins (Id1-4) are cell-cycle regulators that play a key role during development, in cancer and vascular disorders. They contain a conserved helix-loop-helix (HLH) domain that folds into a parallel four-helix bundle upon self- or hetero-association with basic-HLH transcription factors. By using such protein-protein interactions, the Id proteins inhibit cell differentiation and promote cell-cycle progression. Accordingly, their supporting role in cancer has been convincingly demonstrated, which makes these proteins interesting therapeutic targets. Herein we present a short peptide containing an (i,i+4)-lactam bridge and a hydrophobic (Φ) three-residue motif Φ(i)-Φ(i+3)-Φ(i+6), which adopts a helical conformation in water, shows Id protein binding in the low-micromolar range, penetrates into breast (MCF-7 and T47D) and bladder (T24) cancer cells, accumulates in the nucleus, and decreases cell viability to ∼50 %. Thus, this cyclopeptide is a promising scaffold for the development of Id protein binders that impair cancer cell viability.

Impact of the amino acid sequence on the conformation of side chain lactam-bridged octapeptides.

Synthetic helical peptides are valuable scaffolds for the development of modulators of protein-protein interactions involving helical motifs. Backbone-to-side chain or side chain-to-side chain constraints have been and still are intensively exploited to stabilize short α-helices. Very often, these constraints have been combined with backbone modifications induced by Cα-tetrasubstituted, ß-, or γ-amino acids, which facilitate the α-peptide or α/ß/γ-peptide adopting an α-helical conformation. In this work, we investigated the helical character of octapeptides that were cyclized by a Lys-Asp-(i,i + 4)-lactam bridge. We started with two sequences extracted from the helix-loop-helix region of the Id proteins, which are inhibitors of cell differentiation during development and in cancer. Nineteen analogs containing the lactam bridge at different positions and displaying different amino acid core triads (i + 1,2,3) as well as outer residues were prepared by solid-phase methodology. Their conformation in water and water/2,2,2-trifluoroethanol mixtures was investigated by circular dichroism (CD) spectroscopy. The cyclopeptides could be grouped in helix-prone and non-helix-prone structures. Both the amino acid core triad (i + 1,2,3) and the pendant residues positively or negatively affected the formation of a helical structure. Computational studies based on the NMR-derived helical structure of a cyclopeptide containing Aib at position (i + 2) of the triad were generally in agreement with the secondary structure propensity of the cyclopeptides observed by CD spectroscopy. In conclusion, the Lys-Asp-(i,i + 4)-lactam bridge may succeed or fail in the stabilization of short helices, depending on the primary structure. Moreover, computational methods may be valuable tools to discriminate helix-prone from non-helix-prone peptide-based macrolactams. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

QM/MM simulation (B3LYP) of the RNase A cleavage-transesterification reaction supports a triester A(N) + D(N) associative mechanism with an O2' H internal proton transfer.

The mechanism of the backbone cleavage-transesterification step of the RNase A enzyme remains controversial even after 60 years of study. We report quantum mechanics/molecule mechanics (QM/MM) free energy calculations for two optimized reaction paths based on an analysis of all structural data and identified by a search for reaction coordinates using a reliable quantum chemistry method (B3LYP), equilibrated structural optimizations, and free energy estimations. Both paths are initiated by nucleophilic attack of the ribose O2' oxygen on the neighboring diester phosphate bond, and both reach the same product state (PS) (a O3'-O2' cyclic phosphate and a O5' hydroxyl terminated fragment). Path 1, resembles the widely accepted dianionic transition-state (TS) general acid (His119)/base (His12) classical mechanism. However, this path has a barrier (25 kcal/mol) higher than that of the rate-limiting hydrolysis step and a very loose TS. In Path 2, the proton initially coordinating the O2' migrates to the nonbridging O1P in the initial reaction path rather than directly to the general base resulting in a triester (substrate as base) AN + DN mechanism with a monoanionic weakly stable intermediate. The structures in the transition region are associative with low barriers (TS1 10, TS2 7.5 kcal/mol). The Path 2 mechanism is consistent with the many results from enzyme and buffer catalyzed and uncatalyzed analog reactions and leads to a PS consistent with the reactive state for the following hydrolysis step. The differences between the consistently estimated barriers in Path 1 and 2 lead to a 10(11) difference in rate strongly supporting the less accepted triester mechanism.

Theoretical investigation of the enzymatic phosphoryl transfer of ß-phosphoglucomutase: revisiting both steps of the catalytic cycle.

Enzyme catalyzed phosphate transfer is a part of almost all metabolic processes. Such reactions are of central importance for the energy balance in all organisms and play important roles in cellular control at all levels. Mutases transfer a phosphoryl group while nucleases cleave the phosphodiester linkages between two nucleotides. The subject of our present study is the Lactococcus lactis ß-phosphoglucomutase (ß-PGM), which effectively catalyzes the interconversion of ß-D-glucose-1-phosphate (ß-G1P) to ß-D-glucose-6-phosphate (ß-G6P) and vice versa via stabile intermediate ß-D-glucose-1,6-(bis)phosphate (ß-G1,6diP) in the presence of Mg(2+). In this paper we revisited the reaction mechanism of the phosphoryl transfer starting from the bisphosphate ß-G1,6diP in both directions (toward ß-G1P and ß-G6P) combining docking techniques and QM/MM theoretical method at the DFT/PBE0 level of theory. In addition we performed NEB (nudged elastic band) and free energy calculations to optimize the path and to identify the transition states and the energies involved in the catalytic cycle. Our calculations reveal that both steps proceed via dissociative pentacoordinated phosphorane, which is not a stabile intermediate but rather a transition state. In addition to the Mg(2+) ion, Ser114 and Lys145 also play important roles in stabilizing the large negative charge on the phosphate through strong coordination with the phosphate oxygens and guiding the phosphate group throughout the catalytic process. The calculated energy barrier of the reaction for the ß-G1P to ß-G1,6diP step is only slightly higher than for the ß-G1,6diP to ß-G6P step (16.10 kcal mol(-1) versus 15.10 kcal mol(-1)) and is in excellent agreement with experimental findings (14.65 kcal mol(-1)).

Nucleotide docking: prediction of reactant state complexes for ribonuclease enzymes.

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools-MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0-for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.

Atomistic details of the associative phosphodiester cleavage in human ribonuclease H.

During translation of the genetic information of DNA into proteins, mRNA is synthesized by RNA polymerase and after the transcription process degraded by RNase H. The endoribonuclease RNase H is a member of the nucleotidyl-transferase (NT) superfamily and is known to hydrolyze the phosphodiester bonds of RNA which is hybridized to DNA. Retroviral RNase H is part of the viral reverse transcriptase enzyme that is indispensable for the proliferation of retroviruses, such as HIV. Inhibitors of this enzyme could therefore provide new drugs against diseases like AIDS. In our study we investigated the molecular mechanism of RNA cleavage by human RNase H using a comprehensive high level DFT/B3LYP QM/MM theoretical method for the calculation of the stationary points and nudged elastic band (NEB) and free energy calculations to identify the transition state structures, the rate limiting step and the reaction barrier. Our calculations reveal that the catalytic mechanism proceeds in two steps and that the nature of the nucleophile is a water molecule. In the first step, the water attack on the scissile phosphorous is followed by a proton transfer from the water to the O2P oxygen and a trigonal bipyramidal pentacoordinated phosphorane is formed. Subsequently, in the second step the proton is shuttled to the O3' oxygen to generate the product state. During the reaction mechanism two Mg(2+) ions support the formation of a stable associated in-line S(N)2-type phosphorane intermediate. Our calculated energy barrier of 19.3 kcal mol(-1) is in excellent agreement with experimental findings (20.5 kcal mol(-1)). These results may contribute to the clarification and understanding of the RNase H reaction mechanism and of further enzymes from the RNase family.

A dianionic phosphorane intermediate and transition states in an associative A(N)+D(N) mechanism for the ribonucleaseA hydrolysis reaction.

The RNaseA enzyme efficiently cleaves phosphodiester bonds in the RNA backbone. Phosphoryl transfer plays a central role in many biochemical reactions, and this is one of the most studied enzymes. However, there remains considerable controversy about the reaction mechanism. Most of this debate centers around the roles of the conserved residues, structures of the transition state or states, the possibility of a stable intermediate, and the charge and structure of this intermediate. In this communication we report calculations of the mechanism of the hydrolysis step in this reaction using a comprehensive QM/MM theoretical approach that includes a high level calculation of the interactions in the QM region, free energy estimates along an NEB optimized reaction path, and the inclusion of the interaction of the protein surroundings and solvent. Contrary to prior calculations we find a stable pentacoordinated dianionic phosphorane intermediate in the reaction path supporting an A(N)+D(N) reaction mechanism. In the transition state in the path from the reactant to the intermediate state (with barrier of 3.96 kcal/mol and intermediate stability of 2.21 kcal/mol) a proton from the attacking water is partially transferred to the His119 residue and the PO bond only partially formed from the remaining nucleophilic OH(-) species (bond order (BO) 0.11). In passing from the intermediate to the product state (barrier 13.22 kcal/mol) the PO bond on the cyclic phosphorane intermediate is nearly broken (BO 0.28) and the transfer of the proton from the Lys41 is almost complete (Lys41-H BO 0.87). In the product state a proton has been transferred from Lys41 to the O2' position of the sugar. The role of Lys41 as the catalytic acid is a result of the relative positioning of the Lys41 and His12 in the catalytic site. This configuration is supported by calculations and docking studies.

Revision of the absolute configuration of plumericin and isoplumericin from Plumeria rubra.

The absolute configurations of plumericin (1) and isoplumericin (2), isolated from Plumeria rubra, were re-assigned based on a combination of X-ray crystal-structure determination and quantum-mechanical calculations of their circular dichroism (CD) spectra. The experimental CD spectra showed an excellent match with those calculated for the (1S,5R,8R,9R,10R) absolute configuration (corresponding to ent-1 and ent-2, resp.), opposite to that generally accepted and published in the literature. Since the (false) plumericin configuration has been often used to derive the absolute configuration of related iridoids by chemical correlation, their absolute configurations also have to be reconsidered.

Synthesis and structure-activity relationship of antifungal coniothyriomycin analogues.

The structure of the antifungal metabolite coniothyriomycin was systematically modified by changing the acids of the open chain imide, modification of the hydrophobicity, variation in the degree of saturation, replacement of carbons by nitrogen or oxygen, and incorporation of the open chain molecule into cyclic arrangements. Structure-activity studies showed that antifungal activity was retained by replacement of phenylacetic acids by benzoic acids in the imide structure but diminished by hydrogenation of the fumaric ester part.

Structure-activity relationships in allergic contact dermatitis. Part III. The sensitizing capacity of substituted phenanthrenequinones: a quantum-mechanical approach.

BACKGROUND: Nonterpenoid and diterpenoid phenanthrenequinones (PACs) have been found in the plant kingdom. Some of them occur in plants used in traditional Chinese medicine like Tan-Shen whereas others are constituents of orchids that are popular as ornamental plants. OBJECTIVE: Case reports and our own observations in orchid nurseries suggest that some or even all of these PACs possess a distinct sensitizing potency. Occasional exposure (particularly of botanists) to field-grown orchids, as well as occupational contact with sawdust of PAC-containing tropical timbers, caused allergic contact dermatitis. However, experimental studies in guinea pigs to determine the sensitizing capacity of PACs have not been performed so far. METHODS: Guinea pigs were sensitizied by a modified Freund's complete adjuvant method with four naturally occurring and 22 synthetic PACs in order to find out which and how many substituents at the carbons of the three rings of the PAC will influence the sensitizing power of the molecule. Subsequently, the lowest unoccupied molecular orbital (LUMO) coefficients were calculated to show whether a correlation exists between chemical reactivity and sensitizing capacity. RESULTS: Sensitizing capacity was found to be strong in two PACs, moderate in eight PACs, and weak in ten PACs. Five PACs were extremely weak in sensitizing capacity, and one PAC was completely negative. Two substituents on the left-hand carbons C-7 and C-8 of ring C were shown to be responsible for a strong sensitizing capacity. One methoxy group alone or three of them, especially when localized at C-5, decreased the sensitizing capacity to moderate. Substitution with a methoxy group at C-3 and/or at C-2 of the quinonoid ring itself (ring A) led to a weak sensitizing capacity. The ortho-quinones 1,2-PAC and 9,10-PAC were also weakly sensitizing. In fact, LUMO coefficient calculations corroborated a good correlation between chemical reactivity and sensitizing capacity. CONCLUSION: Substitution with methoxy groups at C-7 and/or at C-8 of ring C of 1,4-phenanthrenequinone increases the LUMO coefficients at the 2,3 double bond of ring A and thus facilitates nucleophilic substitution of protein nitrogen or sulfur nucleophiles at this electron-deficient double bond. The four naturally occurring PACs that were investigated--cypripedin, denbinobin, annoquinone-A, and latinone--do not fulfill these criteria and are thus only weak sensitizers. However, as-yet-unstudied phenanthrenequinones occurring in plants or trees and having no substituents at C-2 or C-3 of the quinonoid ring must be considered potentially strong allergens.