Population Dynamics Modeling of Arapaima gigas

Pirarucu (Arapaima gigas) has been of the most important natural fishing resources of the Amazon region. Due to its economic importance, and the necessity to preserve the species hand, field research concerning the habits and behavior of the pirarucu has been increasing for the last 20 years. The aim of this paper is to present a mathematical model for the pirarucu population dynamics considering the species peculiarities, particularly the male parental care over the offspring. The solution of the dynamical systems indicates three possible equilibrium points for the population. The first corresponds to extinction; the third corresponds to a stable population close to the environmental carrying capacity. The second corresponds to an unstable equilibrium located between extinction and full use of the carrying capacity. It is shown that lack of males’ parental care closes the gap between the point corresponding to the unstable equilibrium and the point of stable non-trivial equilibrium. If guarding failure reaches a critical point the two points coincide and the population tends irreversibly to extinction. If some event tends to destabilize the population equilibrium, as for instance inadequate parental care, the model responds in such a way as to restore the trajectory towards the stable equilibrium point avoiding the route to extinction. The parameters introduced to solve the system of equations are partially derived from limited but reliable field data collected at the Mamirauá Sustainable Development Reserve (MSDR) in the Brazilian Amazonian Region.

Neural Mechanisms and Delayed Gastric Emptying of Liquid Induced Through Acute Myocardial Infarction in Rats

Background: In pathological situations, such as acute myocardial infarction, disorders of motility of the proximal gut can trigger symptoms like nausea and vomiting. Acute myocardial infarction delays gastric emptying (GE) of liquid in rats. Objective: Investigate the involvement of the vagus nerve, α 1-adrenoceptors, central nervous system GABAB receptors and also participation of paraventricular nucleus (PVN) of the hypothalamus in GE and gastric compliance (GC) in infarcted rats. Methods: Wistar rats, N = 8-15 in each group, were divided as INF group and sham (SH) group and subdivided. The infarction was performed through ligation of the left anterior descending coronary artery. GC was estimated with pressure-volume curves. Vagotomy was performed by sectioning the dorsal and ventral branches. To verify the action of GABAB receptors, baclofen was injected via icv (intracerebroventricular). Intravenous prazosin was used to produce chemical sympathectomy. The lesion in the PVN of the hypothalamus was performed using a 1mA/10s electrical current and GE was determined by measuring the percentage of gastric retention (% GR) of a saline meal. Results: No significant differences were observed regarding GC between groups; vagotomy significantly reduced % GR in INF group; icv treatment with baclofen significantly reduced %GR. GABAB receptors were not conclusively involved in delaying GE; intravenous treatment with prazosin significantly reduced GR% in INF group. PVN lesion abolished the effect of myocardial infarction on GE. Conclusion: Gastric emptying of liquids induced through acute myocardial infarction in rats showed the involvement of the vagus nerve, alpha1- adrenergic receptors and PVN.

Modeling AIDS vaccines: the cellular level

This paper discuses current strategies for the development of AIDS vaccines wich allow immunzation to disturb the natural course of HIV at different detailed stages of its life cycle. Mathematical models describing the main biological phenomena (i.e. virus and vaccine induced T4 cell growth; virus and vaccine induced activation latently infected T4 cells; incremental changes immune response as infection progress; antibody dependent enhancement and neutralization of infection) and allowing for different vaccination strategies serve as a backgroud for computer simulations. The mathematical models reproduce updated information on the behavior of immune cells, antibody concentrations and free viruses. The results point to some controversial outcomes of an AIDS vaccine such as an early increase in virus concentration among vaccinated when compared to nonvaccinated individuals.

The stabilizing effects of the acquired immunity on the schistosomiasis transmission modeling - the sensitivity Analysis

