Ground-state functional for quantum spin chains with bond defects

We have developed a nonlocal functional of the exchange interaction for the ground-state energy of quantum spin chains described by the Heisenberg Hamiltonian. An alternating chain is used to obtain the correlation energy and a local unit-cell approximation is defined in the context of the density-functional theory. The agreement with our exact numerical data, for small chains, is significantly better than a previous formulation, even for chains with several ferromagnetic or antiferromagnetic bond defects. The results can be particularly relevant in the study of finite spin-1/2 Heisenberg chains, with exchange couplings changing, magnitude, or even sign, from bond-to-bond.

Protein preparation, crystallization and preliminary X-ray analysis of Trypanosoma cruzi nucleoside diphosphate kinase 1

The flagellated protozoan parasite Trypanosoma cruzi is the aetiological agent of Chagas disease. Nucleoside diphosphate kinases (NDPKs) are enzymes that are involved in energy management and nucleoside balance in the cell. T. cruzi TcNDPK1, a canonical isoform, was overexpressed in Escherichia coli as an N-terminally poly-His-tagged fusion protein and crystallized. Crystals grew after 72 h in 0.2 M MgCl(2), 20% PEG 3350. Data were collected to 3.5 angstrom resolution using synchrotron X-ray radiation at the National Synchrotron Light Laboratory (Campinas, Brazil). The crystals belonged to the trigonal space group P3, with unit-cell parameters a = b = 127.84, c = 275.49 angstrom. Structure determination is under way and will provide relevant information that may lead to the first step in rational drug design for the treatment of Chagas disease.

Chemical and in vitro study of the potential of 3-(aryl)-2-sulfanylpropenoic acids and their Zn(II) complexes as protective agents against cadmium toxicity

The free H(2)xspa ligands [xspa = pspa, Clpspa, tspa or fspa where p = 3-(phenyl), Clp = 3-(2-chlorophenyl), t = 3-(2-thienyl), f = 3-(2-furyl) and spa = 2-sulfanylpropenoato], their Zn(II) complexes of formula [HQ](2)[Zn(xspa)(2)] (HQ=diisopropylammonium) and the Cd(II) equivalents were prepared and characterized by elemental analysis and by IR, Raman and NMR ((1)H, (13)C) spectroscopy. X-Ray studies of the crystal structures of [HQ](2)[Zn(pspa)(2)], [HQ](2)[Zn(Clpspa)2], [HQ](2)[Zn(tspa)(2)] and [HQ](2)[Zn(fspa)(2)] show that the zinc atom is coordinated to two O atoms and two S atoms of the ligands in a distorted tetrahedral ZnO(2)S(2) environment. In the structures of [HQ](2)[Cd(pspa)(2)] and [HQ](2)[Cd(Clpspa)(2)] the cadmium atom is coordinated to three S atoms and two carboxylato O atoms of the ligands in a distorted trigonal bipyramidal environment. The interchange of ligands between Zn( II) and Cd( II) was studied by (113)Cd NMR spectroscopy. The in vitro protective effect of H(2)xspa and their Zn( II) complexes against Cd toxicity was investigated using the human hepatocarcinoma HepG2 cell line and the pig renal proximal tubule LLC-PK1 cell line. The incorporation of Zn( II) was found to be relevant in the case of H(2)pspa, with an increase observed in the cell viability of the LCC-PK1 cells with respect to the value for the free ligand.

The electromagnetics of light transmission through subwavelength slits in metallic films

By numerically calculating the relevant electromagnetic fields and charge current densities, we show how local charges and currents near subwavelength structures govern light transmission through subwavelength apertures in a real metal film. The illumination of a single aperture generates surface waves; and in the case of slits, generates them with high efficiency and with a phase close to - pi with respect to a reference standing wave established at the metal film front facet. This phase shift is due to the direction of induced charge currents running within the slit walls. The surface waves on the entrance facet interfere with the standing wave. This interference controls the profile of the transmission through slit pairs as a function of their separation. We compare the calculated transmission profile for a two-slit array to simple interference models and measurements [Phys. Rev. B 77(11), 115411 (2008)]. (C) 2011 Optical Society of America

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

Cross-Talk Between Mitochondria and NADPH Oxidase: Effects of Mild Mitochondrial Dysfunction on Angiotensin II-Mediated Increase in Nox Isoform Expression and Activity in Vascular Smooth Muscle Cells

