Concentration effect of cimetidine with POPC bilayer: a molecular dynamics simulation study

All-atom molecular dynamics simulations have been performed on cimetidine in the presence of a palmitoyloleoylphosphatidylcholine (POPC) bilayer. The free energy profile of a single cimetidine molecule passing across POPC bilayer displays a minimum at the interface of bilayer and water. Ten cimetidine molecules were inserted into POPC bilayer to obtain an 8 mol % drug model, and molecular dynamics results showed that cimetidine molecules reside at the polar region of POPC bilayer with sulphur atoms directing to the hydrophobic region. By comparing the one drug model with 8 mol % drug model, one can see that the central barrier to cross the membrane increases while the free energy in bulk water decreases, indicating that the ability of cimetidine passing across the POPC bilayer weakens at increased concentration. In addition, the free energy minimum shifts closer to the hydrophobic core. Our results indicate that with the increased drug concentration, it is more difficult for cimetidine to enter and pass across POPC bilayer.     

Acyl chain order parameter profiles in phospholipid bilayers: computation from molecular dynamics simulations and comparison with <Superscript>2</Superscript>H NMR experiments

Order parameters from deuterium NMR are often used to validate or calibrate molecular dynamics simulations. This paper gives a short overview of the literature in which experimental order parameters from ²H NMR are compared to those calculated from MD simulations. The different ways in which order parameters from experiment are used to calibrate and validate simulations are reviewed. In the second part of this review, a case study of cholesterol in a DMPC bilayer is presented. It is concluded that the agreement between experimental data and simulation is favorable in the hydrophobic region of the membrane, for both the phospholipids and cholesterol. In the interfacial region the agreement is less satisfactory, probably because of the high polarity of this region which makes the correct computation of the electrostatics more complex. Presented at the joint biannual meeting of the SFB-GEIMM-GRIP, Anglet France, 14–19 October, 2006.     

Potential of mean force of the hepatitis C virus core protein–monoclonal 19D9D6 antibody interaction

Antigen–antibody interactions are critical for understanding antigen–antibody associations in immunology. To shed further light on this question, we studied a dissociation of the 19D9D6-HCV core protein antibody complex structure. However, forced separations in single molecule experiments are difficult, and therefore molecular simulation techniques were applied in our study. The stretching, that is, the distance between the center of mass of the HCV core protein and the 19D9D6 antibody, has been studied using the potential of mean force calculations based on molecular dynamics and the explicit water model. Our simulations indicate that the 7 residues Gly70, Gly72, Gly134, Gly158, Glu219, Gln221 and Tyr314, the interaction region (antibody), and the 14 interprotein molecular hydrogen bonds might play important roles in the antigen–antibody interaction, and this finding may be useful for protein engineering of this antigen–antibody structure. In addition, the 3 residues Gly134, Gly158 and Tyr314 might be more important in the development of bioactive antibody analogs.     

Selectivity and permeation of alkali metal ions in K+-channels

Ion conduction in K⁺-channels is usually described in terms of concerted movements of K⁺ progressing in a single file through a narrow pore. Permeation is driven by an incoming ion knocking on those ions already inside the protein. A fine-tuned balance between high-affinity binding and electrostatic repulsive forces between permeant ions is needed to achieve efficient conduction. While K⁺-channels are known to be highly selective for K⁺ over Na⁺, some K⁺ channels conduct Na⁺ in the absence of K⁺. Other ions are known to permeate K⁺-channels with a more moderate preference and unusual conduction features. We describe an extensive computational study on ion conduction in K⁺-channels rendering free energy profiles for the translocation of three different alkali ions and some of their mixtures. The free energy maps for Rb⁺ translocation show at atomic level why experimental Rb⁺ conductance is slightly lower than that of K⁺. In contrast to K⁺ or Rb⁺, external Na⁺ block K⁺ currents, and the sites where Na⁺ transport is hindered are characterized. Translocation of K⁺/Na⁺ mixtures is energetically unfavorable owing to the absence of equally spaced ion-binding sites for Na⁺, excluding Na⁺ from a channel already loaded with K⁺.     

Stochastic Gating and Drug-Ribosome Interactions

