PubMed Health⌕ Search

Biomedical subjects

Benoît Roux

Publications and source records attributed to Benoît Roux.

At least 19 recordsLinked to original sources
PubMed2007

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.

Bacterial Proteins↗
PubMed2006

Mechanism of intracellular block of the KcsA K+ channel by tetrabutylammonium: insights from X-ray crystallography, electrophysiology and replica-exchange molecular dynamics simulations.

The mechanism of intracellular blockade of the KcsA potassium channel by tetrabutylammonium (TBA) is investigated through functional, structural and computational studies. Using planar-membrane electrophysiological recordings, we characterize the binding kinetics as well as the dependence on the transmembrane voltage and the concentration of the blocker. It is found that the apparent affinity of the complex is significantly greater than that of any of the eukaryotic K(+) channels studied previously, and that the off-rate increases with the applied transmembrane voltage. In addition, we report a crystal structure of the KcsA-TBA complex at 2.9 A resolution, with TBA bound inside the large hydrophobic cavity located at the center of the channel, consistent with the results of previous functional and structural studies. Of particular interest is the observation that the presence of TBA has a negligible effect on the channel structure and on the position of the potassium ions occupying the selectivity filter. Inspection of the electron density corresponding to TBA suggests that the ligand may adopt more than one conformation in the complex, though the moderate resolution of the data precludes a definitive interpretation on the basis of the crystallographic refinement methods alone. To provide a rationale for these observations, we carry out an extensive conformational sampling of an atomic model of TBA bound in the central cavity of KcsA, using the Hamiltonian replica-exchange molecular dynamics simulation method. Comparison of the simulated and experimental density maps indicates that the latter does reflect at least two distinct binding orientations of TBA. The simulations show also that the relative population of these binding modes is dependent on the ion configuration occupying the selectivity filter, thus providing a clue to the nature of the voltage-dependence of the binding kinetics.

Alanine↗
PubMed2006

Modeling the structure of the StART domains of MLN64 and StAR proteins in complex with cholesterol.

Steroidogenic acute regulatory protein-related lipid transfer (StART) domains are ubiquitously involved in intracellular lipid transport and metabolism and other cell-signaling events. In this work, we use a flexible docking algorithm, comparative modeling, and molecular dynamics (MD) simulations to generate plausible three-dimensional atomic models of the StART domains of human metastatic lymph node 64 (MLN64) and steroidogenic acute regulatory protein (StAR) proteins in complex with cholesterol. Our results show that cholesterol can adopt a similar conformation in the binding cavity in both cases and that the main contribution to the protein-ligand interaction energy derives from hydrophobic contacts. However, hydrogen-bonding and water-mediated interactions appear to be important in the fine-tuning of the binding affinity and the position of the ligand. To gain insights into the mechanism of binding, we carried out steered MD simulations in which cholesterol was gradually extracted from within the StAR model. These simulations indicate that a transient opening of loop Omega1 may be sufficient for uptake and release, and they also reveal a pathway of intermediate states involving residues known to be crucial for StAR activity. Based on these observations, we suggest specific mutagenesis targets for binding studies of cholesterol and its derivatives that could improve our understanding of the structural determinants for ligand binding by sterol carrier proteins.

Algorithms↗
PubMed2006

Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials.

The absolute (standard) binding free energy of eight FK506-related ligands to FKBP12 is calculated using free energy perturbation molecular dynamics (FEP/MD) simulations with explicit solvent. A number of features are implemented to improve the accuracy and enhance the convergence of the calculations. First, the absolute binding free energy is decomposed into sequential steps during which the ligand-surrounding interactions as well as various biasing potentials restraining the translation, orientation, and conformation of the ligand are turned "on" and "off." Second, sampling of the ligand conformation is enforced by a restraining potential based on the root mean-square deviation relative to the bound state conformation. The effect of all the restraining potentials is rigorously unbiased, and it is shown explicitly that the final results are independent of all artificial restraints. Third, the repulsive and dispersive free energy contribution arising from the Lennard-Jones interactions of the ligand with its surrounding (protein and solvent) is calculated using the Weeks-Chandler-Andersen separation. This separation also improves convergence of the FEP/MD calculations. Fourth, to decrease the computational cost, only a small number of atoms in the vicinity of the binding site are simulated explicitly, while all the influence of the remaining atoms is incorporated implicitly using the generalized solvent boundary potential (GSBP) method. With GSBP, the size of the simulated FKBP12/ligand systems is significantly reduced, from approximately 25,000 to 2500. The computations are very efficient and the statistical error is small ( approximately 1 kcal/mol). The calculated binding free energies are generally in good agreement with available experimental data and previous calculations (within approximately 2 kcal/mol). The present results indicate that a strategy based on FEP/MD simulations of a reduced GSBP atomic model sampled with conformational, translational, and orientational restraining potentials can be computationally inexpensive and accurate.

