The fully implicit stochastic-alpha method for stiff stochastic differential equations

A fully implicit integration method for stochastic differential equations with significant multiplicative noise and stiffness in both the drift and diffusion coefficients has been constructed, analyzed and illustrated with numerical examples in this work. The method has strong order 1.0 consistency and has user-selectable parameters that allow the user to expand the stability region of the method to cover almost the entire drift-diffusion stability plane. The large stability region enables the method to take computationally efficient time steps. A system of chemical Langevin equations simulated with the method illustrates its computational efficiency.

Structure and function of GDNF receptor alpha splice variants

The growth factors of the glial cell line-derived neurotrophic factor (GDNF) family consisting of GDNF, neurturin (NRTN), artemin (ARTN) and persephin (PSPN), are involved in the development, differentiation and maintenance of many types of neurons. They also have important functions outside the nervous system in the development of kidney, testis and thyroid gland. Each of these GFLs preferentially binds to one of the glycosylphosphatidylinositol (GPI)-anchored GDNF family receptors α (GFRα). GDNF binds to GFRα1, NRTN to GFRα2, ARTN to GFRα3 and PSPN to GFRα4. The GFLs in the complex with their cognate GFRα receptors all bind to and signal through the receptor tyrosine kinase RET. Alternative splicing of the mouse GFRα4 gene yields three splice isoforms. These had been described as putative GPI-anchored, transmembrane and soluble forms. My goal was to characterise the function of the different forms of mouse GFRα4. I firstly found that the putative GPI-anchored GFRα4 (GFRα4-GPI) is glycosylated, membrane-bound, GPI-anchored and interacts with PSPN and RET. We also showed that mouse GFRα4-GPI mediates PSPN-induced phosphorylation of RET, promotes PSPN-dependent neuronal differentiation of the rat pheochromocytoma cell line PC6-3 and PSPN-dependent survival of cerebellar granule neurons (CGN). However, although this receptor can mediate PSPN-signalling and activate RET, GFRα4-GPI does not recruit RET into lipid rafts. The recruitment of RET into lipid rafts has previously been thought to be a crucial event for GDNF- and GFL-mediated signalling via RET. I secondly demonstrated that the putative transmembrane GFRα4 (GFRα4-TM) is indeed a real transmembrane GFRα4 protein. Although it has a weak binding capacity for PSPN, it can not mediate PSPN-dependent phosphorylation of RET, neuronal differentiation or survival. These data show that GFRα4-TM is inactive as a receptor for PSPN. Surprisingly, GFRα4-TM can negatively regulate PSPN-mediated signalling via GFRα4-GPI. GFRα4-TM interacts with GFRα4-GPI and blocks PSPN-induced phosphorylation of RET, neuronal differentiation as well as survival. Taken together, our data show that GFRα4-TM may act as a dominant negative inhibitor of PSPN-mediated signaling. The most exciting part of my work was the finding that the putative soluble GFRα4 (GFRα4-sol) can form homodimers and function as an agonist of the RET receptor. In the absence of PSPN, GFRα4-sol can promote the phosphorylation of RET, trigger the activation of the PI-3K/AKT pathway, induce neuronal differentiation and support the survival of CGN. Our findings are in line with a recent publication showing the GFRα4-sol might contribute to the inherited cancer syndrome multiple endocrine neoplasia type 2. Our data provide an explanation to how GFRα4-sol may cause or modify the disease. Mammalian GFRα4 receptors all lack the first Cys-rich domain which is present in other GFRα receptors. In the final part of my work I have studied the function of this particular domain. I created a truncated GFRα1 construct lacking the first Cys-rich domain. Using binding assays in both cellular and cell-free systems, phosphorylation assays with RET, as well as neurite outgrowth assays, we found that the first Cys-rich domain contributes to an optimal function of GFRα1, by stabilizing the interaction between GDNF and GFRα1.

Kinetic characterization of rat brain type IIA sodium channel alpha-subunit stably expressed in a somatic cell line

