Biomedical subjects

Christophe Chipot

Publications and source records attributed to Christophe Chipot.

17 recordsLinked to original sources

PubMed2006

Conformational equilibrium in alanine-rich peptides probed by reversible stretching simulations.

Reversible stretching of the alanine-rich peptide 3K (Proc. Natl. Acad. Sci. USA 1989, 86, 5286-5290) and its analogue MW (Nature 1992, 359, 653-655) is examined using molecular dynamics simulations in explicit water. In both cases, sampling of the extension pathway is obtained on the 10 ns time scale by applying an adaptive biasing force. The free energy profile reveals a single minimum associated with a contiguous alpha-helix. Short 3(10)-helical motifs are observed in folded as well as extended conformations, in accordance with their proposed role as folding intermediates. The native 3(10)-helical content of both peptides is found, however, to be no higher than a few percent. Difficulties in both the definition and the detection of secondary structure motifs, most notably in relation to bifurcated hydrogen bonds, are proposed to account for the discrepancy between 3(10)-helical propensities reported by several authors, based on experimental and computational results.

Alanine↗

PubMed2006

Molecular dynamics study of the inclusion of cholesterol into cyclodextrins.

The interaction of three cyclodextrins (CDs), viz. beta-CD, heptakis (2,6-di-O-methyl)-beta-CD (DM-beta-CD), and 2-hydroxypropyl-beta-CD (HP-beta-CD), with cholesterol was investigated using molecular dynamics (MD) simulations. The free energy along the reaction pathway delineating the inclusion of cholesterol into each CD was computed using the adaptive biasing force method. The association constant and the corresponding association free energy were derived by integrating the potential of mean force (PMF) over a representative ordering parameter. The results show that the free energy profiles possess two local minima corresponding to roughly equally probable binding modes. Among the three CDs, DM-beta-CD exhibits the highest propensity to associate with cholesterol. Ranking for binding cholesterol, viz. DM-beta-CD > HP-beta-CD > beta-CD, agrees nicely with experiment. Partitioning of the PMF into free energy components illuminates that entering of cholesterol into the CD cavity is driven mainly by electrostatic interactions, whereas deeper inclusion results from van der Waals forces and solvation effects. Additional MD simulations were performed to investigate the structural stability of the host-guest complexes near the free energy minima. The present results demonstrate that association of cholesterol and CDs follows two possible binding modes. Although the latter are thermodynamically favorable for all CDs, one of the two inclusion complexes appears to be preferred kinetically in the case of DM-beta-CD.

Cholesterol↗

PubMed2006

NMR structure and molecular dynamics of the in-plane membrane anchor of nonstructural protein 5A from bovine viral diarrhea virus.

Hepatitis C virus (HCV) nonstructural protein 5A (NS5A) is a monotopic membrane protein anchored to the membrane by an N-terminal in-plane amphipathic alpha-helix. This membrane anchor is essential for the assembly of a functional viral replication complex. Although amino acid sequences differ considerably, putative membrane anchors with amphipathic features were predicted in NS5A from related Flaviviridae family members, in particular bovine viral diarrhea virus (BVDV), the prototype representative of the genus Pestivirus. We report here the NMR structure of the membrane anchor 1-28 of NS5A from BVDV in the presence of different membrane mimetic media. This anchor includes a long amphipathic alpha-helix of 21 residues interacting in-plane with the membrane interface and including a putative flexible region. Molecular dynamic simulation at a water-dodecane interface used to mimic the surface separating a lipid bilayer and an aqueous medium demonstrated the stability of the helix orientation and the location at the hydrophobic-hydrophilic interface. The flexible region of the helix appears to be required to allow the most favorable interaction of hydrophobic and hydrophilic side chain residues with their respective environment at the membrane interface. Despite the lack of amino acid sequence similarity, this amphipathic helix shares common structural features with that of the HCV counterpart, including a stable, hydrophobic N-terminal segment separated from the more hydrophilic C-terminal segment by a local, flexible region. These structural conservations point toward conserved roles of the N-terminal in-plane membrane anchors of NS5A in replication complex formation of HCV, BVDV, and other related viruses.