Computer Simulation↗
PubMed2006

Ion selectivity in potassium channels.

Potassium channels are tetrameric membrane-spanning proteins that provide a selective pore for the conduction of K(+) across the cell membranes. One of the main physiological functions of potassium channels is efficient and very selective transport of K(+) ions through the membrane to the cell. Classical views of ion selectivity are summarized within a historical perspective, and contrasted with the molecular dynamics (MD) simulations free energy perturbation (FEP) performed on the basis of the crystallographic structure of the KcsA phospholipid membrane. The results show that the KcsA channel does not select for K(+) ions by providing a binding site of an appropriate (fixed) cavity size. Rather, selectivity for K(+) arises directly from the intrinsic local physical properties of the ligands coordinating the cation in the binding site, and is a robust feature of a pore symmetrically lined by backbone carbonyl groups. Further analysis reveals that it is the interplay between the attractive ion-ligand (favoring smaller cation) and repulsive ligand-ligand interactions (favoring larger cations) that is the basic element governing Na(+)/K(+) selectivity in flexible protein binding sites. Because the number and the type of ligands coordinating an ion directly modulate such local interactions, this provides a potent molecular mechanism to achieve and maintain a high selectivity in protein binding sites despite a significant conformational flexibility.

Bacterial Proteins↗
PubMed2006

Molecular determinants of gating at the potassium-channel selectivity filter.

We show that in the potassium channel KcsA, proton-dependent activation is followed by an inactivation process similar to C-type inactivation, and this process is suppressed by an E71A mutation in the pore helix. EPR spectroscopy demonstrates that the inner gate opens maximally at low pH regardless of the magnitude of the single-channel-open probability, implying that stationary gating originates mostly from rearrangements at the selectivity filter. Two E71A crystal structures obtained at 2.5 A reveal large structural excursions of the selectivity filter during ion conduction and provide a glimpse of the range of conformations available to this region of the channel during gating. These data establish a mechanistic basis for the role of the selectivity filter during channel activation and inactivation.

Amino Acid Sequence↗
PubMed2006

Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field.

A polarizable potential function for the hydration of alkali and halide ions is developed on the basis of the recent SWM4-DP water model [Lamoureux, G.; MacKerell, A. D., Jr.; Roux, B. J. Chem. Phys. 2003, 119, 5185]. Induced polarization is incorporated using classical Drude oscillators that are treated as auxiliary dynamical degrees of freedom. The ions are represented as polarizable Lennard-Jones centers, whose parameters are optimized to reproduce the binding energies of gas-phase monohydrates and the hydration free energies in the bulk liquid. Systematic exploration of the parameters shows that the monohydrate binding energies can be consistent with a unique hydration free energy scale if the computed hydration free energies incorporate the contribution from the air/water interfacial electrostatic potential (-540 mV for SWM4-DP). The final model, which can satisfyingly reproduce both gas and bulk-phase properties, corresponds to an absolute scale in which the intrinsic hydration free energy of the proton is -247 kcal/mol.

Air↗
PubMed2005

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in micelles and bilayers: conformational analysis of peptide and lipids.

