Enzyme Localization, Crowding, and Buffers Collectively Modulate Diffusion-Influenced Signal Transduction: Insights from Continuum Diffusion Modeling

Biochemical reaction networks consisting of coupled enzymes connect substrate signaling events with biological function. Substrates involved in these reactions can be strongly influenced by diffusion “barriers” arising from impenetrable cellular structures and macromolecules, as well as interactions with biomolecules, especially within crowded environments. For diffusion-influenced reactions, the spatial organization of diffusion barriers arising from intracellular structures, non-specific crowders, and specific-binders (buffers) strongly controls the temporal and spatial reaction kinetics. In this study, we use two prototypical biochemical reactions, a Goodwin oscillator, and a reaction with a periodic source/sink term to examine how a diffusion barrier that partitions substrates controls reaction behavior. Namely, we examine how conditions representative of a densely packed cytosol, including reduced accessible volume fraction, non-specific interactions, and buffers, impede diffusion over nanometer length-scales. We find that diffusion barriers can modulate the frequencies and amplitudes of coupled diffusion-influenced reaction networks, as well as give rise to “compartments” of decoupled reactant populations. These effects appear to be intensified in the presence of buffers localized to the diffusion barrier. These findings have strong implications for the role of the cellular environment in tuning the dynamics of signaling pathways

Dynamics and calcium association to the N-terminal regulatory domain of human cardiac troponin C: a multiscale computational study.

Troponin C (TnC) is an important regulatory molecule in cardiomyocytes. Calcium binding to site II in TnC initiates a series of molecular events that result in muscle contraction. The most direct change upon Ca(2+) binding is an opening motion of the molecule that exposes a hydrophobic patch on the surface allowing for Troponin I to bind. Molecular dynamics simulations were used to elucidate the dynamics of this crucial protein in three different states: apo, Ca(2+)-bound, and Ca(2+)-TnI-bound. Dynamics between the states are compared, and the Ca(2+)-bound system is investigated for opening motions. On the basis of the simulations, NMR chemical shifts and order parameters are calculated and compared with experimental observables. Agreement indicates that the simulations sample the relevant dynamics of the system. Brownian dynamics simulations are used to investigate the calcium association of TnC. We find that calcium binding gives rise to correlative motions involving the EF hand and collective motions conducive of formation of the TnI-binding interface. We furthermore indicate the essential role of electrostatic steering in facilitating diffusion-limited binding of Ca(2+)

Fidelity of Phenylalanyl-tRNA Synthetase in Binding the Natural Amino Acids

Aminoacyl-tRNA synthetases guard the fidelity of cognate amino acid incorporation during protein biosynthesis; for example, phenylalanyl-tRNA synthetase (PheRS) activates and transfers only Phe to its tRNA. Since we are interested in using a computational protocol to identify nonnatural amino acids that are incorporated by wild-type PheRS, it is critical to understand the fidelity of PheRS in binding the 20 natural amino acids. To this end, HierDock, a computational protocol for predicting binding sites and relative binding affinities, was used for testing the natural amino acids in PheRS. Scanning the entire ligand-accessible protein surface for the best binding region, we find that HierDock correctly identifies the active site of Phe in PheRS and predicts Phe within 0.61 Å RMSD of the crystal structure. HierDock also successfully shows PheRS discriminates for Phe, as the noncognate amino acids bind less favorably in the binding site of Phe. However, we find that Met, Cys, and Tyr bind competitively but at positions distant from the Phe binding site. This result corroborates in vitro measurements of aminoacyl adenylate formation, which show Met competes with Phe at the amino acid binding stage. We predict that the binding site of Met would not activate PheRS, as the noncognate amino acid cannot establish suitable hydrogen bonds with the PheRS reaction center. These results validate the use of HierDock in predicting the binding sites of the cognate amino acids in PheRS. The HierDock procedure calculates the discrimination of aminoacyl-tRNA synthetases at the stage of binding the cognate amino acid and offers a molecular level understanding of the mistakes made in protein biosynthesis that are not readily uncovered through experiments. This technique is also useful for predicting the binding of a selected nonnatural amino acid analogue, thereby indicating whether the molecule would be incorporated into a wild-type aminoacyl-tRNA synthetase

Increasing Salt Rejection of Polybenzimidazole Nanofiltration Membranes via the Addition of Immobilized and Aligned Aquaporins