A mathematical model is proposed to analyze the effects of acquired immunity on the transmission of schistosomiasis in the human host. From this model the prevalence curve dependent on four parameters can be obtained. These parameters were estimated fitting the data by the maximum likelihood method. The model showed a good retrieving capacity of real data from two endemic areas of schistosomiasis: Touros, Brazil (Schistosoma mansoni) and Misungwi, Tanzania (S. haematobium). Also, the average worm burden per person and the dispersion of parasite per person in the community can be obtained from the model. In this paper, the stabilizing effects of the acquired immunity assumption in the model are assessed in terms of the epidemiological variables as follows. Regarded to the prevalence curve, we calculate the confidence interval, and related to the average worm burden and the worm dispersion in the community, the sensitivity analysis (the range of the variation) of both variables with respect to their parameters is performed.

Distinction of males of the Lutzomyia intermedia (Lutz & Neiva, 1912) species complex by ratios between dimensions and by an artificial neural network (Diptera: Psychodidae, Phlebotominae)

The females of the two species of the Lutzomyia intermedia complex can be easily distinguished, but the males of each species are quite similar. The ratios between the extra-genital and the genital structures of L. neivai are larger than those of L. intermedia s. s., according to ANOVA. An artificial neural network was trained with a set of 300 examples, randomly taken from a sample of 358 individuals. The input vectors consisted of several ratios between some structures of each insect. The model was tested on the remaining 58 insects, 56 of which (96.6%) were correctly identified. This ratio of success can be considered remarkable if one takes into account the difficulty of attaining comparable results using traditional statistical techniques.

Molecular Modeling Approaches for Determining Gene Function: application to a Putative Poly-A Binding Protein from Leishmania amazonensis (LaPABP)

The great expansion in the number of genome sequencing projects has revealed the importance of computational methods to speed up the characterization of unknown genes. These studies have been improved by the use of three dimensional information from the predicted proteins generated by molecular modeling techniques. In this work, we disclose the structure-function relationship of a gene product from Leishmania amazonensis by applying molecular modeling and bioinformatics techniques. The analyzed sequence encodes a 159 aminoacids polypeptide (estimated 18 kDa) and was denoted LaPABP for its high homology with poly-A binding proteins from trypanosomatids. The domain structure, clustering analysis and a three dimensional model of LaPABP, basically obtained by homology modeling on the structure of the human poly-A binding protein, are described. Based on the analysis of the electrostatic potential mapped on the model's surface and conservation of intramolecular contacts responsible for folding stabilization we hypothesize that this protein may have less avidity to RNA than it's L. major counterpart but still account for a significant functional activity in the parasite. The model obtained will help in the design of mutagenesis experiments aimed to elucidate the mechanism of gene expression in trypanosomatids and serve as a starting point for its exploration as a potential source of targets for a rational chemotherapy.

Uncertainties regarding dengue modeling in Rio de Janeiro, Brazil

Dengue fever is currently the most important arthropod-borne viral disease in Brazil. Mathematical modeling of disease dynamics is a very useful tool for the evaluation of control measures. To be used in decision-making, however, a mathematical model must be carefully parameterized and validated with epidemiological and entomological data. In this work, we developed a simple dengue model to answer three questions: (i) which parameters are worth pursuing in the field in order to develop a dengue transmission model for Brazilian cities; (ii) how vector density spatial heterogeneity influences control efforts; (iii) with a degree of uncertainty, what is the invasion potential of dengue virus type 4 (DEN-4) in Rio de Janeiro city. Our model consists of an expression for the basic reproductive number (R0) that incorporates vector density spatial heterogeneity. To deal with the uncertainty regarding parameter values, we parameterized the model using a priori probability density functions covering a range of plausible values for each parameter. Using the Latin Hypercube Sampling procedure, values for the parameters were generated. We conclude that, even in the presence of vector spatial heterogeneity, the two most important entomological parameters to be estimated in the field are the mortality rate and the extrinsic incubation period. The spatial heterogeneity of the vector population increases the risk of epidemics and makes the control strategies more complex. At last, we conclude that Rio de Janeiro is at risk of a DEN-4 invasion. Finally, we stress the point that epidemiologists, mathematicians, and entomologists need to interact more to find better approaches to the measuring and interpretation of the transmission dynamics of arthropod-borne diseases.

Schistosoma mansoni major egg antigen Smp40: molecular modeling and potential immunoreactivity for anti-pathology vaccine development

