Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
SlideShare a Scribd company logo
DYNAMICS OF ACTIN FILAMENTS IN
THE CONTRACTILE RING
Summer Internship Report
Submitted to:
Bio-Physics Group, AG Karsten Kruse
University of Saarland , Germany
Submitted by:
Prafull Kumar Sharma
2nd year undergraduate B.Tech .(Engineering Physics)
Indian Institute of Technology (Delhi), India
Numerical Analysis of Diffusion
equation
• First I was supposed to learn basic techniques used in
numerical analysis to analyse solutions of a diffusion
equation.We use Numerical Analysis to understand the
behaviour of solution of the equations involved in our
project.To begin with ,I solved 1-D diffusion equation using
backward euler and forwad euler method algorithms to
show on computer.I have used MATLAB to show solutions
for the equation.Here are some screenshots of theoretical
methods and equations involved.In the screenshots,I have
solved diffusion equation analytically with initial condition
as dirac delta function.Solution is a gaussian distribution as
expected theoretically when solved numerically on
MATLAB.
Solution of 1D diffusion equation using delta function
as initial condition
Gaussian Distribution as expected with theoretical analysis
Diffusion Equation and it’s analytical solution with
initial condition (c(x,0)) as “dirac delta function”
Theory Behind Forward Euler Method
Theoretical explanation behind Backward Euler Method
Numerical Analysis of Dynamics of
actin filaments without considering
Bipolar Filament
• After Basic Numerical analysis techniques,I was supposed to
analyse dynamics of actin filaments as described in “Actively
Contracting Bundles of Polar Filaments” by K.Kruse && F. Jullicher,
published in Physical Review Letters Volume 85,Number 8.
• In this Model, we consider (with our theoretical considerations) we
try to capture essential features of the ring dynamics, such as,
filament polarity, interaction between filaments through protein
motors. Here we assume that actin filaments align with perimeter
of ring. We denote the co-ordinate along the ring perimeter by x
and describe the distribution of (polar) actin filaments with respect
to x coordinates by the densities c+ (x,t) for filaments with their
plus-end pointing clockwise and c- (x,t) for filaments of the opposite
orientation. Here are some screenshots of Theoretical explanations
and equations involved.
Actin Dynamics Equation without bipolar filament with
consideration of treadmilling
Understanding Dynamics of Filaments
Adding Preturbations to analyse stability of solutions
Fourier Transform of equation governing
dynamics of filaments
Matrix Elements which constitutes the matrix A and
Real part of eigenvalues of A is crucial for stability for
these steady states.
Graph of alpha_c vs beta for L=5 with no treadmilling
In this graph,we have dimensionalised α and β using length of
filament and diffusion constant.
C0
+ =.3 C0
- =.7,
Graph of alpha_c vs beta for L=5 with v_tr =.05
C0
+ =.3 C0
- =.7,
Graph of alpha_c vs beta for L=10 with no treadmilling
C0
+ =.3 C0
- =.7,
Graph of alpha_c vs beta for L=10 with v_tr =.05
C0
+ =.3 C0
- =.7,
Numerical Analysis of Dynamics of
actin filaments with consideration of
Bipolar Filament
• After Basic Numerical analysis techniques,I was supposed to analyse
dynamics of actin filaments as described in “self-organization and
mechanical properties of active filament bundles” by K.Kruse && F.
Jullicher, published in Physical Review E 67, 051913 (2003)
• In this Model, we consider (with our theoretical considerations) we try to
capture essential features of the ring dynamics, such as, filament polarity,
interaction between filaments through protein motors. Here we assume
that actin filaments align with perimeter of ring. We denote the co-
ordinate along the ring perimeter by x and describe the distribution of
(polar) actin filaments with respect to x coordinates by the densities c+
(x,t) for filaments with their plus-end pointing clockwise and c- (x,t) for
filaments of the opposite orientation. The distribution of bipolar filaments
is denoted by cbp (x,t) giving the density of their centers. In this wd is rate
of breaking of bipolar ones and wc is rate of combination of two filaments.
Here are some screenshots of Theoretical explanations and equations
involved.
Dynamics of actin filaments in the contractile ring
Notions and conditions used in Dynamics
Actin Dynamics Equation with consideration of bipolar
filament
It is sufficient to check for k=2*pi/L for
critical values of α for given β
• As it seemed while observing graphs , we
don’t need all values of k ( wave number
arising from fourier analysis) to check for
critical values for α vs β.As it is evident from
coding(“numerically” or “graphically”) that we
need only k=2*pi/L as maxmum of
eigenvalues of matrix that I got for stability
analysis is always decreasing with respect to k
for given α and β.
Max . of eigenvalue is decreasing w.r.t. k.
That’s why study of k=2*pi/L is sufficient as evident from graph.
alpha_c vs beta without treadmilling with w_c=0
v=0.0; C0
+ =.3 C0
- =.7, L=5,D=1
alpha_c vs beta without treadmilling
v=0.0; C0
+ =.3 C0
- =.7, L=10,D=1
alpha_c vs beta with treadmilling
v=0.5; C0
+ =.3 C0
- =.7, L=10,D=1
alpha_c vs Beta without treadmilling
v=0; C0
+ =.3 C0
- =.7, L=5,D=1
alpha_c vs Beta with treadmilling
v=0.5; C0
+ =.3 C0
- =.7, L=5,D=1
alpha_c vs treadmilling velocity with β=.5
C0
+ =.3 C0
- =.7, L=10,D=1
alpha_c vs treadmilling velocity with β=.5
C0
+ =.3 C0
- =.7, L=5,D=1
alpha_c v/s w_c with observation of the values in
stable region for which steady state will be stationary
with β= 0.5.
For α≤αc , solutions will be stable. All the lines that are inside this region
will tend give non oscillatory stable steady state solutions.
Here in this graph ,0<α<5,-0.8<wc<4 ; In this graph, descretization
number is shown on y-axis and x-axis.
C0
+ =.3 C0
- =.7, L=10,D=1,v=0
alpha_c v/s w_c with observation of the values in
stable region for which steady state will be stationary
with β= 0.5.
For α≤αc , solutions will be stable. All the lines that are inside this region
will tend give non oscillatory stable steady state solutions.
Here in this graph ,0<α<5,-0.8<wc<4 ; In this graph, descretization
number is shown on y-axis and x-axis.
C0
+ =.3 C0
- =.7, L=10,D=1,v=0.5
Numerical Solutions of Dynamics of
actin filaments without considering
Bipolar Filament
• After doing stability analysis, I solved the actin
dynamics equation numerically using first
order upwind scheme with adaptive time
stepping . Here is a snapshot of theoretical
explanations behind it.
Dynamics of actin filaments in the contractile ring
First Order Upwind scheme for actin dynamics scheme
For solutions of actin dynamics equation without
bipolar filament
c2(1:N,1)=0.3*ones(1,N).*(1+rand(1,N));
c1(1:N,1)=0.7*ones(1,N).*(1+rand(1,N));
L=10;a=0.6;b=2;
For solutions of actin dynamics equation with
consideration of bipolar filament
D=1;L=10;a=0.6;b=2;w1=0.3;w2=0.7;
c1(1:N,1)=0.7*ones(1,N).*(1+rand(1,N));
c2(1:N,1)= 0.3*ones(1,N).*(1+rand(1,N));
c3(1:N,1)=0.09*ones(1,N).*(1+rand(1,N));
Results
• It is sufficient to check for k=2*pi/L for
stability analysis.(refer page22)
• Upwind scheme verifies the results obtained
in stability analysis.
• Verification of fact that solution is unstable for
D*dt/(dx)2 >0.5 for forward euler case.
• Stability patterns obtained are in consensus
with expected data as described in paper.
Experiences
• During the project , I got introduced to various
techniques to solve Non linear Equations. I also
got to learn about programming on MATLAB.
• Back to project, so far I have done numerical
analysis of solutions for actin dynamics
equations. As we have discussed, I am currently
working on stress calculations.
I hope I was good during internship!!!
It was my first research internship and I have
learnt a lot from you. Thanks for your guidance!!!

