Statistical Modeling and Application of Gramicidin A Ion Channel by Kai Y i u Luk B. A. Sc., The University of British Columbia, 2005 A THESIS SUBMITTED IN PARTIAL F U L F I L M E N T OF T H E REQUIREMENTS F O R T H E D E G R E E OF Master of Applied Science in The Faculty of Graduate Studies (Electrical and Computer Engineering) The University Of British Columbia July, 2007 © Kai Y i u Luk 2007 Abstract Ion channels are aqueous pores in the cell membrane for selected ions to flow down their electrochemical gradient. These channels play a prominent role i n a variety of biological processes i n the human body. Determining the structure and function of ion channels is of fundamental importance in biology. Also, the selective conductivity and specific gating mechanism of ion channels have attracted much interest in the area of artificial molecular detectors. Ion channel based biosensors are developed to detect molecu- lar species of interest i n medical diagnostics, environmental monitoring and general bio-hazard detection. This thesis is concerned with statistical techniques used to describe ion channel permeation and to develop ion channel based biosensors. Brownian dynamics is a popular technique to simulate ion channel permeation but is too computationally expensive to run when ionic concentration is high. B y fitting binding site statistics of B D simulation to a semi-Markov chain, we obtain a simpler model with conduction properties that are statistically the same as the simulations. This approach enables the use of extrapolation techniques to predict channel conduction when performing the actual simulation is computationally infeasible. Numerical studies on the simulation of gramicidin A channels are presented. In a separate study, we show the use of statistical modeling and detection techniques as ii Abstract part of a sensitive biosensing platform. A nano-scale biosensor is built by incorporating dimeric gramicidin A channels into bilayer membranes of giant unilamellar liposomes. The presence of specific target molecules changes the statistics of the biosensor's conduction. B y capturing the change i n real time, we devise a maximum likelihood detector to detect the presence of target molecules. T h e performance of the biosensor is tested with the addition of various target molecules known to inhibit conduction of gramicidin A channels. Experimental results show that the detection performed well even when the conductance change was difficult to visualize. The detection algorithm provides a sensitive detection system for ongoing development of membrane-based biosensors. iii Table of Contents Abstract ii Table of Contents iv List of Tables vii List of Figures viii Acknowledgements x 1 Introduction 1 1.1 Permeation Model of Ion Channel 3 1.2 Gated Ion Channel Based Biosensor 4 1.3 Organization 5 2 Background 6 2.1 Gramicidin A Channel 7 2.2 Brownian Dynamics 9 2.2.1 Simulation Setup 2.2.2 Solving Langevin Equations for Ion Dynamics 2.2.3 Solving Poisson's Equation for Systematic Forces 10 . . . . 10 . . 11 iv Table of Contents 2.3 3 4 2.2.4 Estimating Brownian Dynamics Current 12 2.2.5 Motivation for Statistical Modeling 13 Previous Work on Ion Channel Biosensor 15 2.3.1 Composition of I C S ™ B i o s e n s o r 15 2.3.3 Motivation for Statistical Detection 17 A Simple Statistical M o d e l to Characterize B r o w n i a n D y namics Permeation 19 3.1 Statistical Modeling of Ion Channel Permeation 20 3.1.1 Finite State Representation 21 3.1.2 Semi Markov Model of Ion Channel Permeation . . . 24 3.2 Extrapolation Technique for Model Prediction 30 3.3 Numerical Results 33 3.3.1 Parameter Estimation for Binding Site Kinetics 3.3.2 Model Prediction at High Concentration . . . 33 34 B i s - g A Channel Biosensor 40 4.1 Experimental Construction of bis-gA Channel Biosensor . . . 42 4.2 Statistical Modeling and Validation of Biosensor Current 44 . . 4.2.1 Parameter Estimation for Biosensor Current 47 4.2.2 Statistical Validation of Biosensor Model 49 4.3 Target Molecule Detection Algorithm 50 4.4 Experimental Results of Biosensor 52 4.4.1 Model Estimation and Validation for Biosensor . . . . 53 4.4.2 Real-Time Detection of Target Molecules 56 Table of Contents 5 Conclusion and Extensions Bibliography 61 63 Appendices A List of A c r o n y m s 69 List of Tables 3.1 Anderson Darling statistics of holding time distributions . . . 29 3.2 Estimated transition matrices at various ionic concentrations. 34 3.3 Extrapolated transition matrices at high ionic concentrations. 35 4.1 M a x i m u m likelihood estimation of H M M parameters 54 vii List of Figures 2.1 Chemical structure of gramicidin A dimers 2.2 Computing time of 1 jus of Brownian dynamics simulation of 8 gramicidin A channel 14 2.3 Components of the I C S ™ b i o s e n s o r 16 2.4 Detection mechanism of I C S ™ b i o s e n s o r 17 3.1 Shape of gramicidin A channel and location of binding sites. . 21 3.2 Sample path of BD-simulated ion movement inside gramicidin A channel 3.3 22 Effect of binding K i o n on the potential energy of gramicidin + A channel 23 3.4 Finite state mapping of ion sample paths 24 3.5 Comparision between E C D F and estimated gamma C D F for Q31 (£) at ionic concentration 400 m M 30 3.6 Example of logistic function 32 3.7 Estimated state holding time distributions 35 3.8 Comparison between estimated and extrapolated state holding time distributions 3.9 36 Comparison between SMC-simulated and BD-estimated current 37 viii List of Figures 3.10 Comparison between SMC-extrapolated results and best-fit Michaelis-Menten curve 39 4.1 Fluorescence and phase-contrast image of biosensor 42 4.2 Block diagram of experimental setup 43 4.3 Power spectral density of biosensor response 45 4.4 Chemical structure of Methylbenzthonium Chloride 51 4.5 Chemical structure of 2-Methyl-4-tert.-octylph.enol 51 4.6 Experimental measurement of biosensor response 53 4.7 M a x i m u m likelihood estimate of biosensor conductance . . . . 55 4.8 Q-statistics of residuals 55 4.9 Biosensor's response to addition of M B C 56 4.10 Filtered likelihoods for estimated models of M B C 57 4.11 M a x i m u m likelihood detection of M B C 58 4.12 Biosensor's response to various concentrations of M T O P . . . 58 4.13 Filtered likelihoods for estimated models of M T O P 59 4.14 M a x i m u m likelihood detection of M T O P 60 Acknowledgements I would like to thank my supervisor, D r . V i k r a m Krishnamurthy, for his guidence. Thanks to my collaborator, Bruce Cornell, for providing figures for this thesis. Thanks to my sister Candy for helping me with the biology and Phoebe for her support and encouragement. Thanks to my parents for their unconditional support. This thesis is dedicated to them. Chapter 1 Introduction Ion channels are aqueous pores in the cell membrane for selected ions to now down their electrochemical gradient. The class of ion channels is very diverse. In many channels, ion conduction is controlled by a gate, which opens and closes i n response to electrical and chemical signals. In addition, ion channels are highly selective, allowing only desired ions to pass through. Thus ion channels are used to regulate membrane potential as well as intracellular and extracellular concentrations of ions. Ion channels play a prominent role in a variety of biological processes in the human body. In the nervous system, voltage-gated sodium and potassium channels open and close co-operatively to propagate nerve impulses, also known as action potentials, along the axon. In cardiac muscles, voltagegated calcium channels release calcium ions to stimulate contraction of the muscles. Not surprisingly, the malfunctioning of ion channels is also the cause of many diseases, collectively known as channelopathies. These diseases include epilepsy, muscular disorders, cystic fibrosis and diabetes. [3] Knowledge of the structure and function of ion channels can help biologists understand the cause of and ultimately treat these diseases. Therefore, the study of ion channels is of fundamental importance in biology. In 2003, R o d - 1 Chapter 1. Introduction erick MacKinnon's work i n the determination of structure and operation of ion channels with crystallographic analyses won him a share of the Nobel Prize i n Chemistry. [17, 18] Ion channels have also attracted much interest i n the area of artificial molecular detectors. Their selective conductivity and specific gating mechanisms make ion channels an ideal biological recognition component of a biosensor. In addition, ion channels are also highly sensitive; the binding of a single target molecule to an ion channel can lead to macroscopic level of conduction. A s a result, ion channel based biosensors are designed to detect molecular species of interest across a wide range of applications that include medical diagnostics, environmental monitoring and general biohazard detection. This thesis is concerned with statistical techniques used to describe ion channel permeation and to develop ion channel based biosensors. First we study the statistics of computer-simulated ion channel permeation and predict the channel behaviour in conditions where the simulation approach is computationally infeasible. In a separate study, we apply a statistical model to describe macroscopic output current recorded from an ion channel membrane-based biosensor. We show that the model, together with our detection algorithm, provides a sensitive and accurate biosensing platform. The reminder of this chapter is organized as follows. Sec. 1.1 introduces the computer simulation approach used to study ion channel permeation. Sec. 1.2 discusses some considerations for using gated ion channels to develop a biosensor. Finally, Sec. 1.3 discusses the organization for the rest of this thesis. 2 Chapter 1. Introduction 1.1 Permeation Model of Ion Channel A permeation model describes the propagation of individual ions through an ion channel. With recent advances in computational power, one popular approach is to study the permeation of ion channels with computer simulations of channel interactions at an atomic level. These techniques include molecular dynamics (MD) and Brownian dynamics (BD). The general philosophy is to observe permeation properties of the channel in the simulations and to match with electrophysiological measurements. In molecular dynamics, each atom in the channel protein, channel water and the environment near the channel mouth are modeled explicitly. At each discrete time step of the simulation, the forces on each atom are calculated and the location of each atom is adjusted accordingly. This approach computes the exact model of the channel, but the computational cost is very high, making it infeasible to use molecular dynamics to simulate ion conductions. However, molecular dynamics remains a useful tool to study different properties of ion channels. [4, 6] Brownian dynamics is similar to molecular dynamics in that the movement of individual ions are calculated; however, it reduces computational complexity by treating water molecules in the channel as a continuum. The interaction between ions and the water molecules can be modeled as a random-force and a frictional force. The charges of channel protein is assumed to be rigid. Also, B D simulations can operate at a much slower timescale than M D simulations. [39] Typical B D simulations are performed in the order of 100 fs per time step while M D simulations operate at 1 fs per time step. As 3 Chapter 1. Introduction a result, B D has computational complexity several magnitudes lower than that of MD and is the preferred technique for simulating conduction events to estimate channel current. Other theoretical approaches are used to study ion channel permeation. Reaction rate theory is based on partitioning of the channel into specific energy wells and estimating conductance from the transition probabilities between the wells. Poisson-Nernst-Planck theory model both ions and water molecules as continua and compute the dynamics using electrodiffusion theory and macroscopic electrostatics. These two techniques are not discussed in this thesis. For a detailed review of these techniques to model ion channel permeation, see [24, 31]. 1.2 Gated Ion Channel Based Biosensor Triggered by recent threats of biological warfare, there is a high demand for fast, economical and portable methods to detect the presence of harzardous biochemicals. Traditional laboratory techniques such as mass spectrometry and polymerase chain reaction (PCR) are not suited for this type of usage. One particular type of detector that is gaining popularity is a receptor-based sensor that performs biological recognition at the molecular level. Receptorbased sensor uses complementary molecules that bind to a specific type of target molecules and then indicate the result in the form of fluorescence [35], surface plasmon resonance [38] or other observable signals. In this part of the thesis, we focus on the ion channel based biosensor. An ion channel based biosensor is typically built by inserting ion chan- 4 Chapter 1. nels into artificial lipid membranes. Introduction The biosensor consists of a receptor element engineered to bind to a specific type of target molecules. The successful binding of the target molecule to the receptor element distorts the channel pore. This results in an alteration in the channel's conductance or gating mechanism. Detection of this change can be acheived by monitoring membrane conductivity, p H levels or intracellular concentration of ions. In addition, ion channels are highly sensitive devices and are suitable tools for signal amplification. The binding of a single target molecule to the receptor system can lead to an open channel current in the order of pico-Amperes. 1.3 Organization The rest of this thesis is organized as follows. Chapter 2 studies some background information i n ion channel permeation model and ion channel based biosensor. In particular, the technique of simulating ion channel conduction with Brownian dynamics is discussed in detail to identify the difficulty in running these simulations at high ionic concentration. Also, we summarize the development efforts in ion channel based biosensing and use a specific example of ion channel based biosensor to illustrate its functional properties. In Chapter 3, we introduce a statistical model to capture the binding site kinetics of an ion channel i n Brownian dynamics simulation. The modeling, combined with extrapolation techniques, are used to predict the permeation behaviour of Brownian dynamics at high ionic concentration. Finally, in Chapter 4, we introduce the use of stochastic modeling and maximum likelihood detection results as part of a sensitive biosensing platform. 5 Chapter 2 Background In this chapter, we study Brownian dynamics, a practical technique used to simulate ion channel conduction events, and investigate previous development efforts in the area of ion channel based biosensor. Through the discussion of theories and examples, we hope to motivate the use of statistical techniques to improve our understanding in ion channel permeation and to assist the development of sensitive and reliable biosensors. Brownian dynamics have been widely used to determine the structurefunction relationship of various ion channels, such as sodium channels [28, 40], potassium channels [26], calcium channels [14] and gramicidin channels [27]. These papers utlize different parameterization and optimization techniques to determine various structural properties of ion channels. In [27], the authors computed the potential of mean force experienced by a cation inside gramicidin channel. In [28], the authors parameterized and adaptively determined the molecular structure of sodium channels. The general approach in these papers is similar. The structural properties are determined by optimizing the fit between experimental current measured patch clamp techniques and simulated current estimated from Brownian dynamics. Ion channels can perform chemical recognition at a molecular level and 6 Chapter 2. Background the opening and closing of channels results in a macroscopic change in the conduction level. Thus the use of artificial gated ion channels for biosensing is a promising strategy. The majority of the development focuses on improving the sensitivity and stability of the sensor. T h a t has prompted the use of a variety of ion channels and membranes to construct artificial biosensors. Multiple studies have shown the use of engineered pores to detect a variety of ions [9], proteins [32] and D N A [21]. The rest of this chapter is organized as follows. First, in Sec. 2.1 we discuss the permeation and conduction properties of gramicidin A (gA) channels. The algorithms presented in C h . 3 and C h . 4 are based on experimental and simulation results of these channels. In Sec. 2.2, we provide an overview of Brownian dynamics simulation of an ion channel. In particular, we describe the simulation setup, dynamic equations to determine ion movement, and estimation of current using Brownian dynamics. In Sec. 2.3, we provide a case study of a biosensor developed using gramicidin A channels. The specific example is used to identify the major components of an ion channel based biosensor and to highlight the development efforts in this area. 2.1 Gramicidin A Channel Gramicidin A channels are used in antibotics to increase cation permeability across the bacterial cell wall, thus destroying the cell's ion gradient. In addition to its pharmaceutical use, gramicidin A is important as a model ion channel and is often used to test molecular theories of ion channel per- 7 Chapter 2. Background meation. [1] For a recent review of g A channels, see [41]. Gramicidin A is a linear polypeptide consisting of fifteen amino acids. In its conducting form, two gramicidin A peptides are linked transiently by six hydrogen bonds to form a channel, which is 4 A i n diameter and 25 A in length. The chemical structure of gramicidin A dimer is shown in Figure 2.1. The gating mechanism of g A channels results from the formation and breakdown of peptide bonds. Inside the g A channel exists six water molecules. During the permeation of an ion through the g A channel, the ion drags the water molecules i n single file. T h e structure of dimerized g A channels in a lipid bilayer has been resolved by N M R studies. [25] Figure 2.1: Chemical structure of gramicidin A dimers 8 Chapter 2. Background Gramicidin A channels are selective for monovalent cations and impermeable to divalent cations and anions. N M R experiments [33] and M D simulations [5] show that there exists two monovalent binding sites at approximately ± 9A from the center of the gA channel. The cation conductivity of gA channels varies linearly as a function of the bias voltage. The single channel current i depends on the cation concentration \S] in the form s of the Michaelis-Menten curve: W[S] ( 2 1 ) [S}+K S where i max and Ks are parameter constants that depend on the cation type [20]. 2.2 Brownian Dynamics In this section, we summarize the use of Brownian dynamics to estimate channel current. Poisson's equations and Langevin's equations are solved iteratively to compute the movement of ions. Channel current can be estimated by counting ion crossing events. At the end of the chapter, we discuss some computational difficulties faced in B D simulations and show a simpler statistical approach to study Brownian dynamics. The actual Brownian dynamics computer program belongs to the Computational Biophysics Group of Research School of Biological Sciences at the Australian National University (http://langevin.anu.edu.au/). For more detailed analysis on Brownian dynamics simulations, see [22, 23]. 9 Chapter 2. 2.2.1 Background Simulation Setup Our simulation setup consists of two reservoirs TZi and TZ connected by 2 an ion channel C. A n external potential is applied across the channel. A t the start of B D simulation, each reservoir contains N K + ions and N C l ~ ions indexed by i = 1,2,...,N and i = N + 1 , N + 2,...,2N respectively. Let A denote the set of different experimental conditions, which represents varying external potentials across the channel and ion concentrations in the € Ti? and reservoirs. Let of ion i at time t. £ Ti? denote the position and velocity = (x^',y^', z^'), where x[ \ and zf l are the 1 (i) Cartesian components of the position vector. The velocity vector ' follows a similar construct. 2.2.2 S o l v i n g L a n g e v i n E q u a t i o n s for I o n D y n a m i c s The position and velocity of each ion evolve according to the Langevin equations. Let X = (xf ^, x p ' ' , . . . , X j 1 t 2 A ^')' denote the positions and V t — (v^', v p ' ' , . . . , v p ^ ' ) ' denote the velocities of all 2N ions in our setup. 4 =4 + t^ds, Jo ] (2.2) ] For i = 1,2, ...,N, m+v (i) t = m+v^ - + [ FJ> (X )ds ] Jo s / Jo + 2m m -y {x®)v®ds + + + 7 + (xf )wf ) ) (2-3) 10 Chapter 2. Background For i — N + 1,N + 2,..., 2N m 7-(xW)vWd. _ !s + 2 m - - ( x f ) w f ) (2.4) : 7 . F ^ ( X t ) is the systematic forces acting on ion i for experimental con- dition A € A. Details about the systematic forces are given in Sec. 2.2.3. The frictional coefficient m 7 ( x i ) = m 7 ± ± t ) ± the diffusion coefficient of K+ is D = ± if x£ l) eKiUR 2 and = 1.96 x 1CT m /s in a bulk solution, + 9 2 while the diffusion coefficient for Cl~ions is D~ = 2.03 x 10~ m /s. The 9 2 (i) process is a 3-dimensional Brownian motion with component-wise in- dependence and can be written as a zero-mean Gaussian random variable with a 3 x 3 diagonal covariance matrix. 2.2.3 S o l v i n g P o i s s o n ' s E q u a t i o n for S y s t e m a t i c Forces The systematic forces in (2.3) and (2.4) can be rewritten as: Fi i ) (X ) t where the scalar valued process = - . « V <_^(X ) t rienced by ion i given the position X t x W $ i ° ( X t (2.5) ) is the total electric potential expe- of the 27V ions. The potential <_^(X ) t experienced by each ion i consists of five components: * « ( X t ) = < ^ ( X t ) + $ e x t ( + c _ f > ( X ) l t x W + } + sjw C/(x«) ( x M) (2.6) Chapter 2. «£ext^ w^ d x eno Background t e s the external potential applied along the z axis of the ion channel. $ ^ ( x ^ ) denotes the ion-wall interaction potential, also known as the c r / r potential. ^' {X.t) 9 tential and ^f ' (X.t) R l l denotes the inter-ion Coulomb po- denotes the short range ion-ion potential. C/(x^) denotes the potential of mean force ( P M F ) , which is a smooth function used to model the Coulomb interaction between the ion and channel protein as well as induction of surface charges at the water-protein interface when the ion is near the protein wall. For the simulations discussed in this thesis, we have used the P M F that has the best fit to experimental results of gramicidin A current measurements. Details of this result can be found in [29]. The systematic forces F ^ ( X j ) are calculated as a function of the ion positions. This is done by solving Poisson's equation. By assuming that the simulation space can be partitioned into several regions with different dielectric constants. Let €j be the dielectric constant of region j and we can write the modified Poisson's equation as: V $,- = 2 2.2.4 (2.7) Estimating Brownian Dynamics Current The computer simulation of ion channel implements a discretized version of the dynamics described in Sec. 2.2.2 and Sec. 2.2.3. Suppose the channel model for experimental condition A is simulated for L seconds. At each discrete time step, the Langevin's equations in (2.2), (2.3) and (2.4) and modified Poisson's equation in (2.7) are solved to compute the movement of ions. Two crossing events are of importance to the estimation of channel 12 Chapter 2. Background conductance: 1) forward crossing event, during which an ion crosses from 7-i to 7_2, and 2) backward crossing event, during which an ion crosses from 7.2 to TZ\. Let and N-ji -ji lt 2 be the number of forward and back- ward crossing events respectively. We can compute the Brownian dynamics estimate of the single channel current: x = l where q + ^ ( is the electric charge of a proton, 1.6 x 10~ 19 ' C. The current estimate in (2.8) is based on the random samples observed in the simulation and therefore cannot be represented by a closed form expression. 2.2.5 M o t i v a t i o n for S t a t i s t i c a l Modeling Brownian dynamics is a very useful tool to study the structure-function relationship of ion channels. The technique simulates ion channel interaction at an atomic level and is computationally tractable to use BD to simulate channel crossing events because of a few simplifying assumptions discussed in Sec. 1.1. As a result, biophysicists often use BD to estimate single channel current at various applied voltages and concentrations and compare with experimental patch-clamp measurements. Although B D is generally computationally tractable, obtaining an accurate current estimates at multiple voltages and concentrations can take a very long time. Increasing ionic concentration involves placing additional ions in both reservoirs. That has a quadratic effect on the computational 13 Chapter 2. Background complexity of the B D simulations. Fig. 2.2 shows the time required to finish 1 [is of BD simulation of gramicidin A channel. The simulation is performed on the supercomputer cluster on Westgrid (glacier.westgrid.ca). 25 0 i 100 i i 200 300 1 1 . 1 400 500 600 700 Reservoir Concentration (mM) . 1 800 900 Figure 2.2: Computing time of 1 [is of Brownian dynamics simulation of gramicidin A channel. It can be seen that the computing time increases dramatically as more ions are included in the simulation. To obtain an accurate estimate of the channel current, typical BD simulations of an ion channel last 10 [is to 100 [is [11]. Generally it is rare to see BD simulations with reservoir concentrations above IM. In Chapter 3, we construct a statistical fully-parameterized model to characterize ion channel permeation observed in Brownian dynamics. Instead of modeling movements of all ions in the system, we focus on the kinetics of ion activities inside the channel. As a result, the model is a much simpler representation of Brownian dynamics and we show that conduction 14 Chapter 2. Background statistics are completely recoverable. One useful application of this model is model prediction. We can apply extrapolation techniques to predict channel current at conditions where BD simulation is computationally intractable. 2.3 Previous Work on Ion Channel Biosensor To better illustrate the biosensing problem, we investigate previous work in the area of ion channel based biosensor. In particular, we describe the composition and detection mechanism of the Ion Channel Switch ( I C S ™ ) biosensor. [2, 13] The gating mechanism of the novel biosensor exploits the association and disassociation probabilities of the gramicidin A dimers. This flexible and adaptive sensor can detect a variety of molecules such as growth factors, glucose, and D N A strands. [11] 2.3.1 Composition of I C S ™ B i o s e n s o r Figure 2.3 shows a diagram of the major components in the ICS™sensor. The ICS™biosensor is built by incorporating gramicidin A monomers into a tethered lipid bilayer. Because of their chemical and structural stability, gA channels are often used in the development of biosensors. The gA monomers in the inner leaflet of the bilayer are immobile while the monomers in the outer leaflet are free to diffuse. As discussed earlier in the chapter, gA molecules can only conduct when two gA monomers dimerize to form a channel. Therefore, the conductance of the lipid bilayer arises from random diffusion of free-moving monomers that dimerize to tethered monomers in the inner leaflet. 15 Chapter 2. Background JKmt * km Figure 2.3: Components of the I C S ™ b i o s e n s o r consists of (a) gold electrode and (b) tethered lipid bilayer. (c) Immobile g A monomer is embedded in the inner leaflet of the bilayer and (d) free-moving g A monomer embedded in the outer leaflet. In this setup, the receptor system, made of (e) streptavidin and (f) biotinylated receptor, is used to detect a specific (g) target molecule. [11] 1G Chapter 2. Background molecules with matching antibodies. (a) (b) Figure 2.4: Detection mechanism of I C S ™ b i o s e n s o r . (a) W h e n no target molecules are present, the g A monomer i n the outer leaflet is free to diffuse and couple with immobile monomers in the inner leaflet, (b) When target molecule is present, the receptor system binds to the target and the attaching g A monomer is no longer mobile. [11] 2.3.3 M o t i v a t i o n for S t a t i s t i c a l D e t e c t i o n The above example of a biosensor illustrates how an artificial ion channel based system can be designed to detect the presence of target molecules. Other ion channel biosensors have employed different gating mechanisms, such as channel blockage by molecules or channel distortion using antichannel antibodies. In general, the majority of the research focuses on achieving large differences in the biosensor conductance with and without the presence of target molecules. Little research has been done on extracting information from the biosensor's response using signal processing techniques. In C h . 4, we propose the use of statistical modeling and maximum likelihood detection algorithms as part of a biosensing platfom. The algorithm 17 Chapter 2. Background offers an automatic approach to analyze the results from the biosensor in real-time. This approach eliminates the majority of electrical and mechanical interfering effects and offers the ability to enhance the specificity and affinity of the biosensor. 18 Chapter 3 A Simple Statistical M o d e l to Characterize Brownian Dynamics Permeation of Gramicidin A Channels Brownian dynamics have been widely used to determine structure and function of various ion channels. The general approach is to optimize the fit between the channel current measured with patch clamp techniques and the one estimated from B D simulations at various external voltages and ionic concentrations. However, as discussed in Sec. 2.2.5, the computational complexity of Brownian dynamics makes it computationally intractable to perform the simulation at high concentrations. We propose a simpler statistical model to characterize ion permeation simulated in Brownian dynamics. In particular, we use the simple gramicidin A channel to demonstrate the feasibility of this approach. We model only the binding site statistics of gA channel with a finite-state semi-Markov process. To evaluate the performance of this model, we compare conduction statistics 19 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation of the proposed model with that obtained from Brownian dynamics simulation. The model captures dwell time and transition statistics with a small number of parameters and enables us to apply extrapolation techniques to predict the behaviour of Brownian Dynamics at high concentrations. The rest of this chapter is organized as follows. In Sec. 3.1, we map the continuous simulation space into finite states of binding site occupation and fit it to a three-state semi-Markov chain (SMC). In Sec. 3.2, we apply extrapolation techniques to predict channel current at high concentrations where performing the actual simulation is computationally intractable. Finally, numerical results of the modeling and prediction of gA channels are presented in Sec. 3.3. 3.1 Statistical Modeling of Ion Channel Permeation In this section, we use statistical modeling techniques to capture ion dynamics inside the channel. The aim is to develop a simple model that describes only the binding site statistics. By mapping the continuous space into discrete states of the channel, we construct and estimate a finite-state semiMarkov chain that is statistically indistinguishable from sample paths of actual BD simulations. Statistical validation methods are carefully applied at various stages of the modeling to verify goodness-of-fit. 20 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation 3.1.1 Finite State Representation A Brownian dynamics simulation is extremely complex, consisting of 2N ions each with its own position and velocity vectors. To reduce the dimensionality of the space, we choose to only model the occupancy of the binding sites in the gramicidin A channel. Figure 3.1 shows the shape of the channel. The shaded regions — z 2 < z < —z\ and z\ < z < z are the left and 2 right binding sites respectively, z = —Z3 and z — z% denote the reservoirchannel boundaries of the simulation space. "Transition regions" are defined in the regions —Z3 < z < —z , —z\ < z < z\ and z < z < Z3 to smooth 2 2 out unimportant events caused by ions bouncing around near the region boundaries. Only the z-coordinates are considered because K ions can only + pass through the narrow gA channel in single-file. -20 -16 -10 -5 0 5 10 Channel Axial Position (Angstroms) 1E 20 Figure 3.1: Shape of gramicidin A channel and location of binding sites. Figure 3.2 shows a sample path of BD simulated ion movements. Ions 21 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation inside the channel spend the majority of the time in the monovalent binding sites centered at around ± 9 A . The existence of the binding sites are discussed in Sec. 2.1. We argue that since the ions spend negligible time elsewhere in the channel, it is sufficient to model the dwell time and transition statistics in the binding site regions. 0 0.1 0.2 0.3 0.4 0.5 Time (^s) Figure 3.2: Sample path of BD-simulated ion movement inside gramicidin A channel. The sequence clearly shows the ions spending most of the time in the binding sites. The color of the trace is used to differentiate between ions currently in the channel and does not belong to the same ion for the entire sequence. Let C 6 {0,1} and 1Z € {0,1} denote the occupancy of the left and right binding site respectively, where 0 indicates an empty binding site and 1 an occupied binding site. When an ion is within a transition region, C and TZ simply assume their previous values. Due to the close proximity of the two binding sites, it is extremely rare to have both sites occupied for a long period of time. For example, in Figure 22 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation 3.2, an ion (in red) enters the left binding site at 0.21 us and both binding sites are occupied concurrently. However, the presence of the left ion forces the right ion (in blue) out of the channel soon after. To further illustrate this scenario, we observe the change in potential energy when a K i o n is at various locations near the left binding site. In + Figure 3.3 we plot only the energy due to the P M F and inter-ion Coulomb potential. The other three components in (2.6) are ignored for simplicity. The presence of an ion at the left binding site reduces the energy well from 6 kT to 2 kT. 12 -6 -8 0 2 4 6 8 10 12 Channel Axial Coordinates (Angstroms) 14 16 Figure 3.3: Effect of binding K i o n on the potential energy of gramicidin A channel. + The energy E is related to the occupation probability p as follows: p oc exp(—E) (3.1) 23 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation Based on the above relationship, the probability of an ion binding to the right site is 55 times smaller when the left site is occupied than when it is empty. We conclude that the state C = TZ = 1 is very rare and treat it a transition state. Let the state space be <S e {00,10,01}, where the first and second digit represent £ and TZ respectively. 3.1.2 Semi M a r k o v M o d e l of Ion C h a n n e l Permeation In Sec. 3.1.1, we showed the mapping from all the ion positions X to three t finite states, representing the occupancy of the binding sites in the channel. In Figure 3.4, we plot the finite state mapping of the ion sample paths observed in Figure 3.2. 01 0.2 0.3 Time (^ts) Figure 3.4: Finite state mapping of ion sample paths observed in Fig. 3.2. The smoothing technique results in crisp transitions that adequately describes the state of the binding sites. 24 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation M o d e l i n g and E s t i m a t i o n of S M C In this section, we introduce the semi-Markov chain, a probabilistic model used to capture the transition probabilities and rates of a finite state model. SMC, a popular technique used in reliability engineering and lifetime analysis, is parameterized by state space <S, transition matrix P and state holding time matrix Q. Here is a summary of the binding site kinetics that result in a state transition: • 00 -> 10 An ion from the left reservoir binds to the left site. • 10 ^ 00 An ion from the left site enters the left reservoir. • 10 -> 01 (a) A n ion from the left site jumps to the right site. (b) An ion from the right reservoir binds to the right site and pushes the ion at the left site out of the channel. • 00 -> 01 An ion from the right reservoir binds to the right site. • 01 -> 00 An ion from the right site enters the right reservoir. • 01 -> 10 (a) An ion from the right site jumps to the left site. 25 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation (b) An ion from the left reservoir binds to the left site and pushes the ion at the right site out of the channel. Note that both 10 —> 01 and 01 —> 10 can result from one of two possible processes. For the rest of this chapter, we refer to process (a) as the jumping process and process (b) as the pushing process. Let /3 23 denote the proba- bility that 10 —> 01 event is due to the jumping process and (1 — /? ) the 23 probability of the pushing process. Similarly, let /3 32 denote the probability that the 01 —> 10 event is due to the jumping process and (1 — /3 ) the 32 probability of the pushing process. To simplify the notations, we rewrite state 00, 10 and 01 to be state 1, 2 and 3 respectively. Let pij denote the transition probability from state i to state j, i, j = 1, 2 and 3. The transition matrix P can be written as: P = 0 Pl2 P21 0 P31 P32 P13 p (3.2) 2 3 0 Similarly we can write Q as a matrix of state holding time distributions. — Q = 921 (t) 931(*) 9i2 (*) — 9 2 (*) 3 qu (t) q£\t) (3.3) — 26 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation where qij (t) denote the probability distribution function (PDF) of the state holding times from state i to state j, i,j = 1, 2 and 3. The distribu- tions qn(t), 9 2 2 ( 0 and 533(t) are undefined because state transitions in a semi-Markov model must have different start and end states. The model is complicated by the fact that the transitions 10 —> 01 and 01 —> 10 can arise from two different processes and therefore possess different PDFs. For the transition 10 —> 01, let 923,jump(*) denote the P D F for the jumping process and <723,push(£) the P D F for the pushing process. We can write q2l {t) as: 3 23, 923 (*) = \ fl N with probability/3 923,jum (i) P 23 ( - ) 3 I 923, ushM with probability (1 - / 3 P 932,jump (t), 932, ushW and q^{t) P 2 3 4 ) can be defined similarly. As a result, Q is a 3 x 3 matrix of probability distributions, of which two distributions are Bernoulli-modulated. Experimentally we found that the gamma distribution provided the best fit to the state holding time distributions. Other PDFs considered include exponential, normal, log-normal and Weibull distribution. The gamma distribution is defined as: * * ' ( | a 6 ) = &M^* ~ a 1 ( r *' ( 3 ' 5 ) where a and b are known as the shape and scale parameter and T(-) is the gamma function. The estimation of the holding times as a gamma distribution allows us to model each state holding time distribution with only two parameters. Suppose n different transition times of a particular transition 27 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation are observed in a B D simulation. Let ti be the i th largest state holding time. T h e ordering of the holding times is irrelevant for the estimation but is important for the statistical validation method to be introduced later. The maximum-likelihood estimates for the gamma distribution can be computed by solving the following set of equations: l n o ^ _ t = iog </>(a)- log (a) (3.6) < - > 3 7 where t is the mean of the random samples and d and b are the maximum likelihood estimated shape and scale parameter. Statistical Validation of S M C The Anderson-Darling statistics is the standard technique used to evaluate the goodness of fit in distribution fitting. T h e Anderson-Darling statistics is part of a large class known as the empirical distribution function ( E D F ) statistics, which measures the difference between the estimated distribution and the empirical distribution of the random samples: A — -n 2 where Zj = 1 Q(ti\a,b) " i J2( ~ !) § Zi + (2n + l- 2i) log (1 - z ) (3.8) 2i lo { is the cumulative distribution function ( C D F ) of the estimated gamma distribution. Table 3.1 shows the Anderson-Darling statis- 28 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation tics computed at concentration 100 mM, 400 mM and 800 mM. (a) - 0.994 0.779 0.154 - 0.669 0.608, 0.916 0.256, 0.343 - (b) - 1.237 2.054 0.250 - 0.506 0.898, 0.513 0.483, 0.219 - (c) - 1.537 0.940 1.277 - 0.381 0.636, 1.717 0.727, 0.456 - Table 3.1: Anderson Darling statistics of holding time distributions at ionic concentration (a) 100 mM, (b) 400 mM and (c) 800 mM. The statistics for qij(t) are shown in cell (i,j) of each table. Multiple entries exist in cell (2,3) and (3,2) because separate statistics are computed for the jumping (left) and pushing (right) process. It can be seen that the Anderson Darling statistics for most state holding time distributions fall below the critical value of 1.092 found in [15]. In fact, 81% of the transition times satisfied the Anderson-Darling critical value for the gamma distribution at 0.01 significance level. Out of all transitions in Table 3.1, the holding time distribution <?3i(£) at concentration of 400 mM showed the biggest deviation from the gamma distribution. The empirical C D F of (j3i (t) and the estimated gamma C D F are compared in Figure. 3.5. It can be seen that the estimated gamma C D F followed the empirical C D F of the holding times closely even though it failed the Anderson Darling test. Thus it is a reasonable approximation to estimate all the holding times as 29 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation gamma distributions. 0.9 • 0.8 0.7 • 0.6 • 0.5 • 0.4 • 0.3 0.2 • o" 0 Figure 3.5: Comparision between E C D F and estimated gamma C D F for Qsi{t) at ionic concentration 400 mM. Once the transition matrix P and the parameters in the state holding time matrix Q are estimated, a Monte-Carlo simulation of the SMC produces an estimate of the channel current. 3.2 Extrapolation Technique for Model Prediction In this section, we devise an extrapolation procedure that can predict the behaviour of Brownian dynamics. In particular, we are interested in predicting the channel current at high ionic concentration because BD is intractable at that region. It is important to choose a reasonable functional form for the extrapolation to reduce the effect of outliers and estimation errors. General 30 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation extrapolation such as cubic spline leads to overfitting. Here we make two reasonable assumptions about the effect of increasing ionic concentrations: 1. Increasing ionic concentration has a monotonic effect on some transition probabilities and dwell times. For example, an empty binding site is likely to be filled quickly when more ions are in the reservoirs. Other transitions may be unaffected by concentration. 2. The effect of increasing ionic concentration saturates at very high concentrations. Linear extrapolation is monotonic but its simple functional form does not model saturation very well. Polynomial extrapolation of second order or higher is much more flexible but does not guarantee monotonicity. To satisfy the above assumptions, we choose the four-parameter logistic function: f{x) = A 0 + ^~ l (3.9) A 1+ ke d where AQ, A^, k and d are the function parameters. Fig. 3.6 plots an example of the logistic function, commonly used in biology and economics. In addition to its monotonicity, the logistic function is bounded from above or below as x —* oo. The function is flexible enough to model increasing, decreasing and constant functions. The effect of increasing ionic concentration on the transition probabilities and state holding times is unknown and it seems reasonable to use the logistic function as the basis to extrapolate these parameters. 31 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation X Figure 3.6: Example of logistic function, with AQ = 0, A^ = 1, k = 1 and d= 1. 32 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation To perform the model prediction, permeation properties of the g A channel at multiple reservoir concentrations are simulated using Brownian dynamics. Each set of simulation result is mapped to the three-state semi- Markov process described in Sec. 3.1.2. T h e estimated parameters in transition probabilities P and holding time matrix Q are then fitted to the logistic function using nonlinear least square methods such as the Gauss-Newton method. Finally, the parameters of the predicted model can be extrapolated from the fitted logistic function. 3.3 Numerical Results In this section, we present some numerical results based on B D simulations of the g A channel and evaluate the performance of our modeling and prediction techniques. We simulated the ion channel at an external voltage of 100 m V . A total of nine independent simuations were performed with the reservoir concentration varying from 100 m M to 900 m M , in increments of 100 m M . The duration of each simulation is 8 ^s. 3.3.1 P a r a m e t e r E s t i m a t i o n for B i n d i n g Site K i n e t i c s Table 3.2 shows the estimated transition probabilities, which nicely captures statistical changes of the the binding site kinetics at various ionic concentrations. For example, P(l,2) and P(l,3) does not vary with increasing concentration. W i t h concentration in both reservoirs being equal, it is reasonable to see that the ratio of ion entry from the left and right reservoir is constant. 33 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation (a) (b) 0 0.811 0.879 0.639 0 0.121 0 0.622 0.616 0.638 0 0.384 0.361 0.189 0 0 0.715 0.723 0.361 0.378 0 0 0.607 0.581 0.639 0 0.277 (c) 0.361 0.285 0 (d) 0.633 0 0.419 0.367 0.393 0 (e) 0 0.582 0.527 0.647 0 0.473 0.353 0.418 0 Table 3.2: Estimated transition matrices at ionic concentration (a) 100 mM, (b) 400 mM, (c) 700 mM, (d) 800 mM and (e) 900 mM. The probability for Pij are shown in cell of each table. Furthermore, as the concentration increases, it is less likely for an ion in the binding site to leave and go back into the reservoir. Consequently, both P(2,l) and P(3,l) are monotonically decreasing. The effect of increasing concentration on several holding time distribution can be seen in Figure 3.7. In particular, the biggest effect of increasing concentration is the decrease in state holding time 912 (0- An e m p t y left binding site is filled much more quickly as the left reservoir becomes more crowded with ions. 3.3.2 M o d e l P r e d i c t i o n at H i g h C o n c e n t r a t i o n To demonstrate the feasibility of our extrapolation technique, we predict the channel current at 800 m M and 900 mM, using simulation results from 34 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation (c) Figure 3.7: (d) Estimated state holding time distributions (a) qn{t), (b) _23,jump(0) ( ) _3i(*) and (d) <732, ush(£) at ionic concentrations of 100 m M , c P 400 m M and 700 m M . (a) (b) 0 0.636 0.364 0 0.636 0.364 0.591 0 0.409 0.568 0 0.432 0.601 0.399 0 0.576 0.424 0 Table 3.3: Extrapolated transition matrices at ionic concentration (a) 800 m M and (b) 900 m M . T h e probability for pij are shown in cell of each table. 35 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation 100 mM to 700 mM. This allows us to match BD-simulated and SMCextrapolated results at 800 mM and 900 mM and evaluate the performance of our model. The first seven BD simulations, at ionic concentration 100 mM to 700 mM, are estimated and the SMC parameters at each concentration are fitted to the logistic function as a function of concentration. The extrapolation is used to determine the SMC parameters at 800 mM and 900 mM. The extraporated transition matrix is shown in Table 3.3. (c) (d) Figure 3.8: Comparison between estimated and extrapolated state holding time distributions (a) qu{t), (b) q23,jum (i), (c) 931 (*) and (d) <732, ushW at ionic concentrations of 800 mM and 900 mM. P P In Figure 3.8, we compare the extrapolated distribution with results 36 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation estimated from BD simulations. With the exception of g23,jump(*)i the ex- trapolated distributions are very similar to the estimated distributions. One possible explanation of the poor performance in extrapolating q23 jump{t) t is the lack of observations of that transition. The 10 —> 01 jumping process is one of the rarest transitions in BD simulations of gA channels; therefore, the estimation of this process may have higher variance than other transitions. The extrapolated SMC model at 800 mM and 900 mM are simulated for 1000 independent runs. During each 8 ps run, the ion movements into and out of the binding sites are simulated not with Poisson's and Langevin's equations, but instead with the extrapolated SMC models. Figure 3.9 shows the result of the simulations. 100 200 300 400 500 600 Ionic Concentration (mM) 700 800 900 Figure 3.9: Comparison between SMC-simulated and B D estimated currents show similar conduction levels. At the concentrations lOOnM to 700nM, it is not surprising that the BDsimulated and SMC-simulated current have similar conduction level. The 37 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation SMC model was estimated directly from random paths observed in Brownian dynamics. The extrapolation at 800 mM and 900 mM also performed extremely well. Clearly the SMC-extrapolated exhibits very similar lev- els as the ones estimated with Brownian dynamics. Using the two-sample Kolmogorov-Smirnov test, we compare the similarities of the conduction distribution simulated in B D and SMC. The Kolmogorov-Smirnov test is another type of E D F testing techniques and can be used to test general distributions. The tradeoff for this generality is that the test is less effective than other E D F testing techniques designed for specific distributions. [10] Out of 1000 independent simuations at the concentration 800 nM, 96.5% of the time the Komogorov-Smirnov test requirement were satisfied, with an average P-value of 0.5295. At 900 nM, 97.8% of the test passed with an average P-value of 0.5069. Thus the hypothesis that the BD-simulated and SMC-simulated conduction events share the same distribution cannot be rejected. At this point, we do not have simulation results beyond ionic concentration of 900 mM and therefore cannot test our modeling and prediction technique at higher concentrations. We conclude this section with a prediction of channel current up to 5 M . These estimates are also obtained by simulating the extrapolated model for 1000 independent runs. The extrapolation results are plotted in Figure 3.10. Since there are no BD results to compare with our extrapolated current, we compare our results with the Michaelis-Menten kinetics described in 2.1. The plot shows that the SMC technique predicts that BD estimated current saturates at a lower conduction level than that of the Michaelis-Menten 38 Chapter 3. A Simple Statistical Model to Characterize Brownian Dynamics Permeation 3 2.5 + + + + + • © © / / .© / / O B D simulated •f S M C extrapolated - Michaelis-Menten extrapolated © Ionic Concentration (M) Figure 3.10: Comparison between SMC-extrapolated results and best-fit Michaelis-Menten curve. The comparison shows that the SMC extrapolation expects the current to saturate at a lower conduction level than that of the Michaelis-Menten curve. curve. 39 Chapter 4 Bis-gA Channel Biosensor In this chapter, we discuss the construction, statistical modeling and detection algorithm of a new type of ion channel based biosensor. The biosensor constructed comprises engineered dimeric gramicidin A (bis-gA) ion channels incorporated into a lipid bilayer membrane, supported over a 1-micron diameter opening of a micropipette, which was excised from a giant lipid vesicle. In a giant lipid vesicle, covalent dimeric gramicidin A ion chan- nels were incorporated by codispersion with the vesicle forming lipids. This type of artificially constructed biosensor mimics the naturally occurring ion transport processes of a living cell. The biosensor is constructed by Bruce Cornell in AMBRI Ltd. and Donald Martin and his students in the Department of Medical and Molecular Biosciences in the University of Technology, Sydney. Having constructed the biosensor, we formulate a stochastic dynamical model to capture its experimental behaviour. The engineered dimeric gA channels provided the conducting pore for the vesicle membrane, but not with the kinetics of the ICS™sensor discussed in Sec. 2.3. Instead, the gating mechanism of the dimeric gA ion channels in this biosensor is thought to arise from random movement of excess lipid lenses in the liposome that 40 Chapter 4. Bis-gA Channel Biosensor diffuse over the membrane surface and block the conducting channels. In experiments we found that a hidden Markov model (HMM), which takes into account of the 1/f noise in the biosensor's response, is an adequate model for the biosensor currents. In the presence of target molecules, the stochastic behaviour of the biosensor current changes. By using a sequential maximum likelihood detector, we show that the biosensor can be used in real-time target molecule detection. We illustrate the use of the biosensor in detecting two types of target molecules, methylbenzthonium chloride (MBC) and 2-methyl-4tert.-octylphenol (MTOP). The experimental data shows that the detection algorithm performed remarkably well even when it was difficult to visually identify the model change. Thus the system we present in this chapter provides a sensitive platform for the development of artificially constructed biosensors that better mimic the function of living cells. The remainder of this chapter is organized as follows. Sec. 4.1 de- scribes the construction of the biosensor. Sec. 4.2 presents a hidden Markov model to describe the dynamical behaviour of the biosensor. We also present model valididation methods to verify the goodness-of-fit of the H M M to the biosensor response. In Sec. 4.3, we discuss sequential detection algorithms for detecting target molecules in real-time; and finally, in Sec. 4.4 we demonstrate the experimental detection of two target moecules, M B C and M T O P , both of which are known to inhibit the conduction of bis-gA channels. 41 Chapter 4. Bis-gA Channel Biosensor 4.1 Experimental Construction of bis-gA Channel Biosensor The biosensor considered in this chapter was constructed by incorporating bis-gA ion channels into the lipid bilayer membrane of giant unilamellar liposomes and then excising small patches (1 ^im in diameter) of the lipid membrane using a patch-clamp micropipette. T h e bis-gA was synthesized at the A m b r i laboratories. Figure 4.1 shows the fluorescence and phase- contrast image of the optical section through the diameter of the biosensor. The solutions and chemicals used for the model biosensor included DL-alphaphosphatidylcholine from soybean, cholesterol, chloroform, sucrose, glucose, sodium chloride, potassium chloride and 4-(2-Hydroxyethyl)piperazine-lethanesulfonic acid. Figure 4.1: (a) Fluorescence image of biosensor's horizontal optical section shows the bis-gA channels labelled using fluorescein and identified by the green color, (b) Phase-contrast image of the same horizontal slice shows the overall shape of the biosensor. 42 Chapter 4. Bis-gA Channel Biosensor Giant unilamellar liposomes were prepared from a standard hydration procedure w i t h minor modifications. Bis-gA channels were incorporated into lipid membrane patches excised from the giant unilamellar liposomes. A patch-clamp pipette with a tip-opening between 0.9 and 1.5 /Ltm was used to perform the excision and to record the ionic currents from the model biosensor that resulted from the permeation of ions through the bis-gA ion channels. The recording from these membrane patches was amplified and filtered at 1 kHz (4-pole Bessel) using an Axopatch 200B amplifier (Axon Instruments) and sampled online at 10 k H z . Figure 4.2 shows a schematic of the experimental setup. 1 um contact area between glass micropipette and liposome Figure 4.2: Block diagram of experimental setup with photo of glass micropipette and liposome 43 Chapter 4. Bis-gA Channel Biosensor 4.2 Statistical Modeling and Validation of Biosensor Current In this section we discuss the use of hidden Markov modeling to describe the biosensor response and introduce estimation techniques to extract the model parameters. Statistical model validation tests are used to verify the goodness-of-fit of the model. Suppose a patch clamp experiment is conducted with (TV — 1) gramicidin A channels in the biosensor. At discrete time k, each bis-gA channel can be either in the "open" or "closed" state and each open channel conducts a fixed current. Thus the total current due to all (TV — 1) ion channels at any given time can take on one of N possible levels { ^ i , . . . ,PN} be modeled as a N-st&te Markov chain. Write as p = (u\,..., and can HN)- Let I k denote the total channel current at discrete time k. Let a ij = P(I k = ij\I _ =p ),i,je{l,...,N} f k 1 (4.1) H denote the transition probabilities of the Markov chain. Also let no(i) = P{I =(H),ie{l,...,N} (4.2) 1 denote the initial distribution of the Markov chain. Write A — [aij] . N> 7i"o = bToWljvxi- N and The measured current from the biosensor is a distorted version of the signal I . k The distortion arises from thermal noise as well as open channel noise that has its power proportional to the inverse of frequency. Thus this is also known as 1/f noise and is discussed in other 44 Chapter 4. Bis-gA Channel Biosensor studies of bis-gA ion channels. [36, 37]. Figure 4.3 shows the power spectral density of a typical sequence of biosensor recordings. 50,— — •40 10000 1000 — — 1 1 100 Frequency (Hz) 10 Figure 4.3: Power spectral density of biosensor response clearly shows the 1/f noise and anti-aliasing effect. A t low frequencies, the 1/f noise process is evident as the power spectrum decreases at a rate of -10 d B / d e c . The effect of the anti-aliasing filter causes a sharp cutoff at approximately 1 k H z . To model this correlated noise process, we can use an auto-regressive ( A R ) Gaussian process that comprises white Gaussian noise process Wk filtered by an all-pole filter. We represent the filter with transfer function H(q~ ), where q~ denotes the unit delay 1 l operator: Y = Ik + k tffcr ) 1 (4.3) 45 Chapter 4. Bis-gA Channel Biosensor or equivalently Hiq-^Yk = H{q- )I + W l k (4.4) k where H(q- ) 1 = 1 + hxq' + ... + h q~ 1 (4.5) M M The model in (4.4) can be rewritten as: h Y = h I + W r where Y = {Y , F/b-i, • • •, Y -M) h = (1, hi,..., k k k k T k (4.6) T f c and I k = (h, h-i, • • •, 4 - M ) T Let / I M ) - The standard algorithm to estimate the model paT rameters in (4.6) is the expectation-maximization (EM) algorithm, which will be explained in detail in Sec. 4.2.1. The key point is that the E M algorithm involves computing the marginal probabilities of the state space P(Ik). The size of the state space increases exponentially with M. This leads to huge complexity issues because 1/f noise is a long memory process and modeling it requires choosing the filter order, M, to be large. To alleviate the computational complexity, we adjust the model so that the size of the state space is indepedent of the filter order M. We observe that in our current setup, the gating mechanism of the bis-gA channels is a much slower process relative to the sampling rate and most of the power spectrum of the channel resides in low frequencies. Therefore, we argue that 46 Chapter 4. Bis-gA Channel Biosensor we can redefine the response of the noisy observation as: Y k = I k + W (47) k and closely approximate the behaviour of the gramicidin channels with our model in (4.7), provided that the filter i J ( g ) has unity D C gain. We apply _1 model validation methods in Sec. 4.2.1 to evaluate the performance of our model. It is convenient to model the noise corrupting the state of the biosensor as state dependent noise - that is the noise variance at any given time instant is dependent on the state of the biosensor at that time instant. Let o\ be the variance of state i, i — (1,..., N). Write o = (cr ,..., o\f). As a result, the 2 2 observations can be formulated as a hidden Markov model sequence. Let 9 = (A, TT, u, a , h) be the H M M that characterizes the output measured 2 current from the biosensor. 4.2.1 P a r a m e t e r E s t i m a t i o n for B i o s e n s o r C u r r e n t Given an observation sequence {Y } of length T, we define L (9) as the k k log-likelihood of our model at discrete time k. The estimation of the model 0 involves processing {Y } through a H M M maximum likelihood estimator k (MLE). The system in (4.7) can be rewritten as: h Y T f c = I k + W k (4.8) This formulation is analogous to a standard H M M , except that the ob47 Chapter 4. Bis-gA Channel Biosensor servation sequence is modulated by a finite impulse response (FIR) filter. As mentioned earlier, the E M algorithm is an iterative procedure that solves for local maximum of the likelihood function. The E-step uses the forwardbackward recursions to compute the forward probability a k and backward probability f3 : k a {i) = P(Y ,Y2,...,Y ,I k 1 Pki}) = P{Y i, k Y k+ = m\0) k ,..., Y , h = m\0) k+2 T and evaluates the log-likelihood k N N + £££Ct(M')logK-) t=l 1=1 (4.9) j=l where the probabilities ^ and ( k k can be computed from the forward and backward probabilities. The M-step maximizes the log-likelihood with respect to the model parameters as follows: Efc=iCfc(».j) ELi7fcW ELi7fc(i)(h Y - ) ELi7fc(«) ELi7fc(Qh Y r k 2 w r k ELi7fc(») 48 Chapter 4. Bis-gA Channel Biosensor The coefficients of the filter H{q~ ) can also be estimated by taking firstl order derivatives of the likelihood function. Alternatively, the estimation can be implemented with the Yule-Walker type equations. For details about H M M estimation and the E M algorithm, please see [12, 16, 30, 34]. 4.2.2 Statistical Validation of Biosensor M o d e l Statistical model validation is of key importance for the biosensor since, once a satisfactory stochastic model is determined, an appropriate molecular detection algorithm can be constructed. In Sec. 4.2.1, we discussed model estimation techniques with the assumption that the N-state hidden Markov model with an A R filter is an approriate model for our biosensor measurements. In this section, we need to statistically validate that assumption. Model validation can be done by analyzing the autocorrelations of the residuals, which are generated via a H M M one-step predictor: N e/c|fc-i = V fc N M a a -^2Y^ iJ k-i{i)^j - ^KY J =l 71=1 1=1 k _ n (4.10) The lag-/ autocorrelation function of the residual is defined as: r(0 _O'-*-^-) Efc=i(e/c-e) (4 . n) 2 where e is the mean of the residual process. The residuals of an adequately fitted model should be uncorrelated and the autocorrelation should approach 0 as T —> oo. Rather than examining the autocorrelation at each lag /, 49 Chapter 4. Bis-gA Channel Biosensor the standard procedure is to compute the Ljung-Box Q-statistic defined in (4.12), which computes the cumulative sum of autocorrelations at the first L lags. T h e Ljung-Box test is used as a portmanteau lack of fit test for model adequacy. [7] (4.12) i=i v ' It is shown in [8] that for an adequate model, the Q-statistics of the residual is approximately distributed as x (^)2 We report on the validity of the model on the biosensor experimental data in Sec. 4.4. 4.3 Target Molecule Detection Algorithm After the design and validation of the dynamical models, in this section we discuss algorithms for the detection of analytes that are known to change the statistics of bis-gA channel conductance. M B C and M T O P are two analytes that are known to inhibit conductance of bis-gA channels. T h e chemical structure of M B C and M T O P are shown in Figure 4.4 and 4.5 respectively. The interaction of M B C and M T O P with g A channels is described in the manuscript [19]. 50 Chapter 4. Bis-gA Channel Biosensor i" i V A 3 1+ 3 Ch^-C-CHj-C^ V O C H , Ch-12 OCH, C H , - N - C H , CH, CH3H C H 3 _j ~\ r y Cl- Figure 4.4: Chemical structure of Methylbenzthonium Chloride H C 3 CH 3 Figure 4.5: Chemical structure of 2-Methyl-4-tert.-octylphenol Given a measured sequence observed in an unknown condition, the detection problem involves the identification of the condition that most likely contributes to the biosensor's response. This is a model classification problem and can be solved by comparing the likelihood of each known model. Let Y = (Y\,..., Yp) be a sequence of observed responses of the biosensor and let 0 = . . . , 6 M) denote the set of model parameters that charac- terizes the biosensor's response for known conditions. It is assumed that at each time point k, the sequence Y behaves according to one of M possible models in 0. The model parameters 0 are estimated and the log-likelihood at each time point Lh can be computed from (4.9). To make the detection more robust to nonstationary disturbances and outliers in the measurements, we apply a 51 Chapter 4. Bis-gA Channel Biosensor geometric moving-aver age filter to the log-likelihood. Let p be the forgetting factor, 0 < p < 1. Define the filtered likelihood of model 9 at time k: S (9) ={ LAB) forfc = l k (4.13) (1 - p)S -! + pL {9) for 2 < k < T k k The filtered likelihood is a weighted sum of the likelihood of the entire sequence {Y\,... ,Y }, with higher weights on the recent observations. The k apriori probabilities of each model is generally unknown, so we devise a maximum-likelihood detector which picks the most likely model at time k given the measured sequence {Yi,..., Y }: k c9 = argmax(5 (0)) fc fc (4.14) where 9 is the maximum likelihood detection of the model at time k. k 4.4 Experimental Results of Biosensor Here we report on the goodness-of-fit of the dynamical model and the performance of the detection algorithm on experimental data. We recorded output from the biosensor by measuring the activity of the bis-gA ion channels that were incorporated into the small lipid membrane patches excised from unilamellar giant liposomes. Bis-gA ion channels were incorporated into the unilamellar giant liposomes at a concentration of 1/100 from a 66 nM stock solution. We used 0.5 M KC1 solution in the recording pipette 52 Chapter 4. Bis-gA Channel Biosensor with the micro-particles suspended in a 0.5 M NaCl solution. 4.4.1 M o d e l E s t i m a t i o n a n d V a l i d a t i o n for B i o s e n s o r We recorded the biosensor response at an applied voltage of 50 mV and an amplifer gain of 200. Figure 4.6 shows a three-second recording of biosensor response. ii 0 1 . 1 0.5 1 1.5 Time (s) , 2 , 2.5 1 3 Figure 4.6: Experimental measurement of biosensor response was recorded at an applied voltage of 50 mV and an amplifer gain of 200. The giant liposome was constructed from P C with 10% cholesterol in chloroform. Notice that the bis-gA current measurements obtained here are more noisy compared to the ones obtained from a solvent sealed black lipid membrane (BLM). The quality is typical of patch clamp recordings in which a membrane patch is captured on the tip of a 1 fim pipette. The seal obtained using the patch clamp technique is generally poorer than that obtained from a B L M , resulting a worse signal-to-noise ratio for the measurements. 53 Chapter 4. Bis-gA Channel Biosensor Transition Probabilities 0.9964 0.0006 Conductance Level (PA) Variance (PA) 5.831 3.057 0.513 0.371 0.0036 0.9994 2 Table 4.1: Maximum Likelihood Estimation of H M M parameters of the biosensor recording in Figure 4.6. A further potential difficulty is the possibility of multiple lipid layers being present in the patch as compared to the B L M geometry which is thinned to a single bilayer. When recording from multiple lipid membrane stacks the signal-to-noise ratio of the individual conduction events is reduced by noise arising from the complex impedance of the series elements. We fitted the sequence with two-state H M M and a 12~ th order A R filter. The model parameters were estimated with the M L E and listed in Table 4.1. The most likely conductance level sequence is extracted from the H M M procedure and plotted in Figure 4.7. To verify that the two-state H M M outlined in Table 4.1 provides sufficient statistics to model the response of the biosensor, we computed the Q-statistics of the residuals, using (4.10) through (4.12). It can be seen in Figure 4.8 that for the first 13 lags, the Q-statistics are below the critical values of the chi-square distribution at 0.05 significance level. Therefore, the hypothesis that the residual is a white process cannot be rejected. It is interesting to note that the estimated state levels, at approximately 3 pA and 6 pA, are rather large compared to single channel measurement of channel current at a low applied voltage. One plausible explanation for this observation is that the activity of some of the bis-gA channels are correlated. In other words, the level switch from 3 pA to 6 pA may be contributed by 54 Chapter 4. Bis-gA Channel Biosensor 5.86 < Q. C ID fc o 3.05 I 0 i 0.5 i 1 , 1.5 Time (s) i 2 i 2.5 I 3 Figure 4.7: Maximum likelihood estimate of biosensor conductance was extracted from the forward-backward procedure in the H M M estimation algorithm. 90 Lags Figure 4.8: Q-statistics of residuals are plotted against the critical values of the Ljung-Box test at significance level of 0.05. For the first 13 lags, the Q-statistics are below the critical values and the hypothesis that the residual is a white process cannot be rejected. Thus the proposed two-state H M M is an adequate model for the biosensor current. 55 Chapter 4. Bis-gA Channel Biosensor several coupled channels opening at the same time. 4.4.2 R e a l - T i m e D e t e c t i o n o f Target M o l e c u l e s In this section we illustrate the performance of the algorithm in detecting the presence of target molecules in real-time. We present experimental results of the biosensor i n detecting two analytes, M B C and M T O P . B o t h compounds are known to inhibit conduction of the bis-gA channels. Patch-clamp experiments are conducted with and without M B C in the bath solution. Let 9\ = w i t h M B C in the bath solution #2 = with no M B C i n the bath solution 110 Time (s) Figure 4.9: Biosensor's response to addition of M B C at fc = 57.828 seconds. To simulate the addition of M B C into the bath solution, we merged together sequences recorded with and without M B C . The merged sequence is plotted in Figure 4.9. Assuming a four-state H M M , the parameters of 9\ 56 Chapter 4. Bis-gA Channel Biosensor and #2 are estimated with the E M algorithm. The filtered likelihood of each model is shown in Figure 4.10. Since the stochastic model successfully captures the dependencies in the biosensor's response, the detection algorithm quickly and accurately estimates the model switching point. The detection trace in Figure 4.11 indicates a switch in the most likely model from 9\ to 62 at k — 59.865 seconds, approximately 0.2 seconds after M B C is added. Time (s) Figure 4.10: Filtered likelihoods for estimated models of M B C , with p = 0.0001. In the second example, we test the algorithm's performance in identifying the concentration of M T O P in the bath solution. Let 01 = no M T O P in the bath solution 0 = 50 uM M T O P in the bath solution 2 0 = 100 uM M T O P in the bath solution 3 We simulated the change in concentration of M T O P by merging channel 57 Chapter 4. Bis-gA Channel Biosensor Figure 4.11: M a x i m u m likelihood detection of M B C is generated using the likelihoods i n Figure 4.10. While it is difficult to visually identity the transition from the recordings, the detection algorithm detected a model change shortly after M B C was added to the solution. 12 Time (s) Figure 4.12: Biosensor's response to various concentrations of M T O P . 50 / / M of M T O P was added to the bath solution every 50 seconds. 58 Chapter 4. Bis-gA Channel Biosensor recordings independent from the training sequences. The merged sequence is plotted in Figure 4.12. The parameters of Q\, 62 and #3 were estimated offline with the E M algorithm. The filtered likelihood and the detection trace are plotted i n Figure 4.13 and Figure 4.14 respectively. The algorithm identified the correct concentration of M T O P 98.1% of the time. 2.4 2.2 50 |iM - ""0 50 100 100 uM - 150 Time (s) Figure 4.13: Filtered likelihoods for estimated models of M T O P , with p — 0.0001. 59 Chapter 4. Bis-gA Channel Biosensor e ,| 0 1 u 50 100 150 Time (s) Figure 4.14: Maximum likelihood detection of M T O P generated using the likelihoods in Figure 4.13, showing the most likely condition. The algorithm detected the correct model 98.1% of the time. 60 Chapter 5 Conclusion and Extensions This thesis is concerned with the statistical modeling and application of gramicidin A ion channels. In particular, we have addressed two important issues in the area of ion channel permeation and ion channel based biosensor. Brownian dynamics is a practical technique to simulate ion channel conduction but it suffers from huge computational complexity at high ionic concentration. We proposed the use of a finite-state semi-Markov model to describe the binding site kinetics of gramicidin A channels. We extrapolated the model parameters with the logistic function and used the results to predict the channel current at high ionic concentrations. Numerical results showed that the technique correctly predicted the channel current and conduction event distribution. In future work, we will use the same modeling and extrapolation technique to study single channel current as a function of applied voltage. In a separate study, we developed stochastic modeling and detection algorithms as part of an ion channel based bionsensing platform that can automatically detect a variety of target molecules. The conduction of an ion channel biosensor can be modeled by an AR-modulated hidden Markov model. The presence of specific target molecules distorts the channels in 61 Chapter 5. Conclusion and Extensions the biosensor, thus changing the statistics of the biosensor's conduction. We devise a real-time maximum likelihood detector to capture the change and detect the presence of target molecules. The algorithm is tested on a biosensor built by incorporating dimeric gramicidin A channels into bilayer membranes of giant unilamellar liposomes. The performance of the biosensor is tested with addition of M B C and M T O P , both known to inhibit conduction of gramicidin A channels. Experimental results show that the detection performed well even when the change in the conduction level was difficult to visualize. The detection algorithm provides the sensitive detection system for ongoing development of membrane-based biosensors. In future work, we will construct a more robust biosensor by incorporating bis-gA channels into a B L M , where we can observe the activity of a single bis-gA channel and investigate in more detail biological properties of the bis-gA channels. 62 Bibliography [1] T . W. Allen, T . Bastug, S. Kuyucak, and S. H. Chung. Gramicidin A channel as a test ground for molecular dynamics force fields. Biophysical Journal, 84:2159-2168, 2003. [2] A. Anastasiadis and F. Separovic. Solid-state N M R structural determination of components in an ion channel switch biosensor. Australian Journal of Chemistry, 56:163-166, 2003. [3] F. M . Ashcroft. Ion Channels and Disease: Channelopathies. Academic Press, 2000. [4] T . Bastug, S. M . Patra, and S. Kuyucak. Molecular dynamics sim- ulations of gramicidin A in a lipid bilayer: From structure-function relations to force fields. Chemistry and Physics of Lipids, 141:197-204, 2006. [5] Turgut Bastug and Serdar Kuyucak. Molecular dynamics simulation of calcium binding in gramicidin. Chemical Physics Letter, 424:82-85, 2006. [6] S. Berneche and B. Roux. Molecular dynamics of the KcsA K + channel in a bilayer membrane. Biophysical Journal, 78:2900-2917, 2000. 63 Bibliography [7] G. E . P. Box and G. M . Jenkins. Time Series Analysis - Forecasting and Control. Holden-Day, 1976. [8] G. E . P. Box and D. A. Pierce. Distribution of residual autocorrelations in autoregressive-integrated moving average time series models. Journal of American Statistical Association, 65:1509-1526, December 1970. [9] O. Braha, L. Gu, L. Zhou, X. Lu, S. Cheley, and H. Bayley. Simultaneous stochastic sensing of divalent metal ions. Nature Biotechnology, 18:1005-1007, 2000. [10] M . Chakravarti, R. G. Laha, and J. Roy. Handbook of Methods of Applied Statistics, volume 1. John Wiley and Sons, 1967. [11] S. H . Chung, O. S. Andersen, and V. Krishnamurthy. Handbook of Ion Channels: Dynamics, Structure, and Applications. Springer Verlag, 2006. [12] S. H. Chung, V. Krishnamurthy, and J. B. Moore. Adaptive processing techniques based on hidden Markov models for characterising very small channel currents buried in noise and deterministic interferences. Proc. Phil. Trans. Roy. Soc. Lond. B, 334:357-384, 1991. [13] B. A. Cornell, V. L. Braach-Maksvytis, L. G. King, P. D. Osman, B. Raguse, L. Wieczorek, and R.J. Pace. A biosensor that uses ion-channel switches. Nature, 387:580-583, 1997. 64 Bibliography [14] B. Corry, T. Allen, S. Kuyucak, and S. H. Chung. Mechanisms of permeation and selectivity in calcium channels. Biophysical Journal, 80:195-214, 2001. [15] R. B. D'Agostino and M . A. Stephens. Goodness-of-fit techniques. CRC, 1986. [16] A. P. Dempster, M . M . Laird, and D. B. Rubin. Maximum likelihood estimation from incomplete data via the E M algorithm. Journal of the Royal Statistical Society B, 39:1-38, 1977. [17] D. A. Doyle, J. M . Cabral, R. A. Pfuetzner, A. Kuo, J. M . Gulbis, S. L. Cohen, B. T. Chait, and R. Mackinnon. The structure of the potassium channel: molecular basis of K + conduction and selectivity. Science, 280:69-77, 1998. [18] R. Dutzler, E . B. Campbell, M . Cadene, B. T. Chait, and R. MacKinnon. X-ray structure of a CIC chloride channel at 3.0 A reveals the molecular basis of anion selectivity. Nature, 415:287-294, 2002. [19] B. A. Cornell et al. Making lipid membranes even tougher. Journal of Material Research, 2007. In Press. [20] B. Hille. Ionic Channels of Excitable Membranes. Sinauer Associates, 2001. [21] S. Howorka, S. Cheley, and H. Bayley. Sequence-specific detection of individual D N A strands using engineered nanopores. Nature Biotechnology, 19:636-639, 2001. 65 Bibliography [22] Matthew Hoyles. Computer simulation of biological ion channels. PhD Thesis, 2000. [23] Matthew Hoyles, Serdar Kuyucak, and S. H. Chung. Computer simulation of ion conductance in membrane channels. Computational Phys. Rev. E. of and Theoretical Nanoscience, 58:3654-3661, 1998. [24] E . Jakobsson. Using theory and simulation to understand permeation and selectivity in ion channels. Methods, 14:342-351, 1998. [25] R. R. Ketchem, W. Hu, and T. A. Cross. High-resolution conformation of gramicidin A in a lipid bilayer by solid-state NMR. Science, 261:1457-1460, 1993. [26] V . Krishnamurthy and S. H. Chung. Brownian dynamics simulation for modeling ion permeation across bionanotubes. IEEE Nanobioscience, Transactions on 4:102-111, 2005. [27] V. Krishnamurthy, M . Hoyles, R. Saab, and S. H. Chung. Permeation in gramicidin ion channels by directly estimating the potential of mean force using Brownian dynamics simulations. Journal of Computational and Theoretical Nanoscience, 3:702-711, 2006. [28] V. Krishnamurthy, T. Vora, and S. H. Chung. Adaptive Brownian dynamics algorithms for estimating the shape of sodium ion channels. Journal of Nanoscience and Nanotechnology, 7:2273-2282, 2007. [29] Vikram Krishnamurthy, Matthew Hoyles, Rayan Saab, and S. H. Chung. Permeation in gramicidin ion channels by directly estimating 66 Bibliography the potential of mean force using Brownian dynamics simulations. Journal of Computational and Theoretical Nanoscience, 3:702-711, 2006. [30] Vikram Krishnamurthy and A. Logothetis. Iterative and recursive estimators for hidden Markov errors-in-variables models. IEEE Trans. Signal Processing, 44:629-639, March 1996. [31] David G. Levitt. Modeling of ion channels. Journal of General Physiology, 113:789-794, 1999. [32] L . Movileanu, S. Howorka, O. Braha, and H. Bayley. Detecting protein analytes that modulate transmembrane movement of a polymer chain within a single protein pore. Nature Biotechnology, 18:1091-1095, 2000. [33] G. A. Olah, H. W. Huang, W. H. Liu, and Y . L Wu. Location of ionbinding sites in the gramicidin channel by x-ray diffraction. Journal of Molecular Biology, 218:847-858, 1991. [34] L . R. Rabiner. A tutorial on hidden Markov models and selected applications in speech recognition. Proc. IEEE, 77:257-285, 1989. [35] K . E . Sapsford, Y . S. Shubin, J. B. Delehanty, J. P. Golden, C. R. Taitt, L. CI Shriver-Lake, and F. S. Ligler. Fluorescence-based array biosensors for detection of biohazards. Journal of Applied Microbiology, 96:47-58, 2004. [36] R. Sauve and E . Bamberg. 1/f noise in black lipid membranes induced by ionic channels formed by chemically dimerized gramicidin A. Journal of Membrane Biology, 43:317-333, November 1978. 67 Bibliography [37] R. Sauve and G. Szabo. Interpretation of 1/f fluctuations in ion conducting membranes. Journal of Theoretical Biology, 113:501-516, April 1985. [38] E . A. Smith and R. M . Corn. Surface plasmon resonance imaging as a tool to monitor biomolecular interactions in an array based format. Appl. Spectroscopy, 57:320-332, 2003. [39] W. F. van Gunsteren and H. J. C. Berendsen. Algorithms for Brownian dynamics. Molecular Physics, 45:637-647, 1982. [40] T . Vora, B. Corry, and S. H. Chung. A model of sodium channels. Biochimica et Biophysica Acta - Biomembranes, 1668:106-116, 2004. [41] B. A. Wallace. Recent advances in the high resolution structures of bacterial channels: Gramicidin A. Journal of Struct. Biol., 121:123141, 1968. 68 Appendix A List of Acronyms AR BD bis-gA BLM CDF EDF EM FIR gA HMM ICS™ MBC MD MLE MTOP NMR PCR PDF PMF SMC Auto-regressive Brownian dynamics dimeric gramicidin A black lipid membrane cumulative distribution function empirical distribution function expectation-maximization finite impulse response gramicidin A hidden Markov model Ion Channel Switch methylbenzthonium chloride Molecular dynamics maximum likelihood estimator 2-methyl-4-tert .-octylphenol nuclear magnetic resonance polymerase chain reaction probability distribution function potential of mean force semi-Markov chain
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Statistical modeling and application of gramicidin...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
Statistical modeling and application of gramicidin A ion channel Luk, Kai Yiu 2007
pdf
Page Metadata
Item Metadata
Title | Statistical modeling and application of gramicidin A ion channel |
Creator |
Luk, Kai Yiu |
Publisher | University of British Columbia |
Date Issued | 2007 |
Description | Ion channels are aqueous pores in the cell membrane for selected ions to flow down their electrochemical gradient. These channels play a prominent role in a variety of biological processes in the human body. Determining the structure and function of ion channels is of fundamental importance in biology. Also, the selective conductivity and specific gating mechanism of ion channels have attracted much interest in the area of artificial molecular detectors. Ion channel based biosensors are developed to detect molecular species of interest in medical diagnostics, environmental monitoring and general bio-hazard detection. This thesis is concerned with statistical techniques used to describe ion channel permeation and to develop ion channel based biosensors. Brownian dynamics is a popular technique to simulate ion channel permeation but is too computationally expensive to run when ionic concentration is high. By fitting binding site statistics of BD simulation to a semi-Markov chain, we obtain a simpler model with conduction properties that are statistically the same as the simulations. This approach enables the use of extrapolation techniques to predict channel conduction when performing the actual simulation is computationally infeasible. Numerical studies on the simulation of gramicidin A channels are presented. In a separate study, we show the use of statistical modeling and detection techniques as part of a sensitive biosensing platform. A nano-scale biosensor is built by incorporating dimeric gramicidin A channels into bilayer membranes of giant unilamellar liposomes. The presence of specific target molecules changes the statistics of the biosensor's conduction. By capturing the change in real time, we devise a maximum likelihood detector to detect the presence of target molecules. The performance of the biosensor is tested with the addition of various target molecules known to inhibit conduction of gramicidin A channels. Experimental results show that the detection performed well even when the conductance change was difficult to visualize. The detection algorithm provides a sensitive detection system for ongoing development of membrane-based biosensors. |
Genre |
Thesis/Dissertation |
Type |
Text |
Language | eng |
Date Available | 2011-03-03 |
Provider | Vancouver : University of British Columbia Library |
Rights | For non-commercial purposes only, such as research, private study and education. Additional conditions apply, see Terms of Use https://open.library.ubc.ca/terms_of_use. |
DOI | 10.14288/1.0101099 |
URI | http://hdl.handle.net/2429/31984 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
Campus |
UBCV |
Scholarly Level | Graduate |
Aggregated Source Repository | DSpace |
Download
- Media
- 831-ubc_2007-0490.pdf [ 3.85MB ]
- Metadata
- JSON: 831-1.0101099.json
- JSON-LD: 831-1.0101099-ld.json
- RDF/XML (Pretty): 831-1.0101099-rdf.xml
- RDF/JSON: 831-1.0101099-rdf.json
- Turtle: 831-1.0101099-turtle.txt
- N-Triples: 831-1.0101099-rdf-ntriples.txt
- Original Record: 831-1.0101099-source.json
- Full Text
- 831-1.0101099-fulltext.txt
- Citation
- 831-1.0101099.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0101099/manifest