Interactions Outside the Proteinase-binding Loop Contribute Significantly to the Inhibition of Activated Coagulation Factor XII by Its Canonical Inhibitor from Corn

Activated factor XII (FXIIa) is selectively inhibited by corn Hageman factor inhibitor (CHFI) among other plasma proteases. CHFI is considered a canonical serine protease inhibitor that interacts with FXIIa through its protease-binding loop. Here we examined whether the protease-binding loop alone is sufficient for the selective inhibition of serine proteases or whether other regions of a canonical inhibitor are involved. Six CHFI mutants lacking different N- and C-terminal portions were generated. CHFI-234, which lacks the first and fifth disulfide bonds and 11 and 19 amino acid residues at the N and C termini, respectively, exhibited no significant changes in FXIIa inhibition (K_i = 3.2 ± 0.4 nm). CHFI-123, which lacks 34 amino acid residues at the C terminus and the fourth and fifth disulfide bridges, inhibited FXIIa with a K_i of 116 ± 16 nm. To exclude interactions outside the FXIIa active site, a synthetic cyclic peptide was tested. The peptide contained residues 20–45 (Protein Data Bank code 1BEA), and a C29D substitution was included to avoid unwanted disulfide bond formation between unpaired cysteines. Surprisingly, the isolated protease-binding loop failed to inhibit FXIIa but retained partial inhibition of trypsin (K_i = 11.7 ± 1.2 μm) and activated factor XI (K_i = 94 ± 11 μm). Full-length CHFI inhibited trypsin with a K_i of 1.3 ± 0.2 nm and activated factor XI with a K_i of 5.4 ± 0.2 μm. Our results suggest that the protease-binding loop is not sufficient for the interaction between FXIIa and CHFI; other regions of the inhibitor also contribute to specific inhibition.     

Effectors of rat-heart hexokinases and the control of rates of glucose phosphorylation in the perfused rat heart

1. The kinetic properties of the soluble and particulate hexokinases from rat heart have been investigated. 2. For both forms of the enzyme, the K_m for glucose was 45μm and the K_m for ATP 0·5mm. Glucose 6-phosphate was a non-competitive inhibitor with respect to glucose (K_i 0·16mm for the soluble and 0·33mm for the particulate enzyme) and a mixed inhibitor with respect to ATP (K_i 80μm for the soluble and 40μm for the particulate enzyme). ADP and AMP were competitive inhibitors with respect to ATP (K_i for ADP was 0·68mm for the soluble and 0·60mm for the particulate enzyme; K_i for AMP was 0·37mm for the soluble and 0·16mm for the particulate enzyme). P_i reversed glucose 6-phosphate inhibition with both forms at 10mm but not at 2mm, with glucose 6-phosphate concentrations of 0·3mm or less for the soluble and 1mm or less for the particulate enzyme. 3. The total activity of hexokinase in normal hearts and in hearts from alloxan-diabetic rats was 21·5μmoles of glucose phosphorylated/min./g. dry wt. of ventricle at 25°. The temperature coefficient Q₁₀ between 22° and 38·5° was 1·93; the ratio of the soluble to the particulate enzyme was 3:7. 4. The kinetic data have been used to predict rates of glucose phosphorylation in the perfused heart at saturating concentrations of glucose from measured concentrations of ATP, glucose 6-phosphate, ADP and AMP. These have been compared with the rates of glucose phosphorylation measured with precision in a small-volume recirculation perfusion apparatus, which is described. The correlation between predicted and measured rates was highly significant and their ratio was 1·07. 5. These findings are consistent with the control of glucose phosphorylation in the perfused heart by glucose 6-phosphate concentration, subject to certain assumptions that are discussed in detail.     

Certhrax Toxin,an Anthrax-related ADP-ribosyltransferase from Bacillus cereus

We identified Certhrax, the first anthrax-like mART toxin from the pathogenic G9241 strain of Bacillus cereus. Certhrax shares 31% sequence identity with anthrax lethal factor from Bacillus anthracis; however, we have shown that the toxicity of Certhrax resides in the mART domain, whereas anthrax uses a metalloprotease mechanism. Like anthrax lethal factor, Certhrax was found to require protective antigen for host cell entry. This two-domain enzyme was shown to be 60-fold more toxic to mammalian cells than anthrax lethal factor. Certhrax localizes to distinct regions within mouse RAW264.7 cells by 10 min postinfection and is extranuclear in its cellular location. Substitution of catalytic residues shows that the mART function is responsible for the toxicity, and it binds NAD⁺ with high affinity (K_D = 52.3 ± 12.2 μm). We report the 2.2 Å Certhrax structure, highlighting its structural similarities and differences with anthrax lethal factor. We also determined the crystal structures of two good inhibitors (P6 (K_D = 1.7 ± 0.2 μm, K_i = 1.8 ± 0.4 μm) and PJ34 (K_D = 5.8 ± 2.6 μm, K_i = 9.6 ± 0.3 μm)) in complex with Certhrax. As with other toxins in this family, the phosphate-nicotinamide loop moves toward the NAD⁺ binding site with bound inhibitor. These results indicate that Certhrax may be important in the pathogenesis of B. cereus.     