Gentamicin is a potent antibiotic that is used in combination therapy for inhalation anthrax disease. The drug is also often used in therapy for methicillin-resistant Staphylococcusaureus. Gentamicin works by flipping a conformational switch on the ribosome, disrupting the reading head (i.e., 16S ribosomal decoding bases 1492-1493) used for decoding messenger RNA. We use explicit solvent all-atom molecular simulation to study the thermodynamics of the ribosomal decoding site and its interaction with gentamicin. The replica exchange molecular dynamics simulations used an aggregate sampling of 15 μs when summed over all replicas, allowing us to explicitly calculate the free-energy landscape, including a rigorous treatment of enthalpic and entropic effects. Here, we show that the decoding bases flip on a timescale faster than that of gentamicin binding, supporting a stochastic gating mechanism for antibiotic binding, rather than an induced-fit model where the bases only flip in the presence of a ligand. The study also allows us to explore the nonspecific binding landscape near the binding site and reveals that, rather than a two-state bound/unbound scenario, drug dissociation entails shuttling between many metastable local minima in the free-energy landscape. Special care is dedicated to validation of the obtained results, both by direct comparison to experiment and by estimation of simulation convergence.     

Structure, force, and energy of a double-stranded DNA oligonucleotide under tensile loads

The end-to-end stretching of a duplex DNA oligonucleotide has been studied using potential of mean force (PMF) calculations based on molecular dynamics (MD) simulations and atomic force microscopy (AFM) experiments. Near quantitative agreement between the calculations and experiments was obtained for both the extension length and forces associated with strand separation. The PMF calculations show that the oligonucleotide extends without a significant energetic barrier from a length shorter than A-DNA to a length 2.4 times the contour length of B-DNA at which the barrier to strand separation is encountered. Calculated forces associated with the barrier are 0.09±0.03 nN, based on assumptions concerning tip and thermal-activated barrier crossing contributions to the forces. Direct AFM measurements show the oligonucleotide strands separating at 2.6±0.8 contour lengths with a force of 0.13±0.05 nN. Analysis of the energies from the MD simulations during extension reveals compensation between increases in the DNA-self energy and decreases in the DNA-solvent interaction energy, allowing for the barrierless extension of DNA beyond the canonical B form. The barrier to strand separation occurs when unfavorable DNA interstrand repulsion cannot be compensated for by favorable DNA-solvent interactions. The present combination of single molecule theoretical and experimental approaches produces a comprehensive picture of the free energy surface of biological macromolecular structural transitions. Received: 2 June 1998 / Revised version: 25 January 1999 / Accepted: 11 February 1999     

Study on the mechanism of the BtuF periplasmic-binding protein for vitamin B12

BtuF is the periplasmic binding protein (PBP) that binds vitamin B₁₂ and delivers it to the periplasmic surface of the ABC transporter BtuCD. PBPs generally exhibit considerable conformational changes during ligand binding process, however, BtuF belongs to a subclass of PBPs that, doesn't show such behavior on the basis of the crystal structures. Employing steered molecular dynamics on the B₁₂-bound BtuF, we investigated the energetics and mechanism of BtuF. A potential of mean force along the postulated vitamin B₁₂ unbinding pathway was constructed through Jarzynski's equality. The large free energy differences of the postulated B₁₂ unbinding process suggests the B₁₂-bound structure is in a stable closed state and some conformation changes may be necessary to the B₁₂ unbinding. From the result of the principal component analysis, we found the BtuF-B₁₂ complex shows clear opening-closing and twisting motion tendencies which may facilitate the unbinding of B₁₂ from the binding pocket. The intrinsic flexibility of BtuF was also explored, and it's suggested the Trp44-Gln45 pair, which is situated at the mouth of the B₁₂ binding pocket, may act as a gate in the B₁₂ binding and unbinding process.     

Potential of mean force for separation of the repeating units in cellulose and hemicellulose

As a first step toward understanding the energetics of removal of cello-oligomers from the cellulose surface, we have performed umbrella sampling calculations to determine the free energy required for separation of repeating units of cellulose and hemicellulose from each other. Molecular dynamics (MD) simulations were performed for both the stacked and non-stacked arrangements of the cellobiose pair system and the xylobiose pair system. These stacked and non-stacked arrangements were taken as representative systems for the crystalline and amorphous domains in cellulose and hemicellulose. In addition, similar calculations were also carried out to determine the energetics involved in the separation of the cellobiose–xylobiose molecule pair in the non-stacked arrangement. The potential of mean force profiles exhibit a single minimum in all cases and are qualitatively similar. Our results show that the location of the minimum as well as the depth of the well can be correlated with the size of the disaccharide molecules.     

Computational studies of gramicidin permeation: An entryway sulfonate enhances cation occupancy at entry sites