1. The rat brain type IIA Na+ channel alpha-subunit was stably expressed in Chinese hamster ovary (CHO) cells. Current through the expressed Na+ channels was studied using the whole-cell configuration of the patch clamp technique. The transient Na+ current was sensitive to TTX and showed a bell-shaped peak current vs. membrane potential relation. 2. Na+ current inactivation was better described by the sum of two exponentials in the potential range -30 to +40 mV, with. a dominating fast component and a small slower component. 3. The steady-state inactivation, h(infinity), was related to potential by a Boltzmann distribution, underlying thr ee states of the inactivation gate. 4. Recovery of the channels from inactivation at different potentials in the range -70 to -120 mV were characterized by al? initial delay which decreased with hyperpolarization. The time course was well fitted by the sum of two exponentials. In this case the slower exponential was the major component, and both time constants decreased with hyperpolarization. 5. For a working description of the Na+ channel inactivation in this preparation, with a minimal deviation from the Hodgkin-Huxley model, a three-state scheme of the form O reversible arrow I-1 reversible arrow I-2 was proposed, replacing the original two-state scheme of the Hodgkin-Huxley model, and the rate constants are reported. 6. The instantaneous current-voltage relationship showed marked deviation from linearity and was satisfactorily fitted by the constant-field equation. 7. The time course of activation was described by an m(x) model. However, the best-fitted value of x varied with the membrane potential and had a mean value of 2. 8. Effective gating charge was determined to be 4.7e from the slope of the activation plot, plotted on a logarithmic scale. 9. The rate constants of activation, alpha(m) and beta(m), were determined. Their functional dependence on the membrane potential was investigated.

Identification of alpha- and beta-Hydroxy Acid Containing Cyclodepsipeptides in Natural Peptide Mixtures Using Negative Ion Mass Spectrometry

Natural peptide libraries often contain cyclodepsipeptides containing alpha or beta hydroxy residues. Extracts of fungal hyphae of Isaria yield a microheterogenous cyclodepsipeptide mixture in which two classes of molecules can be identified by mass spectral fragmentation of negative ions. In the case of isaridins, which contain an alpha-hydroxy residue and a beta-amino acid residue, a characteristic product ion corresponding to a neutral loss of 72 Da is obtained. hi addition, neutral loss of water followed by a 72 Da loss is also observed. Two distinct modes of fragmentation rationalize the observed product ion distribution. The neutral loss of 72 Da has also been obtained for a roseotoxin component, which is also an alpha-hydroxy residue containing cyclodepsipeptide. In the case of isariins, which contain a beta-hydroxy acid residue, ring opening and subsequent loss of the terminal residue as an unsaturated ketene fragment, rationalizes the observed product ion formation. Fragmentation of negative ions provide characteristic neutral losses, which are diagnostic of the presence of alpha-hydroxy or beta-hydroxy residues.

The alpha method for solving differential algebraic inequality (DAI) systems

This paper describes an algorithm for ``direct numerical integration'' of the initial value Differential-Algebraic Inequalities (DAI) in a time stepping fashion using a sequential quadratic programming (SQP) method solver for detecting and satisfying active path constraints at each time step. The activation of a path constraint generally increases the condition number of the active discretized differential algebraic equation's (DAE) Jacobian and this difficulty is addressed by a regularization property of the alpha method. The algorithm is locally stable when index 1 and index 2 active path constraints and bounds are active. Subject to available regularization it is seen to be stable for active index 3 active path constraints in the numerical examples. For the high index active path constraints, the algorithm uses a user-selectable parameter to perturb the smaller singular values of the Jacobian with a view to reducing the condition number so that the simulation can proceed. The algorithm can be used as a relatively cheaper estimation tool for trajectory and control planning and in the context of model predictive control solutions. It can also be used to generate initial guess values of optimization variables used as input to inequality path constrained dynamic optimization problems. The method is illustrated with examples from space vehicle trajectory and robot path planning.

Recognition of Interaction Interface Residues in Low-Resolution Structures of Protein Assemblies Solely from the Positions of C alpha Atoms