Identification of difructose dianhydride I synthase/hydrolase from an oral bacterium establishes a novel glycoside hydrolase family

Fructooligosaccharides and their anhydrides are widely used as health-promoting foods and prebiotics. Various enzymes acting on β-D-fructofuranosyl linkages of natural fructan polymers have been used to produce functional compounds. However, enzymes that hydrolyze and form α-D-fructofuranosyl linkages have been less studied. Here, we identified the BBDE_2040 gene product from Bifidobacterium dentium (α-D-fructofuranosidase and difructose dianhydride I synthase/hydrolase from Bifidobacterium dentium [αFFase1]) as an enzyme with α-D-fructofuranosidase and α-D-arabinofuranosidase activities and an anomer-retaining manner. αFFase1 is not homologous with any known enzymes, suggesting that it is a member of a novel glycoside hydrolase family. When caramelized fructose sugar was incubated with αFFase1, conversions of β-D-Frup-(2→1)-α-D-Fruf to α-D-Fruf-1,2′:2,1′-β-D-Frup (diheterolevulosan II) and β-D-Fruf-(2→1)-α-D-Fruf (inulobiose) to α-D-Fruf-1,2′:2,1′-β-D-Fruf (difructose dianhydride I [DFA I]) were observed. The reaction equilibrium between inulobiose and DFA I was biased toward the latter (1:9) to promote the intramolecular dehydrating condensation reaction. Thus, we named this enzyme DFA I synthase/hydrolase. The crystal structures of αFFase1 in complex with β-D-Fruf and β-D-Araf were determined at the resolutions of up to 1.76 Å. Modeling of a DFA I molecule in the active site and mutational analysis also identified critical residues for catalysis and substrate binding. The hexameric structure of αFFase1 revealed the connection of the catalytic pocket to a large internal cavity via a channel. Molecular dynamics analysis implied stable binding of DFA I and inulobiose to the active site with surrounding water molecules. Taken together, these results establish DFA I synthase/hydrolase as a member of a new glycoside hydrolase family (GH172).     

Sugar transport across the peritubular face of renal cells of the flounder

The transport of some sugars at the antiluminal face of renal cells was studied using teased tubules of flounder (Pseudopleuronectes americanus). The analytical procedure allowed the determination of both free and total (free plus phosphorylated) tissue sugars. The inulin space of the preparation was 0.333 ± 0.017 kg/kg wet wt (7 animals, 33 analyses). The nonmetabolizable α-methyl-D-glucoside entered the cells by a carrier-mediated (phloridzin-sensitive), ouabain-insensitive process. The steady-state tissue/medium ratio was systematically below that for diffusion equilibrium. D-Glucose was a poor inhibitor of α-methyl-glucoside transport, D-galactose was ineffective. The phloridzin-sensitive transport processes of 2-deoxy-D-glucose,D-galactose,and 2-deoxy-D-galactose were associated with considerable phosphorylation. Kinetic evidence suggested that these sugars were transported in free form and subsequently were phosphorylated. 2-Deoxy-D-glucose accumulated in the cells against a slight concentration gradient. This transport was greatly inhibited by D-glucose, whereas α-methyl-glucoside and also D-galactose and its 2-deoxy-derivative were ineffective. D-Galactose and 2-deoxy-D-galactose mutually competed for transport; D-glucose, 2-deoxy-D-glucose, and α-methyl-D-glucoside were ineffective. Studies using various sugars as inhibitors suggest the presence of three carrier-mediated pathways of sugar transport at the antiluminal cell face of the flounder renal tubule: the pathway of α-methyl-D-glucoside (not shared by D-glucose); the pathway commonly shared by 2-deoxy-D-glucose and D-glucose; the pathway shared by D-galactose and 2-deoxy-D-galactose.     

Inhibition of purine phosphoribosyltransferases of Ehrlich ascites-tumour cells by 6-mercaptopurine

