Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A

Effects of Dipole Potential Modifiers on Heterogenic Lipid Bilayers

In this work, we examine the ability of dipole modifiers, flavonoids, and RH dyes to affect the dipole potential (φ d) and phase separation in membranes composed of ternary mixtures of POPC with different sphingolipids and sterols. Changes in the steady-state conductance induced by cation-ionophore complexes have been measured to evaluate the changes in dipole potential of planar lipid bilayers. Confocal fluorescence microscopy has been employed to investigate lipid segregation in giant unilamellar vesicles. The effects of flavonoids on φ d depend on lipid composition and dipole modifier type. The effectiveness of RH dyes to increase φ d depends on sphingolipid type but is not influenced by sterol content. Tested modifiers lead to partial or complete disruption of gel domains in bilayers composed of POPC, sphingomyelin, and cholesterol. Substitution of cholesterol to ergosterol or 7-dehydrocholesterol leads to a loss of fluidizing effects of modifiers except phloretin. This may be due to various compositions of gel domains. The lack of influence of modifiers on phase scenario in vesicles composed of ternary mixtures of POPC, cholesterol, and phytosphingosine or sphinganine is related to an absence of gel-like phase. It was concluded that the membrane lateral heterogeneity affects the dipole-modifying abilities of the agents that influence the magnitude of φ d by intercalation into the bilayer and orientation of its own large dipole moments (phloretin and RH dyes). The efficacy of modifiers that do not penetrate deeply and affect φ d through water adsorption (phlorizin, quercetin, and myricetin) is not influenced by lateral heterogeneity of membrane.

Unraveling the biomolecular snapshots of mitosis in healthy and cancer cells using plasmonically-enhanced Raman spectroscopy

Owing to the dynamic and complex nature of mitosis, precise and timely executions of biomolecular events are critical for high fidelity cell division. In this context, visualization of such complex events at the molecular level can provide vital information on the biomolecular processes in abnormal cells. Here, we explored the plasmonically enhanced light scattering properties of functionalized gold nanocubes (AuNCs) together with surface-enhanced Raman spectroscopy (SERS) to unravel the complex and dynamic biological processes involved in mitosis of healthy and cancerous cells from its molecular perspectives. By monitoring various stages of mitosis using SERS, we noticed that relatively high rate of conversion of mitotic proteins from their α-helix structure to β-sheet conformation is likely in the cancer cells during meta-, ana-, and telophases. Unique biochemical modifications to the lipid and amino acid moieties, associated with the observed protein conformational modifications, were also identified. However, in healthy cells, the existence of proteins in their β conformation was momentary and was largely in the α-helix form. The role of abnormal conformational modifications of mitotic proteins on the development of anomalous mitotic activities was further confirmed by looking at plasmonic nanoparticle-induced cytokinesis failure in cancer cells. Our findings illustrate the vast possibilities of SERS in real-time tracking of complex, subtle, and momentary modifications of biomolecules in live cells, which could provide new insights to the role of protein conformation dynamics during mitosis on the development of cancer and many other diseases.

The membrane interface dictates different anchor roles for "inner pair" and "outer pair" tryptophan indole rings in gramicidin A channels

We investigated the effects of substituting two of the four tryptophans (the "inner pair" Trp(9) and Trp(11) or the "outer pair" Trp(13) and Trp(15)) in gramicidin A (gA) channels. The conformational preferences of the doubly substituted gA analogues were assessed using circular dichroism spectroscopy and size-exclusion chromatography, which show that the inner tryptophans 9 and 11 are critical for the gA's conformational preference in lipid bilayer membranes. [Phe(13,15)]gA largely retains the single-stranded helical channel structure, whereas [Phe(9,11)]gA exists primarily as double-stranded conformers. Within this context, the (2)H NMR spectra from labeled tryptophans were used to examine the changes in average indole ring orientations, induced by the Phe substitutions and by the shift in conformational preference. Using a method for deuterium labeling of already synthesized gAs, we introduced deuterium selectively onto positions C2 and C5 of the remaining tryptophan indole rings in the substituted gA analogues for solid-state (2)H NMR spectroscopy. The (least possible) changes in orientation and overall motion of each indole ring were estimated from the experimental spectra. Regardless of the mixture of backbone folds, the indole ring orientations observed in the analogues are similar to those found previously for gA channels. Both Phe-substituted analogues form single-stranded channels, as judged from the formation of heterodimeric channels with the native gA. [Phe(13,15)]gA channels have Na(+) currents that are ~50% and lifetimes that are ~80% of those of native gA channels. The double-stranded conformer(s) of [Phe(9,11)]gA do not form detectable channels. The minor single-stranded population of [Phe(9,11)]gA forms channels with Na(+) currents that are ~25% and single-channel lifetimes that are ~300% of those of native gA channels. Our results suggest that Trp(9) and Trp(11), when "reaching" for the interface, tend to drive both monomer folding (to "open" a channel) and dimer dissociation (to "close" a channel). Furthermore, the dipoles of Trp(9) and Trp(11) are relatively more important for the single-channel conductance than are the dipoles of Trp(13) and Trp(15).

