Christopher E. Dempsey

Structure, Location, and Lipid Perturbations of Melittin at the Membrane Interface

Melittin is arguably the most widely studied amphipathic, membrane-lytic alpha-helical peptide. Although several lines of evidence suggest an interfacial membrane location at low concentrations, melittins exact position and depth of penetration into the hydrocarbon core are unknown. Furthermore, the structural basis for its lytic action remains largely a matter of conjecture. Using a novel x-ray absolute-scale refinement method, we have now determined the location, orientation, and likely conformation of monomeric melittin in oriented phosphocholine lipid multilayers. Its helical axis is aligned parallel to the bilayer plane at the depth of the glycerol groups, but its average conformation differs from the crystallographic structure. As observed earlier for another amphipathic alpha-helical peptide, the lipid perturbations induced by melittin are remarkably modest. Small bilayer perturbations thus appear to be a general feature of amphipathic helices at low concentrations. In contrast, a dimeric form of melittin causes larger structural perturbations under otherwise identical conditions. These results provide direct structural evidence that self-association of amphipathic helices may be the crucial initial step toward membrane lysis.

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.

Preferential closed channel blockade of HERG potassium currents by chemically synthesised BeKm-1 scorpion toxin

The scorpion toxin peptide BeKm‐1 was synthesised by fluorenylmethoxycarbonyl solid phase chemistry and folded by air oxidation. The peptides effects on heterologous human ether‐a‐go‐go‐related gene potassium current (I HERG) in HEK293 cells were assessed using ‘whole‐cell’ patch clamp. Blockade of I HERG by BeKm‐1 was concentration‐dependent, temperature‐dependent, and rapid in onset and reversibility. Blockade also exhibited inverse voltage dependence, inverse dependence on duration of depolarisation, and reverse use‐ and frequency‐dependence. Blockade by BeKm‐1 and recombinant ergtoxin, another scorpion toxin known to block HERG, differed in their recovery from HERG current inactivation elicited by strong depolarisation and in their ability to block HERG when the channels were already activated. We conclude that synthetic BeKm‐1 toxin blocks HERG preferentially through a closed (resting) state channel blockade mechanism, although some open channel blockade also occurs.

Hydrogen Bonding in Helical Polypeptides from Molecular Dynamics Simulations and Amide Hydrogen Exchange Analysis: Alamethicin and Melittin in Methanol

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NHs of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NHs of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.

Origin of low mammalian cell toxicity in a class of highly active antimicrobial amphipathic helical peptides

We recently described a novel antimicrobial peptide, RTA3, derived from the commensal organism Streptococcus mitis, with strong anti-Gram-negative activity, low salt sensitivity, and minimal mammalian cell toxicity in vitro and in vivo. This peptide conforms to the positively charged, amphipathic helical peptide motif, but has a positively charged amino acid (Arg-5) on the nonpolar face of the helical structure that is induced upon membrane binding. We surmised that disruption of the hydrophobic face with a positively charged residue plays a role in minimizing eukaryotic cell toxicity, and we tested this using a mutant with an R5L substitution. The greatly enhanced toxicity in the mutant peptide correlated with its ability to bind and adopt helical conformations upon interacting with neutral membranes; the wild type peptide RTA3 did not bind to neutral membranes (binding constant reduced by at least 1000-fold). Spectroscopic analysis indicates that disruption of the hydrophobic face of the parent peptide is accommodated in negatively charged membranes without partial peptide unfolding. These observations apply generally to amphipathic helical peptides of this class as we obtained similar results with a peptide and mutant pair (Chen, Y., Mant, C. T., Farmer, S. W., Hancock, R. E., Vasil, M. L., and Hodges, R. S. (2005) J. Biol. Chem. 280, 12316–12329) having similar structural properties. In contrast to previous interpretations, we demonstrate that these peptides simply do not bind well to membranes (like those of eukaryotes) with exclusively neutral lipids in their external bilayer leaflet. We highlight a significant role for tryptophan in promoting binding of amphipathic helical peptides to neutral bilayers, augmenting the arsenal of strategies to reduce mammalian toxicity in antimicrobial peptides.

Specificity of Ion−Protein Interactions: Complementary and Competitive Effects of Tetrapropylammonium, Guanidinium, Sulfate, and Chloride Ions

The interactions of ions with a model peptide (a single melittin alpha-helix) in solutions of tetrapropylammonium sulfate or guanidinium chloride were examined by molecular dynamics simulations. The tetrapropylammonium cation shares the geometrical property of essentially flat faces with the previously examined guanidinium cation, and it was found that that this geometry leads to a strong preference for tetrapropylammonium to interact in a similar stacking-type fashion with flat nonpolar groups such as the indole side chain of tryptophan. In contrast to guanidinium, however, tetrapropylammonium does not exhibit strong ion pairing or clustering with sulfate counterions in the solution. Sulfate was found to interact almost exclusively and strongly with the cationic groups of the peptide, such that, already in a 0.1 m solution of tetrapropylammonium sulfate, the 6+ charge of the peptide is effectively locally neutralized. In combination with previous simulations, neutron scattering studies, and experiments on the conformational stability of model peptides, the present results suggest that the Hofmeister series can be explained in higher detail by splitting ions according to the effect they have on hydrogen bonding, salt bridges, and hydrophobic interactions in the protein and how these effects are altered by the counterion.

