Cis-acting elements and developmental regulators govern toxin gene expression in Bacillus anthracis

Expression of the structural genes for the anthrax toxin proteins is coordinately controlled by host-related signals such as elevated CO2 , and the trans-acting positive regulator, AtxA. No specific binding of AtxA to the toxin gene promoters has been demonstrated and no sequence-based similarities are apparent in the promoter regions of toxin genes. We hypothesized that the toxin genes possess common structural features that are required for positive regulation. To test this hypothesis, I performed an extensive characterization of the toxin gene promoters. I determined the minimal sequences required for atxA-mediated toxin gene expression and compared these sequences for structural similarities. In silico modeling and in vitro experiments indicated significant curvature within these regions. Random mutagenesis revealed that point mutations associated with reduced transcriptional activity, mostly mapped to areas of high curvature. This work enabled the identification of two potential cis-acting elements implicated in AtxA-mediated regulation of the toxin genes. In addition to the growth condition requirements and AtxA, toxin gene expression is under growth phase regulation. The transition state regulator AbrB represses atxA expression to influence toxin synthesis. Here I report that toxin gene expression also requires sigH, a gene encoding the RNA polymerase sigma factor associated with development in B. subtilis. In the well-studied B. subtilis system, σH is part of a feedback control pathway that involves AbrB and the major response regulator of sporulation initiation, Spo0A. My data indicate that in B. anthracis, regulatory relationships exist between these developmental regulators and atxA . Interestingly, during growth in toxin-inducing conditions, sigH and abrB expression deviates from that described for B. subtilis, affecting expression of the atxA gene. These findings, combined with previous observations, suggest that the steady state level of atxA expression is critical for optimal toxin gene transcription. I propose a model whereby, under toxin-inducing conditions, control of toxin gene expression is fine-tuned by the independent effects of the developmental regulators on the expression of atxA . The growth condition-dependent changes in expression of these regulators may be crucial for the correct timing and uninterrupted expression of the toxin genes during infection. ^

Ser88 phosphorylation regulates dynein light chain 1 (DLC1) function in the mouse mammary gland

Dynein light chain 1 (DLC1) is a highly conserved and ubiquitously expressed protein which might have critical cellular function as total loss of DLC1 caused Drosophila embryonic death. Despite many proteins and RNAs interaction with it identified, DLC1's function(s) and regulation are largely unknown. Recently, DLC1 was identified as a physiological substrate of P21-activate kinase 1(Pak1) kinase from a human mammary cDNA library in a yeast-2-hybridization screening assay. Studies in primary human tumors and cell culture implicated that DLC1 could promote mammary cancerous phenotypes, and more importantly, Ser88 phosphorylation of DLC1by Pak1 kinase was found to be essential for DLC1's tumorigenic activities. Based on the above tissue culture studies, we hypothesized that Ser88 phosphorylation regulates DLC1. ^ To test this hypothesis, we generated two transgenic mouse models: MMTV-DLC1 and MMTV-DLC1-S88A mice with mammary specific expression of the DLC1 and DLC1-S88A cDNAs. Both of the transgenic mice mammary glands showed rare tumor incidence which indicated DLC1 alone may not be sufficient for tumorigenesis in vivo. However, these mice showed a significant alteration of mammary development. Mammary glands from the MMTV-DLC1 mice had hyperbranching and alveolar hyperplasia, with elevated cell proliferation. Intriguingly, these phenotypes were not seen in the mammary glands from the MMTV-S88A mice. Furthermore, while MMTV-DLC1 glands were normal during involution, MMTV-S88A mice showed accelerated mammary involution with increase apoptosis and altered expression of involution-associated genes. Further analysis of the MMTV-S88A glands showed they had increased steady state level of Bim protein which might be responsible for the early involution. Finally, our in vitro data showed that Ser88 phosphorylation abolished DLC1 dimer and consequently might disturb its interaction with Bim and destabilize Bim. ^ Collectively, our findings provided in vivo evidence that Ser88 phosphorylation of DLC1 can regulate DLC1's function. In addition, Ser88 phosphorylation might be critical for DLC1 dimer-monomer transition. ^

2-methylthio-N$\sp6$-dimethylallyl modification of adenosine 37 of tRNAs in {\it Escherichia coli\/} K-12