The gramicidin channel ion permeation free-energy profile: direct and indirect effects of CHARMM force field improvements

A revised CHARMM force field for tryptophan residues is studied as well as a new grid-based correction algorithm, called CMAP, using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl or KCl saline solution. The conformational stability of the interfacial side chains is studied, which shows good stability on the 10 ns time scale. The revised force field for the tryptophan side chain produces, in the decomposition, a Na(+) PMF(Trp) profile that is consonant with the prediction from the experimental results, analyzed with rate theory by Durrant et al. (2006), but in stark contrast to the prediction of the original CHARMM force field, version 22. However, the effect is diluted in the PMF profile due to indirect effects mediated by other components of the system (polypeptide, lipid molecules, ions, and water molecules). CMAP corrections to the L-amino acids help reduce the excessive translocation barrier. Decomposition demonstrates that this effect is due to effects on the K(+) PMF(H(2)O) profile rather than on the K(+) PMF(gA) profile. The results have been confirmed to be robust using an alternative umbrella-potential method. Further force field balancing efforts (direct and indirect) are required for future studies to evaluate whether these effects give rise to predictions that are consistent with those observables extracted from real experiments.

Noncontact dipole effects on channel permeation. III. Anomalous proton conductance effects in gramicidin

Proton transport on water wires, of interest for many problems in membrane biology, is analyzed in side-chain analogs of gramicidin A channels. In symmetrical 0.1N HCl solutions, fluorination of channel Trp(11), Trp-(13), or Trp(15) side chains is found to inhibit proton transport, and replacement of one or more Trps with Phe enhances proton transport, the opposite of the effects on K(+) transport in lecithin bilayers. The current-voltage relations are superlinear, indicating that some membrane field-dependent process is rate limiting. The interfacial dipole effects are usually assumed to affect the rate of cation translocation across the channel. For proton conductance, however, water reorientation after proton translocation is anticipated to be rate limiting. We propose that the findings reported here are most readily interpreted as the result of dipole-dipole interactions between channel waters and polar side chains or lipid headgroups. In particular, if reorientation of the water column begins with the water nearest the channel exit, this hypothesis explains the negative impact of fluorination and the positive impact of headgroup dipole on proton conductance.

The gramicidin ion channel: A model membrane protein

The linear peptide gramicidin forms prototypical ion channels specific for monovalent cations and has been extensively used to study the organization, dynamics and function of membrane-spanning channels. In recent times, the availability of crystal structures of complex ion channels has challenged the role of gramicidin as a model membrane protein and ion channel. This review focuses on the suitability of gramicidin as a model membrane protein in general, and the information gained from gramicidin to understand lipid-protein interactions in particular. Special emphasis is given to the role and orientation of tryptophan residues in channel structure and function and recent spectroscopic approaches that have highlighted the organization and dynamics of the channel in membrane and membrane-mimetic media.

Free energy simulations of single and double ion occupancy in gramicidin A

Simultaneous occupancy of the two binding sites in gramicidin A by monovalent cations is a well known property of this channel, but the energetic feasibility of this process in molecular dynamics simulations has not been established so far. Here the authors study the energetics of single and double ion occupancy in gramicidin A by constructing the potential of mean force for single and pair of cations. As representatives of small and large ions, they consider both Na+ and K+ ions in the calculations. Binding constants of ions are estimated from the free energy profiles. Comparisons with the experimental results indicate 3-4 kT discrepancy in the binding energies. They also study the coordination of the ions in their respective binding sites and the dynamic behavior of the channel water during the double ion binding process.

Tryptophan contributions to the empirical free-energy profile in gramicidin A/M heterodimer channels

