Approximation for Cooperative Interactions of a Spatially-Detailed Cardiac Sarcomere Model

We developed a novel ordinary differential equation (ODE) model, which produced results that correlated well with the Monte Carlo (MC) simulation when applied to a spatially-detailed model of the cardiac sarcomere. Configuration of the novel ODE model was based on the Ising model of myofilaments, with the “co-operative activation” effect introduced to incorporate nearest-neighbor interactions. First, a set of parameters was estimated using arbitrary Ca transient data to reproduce the combinational probability for the states of three consecutive regulatory units, using single unit probabilities for central and neighboring units in the MC simulation. The parameter set thus obtained enabled the calculation of the state transition of each unit using the ODE model with reference to the neighboring states. The present ODE model not only provided good agreement with the MC simulation results but was also capable of reproducing a wide range of experimental results under both steady-state and dynamic conditions including shortening twitch. The simulation results suggested that the nearest-neighbor interaction is a reasonable approximation of the cooperativity based on end-to-end interactions. Utilizing the modified ODE model resulted in a reduction in computational costs but maintained spatial integrity and co-operative effects, making it a powerful tool in cardiac modeling

Transient disruption of M1 during response planning impairs subsequent offline consolidation

Transcranial magnetic stimulation (TMS) was used to probe the involvement of the left primary motor cortex (M1) in the consolidation of a sequencing skill. In particular we asked: (1) if M1 is involved in consolidation of planning processes prior to response execution (2) whether movement preparation and movement execution can undergo consolidation independently and (3) whether sequence consolidation can occur in a stimulus specific manner. TMS was applied to left M1 while subjects prepared left hand sequential finger responses for three different movement sequences, presented in an interleaved fashion. Subjects also trained on three control sequences, where no TMS was applied. Disruption of subsequent consolidation was observed, but only for sequences where subjects had been exposed to TMS during training. Further, reduced consolidation was only observed for movement preparation, not movement execution. We conclude that left M1 is causally involved in the consolidation of effective response planning for left hand movements prior to response execution, and mediates consolidation in a sequence specific manner. These results provide important new insights into the role of M1 in sequential memory consolidation and sequence response planning

Simpson's paradox and calculation of number needed to treat from meta-analysis

