Biomedical subjects

Shin-Ho Chung

Publications and source records attributed to Shin-Ho Chung.

18 recordsLinked to original sources

PubMed2006

Brownian dynamics investigation into the conductance state of the MscS channel crystal structure.

We suggest that the crystal structure of the mechanosensitive channel of small conductance is in a minimally conductive state rather than being fully activated. Performing Brownian dynamics simulations on the crystal structure show that no ions pass through it. When simulations are conducted on just the transmembrane domain (excluding the cytoplasmic residues 128 to 280) ions are seen to pass through the channel, but the conductance of approximately 30 pS is well below experimentally measured values. The mutation L109S that replaces a pore lining hydrophobic residue with a polar one is found to have little effect on the conductance of the channel. Widening the hydrophobic region of the pore by 2.5 Angstrom however, increases the channel conductance to over 200 pS suggesting that only a minimal conformational change is required to gate the pore.

Crystallization↗

PubMed2006

Adaptive Brownian dynamics simulation for estimating potential mean force in ion channel permeation.

Ion channels are biological nanotubes formed by large protein molecules in the cell membrane. This paper presents a novel multiparticle simulation methodology, which we call adaptive controlled Brownian dynamics, for estimating the force experienced by a permeating ion at each discrete position along the ion-conducting pathway. The profile of this force, commonly known as the potential of mean force, results from the electrostatic interactions between the ions in the conduit and all the charges carried by atoms forming the channel the protein, as well as the induced charges on the protein wall. The current across the channel is solely determined by the potential of mean force encountered by the permeant ions. The simulation algorithm yields consistent estimates of this profile. The algorithm operates on an angstrom unit spatial scale and femtosecond time scale. Numerical simulations on the gramicidin ion channel show that the algorithm yields the potential of mean force profile that accurately reproduces experimental observations.

Adaptation, Physiological↗

PubMed2005

The fast gating mechanism in ClC-0 channels.

We investigate and then modify the hypothesis that a glutamate side chain acts as the fast gate in ClC-0 channels. We first create a putative open-state configuration of the prokaryotic ClC Cl- channel using its crystallographic structure as a basis. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to ClC-0 by replacing all the nonconserved polar and charged residues. Using this open-state channel model, we carry out molecular dynamics simulations to study how the glutamate side chain can move between open and closed configurations. When the side chain extends toward the extracellular end of the channel, it presents an electrostatic barrier to Cl- conduction. However, external Cl- ions can push the side chain into a more central position where, pressed against the channel wall, it does not impede the motion of Cl- ions. Additionally, a proton from a low-pH external solution can neutralize the extended glutamate side chain, which also removes the barrier to conduction. Finally, we use Brownian dynamics simulations to demonstrate the influence of membrane potential and external Cl- concentration on channel open probability.

Biophysics↗

PubMed2005

Electrostatic basis of valence selectivity in cationic channels.

We examine how a variety of cationic channels discriminate between ions of differing charge. We construct models of the KcsA potassium channel, voltage gated sodium channel and L-type calcium channel, and show that they all conduct monovalent cations, but that only the calcium channel conducts divalent cations. In the KcsA and sodium channels divalent ions block the channel and prevent any further conduction. We demonstrate that in each case, this discrimination and some of the more complex conductance properties of the channels is a consequence of the electrostatic interaction of the ions with the charges in the channel protein. The KcsA and sodium channels bind divalent ions strongly enough that they cannot be displaced by other ions and thereby block the channel. On the other hand, the calcium channel binds them less strongly such that they can be destabilized by the repulsion of another incoming divalent ion, but not by the lesser repulsion from monovalent ions.

Animals↗

PubMed2005

Homology model of the GABAA receptor examined using Brownian dynamics.