1. The formation of adenosine 5′-phosphate, guanosine 5′-phosphate and inosine 5′-phosphate from [8-¹⁴C]adenine, [8-¹⁴C]guanine and [8-¹⁴C]hypoxanthine respectively in the presence of 5-phosphoribosyl pyrophosphate and an extract from Ehrlich ascites-tumour cells was assayed by a method involving liquid-scintillation counting of the radioactive nucleotides on diethylaminoethylcellulose paper. The results obtained with guanine were confirmed by a spectrophotometric assay which was also used to assay the conversion of 6-mercaptopurine and 5-phosphoribosyl pyrophosphate into 6-thioinosine 5′-phosphate in the presence of 6-mercaptopurine phosphoribosyltransferase from these cells. 2. At pH 7·8 and 25° the Michaelis constants for adenine, guanine and hypoxanthine were 0·9 μm, 2·9 μm and 11·0 μm in the assay with radioactive purines; the Michaelis constant for guanine in the spectrophotometric assay was 2·6 μm. At pH 7·9 the Michaelis constant for 6-mercaptopurine was 10·9 μm. 3. 25 μm-6-Mercaptopurine did not inhibit adenine phosphoribosyltransferase. 6-Mercaptopurine is a competitive inhibitor of guanine phosphoribosyltransferase (K_i 4·7 μm) and hypoxanthine phosphoribosyltransferase (K_i 8·3 μm). Hypoxanthine is a competitive inhibitor of guanine phosphoribosyltransferase (K_i 3·4 μm). 4. Differences in kinetic parameters and in the distribution of phosphoribosyltransferase activities after electrophoresis in starch gel indicate that different enzymes are involved in the conversion of adenine, guanine and hypoxanthine into their nucleotides. 5. From the low values of K_i for 6-mercaptopurine, and from published evidence that ascites-tumour cells require supplies of purines from the host tissues, it is likely that inhibition of hypoxanthine and guanine phosphoribosyltransferases by free 6-mercaptopurine is involved in the biological activity of this drug.     

Novel human D-amino acid oxidase inhibitors stabilize an active-site lid-open conformation

The NMDAR (N-methyl-D-aspartate receptor) is a central regulator of synaptic plasticity and learning and memory. hDAAO (human D-amino acid oxidase) indirectly reduces NMDAR activity by degrading the NMDAR co-agonist D-serine. Since NMDAR hypofunction is thought to be a foundational defect in schizophrenia, hDAAO inhibitors have potential as treatments for schizophrenia and other nervous system disorders. Here, we sought to identify novel chemicals that inhibit hDAAO activity. We used computational tools to design a focused, purchasable library of compounds. After screening this library for hDAAO inhibition, we identified the structurally novel compound, ‘compound 2’ [3-(7-hydroxy-2-oxo-4-phenyl-2H-chromen-6-yl)propanoic acid], which displayed low nM hDAAO inhibitory potency (K_i=7 nM). Although the library was expected to enrich for compounds that were competitive for both D-serine and FAD, compound 2 actually was FAD uncompetitive, much like canonical hDAAO inhibitors such as benzoic acid. Compound 2 and an analog were independently co-crystalized with hDAAO. These compounds stabilized a novel conformation of hDAAO in which the active-site lid was in an open position. These results confirm previous hypotheses regarding active-site lid flexibility of mammalian D-amino acid oxidases and could assist in the design of the next generation of hDAAO inhibitors.     

Therapeutic Benefits of Induced Pluripotent Stem Cells in Monocrotaline-Induced Pulmonary Arterial Hypertension

Pulmonary arterial hypertension (PAH) is characterized by progressive increases in vascular resistance and the remodeling of pulmonary arteries. The accumulation of inflammatory cells in the lung and elevated levels of inflammatory cytokines in the bloodstream suggest that inflammation may play a role in PAH. In this study, the benefits of induced pluripotent stem cells (iPSCs) and iPSC-conditioned medium (iPSC CM) were explored in monocrotaline (MCT)-induced PAH rats. We demonstrated that both iPSCs and iPSC CM significantly reduced the right ventricular systolic pressure and ameliorated the hypertrophy of the right ventricle in MCT-induced PAH rats in models of both disease prevention and disease reversal. In the prevention of MCT-induced PAH, iPSC-based therapy led to the decreased accumulation of inflammatory cells and down-regulated the expression of the IL-1β, IL-6, IL-12α, IL-12β, IL-23 and IFNγ genes in lung specimens, which implied that iPSC-based therapy may be involved in the regulation of inflammation. NF-κB signaling is essential to the inflammatory cascade, which is activated via the phosphorylation of the NF-κB molecule. Using the chemical inhibitor specifically blocked the phosphorylation of NF-κB, and in vitro assays of cultured human M1 macrophages implied that the anti-inflammation effect of iPSC-based therapy may contribute to the disturbance of NF-κB activation. Here, we showed that iPSC-based therapy could restore the hemodynamic function of right ventricle with benefits for preventing the ongoing inflammation in the lungs of MCT-induced PAH rats by regulating NF-κB phosphorylation.     

