New insights about HERG blockade obtained from protein modeling, potential energy mapping, and docking studies

We created a homology model of the homo-tetrameric pore domain of HERG using the crystal structure of the bacterial potassium channel, KvAP, as a template. We docked a set of known blockers with well-characterized effects on channel function into the lumen of the pore between the selectivity filter and extracellular entrance using a novel docking and refinement procedure incorporating Glide and Prime. Key aromatic groups of the blockers are predicted to form multiple simultaneous ring stacking and hydrophobic interactions among the eight aromatic residues lining the pore. Furthermore, each blocker can achieve these interactions via multiple docking configurations. To further interpret the docking results, we mapped hydrophobic and hydrophilic potentials within the lumen of each refined docked complex. Hydrophilic iso-potential contours define a 'propeller-shaped' volume at the selectivity filter entrance. Hydrophobic contours define a hollow 'crown-shaped' volume located above the 'propeller', whose hydrophobic 'rim' extends along the pore axis between Tyr652 and Phe656. Blockers adopt conformations/binding orientations that closely mimic the shapes and properties of these contours. Blocker basic groups are localized in the hydrophilic 'propeller', forming electrostatic interactions with Ser624 rather than a generally accepted pi-cation interaction with Tyr652. Terfenadine, cisapride, sertindole, ibutilide, and clofilium adopt similar docked poses, in which their N-substituents bridge radially across the hollow interior of the 'crown' (analogous to the hub and spokes of a wheel), and project aromatic/hydrophobic portions into the hydrophobic 'rim'. MK-499 docks with its longitudinal axis parallel to the axis of the pore and 'crown', and its hydrophobic groups buried within the hydrophobic 'rim'.     

Hydroxyethylene isostere inhibitors of human immunodeficiency virus-1 protease: structure-activity analysis using enzyme kinetics, X-ray crystallography, and infected T-cell assays.

Analogues of peptides ranging in size from three to six amino acids and containing the hydroxyethylene dipeptide isosteres Phe psi Gly, Phe psi Ala, Phe psi NorVal, Phe psi Leu, and Phe psi Phe, where psi denotes replacement of CONH by (S)-CH(OH)CH2, were synthesized and studied as HIV-1 protease inhibitors. Inhibition constants (Ki) with purified HIV-1 protease depend strongly on the isostere in the order Phe psi Gly greater than Phe psi Ala greater than Phe psi NorVal greater than Phe psi Leu greater than Phe psi Phe and decrease with increasing length of the peptide analogue, converging to a value of 0.4 nM. Ki values are progressively less dependent on inhibitor length as the size of the P1' side chain within the isostere increases. The structures of HIV-1 protease complexed with the inhibitors Ala-Ala-X-Val-Val-OMe, where X is Phe psi Gly, Phe psi Ala, Phe psi NorVal, and Phe psi Phe, have been determined by X-ray crystallography (resolution 2.3-3.2 A). The crystals exhibit symmetry consistent with space group P6(1) with strong noncrystallographic 2-fold symmetry, and the inhibitors all exhibit 2-fold disorder. The inhibitors bind in similar conformations, forming conserved hydrogen bonds with the enzyme. The Phe psi Gly inhibitor adopts an altered conformation that places its P3' valine side chain partially in the hydrophobic S1' pocket, thus suggesting an explanation for the greater dependence of the Ki value on inhibitor length in the Phe psi Gly series. From the kinetic and crystallographic data, a minimal inhibitor model for tight-binding inhibition is derived in which the enzyme subsites S2-S2' are optimally occupied. The Ki values for several compounds are compared with their potencies as inhibitors of proteolytic processing in T-cell cultures chronically infected with HIV-1 (MIC values) and as inhibitors of acute infectivity (IC50 values). There is a rank-order correspondence, but a 20-1000-fold difference, between the values of Ki and those of MIC or IC50. IC50 values can approach those of Ki but are highly dependent on the conditions of the acute infectivity assay and are influenced by physiochemical properties of the inhibitors such as solubility.     

Design of inhibitors for human tissue kallikrein using non-natural aromatic and basic amino acids