More Related Content

What's hot

Computational electromagnetics
Computational electromagneticsComputational electromagnetics
Computational electromagnetics
Awaab Fakih
 
Numerical Methods in Mechanical Engineering - Final Project
Numerical Methods in Mechanical Engineering - Final ProjectNumerical Methods in Mechanical Engineering - Final Project
Numerical Methods in Mechanical Engineering - Final Project
Stasik Nemirovsky
 
Statistical Physics Assignment Help
Statistical Physics Assignment HelpStatistical Physics Assignment Help
Statistical Physics Assignment Help
Statistics Assignment Help
 
Chapter 5
Chapter 5Chapter 5
Mom slideshow
Mom slideshowMom slideshow
Mom slideshow
ashusuzie
 
Diffusion Homework Help
Diffusion Homework HelpDiffusion Homework Help
Diffusion Homework Help
Statistics Assignment Help
 
Quantum algorithm for solving linear systems of equations
 Quantum algorithm for solving linear systems of equations Quantum algorithm for solving linear systems of equations
Quantum algorithm for solving linear systems of equations
XequeMateShannon
 
Mechanical Engineering Assignment Help
Mechanical Engineering Assignment HelpMechanical Engineering Assignment Help
Mechanical Engineering Assignment Help
Matlab Assignment Experts
 