The hypermodified, hydrophobic 2-methylthio-N$\sp6$-(dimethylallyl)-adenosine (ms${2{\cdot}6}\atop1$A) residue occurs $3\sp\prime$ to the anticodon in tRNA species that read codons beginning with U. The first step (i$\sp6$A37 formation) of this modification is catalyzed by dimethylallyl diphosphate:tRNA dimethyallyltransferase (EC 2.5.1.8), which is the product of the miaA gene. Subsequent steps were proposed to be catalyzed by MiaB and MiaC enzymes to complete the ms${2{\cdot}6}\atop1$A37 modification. The study of functions of the ms${2{\cdot}6}\atop1$A37 is very important because this modified base is one of the best candidates for a role in global control in response to environmental stress. This dissertation describes the further delineation of functions of the ms${2{\cdot}6}\atop1$A37 modification in E. coli K-12 cells. This work provides significant information on functions of tRNA modifications in E. coli cells to adapt to stressful environmental conditions. Three hypotheses were tested in this work.^ The first hypothesis tested was that non-optimal translation processes cause increased spontaneous mutagenesis by the induction of SOS response in starving cells. To test this hypothesis, I measured spontaneous mutation rates of wild type cells and various mutant strains which are defective in tRNA modification, SOS response, or oxidative damage repair. I found that the miaA mutation acts as a mutator that increased Lac$\sp+$ reversion rates and Trp$\sp+$ reversion frequencies of the wild-type cells in starving conditions. However, the lexA3(Ind)(which abolishes the induction of SOS response) mutation abolished the mutator phenotype of the miaA mutant. The recA430 mutation, not other identified SOS genes, decreased the Lac$\sp+$ reversion to a less extent than that of the lexA3(Ind) mutation. These results suggest that RecA together with another unidentified SOS gene product are responsible for the process.^ The second hypothesis tested was that MiaA protein binds to full-length tRNA$\sp{\rm Phe}$ molecules in form of a protein dimer. To test this hypothesis, three versions of the MiaA protein and seven species of tRNA substrates were purified. Binding studies by gel mobility shift assays, filter binding assays and gel filtration shift assays support the hypothesis that MiaA protein binds to full-length tRNA$\sp{\rm Phe}$ as a protein dimer but as a monomer to the anticodon stem-and-loop. These results were further supported by using steady state enzyme kinetic studies.^ The third hypothesis tested in this work was that the miaB gene in E. coli exists and is clonable. The miaB::Tn10dCm insertion mutation of Salmonella typhimurium was transduced to E. coli K-12 cells by using P$\sb1$ and P$\sb{22}$ bacteriophages. The insertion was confirmed by HPLC analyses of nucleotide profiles of miaB mutants of E. coli. The insertion mutation was cloned and DNA sequences adjacent to the transposon were sequenced. These DNA sequences were 86% identical to the f474 gene at 14.97 min chromosome of E. coli. The f474 gene was then cloned by PCR from the wild-type chromosome of E. coli. The recombinant plasmid complemented the mutant phenotype of the miaB mutant of E. coli. These results support the hypothesis that the miaB gene of E. coli exists and is clonable. In summary, functions of the ms${2{\cdot}6}\atop1$A37 modification in E. coli cells are further delineated in this work in perspectives of adaptation to stressful environmental conditions and protein:tRNA interaction. (Abstract shortened by UMI.) ^

Turbulent flux and meteorological measurements during ARCTEX-2006 campaign

