Programming Cell-Derived Vesicles with Enhanced Immunomodulatory Properties.

Tumor-associated macrophages are the predominant immune cells present in the tumor microenvironment and mostly exhibit a pro-tumoral M2-like phenotype. However, macrophage biology is reversible allowing them to acquire an anti-tumoral M1-like phenotype in response to external stimuli. A potential therapeutic strategy for treating cancer may be achieved by modulating macrophages from an M2 to an M1-like phenotype with the tumor microenvironment. Here, programmed nanovesicles are generated as an immunomodulatory therapeutic platform with the capability to re-polarize M2 macrophages toward a proinflammatory phenotype. Programmed nanovesicles are engineered from cellular membranes to have specific immunomodulatory properties including the capability to bidirectionally modulate immune cell polarization. These programmed nanovesicles decorated with specific membrane-bound ligands can be targeted toward specific cell types including immune cells. Macrophage-derived vesicles are engineered to enhance immune cell reprogramming toward a proinflammatory phenotype.

Inhibitory and stimulatory micropeptides preferentially bind to different conformations of the cardiac calcium pump.

The ATP-dependent ion pump sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) sequesters Ca2+ in the endoplasmic reticulum to establish a reservoir for cell signaling. Because of its central importance in physiology, the activity of this transporter is tightly controlled via direct interactions with tissue-specific regulatory micropeptides that tune SERCA function to match changing physiological conditions. In the heart, the micropeptide phospholamban (PLB) inhibits SERCA, while dwarf open reading frame (DWORF) stimulates SERCA. These competing interactions determine cardiac performance by modulating the amplitude of Ca2+ signals that drive the contraction/relaxation cycle. We hypothesized that the functions of these peptides may relate to their reciprocal preferences for SERCA binding; SERCA binds PLB more avidly at low cytoplasmic [Ca2+] but binds DWORF better when [Ca2+] is high. In the present study, we demonstrated this opposing Ca2+ sensitivity is due to preferential binding of DWORF and PLB to different intermediate states that SERCA samples during the Ca2+ transport cycle. We show PLB binds best to the SERCA E1-ATP state, which prevails at low [Ca2+]. In contrast, DWORF binds most avidly to E1P and E2P states that are more populated when Ca2+ is elevated. Moreover, FRET microscopy revealed dynamic shifts in SERCA-micropeptide binding equilibria during cellular Ca2+ elevations. A computational model showed that DWORF exaggerates changes in PLB-SERCA binding during the cardiac cycle. These results suggest a mechanistic basis for inhibitory versus stimulatory micropeptide function, as well as a new role for DWORF as a modulator of dynamic oscillations of PLB-SERCA regulatory interactions.

Fluorescence lifetime imaging microscopy reveals sodium pump dimers in live cells.

The sodium-potassium ATPase (Na/K-ATPase, NKA) establishes ion gradients that facilitate many physiological functions including action potentials and secondary transport processes. NKA comprises a catalytic subunit (alpha) that interacts closely with an essential subunit (beta) and regulatory transmembrane micropeptides called FXYD proteins. In the heart, a key modulatory partner is the FXYD protein phospholemman (PLM, FXYD1), but the stoichiometry of the alpha-beta-PLM regulatory complex is unknown. Here, we used fluorescence lifetime imaging and spectroscopy to investigate the structure, stoichiometry, and affinity of the NKA-regulatory complex. We observed a concentration-dependent binding of the subunits of NKA-PLM regulatory complex, with avid association of the alpha subunit with the essential beta subunit as well as lower affinity alpha-alpha and alpha-PLM interactions. These data provide the first evidence that, in intact live cells, the regulatory complex is composed of two alpha subunits associated with two beta subunits, decorated with two PLM regulatory subunits. Docking and molecular dynamics (MD) simulations generated a structural model of the complex that is consistent with our experimental observations. We propose that alpha-alpha subunit interactions support conformational coupling of the catalytic subunits, which may enhance NKA turnover rate. These observations provide insight into the pathophysiology of heart failure, wherein low NKA expression may be insufficient to support formation of the complete regulatory complex with the stoichiometry (alpha-beta-PLM)2.