We have developed a homology model of the GABA(A) receptor, using the subunit combination of alpha1beta2gamma2, the most prevalent type in the mammalian brain. The model is produced in two parts: the membrane-embedded channel domain and the extracellular N-terminal domain. The pentameric transmembrane domain model is built by modeling each subunit by homology with the equivalent subunit of the heteropentameric acetylcholine receptor transmembrane domain. This segment is then joined with the extracellular domain built by homology with the acetylcholine binding protein. The all-atom model forms a wide extracellular vestibule that is connected to an oval chamber near the external surface of the membrane. A narrow, cylindrical transmembrane channel links the outer segment of the pore to a shallow intracellular vestibule. The physiological properties of the model so constructed are examined using electrostatic calculations and Brownian dynamics simulations. A deep energy well of approximately 80 kT accommodates three Cl(-) ions in the narrow transmembrane channel and seven Cl(-) ions in the external vestibule. Inward permeation takes place when one of the ions queued in the external vestibule enters the narrow segment and ejects the innermost ion. The model, when incorporated into Brownian dynamics, reproduces key experimental features, such as the single-channel current-voltage-concentration profiles. Finally, we simulate the gamma2 K289M epilepsy inducing mutation and examine Cl(-) ion permeation through the mutant receptor.

Acetylcholine↗

PubMed2005

A model of sodium channels.

We have explored the permeation and blockage of ions in sodium channels, relating the channel structure to function using electrostatic profiles and Brownian dynamics simulations. The model used resembles the KcsA potassium channel with an added external vestibule and a shorter selectivity filter. The electrostatic energy landscape seen by permeating ions is determined by solving Poisson's equation. The two charged amino acid rings of Glu-Glu-Asp-Asp (EEDD) and Asp-Glu-Lys-Ala (DEKA) around the selectivity filter region are seen to play a crucial role in making the channel sodium selective, and strongly binding calcium ions such that they block the channel. Our model closely reproduces a range of experimental data including the current-voltage curves, current-concentration curves and blockage of monovalent ions by divalent ions.

Calcium↗

PubMed2005

Brownian dynamics simulation for modeling ion permeation across bionanotubes.

The principles underlying Brownian dynamics (BD), its statistical consistency, and algorithms for practical implementation are outlined here. The ability to compute current flow across ion channels confers a distinct advantage to BD simulations compared to other simulation techniques. Thus, two obvious applications of BD ion channels are in calculation of the current-voltage and current-concentration curves, which can be directly compared to the physiological measurements to assess the reliability of the model and predictive power of the method. We illustrate how BD simulations are used to unravel the permeation dynamics in two biological ion channels-the KcsA K+ channel and CIC Cl- channel.

Cell Membrane↗

PubMed2004

Influence of protein flexibility on the electrostatic energy landscape in gramicidin A.

We describe an electrostatic model of the gramicidin A channel that allows protein atoms to move in response to the presence of a permeating ion. To do this, molecular dynamics simulations are carried out with a permeating ion at various positions within the channel. Then an ensemble of atomic coordinates taken from the simulations are used to construct energy profiles using macroscopic electrostatic calculations. The energy profiles constructed are compared to experimentally-determined conductance data by inserting them into Brownian dynamics simulations. We find that the energy landscape seen by a permeating ion changes significantly when we allow the protein atoms to move rather than using a rigid protein structure. However, the model developed cannot satisfactorily reproduce all of the experimental data. Thus, even when protein atoms are allowed to move, the dielectric model used in our electrostatic calculations breaks down when modeling the gramicidin channel.

Computer Simulation↗

PubMed2004

Conduction mechanisms of chloride ions in ClC-type channels.

The conduction properties of ClC-0 and ClC-1 chloride channels are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We create an open-state configuration of the prokaryotic ClC Cl(-) channel using its known crystallographic structure as a basis. Two residues that are occluding the channel are slowly pushed outward with molecular dynamics to create a continuous ion-conducting path with the minimum radius of 2.5 A. Then, retaining the same pore shape, the prokaryotic ClC channel is converted to either ClC-0 or ClC-1 by replacing all the nonconserved dipole-containing and charged amino acid residues. Employing open-state ClC-0 and ClC-1 channel models, current-voltage curves consistent with experimental measurements are obtained. We find that conduction in these pores involves three ions. We locate the binding sites, as well as pinpointing the rate-limiting steps in conduction, and make testable predictions about how the single channel current across ClC-0 and ClC-1 will vary as the ionic concentrations are increased. Finally, we demonstrate that a ClC-0 homology model created from an alternative sequence alignment fails to replicate any of the experimental observations.

