Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

Spatial Modeling of a Soil Fertility Index using Visible-Near-Infrared Spectra and Terrain Attributes

Our objective was to develop a methodology to predict soil fertility using visible near-infrared (vis-NIR) diffuse reflectance spectra and terrain attributes derived from a digital elevation model (DEM). Specifically, our aims were to: (i) assemble a minimum data set to develop a soil fertility index for sugarcane (Sarcharum officinarum L.) (SFI-SC) for biofuel production in tropical soils; (ii) construct a model to predict the SFI-SC using soil vis-NIR spectra and terrain attributes; and (iii) produce a soil fertility map for our study area and assess it by comparing it with a green vegetation index (GVI). The study area was 185 ha located in sao Paulo State, Brazil. In total, 184 soil samples were collected and analyzed for a range of soil chemical and physical properties. Their vis-NIR spectra were collected from 400 to 2500 nm. The Shuttle Radar Topographic Mission 3-arcsec (90-m resolution) DEM of the area was used to derive 17 terrain attributes. A minimum data set of soil properties was selected to develop the SFI-SC. The SFI-SC consisted of three classes: Class 1, the highly fertile soils; Class 2, the fertile soils; and Class 3, the least fertile soils. It was derived heuristically with conditionals and using expert knowledge. The index was modeled with the spectra and terrain data using cross-validated decision trees. The cross-validation of the model correctly predicted Class 1 in 75% of cases, Class 2 in 61%, and Class 3 in 65%. A fertility map was derived for the study area and compared with a map of the GVI. Our approach offers a methodology that incorporates expert knowledge to derive the SFI-SC and uses a versatile spectro-spatial methodology that may be implemented for rapid and accurate determination of soil fertility and better exploration of areas suitable for production.

Modeling soil organic carbon dynamics in Oxisols of Ibiruba (Brazil) with the Century Model

introduction of conservation practices in degraded agricultural land will generally recuperate soil quality, especially by increasing soil organic matter. This aspect of soil organic C (SOC) dynamics under distinct cropping and management systems can be conveniently analyzed with ecosystem models such as the Century Model. In this study, Century was used to simulate SOC stocks in farm fields of the Ibiruba region of north central Rio Grande do Sul state in Southern Brazil. The region, where soils are predominantly Oxisols, was originally covered with subtropical woodlands and grasslands. SOC dynamics was simulated with a general scenario developed with historical data on soil management and cropping systems beginning with the onset of agriculture in 1900. From 1993 to 2050, two contrasting scenarios based on no-tillage soil management were established: the status quo scenario, with crops and agricultural inputs as currently practiced in the region and the high biomass scenario with increased frequency of corn in the cropping system, resulting in about 80% higher biomass addition to soils. Century simulations were in close agreement with SOC stocks measured in 2005 in the Oxisols with finer texture surface horizon originally under woodlands. However, simulations in the Oxisols with loamy surface horizon under woodlands and in the grassland soils were not as accurate. SOC stock decreased from 44% to 50% in fields originally under woodland and from 20% to 27% in fields under grasslands with the introduction of intensive annual grain crops with intensive tillage and harrowing operations. The adoption of conservation practices in the 1980s led to a stabilization of SOC stocks followed by a partial recovery of native stocks. Simulations to 2050 indicate that maintaining status quo would allow SOC stocks to recover from 81% to 86% of the native stocks under woodland and from 80% to 91 % of the native stocks under grasslands. Adoption of a high biomass scenario would result in stocks from 75% to 95% of the original stocks under woodlands and from 89% to 102% in the grasslands by 2050. These simulations outcomes underline the importance of cropping system yielding higher biomass to further increase SOC content in these Oxisols. This application of the Century Model could reproduce general trends of SOC loss and recovery in the Oxisols of the Ibiruba region. Additional calibration and validation should be conducted before extensive usage of Century as a support tool for soil carbon sequestration projects in this and other regions can be recommended. (C) 2009 Elsevier B.V. All rights reserved.

Changes in the genetic structure of Bacteria and microbial activity in an agricultural soil amended with tannery sludge