Mitochondria and NADPH oxidase activation are concomitantly involved in pathogenesis of many vascular diseases. However, possible cross-talk between those ROS-generating systems is unclear. We induced mild mitochondrial dysfunction due to mitochondrial DNA damage after 24 h incubation of rabbit aortic smooth muscle (VSMC) with 250 ng/mL ethidium bromide (EtBr). VSMC remained viable and had 29% less oxygen consumption, 16% greater baseline hydrogen peroxide, and unchanged glutathione levels. Serum-stimulated proliferation was unaltered at 24 h. Although PCR amplification of several mtDNA sequences was preserved, D-Loop mtDNA region showed distinct amplification of shorter products after EtBr. Such evidence for DNA damage was further enhanced after angiotensin-II (AngII) incubation. Remarkably, the normally observed increase in VSMC membrane fraction NADPH oxidase activity after AngII was completely abrogated after EtBr, together with failure to upregulate Nox1 mRNA expression. Conversely, basal Nox4 mRNA expression increased 1.6-fold, while being unresponsive to AngII. Similar loss in AngII redox response occurred after 24 h antimycin-A incubation. Enhanced Nox4 expression was unassociated with endoplasmic reticulum stress markers. Protein disulfide isomerase, an NADPH oxidase regulator, exhibited increased expression and inverted pattern of migration to membrane fraction after EtBr. These results unravel functionally relevant cross-talk between mitochondria and NADPH oxidase, which markedly affects redox responses to AngII. Antioxid Redox Signal 11, 1265-1278.

LATENT RANK CHANGE DETECTION FOR ANALYSIS OF SPLICE-JUNCTION MICROARRAYS WITH NONLINEAR EFFECTS

Alternative splicing of gene transcripts greatly expands the functional capacity of the genome, and certain splice isoforms may indicate specific disease states such as cancer. Splice junction microarrays interrogate thousands of splice junctions, but data analysis is difficult and error prone because of the increased complexity compared to differential gene expression analysis. We present Rank Change Detection (RCD) as a method to identify differential splicing events based upon a straightforward probabilistic model comparing the over-or underrepresentation of two or more competing isoforms. RCD has advantages over commonly used methods because it is robust to false positive errors due to nonlinear trends in microarray measurements. Further, RCD does not depend on prior knowledge of splice isoforms, yet it takes advantage of the inherent structure of mutually exclusive junctions, and it is conceptually generalizable to other types of splicing arrays or RNA-Seq. RCD specifically identifies the biologically important cases when a splice junction becomes more or less prevalent compared to other mutually exclusive junctions. The example data is from different cell lines of glioblastoma tumors assayed with Agilent microarrays.

GPIomics: global analysis of glycosylphosphatidylinositol-anchored molecules of Trypanosoma cruzi

Glycosylphosphatidylinositol (GPI) anchoring is a common, relevant posttranslational modification of eukaryotic surface proteins. Here, we developed a fast, simple, and highly sensitive (high attomole-low femtomole range) method that uses liquid chromatography-tandem mass spectrometry (LC-MS(n)) for the first large-scale analysis of GPI-anchored molecules (i.e., the GPIome) of a eukaryote, Trypanosoma cruzi, the etiologic agent of Chagas disease. Our genome-wise prediction analysis revealed that approximately 12% of T. cruzi genes possibly encode GPI-anchored proteins. By analyzing the GPIome of T. cruzi insect-dwelling epimastigote stage using LC-MS(n), we identified 90 GPI species, of which 79 were novel. Moreover, we determined that mucins coded by the T. cruzi small mucin-like gene (TcSMUG S) family are the major GPI-anchored proteins expressed on the epimastigote cell surface. TcSMUG S mucin mature sequences are short (56-85 amino acids) and highly O-glycosylated, and contain few proteolytic sites, therefore, less likely susceptible to proteases of the midgut of the insect vector. We propose that our approach could be used for the high throughput GPIomic analysis of other lower and higher eukaryotes. Molecular Systems Biology 7 April 2009; doi:10.1038/msb.2009.13

Gold is not a Faradaic-Efficient Borohydride Oxidation Electrocatalyst: An Online Electrochemical Mass Spectrometry Study

