Search

Highly pathogenic avian influenza (HPAI) H5N1 has been associated with significant morbidity and mortality. Given recent detections of HPAI H5N1 in dairy cattle and RNA detections in pasteurized retail milk in the United States, we established the Pan-Canadian Milk Network in April 2024. Through this, retail milk was procured longitudinally and sent to a central laboratory to test for the presence of influenza A virus RNA. We tested 109 retail milk samples from all ten Canadian provinces and all samples tested negative. Our independent testing results have aligned with reporting from federal retail milk testing initiatives. Despite no known HPAI infections of dairy cattle in Canada to date, H5N1 poses a significant threat to the health of both humans and other animals. By performing routine surveillance of retail milk on a national scale, we have shown that academic networks and initiatives can rapidly establish nationwide emerging infectious disease surveillance that is cost-effective, standardized, scalable, and easily accessible. Our network can serve as an early detection system to help inform containment and mitigation activities if positive samples are identified and can be readily reactivated should H5N1 or other emerging zoonotic viruses be identified in agricultural or livestock settings.

Cellular molecules interact with one another in a structured manner, defining a regulatory network topology that describes cellular mechanisms. Genetic mutations alter these networks' pathways, generating complex disorders such as autism spectrum disorder (ASD). Boolean models have assisted in understanding biological system dynamics since Kauffman's 1969 discovery, and various analytical tools for regulatory networks have been developed. This study examined the protein-protein interaction network created in our previous publication of four ASD patients using the SPIDDOR R package, a Boolean model-based method. The aim is to examine how patients' genetic variations in INTS6L, USP9X, RSK4, FGF5, FLNA, SUMF1, and IDS affect mTOR and Wnt cell signaling convergence. The Boolean network analysis revealed abnormal activation levels of essential proteins such as β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD. These proteins affect gene expression, translation, cell adhesion, shape, and migration. Patients 1 and 2 showed consistent patterns of increased β-catenin activity and decreased MTORC1, RPS6, and eIF4E activity. However, patient 2 had an independent decrease in Cadherin and SMAD activity due to the FLNA mutation. Patients 3 and 4 have an abnormal activation of the mTOR pathway, which includes the MTORC1, RPS6, and eIF4E genes. The shared mTOR pathway behavior in these patients is explained by a shared mutation in two closely related proteins (SUMF1 and IDS). Diverse activities in β-catenin, MTORC1, RPS6, eIF4E, Cadherin, and SMAD contributed to the reported phenotype in these individuals. Furthermore, it unveiled the potential therapeutic options that could be suggested to these individuals.

Background: Ovarian cancer (OC) is the most lethal gynecological cancer in the United States. Among the different types of OC, serous ovarian cancer (SOC) stands out as the most prevalent. Transcriptomics techniques generate extensive gene expression data, yet only a few of these genes are relevant to clinical diagnosis., Methods: Methods for feature selection (FS) address the challenges of high dimensionality in extensive datasets. This study proposes a computational framework that applies FS techniques to identify genes highly associated with platinum-based chemotherapy response on SOC patients. Using SOC datasets from the Gene Expression Omnibus (GEO) database, LASSO and varSelRF FS methods were employed. Machine learning classification algorithms such as random forest (RF) and support vector machine (SVM) were also used to evaluate the performance of the models., Results: The proposed framework has identified biomarkers panels with 9 and 10 genes that are highly correlated with platinum-paclitaxel and platinum-only response in SOC patients, respectively. The predictive models have been trained using the identified gene signatures and accuracy of above 90% was achieved., Conclusions: In this study, we propose that applying multiple feature selection methods not only effectively reduces the number of identified biomarkers, enhancing their biological relevance, but also corroborates the efficacy of drug response prediction models in cancer treatment., Competing Interests: Conflicts of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Introduction: FOLFOX and FOLFIRI chemotherapy are considered standard first-line treatment options for colorectal cancer (CRC). However, the criteria for selecting the appropriate treatments have not been thoroughly analyzed. Methods: A newly developed machine learning model was applied on several gene expression data from the public repository GEO database to identify molecular signatures predictive of efficacy of 5-FU based combination chemotherapy (FOLFOX and FOLFIRI) in patients with CRC. The model was trained using 5-fold cross validation and multiple feature selection methods including LASSO and VarSelRF methods. Random Forest and support vector machine classifiers were applied to evaluate the performance of the models. Results and Discussion: For the CRC GEO dataset samples from patients who received either FOLFOX or FOLFIRI, validation and test sets were >90% correctly classified (accuracy), with specificity and sensitivity ranging between 85%-95%. In the datasets used from the GEO database, 28.6% of patients who failed the treatment therapy they received are predicted to benefit from the alternative treatment. Analysis of the gene signature suggests the mechanistic difference between colorectal cancers that respond and those that do not respond to FOLFOX and FOLFIRI. Application of this machine learning approach could lead to improvements in treatment outcomes for patients with CRC and other cancers after additional appropriate clinical validation., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Amniouel and Jafri.)

Text mining methods are being developed to assimilate the volume of biomedical textual materials that are continually expanding. Understanding protein-protein interaction (PPI) deficits would assist in explaining the genesis of diseases. In this study, we designed an automated system to extract PPIs from the biomedical literature that uses a deep learning sentence classification model, a pretrained word embedding, and a BiLSTM recurrent neural network with additional layers, a conditional random field (CRF) named entity recognition (NER) model, and shortest-dependency path (SDP) model using the SpaCy library in Python. The automated system ensures that it targets sentences that contain PPIs and not just these proteins mentioned in the framework of disease discovery or other context. Our first model achieved 13% greater precision on the Aimed/BioInfr benchmark corpus than the previous state-of-the-art BiLSTM neural network models. The NER model presented in this study achieved 98% precision on the Aimed/BioInfr corpus over previous models. In order to facilitate the production of an accurate representation of the PPI network, the processes were developed to systematically map the protein interactions in the texts. Overall, evaluating our system through the use of 6027 abstracts pertaining to seven proteins associated with Autism Spectrum Disorder completed the manually curated PPI network for these proteins. When it comes to complicated diseases, these networks would assist in understanding how PPI deficits contribute to disease development while also emphasizing the influence of interactions on protein function and biological processes.

Calcium (Ca 2+ ) sparks are the elementary events of excitation-contraction coupling, yet they are not explicitly represented in human ventricular myocyte models. A stochastic ventricular cardiomyocyte human model that adapts to intracellular Ca 2+ ([Ca 2+ ] i ) dynamics, spark regulation, and frequency-dependent changes in the form of locally controlled Ca 2+ release was developed. The 20,000 CRUs in this model are composed of 9 individual LCCs and 49 RyRs that function as couplons. The simulated action potential duration at 1 Hz steady-state pacing is ~0.280 s similar to human ventricular cell recordings. Rate-dependence experiments reveal that APD shortening mechanisms are largely contributed by the L-type calcium channel inactivation, RyR open fraction, and [Ca 2+ ] myo concentrations. The dynamic slow-rapid-slow pacing protocol shows that RyR open probability during high pacing frequency (2.5 Hz) switches to an adapted "nonconducting" form of Ca 2+ -dependent transition state. The predicted force was also observed to be increased in high pacing, but the SR Ca 2+ fractional release was lower due to the smaller difference between diastolic and systolic [Ca 2+ ] SR . Restitution analysis through the S1S2 protocol and increased LCC Ca 2+ -dependent activation rate show that the duration of LCC opening helps modulate its effects on the APD restitution at different diastolic intervals. Ultimately, a longer duration of calcium sparks was observed in relation to the SR Ca 2+ loading at high pacing rates. Overall, this study demonstrates the spontaneous Ca 2+ release events and ion channel responses throughout various stimuli.