Amino Acid Sequence↗

PubMed2006

Interaction of a peptide nanotube with a water-membrane interface.

Inserting peptide nanotubes into lipid bilayers modulates the permeability properties of the cell wall, thus conferring potential bacteriocidal capability. Interaction of a peptide nanotube formed by eight cyclo[RRKWLWLW] subunits with the surface of a hydrated dimyristoylphosphatidylcholine bilayer is investigated using molecular dynamics simulations. The present sequence of alternated D-L-alpha-amino acids has been shown to yield remarkable antibacterial in vitro activity, and the chosen topoisomer corresponds to the optimum amphipathy of the tubular structure, whereby non-polar and charged side chains are segregated by the aqueous interface. The cohesion of the nanotube is ensured by a scaffold of intermolecular hydrogen bonds between adjacent cyclic peptides, supplemented by favorable like-charged contacts of arginine side chains. It is further reinforced by interactions of charged residues with the lipid head groups and of non-polar residues with the lipid acyl chains. The simulation reveals a partial breaking of the synthetic channel accompanying its early insertion into the lipid bilayer. The latter opens new questions about how peptide nanotubes permeate the membrane, in particular whether or not (i) self-assembly precedes partitioning and (ii) translocation occurs with the complete tubular structure.

Computer Simulation↗

PubMed2005

Exploring the free-energy landscape of a short peptide using an average force.

The reversible folding of deca-alanine is chosen as a test case for characterizing a method that uses an adaptive biasing force (ABF) to escape from the minima and overcome the barriers of the free-energy landscape. This approach relies on the continuous estimation of a biasing force that yields a Hamiltonian in which no average force is exerted along the ordering parameter xi. Optimizing the parameters that control how the ABF is applied, the method is shown to be extremely effective when a nonequivocal ordering parameter can be defined to explore the folding pathway of the peptide. Starting from a beta-turn motif and restraining xi to a region of the conformational space that extends from the alpha-helical state to an ensemble of extended structures, the ABF scheme is successful in folding the peptide chain into a compact alpha helix. Sampling of this conformation is, however, marginal when the range of xi values embraces arrangements of greater compactness, hence demonstrating the inherent limitations of free-energy methods when ambiguous ordering parameters are utilized.

Alanine↗

PubMed2005

Can free energy calculations be fast and accurate at the same time? Binding of low-affinity, non-peptide inhibitors to the SH2 domain of the src protein.

The usefulness of free-energy calculations in non-academic environments, in general, and in the pharmaceutical industry, in particular, is a long-time debated issue, often considered from the angle of cost/performance criteria. In the context of the rational drug design of low-affinity, non-peptide inhibitors to the SH2 domain of the (pp60)src tyrosine kinase, the continuing difficulties encountered in an attempt to obtain accurate free-energy estimates are addressed. free-energy calculations can provide a convincing answer, assuming that two key-requirements are fulfilled: (i) thorough sampling of the configurational space is necessary to minimize the statistical error, hence raising the question: to which extent can we sacrifice the computational effort, yet without jeopardizing the precision of the free-energy calculation? (ii) the sensitivity of binding free-energies to the parameters utilized imposes an appropriate parametrization of the potential energy function, especially for non-peptide molecules that are usually poorly described by multipurpose macromolecular force fields. Employing the free-energy perturbation method, accurate ranking, within +/-0.7 kcal/mol, is obtained in the case of four non-peptide mimes of a sequence recognized by the (pp60)src SH2 domain.

Binding Sites↗

PubMed2005

Probing a model of a GPCR/ligand complex in an explicit membrane environment: the human cholecystokinin-1 receptor.

