E3 ubiquitin ligase rififylin has yin and yang effects on rabbit cardiac transient outward potassium currents (Ito) and corresponding channel proteins.

Genome-wide association studies have reported a correlation between a SNP of the RING finger E3 ubiquitin protein ligase rififylin (RFFL) and QT interval variability in humans (Newton-Cheh et al., 2009). Previously, we have shown that RFFL downregulates expression and function of the human-like ether-a-go-go-related gene potassium channel and corresponding rapidly activating delayed rectifier potassium current (IKr) in adult rabbit ventricular cardiomyocytes. Here, we report that RFFL also affects the transient outward current (Ito), but in a peculiar way. RFFL overexpression in adult rabbit ventricular cardiomyocytes significantly decreases the contribution of its fast component (Ito,f) from 35% to 21% and increases the contribution of its slow component (Ito,s) from 65% to 79%. Since Ito,f in rabbits is mainly conducted by Kv4.3, we investigated the effect of RFFL on Kv4.3 expressed in HEK293A cells. We found that RFFL overexpression reduced Kv4.3 expression and corresponding Ito,f in a RING domain-dependent manner in the presence or absence of its accessory subunit Kv channel-interacting protein 2. On the other hand, RFFL overexpression in Kv1.4-expressing HEK cells leads to an increase in both Kv1.4 expression level and Ito,s, similarly in a RING domain-dependent manner. Our physiologically detailed rabbit ventricular myocyte computational model shows that these yin and yang effects of RFFL overexpression on Ito,f, and Ito,s affect phase 1 of the action potential waveform and slightly decrease its duration in addition to suppressing IKr. Thus, RFFL modifies cardiac repolarization reserve via ubiquitination of multiple proteins that differently affect various potassium channels and cardiac action potential duration.

Role of ephaptic coupling in discordant alternans domain sizes and action potential propagation in the heart.

Discordant alternans, the spatially out-of-phase alternation of the durations of propagating action potentials in the heart, has been linked to the onset of fibrillation, a major cardiac rhythm disorder. The sizes of the regions, or domains, within which these alternations are synchronized are critical in this link. However, computer models employing standard gap junction-based coupling between cells have been unable to reproduce simultaneously the small domain sizes and rapid action potential propagation speeds seen in experiments. Here we use computational methods to show that rapid wave speeds and small domain sizes are possible when a more detailed model of intercellular coupling that accounts for so-called ephaptic effects is used. We provide evidence that the smaller domain sizes are possible, because different coupling strengths can exist on the wavefronts, for which both ephaptic and gap-junction coupling are involved, in contrast to the wavebacks, where only gap-junction coupling plays an active role. The differences in coupling strength are due to the high density of fast-inward (sodium) channels known to localize on the ends of cardiac cells, which are only active (and thus engage ephaptic coupling) during wavefront propagation. Thus, our results suggest that this distribution of fast-inward channels, as well as other factors responsible for the critical involvement of ephaptic coupling in wave propagation, including intercellular cleft spacing, play important roles in increasing the vulnerability of the heart to life-threatening tachyarrhythmias. Our results, combined with the absence of short-wavelength discordant alternans domains in standard gap-junction-dominated coupling models, also provide evidence that both gap-junction and ephaptic coupling are critical in wavefront propagation and waveback dynamics.

Scaling laws for two-dimensional dendritic crystal growth in a narrow channel.

We investigate analytically and computationally the dynamics of two-dimensional needle crystal growth from the melt in a narrow channel. Our analytical theory predicts that, in the low supersaturation limit, the growth velocity V decreases in time t as a power law V∼t^{-2/3}, which we validate by phase-field and dendritic-needle-network simulations. Simulations further reveal that, above a critical channel width Λ≈5l_{D}, where l_{D} is the diffusion length, needle crystals grow with a constant V

Ephaptic Coupling as a Resolution to the Paradox of Action Potential Wave Speed and Discordant Alternans Spatial Scales in the Heart.

