Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange-correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew-Burke-Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew-Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ∼4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. These results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

CoeffNet: predicting activation barriers through a chemically-interpretable, equivariant and physically constrained graph neural network

Activation barriers of elementary reactions are essential to predict molecular reaction mechanisms and kinetics. However, computing these energy barriers by identifying transition states with electronic structure methods (e.g., density functional theory) can be time-consuming and computationally expensive. In this work, we introduce CoeffNet, an equivariant graph neural network that predicts activation barriers using coefficients of any frontier molecular orbital (such as the highest occupied molecular orbital) of reactant and product complexes as graph node features. We show that using coefficients as features offer several advantages, such as chemical interpretability and physical constraints on the network's behaviour and numerical range. Model outputs are either activation barriers or coefficients of the chosen molecular orbital of the transition state; the latter quantity allows us to interpret the results of the neural network through chemical intuition. We test CoeffNet on a dataset of SN2 reactions as a proof-of-concept and show that the activation barriers are predicted with a mean absolute error of less than 0.025 eV. The highest occupied molecular orbital of the transition state is visualized and the distribution of the orbital densities of the transition states is described for a few prototype SN2 reactions.

Ryanodine Receptor Stabilization Therapy Suppresses Ca 2+ -Based Arrhythmias in a Novel Model of Metabolic HFpEF

Heart Failure with preserved ejection fraction (HFpEF) is the most prevalent form of heart failure worldwide and its significant mortality is associated with a high rate of sudden cardiac death (SCD; 30% - 40%). Chronic metabolic stress is an important driver of HFpEF, and clinical data show metabolic stress as a significant risk factor for ventricular arrhythmias in HFpEF patients. The mechanisms of SCD and ventricular arrhythmia in HFpEF remain critically understudied and empirical treatment is ineffective. To address this important knowledge gap, we developed a novel preclinical model of metabolic-stress induced HFpEF using Western diet (High fructose and fat) and hypertension induced by nitric oxide synthase inhibition (with L-NAME) in wildtype C57BL6/J mice. After 5 months, mice display all clinical characteristics of HFpEF and present with stress-induced sustained ventricular tachycardia (VT). Mechanistically, we found a novel pattern of arrhythmogenic intracellular Ca 2+ handling that is distinct from the well-characterized changes pathognomonic for heart failure with reduced ejection fraction. In addition, we show that the transverse tubular system remains intact in HFpEF and that arrhythmogenic, intracellular Ca 2+ mobilization becomes hyper-sensitive to ß- adrenergic activation. Finally, in proof-of-concept experiments we show in vivo that the clinically used intracellular calcium stabilizer dantrolene, which acts on the Ca 2+ release channels of the sarcoplasmic reticulum (SR), the ryanodine receptors, acutely prevents stress-induced VT in HFpEF mice. Therapeutic control of SR Ca 2+ leak may present a novel mechanistic treatment approach in metabolic HFpEF.

Unconventional Error Cancellation Explains the Success of Hartree-Fock Density Functional Theory for Barrier Heights

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree-Fock electron densities, which has been known for ∼30 years. The conventional explanation is that the Hartree-Fock theory yields the more accurate density. This work presents a benchmark Kohn-Sham inversion of accurate coupled-cluster densities for the reaction H2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree-Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532-543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.

Testing the r2SCAN Density Functional for the Thermodynamic Stability of Solids with and without a van der Waals Correction

