Biomedical subjects

Benoît Roux

Publications and source records attributed to Benoît Roux.

At least 37 records · Page 2Linked to original sources

PubMed2005

The N-terminal end of the catalytic domain of SRC kinase Hck is a conformational switch implicated in long-range allosteric regulation.

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.

Allosteric Regulation↗

PubMed2005

Gramicidin channels.

Gramicidin channels are mini-proteins composed of two tryptophan-rich subunits. The conducting channels are formed by the transbilayer dimerization of nonconducting subunits, which are tied to the bilayer/solution interface through hydrogen bonds between the indole NH groups and the phospholipid backbone and water. The channel structure is known at atomic resolution and the channel's permeability characteristics are particularly well defined: gramicidin channels are selective for monovalent cations, with no measurable permeability to anions or polyvalent cations; ions and water move through a pore whose wall is formed by the peptide backbone; and the single-channel conductance and cation selectivity vary when the amino acid sequence is varied, even though the permeating ions make no contact with the amino acid side chains. Given the amount of experimental information that is available--for both the wild-type channels and for channels formed by amino acid-substituted gramicidin analogues--gramicidin channels provide important insights into the microphysics of ion permeation through bilayer-spanning channels. For the same reason, gramicidin channels constitute the system of choice for evaluating computational strategies for obtaining mechanistic insights into ion permeation through the complex channels formed by integral membrane proteins.

Animals↗

PubMed2005

Ion conduction and selectivity in K(+) channels.

Potassium (K(+)) channels are tetrameric membrane-spanning proteins that provide a selective pore for the conductance of K(+) across the cell membranes. These channels are most remarkable in their ability to discriminate K(+) from Na(+) by more than a thousandfold and conduct at a throughput rate near diffusion limit. The recent progress in the structural characterization of K(+) channel provides us with a unique opportunity to understand their function at the atomic level. With their ability to go beyond static structures, molecular dynamics simulations based on atomic models can play an important role in shaping our view of how ion channels carry out their function. The purpose of this review is to summarize the most important findings from experiments and computations and to highlight a number of fundamental mechanistic questions about ion conduction and selectivity that will require further work.

Animals↗

PubMed2004

Rapid intracellular TEA block of the KcsA potassium channel.

Intracellular tetraethylammonium (TEA) inhibition was studied at the single-channel level in the KcsA potassium channel reconstituted in planar lipid bilayers. TEA acts as a fast blocker (resulting in decreased current amplitude) with an affinity in the 75 mM range even at high bandwidth. Studies over a wide voltage range reveal that TEA block has a complex voltage-dependence that also depends on the ionic conditions. These observations are examined in the context of permeation models to extend our understanding of the coupling between permeant ions and TEA blockade.

Bacterial Proteins↗

PubMed2004

Structural determinants of proton blockage in aquaporins.

Aquaporins are an important class of membrane channels selective for water and linear polyols but impermeable to ions, including protons. Recent computational studies have revealed that the relay of protons through the water-conduction pathway of aquaporin channels is opposed by a substantial free energy barrier peaking at the signature NPA motifs. Here, free-energy simulations and continuum electrostatic calculations are combined to examine the nature and the magnitude of the contribution of specific structural elements to proton blockage in the bacterial glycerol uptake facilitator, GlpF. Potential of mean-force profiles for both hop and turn steps of structural diffusion in the narrow pore are obtained for artificial variants of the GlpF channel in which coulombic interactions between the pore contents and conserved residues Asn68 and Asn203 at the NPA signature motifs, Arg206 at the selectivity filter, and the peptidic backbone of the two half-helices M3 and M7, which are arranged in head-to-head fashion around the NPA motifs, are turned off selectively. A comparison of these results with electrostatic energy profiles for the translocation of a probe cation throughout the water permeation pathway indicates that the free-energy profile for proton movement inside the narrow pore is dominated by static effects arising from the distribution of charged and polar groups of the channel, whereas dielectric effects contribute primarily to opposing the access of H+ to the pore mouths (desolvation penalty). The single most effective way to abolish the free-energy gradients opposing the movement of H+ around the NPA motif is to turn off the dipole moments of helices M3 and M7. Mutation of either of the two NPA Asn residues to Asp compensates for charge-dipole and dipole-dipole effects opposing the hop and turn steps of structural diffusion, respectively, and dramatically reduces the free energy barrier of proton translocation, suggesting that these single mutants could leak protons.