We explored the unique substrate specificity of the primary S, subsite of human urinary kallikrein (hK1), which accepts both Phe or Arg synthesizing and assaying peptides derived from Phenylacetyl-Phe-Ser-Arg-EDDnp, a previously described inhibitor with analgesic and anti-inflammatory activities [Emim et al., Br. J. Pharmacol. 130 (2000), 1099-1107]. Phe was substituted by amino acids containing larger aliphatic or aromatic side chains as well as by non-natural basic amino acids, which were designed to combine a large hydrophobic and/or aromatic group with a positively-charged group at their side chains. In general, all peptides with basic amino acids represented better inhibitors than those with hydrophobic amino acids. Furthermore, the S1 subsite specificity proved to be much more selective than the mere distinction between Phe and Arg, for minor differences in the side chains of the non-natural amino acids resulted in major differences in the Ki values. Finally, we present a series of peptides that were assayed as competitive inhibitors for human tissue kallikrein that may lead to the development of novel peptides, which are both more potent and selective.     

Carbonic anhydrase inhibitors: Valdecoxib binds to a different active site region of the human isoform II as compared to the structurally related cyclooxygenase II "selective" inhibitor celecoxib

The high resolution X-ray crystal structure of the adduct of human carbonic anhydrase (CA, EC 4.2.1.1) isoform II (hCA II) with the clinically used painkiller valdecoxib, acting as a potent CA II and cyclooxygenase-2 (COX-2) inhibitor, is reported. The ionized sulfonamide moiety of valdecoxib is coordinated to the catalytic Zn(II) ion with a tetrahedral geometry. The phenyl-isoxazole moiety of the inhibitor fills the active site channel and interacts with the side chains of Gln92, Val121, Leu198, Thr200, and Pro202. Its 3-phenyl group is located into a hydrophobic pocket, simultaneously establishing van der Waals interactions with the aliphatic side chain of various hydrophobic residues (Val135, Ile91, Val121, Leu198, and Leu141) and a strong offset face-to-face stacking interaction with the aromatic ring of Phe131 (the chi1 angle of which is rotated about 90 degrees with respect to what was observed in the structure of the native enzyme and those of other sulfonamide complexes). Celecoxib, a structurally related COX-2 inhibitor for which the X-ray crystal structure was reported earlier, binds in a completely different manner to hCA II as compared to valdecoxib. Celecoxib completely fills the entire CA II active site, with its trifluoromethyl group in the hydrophobic part of the active site and the p-tolyl moiety in the hydrophilic one, not establishing any interaction with Phe131. In contrast to celecoxib, valdecoxib was rotated about 90 degrees around the chemical bond connecting the benzensulfonamide and the substituted isoxazole ring allowing for these multiple favorable interactions. These different binding modes allow for the further drug design of various CA inhibitors belonging to the benzenesulfonamide class.     

Peripheral but crucial: A hydrophobic pocket (Tyr706, Leu337, and Met336) for potent and selective inhibition of neuronal nitric oxide synthase

Selective inhibition of the neuronal isoform of nitric oxide synthase (nNOS) over endothelial nitric oxide synthase (eNOS) and inducible nitric oxide synthase (iNOS) has become a promising strategy for the discovery of new therapeutic agents for neurodegenerative diseases. However, because of the high sequence homology of different isozymes in the substrate binding pocket, developing inhibitors with both potency and excellent isoform selectivity remains a challenging problem. Herein, we report the evaluation of a recently discovered peripheral hydrophobic pocket (Tyr⁷⁰⁶, Leu³³⁷, and Met³³⁶) that opens up upon inhibitor binding and its potential in designing potent and selective nNOS inhibitors using three compounds, 2a, 2b, and 3. Crystal structure results show that inhibitors 2a and 3 adopted the same binding mode as lead compound 1. We also found that hydrophobic interactions between the 4-methyl group of the aminopyridine ring of these compounds with the side chain of Met³³⁶, as well as the π–π stacking interaction between the pyridinyl motif and the side chain of Tyr⁷⁰⁶ are important for the high potency and selectivity of these nNOS inhibitors.     