Interactives Methods
Interactives MethodsInteractives Methods
Interactives Methods
UIS
 
Lecture 9 dim & rank - 4-5 & 4-6
Lecture 9   dim & rank -  4-5 & 4-6Lecture 9   dim & rank -  4-5 & 4-6
Lecture 9 dim & rank - 4-5 & 4-6
njit-ronbrown
 
Numerical Solution of Diffusion Equation by Finite Difference Method
Numerical Solution of Diffusion Equation by Finite Difference MethodNumerical Solution of Diffusion Equation by Finite Difference Method
Numerical Solution of Diffusion Equation by Finite Difference Method
iosrjce
 
Physics Assignment Help
Physics Assignment Help Physics Assignment Help
Physics Assignment Help
Edu Assignment Help
 
Lecture 8 nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
Lecture 8   nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6Lecture 8   nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
Lecture 8 nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
njit-ronbrown
 
S6 l04 analytical and numerical methods of structural analysis
S6 l04 analytical and numerical methods of structural analysisS6 l04 analytical and numerical methods of structural analysis
S6 l04 analytical and numerical methods of structural analysis
Shaikh Mohsin
 
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMSFINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
roymeister007
 
Fourier Transform Assignment Help
Fourier Transform Assignment HelpFourier Transform Assignment Help
Fourier Transform Assignment Help
Matlab Assignment Experts
 
Iterative methods
Iterative methodsIterative methods
Iterative methods
Kt Silva
 
Integral Transform
Integral  TransformIntegral  Transform
Integral Transform
SheharBano31
 
FDM Numerical solution of Laplace Equation using MATLAB
FDM Numerical solution of Laplace Equation using MATLABFDM Numerical solution of Laplace Equation using MATLAB
FDM Numerical solution of Laplace Equation using MATLAB
Aya Zaki
 
Application of Numerical Methods (Finite Difference) in Heat Transfer
Application of Numerical Methods (Finite Difference) in Heat TransferApplication of Numerical Methods (Finite Difference) in Heat Transfer
Application of Numerical Methods (Finite Difference) in Heat Transfer
Shivshambhu Kumar
 

What's hot (20)

Computational electromagnetics
Computational electromagneticsComputational electromagnetics
Computational electromagnetics
 
Numerical Methods in Mechanical Engineering - Final Project
Numerical Methods in Mechanical Engineering - Final ProjectNumerical Methods in Mechanical Engineering - Final Project
Numerical Methods in Mechanical Engineering - Final Project
 
Statistical Physics Assignment Help
Statistical Physics Assignment HelpStatistical Physics Assignment Help
Statistical Physics Assignment Help
 
Chapter 5
Chapter 5Chapter 5
Chapter 5
 
Mom slideshow
Mom slideshowMom slideshow
Mom slideshow
 
Diffusion Homework Help
Diffusion Homework HelpDiffusion Homework Help
Diffusion Homework Help
 
Quantum algorithm for solving linear systems of equations
 Quantum algorithm for solving linear systems of equations Quantum algorithm for solving linear systems of equations
Quantum algorithm for solving linear systems of equations
 
Mechanical Engineering Assignment Help
Mechanical Engineering Assignment HelpMechanical Engineering Assignment Help
Mechanical Engineering Assignment Help
 
Interactives Methods
Interactives MethodsInteractives Methods
Interactives Methods
 
Lecture 9 dim & rank - 4-5 & 4-6
Lecture 9   dim & rank -  4-5 & 4-6Lecture 9   dim & rank -  4-5 & 4-6
Lecture 9 dim & rank - 4-5 & 4-6
 
Numerical Solution of Diffusion Equation by Finite Difference Method
Numerical Solution of Diffusion Equation by Finite Difference MethodNumerical Solution of Diffusion Equation by Finite Difference Method
Numerical Solution of Diffusion Equation by Finite Difference Method
 
Physics Assignment Help
Physics Assignment Help Physics Assignment Help
Physics Assignment Help
 
Lecture 8 nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
Lecture 8   nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6Lecture 8   nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
Lecture 8 nul col bases dim & rank - section 4-2, 4-3, 4-5 & 4-6
 
S6 l04 analytical and numerical methods of structural analysis
S6 l04 analytical and numerical methods of structural analysisS6 l04 analytical and numerical methods of structural analysis
S6 l04 analytical and numerical methods of structural analysis
 
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMSFINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
FINITE DIFFERENCE MODELLING FOR HEAT TRANSFER PROBLEMS
 
