Molecular and Subcellular-Scale Modeling of Nucleotide Diffusion in the Cardiac Myofilament Lattice

AbstractContractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions

Long QT Syndrome Type 2: Emerging Strategies for Correcting Class 2 \u3cem\u3eKCNH2 (hERG)\u3c/em\u3e Mutations and Identifying New Patients

Significant advances in our understanding of the molecular mechanisms that cause congenital long QT syndrome (LQTS) have been made. A wide variety of experimental approaches, including heterologous expression of mutant ion channel proteins and the use of inducible pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) from LQTS patients offer insights into etiology and new therapeutic strategies. This review briefly discusses the major molecular mechanisms underlying LQTS type 2 (LQT2), which is caused by loss-of-function (LOF) mutations in the KCNH2 gene (also known as the human ether-à-go-go-related gene or hERG). Almost half of suspected LQT2-causing mutations are missense mutations, and functional studies suggest that about 90% of these mutations disrupt the intracellular transport, or trafficking, of the KCNH2-encoded Kv11.1 channel protein to the cell surface membrane. In this review, we discuss emerging strategies that improve the trafficking and functional expression of trafficking-deficient LQT2 Kv11.1 channel proteins to the cell surface membrane and how new insights into the structure of the Kv11.1 channel protein will lead to computational approaches that identify which KCNH2 missense variants confer a high-risk for LQT2

Myofilament Calcium Sensitivity: Consequences of the Effective Concentration of Troponin I

Control of calcium binding to and dissociation from cardiac troponin C (TnC) is essential to healthy cardiac muscle contraction/relaxation. There are numerous aberrant post-translational modifications and mutations within a plethora of contractile, and even non-contractile, proteins that appear to imbalance this delicate relationship. The direction and extent of the resulting change in calcium sensitivity is thought to drive the heart toward one type of disease or another. There are a number of molecular mechanisms that may be responsible for the altered calcium binding properties of TnC, potentially the most significant being the ability of the regulatory domain of TnC to bind the switch peptide region of TnI. Considering TnI is essentially tethered to TnC and cannot diffuse away in the absence of calcium, we suggest that the apparent calcium binding properties of TnC are highly dependent upon an “effective concentration” of TnI available to bind TnC. Based on our previous work, TnI peptide binding studies and the calcium binding properties of chimeric TnC-TnI fusion constructs, and building upon the concept of effective concentration, we have developed a mathematical model that can simulate the steady-state and kinetic calcium binding properties of a wide assortment of disease-related and post-translational protein modifications in the isolated troponin complex and reconstituted thin filament. We predict that several TnI and TnT modifications do not alter any of the intrinsic calcium or TnI binding constants of TnC, but rather alter the ability of TnC to “find” TnI in the presence of calcium. These studies demonstrate the apparent consequences of the effective TnI concentration in modulating the calcium binding properties of TnC

Predicting the influence of long-range molecular interactions on macroscopic-scale diffusion by homogenization of the Smoluchowski equation

The macroscopic diffusion constant for a charged diffuser is in part dependent on (1) the volume excluded by solute “obstacles” and (2) long-range interactions between those obstacles and the diffuser. Increasing excluded volume reduces transport of the diffuser, while long-range interactions can either increase or decrease diffusivity, depending on the nature of the potential. We previously demonstrated [P. M. Kekenes-Huskey , Biophys. J. 105, 2130 (2013)] using homogenization theory that the configuration of molecular-scale obstacles can both hinder diffusion and induce diffusional anisotropy for small ions. As the density of molecular obstacles increases, van der Waals (vdW) and electrostatic interactions between obstacle and a diffuser become significant and can strongly influence the latter's diffusivity, which was neglected in our original model. Here, we extend this methodology to include a fixed (time-independent) potential of mean force, through homogenization of the Smoluchowski equation. We consider the diffusion of ions in crowded, hydrophilic environments at physiological ionic strengths and find that electrostatic and vdW interactions can enhance or depress effective diffusion rates for attractive or repulsive forces, respectively. Additionally, we show that the observed diffusion rate may be reduced independent of non-specific electrostatic and vdW interactions by treating obstacles that exhibit specific binding interactions as “buffers” that absorb free diffusers. Finally, we demonstrate that effective diffusion rates are sensitive to distribution of surface charge on a globular protein, Troponin C, suggesting that the use of molecular structures with atomistic-scale resolution can account for electrostatic influences on substrate transport. This approach offers new insight into the influence of molecular-scale, long-range interactions on transport of charged species, particularly for diffusion-influenced signaling events occurring in crowded cellular environments

Understanding Ion Binding Affinity and Selectivity in β‑Parvalbumin Using Molecular Dynamics and Mean Spherical Approximation Theory

Parvalbumin (PV) is a globular calcium (Ca²⁺)-selective protein expressed in a variety of biological tissues. Our computational studies of the rat β-parvalbumin (β-PV) isoform seek to elucidate the molecular thermodynamics of Ca²⁺ versus magnesium (Mg²⁺) binding at the protein’s two EF-hand motifs. Specifically, we have utilized molecular dynamics (MD) simulations and a mean-field electrolyte model (mean spherical approximation (MSA) theory) to delineate how the EF-hand scaffold controls the “local” thermodynamics of Ca²⁺ binding selectivity over Mg²⁺. Our MD simulations provide the probability density of metal-chelating oxygens within the EF-hand scaffolds for both Ca²⁺ and Mg²⁺, as well the conformational strain induced by Mg²⁺ relative to Ca²⁺ binding. MSA theory utilizes the binding domain oxygen and charge distributions to predict the chemical potential of ion binding, as well as their corresponding concentrations within the binding domain. We find that the electrostatic and steric contributions toward ion binding were similar for Mg²⁺ and Ca²⁺, yet the latter was 5.5 kcal/mol lower in enthalpy when internal strain within the EF hand was considered. We therefore speculate that beyond differences in dehydration energies for the Ca²⁺ versus Mg²⁺, strain induced in the β-PV EF hand by cation binding significantly contributes to the nearly 10,000-fold difference in binding affinity reported in the literature. We further complemented our analyses of local factors governing cation binding selectivity with whole-protein (global) contributions, such as interhelical residue–residue contacts and solvent exposure of hydrophobic surface. These contributions were found to be comparable for both Ca²⁺- and Mg²⁺-bound β-PV, which may implicate local factors, EF-hand strain, and dehydration, in providing the primary means of selectivity. We anticipate these methods could be used to estimate metal binding thermodynamics across a broad range of PV sequence homologues and EF-hand-containing, Ca²⁺ binding proteins