Discrete element methods for asphalt concrete : development and application of user-defined microstructural models and a viscoelastic micromechanical model

As an important Civil Engineering material, asphalt concrete (AC) is commonly used to build road surfaces, airports, and parking lots. With traditional laboratory tests and theoretical equations, it is a challenge to fully understand such a random composite material. Based on the discrete element method (DEM), this research seeks to develop and implement computer models as research approaches for improving understandings of AC microstructure-based mechanics. In this research, three categories of approaches were developed or employed to simulate microstructures of AC materials, namely the randomly-generated models, the idealized models, and image-based models. The image-based models were recommended for accurately predicting AC performance, while the other models were recommended as research tools to obtain deep insight into the AC microstructure-based mechanics. A viscoelastic micromechanical model was developed to capture viscoelastic interactions within the AC microstructure. Four types of constitutive models were built to address the four categories of interactions within an AC specimen. Each of the constitutive models consists of three parts which represent three different interaction behaviors: a stiffness model (force-displace relation), a bonding model (shear and tensile strengths), and a slip model (frictional property). Three techniques were developed to reduce the computational time for AC viscoelastic simulations. It was found that the computational time was significantly reduced to days or hours from years or months for typical three-dimensional models. Dynamic modulus and creep stiffness tests were simulated and methodologies were developed to determine the viscoelastic parameters. It was found that the DE models could successfully predict dynamic modulus, phase angles, and creep stiffness in a wide range of frequencies, temperatures, and time spans. Mineral aggregate morphology characteristics (sphericity, orientation, and angularity) were studied to investigate their impacts on AC creep stiffness. It was found that aggregate characteristics significantly impact creep stiffness. Pavement responses and pavement-vehicle interactions were investigated by simulating pavement sections under a rolling wheel. It was found that wheel acceleration, steadily moving, and deceleration significantly impact contact forces. Additionally, summary and recommendations were provided in the last chapter and part of computer programming codes wree provided in the appendixes.

Four-window technique for measuring optical-phase-space-time-frequency tomography

A new approach, the four-window technique, was developed to measure optical phase-space-time-frequency tomography (OPSTFT). The four-window technique is based on balanced heterodyne detection with two local oscillator (LO) fields. This technique can provide independent control of position, momentum, time and frequency resolution. The OPSTFT is a Wigner distribution function of two independent Fourier transform pairs, phase-space and time-frequency. The OPSTFT can be applied for early disease detection.

Pupil phase apodization for achromatic imaging of extra-solar planets

Direct imaging of extra-solar planets in the visible and infrared region has generated great interest among scientists and the general public as well. However, this is a challenging problem. Diffculties of detecting a planet (faint source) are caused, mostly, by two factors: sidelobes caused by starlight diffraction from the edge of the pupil and the randomly scattered starlight caused by the phase errors from the imperfections in the optical system. While the latter diffculty can be corrected by high density active deformable mirrors with advanced phase sensing and control technology, the optimized strategy for suppressing the diffraction sidelobes is still an open question. In this thesis, I present a new approach to the sidelobe reduction problem: pupil phase apodization. It is based on a discovery that an anti-symmetric spatial phase modulation pattern imposed over a pupil or a relay plane causes diffracted starlight suppression sufficient for imaging of extra-solar planets. Numerical simulations with specific square pupil (side D) phase functions, such as ... demonstrate annulling in at least one quadrant of the diffraction plane to the contrast level of better than 10^12 with an inner working angle down to 3.5L/D (with a = 3 and e = 10^3). Furthermore, our computer experiments show that phase apodization remains effective throughout a broad spectrum (60% of the central wavelength) covering the entire visible light range. In addition to the specific phase functions that can yield deep sidelobe reduction on one quadrant, we also found that a modified Gerchberg-Saxton algorithm can help to find small sized (101 x 101 element) discrete phase functions if regional sidelobe reduction is desired. Our simulation shows that a 101x101 segmented but gapless active mirror can also generate a dark region with Inner Working Distance about 2.8L/D in one quadrant. Phase-only modulation has the additional appeal of potential implementation via active segmented or deformable mirrors, thereby combining compensation of random phase aberrations and diffraction halo removal in a single optical element.