Mitochondria in mammalian cardiomyocytes display considerable structural heterogeneity, the significance of which is not currently understood. We use electron microscopic tomography to analyze a dataset of 68 mitochondrial subvolumes to look for correlations among mitochondrial size and shape, crista morphology and membrane density, and organelle location within rat cardiac myocytes. A tomographic analysis guided the definition of four classes of crista morphology: lamellar, tubular, mixed and transitional, the last associated with remodeling between lamellar and tubular cristae. Correlations include an apparent bias for mitochondria with lamellar cristae to be located in the regions between myofibrils and a two-fold larger crista membrane density in mitochondria with lamellar cristae relative to mitochondria with tubular cristae. The examination of individual cristae inside mitochondria reveals local variations in crista topology, such as extent of branching, alignment of fenestrations and progressive changes in membrane morphology and packing density. The findings suggest both a rationale for the interfibrillar location of lamellar mitochondria and a pathway for crista remodeling from lamellar to tubular morphology.

Mutations in the calcium-sensing protein calmodulin ( CaM ) have been linked to two cardiac arrhythmia diseases, Long QT Syndrome 14 (LQT14) and Catecholaminergic Polymorphic Ventricular Tachycardia Type 4 (CPVT4), with varying degrees of severity. Functional characterization of the CaM mutants most strongly associated with LQT14 show a clear disruption of the calcium-dependent inactivation (CDI) of the L-Type calcium channel (LCC). CPVT4 mutants on the other hand are associated with changes in their affinity to the ryanodine receptor. In clinical studies, some variants have been associated with both CPVT4 and LQT15. This study uses simulations in a model for excitation-contraction coupling in the rat ventricular myocytes to understand how LQT14 variant might give the functional phenotype similar to CPVT4. Changing the CaM -dependent transition rate by a factor of 0.75 corresponding to the D96V variant and by a factor of 0.90 corresponding to the F142L or N98S variants, in a physiologically based stochastic model of the LCC prolonger, the action potential duration changed by a small amount in a cardiac myocyte but did not disrupt CICR at 1, 2, and 4 Hz. Under beta-adrenergic simulation abnormal excitation-contraction coupling was observed above 2 Hz pacing for the mutant CaM . The same conditions applied under beta-adrenergic stimulation led to the rapid onset of arrhythmia in the mutant CaM simulations. Simulations with the LQT14 mutations under the conditions of rapid pacing with beta-adrenergic stimulation drives the cardiac myocyte toward an arrhythmic state known as Ca 2+ overload. These simulations provide a mechanistic link to a disease state for LQT14-associated mutations in CaM to yield a CPVT4 phenotype. The results show that small changes to the CaM -regulated inactivation of LCC promote arrhythmia and underscore the significance of CDI in proper heart function.

Calsequestrin Type 2 (CASQ2) is a high-capacity, low-affinity, Ca 2+ -binding protein expressed in the sarcoplasmic reticulum (SR) of the cardiac myocyte. Mutations in CASQ2 have been linked to the arrhythmia catecholaminergic polymorphic ventricular tachycardia (CPVT2) that occurs with acute emotional stress or exercise can result in sudden cardiac death (SCD). CASQ2 G112+5X is a 16 bp (339-354) deletion CASQ2 mutation that prevents the protein expression due to premature stop codon. Understanding the subcellular mechanisms of CPVT2 is experimentally challenging because the occurrence of arrhythmia is rare. To obtain an insight into the characteristics of this rare disease, simulation studies using a local control stochastic computational model of the Guinea pig ventricular myocyte investigated how the mutant CASQ2s may be responsible for the development of an arrhythmogenic episode under the condition of β-adrenergic stimulation or in the slowing of heart rate afterward once β-adrenergic stimulation ceases. Adjustment of the computational model parameters based upon recent experiments explore the functional changes caused by the CASQ2 mutation. In the simulation studies under rapid pacing (6 Hz), electromechanically concordant cellular alternans appeared under β-adrenergic stimulation in the CPVT mutant but not in the wild-type nor in the non-β-stimulated mutant. Similarly, the simulations of accelerating pacing from slow to rapid and back to the slow pacing did not display alternans but did generate early afterdepolarizations (EADs) during the period of second slow pacing subsequent acceleration of rapid pacing.

Protein phosphorylation is a post-translational modification that enables various cellular activities and plays essential roles in protein interactions. Phosphorylation is an important process for the replication of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). To shed more light on the effects of phosphorylation, we used an ensemble of neural networks to predict potential kinases that might phosphorylate SARS-CoV-2 nonstructural proteins (nsps) and molecular dynamics (MD) simulations to investigate the effects of phosphorylation on nsps structure, which could be a potential inhibitory target to attenuate viral replication. Eight target candidate sites were found as top-ranked phosphorylation sites of SARS-CoV-2. During the process of molecular dynamics (MD) simulation, the root-mean-square deviation (RMSD) analysis was used to measure conformational changes in each nsps. Root-mean-square fluctuation (RMSF) was employed to measure the fluctuation in each residue of 36 systems considered, allowing us to evaluate the most flexible regions. These analysis shows that there are significant structural deviations in the residues namely nsp1 THR 72, nsp2 THR 73, nsp3 SER 64, nsp4 SER 81, nsp4 SER 455, nsp5 SER284, nsp6 THR 238, and nsp16 SER 132. The identified list of residues suggests how phosphorylation affects SARS-CoV-2 nsps function and stability. This research also suggests that kinase inhibitors could be a possible component for evaluating drug binding studies, which are crucial in therapeutic discovery research.

Calcium sparks are the elementary Ca 2+ release events in excitation-contraction coupling that underlie the Ca 2+ transient. The frequency-dependent contractile force generated by cardiac myocytes depends upon the characteristics of the Ca 2+ transients. A stochastic computational local control model of a guinea pig ventricular cardiomyocyte was developed, to gain insight into mechanisms of force-frequency relationship (FFR). This required the creation of a new three-state RyR2 model that reproduced the adaptive behavior of RyR2, in which the RyR2 channels transition into a different state when exposed to prolonged elevated subspace [Ca 2+ ]. The model simulations agree with previous experimental and modeling studies on interval-force relations. Unlike previous common pool models, this local control model displayed stable action potential trains at 7 Hz. The duration and the amplitude of the [Ca 2+ ] myo transients increase in pacing rates consistent with the experiments. The [Ca 2+ ] myo transient reaches its peak value at 4 Hz and decreases afterward, consistent with experimental force-frequency curves. The model predicts, in agreement with previous modeling studies of Jafri and co-workers, diastolic sarcoplasmic reticulum, [Ca 2+ ] sr , and RyR2 adaptation increase with the increased stimulation frequency, producing rising, rather than falling, amplitude of the myoplasmic [Ca 2+ ] transients. However, the local control model also suggests that the reduction of the L-type Ca 2+ current, with an increase in pacing frequency due to Ca 2+ -dependent inactivation, also plays a role in the negative slope of the FFR. In the simulations, the peak Ca 2+ transient in the FFR correlated with the highest numbers of SR Ca 2+ sparks: the larger average amplitudes of those sparks, and the longer duration of the Ca 2+ sparks.

Introduction and Objectives: In patients with coronavirus disease 2019 (COVID-19), several abnormal hematological biomarkers have been reported. The current study aimed to find out the association of neutrophil to lymphocyte ratio (NLR) and derived NLR (dNLR) with COVID-19. The objective was to compare the accuracy of both of these markers in predicting the severity of the disease., Materials and Methods: The study was conducted in a single-center having patients with COVID-19 with a considerable hospital stay. NLR is easily calculated by dividing the absolute neutrophil count (ANC) with the absolute lymphocyte count (ALC) {ANC/ALC}, while dNLR is calculated by ANC divided by total leukocyte count minus ANC {ANC/(WBC-ANC)}. Medians and interquartile ranges (IQR) were represented by box plots. Multivariable logistic regression was performed obtaining an odds ratio (OR), 95% CI, and further adjusted to discover the independent predictors and risk factors associated with elevated NLR and dNLR., Results: A total of 1,000 patients with COVID-19 were included. The baseline NLR and dNLR were 5.00 (2.91-10.46) and 4.00 (2.33-6.14), respectively. A cut-off value of 4.23 for NLR and 2.63 for dNLR were set by receiver operating characteristic (ROC) analysis. Significant associations of NLR were obtained by binary logistic regression for dependent outcome variables as ICU stay ( p < 0.001), death ( p < 0.001), and invasive ventilation ( p < 0.001) while that of dNLR with ICU stay ( p = 0.002), death ( p < 0.001), and invasive ventilation ( p = 0.002) on multivariate analysis when adjusted for age, gender, and a wave of pandemics. Moreover, the indices were found correlating with other inflammatory markers such as C-reactive protein (CRP), D-dimer, and procalcitonin (PCT)., Conclusion: Both markers are equally reliable and sensitive for predicting in-hospital outcomes of patients with COVID-19. Early detection and predictive analysis of these markers can allow physicians to risk assessment and prompt management of these patients., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Asghar, Akram, Yasmin, Najeeb, Naeem, Gaddam, Jafri, Tahir, Yasin, Mahmood, Mehmood and Marzo.)

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca 2+ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca 2+ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca 2+ sparks can combine to form Ca 2+ waves. These propagating elevations of (Ca 2+ ) i can activate an inward Na + -Ca 2+ exchanger current (I NCX ) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.