Background: The number of available structures of large multi-protein assemblies is quite small. Such structures provide phenomenal insights on the organization, mechanism of formation and functional properties of the assembly. Hence detailed analysis of such structures is highly rewarding. However, the common problem in such analyses is the low resolution of these structures. In the recent times a number of attempts that combine low resolution cryo-EM data with higher resolution structures determined using X-ray analysis or NMR or generated using comparative modeling have been reported. Even in such attempts the best result one arrives at is the very course idea about the assembly structure in terms of trace of the C alpha atoms which are modeled with modest accuracy. Methodology/Principal Findings: In this paper first we present an objective approach to identify potentially solvent exposed and buried residues solely from the position of C alpha atoms and amino acid sequence using residue type-dependent thresholds for accessible surface areas of C alpha. We extend the method further to recognize potential protein-protein interface residues. Conclusion/Significance: Our approach to identify buried and exposed residues solely from the positions of C alpha atoms resulted in an accuracy of 84%, sensitivity of 83-89% and specificity of 67-94% while recognition of interfacial residues corresponded to an accuracy of 94%, sensitivity of 70-96% and specificity of 58-94%. Interestingly, detailed analysis of cases of mismatch between recognition of interface residues from C alpha positions and all-atom models suggested that, recognition of interfacial residues using C alpha atoms only correspond better with intuitive notion of what is an interfacial residue. Our method should be useful in the objective analysis of structures of protein assemblies when positions of only C alpha positions are available as, for example, in the cases of integration of cryo-EM data and high resolution structures of the components of the assembly.

Role of thrombin in coronary artery bypass grafting

Thrombin is a multifunctional protease, which has a central role in the development and progression of coronary atherosclerotic lesions and it is a possible mediator of myocardial ischemia-reperfusion injury. Its generation and procoagulant activity are greatly upregulated during cardiopulmonary bypass (CPB). On the other hand, activated protein C, a physiologic anticoagulant that is activated by thrombomodulin-bound thrombin, has been beneficial in various models of ischemia-reperfusion. Therefore, our aim in this study was to test whether thrombin generation or protein C activation during coronary artery bypass grafting (CABG) associate with postoperative myocardial damage or hemodynamic changes. To further investigate the regulation of thrombin during CABG, we tested whether preoperative thrombophilic factors associate with increased CPB-related generation of thrombin or its procoagulant activity. We also measured the anticoagulant effects of heparin during CPB with a novel coagulation test, prothrombinase-induced clotting time (PiCT), and compared the performance of this test with the present standard of laboratory-based anticoagulation monitoring. One hundred patients undergoing elective on-pump CABG were studied prospectively. A progressive increase in markers of thrombin generation (F1+2), fibrinolysis (D-dimer), and fibrin formation (soluble fibrin monomer complexes) was observed during CPB, which was further distinctly propagated by reperfusion after myocardial ischemia, and continued to peak after the neutralization of heparin with protamine. Thrombin generation during reperfusion after CABG associated with postoperative myocardial damage and increased pulmonary vascular resistance. Activated protein C levels increased only slightly during CPB before the release of the aortic clamp, but reperfusion and more significantly heparin neutralization caused a massive increase in activated protein C levels. Protein C activation was clearly delayed in relation to both thrombin generation and fibrin formation. Even though activated protein C associated dynamically with postoperative hemodynamic performance, it did not associate with postoperative myocardial damage. Preoperative thrombophilic variables did not associate with perioperative thrombin generation or its procoagulant activity. Therefore, our results do not favor routine thrombophilia screening before CABG. There was poor agreement between PiCT and other measurements of heparin effects in the setting of CPB. However, lower heparin levels during CPB associated with inferior thrombin control and high heparin levels during CPB associated with fewer perioperative transfusions of blood products. Overall, our results suggest that hypercoagulation after CABG, especially during reperfusion, might be clinically important.

H N.M.R. studies of protected alpha-aminoisobutyric acid containing peptides. Chemical shift nonequivalence of benzyloxycarbonyl methylene protons

The benzylic methylene protons in a large number of benzyloxycarbonyl alpha-aminoisobutyric acid (Z-Aib) containing peptides, show chemical shift nonequivalence. The magnitude of the geminal nonequivalence is correlated with the involvement of the urethane carbonyl group, in an intramolecular hydrogen bond. Studies of the model compounds Z-Aib-Aib-Ala-NHMe, and Z-Aib-Aib-Aib-Pro-OMe clearly establish the presence of intramolecular hydrogen bonds, involving the urethane CO group. In both compounds marked anisochrony of the benzylic methylene protons is demonstrated. In Z-Aib-Aib-Pro-OMe, where a 4 leads to 1 hydrogen bonded beta-turn is not possible, the benzylic-CH2-protons appear as a singlet in CDCl3 and have a very small chemical shift difference in (CD3)2SO. The observation of such nonequivalence is of value in establishing whether the amino terminal Aib-Pro beta-turn is retained in large peptide-fragments of alamethicin.