The human telomeric proteome during telomere replication

Telomere shortening can cause detrimental diseases and contribute to aging. It occurs due to the end replication problem in cells lacking telomerase. Furthermore, recent studies revealed that telomere shortening can be attributed to difficulties of the semi-conservative DNA replication machinery to replicate the bulk of telomeric DNA repeats. To investigate telomere replication in a comprehensive manner, we develop QTIP-iPOND - Quantitative Telomeric chromatin Isolation Protocol followed by isolation of Proteins On Nascent DNA - which enables purification of proteins that associate with telomeres specifically during replication. In addition to the core replisome, we identify a large number of proteins that specifically associate with telomere replication forks. Depletion of several of these proteins induces telomere fragility validating their importance for telomere replication. We also find that at telomere replication forks the single strand telomere binding protein POT1 is depleted, whereas histone H1 is enriched. Our work reveals the dynamic changes of the telomeric proteome during replication, providing a valuable resource of telomere replication proteins. To our knowledge, this is the first study that examines the replisome at a specific region of the genome.     

Role for Heat Shock Protein 90α in the Proliferation and Migration of HaCaT Cells and in the Deep Second-Degree Burn Wound Healing in Mice

Inflammation, proliferation, and tissue remodeling are essential steps for wound healing. The hypoxic wound microenvironment promotes cell migration through a hypoxia—heat shock protein 90 alpha (Hsp90α)—low density lipoprotein receptor-related protein-1 (LRP-1) autocrine loop. To elucidate the role of this autocrine loop on burn wound healing, we investigated the expression profile of Hsp90α at the edge of burn wounds and found a transient increase in both mRNA and protein levels. Experiments performed with a human keratinocyte cell line—HaCaT also confirmed above results. 17-dimethylaminoethylamino-17demethoxygeldanamycin hydrochloride (17-DMAG), an Hsp90α inhibitor, was used to further evaluate the function of Hsp90α in wound healing. Consistently, topical application of Hsp90α in the early stage of deep second-degree burn wounds led to reduced inflammation and increased tissue granulation, with a concomitant reduction in the size of the wound at each time point tested (p<0.05). Consequently, epidermal cells at the wound margin progressed more rapidly causing an expedited healing process. In conclusion, these results provided a rationale for the therapeutic effect of Hsp90α on the burn wound management.     

Structure of the Type VI Effector-Immunity Complex (Tae4-Tai4) Provides Novel Insights into the Inhibition Mechanism of the Effector by Its Immunity Protein

The type VI secretion system (T6SS), a multisubunit needle-like apparatus, has recently been found to play a role in interspecies interactions. The Gram-negative bacteria harboring T6SS (donor) deliver the effectors into their neighboring cells (recipient) to kill them. Meanwhile, the cognate immunity proteins were employed to protect the donor cells against the toxic effectors. Tae4 (type VI amidase effector 4) and Tai4 (type VI amidase immunity 4) are newly identified T6SS effector-immunity pairs. Here, we report the crystal structures of Tae4 from Enterobacter cloacae and Tae4-Tai4 complexes from both E. cloacae and Salmonella typhimurium. Tae4 acts as a dl-endopeptidase and displays a typical N1pC/P60 domain. Unlike Tsi1 (type VI secretion immunity 1), Tai4 is an all-helical protein and forms a dimer in solution. The small angle x-ray scattering study combined with the analytical ultracentrifugation reveal that the Tae4-Tai4 complex is a compact heterotetramer that consists of a Tai4 dimer and two Tae4 molecules in solution. Structure-based mutational analysis of the Tae4-Tai4 interface shows that a helix (α3) of one subunit in dimeric Tai4 plays a major role in binding of Tae4, whereas a protruding loop (L4) in the other subunit is mainly responsible for inhibiting Tae4 activity. The inhibition process requires collaboration between the Tai4 dimer. These results reveal a novel and unique inhibition mechanism in effector-immunity pairs and suggest a new strategy to develop antipathogen drugs.     

Transport of glucose and galactose in kidney-cortex cells