Modulation of P2X4 pore closure by magnesium, potassium, and ATP.

The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 µs we have investigated how metal-complexed ATP stabilizes the channel's open state and prevents its closing. We have identified two metal-binding sites, Mg2+ and potassium K+, one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP-binding site (MBS2), similar to those characterized in Gulf Coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 impeded the collapse of the open-state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and the pore collapsed. This collapse was apparent when fewer than two ATPs were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study therefore provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+.

Crowding within synaptic junctions influences the degradation of nucleotides by CD39 and CD73 ectonucleotidases.

Synapsed cells can communicate using exocytosed nucleotides like adenosine triphosphate (ATP). Ectonucleotidases localized to synaptic junctions degrade nucleotides into metabolites like adenosine monophosphate (AMP) or adenosine. Oftentimes nucleotide degradation occurs in a sequential manner, of which ATP degradation by CD39 and CD73 is a representative example. Here, CD39 first converts ATP and adenosine diphosphate (ADP) into AMP, after which AMP is dephosphorylated into adenosine by CD73. Hence, the concerted activity of CD39 and CD73 can help shape cellular responses to extracellular ATP. In a previous study, we demonstrated that coupled CD39 and CD73 activity within synapse-like junctions is strongly controlled by the enzymes' co-localization, their surface charge densities, and the electrostatic potential of the surrounding cell membranes. In this study, we demonstrate that crowders within synaptic junctions, which can include globular proteins like cytokines and membrane-bound proteins, impact coupled CD39 and CD73 ectonucleotidase activity and, in turn, the availability of intrasynapse ATP. Specifically, we developed a spatially explicit, reaction-diffusion model for the coupled conversion of ATP → AMP and AMP → adenosine in a model synaptic junction with crowders that is solved via the finite element method. Our modeling results suggest that the association rate for ATP to CD39 is strongly influenced by the density of intrasynaptic protein crowders, as increasing crowder density generally suppressed ATP association kinetics. Much of this suppression can be rationalized based on a loss of configurational entropy. The surface charges of crowders can further influence the association rate, with the surprising result that favorable crowder-nucleotide electrostatic interactions can yield CD39 association rates that are faster than crowder-free configurations. However, attractive crowder-nucleotide interactions decrease the rate and efficiency of adenosine production, which in turn increases the availability of ATP and AMP within the synapse relative to crowder-free configurations. These findings highlight how CD39 and CD73 ectonucleotidase activity, electrostatics, and crowding within synapses influence the availability of nucleotides for intercellular communication.

Purinoreceptors and ectonucleotidases control ATP-induced calcium waveforms and calcium-dependent responses in microglia: Roles of P2 receptors and CD39 in ATP-stimulated microglia.

