Gramicidin A channel as a test ground for molecular dynamics force fields

We use the well-known structural and functional properties of the gramicidin A channel to test the appropriateness of force fields commonly used in molecular dynamics (MD) simulations of ion channels. For this purpose, the high-resolution structure of the gramicidin A dimer is embedded in a dimyristoylphosphatidylcholine bilayer, and the potential of mean force of a K(+) ion is calculated along the channel axis using the umbrella sampling method. Calculations are performed using two of the most common force fields in MD simulations: CHARMM and GROMACS. Both force fields lead to large central barriers for K(+) ion permeation, that are substantially higher than those deduced from the physiological data by inverse methods. In long MD simulations lasting over 60 ns, several ions are observed to enter the binding site but none of them crossed the channel despite the presence of a large driving field. The present results, taken together with many earlier studies, highlights the shortcomings of the standard force fields used in MD simulations of ion channels and calls for construction of more appropriate force fields for this purpose.     

Potential responses of nutrient membrane of frog's stomach to step changes in external K+ and Cl- concentrations

With an in vitro chamber method experiments were performed to determine the relative ionic conductances of the nutrient membrane (membrane facing muscularis mucosa). The concentration of a given ion in the nutrient bathing solution was changed, and the ensuing time course of the change in transmucosal potential difference (PD) was recorded. Changing K(+) from 4 to 79 mM produced a response in PD which occurred markedly faster than the response for the reverse change, and similar results were obtained by changing the Cl(-) concentration. It was found that these differences were predicted by the analysis of an idealized model consisting of a membrane in series with a diffusion barrier. When both the K(+) and Cl(-) were changed, such that the product of their concentrations remained constant, the time courses of the responses were again similiar to those predicted on the basis of the model. From the magnitudes of the total PD responses it is shown that in the presence of a 4 mM K(+) nutrient solution, the conductivity of the nutrient membrane appears to be entirely due to the K(+) and Cl(-) conductances, the K(+) conductance being about twice that of the Cl(-). It is also shown that with a 79 mM K(+) nutrient solution the parameters of the membrane were changed such that the conductances of the two ions were approximately equal. The time constant for diffusion of KCl or NaCl across the barrier consisting of the submucosa, muscularis mucosa, and lamina propria is about 1 min.     

Molecular dynamics simulations of the gramicidin A-dimyristoylphosphatidylcholine system with an ion in the channel pore region

To investigate the process of ion permeation in an ion channel systematically, we performed molecular dynamics (MD) simulations on a gramicidin A (GA)-phospholipid model system with an ion in the channel pore region. Each of the three types of ions (Ca2+, Na+ Cl-) was placed at five different positions along the channel axis by replacing a water molecule. MD simulations were performed on each system at constant pressure and constant temperature. The MD trajectories showed that the Ca2+ and Na+ ions could stably fluctuate in the pore region, but the Cl- ion was pushed out because of the unfavorable interaction with the channel. This result is consistent with experimental data. It was also found that the conformation of the GA channel underwent a significant change due to the presence of the ion, and the two ends of the GA monomer were more flexible than its middle region. In particular, the dramatic change of local pore radius near the ion indicated this kind of deformation. The strong interaction between the ion and carbonyl oxygen atoms of GA was the major contributor to this change. Furthermore, it was found that the ethanolamine group of the GA molecule was the most flexible group in the GA channel and often observed to block the entrance of GA. These results imply that the deformation of channel structure plays a very important factor in ion permeation, and the ethanolamine group may play a key role in regulating ion entry into the pore. In conclusion, our results indicate that the ion has a dominant influence on the structure of the GA channel and that the flexibility of the ion channel is a crucial factor in the ion permeation process.     

Brownian dynamics simulations of ion transport through the VDAC

It is important to gain a physical understanding of ion transport through the voltage-dependent anion channel (VDAC) because this channel provides primary permeation pathways for metabolites and electrolytes between the cytosol and mitochondria. We performed grand canonical Monte Carlo/Brownian dynamics (GCMC/BD) simulations to explore the ion transport properties of human VDAC isoform 1 (hVDAC1; PDB:2K4T) embedded in an implicit membrane. When the MD-derived, space-dependent diffusion constant was used in the GCMC/BD simulations, the current-voltage characteristics and ion number profiles inside the pore showed excellent agreement with those calculated from all-atom molecular-dynamics (MD) simulations, thereby validating the GCMC/BD approach. Of the 20 NMR models of hVDAC1 currently available, the third one (NMR03) best reproduces both experimental single-channel conductance and ion selectivity (i.e., the reversal potential). In addition, detailed analyses of the ion trajectories, one-dimensional multi-ion potential of mean force, and protein charge distribution reveal that electrostatic interactions play an important role in the channel structure and ion transport relationship. Finally, the GCMC/BD simulations of various mutants based on NMR03 show good agreement with experimental ion selectivity. The difference in ion selectivity between the wild-type and the mutants is the result of altered potential of mean force profiles that are dominated by the electrostatic interactions.     

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.     