1. The aerobic transport of d-glucose and d-galactose in rabbit kidney tissue at 25° was studied. 2. In slices forming glucose from added substrates an accumulation of glucose against its concentration gradient was found. The apparent ratio of intracellular ([S]_i) and extracellular ([S]_o) glucose concentrations was increased by 0·4mm-phlorrhizin and 0·3mm-ouabain. 3. Slices and isolated renal tubules actively accumulated glucose from the saline; the apparent [S]_i/[S]_o fell below 1·0 only at [S]_o higher than 0·5mm. 4. The rate of glucose oxidation by slices was characterized by the following parameters: K_m 1·16mm; V_max. 4·5μmoles/g. wet wt./hr. 5. The active accumulation of glucose from the saline was decreased by 0·1mm-2,4-dinitrophenol, 0·4mm-phlorrhizin and by the absence of external Na⁺. 6. The kinetic parameters of galactose entry into the cells were: K_m 1·5mm; V_max 10μmoles/g. wet wt./hr. 7. The efflux kinetics from slices indicated two intracellular compartments for d-galactose. The galactose efflux was greatly diminished at 0°, was inhibited by 0·4mm-phlorrhizin, but was insensitive to ouabain. 8. The following mechanism of glucose and galactose transport in renal tubular cells is suggested: (a) at the tubular membrane, these sugars are actively transported into the cells by a metabolically- and Na⁺-dependent phlorrhizin-sensitive mechanism; (b) at the basal cell membrane, these sugars are transported in accordance with their concentration gradient by a phlorrhizin-sensitive Na⁺-independent facilitated diffusion. The steady-state intracellular sugar concentration is determined by the kinetic parameters of active entry, passive outflow and intracellular utilization.     

Preferential Utilization of d-Lactate by Shewanella oneidensis

Shewanella oneidensis couples oxidation of lactate to respiration of many substrates. Here we report that llpR (l-lactate-positive regulator, SO_3460) encodes a positive regulator of l-lactate utilization distinct from previously studied regulators. We also demonstrate d-lactate inhibition of l-lactate utilization in S. oneidensis, resulting in preferential utilization of the d isomer.     

Isolation,cDNA Cloning,and Structure-based Functional Characterization of Oryctin,a Hemolymph Protein from the Coconut Rhinoceros Beetle,Oryctes rhinoceros,as a Novel Serine Protease Inhibitor

We isolated oryctin, a 66-residue peptide, from the hemolymph of the coconut rhinoceros beetle Oryctes rhinoceros and cloned its cDNA. Oryctin is dissimilar to any other known peptides in amino acid sequence, and its function has been unknown. To reveal that function, we determined the solution structure of recombinant ¹³C,¹⁵N-labeled oryctin by heteronuclear NMR spectroscopy. Oryctin exhibits a fold similar to that of Kazal-type serine protease inhibitors but has a unique additional C-terminal α-helix. We performed protease inhibition assays of oryctin against several bacterial and eukaryotic proteases. Oryctin does inhibit the following serine proteases: α-chymotrypsin, endopeptidase K, subtilisin Carlsberg, and leukocyte elastase, with K_i values of 3.9 × 10⁻¹⁰ m, 6.2 × 10⁻¹⁰ m, 1.4 × 10⁻⁹ m, and 1.2 × 10⁻⁸ m, respectively. Although the target molecule of oryctin in the beetle hemolymph remains obscure, our results showed that oryctin is a novel single domain Kazal-type inhibitor and could play a key role in protecting against bacterial infections.     

A Principled Approach to Deriving Approximate Conditional Sampling Distributions in Population Genetics Models with Recombination

