Biomedical subjects
Publications and source records attributed to Benoit Roux.
The gating mechanism of K(v) channels is not known. In this issue of Neuron, Soler-Llavina et al. present fascinating results that support the concept of relatively independent voltage-sensing modules. However, they also find that its interactions with the pore domain are rather complex, with specific S4-S5 intersubunit contacts underlying the concerted transition leading to the channel opening.
Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.
We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.
The determination of the structure of several members of the K+ channel and aquaporin family represents a unique opportunity to explain the mechanism of these biomolecular systems. With their ability to go beyond static structures, molecular dynamics simulations offer a unique route for relating functional properties to membrane channel structure. The recent progress in this area is reviewed.
Proteins, including ion channels, often are described in terms of some average structure and pictured as rigid entities immersed in a featureless solvent continuum. This simplified view, which provides for a convenient representation of the protein's overall structure, incurs the risk of deemphasizing important features underlying protein function, such as thermal fluctuations in the atom positions and the discreteness of the solvent molecules. These factors become particularly important in the case of ion movement through narrow pores, where the magnitude of the thermal fluctuations may be comparable to the ion pore atom separations, such that the strength of the ion channel interactions may vary dramatically as a function of the instantaneous configuration of the ion and the surrounding protein and pore water. Descriptions of ion permeation through narrow pores, which employ static protein structures and a macroscopic continuum dielectric solvent, thus face fundamental difficulties. We illustrate this using simple model calculations based on the gramicidin A and KcsA potassium channels, which show that thermal atomic fluctuations lead to energy profiles that vary by tens of kcal/mol. Consequently, within the framework of a rigid pore model, ion-channel energetics is extremely sensitive to the choice of experimental structure and how the space-dependent dielectric constant is assigned. Given these observations, the significance of any description based on a rigid structure appears limited. Creating a conducting channel model from one single structure requires substantial and arbitrary engineering of the model parameters, making it difficult for such approaches to contribute to our understanding of ion permeation at a microscopic level.
Two different high-resolution structures recently have been proposed for the membrane-spanning gramicidin A channel: one based on solid-state NMR experiments in oriented phospholipid bilayers (Ketchem, R. R.; Roux, B.; Cross, T. A. Structure 1997, 5, 1655-1669; Protein Data Bank, PDB:1MAG); and one based on two-dimensional NMR in detergent micelles (Townsley, L. E.; Tucker, W. A.; Sham, S.; Hinton, J. F. Biochemistry 2001, 40, 11676-11686; PDB:1JNO). Despite overall agreement, the two structures differ in peptide backbone pitch and the orientation of several side chains; in particular that of the Trp at position 9. Given the importance of the peptide backbone and Trp side chains for ion permeation, we undertook an investigation of the two structures using molecular dynamics simulation with an explicit lipid bilayer membrane, similar to the system used for the solid-state NMR experiments. Based on 0.1 micros of simulation, both backbone structures converge to a structure with 6.25 residues per turn, in agreement with X-ray scattering, and broad agreement with SS backbone NMR observables. The side chain of Trp 9 is mobile, more so than Trp 11, 13, and 15, and undergoes spontaneous transitions between the orientations in 1JNO and 1MAG. Based on empirical fitting to the NMR results, and umbrella sampling calculations, we conclude that Trp 9 spends 80% of the time in the 1JNO orientation and 20% in the 1MAG orientation. These results underscore the utility of molecular dynamics simulations in the analysis and interpretation of structural information from solid-state NMR.
A recently proposed model for voltage-dependent activation in K+ channels, largely influenced by the KvAP X-ray structure, suggests that S4 is located at the periphery of the channel and moves through the lipid bilayer upon depolarization. To investigate the physical distance between S4 and the pore domain in functional channels in a native membrane environment, we engineered pairs of cysteines, one each in S4 and the pore of Shaker channels, and identified two instances of spontaneous intersubunit disulfide bond formation, between R362C/A419C and R362C/F416C. After reduction, these cysteine pairs bound Cd2+ with high affinity, verifying that the residues are in atomic proximity. Molecular modeling based on the MthK structure revealed a single position for S4 that was consistent with our results and many other experimental constraints. The model predicts that S4 is located in the groove between pore domains from different subunits, rather than at the periphery of the protein.
Nanosecond dynamics simulations for DNA polymerase beta (pol beta)/DNA complexes with three mismatched base-pairs, namely GG, CA, or CC (primer/template) at the DNA polymerase active site, are performed to investigate the mechanism of polymerase opening and how the mispairs may affect the DNA extension step; these trajectories are compared to the behavior of a pol beta/DNA complex with the correct GC base-pair, and assessed with the aid of targeted molecular dynamics (TMD) simulations of all systems from the closed to the open enzyme state. DNA polymerase conformational changes (subdomain closing and opening) have been suggested to play a critical role in DNA synthesis fidelity, since these changes are associated with the formation of the substrate-binding pocket for the nascent base-pair. Here we observe different large C-terminal subdomain (thumb) opening motions in the simulations of pol beta with GC versus GG base-pairs. Whereas the conformation of pol beta in the former approaches the observed open state in the crystal structures, the enzyme in the latter does not. Analyses of the motions of active-site protein/DNA residues help explain these differences. Interestingly, rotation of Arg258 toward Asp192, which coordinates both active-site metal ions in the closed "active" complex, occurs rapidly in the GG simulation. We have previously suggested that this rotation is a key slow step in the closed to open transition. TMD simulations also point to a unique pathway for Arg258 rotation in the GG mispair complex. Simulations of the mismatched systems also reveal distorted geometries in the active site of all the mispair complexes examined. The hierarchy of the distortions (GG>CC>CA) parallels the experimentally deduced inability of pol beta to extend these mispairs. Such local distortions would be expected to cause inefficient DNA extension and polymerase dissociation and thereby might lead to proofreading by an extrinsic exonuclease. Thus, our studies on the dynamics of pol beta opening in mismatch systems provide structural and dynamic insights to explain experimental results regarding inefficient DNA extension following misincorporation; these details shed light on how proofreading may be invoked by the abnormal active-site geometry.
The atomic force microscope (AFM) was used to image native OmpF porin and to detect the electrostatic potential generated by the protein. To this end the OmpF porin trimers from Escherichia coli was reproducibly imaged at a lateral resolution of approximately 0.5 nm and a vertical resolution of approximately 0.1 nm at variable electrolyte concentrations of the buffer solution. At low electrolyte concentrations the charged AFM probe not only contoured structural details of the membrane protein surface but also interacted with local electrostatic potentials. Differences measured between topographs recorded at variable ionic strength allowed mapping of the electrostatic potential of OmpF porin. The potential map acquired by AFM showed qualitative agreement with continuum electrostatic calculations based on the atomic OmpF porin embedded in a lipid bilayer at the same electrolyte concentrations. Numerical simulations of the experimental conditions showed the measurements to be reproduced quantitatively when the AFM probe was included in the calculations. This method opens a novel avenue to determine the electrostatic potential of native protein surfaces at a lateral resolution better than 1 nm and a vertical resolution of approximately 0.1 nm.
Cysteine-scanning mutagenesis (SCAM) and computer-based modeling were used to investigate key structural features of the S6 transmembrane segment of the calcium-activated K(+) channel of intermediate conductance IKCa. Our SCAM results show that the interaction of [2-(trimethylammonium)ethyl] methanethiosulfonate bromide (MTSET) with cysteines engineered at positions 275, 278, and 282 leads to current inhibition. This effect was state dependent as MTSET appeared less effective at inhibiting IKCa in the closed (zero Ca(2+) conditions) than open state configuration. Our results also indicate that the last four residues in S6, from A283 to A286, are entirely exposed to water in open IKCa channels, whereas MTSET can still reach the 283C and 286C residues with IKCa maintained in a closed state configuration. Notably, the internal application of MTSET or sodium (2-sulfonatoethyl) methanethiosulfonate (MTSES) caused a strong Ca(2+)-dependent stimulation of the A283C, V285C, and A286C currents. However, in contrast to the wild-type IKCa, the MTSET-stimulated A283C and A286C currents appeared to be TEA insensitive, indicating that the MTSET binding at positions 283 and 286 impaired the access of TEA to the channel pore. Three-dimensional structural data were next generated through homology modeling using the KcsA structure as template. In accordance with the SCAM results, the three-dimensional models predict that the V275, T278, and V282 residues should be lining the channel pore. However, the pore dimensions derived for the A283-A286 region cannot account for the MTSET effect on the closed A283C and A286 mutants. Our results suggest that the S6 domain extending from V275 to V282 possesses features corresponding to the inner cavity region of KcsA, and that the COOH terminus end of S6, from A283 to A286, is more flexible than predicted on the basis of the closed KcsA crystallographic structure alone. According to this model, closure by the gate should occur at a point located between the T278 and V282 residues.
Voltage-gated K+ channels are transmembrane proteins that control and regulate the flow of K+ ions across cell membranes in response to changes in membrane potential and are essential for the propagation of action potentials in the nervous system. One of the most studied voltage-gated channels is Shaker. Available experimental results clearly provide specific constraints on the structure of the channel, even though the direct translation of the available information into 3D structures is not trivial. The goal of this work is to develop a computational approach to construct and refine 3D models of Shaker by incorporating and integrating available experimental data. Our approach is based on comparative modelization and global conformational optimization using energy restraints extracted from experimental data.