Does using stepwise variable selection to build sequential path analysis models make sense?

Causal inference methods - mainly path analysis and structural equation modeling - offer plant physiologists information about cause-and-effect relationships among plant traits. Recently, an unusual approach to causal inference through stepwise variable selection has been proposed and used in various works on plant physiology. The approach should not be considered correct from a biological point of view. Here, it is explained why stepwise variable selection should not be used for causal inference, and shown what strange conclusions can be drawn based upon the former analysis when one aims to interpret cause-and-effect relationships among plant traits.

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.

EFFECT OF THERMAL TREATMENT ON PHENOLIC COMPOUNDS AND FUNCTIONALITY LINKED TO TYPE 2 DIABETES AND HYPERTENSION MANAGEMENT OF PERUVIAN AND BRAZILIAN BEAN CULTIVARS (PHASEOLUS VULGARIS L.) USING IN VITRO METHODS

The effect of thermal treatment on phenolic compounds and type 2 diabetes functionality linked to alpha-glucosidase and alpha-amylase inhibition and hypertension relevant angiotensin I-converting enzyme (ACE) inhibition were investigated in selected bean (Phaseolus vulgaris L,) cultivars from Peru and Brazil using in vitro models. Thermal processing by autoclaving decreased the total phenolic content in all cultivars, whereas the 1,1-diphenyl-2-picrylhydrazyl radical scavenging activity-linked antioxidant activity increased among Peruvian cultivars, alpha-Amylase and alpha-glucosidase inhibitory activities were reduced significantly after heat treatment (73-94% and 8-52%, respectively), whereas ACE inhibitory activity was enhanced (9-15%). Specific phenolic acids such as chlorogenic and caffeic acid increased moderately following thermal treatment (2-16% and 5-35%, respectively). No correlation was found between phenolic contents and functionality associated to antidiabetes and antihypertension potential, indicating that non phenolic compounds may be involved. Thermally processed bean cultivars are interesting sources of phenolic acids linked to high antioxidant activity and show potential for hypertension prevention.

Antioxidant capacity of the striped sunflower (Helianthus annuus L.) seed extracts evaluated by three in vitro methods

The antioxidant capacity of the striped sunflower seed cotyledon extracts, obtained by sequential extraction with different polarities of solvents, was evaluated by three different in vitro methods: ferric reducing/antioxidant power (FRAP), 2.2-diphenyl-1-picrylhydrazyl radical (DPPH) and oxygen radical absorbance capacity (ORAC) assays. In the three methods, the aqueous extract at 30 mu g/ml showed a higher antioxidant capacity value (FRAP, 45.27 mu mol; DPPH, 50.18%; ORAC, 1.5 Trolox equivalents) than the ethanolic extract (FRAP, 32.17 mu mol; DPPH, 15.21%; ORAC, 0.50 Trolox equivalents). When compared with the synthetic antioxidant butylated hydroxyl toluene, the antioxidant capacity of the aqueous extract varied from 45% to 66%, according to the used method. The high antioxidant capacity observed for the aqueous extract of the studied sunflower seed suggests that the intake of this seed may prevent in vivo oxidative reactions responsible for the development of several diseases, such as cancer.

Methods for detection of anatoxin-a(s) by liquid chromatography coupled to electrospray ionization-tandem mass spectrometry

Anatoxin-a(s) is a potent irreversible inhibitor of the enzyme acetylcholinesterase with a unique N-hydroxyguanidine methylphosphate ester chemical structure. Determination of this toxin in environmental samples is hampered by the lack of specific methods for its detection. Using the toxic strain of Anabaena lemmermani PH-160 B as positive control, the fragmentation characteristics of anatoxin-a(s) under collision-induced dissociation conditions have been investigated and new LC-MS/MS methods proposed. Recommended ion transitions for correct detection of this toxin are 253 > 58, 253 > 159, 235 > 98 and 235 > 96. Chromatographic separation is better achieved under HILIC conditions employing a ZIC-HILIC column. This method was used to confirm for the first time the production of anatoxin-a(s) by strains of Anabaena oumiana ITEP-025 and ITEP-026. Considering no standard solutions are commercially available, our results will be of significant use for the correct identification of this toxin by LC-MS/MS. (C) 2009 Elsevier Ltd. 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.

QSAR Modeling of a Set of Pyrazinoate Esters as Antituberculosis Prodrugs

Tuberculosis is an infection caused mainly by Mycobacterium tuberculosis. A first-line antimycobacterial drug is pyrazinamide (PZA), which acts partially as a prodrug activated by a pyrazinamidase releasing the active agent, pyrazinoic acid (POA). As pyrazinoic acid presents some difficulty to cross the mycobacterial cell wall, and also the pyrazinamide-resistant strains do not express the pyrazinamidase, a set of pyrazinoic acid esters have been evaluated as antimycobacterial agents. In this work, a QSAR approach was applied to a set of forty-three pyrazinoates against M. tuberculosis ATCC 27294, using genetic algorithm function and partial least squares regression (WOLF 5.5 program). The independent variables selected were the Balaban index (I), calculated n-octanol/water partition coefficient (ClogP), van-der-Waals surface area, dipole moment, and stretching-energy contribution. The final QSAR model (N = 32, r(2) = 0.68, q(2) = 0.59, LOF = 0.25, and LSE = 0.19) was fully validated employing leave-N-out cross-validation and y-scrambling techniques. The test set (N = 11) presented an external prediction power of 73%. In conclusion, the QSAR model generated can be used as a valuable tool to optimize the activity of future pyrazinoic acid esters in the designing of new antituberculosis agents.

