Catalytic roles of active-site amino acid residues of coenzyme B12-dependent diol dehydratase: protonation state of histidine and pull effect of glutamate

The hydrogen abstraction and the OH migration processes catalyzed by diol dehydratase are discussed by means of a quantum mechanical/molecular mechanical method. To evaluate the push effect of His143 and the pull effect of Glu170, we considered three kinds of whole-enzyme model, the protonated and two unprotonated His143 models. A calculated activation energy for the hydrogen abstraction by the adenosyl radical is 15.6 (13.6) kcal/mol in the protonated (unprotonated) His143 model. QM/MM calculational results show that the mechanism of the OH migration is significantly changed by the protonation of His143. In the protonated His143 model, the OH group migration triggered by the full proton donation from the imidazolium to the migrating OH group occurs by a stepwise OH abstraction/re-addition process in which the water production reduces the barrier for the C-O bond cleavage. On the other hand, the OH migration in the unprotonated His143 model proceeds in a concerted manner, as we previously proposed using a simple model including only K+ ion and substrate. The latter mechanism seems to be kinetically more favorable from the calculated energy profiles and is consistent with experimental results. The activation barrier of the OH group migration step is only 1.6 kcal/mol reduced by the hydrogen-bonding interaction between the O2 of the substrate and unprotonated His143. Thus, it is predicted that His143 is not protonated, and therefore the main active-site amino acid residue that lowers the energy of the transition state for the OH group migration is determined to be Glu170.     

Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.     

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.     

Protonation of the acidic residues in the transmembrane cation-binding sites of the ca(2+) pump

The ionization states of the acidic residues around the Ca2+-binding sites of sarcoplasmic reticulum Ca2+ ATPase are studied by continuum electrostatic calculations and all-atom molecular dynamics simulations with explicit solvent and phospholipids. The two methods consistently indicate that Glu58 and Glu908 are protonated at neutral pH. The Ca2+ coordination and the H-bonds formed by the protonation of Glu58 and Glu908 are stable in an MD simulation, whereas the H-bonds are disrupted and the Ca2+ coordination geometry is severely altered in another simulation treating these residues unprotonated. The results clearly indicate that the H-bonds formed by protonation of Glu58 and Glu908 provide extra stability for the Ca2+-binding sites of Ca2+ ATPase.     

Simulation of pH-Dependent Conformational Transitions in Membrane Proteins: The CLC-ec1 Cl−/H+ Antiporter

Intracellular transport of chloride by members of the CLC transporter family involves a coupled exchange between a Cl⁻ anion and a proton (H⁺), which makes the transport function dependent on ambient pH. Transport activity peaks at pH 4.5 and stalls at neutral pH. However, a structure of the WT protein at acidic pH is not available, making it difficult to assess the global conformational rearrangements that support a pH-dependent gating mechanism. To enable modeling of the CLC-ec1 dimer at acidic pH, we have applied molecular dynamics simulations (MD) featuring a new force field modification scheme—termed an Equilibrium constant pH approach (ECpH). The ECpH method utilizes linear interpolation between the force field parameters of protonated and deprotonated states of titratable residues to achieve a representation of pH-dependence in a narrow range of physiological pH values. Simulations of the CLC-ec1 dimer at neutral and acidic pH comparing ECpH-MD to canonical MD, in which the pH-dependent protonation is represented by a binary scheme, substantiates the better agreement of the conformational changes and the final model with experimental data from NMR, cross-link and AFM studies, and reveals structural elements that support the gate-opening at pH 4.5, including the key glutamates Glu_in and Glu_ex.     

A gate mechanism indicated in the selectivity filter of the potassium channel KscA