Fourier Transform Assignment Help
Fourier Transform Assignment HelpFourier Transform Assignment Help
Fourier Transform Assignment Help
 
Iterative methods
Iterative methodsIterative methods
Iterative methods
 
Integral Transform
Integral  TransformIntegral  Transform
Integral Transform
 
FDM Numerical solution of Laplace Equation using MATLAB
FDM Numerical solution of Laplace Equation using MATLABFDM Numerical solution of Laplace Equation using MATLAB
FDM Numerical solution of Laplace Equation using MATLAB
 
Application of Numerical Methods (Finite Difference) in Heat Transfer
Application of Numerical Methods (Finite Difference) in Heat TransferApplication of Numerical Methods (Finite Difference) in Heat Transfer
Application of Numerical Methods (Finite Difference) in Heat Transfer
 

Similar to Dynamics of actin filaments in the contractile ring

DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptxDIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
SubhKr1
 
Numerical Solution of Linear algebraic Equation
Numerical Solution of Linear algebraic EquationNumerical Solution of Linear algebraic Equation
Numerical Solution of Linear algebraic Equation
payalpriyadarshinisa1
 
Signals and Systems Assignment Help
Signals and Systems Assignment HelpSignals and Systems Assignment Help
Signals and Systems Assignment Help
Matlab Assignment Experts
 
Scilab for real dummies j.heikell - part 2
Scilab for real dummies j.heikell - part 2Scilab for real dummies j.heikell - part 2
Scilab for real dummies j.heikell - part 2
Scilab
 
Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05
Charlton Inao
 
project report(1)
project report(1)project report(1)
project report(1)
Conor Bradley
 
circuit_modes_v5
circuit_modes_v5circuit_modes_v5
circuit_modes_v5
Olivier Buu
 
Determinant untuk kuliahteknik sipil atau umum
Determinant untuk kuliahteknik sipil atau umumDeterminant untuk kuliahteknik sipil atau umum
Determinant untuk kuliahteknik sipil atau umum
adjie25
 
Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05
Charlton Inao
 
Matrices ppt
Matrices pptMatrices ppt
Matrices ppt
aakashray33
 
Ijetr021210
Ijetr021210Ijetr021210
Ijetr021210
ER Publication.org
 
Physical Chemistry Exam Help
Physical Chemistry Exam HelpPhysical Chemistry Exam Help
Physical Chemistry Exam Help
Live Exam Helper
 
Signal Processing Assignment Help
Signal Processing Assignment HelpSignal Processing Assignment Help
Signal Processing Assignment Help
Matlab Assignment Experts
 
MODULE_05-Matrix Decomposition.pptx
MODULE_05-Matrix Decomposition.pptxMODULE_05-Matrix Decomposition.pptx
MODULE_05-Matrix Decomposition.pptx
AlokSingh205089
 
A some basic rules of tensor calculus
A some basic rules of tensor calculusA some basic rules of tensor calculus
A some basic rules of tensor calculus
Tarun Gehlot
 
Lesson 8: Determinants III
Lesson 8: Determinants IIILesson 8: Determinants III
Lesson 8: Determinants III
Matthew Leingang
 
Mtc ssample05
Mtc ssample05Mtc ssample05
Mtc ssample05
bikram ...
 
Mtc ssample05
Mtc ssample05Mtc ssample05
Mtc ssample05
bikram ...
 
Mat 223_Ch5-Eigenvalues.ppt
Mat 223_Ch5-Eigenvalues.pptMat 223_Ch5-Eigenvalues.ppt
Mat 223_Ch5-Eigenvalues.ppt
SatyabrataMandal7
 
Stability and pole location
Stability and pole locationStability and pole location
Stability and pole location
ssuser5d64bb
 

Similar to Dynamics of actin filaments in the contractile ring (20)

DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptxDIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
DIP_Lecture-35_36_RKJ_Interpolation_Resampling_Unitary_Transformations1.pptx
 
Numerical Solution of Linear algebraic Equation
Numerical Solution of Linear algebraic EquationNumerical Solution of Linear algebraic Equation
Numerical Solution of Linear algebraic Equation
 
Signals and Systems Assignment Help
Signals and Systems Assignment HelpSignals and Systems Assignment Help
Signals and Systems Assignment Help
 
Scilab for real dummies j.heikell - part 2
Scilab for real dummies j.heikell - part 2Scilab for real dummies j.heikell - part 2
Scilab for real dummies j.heikell - part 2
 
Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05
 