A three-dimensional model structure of a complex formed by a G-protein-coupled receptor (GPCR) and an agonist ligand is probed and refined using molecular-dynamics simulations and free energy calculations in a realistic environment. The model of the human receptor of cholecystokinin associated to agonist ligand CCK9 was obtained from a synergistic procedure combining site-directed mutagenesis experiments and in silico modeling. The 31-ns molecular-dynamics simulation in an explicit membrane environment indicates that both the structure of the receptor and its interactions with the ligand are robust. Whereas the secondary structure of the alpha-helix bundle is well preserved, the region of the intracellular loops exhibits a significant flexibility likely to be ascribed to the absence of G-protein subunits in the model. New insight into the structural features of the binding pocket is gained, in particular, the interplay of the ligand with both the receptor and internal water molecules. Water-mediated interactions are shown to participate in the binding, hence, suggesting additional site-directed mutagenesis experiments. Accurate free energy calculations on mutated ligands provide differences in the receptor-ligand binding affinity, thus offering a direct, quantitative comparison to experiment. We propose that this detailed consistency-checking procedure be used as a routine refinement step of in vacuo GPCR models, before further investigation and application to structure-based drug design.

Cell Membrane↗

PubMed2005

Insights into the recognition and association of transmembrane alpha-helices. The free energy of alpha-helix dimerization in glycophorin A.

The free energy of alpha-helix dimerization of the transmembrane (TM) region of glycophorin A was estimated from a 125-ns molecular dynamics (MD) simulation in a membrane mimetic. The free energy profile was obtained by allowing the TM helical segments to diffuse reversibly along the reaction pathway. Partition of the potential of mean force into free energy components illuminates the critical steps of alpha-helix recognition and association. At large separations, the TM segments are pushed together by the solvent, allowing initial, but not necessarily native, interhelical interactions to occur. This early recognition stage precedes the formation of native contacts, which is accompanied by a tilt of the helices, characteristic of the dimeric structure. This step is primarily driven by the van der Waals helix-helix interactions. Free energy perturbation calculations of the L75A and I76A point mutations reveal a disruption in helix-helix association due to a loss of favorable dispersion interactions. Additional MD simulations of the native TM dimer and of a single alpha-helix confirm that, prior to association, individual alpha-helices are independently stable, in agreement with the "two-stage" model of integral membrane protein folding.

Amino Acid Sequence↗

PubMed2005

Scalable molecular dynamics with NAMD.

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Algorithms↗

PubMed2004

Overcoming free energy barriers using unconstrained molecular dynamics simulations.

Association of unconstrained molecular dynamics (MD) and the formalisms of thermodynamic integration and average force [Darve and Pohorille, J. Chem. Phys. 115, 9169 (2001)] have been employed to determine potentials of mean force. When implemented in a general MD code, the additional computational effort, compared to other standard, unconstrained simulations, is marginal. The force acting along a chosen reaction coordinate xi is estimated from the individual forces exerted on the chemical system and accumulated as the simulation progresses. The estimated free energy derivative computed for small intervals of xi is canceled by an adaptive bias to overcome the barriers of the free energy landscape. Evolution of the system along the reaction coordinate is, thus, limited by its sole self-diffusion properties. The illustrative examples of the reversible unfolding of deca-L-alanine, the association of acetate and guanidinium ions in water, the dimerization of methane in water, and its transfer across the water liquid-vapor interface are examined to probe the efficiency of the method.

Journal Article↗

PubMed2004

Coupled motions between pore and voltage-sensor domains: a model for Shaker B, a voltage-gated potassium channel.