Previous computer simulations have suggested that existing models of action potential wave propagation in the heart are not consistent with observed wave propagation behavior. Specifically, computer models cannot simultaneously reproduce the rapid wave speeds and small spatial scales of discordant alternans patterns measured experimentally in the same simulation. The discrepancy is important, because discordant alternans can be a key precursor to the development of abnormal and dangerous rapid rhythms in the heart. In this Letter, we show that this paradox can be resolved by allowing so-called ephaptic coupling to play a primary role in wave front propagation in place of conventional gap-junction coupling. With this modification, physiological wave speeds and small discordant alternans spatial scales both occur with gap-junction resistance values that are more in line with those observed in experiments. Our theory thus also provides support to the hypothesis that ephaptic coupling plays an important role in normal wave propagation.

Hierarchical structure formation by crystal growth-front instabilities during ice templating.

Directional solidification of aqueous solutions and slurries in a temperature gradient is widely used to produce cellular materials through a phase separation of solutes or suspended particles between growing ice lamellae. While this process has analogies to the directional solidification of metallurgical alloys, it forms very different hierarchical structures. The resulting honeycomb-like porosity of freeze-cast materials consists of regularly spaced, lamellar cell walls which frequently exhibit unilateral surface features of morphological complexity reminiscent of living forms, all of which are unknown in metallurgical structures. While the strong anisotropy of ice-crystal growth has been hypothesized to play a role in shaping those structures, the mechanism by which they form has remained elusive. By directionally freezing binary water mixtures containing small solutes obeying Fickian diffusion, and phase-field modeling of those experiments, we reveal how those structures form. We show that the flat side of lamellae forms because of slow faceted ice-crystal growth along the c-axis, while weakly anisotropic fast growth in other directions, including the basal plane, is responsible for the unilateral features. Diffusion-controlled morphological primary instabilities on the solid-liquid interface form a cellular structure on the atomically rough side of the lamellae, which template regularly spaced "ridges" while secondary instabilities of this structure are responsible for the more complex features. Collating the results, we obtain a scaling law for the lamellar spacing,  [Formula: see text] , where [Formula: see text] and [Formula: see text] are the local growth rate and temperature gradient, respectively.

Cell invasion during competitive growth of polycrystalline solidification patterns.

Spatially extended cellular and dendritic array structures forming during solidification processes such as casting, welding, or additive manufacturing are generally polycrystalline. Both the array structure within each grain and the larger scale grain structure determine the performance of many structural alloys. How those two structures coevolve during solidification remains poorly understood. By in situ observations of microgravity alloy solidification experiments onboard the International Space Station, we have discovered that individual cells from one grain can unexpectedly invade a nearby grain of different misorientation, either as a solitary cell or as rows of cells. This invasion process causes grains to interpenetrate each other and hence grain boundaries to adopt highly convoluted shapes. Those observations are reproduced by phase-field simulations further demonstrating that invasion occurs for a wide range of misorientations. Those results fundamentally change the traditional conceptualization of grains as distinct regions embedded in three-dimensional space.

Microstructural Pattern Formation during Far-from-Equilibrium Alloy Solidification.

We introduce a new phase-field formulation of rapid alloy solidification that quantitatively incorporates nonequilibrium effects at the solid-liquid interface over a very wide range of interface velocities. Simulations identify a new dynamical instability of dendrite tip growth driven by solute trapping at velocities approaching the absolute stability limit. They also reproduce the formation of the widely observed banded microstructures, revealing how this instability triggers transitions between dendritic and microsegregation-free solidification. Predicted band spacings agree quantitatively with observations in rapidly solidified Al-Cu thin films.

Topological control of liquid-metal-dealloyed structures.

The past few years have witnessed the rapid development of liquid metal dealloying to fabricate nano-/meso-scale porous and composite structures with ultra-high interfacial area for diverse materials applications. However, this method currently has two important limitations. First, it produces bicontinuous structures with high-genus topologies for a limited range of alloy compositions. Second, structures have a large ligament size due to substantial coarsening during dealloying at high temperature. Here we demonstrate computationally and experimentally that those limitations can be overcome by adding to the metallic melt an element that promotes high-genus topologies by limiting the leakage of the immiscible element during dealloying. We further interpret this finding by showing that bulk diffusive transport of the immiscible element in the liquid melt strongly influences the evolution of the solid fraction and topology of the structure during dealloying. The results shed light on fundamental differences in liquid metal and electrochemical dealloying and establish a new approach to produce liquid-metal-dealloyed structures with desired size and topologies.