Molecular modeling as a promising tool to study dendrimer prodrugs delivery

Molecular modeling methodologies were applied to perform preliminary studies concerning the release of active agents from potentially antichagasic and antileishmanial dendrimer prodrugs. The dendrimer was designed having myo-inositol as a core, L-malic acid as a spacer group, and hydroxymethylnitrofurazone (NFOH), 3-hydroxyflavone or quercetin, as active compounds. Each dendrimer presented a particular behavior concerning to the following investigated properties: spatial hindrance, map of electrostatic potential (MEP), and the lowest unoccupied molecular orbital energy (E(LUMO)). Additionally, the findings suggested that the carbonyl group next to the active agent seems to be the most promising ester breaking point. (C) 2009 Elsevier B.V. All rights reserved.

Development and validation of sensitive methods for determination of sibutramine hydrochloride monohydrate and direct enantiomeric separation on a protein-based chiral stationary phase

Sibutramine hydrochloride monohydrate, chemically 1-(4-chlorophenyl)-N,N-dimethyl-alpha-(2-methylpropyl) hydrochloride monohydrate (SB center dot HCl center dot H2O), was approved by the U.S. Food and Drug Administration for the treatment of obesity. The objective of this study was to develop, validate, and compare methods using UV-derivative spectrophotometry (UVDS) and reversed-phase high-performance liquid chromatography (HPLC) for the determination of SB center dot HCl center dot H2O in pharmaceutical drug products. The UVDS and HPLC methods were found to be rapid, precise, and accurate. Statistically, there was no significant difference between the proposed UVDS and HPLC methods. The enantiomeric separation of SB was obtained on an alpha-1 acid glycoprotein column. The R- and S-sibutramine were eluted in < 5 min with baseline separation of the chromatographic peaks (alpha = 1.9 and resolution = 1.9).

Molecular modeling and QSAR studies of a set of indole and benzimidazole derivatives as H(4) receptor antagonists

Histamine is an important biogenic amine, which acts with a group of four G-protein coupled receptors (GPCRs), namely H(1) to H(4) (H(1)R - H(4)R) receptors. The actions of histamine at H(4)R are related to immunological and inflammatory processes, particularly in pathophysiology of asthma, and H(4)R ligands having antagonistic properties could be helpful as antiinflammatory agents. In this work, molecular modeling and QSAR studies of a set of 30 compounds, indole and benzimidazole derivatives, as H(4)R antagonists were performed. The QSAR models were built and optimized using a genetic algorithm function and partial least squares regression (WOLF 5.5 program). The best QSAR model constructed with training set (N = 25) presented the following statistical measures: r (2) = 0.76, q (2) = 0.62, LOF = 0.15, and LSE = 0.07, and was validated using the LNO and y-randomization techniques. Four of five compounds of test set were well predicted by the selected QSAR model, which presented an external prediction power of 80%. These findings can be quite useful to aid the designing of new anti-H(4) compounds with improved biological response.

Determination of Vecuronium Bromide in Pharmaceuticals: Development, Validation and Comparative Study of HPLC and CZE Analytical Methods

Vecuronium bromide is a neuromuscular blocking agent used for anesthesia to induce skeletal muscle relaxation. HPLC and CZE analytical methods were developed and validated for the quantitative determination of vecuronium bromide. The HPLC method was achieved on an amino column (Luna 150 x 4.6 mm, 5 mu m) using UV detection at 205 nm. The mobile phase was composed of acetonitrile:water containing 25.0 mmol L(-1) of sodium phosphate monobasic (50:50 v/v), pH 4.6 and flow rate of 1.0 mL min(-1). The CZE method was achieved on an uncoated fused-silica capillary (40.0 cm total length, 31.5 cm effective length and 50 mu m i.d.) using indirect UV detection at 230 nm. The electrolyte comprised 1.0 mmol L(-1) of quinine sulfate dihydrate at pH 3.3 and 8.0% of acetonitrile. The results were used to compare both techniques. No significant differences were observed (p > 0.05).

New Rapid Derivative Spectrophotometric and Chromatographic Methods for Assay of Loratadine in Tablets and Syrups

New rapid first-derivative spectrophotometric (UVDS) and a stability-indicating high performance liquid chromatographic (HPLC) methods were developed, validated and successfully applied in the analysis of loratadine (LT) in tablets and syrups. In the UVDS method, 0.1 M HCl was used as solvent. The measurements were made at 312.4 nm in the first order derivative spectra. The HPLC method was carried out on a RP-18 column with a mobile phase composed of methanol-water-tetrahydrofuran (50:30:20, v/v/v). UV detection was made at 247 nm. For HPLC methods the total analysis time was <3min, adequate for routine quality control of tablets and syrups containing loratadine.