Combining data sets of satellite-retrieved products for basin-scale water balance study: 2. Evaluation on the Mississippi Basin and closure correction model

In this study, we applied the integration methodology developed in the companion paper by Aires (2014) by using real satellite observations over the Mississippi Basin. The methodology provides basin-scale estimates of the four water budget components (precipitation P, evapotranspiration E, water storage change Delta S, and runoff R) in a two-step process: the Simple Weighting (SW) integration and a Postprocessing Filtering (PF) that imposes the water budget closure. A comparison with in situ observations of P and E demonstrated that PF improved the estimation of both components. A Closure Correction Model (CCM) has been derived from the integrated product (SW+PF) that allows to correct each observation data set independently, unlike the SW+PF method which requires simultaneous estimates of the four components. The CCM allows to standardize the various data sets for each component and highly decrease the budget residual (P - E - Delta S - R). As a direct application, the CCM was combined with the water budget equation to reconstruct missing values in any component. Results of a Monte Carlo experiment with synthetic gaps demonstrated the good performances of the method, except for the runoff data that has a variability of the same order of magnitude as the budget residual. Similarly, we proposed a reconstruction of Delta S between 1990 and 2002 where no Gravity Recovery and Climate Experiment data are available. Unlike most of the studies dealing with the water budget closure at the basin scale, only satellite observations and in situ runoff measurements are used. Consequently, the integrated data sets are model independent and can be used for model calibration or validation.

Enriching the annotation of Mycobacterium tuberculosis H37Rv proteome using remote homology detection approaches: Insights into structure and function

The availability of the genome sequence of Mycobacterium tuberculosis H37Rv has encouraged determination of large numbers of protein structures and detailed definition of the biological information encoded therein; yet, the functions of many proteins in M. tuberculosis remain unknown. The emergence of multidrug resistant strains makes it a priority to exploit recent advances in homology recognition and structure prediction to re-analyse its gene products. Here we report the structural and functional characterization of gene products encoded in the M. tuberculosis genome, with the help of sensitive profile-based remote homology search and fold recognition algorithms resulting in an enhanced annotation of the proteome where 95% of the M. tuberculosis proteins were identified wholly or partly with information on structure or function. New information includes association of 244 proteins with 205 domain families and a separate set of new association of folds to 64 proteins. Extending structural information across uncharacterized protein families represented in the M. tuberculosis proteome, by determining superfamily relationships between families of known and unknown structures, has contributed to an enhancement in the knowledge of structural content. In retrospect, such superfamily relationships have facilitated recognition of probable structure and/or function for several uncharacterized protein families, eventually aiding recognition of probable functions for homologous proteins corresponding to such families. Gene products unique to mycobacteria for which no functions could be identified are 183. Of these 18 were determined to be M. tuberculosis specific. Such pathogen-specific proteins are speculated to harbour virulence factors required for pathogenesis. A re-annotated proteome of M. tuberculosis, with greater completeness of annotated proteins and domain assigned regions, provides a valuable basis for experimental endeavours designed to obtain a better understanding of pathogenesis and to accelerate the process of drug target discovery. (C) 2014 Elsevier Ltd. All rights reserved.

Paralogous cAMP Receptor Proteins in Mycobacterium smegmatis Show Biochemical and Functional Divergence

The cyclic AMP receptor protein (CRP) family of transcription factors consists of global regulators of bacterial gene expression. Here, we identify two paralogous CRPs in the genome of Mycobacterium smegmatis that have 78% identical sequences and characterize them biochemically and functionally. The two proteins (MSMEG_0539 and MSMEG_6189) show differences in cAMP binding affinity, trypsin sensitivity, and binding to a CRP site that we have identified upstream of the msmeg_3781 gene. MSMEG_6189 binds to the CRP site readily in the absence of cAMP, while MSMEG_0539 binds in the presence of cAMP, albeit weakly. msmeg_6189 appears to be an essential gene, while the ?msmeg_0539 strain was readily obtained. Using promoter-reporter constructs, we show that msmeg_3781 is regulated by CRP binding, and its transcription is repressed by MSMEG_6189. Our results are the first to characterize two paralogous and functional CRPs in a single bacterial genome. This gene duplication event has subsequently led to the evolution of two proteins whose biochemical differences translate to differential gene regulation, thus catering to the specific needs of the organism.

Dendritic atrophy constricts functional maps in resonance and impedance properties of hippocampal model neurons