BACKGROUND: Calculation of numbers needed to treat (NNT) is more complex from meta-analysis than from single trials. Treating the data as if it all came from one trial may lead to misleading results when the trial arms are imbalanced. DISCUSSION: An example is shown from a published Cochrane review in which the benefit of nursing intervention for smoking cessation is shown by formal meta-analysis of the individual trial results. However if these patients were added together as if they all came from one trial the direction of the effect appears to be reversed (due to Simpson's paradox). Whilst NNT from meta-analysis can be calculated from pooled Risk Differences, this is unlikely to be a stable method unless the event rates in the control groups are very similar. Since in practice event rates vary considerably, the use a relative measure, such as Odds Ratio or Relative Risk is advocated. These can be applied to different levels of baseline risk to generate a risk specific NNT for the treatment. SUMMARY: The method used to calculate NNT from meta-analysis should be clearly stated, and adding the patients from separate trials as if they all came from one trial should be avoided

Human Platelet-Rich Plasma- and Extracellular Matrix-Derived Peptides Promote Impaired Cutaneous Wound Healing In Vivo

Previous work in our laboratory has described several pro-angiogenic short peptides derived from endothelial extracellular matrices degraded by bacterial collagenase. Here we tested whether these peptides could stimulate wound healing in vivo. Our experiments demonstrated that a peptide created as combination of fragments of tenascin X and fibrillin 1 (comb1) applied into cranial dermal wounds created in mice treated with cyclophosphamide to impair wound healing, can improve the rate of wound closure. Furthermore, we identify and characterize a novel peptide (UN3) created and modified from two naturally-occurring peptides, which are present in human platelet-rich plasma. In vitro testing of UN3 demonstrates that it causes a 50% increase in endothelial proliferation, 250% increase in angiogenic response and a tripling of epithelial cell migration in response to injury. Results of in vivo experiments where comb1 and UN3 peptides were added together to cranial wounds in cyclophosphamide-treated mice leads to improvement of wound vascularization as shown by an increase of the number of blood vessels present in the wound beds. Application of the peptides markedly promotes cellular responses to injury and essentially restores wound healing dynamics to those of normal, acute wounds in the absence of cyclophosphamide impairment. Our current work is aimed at understanding the mechanisms underlying the stimulatory effects of these peptides as well as identification of the cellular receptors mediating these effects.National Institutes of Health (U.S.) (Grant EY15125)National Institutes of Health (U.S.) (Grant EY19533)Wound Care Partners, LL

Congenital hereditary endothelial dystrophy caused by SLC4A11 mutations progresses to Harboyan syndrome.

Purpose: Homozygous mutations in SLC4A11 cause 2 rare recessive conditions: congenital hereditary endothelial dystrophy (CHED), affecting the cornea alone, and Harboyan syndrome consisting of corneal dystrophy and sensorineural hearing loss. In addition, adult-onset Fuchs endothelial corneal dystrophy (FECD) is associated with dominant mutations in SLC4A11. In this report, we investigate whether patients with CHED go on to develop hearing loss and whether their parents, who are carriers of an SLC4A11 mutation, show signs of having FECD. Methods: Patients with CHED were screened for mutations in the SLC4A11 gene and underwent audiometric testing. The patients and their parents underwent a clinical examination and specular microscopy. Results: Molecular analyses confirmed SLC4A11 mutations in 4 affected individuals from 3 families. All the patients were found to have varying degrees of sensorineural hearing loss at a higher frequency range. Guttate lesions were seen in 2 of the 4 parents who were available for examination. Conclusions: Our observations suggest that CHED caused by homozygous SLC4A11 mutations progresses to Harboyan syndrome, but the severity of this may vary considerably. Patients with CHED should therefore be monitored for progressive hearing loss. We could not determine conclusively whether the parents of the patients with CHED were at increased risk of developing late-onset FECD

Ischemia reperfusion dysfunction changes model-estimated kinetics of myofilament interaction due to inotropic drugs in isolated hearts

BACKGROUND: The phase-space relationship between simultaneously measured myoplasmic [Ca(2+)] and isovolumetric left ventricular pressure (LVP) in guinea pig intact hearts is altered by ischemic and inotropic interventions. Our objective was to mathematically model this phase-space relationship between [Ca(2+)] and LVP with a focus on the changes in cross-bridge kinetics and myofilament Ca(2+ )sensitivity responsible for alterations in Ca(2+)-contraction coupling due to inotropic drugs in the presence and absence of ischemia reperfusion (IR) injury. METHODS: We used a four state computational model to predict LVP using experimentally measured, averaged myoplasmic [Ca(2+)] transients from unpaced, isolated guinea pig hearts as the model input. Values of model parameters were estimated by minimizing the error between experimentally measured LVP and model-predicted LVP. RESULTS: We found that IR injury resulted in reduced myofilament Ca(2+ )sensitivity, and decreased cross-bridge association and dissociation rates. Dopamine (8 μM) reduced myofilament Ca(2+ )sensitivity before, but enhanced it after ischemia while improving cross-bridge kinetics before and after IR injury. Dobutamine (4 μM) reduced myofilament Ca(2+ )sensitivity while improving cross-bridge kinetics before and after ischemia. Digoxin (1 μM) increased myofilament Ca(2+ )sensitivity and cross-bridge kinetics after but not before ischemia. Levosimendan (1 μM) enhanced myofilament Ca(2+ )affinity and cross-bridge kinetics only after ischemia. CONCLUSION: Estimated model parameters reveal mechanistic changes in Ca(2+)-contraction coupling due to IR injury, specifically the inefficient utilization of Ca(2+ )for contractile function with diastolic contracture (increase in resting diastolic LVP). The model parameters also reveal drug-induced improvements in Ca(2+)-contraction coupling before and after IR injury

Estimating Genetic Variability in Non-Model Taxa: A General Procedure for Discriminating Sequence Errors from Actual Variation

Genetic variation is the driving force of evolution and as such is of central interest for biologists. However, inadequate discrimination of errors from true genetic variation could lead to incorrect estimates of gene copy number, population genetic parameters, phylogenetic relationships and the deposition of gene and protein sequences in databases that are not actually present in any organism. Misincorporation errors in multi-template PCR cloning methods, still commonly used for obtaining novel gene sequences in non-model species, are difficult to detect, as no previous information may be available about the number of expected copies of genes belonging to multi-gene families. However, studies employing these techniques rarely describe in any great detail how errors arising in the amplification process were detected and accounted for. Here, we estimated the rate of base misincorporation of a widely-used PCR-cloning method, using a single copy mitochondrial gene from a single individual to minimise variation in the template DNA, as 1.62×10−3 errors per site, or 9.26×10−5 per site per duplication. The distribution of errors among sequences closely matched that predicted by a binomial distribution function. The empirically estimated error rate was applied to data, obtained using the same methods, from the Phospholipase A2 toxin family from the pitviper Ovophis monticola. The distribution of differences detected closely matched the expected distribution of errors and we conclude that, when undertaking gene discovery or assessment of genetic diversity using this error-prone method, it will be informative to empirically determine the rate of base misincorporation

Model-Derived Dispersal Pathways from Multiple Source Populations Explain Variability of Invertebrate Larval Supply

Background: Predicting the spatial and temporal patterns of marine larval dispersal and supply is a challenging task due to the small size of the larvae and the variability of oceanographic processes. Addressing this problem requires the use of novel approaches capable of capturing the inherent variability in the mechanisms involved. Methodology/Principal Findings: In this study we test whether dispersal and connectivity patterns generated from a biophysical model of larval dispersal of the crab Carcinus maenas, along the west coast of the Iberian Peninsula, can predict the highly variable daily pattern of wind-driven larval supply to an estuary observed during the peak reproductive season (March–June) in 2006 and 2007. Cross-correlations between observed and predicted supply were significant (p,0.05) and strong, ranging from 0.34 to 0.81 at time lags of 26 to+5 d. Importantly, the model correctly predicted observed cross-shelf distributions (Pearson r = 0.82, p,0.001, and r = 0.79, p,0.01, in 2006 and 2007) and indicated that all supply events were comprised of larvae that had been retained within the inner shelf; larvae transported to the outer shelf and beyond never recruited. Estimated average dispersal distances ranged from 57 to 198 km and were only marginally affected by mortality. Conclusions/Significance: The high degree of predicted demographic connectivity over relatively large geographic scales is consistent with the lack of genetic structuring in C. maenas along the Iberian Peninsula. These findings indicate that the dynamic nature of larval dispersal can be captured by mechanistic biophysical models, which can be used to provid