Characterization of adjacent breast tumors using oligonucleotide microarrays.

BACKGROUND: Current methodology often cannot distinguish second primary breast cancers from multifocal disease, a potentially important distinction for clinical management. In the present study we evaluated the use of oligonucleotide-based microarray analysis in determining the clonality of tumors by comparing gene expression profiles. METHOD: Total RNA was extracted from two tumors with no apparent physical connection that were located in the right breast of an 87-year-old woman diagnosed with invasive ductal carcinoma (IDC). The RNA was hybridized to the Affymetrix Human Genome U95A Gene Chip (12,500 known human genes) and analyzed using the Gene Chip Analysis Suite 3.3 (Affymetrix, Inc, Santa Clara, CA, USA) and JMPIN 3.2.6 (SAS Institute, Inc, Cary, NC, USA). Gene expression profiles of tumors from five additional patients were compared in order to evaluate the heterogeneity in gene expression between tumors with similar clinical characteristics. RESULTS: The adjacent breast tumors had a pairwise correlation coefficient of 0.987, and were essentially indistinguishable by microarray analysis. Analysis of gene expression profiles from different individuals, however, generated a pairwise correlation coefficient of 0.710. CONCLUSION: Transcriptional profiling may be a useful diagnostic tool for determining tumor clonality and heterogeneity, and may ultimately impact on therapeutic decision making.

A new, non-destructive method for analysis of clinical samples with FT-IR microspectroscopy. Breast cancer tissue as an example.

A new method for infrared analysis of tissues and cells is presented. The method is based on Fourier transform infrared microspectroscopy coupled with attenuated total reflectance. The technique allows spectroscopic measurements on the same samples used by pathologists for histopathological evaluation, e.g. stained samples on plain glass slides. Since the same specimen can be used as for histopathology, the method does not require sample preparation or modification. Significantly, the sample is not damaged. Glass absorbs in the infrared and thus has not been used previously in infrared analysis of tissues and cells. Conventional infrared techniques utilize expensive substrates, such as BaF2 windows and gold coated slides which do not absorb infrared radiation. However, these measurements require special preparation and result in the destruction of the sample. Breast cancer tissues were examined to demonstrate the feasibility and reproducibility of the new method. Linear discriminant analysis was used to discriminate and classify three types of cells: benign, atypical hyperplasia and malignant. It was demonstrated that benign vs. malignant cells were discriminated with 100% accuracy, benign vs. atypical hyperplasia were discriminated with 100% accuracy and malignant vs. atypical hyperplasia were discriminated with an accuracy of 90% and higher.

Qualitative modeling of normal blood coagulation and its pathological states using stochastic activity networks.

We have developed a method for representing biological pathways and simulating their behavior based on the use of stochastic activity networks (SANs). SANs, an extension of the original Petri net, have been used traditionally to model flow systems including data-communications networks and manufacturing processes. We apply the methodology to the blood coagulation cascade, a biological flow system, and present the representation method as well as results of simulation studies based on published experimental data. In addition to describing the dynamic model, we also present the results of its utilization to perform simulations of clinical states including hemophilia's A and B as well as sensitivity analysis of individual factors and their impact on thrombin production.

Qualitative analysis of biochemical reaction systems.

The qualitative analysis of biochemical reaction systems is presented. A discrete event systems approach is used to represent and analyze bioreaction pathways. The approach is based on Petri nets, which are particularly suited to modeling stoichiometric transformations, i.e. the inter-conversion of metabolites in fixed proportions. The properties and methods for the analysis of Petri nets, along with their interpretation for biochemical systems, are presented. As an example, the combined glycolytic and pentose phosphate pathway of the erythrocyte cell is presented to illustrate the concepts of the methodology.

Infrared spectroscopic studies of calcium binding to inhibited beta-trypsins.

Fourier transform infrared spectroscopy was used to examine the effect of calcium binding on the secondary structure of two inhibited bovine beta-trypsins. Neither the diisopropyl fluorophosphate- nor benzamidine-inhibited forms showed detectable secondary structure perturbation upon calcium binding at pD 6.9 and 5.0, respectively. Considered in light of the recent assignment of an amide I' band to the autolysis loop of bovine beta-trypsin, these results contradict the generally held hypothesis that calcium slows trypsin autolysis by induction of a conformational change at this site and support the recent contention that the mechanism of action has a specific electrostatic origin. In addition, the appearance of a band at 1699 cm-1 in the benzamidine-inhibited form can be interpreted as resulting from the NC-N stretching vibrations of the amidinium moiety, which the observed crystal structure indicates is hydrogen-bonded to the carboxyl group of active-site Asp-189.

Petri net representations in metabolic pathways.

The present methods for representing metabolic pathways are limited in their ability to handle complex systems, incorporate new information, and to provide for drawing qualitative conclusions from the structure of pathways. The theory of Petri nets is introduced as a tool for computer-implementable representation of pathways. Petri nets have the potential to overcome the present limitations, and through a multitude of properties, enable the preliminary qualitative analysis of pathways.

