FigureS9Heijman.m

function datatable = FigureS9Heijman
%% Figure S9 Heijman: Inducing proarrhythmic behavior through constant inward current 
%% inject in Heijman model. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%--- "Slow delayed rectifier current protects ventricular myocytes from
% arrhythmic dynamics across multiple species: a computational study" ---%

% By: Varshneya,Devenyi,Sobie 
% For questions, please contact Dr.Eric A Sobie -> eric.sobie@mssm.edu 
% or put in a pull request or open an issue on the github repository:
% https://github.com/meeravarshneya1234/IKs_stabilizes_APs.git. 

%--- Note:
% Results displayed in manuscript were run using MATLAB 2016a on a 64bit
% Intel Processor. For exact replication of figures it is best to use these
% settings.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%--------------------------------------------------------------------------
                            %% -- FigureS9Heijman.m -- %%
% Description: Injecting a constant inward current of 0.1 A/F while voltage
% is greater than -60 mV to induce proarrhythmic behavior in each of the
% models. Tetsing the percent change in APD from baseline.

% Outputs:
% --> X_Heijman - struct that outputs the APDs, time, voltage, and state variables 

%---: Functions required to run this script :---%
% mainHRdBA_Inject.m - runs injection simulation for Heijman only 
%--------------------------------------------------------------------------
%%
% settings
settings.bcl = 1000;% Interval bewteen stimuli,[ms]
settings.freq = 100; %number of beats to stimulate first EAD 
settings.storeLast = 1; % Determine how many beats to keep. 1 = last beat, 2 = last two beats 
settings.stimdur = 2;% Stimulus duration
settings.Istim = -36.7;%Stimulus amplitude for each model, see Table S1 for details
settings.ISO = 0; %concentration of ISO to add; 0 = no ISO 
flags.SS = 1; % 1 - run steady state conditions 0 - do not run steady state conditions 
settings.showProgress = 0;

% option to block PKA targets: no block = 1; 100% block = 0 
flags.ICaL = 1; flags.IKs = 1; flags.PLB = 1; flags.TnI = 1; flags.INa = 1;
flags.INaK = 1; flags.RyR = 1; flags.IKur = 1;

Injects = [0 0.1]; %amount of current to inject into cell [A/F]

disp('Heijman Model')
for ii = 1:length(Injects)
    
    settings.Inject = Injects(ii);
    [currents,State,Ti,APDs,settings]=mainHRdBA_Inject(settings,flags);
       
    datatable.APDs(:,ii) = APDs;
    datatable.times{:,ii} = Ti;
    datatable.V{:,ii} =  State(:,1);
    datatable.statevars{:,ii} =  State;
    datatable.currents{:,ii} = currents;       
   
    disp(['Current Inject = ' num2str(Injects(ii))])

end