A Crucial Role for Side-Chain Conformation in the Versatile Charge Selectivity of Cys-Loop Receptors.

Whether synaptic transmission is excitatory or inhibitory depends, to a large extent, on whether the ion channels that open upon binding the released neurotransmitter conduct cations or anions. The mechanistic basis of the opposite charge selectivities of Cys-loop receptors has only recently begun to emerge. It is now clear that ionized side chains-whether pore-facing or buried-in the first α-helical turn of the second transmembrane segments underlie this phenomenon and that the electrostatics of backbone atoms are not critically involved. Moreover, on the basis of electrophysiological observations, it has recently been suggested that not only the sign of charged side chains but also their conformation are crucial determinants of cation-anion selectivity. To challenge these ideas with the chemical and structural rigor that electrophysiological observations naturally lack, we performed molecular dynamics, Brownian dynamics, and electrostatics calculations of ion permeation. To this end, we used structural models of the open-channel conformation of the α1 glutamate-gated Cl- channel and the α1 glycine receptor. Our results provided full support to the notion that the conformation of charged sides chains matters for charge selectivity. Indeed, whereas some rotamers of the buried arginines at position 0' conferred high selectivity for anions, others supported the permeation of cations and anions at similar rates or even allowed the faster permeation of cations. Furthermore, we found that modeling glutamates at position -1' of the anion-selective α1 glycine receptor open-state structure-instead of the five native alanines-switches charge selectivity also in a conformation-dependent manner, with some glutamate rotamers being much more effective at conferring selectivity for cations than others. Regarding pore size, we found that the mere expansion of the pore has only a minimal impact on cation-anion selectivity. Overall, these results bring to light the previously unappreciated impact of side-chain conformation on charge selectivity in Cys-loop receptors.

Gating interaction maps reveal a noncanonical electromechanical coupling mode in the Shaker K+ channel.

Membrane potential regulates the activity of voltage-dependent ion channels via specialized voltage-sensing modules, but the mechanisms involved in coupling voltage-sensor movement to pore opening remain unclear owing to a lack of resting state structures and robust methods to identify allosteric pathways. Here, using a newly developed interaction-energy analysis, we probe the interfaces of the voltage-sensing and pore modules in the Drosophila Shaker K+ channel. Our measurements reveal unexpectedly strong equilibrium gating interactions between contacts at the S4 and S5 helices in addition to those between S6 and the S4-S5 linker. Network analysis of MD trajectories shows that the voltage-sensor and pore motions are linked by two distinct pathways: a canonical pathway through the S4-S5 linker and a hitherto unknown pathway akin to rack-and-pinion coupling involving the S4 and S5 helices. Our findings highlight the central role of the S5 helix in electromechanical transduction in the voltage-gated ion channel (VGIC) superfamily.

Inference of Calmodulin's Ca2+-Dependent Free Energy Landscapes via Gaussian Mixture Model Validation.

A free energy landscape estimation method based on the well-known Gaussian mixture model (GMM) is used to compare the efficiencies of thermally enhanced sampling methods with respect to regular molecular dynamics. The simulations are carried out on two binding states of calmodulin, and the free energy estimation method is compared with other estimators using a toy model. We show that GMM with cross-validation provides a robust estimate that is not subject to overfitting. The continuous nature of Gaussians provides better estimates on sparse data than canonical histogramming. We find that diffusion properties determine the sampling method effectiveness, such that diffusion-dominated apo calmodulin is most efficiently sampled by regular molecular dynamics, while holo calmodulin, with its rugged free energy landscape, is better sampled by enhanced sampling methods.

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations.

The expansion of computational power, better parameterization of force fields, and the development of novel algorithms to enhance the sampling of the free energy landscapes of proteins have allowed molecular dynamics (MD) simulations to become an indispensable tool to understand the function of biomolecules. The temporal and spatial resolution of MD simulations allows for the study of a vast number of processes of interest. Here, we review the computational efforts to uncover the conformational free energy landscapes of a subset of membrane proteins: ion channels, transporters and G-protein coupled receptors. We focus on the various enhanced sampling techniques used to study these questions, how the conclusions come together to build a coherent picture, and the relationship between simulation outcomes and experimental observables.

Side-chain conformation at the selectivity filter shapes the permeation free-energy landscape of an ion channel.

On the basis of single-channel currents recorded from the muscle nicotinic acetylcholine receptor (AChR), we have recently hypothesized that the conformation adopted by the glutamate side chains at the first turn of the pore-lining α-helices is a key determinant of the rate of ion permeation. In this paper, we set out to test these ideas within a framework of atomic detail and stereochemical rigor by conducting all-atom molecular dynamics and Brownian dynamics simulations on an extensively validated model of the open-channel muscle AChR. Our simulations provided ample support to the notion that the different rotamers of these glutamates partition into two classes that differ markedly in their ability to catalyze ion conduction, and that the conformations of the four wild-type glutamates are such that two of them "fall" in each rotamer class. Moreover, the simulations allowed us to identify the mm (χ(1) ≅ -60°; χ(2) ≅ -60°) and tp (χ(1) ≅ 180°; χ(2) ≅ +60°) rotamers as the likely conduction-catalyzing conformations of the AChR's selectivity-filter glutamates. More generally, our work shows an example of how experimental benchmarks can guide molecular simulations into providing a type of structural and mechanistic insight that seems otherwise unattainable.