Direct borohydride fuel cells are promising high energy density portable generators. However, their development remains limited by the complexity of the anodic reaction: The borohydride oxidation reaction (BOR) kinetics is slow and occurs at high overvoltages, while it may compete with the heterogeneous hydrolysis of BH(4)(-). Nevertheless, one usually admits that gold is rather inactive toward the heterogeneous hydrolysis of BH(4)(-) and presents some activity regarding the BOR, therefore yielding to the complete eight-electron BOR. In the present paper, by coupling online mass spectrometry to electrochemistry, we in situ monitored the H(2) yield during BOR experiments on sputtered gold electrodes. Our results show non-negligible H(2) generation on Au on the whole BOR potential range (0-0.8 V vs reversible hydrogen electrode), thus revealing that gold cannot be considered as a faradaic-efficient BOR electrocatalyst. We further propose a relevant reaction pathway for the BOR on gold that accounts for these findings.

In situ infrared (FTIR) study of the mechanism of the borohydride oxidation reaction

Early reports stated that Au was a catalyst of choice for the BOR because it would yield a near complete faradaic efficiency. However, it has recently been suggested that gold could yield to some extent the heterogeneous hydrolysis of BH(4)(-),therefore lowering the electron count per BH(4)(-), especially at low potential. Actually, the blur will exist regarding the BOR mechanism on Au as long as no physical proof regarding the reaction intermediates is not put forward. In that frame, in situ physical techniques like FTIR exhibit some interest to study the BOR. Consequently, in situ infrared reflectance spectroscopy measurements (SPAIRS technique) have been performed in 1 M NaOH/1 M NaBH(4) on a gold electrode with the aim to detect the intermediate species. We monitored several bands in B-H ((nu) over bar similar to 1180,1080 and 972 cm(-1)) and B-O bond regions ((nu) over bar =1325 and similar to 1425cm(-1)), which appear sequentially as a function of the electrode polarization. These absorption bands are assigned to BH(3), BH(2) and BO(2)(-) species. At the light of the experimental results, possible initial elementary steps of the BOR on gold electrode have been proposed and discussed according to the relevant literature data.

Self-organized distribution of periodicity and chaos in an electrochemical oscillator

We report a detailed numerical investigation of a prototype electrochemical oscillator, in terms of high-resolution phase diagrams for an experimentally relevant section of the control (parameter) space. The prototype model consists of a set of three autonomous ordinary differential equations which captures the general features of electrochemical oscillators characterized by a partially hidden negative differential resistance in an N-shaped current-voltage stationary curve. By computing Lyapunov exponents, we provide a detailed discrimination between chaotic and periodic phases of the electrochemical oscillator. Such phases reveal the existence of an intricate structure of domains of periodicity self-organized into a chaotic background. Shrimp-like periodic regions previously observed in other discrete and continuous systems were also observed here, which corroborate the universal nature of the occurrence of such structures. In addition, we have also found a structured period distribution within the order region. Finally we discuss the possible experimental realization of comparable phase diagrams.

Dynamic light scattering and optical absorption spectroscopy study of pH and temperature stabilities of the extracellular hemoglobin of Glossoscolex paulistus

The extracellular hemoglobin of Glossoscolex paulistus (HbGp) is constituted of subunits containing heme groups, monomers and trimers, and nonheme structures, called linkers, and the whole protein has a minimum molecular mass near 3.1 x 10(6) Da. This and other proteins of the same family are useful model systems for developing blood substitutes due to their extracellular nature, large size, and resistance to oxidation. HbGp samples were studied by dynamic light scattering (DLS). In the pH range 6.0-8.0, HbGp is stable and has a monodisperse size distribution with a z-average hydrodynamic diameter (D-h) of 27 +/- 1 nm. A more alkaline pH induced an irreversible dissociation process, resulting in a smaller D-h of 10 +/- 1 nm. The decrease in D-h suggests a complete hemoglobin dissociation. Gel filtration chromatography was used to show unequivocally the oligomeric dissociation observed at alkaline pH. At pH 9.0, the dissociation kinetics is slow, taking a minimum of 24 h to be completed. Dissociation rate constants progressively increase at higher pH, becoming, at pH 10.5, not detectable by DILS. Protein temperature stability was also pH-dependent. Melting curves for HbGp showed oligomeric dissociation and protein denaturation as a function of pH. Dissociation temperatures were lower at higher pH. Kinetic studies were also performed using ultraviolet-visible absorption at the Soret band. Optical absorption monitors the hemoglobin autoxidation while DLS gives information regarding particle size changes in the process of protein dissociation. Absorption was analyzed at different pH values in the range 9.0-9.8 and at two temperatures, 25 degrees C and 38 degrees C. At 25 degrees C, for pH 9.0 and 9.3, the kinetics monitored by ultraviolet-visible absorption presents a monoexponential behavior, whereas for pH 9.6 and 9.8, a biexponential behavior was observed, consistent with heme heterogeneity at more alkaline pH. The kinetics at 38 degrees C is faster than that at 25 degrees C and is biexponential in the whole pH range. DLS dissociation rates are faster than the autoxidation dissociation rates at 25 degrees C. Autoxiclation and dissociation processes are intimately related, so that oligomeric protein dissociation promotes the increase of autoxidation rate and vice versa. The effect of dissociation is to change the kinetic character of the autoxidation of hemes from monoexponential to biexponential, whereas the reverse change is not as effective. This work shows that DLS can be used to follow, quantitatively and in real time, the kinetics of changes in the oligomerization of biologic complex supramolecular systems. Such information is relevant for the development of mimetic systems to be used as blood substitutes.