Generation of a substructure library for the description and classification of protein secondary structure. I. Overview of the methods and results.

Protein secondary structure has been typically classified into four major classes--alpha-helices, extended strands, reverse turns, and loops. Available methods for secondary structure analysis utilize predefined structure templates to search for structural matches among proteins. By this approach a significant portion of a proteins backbone conformation is assigned to one of a limited number of conformations or, if unassigned, to random coil. To expand our ability to describe protein secondary structure, we have developed an algorithm that operates independently of a predefined structure template. The procedure uses two geometric descriptors, the linear distance and the backbone dihedral angle, to represent the conformation form the alpha-carbon coordinates. The algorithm functions by searching for conformationally equivalent, contiguous fragments without regard to secondary structural classification and is thus independent of the complexity of the backbone fold. The result is a library of conformationally equivalent structure fragments that exhibit some novel characteristics. The library contains features that reproduce the major secondary structure classes as well as defining conformations previously described only as random or undefined conformations. Additionally, the library defines several subclassifications of beta-strands. We present here a validation of this method and a presentation and discussion of the most significant results. In a second study, we report the results of application of this method to spectra-structure correlations in Fourier transform infrared spectroscopy.

Generation of a substructure library for the description and classification of protein secondary structure. II. Application to spectra-structure correlations in Fourier transform infrared spectroscopy.

Fourier transform infrared spectroscopy has become well known as a sensitive and informative tool for studying secondary structure in proteins. Present analysis of the conformation-sensitive amide I region in protein infrared spectra, when combined with band narrowing techniques, provides more information concerning protein secondary structure than can be meaningfully interpreted. This is due in part to limited models for secondary structure. Using the algorithm described in the previous paper of this series, we have generated a library of substructures for several trypsin-like serine proteases. This library was used as a basis for spectra-structure correlations with infrared spectra in the amide I' region, for five homologous proteins for which spectra were collected. Use of the substructure library has allowed correlations not previously possible with template-based methods of protein conformational analysis.

Comparison of various molecular forms of bovine trypsin: correlation of infrared spectra with X-ray crystal structures.

Fourier-transform infrared spectroscopy is a valuable method for the study of protein conformation in solution primarily because of the sensitivity to conformation of the amide I band (1700-1620 cm-1) which arises from the backbone C = O stretching vibration. Combined with resolution-enhancement techniques such as derivative spectroscopy and self-deconvolution, plus the application of iterative curve-fitting techniques, this method provides a wealth of information concerning protein secondary structure. Further extraction of conformational information from the amide I band is dependent upon discerning the correlations between specific conformational types and component bands in the amide I region. In this paper, we report spectra-structure correlations derived from conformational perturbations in bovine trypsin which arise from autolytic processing, zymogen activation, and active-site inhibition. IR spectra were collected for the single-chain (beta-trypsin) and once-cleaved, double-chain (alpha-trypsin) forms as well as at various times during the course of autolysis and also for zymogen, trypsinogen, and beta-trypsin inhibited with diisopropyl fluorophosphate. Spectral differences among the various molecular forms were interpreted in light of previous biochemical studies of autolysis and the known three-dimensional structures of the zymogen, the active enzyme, and the DIP-inhibited form. Our spectroscopic results from these proteins in D2O imply that certain loop structures may absorb in the region of 1655 cm-1. Previously, amide I' infrared bands near 1655 cm-1 have been interpreted as arising solely from alpha-helices. These new data suggest caution in interpreting this band. We have also proposed that regions of protein molecules which are known from crystallographic experiments to be disordered absorb in the 1645 cm-1 region and that type II beta-turns absorb in the region of 1672-1685 cm-1. Our results also corroborate assignment of the low-frequency component of extended strands to bands below 1636 cm-1. Additionally, the results of multiple measurements have allowed us to estimate the variability present in component band areas calculated by curve fitting the resolution-enhanced IR spectra. We estimate that this approach to data analysis and interpretation is sensitive to changes of 0.01 unit or less in the relative integrated intensities of component bands in spectra whose peaks are well resolved.

Domain interaction in rabbit muscle pyruvate kinase. II. Small angle neutron scattering and computer simulation.