Ion permeation through a narrow channel: using gramicidin to ascertain all-atom molecular dynamics potential of mean force methodology and biomolecular force fields

We investigate methods for extracting the potential of mean force (PMF) governing ion permeation from molecular dynamics simulations (MD) using gramicidin A as a prototypical narrow ion channel. It is possible to obtain well-converged meaningful PMFs using all-atom MD, which predict experimental observables within order-of-magnitude agreement with experimental results. This was possible by careful attention to issues of statistical convergence of the PMF, finite size effects, and lipid hydrocarbon chain polarizability. When comparing the modern all-atom force fields of CHARMM27 and AMBER94, we found that a fairly consistent picture emerges, and that both AMBER94 and CHARMM27 predict observables that are in semiquantitative agreement with both the experimental conductance and dissociation coefficient. Even small changes in the force field, however, result in significant changes in permeation energetics. Furthermore, the full two-dimensional free-energy surface describing permeation reveals the location and magnitude of the central barrier and the location of two binding sites for K(+) ion permeation near the channel entrance--i.e., an inner site on-axis and an outer site off-axis. We conclude that the MD-PMF approach is a powerful tool for understanding and predicting the function of narrow ion channels in a manner that is consistent with the atomic and thermally fluctuating nature of proteins.     

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.     

Gating kinetics of potassium channel in rat dorsal root ganglion neurons analyzed with fractal model

The kinetics of ion channels have been widely modeled as a Markov process. In these models it is assumed that the channel protein has a small number of discrete conformational states and kinetic rate constants connecting these states are constant. To study the gating kinetics of voltage-dependent K(+) channel in rat dorsal root ganglion neurons, K(+) channel current were recorded using cell-attached patch-clamp technique. The K(+) channel characteristic of kinetics were found to be statistically self-similar at different time scales as predicted by the fractal model. The fractal dimension D for the closed times and for the open times depend on the pipette potential. For the open and closed times of kinetic setpoint, it was found dependent on the applied pipette potential, which indicated that the ion channel gating kinetics had nonlinear kinetic properties. Thus, the open and closed durations, which had the voltage dependence of the gating of this ion channel, were well described by the fractal model.     

A Stochastic Analysis of a Brownian Ratchet Model for Actin-Based Motility

In recent single-particle tracking (SPT) measurements on Listeria monocytogenes motility in cells [Kuo and McGrath (2000)], the actin-based stochastic dynamics of the bacterium movement has been analyzed statistically in terms of the mean-square displacement (MSD) of the trajectory. We present a stochastic analysis of a simplified polymerization Brownian ratchet (BR) model in which motions are limited by the bacterium movement. Analytical results are obtained and statistical data analyses are investigated. It is shown that the MSD of the stochastic bacterium movement is a monotonic quadratic function while the MSD for detrended trajectories is linear. Both the short-time relaxation and the long-time kinetics in terms the mean velocity and effective diffusion constant of the propelled bacterium are obtained from the MSD analysis. The MSD of the gap between actin tip and the bacterium exhibits an oscillatory behavior when there is a large resistant force from the bacterium. For comparison, a continuous diffusion formalism of the BR model with great analytical simplicity is also studied.     

Molecular Dynamics Studies of Diffusion in Model Cylindrical Pores at Very Low Densities

Abstract

Molecular dynamics simulation has been used to study diffusion of methane at ambient temperature in cylindrical pores at very low densities. The cylinders were modelled as a continuum solid which interacts with the methane in the radial direction only. At the lowest densities, the VACF method does not yield reliable values of the self diffusion coefficient, D_s , but a suitable choice of time step and run length enables values of D_s to be found from MSD plots that are below the classical Knudsen diffusion coefficients. When density is increased, D_s passes through a maximum although the adsorption isotherm remains inside the Henry law region. Maxima are found for two cylinder radii and for two adsorbent field strengths. The existence of a maximum is attributed to transient intermolecular interactions. Analysis of a molecular trajectory demonstrates that long diffusion paths can be triggered by the rare event of an intermolecular encounter which forces a molecule into the repulsive part of the wall potential. At sufficiently high density, subsequent collisions quench the tendency towards long paths, and D_s decreases again. The issue of simulation artefact as a source of these observations is discussed.     

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

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

Ion transport across a phospholipid membrane mediated by the peptide trichogin GA IV