A gradient in the density of hyperpolarization-activated cyclic-nucleotide gated (HCN) channels is necessary for the emergence of several functional maps within hippocampal pyramidal neurons. Here, we systematically analyzed the impact of dendritic atrophy on nine such functional maps, related to input resistance and local/transfer impedance properties, using conductance-based models of hippocampal pyramidal neurons. We introduced progressive dendritic atrophy in a CA1 pyramidal neuron reconstruction through a pruning algorithm, measured all functional maps in each pruned reconstruction, and arrived at functional forms for the dependence of underlying measurements on dendritic length. We found that, across frequencies, atrophied neurons responded with higher efficiency to incoming inputs, and the transfer of signals across the dendritic tree was more effective in an atrophied reconstruction. Importantly, despite the presence of identical HCN-channel density gradients, spatial gradients in input resistance, local/transfer resonance frequencies and impedance profiles were significantly constricted in reconstructions with dendrite atrophy, where these physiological measurements across dendritic locations converged to similar values. These results revealed that, in atrophied dendritic structures, the presence of an ion channel density gradient alone was insufficient to sustain homologous functional maps along the same neuronal topograph. We assessed the biophysical basis for these conclusions and found that this atrophy-induced constriction of functional maps was mediated by an enhanced spatial spread of the influence of an HCN-channel cluster in atrophied trees. These results demonstrated that the influence fields of ion channel conductances need to be localized for channel gradients to express themselves as homologous functional maps, suggesting that ion channel gradients are necessary but not sufficient for the emergence of functional maps within single neurons.

Study of structures, energies and vibrational frequencies of (O-2)(n)(+) (n=2-5) clusters by GGA and meta-GGA density functional methods

Using Generalized Gradient Approximation (GGA) and meta-GGA density functional methods, structures, binding energies and harmonic vibrational frequencies for the clusters O-4(+), O-6(+), O-8(+) and O-10(+) have been calculated. The stable structures of O-4(+), O-6(+), O-8(+) and O-10(+) have point groups D-2h, D-3h, D-4h, and D-5h optimized on the quartet, sextet, octet and dectet potential energy surfaces, respectively. Rectangular (D-2h) O-4(+) has been found to be more stable compared to trans-planar (C-2h) on the quartet potential energy surface. Cyclic structure (D-3h) of CA cluster ion has been calculated to be more stable than other structures. Binding energy (B.E.) of the cyclic O-6(+) is in good agreement with experimental measurement. The zero-point corrected B.E. of O-8(+) with D4h symmetry on the octet potential energy surface and zero-point corrected B.E. of O-10(+) with D-5h symmetry on the dectet potential energy surface are also in good agreement with experimental values. The B.E. value for O-4(+) is close to the experimental value when single point energy is calculated by Brueckner coupled-cluster method, BD(T). (C) 2014 Elsevier B.V. All rights reserved.

Structural and functional analysis of two universal stress proteins YdaA and YnaF from Salmonella typhimurium: possible roles in microbial stress tolerance

In many organisms ``Universal Stress Proteins'' CUSPS) are induced in response to a variety of environmental stresses. Here we report the structures of two USPs, YnaF and YdaA from Salmonella typhimurium determined at 1.8 angstrom and 2.4 angstrom resolutions, respectively. YnaF consists of a single USP domain and forms a tetrameric organization stabilized by interactions mediated through chloride ions. YdaA is a larger protein consisting of two tandem USP domains. Two protomers of YdaA associate to form a structure similar to the YnaF tetramer. YdaA showed ATPase activity and an ATP binding motif G-2X-G-9X-G(S/T/N) was found in its C-terminal domain. The residues corresponding to this motif were not conserved in YnaF although YnaF could bind ATP. However, unlike YdaA, YnaF did not hydrolyse ATP in vitro. Disruption of interactions mediated through chloride ions by selected mutations converted YnaF into an ATPase. Residues that might be important for ATP hydrolysis could be identified by comparing the active sites of native and mutant structures. Only the C-terminal domain of YdaA appears to be involved in ATP hydrolysis. The structurally similar N-terminal domain was found to bind a zinc ion near the segment equivalent to the phosphate binding loop of the C-terminal domain. Mass spectrometric analysis showed that YdaA might bind a ligand of approximate molecular weight 800 daltons. Structural comparisons suggest that the ligand, probably related to an intermediate in lipid A biosynthesis, might bind at a site close to the zinc ion. Therefore, the N-terminal domain of YdaA binds zinc and might play a role in lipid metabolism. Thus, USPs appear to perform several distinct functions such as ATP hydrolysis, altering membrane properties and chloride sensing. (C) 2015 Elsevier Inc. All rights reserved.

Metallophosphoesterases: structural fidelity with functional promiscuity