project report(1)
project report(1)project report(1)
project report(1)
 
circuit_modes_v5
circuit_modes_v5circuit_modes_v5
circuit_modes_v5
 
Determinant untuk kuliahteknik sipil atau umum
Determinant untuk kuliahteknik sipil atau umumDeterminant untuk kuliahteknik sipil atau umum
Determinant untuk kuliahteknik sipil atau umum
 
Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05Wk 6 part 2 non linearites and non linearization april 05
Wk 6 part 2 non linearites and non linearization april 05
 
Matrices ppt
Matrices pptMatrices ppt
Matrices ppt
 
Ijetr021210
Ijetr021210Ijetr021210
Ijetr021210
 
Physical Chemistry Exam Help
Physical Chemistry Exam HelpPhysical Chemistry Exam Help
Physical Chemistry Exam Help
 
Signal Processing Assignment Help
Signal Processing Assignment HelpSignal Processing Assignment Help
Signal Processing Assignment Help
 
MODULE_05-Matrix Decomposition.pptx
MODULE_05-Matrix Decomposition.pptxMODULE_05-Matrix Decomposition.pptx
MODULE_05-Matrix Decomposition.pptx
 
A some basic rules of tensor calculus
A some basic rules of tensor calculusA some basic rules of tensor calculus
A some basic rules of tensor calculus
 
Lesson 8: Determinants III
Lesson 8: Determinants IIILesson 8: Determinants III
Lesson 8: Determinants III
 
Mtc ssample05
Mtc ssample05Mtc ssample05
Mtc ssample05
 
Mtc ssample05
Mtc ssample05Mtc ssample05
Mtc ssample05
 
Mat 223_Ch5-Eigenvalues.ppt
Mat 223_Ch5-Eigenvalues.pptMat 223_Ch5-Eigenvalues.ppt
Mat 223_Ch5-Eigenvalues.ppt
 
Stability and pole location
Stability and pole locationStability and pole location
Stability and pole location
 

Recently uploaded

Founders Of Modern Science 16th Century to the 21st Century.pdf
Founders Of Modern Science 16th Century to the 21st Century.pdfFounders Of Modern Science 16th Century to the 21st Century.pdf
Founders Of Modern Science 16th Century to the 21st Century.pdf
Steven Camilleri
 
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
Sérgio Sacani
 
Principles of Colorimetry - Mastering the Art of Colour Measurement
Principles of Colorimetry - Mastering the Art of Colour MeasurementPrinciples of Colorimetry - Mastering the Art of Colour Measurement
Principles of Colorimetry - Mastering the Art of Colour Measurement
Colours Guide
 
Cause and solution of Water hyacinth (Terror of Bengal)
Cause and solution of Water hyacinth (Terror of Bengal)Cause and solution of Water hyacinth (Terror of Bengal)
Cause and solution of Water hyacinth (Terror of Bengal)
saloniswain225
 
Science grade 09 Lesson1-2 NLC-pptx.pptx
Science grade 09 Lesson1-2 NLC-pptx.pptxScience grade 09 Lesson1-2 NLC-pptx.pptx
Science grade 09 Lesson1-2 NLC-pptx.pptx
JoanaBanasen1
 
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
ThrombUS+ Project
 
Mathematics 7th class English medium text book
Mathematics 7th class English medium text bookMathematics 7th class English medium text book
Mathematics 7th class English medium text book
AnilKumar210119
 
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDYSUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
Dr Kirpa Ram Jangra
 
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
Steffi Friedrichs
 
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
Sérgio Sacani
 
What is data analytics ,Data science,Data processing chain,regression,decisio...
What is data analytics ,Data science,Data processing chain,regression,decisio...What is data analytics ,Data science,Data processing chain,regression,decisio...
What is data analytics ,Data science,Data processing chain,regression,decisio...
hannahroseline2
 
@Call @Girls in Asansol 🐱‍🐉 XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
 @Call @Girls in Asansol 🐱‍🐉  XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A... @Call @Girls in Asansol 🐱‍🐉  XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
@Call @Girls in Asansol 🐱‍🐉 XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
JoyaBansal
 
Comparison of RNA Viruses. docx
Comparison of RNA Viruses.          docxComparison of RNA Viruses.          docx
Comparison of RNA Viruses. docx
Dr Sumitha Jagadibabu
 
Gijubhai Badheka bed 1st year pppt presentation
Gijubhai Badheka bed 1st year pppt presentationGijubhai Badheka bed 1st year pppt presentation
Gijubhai Badheka bed 1st year pppt presentation
PRITIKUMARI117
 
