THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Molecular evolution of primate immunodeficiency viruses and hepatitis delta virus

Primate immunodeficiency viruses, or lentiviruses (HIV-1, HIV-2, and SIV), and hepatitis delta virus (HDV) are RNA viruses characterized by rapid evolution. Infection by primate immunodeficiency viruses usually results in the development of acquired immunodeficiency syndrome (AIDS) in humans and AIDS-like illnesses in Asian macaques. Similarly, hepatitis delta virus infection causes hepatitis and liver cancer in humans. These viruses are heterogeneous within an infected patient and among individuals. Substitution rates in the virus genomes are high and vary in different lineages and among sites. Methods of phylogenetic analysis were applied to study the evolution of primate lentiviruses and the hepatitis delta virus. The following results have been obtained: (1) The substitution rate varies among sites of primate lentivirus genes according to the two parameter gamma distribution, with the shape parameter $\alpha$ being close to 1. (2) Primate immunodeficiency viruses fall into species-specific lineages. Therefore, viral transmissions across primate species are not as frequent as suggested by previous authors. (3) Primate lentiviruses have acquired or lost their pathogenicity several times in the course of evolution. (4) Evidence was provided for multiple infections of a North American patient by distinct HIV-1 strains of the B subtype. (5) Computer simulations indicate that the probability of committing an error in testing HIV transmission depends on the number of virus sequences and their length, the divergence times among sequences, and the model of nucleotide substitution. (6) For future investigations of HIV-1 transmissions, using longer virus sequences and avoiding the use of distant outgroups is recommended. (7) Hepatitis delta virus strains are usually related according to the geographic region of isolation. (8) Evolution of HDV is characterized by the rate of synonymous substitution being lower than the nonsynonymous substitution rate and the rate of evolution of the noncoding region. (9) There is a strong preference for G and C nucleotides at the third codon positions of the HDV coding region. ^

STUDIES INTO THE MOLECULAR MECHANISMS REGULATING MUSCULAR ATROPHY RESULTING FROM HINDLIMB IMMOBILIZATION IN THE RAT (PROTEIN, SYNTHESIS, RNA)

The purpose of the work performed in this dissertation was to examine some of the possible regulatory mechanisms involved in the initiation of muscular atrophy during periods of decreased muscle utilization resulting from hindlimb immobilization in the rat. A 37% decrease in the rate of total muscle protein synthesis which has been observed to occur in the first 6 h of immobilization contributes significantly to the observed loss of protein during immobilization.^ The rates of cytochrome c and actin synthesis were determined in adult rat red vastus lateralis and gastrocnemius muscles, respectively, by the constant infusion and incorporation of ('3)H-tyrosine into protein. The fractional synthesis rates of both actin and cytochrome c were significantly decreased (P < 0.05) in the 6th h of hindlimb immobilization.^ RHA was extracted from adult rat gastrocnemius muscle by modification of the phenol: chloroform: SDS extraction procedures commonly used for preparation of RNA for hybridization analysis from other mammalian tissues. RNA content of rat gastrocnemius muscle, as determined by this method of extraction and its subsequent quantification by UV absorbance and orcinol assay, was significantly greater than the RNA content previously determined for adult rat gastrocnemius by other commonly employed methods.^ RNA extracted by this method from gastrocnemius muscles of control and 6h immobilized rats was subjected to "dot blot" hybridization to ('32)P-labelled probe from plasmid p749, containing a cDNA sequence complementary to (alpha)-actin mRNA and from rat skeletal muscle. (alpha)-Actin specific mRNA content as estimated by this procedure is not significantly decreased in rat gastrocnemius following 6h or hindlimb immobilization. However, (alpha)-actin specific mRNA content is significantly decreased (P < 0.05) in adult rat gastrocnemius (alpha)-actin specific mRNA is not decreased in adult rat gastrocnemius muscle following 6h of immobilization, a time when actin synthesis is significantly decreased, it is concluded that a change in (alpha)-actin specific mRNA content is not the initiating event responsible for the early decrease in actin synthesis observed in the 6th h of immobilization. ^

