An (O)over-bar(m(2)n) randomized algorithm to compute a minimum cycle basis of a directed graph

We consider the problem of computing a minimum cycle basis in a directed graph G. The input to this problem is a directed graph whose arcs have positive weights. In this problem a {- 1, 0, 1} incidence vector is associated with each cycle and the vector space over Q generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of weights of the cycles is minimum is called a minimum cycle basis of G. The current fastest algorithm for computing a minimum cycle basis in a directed graph with m arcs and n vertices runs in O(m(w+1)n) time (where w < 2.376 is the exponent of matrix multiplication). If one allows randomization, then an (O) over tilde (m(3)n) algorithm is known for this problem. In this paper we present a simple (O) over tilde (m(2)n) randomized algorithm for this problem. The problem of computing a minimum cycle basis in an undirected graph has been well-studied. In this problem a {0, 1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of the graph. The fastest known algorithm for computing a minimum cycle basis in an undirected graph runs in O(m(2)n + mn(2) logn) time and our randomized algorithm for directed graphs almost matches this running time.

Structural and molecular basis of interaction of HCV non-structural protein 5A with human casein kinase 1 alpha and PKR

Background: Interaction of non-structural protein 5A (NS5A) of Hepatitis C virus (HCV) with human kinases namely, casein kinase 1 alpha (ck1 alpha) and protein kinase R (PKR) have different functional implications such as regulation of viral replication and evasion of interferon induced immune response respectively. Understanding the structural and molecular basis of interactions of the viral protein with two different human kinases can be useful in developing strategies for treatment against HCV. Results: Serine 232 of NS5A is known to be phosphorylated by human ck1 alpha. A structural model of NS5A peptide containing phosphoacceptor residue Serine 232 bound to ck1 alpha has been generated using the known 3-D structures of kinase-peptide complexes. The substrate interacting residues in ck1 alpha has been identified from the model and these are found to be conserved well in the ck1 family. ck1 alpha - substrate peptide complex has also been used to understand the structural basis of association between ck1 alpha and its other viral stress induced substrate, tumour suppressor p53 transactivation domain which has a crystal structure available. Interaction of NS5A with another human kinase PKR is primarily genotype specific. NS5A from genotype 1b has been shown to interact and inhibit PKR whereas NS5A from genotype 2a/3a are unable to bind and inhibit PKR efficiently. This is one of the main reasons for the varied response to interferon therapy in HCV patients across different genotypes. Using PKR crystal structure, sequence alignment and evolutionary trace analysis some of the critical residues responsible for the interaction of NS5A 1b with PKR have been identified. Conclusions: The substrate interacting residues in ck1 alpha have been identified using the structural model of kinase substrate peptide. The PKR interacting NS5A 1b residues have also been predicted using PKR crystal structure, NS5A sequence analysis along with known experimental results. Functional significance and nature of interaction of interferon sensitivity determining region and variable region 3 of NS5A in different genotypes with PKR which was experimentally shown are also supported by the findings of evolutionary trace analysis. Designing inhibitors to prevent this interaction could enable the HCV genotype 1 infected patients respond well to interferon therapy.

Solution nuclear magnetic resonance structure of the GATase subunit and structural Basis of the interaction between GATase and ATPPase subunits in a two-subunit-type GMPS from methanocaldococcus jannaschii

The solution structure of the monomeric glutamine amidotransferase (GATase) subunit of the Methanocaldococcus janaschii (Mj) guanosine monophosphate synthetase (GMPS) has been determined using high-resolution nuclear magnetic resonance methods. Gel filtration chromatography and N-15 backbone relaxation studies have shown that the Mj GATase subunit is present in solution as a 21 kDa (188-residue) monomer. The ensemble of 20 lowest-energy structures showed root-mean-square deviations of 0.35 +/- 0.06 angstrom for backbone atoms and 0.8 +/- 0.06 angstrom for all heavy atoms. Furthermore, 99.4% of the backbone dihedral angles are present in the allowed region of the Ramachandran map, indicating the stereochemical quality of the structure. The core of the tertiary structure of the GATase is composed of a seven-stranded mixed beta-sheet that is fenced by five alpha-helices. The Mj GATase is similar in structure to the Pyrococcus horikoshi (Ph) GATase subunit. Nuclear magnetic resonance (NMR) chemical shift perturbations and changes in line width were monitored to identify residues on GATase that were responsible for interaction with magnesium and the ATPPase subunit, respectively. These interaction studies showed that a common surface exists for the metal ion binding as well as for the protein-protein interaction. The dissociation constant for the GATase-Mg2+ interaction has been found to be similar to 1 mM, which implies that interaction is very weak and falls in the fast chemical exchange regime. The GATase-ATPPase interaction, on the other hand, falls in the intermediate chemical exchange regime on the NMR time scale. The implication of this interaction in terms of the regulation of the GATase activity of holo GMPS is discussed.

Structural studies on a non-toxic homologue of type II RIPs from bitter gourd: Molecular basis of non-toxicity, conformational selection and glycan structure

The structures of nine independent crystals of bitter gourd seed lectin (BGSL), a non-toxic homologue of type II RIPs, and its sugar complexes have been determined. The four-chain, two-fold symmetric, protein is made up of two identical two-chain modules, each consisting of a catalytic chain and a lectin chain, connected by a disulphide bridge. The lectin chain is made up of two domains. Each domain carries a carbohydrate binding site in type II RIPs of known structure. BGSL has a sugar binding site only on one domain, thus impairing its interaction at the cell surface. The adenine binding site in the catalytic chain is defective. Thus, defects in sugar binding as well as adenine binding appear to contribute to the non-toxicity of the lectin. The plasticity of the molecule is mainly caused by the presence of two possible well defined conformations of a surface loop in the lectin chain. One of them is chosen in the sugar complexes, in a case of conformational selection, as the chosen conformation facilitates an additional interaction with the sugar, involving an arginyl residue in the loop. The N-glycosylation of the lectin involves a plant-specific glycan while that in toxic type II RIPs of known structure involves a glycan which is animal as well as plant specific.

Prediction of Indian rainfall during the summer monsoon season on the basis of links with equatorial Pacific and Indian Ocean climate indices

Interannual variation of Indian summer monsoon rainfall (ISMR) is linked to El Nino-Southern oscillation (ENSO) as well as the Equatorial Indian Ocean oscillation (EQUINOO) with the link with the seasonal value of the ENSO index being stronger than that with the EQUINOO index. We show that the variation of a composite index determined through bivariate analysis, explains 54% of ISMR variance, suggesting a strong dependence of the skill of monsoon prediction on the skill of prediction of ENSO and EQUINOO. We explored the possibility of prediction of the Indian rainfall during the summer monsoon season on the basis of prior values of the indices. We find that such predictions are possible for July-September rainfall on the basis of June indices and for August-September rainfall based on the July indices. This will be a useful input for second and later stage forecasts made after the commencement of the monsoon season.

The ecological basis of fishery yield of the Puerto Rico-Virgin Islands Insular Shelf: 1987 Assessment

A literature review was conducted to locate information on the flow of energy from primary producers to the fishery stocks of the Puerto Rican-Virgin Islands insular shelf. This report uses site-specific information to describe the major ecological subsystems, or habitats, of the region, to identify the more common species and the subsystems in which they occur, to quantify productivity and biomass, and to outline trophic relationships. Discussions on each topic and subsystem vary in substance and detail, being limited by the availability and accessibility of information. (PDF contains 189 pages) Seven distinct subsystems are described: mangrove estuary, seagrass bed, coral reef, algal plain, sand/mud bottom, shelf break, and overlying pelagic. Over 50 tables provide lists of species found in each habitat on various surveys dating back to 1956. Estimates of density, relative abundance, and productivity are provided when possible. We evaluated whether sufficient information exists to support an analysis of the energy basis of fishery production in the area, beginning with the design and development of an ecosystem model. Data needs in three categories - species lists, biomass, and trophic relations - were examined for each subsystem and for each of three species groups - primary producers, invertebrates, and fish. We concluded that adequate data, sufficient for modeling purposes, are available in 16 (25%) of 64 categories; limited data, those requiring greater extrapolation, are available in 35 (55%) categories; and no data are available in 13 (20%) categories. The best-studied subsystems are seagrass beds and coral reefs, with at least limited data in all categories. Invertebrates, the intermediate link in the food web between primary producers and fishes, are the least quantified group in the region. Primary production and fishes, however, are relatively well-studied, providing sufficient data to support an ecosystem-level analysis and to initiate a modeling effort.

The Molecular Basis of Lysine Acetylation: Addition, Removal, and Recognition

Acetyltransferases and deacetylases catalyze the addition and removal, respectively, of acetyl groups to the epsilon-amino group of protein lysine residues. This modification can affect the function of a protein through several means, including the recruitment of specific binding partners called acetyl-lysine readers. Acetyltransferases, deacetylases, and acetyl-lysine readers have emerged as crucial regulators of biological processes and prominent targets for the treatment of human disease. This work describes a combination of structural, biochemical, biophysical, cell-biological, and organismal studies undertaken on a set of proteins that cumulatively include all steps of the acetylation process: the acetyltransferase MEC-17, the deacetylase SIRT1, and the acetyl-lysine reader DPF2. Tubulin acetylation by MEC-17 is associated with stable, long-lived microtubule structures. We determined the crystal structure of the catalytic domain of human MEC-17 in complex with the cofactor acetyl-CoA. The structure in combination with an extensive enzymatic analysis of MEC-17 mutants identified residues for cofactor and substrate recognition and activity. A large, evolutionarily conserved hydrophobic surface patch distal to the active site was shown to be necessary for catalysis, suggesting that specificity is achieved by interactions with the alpha-tubulin substrate that extend outside of the modified surface loop. Experiments in C. elegans showed that while MEC-17 is required for touch sensitivity, MEC-17 enzymatic activity is dispensible for this behavior. SIRT1 deacetylates a wide range of substrates, including p53, NF-kappaB, FOXO transcription factors, and PGC-1-alpha, with roles in cellular processes ranging from energy metabolism to cell survival. SIRT1 activity is uniquely controlled by a C-terminal regulatory segment (CTR). Here we present crystal structures of the catalytic domain of human SIRT1 in complex with the CTR in an apo form and in complex with a cofactor and a pseudo-substrate peptide. The catalytic domain adopts the canonical sirtuin fold. The CTR forms a beta-hairpin structure that complements the beta-sheet of the NAD^+-binding domain, covering an essentially invariant, hydrophobic surface. A comparison of the apo and cofactor bound structures revealed conformational changes throughout catalysis, including a rotation of a smaller subdomain with respect to the larger NAD^+-binding subdomain. A biochemical analysis identified key residues in the active site, an inhibitory role for the CTR, and distinct structural features of the CTR that mediate binding and inhibition of the SIRT1 catalytic domain. DPF2 represses myeloid differentiation in acute myelogenous leukemia. Finally, we solved the crystal structure of the tandem PHD domain of human DPF2. We showed that DPF2 preferentially binds H3 tail peptides acetylated at Lys14, and binds H4 tail peptides with no preference for acetylation state. Through a structural and mutational analysis we identify the molecular basis of histone recognition. We propose a model for the role of DPF2 in AML and identify the DPF2 tandem PHD finger domain as a promising novel target for anti-leukemia therapeutics.

The ionic basis of frequency selectivity in hair cells of the bullfrog's sacculus

Hair cells from the bull frog's sacculus, a vestibular organ responding to substrate-borne vibration, possess electrically resonant membrane properties which maximize the sensitivity of each cell to a particular frequency of mechanical input. The electrical resonance of these cells and its underlying ionic basis were studied by applying gigohm-seal recording techniques to solitary hair cells enzymatically dissociated from the sacculus. The contribution of electrical resonance to frequency selectivity was assessed from microelectrode recordings from hair cells in an excised preparation of the sacculus.

Electrical resonance in the hair cell is demonstrated by damped membrane-potential oscillations in response to extrinsic current pulses applied through the recording pipette. This response is analyzed as that of a damped harmonic oscillator. Oscillation frequency rises with membrane depolarization, from 80-160 Hz at resting potential to asymptotic values of 200-250 Hz. The sharpness of electrical tuning, denoted by the electrical quality factor, Q_e, is a bell-shaped function of membrane voltage, reaching a maximum value around eight at a membrane potential slightly positive to the resting potential.

In whole cells, three time-variant ionic currents are activated at voltages more positive than -60 to -50 mV; these are identified as a voltage-dependent, non-inactivating Ca current (I_ca), a voltage-dependent, transient K current (I_a), and a Ca-dependent K current (I_c). The C channel is identified in excised, inside-out membrane patches on the basis of its large conductance (130-200 pS), its selective permeability to Kover Na or Cl, and its activation by internal Ca ions and membrane depolarization. Analysis of open- and closed-lifetime distributions suggests that the C channel can assume at least two open and three closed kinetic states.

Exposing hair cells to external solutions that inhibit the Ca or C conductances degrades the electrical resonance properties measured under current-clamp conditions, while blocking the A conductance has no significant effect, providing evidence that only the Ca and C conductances participate in the resonance mechanism. To test the sufficiency of these two conductances to account for electrical resonance, a mathematical model is developed that describes I_ca, I_c, and intracellular Ca concentration during voltage-clamp steps. I_ca activation is approximated by a third-order Hodgkin-Huxley kinetic scheme. Ca entering the cell is assumed to be confined to a small submembrane compartment which contains an excess of Ca buffer; Ca leaves this space with first-order kinetics. The Ca- and voltage-dependent activation of C channels is described by a five-state kinetic scheme suggested by the results of single-channel observations. Parameter values in the model are adjusted to fit the waveforms of I_ca and I_c evoked by a series of voltage-clamp steps in a single cell. Having been thus constrained, the model correctly predicts the character of voltage oscillations produced by current-clamp steps, including the dependencies of oscillation frequency and Q_e on membrane voltage. The model shows quantitatively how the Ca and C conductances interact, via changes in intracellular Ca concentration, to produce electrical resonance in a vertebrate hair cell.

Leukemia-Associated Mutations in Nucleophosmin Alter Recognition by CRM1: Molecular Basis of Aberrant Transport

Nucleophosmin (NPM) is a nucleocytoplasmic shuttling protein, normally enriched in nucleoli, that performs several activities related to cell growth. NPM mutations are characteristic of a subtype of acute myeloid leukemia (AML), where mutant NPM seems to play an oncogenic role. AML-associated NPM mutants exhibit altered subcellular traffic, being aberrantly located in the cytoplasm of leukoblasts. Exacerbated export of AML variants of NPM is mediated by the nuclear export receptor CRM1, and due, in part, to a mutationally acquired novel nuclear export signal (NES). To gain insight on the molecular basis of NPM transport in physiological and pathological conditions, we have evaluated the export efficiency of NPM in cells, and present new data indicating that, in normal conditions, wild type NPM is weakly exported by CRM1. On the other hand, we have found that AML-associated NPM mutants efficiently form complexes with CRM1HA (a mutant CRM1 with higher affinity for NESs), and we have quantitatively analyzed CRM1HA interaction with the NES motifs of these mutants, using fluorescence anisotropy and isothermal titration calorimetry. We have observed that the affinity of CRM1HA for these NESs is similar, which may help to explain the transport properties of the mutants. We also describe NPM recognition by the import machinery. Our combined cellular and biophysical studies shed further light on the determinants of NPM traffic, and how it is dramatically altered by AML-related mutations.