This dataset present result from the DFG- funded Arctic-Turbulence-Experiment (ARCTEX-2006) performed by the University of Bayreuth on the island of Svalbard, Norway, during the winter/spring transition 2006. From May 5 to May 19, 2006 turbulent flux and meteorological measurements were performed on the monitoring field near Ny-Ålesund, at 78°55'24'' N, 11°55'15'' E Kongsfjord, Svalbard (Spitsbergen), Norway. The ARCTEX-2006 campaign site was located about 200 m southeast of the settlement on flat snow covered tundra, 11 m to 14 m above sea level. The permanent sites used for this study consisted of the 10 m meteorological tower of the Alfred Wegener Institute for Polar- and Marine Research (AWI), the international standardized radiation measurement site of the Baseline Surface Radiation Network (BSRN), the radiosonde launch site and the AWI tethered balloon launch sites. The temporary sites - set up by the University of Bayreuth - were a 6 m meteorological gradient tower, an eddy-flux measurement complex (EF), and a laser-scintillometer section (SLS). A quality assessment and data correction was applied to detect and eliminate specific measurement errors common at a high arctic landscape. In addition, the quality checked sensible heat flux measurements are compared with bulk aerodynamic formulas that are widely used in atmosphere-ocean/land-ice models for polar regions as described in Ebert and Curry (1993, doi:10.1029/93JC00656) and Launiainen and Cheng (1995). These parameterization approaches easily allow estimation of the turbulent surface fluxes from routine meteorological measurements. The data show: - the role of the intermittency of the turbulent atmospheric fluctuation of momentum and scalars, - the existence of a disturbed vertical temperature profile (sharp inversion layer) close to the surface, - the relevance of possible free convection events for the snow or ice melt in the Arctic spring at Svalbard, and - the relevance of meso-scale atmospheric circulation pattern and air-mass advection for the near-surface turbulent heat exchange in the Arctic spring at Svalbard. Recommendations and improvements regarding the interpretation of eddy-flux and laser-scintillometer data as well as the arrangement of the instrumentation under polar distinct exchange conditions and (extreme) weather situations could be derived.

Eruption of the Håkon Mosby mud volcano recorded by the long-term observatory on mud-volcano eruptions (LOOME) between 2009 and 2010

Submarine mud volcanoes are considered an important source of methane to the water column. However, the temporal variability of their fluid transport including mud and methane emissions is largely unknown. Assuming that this transport was continuous and at steady state, methane emissions were previously proposed to result from a dynamic equilibrium between upward migration and consumption at the seabed by methane-consuming microbes. Here we have investigated non-steady state situations of vigorous mud movements and their reflection in fluid flow, seabed temperature and bathymetry. Time series of pressure, temperature, pH and seafloor photography were collected by a benthic observatory (LOOME) for 431 days at the active Håkon Mosby mud volcano. These new data document eruptions, which were accompanied by pulses of hot subsurface fluids and triggered rapid sediment uplift and lateral movement, as well as emissions of free gas.

Solid phase analysis of sediment core GeoB9309-1 from the Zambezi (or Mozambique) deep-sea fan

The Zambezi deep-sea fan, the largest of its kind along the east African continental margin, is poorly studied to date, despite its potential to record marine and terrestrial climate signals in the southwest Indian Ocean. Therefore, gravity core GeoB 9309-1, retrieved from 1219 m water depth, was investigated for various geophysical (magnetic susceptibility, porosity, colour reflectance) and geochemical (pore water and sediment geochemistry, Fe and P speciation) properties. Onboard and onshore data documented a sulphate/methane transition (SMT) zone at ~ 450-530 cm sediment depth, where the simultaneous consumption of pore water sulphate and methane liberates hydrogen sulphide and bi-carbonate into the pore space. This leads to characteristic changes in the sediment and pore water chemistry, as the reduction of primary Fe (oxyhydr)oxides, the precipitation of Fe sulphides, and the mobilization of Fe (oxyhydr)oxide-bound P. These chemical processes also lead to a marked decrease in magnetic susceptibility. Below the SMT, we find a reduction of porosity, possibly due to pore space cementation by authigenic minerals. Formation of the observed geochemical, magnetic and mineralogical patterns requires a fixation of the SMT at this distinct sediment depth for a considerable time-which we calculated to be ~ 10 000 years assuming steady-state conditions-following a period of rapid upward migration towards this interval. We postulate that the worldwide sea-level rise at the last glacial/interglacial transition (~ 10 000 years B.P.) most probably caused the fixation of the SMT at its present position, through drastically reduced sediment delivery to the deep-sea fan. In addition, we report an internal redistribution of P occurring around the SMT, closely linked to the (de)coupling of sedimentary Fe and P, and leaving a characteristic pattern in the solid P record. By phosphate re-adsorption onto Fe (oxyhydr)oxides above, and formation of authigenic P minerals (e.g. vivianite) below the SMT, deep-sea fan deposits may potentially act as long-term sinks for P.