Glutamic acid analogues as potent dipeptidyl peptidase IV and 8 inhibitors

To find potent and selective inhibitors of dipeptidyl peptidase IV (DPP-IV), we synthesized a series of 2-cyanopyrrolidine with P2-site 4-substituted glutamic acid derivatives and tested their activities against DPP-IV, DPP8, and DPP-II. Analogues that incorporated a bulky substituent at the first carbon position of benzylamine or isoquinoline showed over 30-fold selectivity for DPP-IV over both DPP8 and DPP-II. From structure-activity relationship studies, we speculate that the S2 site of DPP8 might be similar to that of DPP-IV, while DPP-IV inhibitor with N-substituted glycine in the P2 site and/or with a moiety involving in hydrophobic interaction with the side chain of Phe357 might provide a better selectivity for DPP-IV over DPP8.     

Mutational analysis of block and facilitation of HERG current by a class III anti-arrhythmic agent, nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

Mutational Analysis of Block and Facilitation of HERG Current by A Class III Anti-Arrhythmic Agent,Nifekalant

Chemicals and toxins are useful tools to elucidate the structure-function relationship of various proteins including ion channels. The HERG channel is blocked by many compounds and this may cause life-threatening cardiac arrhythmia. Besides block, some chemicals such as the class III anti-arrhythmic agent nifekalant stimulate HERG at low potentials by shifting its activation curve towards hyperpolarizing voltages. This is called "facilitation". Here, we report mutations and simulations analyzing the association between nifekalant and channel pore residues for block and facilitation. Alanine-scanning mutagenesis was performed in the pore region of HERG. The mutations at the base of the pore helix (T623A), the selectivity filter (V625A) and the S6 helix (G648A, Y652A and F656A) abolished and S624A attenuated both block and facilitation induced by the drug. On the other hand, the mutation of other residues caused either an increase or a decrease in nifekalant-induced facilitation without affecting block. An open-state homology model of the HERG pore suggested that T623, S624, Y652 and F656 faced the central cavity, and were positioned within geometrical range for the drug to be able to interact with all of them at the same time. Of these, S649 was the only polar residue located within possible interaction distance from the drug held in its blocking position. Further mutations and flexible-docking simulations suggest that the size, but not the polarity, of the side chain at S649 is critical for drug induced facilitation.     

NMR studies of interactions between inhibitors and porcine pancreatic phospholipase A2.

Two-dimensional NMR studies were performed on the complexes of porcine pancreatic phospholipase A2, bound to a micellar lipid-water interface of fully deuterated dodecylphosphocholine, with competitive inhibitors derived from the following general structure: [formula: see text] X and Y are alkyl chains with various 'reporter groups'. The interactions between the inhibitor and the enzyme were localized by comparison of 2-D nuclear Overhauser effect spectra using protonated and selectively deuterated inhibitors, and inhibitors with groups having easily identifiable chemical shifts. These experiments led us to the following conclusions for the phospholipase A2/inhibitor/micelle complex: i) the His48 C2 ring proton is in close proximity to both the amide proton and the methylene protons at the sn-1 position of the glycerol skeleton of the inhibitor, ii) the acyl chain of the inhibitor at the sn-2 position makes hydrophobic contacts near Phe5, Ile9, Phe22 and Phe106; iii) no interactions between the acyl chain at the sn-1 position and the protein could be identified. Comparison of our results on the enzyme/inhibitor/micelle ternary complex with the crystal structure of the enzyme-inhibitor complex shows that the mode of inhibitor binding is similar. However, in several cases we found indications that the hydrophobic chains of the inhibitors can have multiple conformations.     

Investigation of the nature of the methionine-pi interaction in beta-hairpin peptide model systems