Classical molecular dynamics (MD) and non-equilibrium steered molecular dynamics (SMD) simulations were performed on the molecular structure of the potassium channel KcsA using the GROMOS 87 force fields. Our simulations focused on mechanistic and dynamic properties of the permeation of potassium ions through the selectivity filter of the channel. According to the SMD simulations a concerted movement of ions inside the selectivity filter from the cavity to extracellular side depends on the conformation of the peptide linkage between Val76 and Gly77 residues in one subunit of the channel. In SMD simulations, if the carbonyl oxygen of Val76 is positioned toward the ion bound at the S3 site (gate-opened conformation) the net flux of ions through the filter is observed. When the carbonyl oxygen leaped out from the filter (gate-closed conformation), ions were blocked at the S3 site and no flux occurred. A reorientation of the Thr75-Val76 linkage indicated by the CHARMM-based MD simulations performed Berneche and Roux [(2005) Structure 13:591–600; (2000) Biophys J 78:2900–2917] as a concomitant process of the Val76-Gly77 conformational interconversion was not observed in our GROMOS-based MD simulations.     

Solid-state 19F NMR spectroscopy reveals that Trp41 participates in the gating mechanism of the M2 proton channel of influenza A virus

The integral membrane protein M2 of influenza A virus assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH. The side chain of His37 in the transmembrane alpha-helix is known to play an important role in the pH activation of the proton channel. It has also been suggested that Trp41, which is located in an adjacent turn of the helix, forms part of the gating mechanism. Here, a synthetic 25-residue peptide containing the M2 transmembrane domain was labeled with 6F-Trp41 and studied in lipid membranes by solid-state 19F NMR. We monitored the pH-dependent differences in the 19F dipolar couplings and motionally narrowed chemical shift anisotropies of this 6F-Trp41 residue, and we discuss the pH activation mechanism of the H+ channel. At pH 8.0, the structural parameters implicate an inactivated state, while at pH 5.3 the tryptophan conformation represents the activated state. With the aid of COSMOS force field simulations, we have obtained new side-chain torsion angles for Trp41 in the inactivated state (chi1 = -100 degrees +/- 10 degrees , chi2 = +110 degrees +/- 10 degrees ), and we predict a most probable activated state with chi1 = -50 degrees +/- 10 degrees and chi2 = +115 degrees +/- 10 degrees . We have also validated the torsion angles of His37 in the inactivated state as chi1 = -175 degrees +/- 10 degrees and chi2 = -170 degrees +/- 10 degrees .     

Protein engineering of nitrile hydratase activity of papain: molecular dynamics study of a mutant and wild-type enzyme

The mechanism of hydrolysis of the nitrile (N-acetyl-phenylalanyl-2-amino-propionitrile, I) catalyzed by Gln19Glu mutant of papain has been studied by nanosecond molecular dynamics (MD) simulations. MD simulations of the complex of mutant enzyme with I and of mutant enzyme covalently attached to both neutral (II) and protonated (III) thioimidate intermediates were performed. An MD simulation with the wild-type enzyme.I complex was undertaken as a reference. The ion pair between protonated His159 and thiolate of Cys25 is coplanar, and the hydrogen bonding interaction S(-)(25).HD1-ND1(159) is observed throughout MD simulation of the mutant enzyme.I complex. Such a sustained hydrogen bond is absent in nitrile-bound wild-type papain due to the flexibility of the imidazole ring of His159. The nature of the residue at position 19 plays a critical role in the hydrolysis of the covalent thioimidate intermediate. When position 19 represents Glu, the imidazolium ion of His159-ND1(+).Cys25-S(-) ion pair is distant, on average, from the nitrile nitrogen of substrate I. Near attack conformers (NACs) have been identified in which His159-ImH(+) is positioned to initiate a general acid-catalyzed addition of Cys-S(-) to nitrile. Though Glu19-CO(2)H is distant from nitrile nitrogen in the mutant.I structure, MD simulations of the mutant.II covalent adduct finds Glu19-CO(2)H hydrogen bonded to the thioimide nitrogen of II. This hydrogen bonded species is much less stable than the hydrogen bonded Glu19-CO(2)(-) with mutant-bound protonated thioimidate (III). This observation supports Glu19-CO(2)H general acid catalysis of the formation of mutant.III. This is the commitment step in the Gln19Glu mutant catalysis of nitrile hydrolysis.     

Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics

pH is an important parameter in condensed-phase systems, because it determines the protonation state of titratable groups and thus influences the structure, dynamics, and function of molecules in solution. In most force field simulation protocols, however, the protonation state of a system (rather than its pH) is kept fixed and cannot adapt to changes of the local environment. Here, we present a method, implemented within the MD package GROMACS, for constant pH molecular dynamics simulations in explicit solvent that is based on the λ-dynamics approach. In the latter, the dynamics of the titration coordinate λ, which interpolates between the protonated and deprotonated states, is driven by generalized forces between the protonated and deprotonated states. The hydration free energy, as a function of pH, is included to facilitate constant pH simulations. The protonation states of titratable groups are allowed to change dynamically during a simulation, thus reproducing average protonation probabilities at a certain pH. The accuracy of the method is tested against titration curves of single amino acids and a dipeptide in explicit solvent.     

Molecular dynamics and brownian dynamics investigation of ion permeation and anesthetic halothane effects on a proton-gated ion channel

Bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) is activated to cation permeation upon lowering the solution pH. Its function can be modulated by anesthetic halothane. In the present work, we integrate molecular dynamics (MD) and Brownian dynamics (BD) simulations to elucidate the ion conduction, charge selectivity, and halothane modulation mechanisms in GLIC, based on recently resolved X-ray crystal structures of the open-channel GLIC. MD calculations of the potential of mean force (PMF) for a Na(+) revealed two energy barriers in the extracellular domain (R109 and K38) and at the hydrophobic gate of transmembrane domain (I233), respectively. An energy well for Na(+) was near the intracellular entrance: the depth of this energy well was modulated strongly by the protonation state of E222. The energy barrier for Cl(-) was found to be 3-4 times higher than that for Na(+). Ion permeation characteristics were determined through BD simulations using a hybrid MD/continuum electrostatics approach to evaluate the energy profiles governing the ion movement. The resultant channel conductance and a near-zero permeability ratio (P(Cl)/P(Na)) were comparable to experimental data. On the basis of these calculations, we suggest that a ring of five E222 residues may act as an electrostatic gate. In addition, the hydrophobic gate region may play a role in charge selectivity due to a higher dehydration energy barrier for Cl(-) ions. The effect of halothane on the Na(+) PMF was also evaluated. Halothane was found to perturb salt bridges in GLIC that may be crucial for channel gating and open-channel stability, but had no significant impact on the single ion PMF profiles.     

Multiple Proton Confinement in the M2 Channel from the Influenza A Virus

The tetrameric M2 protein bundle of the influenza A virus is the proton channel responsible for the acidification of the viral interior, a key step in the infection cycle. Selective proton transport is achieved by successive protonation of the conserved histidine amino acids at position 37. A recent X-ray structure of the tetrameric transmembrane (TM) domain of the protein (residues 22-46) resolved several water clusters in the channel lumen, which suggest possible proton pathways to the His37 residues. To explore this hypothesis, we have carried out molecular dynamics (MD) simulations of a proton traveling towards the His37 side chains using MD with classical and quantum force fields. Diffusion through the first half of the channel to the "entry" water cluster near His37 may be hampered by significant kinetic barriers due to electrostatic repulsion. However, once in the entry cluster, a proton can move to one of the acceptor His37 in a nearly barrierless fashion, as evidenced both by MD simulations and a scan of the potential energy surface (PES). Water molecules of the entry cluster, although confined in the M2 pore and restricted in their motions, can conduct protons with a rate very similar to that of bulk water.     

NMR detection of pH-dependent histidine-water proton exchange reveals the conduction mechanism of a transmembrane proton channel