Temperature- and strain-rate-dependent microfabric evolution in monomineralic mylonite: evidence from in situ deformation of norcamphor

Norcamphor (C7H10O) was subjected to plane strain simple shear in a see-through deformation rig at four different strain rate and temperature conditions. Two transient stages in the microfabric evolution to steady state are distinguished. The grain scale mechanisms associated with the microstructural and textural evolution vary with the applied temperature, strain rate and strain. In high-temperature-low-strain-rate experiments, computer integrated polarization microscopy reveals that the texture evolution is closely related to the crystallographic rotation paths and rotation rates of individual grains. High c-axis rotation rates at low to intermediate shear strains are related to the development of a symmetrical c-axis cross girdle by the end of the first transient stage (γ = 1.5 to 2). During the second transient stage (γ = 1.5 to 6), the cross girdle yields to an oblique c-axis single girdle as c-axis rotation rates decrease and the relative activity of grain boundary migration recrystallization increases. Steady state (γ > 8) is characterized by a stable end orientation of the sample texture and the cyclic growth, rotation and consumption of individual grains within the aggregate.

The evolution of high-temperature mylonitic microfabrics: evidence from simple shearing of a quartz analogue (norcamphor)

Plane strain simple shearing of norcamphor (C7H10O) in a see-through deformation rig to a shear strain of γ = 10.5 at a homologous temperature of Th = 0.81 yields a microfabric similar to that of quartz in amphibolite facies mylonite. Synkinematic analysis of the norcamphor microfabric reveals that the development of a steady-state texture is linked to changes in the relative activities of several grain-scale mechanisms. Three stages of textural and microstructural evolution are distinguished: (1) rotation and shearing of the intracrystalline glide planes are accommodated by localized deformation along three sets of anastomozing microshears. A symmetrical c-axis girdle reflects localized pure shear extension along the main microshear set (Sa) oblique to the bulk shear zone boundary (abbreviated as SZB); (2) progressive rotation of the microshears into parallelism with the SZB increases the component of simple shear on the Sa microshears. Grain-boundary migration recrystallization favours the survival of grains with slip systems oriented for easy glide. This is associated with a textural transition towards two stable c-axis point maxima whose skeletal outline is oblique with respect to the Sa microshears and the SZB; and (3) at high shear strains (γ > 8), the microstructure, texture and mechanism assemblage are strain invariant, but strain continues to partition into rotating sets of microshears. Steady state is therefore a dynamic, heterogeneous condition involving the cyclic nucleation, growth and consumption of grains.

Steps Ahead in the Few-Group CRoss-Section Library Generation at the Pin Level

There exists an interest in performing pin-by-pin calculations coupled with thermal hydraulics so as to improve the accuracy of nuclear reactor analysis. In the framework of the EU NURISP project, INRNE and UPM have generated an experimental version of a few group diffusion cross sections library with discontinuity factors intended for VVER analysis at the pin level with the COBAYA3 code. The transport code APOLLO2 was used to perform the branching calculations. As a first proof of principle the library was created for fresh fuel and covers almost the full parameter space of steady state and transient conditions. The main objective is to test the calculation schemes and post-processing procedures, including multi-pin branching calculations. Two library options are being studied: one based on linear table interpolation and another one using a functional fitting of the cross sections. The libraries generated with APOLLO2 have been tested with the pin-by-pin diffusion model in COBAYA3 including discontinuity factors; first comparing 2D results against the APOLLO2 reference solutions and afterwards using the libraries to compute a 3D assembly problem coupled with a simplified thermal-hydraulic model.

Assessment of CTF boiling transition and critical heat flux modeling capabilities using the OECD/NRC BFBT and PSBT benchmark databases