Delineating the molecular basis of human genetic diseases: Epigenetic and functional studies

Identifying and characterizing the genes responsible for inherited human diseases will ultimately lead to a more holistic understanding of disease pathogenesis, catalyze new diagnostic and treatment modalities, and provide insights into basic biological processes. This dissertation presents research aimed at delineating the genetic and molecular basis of human diseases through epigenetic and functional studies and can be divided into two independent areas of research. The first area of research describes the development of two high-throughput melting curve based methods to assay DNA methylation, referred to as McMSP and McCOBRA. The goal of this project was to develop DNA methylation methods that can be used to rapidly determine the DNA methylation status at a specific locus in a large number of samples. McMSP and McCOBRA provide several advantages over existing methods, as they are simple, accurate, robust, and high-throughput making them applicable to large-scale DNA methylation studies. McMSP and McCOBRA were then used in an epigenetic study of the complex disease Ankylosing spondylitis (AS). Specifically, I tested the hypothesis that aberrant patterns of DNA methylation in five AS candidate genes contribute to disease susceptibility. While no statistically significant methylation differences were observed between cases and controls, this is the first study to investigate the hypothesis that epigenetic variation contributes to AS susceptibility and therefore provides the conceptual framework for future studies. ^ In the second area of research, I performed experiments to better delimit the function of aryl hydrocarbon receptor-interacting protein-like 1 (AIPL1), which when mutated causes various forms of inherited blindness such as Leber congenital amaurosis. A yeast two-hybrid screen was performed to identify putative AIPL1-interacting proteins. After screening 2 × 106 bovine retinal cDNA library clones, 6 unique putative AIPL1-interacting proteins were identified. While these 6 AIPL1 protein-protein interactions must be confirmed, their identification is an important step in understanding the functional role of AIPL1 within the retina and will provide insight into the molecular mechanisms underlying inherited blindness. ^

The dynamics of calmodulin signaling revealed using optical methods

The Ca2+-binding protein calmodulin (CaM) is a key transducer of Ca2+ oscillations by virtue of its ability to bind Ca 2+ selectively and then interact specifically with a large number of downstream enzymes and proteins. It remains unclear whether Ca2+ -dependent signaling alone can activate the full range of Ca 2+/CaM regulated processes or whether other regulatory schemes in the cell exist that allow specific targeting of CaM to subsets of Ca 2+/CaM binding sites or regions of the cell. Here we investigate the possibility that alterations of the availability of CaM may serve as a potential cellular mechanism for regulating the activation of CaM-dependent targets. By utilizing sensitive optical techniques with high spatial and temporal resolution, we examine the intracellular dynamics of CaM signaling at a resolution previously unattainable. After optimizing and characterizing both the optical methods and fluorescently labeled probes for intracellular measurements, the diffusion of CaM in the cytoplasm of HEK293 cells was analyzed. It was discovered that the diffusion characteristics of CaM are similar to that of a comparably sized inert molecule. Independent manipulation of experimental parameters, including increases in total concentrations of CaM and intracellular Ca2+ levels, did not change the diffusion of CaM in the cytoplasm. However, changes in diffusion were seen when the concentration of Ca2+/CaM-binding targets was increased in conjunction with elevated Ca2+. This indicates that CaM is not normally limiting for the activation of Ca 2+/CaM-dependent enzymes in HEK293 cells but reveals that the ratio of CaM to CaM-dependent targets is a potential mechanism for changing CaM availability. Next we considered whether cellular compartmentalization may act to regulate concentrations of available Ca2+/CaM in hippocampal neurons. We discovered changes in diffusion parameters of CaM under elevated Ca2+ conditions in the soma, neurite and nucleus which suggest that either the composition of cytoplasm is different in these compartments and/or they are composed of unique families of CaM-binding proteins. Finally, we return to the HEK293 cell and for the first time directly show the intracellular binding of CaM and CaMKII, an important target for CaM critical for neuronal function and plasticity. Furthermore, we analyzed the complex binding stoichiometry of this molecular interaction in the basal, activated and autophosphorylated states of CaMKII and determined the impact of this binding on CaM availability in the cell. Overall these results demonstrate that regulation of CaM availability is a viable cellular mechanism for regulating the output of CaM-dependent processes and that this process is tuned to the specific functional needs of a particular cell type and subcellular compartment. ^