A central aim of materials discovery is an accurate and numerically reliable description of thermodynamic properties, such as the enthalpies of formation and decomposition. The r2SCAN revision of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (meta-GGA) balances numerical stability with high general accuracy. To assess the r2SCAN description of solid-state thermodynamics, we evaluate the formation and decomposition enthalpies, equilibrium volumes, and fundamental band gaps of more than 1000 solids using r2SCAN, SCAN, and PBE, as well as two dispersion-corrected variants, SCAN+rVV10 and r2SCAN+rVV10. We show that r2SCAN achieves accuracy comparable to SCAN and often improves upon SCAN's already excellent accuracy. Although SCAN+rVV10 is often observed to worsen the formation enthalpies of SCAN and makes no substantial correction to SCAN's cell volume predictions, r2SCAN+rVV10 predicts marginally less accurate formation enthalpies than r2SCAN, and slightly more accurate cell volumes than r2SCAN. The average absolute errors in predicted formation enthalpies are found to decrease by a factor of 1.5 to 2.5 from the GGA level to the meta-GGA level. Smaller decreases in error are observed for decomposition enthalpies. For formation enthalpies r2SCAN improves over SCAN for intermetallic systems. For a few classes of systems-transition metals, intermetallics, weakly bound solids, and enthalpies of decomposition into compounds-GGAs are comparable to meta-GGAs. In total, r2SCAN and r2SCAN+rVV10 can be recommended as stable, general-purpose meta-GGAs for materials discovery.

The Predictive Power of Exact Constraints and Appropriate Norms in Density Functional Theory

Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About 20 mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that (a) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems and (b) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss both how the class of well-described systems has grown along with constraint satisfaction and the possibilities for future functional development. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 74 is April 2023. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C₂ singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B₂, C₂, O₂, and F₂ in the local spin density approximation (LSDA), the Perdew-Burke-Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C₂. This work shows that spin symmetry breaking in singlet C₂, which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

Understanding Density-Driven Errors for Reaction Barrier Heights

Delocalization errors, such as charge-transfer and some self-interaction errors, plague computationally efficient and otherwise accurate density functional approximations (DFAs). Evaluating a semilocal DFA non-self-consistently on the Hartree-Fock (HF) density is often recommended as a computationally inexpensive remedy for delocalization errors. For sophisticated meta-GGAs like SCAN, this approach can achieve remarkable accuracy. This HF-DFT (also known as DFA@HF) is often presumed to work, when it significantly improves over the DFA, because the HF density is more accurate than the self-consistent DFA density in those cases. By applying the metrics of density-corrected density functional theory (DFT), we show that HF-DFT works for barrier heights by making a localizing charge-transfer error or density overcorrection, thereby producing a somewhat reliable cancellation of density- and functional-driven errors for the energy. A quantitative analysis of the charge-transfer errors in a few randomly selected transition states confirms this trend. We do not have the exact functional and electron densities that would be needed to evaluate the exact density- and functional-driven errors for the large BH76 database of barrier heights. Instead, we have identified and employed three fully nonlocal proxy functionals (SCAN 50% global hybrid, range-separated hybrid LC-ωPBE, and SCAN-FLOSIC) and their self-consistent proxy densities. These functionals are chosen because they yield reasonably accurate self-consistent barrier heights and because their self-consistent total energies are nearly piecewise linear in fractional electron number─two important points of similarity to the exact functional. We argue that density-driven errors of the energy in a self-consistent density functional calculation are second order in the density error and that large density-driven errors arise primarily from incorrect electron transfers over length scales larger than the diameter of an atom.

Calcium Signaling Silencing in Atrial Fibrillation: Implications for Atrial Sodium Homeostasis

Atrial fibrillation (AF) is the most common type of cardiac arrhythmia, affecting more than 33 million people worldwide. Despite important advances in therapy, AF's incidence remains high, and treatment often results in recurrence of the arrhythmia. A better understanding of the cellular and molecular changes that (1) trigger AF and (2) occur after the onset of AF will help to identify novel therapeutic targets. Over the past 20 years, a large body of research has shown that intracellular Ca2+ handling is dramatically altered in AF. While some of these changes are arrhythmogenic, other changes counteract cellular arrhythmogenic mechanisms (Calcium Signaling Silencing). The intracellular Na+ concentration ([Na+])i is a key regulator of intracellular Ca2+ handling in cardiac myocytes. Despite its importance in the regulation of intracellular Ca2+ handling, little is known about [Na+]i, its regulation, and how it might be changed in AF. Previous work suggests that there might be increases in the late component of the atrial Na+ current (INa,L) in AF, suggesting that [Na+]i levels might be high in AF. Indeed, a pharmacological blockade of INa,L has been suggested as a treatment for AF. Here, we review calcium signaling silencing and changes in intracellular Na+ homeostasis during AF. We summarize the proposed arrhythmogenic mechanisms associated with increases in INa,L during AF and discuss the evidence from clinical trials that have tested the pharmacological INa,L blocker ranolazine in the treatment of AF.

