Journal Club: How Anti-Arrhythmic Drugs Increase the Rate of Sudden Cardiac Death

We had the first meeting of our cardiac modelling journal club yesterday, and we looked at this paper by C. Frank Starmer from 2002.

The paper was investigating "Sanderson's paradox", whereby drugs intended to treat cardiac arrhythmia actually increased the incidence of sudden cardiac death in patients, even though they displayed what were thought to be anti-arrhythmic properties in single cells. (Incidentally, I think "paradox" is an overly dramatic word for "we discovered that we were wrong about something"). He was looking specifically at sodium blockers, because they were common anti-arrhythmic drugs.

He first describes a method for modelling drug block: it has been observed that sodium blockers are more effective when the heart is paced faster. This would suggest that the drug might bind only to the open ion channel. He formulated this as an extra gate on the sodium channel in addition to the h, m, and j gates.

He next talks about the refractory interval, which is the length of time a cell spends being unable to depolarise again after the initial stimulus. It was generally considered that increasing the refractory interval would be anti-arrhythmic, because it would make the cell less vulnerable to being stimulated by some random event.
So the thought was that you'd start with this:



And then your drug would change it to this:



Judging by the results of these drugs in humans, however, it looks like this isn't what happens on a larger scale. Starmer then looked at a group of cells arranged in a 1-dimensional cable, and talked about finding a different definition for the vulnerable period, based on the properties of the whole cable.