The acid-activated proton channel formed by the influenza M2 protein is important for the life cycle of the virus. A single histidine, His37, in the M2 transmembrane domain (M2TM) is responsible for pH activation and proton selectivity of the channel. Recent studies suggested three models for how His37 mediates proton transport: a shuttle mechanism involving His37 protonation and deprotonation, a H-bonded imidazole-imidazolium dimer model, and a transporter model involving large protein conformational changes in synchrony with proton conduction. Using magic-angle-spinning (MAS) solid-state NMR spectroscopy, we examined the proton exchange and backbone conformational dynamics of M2TM in a virus-envelope-mimetic membrane. At physiological temperature and pH, (15)N NMR spectra show fast exchange of the imidazole (15)N between protonated and unprotonated states. To quantify the proton exchange rates, we measured the (15)N T(2) relaxation times and simulated them for chemical-shift exchange and fluctuating N-H dipolar fields under (1)H decoupling and MAS. The exchange rate is 4.5 × 10(5) s(-1) for Nδ1 and 1.0 × 10(5) s(-1) for Nε2, which are approximately synchronized with the recently reported imidazole reorientation. Binding of the antiviral drug amantadine suppressed both proton exchange and ring motion, thus interfering with the proton transfer mechanism. By measuring the relative concentrations of neutral and cationic His as a function of pH, we determined the four pK(a) values of the His37 tetrad in the viral membrane. Fitting the proton current curve using the charge-state populations from these pK(a)'s, we obtained the relative conductance of the five charge states, which showed that the +3 channel has the highest time-averaged unitary conductance. At physiologically relevant pH, 2D correlation spectra indicated that the neutral and cationic histidines do not have close contacts, ruling out the H-bonded dimer model. Moreover, a narrowly distributed nonideal helical structure coexists with a broadly distributed ideal helical conformation without interchange on the sub-10 ms time scale, thus excluding the transporter model in the viral membrane. These data support the shuttle mechanism of proton conduction, whose essential steps involve His-water proton exchange facilitated by imidazole ring reorientations.     

Theoretical studies on water–tetracaine interaction

The action of local anesthetics (LA) is controversial. There is experimental evidence that the unprotonated form of LA penetrates the axon, while the charged form acts in the intracellular phase. To obtain some insight on the structure of the local anesthetics tetracaine and its pharmacological action, we made calculations using the density functional theory (DFT) method. After those calculations, we performed molecular dynamics (MD) simulations in a p, N, T ensemble, in an aqueous environment, on both unprotonated and protonated forms of the molecule. The radial distribution function was used to study water solvent effects, through the characterization of the affinity of tetracaine to water. The results indicate that the molecule has regions with different degree of hydrophobicity, and the N‐terminal of the anesthetic was primarily affected by changes in the protonation state of the anesthetic. The pH‐dependent activity of TTC should then be analyzed in view of local changes in different regions of the molecule, rather than in terms of general effects on the hydrophobicity of the molecule as a whole. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