Fracture toughness of bone at the microscale.

Bone's hierarchical arrangement of collagen and mineral generates a confluence of toughening mechanisms acting at every length scale from the molecular to the macroscopic level. Molecular defects, disease, and age alter bone structure at different levels and diminish its fracture resistance. However, the inability to isolate and quantify the influence of specific features hampers our understanding and the development of new therapies. Here, we combine in situ micromechanical testing, transmission electron microscopy and phase-field modelling to quantify intrinsic deformation and toughening at the fibrillar level and unveil the critical role of fibril orientation on crack deflection. At this level dry bone is highly anisotropic, with fracture energies ranging between 5 and 30 J/m2 depending on the direction of crack propagation. These values are lower than previously calculated for dehydrated samples from large-scale tests. However, they still suggest a significant amount of energy dissipation. This approach provides a new tool to uncouple and quantify, from the bottom up, the roles played by the structural features and constituents of bone on fracture and how can they be affected by different pathologies. The methodology can be extended to support the rational development of new structural composites.

The endosomal trafficking regulator LITAF controls the cardiac Nav1.5 channel via the ubiquitin ligase NEDD4-2.

The QT interval is a recording of cardiac electrical activity. Previous genome-wide association studies identified genetic variants that modify the QT interval upstream of LITAF (lipopolysaccharide-induced tumor necrosis factor-α factor), a protein encoding a regulator of endosomal trafficking. However, it was not clear how LITAF might impact cardiac excitation. We investigated the effect of LITAF on the voltage-gated sodium channel Nav1.5, which is critical for cardiac depolarization. We show that overexpressed LITAF resulted in a significant increase in the density of Nav1.5-generated voltage-gated sodium current INa and Nav1.5 surface protein levels in rabbit cardiomyocytes and in HEK cells stably expressing Nav1.5. Proximity ligation assays showed co-localization of endogenous LITAF and Nav1.5 in cardiomyocytes, whereas co-immunoprecipitations confirmed they are in the same complex when overexpressed in HEK cells. In vitro data suggest that LITAF interacts with the ubiquitin ligase NEDD4-2, a regulator of Nav1.5. LITAF overexpression down-regulated NEDD4-2 in cardiomyocytes and HEK cells. In HEK cells, LITAF increased ubiquitination and proteasomal degradation of co-expressed NEDD4-2 and significantly blunted the negative effect of NEDD4-2 on INa We conclude that LITAF controls cardiac excitability by promoting degradation of NEDD4-2, which is essential for removal of surface Nav1.5. LITAF-knockout zebrafish showed increased variation in and a nonsignificant 15% prolongation of action potential duration. Computer simulations using a rabbit-cardiomyocyte model demonstrated that changes in Ca2+ and Na+ homeostasis are responsible for the surprisingly modest action potential duration shortening. These computational data thus corroborate findings from several genome-wide association studies that associated LITAF with QT interval variation.

Oscillatory-nonoscillatory transitions for inclined cellular patterns in three-dimensional directional solidification.

The oscillatory behavior of cellular patterns produced by directional solidification of a transparent alloy under microgravity conditions was recently observed to depend on the misorientation of the main crystal axis with respect to the direction of the imposed thermal gradient [Pereda et al., Phys. Rev. E 95, 012803 (2017)2470-004510.1103/PhysRevE.95.012803]. To characterize the oscillatory-nonoscillatory transition resulting from the variations of the crystal misorientation, new experiments performed in DECLIC-DSI onboard the International Space Station and phase-field simulations are analyzed and combined in the present study. Experimental results are extracted from movies showing regions that extend on both sides of a boundary between two grains with respective misorientations of roughly 3 and 7 degrees. A set of tools are developed to analyze the experimental data and the same analysis is reproduced for the numerical data. A number of points are addressed in the simulations, like the effects of the system dimensions. The oscillatory state is found to be favored by the increase of the geometrical degrees of freedom. In bulk samples, a good agreement is found between the experimental and the numerical oscillatory-nonoscillatory threshold given by the ratio of the drift time to the oscillation period at the transition. The existence and the origin of bursts of localized groups of oscillating cells within a globally nonoscillatory pattern are characterized. A qualitative description of the physical mechanism that governs the oscillatory-nonoscillatory transition is provided.