The application of tannery sludge to soils is a form of recycling; however, few studies have examined the impacts of this practice on soil microbial properties. We studied effects of two applications (2006 and 2007) of tannery sludge (with a low chromium content) on the structure of the bacterial community and on the microbial activity of soils. We fertilized an agricultural area in Rolandia, Parana state, Brazil with different doses of sludge based on total N content, which ranged from 0 to 1200 kg N ha(-1). Sludge remained on the soil surface for three months before being plowed. Soils were sampled seven times during the experiment. Bacterial community structure, assessed by denaturing gradient gel electrophoresis (DGGE), was modified by the application of tannery sludge. Soon after the first application, there was clear separation between the bacterial communities in different treatments, such that each dose of sludge was associated with a specific community. These differences remained until 300 days after application and also after the second sludge application, but 666 days after the beginning of the experiment no differences were found in the bacterial communities of the lowest doses and the control. The principal response curve (PRC) analysis showed that the first sludge application strongly stimulated biological activity even 300 days after application. The second application also stimulated activity, but at a lower magnitude and for a shorter time, given that 260 days after the second application there was no difference in biological activity among treatments. PRC also showed that the properties most influenced by the application of tannery sludge were enzymatic activities related to N cycling (asparaginase and urease). The redundancy analysis (RDA) showed that tannery sludge`s influence on microbial activity is mainly related to increases in inorganic N and soil pH. Results showed that changes in the structure of the bacterial community in the studied soils were directly related to changes of their biological activity. (C) 2010 Elsevier Ltd. All rights reserved.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Relating MODIS vegetation index time-series with structure, light absorption and stem production of fast-growing Eucalyptus plantations

By allowing the estimation of forest structural and biophysical characteristics at different temporal and spatial scales, remote sensing may contribute to our understanding and monitoring of planted forests. Here, we studied 9-year time-series of the Normalized Difference Vegetation Index (NDVI) from the Moderate Resolution Imaging Spectroradiometer (MODIS) on a network of 16 stands in fast-growing Eucalyptus plantations in Sao Paulo State, Brazil. We aimed to examine the relationships between NDVI time-series spanning entire rotations and stand structural characteristics (volume, dominant height, mean annual increment) in these simple forest ecosystems. Our second objective was to examine spatial and temporal variations of light use efficiency for wood production, by comparing time-series of Absorbed Photosynthetically Active Radiation (APAR) with inventory data. Relationships were calibrated between the NDVI and the fractions of intercepted diffuse and direct radiation, using hemispherical photographs taken on the studied stands at two seasons. APAR was calculated from the NDVI time-series using these relationships. Stem volume and dominant height were strongly correlated with summed NDVI values between planting date and inventory date. Stand productivity was correlated with mean NDVI values. APAR during the first 2 years of growth was variable between stands and was well correlated with stem wood production (r(2) = 0.78). In contrast, APAR during the following years was less variable and not significantly correlated with stem biomass increments. Production of wood per unit of absorbed light varied with stand age and with site index. In our study, a better site index was accompanied both by increased APAR during the first 2 years of growth and by higher light use efficiency for stem wood production during the whole rotation. Implications for simple process-based modelling are discussed. (C) 2009 Elsevier B.V. All rights reserved.

HFE gene mutations in patients with primary iron overload: Is there a significant improvement in molecular diagnosis yield with HFE sequencing?

Rare HFE variants have been shown to be associated with hereditary hemochromatosis (HH), an iron overload disease. The low frequency of the HFE p.C282Y mutation in HH-affected Brazilian patients may suggest that other HFE-related mutations may also be implicated in the pathogenesis of HH in this population. The main aim was to screen for new HFE mutations in Brazilian individuals with primary iron overload and to investigate their relationship with HH. Fifty Brazilian patients with primary iron overload (transferrin saturation >50% in females and 60% in males) were selected. Subsequent bidirectional sequencing for each HFE exon was performed. The effect of HFE mutations on protein structure were analyzed by molecular dynamics simulation and free binding energy calculations. p.C282Y in homozygosis or in heterozygosis with p.H63D were the most frequent genotypic combinations associated with HH in our sample population (present in 17 individuals, 34%). Thirty-six (72.0%) out of the 50 individuals presented at least one HFE mutation. The most frequent genotype associated with HH was the homozygous p.C282Y mutation (n = 11, 22.0%). One novel mutation (p.V256I) was indentified in heterozygosis with the p.H63D mutation. In silico modeling analysis of protein behavior indicated that the p.V256I mutation does not reduce the binding affinity between HFE and beta 2-microglobulin ((beta 2M) in the same way the p.C282Y mutation does compared with the native HFE protein. In conclusion, screening of HFE through direct sequencing, as compared to p.C282Y/p.H63D genotyping, was not able to increase the molecular diagnosis yield of HH. The novel p.V256I mutation could not be implicated in the molecular basis of the HH phenotype, although its role cannot be completely excluded in HH-phenotype development. Our molecular modeling analysis can help in the analysis of novel, previously undescribed, HFE mutations. (C) 2010 Elsevier Inc. All rights reserved.

Molecular Modeling Suggests Cruzain Specificity for Peptide Primaquine Prodrugs

Chagas` disease, infection caused by the protozoan Trypanosoma cruzi, is an important, social and medical ailment in the Latin America. This disease is endemic in 21 countries, mostly Latin America countries, with more than 300,000 new cases every year and about 16-18 million infected people. Current therapy is not effective in the chronic phase of the disease. Thus, new and better drugs are urgently needed. In this sense, the in vitro activity of primaquine (PQ) was reported. Based on this, peptide prodrugs of primaquine containing dipeptides - lysine-arginine (LysArg), phenylalanine-alanine (PheAla) and phenylalanine-arginine (PheArg) -- as carriers, were designed to be selectively cleaved by cruzain, a specific cysteine protease of T. cruzi. The prodrugs have shown to be active against tripomastigote forms according to this order: LysArg-PQ> PheAla-PQ> PheArg-PQ. The molecular mechanism of action considered a probable nucleophilic attack of the catalytic residue of cruzain (Cys25) on the respective prodrug amide carbonyl carbon, releasing PQ. In order to test this hypothesis, molecular modeling studies were performed, physicochemical parameters and stereoelectronic features calculated by using the AM1 semi-empirical method suggest that the amide carbonyl carbon is favorable for cleavage, where the LysArg showed the most electronic reactive and sterically disposable, leading to the prodrug release and action. In addition, the docking study indicates the occurrence of specific interactions between prodrugs and the pockets S1 and S2 of cruzain through the dipeptides carriers, being the distance between cruzain Cys25 and the amide carbonyl group related to the biological activity of the prodrugs.