Aquaporins↗

PubMed2004

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.

Bacterial Proteins↗

PubMed2004

Grand canonical Monte Carlo simulations of water in protein environments.

The grand canonical simulation algorithm is considered as a general methodology to sample the configuration of water molecules confined within protein environments. First, the probability distribution of the number of water molecules and their configuration in a region of interest for biochemical simulations, such as the active site of a protein, is derived by considering a finite subvolume in open equilibrium with a large system serving as a bulk reservoir. It is shown that the influence of the bulk reservoir can be represented as a many-body potential of mean force acting on the atoms located inside the subvolume. The grand canonical Monte Carlo (GCMC) algorithm, augmented by a number of technical advances to increase the acceptance of insertion attempts, is implemented, and tested for simple systems. In particular, the method is illustrated in the case of a pure water box with periodic boundary conditions. In addition, finite spherical systems of pure water and containing a dialanine peptide, are simulated with GCMC while the influence of the surrounding infinite bulk is incorporated using the generalized solvent boundary potential [W. Im, S. Berneche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)]. As a last illustration of water confined in the interior of a protein, the hydration of the central cavity of the KcsA potassium channel is simulated.

Algorithms↗

PubMed2004

Electrostatics of ion stabilization in a ClC chloride channel homologue from Escherichia coli.

The structural determinants of electrostatics of ion stabilization within EcClC, a ClC-type chloride channel homologue from Escherichia coli, are studied using a continuum dielectric approximation. Specifically, the ion occupancy is investigated in the wild-type protein and a mutant thereof, and the contribution to the electrostatic binding free energy of local and non-local interactions is characterized at the single-residue level. This analysis shows that, in spite of the desolvation cost and the strong ion-ion repulsion, all previously reported binding sites can be occupied simultaneously. The stabilizing effect of the protein arises from hydrogen bonding as well as from longer-range favorable interactions, such as with the strictly conserved Lys131 side-chain. The latter is involved in the stabilization of the conserved GSGIP motif that delimits two of the binding sites. Interestingly, an additional low-affinity binding site, mediated by a structurally analogous motif including the side-chain of Arg340, can be identified on the extracellular side of the permeation pathway. Finally, it is shown that, in contrast to K-channels, and in analogy to the SBP/PBP sulfate/phosphate-binding proteins, the contribution of helix macrodipoles to chloride binding in EcClC is only marginal.

Amino Acid Sequence↗

PubMed2004

Critical assessment of a proposed model of Shaker.

Detailed three-dimensional structures at atomic resolution are essential to understand how voltage-activated K(+) channels function. The X-ray crystallographic structure of the KvAP channel has offered the first view at atomic resolution of the molecular architecture of a voltage-activated K(+) channel. In the crystal, the voltage sensors are bound by monoclonal Fab fragments, which apparently induce a non-native conformation of the tetrameric channel. Thus, despite this significant advance our knowledge of the native conformation of a Kv channel in a membrane remains incomplete. Numerous results from different experimental approaches provide very specific constraints on the structure of K(+) channels in functional conformations. These results can be used to go further in trying to picture the native conformation of voltage-gated K(+) channels. However, the direct translation of all the available information into three-dimensional models is not straightforward and many questions about the structure of voltage-activated K(+) channels are still unanswered. Our aim in this review is to summarize the most important pieces of information currently available and to provide a critical assessment of the model of Shaker recently proposed by Lainé et al.

Animals↗

PubMed2004

Molecular basis of proton blockage in aquaporins.