Usually when you apply an electrical stimulus (S1) (like you'd get during a heartbeat) to the end of a 1D cable, like so:




You get a wave propagating away from the stimulus, as expected:



This is all fine and normal, but what if you had a second stimulus (S2) outside of the usual time and place of stimulation?



You'd either get a wave propagating backwards along the cable:



or waves going in both directions:



In 1D, having an aberrant wave that propagates in one direction only is a pro-arrhythmic marker, because it could lead to a re-entrant arrhythmia in the whole heart. It's important, therefore, to work out what amount of drug block is going to cause unidirectional propagation and what amount bidirectional, and to redefine the "vulnerable period" as a property of the propagating wave. The vulnerable period depends on such things as the distance and time between S1 and S2, the speed of the wave, and the sodium channel availability.


If I have understood this paper correctly, the vulnerable region is the time between there being enough sodium channel availability at the location of S2 to create a backwards wave and there being enough availabiltity to create a forwards wave at the same location.

This was a really interesting thing to do, and I'm looking forward to having more journal club sessions, because it helped me to pick out the parts of the paper that I thought I understood but didn't, and brought together a lot of clever people in a room to explain them to me.

...

NC3Rs workshop

The main motivation behind my research is coming up with new and better methods for testing drugs and other substances to see if they are effective and safe. I've written previously about some of the promising animal testing alternatives that are being created.

Last week I went to a workshop run by the national centre for the replacement/refinement/reduction of animals in research (NC3Rs) and POEMS. Biologists brought 3Rs problems, and a group of computational and mathematical modellers tried to solve them.

There were five different projects presented. One was about using fruit flies to model diseases like Alzheimer's and Parkinson's by seeing how their courting behaviour is affected, another was modelling the way that nerves respond to bladder filling, one was related to immunoglobulin, and another was modelling retina damage in shaken baby syndrome (so much respect for the people who picked that project. It's really important but heartbreaking). The project I chose was about testing for neurodevelopmental safety.

The problem

If you affect the amount of thyroid hormones in the body of a pregnant woman, it can have huge effects on the developing foetus. A lot of potential drugs, potential industrial chemicals and potential pesticides turn out to significantly reduce the amount of thyroid hormones in the blood, meaning that they aren't safe to use around people.

Currently if you want to bring a new chemical to market, you have to test for thyroid effects using about 40 pregnant rats, each of which could be pregnant with up to 10 pups. We wanted to see if it was possible to use data from earlier studies on a few adult male rats to predict these effects, without having to use any extra animals. This method would still use data from animal tests, but would reduce the number of animals used, and could pave the way towards a non-animal method later on.

The picture below shows a simplified version of the pathway that thyroid hormones take through the body. T4 is produced in the mother's thyroid gland and then secreted into her blood. From here, it can either be removed by the liver, or carried into the blood of her foetus. It then travels from the foetus' blood to its developing brain cells. Inside the brain cells, T4 is converted into T3, which is essential for the brain to develop correctly (for a review, see Williams, 2008).

Diagram of thyroid hormone pathway
We concentrated on modelling the last two steps of this pathway: the transport of T4 from the outside of the cell to the inside, and its conversion to T3. The amount of T3 that ends up bound to its receptor in the nucleus is the most important variable for brain development. We wanted to see how the concentration of T4 in the foetus' blood affected the concentration of receptor-bound T3 in the nucleus.

The process

First, we all sat down together and came up with a set of equations to describe the way that all the components of the regulatory network interacted with each other. We used simple ordinary differential equations, which you might recognise from chemistry lessons.

For example, if a molecule A turns into a different molecule B at a rate of c molecules per second, and B turns into A at a rate of d molecules per second, you can say that


and

This way, given starting concentrations for A and B, you can work out how many molecules of each you have at a given time. We used this method to write down equations to describe the whole network.

After this, we each followed a different path of investigation. One researcher tried to derive analytically an equation to describe what steady state the system will reach, while others performed sensitivity analysis and parameter fitting. I was in a group that were looking at Petri nets, which are a different way of describing reaction networks, with some useful tools for analysis.

An animation of the petri net for this system is shown below. The circles (known as "places"), represent possible states that molecules can be in, the numbers inside the places represent the number of molecules in that state. The squares are transitions - usually representing some kind of reaction or binding event - and the red dots moving between places and transitions are the molecules as they change from place to place.

Petri nets seem like a really useful way of modelling biological systems for people without much experience of ODEs or programming, because they are a direct analogue for a system of ODEs with all the rate constants set to 1. The next stage in our analysis will be to parameterise this model, and use some graph theory tools to find out new things about the system. I haven't done graph theory since A-level, so this should be quite interesting!

...

Lessons from talking to teenagers

My new anatomically correct friend

I've had the opportunity twice recently to talk about my research to a group of 14-17 year old secondary school students as part of the Women in Computer Science schools outreach programme in my department. This has been strange and interesting. The last time I had to address such a group, I was the same age as them, and they were people I had known for years. The dynamic feels very different from the other side.

I felt like my first talk was more successful than my second. The first time I made a point of talking about the different types of computational modelling that went on in the department - aided, beautifully, by the lovely detachably-organned anatomical model pictured above. I took out her organs and passed them around the class, explaining that yes, people did model all of these things. Yes, even the naughty bits (seriously, I have a friend whose PhD project is modelling the C. elegans gonad). I think I planned this talk better and aimed it more clearly at a young audience, and I managed make the experience quite interactive with questions and answers.

They were a small group, and there was one student who was clearly more engaged and interested than the rest. He answered my every question eagerly and accurately. I hadn't realised before just how much of an impediment this can be if you're trying to talk to a group! I was always that kid at school, and it drove me crazy that the teacher would wait for someone else to respond when I knew the answer and I wanted to get to the next part of the lesson already.

The second time I gave a talk to school children (today!), there were about a hundred of 'em and I was concentrating more on describing my own research, because they'd heard all about computational medicine that morning. There were fewer jokes and I spent more time trying to get them to understand the material. I think I made it a bit too dry and technical, when it needed to be more of a fun overview of the subject. However, they totally laughed at all two of my jokes. I now understand how my dad managed to maintain his stock of terrible jokes over the course of his teaching career.

I was following a brilliantly charismatic speaker (this guy!) who really managed to get the kids' attention and keep them engaged. He really made sure to keep the content appropriate for the audience (social networking and selfies came up a lot!) and gave out practically every piece of information in the form of a question and answer, forcing everyone to stay alert. He made jokes all the way through and kept them laughing and happy. I'm more used to a straight lecture-and-questions format from years of university education, and I'd almost forgotten the give and take that happens with school classes.

I really enjoy talking to young people and I hope to do lots more schools outreach in the future. The next time I give a talk I'm going to make it a lot more fun and a lot more interactive!

...