Ion selectivity of the anthrax toxin channel and its effect on protein translocation.

Anthrax toxin consists of three ∼ 85-kD proteins: lethal factor (LF), edema factor (EF), and protective antigen (PA). PA63 (the 63-kD, C-terminal portion of PA) forms heptameric channels ((PA63)7) in planar phospholipid bilayer membranes that enable the translocation of LF and EF across the membrane. These mushroom-shaped channels consist of a globular cap domain and a 14-stranded ß-barrel stem domain, with six anionic residues lining the interior of the stem to form rings of negative charges. (PA63)7 channels are highly cation selective, and, here, we investigate the effects on both cation selectivity and protein translocation of mutating each of these anionic residues to a serine. We find that although some of these mutations reduce cation selectivity, selectivity alone does not directly predict the rate of protein translocation; local changes in electrostatic forces must be considered as well.

Ion conductance of the stem of the anthrax toxin channel during lethal factor translocation.

The tripartite anthrax toxin consists of protective antigen, lethal factor (LF), and edema factor. PA63 (the 63-kDa, C-terminal part of protective antigen) forms heptameric channels in cell membranes that allow for the transport of LF and edema factor into the cytosol. These channels are mushroom shaped, with a ring of seven phenylalanine residues (known as the phenylalanine clamp) lining the junction between the cap and the stem. It is known that when LF is translocated through the channel, the phenylalanine clamp creates a seal that causes an essentially complete block of conduction. In order to examine ion conductance in the stem of the channel, we used Venus yellow fluorescent protein as a molecular stopper to trap LFN (the 30-kDa, 263-residue N-terminal segment of LF), as well as various truncated constructs of LFN, in mutant channels in which the phenylalanine clamp residues were mutated to alanines. Here we present evidence that ion movement occurs within the channel stem (but is stopped, of course, at the phenylalanine clamp) during protein translocation. Furthermore, we also propose that the lower region of the stem plays an important role in securing peptide chains during translocation.

A molecular mechanics approach to modeling protein-ligand interactions: relative binding affinities in congeneric series.

We introduce the "Prime-ligand" method for ranking ligands in congeneric series. The method employs a single scoring function, the OPLS-AA/GBSA molecular mechanics/implicit solvent model, for all stages of sampling and scoring. We evaluate the method using 12 test sets of congeneric series for which experimental binding data is available in the literature, as well as the structure of one member of the series bound to the protein. Ligands are "docked" by superimposing a common stem fragment among the compounds in the series using a crystal complex from the Protein Data Bank and sampling the conformational space of the variable region. Our results show good correlation between our predicted rankings and the experimental data for cases in which binding affinities differ by at least 1 order of magnitude. For 11 out of 12 cases, >90% of such ligand pairs could be correctly ranked, while for the remaining case, Factor Xa, 76% of such pairs were correctly ranked. A small number of compounds could not be docked using the current protocol because of the large size of functional groups that could not be accommodated by a rigid receptor. CPU requirements for the method, involving CPU minutes per ligand, are modest compared with more rigorous methods that use similar force fields, such as free energy perturbation. We also benchmark the scoring function using series of ligands bound to the same protein within the CSAR data set. We demonstrate that energy minimization of ligands in the crystal structures is critical to obtain any correlation with experimentally determined binding affinities.