The multilocus conditional sampling distribution (CSD) describes the probability that an additionally sampled DNA sequence is of a certain type, given that a collection of sequences has already been observed. The CSD has a wide range of applications in both computational biology and population genomics analysis, including phasing genotype data into haplotype data, imputing missing data, estimating recombination rates, inferring local ancestry in admixed populations, and importance sampling of coalescent genealogies. Unfortunately, the true CSD under the coalescent with recombination is not known, so approximations, formulated as hidden Markov models, have been proposed in the past. These approximations have led to a number of useful statistical tools, but it is important to recognize that they were not derived from, though were certainly motivated by, principles underlying the coalescent process. The goal of this article is to develop a principled approach to derive improved CSDs directly from the underlying population genetics model. Our approach is based on the diffusion process approximation and the resulting mathematical expressions admit intuitive genealogical interpretations, which we utilize to introduce further approximations and make our method scalable in the number of loci. The general algorithm presented here applies to an arbitrary number of loci and an arbitrary finite-alleles recurrent mutation model. Empirical results are provided to demonstrate that our new CSDs are in general substantially more accurate than previously proposed approximations.THE probability of observing a sample of DNA sequences under a given population genetics model—which is referred to as the sampling probability or likelihood—plays an important role in a wide range of problems in a genetic variation study. When recombination is involved, however, obtaining an analytic formula for the sampling probability has hitherto remained a challenging open problem (see Jenkins and Song 2009, 2010 for recent progress on this problem). As such, much research (Griffiths and Marjoram 1996; Kuhner et al. 2000; Nielsen 2000; Stephens and Donnelly 2000; Fearnhead and Donnelly 2001; De Iorio and Griffiths 2004a,b; Fearnhead and Smith 2005; Griffiths et al. 2008; Wang and Rannala 2008) has focused on developing Monte Carlo methods on the basis of the coalescent with recombination (Griffiths 1981; Kingman 1982a,b; Hudson 1983), a well-established mathematical framework that models the genealogical history of sample chromosomes. These Monte Carlo-based full-likelihood methods mark an important development in population genetics analysis, but a well-known obstacle to their utility is that they tend to be computationally intensive. For a whole-genome variation study, approximations are often unavoidable, and it is therefore important to think of ways to minimize the trade-off between scalability and accuracy.A popular likelihood-based approximation method that has had a significant impact on population genetics analysis is the following approach introduced by Li and Stephens (2003): Given a set Φ of model parameters (e.g., mutation rate, recombination rate, etc.), the joint probability p(h₁, … , h_n | Φ) of observing a set {h₁, … , h_n} of haplotypes sampled from a population can be decomposed as a product of conditional sampling distributions (CSDs), denoted by π,(1)where π(h_k+1|h₁, …, h_k, Φ) is the probability of an additionally sampled haplotype being of type h_k+1, given a set of already observed haplotypes h₁, …, h_k. In the presence of recombination, the true CSD π is unknown, so Li and Stephens proposed using an approximate CSD in place of π, thus obtaining the following approximation of the joint probability:(2)Li and Stephens referred to this approximation as the product of approximate conditionals (PAC) model. In general, the closer is to the true CSD π, the more accurate the PAC model becomes. Notable applications and extensions of this framework include estimating crossover rates (Li and Stephens 2003; Crawford et al. 2004) and gene conversion parameters (Gay et al. 2007; Yin et al. 2009), phasing genotype data into haplotype data (Stephens and Scheet 2005; Scheet and Stephens 2006), imputing missing data to improve power in association mapping (Stephens and Scheet 2005; Li and Abecasis 2006; Marchini et al. 2007; Howie et al. 2009), inferring local ancestry in admixed populations (Price et al. 2009), inferring human colonization history (Hellenthal et al. 2008), inferring demography (Davison et al. 2009), and so on.Another problem in which the CSD plays a fundamental role is importance sampling of genealogies under the coalescent process (Stephens and Donnelly 2000; Fearnhead and Donnelly 2001; De Iorio and Griffiths 2004a,b; Fearnhead and Smith 2005; Griffiths et al. 2008). In this context, the optimal proposal distribution can be written in terms of the CSD π (Stephens and Donnelly 2000), and as in the PAC model, an approximate CSD may be used in place of π. The performance of an importance sampling scheme depends critically on the proposal distribution and therefore on the accuracy of the approximation . Often in conjunction with composite-likelihood frameworks (Hudson 2001; Fearnhead and Donnelly 2002), importance sampling has been used in estimating fine-scale recombination rates (McVean et al. 2004; Fearnhead and Smith 2005; Johnson and Slatkin 2009).So far, a significant scope of intuition has gone into choosing the approximate CSDs used in these problems (Marjoram and Tavaré 2006). In the case of completely linked loci, Stephens and Donnelly (2000) suggested constructing an approximation by assuming that the additional haplotype h_k+1 is an imperfect copy of one of the first k haplotypes, with copying errors corresponding to mutation. Fearnhead and Donnelly (2001) generalized this construction to include crossover recombination, assuming that the haplotype h_k+1 is an imperfect mosaic of the first k haplotypes (i.e., h_k+1 is obtained by copying segments from h₁, …, h_k, where crossover recombination can change the haplotype from which copying is performed). The associated CSD, which we denote by , can be interpreted as a hidden Markov model and so admits an efficient dynamic programming solution. Finally, Li and Stephens (2003) proposed a modification to Fearnhead and Donnelly''s model that limits the hidden state space, thereby providing a computational simplification; we denote the corresponding approximate CSD by .Although these approaches are computationally appealing, it is important to note that they are not derived from, though are certainly motivated by, principles underlying typical population genetics models, in particular the coalescent process (Griffiths 1981; Kingman 1982a,b; Hudson 1983). The main objective of this article is to develop a principled technique to derive an improved CSD directly from the underlying population genetics model. Rather than relying on intuition, we base our work on mathematical foundation. The theoretical framework we employ is the diffusion process. De Iorio and Griffiths (2004a,b) first introduced the diffusion-generator approximation technique to obtain an approximate CSD in the case of a single locus (i.e., no recombination). Griffiths et al. (2008) later extended the approach to two loci to include crossover recombination, assuming a parent-independent mutation model at each locus. In this article, we extend the framework to develop a general algorithm that applies to an arbitrary number of loci and an arbitrary finite-alleles recurrent mutation model.Our work can be summarized as follows. Using the diffusion-generator approximation technique, we derive a recursion relation satisfied by an approximate CSD. This recursion can be used to construct a closed system of coupled linear equations, in which the conditional sampling probability of interest appears as one of the unknown variables. The system of equations can be solved using standard numerical analysis techniques. However, the size of the system grows superexponentially with the number of loci and, consequently, so does the running time. To remedy this drawback, we introduce additional approximations to make our approach scalable in the number of loci. Specifically, the recursion admits an intuitive genealogical interpretation, and, on the basis of this interpretation, we propose modifications to the recursion, which then can be easily solved using dynamic programming. The computational complexity of the modified algorithm is polynomial in the number of loci, and, importantly, the resulting CSD has little loss of accuracy compared to that following from the full recursion.The accuracy of approximate CSDs has not been discussed much in the literature, except in the application-specific context for which they are being employed. In this article, we carry out an empirical study to explicitly test the accuracy of various CSDs and demonstrate that our new CSDs are in general substantially more accurate than previously proposed approximations. We also consider the PAC framework and show that our approximations also produce more accurate PAC-likelihood estimates. We note that for the maximum-likelihood estimation of recombination rates, the actual value of the likelihood may not be so important, as long as it is maximized near the true recombination rate. However, in many other applications—e.g., phasing genotype data into haplotype data, imputing missing data, importance sampling, and so on—the accuracy of the CSD and PAC-likelihood function over a wide range of parameter values may be important. Thus, we believe that the theoretical work presented here will have several practical implications; our method can be applied in a wide range of statistical tools that use CSDs, improving their accuracy.The remainder of this article is organized as follows. To provide intuition for the ensuing mathematics, we first describe a genealogical process that gives rise to our CSD. Using our genealogical interpretation, we consider two additional approximations and relate these to previously proposed CSDs. Then, in the following section, we derive our CSD using the diffusion-generator approach and provide mathematical statements for the additional approximations; some interesting limiting behavior is also described there. This section is self-contained and may be skipped by the reader uninterested in mathematical details. Finally, in the subsequent section, we carry out a simulation study to compare the accuracy of various approximate CSDs and demonstrate that ours are generally the most accurate.     