Molecular dynamics simulations of the influenza hemagglutinin fusion peptide in two differently sized dodecylphosphocholine micelles and a palmitoyl oleoyl phosphatidylcholine bilayer were generated to analyze the influence of the environment. Four independent trajectories (5 ns each for the bilayer, and 2 ns each for the micelles) were generated for each system. The peptide lies at the surface of the micelles, while its N-terminal region inserts deeply in the bilayer. This leads to a substantial increase of the solvation and rigidity of the peptide in micelles as compared to the bilayer. The average structures, nevertheless, are similar in all three systems and agree reasonably with micelle-based NMR structures. When in the bilayer, the peptide increases the chain gauche population and area of adjacent lipids in the same binding leaflet, while it has the opposite effect for the nearby lipids of the other leaflet. These changes, which occur spontaneously to fill voids and defects, cause a decrease in the thickness of the membrane in the neighborhood of the peptide. They would be expected to promote positive curvature, as consistent with the formation of the convex bulge, or "nipple", in the initial stage of membrane fusion. An extension of the classical surfactant theory of Israelachvili based on shapes is proposed to introduce the concept of a "dynamically induced shape" of the membrane lipids by the peptide.

Computer Simulation↗
PubMed2005

Electrostatics of the intracellular vestibule of K+ channels.

Previous calculations using continuum electrostatic calculations showed that a fully hydrated monovalent cation is electrostatically stabilized at the center of the cavity of the KcsA potassium channel. Further analysis demonstrated that this cavity stabilization was controlled by a balance between the unfavorable reaction field due to the finite size of the cavity and the favorable electrostatic field arising from the pore helices. In the present study, continuum electrostatic calculations are used to investigate how the stability of an ion in the intracellular vestibular cavity common to known potassium channels is affected as the inner channel gate opens and the cavity becomes larger and contiguous with the intracellular solution. The X-ray structure of the calcium-activated potassium channel MthK, which was crystallized in the open state, is used to construct models of the KcsA channel in the open state. It is found that, as the channel opens, the barrier at the helix bundle crossing decreases to approximately 0 kcal/mol, but that the ion in the cavity is also significantly destabilized. The results are compared and contrasted with additional calculations performed on the KvAP (voltage-activated) and KirBac1.1 (inward rectifier) channels, as well as models of the pore domain of Shaker in the open and closed state. In conclusion, electrostatic factors give rise to energetic constraints on ion permeation that have important functional consequences on the various K+ channels, and partly explain the presence or absence of charged residues near the inner vestibular entry.

Amino Acid Sequence↗
PubMed2005

The art of dissecting the function of a potassium channel.

The essential function of potassium channels is to selectively facilitate the passage of K(+) ions across the cell membrane by lowering the energy barrier that opposes this process. In this issue of Neuron, Chatelain et al. report the results from an elegant study that sheds new light on the microscopic factors that modulate the function of these important membrane channels.

Animals↗
PubMed2005

Gating charge displacement in voltage-gated ion channels involves limited transmembrane movement.

Voltage-gated ion channels are responsible for generating electrical impulses in nerves and other excitable cells. The fourth transmembrane helix (S4) in voltage-gated channels is the primary voltage-sensing unit that mediates the response to a changing membrane electric field. The molecular mechanism of voltage sensing, particularly with respect to the magnitude of the transmembrane movement of S4, remains controversial. To determine the extent of this transmembrane movement, we use fluorescent resonance energy transfer between the S4 domain and a reference point in the lipid bilayer. The lipophilic ion dipicrylamine distributes on either side of the lipid bilayer depending on the membrane potential, and is used here as a resonance-energy-transfer acceptor from donor molecules attached to several positions in the Shaker K+ channel. A voltage-driven transmembrane movement of the donor should produce a transient fluorescence change because the acceptor also translocates as a function of voltage. In Shaker K+ channels no such transient fluorescence is observed, indicating that the S4 segment does not translocate across the lipid bilayer. Based on these observations, we propose a molecular model of voltage gating that can account for the observed 13e gating charge with limited transmembrane S4 movement.

Cell Membrane↗
PubMed2005

A variable residue in the pore of Kv1 channels is critical for the high affinity of blockers from sea anemones and scorpions.

Animal toxins are associated with well defined selectivity profiles; however the molecular basis for this property is not understood. To address this issue we refined our previous three-dimensional models of the complex between the sea anemone toxin BgK and the S5-S6 region of Kv1.1 (Gilquin, B., Racape, J., Wrisch, A., Visan, V., Lecoq, A., Grissmer, S., Ménez, A., and Gasparini, S. (2002) J. Biol. Chem. 277, 37406-37413) using a docking procedure that scores and ranks the structures by comparing experimental and back-calculated values of coupling free energies DeltaDeltaGint obtained from double-mutant cycles. These models further highlight the interaction between residue 379 of Kv1.1 and the conserved dyad tyrosine residue of BgK. Because the nature of the residue at position 379 varies from one channel subtype to another, we explored how these natural mutations influence the sensitivity of Kv1 channel subtypes to BgK using binding and electrophysiology experiments. We demonstrated that mutations at this single position indeed suffice to abolish or enhance the sensitivity of Kv1 channels for BgK and other sea anemone and scorpion toxins. Altogether, our data suggest that the residue at position 379 of Kv1 channels controls the affinity of a number of blocking toxins.