The need to refine models for best-estimate calculations, based on good-quality experimental data, has been expressed in many recent meetings in the field of nuclear applications. The modeling needs arising in this respect should not be limited to the currently available macroscopic methods but should be extended to next-generation analysis techniques that focus on more microscopic processes. One of the most valuable databases identified for the thermalhydraulics modeling was developed by the Nuclear Power Engineering Corporation (NUPEC), Japan. From 1987 to 1995, NUPEC performed steady-state and transient critical power and departure from nucleate boiling (DNB) test series based on the equivalent full-size mock-ups. Considering the reliability not only of the measured data, but also other relevant parameters such as the system pressure, inlet sub-cooling and rod surface temperature, these test series supplied the first substantial database for the development of truly mechanistic and consistent models for boiling transition and critical heat flux. Over the last few years the Pennsylvania State University (PSU) under the sponsorship of the U.S. Nuclear Regulatory Commission (NRC) has prepared, organized, conducted and summarized the OECD/NRC Full-size Fine-mesh Bundle Tests (BFBT) Benchmark. The international benchmark activities have been conducted in cooperation with the Nuclear Energy Agency/Organization for Economic Co-operation and Development (NEA/OECD) and Japan Nuclear Energy Safety (JNES) organization, Japan. Consequently, the JNES has made available the Boiling Water Reactor (BWR) NUPEC database for the purposes of the benchmark. Based on the success of the OECD/NRC BFBT benchmark the JNES has decided to release also the data based on the NUPEC Pressurized Water Reactor (PWR) subchannel and bundle tests for another follow-up international benchmark entitled OECD/NRC PWR Subchannel and Bundle Tests (PSBT) benchmark. This paper presents an application of the joint Penn State University/Technical University of Madrid (UPM) version of the well-known subchannel code COBRA-TF, namely CTF, to the critical power and departure from nucleate boiling (DNB) exercises of the OECD/NRC BFBT and PSBT benchmarks

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography. The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.

Connexin-43 in the osteogenic BM niche regulates its cellular composition and the bidirectional traffic of hematopoietic stem cells and progenitors

Connexin-43 (Cx43), a gap junction protein involved in control of cell proliferation, differentiation and migration, has been suggested to have a role in hematopoiesis. Cx43 is highly expressed in osteoblasts and osteogenic progenitors (OB/P). To elucidate the biologic function of Cx43 in the hematopoietic microenvironment (HM) and its influence in hematopoietic stem cell (HSC) activity, we studied the hematopoietic function in an in vivo model of constitutive deficiency of Cx43 in OB/P. The deficiency of Cx43 in OB/P cells does not impair the steady state hematopoiesis, but disrupts the directional trafficking of HSC/progenitors (Ps) between the bone marrow (BM) and peripheral blood (PB). OB/P Cx43 is a crucial positive regulator of transstromal migration and homing of both HSCs and progenitors in an irradiated microenvironment. However, OB/P Cx43 deficiency in nonmyeloablated animals does not result in a homing defect but induces increased endosteal lodging and decreased mobilization of HSC/Ps associated with proliferation and expansion of Cxcl12-secreting mesenchymal/osteolineage cells in the BM HM in vivo. Cx43 controls the cellular content of the BM osteogenic microenvironment and is required for homing of HSC/Ps in myeloablated animals

Assessment of CTF boiling transition and critical heat flux modeling capabilities using the OECD/NRC BFBT and PSBT benchmark databases

Over the last few years, the Pennsylvania State University (PSU) under the sponsorship of the US Nuclear Regulatory Commission (NRC) has prepared, organized, conducted, and summarized two international benchmarks based on the NUPEC data—the OECD/NRC Full-Size Fine-Mesh Bundle Test (BFBT) Benchmark and the OECD/NRC PWR Sub-Channel and Bundle Test (PSBT) Benchmark. The benchmarks’ activities have been conducted in cooperation with the Nuclear Energy Agency/Organization for Economic Co-operation and Development (NEA/OECD) and the Japan Nuclear Energy Safety (JNES) Organization. This paper presents an application of the joint Penn State University/Technical University of Madrid (UPM) version of the well-known sub-channel code COBRA-TF (Coolant Boiling in Rod Array-Two Fluid), namely, CTF, to the steady state critical power and departure from nucleate boiling (DNB) exercises of the OECD/NRC BFBT and PSBT benchmarks. The goal is two-fold: firstly, to assess these models and to examine their strengths and weaknesses; and secondly, to identify the areas for improvement.

Space-state robust control of a Buck converter with amorphous core coil and variable load