The mitochondrial permeability transition pore (mPTP) is a non-selective pore in the inner mitochondrial membrane (IMM) which causes depolarization when it opens under conditions of oxidative stress and high concentrations of Ca 2+ . In this study, a stochastic computational model was developed to better understand the dynamics of mPTP opening and closing associated with elevated reactive oxygen species (ROS) in cardiomyocytes. The data modeled are from "photon stress" experiments in which the fluorescent dye TMRM (tetramethylrhodamine methyl ester) is both the source of ROS (induced by laser light) and sensor of the electrical potential difference across the IMM. Monte Carlo methods were applied to describe opening and closing of the pore along with the Hill Equation to account for the effect of ROS levels on the transition probabilities. The amplitude distribution of transient mPTP opening events, the number of transient mPTP opening events per minute in a cell, the time it takes for recovery after transient depolarizations in the mitochondria, and the change in TMRM fluorescence during the transition from transient to permanent mPTP opening events were analyzed. The model suggests that mPTP transient open times have an exponential distribution that are reflected in TMRM fluorescence. A second multiple pore model in which individual channels have no permanent open state suggests that 5-10 mPTP per mitochondria would be needed for sustained mitochondrial depolarization at elevated ROS with at least 1 mPTP in the transient open state.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) nucleocapsid protein (N-protein) is responsible for viral replication by assisting in viral RNA synthesis and attaching the viral genome to the replicase-transcriptase complex (RTC). Numerous studies suggested the N-protein as a drug target. However, the specific N-protein active sites for SARS-CoV-2 drug treatments are yet to be discovered. The purpose of this study was to determine active sites of the SARS-CoV-2 N-protein by identifying torsion angle classifiers for N-protein structural changes that correlated with the respective angle differences between the active and inactive N-protein. In the study, classifiers with a minimum accuracy of 80% determined from molecular simulation data were analyzed by Principal Component Analysis and cross-validated by Logistic Regression, Support Vector Machine, and Random Forest Classification. The ability of torsion angles ψ252 and φ375 to differentiate between phosphorylated and unphosphorylated structures suggested that residues 252 and 375 in the RNA binding domain might be important in N-protein activation. Furthermore, the φ and ψ angles of residue S189 correlated to a 90.7% structural determination accuracy. The key residues involved in the structural changes identified here might suggest possible important functional sites on the N-protein that could be the focus of further study to understand their potential as drug targets., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Authors.)

Coronaviruses, including the recent pandemic strain SARS-Cov-2, use a multifunctional 2'-O-methyltransferase (2'-O-MTase) to restrict the host defense mechanism and to methylate RNA. The nonstructural protein 16 2'-O-MTase (nsp16) becomes active when nonstructural protein 10 (nsp10) and nsp16 interact. Novel peptide drugs have shown promise in the treatment of numerous diseases and new research has established that nsp10 derived peptides can disrupt viral methyltransferase activity via interaction of nsp16. This study had the goal of optimizing new analogous nsp10 peptides that have the ability to bind nsp16 with equal to or higher affinity than those naturally occurring. The following research demonstrates that in silico molecular simulations can shed light on peptide structures and predict the potential of new peptides to interrupt methyltransferase activity via the nsp10/nsp16 interface. The simulations suggest that misalignments at residues F68, H80, I81, D94, and Y96 or rotation at H80 abrogate MTase function. We develop a new set of peptides based on conserved regions of the nsp10 protein in the Coronaviridae species and test these to known MTase variant values. This results in the prediction that the H80R variant is a solid new candidate for potential new testing. We envision that this new lead is the beginning of a reputable foundation of a new computational method that combats coronaviruses and that is beneficial for new peptide drug development., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Authors.)

