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

Natural and Synthetic Halogenated Amino Acids-Structural and Bioactive Features in Antimicrobial Peptides and Peptidomimetics

The 3D structure and surface characteristics of proteins and peptides are crucial for interactions with receptors or ligands and can be modified to some extent to modulate their biological roles and pharmacological activities. The introduction of halogen atoms on the side-chains of amino acids is a powerful tool for effecting this type of tuning, influencing both the physico-chemical and structural properties of the modified polypeptides, helping to first dissect and then rationally modify features that affect their mode of action. This review provides examples of the influence of different types of halogenation in amino acids that replace native residues in proteins and peptides. Examples of synthetic strategies for obtaining halogenated amino acids are also provided, focusing on some representative compounds and their biological effects. The role of halogenation in native and designed antimicrobial peptides (AMPs) and their mimetics is then discussed. These are in the spotlight for the development of new antimicrobial drugs to counter the rise of antibiotic-resistant pathogens. AMPs represent an interesting model to study the role that natural halogenation has on their mode of action and also to understand how artificially halogenated residues can be used to rationally modify and optimize AMPs for pharmaceutical purposes.

Recombinant expression, antimicrobial activity and mechanism of action of tritrpticin analogs containing fluoro-tryptophan residues

The increase in antibiotic-resistant bacterial infections has prompted significant academic research into new therapeutic agents targeted against these pathogens. Antimicrobial peptides (AMPs) appear as promising candidates, due their potent antimicrobial activity and their ubiquitous presence in almost all organisms. Tritrpticin is a member of this family of peptides and has been shown to exert a strong antimicrobial activity against several bacterial strains. Tritrpticin's main structural characteristic is the presence of three consecutive Trp residues at the center of the peptide. These residues play an important role in the activity of tritrpticin against Escherichia coli. In this work, a recombinant version of tritrpticin was produced in E. coli using calmodulin as a fusion protein expression tag to overcome the toxicity of the peptide. When used in combination with glyphosate, an inhibitor of the endogenous synthesis of aromatic amino acids, this expression system allowed for the incorporation of fluorinated Trp analogs at very high levels (>90%). The antimicrobial activity of the 4-, 5- and 6-fluoro-Trp-containing tritrpticins against E. coli was as strong as the activity of the native peptide. Similarly, the tritrpticin analogs exhibited comparable abilities to perturb and permeabilize synthetic lipid bilayers as well as the outer and inner membrane of E. coli. Furthermore, the use of 19F NMR spectroscopy established that each individual fluoro-Trp residue interacts differently with SDS micelles, supporting the idea that each Trp in the original tritrpticin plays a different role in the perturbing/permeabilizing activity of the peptide. Moreover, our work demonstrates that the use of fluoro-Trp in solvent perturbation 19F NMR experiments provides detailed site-specific information on the insertion of the Trp residues in biological membrane mimetics. This article is part of a Special Issue entitled: Antimicrobial peptides edited by Karl Lohner and Kai Hilpert.

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.

Monitoring membrane protein conformational heterogeneity by fluorescence lifetime distribution analysis using the maximum entropy method

Due to the inherent difficulty in crystallizing membrane proteins, approaches based on fluorescence spectroscopy have proved useful in elucidating their conformational characteristics. The ion channel peptide gramicidin serves as an excellent prototype for monitoring membrane protein conformation and dynamics due to a number of reasons. We have analyzed conformational heterogeneity in membrane-bound gramicidin using fluorescence lifetime distribution analysis of tryptophan residues by the maximum entropy method (MEM). MEM represents a model-free and robust approach for analyzing fluorescence lifetime distribution. In this paper, we show for the first time, that fluorescence lifetime distribution analysis using MEM could be a convenient approach to monitor conformational heterogeneity in membrane-bound gramicidin in particular and membrane proteins in general. Lifetime distribution analysis by MEM therefore provides a novel window to monitor conformational transitions in membrane proteins.

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.

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.

Monitoring ion channel conformations in membranes utilizing a novel dual fluorescence quenching approach

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. We have analyzed the localization of the functionally important tryptophan residues of the membrane-bound channel and non-channel conformations of gramicidin utilizing a novel dual fluorescence quenching approach [G.A. Caputo, E. London, Biochemistry 42 (2003) 3265-3274]. In this paper, we show for the first time that the dual quenching approach is applicable to multiple tryptophan containing functional ion channel peptides such as gramicidin. Importantly, dual quenching is found to be sensitive to the membrane-bound conformations of this important model ion channel.

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.