Pulse-width modulation is widely used to control electronic converters. One of the most frequently used topologies for high DC voltage/low DC voltage conversion is the Buck converter. These converters are described by a second order system with an LC filter between the switching subsystem and the load. The use of a coil with an amorphous magnetic material core rather than an air core permits the design of smaller converters. If high switching frequencies are used to obtain high quality voltage output, then the value of the auto inductance L is reduced over time. Robust controllers are thus needed if the accuracy of the converter response must be preserved under auto inductance and payload variations. This paper presents a robust controller for a Buck converter based on a state space feedback control system combined with an additional virtual space variable which minimizes the effects of the inductance and load variations when a switching frequency that is not too high is applied. The system exhibits a null steady-state average error response for the entire range of parameter variations. Simulation results and a comparison with a standard PID controller are also presented.

Desarrollo de métodos numérico-interpretativos para la realización de ensayos de bombeo

Los ensayos de bombeo son, sin lugar a dudas, una de las pruebas más fiables y de mayor interés que se hacen en el medio físico. No son pruebas estrictamente puntuales, dado que el bombeo atrae flujo desde distancias lejanas al pozo, la prueba tiene una excelente representatividad espacial. Los métodos de interpretación mediante ensayos de bombeo se empezaron a plantear en la primera mitad del pasado siglo. Con los ensayos de bombeo se puede calcular la transmisividad y coeficiente de almacenamiento de las formaciones acuíferas y suministran información sobre el tipo de acuífero, la calidad constructiva del pozo de extracción, la existencia de barreras impermeable o bordes de recarga próximos, e incluso en algunas circunstancias permiten el cálculo del área de embalse subterráneo. Desde mediados del siglo 20 existe una eficaz y abundante gama de métodos analítico-interpretativos de ensayos de bombeo, tanto en régimen permanente como transitorio. Estos métodos son ampliamente conocidos y están muy experimentados a lo largo de muchos países, sin embargo, hoy día, podrían utilizarse modelos de flujo para la interpretación, logrando la misma fiabilidad e incluso mejores posibilidades de análisis. Muchos ensayos que no pueden interpretarse porque las configuraciones del medio son demasiado complejas y no están disponibles, o no es posible, el desarrollo de métodos analíticos, tienen buena adaptación y en ocasiones muy fácil solución haciendo uso de los métodos numéricos de simulación del flujo. En esta tesis se ha buscado una vía de interpretar ensayos de bombeo haciendo uso de modelos de simulación del flujo. Se utiliza el modelo universal MODFLOW del United States Geological Survey, en el cual se configura una celda de simulación y mallado particularmente adecuados para el problema a tratar, se valida con los métodos analíticos existentes. Con la célula convenientemente validada se simulan otros casos en los que no existen métodos analíticos desarrollados dada la complejidad del medio físico a tratar y se sacan las oportunas conclusiones. Por último se desarrolla un modelo específico y la correspondiente aplicación de uso general para la interpretación numérica de ensayos de bombeo tanto con las configuraciones normales como con configuraciones complejas del medio físico. ABSTRACT Pumping tests are, without doubt, one of the most reliable and most interesting tests done in the physical environment. They are not strictly anecdotal evidence, since pumping flow attracts from far distances to the well, the test has excellent spatial representation. Methods of interpretation by pumping tests began to arise in the first half of last century. With pumping tests, can be calculated transmissivity and storage coefficient of the aquifer formations, and provide information on the type of aquifer, the construction quality of the well, the existence of waterproof barriers or borders next recharge, and even in some circumstances allow calculating the area of underground reservoir. Since the mid-20th century there is effective and abundant range of analytical interpretative pumping tests, both in steady state and transient methods. These methods are very widely known and experienced over many countries, however, nowadays, may flow models used for interpretation, obtaining equally reliable or even better possibilities for analysis. Many trials cannot be interpreted as environmental settings are too complex and are not available, or not possible, the development of analytical methods, have good adaptation and sometimes very easily solved using numerical flow simulation methods. This thesis has sought a way to interpret pumping tests using flow simulation models. MODFLOW universal model of United States Geological Survey, in which a simulation cell and meshing particularly suitable for the problem to be treated, is validated with existing analytical methods used is set. With suitably validated cell other cases where there are no analytical methods developed given the complexity of the physical environment to try and draw appropriate conclusions are simulated. Finally, a specific model and the corresponding application commonly used for numerical interpretation of pumping tests both with normal settings as complex configurations of the physical environment is developed.