ThrombUS+ Project Overview - January 2024
ThrombUS+ Project Overview - January 2024ThrombUS+ Project Overview - January 2024
ThrombUS+ Project Overview - January 2024
ThrombUS+ Project
 
Comparative Accounts of Selaginella And Isoetes
Comparative Accounts of Selaginella And IsoetesComparative Accounts of Selaginella And Isoetes
Comparative Accounts of Selaginella And Isoetes
saloniswain225
 
Capparidaceae ( Caper family) from Ahmedabad
Capparidaceae ( Caper family) from AhmedabadCapparidaceae ( Caper family) from Ahmedabad
Capparidaceae ( Caper family) from Ahmedabad
saloniswain225
 
ThrombUS+ Gender Balance Awareness - May 2024
ThrombUS+ Gender Balance Awareness - May 2024ThrombUS+ Gender Balance Awareness - May 2024
ThrombUS+ Gender Balance Awareness - May 2024
ThrombUS+ Project
 
ThrombUS+ Project Presentation - June 2024
ThrombUS+ Project Presentation - June 2024ThrombUS+ Project Presentation - June 2024
ThrombUS+ Project Presentation - June 2024
elenikaldoudi1
 
SPERM FUNCTION TEST IN EMBRYOLOGY .pptx
SPERM FUNCTION TEST  IN EMBRYOLOGY .pptxSPERM FUNCTION TEST  IN EMBRYOLOGY .pptx
SPERM FUNCTION TEST IN EMBRYOLOGY .pptx
SRI AUROBINDO UNIVERSITY
 

Recently uploaded (20)

Founders Of Modern Science 16th Century to the 21st Century.pdf
Founders Of Modern Science 16th Century to the 21st Century.pdfFounders Of Modern Science 16th Century to the 21st Century.pdf
Founders Of Modern Science 16th Century to the 21st Century.pdf
 
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
The X‐Pattern Merging of the Equatorial IonizationAnomaly Crests During Geoma...
 
Principles of Colorimetry - Mastering the Art of Colour Measurement
Principles of Colorimetry - Mastering the Art of Colour MeasurementPrinciples of Colorimetry - Mastering the Art of Colour Measurement
Principles of Colorimetry - Mastering the Art of Colour Measurement
 
Cause and solution of Water hyacinth (Terror of Bengal)
Cause and solution of Water hyacinth (Terror of Bengal)Cause and solution of Water hyacinth (Terror of Bengal)
Cause and solution of Water hyacinth (Terror of Bengal)
 
Science grade 09 Lesson1-2 NLC-pptx.pptx
Science grade 09 Lesson1-2 NLC-pptx.pptxScience grade 09 Lesson1-2 NLC-pptx.pptx
Science grade 09 Lesson1-2 NLC-pptx.pptx
 
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
Towards Wearable Continuous Point-of-Care Monitoring for Deep Vein Thrombosis...
 