The effects of ligands on the structure of rabbit muscle pyruvate kinase were studied by small angle neutron scattering. The radius of gyration, RG, decreases by about 1 A in the presence of the substrate phosphoenolpyruvate, but increases by about the same magnitude in the presence of the allosteric inhibitor phenylalanine. With increasing pH or in the absence of Mg2+ and K+, the RG of pyruvate kinase increases. Hence, there is a 2-A difference in RG between two alternative conformations. Length distribution analysis indicates that, under all experimental conditions which increase the radius of gyration, there is a pronounced increase observed in the probability for interatomic distance between 80 and 110 A. These small angle neutron scattering results indicate a "contraction" and "expansion" of the enzyme when it transforms between its active and inactive forms. Using the alpha-carbon coordinates of crystalline cat muscle pyruvate kinase, a length distribution profile was calculated, and it matches the scattering profile of the inactive form. These observations are expected since the crystals were grown in the absence of divalent cations (Stuart, D. I., Levine, M., Muirhead, H., and Stammers, D. K. (1979) J. Mol. Biol. 134, 109-142). Hence, results from neutron scattering, x-ray crystallographic, and sedimentation studies (Oberfelder, R. W., Lee, L. L.-Y., and Lee, J.C. (1984) Biochemistry 23, 3813-3821) are totally consistent with each other. With the aid of computer modeling, the crystal structure has been manipulated in order to effect changes that are consistent with the conformational change described by the solution scattering data. The structural manipulation involves the rotation of the B domain relative to the A domain, leading to the closure of the cleft between these domains. These manipulations resulted in the generation of new sets of atomic (C-alpha) coordinates, which were utilized in calculations, the result of which compared favorably with the solution data.

Molecular modeling of protein structure and function: a bioinformatic approach.

This paper reports on the data/information structure of macromolecules as it extends beyond the three-dimensional conformation to include functional descriptors of biochemical (in vitro) and biological (in vivo) characteristics and as it contrasts with the limitations imposed by the data reduction and data classification techniques of traditional molecular modeling. Methodologies for structure-function representation are presented which are being incorporated within a knowledge-acquisition expert system. Examples of the bioinformatic approach are presented concerning macromolecular recognition by serine proteases and the use of Fourier transform-infrared (FT-IR) spectroscopy for structural assignment and analysis by a novel structure-perturbation approach.

Structural organization in the serine proteases. I. Macromolecular specificity in limited proteolysis.

We have been developing computational approaches to increase our ability to analyze the growing body of three-dimensional structural data with applications centered about the serine proteases. The emphasis of these approaches is to compare and contrast macromolecules at the separate levels of secondary, tertiary, and quaternary structure. Our assumption is that in functionally related molecules, regions of structural and/or physicochemical similarity will exhibit functional similarity; regions that are different in structure and/or physicochemical properties will function differently and, therefore, be the source of specificity. Based on this assumption, the independent observations from studies of these enzymes in solution and in biological systems are combined with the structural observations from X-ray crystallographic analysis. A goal of the present research effort is to probe the biomolecular architecture of the serine proteases to evaluate the role of the three-dimensional structure beyond that of the active site in determining recognition and reactivity determinants for these enzymes, and to determine those principles that might be applied successfully to other enzyme systems as well. Of particular note has been our observation of a macromolecular recognition surface which occurs as a topographical feature outside of the active site and under evolutionary control to produce specificity towards macromolecular substrates and inhibitors. In addition we have established the important role of conformational changes that occur beyond the active site of an enzyme and differentiate between natural and synthetic inhibitor-enzyme interactions. This suggests that the specificity and reactivity determinants of a macromolecule are derived from its architecture and structural organization.

Determinants of molecular reactivity as criteria for predicting toxicity: problems and approaches.

We discuss the physicochemical basis for mechanisms of action of toxic chemicals and theoretical methods that can be used to understand the relation to the structure of these chemicals. Molecular properties that determine the chemical reactivity of the compounds are proposed as parameters in the analysis of such structure-activity relationships and as criteria for predicting potential toxicity. The theoretical approaches include quantitative methods for structural superposition of molecules and for superposition of their reactivity characteristics. Applications to polychlorinated hydrocarbons are used to illustrate both rigid superposition methods, and methods that take advantage of structural flexibility. These approaches and their results are discussed and compared with methods that afford quantitative structural comparisons without direct superposition, with special emphasis on the need for efficient automated methods suitable for rapid scans of large structural data bases. Quantum mechanical methods for the calculation of molecular properties that can serve as reactivity criteria are presented and illustrated. Special attention is given to the electrostatic properties of the molecules such as the molecular electrostatic potential, the electric fields, and the polarizability terms calculated from perturbation expansions. The practical considerations related to the rapid calculation of these properties on relevant molecular surfaces (e.g., solvent- or reagent-accessible surfaces) are discussed and exemplified, stressing the special problems posed by the structural variety of toxic substances and the paucity of information on their mechanisms of action. The discussion leads to a rationale for the use of the combination of theoretical methods to reveal discriminant criteria for toxicity and to analyze the initial steps in the metabolic processes that could yield toxic products.

Macro-structural organization of phosphoglycerate mutase.

The three-dimensional structure of phosphoglycerate mutase has been analyzed using a contoured distance matrix and by visual inspection using three-dimensional computer graphics. Three folding lobes have been identified and their internal structure tentatively characterized. The active site is located at a lobe interface with a channel providing possible access from above and below. The arrangement of active site residues on two lobes suggests that the active site might be conformationally flexible. The remaining interface not associated with the active site channel appears to be predominately hydrophobic and thus may contribute to inter-lobe stability.