Trichogin GA IV is a special member of a class of peptaibols that are linear peptide antibiotics of fungal origin, characterised by the presence of a variable number of alpha-aminoisobutyric acid residues, an acyl group at the N-terminus and a 1,2-amino alcohol at the C-terminus. Most of the peptaibols display ion-channel-forming or at least membrane-modifying properties. The 11-residue-long trichogin GA IV is not only one of shortest peptaibols, but it is also unique for its n-octanoyl group instead of the more common found acetyl group at the N-terminus. For the first time we have found that this lipopeptaibol is able to enhance conduction of monovalent cations through membranes of large unilamellar vesicles (LUVs). The influence of the [Leu-OMe]trichogin GA IV analogue (TRI) on ion permeation was studied under a variety of conditions (lipid composition, lipid-to-peptide ratio and a transmembrane potential). Parallel experiments were performed with the 16-residue long, channel-forming peptaibol, zervamicin (ZER). For the two peptides, the permeability between K(+) and Na(+) was found to be different. In addition, the ion diffusion rate dependencies on the peptide concentration are observed to be different. This might indicate that a different number of aggregated molecules are involved in the rate-limiting step, i.e. 3-4 (TRI) and 4-7 (ZER). In the presence of TRI, dissipation of the transmembrane potential, Delta psi, was observed with a rate to be dependent on the magnitude of both initial Delta psi and peptide concentration. Both peptides were activated by a cis-positive but not by cis-negative Delta psi. Under identical conditions the ion-conducting efficiency of zervamicin was 100-200 times higher than that of trichogin. Our results show that, unlike for zervamicin, the membrane-modifying activity of trichogin is not associated with a channel mechanism.     

Computational modelling of the open-state Kv 1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.     

Water and deuterium oxide permeability through aquaporin 1: MD predictions and experimental verification

Determining the mechanisms of flux through protein channels requires a combination of structural data, permeability measurement, and molecular dynamics (MD) simulations. To further clarify the mechanism of flux through aquaporin 1 (AQP1), osmotic p(f) (cm(3)/s/pore) and diffusion p(d) (cm(3)/s/pore) permeability coefficients per pore of H(2)O and D(2)O in AQP1 were calculated using MD simulations. We then compared the simulation results with experimental measurements of the osmotic AQP1 permeabilities of H(2)O and D(2)O. In this manner we evaluated the ability of MD simulations to predict actual flux results. For the MD simulations, the force field parameters of the D(2)O model were reparameterized from the TIP3P water model to reproduce the experimentally observed difference in the bulk self diffusion constants of H(2)O vs. D(2)O. Two MD systems (one for each solvent) were constructed, each containing explicit palmitoyl-oleoyl-phosphatidyl-ethanolamine (POPE) phospholipid molecules, solvent, and AQP1. It was found that the calculated value of p(f) for D(2)O is approximately 15% smaller than for H(2)O. Bovine AQP1 was reconstituted into palmitoyl-oleoyl-phosphatidylcholine (POPC) liposomes, and it was found that the measured macroscopic osmotic permeability coefficient P(f) (cm/s) of D(2)O is approximately 21% lower than for H(2)O. The combined computational and experimental results suggest that deuterium oxide permeability through AQP1 is similar to that of water. The slightly lower observed osmotic permeability of D(2)O compared to H(2)O in AQP1 is most likely due to the lower self diffusion constant of D(2)O.     

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.     

獐茅盐腺形态结构及其泌盐性

显微观察结果显示,獐茅盐腺为典型的双细胞结构,主要分布在叶脉附近,这样有利于它快速收集来自根部的盐离子。同时,X-ray微区分析结果表明,獐茅盐腺可以有效地将Na^＋从表皮细胞和叶肉细胞转运到基细胞,由帽细胞分泌到体外,从而降低植物体内盐分水平以适应盐渍环境。用不同盐类处理材料,獐茅对Na^＋、K^＋和Ca^2＋的吸收和分泌均表现为具有不同的选择性,分泌量的顺序为Na^＋〉K^＋〉Ca^2＋,并且仅仅24h内,盐分的分泌量就已经超过了叶片内的含量,结果植株体内的总离子含量水平几乎不变,说明獐茅具有较强的泌盐能力。     

On the prediction of transport properties of ionic liquid using 1-n-butylmethylpyridinium tetrafluoroborate as an example

