In mammalian skeletal muscle, the propagation of surface membrane depolarization into the interior of the muscle fibre along the transverse (T) tubular network is essential for the synchronized release of calcium from the sarcoplasmic reticulum (SR) via ryanodine receptors (RyRs) in response to the conformational change in the voltage-sensor dihydropyridine receptors. Deficiency in 3-phosphoinositide phosphatase myotubularin (MTM1) has been reported to disrupt T-tubules, resulting in impaired SR calcium release. Here confocal calcium transients recorded in muscle fibres of MTM1-deficient mice were compared with the results from a model where propagation of the depolarization along the T-tubules was modelled mathematically with disruptions in the network assumed to modify the access and transmembrane resistance as well as the capacitance. If, in simulations, T-tubules were assumed to be partially or completely inaccessible to the depolarization and RyRs at these points to be prime for calcium-induced calcium release, all the features of measured SR calcium release could be reproduced. We conclude that the inappropriate propagation of the depolarization into the fibre interior is the initial critical cause of severely impaired SR calcium release in MTM1 deficiency, while the Ca2+ -triggered opening of RyRs provides an alleviating support to the diseased process. KEY POINTS: Myotubular myopathy is a fatal disease due to genetic deficiency in the phosphoinositide phosphatase MTM1. Although the causes are known and corresponding gene therapy strategies are being developed, there is no mechanistic understanding of the disease-associated muscle function failure. Resolving this issue is of primary interest not only for a fundamental understanding of how MTM1 is critical for healthy muscle function, but also for establishing the related cellular mechanisms most primarily or stringently affected by the disease, which are thus of potential interest as therapy targets. The mathematical modelling approach used in the present work proves that the disease-associated alteration of the plasma membrane invagination network is sufficient to explain the dysfunctions of excitation-contraction coupling, providing the first integrated quantitative framework that explains the associated contraction failure.