Interpretations of ground-state symmetry breaking and strong correlation in wavefunction and density functional theories

Strong correlations within a symmetry-unbroken ground-state wavefunction can show up in approximate density functional theory as symmetry-broken spin densities or total densities, which are sometimes observable. They can arise from soft modes of fluctuations (sometimes collective excitations) such as spin-density or charge-density waves at nonzero wavevector. In this sense, an approximate density functional for exchange and correlation that breaks symmetry can be more revealing (albeit less accurate) than an exact functional that does not. The examples discussed here include the stretched H2 molecule, antiferromagnetic solids, and the static charge-density wave/Wigner crystal phase of a low-density jellium. Time-dependent density functional theory is used to show quantitatively that the static charge-density wave is a soft plasmon. More precisely, the frequency of a related density fluctuation drops to zero, as found from the frequency moments of the spectral function, calculated from a recent constraint-based wavevector- and frequency-dependent jellium exchange-correlation kernel.

Accurate and Numerically Efficient r2SCAN Meta-Generalized Gradient Approximation

The recently proposed rSCAN functional [ J. Chem. Phys. 2019 150, 161101] is a regularized form of the SCAN functional [ Phys. Rev. Lett. 2015 115, 036402] that improves SCAN's numerical performance at the expense of breaking constraints known from the exact exchange-correlation functional. We construct a new meta-generalized gradient approximation by restoring exact constraint adherence to rSCAN. The resulting functional maintains rSCAN's numerical performance while restoring the transferable accuracy of SCAN.

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to the leading order (-0.221Z^5/3 and -0.021Z ln Z, respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E_x(N, Z) ≈ -0.354N^2/3Z (as known before only for Z ≫ N ≫ 1) and E_c ≈ -0.02N ln N. These asymptotic estimates are most correct for atomic ions with large N and Z ≫ N, but we find that they are qualitatively and semi-quantitatively correct even for small N and N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

Identification of an INa-dependent and Ito-mediated proarrhythmic mechanism in cardiomyocytes derived from pluripotent stem cells of a Brugada syndrome patient

Brugada syndrome (BrS) is an inherited cardiac arrhythmia commonly associated with SCN5A mutations, yet its ionic mechanisms remain unclear due to a lack of cellular models. Here, we used human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) from a BrS patient (BrS1) to evaluate the roles of Na⁺ currents (I_Na) and transient outward K⁺ currents (I_to) in BrS induced action potential (AP) changes. To understand the role of these current changes in repolarization we employed dynamic clamp to “electronically express” I_K1 and restore normal resting membrane potentials and allow normal recovery of the inactivating currents, I_Na, I_Ca and I_to. HiPSC-CMs were generated from BrS1 with a compound SCN5A mutation (p. A226V & p. R1629X) and a healthy sibling control (CON1). Genome edited hiPSC-CMs (BrS2) with a milder p. T1620M mutation and a commercial control (CON2) were also studied. CON1, CON2 and BrS2, had unaltered peak I_Na amplitudes, and normal APs whereas BrS1, with over 75% loss of I_Na, displayed a loss-of-I_Na basal AP morphology (at 1.0 Hz) manifested by a reduced maximum upstroke velocity (by ~80%, p < 0.001) and AP amplitude (p < 0.001), and an increased phase-1 repolarization pro-arrhythmic AP morphology (at 0.1 Hz) in ~25% of cells characterized by marked APD shortening (~65% shortening, p < 0.001). Moreover, I_to densities of BrS1 and CON1 were comparable and increased from 1.0 Hz to 0.1 Hz by ~ 100%. These data indicate that a repolarization deficit could be a mechanism underlying BrS.