Aquaporins are water channel proteins in cell membrane, highly specific for water molecules while restricting the passage of contaminants and small molecules, such as urea and boric acid. Cysteine functional groups were installed on aquaporin Z for covalent attachment to the polymer membrane matrix so that the proteins could be immobilized to the membranes and aligned in the direction of the flow. Depth profiling using x-ray photoelectron spectrometer (XPS) analysis showed the presence of functional groups corresponding to aquaporin Z modified with cysteine (Aqp-SH). Aqp-SH modified membranes showed a higher salt rejection as compared to unmodified membranes. For 2 M NaCl and CaCl2 solutions, the rejection obtained from Aqp-SH membranes was 49.3 ± 7.5% and 59.1 ± 5.1%. On the other hand, the rejections obtained for 2 M NaCl and CaCl2 solutions from unmodified membranes were 0.8 ± 0.4% and 1.3 ± 0.2% respectively. Furthermore, Aqp-SH membranes did not show a significant decrease in salt rejection with increasing feed concentrations, as was observed with other membranes. Through simulation studies, it was determined that there was approximately 24% capping of membrane pores by dispersed aquaporins

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

Modeling effects of L-type ca(2+) current and na(+)-ca(2+) exchanger on ca(2+) trigger flux in rabbit myocytes with realistic T-tubule geometries.

The transverse tubular system of rabbit ventricular myocytes consists of cell membrane invaginations (t-tubules) that are essential for efficient cardiac excitation-contraction coupling. In this study, we investigate how t-tubule micro-anatomy, L-type Ca(2+) channel (LCC) clustering, and allosteric activation of Na(+)/Ca(2+) exchanger by L-type Ca(2+) current affects intracellular Ca(2+) dynamics. Our model includes a realistic 3D geometry of a single t-tubule and its surrounding half-sarcomeres for rabbit ventricular myocytes. The effects of spatially distributed membrane ion-transporters (LCC, Na(+)/Ca(2+) exchanger, sarcolemmal Ca(2+) pump, and sarcolemmal Ca(2+) leak), and stationary and mobile Ca(2+) buffers (troponin C, ATP, calmodulin, and Fluo-3) are also considered. We used a coupled reaction-diffusion system to describe the spatio-temporal concentration profiles of free and buffered intracellular Ca(2+). We obtained parameters from voltage-clamp protocols of L-type Ca(2+) current and line-scan recordings of Ca(2+) concentration profiles in rabbit cells, in which the sarcoplasmic reticulum is disabled. Our model results agree with experimental measurements of global Ca(2+) transient in myocytes loaded with 50 μM Fluo-3. We found that local Ca(2+) concentrations within the cytosol and sub-sarcolemma, as well as the local trigger fluxes of Ca(2+) crossing the cell membrane, are sensitive to details of t-tubule micro-structure and membrane Ca(2+) flux distribution. The model additionally predicts that local Ca(2+) trigger fluxes are at least threefold to eightfold higher than the whole-cell Ca(2+) trigger flux. We found also that the activation of allosteric Ca(2+)-binding sites on the Na(+)/Ca(2+) exchanger could provide a mechanism for regulating global and local Ca(2+) trigger fluxes in vivo. Our studies indicate that improved structural and functional models could improve our understanding of the contributions of L-type and Na(+)/Ca(2+) exchanger fluxes to intracellular Ca(2+) dynamics

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

3D dSTORM imaging reveals novel detail of ryanodine receptor localization in rat cardiac myocytes

Cardiomyocyte contraction is dependent on Ca2+ release from ryanodine receptors (RyRs). However, the precise localization of RyRs remains unknown, due to shortcomings of imaging techniques which are diffraction limited or restricted to 2D. We aimed to determine the 3D nanoscale organization of RyRs in rat cardiomyocytes by employing direct stochastic optical reconstruction microscopy (dSTORM) with phase ramp technology. Initial observations at the cell surface showed an undulating organization of RyR clusters, resulting in their frequent overlap in the z‐axis and obscured detection by 2D techniques. Non‐overlapping clusters were imaged to create a calibration curve for estimating RyR number based on recorded fluorescence blinks. Employing this method at the cell surface and interior revealed smaller RyR clusters than 2D estimates, as erroneous merging of axially aligned RyRs was circumvented. Functional groupings of RyR clusters (Ca2+ release units, CRUs), contained an average of 18 and 23 RyRs at the surface and interior, respectively, although half of all CRUs contained only a single ‘rogue’ RyR. Internal CRUs were more tightly packed along z‐lines than surface CRUs, contained larger and more numerous RyR clusters, and constituted ∼75% of the roughly 1 million RyRs present in an average cardiomyocyte. This complex internal 3D geometry was underscored by correlative imaging of RyRs and t‐tubules, which enabled quantification of dyadic and non‐dyadic RyR populations. Mirroring differences in CRU size and complexity, Ca2+ sparks originating from internal CRUs were of longer duration than those at the surface. These data provide novel, nanoscale insight into RyR organization and function across cardiomyocytes