Rational Design and 3D-Pharmacophore Mapping of 5 `-Thiourea-Substituted alpha-Thymidine Analogues as Mycobacterial TMPK Inhibitors

Thymidine monophosphate kinase (TMPK) has emerged as an attractive target for developing inhibitors of Mycobacterium tuberculosis growth. In this study the receptor-independent (RI) 4D-QSAR formalism has been used to develop QSAR models and corresponding 3D-pharmacophores for a set of 5`-thiourea-substituted alpha-thymidine inhibitors. Models were developed for the entire training set and for a subset of the training set consisting of the most potent inhibitors. The optimized (RI) 4D-QSAR models are statistically significant (r(2) = 0.90, q(2) = 0.83 entire set, r(2) = 0.86, q(2) = 0.80 high potency subset) and also possess good predictivity based on test set predictions. The most and least potent inhibitors, in their respective postulated active conformations derived from the models, were docked in the active site of the TMPK crystallographic structure. There is a solid consistency between the 3D-pharmacophore sites defined by the QSAR models and interactions with binding site residues. This model identifies new regions of the inhibitors that contain pharmacophore sites, such as the sugar-pyrimidine ring structure and the region of the 5`-arylthiourea moiety. These new regions of the ligands can be further explored and possibly exploited to identify new, novel, and, perhaps, better antituberculosis inhibitors of TMPKmt. Furthermore, the 3D-pharmacophores defined by these models can be used as a starting point for future receptor-dependent antituberculosis drug design as well as to elucidate candidate sites for substituent addition to optimize ADMET properties of analog inhibitors.

Host-guest complexation of a nitroheterocyclic compound with cyclodextrins: a spectrofluorimetric and molecular modeling study

Nitroheterocyclic compounds (NC) were candidate drugs proposed for Chagas disease chemotherapy. In this study, we investigated the complexation of hydroxymethylnitrofurazone (NFOH), a potential antichagasic compound, with alpha-cyclodextrin (alpha-CD), beta-cyclodextrin (beta-CD), Hydroxypropyl-beta-cyclodextrin (HP-beta-CD), Dimethyl-beta-cyclodextrin (DM-beta-CD) and gamma-cyclodextrin (gamma-CD) by fluorescence spectroscopy and molecular modeling studies. Hildebrand-Benesi equation was used to calculate the formation constants of NFOH with cyclodextrins based on the fluorescence differences in the CDs solution. The complexing capacity of NFOH with different CDs was compared through the results of association constant according to the following order: DM-beta-CD > beta-CD > alpha-CD > HP-beta-CD > gamma-CD. Molecular modeling studies give support for the experimental assignments, in favor of the formation of an inclusion complex between cyclodextrins with NFOH. This is an important study to investigate the effects of different kinds of cyclodextrins on the inclusion complex formation with NFOH and to better characterize a potential formulations to be used as therapeutic options for the oral treatment of Chagas disease.

Design, synthesis, antimicrobial activity and molecular modeling studies of novel benzofuroxan derivatives against Staphylococcus aureus

Molecular modi. cation is a quite promising strategy in the design and development of drug analogs with better bioavailability, higher intrinsic activity and less toxicity. In the search of new leads with potential antimicrobial activity, a new series of 14 4-substituted [N`-(benzofuroxan-5-yl) methylene] benzohydrazides, nifuroxazide derivatives, were synthesized and tested against standard and multidrug-resistant Staphylococcus aureus strains. The selection of the substituent groups was based on physicochemical properties, such as hydrophobicity and electronic effect. These properties were also evaluated through the lipophilic and electrostatic potential maps, respectively, considering the compounds with better biological pro. le. Twelve compounds exhibited similar bacteriostatic activity against standard and multidrug-resistant strains. The most active compound was the 4-CF(3) substituted derivative, which presented a minimum inhibitory concentration (MIC) value of 14.6-13.1 mu g/mL, and a ClogP value of 1.87. The results highlight the benzofuroxan derivatives as potential leads for designing new future antimicrobial drug candidates. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular modeling studies and in vitro bioactivity evaluation of a set of novel 5-nitro-heterocyclic derivatives as anti-T. cruzi agents