The Toxoplasma Dense Granule Proteins GRA17 and GRA23 Mediate the Movement of Small Molecules between the Host and the Parasitophorous Vacuole

Toxoplasma gondii is a protozoan pathogen in the phylum Apicomplexa that resides within an intracellular parasitophorous vacuole (PV) that is selectively permeable to small molecules through unidentified mechanisms. We have identified GRA17 as a Toxoplasma-secreted protein that localizes to the parasitophorous vacuole membrane (PVM) and mediates passive transport of small molecules across the PVM. GRA17 is related to the putative Plasmodium translocon protein EXP2 and conserved across PV-residing Apicomplexa. The PVs of GRA17-deficient parasites have aberrant morphology, reduced permeability to small molecules, and structural instability. GRA17-deficient parasites proliferate slowly and are avirulent in mice. These GRA17-deficient phenotypes are rescued by complementation with Plasmodium EXP2. GRA17 functions synergistically with a related protein, GRA23. Exogenous expression of GRA17 or GRA23 alters the membrane conductance properties of Xenopus oocytes in a manner consistent with a large non-selective pore. Thus, GRA17 and GRA23 provide a molecular basis for PVM permeability and nutrient access.

Regional variation of the inwardly rectifying potassium current in the canine heart and the contributions to differences in action potential repolarization

The inward rectifier potassium current, IK1, contributes to the terminal phase of repolarization of the action potential (AP), as well as the value and stability of the resting membrane potential. Regional variation in IK1 has been noted in the canine heart, but the biophysical properties have not been directly compared. We examined the properties and functional contribution of IK1 in isolated myocytes from ventricular, atrial and Purkinje tissue. APs were recorded from canine left ventricular midmyocardium, left atrial and Purkinje tissue. The terminal rate of repolarization of the AP in ventricle, but not in Purkinje, depended on changes in external K(+) ([K(+)]o). Isolated ventricular myocytes had the greatest density of IK1 while atrial myocytes had the lowest. Furthermore, the outward component of IK1 in ventricular cells exhibited a prominent outward component and steep negative slope conductance, which was also enhanced in 10 mM [K(+)]o. In contrast, both Purkinje and atrial cells exhibited little outward IK1, even in the presence of 10 mM [K(+)]o, and both cell types showed more persistent current at positive potentials. Expression of Kir2.1 in the ventricle was 76.9-fold higher than that of atria and 5.8-fold higher than that of Purkinje, whereas the expression of Kir2.2 and Kir2.3 subunits was more evenly distributed in Purkinje and atria. Finally, AP clamp data showed distinct contributions of IK1 for each cell type. IK1 and Kir2 subunit expression varies dramatically in regions of the canine heart and these regional differences in Kir2 expression likely underlie regional distinctions in IK1 characteristics, contributing to variations in repolarization in response to in [K(+)]o changes.

Modeling and study of the mechanism of dilated cardiomyopathy using induced pluripotent stem cells derived from individuals with Duchenne muscular dystrophy

Duchenne muscular dystrophy (DMD) is caused by mutations in the dystrophin gene (DMD), and is characterized by progressive weakness in skeletal and cardiac muscles. Currently, dilated cardiomyopathy due to cardiac muscle loss is one of the major causes of lethality in late-stage DMD patients. To study the molecular mechanisms underlying dilated cardiomyopathy in DMD heart, we generated cardiomyocytes (CMs) from DMD and healthy control induced pluripotent stem cells (iPSCs). DMD iPSC-derived CMs (iPSC-CMs) displayed dystrophin deficiency, as well as the elevated levels of resting Ca²⁺, mitochondrial damage and cell apoptosis. Additionally, we found an activated mitochondria-mediated signaling network underlying the enhanced apoptosis in DMD iPSC-CMs. Furthermore, when we treated DMD iPSC-CMs with the membrane sealant Poloxamer 188, it significantly decreased the resting cytosolic Ca²⁺ level, repressed caspase-3 (CASP3) activation and consequently suppressed apoptosis in DMD iPSC-CMs. Taken together, using DMD patient-derived iPSC-CMs, we established an in vitro model that manifests the major phenotypes of dilated cardiomyopathy in DMD patients, and uncovered a potential new disease mechanism. Our model could be used for the mechanistic study of human muscular dystrophy, as well as future preclinical testing of novel therapeutic compounds for dilated cardiomyopathy in DMD patients.

