The interaction energy of charged filaments in an electrolyte: Results for all filament spacings

Muscle active force-length curve explained by an electrophysical model of interfilament spacing

The active isometric force-length relation (FLR) of striated muscle sarcomeres is central to understanding and modeling muscle function. The mechanistic basis of the descending arm of the FLR is well explained by the decreasing thin:thick filament overlap that occurs at long sarcomere lengths. The mechanistic basis of the ascending arm of the FLR (the decrease in force that occurs at short sarcomere lengths), alternatively, has never been well explained. Because muscle is a constant-volume system, interfilament lattice distances must increase as sarcomere length shortens. This increase would decrease thin and thick-filament electrostatic interactions independently of thin:thick filament overlap. To examine this effect, we present here a fundamental, physics-based model of the sarcomere that includes filament molecular properties, calcium binding, sarcomere geometry including both thin:thick filament overlap and interfilament radial distance, and electrostatics. The model gives extremely good fits to existing FLR data from a large number of different muscles across their entire range of measured activity levels, with the optimized parameter values in all cases lying within anatomically and physically reasonable ranges. A local first-order sensitivity analysis (varying individual parameters while holding the values of all others constant) shows that model output is most sensitive to a subset of model parameters, most of which are related to sarcomere geometry, with model output being most sensitive to interfilament radial distance. This conclusion is supported by re-running the fits with only this parameter subset being allowed to vary, which increases fit errors only moderately. These results show that the model well reproduces existing experimental data, and indicate that changes in interfilament spacing play as central a role as changes in filament overlap in determining the FLR, particularly on its ascending arm.

An electrostatic model with weak actin-myosin attachment resolves problems with the lattice stability of skeletal muscle

The stability of the filament lattice in relaxed striated muscle can be viewed as a balance of electrostatic and van der Waals forces. The simplest electrostatic model, where actin and myosin filaments are treated as charged cylinders, generates reasonable lattice spacings for skinned fibers. However, this model predicts excessive radial stiffness under osmotic pressure and cannot account for the initial pressure (∼1 kPa) required for significant compression. Good agreement with frog compression data is obtained with an extended model, in which S1 heads are weakly attached to actin when the lattice spacing is reduced below a critical value; further compression moves fixed negative charges on the heads closer to the myofilament backbone as they attach at a more acute angle to actin. The model predicts pH data in which the lattice shrinks as pH is lowered and protons bind to filaments. Electrostatic screening implies that the lattice shrinks with increasing ionic strength, but the observed expansion of the frog lattice at ionic strengths above 0.1 M with KCl might be explained if Cl(-) binds to sites on the motor domain of S1. With myosin-myosin and actin-actin interactions, the predicted lattice spacing decreases slightly with sarcomere length, with a more rapid decrease when actin-myosin filament overlap is very small.