Uptake of monosaccharides by guinea-pig cerebral-cortex slices

By the use of 1mm-iodoacetate to inhibit glycolysis in guinea-pig cerebral tissue slices, the kinetics of the uptake of monosaccharides on transfer of tissue from 0° to 37° were studied. d-Ribose, d-galactose, d-mannose, l-sorbose, and d-fructose showed diffusion kinetics, whereas 2-deoxy-d-glucose, d-glucose, d-arabinose and d-xylose showed saturation kinetics.     

Structural basis for the sequence-dependent effects of platinum–DNA adducts

The differences in efficacy and molecular mechanisms of platinum based anti-cancer drugs cisplatin (CP) and oxaliplatin (OX) have been hypothesized to be in part due to the differential binding affinity of cellular and damage recognition proteins to CP and OX adducts formed on adjacent guanines in genomic DNA. HMGB1a in particular exhibits higher binding affinity to CP-GG adducts, and the extent of discrimination between CP- and OX-GG adducts is dependent on the bases flanking the adducts. However, the structural basis for this differential binding is not known. Here, we show that the conformational dynamics of CP- and OX-GG adducts are distinct and depend on the sequence context of the adduct. Molecular dynamics simulations of the Pt-GG adducts in the TGGA sequence context revealed that even though the major conformations of CP- and OX-GG adducts were similar, the minor conformations were distinct. Using the pattern of hydrogen bond formation between the Pt–ammines and the adjacent DNA bases, we identified the major and minor conformations sampled by Pt–DNA. We found that the minor conformations sampled exclusively by the CP-GG adduct exhibit structural properties that favor binding by HMGB1a, which may explain its higher binding affinity to CP-GG adducts, while these conformations are not sampled by OX-GG adducts because of the constraints imposed by its cyclohexane ring, which may explain the negligible binding affinity of HMGB1a for OX-GG adducts in the TGGA sequence context. Based on these results, we postulate that the constraints imposed by the cyclohexane ring of OX affect the DNA conformations explored by OX-GG adduct compared to those of CP-GG adduct, which may influence the binding affinities of HMG-domain proteins for Pt-GG adducts, and that these conformations are further influenced by the DNA sequence context of the Pt-GG adduct.     

