Boris S. Zhorov

Kinked-Helices Model of the Nicotinic Acetylcholine Receptor Ion Channel and Its Complexes with Blockers: Simulation by the Monte Carlo Minimization Method

A model of the nicotinic acetylcholine receptor ion channel was elaborated based on the data from electron microscopy, affinity labeling, cysteine scanning, mutagenesis studies, and channel blockade. A restrained Monte Carlo minimization method was used for the calculations. Five identical M2 segments (the sequence EKMTLSISVL10LALTVFLLVI20V) were arranged in five-helix bundles with various geometrical profiles of the pore. For each bundle, energy profiles for chlorpromazine, QX-222, pentamethonium, and other blocking drugs pulled through the pore were calculated. An optimal model obtained allows all of the blockers free access to the pore, but retards them at the rings of residues known to contribute to the corresponding binding sites. In this model, M2 helices are necessarily kinked. They come into contact with each other at the cytoplasmic end but diverge at the synaptic end, where N-termini of M1 segments may contribute to the pore. The kinks disengage alpha-helical H-bonds between Ala12 and Ser8. The uncoupled lone electron pairs of Ser8 carbonyl oxygens protrude into the pore, forming a hydrophilic ring that may be important for the permeation of cations. A split network of H-bonds provides a flexibility to the chains Val9-Ala12, the numerous conformations of which form only two or three intrasegment H-bonds. The cross-ectional dimensions of the interface between the flexible chains vary essentially at the level of Leu11. We suggest that conformational transitions in the chains Val9-Ala12 are responsible for the channel gating, whereas rotations of more stable alpha-helical parts of M2 segments may be necessary to transfer the channel in the desensitized state.

Intersegment hydrogen bonds as possible structural determinants of the N/Q/R site in glutamate receptors.

Specific electrophysiological and pharmacological properties of ionic channels in NMDA, AMPA, and kainate subtypes of ionotropic glutamate receptors (GluRs) are determined by the Asn (N), Gln (Q), and Arg (R) residues located at homologous positions of the pore-lining M2 segments (the N/Q/R site). Presumably, the N/Q/R site is located at the apex of the reentrant membrane loop and forms the narrowest constriction of the pore. Although the shorter Asn residues are expected to protrude in the pore to a lesser extent than the longer Gln residues, the effective dimension of the NMDA channel (corresponding to the size of the largest permeant organic cation) is, surprisingly, smaller than that of the AMPA channel. To explain this paradox, we propose that the N/Q/R residues form macrocyclic structures (rings) stabilized by H-bonds between a NH(2) group in the side chain of a given M2 segment and a C==O group of the main chain in the adjacent M2 segment. Using Monte Carlo minimization, we have explored conformational properties of the rings. In the Asn, but not in the Gln ring, the side-chain oxygens protruding into the pore may facilitate ion permeation and accept H-bonds from the blocking drugs. In this way, the model explains different electrophysiological and pharmacological properties of NMDA and non-NMDA GluR channels. The ring of H-bonded polar residues at the pore narrowing resembles the ring of four Thr(75) residues observed in the crystallographic structure of the KcsA K(+) channel.

Architecture of the Neuronal Nicotinic Acetylcholine Receptor Ion Channel at the Binding Site of bis-Ammonium Blockers

Abstract. Structure-activity relationships of 56 pentamethylenbis-ammonium compounds, the blockers of the neuronal nicotinic acetylcholine receptor (nAChR) ion channel, have been studied to estimate the cross-sectional dimensions of the channel pore. The cat superior cervical sympathetic ganglion in situ and isolated guinea pig ileum were used to evaluate the potency of the compounds to block ganglionic transmission. Minimum-energy conformations of each compound were calculated by the molecular mechanics method. A topographic model of the binding site of the blockers was proposed. It incorporates two narrowings, a large and a small one. The small narrowing is located between the large one and the cytoplasmic end of the pore. The cross-sectional dimensions of the large and small narrowings estimated from the dimensions of the blockers are 6.1 × 8.3 Å and 5.5 × 6.4 Å, respectively, the distance between the narrowings along the pore being approximately 7 Å. Most potent blockers would occlude the pore via binding to the channel at the levels of both narrowings. Less potent blockers are either too large or too small to bind to both narrowings simultaneously: large blockers would occlude the pore at the level of large narrowing, while small blockers would pass the large narrowing and occlude the pore at the level of small narrowing only. A comparison of the topographic model with a molecular five-helix bundle model of nAChR pore predicts Serine and Threonine rings to be the most probable candidates for the large and small narrowings, respectively.