The Green–Kubo and Einstein–Helfand approaches are examined for calculating diffusion coefficient, electronic conductivity and shear viscosity of ionic liquid using 1-n-butylmethylpyridinium tetrafluoroborate [C₄PY][BF₄] as an example. Both methods suffer numerical errors accumulated at long simulation time, resulting divergences in the integrated autocorrelation time (IAT) and nonlinearity in the mean square displacement (MSD). The numerical errors can be reduced using smaller time step in the simulation. By identifying a converged plateau in IAT and a linear segment in MSD both approaches yield consistent predictions. Using a validated force field, the predicted diffusion coefficient and electrical conductivity agree reasonably well with the experimental data. However, the shear viscosity is significantly underestimated. Analysis of the simulation data indicates that a much slow relaxation in the pressure tensor must be considered, which is unfortunately infeasible due to the accumulated numerical error. Alternatively, the non-equilibrium periodic perturbation method shows promising improvement in the prediction.     

Effects of the gap junction blocker glycyrrhetinic acid on gastrointestinal smooth muscle cells

In the tunica muscularis of the gastrointestinal (GI) tract, gap junctions form low-resistance pathways between pacemaker cells known as interstitial cells of Cajal (ICCs) and between ICC and smooth muscle cells. Coupling via these junctions facilitates electrical slow-wave propagation and responses of smooth muscle to enteric motor nerves. Glycyrrhetinic acid (GA) has been shown to uncouple gap junctions, but previous studies have shown apparent nonspecific effects of GA in a variety of tissues. We tested the effects of GA using isometric force measurements, intracellular microelectrode recordings, the patch-clamp technique, and the spread of Lucifer yellow within cultured ICC networks. In murine small intestinal muscles, beta-GA (10 muM) decreased phasic contractions and depolarized resting membrane potential. Preincubation of GA inhibited the spread of Lucifer yellow, increased input resistance, and decreased cell capacitance in ICC networks, suggesting that GA uncoupled ICCs. In patch-clamp experiments of isolated jejunal myocytes, GA significantly decreased L-type Ca(2+) current in a dose-dependent manner without affecting the voltage dependence of this current. The IC(50) for Ca(2+) currents was 1.9 muM, which is lower than the concentrations used to block gap junctions. GA also significantly increased large-conductance Ca(2+)-activated K(+) currents but decreased net delayed rectifier K(+) currents, including 4-aminopyridine and tetraethylammonium-resistant currents. In conclusion, the reduction of phasic contractile activity of GI muscles by GA is likely a consequence of its inhibitory effects on gap junctions and voltage-dependent Ca(2+) currents. Membrane depolarization may be a consequence of uncoupling effects of GA on gap junctions between ICCs and smooth muscles and inhibition of K(+) conductances in smooth muscle cells.     

Roles of Extensibility and Turgor in Gibberellin- and Dark-stimulated Growth

The elongation response elicited by incubating excised hypocotyl sections of lettuce (Lactuca sativa L.) in light in gibberellin (GA) can be enhanced by the addition of Cl(-), Br(-), and NO(3) (-) salts of K(+) and Na(+). Sections incubated in light in the absence of GA do not elongate in response to the addition of salts. In contrast, excised hypocotyls incubated in darkness elongate equally in both GA and water, and their elongation can also be enhanced by KCl treatment. Growth stimulation by the salts of K(+) and Na(+) occurs optimally at 10 mm and the magnitude of the response is proportional to the duration of salt treatment. Although the growth of sections incubated in light in the absence of GA is not enhanced by various salts of K(+) and Na(+), the concentration of these cations exceeds that in GA-treated sections. In dark-grown tissue, uptake of K(+) also occurs in both GA- and H(2)O-treated sections incubated in 10 mm KCl. Since increased osmotic potential resulting from cation uptake does not correlate with growth stimulation resulting from salt treatments, we conclude that increased cell turgor is not the principal driving force for growth in hypocotyl sections. Changes in the extensibility of GA-treated, light-grown tissue and dark-grown tissue incubated with and without GA correlate with the increased growth rate of these sections. Incubation of sections in KCl results only in changes in water potential of sections without having a significant effect on extensibility. When changes in water potential are accompanied by increased extensibility, however, a marked increase in growth rate is observed.     

短杆菌肽A-DMPC通道内离子输运的分子动力学模拟

用最近提出的构建膜体系初始构象的有效方法 ,构建了在DMPC脂膜环境下短杆菌肽A通道模型 (GA -DMPC)。通过对Na 、Ca2 、Cl-三种不同离子在GA -DMPC通道内不同位置的分子动力学模拟 ,研究离子在通道内输运过程中与通道及通道内水分子的相互作用 ,从分子动力学的角度阐明离子在通道内的输运机制。主要计算结果表明 :(1)离子在通道内的输运使GA的构象发生变化 ,GA的柔性是离子在通道内通透的重要因素 ;(2)Cl- 离子可扩大通道半径 ,Na 离子和Ca2 离子则减小通道半径。Cl-离子不能在GA通道内通透 ;(3)离子的出现使通道内水分子的偶极方向发生变化。上述结果均与实验相符。