High-resolution mining of the SARS-CoV-2 main protease conformational space: supercomputer-driven unsupervised adaptive sampling

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs). The global exploration problem is decomposed into a set of separate MD trajectories that can be restarted within a selective process to achieve sufficient phase-space sampling. Accurate statistical properties can be obtained through reweighting. Within this highly parallel setup, the Tinker-HP package can be powered by an arbitrary large number of GPUs on supercomputers, reducing exploration time from years to days. This approach is used to tackle the urgent modeling problem of the SARS-CoV-2 Main Protease (M^pro) producing more than 38 μs of all-atom simulations of its apo (ligand-free) dimer using the high-resolution AMOEBA PFF. The first 15.14 μs simulation (physiological pH) is compared to available non-PFF long-timescale simulation data. A detailed clustering analysis exhibits striking differences between FFs, with AMOEBA showing a richer conformational space. Focusing on key structural markers related to the oxyanion hole stability, we observe an asymmetry between protomers. One of them appears less structured resembling the experimentally inactive monomer for which a 6 μs simulation was performed as a basis for comparison. Results highlight the plasticity of the M^pro active site. The C-terminal end of its less structured protomer is shown to oscillate between several states, being able to interact with the other protomer, potentially modulating its activity. Active and distal site volumes are found to be larger in the most active protomer within our AMOEBA simulations compared to non-PFFs as additional cryptic pockets are uncovered. A second 17 μs AMOEBA simulation is performed with protonated His172 residues mimicking lower pH. Data show the protonation impact on the destructuring of the oxyanion loop. We finally analyze the solvation patterns around key histidine residues. The confined AMOEBA polarizable water molecules are able to explore a wide range of dipole moments, going beyond bulk values, leading to a water molecule count consistent with experimental data. Results suggest that the use of PFFs could be critical in drug discovery to accurately model the complexity of the molecular interactions structuring M^pro.

We provide an unsupervised adaptive sampling strategy capable of producing μs-timescale molecular dynamics (MD) simulations of large biosystems using many-body polarizable force fields (PFFs).     

Nonequilibrium dynamics simulations of nitric oxide release: comparative study of nitrophorin and myoglobin

Nitrophorin 4 (NP4) is a heme protein that reversibly binds nitric oxide (NO), with release rates modulated by pH change. High-resolution structures of NP4 revealed that pH changes and NO binding induce a large conformational rearrangement in two loops that serve to protect the heme-bound NO molecule from solvent. We used extended (110 ns) molecular dynamics simulations of NP4 at pH 5 and pH 7, modeled by selective deprotonation of acidic groups. Conformational and dynamic changes were observed, consistent with those found in the crystal. Further, major solvent movement and NO escape were observed at pH 7, while the ligand remained in the heme binding pocket at pH 5. As a control, we also performed molecular dynamics (MD) simulations of sperm whale myoglobin, where NO migration into the interior cavities of the protein was observed, consistent with previous reports. We constructed a kinetic model of ligand escape to quantitatively relate the microscopic rate constants to the observed rates, and tested the predictions against the experimental data. The results suggest that release rates of diatomic molecules from heme proteins can be varied by several orders of magnitude through modest adjustments in geminate rebinding and gating behavior.     

A FLUORESCENCE STUDY OF TRYPTOPHAN-HISTIDINE INTERACTIONS IN THE PEPTIDE ANANTIN AND IN SOLUTION

Abstract—Anantin is a heptadecapeptide in which the C-terminal peptide chain pierces the covalently cyclized peptide ring formed by an amide link between the α-NH₂ end group and the β-carboxyl group of Asp(8). It contains a tryptophan and a histidine at positions 5 and 12 , respectively. Des-Phe(17)-anantin lacks the C-terminal phenylalanine. Fluorescence emission intensity as a function of pH follows the ionization of a single residue. The pK_a amounts to 7.23 ± 0.03 for anantin and is attributed to His(12). At pH 9 the quantum yield is 0.12 ± 0.01 for anantin, whereas at pH 4.5 the quantum yield decreases more than two-fold (0.05 2 0.01). Practically identical parameters are observed for des-Phe(17)-anantin. This pH dependency reveals intramolecular quenching of the excited indole ring of Trp(5) by the imidazole of His(12), which results in a marked decrease of the tryptophan fluorescence at low pH. In a multifrequency phase fluorometric study the fluorescence lifetimes for both peptides at pH 4.5 and pH 9 are determined. At both, pH fluorescence decay is well described by a sum of two exponentials. For anantin at pH 4.5 the lifetimes are 0.72 ± 0.07 ns and 1.67 ± 0.07 ns. At pH 9 the lifetimes are 1.11 ±0.12 ns and 2.55 ± 0.03 ns. In methanol we find two lifetimes for anantin: 0.68 ± 0.01 ns and 2.57 ± 0.01 ns. The lifetimes are found to be slightly dependent upon emission wavelength. For des-Phe(17)-anantin practically the same values are observed. The quenching of tryptophan fluorescence by histidine is further studied in solution using N-acetyl-tryptophanamide in the presence of increasing concentrations of imidazole in the protonated (pH 4.5) and unprotonated (pH 9) state and in methanol. At both pH values and in methanol, a linear increase in both the inverse of the steady-state fluorescence F_o/F and the inverse of the lifetime 1/τ with increasing imidazole concentration indicates that a collisional mechanism is at the root of the observed quenching. The quenching efficiency values, γ, are calculated and amount to about 0.32 at pH 4. 5 , 0.02 at pH 9 and 0.002 in methanol, showing that protonated imidazole is a better quencher than the unprotonated form, and that the nature of the solvent is involved even in the quenching by unprotonated imidazole. Tryptophan-histidine interactions in solution and in the peptide are compared.     

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.     