Conformational analysis of d-tubocurarine: implications for minimal dimensions of its binding site within ion channels.

All the minimum-energy conformations of d-tubocurarine were calculated by the method of molecular mechanics. The energy was minimized from 413 closed forms of the 18-member ring. The set of minimum-energy conformations includes 10 forms with energies less than 6 kcal/mol from the most stable one. Among the four lowest minimum-energy conformations, two forms correspond to those known from X-ray studies, whereas two conformations were not detected experimentally earlier. The flexibility of d-tubocurarine was estimated by calculating six paths of interconversion between the four lowest minimum-energy conformations. Using a molecular graphics technique, it was found that the most extended minimum-energy conformation of d-tubocurarine may fit in an ion channel of a rectangular profile of 8.7 × 11.2 Å, while one tetrahydroisoquinoline head may fit a profile as small as 6.9 × 11.0 Å. A possible model of d-tubocurarine location within the ion channel of the neuronal nicotinic acetylcholine receptor is suggested.

Relationships Between Activity of Organophosphorus Inhibitors of Acetylcholinesterase and Accessibility of Phosphorus Atom as Estimated by Molecular Mechanics Calculations

Irreversible inhibition of acetylcholinesterase (AChE) by organophosphorus compounds is believed to be due to phosphorylation of serine hydroxyl in the enzyme active center (Cohen, Oosterbaan, 1963). The attack by the hydroxyl may occur from the different faces of phosphorus atom (Berman and Decker, 1989). The efficiency of these reactions should depend on sterical accessibility of the phosphorus atom for the attacking agent. In particular, low activity of some inhibitors may be due to low accessibility of their phosphorus atom in the reaction with the enzyme. To obtain quantitative estimates of the accessibility of different faces of the phosphorus atom, conformational flexibility of inhibitors should be taken into account. In this work, using molecular mechanics method, we have calculated all minimum-energy conformations (conformers) of 15 organophosphorus inhibitors with nitro-phenyl and S-ethylthioethyl leaving groups and different structure of phosphoryl moiety, determined the accessibility of the phosphorus atom in every conformer, and compared anti-AChE efficiency of the compounds with the total population of those conformers which have accessible phosphorus atom.

Interaction of the BKCa channel gating ring with dendrotoxins

Two classes of small homologous basic proteins, mamba snake dendrotoxins (DTX) and bovine pancreatic trypsin inhibitor (BPTI), block the large conductance Ca2+-activated K+ channel (BKCa, KCa1.1) by production of discrete subconductance events when added to the intracellular side of the membrane. This toxin-channel interaction is unlikely to be pharmacologically relevant to the action of mamba venom, but as a fortuitous ligand-protein interaction, it has certain biophysical implications for the mechanism of BKCa channel gating. In this work we examined the subconductance behavior of 9 natural dendrotoxin homologs and 6 charge neutralization mutants of δ-dendrotoxin in the context of current structural information on the intracellular gating ring domain of the BKCa channel. Calculation of an electrostatic surface map of the BKCa gating ring based on the Poisson-Boltzmann equation reveals a predominantly electronegative surface due to an abundance of solvent-accessible side chains of negatively charged amino acids. Available structure-activity information suggests that cationic DTX/BPTI molecules bind by electrostatic attraction to site(s) on the gating ring located in or near the cytoplasmic side portals where the inactivation ball peptide of the β2 subunit enters to block the channel. Such an interaction may decrease the apparent unitary conductance by altering the dynamic balance of open versus closed states of BKCa channel activation gating.