A high-resolution crystal structure of KvAP, an archeabacterial voltage-gated potassium (Kv) channel, complexed with a monoclonal Fab fragment has been recently determined. Based on this structure, a mechanism for the activation (opening) of Kv channels has been put forward. This mechanism has since been criticized, suggesting that the resolved structure is not representative of the family of voltage-gated potassium channels. Here, we propose a model of the transmembrane domain of Shaker B, a well-characterized Kv channel, built by homology modeling and docking calculations. In this model, the positively charged S4 helices are oriented perpendicular to the membrane and localized in the groove between segments S5 and S6 of adjacent subunits. The structure and the dynamics of the full atomistic model embedded in a hydrated lipid bilayer were investigated by means of two large-scale molecular dynamics simulations under transmembrane-voltage conditions known to induce, respectively, the resting state (closed) and the activation (opening) of voltage-gated channels. Upon activation, the model undergoes conformational changes that lead to an increase of the hydration of the charged S4 helices, correlated with an upward translation and a tilting of the latter, concurrently with movements of the S5 helices and the activation gate. Although small, these conformational changes ultimately result in an alteration of the ion-conduction pathway. Our findings support the transporter model devised by Bezanilla and collaborators, and further underline the crucial role played by internal hydration in the activation of the channel.

Computer Simulation↗

PubMed2003

Non-arrhenius behavior in the unfolding of a short, hydrophobic alpha-helix. Complementarity of molecular dynamics and lattice model simulations.

The unfolding of the last, C-terminal residue of AcNH(2)-(l-Leu)(11)-NHMe in its alpha-helical form has been investigated by measuring the variation of free energy involved in the alpha(R) to beta conformational transition. These calculations were performed using large-scale molecular dynamics simulations in conjunction with the umbrella sampling method. For different temperatures ranging from 280 to 370 K, the free energy of activation was estimated. Concurrently, unfolding simulations of a homopolypeptide formed by twelve hydrophobic residues were carried out, employing a three-dimensional lattice model description of the peptide, with a temperature-dependent interaction potential. Using a Monte Carlo approach, the lowest free energy conformation, an analogue of a right-handed alpha-helix, was determined in the region where the peptide chain is well ordered. The free energy barrier separating this state from a distinct, compact conformation, analogue to a beta-strand, was determined over a large enough range of temperatures. The results of these molecular dynamics and lattice model simulations are consistent and indicate that the kinetics of the unfolding of a hydrophobic peptide exhibits a non-Arrhenius behavior closely related to the temperature dependence of the hydrophobic effect. These results further illuminate the necessity to include a temperature dependence in potential energy functions designed for coarse-grained models of proteins.

Computer Simulation↗

PubMed2003

Rational determination of charge distributions for free energy calculations.

Point charges derived from RHF/6-31G* electrostatic potentials are attractive because they tend to exaggerate the polarity of solvated molecules, thereby compensating in an average fashion missing induction effects. In the context of free energy calculations, wherein the molecule is transferred from a polar environment to a nonpolar one, we propose a more rational approach based on a self-consistent reaction field computation at a higher level of theory, supplemented by an estimation of the corresponding distortion energy to account for the change of polarity of the surroundings. Application of this method to the test cases acetamide, acetic acid, methyl acetate and phenol, using multinanosecond molecular dynamics/"umbrella sampling" simulations, yields consistent hydration free energies in reasonably good agreement with experiment.

Journal Article↗

PubMed2003

OPEP: a tool for the optimal partitioning of electric properties.

OPEP is a suite of FORTRAN programs targeted at the optimal partitioning of molecular electric properties. It includes an interactive module for the construction of Cartesian grids of points, on which either the molecular electrostatic potential or the induction energy is mapped. The generation of distributed multipoles and polarizabilities is achieved using either the formalism of the normal equations of the least-squares problem, which restates the fitting procedure in terms of simple matrix operations, or a statistical approach, which provides a pictorial description of the distributed models of multipoles and polarizabilities, thereby allowing the pinpointing of pathological cases. Molecular symmetry is accounted for by means of local atomic frames, which are generated in an automated fashion. A Tcl/Tk graphical user interface wraps the suite of programs, thereby making OPEP a user-friendly package for building models of distributed multipoles and polarizabilities. OPEP is an open-source suite of programs distributed free of charge under the GNU general public license (GPL) at http://www.lctn.uhp-nancy.fr/Opep.

