Sergei Y. Noskov

Control of ion selectivity in potassium channels by electrostatic and dynamic properties of carbonyl ligands.

Potassium channels are essential for maintaining a normal ionic balance across cell membranes. Central to this function is the ability of such channels to support transmembrane ion conduction at nearly diffusion-limited rates while discriminating for K+ over Na+ by more than a thousand-fold. This selectivity arises because the transfer of the K+ ion into the channel pore is energetically favoured, a feature commonly attributed to a structurally precise fit between the K+ ion and carbonyl groups lining the rigid and narrow pore. But proteins are relatively flexible structures that undergo rapid thermal atomic fluctuations larger than the small difference in ionic radius between K+ and Na+. Here we present molecular dynamics simulations for the potassium channel KcsA, which show that the carbonyl groups coordinating the ion in the narrow pore are indeed very dynamic (‘liquid-like’) and that their intrinsic electrostatic properties control ion selectivity. This finding highlights the importance of the classical concept of field strength. Selectivity for K+ is seen to emerge as a robust feature of a flexible fluctuating pore lined by carbonyl groups.

Importance of Hydration and Dynamics on the Selectivity of the KcsA and NaK Channels

Fundamental concepts governing ion selectivity in narrow pores are reviewed and the microscopic factors responsible for the lack of selectivity of the NaK channel, which is structurally similar to the K+-selective KcsA channel, are elucidated on the basis of all-atom molecular dynamics free energy simulations. The results on NaK are contrasted and compared with previous studies of the KcsA channel. Analysis indicates that differences in hydration of the cation in the pore of NaK is at the origin of the lack of selectivity of NaK.

Ion selectivity in channels and transporters

A multitude of biological processes requires the participation of specific cations, such as H+, Na+, K+, Ca2+, and Mg2+. Many of these processes can take place only when proteins have the ability to discriminate between different ions with a very high fidelity. How this is possible is a fundamental

Combined Receptor and Ligand-Based Approach to the Universal Pharmacophore Model Development for Studies of Drug Blockade to the hERG1 Pore Domain

Long QT syndrome, LQTS, results in serious cardiovascular disorders, such as tachyarrhythmia and sudden cardiac death. A promiscuous binding of different drugs to the intracavitary binding site in the pore domain (PD) of human ether-a-go-go related gene (hERG) channels leads to a similar dysfunction, known as a drug-induced LQTS. Therefore, an assessment of the blocking ability for potent drugs is of great pragmatic value for molecular pharmacology and medicinal chemistry of hERGs. Thus, we attempted to create an in silico model aimed at blinded drug screening for their blocking ability to the hERG1 PD. Two distinct approaches to the drug blockage, ligand-based QSAR and receptor-based molecular docking methods, are combined for development of a universal pharmacophore model, which provides rapid assessment of drug blocking ability to the hERG1 channel. The best 3D-QSAR model (AAADR.7) from PHASE modeling was selected from a pool consisting of 44 initial candidates. The constructed model using 31 hERG blockers was validated with 9 test set compounds. The resulting model correctly predicted the pIC(50) values of test set compounds as true unknowns. To further evaluate the pharmacophore model, 14 hERG blockers with diverse hERG blocking potencies were selected from literature and they were used as additional external blind test sets. The resulting average deviation between in vitro and predicted pIC(50) values of external test set blockers is found as 0.29 suggesting that the model is able to accuretely predict the pIC(50) values as true unknowns. These pharmacophore models were merged with a previously developed atomistic receptor model for the hERG1 PD and exhibited a high consistency between ligand-based and receptor-based models. Therefore, the developed 3D-QSAR model provides a predictive tool for profiling candidate compounds before their synthesis. This model also indicated the key functional groups determining a high-affinity blockade of the hERG1 channel. To cross-validate consistency between the constructed hERG1 pore domain and the pharmacophore models, we performed docking studies using the homology model of hERG1. To understand how polar or nonpolar moieties of inhibitors stimulate channel inhibition, critical amino acid replacement (i.e., T623, S624, S649, Y652 and F656) at the hERG cavity was examined by in silico mutagenesis. The average docking score differences between wild type and mutated hERG channels was found to have the following order: F656A > Y652A > S624A > T623A > S649A. These results are in agreement with experimental data.

Evidence for a third sodium-binding site in glutamate transporters suggests an ion/substrate coupling model.

Excitatory amino acid transporters (EAATs) remove glutamate from synapses. They maintain an efficient synaptic transmission and prevent glutamate from reaching neurotoxic levels. Glutamate transporters couple the uptake of one glutamate to the cotransport of three sodium ions and one proton and the countertransport of one potassium ion. The molecular mechanism for this coupled uptake of glutamate and its co- and counter-transported ions is not known. In a crystal structure of the bacterial glutamate transporter homolog, GltPh, only two cations are bound to the transporter, and there is no indication of the location of the third sodium site. In experiments using voltage clamp fluorometry and simulations based on molecular dynamics combined with grand canonical Monte Carlo and free energy simulations performed on different isoforms of GltPh as well on a homology model of EAAT3, we sought to locate the third sodium-binding site in EAAT3. Both experiments and computer simulations suggest that T370 and N451 (T314 and N401 in GltPh) form part of the third sodium-binding site. Interestingly, the sodium bound at T370 forms part of the binding site for the amino acid substrate, perhaps explaining both the strict coupling of sodium transport to uptake of glutamate and the ion selectivity of the affinity for the transported amino acid in EAATs.