Amino Acids↗

PubMed2003

A model of the glycine receptor deduced from Brownian dynamics studies.

We have developed a three-dimensional model of the alpha1 homomeric glycine receptor by using Brownian dynamics simulations to account for its observed physiological properties. The model channel contains a large external vestibule and a shallow internal vestibule, connected by a narrow, cylindrical selectivity filter. Three rings of charged residues from the pore-lining M2 domain are modeled as point charges in the protein. Our simulations reproduce many of the key features of the channel, such as the current-voltage profiles, permeability ratios, and ion selectivity. When we replace the ring of alanine residues lining the selectivity filter with glutamates, the mutant model channel becomes permeable to cations, as observed experimentally. In this mutation, anions act as chaperones for sodium ions in the extracellular vestibule, and together they penetrate deep inside the channel against a steep energy barrier encountered by unaccompanied ions. Two subsequent amino acid mutations increase the cation permeability, enabling monovalent cations to permeate through the channel unaided and divalent cations to permeate when chaperoned by anions. These results illustrate the key structural features and underlying mechanism for charge selectivity in the glycine receptor.

Alanine↗

PubMed2003

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.

Binding Sites↗

PubMed2003

Dielectric self-energy in Poisson-Boltzmann and Poisson-Nernst-Planck models of ion channels.

We demonstrated previously that the two continuum theories widely used in modeling biological ion channels give unreliable results when the radius of the conduit is less than two Debye lengths. The reason for this failure is the neglect of surface charges on the protein wall induced by permeating ions. Here we attempt to improve the accuracy of the Poisson-Boltzmann and Poisson-Nernst-Planck theories, when applied to channel-like environments, by including a specific dielectric self-energy term to overcome spurious shielding effects inherent in these theories. By comparing results with Brownian dynamics simulations, we show that the inclusion of an additional term in the equations yields significant qualitative improvements. The modified theories perform well in very wide and very narrow channels, but are less successful at intermediate sizes. The situation is worse in multi-ion channels because of the inability of the continuum theories to handle the ion-to-ion interactions correctly. Thus, further work is required if these continuum theories are to be reliably salvaged for quantitative studies of biological ion channels in all situations.

Calcium↗

PubMed2003

Adaptive learning algorithms for nernst potential and I-V curves in nerve cell membrane ion channels modeled as hidden Markov models.

We present discrete stochastic optimization algorithms that adaptively learn the Nernst potential in membrane ion channels. The proposed algorithms dynamically control both the ion channel experiment and the resulting hidden Markov model signal processor and can adapt to time-varying behavior of ion channels. One of the most important properties of the proposed algorithms is their self-learning capability--they spend most of the computational effort at the global optimizer (Nernst potential). Numerical examples illustrate the performance of the algorithms on computer-generated synthetic data.

Algorithms↗

PubMed2002

Ion channels: recent progress and prospects.

Determination of the crystal structure of the KcsA potassium channel and its subsequent refinement at 2 A resolution have stimulated much interest in modelling of ion channels. Here we review the recent developments in ion channels research, focusing especially on the question of structure-function relationships, and discuss how permeation models based on Brownian and molecular dynamics simulations can be used fruitfully in this endeavour.

Calcium Channels↗

PubMed2002

Continuum electrostatics fails to describe ion permeation in the gramicidin channel.

