Structural modeling of calcium binding in the selectivity filter of the L-type calcium channel.

Calcium channels play crucial physiological roles. In the absence of high-resolution structures of the channels, the mechanism of ion permeation is unknown. Here we used a method proposed in an accompanying paper (Cheng and Zhorov in Eur Biophys J, 2009) to predict possible chelation patterns of calcium ions in a structural model of the L-type calcium channel. We compared three models in which two or three calcium ions interact with the four selectivity filter glutamates and a conserved aspartate adjacent to the glutamate in repeat II. Monte Carlo energy minimizations yielded many complexes with calcium ions bound to at least two selectivity filter carboxylates. In these complexes calcium-carboxylate attractions are counterbalanced by calcium-calcium and carboxylate-carboxylate repulsions. Superposition of the complexes suggests a high degree of mobility of calcium ions and carboxylate groups of the glutamates. We used the predicted complexes to propose a permeation mechanism that involves single-file movement of calcium ions. The key feature of this mechanism is the presence of bridging glutamates that coordinate two calcium ions and enable their transitions between different chelating patterns involving four to six oxygen atoms from the channel protein. The conserved aspartate is proposed to coordinate a calcium ion incoming to the selectivity filter from the extracellular side. Glutamates in repeats III and IV, which are most distant from the repeat II aspartate, are proposed to coordinate the calcium ion that leaves the selectivity filter to the inner pore. Published experimental data and earlier proposed permeation models are discussed in view of our model.

Docking of calcium ions in proteins with flexible side chains and deformable backbones.

A method of docking Ca(2+) ions in proteins with flexible side chains and deformable backbones is proposed. The energy was calculated with the AMBER force field, implicit solvent, and solvent exposure-dependent and distance-dependent dielectric function. Starting structures were generated with Ca(2+) coordinates and side-chain torsions sampled in 1000 A(3) cubes centered at the experimental Ca(2+) positions. The energy was Monte Carlo-minimized. The method was tested on fourteen Ca(2+)-binding sites. For twelve Ca(2+)-binding sites the root mean square (RMS) deviation of the apparent global minimum from the experimental structure was below 1.3 and 1.7 A for Ca(2+) ions and side-chain heavy atoms, respectively. Energies of multiple local minima correlate with the RMS deviations from the X-ray structures. Two Ca(2+)-binding sites at the surface of proteinase K were not predicted, because of underestimation of Ca(2+) hydration energy by the implicit-solvent method.

Structural model for phenylalkylamine binding to L-type calcium channels.

Phenylalkylamines (PAAs), a major class of L-type calcium channel (LTCC) blockers, have two aromatic rings connected by a flexible chain with a nitrile substituent. Structural aspects of ligand-channel interactions remain unclear. We have built a KvAP-based model of LTCC and used Monte Carlo energy minimizations to dock devapamil, verapamil, gallopamil, and other PAAs. The PAA-LTCC models have the following common features: (i) the meta-methoxy group in ring A, which is proximal to the nitrile group, accepts an H-bond from a PAA-sensing Tyr_IIIS6; (ii) the meta-methoxy group in ring B accepts an H-bond from a PAA-sensing Tyr_IVS6; (iii) the ammonium group is stabilized at the focus of P-helices; and (iv) the nitrile group binds to a Ca(2+) ion coordinated by the selectivity filter glutamates in repeats III and IV. The latter feature can explain Ca(2+) potentiation of PAA action and the presence of an electronegative atom at a similar position of potent PAA analogs. Tyr substitution of a Thr in IIIS5 is known to enhance action of devapamil and verapamil. Our models predict that the para-methoxy group in ring A of devapamil and verapamil accepts an H-bond from this engineered Tyr. The model explains structure-activity relationships of PAAs, effects of LTCC mutations on PAA potency, data on PAA access to LTCC, and Ca(2+) potentiation of PAA action. Common and class-specific aspects of action of PAAs, dihydropyridines, and benzothiazepines are discussed in view of the repeat interface concept.