Alamethicin channels in a membrane: molecular dynamics simulations

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.     

Realistic simulations of proton transport along the gramicidin channel: demonstrating the importance of solvation effects

The nature of proton transduction (PTR) through a file of water molecules, along the gramicidin A (gA) channel, has long been considered as being highly relevant to PTR in biological systems. Previous attempts to model this process implied that the so-called Grotthuss mechanism and the corresponding orientation of the water file plays a major role. The present work reexamines the PTR in gA by combining a fully microscopic empirical valence bond (EVB) model and a recently developed simplified EVB-based model with Langevin dynamics (LD) simulations. The full model is used first to evaluate the free energy profile for a stepwise PTR process. The corresponding results are then used to construct the effective potential of the simplified EVB. This later model is then used in Langevin dynamics simulations, taking into account the correct physics of possible concerted motions and the effect of the solvent reorganization. The simulations reproduce the observed experimental trend and lead to a picture that is quite different from that assumed previously. It is found that the PTR in gA is controlled by the change in solvation energy of the transferred proton along the channel axis. Although the time dependent electrostatic fluctuations of the channel and water dipoles play their usual role in modulating the proton-transfer process (Proc. Natl. Acad. Sci. U.S.A. 1984, 81, 444), the PTR rate is mainly determined by the free energy profile. Furthermore, the energetics of the reorientation of the unprotonated water file do not appear to provide a consistent way of assessing the activation barrier for the PTR process. It seems to us that in the case of gA, and probably other systems with significant electrostatic barriers for the transfer of the proton charge, the PTR rate is controlled by the electrostatic barrier. This finding has clear consequences with regards to PTR processes in biological systems.     

Rectification of the current in alpha-hemolysin pore depends on the cation type: the alkali series probed by MD simulations and experiments

A striking feature of the alpha-hemolysin channel-a prime candidate for biotechnological applications-is the dependence of its ionic conductance on the magnitude and direction of the applied bias. Through a combination of lipid bilayer single-channel recording and molecular dynamics (MD) simulations, we characterized the current-voltage relationship of alpha-hemolysin for all alkali chloride salts at neutral pH. The rectification of the ionic current was found to depend on the type of cations and increase from Li(+) to Cs(+). Analysis of the MD trajectories yielded a simple quantitative model that related the ionic current to the electrostatic potential, the concentration and effective mobility of ions in the channel. MD simulations reveal that the major contribution to the current asymmetry and rectification properties originates from the cationic contribution to the current that is significantly reduced in a cationic dependent way when the membrane polarity is reversed. The variation of chloride current was found to be less important. We report that the differential affinity of cations for the charged residues positioned at the channel's end modulates the number of ions inside the channel stem thus affecting the current properties. Through direct comparison of simulation and experiment, this study evaluates the accuracy of the MD method for prediction of the asymmetric, voltage dependent conductances of a membrane channel.