Water transport channels in membrane proteins of the aquaporin superfamily are impermeable to ions, including H+ and OH-. We examine the molecular basis for the blockage of proton translocation through the single-file water chain in the pore of a bacterial aquaporin, GlpF. We compute the reversible thermodynamic work for the two complementary steps of the Grotthuss "hop-and-turn" relay mechanism: consecutive transfers of H+ along the hydrogen-bonded chain (hop) and conformational reorganization of the chain (turn). In the absence of H+, the strong preference for the bipolar orientation of water around the two Asn-Pro-Ala (NPA) motifs lining the pore over both unidirectional polarization states of the chain precludes the reorganization of the hydrogen-bonded network. Inversely, translocation of an excess proton in either direction is opposed by a free-energy barrier centered at the NPA region. Both hop and turn steps of proton translocation are opposed by the electrostatic field of the channel.

Aquaporins↗

PubMed2004

Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory.

Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions.

Bacterial Toxins↗

PubMed2003

Energetics of ion conduction through the gramicidin channel.

The free energy governing K(+) conduction through gramicidin A channels is characterized by using over 0.1 micros of all-atom molecular dynamics simulations with explicit solvent and membrane. The results provide encouraging agreement with experiments and insights into the permeation mechanism. The free energy surface of K(+), as a function of both axial and radial coordinates, is calculated. Correcting for simulation artifacts due to periodicity and the lack of hydrocarbon polarizability, the calculated single-channel conductance for K(+) ions is 0.8 pS, closer to experiment than any previous calculation. In addition, the estimated single ion dissociation constants are within the range of experimental determinations. The relatively small free energy barrier to ion translocation arises from a balance of large opposing contributions from protein, single-file water, bulk electrolyte, and membrane. Mean force decomposition reveals a remarkable ability of the single-file water molecules to stabilize K(+) by -40 kcal/mol, roughly half the bulk solvation free energy. The importance of the single-file water confirms the conjecture of Mackay et al. [Mackay, D. H. J., Berens, P. H., Wilson, K. R. & Hagler, A. T. (1984) Biophys. J. 46, 229-248]. Ion association with the channel involves gradual dehydration from approximately six to seven water molecules in the first shell, to just two inside the narrow pore. Ion permeation is influenced by the orientation of the single-file water column, which can present a barrier to conduction and give rise to long-range coupling of ions on either side of the pore. Small changes in the potential function, including contributions from electronic polarization, are likely to be sufficient to obtain quantitative agreement with experiments.

Amino Acid Sequence↗

PubMed2003

A microscopic view of ion conduction through the K+ channel.

Recent results from x-ray crystallography and molecular dynamics free-energy simulations have revealed the existence of a number of specific cation-binding sites disposed along the narrow pore of the K+ channel from Streptomyces lividans (KcsA), suggesting that K+ ions might literally "hop" in single file from one binding site to the next as permeation proceeds. In support of this view, it was found that the ion configurations correspond to energy wells of similar depth and that ion translocation is opposed only by small energy barriers. Although such features of the multiion potential energy surface are certainly essential for achieving a high throughput rate, diffusional and dissipative dynamical factors must also be taken into consideration to understand how rapid conduction of K+ is possible. To elucidate the mechanism of ion conduction, we established a framework theory enabling the direct simulation of nonequilibrium fluxes by extending the results of molecular dynamics over macroscopically long times. In good accord with experimental measurements, the simulated maximum conductance of the channel at saturating concentration is on the order of 550 and 360 pS for outward and inward ions flux, respectively, with a unidirectional flux-ratio exponent of 3. Analysis of the ion-conduction process reveals a lack of equivalence between the cation-binding sites in the selectivity filter.

Algorithms↗

PubMed2003

Structural basis of two-stage voltage-dependent activation in K+ channels.

The structure of the voltage sensor and the detailed physical basis of voltage-dependent activation in ion channels have not been determined. We now have identified conserved molecular rearrangements underlying two major voltage-dependent conformational changes during activation of divergent K(+) channels, ether-à-go-go (eag) and Shaker. Two conserved arginines of the S4 voltage sensor move sequentially into an extracellular gating pocket, where they interact with an acidic residue in S2. In eag, these transitions are modulated by a divalent ion that binds in the gating pocket. Conservation of key molecular details in the activation mechanism confirms that voltage sensors in divergent K(+) channels share a common structure. Molecular modeling reveals that structural constraints derived from eag and Shaker specify the unique packing arrangement of transmembrane segments S2, S3, and S4 within the voltage sensor.

Amino Acid Sequence↗