Mitochondria are complex organelles that use catabolic metabolism to produce ATP which is the critical energy source for cell function. Oxidative phosphorylation by the electron transport chain, which receives reducing equivalents (NADH and FADH 2 ) from the tricarboxylic acid cycle, also produces reactive oxygen species (ROS) as a by-product at complex I and III. ROS play a significant role in health and disease. In order to better understand this process, a computational model of mitochondrial energy metabolism and the production of ROS has been developed. The model demonstrates the process regulating ROS production and removal and how different energy substrates can affect ROS production., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Venetoclax is a BH3 (BCL-2 Homology 3) mimetic used to treat leukemia and lymphoma by inhibiting the anti-apoptotic BCL-2 protein thereby promoting apoptosis of cancerous cells. Acquired resistance to Venetoclax via specific variants in BCL-2 is a major problem for the successful treatment of cancer patients. Replica exchange molecular dynamics (REMD) simulations combined with machine learning were used to define the average structure of variants in aqueous solution to predict changes in drug and ligand binding in BCL-2 variants. The variant structures all show shifts in residue positions that occlude the binding groove, and these are the primary contributors to drug resistance. Correspondingly, we established a method that can predict the severity of a variant as measured by the inhibitory constant (K i ) of Venetoclax by measuring the structure deviations to the binding cleft. In addition, we also applied machine learning to the phi and psi angles of the amino acid backbone to the ensemble of conformations that demonstrated a generalizable method for drug resistant predictions of BCL-2 proteins that elucidates changes where detailed understanding of the structure-function relationship is less clear., (Copyright © 2021 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Calcium (Ca 2+ ) plays a central role in the excitation and contraction of cardiac myocytes. Experiments have indicated that calcium release is stochastic and regulated locally suggesting the possibility of spatially heterogeneous calcium levels in the cells. This spatial heterogeneity might be important in mediating different signaling pathways. During more than 50 years of computational cell biology, the computational models have been advanced to incorporate more ionic currents, going from deterministic models to stochastic models. While periodic increases in cytoplasmic Ca 2+ concentration drive cardiac contraction, aberrant Ca 2+ release can underly cardiac arrhythmia. However, the study of the spatial role of calcium ions has been limited due to the computational expense of using a three-dimensional stochastic computational model. In this paper, we introduce a three-dimensional stochastic computational model for rat ventricular myocytes at the whole-cell level that incorporate detailed calcium dynamics, with (1) non-uniform release site placement, (2) non-uniform membrane ionic currents and membrane buffers, (3) stochastic calcium-leak dynamics and (4) non-junctional or rogue ryanodine receptors. The model simulates spark-induced spark activation and spark-induced Ca 2+ wave initiation and propagation that occur under conditions of calcium overload at the closed-cell condition, but not when Ca 2+ levels are normal. This is considered important since the presence of Ca 2+ waves contribute to the activation of arrhythmogenic currents.

Dengue virus infection is a global health problem for which there have been challenges to obtaining a cure. Current vaccines and anti-viral drugs can only be narrowly applied in ongoing clinical trials. We employed computational methods based on structure-function relationships between human host kinases and viral nonstructural protein 3 (NS3) to understand viral replication inhibitors' therapeutic effect. Phosphorylation at each of the two most evolutionarily conserved sites of NS3, serine 137 and threonine 189, compared to the unphosphorylated state were studied with molecular dynamics and docking simulations. The simulations suggested that phosphorylation at serine 137 caused a more remarkable structural change than phosphorylation at threonine 189, specifically located at amino acid residues 49-95. Docking studies supported the idea that phosphorylation at serine 137 increased the binding affinity between NS3 and nonstructural Protein 5 (NS5), whereas phosphorylation at threonine 189 decreased it. The interaction between NS3 and NS5 is essential for viral replication. Docking studies with the antiviral plant flavonoid Quercetin with NS3 indicated that Quercetin physically occluded the serine 137 phosphorylation site. Taken together, these findings suggested a specific site and mechanism by which Quercetin inhibits dengue and possible other flaviviruses.

Alzheimer's disease, the most common form of dementia, currently has no cure. There are only temporary treatments that reduce symptoms and the progression of the disease. Alzheimer's disease is characterized by the prevalence of plaques of aggregated amyloid β (Aβ) peptide. Recent treatments to prevent plaque formation have provided little to relieve disease symptoms. Although there have been numerous molecular simulation studies on the mechanisms of Aβ aggregation, the signaling role has been less studied. In this study, a total of over 38,000 simulated structures, generated from molecular dynamics (MD) simulations, exploring different conformations of the Aβ42 mutants and wild-type peptides were used to examine the relationship between Aβ torsion angles and disease measures. Unique methods characterized the data set and pinpointed residues that were associated in aggregation and others associated with signaling. Machine learning techniques were applied to characterize the molecular simulation data and classify how much each residue influenced the predicted variant of Alzheimer's Disease. Orange3 data mining software provided the ability to use these techniques to generate tables and rank the data. The test and score module coupled with the confusion matrix module analyzed data with calculations of specificity and sensitivity. These methods evaluating frequency and rank allowed us to analyze and predict important residues associated with different phenotypic measures. This research has the potential to help understand which specific residues of Aβ should be targeted for drug development.

Cardiac alternans is characterized by alternating weak and strong beats of the heart. This signaling at the cellular level may appear as alternating long and short action potentials (APs) that occur in synchrony with alternating large and small calcium transients, respectively. Previous studies have suggested that alternans manifests itself through either a voltage dependent mechanism based upon action potential restitution or as a calcium dependent mechanism based on refractoriness of calcium release. We use a novel model of cardiac excitation-contraction (EC) coupling in the rat ventricular myocyte that includes 20,000 calcium release units (CRU) each with 49 ryanodine receptors (RyR2s) and 7 L-type calcium channels that are all stochastically gated. The model suggests that at the cellular level in the case of alternans produced by rapid pacing, the mechanism requires a synergy of voltage- and calcium-dependent mechanisms. The rapid pacing reduces AP duration and magnitude reducing the number of L-type calcium channels activating individual CRUs during each AP and thus increases the population of CRUs that can be recruited stochastically. Elevated myoplasmic and sarcoplasmic reticulum (SR) calcium, [Ca 2+ ] myo and [Ca 2+ ] SR respectively, increases ryanodine receptor open probability (P o ) according to our model used in this simulation and this increased the probability of activating additional CRUs. A CRU that opens in one beat is less likely to open the subsequent beat due to refractoriness caused by incomplete refilling of the junctional sarcoplasmic reticulum (jSR). Furthermore, the model includes estimates of changes in Na + fluxes and [Na + ] i and thus provides insight into how changes in electrical activity, [Na + ] i and sodium-calcium exchanger activity can modulate alternans. The model thus tracks critical elements that can account for rate-dependent changes in [Na + ] i and [Ca 2+ ] myo and how they contribute to the generation of Ca 2+ signaling alternans in the heart.

The stretching of a cardiomyocyte leads to the increased production of reactive oxygen species that increases ryanodine receptor open probability through a process termed X-ROS signaling. The stretching of the myocyte also increases the calcium affinity of myofilament Troponin C, which increases its calcium buffering capacity. Here, an integrative experimental and modeling study is pursued to explain the interplay of length-dependent changes in calcium buffering by troponin and stretch-activated X-ROS calcium signaling. Using this combination, we show that the troponin C-dependent increase in myoplasmic calcium buffering during myocyte stretching largely offsets the X-ROS-dependent increase in calcium release from the sarcoplasmic reticulum. The combination of modeling and experiment are further informed by the elimination of length-dependent changes to troponin C calcium binding in the presence of blebbistatin. Here, the model suggests that it is the X-ROS signaling-dependent Ca 2+ release increase that serves to maintain free myoplasmic calcium concentrations during a change in myocyte length. Together, our experimental and modeling approaches have further defined the relative contributions of X-ROS signaling and the length-dependent calcium buffering by troponin in shaping the myoplasmic calcium transient.

Biological processes are incredibly complex-integrating molecular signaling networks involved in multicellular communication and function, thus maintaining homeostasis. Dysfunction of these processes can result in the disruption of homeostasis, leading to the development of several disease processes including atherosclerosis. We have significantly advanced our understanding of bioprocesses in atherosclerosis, and in doing so, we are beginning to appreciate the complexities, intricacies, and heterogeneity atherosclerosi. We are also now better equipped to acquire, store, and process the vast amount of biological data needed to shed light on the biological circuitry involved. Such data can be analyzed within machine learning frameworks to better tease out such complex relationships. Indeed, there has been an increasing number of studies applying machine learning methods for patient risk stratification based on comorbidities, multi-modality image processing, and biomarker discovery pertaining to atherosclerotic plaque formation. Here, we focus on current applications of machine learning to provide insight into atherosclerotic plaque formation and better understand atherosclerotic plaque progression in patients with cardiovascular disease., (© 2021 Author(s).)

Folding of the mitochondrial inner membrane (IM) into cristae greatly increases the ATP-generating surface area, S , per unit volume but also creates diffusional bottlenecks that could limit reaction rates inside mitochondria. This study explores possible effects of inner membrane folding on mitochondrial ATP output, using a mathematical model for energy metabolism developed by the Jafri group and two- and three-dimensional spatial models for mitochondria, implemented on the Virtual Cell platform. Simulations demonstrate that cristae are micro-compartments functionally distinct from the cytosol. At physiological steady states, standing gradients of ADP form inside cristae that depend on the size and shape of the compartments, and reduce local flux (rate per unit area) of the adenine nucleotide translocase. This causes matrix ADP levels to drop, which in turn reduces the flux of ATP synthase. The adverse effects of membrane folding on reaction fluxes increase with crista length and are greater for lamellar than tubular crista. However, total ATP output per mitochondrion is the product of flux of ATP synthase and IM , per unit volume but also creates diffusional bottlenecks that could limit reaction rates inside mitochondria. This study explores possible effects of inner membrane folding on mitochondrial ATP output, using a mathematical model for energy metabolism developed by the Jafri group and two- and three-dimensional spatial models for mitochondria, implemented on the Virtual Cell platform. Simulations demonstrate that cristae are micro-compartments functionally distinct from the cytosol. At physiological steady states, standing gradients of ADP form inside cristae that depend on the size and shape of the compartments, and reduce local flux (rate per unit area) of the adenine nucleotide translocase. This causes matrix ADP levels to drop, which in turn reduces the flux of ATP synthase. The adverse effects of membrane folding on reaction fluxes increase with crista length and are greater for lamellar than tubular crista. However, total ATP output per mitochondrion is the product of flux of ATP synthase and S IM which can be two-fold greater for mitochondria with lamellar than tubular cristae, resulting in greater ATP output for the former. The simulations also demonstrate the crucial role played by intracristal kinases (adenylate kinase, creatine kinase) in maintaining the energy advantage of IM folding., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Gain-of-function RyR1-p.R163C mutation in ryanodine receptors type 1 (RyR1) deregulates Ca 2+ signaling and mitochondrial function in skeletal muscle and causes malignant hyperthermia in humans and mice under triggering conditions. We investigated whether T lymphocytes from heterozygous RyR1-p.R163C knock-in mutant mice (HET T cells) display measurable aberrations in resting cytosolic Ca 2+ concentration ([Ca 2+ ] i ), Ca 2+ release from the store, store-operated Ca 2+ entry (SOCE), and mitochondrial inner membrane potential (ΔΨ m ) compared with T lymphocytes from wild-type mice (WT T cells). We explored whether these variables can be used to distinguish between T cells with normal and altered RyR1 genotype. HET and WT T cells were isolated from spleen and lymph nodes and activated in vitro using phytohemagglutinin P. [Ca 2+ ] i and ΔΨ m dynamics were examined using Fura 2 and tetramethylrhodamine methyl ester fluorescent dyes, respectively. Activated HET T cells displayed elevated resting [Ca 2+ ] i , diminished responses to Ca 2+ mobilization with thapsigargin, and decreased rate of [Ca 2+ ] i elevation in response to SOCE compared with WT T cells. Pretreatment of HET T cells with ryanodine or dantrolene sodium reduced disparities in the resting [Ca 2+ ] i and ability of thapsigargin to mobilize Ca 2+ between HET and WT T cells. While SOCE elicited dissipation of the ΔΨ m in WT T cells, it produced ΔΨ m hyperpolarization in HET T cells. When used as the classification variable, the amplitude of thapsigargin-induced Ca 2+ transient showed the best promise in predicting the presence of RyR1-p.R163C mutation. Other significant variables identified by machine learning analysis were the ratio of resting cytosolic Ca 2+ level to the amplitude of thapsigargin-induced Ca 2+ transient and an integral of changes in ΔΨ m in response to SOCE. Our study demonstrated that gain-of-function mutation in RyR1 significantly affects Ca 2+ signaling and mitochondrial fiction in T lymphocytes, which suggests that this mutation may cause altered immune responses in its carrier. Our data link the RyR1-p.R163C mutation, which causes inherited skeletal muscle diseases, to deregulation of Ca 2+ signaling and mitochondrial function in immune T cells and establish proof-of-principle for in vitro T cell-based diagnostic assay for hereditary RyR1 hyperfunction., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Neuropeptide S (NPS) is a naturally occurring appetite stimulant, associated with anxiety, stress, and excitement regulation. Neuropeptide S serves as a hypothalamic energy regulator that enhances food intake with a reduced level of satiety. NPS activates fat angiogenesis and the proliferation of new adipocytes in obesity. NPS has an established role in energy regulation by many pre-clinical investigations; however we have limited data available to support this notion in humans. We found significant association of Neuropeptide S receptor (NPSR1) Asn107Ile (rs324981, A>T) polymorphism with obese male participants. The current investigation carried out genotype screening of NPSR1 allele to assess the spectrum of the Asn107Ile polymorphism in obese and healthy Pakistani individuals. We revealed a significant (p = 0.04) difference between AA vs TT + AT genotype distribution of NPSR1 (SNP rs324981,) between obese and healthy individuals (p = 0.04). In this genotype analysis of (SNP rs324981) of the NPSR1 gene, T allele was marked as risk allele with higher frequency in the obese (38%) compared to its frequency in the controls (25%). Single Nucleotide Polymorphism (SNP, rs324981) Asn107Ile of NPSR1gene, that switches an amino acid from Asn to Ile, has been found associated with increased susceptibility to obesity in Pakistani individuals. Furthermore, molecular simulation studies predicted a lower binding affinity of NPSR1 Asn107Ile variant to NPS than the wild-type consistent with the genotype studies. These molecular simulation studies predict a possible molecular mechanism of this interaction by defining the key amino acid residues. However, a significantly (p<0.0001) lower concentration of NPS was recorded independent of genotype frequencies in obese subjects compared to healthy controls. We believe that large scale polymorphism data of population for important gene players including NPSR1 will be more useful to understand obesity and its associated risk factors., Competing Interests: The authors declare no conflict of interest.

Background: Psoriasis is associated with elevated risk of heart attack and increased accumulation of subclinical noncalcified coronary burden by coronary computed tomography angiography (CCTA). Machine learning algorithms have been shown to effectively analyze well-characterized data sets., Objective: In this study, we used machine learning algorithms to determine the top predictors of noncalcified coronary burden by CCTA in psoriasis., Methods: The analysis included 263 consecutive patients with 63 available variables from the Psoriasis Atherosclerosis Cardiometabolic Initiative. The random forest algorithm was used to determine the top predictors of noncalcified coronary burden by CCTA. We evaluated our results using linear regression models., Results: Using the random forest algorithm, we found that the top 10 predictors of noncalcified coronary burden were body mass index, visceral adiposity, total adiposity, apolipoprotein A1, high-density lipoprotein, erythrocyte sedimentation rate, subcutaneous adiposity, small low-density lipoprotein particle, cholesterol efflux capacity and the absolute granulocyte count. Linear regression of noncalcified coronary burden yielded results consistent with our machine learning output., Limitation: We were unable to provide external validation and did not study cardiovascular events., Conclusion: Machine learning methods identified the top predictors of noncalcified coronary burden in psoriasis. These factors were related to obesity, dyslipidemia, and inflammation, showing that these are important targets when treating comorbidities in psoriasis., (Published by Elsevier Inc.)

Long non-coding RNA Knowledgebase (lncRNAKB) is an integrated resource for exploring lncRNA biology in the context of tissue-specificity and disease association. A systematic integration of annotations from six independent databases resulted in 77,199 human lncRNA (224,286 transcripts). The user-friendly knowledgebase covers a comprehensive breadth and depth of lncRNA annotation. lncRNAKB is a compendium of expression patterns, derived from analysis of RNA-seq data in thousands of samples across 31 solid human normal tissues (GTEx). Thousands of co-expression modules identified via network analysis and pathway enrichment to delineate lncRNA function are also accessible. Millions of expression quantitative trait loci (cis-eQTL) computed using whole genome sequence genotype data (GTEx) can be downloaded at lncRNAKB that also includes tissue-specificity, phylogenetic conservation and coding potential scores. Tissue-specific lncRNA-trait associations encompassing 323 GWAS (UK Biobank) are also provided. LncRNAKB is accessible at http://www.lncrnakb.org/ , and the data are freely available through Open Science Framework ( https://doi.org/10.17605/OSF.IO/RU4D2 ).

T cell receptor stimulation initiates a cascade of reactions that cause an increase in intracellular calcium (Ca 2+ ) concentration mediated through inositol 1,4,5-trisphosphate (IP 3 ). To understand the basic mechanisms by which the immune response in T cells is activated, it is useful to understand the signaling pathways that contain important targets for drugs in a quantitative fashion. A computational model helps us to understand how the selected elements in the pathways interact with each other, and which component plays the crucial role in systems. We have developed a mathematical model to explore the mechanism for controlling transcription factor activity, which regulates gene expression, by the modulation of calcium signaling triggered during T cell activation. The model simulates the activation and modulation of Ca 2+ release-activated Ca 2+ (CRAC) channels by mitochondrial dynamics and depletion of endoplasmic reticulum (ER) store, and also includes membrane potential in T-cells. The model simulates the experimental finding that increases in Ca 2+ current enhances the activation of transcription factors and the Ca 2+ influx through CRAC is also essential for the NFAT and NFκB activation. The model also suggests that plasma membrane Ca 2+ -ATPase (PMCA) controls a majority of the extrusion of Ca 2+ and modulates the activation of CRAC channels. Furthermore, the model simulations explain how the complex interaction of the endoplasmic reticulum, membrane potential, mitochondria, and ion channels such as CRAC channels control T cell activation., (© 2020 Published by Elsevier Ltd.)

Background: Molecular simulations are used to provide insight into protein structure and dynamics, and have the potential to provide important context when predicting the impact of sequence variation on protein function. In addition to understanding molecular mechanisms and interactions on the atomic scale, translational applications of those approaches include drug screening, development of novel molecular therapies, and targeted treatment planning. Supporting the continued development of these applications, we have developed the SNP2SIM workflow that generates reproducible molecular dynamics and molecular docking simulations for downstream functional variant analysis. The Python workflow utilizes molecular dynamics software (NAMD (Phillips et al., J Comput Chem 26(16):1781-802, 2005), VMD (Humphrey et al., J Mol Graph 14(1):33-8, 27-8, 1996)) to generate variant specific scaffolds for simulated small molecule docking (AutoDock Vina (Trott and Olson, J Comput Chem 31(2):455-61, 2010))., Results: SNP2SIM is composed of three independent modules that can be used sequentially to generate the variant scaffolds of missense protein variants from the wildtype protein structure. The workflow first generates the mutant structure and configuration files required to execute molecular dynamics simulations of solvated protein variant structures. The resulting trajectories are clustered based on the structural diversity of residues involved in ligand binding to produce one or more variant scaffolds of the protein structure. Finally, these unique structural conformations are bound to small molecule ligand libraries to predict variant induced changes to drug binding relative to the wildtype protein structure., Conclusions: SNP2SIM provides a platform to apply molecular simulation based functional analysis of sequence variation in the protein targets of small molecule therapies. In addition to simplifying the simulation of variant specific drug interactions, the workflow enables large scale computational mutagenesis by controlling the parameterization of molecular simulations across multiple users or distributed computing infrastructures. This enables the parallelization of the computationally intensive molecular simulations to be aggregated for downstream functional analysis, and facilitates comparing various simulation options, such as the specific residues used to define structural variant clusters. The Python scripts that implement the SNP2SIM workflow are available (SNP2SIM Repository. https://github.com/mccoymd/SNP2SIM , Accessed 2019 February ), and individual SNP2SIM modules are available as apps on the Seven Bridges Cancer Genomics Cloud (Lau et al., Cancer Res 77(21):e3-e6, 2017; Cancer Genomics Cloud [ www.cancergenomicscloud.org ; Accessed 2018 November]).

The pharmacological targeting of the α 7 nicotinic acetylcholine receptor ( α 7) is a promising strategy in the development of new drugs for neurologic diseases. Because α 7 receptors regulate cellular calcium, we investigated how the prototypical type II-positive allosteric modulator PNU120596 affects α 7-mediated calcium signaling. Live imaging experiments show that PNU120596 augments ryanodine receptor-driven calcium-induced calcium release (CICR), inositol-induced calcium release (IICR), and phospholipase C activation by the α 7 receptor. Both influx of calcium through the α 7 nicotinic acetylcholine receptor (nAChR) channel as well as the binding of intracellular G proteins were involved in the effect of PNU120596 on intracellular calcium. This is evidenced by the findings that chelation of extracellular calcium, expression of α 7 D44A or α 7 345-348A mutant subunits, or blockade of calcium store release compromised the ability of PNU120596 to increase intracellular calcium transients generated by α 7 ligand activation. Spatiotemporal stochastic modeling of calcium transient responses corroborates these results and indicates that α 7 receptor activation enables calcium microdomains locally and to lesser extent in the distant cytosol. From the model, allosteric modulation of the receptor activates CICR locally via ryanodine receptors and augments IICR through enhanced calcium influx due to prolonged α 7 nAChR opening. These findings provide a new mechanistic framework for understanding the effect of α 7 receptor allosteric modulation on both local and global calcium dynamics., (Copyright © 2018 by The American Society for Pharmacology and Experimental Therapeutics.)

Calcium-induced calcium release is the principal mechanism that triggers the cell-wide [Ca(2+)]i transient that activates muscle contraction during cardiac excitation-contraction coupling (ECC). Here, we characterize this process in mouse cardiac myocytes with a novel mathematical action potential (AP) model that incorporates realistic stochastic gating of voltage-dependent L-type calcium (Ca(2+)) channels (LCCs) and sarcoplasmic reticulum (SR) Ca(2+) release channels (the ryanodine receptors, RyR2s). Depolarization of the sarcolemma during an AP stochastically activates the LCCs elevating subspace [Ca(2+)] within each of the cell's 20,000 independent calcium release units (CRUs) to trigger local RyR2 opening and initiate Ca(2+) sparks, the fundamental unit of triggered Ca(2+) release. Synchronization of Ca(2+) sparks during systole depends on the nearly uniform cellular activation of LCCs and the likelihood of local LCC openings triggering local Ca(2+) sparks (ECC fidelity). The detailed design and true SR Ca(2+) pump/leak balance displayed by our model permits investigation of ECC fidelity and Ca(2+) spark fidelity, the balance between visible (Ca(2+) spark) and invisible (Ca(2+) quark/sub-spark) SR Ca(2+) release events. Excess SR Ca(2+) leak is examined as a disease mechanism in the context of "catecholaminergic polymorphic ventricular tachycardia (CPVT)", a Ca(2+)-dependent arrhythmia. We find that that RyR2s (and therefore Ca(2+) sparks) are relatively insensitive to LCC openings across a wide range of membrane potentials; and that key differences exist between Ca(2+) sparks evoked during quiescence, diastole, and systole. The enhanced RyR2 [Ca(2+)]i sensitivity during CPVT leads to increased Ca(2+) spark fidelity resulting in asynchronous systolic Ca(2+) spark activity. It also produces increased diastolic SR Ca(2+) leak with some prolonged Ca(2+) sparks that at times become "metastable" and fail to efficiently terminate. There is a huge margin of safety for stable Ca(2+) handling within the cell and this novel mechanistic model provides insight into the molecular signaling characteristics that help maintain overall Ca(2+) stability even under the conditions of high SR Ca(2+) leak during CPVT. Finally, this model should provide tools for investigators to examine normal and pathological Ca(2+) signaling characteristics in the heart., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Aims: In atrial fibrillation (AF), abnormalities in Ca(2+) release contribute to arrhythmia generation and contractile dysfunction. We explore whether ryanodine receptor (RyR) cluster ultrastructure is altered and is associated with functional abnormalities in AF., Methods and Results: Using high-resolution confocal microscopy (STED), we examined RyR cluster morphology in fixed atrial myocytes from sheep with persistent AF (N = 6) and control (Ctrl; N = 6) animals. RyR clusters on average contained 15 contiguous RyRs; this did not differ between AF and Ctrl. However, the distance between clusters was significantly reduced in AF (288 ± 12 vs. 376 ± 17 nm). When RyR clusters were grouped into Ca(2+) release units (CRUs), i.e. clusters separated by <150 nm, CRUs in AF had more clusters (3.43 ± 0.10 vs. 2.95 ± 0.02 in Ctrl), which were more dispersed. Furthermore, in AF cells, more RyR clusters were found between Z lines. In parallel experiments, Ca(2+) sparks were monitored in live permeabilized myocytes. In AF, myocytes had >50% higher spark frequency with increased spark time to peak (TTP) and duration, and a higher incidence of macrosparks. A computational model of the CRU was used to simulate the morphological alterations observed in AF cells. Increasing cluster fragmentation to the level observed in AF cells caused the observed changes, i.e. higher spark frequency, increased TTP and duration; RyR clusters dispersed between Z-lines increased the occurrence of macrosparks., Conclusion: In persistent AF, ultrastructural reorganization of RyR clusters within CRUs is associated with overactive Ca(2+) release, increasing the likelihood of propagating Ca(2+) release., (© The Author 2015. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Activated phospholipase C (PLC*) generates 1,4,5-triphosphate (IP 3 ) and diacylglycerol (DAG) from phosphatidyl inositol (PIP 2 ). The DAG remains in the plasma membrane and co-activates conventional protein kinase C (PKC) with Ca 2+ . We have developed a mathematical model for the activation of the Ca 2+ -dependent PKC and its negative feedback on phospholipase C (PLC) and coupled it to the De Young-Keizer model for IP 3 mediated Ca 2+ oscillations. The model describes the cascade of reactions for the translocation of PKC to plasma membrane, and simulates activation of Ca 2+ and diacylglycerol (DAG) oscillations. The model demonstrates that oscillations in Ca 2+ and DAG are possible with or without a positive Ca 2+ feedback on phospholipase C consistent with experiment. In many experimental studies, the timing of the peaks of the Ca 2+ and IP 3 oscillations have been used to suggest causality, i.e. that the IP 3 oscillations cause the Ca 2+ oscillations. The model is used to explore this question. To this end, the positive and negative feedback between Ca 2+ and IP 3 production are modulated, resulting in changes to the phase lag between the peaks in [Ca 2+ ] cyt and [IP] cyt . The model simulates a possible experimental protocol that can be used to differentiate whether or not the positive feedback of Ca 2+ on PLC is needed for the oscillations.

Stretching single ventricular cardiac myocytes has been shown experimentally to activate transmembrane nicotinamide adenine dinucleotide phosphate oxidase type 2 to produce reactive oxygen species (ROS) and increase the Ca2+ spark rate in a process called X-ROS signaling. The increase in Ca2+ spark rate is thought to be due to an increase in ryanodine receptor type 2 (RyR2) open probability by direct oxidation of the RyR2 protein complex. In this article, a computational model is used to examine the regulation of ROS and calcium homeostasis by local, subcellular X-ROS signaling and its role in cardiac excitation-contraction coupling. To this end, a four-state RyR2 model was developed that includes an X-ROS-dependent RyR2 mode switch. When activated, [Ca2+]i-sensitive RyR2 open probability increases, and the Ca2+ spark rate changes in a manner consistent with experimental observations. This, to our knowledge, new model is used to study the transient effects of diastolic stretching and subsequent ROS production on RyR2 open probability, Ca2+ sparks, and the myoplasmic calcium concentration ([Ca2+]i) during excitation-contraction coupling. The model yields several predictions: 1) [ROS] is produced locally near the RyR2 complex during X-ROS signaling and increases by an order of magnitude more than the global ROS signal during myocyte stretching; 2) X-ROS activation just before the action potential, corresponding to ventricular filling during diastole, increases the magnitude of the Ca2+ transient; 3) during prolonged stretching, the X-ROS-induced increase in Ca2+ spark rate is transient, so that long-sustained stretching does not significantly increase sarcoplasmic reticulum Ca2+ leak; and 4) when the chemical reducing capacity of the cell is decreased, activation of X-ROS signaling increases sarcoplasmic reticulum Ca2+ leak and contributes to global oxidative stress, thereby increases the possibility of arrhythmia. The model provides quantitative information not currently obtainable through experimental means and thus provides a framework for future X-ROS signaling experiments., (Copyright © 2015 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Stable calcium-induced calcium release (CICR) is critical for maintaining normal cellular contraction during cardiac excitation-contraction coupling. The fundamental element of CICR in the heart is the calcium (Ca(2+)) spark, which arises from a cluster of ryanodine receptors (RyR). Opening of these RyR clusters is triggered to produce a local, regenerative release of Ca(2+) from the sarcoplasmic reticulum (SR). The Ca(2+) leak out of the SR is an important process for cellular Ca(2+) management, and it is critically influenced by spark fidelity, i.e., the probability that a spontaneous RyR opening triggers a Ca(2+) spark. Here, we present a detailed, three-dimensional model of a cardiac Ca(2+) release unit that incorporates diffusion, intracellular buffering systems, and stochastically gated ion channels. The model exhibits realistic Ca(2+) sparks and robust Ca(2+) spark termination across a wide range of geometries and conditions. Furthermore, the model captures the details of Ca(2+) spark and nonspark-based SR Ca(2+) leak, and it produces normal excitation-contraction coupling gain. We show that SR luminal Ca(2+)-dependent regulation of the RyR is not critical for spark termination, but it can explain the exponential rise in the SR Ca(2+) leak-load relationship demonstrated in previous experimental work. Perturbations to subspace dimensions, which have been observed in experimental models of disease, strongly alter Ca(2+) spark dynamics. In addition, we find that the structure of RyR clusters also influences Ca(2+) release properties due to variations in inter-RyR coupling via local subspace Ca(2+) concentration ([Ca(2+)]ss). These results are illustrated for RyR clusters based on super-resolution stimulated emission depletion microscopy. Finally, we present a believed-novel approach by which the spark fidelity of a RyR cluster can be predicted from structural information of the cluster using the maximum eigenvalue of its adjacency matrix. These results provide critical insights into CICR dynamics in heart, under normal and pathological conditions., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Wilbrand and Saenger(1) studied optic chiasms after unilateral enucleation, noting inferonasal crossing fibers curved anteriorly into the contralateral optic nerve (Wilbrand knee; figure, A). This explains contralateral superotemporal visual field defects (junctional scotomas) with optic nerve lesions at the chiasmal junction. However, Wilbrand knee may be an enucleation artifact.(2) The anisotropic light-reflecting properties of myelinated axons permitted imaging of normal human chiasms. Thin sections (25 µm) were illuminated and digitally imaged from 3 incident angles. Each of the images was pseudocolored (red, green, or blue) and merged, revealing an anomalously oriented fiber tract (appearing white) that reversed direction at the optic nerve-chiasm junction, found in inferior (figure, C) but not in superior sections (figure, B), consistent with Wilbrand and Saenger's original description.

Many neurodegenerative diseases involve defects in mitochondrial function. These defects often arise from mutations carried in the genes that code mitochondrial proteins. Many of these defective proteins are involved in mitochondrial energy metabolism which is one of the primary roles of the mitochondria. However, others proteins have other roles that are related to signaling or mitochondrial structure. The interaction of the mitochondrial proteins is complex. Understanding these complex dynamics requires the use of computational models. Studies have started to exploit such models, but much further work is necessary to understand mitochondrial function and its role in neurodegenerative disease., (© 2014 Elsevier Inc. All rights reserved.)

Myofascial pain syndrome is an important health problem. It affects a majority of the general population, impairs mobility, causes pain, and reduces the overall sense of well-being. Underlying this syndrome is the existence of painful taut bands of muscle that contain discrete, hypersensitive foci called myofascial trigger points. In spite of the significant impact on public health, a clear mechanistic understanding of the disorder does not exist. This is likely due to the complex nature of the disorder which involves the integration of cellular signaling, excitation-contraction coupling, neuromuscular inputs, local circulation, and energy metabolism. The difficulties are further exacerbated by the lack of an animal model for myofascial pain to test mechanistic hypothesis. In this review, current theories for myofascial pain are presented and their relative strengths and weaknesses are discussed. Based on new findings linking mechanoactivation of reactive oxygen species signaling to destabilized calcium signaling, we put forth a novel mechanistic hypothesis for the initiation and maintenance of myofascial trigger points. It is hoped that this lays a new foundation for understanding myofascial pain syndrome and how current therapies work, and gives key insights that will lead to the improvement of therapies for its treatment.

Introduction: Pegylated interferon and ribavirin (PEGIFN/RBV) therapy for recurrent hepatitis C after liver transplantation (LT) is associated with a lower sustained virological response (SVR) rate as well as more frequent side effects compared to non-transplant patients. We aimed to determine the incidence and clinical characteristics of LT recipients with recurrent hepatitis C who developed immunological dysfunction (ID) during or after PEG-IFN/RBV therapy and to assess its impact on patient and graft survival., Methods: Seventy-four deceased donor LT recipients with histological recurrence of hepatitis C were treated with PEG-IFN/RBV from 1/00 to 12/08. ID was defined as biopsy-proven rejection or moderate plasma cell hepatitis. Patients were followed up until death, re-LT or 30 September 2011., Results: Twelve patients (16 %) had ID, 8 (10.7 %) had cholestasis without ID, while 54 had no ID/cholestasis during or after discontinuation of PEG-IFN/RBV therapy. Biopsy-proven acute cellular rejection prior to (hazard ratio = 4.87, p = 0.009) and type of immunosuppression at the time of initiation of PEG-IFN/RBV were the only independent predictors of ID. Patients who were on tacrolimus at the time of initiation of PEG-IFN/RBV had a significantly lower risk of ID compared to those who were on cyclosporine (HR 0.254, p = 0.023). Patients with ID had a trend toward a lower SVR rate (25 vs. 54 %, p = 0.18) and a significantly higher rate of graft failure (33 vs. 4 %, p = 0.004) compared to patients with no ID/cholestasis., Conclusions: ID is common during or after PEG-IFN/RBV therapy for recurrent hepatitis C and frequently associated with decreased graft survival, trending toward low rates of SVR. Careful monitoring of liver biochemistries during or after PEG-IFN/RBV therapy with a low threshold to biopsy patients and particularly those receiving cyclosporine-based immunosuppression may improve outcomes in these patients.

The mitochondrial inner membrane has a complex and dynamic structure that plays an important role in the function of this organelle. The internal compartments called cristae are created by processes that are just beginning to be understood. Crista size and morphology influence the internal diffusion of solutes and the surface area of the inner membrane, which is home to critical membrane proteins including ATP synthase and electron transport chain complexes; metabolite and ion transporters including the adenine nucleotide translocase, the calcium uniporter (MCU), and the sodium/calcium exchanger (NCLX); and many more. Here we provide a brief overview of what is known about crista structure and formation, and discuss mitochondrial function in the context of that structure. We also suggest that mathematical modeling of mitochondria that incorporates accurate information about the organelle's internal architecture can lead to a better understanding of its diverse functions. This article is part of a Special Issue entitled 'Calcium Signalling in Heart'., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Miniature optical sensors that can detect blood vessels in front of advancing instruments will significantly benefit many interventional procedures. Towards this end, we developed a thin and flexible coherence-gated Doppler (CGD) fiber probe (O.D. = 0.125 mm) that can be integrated with minimally-invasive tools to provide real-time audio feedback of blood flow at precise locations in front of the probe. Coherence-gated Doppler (CGD) is a hybrid technology with features of laser Doppler flowmetry (LDF) and Doppler optical coherence tomography (DOCT). Because of its confocal optical design and coherence-gating capabilities, CGD provides higher spatial resolution than LDF. And compared to DOCT imaging systems, CGD is simpler and less costly to produce. In vivo studies of rat femoral vessels using CGD demonstrate its ability to distinguish between artery, vein and bulk movement of the surrounding soft tissue. Finally, by placing the CGD probe inside a 30-gauge needle and advancing it into the brain of an anesthetized sheep, we demonstrate that it is capable of detecting vessels in front of advancing probes during simulated stereotactic neurosurgical procedures. Using simultaneous ultrasound (US) monitoring from the surface of the brain we show that CGD can detect at-risk blood vessels up to 3 mm in front of the advancing probe. The improved spatial resolution afforded by coherence gating combined with the simplicity, minute size and robustness of the CGD probe suggest it may benefit many minimally invasive procedures and enable it to be embedded into a variety of surgical instruments.

Rationale: Transverse tubules (TTs) couple electric surface signals to remote intracellular Ca(2+) release units (CRUs). Diffraction-limited imaging studies have proposed loss of TT components as disease mechanism in heart failure (HF)., Objectives: Objectives were to develop quantitative super-resolution strategies for live-cell imaging of TT membranes in intact cardiomyocytes and to show that TT structures are progressively remodeled during HF development, causing early CRU dysfunction., Methods and Results: Using stimulated emission depletion (STED) microscopy, we characterized individual TTs with nanometric resolution as direct readout of local membrane morphology 4 and 8 weeks after myocardial infarction (4pMI and 8pMI). Both individual and network TT properties were investigated by quantitative image analysis. The mean area of TT cross sections increased progressively from 4pMI to 8pMI. Unexpectedly, intact TT networks showed differential changes. Longitudinal and oblique TTs were significantly increased at 4pMI, whereas transversal components appeared decreased. Expression of TT-associated proteins junctophilin-2 and caveolin-3 was significantly changed, correlating with network component remodeling. Computational modeling of spatial changes in HF through heterogeneous TT reorganization and RyR2 orphaning (5000 of 20 000 CRUs) uncovered a local mechanism of delayed subcellular Ca(2+) release and action potential prolongation., Conclusions: This study introduces STED nanoscopy for live mapping of TT membrane structures. During early HF development, the local TT morphology and associated proteins were significantly altered, leading to differential network remodeling and Ca(2+) release dyssynchrony. Our data suggest that TT remodeling during HF development involves proliferative membrane changes, early excitation-contraction uncoupling, and network fracturing.

Excitation-contraction coupling describes the processes relating to electrical excitation through force generation and contraction in the heart. It occurs at multiple levels from the whole heart, to single myocytes and down to the sarcomere. A central process that links electrical excitation to contraction is calcium mobilization. Computational models that are well grounded in experimental data have been an effective tool to understand the complex dynamics of the processes involved in excitation-contraction coupling. Presented here is a summary of some computational models that have added to the understanding of the cellular and subcellular mechanisms that control ventricular myocyte calcium dynamics. Models of cardiac ventricular myocytes that have given insight into termination of calcium release and interval-force relations are discussed in this manuscript. Computational modeling of calcium sparks, the elementary events in cardiac excitation-contraction coupling, has given insight into mechanism governing their dynamics and termination as well as their role in excitation-contraction coupling and is described herein.

A forward-imaging needle-type optical coherence tomography (OCT) probe with Doppler OCT (DOCT) capability has the potential to solve critical challenges in interventional procedures. A case in point is stereotactic neurosurgery where probes are advanced into the brain based on predetermined coordinates. Laceration of blood vessels in front of the advancing probe is an unavoidable complication with current methods. Moreover, cerebrospinal fluid (CSF) leakage during surgery can shift the brain rendering the predetermined coordinates unreliable. In order to address these challenges, we developed a forward-imaging OCT probe (740 μm O.D.) using a gradient-index (GRIN) rod lens that can provide real-time imaging feedback for avoiding at-risk vessels (8 frames/s with 1024 A-scans per frame for OCT/DOCT dual imaging) and guiding the instrument to specific targets with 12 μm axial resolution (100 frames/s with 160 A-scans per frame for OCT imaging only). The high signal-to-background characteristic of DOCT provides exceptional sensitivity in detecting and quantifying the blood flow within the sheep brain parenchyma in real time. The OCT/DOCT dual imaging also demonstrated its capability to differentiate the vessel type (artery/vein) on rat's femoral vessels. We also demonstrated in ex vivo human brain that the location of the tip of the OCT probe can be inferred from micro-anatomical landmarks in OCT images. These findings demonstrate the suitability of OCT guidance during stereotactic procedures in the brain and its potential for reducing the risk of cerebral hemorrhage.

We present what we believe to be a new mathematical model of Ca(2+) leak from the sarcoplasmic reticulum (SR) in the heart. To our knowledge, it is the first to incorporate a realistic number of Ca(2+)-release units, each containing a cluster of stochastically gating Ca(2+) channels (RyRs), whose biophysical properties (e.g., Ca(2+) sensitivity and allosteric interactions) are informed by the latest molecular investigations. This realistic model allows for the detailed characterization of RyR Ca(2+)-release properties, and shows how this balances reuptake by the SR Ca(2+) pump. Simulations reveal that SR Ca(2+) leak consists of brief but frequent single RyR openings (~3000 cell(-1) s(-1)) that are likely to be experimentally undetectable, and are, therefore, "invisible". We also observe that these single RyR openings can recruit additional RyRs to open, due to elevated local (Ca(2+)), and occasionally lead to the generation of Ca(2+) sparks (~130 cell(-1) s(-1)). Furthermore, this physiological formulation of "invisible" leak allows for the removal of the ad hoc, non-RyR mediated Ca(2+) leak terms present in prior models. Finally, our model shows how Ca(2+) sparks can be robustly triggered and terminated under both normal and pathological conditions. Together, these discoveries profoundly influence how we interpret and understand diverse experimental and clinical results from both normal and diseased hearts., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

A three dimensional model of calcium dynamics in the rat ventricular myocyte was developed to study the mechanism of calcium homeostasis and pathological calcium dynamics during calcium overload. The model contains 20,000 calcium release units (CRUs) each containing 49 ryanodine receptors. The model simulates calcium sparks with a realistic spontaneous calcium spark rate. It suggests that in addition to the calcium spark-based leak, there is an invisible calcium leak caused by the stochastic opening of a small number of ryanodine receptors in each CRU without triggering a calcium spark. The model also explores the mechanism of calcium wave propagation between release sites under the conditions of calcium overload.