Late INa Blocker GS967 Supresses Polymorphic Ventricular Tachycardia in a Transgenic Rabbit Model of Long QT Type 2.

BACKGROUND: Long QT syndrome has been associated with sudden cardiac death likely caused by early afterdepolarizations (EADs) and polymorphic ventricular tachycardias (PVTs). Suppressing the late sodium current (INaL) may counterbalance the reduced repolarization reserve in long QT syndrome and prevent EADs and PVTs. METHODS: We tested the effects of the selective INaL blocker GS967 on PVT induction in a transgenic rabbit model of long QT syndrome type 2 using intact heart optical mapping, cellular electrophysiology and confocal Ca2+ imaging, and computer modeling. RESULTS: GS967 reduced ventricular fibrillation induction under a rapid pacing protocol (n=7/14 hearts in control versus 1/14 hearts at 100 nmol/L) without altering action potential duration or restitution and dispersion. GS967 suppressed PVT incidences by reducing Ca2+-mediated EADs and focal activity during isoproterenol perfusion (at 30 nmol/L, n=7/12 and 100 nmol/L n=8/12 hearts without EADs and PVTs). Confocal Ca2+ imaging of long QT syndrome type 2 myocytes revealed that GS967 shortened Ca2+ transient duration via accelerating Na+/Ca2+ exchanger (INCX)-mediated Ca2+ efflux from cytosol, thereby reducing EADs. Computer modeling revealed that INaL potentiates EADs in the long QT syndrome type 2 setting through (1) providing additional depolarizing currents during action potential plateau phase, (2) increasing intracellular Na+ (Nai) that decreases the depolarizing INCX thereby suppressing the action potential plateau and delaying the activation of slowly activating delayed rectifier K+ channels (IKs), suggesting important roles of INaL in regulating Nai. CONCLUSIONS: Selective INaL blockade by GS967 prevents EADs and abolishes PVT in long QT syndrome type 2 rabbits by counterbalancing the reduced repolarization reserve and normalizing Nai. Graphic Abstract: A graphic abstract is available for this article.

LITAF (Lipopolysaccharide-Induced Tumor Necrosis Factor) Regulates Cardiac L-Type Calcium Channels by Modulating NEDD (Neural Precursor Cell Expressed Developmentally Downregulated Protein) 4-1 Ubiquitin Ligase.

BACKGROUND: The turnover of cardiac ion channels underlying action potential duration is regulated by ubiquitination. Genome-wide association studies of QT interval identified several single-nucleotide polymorphisms located in or near genes involved in protein ubiquitination. A genetic variant upstream of LITAF (lipopolysaccharide-induced tumor necrosis factor) gene prompted us to determine its role in modulating cardiac excitation. METHODS: Optical mapping was performed in zebrafish hearts to determine Ca2+ transients. Live-cell confocal calcium imaging was performed on adult rabbit cardiomyocytes to determine intracellular Ca2+handling. L-type calcium channel (LTCC) current (ICa,L) was measured using whole-cell recording. To study the effect of LITAF on Cav1.2 (L-type voltage-gated calcium channel 1.2) channel expression, surface biotinylation, and Westerns were performed. LITAF interactions were studied using coimmunoprecipitation and in situ proximity ligation assay. RESULTS: LITAF knockdown in zebrafish resulted in a robust increase in calcium transients. Overexpressed LITAF in 3-week-old rabbit cardiomyocytes resulted in a decrease in ICa,L and Cavα1c abundance, whereas LITAF knockdown increased ICa,L and Cavα1c protein. LITAF-overexpressing decreases calcium transients in adult rabbit cardiomyocytes, which was associated with lower Cavα1c levels. In tsA201 cells, overexpressed LITAF downregulated total and surface pools of Cavα1c via increased Cavα1c ubiquitination and its subsequent lysosomal degradation. We observed colocalization between LITAF and LTCC in tsA201 and cardiomyocytes. In tsA201, NEDD (neural precursor cell expressed developmentally downregulated protein) 4-1, but not its catalytically inactive form NEDD4-1-C867A, increased Cavα1c ubiquitination. Cavα1c ubiquitination was further increased by coexpressed LITAF and NEDD4-1 but not NEDD4-1-C867A. NEDD4-1 knockdown abolished the negative effect of LITAF on ICa,L and Cavα1c levels in 3-week-old rabbit cardiomyocytes. Computer simulations demonstrated that a decrease of ICa,L current associated with LITAF overexpression simultaneously shortened action potential duration and decreased calcium transients in rabbit cardiomyocytes. CONCLUSIONS: LITAF acts as an adaptor protein promoting NEDD4-1-mediated ubiquitination and subsequent degradation of LTCC, thereby controlling LTCC membrane levels and function and thus cardiac excitation.