CD4+ alpha beta T cell and gamma delta T cell responses to Mycobacterium tuberculosis. Similarities and differences in Ag recognition, cytotoxic effector function, and cytokine production.

CD4+ and gamma delta T cells are activated readily by Mycobacterium tuberculosis. To examine their role in the human immune response to M. tuberculosis, CD4+ and gamma delta T cells from healthy tuberculin-positive donor were studied for patterns of Ag recognition, cytotoxicity, and cytokine production in response to M. tuberculosis-infected mononuclear phagocytes. Both T cell subsets responded to intact M. tuberculosis and its cytosolic Ags. However, CD4+ and gamma delta T cells differed in the range of cytosolic Ags recognized: reactivity to a wide m.w. range of Ags for CD4+ T cells, and a restricted pattern for gamma delta T cells, with dominance of Ags of 10 to 15 kDa. Both T cell subsets were equally cytotoxic for M. tuberculosis-infected monocytes. Furthermore, both CD4+ and gamma delta T cells produced large amounts of IFN-gamma: mean pg/ml of IFN-gamma in supernatants was 2458 +/- 213 for CD4+ and 2349 +/- 245 for gamma delta T cells. By filter-spot ELISA (ELISPOT), the frequency of IFN-gamma-secreting gamma delta T cells was one-half of that of CD4+ T cells in response to M. tuberculosis, suggesting that gamma delta T cells on a per cell basis were more efficient producers of IFN-gamma than CD4+ T cells. In contrast, CD4+ T cells produced more IL-2 than gamma delta T cells, which correlated with diminished T cell proliferation of gamma delta T cells compared with CD4+ T cells. These results indicate that CD4+ and gamma delta T cell subsets have similar effector functions (cytotoxicity, IFN-gamma production) in response to M. tuberculosis-infected macrophages, despite differences in the Ags recognized, IL-2 production, and efficiency of IFN-gamma production.

Expanding the Peptide beta-Turn in alpha gamma Hybrid Sequences: 12 Atom Hydrogen Bonded Helical and Hairpin Turns

Hybrid peptide segments containing contiguous alpha and gamma amino acid residues can form C-12 hydrogen bonded turns which may be considered as backbone expanded analogues of C-10 beta-turns) found in alpha alpha segments. Exploration of the regular hydrogen bonded conformations accessible for hybrid alpha gamma sequences is facilitated by the use of a stereochemically constrained gamma amino acid residue gabapentin (1-aminomethylcyclohexaneacetic acid, Gpn), in which the two torsion angles about C-gamma-C-beta (theta(1)) and C-beta-C-alpha (theta(2)) are predominantly restricted to gauche conformations. The crystal structures of the octapeptides Boc-Gpn-Aib-Gpn-Aib-Gpn-Aib-Gpn-Aib-OMe (1) and Boc-Leu-Phe-Val-Aib-Gpn-Leu-Phe-Val-OMe (2) reveal two distinct conformations for the Aib-Gpn segment. Peptide 1 forms a continuous helix over the Aib(2)-Aib(6) segment, while the peptide 2 forms beta-hairpin structure stabilized by four cross-strand hydrogen bonds with the Aib-Gpn segment forming a nonhelical C-12 turn. The robustness of the helix in peptide 1 in solution is demonstrated by NMR methods. Peptide 2 is conformationally fragile in solution with evidence of beta-hairpin conformations being obtained in methanol. Theoretical calculations permit delineation of the various C-12 hydrogen bonded structures which are energetically feasible in alpha gamma and gamma alpha sequences.

Hydrothermal synthesis and structural characterization of novel alpha-Keggin unit-supported Cu(II)- and Mn(II)-bipyridine complexes from a tri-lacunary precursor