TIME LAPSE SEISMIC OBSERVATIONS AND EFFECTS OF RESERVOIR COMPRESSIBILITY AT TEAL SOUTH OIL FIELD

One of the original ocean-bottom time-lapse seismic studies was performed at the Teal South oil field in the Gulf of Mexico during the late 1990’s. This work reexamines some aspects of previous work using modern analysis techniques to provide improved quantitative interpretations. Using three-dimensional volume visualization of legacy data and the two phases of post-production time-lapse data, I provide additional insight into the fluid migration pathways and the pressure communication between different reservoirs, separated by faults. This work supports a conclusion from previous studies that production from one reservoir caused regional pressure decline that in turn resulted in liberation of gas from multiple surrounding unproduced reservoirs. I also provide an explanation for unusual time-lapse changes in amplitude-versus-offset (AVO) data related to the compaction of the producing reservoir which, in turn, changed an isotropic medium to an anisotropic medium. In the first part of this work, I examine regional changes in seismic response due to the production of oil and gas from one reservoir. The previous studies primarily used two post-production ocean-bottom surveys (Phase I and Phase II), and not the legacy streamer data, due to the unavailability of legacy prestack data and very different acquisition parameters. In order to incorporate the legacy data in the present study, all three poststack data sets were cross-equalized and examined using instantaneous amplitude and energy volumes. This approach appears quite effective and helps to suppress changes unrelated to production while emphasizing those large-amplitude changes that are related to production in this noisy (by current standards) suite of data. I examine the multiple data sets first by using the instantaneous amplitude and energy attributes, and then also examine specific apparent time-lapse changes through direct comparisons of seismic traces. In so doing, I identify time-delays that, when corrected for, indicate water encroachment at the base of the producing reservoir. I also identify specific sites of leakage from various unproduced reservoirs, the result of regional pressure blowdown as explained in previous studies; those earlier studies, however, were unable to identify direct evidence of fluid movement. Of particular interest is the identification of one site where oil apparently leaked from one reservoir into a “new” reservoir that did not originally contain oil, but was ideally suited as a trap for fluids leaking from the neighboring spill-point. With continued pressure drop, oil in the new reservoir increased as more oil entered into the reservoir and expanded, liberating gas from solution. Because of the limited volume available for oil and gas in that temporary trap, oil and gas also escaped from it into the surrounding formation. I also note that some of the reservoirs demonstrate time-lapse changes only in the “gas cap” and not in the oil zone, even though gas must be coming out of solution everywhere in the reservoir. This is explained by interplay between pore-fluid modulus reduction by gas saturation decrease and dry-frame modulus increase by frame stiffening. In the second part of this work, I examine various rock-physics models in an attempt to quantitatively account for frame-stiffening that results from reduced pore-fluid pressure in the producing reservoir, searching for a model that would predict the unusual AVO features observed in the time-lapse prestack and stacked data at Teal South. While several rock-physics models are successful at predicting the time-lapse response for initial production, most fail to match the observations for continued production between Phase I and Phase II. Because the reservoir was initially overpressured and unconsolidated, reservoir compaction was likely significant, and is probably accomplished largely by uniaxial strain in the vertical direction; this implies that an anisotropic model may be required. Using Walton’s model for anisotropic unconsolidated sand, I successfully model the time-lapse changes for all phases of production. This observation may be of interest for application to other unconsolidated overpressured reservoirs under production.

DEVELOPMENT OF CONFOCAL IMAGING TECHNIQUES FOR PROBING INTERFACIAL DYNAMICS IN MICROSCALE, GAS-LIQUID, TWO-PHASE FLOW