The impact on the cation-transport free-energy profile of replacing the C-terminal ethanolamine in the gramicidin A channel with a taurine residue is studied using molecular dynamics simulations of gramicidin A (1JNO) embedded in a lipid bilayer (DMPC) with 1 mol/kg NaCl saline solution. The potential of mean force for ion transport is obtained by umbrella sampling. The presence of a negatively charged sulfonate group at the entrance of the gramicidin channel affects the depth and the location of the binding sites, producing a strong attraction for the cations in the bulk. The potential of mean force by the sulfonate acting directly through electrostatics and van der Waals interactions on the test ion is highly modulated by indirect effects (i.e., sulfonate effects on other components of the system that, in turn, affect the ion free-energy profile in the channel). Because the “entry” sites are located symmetrically at both entry and exit of the channel, the deeper free-energy wells should inhibit exit. Given that the channel has increased conductance experimentally, the simulation results suggest that the channel conductance is normally entry limited.     

Free energy analysis of conductivity and charge selectivity of M2GlyR-derived synthetic channels

Significant progresses have been made in the design, synthesis, modeling and in vitro testing of channel-forming peptides derived from the second transmembrane domain of the α-subunit of the glycine receptor (GlyR). The latest designs, including p22 (KKKKP ARVGL GITTV LTMTT QS), are highly soluble in water with minimal aggregation propensity and insert efficiently into cell membranes to form highly conductive ion channels. The last obstacle to a potential lead sequence for channel replacement treatment of CF patients is achieving adequate chloride selectivity. We have performed free energy simulation to analyze the conductance and charge selectivity of M2GlyR-derived synthetic channels. The results reveal that the pentameric p22 pore is non-selective. Moderate barriers for permeation of both K⁺ and Cl⁻ are dominated by the desolvation cost. Despite previous evidence suggesting a potential role of threonine side chains in anion selectivity, the hydroxyl group is not a good surrogate of water for coordinating these ions. We have also tested initial ideas of introducing additional rings of positive changes to various positions along the pore to increase anion selectivity. The results support the feasibility of achieving anion selectivity by modifying the electrostatic properties of the pore, but at the same time suggest that the peptide assembly and pore topology may also be dramatically modified, which could abolish the effects of modified electrostatics on anion selectivity. This was confirmed by subsequent two-electrode voltage clamp measurements showing that none of the tested mono-, di- and tri-Dap substituted sequences was selective. The current study thus highlights the importance of controlling channel topology besides modifying pore electrostatics for achieving anion selectivity. Several strategies are now being explored in our continued efforts to design an anion selective peptide channel with suitable biophysical, physiological and pharmacological properties as a potential treatment modality for channel replacement therapy. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

1-Alkanols and membranes: A story of attraction

Although 1-alkanols have long been known to act as penetration enhancers and anesthetics, the mode of operation is not yet understood. In this study, long-time molecular dynamics simulations have been performed to investigate the effect of 1-alkanols of various carbon chain lengths onto the structure and dynamics of dimyristoylphosphatidylcholine bilayers. The simulations were complemented by microcalorimetry, continuous bleaching and film balance experiments. In the simulations, all investigated 1-alkanols assembled inside the lipid bilayer within tens of nanoseconds. Their hydroxyl groups bound preferentially to the lipid carbonyl group and the hydrocarbon chains stretched into the hydrophobic core of the bilayer. Both molecular dynamics simulations and experiments showed that all 1-alkanols drastically affected the bilayer properties. Insertion of long-chain 1-alkanols decreased the area per lipid while increasing the thickness of the bilayer and the order of the lipids. The bilayer elasticity was reduced and the diffusive motion of the lipids within the bilayer plane was suppressed. On the other hand, integration of ethanol into the bilayer enlarged the area per lipid. The bilayer became softer and lipid diffusion was enhanced.     

The partition and transport behavior of cytotoxic ionic liquids (ILs) through the DPPC bilayer: Insights from molecular dynamics simulation

A molecular dynamics (MD) simulation with atomistic details was performed to examine the partitioning and transport behavior of moderately cytotoxic ionic liquids (ILs), namely choline bis(2-ethylhexyl) phosphate (CBEH), choline bis(2,4,4-trimethylpentyl) phosphinate (CTMP) and choline O,O-diethyl dithiophosphate (CDEP) in a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer in the fluid phase at 323?K. The structure of ILs was so selected to understand if the role of dipole and dispersion forces in the ILs distribution in the membrane can be possible. Several analyses including mass density, electrostatic potential, order parameter, diffusion coefficients and hydrogen bond formation, was carried out to determine the precise location of the anionic species inside the membrane. Moreover, the potential of the mean force (PMF) method was used to calculate free energy profile for transferring anionic species from the DPPC membrane into the bulk water. While less cytotoxic DEP is located within the bulk water, more cytotoxic TMP and BEH ILs were found to remain in the membrane and the energy barrier for crossing through the bilayer center of BEH was higher. Various ILs have no significant effect on P–N vector. The thickness of lipid bilayer decreased in all systems comprising ILs, while area per lipid increased.     

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

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.     