Analysis of Mycobacterium tuberculosis in the state of Texas for rifampin resistance using molecular beacons

Mycobacterium tuberculosis, a bacillus known to cause disease in humans since ancient times, is the etiological agent of tuberculosis (TB). The infection is primarily pulmonary, although other organs may also be affected. The prevalence of pulmonary TB disease in the US is highest along the US-Mexico border, and of the four US states bordering Mexico, Texas had the second highest percentage of cases of TB disease among Mexico-born individuals in 1999 (CDC, 2001). Between the years of 1993 and 1998, the prevalence of drug-resistant (DR) TB was 9.1% among Mexican-born individuals and 4.4% among US-born individuals (CDC, 2001). In the same time period, the prevalence of multi-drug resistant (MDR) TB was 1.4% among Mexican-born individuals and 0.6% among US-born individuals (CDC, 2001). There is a renewed urgency in the quest for faster and more effective screening, diagnosis, and treatment methods for TB due to the resurgence of tuberculosis in the US during the mid-1980s and early 1990s (CDC, 2007a), and the emergence of drug-resistant, multidrug-resistant, and extremely drug-resistant tuberculosis worldwide. Failure to identify DR and MDR-TB quickly leads to poorer treatment outcomes (CDC, 2007b). The recent rise in TB/HIV comorbidity further complicates TB control efforts. The gold standard for identification of DR-TB requires mycobacterial growth in culture, a technique taking up to three weeks, during which time DR/MDR-TB individuals harboring resistant organisms may be receiving inappropriate treatment. The goal of this study was to determine the sensitivity and specificity of real-time quantitative polymerase chain reaction (qPCR) using molecular beacons in the Texas population. qPCR using molecular beacons is a novel approach to detect mycobacterial mutations conferring drug resistance. This technique is time-efficient and has been shown to have high sensitivity and specificity in several populations worldwide. Rifampin (RIF) susceptibility was chosen as the test parameter because strains of M. tuberculosis which are resistant to RIF are likely to also be MDR. Due to its status as a point of entry for many immigrants into the US, control efforts against TB and drug-resistant TB in Texas is a vital component of prevention efforts in the US as a whole. We show that qPCR using molecular beacons has high sensitivity and specificity when compared with culture (94% and 87%, respectively) and DNA sequencing (90% and 96%, respectively). We also used receiver operator curve analysis to calculate cutoff values for the objective determination of results obtained by qPCR using molecular beacons. ^

Macromolecular interaction studies on FF1-3 of human CA150 and STAT1-importin alpha5 using biophysical and biochemical methods

Macromolecular interactions, such as protein-protein interactions and protein-DNA interactions, play important roles in executing biological functions in cells. However the complexity of such interactions often makes it very challenging to elucidate the structural details of these subjects. In this thesis, two different research strategies were applied on two different two macromolecular systems: X-ray crystallography on three tandem FF domains of transcription regulator CA150 and electron microscopy on STAT1-importin α5 complex. The results from these studies provide novel insights into the function-structure relationships of transcription coupled RNA splicing mediated by CA150 and the nuclear import process of the JAK-STAT signaling pathway. ^ The first project aimed at the protein-protein interaction module FF domain, which often occurs as tandem repeats. Crystallographic structure of the first three FF domains of human CA150 was determined to 2.7 Å resolution. This is the only crystal structure of an FF domain and the only structure on tandem FF domains to date. It revealed a striking connectivity between an FF domain and the next. Peptide binding assay with the potential binding ligand of FF domains was performed using fluorescence polarization. Furthermore, for the first time, FF domains were found to potentially interact with DNA. DNA binding assays were also performed and the results were supportive to this newly proposed functionality of an FF domain. ^ The second project aimed at understanding the molecular mechanism of the nuclear import process of transcription factor STAT1. The first structural model of pSTAT1-importin α5 complex in solution was built from the images of negative staining electron microscopy. Two STAT1 molecules were observed to interact with one molecule of importin α5 in an asymmetric manner. This seems to imply that STAT1 interacts with importin α5 with a novel mechanism that is different from canonical importin α-cargo interactions. Further in vitro binding assays were performed to obtain more details on the pSTAT1-importin α5 interaction. ^

