BROWNIAN DYNAMICS SIMULATION OF A FIVE-SITE MODEL FOR A MOTOR PROTEIN ON A BEAD-SPRING SUBSTRATE. A Thesis. Presented to
|
|
- Joy Jacobs
- 5 years ago
- Views:
Transcription
1 BROWNIAN DYNAMICS SIMULATION OF A FIVE-SITE MODEL FOR A MOTOR PROTEIN ON A BEAD-SPRING SUBSTRATE A Thesis Presented to The Graduate Faculty of The University of Akron In Partial Fulfillment of the Requirements for the Degree Master of Science Nabina Paudyal December, 2014
2 BROWNIAN DYNAMICS SIMULATION OF A FIVE-SITE MODEL FOR A MOTOR PROTEIN ON A BEAD-SPRING SUBSTRATE Nabina Paudyal Thesis Approved: Accepted: Advisor Dr. Jutta Luettmer-Strathmann Dean of the College Dr. Chand Midha Faculty Reader Dr. Yu-Kuang Hu Dean of the Graduate School Dr. Rex D. Ramsier Faculty Reader Dr. Alper Buldum Date Department Chair Dr. David N. Steer ii
3 ABSTRACT Motor proteins play an important role in many biological processes. For example, kinesin molecules are responsible for the transport of vesicles in nerve cells and their malfunction has been linked to neurodegenerative diseases. Unfortunately, the complexity of motor proteins and their environment makes it difficult to model the detailed dynamics of molecular motors over long time scales. In this work, we develop a simple coarse-grained model for a motor protein on a bead-spring substrate under tension. In our model, different pair potentials describe interactions between substrate and motor, motor components and substrate components. The movement of motor proteins entails ATP hydrolysis, which is modeled in terms of mechano-chemical states that couple positional and chemical degrees of freedom. We have identified model parameters corresponding to protein and chain conformations of a walking protein as well as transition rates for the states and apply our model to simulate a motor protein walker with cargo. iii
4 ACKNOWLEDGEMENTS I would like to express my sincere gratitude to my advisor Dr. Jutta Luettmer Strathmann for her valuable support and guidance throughout this journey. I would also like to thank my dear friend Maral Adeli Koudehi for her priceless help and company. I am thankful to Dr. Hu and Dr. Buldum for being on my committee and providing me valuable advice and encouragement. Special thanks to my research group and my friends Binod Dhakal, Kiran Khanal, Shavyata Acharya, Dipika Parajuli, Chinta Mani Aryal, and Sajeevi Withanage for their valuable guidance. I am grateful to the Department of Physics and all my colleagues for their direct and indirect support during my study and research. Finally I would like to dedicate all my works and my thesis to my parents and my brother Nabin who have always believed in me and have always been there for me to lift me up. iv
5 TABLE OF CONTENTS Page LIST OF TABLES LIST OF FIGURES vii viii CHAPTER I. INTRODUCTION Outline of thesis II. MODEL Conservative Forces Mechanical State Chemical Steps III. METHOD Brownian Dynamics Simulation Behavior of an isolated protein Conformation of the protein in contact with the substrate Kinetic energy IV. RESULTS Force-extension relation of the substrate Motor protein taking steps on the substrate without cargo v
6 4.3 Effect of changing the interaction of the attached head of a motor protein walking without cargo Motion of protein with cargo Five steps for a protein with cargo V. CONCLUSIONS BIBLIOGRAPHY vi
7 LIST OF TABLES Table Page 2.1 Model units of physical quantities Interaction parameters Interaction parameters for cargo Parameters for the equilibrium states Assigned parameters for the change of states to maintain the desired amplitudes vii
8 LIST OF FIGURES Figure Page 2.1 Coarse grained model of a motor protein on a bead spring substrate. Sites 1, 2, 3 and 4 represents protein heads and 5 represents the central site, which is bonded to the light chain domain for the attachment of cargo structures Graphs of spring potentials as a function of distance for three sets of parameters that describe bonds between different types of particles: BB for bead-bead interactions along the substrate chain, HH for bonds between protein sites in the same head domain, HC for bonds between protein head sites and the central protein sites; values for the parameters are presented in Table Graphs of four interparticle potentials as a function of distance between the particles. The solid lines represent the bead-bead harmonic spring (magenta), soft repulsive Lennard-Jones (black), Lennard-Jones (blue), and Morse (red) potentials, respectively. The parameters of the potentials are shown in Table Illustration of the generalized dihedral planes of our motor protein model. The dihedral angle α is the angle between the normal unit vectors of the planes Graph of the dihedral potential as a function of the angle α for three different dihedral amplitudes V Illustration of the bond angle θ formed by two unit vectors â and ˆb along the bonds between the central protein site and head sites 2 and Graph of the bond angle potential as a function of the bond angle θ. The equilibrium bond angle is θ 0 = in our model Illustration of the interactions between the components of the protein.. 16 viii
9 2.9 Illustration of the interactions between protein and substrate sites. Solid red arrows represent the Morse interaction between protein heads and the substrate, whereas dotted black arrows represent the soft repulsive Lennard-Jones interaction between the central protein and the substrate Illustration of the interactions between substrate sites. The heavy lines connecting the beads represent harmonic springs; the curved arrows indicate Lennard-Jonnes interactions between non bonded beads. An external tension force, represented by F pull, is applied to the end beads of the chain with equal magnitude and opposite directions Chemical state as a function of time for a two-step simulation. The initial state, S 1 = 1, S 2 = S 3 = S 4 = 0 is maintained during an equilibration period at the start of the simulation. The probability S 1 to be in state 1 (red line) falls off rapidly as S 2 (blue line) first increases. As S 2 decreases S 3 (purple line) first increases and then decreases as S 4 starts to grow. At the end of the first step, S 1 =S2 =S2 = 0 and S1 = 1. In the second step the system goes from S 4 = 1 to S 1 through the states S 2 and S Effect of K θ and V 0 on conformations of an isolated protein Illustration of the mechano-chemical states (1, 2, and 3) of a protein walker on a substrate. The potential parameters of the states are presented in Table Instantaneous (red) and segment-averaged (blue) kinetic energy as a function of time for simulations of the protein-vesicle-chain system, which has a total of N = 26 particles. Panel A and B show results for fixed and time-varying interaction parameters, respectively. The length of one time segment in the averaging process is steps, corresponding to a simulation time of 500. In both simulations, the average kinetic energy has a value of 39, in agreement with the equipartition theorem and a reduced temperature of one Force extension results for an isolated chain. The symbols represent simulation results for the extension, defined in equation (4.1), as a function of the pulling force for a chain of length N c = 20. The error bars indicate the standard deviation of block averages (see equation (3.16)). Two chain conformations for intermediate (left) and high (right) extensions are also shown ix
10 4.2 Two steps of a motor protein walker on a chain substrate. Panel A shows the smallest distance between the chain and head site 1 on head one (d 1, red) and the chain and head site 3 on head two (d 3, pink) together with the indices of the beads closest to the head sites (C 1, blue and C 3, cyan). Panel B shows the value of the Morse potential parameters D e1 and D e3 as a function of time Two steps of a motor protein on a substrate with increased attraction for the attached head site. Panel A shows the smallest distance between the chain and head site 1 on head one (d 1, red) and the chain and head site 3 on head two (d 3, pink) together with bead indices closest to head site 1 (C 1, blue) and head site 3 (C 3, cyan). Panel B shows the value of the Morse potential parameters D e1 and D e3 as a function of time Two steps of a motor protein walking with cargo on a substrate. The graph shows the smallest distance between the chain and head site 1 (d 1, red), head site 3 (d 3, pink) and the cargo (ves, orange). The change in Morse parameters is as shown in panel B of graph 4.2. The closest beads are C 1 (blue), C 3 (cyan) and C ves (yellow) for head sites 1, 3, and the cargo, respectively. In these simulations, the Morse parameters of the non-released head remained constant, as shown in panel B of figure Two steps of a motor protein walking with cargo. The lines represent the same quantities as in figure 4.4; the Morse amplitudes for this simulation are shown in figure 4.3 B Results for the same system as in figure 4.5 from a separate simulation step simulation with the new protocol. The lines represent the same quantities as in figure x
11 CHAPTER I INTRODUCTION A machine is called a motor if its output involves mechanical movement [1, 2]. Molecular motors play a role in complex processes such as cell motility, cell division, cell development and also the transport of biomolecules [3]. Axonal transport is the transport of independent cargo structures with the help of motor proteins in the axon of nerve cells [4]. To ensure proper functioning of the nerve cells, cargoes should be delivered to the correct destination in appropriate quantity. Failure of such transport has been linked to several neurodegenerative diseases [4]. The cargo structures most often detected are vesicles, mitochondria, proteins and other cellular structures. In general, cargoes may be membranous or non membranous. However, the ones found in axons are typically vesicles, which consist of a lipid membrane enclosing the molecules to be transported. Microtubules and microfilaments are the most common parts of the cytoskeleton present in axons of nerve cells. They act as tracks for axonal transport. Motor proteins such as kinesin and dynein are found on microtubules and microfilaments, while the motor protein myosin interacts with actin filaments in muscles [4, 5, 6]. All motor proteins consist of heavy and light chain domains. The heavy chain domains have a globular head that interacts with the substrate, whereas the light chain domains attach to cargo and other structures. Although motor pro- 1
12 teins are very small, they are very complex structures and interactions between motor proteins and their substrate involve multiple degrees of freedom [7]. Motors operate by adjusting their conformation in response to chemical changes that fuel the motor proteins. The movement of motor proteins is fueled by ATP (adenosine triphosphate) hydrolysis. To describe how a chemical reaction generates mechanical motion of a motor protein, a number of theoretical models have been developed [7, 8, 9, 10, 11, 12]. Except for Atzberger and Peskin [12], who developed a three-dimensional model for kinesin consisting of two sites in contact with a hard, patchy surface, these are onedimensional mechano-chemical models. In our modeling of a motor protein we consider four issues: (i) degrees of freedom, (ii) forms of interaction, (iii) dynamical equations and (iv) methods of solutions. In this work, we develop a coarse-grained three-dimensional model of a motor protein that transports cargo by walking along a bead spring substrate. In our model, different pair potentials describe interactions between substrate and motor, motor components and substrate components. To couple the mechanical degrees of freedom to chemical states, we learn from the one-dimensional model of Lan and Sun [11]. Lan and Sun investigated the dynamics of the motor protein myosin and the simultaneous muscle contraction [11]. They used a simple one-dimensional model of a single myosin structure. The conformation of myosin was described by the angle between the light chain domain and the motor domain. Mathematically, the myosin 2
13 dynamics is described by the following coupled Langevin equations η θ = E(θ, S) F + f b (t) (1.1) θ s = K(θ) t S (1.2) where η is the friction coefficient, θ is the angle between light chain domain and motor domain, E(θ, S) is the elastic energy of myosin, F is an external force, f b (t) is the Brownian random force, S is the chemical state of the motor and K(θ) is a matrix of chemical transition rates. The chemical states, S i, are identified as ATP hydrolysis, phosphate release, power-stroke, and ADP (adenosine diphosphate) release along with binding and detachment steps of myosin. ATP and ADP release are dependent on the conformation of myosin, therefore the transition rates for these states depend on θ. Lan and Sun determined values of the chemical transition rates from experimental data. In this work, we apply the concept of coupled equations as described in equation (1.1) and (1.2) to our three dimensional, coarse-grained model of a molecular motor on a substrate. We combine Brownian dynamics simulations for the positional degrees of freedom with a numeric solution of the equations of motion for the chemical variables to simulate a protein taking steps on a substrate. The cargo is attached to the motor protein with a harmonic spring. 3
14 1.1 Outline of thesis This thesis is organized as follows: In chapter II we present the model for the protein, the cargo, and the substrate and we define the mechano-chemical states of our model. The simulation method is described in chapter III and results are presented in chapter IV. We give a brief summary and conclusions in chapter V. 4
15 CHAPTER II MODEL We develop a coarse grained model of a motor protein walking on a bead spring substrate under tension. The walker has a total of five sites, with two binding sites in each of two head domains and a central site. The two head sites and the central sites are bonded to each other with harmonic springs. Equilibrium conformation and motion of the motor protein and a substrate are determined by interactions between the constituents. The surrounding medium is treated as an implicit solvent and provides a random force that we discuss in chapter 3. Figure 2.1: Coarse grained model of a motor protein on a bead spring substrate. Sites 1, 2, 3 and 4 represents protein heads and 5 represents the central site, which is bonded to the light chain domain for the attachment of cargo structures. 5
16 We use Lennard Jones potential parameters as units for the energy and length in our model. To complete the units, we set the mass m of a bead of the chain to one. With these choices, the units for temperature, time, etc are also fixed. Table 2.1 shows model units and SI units. Table 2.1: Model units of physical quantities Parameters Model SI Length [L] σ m Energy [E] ɛ J Temperature [T] ɛ k B K Mass [M] m kg Time [t] mσ 2 ɛ s Viscosity [η] m tσ Pa.s To make the connection to biological systems, we assume a temperature of 300 K, which implies ɛ = J. If we set σ = 5 nm, about times the size of a kinesin head domain, and set the density to that of water, ρ = 1000 kg m 3, we find m = ρπ σ3 6 = kg, about one third of the mass of a kinesin head [5, 6]. This gives s for the model unit of time. 6
17 2.1 Conservative Forces Interactions between the sites of the model are described by conservative forces. The total conservative force on particle i at position r i is given by F conservative ( r i ) = j { F hs ( r ij ) + F lj ( r ij ) + F m ( r ij )+ (2.1) F dih ( r ij ) + F bond ( r ij ) + F pull ( r ij )}, where F hs ( r ij ) represents harmonic spring forces, Flj ( r ij ) and F m ( r ij ) are due to Lennard-Jones and Morse potentials, respectively, and F dih ( r ij ) and F bond ( r ij ) are the forces related to the dihedral and bond angle, respectively. Fpull ( r ij ) represents an external tension applied in equal amounts but opposite directions to the end beads of the substrate chain to keep the substrate extended. The interaction forces are calculated from the negative gradients with respect to the particle s coordinates of the corresponding potential energy contributions. Let, U α be a potential energy contribution, then its corresponding force is calculated as, F α ( r i ) = ri U α. (2.2) The potentials describing interactions between the particles are now discussed in detail. Harmonic Spring: Harmonic springs represent bonds between two sites. With the compression or expansion of a bond, the bond length deviates from its equilibrium value. The potential 7
18 energy of the spring is the harmonic potential U hs (r ij ), [13] U hs (r ij ) = 1 2 K s(r ij B 0 ) 2, (2.3) where K s is the spring constant, r ij is the distance between the sites, and B 0 is the equilibrium bond length. Our model has three different sets of parameters to describe bonds between different types of particles. Section 2.2 explains which sites are connected by harmonic springs and Table 2.2 contains the parameter values. In figure 2.2 we represent graphs of the spring potentials used in this work. The force on particle i due to a harmonic spring is calculated as, F hs,i = K s (r ij B 0 )ˆr ij, (2.4) where ˆr ij = ( r j r i ) r ij is the unit vector pointing from particle i to particle j. Lennard-Jones Potential: The Lenard-Jones (LJ) potential is generally used to describe interactions between non-bonded pairs of neutral particles, [14] [ ( ) 12 ( ) ] 6 σ σ U lj (r ij ) = 4ɛ, (2.5) r ij where ɛ is the potential depth, r ij is the distance between the sites and σ sets the length scale. A graph of the LJ potential is included in Fig 2.3. The force on particle i due to the LJ potential is given by r ij [ ( F lj,i = 48ɛ 1 ) 12 σ 1 ( ) ] 6 σ ˆr ij. (2.6) r ij r ij 2 r ij 8
19 Figure 2.2: Graphs of spring potentials as a function of distance for three sets of parameters that describe bonds between different types of particles: BB for beadbead interactions along the substrate chain, HH for bonds between protein sites in the same head domain, HC for bonds between protein head sites and the central protein sites; values for the parameters are presented in Table
20 The LJ potential is a one length-scale potential with σ, the distance at which the potential value is zero, as its adjustable length. The attractive part of the potential is given by the r 12 term and the repulsive part is given by the r 6 term [14]. The cut-off value for the Lennard-Jones potential in our model is 3σ and the potential is shifted to make sure that its value is zero at the cut-off distance. Soft repulsive Lennard-Jones Potential: The soft repulsive Lennard-Jones (SLJ) potential is a purely repulsive Lennard-Jones potential, (U lj + ɛ) for r lj R c, U lj(r ij ) = 0 for r lj > R c, (2.7) where R c = 2 1/6 σ is the cut-off distance. A graph of the SLJ potential is included in Fig Morse Potential: A Morse potential describes interactions from the bond equilibrium to the detachment state, [13] U m (r ij ) = D e [ (1 e κ(r ij B m) ) 2 1 ], (2.8) where D e is the depth of the potential energy minimum, κ is a positive number which has the units of inverse of length and is used to control the decay length of the exponential term, r ij is the distance between the two sites and B m is the equilibrium bond length. The negative differentiation of the Morse potential with respect to distance yields the Morse force on particle i, 10
21 F m,i = 2κD e ( e κ(r ij B m) ) [ 1 e κ(r ij B m) ] ˆr ij. (2.9) The Morse potential has two length scales, B m and κ, which control the distance where the force becomes repulsive and the width of the attractive well, respectively. A graph of a Morse potential is included in Fig. 2.3; potential parameters are presented in Table 2.2 Figure 2.3: Graphs of four interparticle potentials as a function of distance between the particles. The solid lines represent the bead-bead harmonic spring (magenta), soft repulsive Lennard-Jones (black), Lennard-Jones (blue), and Morse (red) potentials, respectively. The parameters of the potentials are shown in Table 2.2. Dihedral Potential: A dihedral potential models the change in energy due to a change in the dihedral angle α, which is caused by a rotation of the dihedral planes [13, 14]. We use a generalized dihedral potential to control the conformation of the walker on the bead 11
22 spring substrate. Fig 2.4 shows the dihedral planes of our walker model. The dihedral potential is given by, U dih (r ij ) = V 0 (1 cos 2 α), (2.10) where V 0 is the dihedral amplitude and α is the angle between the dihedral planes. Fig 2.5 shows a graph of the dihedral potential; values for the parameters are included in Table 2.2. A dihedral force acts on the protein head and the central sites and is calculated from F dih,i = V 0 sin 2α j α r ij ˆr ij. (2.11) Figure 2.4: Illustration of the generalized dihedral planes of our motor protein model. The dihedral angle α is the angle between the normal unit vectors of the planes. Bond Angle Potential: Just as the harmonic spring potential, the bond angle potential is based on Hooke s law but the deviation from the equilibrium configuration is due to a change in the bond angle θ rather than a distance [13]. Fig 2.6 shows the angle θ in our walker 12
23 Figure 2.5: Graph of the dihedral potential as a function of the angle α for three different dihedral amplitudes V 0. model that we control with the bond angle potential, U bond (θ) = 1 2 K θ(θ θ 0 ) 2, (2.12) where K θ is the force constant, θ is the bond angle, and θ 0 is the equilibrium bond angle. Fig 2.7 shows a graph of the potential. The force is calculated from F bond,i = K θ (θ θ 0 ) j θ r ij ˆr ij. (2.13) 13
24 Figure 2.6: Illustration of the bond angle θ formed by two unit vectors â and ˆb along the bonds between the central protein site and head sites 2 and 3. Figure 2.7: Graph of the bond angle potential as a function of the bond angle θ. The equilibrium bond angle is θ 0 = in our model. 14
25 2.2 Mechanical State The mechanical state of the system is described by the positions of all the particles in the model and controlled by the potentials introduced in Section 2.1. Here we describe how we apply the potentials to define the interactions between the constituents of our model Intra protein interactions The sites of the protein interact with both bonded and non-bonded interactions as shown in Fig 2.8. Bonded interactions are represented by harmonic springs, which connect the two sites of each head and also the head sites to the central site. Non-bonded interactions between sites on different heads are described by the soft-repulsive Lennard-Jones potential. The bond angle interaction between two nonbonded head sites is used to control the distance between the two heads. The dihedral interaction controls the orientation of the planes. Figure 2.8 illustrates all the intra protein interactions Protein-substrate interactions Interactions between protein sites and substrate sites are described by two kinds of potentials, Morse and soft-repulsive Lennard-Jones. The Morse interaction between head sites of the protein and chain beads sets the attraction of the head to the substrate and allows us to control the release and attachment of the head to the substrate. The soft repulsive Lennard-Jones potential between the central site and the 15
26 Figure 2.8: Illustration of the interactions between the components of the protein. chain beads prevents overlap of the central site with the chain. Figure 2.9 illustrates the protein substrate interactions Substrate-substrate interactions Bonded interactions between the chain beads are represented by harmonic springs, while non-bonded interactions are described by a Lennard-Jones potential. An external tensile force, given by F pull, is introduced to keep the chain elongated for the protein walk. Figure 2.10 illustrates the substrate-substrate interactions. After defining the form of the interaction forces a parametrization process is required to complete the model [13]. We determined the interaction parameters of our model by evaluating trajectories and energy values for different trial parameters. The potential parameters for our model are presented in table
27 Figure 2.9: Illustration of the interactions between protein and substrate sites. Solid red arrows represent the Morse interaction between protein heads and the substrate, whereas dotted black arrows represent the soft repulsive Lennard-Jones interaction between the central protein and the substrate Figure 2.10: Illustration of the interactions between substrate sites. The heavy lines connecting the beads represent harmonic springs; the curved arrows indicate Lennard-Jonnes interactions between non bonded beads. An external tension force, represented by F pull, is applied to the end beads of the chain with equal magnitude and opposite directions. 17
28 Table 2.2: Interaction parameters. Protein Bonded Head-Head Harmonic Spring K s = 80ɛ/σ 2, B 0 = 0.9 σ Head-Central Harmonic Spring K s = 80ɛ/σ 2, B 0 = 1.8σ σ = 1, ɛ = 1, Non-Bonded Head-Central Soft Repulsive LJ R c = σ σ = 1, ɛ = 1, Head-Head Soft Repulsive LJ R c = σ Dihedral Planes Dihedral V 0 = 100 Head-Head Bond Angle θ 0 = Chain-Chain Bonded Harmonic Spring K s = 320ɛ/σ 2, B 0 = σ Non-Bonded Lennard Jones σ = 1, ɛ = 1, R c = 3σ External Tension F pull = 40ɛ/σ D e = 4ɛ, κ = 4/σ, Protein-Chain Chain-Head Morse B 0 = 0.9σ, B m = B 0 + ln(2)/κ σ = 1, ɛ = 1, Chain-Central Soft Repulsive LJ R c = σ 18
29 To simulate the transport of a cargo, like a vesicle, we add one more site to the model. This cargo site is bonded to the central site of the motor protein with a harmonic spring and interacts with non-bonded sites through the soft repulsive Lennard-Jones potential. Fig 2.1 includes an illustration of the motor protein with cargo. The cargo particle has four times the size and 64 times the mass of a protein or chain bead site. The parameters used to define the cargo interactions are given in table 2.3. Table 2.3: Interaction parameters for cargo. Cargo Head and Chain Soft Repulsive LJ Size = 4σ, σ c = 1 ( Size +σ), ɛ = 1, 2 R c = (σ c + σ) Central Harmonic Spring K s = 80ɛ/σ 2, B 0 = 4σ 19
30 2.3 Chemical Steps In the biological system, the chemical steps regulating the dynamics of a motor protein include ATP hydrolysis, phosphate release, powerstroke and ADP release [11]. In the Introduction, we discussed how Lan and Sun [11] mapped the biological system onto a one-dimensional model for a walker to obtain a mathematical description of the walking process for myosin. They used a model with eight states and experimental values for the transition rates. For our three dimensional model of a walker, we define n s = 4 chemical states. They represent the equilibrium protein conformations that serve as intermediate steps along the walk. The states are associated with different interaction parameter values during the simulations. Mathematically, the chemical state of the system is represented by a vector, S, with n s components. The change of the chemical state of the motor protein is defined by the first order differential equation S t = K S (2.14) where, K is an n s n s matrix of transition rates. States 1 and 4 in our model are identical and represent a conformation where both heads are in contact with the chain. To simulate a motor that can make a step with either head one or head two, 20
31 we define two transition matrices, K and K1. k fast k fast k slow 0 0 K =, (2.15) 0 k slow k slow k slow k slow 0 0 k slow 0 k fast K1 =, (2.16) 0 k slow k slow k fast where k fast and k slow represent fast and slow processes, respectively. They have the values k fast = and k slow = in our model. According to equation (2.14), the set of coupled differential equations for chemical state n is given by S n t = j K nj S j, (2.17) where S n is the probability to find the system in state n so that S n = 1. (2.18) n The chemical state is coupled to the mechanical state through the values of the potential parameters, Q( S) = n Q n S n, (2.19) where Q is any of the parameters (D e, K θ, V 0 ) that change during the simulation. This is further discussed in chapter 4. The system of equations (2.17), (2.18), (2.19) 21
32 is used to advance the chemical state during the simulations. Figure 2.11 shows the evolution of the states. Figure 2.11: Chemical state as a function of time for a two-step simulation. The initial state, S 1 = 1, S 2 = S 3 = S 4 = 0 is maintained during an equilibration period at the start of the simulation. The probability S 1 to be in state 1 (red line) falls off rapidly as S 2 (blue line) first increases. As S 2 decreases S 3 (purple line) first increases and then decreases as S 4 starts to grow. At the end of the first step, S 1 =S2 =S2 = 0 and S 1 = 1. In the second step the system goes from S4 = 1 to S 1 through the states S 2 and S 3. 22
33 CHAPTER III METHOD 3.1 Brownian Dynamics Simulation Since motor proteins are small particles in a liquid environment they undergo Brownian motion when they are not actively changing conformation. Brownian motion is the result of collisions with molecules in the surrounding solvent [15]. Since our motor protein model has an implicit solvent we use Brownian dynamics simulations to describe the motion. The net force on particle i of our model is given by F net,i = F conservative,i + F frictional,i + F random,i. (3.1) where F conservative,i is the conservative force defined in equation (2.1), F frictional,i is the frictional force due to the viscous medium, F frictional,i = ζ i v i. (3.2) where ζ i is the friction coefficient and v i is the velocity of the i th particle. Frandom,i is the random force that describes Brownian motion. From Newton s second law, F net,i = m i a i, (3.3) where m i is the mass of the i th particle and a i its acceleration, we find, m i a i = F conservative,i + F frictional,i + F random,i. (3.4) 23
34 As the inertial term is very small [1], we ignore the left-hand-side of equation (3.4) and, substituting F frictional,i from equation (3.2), we obtain ζ i v i = F conservative,i + F random,i (3.5) v i = v D i + v B i, (3.6) where v D i (t) is the deterministic velocity given by, v D i (t) = 1 ζ i Fconservative,i (t), (3.7) and v B i (t) is the random velocity v B i (t) = 1 ζ i Frandom,i. (3.8) The random force has an average value of zero, F random,i = 0. Its correlations are related to the friction coefficient F random,iα (t)f random,iβ (t ) = γ i δ(t t ), (3.9) where γ i = 2k B T ζ i and δ is the Dirac delta function. In thermal equilibrium, the theorem of equipartition of energy gives, K.E. = 3 2 Nk BT, (3.10) where K.E. is the average kinetic energy of N particles in three dimensions. In a Brownian dynamics simulation with time step δt, the random velocity is calculated from, v α B (t) = 2T δtζ n rα, (3.11) 24
35 where T is the reduced temperature and n r is a random number in Gaussian distribution, n rα = 0, n rα n rβ = δ αβ. The friction coefficient of a spherical particle ζ is given by, ζ = 6πηa, (3.12) where η is the shear viscosity and a is the radius of the particle. Since the random velocity determines the temperature, equations (3.11) and (3.12) help us to find a useful time step for Brownian dynamics simulations, δt = 2m 6πηa, (3.13) where m is the mass of the particle [15]. In our simulations, we use m = 1, a = σ/2 and δt = with η = to satisfy equation (3.13). Equations (3.7) and (3.11) with (3.12) are valid when all particles have the same mass. To simulate the cargo, which has larger mass, we set ζ c = ( mc ) ( ) σ m σ c 3πησ c, (3.14) where m c and σ c are the mass and diameter of the cargo, while m and σ are the mass and diameter of all the other particles. During the simulation, we collect data on quantities such as the kinetic energy, the potential energy and distances between particles. To obtain error estimates, we calculate the average of a quantity ξ i in a block of simulation steps and calculate the average according to equation ξ = 1 N b 25 N b i=1 ξ i, (3.15)
36 where N b is the number of blocks. The standard deviation of the block values is an estimate for the statistical uncertainty and calculated from σ(ξ) = 1 N b ξ i ξ N b 1 2. (3.16) i=1 3.2 Behavior of an isolated protein. To learn to control its conformations, we studied the behavior of the isolated protein for different values of the potential parameters K θ and V 0. Figure 3.1 shows protein conformations for relevant K θ and V 0 values. The figures illustrate that the distance Figure 3.1: Effect of K θ and V 0 on conformations of an isolated protein. between the heads decreases when the value of K θ is changed from 0 to 20. For different V 0 values the dihedral planes have different orientations relative to each 26
37 other. When V 0 = 100, the planes are parallel to each other, for V 0 = 100, on the other hand, the planes are perpendicular in orientation. 3.3 Conformation of the protein in contact with the substrate. We investigated the effect of the model parameters on the protein and chain system to identify the intermediate states for a protein taking a step. Figure 3.2 shows conformations of the chain and walker for the three identified states. The parameters Figure 3.2: Illustration of the mechano-chemical states (1, 2, and 3) of a protein walker on a substrate. The potential parameters of the states are presented in Table
38 of each state described in figure 3.2 are listed in table 3.1. State 1 is the initial state of the protein walk, where both heads are attached to the substrate. The dihedral planes are parallel in orientation and the bond-angle potential is zero. State 2 is the second state in a step; the Morse potential for one head is turned off, the dihedral planes are at right angles, and K θ is zero to help the protein head move round. In state 3 both heads are attached to the substrate and the bond angle interaction is turned on to achieve a small distance between head sites. We used these states to construct a 4-step model for the walk, where state 4 is the same as state 1. During the walk, when head one completes one step, the fourth state becomes the initial state for head 2 for the next step and so on. For head two to take a step the values of D e1, D e2, D e3 and D e4 in state 2 becomes 4, 4, 0 and 0, respectively. Table 3.1: Parameters for the equilibrium states Parameters D e1 D e2 D e3 D e4 K θ V 0 State State State State In our simulations the probabilities S 1, S 2, S 3 and S 4 to be in the states evolve in time as shown in figure Since S 2 and S 3 never come close to a value of one, we increase the range of the potential parameters when we calculate their 28
39 time-dependent values with equation (2.17). In table 3.2, we express the potential Table 3.2: Assigned parameters for the change of states to maintain the desired amplitudes. Parameters D e1 D e2 D e3 D e4 K θ V 0 State 1 State 2 State 3 State 4 D e,max D e,max D e,max D e,max K θ,min V pos D e,min D e,min D e,max D e,max K θ,min V neg D e,max D e,max D e,max D e,max K θ,max V pos D e,max D e,max D e,max D e,max K θ,min V pos parameters in terms of minimum and maximum values that, when used in equation (2.17), reproduce the values in Table 3.1. For the state evolution shown in figure 2.11, we use the values D e,max = 4ɛ, D e,min = -0.95ɛ, K θ,max = 53.5, K θ,min = 0, V pos = 100ɛ and V neg = 148ɛ. 3.4 Kinetic energy The instantaneous kinetic energy of the system is calculated from K.E. = 1 2 N m i v 2 i, (3.17) i=1 where N is the total number of particles and m i and v i are their masses and speeds, respectively. The instantaneous kinetic energy fluctuates rapidly due to the random forces. To examine the kinetic energy over larger time scales, we average its values 29
40 over segments of fixed length. In figure 3.3 we present graphs of the instantaneous and the segment-averaged kinetic energy for two simulations. In both cases, the system consists of the chain, the protein, and the cargo with a total number of N = 26 particles. For the reduced temperature k B T = 1.0 used in the simulations, we expect an average kinetic energy of K.E. = 39 according to the principle of equipartition of energy, Eq (3.10). This is exactly what we observe. Panel A of figure 3.3 shows results for fixed interaction parameters, whereas panel B shows results for changing parameters. During the time period in the figure, head one was detached from the chain substrate. This period was chosen to see if any variation in the segment-averaged kinetic energy occurs due to the change in the parameters. This is not the case, which shows that the system remains in thermal equilibrium. We expected this result because, in our simulations, the kinetic transition rates are very slow compared to the time scale of thermal fluctuations. 30
41 Figure 3.3: Instantaneous (red) and segment-averaged (blue) kinetic energy as a function of time for simulations of the protein-vesicle-chain system, which has a total of N = 26 particles. Panel A and B show results for fixed and time-varying interaction parameters, respectively. The length of one time segment in the averaging process is steps, corresponding to a simulation time of 500. In both simulations, the average kinetic energy has a value of 39, in agreement with the equipartition theorem and a reduced temperature of one. 31
42 CHAPTER IV RESULTS 4.1 Force-extension relation of the substrate The extension of the substrate chain is the distance between the first and last chain bead along the direction of the applied force. Since we apply the external force in the x-direction the extension X is given by X = (X Nc X 1 ), (4.1) where N c is the number of chain beads, N c = 20, and X Nc and X 1 are the x-coordinates of the last and first bead, respectively. In figure 4.1 we present results of the extension as a function of the pulling force. Due to the strong attractive interactions between the chain beads, the chain is found to be in a collapsed state in the absence of external tension. With the increase in the pulling force the extension of the substrate increases. At a particular pulling force, there is a rapid change in the extension as the chain is stretched from its globule state. Then the extension increases more slowly with the force. In our simulations with the protein walker we apply a tension force of 40 so that the chain is highly stretched. 32
43 Figure 4.1: Force extension results for an isolated chain. The symbols represent simulation results for the extension, defined in equation (4.1), as a function of the pulling force for a chain of length N c = 20. The error bars indicate the standard deviation of block averages (see equation (3.16)). Two chain conformations for intermediate (left) and high (right) extensions are also shown. 33
44 4.2 Motor protein taking steps on the substrate without cargo. We have performed simulations for a motor protein without cargo for steps corresponding to a simulation time of We also added an equilibration time of 10 7 steps so that the protein comes to thermal equilibrium before the chemical variables start to change. To investigate if the change in parameter values has the desired effect, we calculate the distance d 1 between site 1 on head one and the closest bead on the chain and the corresponding distance d 3 for site 3 on head two. In figure 4.2 we present simulation results for d 1 and d 3 in panel A together with the Morse potential amplitudes D e1 and D e3 in panel B. At time the value of the Morse potential amplitude D e1 starts decreasing to zero. Almost immediately, the distance d 1 between head one and the chain increases. It reaches a maximum value of about 1.8 before it decreases again to about one when D e1 goes back to its original, attractive value. The distance d 3 is nearly constant during this step. C 1 and C 3 represent the closest chain beads for head sites 1 and 3 respectively. During the time when head site 1 is released the bead index C 1 increases, showing that head one moves forward along the chain. At the same time, head two slides backward along the chain. When head one is attached again, both C 1 and C 3 are slowly sliding backward. At time the Morse amplitude D e3 starts decreasing and head two is released and reaches a maximum distance of about 1.6. The bead index C 3 at that time increases which suggests that the protein is taking a step forward. But at the same time the increase of the bead index C 1 suggest that the protein glides along the 34
45 chain. Since the protein moves along the chain even when both heads are attached, we decided to adjust the interaction parameters. 35
46 Figure 4.2: Two steps of a motor protein walker on a chain substrate. Panel A shows the smallest distance between the chain and head site 1 on head one (d 1, red) and the chain and head site 3 on head two (d 3, pink) together with the indices of the beads closest to the head sites (C 1, blue and C 3, cyan). Panel B shows the value of the Morse potential parameters D e1 and D e3 as a function of time. 36
47 4.3 Effect of changing the interaction of the attached head of a motor protein walking without cargo. In figure 4.2 we saw that both heads are moving when one head is released. keep the non-released head fixed, we change the Morse amplitudes D e3 and D e4 To to 2 D e,max when head one is released and we change D e1 and D e2 to 2 D e,max when head two is released. Figure 4.3 B shows the Morse amplitudes as a function of time while results of a simulation with these parameters are shown in figure 4.3 A. The results for C 1 and C 3 in panel A show that immediately after the attachment of head one and head two, the protein glides backward. 37
48 Figure 4.3: Two steps of a motor protein on a substrate with increased attraction for the attached head site. Panel A shows the smallest distance between the chain and head site 1 on head one (d 1, red) and the chain and head site 3 on head two (d 3, pink) together with bead indices closest to head site 1 (C 1, blue) and head site 3 (C 3, cyan). Panel B shows the value of the Morse potential parameters D e1 and D e3 as a function of time. 38
49 4.4 Motion of protein with cargo. We also studied a walking motor with attached cargo. In figure 4.4 we show results for the distances to the chain and nearest bead indices for head sites and cargo. The Morse amplitudes for this case are the same as in panel B of figure 4.2. The results in figure 4.4 show that when head one reattaches the protein stays in place for a while and then slides forward. After the attachment of head two, the protein starts moving backward along the chain even though both heads are attached to the chain. Figure 4.4: Two steps of a motor protein walking with cargo on a substrate. The graph shows the smallest distance between the chain and head site 1 (d 1, red), head site 3 (d 3, pink) and the cargo (ves, orange). The change in Morse parameters is as shown in panel B of graph 4.2. The closest beads are C 1 (blue), C 3 (cyan) and C ves (yellow) for head sites 1, 3, and the cargo, respectively. In these simulations, the Morse parameters of the non-released head remained constant, as shown in panel B of figure
50 In figure 4.5 we present simulation results for the case that the non-released head is more strongly attached as shown in figure 4.3 B. To find out how reproducible the results are, we performed simulations with identical parameters and initial configurations as in figure 4.5 but a different seed of the random number generator. The results shown in figure 4.6 have qualitatively the same behavior as those in figure 4.5. Figure 4.5: Two steps of a motor protein walking with cargo. The lines represent the same quantities as in figure 4.4; the Morse amplitudes for this simulation are shown in figure 4.3 B. 40
51 Figure 4.6: Results for the same system as in figure 4.5 from a separate simulation. 41
52 4.5 Five steps for a protein with cargo. Considering the results presented in figures , we concluded that two issues prevented the protein from walking effectively. For one, the protein sometimes moves the wrong foot, for example when the heads switched positions during the wait period. For two, the protein tends to slide while waiting for the next step to start. To address the first issue, we consider the position of the heads before a step is taken: the head farthest from the target chain end is determined and that head is released in the next step. To address the second issue, the wait time before the next step is reduced. To illustrate the protocol, we show the Morse amplitudes as a function of time for a five step walk in figure 4.7 B. The results for closest chain beads in panel A show that the motor protein is moving with its cargo in the forward direction. For the first three steps the movement of the motor protein is as intended. After, the third step, the protein is close to the end of the chain and the position of the head sites and the cargo fluctuate about a final position. 42
53 Figure 4.7: 5-step simulation with the new protocol. The lines represent the same quantities as in figure
54 CHAPTER V CONCLUSIONS In this work we have developed a five-site model for a motor protein in contact with a bead-spring substrate and performed Brownian dynamics simulations of the model. We have investigated the effect of model parameters on protein and chain conformations and used this information to determine suitable interaction parameters for our model. To describe a motor taking steps on the substrate, we have identified equilibrium states that serve as chemical states for the walking protein and evolve in time according to coupled rate equations. Since motor proteins are responsible for carrying cargo from one part of the cell to another, we have added a cargo structure in the form of a large sphere that is bonded by a harmonic spring to the central site of the motor protein. To investigate mechanical properties of the substrate and validate our code, we performed a series of simulations for an isolated chain under tension. A force extension curve for the isolated chain shows the expected behavior for stretching a chain from a globule to an elongated state. We have investigated the motor protein taking steps and found that we were able to control the attachment of the protein to the substrate. To investigate the effect of cargo on the protein s walk we have performed simulations with cargo and found that a large sphere is dragged very 44
55 slowly. Since thermal fluctuations play a large role in the system, we tested the reproducibility of our results by repeating simulations with identical conditions but different seeds for the random number generator. We found that the results produced agree qualitatively. Since the protein has a tendency to step back and forth we developed a protocol that allows the walker to keep going in the same direction on average. Our results show that the model and simulations developed in this work can be used as a starting point for a more detailed investigation of cargo-carrying motor protein walkers. 45
56 BIBLIOGRAPHY [1] Debashish Chowdhury. Stochastic mechano-chemical kinetics of molecular motors: A multidisciplinary enterprise from a physicist s perspective. Physics Reports, 529:1 197, [2] Debashish Chowdhury. Resource letter PBM-1: Physics of biomolecular machines. American Journal of Physics, 77: , [3] Philip Nelson. Biological Physics: Energy, Information, Life. W. H. Freeman, New York, [4] Anthony Brown. Axonal transport. In Donald W. Pfaff, editor, Neuroscience in the 21st Century From Basic to Clinical, chapter 10, pages Springer, New York, [5] H. Lodish, A. Berk, S. L. Zipursky, P. Matsudaira, D. Baltimore, and J. Darnell. Molecular Cell Biology. W. H. Freeman, New York, 4th edition, [6] R. Mallik and S. P. Gross. Molecular motors: Strategies to get along. Current Biology, 14:R971 R982, [7] D. Keller and C. Bustamante. The mechanochemistry of molecular motors. Biophys. J., 78: , [8] I. Derényi and T. Vicsek. The kinesin walk: A dynamic model with elastically coupled heads. Proc. Natl. Acad. Sci. USA, 93: , [9] C. Bustamante, D. Keller, and G. Oster. The physics of molecular motors. Acc. Chem. Res, 34: , [10] A. B. Kolomeisky and H. Phillips III. Dynamic properties of motor proteins with two subunits. J. Phys.: Condens. Matter, 17:S3887 S3899, [11] G. Lan and X. Sun. Dynamics of myosin-driven skeletal muscle contraction: I. Steady-state force generation. Biophys. J., 88: ,
For slowly varying probabilities, the continuum form of these equations is. = (r + d)p T (x) (u + l)p D (x) ar x p T(x, t) + a2 r
3.2 Molecular Motors A variety of cellular processes requiring mechanical work, such as movement, transport and packaging material, are performed with the aid of protein motors. These molecules consume
More informationAnatoly B. Kolomeisky. Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS?
Anatoly B. Kolomeisky Department of Chemistry CAN WE UNDERSTAND THE COMPLEX DYNAMICS OF MOTOR PROTEINS USING SIMPLE STOCHASTIC MODELS? Motor Proteins Enzymes that convert the chemical energy into mechanical
More informationChapter 16. Cellular Movement: Motility and Contractility. Lectures by Kathleen Fitzpatrick Simon Fraser University Pearson Education, Inc.
Chapter 16 Cellular Movement: Motility and Contractility Lectures by Kathleen Fitzpatrick Simon Fraser University Two eukaryotic motility systems 1. Interactions between motor proteins and microtubules
More informationTransport of single molecules along the periodic parallel lattices with coupling
THE JOURNAL OF CHEMICAL PHYSICS 124 204901 2006 Transport of single molecules along the periodic parallel lattices with coupling Evgeny B. Stukalin The James Franck Institute The University of Chicago
More informationThe biological motors
Motor proteins The definition of motor proteins Miklós Nyitrai, November 30, 2016 Molecular machines key to understand biological processes machines in the micro/nano-world (unidirectional steps): nm,
More informationBME Engineering Molecular Cell Biology. Basics of the Diffusion Theory. The Cytoskeleton (I)
BME 42-620 Engineering Molecular Cell Biology Lecture 07: Basics of the Diffusion Theory The Cytoskeleton (I) BME42-620 Lecture 07, September 20, 2011 1 Outline Diffusion: microscopic theory Diffusion:
More informationRegularity and synchrony in motor proteins
Regularity and synchrony in motor proteins R. E. Lee DeVille and Eric Vanden-Eijnden Courant Institute of Mathematical Sciences, New York University, New York, NY 10012 1 Corresponding author: R. E. Lee
More informationSupplementary Information
Supplementary Information Switching of myosin-v motion between the lever-arm swing and Brownian search-and-catch Keisuke Fujita 1*, Mitsuhiro Iwaki 2,3*, Atsuko H. Iwane 1, Lorenzo Marcucci 1 & Toshio
More informationHSC PHYSICS ONLINE B F BA. repulsion between two negatively charged objects. attraction between a negative charge and a positive charge
HSC PHYSICS ONLINE DYNAMICS TYPES O ORCES Electrostatic force (force mediated by a field - long range: action at a distance) the attractive or repulsion between two stationary charged objects. AB A B BA
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature09450 Supplementary Table 1 Summary of kinetic parameters. Kinetic parameters were V = V / 1 K / ATP and obtained using the relationships max ( + m [ ]) V d s /( 1/ k [ ATP] + 1 k ) =,
More informationAccording to the diagram, which of the following is NOT true?
Instructions: Review Chapter 44 on muscular-skeletal systems and locomotion, and then complete the following Blackboard activity. This activity will introduce topics that will be covered in the next few
More informationIntracellular transport
Transport in cells Intracellular transport introduction: transport in cells, molecular players etc. cooperation of motors, forces good and bad transport regulation, traffic issues, Stefan Klumpp image
More informationNeurite formation & neuronal polarization
Neurite formation & neuronal polarization Paul Letourneau letou001@umn.edu Chapter 16; The Cytoskeleton; Molecular Biology of the Cell, Alberts et al. 1 An immature neuron in cell culture first sprouts
More informationChapter 1. Introduction
Chapter 1. Introduction 1a) This is an e-book about the constructive effects of thermal energy in biology at the subcellular level. Thermal energy manifests itself as thermal fluctuations, sometimes referred
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationBiophysik der Moleküle!
Biophysik der Moleküle!!"#$%&'()*+,-$./0()'$12$34!4! Molecular Motors:! - linear motors" 6. Dec. 2010! Muscle Motors and Cargo Transporting Motors! There are striking structural similarities but functional
More informationSupplementary Information. SI Text 1: Derivation and assumptions of the effective temperature model
Supplementary Information SI Text 1: Derivation and assumptions of the effective temperature model We assume that the displacements of intracellular particles are due to passive thermal activity and active
More informationLecture 7 : Molecular Motors. Dr Eileen Nugent
Lecture 7 : Molecular Motors Dr Eileen Nugent Molecular Motors Energy Sources: Protonmotive Force, ATP Single Molecule Biophysical Techniques : Optical Tweezers, Atomic Force Microscopy, Single Molecule
More informationLinear Motors. Nanostrukturphysik II, Manuel Bastuck
Molecular l Motors I: Linear Motors Nanostrukturphysik II, Manuel Bastuck Why can he run so fast? Usain Bolt 100 m / 9,58 s Usain Bolt: http://www.wallpaperdev.com/stock/fantesty-usain-bolt.jpg muscle:
More informationLecture 1: Introduction. Keywords: Mathematics as a language, Need of learning mathematics, Applications of mathematics in BIology
NPTEL Syllabus Biomathematics - Video course COURSE OUTLINE Graphs and functions, Derivative of a function, Techniques of differentiation Differentiation and its application in Biology, Finding maxima,
More informationNeurite formation & neuronal polarization. The cytoskeletal components of neurons have characteristic distributions and associations
Mechanisms of neuronal migration & Neurite formation & neuronal polarization Paul Letourneau letou001@umn.edu Chapter 16; The Cytoskeleton; Molecular Biology of the Cell, Alberts et al. 1 The cytoskeletal
More informationBiophysics Biological soft matter
Biophysics Biological soft matter!"#$%&'(&)%*+,-.& /"#$%("%*+,-.0."122,13$(%4(5+& Biophysics lectures outline Biological soft matter 1. Biopolymers 2. Molecular motors 3. The cytoskeleton Biophysics 1.
More informationPhysics of Cellular materials: Filaments
Physics of Cellular materials: Filaments Tom Chou Dept. of Biomathematics, UCLA, Los Angeles, CA 995-766 (Dated: December 6, ) The basic filamentary structures in a cell are reviewed. Their basic structures
More informationChapter 13. Simple Harmonic Motion
Chapter 13 Simple Harmonic Motion Hooke s Law F s = - k x F s is the spring force k is the spring constant It is a measure of the stiffness of the spring A large k indicates a stiff spring and a small
More informationMetabolism: Energy and Enzymes. February 24 th, 2012
Metabolism: Energy and Enzymes February 24 th, 2012 1 Outline Forms of Energy Laws of Thermodynamics Metabolic Reactions ATP Metabolic Pathways Energy of Activation Enzymes Photosynthesis Cellular Respiration
More informationPHYSICS 220 LAB #5: WORK AND ENERGY
Lab Section / 33 pts Name: Partners: PHYSICS 0 LAB #5: WORK AND ENERGY OBJECTIVES 1. To get practice calculating work.. To understand the concept of kinetic energy and its relationship to the net work
More information(1) (3)
1. This question is about momentum, energy and power. (a) In his Principia Mathematica Newton expressed his third law of motion as to every action there is always opposed an equal reaction. State what
More informationMolecular Motors. Structural and Mechanistic Overview! Kimberly Nguyen - December 6, 2013! MOLECULAR MOTORS - KIMBERLY NGUYEN
Molecular Motors Structural and Mechanistic Overview!! Kimberly Nguyen - December 6, 2013!! 1 Molecular Motors: A Structure and Mechanism Overview! Introduction! Molecular motors are fundamental agents
More informationPhysics A - PHY 2048C
Kinetic Mechanical Physics A - PHY 2048C and 11/01/2017 My Office Hours: Thursday 2:00-3:00 PM 212 Keen Building Warm-up Questions Kinetic Mechanical 1 How do you determine the direction of kinetic energy
More informationSection 5: Conservation of Energy and Gravitation
Section 5: Conservation of Energy and Gravitation 5.01 Work 5.02 Kinetic and Gravitational Potential Energy Physics (6)(B) Physics (6)(C) 5.03 Conservation of Energy Physics (6)(A) 5.04 Elastic Potential
More informationDynamics of an inchworm nano walker
Dynamics of an inchworm nano walker A. Ciudad a, J.M. Sancho a A.M. Lacasta b a Departament d Estructura i Constituents de la Matèria, Facultat de Física, Universitat de Barcelona, Diagonal 647, E-08028
More informationPhysics 20 Homework 3 SIMS 2016
Physics 20 Homework 3 SIMS 2016 Due: Thursday, August 25 th Special thanks to Sebastian Fischetti for problems 1, 5, and 6. Edits in red made by Keith Fratus. 1. The ballistic pendulum is a device used
More informationSupplemental Data. Stepping, Strain Gating, and an Unexpected. Force-Velocity Curve. for Multiple-Motor-Based Transport
Supplemental Data Stepping, Strain Gating, and an Unexpected Force-Velocity Curve for Multiple-Motor-Based Transport Ambarish Kunwar, Michael Vershinin, Jing Xu, and Steven P. Gross Supplemental Experimental
More informationActo-myosin: from muscles to single molecules. Justin Molloy MRC National Institute for Medical Research LONDON
Acto-myosin: from muscles to single molecules. Justin Molloy MRC National Institute for Medical Research LONDON Energy in Biological systems: 1 Photon = 400 pn.nm 1 ATP = 100 pn.nm 1 Ion moving across
More informationDiffusion in biological systems and Life at low Reynolds number
Diffusion in biological systems and Life at low Reynolds number Aleksandra Radenovic EPFL Ecole Polytechnique Federale de Lausanne Bioengineering Institute Laboratory of Nanoscale Biology Lausanne September
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationThe graph shows how an external force applied to an object of mass 2.0 kg varies with time. The object is initially at rest.
T2-2 [195 marks] 1. The graph shows how an external force applied to an object of mass 2.0 kg varies with time. The object is initially at rest. What is the speed of the object after 0.60 s? A. 7.0 ms
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationPhysics 218: Waves and Thermodynamics Fall 2002, James P. Sethna Homework 12, due Wednesday Dec. 4 Latest revision: November 26, 2002, 10:45 am
Physics 218: Waves and Thermodynamics Fall 2002, James P. Sethna Homework 12, due Wednesday Dec. 4 Latest revision: November 26, 2002, 10:45 am Reading Feynman, I.44 Laws of Thermodynamics, I.45 Illustrations
More informationLecture 18. Newton s Laws
Agenda: l Review for exam Lecture 18 l Assignment: For Monday, Read chapter 14 Physics 207: Lecture 18, Pg 1 Newton s Laws Three blocks are connected on the table as shown. The table has a coefficient
More informationPHYSICS. Chapter 9 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 9 Lecture RANDALL D. KNIGHT Chapter 9 Work and Kinetic Energy IN THIS CHAPTER, you will begin your study of how energy is transferred
More information3/10/2019. What Is a Force? What Is a Force? Tactics: Drawing Force Vectors
What Is a Force? A force acts on an object. A force requires an agent, something that acts on the object. If you throw a ball, your hand is the agent or cause of the force exerted on the ball. A force
More informationPotential energy and conservation of energy
Chapter 8 Potential energy and conservation of energy Copyright 8.1_2 Potential Energy and Work Potential energy U is energy that can be associated with the configuration (arrangement) of a system of objects
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More informationIntroduction: actin and myosin
Introduction: actin and myosin Actin Myosin Myosin V and actin 375 residues Found in all eukaryotes Polymeric Forms track for myosin Many other cellular functions 36 nm pseudo-helical repeat Catalytic
More informationDirected Assembly of Functionalized Nanoparticles with Amphiphilic Diblock Copolymers. Contents
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2017 Electronic Supplementary Information for Directed Assembly of Functionalized Nanoparticles
More informationJoy of Science Experience the evolution of the Universe, Earth and Life
Joy of Science Experience the evolution of the Universe, Earth and Life Review Introduction Main contents Group discussions Unless otherwise noted, pictures are taken from wikipedia.org Review 1 Organic
More informationENERGY. Conservative Forces Non-Conservative Forces Conservation of Mechanical Energy Power
ENERGY Conservative Forces Non-Conservative Forces Conservation of Mechanical Energy Power Conservative Forces A force is conservative if the work it does on an object moving between two points is independent
More informationWhat Is a Force? Slide Pearson Education, Inc.
What Is a Force? A force acts on an object. A force requires an agent, something that acts on the object. If you throw a ball, your hand is the agent or cause of the force exerted on the ball. A force
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationPSI AP Physics B Dynamics
PSI AP Physics B Dynamics Multiple-Choice questions 1. After firing a cannon ball, the cannon moves in the opposite direction from the ball. This an example of: A. Newton s First Law B. Newton s Second
More informationLectures Chapter 6 (Cutnell & Johnson, Physics 7 th edition)
PH 201-4A spring 2007 Work and Energy Lectures 16-17 Chapter 6 (Cutnell & Johnson, Physics 7 th edition) 1 Work and Energy: Work done by a constant force Constant pushing force F pointing in the same direction
More informationFREEMAN MEDIA INTEGRATION GUIDE Chapter 7: Inside the Cell
FREEMAN MEDIA INTEGRATION GUIDE Chapter 7: Inside the Cell All media is on the Instructors Resource CD/DVD JPEG Resources Figures, Photos, and Tables PowerPoint Resources Chapter Outline with Figures Lecture
More informationMuscle regulation and Actin Topics: Tropomyosin and Troponin, Actin Assembly, Actin-dependent Movement
1 Muscle regulation and Actin Topics: Tropomyosin and Troponin, Actin Assembly, Actin-dependent Movement In the last lecture, we saw that a repeating alternation between chemical (ATP hydrolysis) and vectorial
More informationPhysics General Physics. Lecture 3 Newtonian Mechanics. Fall 2016 Semester. Prof. Matthew Jones
Physics 22000 General Physics Lecture 3 Newtonian Mechanics Fall 2016 Semester Prof. Matthew Jones 1 Review of Lectures 1 and 2 In the previous lectures we learned how to describe some special types of
More informationPHYSICS. Chapter 5 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 5 Lecture RANDALL D. KNIGHT Chapter 5 Force and Motion IN THIS CHAPTER, you will learn about the connection between force and motion.
More informationPhysics 18 Spring 2010 Midterm 1
Physics 18 Spring 2010 Midterm 1 For the midterm, you may use one sheet of notes with whatever you want to put on it, front and back. Please sit every other seat, and please don t cheat! If something isn
More informationLecture Outline Chapter 6. Physics, 4 th Edition James S. Walker. Copyright 2010 Pearson Education, Inc.
Lecture Outline Chapter 6 Physics, 4 th Edition James S. Walker Chapter 6 Applications of Newton s Laws Units of Chapter 6 Frictional Forces Strings and Springs Translational Equilibrium Connected Objects
More informationWiley Plus. Final Assignment (5) Is Due Today: Before 11 pm!
Wiley Plus Final Assignment (5) Is Due Today: Before 11 pm! Final Exam Review December 9, 009 3 What about vector subtraction? Suppose you are given the vector relation A B C RULE: The resultant vector
More informationSection /07/2013. PHY131H1F University of Toronto Class 9 Preclass Video by Jason Harlow. Based on Knight 3 rd edition Ch. 5, pgs.
PHY131H1F University of Toronto Class 9 Preclass Video by Jason Harlow Based on Knight 3 rd edition Ch. 5, pgs. 116-133 Section 5.1 A force is a push or a pull What is a force? What is a force? A force
More informationChapter Four Holt Physics. Forces and the Laws of Motion
Chapter Four Holt Physics Forces and the Laws of Motion Physics Force and the study of dynamics 1.Forces - a. Force - a push or a pull. It can change the motion of an object; start or stop movement; and,
More informationTranslational Motion Rotational Motion Equations Sheet
PHYSICS 01 Translational Motion Rotational Motion Equations Sheet LINEAR ANGULAR Time t t Displacement x; (x = rθ) θ Velocity v = Δx/Δt; (v = rω) ω = Δθ/Δt Acceleration a = Δv/Δt; (a = rα) α = Δω/Δt (
More informationMolecular Motors. Dave Wee 24 Sept Mathematical & Theoretical Biology Seminar
Molecular Motors Dave Wee 24 Sept 2003 Mathematical & Theoretical Biology Seminar Overview Types of motors and their working mechanisms. Illustrate the importance of motors using the example of : ATP-Synthase
More informationSupplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers
Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department
More informationBell Ringer: What is Newton s 3 rd Law? Which force acts downward? Which force acts upward when two bodies are in contact?
Bell Ringer: What is Newton s 3 rd Law? Which force acts downward? Which force acts upward when two bodies are in contact? Does the moon attract the Earth with the same force that the Earth attracts the
More informationAuthor(s) Igarashi, A; Tsukamoto, S; Goko, H. Citation PHYSICAL REVIEW E (2001), 64(5)
Title Transport properties and efficiency Brownian motors Author(s) Igarashi, A; Tsukamoto, S; Goko, H Citation PHYSICAL REVIEW E (2001), 64(5) Issue Date 2001-11 URL http://hdl.handle.net/2433/50297 RightCopyright
More informationImportance of Hydrodynamic Interactions in the Stepping Kinetics of Kinesin
pubs.acs.org/jpcb Importance of Hydrodynamic Interactions in the Stepping Kinetics of Kinesin Yonathan Goldtzvik, Zhechun Zhang, and D. Thirumalai*,, Biophysics Program, Institute for Physical Science
More informationChapter 8: An Introduction to Metabolism. 1. Energy & Chemical Reactions 2. ATP 3. Enzymes & Metabolic Pathways
Chapter 8: An Introduction to Metabolism 1. Energy & Chemical Reactions 2. ATP 3. Enzymes & Metabolic Pathways 1. Energy & Chemical Reactions 2 Basic Forms of Energy Kinetic Energy (KE) energy in motion
More informationChapter 5 Force and Motion
Chapter 5 Force and Motion Chapter Goal: To establish a connection between force and motion. Slide 5-2 Chapter 5 Preview Slide 5-3 Chapter 5 Preview Slide 5-4 Chapter 5 Preview Slide 5-5 Chapter 5 Preview
More informationStudy Guide Solutions
Study Guide Solutions Table of Contents Chapter 1 A Physics Toolkit... 3 Vocabulary Review... 3 Section 1.1: Mathematics and Physics... 3 Section 1.2: Measurement... 3 Section 1.3: Graphing Data... 4 Chapter
More informationPre-Comp Review Questions- 8 th Grade
Pre-Comp Review Questions- 8 th Grade Section 1- Units 1. Fill in the missing SI and English Units Measurement SI Unit SI Symbol English Unit English Symbol Time second s. Temperature K Fahrenheit Length
More informationSMASIS A STOCHASTIC MECHANO-CHEMICAL MODEL FOR COOPERATIVE MOTOR PROTEIN DYNAMICS
Proceedings of SMASIS 8 8 ASME Conference on Smart Materials, Adaptive Structures and Intelligent Systems October 8-3, 8, Ellicott City, Maryland, USA SMASIS8-585 A STOCHASTIC MECHANO-CHEMICAL MODEL FOR
More informationGood Vibes: Introduction to Oscillations
Chapter 14 Solutions Good Vibes: Introduction to Oscillations Description: Several conceptual and qualitative questions related to main characteristics of simple harmonic motion: amplitude, displacement,
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationMechanics of Motor Proteins and the Cytoskeleton Jonathon Howard Chapter 10 Force generation 2 nd part. Andrea and Yinyun April 4 th,2012
Mechanics of Motor Proteins and the Cytoskeleton Jonathon Howard Chapter 10 Force generation 2 nd part Andrea and Yinyun April 4 th,2012 I. Equilibrium Force Reminder: http://www.youtube.com/watch?v=yt59kx_z6xm
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationPhys101 Lectures 9 and 10 Conservation of Mechanical Energy
Phys101 Lectures 9 and 10 Conservation of Mechanical Energy Key points: Conservative and Nonconservative Forces Potential Energy Generalized work-energy principle Mechanical Energy and Its Conservation
More informationTitleElastically coupled two-dimensional. Citation PHYSICAL REVIEW E (2005), 71(6)
TitleElastically coupled two-dimensional Author(s) Goko, H; Igarashi, A Citation PHYSICAL REVIEW E (2005), 71(6) Issue Date 2005-06 URL http://hdl.handle.net/2433/50298 RightCopyright 2005 American Physical
More informationChapter 6- An Introduction to Metabolism*
Chapter 6- An Introduction to Metabolism* *Lecture notes are to be used as a study guide only and do not represent the comprehensive information you will need to know for the exams. The Energy of Life
More informationThe ability to determine protein structure has made great advances in recent years
Weighted-Ensemble Brownian Dynamics Simulation Amanda McCoy, Harvard University Daniel M. Zuckerman, Ph.D., University of Pittsburgh, Department of Computational Biology I. Introduction The ability to
More informationPHYSIOLOGY CHAPTER 9 MUSCLE TISSUE Fall 2016
PHYSIOLOGY CHAPTER 9 MUSCLE TISSUE Fall 2016 2 Chapter 9 Muscles and Muscle Tissue Overview of Muscle Tissue types of muscle: are all prefixes for muscle Contractility all muscles cells can Smooth & skeletal
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationSECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C7: BIOLOGICAL PHYSICS
2757 SECOND PUBLIC EXAMINATION Honour School of Physics Part C: 4 Year Course Honour School of Physics and Philosophy Part C C7: BIOLOGICAL PHYSICS TRINITY TERM 2013 Monday, 17 June, 2.30 pm 5.45 pm 15
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationThe... of a particle is defined as its change in position in some time interval.
Distance is the. of a path followed by a particle. Distance is a quantity. The... of a particle is defined as its change in position in some time interval. Displacement is a.. quantity. The... of a particle
More informationAll questions are of equal value. No marks are subtracted for wrong answers.
(1:30 PM 4:30 PM) Page 1 of 6 All questions are of equal value. No marks are subtracted for wrong answers. Record all answers on the computer score sheet provided. USE PENCIL ONLY! Black pen will look
More informationName Lesson 7. Homework Work and Energy Problem Solving Outcomes
Physics 1 Name Lesson 7. Homework Work and Energy Problem Solving Outcomes Date 1. Define work. 2. Define energy. 3. Determine the work done by a constant force. Period 4. Determine the work done by a
More informationAP Physics C Mechanics
1 AP Physics C Mechanics Simple Harmonic Motion 2015 12 05 www.njctl.org 2 Table of Contents Click on the topic to go to that section Spring and a Block Energy of SHM SHM and UCM Simple and Physical Pendulums
More informationBioelectricity Prof. Mainak Das Department of Biological Sciences, and Bioengineering Indian Institute of Technology, Kanpur.
Bioelectricity Prof. Mainak Das Department of Biological Sciences, and Bioengineering Indian Institute of Technology, Kanpur Lecture 17 Welcome back to the bioelectricity lecture, series. So, in the last
More informationMODEL PAPER CLASS XI PHYSICS (GROUP 1) BLUEPRINT Name of chapter (1)
sr. no. MODEL PAPER CLASS XI PHYSICS (GROUP ) BLUEPRINT Name of chapter VSAQ () SA-I (2) SA-II (3) Value based (4) LA(5) Total 70 Physical world and measurement 3 2 Kinematics 2 3,3 5 3 Laws of motion
More informationChapters 10 & 11: Energy
Chapters 10 & 11: Energy Power: Sources of Energy Tidal Power SF Bay Tidal Power Project Main Ideas (Encyclopedia of Physics) Energy is an abstract quantity that an object is said to possess. It is not
More informationPHYSICS. Chapter 5 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 5 Lecture RANDALL D. KNIGHT Chapter 5 Force and Motion IN THIS CHAPTER, you will learn about the connection between force and motion.
More informationWorksheet #05 Kinetic Energy-Work Theorem
Physics Summer 08 Worksheet #05 June. 8, 08. A 0-kg crate is pulled 5 m up along a frictionless incline as shown in the figure below. The crate starts at rest and has a final speed of 6.0 m/s. (a) Draw
More informationPRACTICE QUESTION PAPER WITH SOLUTION CLASS XI PHYSICS
PRACTICE QUESTION PAPER WITH SOLUTION CLASS XI PHYSICS. A given quantity has both magnitude and direction. Is it necessarily a vector? Justify your answer.. What is the rotational analogue of the force?.
More informationSupramolecular Chemistry of Nanomaterials
Supramolecular Chemistry of Nanomaterials Joachim Steinke Ramon Vilar Lecture 6 Towards the Development of Molecular Machines Department of Chemistry Imperial College of Science, Technology and Medicine
More informationPower: Sources of Energy
Chapter 5 Energy Power: Sources of Energy Tidal Power SF Bay Tidal Power Project Main Ideas (Encyclopedia of Physics) Energy is an abstract quantity that an object is said to possess. It is not something
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
More informationLecture 4: viscoelasticity and cell mechanics
Teaser movie: flexible robots! R. Shepherd, Whitesides group, Harvard 1 Lecture 4: viscoelasticity and cell mechanics S-RSI Physics Lectures: Soft Condensed Matter Physics Jacinta C. Conrad University
More informationPhysics 9 Spring 2012 Midterm 1 Solutions
Physics 9 Spring 22 NAME: TA: Physics 9 Spring 22 Midterm s For the midterm, you may use one sheet of notes with whatever you want to put on it, front and back. Please sit every other seat, and please
More informationAP Physics Free Response Practice Oscillations
AP Physics Free Response Practice Oscillations 1975B7. A pendulum consists of a small object of mass m fastened to the end of an inextensible cord of length L. Initially, the pendulum is drawn aside through
More information