Surface-induced phase separation of a sphingomyelin/cholesterol/ganglioside GM1-planar bilayer on mica surfaces and microdomain molecular conformation that accelerates Aβ oligomerization

Ganglioside GM1 mediates the amyloid beta (Aβ) aggregation that is the hallmark of Alzheimer's disease (AD). To investigate how ganglioside-containing lipid bilayers interact with Aβ, we examined the interaction between Aβ40 and supported planar lipid bilayers (SPBs) on mica and SiO₂ substrates by using atomic force microscopy, fluorescence microscopy, and molecular dynamics computer simulations. These SPBs contained several compositions of sphingomyelin, cholesterol, and GM1 and were treated at physiological salt concentrations. Surprisingly high-speed Aβ aggregation of fibril formations occurred at all GM1 concentrations examined on the mica surface, but on the SiO₂ surface, only globular agglomerates formed and they formed slowly. At a GM1 concentration of 20 mol%, unique triangular regions formed on the mica surface and the rapidly formed Aβ aggregations were observed only outside these regions. We have found that some unique surface-induced phase separations are induced by the GM1 clustering effects and the strong interactions between the GM1 head group and the water layer adsorbed in the ditrigonal cavities on the mica surface. The speed of Aβ40 aggregation and the shape of the agglomerates depend on the molecular conformation of GM1, which varies depending on the substrate materials. We identified the conformation that significantly accelerates Aβ40 aggregation, and we think that the detailed knowledge about the GM1 molecular conformation obtained in this work will be useful to those investigating Aβ-GM1 interactions.     

Mechanism of tetrodotoxin block and resistance in sodium channels

Tetrodotoxin (TTX) has been used for many decades to characterize the structure and function of biological ion channels. Yet, the precise mechanism by which TTX blocks voltage-gated sodium (Na_V) channels is not fully understood. Here molecular dynamics simulations are used to elucidate how TTX blocks mammalian voltage-gated sodium (Nav) channels and why it fails to be effective for the bacterial sodium channel, Na_VAb. We find that, in Na_VAb, a sodium ion competes with TTX for the binding site at the extracellular end of the filter, thus reducing the blocking efficacy of TTX. Using a model of the skeletal muscle channel, Na_V1.4, we show that the conduction properties of the channel observed experimentally are faithfully reproduced. We find that TTX occludes the entrance of Na_V1.4 by forming a network of hydrogen-bonds at the outer lumen of the selectivity filter. The guanidine group of TTX adopts a lateral orientation, rather than pointing into the filter as proposed previously. The acidic residues just above the selectivity filter are important in stabilizing the hydrogen-bond network between TTX and Na_V1.4. The effect of two single mutations of a critical tyrosine residue in the filter of Na_V1.4 on TTX binding observed experimentally is reproduced using computational mutagenesis.     

Dynamics, Energetics, and Selectivity of the Low-K KcsA Channel Structure

Potassium channels are a diverse family of integral membrane proteins through which K⁺ can pass selectively. There is ongoing debate about the nature of conformational changes associated with the opening/closing and conductive/nonconductive states of potassium channels. The channels partly exert their function by varying their conductance through a mechanism known as C-type inactivation. Shortly after the activation of K⁺ channels, their selectivity filter stops conducting ions at a rate that depends on various stimuli. The molecular mechanism of C-type inactivation has not been fully understood yet. However, the X-ray structure of the KcsA channel obtained in the presence of low K⁺ concentration is thought to be representative of a K⁺ channel in the C-type inactivated state. Here, extensive, fully atomistic molecular dynamics and free-energy simulations of the low-K⁺ KcsA structure in an explicit lipid bilayer are performed to evaluate the stability of this structure and the selectivity of its binding sites. We find that the low-K⁺ KcsA structure is stable on the timescale of the molecular dynamics simulations performed, and that ions preferably remain in S1 and S4. In the absence of ions, the selectivity filter evolves toward an asymmetric architecture, as already observed in other computations of the high-K⁺ structure of KcsA and KirBac. The low-K⁺ KcsA structure is not permeable by Na⁺, K⁺, or Rb⁺, and the selectivity of its binding sites is different from that of the high-K⁺ structure.