Interpreting equatorial X-ray diffraction patterns from striated muscle with multiscale modeling.

Abstract

Small-angle X-ray diffraction is the tool of choice to obtain nm-scale structural information from striated muscle under near-physiological conditions. A critical barrier to progress has been the lack of generally applicable computational tools for modeling these X-ray diffraction patterns, limiting our ability to extract all the structural information potentially available. We addressed this need by extending the spatially explicit simulation platform MUSICO to predict equatorial fiber X-ray diffraction patterns from skeletal muscle. Electron densities for thin filaments with bound crossbridges were calculated using PDB structures of myosin S1, actin, tropomyosin, troponin and analogous nebulin structure derived from tropomyosin PDB. Thick filament backbone densities were constructed based on the Squire 1973 model incorporating “ribbon motifs” from a high-resolution cryo-EM structure of the Lethocerus thick filament. Titin and myosin S2 regions were constructed as cylindrical rings of density shells. Myosin heads in an ordered “parked state” were positioned in a shell close to the thick filament backbone with any disordered myosin heads distributed in a shell extending out to actin. X-ray diffraction patterns were calculated from sub-system motifs consisting of 1 thick filament and two thin filaments by evaluating the Fourier transform of this object at the reciprocal hexagonal lattice positions. Intensities were averaged over the 500 thick filaments and 1000 thin filaments used in simulation. Using this approach, we found that we could recapitulate 7 independent experimental intensities (out to the 4,0 reflection) from relaxed rat soleus and EDL muscle by varying only the proportion of heads in the parked state and a temperature factor type disorder term. Fits to contracting data could be achieved by adjusting only the temperature factor term and the number of force-producing crossbridges predicted by MUSICO to match the force levels observed experimentally.