Calcineurin-like metallophosphoesterases (MPEs) form a large superfamily of binuclear metal-ion-centre-containing enzymes that hydrolyse phosphomono-, phosphodi-or phosphotri-esters in a metal-dependent manner. The MPE domain is found in Mre11/SbcD DNA-repair enzymes, mammalian phosphoprotein phosphatases, acid sphingomyelinases, purple acid phosphatases, nucleotidases and bacterial cyclic nucleotide phosphodiesterases. Despite this functional diversity, MPEs show a remarkably similar structural fold and active-site architecture. In the present review, we summarize the available structural, biochemical and functional information on these proteins. We also describe how diversification and specialization of the core MPE fold in various MPEs is achieved by amino acid substitution in their active sites, metal ions and regulatory effects of accessory domains. Finally, we discuss emerging roles of these proteins as non-catalytic protein-interaction scaffolds. Thus we view the MPE superfamily as a set of proteins with a highly conserved structural core that allows embellishment to result in dramatic and niche-specific diversification of function.

Molecular and Functional Characterization of RecD, a Novel Member of the SF1 Family of Helicases, from Mycobacterium tuberculosis

Background: The heterotrimeric M. tuberculosis RecBCD complex, or each of its individual subunits, remains uncharacterized. Results: MtRecD exists as a homodimer in solution, catalyzes ssDNA-dependent ATP hydrolysis, unwinding of DNA replication/recombination intermediates, and interacts with RecA. Conclusion: MtRecD possesses strong 5 3- and weak 3 5-helicase activities. Significance: These findings provide insights into the mechanism underlying DSB repair and homologous recombination in mycobacteria. The annotated whole-genome sequence of Mycobacterium tuberculosis revealed the presence of a putative recD gene; however, the biochemical characteristics of its encoded protein product (MtRecD) remain largely unknown. Here, we show that MtRecD exists in solution as a stable homodimer. Protein-DNA binding assays revealed that MtRecD binds efficiently to single-stranded DNA and linear duplexes containing 5 overhangs relative to the 3 overhangs but not to blunt-ended duplex. Furthermore, MtRecD bound more robustly to a variety of Y-shaped DNA structures having 18-nucleotide overhangs but not to a similar substrate containing 5-nucleotide overhangs. MtRecD formed more salt-tolerant complexes with Y-shaped structures compared with linear duplex having 3 overhangs. The intrinsic ATPase activity of MtRecD was stimulated by single-stranded DNA. Site-specific mutagenesis of Lys-179 in motif I abolished the ATPase activity of MtRecD. Interestingly, although MtRecD-catalyzed unwinding showed a markedly higher preference for duplex substrates with 5 overhangs, it could also catalyze significant unwinding of substrates containing 3 overhangs. These results support the notion that MtRecD is a bipolar helicase with strong 5 3 and weak 3 5 unwinding activities. The extent of unwinding of Y-shaped DNA structures was approximate to 3-fold lower compared with duplexes with 5 overhangs. Notably, direct interaction between MtRecD and its cognate RecA led to inhibition of DNA strand exchange promoted by RecA. Altogether, these studies provide the first detailed characterization of MtRecD and present important insights into the type of DNA structure the enzyme is likely to act upon during the processes of DNA repair or homologous recombination.

Functional characterization of heat-shock protein 90 from Oryza sativa and crystal structure of its N-terminal domain

Heat-shock protein 90 (Hsp90) is an ATP-dependent molecular chaperone that is essential for the normal functioning of eukaryotic cells. It plays crucial roles in cell signalling, cell-cycle control and in maintaining proteome integrity and protein homeostasis. In plants, Hsp90s are required for normal plant growth and development. Hsp90s are observed to be upregulated in response to various abiotic and biotic stresses and are also involved in immune responses in plants. Although there are several studies elucidating the physiological role of Hsp90s in plants, their molecular mechanism of action is still unclear. In this study, biochemical characterization of an Hsp90 protein from rice (Oryza sativa; OsHsp90) has been performed and the crystal structure of its N-terminal domain (OsHsp90-NTD) was determined. The binding of OsHsp90 to its substrate ATP and the inhibitor 17-AAG was studied by fluorescence spectroscopy. The protein also exhibited a weak ATPase activity. The crystal structure of OsHsp90-NTD was solved in complex with the nonhydrolyzable ATP analogue AMPPCP at 3.1 angstrom resolution. The domain was crystallized by cross-seeding with crystals of the N-terminal domain of Hsp90 from Dictyostelium discoideum, which shares 70% sequence identity with OsHsp90-NTD. This is the second reported structure of a domain of Hsp90 from a plant source.

Structural plasticity in Mycobacterium tuberculosis uracil-DNA glycosylase (MtUng) and its functional implications

17 independent crystal structures of family I uracil-DNA glycosylase from Mycobacterium tuberculosis (MtUng) and its complexes with uracil and its derivatives, distributed among five distinct crystal forms, have been determined. Thermodynamic parameters of binding in the complexes have been measured using isothermal titration calorimetry. The two-domain protein exhibits open and closed conformations, suggesting that the closure of the domain on DNA binding involves conformational selection. Segmental mobility in the enzyme molecule is confined to a 32-residue stretch which plays a major role in DNA binding. Uracil and its derivatives can bind to the protein in two possible orientations. Only one of them is possible when there is a bulky substituent at the 50 position. The crystal structures of the complexes provide a reasonable rationale for the observed thermodynamic parameters. In addition to providing fresh insights into the structure, plasticity and interactions of the protein molecule, the results of the present investigation provide a platform for structure-based inhibitor design.