PubMed2002

Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics, Brownian dynamics, and continuum electrodiffusion theory.

Three different theoretical approaches are used and compared to refine our understanding of ion permeation through the channel formed by OmpF porin from Escherichia coli. Those approaches are all-atom molecular dynamics (MD) in which ions, solvent, and lipids are represented explicitly, Brownian dynamics (BD) in which ions are represented explicitly, while solvent and lipids are represented as featureless dielectrics, and Poisson-Nernst-Planck (PNP) electrodiffusion theory in which both solvent and local ion concentrations are represented as a continuum. First, the ability of the different theoretical approaches in reproducing the equilibrium average ion density distribution in OmpF porin bathed by a 1M KCl symmetric salt solution is examined. Under those conditions the PNP theory is equivalent to the non-linear Poisson-Boltzmann (PB) theory. Analysis shows that all the three approaches are able to capture the important electrostatic interactions between ions and the charge distribution of the channel that govern ion permeation and selectivity in OmpF. The K(+) and Cl(-) density distributions obtained from the three approaches are very consistent with one another, which suggests that a treatment on the basis of a rigid protein and continuum dielectric solvent is valid in the case of OmpF. Interestingly, both BD and continuum electrostatics reproduce the distinct left-handed twisted ion pathways for K(+) and Cl(-) extending over the length of the pore which were observed previously in MD. Equilibrium BD simulations in the grand canonical ensemble indicate that the channel is very attractive for cations, particularly at low salt concentration. On an average there is 1.55 K(+) inside the pore in 10mM KCl. Remarkably, there is still 0.17 K(+) on average inside the pore even at a concentration as low as 1microM KCl. Secondly, non-equilibrium ion flow through OmpF is calculated using BD and PNP and compared with experimental data. The channel conductance in 0.2M and 1M KCl calculated using BD is in excellent accord with the experimental data. The calculations reproduce the experimentally well-known conductance-concentration relation and also reveal an asymmetry in the channel conductance (a larger conductance is observed under a positive transmembrane potential). Calculations of the channel conductance for three mutants (R168A, R132A, and K16A) in 1M KCl suggest that the asymmetry in the channel conductance arises mostly from the permanent charge distribution of the channel rather than the shape of the pore itself. Lastly, the calculated reversal potential in a tenfold salt gradient (0.1:1M KCl) is 27.4(+/-1.3)mV (BD) and 22.1(+/-0.6)mV (PNP), in excellent accord with the experimental value of 24.3mV. Although most of the results from PNP are qualitatively reasonable, the calculated channel conductance is about 50% higher than that calculated from BD probably because of a lack of some dynamical ion-ion correlations.

Computer Simulation↗

PubMed2002

Ions and counterions in a biological channel: a molecular dynamics simulation of OmpF porin from Escherichia coli in an explicit membrane with 1 M KCl aqueous salt solution.

A 5 ns all-atom molecular dynamics trajectory of Escherichia coli OmpF porin embedded in an explicit dimyristoyl-phosphatidylcholine (DMPC) bilayer bathed by a 1 M [KCl] aqueous salt solution is generated to explore the microscopic details of the mechanism of ion permeation. The atomic model includes the OmpF trimer, 124 DMPC, 13470 water molecules as well as 231 K+ and 201 Cl-, for a total of 70,693 atoms. The structural and dynamical results are in excellent agreement with the X-ray data. The global root-mean-square deviation of the backbone atoms relative to the X-ray structure is 1.4 A. A cluster of three fully charged arginine (Arg42, Arg82, and Arg132) facing two acidic residues (Asp113 and Glu117) on L3 in the narrowest part of the aqueous pore is observed to be very stable in the crystallographic conformation. In this region of the pore, the water molecules are markedly oriented perpendicular to the channel axis due to the strong transversal electrostatic field arising from those residues. On average the size of the pore is smaller during the simulation than in the X-ray structure, undergoing small fluctuations. No large movements of loop L3 leading to a gating of the pore are observed. Remarkably, it is observed that K+ and Cl- follow two well-separated average pathways spanning over nearly 40 A along the axis of the pore. In the center of the monomer, the two screw-like pathways have a left-handed twist, undergoing a counter-clockwise rotation of 180 degrees from the extracellular vestibule to the pore periplasmic side. In the pore, the dynamical diffusion constants of the ions are reduced by about 50% relative to their value in bulk solvent. Analysis of ion solvation across the channel reveals that the contributions from the water and the protein are complementary, keeping the total solvation number of both ions nearly constant. Unsurprisingly, K+ have a higher propensity to occupy the aqueous pore than Cl-, consistent with the cation selectivity of the channel. However, further analysis suggests that ion-ion pairs play an important role. In particular, it is observed that the passage of Cl- occurs only in the presence of K+ counterions, and isolated K+ can move through the channel and permeate on their own. The presence of K+ in the pore screens the negative electrostatic potential arising from OmpF to help the translocation of Cl- by formation of ion pairs.