Comparison of multiple comparison methods for identifying differential gene expression in simulated and real papillary thyroid cancer microarray data

The difficulty of detecting differential gene expression in microarray data has existed for many years. Several correction procedures try to avoid the family-wise error rate in multiple comparison process, including the Bonferroni and Sidak single-step p-value adjustments, Holm's step-down correction method, and Benjamini and Hochberg's false discovery rate (FDR) correction procedure. Each multiple comparison technique has its advantages and weaknesses. We studied each multiple comparison method through numerical studies (simulations) and applied the methods to the real exploratory DNA microarray data, which detect of molecular signatures in papillary thyroid cancer (PTC) patients. According to our results of simulation studies, Benjamini and Hochberg step-up FDR controlling procedure is the best process among these multiple comparison methods and we discovered 1277 potential biomarkers among 54675 probe sets after applying the Benjamini and Hochberg's method to PTC microarray data.^

Systematic interpretation of molecular beacon polymerase chain reaction for identifying rpoB mutations in Mycobacterium tuberculosis isolates with mixed resistant and susceptible bacteria

Detection of multidrug-resistant tuberculosis (MDR-TB), a frequent cause of treatment failure, takes 2 or more weeks to identify by culture. RIF-resistance is a hallmark of MDR-TB, and detection of mutations in the rpoB gene of Mycobacterium tuberculosis using molecular beacon probes with real-time quantitative polymerase chain reaction (qPCR) is a novel approach that takes ≤2 days. However, qPCR identification of resistant isolates, particularly for isolates with mixed RIF-susceptible and RIF-resistant bacteria, is reader dependent and limits its clinical use. The aim of this study was to develop an objective, reader-independent method to define rpoB mutants using beacon qPCR. This would facilitate the transition from a research protocol to the clinical setting, where high-throughput methods with objective interpretation are required. For this, DNAs from 107 M. tuberculosis clinical isolates with known susceptibility to RIF by culture-based methods were obtained from 2 regions where isolates have not previously been subjected to evaluation using molecular beacon qPCR: the Texas–Mexico border and Colombia. Using coded DNA specimens, mutations within an 81-bp hot spot region of rpoB were established by qPCR with 5 beacons spanning this region. Visual and mathematical approaches were used to establish whether the qPCR cycle threshold of the experimental isolate was significantly higher (mutant) compared to a reference wild-type isolate. Visual classification of the beacon qPCR required reader training for strains with a mixture of RIF-susceptible and RIF-resistant bacteria. Only then had the visual interpretation by an experienced reader had 100% sensitivity and 94.6% specificity versus RIF-resistance by culture phenotype and 98.1% sensitivity and 100% specificity versus mutations based on DNA sequence. The mathematical approach was 98% sensitive and 94.5% specific versus culture and 96.2% sensitive and 100% specific versus DNA sequence. Our findings indicate the mathematical approach has advantages over the visual reading, in that it uses a Microsoft Excel template to eliminate reader bias or inexperience, and allows objective interpretation from high-throughput analyses even in the presence of a mixture of RIF-resistant and RIF-susceptible isolates without the need for reader training.^

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

Genetic epidemiologic methods in risk determination in Li-Fraumeni syndrome