Two mechanisms of ion selectivity in protein binding sites

A theoretical framework is presented to clarify the molecular determinants of ion selectivity in protein binding sites. The relative free energy of a bound ion is expressed in terms of the main coordinating ligands coupled to an effective potential of mean force representing the influence of the rest of the protein. The latter is separated into two main contributions. The first includes all the forces keeping the ion and the coordinating ligands confined to a microscopic subvolume but does not prevent the ligands from adapting to a smaller or larger ion. The second regroups all the remaining forces that control the precise geometry of the coordinating ligands best adapted to a given ion. The theoretical framework makes it possible to delineate two important limiting cases. In the limit where the geometric forces are dominant (rigid binding site), ion selectivity is controlled by the ion-ligand interactions within the matching cavity size according to the familiar “snug-fit” mechanism of host-guest chemistry. In the limit where the geometric forces are negligible, the ion and ligands behave as a “confined microdroplet” that is free to fluctuate and adapt to ions of different sizes. In this case, ion selectivity is set by the interplay between ion-ligand and ligand-ligand interactions and is controlled by the number and the chemical type of ion-coordinating ligands. The framework is illustrated by considering the ion-selective binding sites in the KcsA channel and the LeuT transporter.

A Three-Dimensional Model of a Group II Intron RNA and Its Interaction with the Intron-Encoded Reverse Transcriptase

Group II introns are self-splicing ribozymes believed to be the ancestors of spliceosomal introns. Many group II introns encode reverse transcriptases that promote both RNA splicing and intron mobility to new genomic sites. Here we used a circular permutation and crosslinking method to establish 16 intramolecular distance relationships within the mobile Lactococcus lactis Ll.LtrB-DeltaORF intron. Using these new constraints together with 13 established tertiary interactions and eight published crosslinks, we modeled a complete three-dimensional structure of the intron. We also used the circular permutation strategy to map RNA-protein interaction sites through fluorescence quenching and crosslinking assays. Our model provides a comprehensive structural framework for understanding the function of group II ribozymes, their natural structural variations, and the mechanisms by which the intron-encoded protein promotes RNA splicing and intron mobility. The model also suggests an arrangement of active site elements that may be conserved in the spliceosome.

Modeling of open, closed, and open-inactivated states of the hERG1 channel: structural mechanisms of the state-dependent drug binding.

The human ether-a-go-go related gene 1 (hERG1) K ion channel is a key element for the rapid component of the delayed rectified potassium current in cardiac myocytes. Since there are no crystal structures for hERG channels, creation and validation of its reliable atomistic models have been key targets in molecular cardiology for the past decade. In this study, we developed and vigorously validated models for open, closed, and open-inactivated states of hERG1 using a multistep protocol. The conserved elements were derived using multiple-template homology modeling utilizing available structures for Kv1.2, Kv1.2/2.1 chimera, and KcsA channels. Then missing elements were modeled with the ROSETTA De Novo protein-designing suite and further refined with all-atom molecular dynamics simulations. The final ensemble of models was evaluated for consistency to the reported experimental data from biochemical, biophysical, and electrophysiological studies. The closed state models were cross-validated against available experimental data on toxin footprinting with protein-protein docking using hERG state-selective toxin BeKm-1. Poisson-Boltzmann calculations were performed to determine gating charge and compare it to electrophysiological measurements. The validated structures offered us a unique chance to assess molecular mechanisms of state-dependent drug binding in three different states of the channel.

Representation of Ion–Protein Interactions Using the Drude Polarizable Force-Field

Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K+, Na+, Ca2+, and Cl–) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca2+ binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins.

Molecular Mechanism of Ion-Ion and Ion-Substrate Coupling in the Na+-Dependent Leucine Transporter LeuT

Ion-coupled transport of neurotransmitter molecules by neurotransmitter:sodium symporters (NSS) play an important role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling and the corresponding ion-substrate stoichiometry in NSS transporters has yet to be understood. The recent determination of a high-resolution structure for a bacterial homolog of Na(+)/Cl(-)-dependent neurotransmitter transporters, LeuT, offers a unique opportunity to analyze the functional roles of the multi-ion binding sites within the binding pocket. The binding pocket of LeuT contains two metal binding sites. The first ion in site NA1 is directly coupled to the bound substrate (Leu) with the second ion in the neighboring site (NA2) only approximately 7 A away. Extensive, fully atomistic, molecular dynamics, and free energy simulations of LeuT in an explicit lipid bilayer are performed to evaluate substrate-binding affinity as a function of the ion load (single versus double occupancy) and occupancy by specific monovalent cations. It was shown that double ion occupancy of the binding pocket is required to ensure substrate coupling to Na(+) and not to Li(+) or K(+) cations. Furthermore, it was found that presence of the ion in site NA2 is required for structural stability of the binding pocket as well as amplified selectivity for Na(+) in the case of double ion occupancy.