Amino ketone formation and aminopropanol-dehydrogenase activity in rat-liver preparations

1. Rat tissue homogenates convert dl-1-aminopropan-2-ol into aminoacetone. Liver homogenates have relatively high aminopropanol-dehydrogenase activity compared with kidney, heart, spleen and muscle preparations. 2. Maximum activity of liver homogenates is exhibited at pH9·8. The K_m for aminopropanol is approx. 15mm, calculated for a single enantiomorph, and the maximum activity is approx. 9mμmoles of aminoacetone formed/mg. wet wt. of liver/hr.at 37°. Aminoacetone is also formed from l-threonine, but less rapidly. An unidentified amino ketone is formed from dl-4-amino-3-hydroxybutyrate, the K_m for which is approx. 200mm at pH9·8. 3. Aminopropanol-dehydrogenase activity in homogenates is inhibited non-competitively by dl-3-hydroxybutyrate, the K_i being approx. 200mm. EDTA and other chelating agents are weakly inhibitory, and whereas potassium chloride activates slightly at low concentrations, inhibition occurs at 50–100mm. 4. It is concluded that aminopropanol-dehydrogenase is located in mitochondria, and in contrast with l-threonine dehydrogenase can be readily solubilized from mitochondrial preparations by ultrasonic treatment. 5. Soluble extracts of disintegrated mitochondria exhibit maximum aminopropanol-dehydrogenase activity at pH9·1 At this pH, K_m values for the amino alcohol and NAD⁺ are approx. 200 and 1·3mm respectively. Under optimum conditions the maximum velocity is approx. 70mμmoles of aminoacetone formed/mg. of protein/hr. at 37°. Chelating agents and thiol reagents appear to have little effect on enzyme activity, but potassium chloride inhibits at all concentrations tested up to 80mm. dl-3-Hydroxybutyrate is only slightly inhibitory. 6. Dehydrogenase activities for l-threonine and dl-4-amino-3-hydroxybutyrate appear to be distinct from that for aminopropanol. 7. Intraperitoneal injection of aminopropanol into rats leads to excretion of aminoacetone in the urine. Aminoacetone excretion proportional to the amount of the amino alcohol administered, is complete within 24hr., but represents less than 0·1% of the dose given. 8. The possible metabolic role of amino alcohol dehydrogenases is discussed.     

siah-1 Protein Is Necessary for High Glucose-induced Glyceraldehyde-3-phosphate Dehydrogenase Nuclear Accumulation and Cell Death in M��ller Cells

The translocation and accumulation of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) in the nucleus has closely been associated with cell death induction. However, the mechanism of this process has not been completely understood. The E3 ubiquitin ligase siah-1 (seven in absentia homolog 1) has recently been identified as a potential shuttle protein to transport GAPDH from the cytosol to the nucleus. Previously, we have demonstrated that elevated glucose levels induce GAPDH nuclear accumulation in retinal Müller cells. Therefore, this study investigated the role of siah-1 in high glucose-induced GAPDH nuclear translocation and subsequent cell death in retinal Müller cells. High glucose significantly increased siah-1 expression within 12 h. Under hyperglycemic conditions, siah-1 formed a complex with GAPDH and was predominantly localized in the nucleus of Müller cells. siah-1 knockdown using 50 nm siah-1 small interfering RNA significantly decreased high glucose-induced GAPDH nuclear accumulation at 24 h by 43.8 ± 4.0%. Further, knockdown of siah-1 prevented high glucose-induced cell death of Müller cells potentially by inhibiting p53 phosphorylation consistent with previous observations, indicating that nuclear GAPDH induces cell death via p53 activation. Therefore, inhibition of GAPDH nuclear translocation and accumulation by targeting siah-1 promotes Müller cell survival under hyperglycemic conditions.