Micro-scale, two-phase flow is found in a variety of devices such as Lab-on-a-chip, bio-chips, micro-heat exchangers, and fuel cells. Knowledge of the fluid behavior near the dynamic gas-liquid interface is required for developing accurate predictive models. Light is distorted near a curved gas-liquid interface preventing accurate measurement of interfacial shape and internal liquid velocities. This research focused on the development of experimental methods designed to isolate and probe dynamic liquid films and measure velocity fields near a moving gas-liquid interface. A high-speed, reflectance, swept-field confocal (RSFC) imaging system was developed for imaging near curved surfaces. Experimental studies of dynamic gas-liquid interface of micro-scale, two-phase flow were conducted in three phases. Dynamic liquid film thicknesses of segmented, two-phase flow were measured using the RSFC and compared to a classic film thickness deposition model. Flow fields near a steadily moving meniscus were measured using RSFC and particle tracking velocimetry. The RSFC provided high speed imaging near the menisci without distortion caused the gas-liquid interface. Finally, interfacial morphology for internal two-phase flow and droplet evaporation were measured using interferograms produced by the RSFC imaging technique. Each technique can be used independently or simultaneously when.

Ovarian stimulation to cryopreserve fertilized oocytes in cancer patients can be started in the luteal phase

OBJECTIVE: To analyze if oocytes can be obtained in all patients before cancer treatment within 2 weeks by initiating ovarian stimulation during the follicular or luteal phase. DESIGN: Prospective controlled multicenter trial. SETTING: Four university-based centers. PATIENT(S): Forty cancer patients before chemotherapy. INTERVENTION(S): Twenty-eight patients were stimulated with gonadotropins in the follicular phase (group I). In 12 patients (group II), ovarian stimulation was initiated in the luteal phase, and these received GnRH antagonists and recombinant FSH. In 14 patients, 143 oocytes were further processed for fertilization by intracytoplasmic sperm injection (ICSI). MAIN OUTCOME MEASURE(S): Number of oocytes aspirated after ovarian stimulation, cumulative FSH/hMG dosage, viability and maturity of oocytes, and fertilization rate by ICSI. RESULT(S): Patients in group I (age 27.6 +/- 4.9 yrs) were stimulated on average for 10.6 days, and patients in group II (age 31.2 +/- 5.7 yrs) for 11.4 days. Total amount of FSH was on average 2,255 IU (I) and 2,720 IU (II) per patient. Average and median numbers of aspirated oocytes were, respectively, 13.1 and 11.5 (I) versus 10.0 and 8.5 (II); 83.7% (I) and 80.4% (II) of the oocytes were mature and viable and could be treated by ICSI. Fertilization rate was 61.0% (I) versus 75.6% (II). CONCLUSION(S): This pilot study suggests that oocytes can be obtained before cancer treatment efficiently irrespective of the phase of the menstrual cycle.

Gingival labial recessions in orthodontically treated and untreated individuals: a case - control study

OBJECTIVES To evaluate the long-term development of labial gingival recessions during orthodontic treatment and retention phase. MATERIAL AND METHODS In this retrospective case-control study, the presence of gingival recession was scored (Yes or No) on plaster models of 100 orthodontic patients (cases) and 120 controls at the age of 12 (T12 ), 15 (T15 ), 18 (T18 ), and 21 (T21 ) years. In the treated group, T12 reflected the start of orthodontic treatment and T15 - the end of active treatment and the start of retention phase with bonded retainers. Independent t-tests, Fisher's exact tests and a fitted two-part "hurdle" model were used to identify the effect of orthodontic treatment/retention on recessions. RESULTS The proportion of subjects with recessions was consistently higher in cases than controls. Overall, the odds ratio for orthodontic patients as compared with controls to have recessions is 4.48 (p < 0.001; 95% CI: 2.61-7.70). CONCLUSIONS Within the limits of the present research design, orthodontic treatment and/or the retention phase may be risk factors for the development of labial gingival recessions. In orthodontically treated subjects, mandibular incisors seem to be the most vulnerable to the development of gingival recessions.

2-methiopropamine, a thiopene analogue of metamphetamine: studies on its metabolism abt detectability in the rat and human using