Highlighted Article: Patient-derived induced pluripotent stem cells are used to establish an in vitro model of DMD-associated cardiomyopathy that could be used for future preclinical testing.

Mechanism of automaticity in cardiomyocytes derived from human induced pluripotent stem cells

BACKGROUND AND OBJECTIVES

The creation of cardiomyocytes derived from human induced pluripotent stem cells (hiPS-CMs) has spawned broad excitement borne out of the prospects to diagnose and treat cardiovascular diseases based on personalized medicine. A common feature of hiPS-CMs is their spontaneous contractions but the mechanism(s) remain uncertain.

METHODS

Intrinsic activity was investigated by the voltage-clamp technique, optical mapping of action potentials (APs) and intracellular Ca(2+) (Cai) transients (CaiT) at subcellular-resolution and pharmacological interventions.

RESULTS

The frequency of spontaneous CaiT (sCaiT) in monolayers of hiPS-CMs was not altered by ivabradine, an inhibitor of the pacemaker current, If despite high levels of HCN transcripts (1-4). HiPS-CMs had negligible If and IK1 (inwardly-rectifying K(+)-current) and a minimum diastolic potential of -59.1±3.3mV (n=18). APs upstrokes were preceded by a depolarizing-foot coincident with a rise of Cai. Subcellular Cai wavelets varied in amplitude, propagated and died-off; larger Cai-waves triggered cellular sCaTs and APs. SCaiTs increased in frequency with [Ca(2+)]out (0.05-to-1.8mM), isoproterenol (1μM) or caffeine (100μM) (n≥5, p<0.05). HiPS-CMs became quiescent with ryanodine receptor stabilizers (K201=2μM); tetracaine; Na-Ca exchange (NCX) inhibition (SEA0400=2μM); higher [K(+)]out (5→8mM), and thiol-reducing agents but could still be electrically stimulated to elicit CaiTs. Cell-cell coupling of hiPS-CM in monolayers was evident from connexin-43 expression and CaiT propagation. SCaiTs from an ensemble of dispersed hiPS-CMs were out-of-phase but became synchronous through the outgrowth of inter-connecting microtubules.

CONCLUSIONS

Automaticity in hiPS-CMs originates from a Ca(2+)-clock mechanism involving Ca(2+) cycling across the sarcoplasmic reticulum linked to NCX to trigger APs.

Study familial hypertrophic cardiomyopathy using patient-specific induced pluripotent stem cells

Aims

Familial hypertrophic cardiomyopathy (HCM) is one the most common heart disorders, with gene mutations in the cardiac sarcomere. Studying HCM with patient-specific induced pluripotent stem-cell (iPSC)-derived cardiomyocytes (CMs) would benefit the understanding of HCM mechanism, as well as the development of personalized therapeutic strategies.

Methods and results

To investigate the molecular mechanism underlying the abnormal CM functions in HCM, we derived iPSCs from an HCM patient with a single missense mutation (Arginine442Glycine) in the MYH7 gene. CMs were next enriched from HCM and healthy iPSCs, followed with whole transcriptome sequencing and pathway enrichment analysis. A widespread increase of genes responsible for ‘Cell Proliferation’ was observed in HCM iPSC-CMs when compared with control iPSC-CMs. Additionally, HCM iPSC-CMs exhibited disorganized sarcomeres and electrophysiological irregularities. Furthermore, disease phenotypes of HCM iPSC-CMs were attenuated with pharmaceutical treatments.

Conclusion