Computer Simulation↗

PubMed2002

Modeling the structure of agitoxin in complex with the Shaker K+ channel: a computational approach based on experimental distance restraints extracted from thermodynamic mutant cycles.

Computational methods are used to determine the three-dimensional structure of the Agitoxin (AgTx2)-Shaker complex. In a first stage, a large number of models of the complex are generated using high temperature molecular dynamics, accounting for side chain flexibility with distance restraints deduced from thermodynamic analysis of double mutant cycles. Four plausible binding mode candidates are found using this procedure. In a second stage, the quality and validity of the resulting complexes is assessed by examining the stability of the binding modes during molecular dynamics simulations with explicit water molecules and by calculating the binding free energies of mutant proteins using a continuum solvent representation and comparing with experimental data. The docking protocol and the continuum solvent model are validated using the Barstar-Barnase and the lysozyme-antibody D1.2 complexes, for which there are high-resolution structures as well as double mutant data. This combination of computational methods permits the identification of two possible structural models of AgTx2 in complex with the Shaker K+ channel, additional structural analysis providing further evidence in favor of a single model. In this final complex, the toxin is bound to the extracellular entrance of the channel along the pore axis via a combination of hydrophobic, hydrogen bonding, and electrostatic interactions. The magnitude of the buried solvent accessible area corresponding to the protein-protein contact is on the order of 1000 A with roughly similar contributions from each of the four subunits. Some side chains of the toxin adopt different conformation than in the experimental solution structure, indicating the importance of an induced-fit upon the formation of the complex. In particular, the side chain of Lys-27, a residue highly conserved among scorpion toxins, points deep into the pore with its positively charge amino group positioned at the outer binding site for K+. Specific site-directed mutagenesis experiments are suggested to verify and confirm the structure of the toxin-channel complex.

Animals↗

PubMed2002

The ionization state and the conformation of Glu-71 in the KcsA K(+) channel.

The side chain of Glu-71 of the KcsA K(+) channel, an important residue in the vicinity of the selectivity filter, was not resolved in the crystallographic structure of Doyle et al. (Doyle, D. A., J. M. Cabral, R. A. Pfuetzner, A. Kuo, J. M. Gulbis, S. L. Cohen, B. T. Chait, and R. MacKinnon. 1998. Science. 280:69-77). Its atomic coordinates are undetermined and its ionization state is unknown. For meaningful theoretical and computational studies of the KcsA K(+) channel, it is essential to address questions about the conformation and the ionization state of this residue in detail. In previous MD simulations in which the side chain of Glu-71 is protonated and forming a strong hydrogen bond with Asp-80 it was observed that the channel did not deviate significantly from the crystallographic structure (Bernèche, S., and B. Roux. 2000. Biophys. J. 78:2900-2917). In contrast, we show here that the structure of the selectivity filter of the KcsA channel is significantly disrupted when these side chains are fully ionized on each of the four monomers. To further resolve questions about the ionization state of Glu-71 we calculated the pK(a) value of this residue using molecular dynamics free energy simulations (MD/FES) with a fully flexible system including explicit solvent and membrane and finite-difference Poisson-Boltzmann (PB) continuum electrostatics. It is found that the pK(a) of Glu-71 is shifted by approximately +10 pK(a) units. These results strongly suggest that Glu-71 is protonated under normal conditions.

Bacterial Proteins↗