2-Methiopropamine [1-(thiophen-2-yl)-2-methylaminopropane, 2-MPA], a thiophene analogue of methamphetamine, is available from online vendors selling "Research chemicals." The first samples were seized by the German police in 2011. As it is a recreational stimulant, its inclusion in routine drug screening protocols should be required. The aims of this study were to identify the phase I and II metabolites of 2-MPA in rat and human urine and to identify the human cytochrome-P450 (CYP) isoenzymes involved in its phase I metabolism. In addition, the detectability of 2-MPA in urine samples using the authors' well-established gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-linear ion trap-mass spectrometry (LC-MS(n)) screening protocols was also evaluated. The metabolites were isolated from rat and human urine samples by solid-Phase extraction without or following enzymatic cleavage of conjugates. The phase I metabolites, following acetylation, were separated and identified by GC-MS and/or liquid chromatography-high-resolution linear ion trap mass spectrometry (LC-HR-MS(n)) and the phase II metabolites by LC-HR-MS(n). The following Major metabolic pathways were proposed: N-demethylation, hydroxylation at the side chain and at the thiophene ring, and combination of these transformations followed by glucuronidation and/or sulfation. CYP1A2, CYP2C19, CYP2D6, and CYP3A4 were identified as the major phase I metabolizing enzymes. They were also involved in the N-demethylation of the analogue methamphetamine and CYP2C19, CYP2D6, and CYP3A4 in its ring hydroxylation. Following the administration of a typical user's dose, 2-MPA and its metabolites were identified in rat urine using the authors' GC-MS and the LC-MS(n) screening approaches. Ingestion of 2-MPA could also be detected by both protocols in an authentic human urine sample.

Xenobiotic Metabolism Genes and Clubfoot

Idiopathic or isolated clubfoot is a common orthopedic birth defect that affects approximately 135,000 children worldwide. It is characterized by equinus, varus and adductus deformities of the ankle and foot. Correction of clubfoot involves months of serial manipulations, castings and bracing, with surgical correction needed in forty percent of cases. Multifactorial etiology has been suggested in numerous studies with both environmental and genetic factors playing an etiologic role. Maternal smoking during pregnancy is the only common environmental factor that has consistently been shown to increase the risk for clubfoot. Moreover, a positive family history of clubfoot and maternal smoking increases the risk of clubfoot twenty fold. These findings suggest that genetic variation in smoking metabolism genes may increase susceptibility to clubfoot. Based on this reasoning, we interrogated eight candidate genes, chosen based on their involvement in phase 1 and 2 cigarette smoke metabolism. Twenty-two SNPs and two null alleles in eight genes (CYP1A1, CYP1A2, CYP1B1, CYP2A6, EPHX1, NAT2, GSTM1 and GSTT1) were genotyped in a dataset composed of nonHispanic white and Hispanic multiplex and simplex families. Only one SNP in CYP1A1, rs1048943, had significantly altered transmission in the aggregate and multiplex NHW datasets (p=0.003 and p=0.009). Perturbation of CYP1A1 by rs1048943 polymorphism causes an increase in the amount of harmful, adduct forming metabolic intermediates. A significant gene interaction between EPHX1 and NAT2 was also found (p=0.007). This interaction may affect the metabolism of harmful metabolic intermediates. Additionally, marginal interactions were found for other xenobiotic genes and these interactions may play a contributory role in clubfoot. Importantly, for CYP1A2, significant maternal (p=0.03; RR=1.24; 95% CI: 1.04-1.44) and fetal (p=0.01; RR=1.33; 95% CI: 1.13-1.54) genotypic effects were identified suggesting that both maternal and fetal genotypes impact normal limb development. No association was found for maternal smoking status and tobacco metabolism genes. Together, these results suggest that xenobiotic metabolism genes may play a contributory role in the etiology of clubfoot regardless of maternal smoking status and may impact foot development through perturbation of tobacco metabolic pathways.

MODELING REACTIVE GAS UPTAKE, TRANSPORT, AND TRANSFORMATION IN AGGREGATED SOILS