Gramicidin A/gramicidin M heterodimer conductances were measured in planar lipid bilayers and found to form two distinguishable populations about halfway between the gramicidin A and gramicidin M homodimer conductances. This implies that the principle difference in the gramicidin A and gramicidin M transport free-energy profiles occurs at the channel center, where it would produce similar effects on the rate-limiting barrier for the two heterodimers. Kinetic analysis based on this and nearly all previously published homodimer conductance data for both gramicidin A and gramicidin M channels confirms this conclusion, indicating that the translocation step is approximately 100-fold slower in gramicidin M homodimers than in gramicidin A homodimers and that first- and second-ion exit-rate constants are higher by factors of 24 and 10, respectively. Assuming that the ratios of rate constants are related to the free-energy difference between gramicidin A and gramicidin M, we construct an effective ion-Trp free-energy interaction profile that has a minimum at the channel center.

Water alignment, dipolar interactions, and multiple proton occupancy during water-wire proton transport

A discrete multistate kinetic model for water-wire proton transport is constructed and analyzed using Monte Carlo simulations. In the model, each water molecule can be in one of three states: oxygen lone-pairs pointing leftward, pointing rightward, or protonated (H(3)O(+)). Specific rules for transitions among these states are defined as protons hop across successive water oxygens. Our model also includes water-channel interactions that preferentially align the water dipoles, nearest-neighbor dipolar coupling interactions, and Coulombic repulsion. Extensive Monte Carlo simulations were performed and the observed qualitative physical behaviors discussed. We find the parameters that allow the model to exhibit superlinear and sublinear current-voltage relationships, and show why alignment fields, whether generated by interactions with the pore interior or by membrane potentials, always decrease the proton current. The simulations also reveal a "lubrication" mechanism that suppresses water dipole interactions when the channel is multiply occupied by protons. This effect can account for an observed sublinear-to-superlinear transition in the current-voltage relationship.

Effects of phenylalanine substitutions in gramicidin A on the kinetics of channel formation in vesicles and channel structure in SDS micelles

The common occurrence of Trp residues at the aqueous-lipid interface region of transmembrane channels is thought to be indicative of its importance for insertion and stabilization of the channel in membranes. To further investigate the effects of Trp-->Phe substitution on the structure and function of the gramicidin channel, four analogs of gramicidin A have been synthesized in which the tryptophan residues at positions 9, 11, 13, and 15 are sequentially replaced with phenylalanine. The three-dimensional structure of each viable analog has been determined using a combination of two-dimensional NMR techniques and distance geometry-simulated annealing structure calculations. These phenylalanine analogs adopt a homodimer motif, consisting of two beta6.3 helices joined by six hydrogen bonds at their NH2-termini. The replacement of the tryptophan residues does not have a significant effect on the backbone structure of the channels when compared to native gramicidin A, and only small effects are seen on side-chain conformations. Single-channel conductance measurements have shown that the conductance and lifetime of the channels are significantly affected by the replacement of the tryptophan residues (Wallace, 2000; Becker et al., 1991). The variation in conductance appears to be caused by the sequential removal of a tryptophan dipole, thereby removing the ion-dipole interaction at the channel entrance and at the ion binding site. Channel lifetime variations appear to be related to changing side chain-lipid interactions. This is supported by data relating to transport and incorporation kinetics.

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).

Gramicidin A channel as a test ground for molecular dynamics force fields

We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K(+) ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K(+) ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60 ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose.

Solid-state 19F-NMR analysis of 19F-labeled tryptophan in gramicidin A in oriented membranes

The response of membrane-associated peptides toward the lipid environment or other binding partners can be monitored by solid-state NMR of suitably labeled side chains. Tryptophan is a prominent amino acid in transmembrane helices, and its (19)F-labeled analogues are generally biocompatible and cause little structural perturbation. Hence, we use 5F-Trp as a highly sensitive NMR probe to monitor the conformation and dynamics of the indole ring. To establish this (19)F-NMR strategy, gramicidin A was labeled with 5F-Trp in position 13 or 15, whose chi(1)/chi(2) torsion angles are known from previous (2)H-NMR studies. First, the alignment of the (19)F chemical shift anisotropy tensor within the membrane was deduced by lineshape analysis of oriented samples. Next, the three principal axes of the (19)F chemical shift anisotropy tensor were assigned within the molecular frame of the indole ring. Finally, determination of chi(1)/chi(2) for 5F-Trp in the lipid gel phase showed that the side chain alignment differs by up to 20 degrees from its known conformation in the liquid crystalline state. The sensitivity gain of (19)F-NMR and the reduction in the amount of material was at least 10-fold compared with previous (2)H-NMR studies on the same system and 100-fold compared with (15)N-NMR.

The role of Trp side chains in tuning single proton conduction through gramicidin channels