There are frequent contacts between aromatic rings and sulfur atoms in proteins. However, it is unclear to what degree this putative interaction is stabilizing and what the nature of the interaction is. We have investigated the aryl-sulfur interaction by placing a methionine residue diagonal to an aromatic ring on the same face of a beta-hairpin, which places the methionine side chain in close proximity to the aryl side chain. The methionine (Met)-aryl interaction was compared with an equivalent hydrophobic and cation-pi interaction in the context of the beta-hairpin. The interaction between phenylalanine (Phe), tryptophan (Trp), or cyclohexylalanine (Cha) and Met stabilized the beta-hairpin by -0.3 to -0.5 kcal mole(-1), as determined by double-mutant cycles. The peptides were subjected to thermal denaturations that suggest a hydrophobic driving force for the interactions between Met and Trp or Cha. The observed interaction of Met or norleucine (Nle) with Trp or Cha are quite similar, implying a hydrophobic driving force for the Met-pi interaction. However, the thermodynamic data suggest that there may be some differences between the interaction of Met with Trp and Phe and that there may be a small thermodynamic component to the Met...Phe interaction.     

1,4-Dihydroindeno[1,2-c]pyrazoles as novel multitargeted receptor tyrosine kinase inhibitors

A series of 1,4-dihydroindeno[1,2-c]pyrazoles with a 3-thiophene substituent carrying a urea-type side chain were identified as potent multitargeted (VEGFR and PDGFR families) receptor tyrosine kinase inhibitors. A KDR homology model suggested that the urea moiety is able to interact with a recognition motif in the hydrophobic specificity pocket of the enzyme.     

Thermodynamics of interactions between amino acid side chains: experimental differentiation of aromatic-aromatic, aromatic-aliphatic, and aliphatic-aliphatic side-chain interactions in water

A stationary phase for high-pressure liquid chromatography has been prepared by derivatizing microparticulate silica gel with functionality mimicking the side chain of isoleucine. The chromatographic retentions of a series of hydrophobic and amphiphilic amino acid analytes on this stationary phase (Ile MSP) using an aqueous mobile phase were measured as a function of temperature from 273 K to 323 K. Observed temperature dependencies are consistent with a constant change in heat capacity, DeltaC degrees P, upon binding of the analyte to the stationary phase. The curvatures of plots of retention data versus temperature (related to the magnitude of DeltaC degrees P) are distinctly different for retention of aromatic and aliphatic analytes, with retention of aliphatic analytes Val, Ile, and Leu exhibiting the characteristic signature of the hydrophobic effect, i.e., a large negative DeltaC degrees P upon desolvation from water and a maximum of retention around room temperature. Retention of aromatic analytes (Trp, Phe, and Tyr) involves smaller heat capacity changes and pronounced negative enthalpies of interaction with the stationary phase. Estimates of DeltaC degrees P for the interactions of analyte side chains with the Ile side chain were obtained by fitting the temperature dependence of retention to an expression derived from thermodynamic considerations and chromatographic theory. Similar estimates were made for interactions with the Phe side chain, using previously published data for a phenylalanine mimic stationary phase (Phe MSP) (. Protein Sci. 1:786-795). As with the Ile MSP, the retentions of aliphatic analytes show temperature dependencies markedly different from those of aromatic analytes. Data from both phases indicate that a realistic differentiation can be made between the interactions of various types of amino acid side chains tested (i.e., aliphatic/aliphatic, aliphatic/aromatic, and aromatic/aromatic) by comparison of the corresponding thermodynamic functions for pairwise interactions. The retention of leucine on the Phe MSP and that of phenylalanine on the Ile MSP showed similar DeltaC degrees P values, suggesting that the aromatic-aliphatic interaction is reasonably independent of the residue attached to the stationary phase. This result is consistent with a one-to-one interaction and suggests a simple way to estimate the column-dependent phase factor, making it possible to compare entropies and free energies of interaction obtained using different MSPs. The possibilities for using MSP-derived interaction potentials in folding simulations are discussed.     

In search of AKT kinase inhibitors as anticancer agents: structure-based design,docking, and molecular dynamics studies of 2,4,6-trisubstituted pyridines