In this study, in vitro anti-T. cruzi activity assays of nifuroxazide (NX) analogues, such as 5-nitro-2-furfuryliden and 5-nitro-2-theniliden derivatives, were performed. A molecular modeling approach was also carried out to relate the lipophilicity potential ( LP) property and biological activity data. The majority of the NX derivatives showed increased anti-T. cruzi activity in comparison to the reference drug, benznidazole (BZN). Additionally, the 5-nitro-2-furfuryliden derivatives presented better pharmacological profile than the 5-nitro-2-theniliden analogues. The LP maps and corresponding ClogP values indicate that there is an optimum lipophilicity value, which must be observed in the design of new potential anti-T. cruzi agents. (c) 2009 Elsevier Ltd. All rights reserved.

Crystal structure of 1,2-bis(4-fluoro-2-methylphenyl)ditelluride, [Te(C7H6F)](2)

C14H12F2Te2, monoclinic, C12/c1 (no. 15), a = 23.650(2) angstrom, b = 7.792(1) angstrom, c = 7.556(1) angstrom, beta = 106.47(1)degrees, V = 1335.2 angstrom(3), Z = 4, R-gt(F) = 0.041, wR(ref)(F-2) = 0.123, T= 98 K.

Crystal structure of Trypanosoma cruzi dihydroorotate dehydrogenase from Y strain

Trypanosoma cruzi is the etiological agent of Chagas` disease, a pathogenesis that affects millions of people in Latin America. Here, we report the crystal structure of dihydroorotate dehydrogenase (DHODH) from T cruzi strain Y solved at 2.2 angstrom resolution. DHODH is a flavin mononucleotide containing enzyme, which catalyses the oxidation Of L-dihydroorotate to orotate, the fourth step and only redox reaction in the de novo biosynthesis of pyrimidine nucleotides. Genetic studies have shown that DHODH is essential for T cruzi survival, validating the idea that this enzyme can be considered an attractive target for the development of antichagasic drugs. In our work, a detailed analysis of T cruzi DHODH crystal structure has allowed us to suggest potential sites to be further exploited for the design of highly specific inhibitors through the technology of structure-based drug design. (c) 2008 Elsevier Inc. All rights reserved.

HPLC separation, NMR and QTOF/MS/MS structure elucidation of a prominent nitric oxide donor agent based on an isomeric composition of a nitrosyl ruthenium complex

A new and promising nitrosyl ruthenium complex, [Ru(NO)(bdqi-COOH)(terpy)](PF(6))(3), bdqi-COOH is 3,4-diiminebenzoic acid and terpy is 2,2`-terpyridine, has been synthesized as a NO donor agent. The procedure used for [Ru(NO)(bdqi-COOH)(terpy)](PF(6))(3) synthesis has, apparently, yielded the formation of two isomers in which the ligand bdqi-COOH appears to be coordinated in its reduced form (bdcat-COOH), which could have differences in their pharmacological properties. Therefore, it was intended to separate the two possible isomers by high-performance liquid chromatography (HPLC) and to characterize them by high resolution mass spectrometry (QTOF MS) and by magnetic nuclear resonance spectroscopy (NMR). The results obtained by MS showed that the ESI-MS mass spectra of both HPLC column fractions, e.g. peak 1 and peak 2, are essentially equal, showing that both isomers display nearly identical gas-phase behavior with clusters of isotopologue ions centered at m/z 573, m/z 543 and m/z 513. Regarding the NMR analysis, the results showed that the positional isomerism is located in the bdqi-COOH ligand. From the observed results it can be concluded that the synthesis procedure that has been used results in the formation of two [Ru(terpy)(bdqi-COOH)NO](PF(6))(3) isomers. (c) 2009 Elsevier B.V. All rights reserved.