Location: Signalling Principled Model Example @ ad404d66f734 / Calcium_Handling.cellml

Author:
Michael Cooling <mcoo001@UoA.auckland.ac.nz>
Date:
2016-04-11 12:30:41+12:00
Desc:
Adding .zip of complete model
Permanent Source URI:
https://models.physiomeproject.org/workspace/2c4/rawfile/ad404d66f734c5e175dca6c4601b7a2f00d8a1a3/Calcium_Handling.cellml

<?xml version="1.0" encoding="utf-8"?>

<model xmlns:xlink="http://www.w3.org/1999/xlink" name="Calcium_Handling" cmeta:id="Calcium_Handling" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Units_Cellular_Systems.cellml">
		<units name="per_s" units_ref="per_s"/>
		<units name="per_m" units_ref="per_m"/>
		<units name="per_m2" units_ref="per_m2"/>
		<units name="N_per_m2" units_ref="N_per_m2"/>
		<units name="m2_per_N" units_ref="m2_per_N"/>
		<units name="N_per_m" units_ref="N_per_m"/>
		<units name="N_m_per_K" units_ref="N_m_per_K"/>
		<units name="uM" units_ref="uM"/>
		<units name="uM_per_s" units_ref="uM_per_s"/>
		<units name="per_uM_per_s" units_ref="per_uM_per_s"/>
	</import>
	
	<group>
	<relationship_ref relationship="encapsulation" />
		<component_ref component="Calcium_Handling"> 
			<component_ref component="Ca_st" />
			<component_ref component="Ca_st_delta" />
			<component_ref component="Ca_if" />
			<component_ref component="Ca_if_delta" />
			<component_ref component="CaStore_Release" />
			<component_ref component="CaStore_Resequestration" />
			<component_ref component="Ca_Extrusions">
				<component_ref component="Ca_Extrusion_Basal" />
				<component_ref component="Ca_Extrusion_Pump" />
				<component_ref component="Ca_Extrusion_Exch" />
			</component_ref>
			<component_ref component="Ca_Intrusions" >
				<component_ref component="Ca_Intrusion_StressSensitiveChannels" />
				<component_ref component="Ca_Intrusion_CCE" />
			</component_ref>
		</component_ref>
	</group> 
	
 <component name="Calcium_Handling">
	<variable name="time" units="second" public_interface="in" private_interface="out"/>
	<variable name="Ca_st_init" units="uM" public_interface="in" private_interface="out"/>
	<variable name="V_ce" units="dimensionless" public_interface="in" private_interface="out"/> 

	<variable name="J_Ca_st_External" units="uM_per_s" public_interface="in" private_interface="out"/> 

	<variable name="J_Ca_if_External" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="Ca_if_init" units="uM" public_interface="in" private_interface="out"/>  

	<variable name="qrel_k2" units="per_s" public_interface="in" private_interface="out"/>
	<variable name="IP3" units="uM" public_interface="in" private_interface="out"/>
	<variable name="qrel_k1" units="uM" public_interface="in"  private_interface="out"/>

	<variable name="qres_k1" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qres_k2" units="uM" public_interface="in"  private_interface="out"/>
	<variable name="qres_k3" units="per_uM_per_s" public_interface="in"  private_interface="out"/>

	<variable name="Ca_if_rest" units="uM" public_interface="in" private_interface="out"/>
	<variable name="qout_basal_vmax" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qout_basal_Km" units="uM" public_interface="in" private_interface="out"/>
	<variable name="qout_basal_n" units="dimensionless" public_interface="in" private_interface="out"/>
	<variable name="qout_pump_vmax" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qout_pump_Km" units="uM" public_interface="in" private_interface="out"/>
	<variable name="qout_pump_n" units="dimensionless" public_interface="in" private_interface="out"/>
	<variable name="qout_exch_vmax" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qout_exch_Km" units="uM" public_interface="in" private_interface="out"/>  
	<variable name="qout_exch_n" units="dimensionless" public_interface="in" private_interface="out"/>

	<variable name="Ca_ex" units="uM" public_interface="in" private_interface="out"/>
	<variable name="l" units="meter" public_interface="in" private_interface="out"/>
	<variable name="T" units="kelvin" public_interface="in" private_interface="out"/>
	<variable name="tau" units="N_per_m2" public_interface="in" private_interface="out"/>
	<variable name="Ca_st_rest" units="uM" public_interface="in" private_interface="out"/>
	<variable name="qin_kcce" units="per_uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qin_max" units="uM_per_s" public_interface="in" private_interface="out"/>
	<variable name="qin_alpha" units="dimensionless" public_interface="in" private_interface="out"/>
	<variable name="qinw_f" units="dimensionless" public_interface="in" private_interface="out"/>
	<variable name="epsilon" units="dimensionless" public_interface="in" private_interface="out"/>
	<variable name="k" units="N_m_per_K" public_interface="in" private_interface="out"/>     
	<variable name="cac" units="per_m2" public_interface="in" private_interface="out"/>
	<variable name="qinw_delta" units="N_per_m" public_interface="in" private_interface="out"/>    

	<variable name="Ca_st" units="uM"  private_interface="in" public_interface="out"/>
	<variable name="J_Ca_st" units="uM_per_s" private_interface="in"  public_interface="out" /> 
	<variable name="Ca_if" units="uM" private_interface="in" public_interface="out" />
	<variable name="J_Ca_if" units="uM_per_s" private_interface="in"  public_interface="out" /> 
   </component>
	
   	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Template_Species_uM.cellml">
		<component name="Ca_st" component_ref="Template_Species_uM"/>
	</import>   
	
	<connection>
		<map_components component_1="Ca_st" component_2="Calcium_Handling"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="concentrationInitialValue" variable_2="Ca_st_init"/>
		<map_variables variable_1="concentration" variable_2="Ca_st"/>
	</connection> 
	
	<connection>
		<map_components component_1="Ca_st" component_2="Ca_st_delta"/>
			<map_variables variable_1="JGain" variable_2="JNet"/>
	</connection> 	
	
	<component name="Ca_st_delta">
	<variable name="JPlusCaStore_Resequestration" units="uM_per_s" public_interface="in" />
	<variable name="JMinusCaStore_Release" units="uM_per_s" public_interface="in" />
	<variable name="JExternal" units="uM_per_s" public_interface="in" />
	<variable name="V_ce" units="dimensionless" public_interface="in" /> 
	<variable name="JNet"  units="uM_per_s" public_interface="out" />
	<math xmlns="http://www.w3.org/1998/Math/MathML">
		<apply>
			<eq/>
			<ci>JNet</ci>
			<apply>
				<plus/>
				<apply>
					<times/>
					<ci>V_ce</ci>
					<apply>
						<minus/>
							<ci>JPlusCaStore_Resequestration</ci>
							<ci>JMinusCaStore_Release</ci>
					</apply>
				</apply>
				<ci>JExternal</ci>
			</apply>
		</apply>
	</math>
	</component>
	
	<connection>
		<map_components component_1="Ca_st_delta" component_2="Calcium_Handling"/>
			<map_variables variable_1="V_ce" variable_2="V_ce"/>
			<map_variables variable_1="JNet" variable_2="J_Ca_st"/>
			<map_variables variable_1="JExternal" variable_2="J_Ca_st_External"/>
	</connection> 	
	
	<connection>
		<map_components component_1="Ca_st_delta" component_2="CaStore_Release"/>
		<map_variables variable_1="JMinusCaStore_Release" variable_2="J_rel"/>
	</connection>
	
	<connection>
		<map_components component_1="Ca_st_delta" component_2="CaStore_Resequestration"/>
		<map_variables variable_1="JPlusCaStore_Resequestration" variable_2="J_res"/>
	</connection>
	
	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Template_Species_uM.cellml">
		<component name="Ca_if" component_ref="Template_Species_uM"/>
	</import>   
	
	<connection>
		<map_components component_1="Ca_if" component_2="Calcium_Handling"/>
		<map_variables variable_1="time" variable_2="time"/>
		<map_variables variable_1="concentrationInitialValue" variable_2="Ca_if_init"/>
		<map_variables variable_1="concentration" variable_2="Ca_if"/>
	</connection> 
	
	<connection>
		<map_components component_1="Ca_if" component_2="Ca_if_delta"/>
			<map_variables variable_1="JGain" variable_2="JNet"/>
	</connection> 	
	
	<component name="Ca_if_delta">
	<variable name="JPlusCaStore_Release" units="uM_per_s" public_interface="in" />
	<variable name="JMinusCaStore_Resequestration" units="uM_per_s" public_interface="in" />
	<variable name="JExternal" units="uM_per_s" public_interface="in" />
	<variable name="JMinusCa_if_Out" units="uM_per_s" public_interface="in" /> 
	<variable name="JPlusCa_if_In"  units="uM_per_s" public_interface="in" />
	<variable name="JNet"  units="uM_per_s" public_interface="out" />
	<math xmlns="http://www.w3.org/1998/Math/MathML">
	<apply>
	     <eq/>
	     <ci>JNet</ci>
	     <apply>
		<plus/>
		<apply>
		   <minus/>
		   <apply>
		      <plus/>
		      <apply>
			 <minus/>
			 <ci>JPlusCaStore_Release</ci>
			 <ci>JMinusCaStore_Resequestration</ci>
		      </apply>
		      <ci>JExternal</ci>
		   </apply>
		   <ci>JMinusCa_if_Out</ci>
		</apply>
		<ci>JPlusCa_if_In</ci>
	     </apply>
         </apply>	
	</math>
	</component>
	
	<connection>
		<map_components component_1="Ca_if_delta" component_2="Calcium_Handling"/>		
			<map_variables variable_1="JNet" variable_2="J_Ca_if"/>
			<map_variables variable_1="JExternal" variable_2="J_Ca_if_External"/>
	</connection> 
	
	<connection>
		<map_components component_1="Ca_if_delta" component_2="Ca_Intrusions"/>
			<map_variables variable_1="JPlusCa_if_In" variable_2="J_in"/>	
	</connection> 
	
	<connection>
	<map_components component_1="Ca_if_delta" component_2="Ca_Extrusions"/>
			<map_variables variable_1="JMinusCa_if_Out" variable_2="J_out"/>
	</connection>
	
	<connection>
		<map_components component_1="Ca_if_delta" component_2="CaStore_Release"/>
		<map_variables variable_1="JPlusCaStore_Release" variable_2="J_rel"/>
	</connection>
	
	<connection>
		<map_components component_1="Ca_if_delta" component_2="CaStore_Resequestration"/>
		<map_variables variable_1="JMinusCaStore_Resequestration" variable_2="J_res"/>
	</connection>
	
	<component name="CaStore_Release">
	<variable name="qrel_k2" units="per_s" public_interface="in"/>
	<variable name="IP3" units="uM" public_interface="in" />
	<variable name="qrel_k1" units="uM" public_interface="in"/>
	<variable name="Ca_st" units="uM" public_interface="in" />
	
	<variable name="J_rel" units="uM_per_s" public_interface="out"/>
	<math xmlns="http://www.w3.org/1998/Math/MathML">
	<apply>
	                  <eq/>
	                  <ci>J_rel</ci>
	                  <apply>
	                     <times/>
	                     <ci>qrel_k2</ci>
	                     <apply>
	                        <power/>
	                        <apply>
	                           <divide/>
	                           <ci>IP3</ci>
	                           <apply>
	                              <plus/>
	                              <ci>IP3</ci>
	                              <ci>qrel_k1</ci>
	                           </apply>
	                        </apply>
	                        <cn cellml:units="dimensionless">3</cn>
	                     </apply>
	                     <ci>Ca_st</ci>
	                  </apply>
	               </apply>           
         </math>	
	</component>
	
	<connection>
		<map_components component_1="CaStore_Release" component_2="Calcium_Handling"/>
			<map_variables variable_1="qrel_k2" variable_2="qrel_k2"/>
			<map_variables variable_1="IP3" variable_2="IP3"/>
			<map_variables variable_1="qrel_k1" variable_2="qrel_k1"/>			
	</connection> 
	
	<connection>
		<map_components component_1="CaStore_Release" component_2="Ca_st"/>
			<map_variables variable_1="Ca_st" variable_2="concentration"/>			
	</connection> 
	
	<component name="CaStore_Resequestration">
		<variable name="qres_k1" units="uM_per_s" public_interface="in"/>
		<variable name="Ca_if" units="uM" public_interface="in" />
		<variable name="qres_k2" units="uM" public_interface="in" />
		<variable name="qres_k3" units="per_uM_per_s" public_interface="in"/>
		<variable name="Ca_st" units="uM"  public_interface="in"/>
		<variable name="J_res" units="uM_per_s" public_interface="out"/>
	<math xmlns="http://www.w3.org/1998/Math/MathML">
	 <apply>
	  	<eq/>
	  	<ci>J_res</ci>
	  	<apply>
	    		 <minus/>
	     		<apply>
			<times/>
			<ci>qres_k1</ci>
				<apply>
					<power/>
					<apply>
						<divide/>
						<ci>Ca_if</ci>
						<apply>
							 <plus/>
							 <ci>Ca_if</ci>
							 <ci>qres_k2</ci>
						</apply>
					</apply>
				   <cn cellml:units="dimensionless">2</cn>
				</apply>
	     		</apply>
			     <apply>
				<times/>
				<ci>qres_k3</ci>
				<apply>
				   <power/>
				   <ci>Ca_st</ci>
				   <cn cellml:units="dimensionless">2</cn>
				</apply>
			     </apply>
	  	</apply>
         </apply>              
         </math>	
	</component>
	
	<connection>
		<map_components component_1="CaStore_Resequestration" component_2="Calcium_Handling"/>
			<map_variables variable_1="qres_k2" variable_2="qres_k2"/>
			<map_variables variable_1="qres_k1" variable_2="qres_k1"/>
			<map_variables variable_1="qres_k3" variable_2="qres_k3"/>			
	</connection> 
	
	<connection>
		<map_components component_1="CaStore_Resequestration" component_2="Ca_st"/>
			<map_variables variable_1="Ca_st" variable_2="concentration"/>			
	</connection> 
	
	<connection>
		<map_components component_1="CaStore_Resequestration" component_2="Ca_if"/>
			<map_variables variable_1="Ca_if" variable_2="concentration"/>			
	</connection> 
	
	<component name="Ca_Extrusions">
		<variable name="qout_basal_vmax" units="uM_per_s" public_interface="in" private_interface="out"/>
		<variable name="qout_basal_Km" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="qout_basal_n" units="dimensionless" public_interface="in"  private_interface="out"/>
		<variable name="qout_basal" units="uM_per_s" private_interface="in" />
		
		<variable name="qout_pump_vmax" units="uM_per_s" public_interface="in"  private_interface="out"/>
		<variable name="qout_pump_Km" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="qout_pump_n" units="dimensionless" public_interface="in"  private_interface="out"/>
		<variable name="qout_pump" units="uM_per_s" private_interface="in" />
		
		<variable name="qout_exch_vmax" units="uM_per_s" public_interface="in"  private_interface="out"/>
		<variable name="qout_exch_Km" units="uM" public_interface="in"  private_interface="out"/>  
		<variable name="qout_exch_n" units="dimensionless" public_interface="in"  private_interface="out"/>
		<variable name="qout_exch" units="uM_per_s" private_interface="in" />
		
		<variable name="Ca_if_rest" units="uM" public_interface="in"/>
		<variable name="Ca_if" units="uM" public_interface="in"  private_interface="out"/>
	
		<variable name="J_out" units="uM_per_s" public_interface="out"/>
		<math xmlns="http://www.w3.org/1998/Math/MathML">
		  <apply>
		            <eq/>
		            <ci>J_out</ci>
		            <piecewise>
		               <piece>
		                  <apply>
		                     <plus/>
		                     <ci>qout_basal</ci>
		                     <ci>qout_pump</ci>
		                     <ci>qout_exch</ci>
		                  </apply>
		                  <apply>
		                     <gt/>
		                     <ci>Ca_if</ci>
		                     <ci>Ca_if_rest</ci>
		                  </apply>
		               </piece>
		               <otherwise>
		                  <cn cellml:units="uM_per_s">0</cn>
		               </otherwise>
		            </piecewise>
        	 </apply>      
		
		</math>
	</component>
	
	<connection>
		<map_components component_1="Ca_Extrusions" component_2="Calcium_Handling"/>
			<map_variables variable_1="qout_basal_vmax" variable_2="qout_basal_vmax"/>
			<map_variables variable_1="qout_basal_Km"  variable_2="qout_basal_Km" />
			<map_variables variable_1="qout_basal_n" variable_2="qout_basal_n" />
			<map_variables variable_1="qout_pump_vmax" variable_2="qout_pump_vmax"/>
			<map_variables variable_1="qout_pump_Km" variable_2="qout_pump_Km" />
			<map_variables variable_1="qout_pump_n" variable_2="qout_pump_n" />
			<map_variables variable_1="qout_exch_vmax" variable_2="qout_exch_vmax"/>
			<map_variables variable_1="qout_exch_Km" variable_2="qout_exch_Km" />	
			<map_variables variable_1="qout_exch_n" variable_2="qout_exch_n" />	
			<map_variables variable_1="Ca_if_rest" variable_2="Ca_if_rest" />		
	</connection> 
	
	<connection>
		<map_components component_1="Ca_Extrusions" component_2="Ca_if"/>
			<map_variables variable_1="Ca_if" variable_2="concentration"/>			
	</connection> 
	
   	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Template_HillFunction_uM.cellml">
		<component name="Ca_Extrusion_Basal" component_ref="Template_HillFunction_uM"/>
	</import> 

	<connection>
		<map_components component_1="Ca_Extrusions" component_2="Ca_Extrusion_Basal"/>
			<map_variables variable_1="qout_basal_vmax" variable_2="Vmax"/>	
			<map_variables variable_1="qout_basal_Km" variable_2="Km"/>	
			<map_variables variable_1="qout_basal_n" variable_2="n"/>	
			<map_variables variable_1="qout_basal" variable_2="v"/>	
			<map_variables variable_1="Ca_if" variable_2="S"/>
	</connection>
	
   	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Template_HillFunction_uM.cellml">
		<component name="Ca_Extrusion_Pump" component_ref="Template_HillFunction_uM"/>
	</import> 

	<connection>
		<map_components component_1="Ca_Extrusions" component_2="Ca_Extrusion_Pump"/>
			<map_variables variable_1="qout_pump_vmax" variable_2="Vmax"/>	
			<map_variables variable_1="qout_pump_Km" variable_2="Km"/>	
			<map_variables variable_1="qout_pump_n" variable_2="n"/>	
			<map_variables variable_1="qout_pump" variable_2="v"/>	
			<map_variables variable_1="Ca_if" variable_2="S"/>
	</connection>
	
   	<import xmlns="http://www.cellml.org/cellml/1.1#" xlink:href="Principled_Base/Template_HillFunction_uM.cellml">
		<component name="Ca_Extrusion_Exch" component_ref="Template_HillFunction_uM"/>
	</import> 

	<connection>
		<map_components component_1="Ca_Extrusions" component_2="Ca_Extrusion_Exch"/>
			<map_variables variable_1="qout_exch_vmax" variable_2="Vmax"/>	
			<map_variables variable_1="qout_exch_Km" variable_2="Km"/>	
			<map_variables variable_1="qout_exch_n" variable_2="n"/>	
			<map_variables variable_1="qout_exch" variable_2="v"/>	
			<map_variables variable_1="Ca_if" variable_2="S"/>
	</connection>
	
	<component name="Ca_Intrusions">
		<variable name="epsilon" units="dimensionless" public_interface="in" private_interface="out"/>
		<variable name="tau" units="N_per_m2" public_interface="in"  private_interface="out"/>
		<variable name="l" units="meter" public_interface="in"  private_interface="out"/>
		<variable name="qinw_delta" units="N_per_m" public_interface="in"  private_interface="out"/>
		<variable name="qinw_f" units="dimensionless" public_interface="in"  private_interface="out"/>
		<variable name="k" units="N_m_per_K" public_interface="in"  private_interface="out"/> 
		<variable name="T" units="kelvin" public_interface="in"  private_interface="out"/>
		<variable name="cac" units="per_m2" public_interface="in"  private_interface="out"/>
		<variable name="qin_max" units="uM_per_s" public_interface="in"  private_interface="out"/>
		<variable name="qin_alpha" units="dimensionless" public_interface="in"  private_interface="out"/>
		<variable name="J_chan" units="uM_per_s" private_interface="in"/>
		
		<variable name="qin_kcce" units="per_uM_per_s" public_interface="in"  private_interface="out" />
		<variable name="Ca_st_rest" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="Ca_st" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="Ca_ex" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="Ca_if" units="uM" public_interface="in"  private_interface="out"/>
		<variable name="J_cce" units="uM_per_s" private_interface="in"/>
		
		<variable name="J_in" units="uM_per_s" public_interface="out" />
		
		 <math xmlns="http://www.w3.org/1998/Math/MathML">   
				  <apply>
				     <eq/>
				     <ci>J_in</ci>
				     <apply>
					<plus/>
					<ci>J_cce</ci>
					<ci>J_chan</ci>
				     </apply>
				 </apply> 
	      </math>
	</component>
	
	<connection>
		<map_components component_1="Ca_Intrusions" component_2="Calcium_Handling"/>
			<map_variables variable_1="epsilon" variable_2="epsilon"/>
			<map_variables variable_1="tau" variable_2="tau"/>
			<map_variables variable_1="l" variable_2="l"/>
			<map_variables variable_1="qinw_delta" variable_2="qinw_delta"/>
			<map_variables variable_1="qinw_f" variable_2="qinw_f"/>
			<map_variables variable_1="k" variable_2="k"/>
			<map_variables variable_1="T" variable_2="T"/>
			<map_variables variable_1="cac" variable_2="cac"/>
			<map_variables variable_1="qin_max" variable_2="qin_max"/>
			<map_variables variable_1="qin_alpha" variable_2="qin_alpha"/>
			<map_variables variable_1="qin_kcce" variable_2="qin_kcce"/>
			<map_variables variable_1="Ca_st_rest" variable_2="Ca_st_rest"/>
			
			<map_variables variable_1="Ca_ex" variable_2="Ca_ex"/>
	</connection> 
	<connection>
		<map_components component_1="Ca_Intrusions" component_2="Ca_if"/>
				<map_variables variable_1="Ca_if" variable_2="concentration"/>	
	</connection>
	<connection>
		<map_components component_1="Ca_Intrusions" component_2="Ca_st"/>
				<map_variables variable_1="Ca_st" variable_2="concentration"/>	
	</connection>	
	
	<component name="Ca_Intrusion_StressSensitiveChannels">
		<variable name="epsilon" units="dimensionless" public_interface="in"/>
		<variable name="tau" units="N_per_m2" public_interface="in" />
		<variable name="l" units="meter" public_interface="in" />
		<variable name="qinw_delta" units="N_per_m" public_interface="in" />
		<variable name="qinw_f" units="dimensionless" public_interface="in" />
		<variable name="k" units="N_m_per_K" public_interface="in" /> 
		<variable name="T" units="kelvin" public_interface="in" />
		<variable name="cac" units="per_m2" public_interface="in" />
		<variable name="qin_max" units="uM_per_s" public_interface="in" />
		<variable name="qin_alpha" units="dimensionless" public_interface="in" />
		<variable name="qin_w" units="dimensionless"/>
		<variable name="qinw_a" units="N_per_m"/>
		<variable name="J_chan" units="uM_per_s" public_interface="out"/>
		
	 	<math xmlns="http://www.w3.org/1998/Math/MathML">  
		       <apply>
			    <eq/>
			    <ci>qinw_a</ci>
			    <apply>
			       <plus/>
			       <apply>
				  <times/>
				  <ci>epsilon</ci>
				  <ci>tau</ci>
				  <ci>l</ci>
			       </apply>
			       <apply>
				  <root/>
				  <apply>
				     <plus/>
				     <apply>
					<times/>
					<cn cellml:units="dimensionless">16</cn>
					<apply>
					   <power/>
					   <ci>qinw_delta</ci>
					   <cn cellml:units="dimensionless">2</cn>
					</apply>
				     </apply>
				     <apply>
					<times/>
					<apply>
					   <power/>
					   <ci>epsilon</ci>
					   <cn cellml:units="dimensionless">2</cn>
					</apply>
					<apply>
					   <power/>
					   <ci>tau</ci>
					   <cn cellml:units="dimensionless">2</cn>
					</apply>
					<apply>
					   <power/>
					   <ci>l</ci>
					   <cn cellml:units="dimensionless">2</cn>
					</apply>
				     </apply>
				  </apply>
			       </apply>
			    </apply>
			 </apply>         

			 <apply>
			    <eq/>
			    <ci>qin_w</ci>
			    <apply>
			       <divide/>
			       <apply>
				  <times/>
				  <ci>qinw_f</ci>
				  <apply>
				     <power/>
				     <apply>
					<minus/>
					<ci>qinw_a</ci>
					<apply>
					   <times/>
					   <cn cellml:units="dimensionless">4</cn>
					   <ci>qinw_delta</ci>
					</apply>
				     </apply>
				     <cn cellml:units="dimensionless">2</cn>
				  </apply>
			       </apply>
			       <apply>
				  <times/>
				  <cn cellml:units="dimensionless">8</cn>
				  <ci>k</ci>
				  <ci>T</ci>
				  <ci>cac</ci>
				  <ci>qinw_a</ci>
			       </apply>
			    </apply>
			 </apply>

			 <apply>
			    <eq/>
			    <ci>J_chan</ci>
			    <apply>
			       <divide/>
			       <ci>qin_max</ci>
			       <apply>
				  <plus/>
				  <cn cellml:units="dimensionless">1</cn>
				  <apply>
				     <times/>
				     <ci>qin_alpha</ci>
				     <apply>
					<power/>
					<exponentiale/>
					<ci>qin_w</ci>
				     </apply>
				  </apply>
			       </apply>
			    </apply>
		 </apply>
	 	</math>
	</component>
	
	<connection>
		<map_components component_1="Ca_Intrusions" component_2="Ca_Intrusion_StressSensitiveChannels"/>
			<map_variables variable_1="epsilon" variable_2="epsilon"/>
			<map_variables variable_1="tau" variable_2="tau"/>
			<map_variables variable_1="l" variable_2="l"/>
			<map_variables variable_1="qinw_delta" variable_2="qinw_delta"/>
			<map_variables variable_1="qinw_f" variable_2="qinw_f"/>
			<map_variables variable_1="k" variable_2="k"/>
			<map_variables variable_1="T" variable_2="T"/>
			<map_variables variable_1="cac" variable_2="cac"/>
			<map_variables variable_1="qin_max" variable_2="qin_max"/>
			<map_variables variable_1="qin_alpha" variable_2="qin_alpha"/>
			<map_variables variable_1="J_chan" variable_2="J_chan"/>
	</connection> 
	
	<component name="Ca_Intrusion_CCE">
		<variable name="qin_kcce" units="per_uM_per_s" public_interface="in"  />
		<variable name="Ca_st_rest" units="uM" public_interface="in"  />
		<variable name="Ca_st" units="uM" public_interface="in" />
		<variable name="Ca_ex" units="uM" public_interface="in" />
		<variable name="Ca_if" units="uM" public_interface="in" />
		<variable name="J_cce" units="uM_per_s" public_interface="out"  />
	 	<math xmlns="http://www.w3.org/1998/Math/MathML"> 
			 <apply>
			    <eq/>
			    <ci>J_cce</ci>
			    <apply>
			       <times/>
			       <ci>qin_kcce</ci>
			       <apply>
				  <minus/>
				  <ci>Ca_st_rest</ci>
				  <ci>Ca_st</ci>
			       </apply>
			       <apply>
				  <minus/>
				  <ci>Ca_ex</ci>
				  <ci>Ca_if</ci>
			       </apply>
			    </apply>
			</apply>
	 	</math>
	</component>
	
	<connection>
		<map_components component_1="Ca_Intrusions" component_2="Ca_Intrusion_CCE"/>
			<map_variables variable_1="qin_kcce" variable_2="qin_kcce"/>
			<map_variables variable_1="Ca_st_rest" variable_2="Ca_st_rest"/>
			<map_variables variable_1="Ca_st" variable_2="Ca_st"/>
			<map_variables variable_1="Ca_ex" variable_2="Ca_ex"/>
			<map_variables variable_1="Ca_if" variable_2="Ca_if"/>
			<map_variables variable_1="J_cce" variable_2="J_cce"/>
	</connection>
	
</model>