Adenosine triphosphate (ATP) and its metabolites drive microglia migration and cytokine production by activating P2X- and P2Y- class purinergic receptors. Purinergic receptor activation gives rise to diverse intracellular calcium (Ca2+ signals, or waveforms, that differ in amplitude, duration, and frequency. Whether and how these characteristics of diverse waveforms influence microglia function is not well-established. We developed a computational model trained with data from published primary murine microglia studies. We simulate how purinoreceptors influence Ca2+ signaling and migration, as well as, how purinoreceptor expression modifies these processes. Our simulation confirmed that P2 receptors encode the amplitude and duration of the ATP-induced Ca2+ waveforms. Our simulations also implicate CD39, an ectonucleotidase that rapidly degrades ATP, as a regulator of purinergic receptor-induced Ca2+ responses. Namely, it was necessary to account for CD39 metabolism of ATP to align the model's predicted purinoreceptor responses with published experimental data. In addition, our modeling results indicate that small Ca2+ transients accompany migration, while large and sustained transients are needed for cytokine responses. Lastly, as a proof-of-principal, we predict Ca2+ transients and cell membrane displacements in a BV2 microglia cell line using published P2 receptor mRNA data to illustrate how our computer model may be extrapolated to other microglia subtypes. These findings provide important insights into how differences in purinergic receptor expression influence microglial responses to ATP.

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.

A Matched-Filter-Based Algorithm for Subcellular Classification of T-System in Cardiac Tissues.

In mammalian ventricular cardiomyocytes, invaginations of the surface membrane form the transverse tubular system (T-system), which consists of transverse tubules (TTs) that align with sarcomeres and Z-lines as well as longitudinal tubules (LTs) that are present between Z-lines in some species. In many cardiac disease etiologies, the T-system is perturbed, which is believed to promote spatially heterogeneous, dyssynchronous Ca2+ release and inefficient contraction. In general, T-system characterization approaches have been directed primarily at isolated cells and do not detect subcellular T-system heterogeneity. Here, we present MatchedMyo, a matched-filter-based algorithm for subcellular T-system characterization in isolated cardiomyocytes and millimeter-scale myocardial sections. The algorithm utilizes "filters" representative of TTs, LTs, and T-system absence. Application of the algorithm to cardiomyocytes isolated from rat disease models of myocardial infarction (MI), dilated cardiomyopathy induced via aortic banding, and sham surgery confirmed and quantified heterogeneous T-system structure and remodeling. Cardiomyocytes from post-MI hearts exhibited increasing T-system disarray as proximity to the infarct increased. We found significant (p < 0.05, Welch's t-test) increases in LT density within cardiomyocytes proximal to the infarct (12 ± 3%, data reported as mean ± SD, n = 3) versus sham (4 ± 2%, n = 5), but not distal to the infarct (7 ± 1%, n = 3). The algorithm also detected decreases in TTs within 5° of the myocyte minor axis for isolated aortic banding (36 ± 9%, n = 3) and MI cardiomyocytes located intermediate (37 ± 4%, n = 3) and proximal (34 ± 4%, n = 3) to the infarct versus sham (57 ± 12%, n = 5). Application of bootstrapping to rabbit MI tissue revealed distal sections comprised 18.9 ± 1.0% TTs, whereas proximal sections comprised 10.1 ± 0.8% TTs (p < 0.05), a 46.6% decrease. The matched-filter approach therefore provides a robust and scalable technique for T-system characterization from isolated cells through millimeter-scale myocardial sections.

Simulation of P2X-mediated calcium signalling in microglia.

KEY POINTS: A computational model of P2X channel activation in microglia was developed that includes downfield Ca2+ -dependent signalling pathways. This model provides quantitative insights into how diverse signalling pathways in microglia converge to control microglial function. ABSTRACT: Microglia function is orchestrated through highly coupled signalling pathways that depend on calcium (Ca2+ ). In response to extracellular ATP, transient increases in intracellular Ca2+ driven through the activation of purinergic receptors, P2X and P2Y, are sufficient to promote cytokine synthesis. Although the steps comprising the pathways bridging purinergic receptor activation with transcriptional responses have been probed in great detail, a quantitative model for how these steps collectively control cytokine production has not been established. Here we developed a minimal computational model that quantitatively links extracellular stimulation of two prominent ionotropic purinergic receptors, P2X4 and P2X7, with the graded production of a gene product, namely the tumour necrosis factor α (TNFα) cytokine. In addition to Ca2+ handling mechanisms common to eukaryotic cells, our model includes microglia-specific processes including ATP-dependent P2X4 and P2X7 activation, activation of nuclear factor of activated T-cells (NFAT) transcription factors, and TNFα production. Parameters for this model were optimized to reproduce published data for these processes, where available. With this model, we determined the propensity for TNFα production in microglia, subject to a wide range of ATP exposure amplitudes, frequencies and durations that the cells could encounter in vivo. Furthermore, we have investigated the extent to which modulation of the signal transduction pathways influence TNFα production. Our results suggest that pulsatile stimulation of P2X4 via micromolar ATP may be sufficient to promote TNFα production, whereas high-amplitude ATP exposure is necessary for production via P2X7. Furthermore, under conditions that increase P2X4 expression, for instance, following activation by pathogen-associated molecular factors, P2X4-associated TNFα production is greatly enhanced. Given that Ca2+ homeostasis in microglia is profoundly important to its function, this computational model provides a quantitative framework to explore hypotheses pertaining to microglial physiology.

Computational modeling of amylin-induced calcium dysregulation in rat ventricular cardiomyocytes.

Hyperamylinemia is a condition that accompanies obesity and precedes type II diabetes, and it is characterized by above-normal blood levels of amylin, the pancreas-derived peptide. Human amylin oligomerizes easily and can deposit in the pancreas [1], brain [2], and heart [3], where they have been associated with calcium dysregulation. In the heart, accumulating evidence suggests that human amylin oligomers form moderately cation-selective [4,5] channels that embed in the cell sarcolemma (SL). The oligomers increase membrane conductance in a concentration-dependent manner [5], which is correlated with elevated cytosolic Ca2+. These findings motivate our core hypothesis that non-selective inward Ca2+ conduction afforded by human amylin oligomers increase cytosolic and sarcoplasmic reticulum (SR) Ca2+ load, which thereby magnifies intracellular Ca2+ transients. Questions remain however regarding the mechanism of amylin-induced Ca2+ dysregulation, including whether enhanced SL Ca2+ influx is sufficient to elevate cytosolic Ca2+ load [6], and if so, how might amplified Ca2+ transients perturb Ca2+-dependent cardiac pathways. To investigate these questions, we modified a computational model of cardiomyocytes Ca2+ signaling to reflect experimentally-measured changes in SL membrane permeation and decreased sarcoplasmic/endoplasmic reticulum calcium ATPase (SERCA) function stemming from acute and transgenic human amylin peptide exposure. With this model, we confirmed the hypothesis that increasing SL permeation alone was sufficient to enhance Ca2+ transient amplitudes. Our model indicated that amplified cytosolic transients are driven by increased Ca2+ loading of the SR and that greater fractional release may contribute to the Ca2+-dependent activation of calmodulin, which could prime the activation of myocyte remodeling pathways. Importantly, elevated Ca2+ in the SR and dyadic space collectively drive greater fractional SR Ca2+ release for human amylin expressing rats (HIP) and acute amylin-exposed rats (+Amylin) mice, which contributes to the inotropic rise in cytosolic Ca2+ transients. These findings suggest that increased membrane permeation induced by oligomeratization of amylin peptide in cell sarcolemma contributes to Ca2+ dysregulation in pre-diabetes.

Molecular and subcellular-scale modeling of nucleotide diffusion in the cardiac myofilament lattice.

Contractile 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.

A molecular docking study of estrogenically active compounds with 1,2-diarylethane and 1,2-diarylethene pharmacophores.

Numerous selective estrogen receptor modulators (SERMs) have been synthesized and assayed in recent years. The focus of this study is to apply coarse-grain molecular docking procedures coupled with fine-grain all-atom force field optimization strategies to shed light on the binding mechanisms of currently available estrogen receptor-active compounds. Although the mechanics of ligand binding in estrogen receptors is generally well understood, there is room for surprises. In this paper computational evidence corroborating the experimentally observed type I agonistic binding mode for estradiol (E2) and diethylstilbesterol (DES) and the type II antagonistic binding mode for 4-hydroxytamoxifen and raloxifen is presented. Included in this type I agonistic mode are the DES derivatives, transstilbene and 1,2-diaryldiaminoethane. In addition, a novel 'type II agonistic' binding mode for 2,3-diarylimidazolines, 4,5-diarylimidazoles, 2,3-diarylpiperazines is introduced. This mode is stabilized by suggesting alternative hydrogen bond anchor points in the ligand binding domain as potential leads for future drug design.