Overall, this study explored the possible patient-specific and mutation-specific disease mechanism of HCM, and demonstrates the potential of using HCM iPSC-CMs for future development of therapeutic strategies. Additionally, the whole methodology established in this study could be utilized to study mechanisms of other human-inherited heart diseases.

Relaxin suppresses atrial fibrillation by reversing fibrosis and myocyte hypertrophy and increasing conduction velocity and sodium current in spontaneously hypertensive rat hearts

RATIONALE

Atrial fibrillation (AF) contributes significantly to morbidity and mortality in elderly and hypertensive patients and has been correlated to enhanced atrial fibrosis. Despite a lack of direct evidence that fibrosis causes AF, reversal of fibrosis is considered a plausible therapy.

OBJECTIVE

To evaluate the efficacy of the antifibrotic hormone relaxin (RLX) in suppressing AF in spontaneously hypertensive rats (SHR).

METHODS AND RESULTS

Normotensive Wistar-Kyoto (WKY) and SHR were treated for 2 weeks with vehicle (WKY+V and SHR+V) or RLX (0.4 mg/kg per day, SHR+RLX) using implantable mini-pumps. Hearts were perfused, mapped optically to analyze action potential durations, intracellular Ca²⁺ transients, and restitution kinetics, and tested for AF vulnerability. SHR hearts had slower conduction velocity (CV; P<0.01 versus WKY), steeper CV restitution kinetics, greater collagen deposition, higher levels of transcripts for transforming growth factor-β, metalloproteinase-2, metalloproteinase-9, collagen I/III, and reduced connexin 43 phosphorylation (P<0.05 versus WKY). Programmed stimulation triggered sustained AF in SHR (n=5/5) and SHR+V (n=4/4), but not in WKY (n=0/5) and SHR+RLX (n=1/8; P<0.01). RLX treatment reversed the transcripts for fibrosis, flattened CV restitution kinetics, reduced action potential duration at 90% recovery to baseline, increased CV (P<0.01), and reversed atrial hypertrophy (P<0.05). Independent of antifibrotic actions, RLX (0.1 µmol/L) increased Na⁺ current density, INa (≈2-fold in 48 hours) in human cardiomyocytes derived from inducible pluripotent stem cells (n=18/18; P<0.01).

CONCLUSIONS

RLX treatment suppressed AF in SHR hearts by increasing CV from a combination of reversal of fibrosis and hypertrophy and by increasing INa. The study provides compelling evidence that RLX may provide a novel therapy to manage AF in humans by reversing fibrosis and hypertrophy and by modulating cardiac ionic currents.

Patient-specific induced pluripotent stem-cell-derived models of LEOPARD syndrome

Generation of reprogrammed induced pluripotent stem cells (iPSC) from patients with defined genetic disorders promises important avenues to understand the etiologies of complex diseases, and the development of novel therapeutic interventions. We have generated iPSC from patients with LEOPARD syndrome (LS; acronym of its main features: Lentigines, Electrocardiographic abnormalities, Ocular hypertelorism, Pulmonary valve stenosis, Abnormal genitalia, Retardation of growth and Deafness), an autosomal dominant developmental disorder belonging to a relatively prevalent class of inherited RAS-MAPK signaling diseases, which also includes Noonan syndrome (NS), with pleiomorphic effects on several tissues and organ systems1,2. The patient-derived cells have a mutation in the PTPN11 gene, which encodes the SHP2 phosphatase. The iPSC have been extensively characterized and produce multiple differentiated cell lineages. A major disease phenotype in patients with LEOPARD syndrome is hypertrophic cardiomyopathy. We show that in vitro-derived cardiomyocytes from LS-iPSC are larger, have a higher degree of sarcomeric organization and preferential localization of NFATc4 in the nucleus when compared to cardiomyocytes derived from human embryonic stem cells (HESC) or wild type (wt) iPSC derived from a healthy brother of one of the LS patients. These features correlate with a potential hypertrophic state. We also provide molecular insights into signaling pathways that may promote the disease phenotype.

The cardiomyocyte lineage is critical for optimization of stem cell therapy in a mouse model of myocardial infarction