Insights into the Functional Roles of N-Terminal and C-Terminal Domains of Helicobacter pylori DprA

DNA processing protein A (DprA) plays a crucial role in the process of natural transformation. This is accomplished through binding and subsequent protection of incoming foreign DNA during the process of internalization. DprA along with Single stranded DNA binding protein A (SsbA) acts as an accessory factor for RecA mediated DNA strand exchange. H. pylori DprA (HpDprA) is divided into an N-terminal domain and a C-terminal domain. In the present study, individual domains of HpDprA have been characterized for their ability to bind single stranded (ssDNA) and double stranded DNA (dsDNA). Oligomeric studies revealed that HpDprA possesses two sites for dimerization which enables HpDprA to form large and tightly packed complexes with ss and dsDNA. While the N-terminal domain was found to be sufficient for binding with ss or ds DNA, C-terminal domain has an important role in the assembly of poly-nucleoprotein complex. Using site directed mutagenesis approach, we show that a pocket comprising positively charged amino acids in the N-terminal domain has an important role in the binding of ss and dsDNA. Together, a functional cross talk between the two domains of HpDprA facilitating the binding and formation of higher order complex with DNA is discussed.

Solid-state optical absorption from optimally tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for the accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce the known band gap, this method accurately predicts band structures and optical spectra of silicon and lithium fluoride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

Morphological and functional characterisation of the burrow systems of six earthworm species (Lumbricidae)

Earthworm burrow systems are generally described based on postulated behaviours associated with the three ecological types. In this study, we used X-ray tomography to obtain 3D information on the burrowing behaviour of six very common anecic (Aporrectodea nocturna and Lumbricus terrestris) and endogeic (Aporrectodea rosea, Allolobophora chlorotica, Aporrectodea caliginosa, Aporrectodea icterica) earthworm species, introduced into repacked soil cores for 6 weeks. A simple water infiltration test, the Beerkan method, was also used to assess some functional properties of these burrow systems. Endogeic worms make larger burrow systems, which are more highly branched, less continuous and of smaller diameter, than those of anecic worms. Among the anecic species, L. terrestris burrow systems are shorter (9.2 vs 21.2 m) with a higher number (14.5 vs 23.5) of less branched burrows (12.2 vs 20.2 branches m(-1)), which are also wider (7.78 vs 5.16 mm) than those of A. nocturna. In comparison, the burrow systems made by endogeic species appeared similar to each other. However, A. rosea burrows were short and narrow, whereas A. icterica had a longer burrow system (15.7 m), more intense bioturbation intensity (refilled macropores or soil lateral compaction around them) and thus a greater number of burrows. Regarding water infiltration, anecic burrow systems were far more efficient due to open burrows linking the top and bottom of the cores. For endogeic species, we observed a linear relationship between burrow length and the water infiltration rate (R (2) = 0.49, p < 0.01). Overall, the three main characteristics significantly influencing water infiltration were burrow length, burrow number and bioturbation volume. This last characteristic highlighted the effect of burrow refilling by casts.

Evaluation of Precipitation Retrievals From Orbital Data Products of TRMM Over a Subtropical Basin in India

The spatial error structure of daily precipitation derived from the latest version 7 (v7) tropical rainfall measuring mission (TRMM) level 2 data products are studied through comparison with the Asian precipitation highly resolved observational data integration toward evaluation of the water resources (APHRODITE) data over a subtropical region of the Indian subcontinent for the seasonal rainfall over 6 years from June 2002 to September 2007. The data products examined include v7 data from the TRMM radiometer Microwave Imager (TMI) and radar precipitation radar (PR), namely, 2A12, 2A25, and 2B31 (combined data from PR and TMI). The spatial distribution of uncertainty from these data products were quantified based on performance metrics derived from the contingency table. For the seasonal daily precipitation over a subtropical basin in India, the data product of 2A12 showed greater skill in detecting and quantifying the volume of rainfall when compared with the 2A25 and 2B31 data products. Error characterization using various error models revealed that random errors from multiplicative error models were homoscedastic and that they better represented rainfall estimates from 2A12 algorithm. Error decomposition techniques performed to disentangle systematic and random errors verify that the multiplicative error model representing rainfall from 2A12 algorithm successfully estimated a greater percentage of systematic error than 2A25 or 2B31 algorithms. Results verify that although the radiometer derived 2A12 rainfall data is known to suffer from many sources of uncertainties, spatial analysis over the case study region of India testifies that the 2A12 rainfall estimates are in a very good agreement with the reference estimates for the data period considered.