Universality and Stability Phase Diagram of Two-Dimensional Brittle Fracture.

The two-dimensional oscillatory crack instability, experimentally observed in a class of brittle materials under strongly dynamic conditions, has been recently reproduced by a nonlinear phase-field fracture theory. Here, we highlight the universal character of this instability by showing that it is present in materials exhibiting widely different near crack tip elastic nonlinearity, and by demonstrating that the oscillations wavelength follows a universal master curve in terms of dissipation-related and nonlinear elastic intrinsic length scales. Moreover, we show that upon increasing the driving force for fracture, a high-velocity tip-splitting instability emerges, as experimentally demonstrated. The analysis culminates in a comprehensive stability phase diagram of two-dimensional brittle fracture, whose salient properties and topology are independent of the form of near tip nonlinearity.

NCX-Mediated Subcellular Ca2+ Dynamics Underlying Early Afterdepolarizations in LQT2 Cardiomyocytes.

Long QT syndrome type 2 (LQT2) is a congenital disease characterized by loss of function mutations in hERG potassium channels (IKr). LQT2 is associated with fatal ventricular arrhythmias promoted by triggered activity in the form of early afterdepolarizations (EADs). We previously demonstrated that intracellular Ca2+ handling is remodeled in LQT2 myocytes. Remodeling leads to aberrant late RyR-mediated Ca2+ releases that drive forward-mode Na+-Ca2+ exchanger (NCX) current and slow repolarization to promote reopening of L-type calcium channels and EADs. Forward-mode NCX was found to be enhanced despite the fact that these late releases do not significantly alter the whole-cell cytosolic calcium concentration during a vulnerable period of phase 2 of the action potential corresponding to the onset of EADs. Here, we use a multiscale ventricular myocyte model to explain this finding. We show that because the local NCX current is a saturating nonlinear function of the local submembrane calcium concentration, a larger number of smaller-amplitude discrete Ca2+ release events can produce a large increase in whole-cell forward-mode NCX current without increasing significantly the whole-cell cytosolic calcium concentration. Furthermore, we develop novel insights, to our knowledge, into how alterations of stochastic RyR activity at the single-channel level cause late aberrant Ca2+ release events. Experimental measurements in transgenic LTQ2 rabbits confirm the critical arrhythmogenic role of NCX and identify this current as a potential target for antiarrhythmic therapies in LQT2.

The Ca2+ transient as a feedback sensor controlling cardiomyocyte ionic conductances in mouse populations.

Conductances of ion channels and transporters controlling cardiac excitation may vary in a population of subjects with different cardiac gene expression patterns. However, the amount of variability and its origin are not quantitatively known. We propose a new conceptual approach to predict this variability that consists of finding combinations of conductances generating a normal intracellular Ca2+ transient without any constraint on the action potential. Furthermore, we validate experimentally its predictions using the Hybrid Mouse Diversity Panel, a model system of genetically diverse mouse strains that allows us to quantify inter-subject versus intra-subject variability. The method predicts that conductances of inward Ca2+ and outward K+ currents compensate each other to generate a normal Ca2+ transient in good quantitative agreement with current measurements in ventricular myocytes from hearts of different isogenic strains. Our results suggest that a feedback mechanism sensing the aggregate Ca2+ transient of the heart suffices to regulate ionic conductances.

Transient Outward K+ Current (Ito) Underlies the Right Ventricular Initiation of Polymorphic Ventricular Tachycardia in a Transgenic Rabbit Model of Long-QT Syndrome Type 1.