Gas diffusion research in soils covers, to a large extent, the transport behavior of practically insoluble gases. We extend the mathematical description of gas transport to include reactive gaseous components that hydrolyze in water such as SO2 and CO2. The path between the free atmosphere and the microporous niches is modeled by assuming penetration through gas-filled macropores, air-water phase transfer, and diffusion and speciation in the liquid phase. For hydrolyzable gases, the rate of mass transfer into and the total absorption capacity of the soil solution may be high. Both the capacity and the transfer rate are influenced by the soil-solution pH; for high pH, they become extremely high for SO2. The soil absorption of such gases is also influenced by soil structure. Well-aerated, near-neutral soils are a potentially important sink for SO2.

Historical and idealized climate model experiments: an intercomparison of Earth system models of intermediate complexity

Both historical and idealized climate model experiments are performed with a variety of Earth system models of intermediate complexity (EMICs) as part of a community contribution to the Intergovernmental Panel on Climate Change Fifth Assessment Report. Historical simulations start at 850 CE and continue through to 2005. The standard simulations include changes in forcing from solar luminosity, Earth's orbital configuration, CO2, additional greenhouse gases, land use, and sulphate and volcanic aerosols. In spite of very different modelled pre-industrial global surface air temperatures, overall 20th century trends in surface air temperature and carbon uptake are reasonably well simulated when compared to observed trends. Land carbon fluxes show much more variation between models than ocean carbon fluxes, and recent land fluxes appear to be slightly underestimated. It is possible that recent modelled climate trends or climate–carbon feedbacks are overestimated resulting in too much land carbon loss or that carbon uptake due to CO2 and/or nitrogen fertilization is underestimated. Several one thousand year long, idealized, 2 × and 4 × CO2 experiments are used to quantify standard model characteristics, including transient and equilibrium climate sensitivities, and climate–carbon feedbacks. The values from EMICs generally fall within the range given by general circulation models. Seven additional historical simulations, each including a single specified forcing, are used to assess the contributions of different climate forcings to the overall climate and carbon cycle response. The response of surface air temperature is the linear sum of the individual forcings, while the carbon cycle response shows a non-linear interaction between land-use change and CO2 forcings for some models. Finally, the preindustrial portions of the last millennium simulations are used to assess historical model carbon-climate feedbacks. Given the specified forcing, there is a tendency for the EMICs to underestimate the drop in surface air temperature and CO2 between the Medieval Climate Anomaly and the Little Ice Age estimated from palaeoclimate reconstructions. This in turn could be a result of unforced variability within the climate system, uncertainty in the reconstructions of temperature and CO2, errors in the reconstructions of forcing used to drive the models, or the incomplete representation of certain processes within the models. Given the forcing datasets used in this study, the models calculate significant land-use emissions over the pre-industrial period. This implies that land-use emissions might need to be taken into account, when making estimates of climate–carbon feedbacks from palaeoclimate reconstructions.

The functional and structural interactions between v-{\it mos\/} proteins and cell cycle control elements

The v-mos gene of Moloney murine sarcoma virus (Mo-MuSv) encodes a serine/threonine protein kinase capable of inducing cellular transformation. The c-mos protein is an important cell cycle regulator that functions during meiotic cell division cycles in germ cells. The overall function of c-mos in controlling meiosis is becoming better understood but the role of v-mos in malignant transformation of cells is largely unknown.^ In this study, v-mos protein was shown to be phosphorylated by M phase kinase in vitro and in vivo. The kinase activity and neoplastic transforming ability of v-mos is positively regulated by the phosphorylation. Together with the earlier finding of activation of M phase kinase by c-mos, these results raise the possibility of mutual regulation between M phase kinase and mos kinases.^ In addition to its functional interaction with the M phase kinase, the v-mos protein was shown to be present in the same protein complex with a cyclin-dependent kinase (cdk). In addition, an antibody that recognizes the cdk proteins was shown to co-precipitate the v-mos proteins in the interphase and mitotic cells transformed by p85$\sp{\rm gag-mos}$. Cdk proteins have been shown to be associated with nonmitotic cyclins which are potential oncogenes. The perturbation of cdk kinase or the activation of non-mitotic cyclins as oncogenes by v-mos could contribute directly to v-mos induced cellular transformation. v-mos proteins were also shown to interact with tubulin and vimentin, the essential components of microtubules and type IV intermediate filaments, respectively. The organizations of both microtubules and intermediate filaments are cell cycle-regulated. These results suggest that the v-mos kinase could be directly involved in inducing morphological changes typically seen in transformed cells.^ The interactions between the v-mos protein and these cell cycle control elements in regards to v-mos induced neoplastic transformation are discussed in detail in the text. ^

Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