Direct 1H NMR evidence for conversion of β-d-cellobiose to cellobionolactone by cellobiose dehydrogenase from Phanerochaete chrysosporium

The α‐ and β‐anomers of d‐cellobiose were resolved by 1H NMR spectroscopy. Addition of cellobiose dehydrogenase purified from the white‐rot P. chrysosporium led to selective conversion of β‐d‐cellobiose. The product was identical to cellobionolactone as synthesized from Ca‐cellobionate. Overnight incubation of the product led to an altered NMR spectrum, which was also obtained by incubation of cellobionolactone. The new spectrum matched that for Ca‐cellobionate. The instability of cellobionolactone explains the detection of cellobionic acid as product in earlier studies.

A combined X-ray and neutron diffraction study of selectively deuterated melittin in phospholipid bilayers: effect of pH

In order to study consequences of protonation of the N-terminus upon the interaction of the bee venom melittin with phospholipid bilayers, analogues of melittin, some of which were specifically deuterated at either Ala-12 or 15, were synthesized. These peptides were incorporated into bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine at either low pH (N-terminus protonated) or high pH (N-terminus unprotonated). X-ray and neutron diffraction data were collected from ordered stacks of these bilayers and from peptide-free controls. Phase determination was carried out using the swelling series (X-ray) and isomorphous derivative (neutron) methods. The water distribution between adjacent bilayers in the stacks may be described by a pair of Gaussians whose position and width change with the protonation state of the melittin. Difference Fourier profiles reveal that the melittin largely incorporates into the phospholipid bilayers. Changes in the water, melittin and deuterium label distributions fit a model in which the melittin lies both at the surface and close to the centre of the bilayer, the distribution of peptide between these locations being pH-dependent, with a larger population of surface melittin when the N-terminus is unprotonated.

hERG potassium channel blockade by the HCN channel inhibitor bradycardic agent ivabradine.

Background Ivabradine is a specific bradycardic agent used in coronary artery disease and heart failure, lowering heart rate through inhibition of sinoatrial nodal HCN‐channels. This study investigated the propensity of ivabradine to interact with KCNH2‐encoded human Ether‐à‐go‐go–Related Gene (hERG) potassium channels, which strongly influence ventricular repolarization and susceptibility to torsades de pointes arrhythmia. Methods and Results Patch clamp recordings of hERG current (IhERG) were made from hERG expressing cells at 37°C. IhERG was inhibited with an IC50 of 2.07 μmol/L for the hERG 1a isoform and 3.31 μmol/L for coexpressed hERG 1a/1b. The voltage and time‐dependent characteristics of IhERG block were consistent with preferential gated‐state‐dependent channel block. Inhibition was partially attenuated by the N588K inactivation‐mutant and the S624A pore‐helix mutant and was strongly reduced by the Y652A and F656A S6 helix mutants. In docking simulations to a MthK‐based homology model of hERG, the 2 aromatic rings of the drug could form multiple π‐π interactions with the aromatic side chains of both Y652 and F656. In monophasic action potential (MAP) recordings from guinea‐pig Langendorff‐perfused hearts, ivabradine delayed ventricular repolarization and produced a steepening of the MAPD90 restitution curve. Conclusions Ivabradine prolongs ventricular repolarization and alters electrical restitution properties at concentrations relevant to the upper therapeutic range. In absolute terms ivabradine does not discriminate between hERG and HCN channels: it inhibits IhERG with similar potency to that reported for native If and HCN channels, with S6 binding determinants resembling those observed for HCN4. These findings may have important implications both clinically and for future bradycardic drug design.

Assessing hERG Pore Models As Templates for Drug Docking Using Published Experimental Constraints: The Inactivated State in the Context of Drug Block

Many structurally and therapeutically diverse drugs interact with the human heart K+ channel hERG by binding within the K+ permeation pathway of the open channel, leading to drug-induced ‘long QT syndrome’. Drug binding to hERG is often stabilized by inactivation gating. In the absence of a crystal structure, hERG pore homology models have been used to characterize drug interactions. Here we assess potentially inactivated states of the bacterial K+ channel, KcsA, as templates for inactivated state hERG pore models in the context of drug binding using computational docking. Although Flexidock and GOLD docking produced low energy score poses in the models tested, each method selected a MthK K+ channel-based model over models based on the putative inactivated state KcsA structures for each of the 9 drugs tested. The variety of docking poses found indicates that an optimal arrangement for drug binding of aromatic side chains in the hERG pore can be achieved in several different configurations. This plasticity of the drug “binding site” is likely to be a feature of the hERG inactivated state. The results demonstrate that experimental data on specific drug interactions can be used as structural constraints to assess and refine hERG homology models.