The pathogenesis of Schistosoma mansoni infection is largely determined by host T-cell mediated immune responses such as the granulomatous response to tissue deposited eggs and subsequent fibrosis. The major egg antigens have a valuable role in desensitizing the CD4+ Th cells that mediate granuloma formation, which may prevent or ameliorate clinical signs of schistosomiasis.S. mansoni major egg antigen Smp40 was expressed and completely purified. It was found that the expressed Smp40 reacts specifically with anti-Smp40 monoclonal antibody in Western blotting. Three-dimensional structure was elucidated based on the similarity of Smp40 with the small heat shock protein coded in the protein database as 1SHS as a template in the molecular modeling. It was figured out that the C-terminal of the Smp40 protein (residues 130 onward) contains two alpha crystallin domains. The fold consists of eight beta strands sandwiched in two sheets forming Greek key. The purified Smp40 was used for in vitro stimulation of peripheral blood mononuclear cells from patients infected with S. mansoni using phytohemagglutinin mitogen as a positive control. The obtained results showed that there is no statistical difference in interferon-g, interleukin (IL)-4 and IL-13 levels obtained with Smp40 stimulation compared with the control group (P > 0.05 for each). On the other hand, there were significant differences after Smp40 stimulation in IL-5 (P = 0.006) and IL-10 levels (P < 0.001) compared with the control group. Gaining the knowledge by reviewing the literature, it was found that the overall pattern of cytokine profile obtained with Smp40 stimulation is reported to be associated with reduced collagen deposition, decreased fibrosis, and granuloma formation inhibition. This may reflect its future prospect as a leading anti-pathology schistosomal vaccine candidate.

Histopathological examination of nerve samples from pure neural leprosy patients: obtaining maximum information to improve diagnostic efficiency

Nerve biopsy examination is an important auxiliary procedure for diagnosing pure neural leprosy (PNL). When acid-fast bacilli (AFB) are not detected in the nerve sample, the value of other nonspecific histological alterations should be considered along with pertinent clinical, electroneuromyographical and laboratory data (the detection of Mycobacterium leprae DNA with polymerase chain reaction and the detection of serum anti-phenolic glycolipid 1 antibodies) to support a possible or probable PNL diagnosis. Three hundred forty nerve samples [144 from PNL patients and 196 from patients with non-leprosy peripheral neuropathies (NLN)] were examined. Both AFB-negative and AFB-positive PNL samples had more frequent histopathological alterations (epithelioid granulomas, mononuclear infiltrates, fibrosis, perineurial and subperineurial oedema and decreased numbers of myelinated fibres) than the NLN group. Multivariate analysis revealed that independently, mononuclear infiltrate and perineurial fibrosis were more common in the PNL group and were able to correctly classify AFB-negative PNL samples. These results indicate that even in the absence of AFB, these histopathological nerve alterations may justify a PNL diagnosis when observed in conjunction with pertinent clinical, epidemiological and laboratory data.

Desenvolvimento e análise de uma rede neural artificial para estimativa da erosividade da chuva para o Estado de São Paulo

O conhecimento do valor da erosividade da chuva (R) de determinada localidade é fundamental para a estimativa das perdas de solo feitas a partir da Equação Universal de Perdas de Solo, sendo, portanto, de grande importância no planejamento conservacionista. A fim de obter estimativas do valor de R para localidades onde este é desconhecido, desenvolveu-se uma rede neural artificial (RNA) e analisou-se a acurácia desta com o método de interpolação "Inverso de uma Potência da Distância" (ID). Comparando a RNA desenvolvida com o método de interpolação ID, verificou-se que a primeira apresentou menor erro relativo médio na estimativa de R e melhor índice de confiança, classificado como "Ótimo", podendo, portanto, ser utilizada no planejamento de uso, manejo e conservação do solo no Estado de São Paulo.

Modeling water movement in horizontal columns using fractal theory

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.

Remaining phosphorus estimated by pedotransfer function