We present an extensive set of measurements of proton conduction through gramicidin A (gA), B (gB), and M (gM) homodimer channels which have 4, 3, or 0 Trp residues at each end of the channel, respectively. In gA we find a shoulder separating two domains of conductance increasing with concentration, confirming the results of Eisenman, G., B. Enos, J. Hagglund, and J. Sandblom. 1980. Ann. NY. Acad. Sci. 339:8-20. In gB, the shoulder is shifted by approximately 1/2 pH unit to higher H(+) concentrations and is very sharply defined. No shoulder appears in the gM data, but an associated transition from sublinear to superlinear I-V values occurs at a 100-fold higher [H(+)] in gM than in gA. The data in the low concentration domain are analyzed using a configuration space model of single-proton conduction, assuming that the difference in the proton potential of mean force (PMF) between gA and its analogs is constant, similar to the results of Anderson, D., R. B. Shirts, T. A. Cross, and D. D. Busath. 2001. Biophys. J. 81:1255-1264. Our results suggest that the average amplitudes of the calculated proton PMFs are nearly correct, but that the water reorientation barrier calculated for gA by molecular dynamics using the PM6 water model (Pomès, R., and B. Roux. 1997. Biophys. J. 72:246a) must be reduced in amplitude by 1.5 kcal/mol or more, and is not rate-limiting for gA.

Noncontact dipole effects on channel permeation. IV. Kinetic model of 5F-Trp(13) gramicidin A currents

Nonlinear least squares fitting was used to assign rate constants for the three-barrier, two-site, double-occupancy, single-filing kinetic model for previously reported current-voltage relations of (5F-Indole)Trp(13) gramicidin A and gramicidin A channels (, 75:2830-2844). By judicious coupling of parameters, it was possible to reduce the parameter space from 64 parameters to 24, and a reasonable fit consistent with other experimental data was obtained. The main features of the fit were that fluorination increased the rate constant for translocation by a factor of 2.33, consistent with a free energy change in the translocation barrier of -0.50 kcal/mol, and increased first-ion binding affinity by a factor of 1.13, primarily by decreasing the first-ion exit rate constant. The translocation rate constant was 5.62 times slower in diphytanoyl phosphatidylcholine (DPhPC) bilayers than in monoolein (GMO) bilayers (coupled for the four combinations of peptide and salt), suggesting a 44.2-mV difference in the projection of the interfacial dipole into the channel. Thus fluorination caused increased currents in DPhPC bilayers, where a high interfacial dipole potential makes translocation more rate limiting because the translocation barrier was reduced, and decreased currents in GMO bilayers, where ion exit or entry is rate limiting because these barriers were increased.

Noncontact dipole effects on channel permeation. V. Computed potentials for fluorinated gramicidin

Experimental and theoretical calculations indicate that the dipole moment of the four Trp side chains in gramicidin A (gA) channels modify channel conductance through long-range electrostatic interactions. Electrostatic ion/side-chain interaction energies along the channel were computed with CHARMM using ab initio atom charges for native and 4-, 5-, or 6-fluorinated Trp side chains. The bulk water reaction to the polar side chains was included using the method of images as implemented by, and channel waters in idealized structures were included. Ion/Trp interaction energies were approximately -0.6 kcal/mol throughout the channel for all four of the native Trp pairs. Channel waters produced a modest reduction in the magnitude of interactions, essentially offsetting images representing the bulk water outside the channel. The effects of side-chain fluorination depended on ring position and, to a lesser extent, residue number. Compared with native Trp, 5-fluorination reduces the translocation barrier with minor effects on the exit barrier. In contrast, 6-fluorination primarily reduces exit barrier. 4-Fluorination produces a more complex double-well energy profile. Effects of measured side-chain movements resulting from fluorination or change in lipid bilayer were negligible whereas thermal side chain librations cause large effects, especially in the region of the ion-binding sites.

Comparison of gramicidin A and gramicidin M channel conductance dispersities

To explore the possible role of Trp side chains in gramicidin channel conductance dispersity, we studied the dispersity of gramicidin M (gM), a gramicidin variant in which all four tryptophan residues are replaced with phenylalanine residues, and its enantiomer, gramicidin M(-) (gM(-)), and compared them to that of gramicidin A (gA). The conductances of highly purified gM and gM(-) were studied in alkali metal solutions at a variety of concentrations and voltages, in seven different types of lipid, and in the presence of detergent. Like gA channels, the most common gM channel conductance forms a narrow band. However, unlike gA channels, where the remaining 5-30% of channel conductances are broadly distributed below (and slightly above) the main band, in gM there is a narrow secondary band with <50% of the main peak conductance. This secondary peak was prominent in NaCl and KCl, but significantly diminished in CsCl and RbCl. Under some conditions, minor components can be observed with conductances yet lower than the secondary peak. Interconversions between the primary conductance state and these yet lower conductance states were observed. The current-voltage relations for both primary and secondary gM channel types have about the same curvature. The mean lifetime of the secondary channel type is below one third that of the primary type. The variants represent state deviations in the peptide or adjacent lipid structure.