Journal Article↗

PubMed2003

Molecular dynamics investigation of an oriented cyclic peptide nanotube in DMPC bilayers.

Nanotubes resulting from the self-assembly of cyclic peptides formed by eight alpha-amino acids and inserted into lipid bilayers have been shown to function as synthetic, integral transmembrane channels. A nanotube consisting of eight cyclo[(L-Trp-D-Leu)(3)-L-Gln-D-Leu] subunits, organized in an antiparallel, beta-sheetlike channel embedded in a hydrated dimyristoylphosphatidylcholine bilayer was investigated in an 8-ns molecular dynamics trajectory. This large-scale statistical simulation brings to light not only the atomic-level structural features of the synthetic channel, but also its dynamical properties. Overall, the nanotube conserves its hollow tubular structure. The calculation reproduces the tilt of the channel with respect to the normal of the bilayer, in reasonable agreement with experiment. The results show a dislocation of the nanotube indicative of a possible disassembly process that may influence the channel conduction. The dynamics of the water in the hollow tubular structure has been characterized, and the conductance of the channel has been estimated. Transport properties of the peptide nanotube are discussed in comparison with other transporters.

Computer Simulation↗

PubMed2003

Membrane peptides and their role in protobiological evolution.

How simple membrane peptides performed such essential protocellular functions as transport of ions and organic matter across membranes separating the interior of the cell from the environment, capture and utilization of energy, and transduction of environmental signals, is a key question in protobiological evolution. On the basis of detailed, molecular-level computer simulations we explain how these peptides fold at water-membrane interfaces, insert into membranes, self-assemble into higher-order structures and acquire functions. We have investigated the interfacial behavior and folding of several peptides built of leucine and glutamine residues and have demonstrated that many of them tend to adopt ordered structures. Further, we have studied the insertion of an alpha-helical peptide containing leucine (L) and serine (S) of the form (LSLLLSL)3 into a model membrane. The transmembrane state is metastable, and approximately 15 kcal mol(-1) is required to insert the peptide into the membrane. Investigations of dimers formed by (LSLLLSL)3 and glycophorin A demonstrate how the favorable free energy of helix association can offset the unfavorable free energy of insertion, leading to self-assembly of peptide helices in the membrane. An example of a self-assembled structure is the tetrameric transmembrane pore of the influenza virus M2 protein, which is an efficient and selective voltage-gated proton channel. Our simulations explain the gating mechanism and provide guidelines how to re-engineer the channel to act as a simple proton pump. In general, emergence of integral membrane proteins appears to be quite feasible and may be easier to envision than the emergence of water-soluble proteins.

Amino Acids↗

PubMed2002

Rational determination of transfer free energies of small drugs across the water-oil interface.

The application of statistical simulations to the estimation of transfer free energies of pharmacologically relevant organic molecules is reported. Large-scale molecular dynamics simulations have been carried out on a series of four solutes, viz. antipyrine, caffeine, ganciclovir, and alpha-D-glucose, at the water-dodecane interface as a model of a biological water-membrane interfacial system. Agreement with experimentally determined partition coefficients is remarkable, demonstrating that free energy calculations, when executed with appropriate protocols, have reached the maturity to predict thermodynamic quantities of interest to the pharmaceutical world. The computational effort that warrants accurate, converged free energies remains, however, in large measure, incompatible with the high-throughput exploration of large sets of pharmacologically active drugs sought by industrial settings. Compared to the cost-effective, fast estimation of simple partition coefficients, the present free energy calculations, nevertheless, offer a far more detailed information about the underlying energetics of the system when the solute is translocated across the water-dodecane interface, which can be valuable in the context of de novo drug design.

Alkanes↗