Although the determination of remaining phosphorus (Prem) is simple, accurate values could also be estimated with a pedotransfer function (PTF) aiming at the additional use of soil analysis data and/or Prem replacement by an even simpler determination. The purpose of this paper was to develop a pedotransfer function to estimate Prem values of soils of the State of São Paulo based on properties with easier or routine laboratory determination. A pedotransfer function was developed by artificial neural networks (ANN) from a database of Prem values, pH values measured in 1 mol L-1 NaF solution (pH NaF) and soil chemical and physical properties of samples collected during soil classification activities carried out in the State of São Paulo by the Agronomic Institute of Campinas (IAC). Furthermore, a pedotransfer function was developed by regressing Prem values against the same predictor variables of the ANN-based PTF. Results showed that Prem values can be calculated more accurately with the ANN-based pedotransfer function with the input variables pH NaF values along with the sum of exchangeable bases (SB) and the exchangeable aluminum (Al3+) soil content. In addition, the accuracy of the Prem estimates by ANN-based PTF were more sensitive to increases in the experimental database size. Although the database used in this study was not comprehensive enough for the establishment of a definitive pedotrasnfer function for Prem estimation, results indicated the inclusion of Prem and pH NaF measurements among the soil testing evaluations as promising ind order to provide a greater database for the development of an ANN-based pedotransfer function for accurate Prem estimates from pH NaF, SB, and Al3+ values.

Numerical modeling of the effect of variation of boundary conditions on vadose zone hydraulic properties

An accurate estimation of hydraulic fluxes in the vadose zone is essential for the prediction of water, nutrient and contaminant transport in natural systems. The objective of this study was to simulate the effect of variation of boundary conditions on the estimation of hydraulic properties (i.e. water content, effective unsaturated hydraulic conductivity and hydraulic flux) in a one-dimensional unsaturated flow model domain. Unsaturated one-dimensional vertical water flow was simulated in a pure phase clay loam profile and in clay loam interlayered with silt loam distributed according to the third iteration of the Cantor Bar fractal object Simulations were performed using the numerical model Hydrus 1D. The upper and lower pressure heads were varied around average values of -55 cm for the near-saturation range. This resulted in combinations for the upper and lower constant head boundary conditions, respectively, of -50 and -60 cm, -40 and -70 cm, -30 and -80 cm, -20 and -90 cm, and -10 and -100 cm. For the drier range the average head between the upper and lower boundary conditions was set to -550 cm, resulting in the combinations -500 and -600 cm, -400 and -700 cm, -300 and -800 cm, -200 and -900 cm, and -100 and -1,000 cm, for upper and lower boundary conditions, respectively. There was an increase in water contents, fluxes and hydraulic conductivities with the increase in head difference between boundary conditions. Variation in boundary conditions in the pure phase and interlayered one-dimensional profiles caused significant deviations in fluxes, water contents and hydraulic conductivities compared to the simplest case (a head difference between the upper and lower constant head boundaries of 10 cm in the wetter range and 100 cm in the drier range).

Structural changes in latosols of the cerrado region: II - soil compressive behavior and modeling of additional compaction

Currently in Brazil, as in other parts of the world, the concern is great with the increase of degraded agricultural soil, which is mostly related to the occurrence of soil compaction. Although soil texture is recognized as a very important component in the soil compressive behaviors, there are few studies that quantify its influence on the structural changes of Latosols in the Brazilian Cerrado region. This study aimed to evaluate structural changes and the compressive behavior of Latosols in Rio Verde, Goiás, through the modeling of additional soil compaction. The study was carried out using five Latosols with very different textures, under different soil compaction levels. Water retention and soil compression curves, and bearing capacity models were determined from undisturbed samples collected on the B horizons. Results indicated that clayey and very clayey Latosols were more susceptible to compression than medium-textured soils. Soil compression curves at density values associate with edaphic functions were used to determine the beneficial pressure (σ b) , i.e., pressure with optimal water retention, and critical pressure (σcrMAC), i.e., pressure with macroporosity below critical levels. These pressure values were higher than the preconsolidation pressure (σp), and therefore characterized as additional compaction. Based on the compressive behavior of these Latosols, it can be concluded that the combined preconsolidation pressure, beneficial pressure and critical pressure allow a better understanding of compression processes of Latosols.