Changes in cerebral hemodynamics and oxygenation induced by different speech tasks – an assessment by combining functional near‐infrared spectroscopy and capnography

Introduction In several studies, we found that during guided rhythmic speech exercises, a decrease in cerebral hemodynamics and oxygenation occurred as the result of a decrease in the partial pressure of carbon dioxide in the arterial blood (PaCO2) during speaking. To further explore the effect of PaCO2 variations on cerebral hemodynamics and oxygenation, the aim of the present study was to investigate the impact of spoken, inner and heard speech tasks on these parameters. Material and Methods Speech tasks included recitation or inner recitation or listening to hexameter, alliteration, prose, or performing mental arithmetic. The following physiological parameters were measured: tissue oxygen saturation (StO2) and absolute concentrations of oxyhemoglobin, deoxyhemoglobin, total hemoglobin (over the left and right anterior prefrontal cortex, using an ISS OxiplexTS frequency domain near-infrared spectrometer) and end-tidal CO2 (PETCO2; using Nellcor N1000 and Datex NORMOCAP capnographs). Statistical analysis was applied to the differences between baseline, 2 tasks, and 3 post-baseline periods. Data of 3 studies with 24, 7 and 29 healthy subjects, respectively, were combined, and linear regression analyses were calculated. Results Linear regression analyses revealed significant relations between changes in oxyhemoglobin, deoxyhemoglobin, total hemoglobin or StO2 and the participants’ age, the baseline PETCO2 or certain speech tasks. While hexameter verses affected changes during the tasks, alliteration verses only affected changes during the recovery phase. Discussion and Conclusion The observed effects in hemodynamics and oxygenation indicate a combination of neurovascular coupling (increased neuronal activity leading to an increase in the cerebral metabolic rate of oxygen resulting in an increase in cerebral flood flow/volume) and CO2 reactivity (increased breathing during speech tasks causing a decrease in PaCO2 leading to vasoconstriction and decrease in cerebral blood flow). The neurovascular coupling characteristics are task-dependent. References Scholkmann F, Gerber U, Wolf M, Wolf U. End-tidal CO2: An important parameter for a correct interpretation in functional brain studies using speech tasks. Neuroimage 2013;66:71-79. Scholkmann F, Wolf M, Wolf U. The effect of inner speech on arterial CO2, cerebral hemodynamics and oxygenation – A functional NIRS study. Adv Exp Med Biol 2013;789:81-87.

Metabolomics and its potential in drug development

Metabolomics is the global and unbiased survey of the complement of small molecules (say, <1 kDa) in a biofluid, tissue, organ or organism and measures the end-products of the cellular metabolism of both endogenous and exogenous substrates. Many drug candidates fail during Phase II and III clinical trials at an enormous cost to the pharmaceutical industry in terms of both time lost and of financial resources. The constantly evolving model of drug development now dictates that biomarkers should be employed in preclinical development for the early detection of likely-to-fail candidates. Biomarkers may also be useful in the preselection of patients and through the subclassification of diseases in clinical drug development. Here we show with examples how metabolomics can assist in the preclinical development phases of discovery, pharmacology, toxicology, and ADME. Although not yet established as a clinical trial patient prescreening procedure, metabolomics shows considerable promise in this regard. We can be certain that metabolomics will join genomics and transcriptomics in lubricating the wheels of clinical drug development in the near future.