Mathematics 7th class English medium text book
Mathematics 7th class English medium text bookMathematics 7th class English medium text book
Mathematics 7th class English medium text book
 
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDYSUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
SUBJECT SPECIFIC ETHICAL ISSUES IN STUDY
 
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
MACRAMÉ-ChiPs: Patchwork Project Family & Sibling Projects (24th Meeting of t...
 
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
The extremotolerant desert moss Syntrichia caninervis is a promising pioneer ...
 
What is data analytics ,Data science,Data processing chain,regression,decisio...
What is data analytics ,Data science,Data processing chain,regression,decisio...What is data analytics ,Data science,Data processing chain,regression,decisio...
What is data analytics ,Data science,Data processing chain,regression,decisio...
 
@Call @Girls in Asansol 🐱‍🐉 XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
 @Call @Girls in Asansol 🐱‍🐉  XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A... @Call @Girls in Asansol 🐱‍🐉  XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
@Call @Girls in Asansol 🐱‍🐉 XXXXXXXXXX 🐱‍🐉 Tanisha Sharma Best High Class A...
 
Comparison of RNA Viruses. docx
Comparison of RNA Viruses.          docxComparison of RNA Viruses.          docx
Comparison of RNA Viruses. docx
 
Gijubhai Badheka bed 1st year pppt presentation
Gijubhai Badheka bed 1st year pppt presentationGijubhai Badheka bed 1st year pppt presentation
Gijubhai Badheka bed 1st year pppt presentation
 
ThrombUS+ Project Overview - January 2024
ThrombUS+ Project Overview - January 2024ThrombUS+ Project Overview - January 2024
ThrombUS+ Project Overview - January 2024
 
Comparative Accounts of Selaginella And Isoetes
Comparative Accounts of Selaginella And IsoetesComparative Accounts of Selaginella And Isoetes
Comparative Accounts of Selaginella And Isoetes
 
Capparidaceae ( Caper family) from Ahmedabad
Capparidaceae ( Caper family) from AhmedabadCapparidaceae ( Caper family) from Ahmedabad
Capparidaceae ( Caper family) from Ahmedabad
 
ThrombUS+ Gender Balance Awareness - May 2024
ThrombUS+ Gender Balance Awareness - May 2024ThrombUS+ Gender Balance Awareness - May 2024
ThrombUS+ Gender Balance Awareness - May 2024
 
ThrombUS+ Project Presentation - June 2024
ThrombUS+ Project Presentation - June 2024ThrombUS+ Project Presentation - June 2024
ThrombUS+ Project Presentation - June 2024
 
SPERM FUNCTION TEST IN EMBRYOLOGY .pptx
SPERM FUNCTION TEST  IN EMBRYOLOGY .pptxSPERM FUNCTION TEST  IN EMBRYOLOGY .pptx
SPERM FUNCTION TEST IN EMBRYOLOGY .pptx
 

Dynamics of actin filaments in the contractile ring

  • 1. DYNAMICS OF ACTIN FILAMENTS IN THE CONTRACTILE RING Summer Internship Report Submitted to: Bio-Physics Group, AG Karsten Kruse University of Saarland , Germany Submitted by: Prafull Kumar Sharma 2nd year undergraduate B.Tech .(Engineering Physics) Indian Institute of Technology (Delhi), India
  • 2. Numerical Analysis of Diffusion equation • First I was supposed to learn basic techniques used in numerical analysis to analyse solutions of a diffusion equation.We use Numerical Analysis to understand the behaviour of solution of the equations involved in our project.To begin with ,I solved 1-D diffusion equation using backward euler and forwad euler method algorithms to show on computer.I have used MATLAB to show solutions for the equation.Here are some screenshots of theoretical methods and equations involved.In the screenshots,I have solved diffusion equation analytically with initial condition as dirac delta function.Solution is a gaussian distribution as expected theoretically when solved numerically on MATLAB.
  • 3. Solution of 1D diffusion equation using delta function as initial condition Gaussian Distribution as expected with theoretical analysis
  • 4. Diffusion Equation and it’s analytical solution with initial condition (c(x,0)) as “dirac delta function”
  • 5. Theory Behind Forward Euler Method
  • 6. Theoretical explanation behind Backward Euler Method
  • 7. Numerical Analysis of Dynamics of actin filaments without considering Bipolar Filament • After Basic Numerical analysis techniques,I was supposed to analyse dynamics of actin filaments as described in “Actively Contracting Bundles of Polar Filaments” by K.Kruse && F. Jullicher, published in Physical Review Letters Volume 85,Number 8. • In this Model, we consider (with our theoretical considerations) we try to capture essential features of the ring dynamics, such as, filament polarity, interaction between filaments through protein motors. Here we assume that actin filaments align with perimeter of ring. We denote the co-ordinate along the ring perimeter by x and describe the distribution of (polar) actin filaments with respect to x coordinates by the densities c+ (x,t) for filaments with their plus-end pointing clockwise and c- (x,t) for filaments of the opposite orientation. Here are some screenshots of Theoretical explanations and equations involved.
  • 8. Actin Dynamics Equation without bipolar filament with consideration of treadmilling
  • 10. Adding Preturbations to analyse stability of solutions
  • 11. Fourier Transform of equation governing dynamics of filaments
  • 12. Matrix Elements which constitutes the matrix A and Real part of eigenvalues of A is crucial for stability for these steady states.
  • 13. Graph of alpha_c vs beta for L=5 with no treadmilling In this graph,we have dimensionalised α and β using length of filament and diffusion constant. C0 + =.3 C0 - =.7,
  • 14. Graph of alpha_c vs beta for L=5 with v_tr =.05 C0 + =.3 C0 - =.7,
  • 15. Graph of alpha_c vs beta for L=10 with no treadmilling C0 + =.3 C0 - =.7,
  • 16. Graph of alpha_c vs beta for L=10 with v_tr =.05 C0 + =.3 C0 - =.7,
  • 17. Numerical Analysis of Dynamics of actin filaments with consideration of Bipolar Filament • After Basic Numerical analysis techniques,I was supposed to analyse dynamics of actin filaments as described in “self-organization and mechanical properties of active filament bundles” by K.Kruse && F. Jullicher, published in Physical Review E 67, 051913 (2003) • In this Model, we consider (with our theoretical considerations) we try to capture essential features of the ring dynamics, such as, filament polarity, interaction between filaments through protein motors. Here we assume that actin filaments align with perimeter of ring. We denote the co- ordinate along the ring perimeter by x and describe the distribution of (polar) actin filaments with respect to x coordinates by the densities c+ (x,t) for filaments with their plus-end pointing clockwise and c- (x,t) for filaments of the opposite orientation. The distribution of bipolar filaments is denoted by cbp (x,t) giving the density of their centers. In this wd is rate of breaking of bipolar ones and wc is rate of combination of two filaments. Here are some screenshots of Theoretical explanations and equations involved.
  • 19. Notions and conditions used in Dynamics
  • 20. Actin Dynamics Equation with consideration of bipolar filament
  • 21. It is sufficient to check for k=2*pi/L for critical values of α for given β • As it seemed while observing graphs , we don’t need all values of k ( wave number arising from fourier analysis) to check for critical values for α vs β.As it is evident from coding(“numerically” or “graphically”) that we need only k=2*pi/L as maxmum of eigenvalues of matrix that I got for stability analysis is always decreasing with respect to k for given α and β.
  • 22. Max . of eigenvalue is decreasing w.r.t. k. That’s why study of k=2*pi/L is sufficient as evident from graph.
  • 23. alpha_c vs beta without treadmilling with w_c=0 v=0.0; C0 + =.3 C0 - =.7, L=5,D=1
  • 24. alpha_c vs beta without treadmilling v=0.0; C0 + =.3 C0 - =.7, L=10,D=1
  • 25. alpha_c vs beta with treadmilling v=0.5; C0 + =.3 C0 - =.7, L=10,D=1
  • 26. alpha_c vs Beta without treadmilling v=0; C0 + =.3 C0 - =.7, L=5,D=1
  • 27. alpha_c vs Beta with treadmilling v=0.5; C0 + =.3 C0 - =.7, L=5,D=1
  • 28. alpha_c vs treadmilling velocity with β=.5 C0 + =.3 C0 - =.7, L=10,D=1
  • 29. alpha_c vs treadmilling velocity with β=.5 C0 + =.3 C0 - =.7, L=5,D=1
  • 30. alpha_c v/s w_c with observation of the values in stable region for which steady state will be stationary with β= 0.5. For α≤αc , solutions will be stable. All the lines that are inside this region will tend give non oscillatory stable steady state solutions. Here in this graph ,0<α<5,-0.8<wc<4 ; In this graph, descretization number is shown on y-axis and x-axis. C0 + =.3 C0 - =.7, L=10,D=1,v=0
  • 31. alpha_c v/s w_c with observation of the values in stable region for which steady state will be stationary with β= 0.5. For α≤αc , solutions will be stable. All the lines that are inside this region will tend give non oscillatory stable steady state solutions. Here in this graph ,0<α<5,-0.8<wc<4 ; In this graph, descretization number is shown on y-axis and x-axis. C0 + =.3 C0 - =.7, L=10,D=1,v=0.5
  • 32. Numerical Solutions of Dynamics of actin filaments without considering Bipolar Filament • After doing stability analysis, I solved the actin dynamics equation numerically using first order upwind scheme with adaptive time stepping . Here is a snapshot of theoretical explanations behind it.
  • 34. First Order Upwind scheme for actin dynamics scheme
  • 35. For solutions of actin dynamics equation without bipolar filament c2(1:N,1)=0.3*ones(1,N).*(1+rand(1,N)); c1(1:N,1)=0.7*ones(1,N).*(1+rand(1,N)); L=10;a=0.6;b=2;
  • 36. For solutions of actin dynamics equation with consideration of bipolar filament D=1;L=10;a=0.6;b=2;w1=0.3;w2=0.7; c1(1:N,1)=0.7*ones(1,N).*(1+rand(1,N)); c2(1:N,1)= 0.3*ones(1,N).*(1+rand(1,N)); c3(1:N,1)=0.09*ones(1,N).*(1+rand(1,N));
  • 37. Results • It is sufficient to check for k=2*pi/L for stability analysis.(refer page22) • Upwind scheme verifies the results obtained in stability analysis. • Verification of fact that solution is unstable for D*dt/(dx)2 >0.5 for forward euler case. • Stability patterns obtained are in consensus with expected data as described in paper.
  • 38. Experiences • During the project , I got introduced to various techniques to solve Non linear Equations. I also got to learn about programming on MATLAB. • Back to project, so far I have done numerical analysis of solutions for actin dynamics equations. As we have discussed, I am currently working on stress calculations. I hope I was good during internship!!! It was my first research internship and I have learnt a lot from you. Thanks for your guidance!!!