We recently described a murine embryonic stem cell (ESC) line engineered to express the activated Notch 4 receptor in a tetracycline (doxcycline; Dox) regulated fashion (tet-notch4 ESCs). Notch 4 induction in Flk1(+) hematopoietic and vascular progenitors from this line respecified them to a cardiovascular fate. We reasoned that these cells would be ideal for evaluating the contribution of the cardiomyocyte and vascular lineages to the functional improvement noted following stem cell transplantation in infarcted hearts. Flk-1(+) Tet-notch4 cells from d 3 embryoid bodies exposed to doxycycline (Dox(+)) were compared to uninduced (Dox(-)) Flk-1(+) cells. Mice underwent transplantation of 5 x 10(5) Dox(+) cells, Dox(-)cells, or an equal volume of serum-free medium after surgically induced myocardial infarction. The mean ejection fraction was 59 + or - 15, 46 + or - 17, and 39 + or - 13% in the Dox(+), Dox(-), and serum-free medium groups, respectively (P<0.05 for the differences among all 3 groups). Immunohistochemistry of hearts injected with Dox(+) grafts expressed myocardial and vascular markers, whereas grafts of Dox(-) cells expressed primarily vascular markers. We conclude that cardiovascular progenitors are more effective than vascular progenitors in improving function after myocardial infarction. The transplantation of appropriate cell types is critical for maximizing the benefit of cardiovascular cell therapy.-Adler, E. D., Chen, V. C., Bystrup, A., Kaplan, A. D., Giovannone, S., Briley-Saebo, K., Young, W., Kattman, S., Mani, V., Laflamme, M., Zhu, W.-Z., Fayad, Z., Keller, G. The cardiomyocyte lineage is critical for optimization of stem cell therapy in a mouse model of myocardial infarction.

Extracellular nucleotides act through P2U purinoceptors to elevate [Ca2+]i and enhance basic fibroblast growth factor-induced proliferation in sheep chondrocytes

Extracellular nucleotides interact with specific cell surface receptors to mediate a variety of biological responses, including elevation of the cytosolic free Ca2+ concentration ([Ca2+]i) in a number of cell types. Although extracellular ATP has been shown to affect chondrocyte function, the underlying mechanisms are poorly understood. In the present study, we investigated whether Ca2+-mobilizing purinoceptors are present on sheep chondrocytes. Chondrocytes were isolated from the proximal tibial growth plate of day 120-130 sheep fetuses. Early passage cells were loaded with indo-1 or fluo-3, and [Ca2+]i was monitored by fluorescence spectrophotometry. ATP (0.3-100 microM) induced transient elevation of [Ca2+]i, lasting approximately 1 min. Half-maximal elevation of [Ca2+]i was observed at an ATP concentration of 5.0 +/- 0.2 microM. Responses were still observed in the absence of extracellular Ca2+, and were abolished by pretreatment with thapsigargin, consistent with the release of Ca2+ from intracellular stores. Several nucleotides were tested for their ability to elevate [Ca2+]i. In order of potency, these were UTP approximately ATP >> ADP approximately 2-methylthio-ATP. No responses were elicited by benzoylbenzoic-ATP, a P2Z-selective agonist; alpha,beta-methylene-ATP, an agonist selective for certain P2X purinoceptors; AMP; adenosine; or pyrophosphate (all at 100 microM), demonstrating specificity. Taken together, these data indicate that nucleotides elevate [Ca2+]i in chondrocytes through interaction with the P2U purinoceptor subtype. Although pretreatment with pertussis toxin virtually abolished the Ca2+ response to lysophosphatidic acid, the response to UTP was relatively insensitive, suggesting that P2U purinoceptors are not linked to a pertussis toxin-sensitive G protein in chondrocytes. In contrast, the Ca2+ response to UTP was markedly inhibited by the biologically active phorbol ester 12-O-tetradecanoyl-beta-phorbol 13-acetate, but not by the inactive control compound 4 alpha-phorbol 12,13-didecanoate, suggesting that a 12-O-tetradecanoyl-beta-phorbol 13-acetate-sensitive isoform of protein kinase C regulates P2U purinoceptor signaling in these cells. UTP (10 microM) enhanced the proliferative response to basic fibroblast growth factor. The response to basic fibroblast growth factor was also enhanced by ATP, but not by 2-methylthio-ATP, consistent with involvement of P2U purinoceptors. Nucleotides released during trauma, inflammation, or cell death may act through P2U purinoceptors to regulate chondrocyte function in an autocrine or paracrine manner.