Amino Acid Substitution↗
PubMed2005

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution.

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.

Computer Simulation↗
PubMed2005

Calculation of absolute protein-ligand binding free energy from computer simulations.

A general methodology for calculating the equilibrium binding constant of a flexible ligand to a protein receptor is formulated on the basis of potentials of mean force. The overall process is decomposed into several stages that can be computed separately: the free ligand in the bulk is first restrained into the conformation it adopts in the bound state, position, and orientation by applying biasing potentials, then it is translated into the binding site, where it is released completely. The conformational restraining potential is based on the root-mean-square deviation of the peptide coordinates relative to its average conformation in the bound complex. Free energy contributions from each stage are calculated by means of free energy perturbation potential of mean force techniques by using appropriate order parameters. The present approach avoids the need to decouple the ligand from its surrounding (bulk solvent and receptor protein) as is traditionally performed in the double-decoupling scheme. It is believed that the present formulation will be particularly useful when the solvation free energy of the ligand is very large. As an application, the equilibrium binding constant of the phosphotyrosine peptide pYEEI to the Src homology 2 domain of human Lck has been calculated. The results are in good agreement with experimental values.

Binding Sites↗
PubMed2005

Molecular dynamics study of hydration in ethanol-water mixtures using a polarizable force field.

The abnormal physicochemical characteristics of ethanol solvation in water are commonly attributed to the phenomenon of hydrophobic hydration. To investigate the structural organization of hydrophobic hydration in water-ethanol mixtures, we use molecular dynamics simulations based on detailed atomic models. Induced polarization is incorporated into the potential function on the basis of the classical Drude oscillator model. Water-ethanol mixtures are simulated at 11 ethanol molar fractions, from 0.05 to 0.9. Although the water and ethanol models are parametrized separately to reproduce the vaporization enthalpy, static dielectric constant, and self-diffusion constant of neat liquids at ambient conditions, they also reproduce the energetic and dynamical properties of the mixtures accurately. Furthermore, the calculated dielectric constant for the various water-alcohol mixtures is in excellent agreement with experimental data. The simulations provide a detailed structural characterization of the mixtures. A depletion of water-water hydrogen bonding in the first hydration shell of ethanol is compensated by an enhancement in the second hydration shell. The structuring effect from the second solvation shell gives rise to a net positive hydrogen-bonding excess for ethanol molar fractions up to approximately 0.5. For larger molar fractions, the second hydration shell is not sufficiently populated to overcome the net H-bond depletion from the first shell.

Alcohols↗
PubMed2005

A gate in the selectivity filter of potassium channels.

The selectivity filter of potassium channels is the structural element directly responsible for the selective and rapid conduction of K+, whereas other parts of the protein are thought to function as a molecular gate that either permits or blocks the passage of ions. However, whether the selectivity filter itself also possesses the ability to play the role of a gate is an unresolved question. Using free energy molecular dynamics simulations, it is shown that the reorientation of two peptide linkages in the selectivity filter of the KcsA K+ channel can lead to a stable nonconducting conformational state. Two microscopic factors influence the transition toward such a conformational state: the occupancy of one specific cation binding site in the selectivity filter (S2), and the strength of intersubunit interactions involving the GYG signature sequence. These results suggest that such conformational transitions occurring in the selectivity filter might be related to different K+ channel gating events, including C-type (slow) inactivation.

Ion Channel Gating↗
PubMed2005

One channel: open and closed.

Structural information about the prokaryotic KirBac3.1 inward rectifier family K(+) channel from Magnetospirillum magnetotacticum is reported. These results from two-dimensional electron cryomicroscopy (EM) shed light on the gating mechanism of members of the Kir channel family.

Cryoelectron Microscopy↗