Li-Fraumeni syndrome (LFS) is characterized by a variety of neoplasms occurring at a young age with an apparent autosomal dominant transmission. Individuals in pedigrees with LFS have high incidence of second malignancies. Recently LFS has been found to be associated with germline mutations of a tumor-suppressor gene, p53. Because LFS is rare and indeed not a clear-cut disease, it is not known whether all cases of LFS are attributable to p53 germline mutations and how p53 plays in cancer occurrence in such cancer syndrome families. In the present study, DNAs from constitutive cells of two-hundred and thirty-three family members from ten extended pedigrees were screened for p53 mutations. Six out of the ten LFS families had germline mutations at the p53 locus, including point and deletion mutations. In these six families, 55 out of 146 members were carriers of p53 mutations. Except one, all mutations occurred in exons 5 to 8 (i.e., the "hot spot" region) of the p53 gene. The age-specific penetrance of cancer was estimated after the genotype for each family member at risk was determined. The penetrance was 0.15, 0.29, 0.35, 0.77, and 0.91 by 20, 30, 40, 50 and 60 year-old, respectively, in male carriers; 0.19, 0.44, 0.76, and 0.90 by 20, 30, 40, and 50 year-old, respectively, in female carriers. These results indicated that one cannot escape from tumorigenesis if one inherits a p53 mutant allele; at least ninety percent of p53 carriers will develop cancer by the age of 60. To evaluate the possible bias due to the unexamined blood-relatives in LFS families, I performed a simulation analysis in which a p53 genotype was assigned to each unexamined person based on his cancer status and liability to cancer. The results showed that the penetrance estimates were not biased by the unexamined relatives. I also determined the sex, site, and age-specific penetrance of breast cancer in female carriers and lung cancer in male carriers. The penetrance of breast cancer in female carriers was 0.81 by age 45; the penetrance of lung cancer in male carriers was 0.78 by age 60, indicating that p53 play a key role for tumorigenesis in common cancers. ^

Characterization of fecal contamination targeting 16S rRNA bacteroides molecular markers in the Santa Cruz River, Arizona

Understanding the origins, transport and fate of contamination is essential to effective management of water resources and public health. Individuals and organizations with management responsibilities need to understand the risks to ecosystems and to humans from contact with contamination. Managers also need to understand how key contaminants vary over time and space in order to design and prioritize mitigation strategies. Tumacacori National Historic Park (NHP) is responsible for management of its water resources for the benefit of the park and for the health of its visitors. The existence of microbial contaminants in the park poses risks that must be considered in park planning and operations. The water quality laboratory at the Maricopa Agricultural Center (in collaboration with stakeholder groups and individuals located in the ADEQ-targeted watersheds) identified biological changes in surface water quality in impaired reaches rivers to determine the sources of Escherichia coli (E. coli); bacteria utilizing innovative water quality microbial/bacterial source tracking methods. The end goal was to support targeted watershed groups and ADEQ towards E. coli reductions. In the field monitoring was conducted by the selected targeted watershed groups in conjunction with The University of Arizona Maricopa Agricultural Center Water Quality Laboratory. This consisted of collecting samples for Bacteroides testing from multiple locations on select impaired reaches, to determine contamination resulting from cattle, human recreation, and other contributions. Such testing was performed in conjunction with high flow and base flow conditions in order to accurately portray water quality conditions and variations. Microbial monitoring was conducted by The University of Arizona Water Quality Laboratory at the Maricopa Agricultural Center using genetic typing to differentiate among two categories of Bacteroides: human and all (total). Testing used microbial detection methodologies and molecular source tracking techniques.^

Pathway Semantics: An Algebraic Data Driven Algorithm to Generate Hypotheses about Molecular Patterns Underlying Disease Progression

The overarching goal of the Pathway Semantics Algorithm (PSA) is to improve the in silico identification of clinically useful hypotheses about molecular patterns in disease progression. By framing biomedical questions within a variety of matrix representations, PSA has the flexibility to analyze combined quantitative and qualitative data over a wide range of stratifications. The resulting hypothetical answers can then move to in vitro and in vivo verification, research assay optimization, clinical validation, and commercialization. Herein PSA is shown to generate novel hypotheses about the significant biological pathways in two disease domains: shock / trauma and hemophilia A, and validated experimentally in the latter. The PSA matrix algebra approach identified differential molecular patterns in biological networks over time and outcome that would not be easily found through direct assays, literature or database searches. In this dissertation, Chapter 1 provides a broad overview of the background and motivation for the study, followed by Chapter 2 with a literature review of relevant computational methods. Chapters 3 and 4 describe PSA for node and edge analysis respectively, and apply the method to disease progression in shock / trauma. Chapter 5 demonstrates the application of PSA to hemophilia A and the validation with experimental results. The work is summarized in Chapter 6, followed by extensive references and an Appendix with additional material.