The structure, cation binding, transport, and conductance of Gly15-gramicidin A incorporated into SDS micelles and PC/PG vesicles

To further investigate the effect of single amino acid substitution on the structure and function of the gramicidin channel, an analogue of gramicidin A (GA) has been synthesized in which Trp(15) is replaced by Gly in the critical aqueous interface and cation binding region. The structure of Gly(15)-GA incorporated into SDS micelles has been determined using a combination of 2D-NMR spectroscopy and molecular modeling. Like the parent GA, Gly(15)-GA forms a dimeric channel composed of two single-stranded, right-handed beta(6.3)-helices joined by hydrogen bonds between their N-termini. The replacement of Trp(15) by Gly does not have a significant effect on backbone structure or side chain conformations with the exception of Trp(11) in which the indole ring is rotated away from the channel axis. Measurement of the equilibrium binding constants and Delta G for the binding of monovalent cations to GA and Gly(15)-GA channels incorporated into PC vesicles using (205)Tl NMR spectroscopy shows that monovalent cations bind much more weakly to the Gly(15)-GA channel entrance than to GA channels. Utilizing the magnetization inversion transfer NMR technique, the transport of Na(+) ions through GA and Gly(15)-GA channels incorporated into PC/PG vesicles has been investigated. The Gly(15) substitution produces an increase in the activation enthalpy of transport and thus a significant decrease in the transport rate of the Na(+) ion is observed. The single-channel appearances show that the conducting channels have a single, well-defined structure. Consistent with the NMR results, the single-channel conductances are reduced by 30% and the lifetimes by 70%. It is concluded that the decrease in cation binding, transport, and conductance in Gly(15)-GA results from the removal of the Trp(15) dipole and, to a lesser extent, the change in orientation of Trp(11).

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.

Noncontact dipole effects on channel permeation. VI. 5F- and 6F-Trp gramicidin channel currents

Fluorination of peptide side chains has been shown to perturb gramicidin channel conductance without significantly changing the average side chain structure, which, it is hoped, will allow detailed analysis of electrostatic modulation of current flow. Here we report a 1312-point potassium current-voltage-concentration data set for homodimeric channels formed from gramicidin A (gA) or any of eight fluorinated Trp analogs in both lecithin and monoglyceride bilayers. We fit the data with a three-barrier, two-site, two-ion (3B2S) kinetic model. The fluorination-induced changes in the rate constants were constrained by the same factor in both lipids. The rate constant changes were converted to transition-state free-energy differences for comparison with previous electrostatic potential energy differences based on an ab initio force field. The model allowed a reasonably good fit (chi = 8.29 with 1271 degrees of freedom). The measured changes were subtle. Nevertheless, the fitted energy perturbations agree well with electrostatic predictions for five of the eight peptides. For the other three analogs, the fitted changes suggested a reduced translocation barrier rather than the reduced exit barrier as predicted by electrostatics.

Continuum electrostatics fails to describe ion permeation in the gramicidin channel

We investigate the validity of continuum electrostatics in the gramicidin A channel using a recently determined high-resolution structure. The potential and electric field acting on ions in and around the channel are computed by solving Poisson's equation. These are then used in Brownian dynamics simulations to obtain concentration profiles and the current passing through the channel. We show that regardless of the effective dielectric constant used for water in the channel or the channel protein, it is not possible to reproduce all the experimental data on gramicidin A; thus, continuum electrostatics cannot provide a valid framework for the description of ion dynamics in gramicidin channels. Using experimental data and molecular dynamics simulations as guides, we have constructed potential energy profiles that can satisfactorily describe the available physiological data. These profiles provide useful benchmarks for future potential of mean force calculations of permeating ions from molecular dynamics simulations of gramicidin A. They also offer a convenient starting point for studying structure-function relationships in modified gramicidin channels.

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.

Model channel ion currents in NaCl-extended simple point charge water solution with applied-field molecular dynamics

Using periodic boundary conditions and a constant applied field, we have simulated current flow through an 8.125-A internal diameter, rigid, atomistic channel with polar walls in a rigid membrane using explicit ions and extended simple point charge water. Channel and bath currents were computed from 10 10-ns trajectories for each of 10 different conditions of concentration and applied voltage. An electric field was applied uniformly throughout the system to all mobile atoms. On average, the resultant net electric field falls primarily across the membrane channel, as expected for two conductive baths separated by a membrane capacitance. The channel is rarely occupied by more than one ion. Current-voltage relations are concentration dependent and superlinear at high concentrations.

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.