Source populations in coastal crabs: parameters affecting egg production

Benthic marine invertebrates may form metapopulations connected via propagule dispersal. Conservation efforts often target potential source coastlines to indirectly benefit areas depending on allochthonous offspring production. Besides population density, adult size structure, sex ratio, brooding frequency and the proportion of breeding individuals may significantly influence the reproductive output of benthic populations, but these effects have seldom been tested. We used rocky shore crabs to assess the spatial variability of such parameters at relevant scales for conservation purposes and to test their consistency over 2 consecutive years; we then used the data to address whether bottom-up processes or biological interactions might explain the patterns observed. We decomposed egg production rates into their components for the 2 most abundant brachyuran species inhabiting the intertidal rocky habitat. Adult density and brooding frequency varied consistently among shores for both species and largely explained the overall spatial trends of egg production. Temporally consistent patterns also included among-shore differences in the size of ovigerous females of the grapsid Pachygrapsus transversus and between-bay differences in the fecundity of the spider crab Epialtus brasiliensis. Sex ratio was remarkably constant in both. We found no positive or negative correlations between adult density and brooding frequency to support either the existence of a component Allee effect (lack of mate encounters) or an effect of intra-specific competition. Likewise, shore-specific potential growth in P. transversus does not negatively correlate with frequency of ovigerous individuals, as would be expected under a critical balance between these 2 processes. The patterns observed suggest that bottom-up drivers may best explain spatial trends in the reproductive output of these species.

Simulation of sugarcane residue decomposition and aboveground growth

Due to the worldwide increase in demand for biofuels, the area cultivated with sugarcane is expected to increase. For environmental and economic reasons, an increasing proportion of the areas are being harvested without burning, leaving the residues on the soil surface. This periodical input of residues affects soil physical, chemical and biological properties, as well as plant growth and nutrition. Modeling can be a useful tool in the study of the complex interactions between the climate, residue quality, and the biological factors controlling plant growth and residue decomposition. The approach taken in this work was to parameterize the CENTURY model for the sugarcane crop, to simulate the temporal dynamics of aboveground phytomass and litter decomposition, and to validate the model through field experiment data. When studying aboveground growth, burned and unburned harvest systems were compared, as well as the effect of mineral fertilizer and organic residue applications. The simulations were performed with data from experiments with different durations, from 12 months to 60 years, in Goiana, TimbaA(0)ba and Pradpolis, Brazil; Harwood, Mackay and Tully, Australia; and Mount Edgecombe, South Africa. The differentiation of two pools in the litter, with different decomposition rates, was found to be a relevant factor in the simulations made. Originally, the model had a basically unlimited layer of mulch directly available for decomposition, 5,000 g m(-2). Through a parameter optimization process, the thickness of the mulch layer closer to the soil, more vulnerable to decomposition, was set as 110 g m(-2). By changing the layer of mulch at any given time available for decomposition, the sugarcane residues decomposition simulations where close to measured values (R (2) = 0.93), contributing to making the CENTURY model a tool for the study of sugarcane litter decomposition patterns. The CENTURY model accurately simulated aboveground carbon stalk values (R (2) = 0.76), considering burned and unburned harvest systems, plots with and without nitrogen fertilizer and organic amendment applications, in different climates and soil conditions.