New methods for quantification and analysis of quantitative real-time polymerase chain reaction data

Quantitative real-time polymerase chain reaction (qPCR) is a sensitive gene quantitation method that has been widely used in the biological and biomedical fields. The currently used methods for PCR data analysis, including the threshold cycle (CT) method, linear and non-linear model fitting methods, all require subtracting background fluorescence. However, the removal of background fluorescence is usually inaccurate, and therefore can distort results. Here, we propose a new method, the taking-difference linear regression method, to overcome this limitation. Briefly, for each two consecutive PCR cycles, we subtracted the fluorescence in the former cycle from that in the later cycle, transforming the n cycle raw data into n-1 cycle data. Then linear regression was applied to the natural logarithm of the transformed data. Finally, amplification efficiencies and the initial DNA molecular numbers were calculated for each PCR run. To evaluate this new method, we compared it in terms of accuracy and precision with the original linear regression method with three background corrections, being the mean of cycles 1-3, the mean of cycles 3-7, and the minimum. Three criteria, including threshold identification, max R2, and max slope, were employed to search for target data points. Considering that PCR data are time series data, we also applied linear mixed models. Collectively, when the threshold identification criterion was applied and when the linear mixed model was adopted, the taking-difference linear regression method was superior as it gave an accurate estimation of initial DNA amount and a reasonable estimation of PCR amplification efficiencies. When the criteria of max R2 and max slope were used, the original linear regression method gave an accurate estimation of initial DNA amount. Overall, the taking-difference linear regression method avoids the error in subtracting an unknown background and thus it is theoretically more accurate and reliable. This method is easy to perform and the taking-difference strategy can be extended to all current methods for qPCR data analysis.^

Development of combined micro -preparation, separation, and mass spectrometry methods and application in the microdialysis study of neuropeptide metabolism

Two sets of mass spectrometry-based methods were developed specifically for the in vivo study of extracellular neuropeptide biochemistry. First, an integrated micro-concentration/desalting/matrix-addition device was constructed for matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) to achieve attomole sensitivity for microdialysis samples. Second, capillary electrophoresis (CE) was incorporated into the above micro-liquid chromatography (LC) and MALDI MS system to provide two-dimensional separation and identification (i.e. electrophoretic mobility and molecular mass) for the analysis of complex mixtures. The latter technique includes two parts of instrumentation: (1) the coupling of a preconcentration LC column to the inlet of a CE capillary, and (2) the utilization of a matrix-precoated membrane target for continuous CE effluent deposition and for automatic MALDI MS analysis (imaging) of the CE track.^ Initial in vivo data reveals a carboxypeptidase A (CPA) activity in rat brain involved in extracellular neurotensin metabolism. Benzylsuccinic acid, a CPA inhibitor, inhibited neurotensin metabolite NT1-12 formation by 70%, while inhibitors of other major extracellular peptide metabolizing enzymes increased NT1-12 formation. CPA activity has not been observed in previous in vitro experiments. Next, the validity of the methodology was demonstrated in the detection and structural elucidation of an endogenous neuropeptide, (L)VV-hemorphin-7, in rat brain upon ATP stimulation. Finally, the combined micro-LC/CE/MALDI MS was used in the in vivo metabolic study of peptide E, a mu-selective opioid peptide with 25 amino acid residues. Profiles of 88 metabolites were obtained, their identity being determined by their mass-to-charge ratio and electrophoretic mobility. The results indicate that there are several primary cleavage sites in vivo for peptide E in the release of its enkephalin-containing fragments. ^