The AKT isoforms are a group of key kinases that play a critical role in tumorigenesis. These enzymes are overexpressed in different types of cancers, such as breast, colon, prostate, ovarian, and lung. Because of its relevance the AKT isoforms are attractive targets for the design of anticancer molecules. However, it has been found that AKT1 and AKT3 isoforms have a main role in tumor progression and metastasis; thus, the identification of AKT isoforms specific inhibitors seems to be a challenge. Previously, we identified an ATP binding pocket pan-AKT inhibitor, this compound is a 2,4,6-trisubstituted pyridine (compound 11), which represents a new interesting scaffold for the developing of AKT inhibitors. Starting from the 2,4,6-trisubstituted pyridine scaffold, and guided by structure-based design technique, 42 new inhibitors were designed and further evaluated in the three AKT isoforms by multiple docking approach and molecular dynamics. Results showed that seven compounds presented binding selectivity for AKT1 and AKT3, better than for AKT2. The binding affinities of these seven compounds on AKT1 and AKT3 isoforms were mainly determined by hydrophobic contributions between the aromatic portion at position 4 of the pyridine ring with residues Phe236/234, Phe237/235, Phe438/435, and Phe442/439 in the ATP binding pocket. Results presented in this work provide an addition knowledge leading to promising selective AKT inhibitors.     

Solution structure of the catalytic domain of human stromelysin complexed with a hydrophobic inhibitor.

Stromelysin, a representative matrix metalloproteinase and target of drug development efforts, plays a prominent role in the pathological proteolysis associated with arthritis and secondarily in that of cancer metastasis and invasion. To provide a structural template to aid the development of therapeutic inhibitors, we have determined a medium-resolution structure of a 20-kDa complex of human stromelysin's catalytic domain with a hydrophobic peptidic inhibitor using multinuclear, multidimensional NMR spectroscopy. This domain of this zinc hydrolase contains a mixed beta-sheet comprising one antiparallel strand and four parallel strands, three helices, and a methionine-containing turn near the catalytic center. The ensemble of 20 structures was calculated using, on average, 8 interresidue NOE restraints per residue for the 166-residue protein fragment complexed with a 4-residue substrate analogue. The mean RMS deviation (RMSD) to the average structure for backbone heavy atoms is 0.91 A and for all heavy atoms is 1.42 A. The structure has good stereochemical properties, including its backbone torsion angles. The beta-sheet and alpha-helices of the catalytic domains of human stromelysin (NMR model) and human fibroblast collagenase (X-ray crystallographic model of Lovejoy B et al., 1994b, Biochemistry 33:8207-8217) superimpose well, having a pairwise RMSD for backbone heavy atoms of 2.28 A when three loop segments are disregarded. The hydroxamate-substituted inhibitor binds across the hydrophobic active site of stromelysin in an extended conformation. The first hydrophobic side chain is deeply buried in the principal S'1 subsite, the second hydrophobic side chain is located on the opposite side of the inhibitor backbone in the hydrophobic S'2 surface subsite, and a third hydrophobic side chain (P'3) lies at the surface.     