Experimental autoimmune anterior uveitis. Induction with melanin-associated antigen from the iris and ciliary body

PURPOSE

This study was designed to investigate an animal model of uveitis that resembles anterior uveitis in humans after immunization with iris-ciliary body antigen.

METHODS

Male Lewis rats 6 to 8 weeks of age were immunized with the buffer- and detergent-insoluble bovine iris-ciliary body antigen mixed with complete Freund's adjuvant and pertussis toxin. Antigen was digested with various proteolytic enzymes and tested in different rodent strains for a uveitogenic response.

RESULTS

Acute iridocyclitis developed in both eyes of the Lewis rat during the second week after immunization, and the pattern of inflammation was similar to acute anterior uveitis in humans, with sudden onset, localization to the anterior uvea, and spontaneous resolution. Among the strains tested, F344 rats were susceptible to experimental autoimmune anterior uveitis but Long-Evans rats were not. Experimental autoimmune anterior uveitis did not develop in any of the mice studied, nor was it induced by immunization with synthetic melanin, amelanotic bovine tissues, pigmented bovine skin, or pigmented rat and rabbit iris-ciliary body. A soluble fraction derived from bovine melanin-associated antigen (BMAA) after digestion with the proteolytic enzyme V8 protease resulted in a disease similar to that observed with intact BMAA.

CONCLUSIONS

A model of anterior uveitis has been induced in the Lewis rat after immunization with bovine uveal antigen, and it resembles the acute iridocyclitis observed in humans. These results suggest that the pathogenic antigen is a melanin-associated protein(s) present within the iris-ciliary body.

Extracellular nucleotides potentiate the cytosolic Ca2+, but not cyclic adenosine 3', 5'-monophosphate response to parathyroid hormone in rat osteoblastic cells

Binding to PTH to its cell surface receptor activates both adenylyl cyclase and phospholipase-C, leading to elevation of cytosolic cAMP and free Ca2+. We have shown previously that extracellular nucleotides interact with P2U and P2Y subtypes of purinoceptor on osteoblastic cells, both linked to Ca2+ mobilization. In the present study, we investigated possible interactions between nucleotide and PTH signaling pathways in osteoblastic cells. The cytosolic free Ca2+ concentration ([Ca2+]i) of UMR-106 osteoblastic cells was monitored by fluorescence spectrophotometry. PTH (0.01-1 microM; bovine 1-84 or human 1-34) induced a small transient elevation of [Ca2+]i, lasting less than 1 min. A number of nucleotides, including ATP, UTP, and UDP, induced transient elevation of [Ca2+]i and potentiated the subsequent Ca2+ response to PTH. Of the nucleotides tested, UDP was the most effective at potentiating the PTH-induced Ca2+ transient. Treatment of cells with UDP (100 microM for 2.5 min), but not inorganic phosphate or uridine, reversibly potentiated the Ca2+ response to PTH (0.1 microM) by 11 +/- 2-fold (mean +/- SEM; n = 39). In contrast, UDP did not affect the cAMP response to PTH, indicating a selective action on Ca2+ signaling. Potentiation of the Ca2+ signal was still observed in the absence of extracellular Ca2+, establishing that nucleotides enhance PTH-induced release of Ca2+ from intracellular stores. Studies using selective purinoceptor agonists suggest that potentiation of PTH signaling is mediated by the P2U receptor subtype. In vivo, nucleotides released during trauma or inflammation may modulate PTH-induced Ca2+ signaling in osteoblasts.