Blue [{Cu(2,2'-bipy)(2)}(2){alpha-SiW12O40}] (bipy = bipyridyl) (1) and pale yellow [Mn(2,2'-bipy)(3)](2)[alpha-SiW12O40] (2) have been synthesized hydrothermally and characterized by IR spectroscopy and single crystal X-ray structure analysis. In 1, the [alpha-SiW12O40](4-) ion acts as a bridge between the two [{Cu(2,2'-bipy)(2)](2+) moieties via coordination through the terminal oxygen atoms, while in 2, the [Mn(2,2'-bipy)(3)](2+) ion balances the charge on the polyoxo anion without forming any covalent bond. To the best of our knowledge, this is the first example of transition metal-mediated transformation of [alpha-SiW9O34](10-) to [alpha-SiW12O40](4-).

Inhibition of Thrombin in Cardiac Surgery : experiments in a porcine model

Cardiac surgery involving cardiopulmonary bypass (CPB) induces activation of inflammation and coagulation systems and is associated with ischemia-reperfusion injury (I/R injury)in various organs including the myocardium, lungs, and intestine. I/R injury is manifested as organ dysfunction. Thrombin, the key enzyme of coagulation , plays a cenral role also in inflammation and contributes to regulation of apoptosis as well. The general aim of this thesis was to evaluate the potential of thrombin inhibition in reducing the adverse effects of I/R injury in myocardium, lungs, and intestine associated with the use of CPB and cardiac surgery. Forty five pigs were used for the studies. Two randomized blinded studies were performed. Animals underwent 75 min of normothermic CPB, 60 min of aortic clamping, and 120 min of reperfusion period. Twenty animals received iv. recombinant hirudin, a selective and effective inbitor of thrombin, or placebo. In a similar setting, twenty animals received an iv-bolus (250 IU/kg) of antithrombin (AT) or placebo. An additional group of 5 animals received 500 IU/kg in an open label setting to test dose response. Generation of thrombin (TAT), coagulation status (ACT), and hemodynamics were measured. Intramucosal pH and pCO2 were measured from the luminal surface of ileum using tonometry simultaneusly with arterial gas analysis. In addition, myocardial, lung, and intestinal biopsies were taken to quantitate leukocyte infiltration (MPO), for histological evaluation, and detection of apoptosis (TUNEL, caspase 3). In conclusion, our data suggest that r-hirudin may be an effective inhibitor of reperfusion induced thrombin generation in addition to being a direct inhibitor of preformed thrombin. Overall, the results suggest that inhibition of thrombin, beyond what is needed for efficient anticoagulation by heparin, has beneficial effects on myocardial I/R injury and hemodynamics during cardiac surgery and CPB. We showed that infusion of the thrombin inhibitor r-hirudin during reperfusion was associated with attenuated post ischemia left ventricular dysfunction and decreased systemic vascular resistance. Consequently microvascular flow was improved during ischemia-reperfusion injury. Improved recovery of myocardium during the post-ischemic reperfusion period was associated with significantly less cardiomyocyte apoptosis and with a trend in anti-inflammatory effects. Thus, inhibition of reperfusion induced thrombin may offer beneficial effects by mechanisms other than direct anticoagulant effects. AT, in doses with a significant anticoagulant effect, did not alleviate myocardial I/R injury in terms of myocardial recovery, histological inflammatory changes or post-ischemic troponin T release. Instead, AT attenuated reperfusion induced increase in pulmonary pressure after CPB. Taken the clinical significance of postoperative pulmonary hemodynamics in patients undergoing cardiopulmonary bypass, the potential positive regulatory role of AT and clinical implications needs to be studied further. Inflammatory response in the gut wall proved to be poorly associated with perturbed mucosal perfusion and the animals with the least neutrophil tissue sequestration and I/R related histological alterations tended to have the most progressive mucosal hypoperfusion. Thus, mechanisms of low-flow reperfusion injury during CPB can differ from the mechanisms seen in total ischemia reperfusion injury.

Characterization of mouse neuronal Ca2+/calmodulin kinase II inhibitor alpha

We have overexpressed an 8.5-kDa mouse Ca2+/calmodulin kinase II inhibitor a protein (mCaMKIIN alpha) in Escherichia coli and demonstrate that the recombinant protein is a potent inhibitor of Ca2+/calmodulin kinase 11 (CaMKII) in vitro. However, antibodies raised against recombinant mCaMKIIN alpha. react with an similar to 37-kDa protein present in mouse brain. The pattern of expression of the similar to 37-kDa protein is similar to that of mCaMKIIN alpha mRNA as both are expressed in normal but not Japanese encephalitis virus (JEV)-infected mouse brain. Subcellular localization studies indicate that the similar to 37-kDa protein is present in the post-synaptic density (PSD) where mCaMKII alpha is known to perform key regulatory functions. We conclude that the similar to 37-kDa protein identified in this study is mCaMKIIN alpha. and its localization in the PSD indicates a novel role for this protein in the regulation of neuronal CaMKII alpha. (c) 2007 Elsevier B.V. All rights reserved.

Thrombophilia and direct thrombin inhibitor lepirudin clinical and monitoring aspects

Thrombophilia (TF) predisposes both to venous and arterial thrombosis at a young age. TF may also impact the thrombosis or stenosis of hemodialysis (HD) vascular access in patients with end-stage renal disease (ESRD). When involved in severe thrombosis TF may associate with inappropriate response to anticoagulation. Lepirudin, a potent direct thrombin inhibitor (DTI), indicated for heparin-induced thrombocytopenia-related thrombosis, could offer a treatment alternative in TF. Monitoring of narrow-ranged lepirudin demands new insights also in laboratory. The above issues constitute the targets in this thesis. We evaluated the prevalence of TF in patients with ESRD and its impact upon thrombosis- or stenosis-free survival of the vascular access. Altogether 237 ESRD patients were prospectively screened for TF and thrombogenic risk factors prior to HD access surgery in 2002-2004 (mean follow-up of 3.6 years). TF was evident in 43 (18%) of the ESRD patients, more often in males (23 vs. 9%, p=0.009). Known gene mutations of FV Leiden and FII G20210A occurred in 4%. Vascular access sufficiently matured in 226 (95%). The 1-year thrombosis- and stenosis-free access survival was 72%. Female gender (hazards ratio, HR, 2.5; 95% CI 1.6-3.9) and TF (HR 1.9, 95% CI 1.1-3.3) were independent risk factors for the shortened thrombosis- and stenosis-free survival. Additionally, TF or thrombogenic background was found in relatively young patients having severe thrombosis either in hepatic veins (Budd-Chiari syndrome, BCS, one patient) or inoperable critical limb ischemia (CLI, six patients). Lepirudin was evaluated in an off-label setting in the severe thrombosis after inefficacious traditional anticoagulation without other treatment options except severe invasive procedures, such as lower extremity amputation. Lepirudin treatments were repeatedly monitored clinically and with laboratory assessments (e.g. activated partial thromboplastin time, APTT). Our preliminary studies with lepirudin in thrombotic calamities appeared safe, and no bleeds occurred. An effective DTI lepirudin calmed thrombosis as all patients gradually recovered. Only one limb amputation was performed 3 years later during the follow-up (mean 4 years). Furthermore, we aimed to overcome the limitations of APTT and confounding effects of warfarin (INR of 1.5-3.9) and lupus anticoagulant (LA). Lepirudin responses were assessed in vitro by five specific laboratory methods. Ecarin chromogenic assay (ECA) or anti-Factor IIa (anti-FIIa) correlated precisely (r=0.99) with each other and with spiked lepirudin in all plasma pools: normal, warfarin, and LA-containing plasma. In contrast, in the presence of warfarin and LA both APTT and prothrombinase-induced clotting time (PiCT®) were limited by non-linear and imprecise dose responses. As a global coagulation test APTT is useful in parallel to the precise chromogenic methods ECA or Anti-FIIa in challenging clinical situations. Lepirudin treatment requires multidisciplinary approach to ensure appropriate patient selection, interpretation of laboratory monitoring, and treatment safety. TF seemed to be associated with complicated thrombotic events, in venous (BCS), arterial (CLI), and vascular access systems. TF screening should be aimed to patients with repeated access complications or prior unprovoked thromboembolic events. Lepirudin inhibits free and clot-bound thrombin which heparin fails to inhibit. Lepirudin seems to offer a potent and safe option for treatment of severe thrombosis. Multi-centered randomized trials are necessary to assess the possible management of complicated thrombotic events with DTIs like lepirudin and seek prevention options against access complications.