We investigate the validity of continuum electrostatics in the gramicidin A channel using a recently determined high-resolution structure. The potential and electric field acting on ions in and around the channel are computed by solving Poisson's equation. These are then used in Brownian dynamics simulations to obtain concentration profiles and the current passing through the channel. We show that regardless of the effective dielectric constant used for water in the channel or the channel protein, it is not possible to reproduce all the experimental data on gramicidin A; thus, continuum electrostatics cannot provide a valid framework for the description of ion dynamics in gramicidin channels. Using experimental data and molecular dynamics simulations as guides, we have constructed potential energy profiles that can satisfactorily describe the available physiological data. These profiles provide useful benchmarks for future potential of mean force calculations of permeating ions from molecular dynamics simulations of gramicidin A. They also offer a convenient starting point for studying structure-function relationships in modified gramicidin channels.

Biophysical Phenomena↗

PubMed2002

Modeling diverse range of potassium channels with Brownian dynamics.

Using the experimentally determined KcsA structure as a template, we propose a plausible explanation for the diversity of potassium channels seen in nature. A simplified model of KcsA is constructed from its atomic resolution structure by smoothing out the protein-water boundary and representing the atoms forming the channel protein as a homogeneous, low dielectric medium. The properties of the simplified and atomic-detail models, deduced from electrostatic calculations and Brownian dynamics simulations, are shown to be qualitatively similar. We then study how the current flowing across the simplified model channel changes as the shape of the intrapore region is modified. This is achieved by increasing the radius of the intracellular pore systematically from 1.5 to 5 A while leaving the dimensions of the selectivity filter and inner chamber unaltered. The strengths of the dipoles located near the entrances of the channel, the carbonyl groups lining the selectivity filter, and the helix macrodipoles are kept constant. The channel conductance increases steadily as the radius of the intracellular pore is increased. The rate-limiting step for both the outward and inward current is the time it takes for an ion to cross the residual energy barrier located in the intrapore region. The current-voltage relationship obtained with symmetrical solutions is linear when the applied potential is less than approximately 100 mV but deviates slightly from Ohm's law at higher applied potentials. The nonlinearity in the current-voltage curve becomes less pronounced as the radius of the intracellular pore is increased. When the strengths of the dipoles near the intracellular entrance are reduced, the channel shows a pronounced inward rectification. Finally, the conductance exhibits the saturation property observed experimentally. We discuss the implications of these findings on the transport of ions across the potassium channels and membrane channels in general.

Biophysical Phenomena↗

PubMed2002

Conducting-state properties of the KcsA potassium channel from molecular and Brownian dynamics simulations.

The mechanisms underlying transport of ions across the potassium channel are examined using electrostatic calculations and three-dimensional Brownian dynamics simulations. We first build open-state configurations of the channel with molecular dynamics simulations, by pulling the transmembrane helices outward until the channel attains the desired interior radius. To gain insights into ion permeation, we construct potential energy profiles experienced by an ion traversing the channel in the presence of other resident ions. These profiles reveal that in the absence of an applied field the channel accommodates three potassium ions in a stable equilibrium, two in the selectivity filter and one in the central cavity. In the presence of a driving potential, this three-ion state becomes unstable, and ion permeation across the channel is observed. These qualitative explanations are confirmed by the results of three-dimensional Brownian dynamics simulations. We find that the channel conducts when the ionizable residues near the extracellular entrance are fully charged and those near the intracellular side are partially charged. The conductance increases steeply as the radius of the intracellular mouth of the channel is increased from 2 A to 5 A. Our simulation results reproduce several experimental observations, including the current-voltage curves, conductance-concentration relationships, and outward rectification of currents.

Bacterial Proteins↗

PubMed2002

Reservoir boundaries in Brownian dynamics simulations of ion channels.

Brownian dynamics (BD) simulations provide a practical method for the calculation of ion channel conductance from a given structure. There has been much debate about the implementation of reservoir boundaries in BD simulations in recent years, with claims that the use of improper boundaries could have large effects on the calculated conductance values. Here we compare the simple stochastic boundary that we have been using in our BD simulations with the recently proposed grand canonical Monte Carlo method. We also compare different methods of creating transmembrane potentials. Our results confirm that the treatment of the reservoir boundaries is mostly irrelevant to the conductance properties of an ion channel as long as the reservoirs are large enough.

Biophysical Phenomena↗