Molecular determinants of electrical rectification of single channel conductance in gap junctions formed by connexins 26 and 32

The fully open state of heterotypic gap junction channels formed by pairing cells expressing connexin 32 (Cx32) with those expressing connexin 26 (Cx26) rectifies in a way that cannot be predicted from the current-voltage (I-V) relation of either homotypic channel. Using a molecular genetic analysis, we demonstrate that charged amino acids positioned in the amino terminus (M1 and D2) and first extracellular loop (E42) are major determinants of the current-voltage relation of the fully open state of homotypic and heterotypic channels formed by Cx26 and Cx32. The observed I-V relations of wild-type and mutant channels were closely approximated by those obtained with the electrodiffusive model of Chen and Eisenberg (Chen, D., and R. Eisenberg. 1993. Biophys. J. 64:1405-1421), which solves the Poisson-Nernst-Plank equations in one dimension using charge distribution models inferred from the molecular analyses. The rectification of the Cx32/Cx26 heterotypic channel results from the asymmetry in the number and position of charged residues. The model required the incorporation of a partial charge located near the channel surface to approximate the linear I-V relation observed for the Cx32*Cx26E1 homotypic channel. The best candidate amino acid providing this partial charge is the conserved tryptophan residue (W3). Incorporation of the partial charge of residue W3 and the negative charge of the Cx32E41 residue into the charge profile used in the Poisson-Nernst-Plank model of homotypic Cx32 and heterotypic Cx26/Cx32 channels resulted in I-V relations that closely resembled the observed I-V relations of these channels. We further demonstrate that some channel substates rectify. We suggest that the conformational changes associated with transjunctional voltage (V(j))-dependent gating to these substates involves a narrowing of the cytoplasmic entry of the channel that increases the electrostatic effect of charges in the amino terminus. The rectification that is observed in the Cx32/Cx26 heterotypic channel is similar although less steep than that reported for some rectifying electrical synapses. We propose that a similar electrostatic mechanism, which results in rectification through the open and substates of heterotypic channels, is sufficient to explain the properties of steeply rectifying electrical synapses.

Modulating dipoles for structure-function correlations in the gramicidin A channel

Dipoles of the tryptophan indole side chains have a direct impact on ion conductance in the gramicidin channel. Here, fluorination of the indoles (both 5- and 6-fluoro) is used to manipulate both the orientations and the magnitudes of the dipoles. The orientations and positions with respect to the channel axis were determined using (2)H solid state NMR of uniformly aligned lipid bilayer preparations. By exchange of the remaining four protons in the indole ring for deuterium, comparison could be made to d(5)-indole spectra that have previously been recorded for each of the four indoles of gramicidin A. After making the assignments which were aided by the observation of (19)F-(2)H dipolar interactions, we found that fluorination caused only minor changes in side chain conformation. With the high-resolution structural characterization of the fluorinated indoles in position 11, 13, and 15, the electrostatic interactions with a cation at the channel and bilayer center can be predicted and the influence of the modified dipoles on ion conductance estimated. The importance of the long-range electrostatic interaction was recently documented with the observation of alpha-helical dipoles oriented toward the bilayer center on the ion conductance pathway for the Streptomyces K(+) channel. We present direct measurements of the orientation of gramicidin channel F-Trp positions for use in analysis of dipole effects on channel permeation.

Lorentzian noise in single gramicidin A channel formamidinium currents

Seoh & Busath (1995) showed that in the presence of formamidinium, single gramicidin A channels were lengthened, had uniformly noisy currents at low voltages and had superlinear current-voltage relationships, all three properties being absent in gramicidin M- channels in which the interfacial tryptophan residues in gramicidin A are all replaced by phenylalanine. We measured the single channel noise power spectra (PSDs) in small monoolein (GMO) bilayers with formamidinium chloride solutions to help identify the mechanism of noise process. PSDs were Lorentzian with characteristic frequencies of 0.1-1.0 kHz in 0.1 and 0.3 M formamidinium chloride solutions, and from. 1-6 kHz in 1 M solution. Si(0), where measurable, ranged from approximately 50-200 fA2/Hz. The time course of the noise process could not be detected in these experiments. The low fc suggests slow motions or rare states of the blocking 'gates' which, judging from the result with gramicidin M-, must be equal to or related to the Trp residues.