New insights into DHFR interactions: analysis of Pneumocystis carinii and mouse DHFR complexes with NADPH and two highly potent 5-(omega-carboxy(alkyloxy) trimethoprim derivatives reveals conformational correlations with activity and novel parallel ring stacking interactions

Structural data are reported for two highly potent antifolates, 2,4-diamino-5-[3',4'-dimethoxy-5'-(5-carboxy-1-pentynyl)]benzylpyrimidine (PY1011), with 5000-fold selectivity for Pneumocystis carinii dihydrofolate reductase (pcDHFR), relative to rat liver DHFR, and 2,4-diamino-5-[2-methoxy-5-(4-carboxybutyloxy)benzyl]pyrimidine (PY957), that has 80-fold selectivity for pcDHFR. Crystal structures are reported for NADPH ternary complexes with PY957 and pcDHFR, refined to 2.2 A resolution; with PY1011 and pcDHFR, refined to 2.0 A resolution; and with PY1011 and mouse DHFR (mDHFR), refined to 2.2 A resolution. These results reveal that the carboxylate of the omega-carboxyalkyloxy side chain of these inhibitors form ionic interactions with the conserved Arg in the substrate binding pocket of DHFR. These data suggest that the enhanced inhibitory activity of PY1011 compared with PY957 is, in part, due to the favorable contacts with Phe69 of pcDHFR by the methylene carbons of the inhibitor side chain that are oriented by the triple bond of the 1-pentynyl side chain. These contacts are not present in the PY957 pcDHFR complex, or in the PY1011 mDHFR complex. In the structure of mDHFR the site of Phe69 in pcDHFR is occupied by Asn64. These data also revealed a preference for an unusual parallel ring stacking interaction between Tyr35 of the active site helix and Phe199 of the C-terminal beta sheet in pcDHFR and by Tyr33 and Phe179 in mDHFR that is independent of bound ligand. A unique His174-His187 parallel ring stacking interaction was also observed only in the structure of pcDHFR. These ring stacking interactions are rarely found in any other protein families and may serve to enhance protein stability.     

Insight into the mechanism of inactivation and pH sensitivity in potassium channels from molecular dynamics simulations

Potassium (K (+)) channels can regulate ionic conduction through their pore by a mechanism, involving the selectivity filter, known as C-type inactivation. This process is rapid in the hERG K (+) channel and is fundamental to its physiological role. Although mutations within hERG are known to remove this process, a structural basis for the inactivation mechanism has yet to be characterized. Using MD simulations based on homology modeling, we observe that the carbonyl of the filter aromatic, Phe627, forming the S 0 K (+) binding site, swiftly rotates away from the conduction axis in the wild-type channel. In contrast, in well-characterized non-inactivating mutant channels, this conformational change occurs less frequently. In the non-inactivating channels, interactions with a water molecule located behind the selectivity filter are critical to the enhanced stability of the conducting state. We observe comparable conformational changes in the acid sensitive TASK-1 channel and propose a common mechanism in these channels for regulating efflux of K (+) ions through the selectivity filter.     

Probing the basis of domain-dependent inhibition using novel ketone inhibitors of Angiotensin-converting enzyme

Human angiotensin-converting enzyme (ACE) has two homologous domains, the N and C domains, with differing substrate preferences. X-ray crystal structures of the C and N domains complexed with various inhibitors have allowed identification of active site residues that might be important for the molecular basis of this selectivity. However, it is unclear to what extent the different residues contribute to substrate domain selectivity. Here, cocrystal structures of human testis ACE, equivalent to the C domain, have been determined with two novel C domain-selective ketomethylene inhibitors, (5 S)-5-[( N-benzoyl)amino]-4-oxo-6-phenylhexanoyl- l-tryptophan (kAW) and (5 S)-5-[( N-benzoyl)amino]-4-oxo-6-phenylhexanoyl- l-phenylalanine (kAF). The ketone groups of both inhibitors bind to the zinc ion as a hydrated geminal diolate, demonstrating the ability of the active site to catalyze the formation of the transition state. Moreover, active site residues involved in inhibitor binding have been mutated to their N domain counterparts, and the effect of the mutations on inhibitor binding has been determined. The C domain selectivity of these inhibitors was found to result from interactions between bulky hydrophobic side chain moieties and C domain-specific residues F391, V518, E376, and V380 (numbering of testis ACE). Mutation of these residues decreased the affinity for the inhibitors 4-20-fold. T282, V379, E403, D453, and S516 did not contribute individually to C domain-selective inhibitor binding. Further domain-selective inhibitor design should focus on increasing both the affinity and selectivity of the side chain moieties.     

A cation-pi interaction discriminates among sodium channels that are either sensitive or resistant to tetrodotoxin block

Voltage-gated sodium channels control the upstroke of the action potential in excitable cells of nerve and muscle tissue, making them ideal targets for exogenous toxins that aim to squelch electrical excitability. One such toxin, tetrodotoxin (TTX), blocks sodium channels with nanomolar affinity only when an aromatic Phe or Tyr residue is present at a specific location in the external vestibule of the ion-conducting pore. To test whether TTX is attracted to Tyr401 of NaV1.4 through a cation-pi interaction, this aromatic residue was replaced with fluorinated derivatives of Phe using in vivo nonsense suppression. Consistent with a cation-pi interaction, increased fluorination of Phe401, which reduces the negative electrostatic potential on the aromatic face, caused a monotonic increase in the inhibitory constant for block. Trifluorination of the aromatic ring decreased TTX affinity by approximately 50-fold, a reduction similar to that caused by replacement with the comparably hydrophobic residue Leu. Furthermore, we show that an energetically equivalent cation-pi interaction underlies both use-dependent and tonic block by TTX. Our results are supported by high level ab initio quantum mechanical calculations applied to a model of TTX binding to benzene. Our analysis suggests that the aromatic side chain faces the permeation pathway where it orients TTX optimally and interacts with permeant ions. These results are the first of their kind to show the incorporation of unnatural amino acids into a voltage-gated sodium channel and demonstrate that a cation-pi interaction is responsible for the obligate nature of an aromatic at this position in TTX-sensitive sodium channels.     

Structural characterizations of nonpeptidic thiadiazole inhibitors of matrix metalloproteinases reveal the basis for stromelysin selectivity.

The binding of two 5-substituted-1,3,4-thiadiazole-2-thione inhibitors to the matrix metalloproteinase stromelysin (MMP-3) have been characterized by protein crystallography. Both inhibitors coordinate to the catalytic zinc cation via an exocyclic sulfur and lay in an unusual position across the unprimed (P1-P3) side of the proteinase active site. Nitrogen atoms in the thiadiazole moiety make specific hydrogen bond interactions with enzyme structural elements that are conserved across all enzymes in the matrix metalloproteinase class. Strong hydrophobic interactions between the inhibitors and the side chain of tyrosine-155 appear to be responsible for the very high selectivity of these inhibitors for stromelysin. In these enzyme/inhibitor complexes, the S1'' enzyme subsite is unoccupied. A conformational rearrangement of the catalytic domain occurs that reveals an inherent flexibility of the substrate binding region leading to speculation about a possible mechanism for modulation of stromelysin activity and selectivity.     

Amino acids defining the acyl pocket of an invertebrate cholinesterase

Amphioxus (Branchiostoma floridae) cholinesterase 2 (ChE2) hydrolyzes acetylthiocholine (AsCh) almost exclusively. We constructed a homology model of ChE2 on the basis of Torpedo californica acetylcholinesterase (AChE) and found that the acyl pocket of the enzyme resembles that of Drosophila melanogaster AChE, which is proposed to be comprised of Phe330 (Phe290 in T. californica AChE) and Phe440 (Val400), rather than Leu328 (Phe288) and Phe330 (Phe290), as in vertebrate AChE. In ChE2, the homologous amino acids are Phe312 (Phe290) and Phe422 (Val400). To determine if these amino acids define the acyl pocket of ChE2 and its substrate specificity, and to obtain information about the hydrophobic subsite, partially comprised of Tyr352 (Phe330) and Phe353 (Phe331), we performed site-directed mutagenesis and in vitro expression. The aliphatic substitution mutant F312I ChE2 hydrolyzes AsCh preferentially but also butyrylthiocholine (BsCh), and the change in substrate specificity is due primarily to an increase in k_cat for BsCh; K_m and K_ss are also altered. F422L and F422V produce enzymes that hydrolyze BsCh and AsCh equally due to an increase in k_cat for BsCh and a decrease in k_cat for AsCh. Our data suggest that Phe312 and Phe422 define the acyl pocket. We also screened mutants for changes in sensitivity to various inhibitors. Y352A increases the sensitivity of ChE2 to the bulky inhibitor ethopropazine. Y352A decreases inhibition by BW284c51, consistent with its role as part of the choline-binding site. Aliphatic replacement mutations produce enzymes that are more sensitive to inhibition by iso-OMPA, presumably by increasing access to the active site serine. Y352A, F353A and F353V make ChE2 considerably more resistant to inhibition by eserine and neostigmine, suggesting that binding of these aromatic inhibitors is mediated by π–π or cation–π interactions at the hydrophobic site. Our results also provide information about the aromatic trapping of the active site histidine and the inactivation of ChE2 by sulfhydryl reagents.