BACKGROUND: Sudden death in long-QT syndrome type 1 (LQT1), an inherited disease caused by loss-of-function mutations in KCNQ1, is triggered by early afterdepolarizations (EADs) that initiate polymorphic ventricular tachycardia (pVT). We investigated ionic mechanisms that underlie pVT in LQT1 using a transgenic rabbit model of LQT1. METHODS: Optical mapping, cellular patch clamping, and computer modeling were used to elucidate the mechanisms of EADs in transgenic LQT1 rabbits. RESULTS: The results showed that shorter action potential duration in the right ventricle (RV) was associated with focal activity during pVT initiation. RV cardiomyocytes demonstrated higher incidence of EADs under 50 nmol/L isoproterenol. Voltage-clamp studies revealed that the transient outward potassium current (Ito) magnitude was 28% greater in RV associated with KChiP2 but with no differences in terms of calcium-cycling kinetics and other sarcolemmal currents. Perfusing with the Ito blocker 4-aminopyridine changed the initial focal sites of pVT from the RV to the left ventricle, corroborating the role of Ito in pVT initiation. Computer modeling showed that EADs occur preferentially in the RV because of the larger conductance of the slow-inactivating component of Ito, which repolarizes the membrane potential sufficiently rapidly to allow reactivation of ICa,L before IKr has had sufficient time to activate. CONCLUSIONS: Ito heterogeneity creates both triggers and an arrhythmogenic substrate in LQT1. In the absence of IKs, Ito interactions with ICa,L and IKr promote EADs in the RV while prolonging action potential duration in the left ventricle. This heterogeneity of action potential enhances dispersion of refractoriness and facilitates conduction blocks that initiate pVTs.

A personalized, multiomics approach identifies genes involved in cardiac hypertrophy and heart failure.

A traditional approach to investigate the genetic basis of complex diseases is to identify genes with a global change in expression between diseased and healthy individuals. However, population heterogeneity may undermine the effort to uncover genes with significant but individual contribution to the spectrum of disease phenotypes within a population. Here we investigate individual changes of gene expression when inducing hypertrophy and heart failure in 100 + strains of genetically distinct mice from the Hybrid Mouse Diversity Panel (HMDP). We find that genes whose expression fold-change correlates in a statistically significant way with the severity of the disease are either up or down-regulated across strains, and therefore missed by a traditional population-wide analysis of differential gene expression. Furthermore, those "fold-change" genes are enriched in human cardiac disease genes and form a dense co-regulated module strongly interacting with the cardiac hypertrophic signaling network in the human interactome. We validate our approach by showing that the knockdown of Hes1, predicted as a strong candidate, induces a dramatic reduction of hypertrophy by 80-90% in neonatal rat ventricular myocytes. Our results demonstrate that individualized approaches are crucial to identify genes underlying complex diseases as well as to develop personalized therapies.

Stochastic initiation and termination of calcium-mediated triggered activity in cardiac myocytes.

Cardiac myocytes normally initiate action potentials in response to a current stimulus that depolarizes the membrane above an excitation threshold. Aberrant excitation can also occur due to spontaneous calcium (Ca2+) release (SCR) from intracellular stores after the end of a preceding action potential. SCR drives the Na+/Ca2+ exchange current inducing a "delayed afterdepolarization" that can in turn trigger an action potential if the excitation threshold is reached. This "triggered activity" is known to cause arrhythmias, but how it is initiated and terminated is not understood. Using computer simulations of a ventricular myocyte model, we show that initiation and termination are inherently random events. We determine the probability of those events from statistical measurements of the number of beats before initiation and before termination, respectively, which follow geometric distributions. Moreover, we elucidate the origin of randomness by a statistical analysis of SCR events, which do not follow a Poisson process observed in other eukaryotic cells. Due to synchronization of Ca2+ releases during the action potential upstroke, waiting times of SCR events after the upstroke are narrowly distributed, whereas SCR amplitudes follow a broad normal distribution with a width determined by fluctuations in the number of independent Ca2+ wave foci. This distribution enables us to compute the probabilities of initiation and termination of bursts of triggered activity that are maintained by a positive feedback between the action potential upstroke and SCR. Our results establish a theoretical framework for interpreting complex and varied manifestations of triggered activity relevant to cardiac arrhythmias.