PHYTOMICS: A GENOMIC-BASED APPROACH TO HERBAL COMPOSITIONS

Abstract: The present invention provides the tools and methodologies necessary to guide the standardization of herbal compositions, to determine which specific components of an herbal composition are responsible for any particular biological activity, to predict the biological activities of a particular herbal composition, to determine the relatedness of herbal compositions, and for the development of improved herbal therapeutics. This invention provides the tools and methodologies for creating, maintaining, improving and utilizing Herbal BioResponse Arrays (HBR Arrays), wherein the HBR Arrays constitute data sets associated with particular herbal compositions. The HBR Arrays of the present invention contain gene expression profiles and may also include information on the plant-related parameters of the herbal constituents, marker information collected following the exposure of a biosystem to the herbal composition, and biological response information collected following the exposure of a biosystem to the herbal composition.

Description: PHYTOMICS: A GENOMIC-BASED APPROACH TO HERBAL COMPOSITIONS FIELD OF THE INVENTION This invention relates to herbal compositions. More specifically, this invention provides tools and methodologies for improving the selection, testing, quality control and manufacture of herbal compositions, and to help guide the development of new herbal compositions and identify novel uses of existing herbal compositions. BACKGROUND OF THE INVENTION All publications and patent applications herein are incorporated by reference to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference. Herbal medicine has been in use for centuries by people of Asia and Europe. In the United States (US), herbs have become commercially valuable in the dietary supplement industry as well as in holistic, medicine. Approximately one third ofthe US population has tried some form of alternative medicine at least once (Eisenberg et al, 1993, N. Engl. J. Med. 328:246-252). Botanicals, including herbs,' have also become a focal point for the identification of new active agents to treat diseases. Active compounds, derived from plant extracts, are of continuing interest to the pharmaceutical industry. For example, taxol is an antineoplastic drug obtained from the baric of he western yew tree. It is estimated that approximately 50 percent ofthe thousands of drugs commonly used and prescribed today are either derived from a plant source or contain chemical imitations of a plant compound (Mindell, E.R., 1992, Earl Mindell's Herb Bible, A Fireside Book). Currently, a number of medicinal formulations, food supplements, dietary supplements and the like contain herbal components or extracts from herbs. Herbal medicines have been used for treating various diseases of humans and animals in many different countries for a very long period of time (see, e.g., LA. Ross, 1999, Medicinal Plants ofthe World, Chemical Constituents, Traditional and Modern Medicinal Uses, Humana Press; D. Mplony,,1998, The American Association of Oriental Medicine's Complete Guide to Chinese Herbal Medicine, Berkley Books; Kessler et al., 1996, The Doctor's Complete Guide to Healing Medicines, Berkley Health/Reference Books); Mindell, supra). Herbal Medicines. There are many branches of herbal medicine around the world, such as Ayuryeda, Unani, Sida and Traditional Chinese Medicine (TCM). While modern Western medicine typically consists of administering a single chemical entity capable of intervening, specific biochemical pathways, each formula of TCM typically contains hundreds of chemical entities from several herbs which are designed to interact with multiple targets in the body in a coordinated manner. Although empirical practice contributed in a significant way to the herbal compositipri and prescription of these ancient herbal medicines, they are also supported, to a varying degree, by a set of theories which all are distinct from that of modern Western medicine in tόrms of an'atp'my,. pharmacology, pathology, diagnosis treatment, etc. Among the different herba ήiedicine fields, TCM has developed a more complete set pf theories over several'centύries which have befen well documented and practiced by local physicians caring for a huge population (>1.3 billion people) in greater China and in East Asia including Korea and Japan ... Western medicine' generally uses purified compounds, either natural or synthetic, mostly directed towards a single physiological target. However, the compositions used in TCM are usually composed of multiple herbs and compounds which are aimed at multiple targets in the body based oh unique and holistic concepts. TCM mainly used processed crude natural products, with ' yiarious -.combinations and formulations, to treat different conformations resulting in fewer side effects: . The great potential of TCM has yet to be realized for the majority ofthe world's people-.; ' ,' ' \' ' The herbs in a typical TCM prescription are assigned roles as the principal herb and the secondary herbs, including assistant, adjuvant and guiding herbs. The principal herb produces the leading effects in treating' trie cause or the main symptom of a disease. An assistant herb helps to strengthen the effect ofthe principal herb and produces leading effects in the treatment ofthe accompanying symptoms. There are three types of adjuvant herbs: 1) those that enhance the therapeutic effects of the principal and assistant herbs or treat tertiary symptoms, 2) those that reduce or eliminate the toxicity arid other side effects ofthe principal and the assistant herbs and 3) those which act bή complementary target tissues not specifically affected by the principal herb. A guiding herb directs the effect of other herbs to the affected site and/or coordinates and mediates the " effects of the other herbs in the prescription or formulation. In contrast to most ofthe herbal medicines or supplements that consist of one or more parts of a single plant, the intended effects of TCM, are directed at multiple tissues. For example, a wellrknown TCM recipe, "Ephedra Decoction" used for treating asthma is composed of ephedra, cinnampn twig, bitter apricot kernel and licorice. Ephedra, is the principal herb, which expels cold, induces diaphoresis and facilitates the flow ofthe Lung Qi to relieve asthma, the main symptom. Cinnamon twig, as the assistant herb, enhances ephedra's induction of diaphoresis and warms the Channels to ensure the flow of Yang Qi for reducing headache and pantalgia. Bitter apricot kernel, as the adjuvant herb, facilitates the adverse flow ofthe Lung Qi and strengthens the asthma relief by ephedra. Licorice as the guiding herb moderates the effects of both ephedra and cinnamon to ensure a homeostasis of the vital Qi. While each ofthe four herbs clearly exhibits its respective activity, they complement as well as supplement each other when they are combined. In practice, the principal herb can be prescribed with one pr more secondary herbs, depending on the symptoms at a patient's presentation (Prescriptions of Traditional Chinese Medicine, Chapter One, pplO-16, E. Zhang, editόrin Chief, Publishing House, Shanghai University of Traditional Chinese Medicine, 1998). . , , , ' \ The main theories ' of TCM that guide the treatment of sickness with herbal medicine and other means, such as acupuncture, are 1) the theory of Yin and Yang, 2) the theory of Five Elements, 3) the theory of Viscera and Bowels, 4) the theory of Qi, Blood and Body Fluid, and 5) the theory of Channels and Collaterals. In TCM, the first important aspect of making the proper diagnosis is to ascertain whether the disease is Yiή or Yang. For example, those patients who have a fever, are thirsty, constipated or have a rapid pulse condition are of Yang character. Those individuals who have an aversion to cold, are not thirsty, and diarrhea and a slow pulse condition are of Yin character. The property, flavpt and function of herbs can also be classified according to Ying and Yang theory. For example, herbs of cold and cool nature belong to Ying, while herbs which are warm and hot in, nature belong to Yang. Herbs with sour, bitter and salty flavor belong to Ying, while herbs with pungent, sweet and bland flavor belong to Yang. Herbs with astringent and subsiding function belong to Yin, while herbs with dispersing, ascending and floating function belong to ang. In TCM, the principles of treatment are based on the predominance or weakness of Yin and Yang. Herbs are prescribed according to their property of Ying and Yang and their function for restoring the imbalance ofthe Ying and Yang. In so doing, the benefit of treatment is achieved. According to the theory of Five Elements there are five basic substances that constitute the material world (i.e., wood, fire, earth, metal and water). In TCM, this theory has been used to explain the physiology and pathology ofthe human body and to guide clinical diagnpsis and treatment. Herbal physicians have applied the laws of generation, restriction, subjugation and reverse restriction ofthe five elements to work out many effective and specific treatment regimens, such as reinforcing earth to generate metal (strengthening the function ofthe spleen to benefit the lung), replenishing water to nourish wood (nourishing the essence ofthe kidney to benefit the liver), supporting earth to restrict the wood (supplementing the function ofthe spleen to treat the hyperactivity ofthe liver), and strengthening water to control fire (replenishing the essence ofthe kidney to treat hyperactivity ofthe heart). Specifically, the property of some herbs is assigned to each ofthe five Elements for the purposes of guiding the prescription of a TCM recipe. In TCM, the internal organs ofthe human body are divided into three groups: five Viscera (the Heart, the LiVe , the Spleen, the Lung and the Kidney), Six Bowels (the Gall Bladder, the Stomach, the/large Intestine, the Small Intestine, the Urinary Bladder, and the Triple Warmer), the Extraordinary Organs (the Brain, the Medulla, the Bone, the Blood Vessel, the Gall Bladder, and; the Uterus). In TCM, the Viscera or the Bowel are not only anatomic units, but are also, concepts of physiology and pathology about interactions between different organs. For exarhpϊe, the heart also refers to some ofthe mental functions and influence functions of blood, hair, tongue and skin. Ying- Yang and the Five Elements influence the interactions among these Viscera, Bowels and Organs. The complexity of interplay ofthe theories is used .0 explain the pathology of diseases to which herbs are prescribed, as discussed below. The prescriptiori of herbal medicine in TCM starts with the diagnosis, which consists of four main items: interrogation, inspection, auscμltation and olfaction, pulse taking and palpation. During the interrogation phase, much information is gathered, including the characteristics of the main symptoms. For instance, if the main symptom is characterized by dull pain of epigastric region, ' hich may be relieved by warming and pressing, this suggests the insufficiency ofthe Spleen- Yang. Soreness and weakness ofthe loins and knees, intolerance of coldness with cold extremities manifests a weakness ofthe Kidney- Yang. During inspection, observation ' s- are made for vitality, skin color and the general appearance and the condition ofthe tongue. For example, a pale complexion corresponds internally to the Lung of autumn, whose Qi is dry. This may occur when Yang Qi is lacking and the circulation of Qi and blood is impeded, or when the coldness in the channels and collaterals causes them to contract. in TCM, it is from Qi, blood and body fluid that come energy needed by the Viscera and Bowels, Channels and Collaterals, tissues and other organs for carrying-out their physiological functions; and on which the formation and metabolism of Qi, blood and body fluid depend. Prescriptions of TCM consider the herbal effects on Qi and blood for treatments. TCM holds that Channels, Collaterals and their subsidiary parts are distributed over the entire body. It is through them that herbs exert influence on pathological targets and achieve the improvement of sickness. For example, ephedra acts on the Channels ofthe Lung and Urinary Bladder so as to induce sweat for relieving asthma and promoting diuresis. As noted above, clinical applications of acupuncture are also guided by the theory of Channels and Collaterals. ' / ' . ' In summary, while.the nature or property of each herb in TCM may be assigned as Yin or Yang, and to one ofthe FiVe Elements, they act through Channels and Collaterals and are mediated via Qi, Blood and Fluid to yield therapeutic effects on targets, such as Viscera and Bowels. Pathogenic factors may be. disguised as decoy through the very same systems of Channels and Collaterals to adversely affect the functions of Viscera and Bowels and thus cause sickness. • '. ;'/ > . • . From the foregoing discussion, it is clear that the TCM terminology is as much of a philosophical concept as an anatomical one. For example, the Heart represents a host of tissues, organs or systems in the body that contribute to a function described in TCM. Thus, the concept ofthe Heart requires a multiple dimension data set to describe each concept of TCM. Once this is accomplished, a molecular holistic medicine can be developed. U.S. Regulatory Process. In the US, dietary supplements (such as botanical products, vitamins and minerals, amino' acids and tissue extracts) are regulated under the Dietary Supplement Health and Education Act of 1994 (the DSHE Act). This Act removed the ingredients of dietary supplements from regulation as food additives under the Federal Food, Drag, and Cosmetic Act. ; In addition, the DSHE Act requires that The Food and Drug Administration (FDA) bear.the burden of proof that a marketed dietary supplement presents a serious or unreasonable risk ' under the conditions of use on the label or as commonly consumed. Thus, there are currently no federal regulations that establish specific criteria for purity, identification and manufacturing procedures for dietary supplements. In addition, few published papers on herbal quality have resulted from the establishment ofthe Office pf Alternative Medicine by Congress in 1992 (Angell et al, 1998, N. Engl. J. Med. 339:839-841). At the present time, the FD must approve each one ofthe chemical entities in a drag composition or cocktail, and then clinical trials must be undertaken so as to obtain separate FDA approval for marketing the drag. This process is extremely tedious and costly. A molecular holistic medicine may require a less arduous evaluation since the previous use of a particular herbal composition as a,botanical drag permits clinical trials with multiple chemicals at the outset (i.e., clinical trials using the herbal composition or specific components ofthe herbal composition). Recently, the FDA has approved the testing of some herbal medicines in clinical trials as botanical drags (FDA Guidance on Botanical Drags, April, 1997). While these events represent a positive development for health care in general, it also raises important issues regarding the formulation, manufacturing and quality control of herbal medicines and dietary supplements, iricluding the traditional Chinese medicines. practices (see, e.g., Angell et al, supra). The need to apply scientific testing to the preparation and administration of herbal medicines, and food supplements has been highlighted by several recent reports of toxicity Resulting; from ingesting herb-based formulations. For example, one patient who took an hefbal : based dietary supplement experienced digitalis toxicity (Slifinan et al, 1998, N. Engl. J. Med. '339:806-8111 It was subsequently determined that the herb ingredient labeled as plantain in the supplement was actually contaminated with Digitalis lanata, an herb known to contain at'.least 60 cardiac glycosides. In another instance, an herbal preparation was found to be the cause of chronic lead intoxication in a patient (Beigel et al, 1998, N. Engl. J. Med. 339:827-830). This is not a completely unexpected occurrence since contamination of traditional Asian herbal remedies by lead and other heavy metals is well documented (Woolf et α/., 1994, Ann. Intern. Med. 121:729-735). Characterization of Botanicals. It is well known that the genetic identity (e.g., genera, species, cultivar, variety, clone), age of herbal growth, harvest time, the specific plant part utilized, processing method geog aphical origin, soil type, weather patterns, type and rate of fertilizer, and other growth factors have a great impact on the particular chemical composition of any particular herb "harvested" from any particular area. Increasing numbers of various types of tests have been instituted to assure the consistent quality of herbs used in medicine and as dietary supplements; including inspections at the macro- and microscopic levels as well as a variety of chemical analyses. Recently, high performance liquid chromatography (HPLC) profile of marker molecules in an herbal extract has become one reference standard. However, there are problems with this approach, including that some ofthe bioactive molecules may not adsorb UV or the visible lights for HPLC detection, and the amount of a chemical is not necessarily proportional to its biological potency. For these reasons, herbal manufacturers resort to a practice of mixing raw herbs from different sources to minimize chemical variations. Mass spectrometry (MS) is an analytical method for determining the relative masses and relative abundances of components of a beam of ionized molecules or molecular fragments produced from a sample in a high vacuum. MS, unlike HPLC, is not optical density- dependent. In practice it is iϊsed in conjunction with HPLC or capillary electrophoresis (CE): the HPLC separates the chemicals and the MS then can be used to identify what they are. Commercial systems are available which integrate MS and HPLC for biological uses. Mass spectrometry is limited to samples that are gaseous or volatile at low pressure, or that can be so rendered by derivatization. : . These steps are no longer adequate. Recent publications report a greater variation in the quality of herbs by specific suppliers, and the difficulty of providing biological equivalence of herbal extracts. Furthermore, the correlation between safety and efficacy and chemicals in an herb is not well defined in most cases. Recently, in response to complaints from consumer groups and regulatory agencies (Federal Register, February 6, 1997, Volume 62, No. 25, Docket No. 96M-0417, cGMP in Manufacturing, Packing or Holding Dietary Supplements, Proposed Rules), some herbal rnkriufacturers have begun to implement Good Manufacturing Practice (GMP) which requires stringent controls at all levels. Chemical and spectroscopic methods have been used to characterize the components of herbal medicines and food supplements: For example, three new hederagenin-based acetylated saponins were isolated from the fruits of Gliricidia sepium using these two methods (Kojima et al, 1998, Phytochemistry 48J5):885-888). The botanical sources of Chinese herbal drags in a number of commercial samples were inferred by comparing the contents of some characteristic constituents which were analyzed with high-performance chromatography (HPLC) or capillary electrophoresis (CE) (Shuerin-Jyi Sheu, 1997, Journal of Food and Drag Analysis 5(4):285- 294). For example, the ratio of ephedrine/pseudoephedrine was used as a marker to differentiate Ephedra intermedia from other species; total alkaloid contents were used to distinguish between species of Phellodendron; and the contents of ginsenosides were used to differentiate between species oϊPanax. However, these methods do not provide a direct measurement ofthe effect of 1 the various herbs on the molecular, physiological or morphological responses following human treatment with the herbs. Using gas chromatography-mass spectrometry and atomic-absorption methods, the California Department of Health Sciences, Food and Drag Branch, recently tested Asian medicines obtained from herbal stores for contaminants (R. J. Ko, 1998, N. Engl. J. Med. 339:847). Ofthe 260 products they tested, at least 83 (32 percent) contained undeclared pharmaceuticals or heavy metals, and 23 had more than one adulterant. Using high- performance liquid chromatography, gas chromatography, and mass spectrometry, a commercially available combination of eight herbs (PC-SPES), was found to contain estrogenic organic compounds (DiPaola et al, 1998. N. Engl. J. Med. 339:785-791). The researchers concluded that PC-SPES has potent estrogenic activity and that prostate cancer patients that took PC-SPES, could confound the results of standard therapies and may experience clinically significant adverse effects. Gas chromatography data was also collected for different samples ofthe traditional Chinese medicine 'wei ling xian' and correlated to the antiinflammatory activity of the samples (Wei et al., Study of chemical pattern recognition as applied to quality assessment ofthe traditional Chinese medicine "wei ling xian," Yao Hsueh Pao 26(10): 772-772 (1991)): This study did not provide relevant HBR Array data, such as time course, dose dependent response, control samples to substantiate the differential power of the biomarkers, nor it utilize a reiterative type, of data construction process to establish a comprehensive database for characterizing effects ofthe herbal composition. Changes in protein levels have also been used to characterize the effects of herbal compositions or specific components of herbs. ?For example, the production of granulocyte colony-stimμlating factor (G-CSF) from peripheral blood mononuclear cells was found to vary depending on which specific Chinese herb was added to the culture (Yamashiki et al, 1992, 1 Clin. Lab. Immunol. 37(2): 83-90). Expression of interleukin-1 alpha receptors was markedly up regulated in cultured human- epidermal keratinocytes treated with Sho-saiko-to, the most commonly used herbal medicine in Japan (Matsumoto et al, 1997, Jpn. J. Pharmacol. 73(4):333-336). The expression of Fc gamma 11/111 receptors and complement receptor 3 of macrophages were increased by treatment with Toki-shakuyakusan (TSS) (J. C. Cyong, 1997, Nippon Yakurigaku Za ' sshi 110(Suppl. l):87-92). Tetrandrine, an alkaloid isolated from a natural Chinese herbal medicine, inhibited signal-induced NF-kappa B activation in rat alveolar macrophages (Chen et dl, 1997, Biochem. Biophys. Res. Commun. 231(1):99-102). The herbs Sairei-to, alismatis rhizoma ' (Japanese name "Takusha") and hoelen (Japanese name "Bukuryou'') inhibited the syijthesis and expression of endσthelin-1 in rats with anti-glomeralar basement membrane nephritis (Hattori et al., 1997, Nippon Jinzo Gakkai Shi 39(2):121-128). The increase or decrease in mRNA levels has also been used as an indicator ofthe effect of various herbs and herbal components. Intraperitoneal injection of Qingyangshen (QYS), a traditional Chinese medicine with antiepileptic properties, and diphenylhydantoin sodium reduced alpha- and beta-tublin mRNAs and hippocampal c-fos mRNA induction during kainic acid-induced chronic seizures in rats (Guo et al., 1993, J. Tradit. Chin. Med. 13(4):281-286; Guo et al,.1995, J. Tradit. Chin. Med. 15(4):292-296; Guo et al, 1996, L Tradit. Chin. Med. 16(l):48.-51)i Treatment of cultured human umbilical vein endothelial cells (HUVECs) with the sappnin astragaloside IV, a component purified from Astragalus membranaceus, decreased plasriiiύogen activator inhibitor type I (PAI-1) specific mRNA expression and increased tissue-type plasminogen activator (t-PA) specific mRNA (Zhang et al, 1997, J. Vase. Res. ,34(4):273-280). One component isolated from the root of Panax ginseng was found to be a potent inducer of interleukin-8 (IL-8) production by human monocytes and by the human riionocytic cell line THP-1, with this induction being accompanied by increased IL-8 mRNA expression (Sonoda et al, 1998, Immunopharmacolo y 38:287-294). ■ . , Recent advances in nucleic acid microarray technology enable massive parallel mining of information on gene expression; This process has been used to study cell cycles, biochemical pathways, genome-wide expression in yeast, cell growth, cellular differentiation, cellular responses to a single chemical, compound, and genetic diseases, including the onset and progression ofthe diseases (M. Schena et α/., 1998. TIBTECH 16:301). Because cells respond to the micro-envirόrirhent, changes by changing the expression level of specific genes, the identities of genes expressed in a cell determine what the cell is derived of and what biochemical arid regulatory systems are involved, among other things (Brown et al 1999, Nature genet. 21 (1) supplemerit:33). Thus, cellular gene expression profiles portray the origin, the present differentiation ofthe cell, and the cellular responses to external stimulants. No researchers to date, if any, have attempted to apply these new technologies to study the molecular effects of whole ' herbal treatments and supplements. Some researchers have, attempted to characterize the effects ofthe major active constituents isolated from selected herbs. For example, treatment of HUVECs with notoginsenoside Rl (NR1), purified from Panax notoginseng, resulted in a dose- and time- dependent increase in TPA synthesis (Zhang et al, 1994, Arteriosclerosis and Thromobosis 14(7): 1040-1046). Treatment with NR1 did not change urokinase-type plasminogen activator and PAI-1 antigen synthesis, nor did it effect the deposition of PAI-1 in the extracellular matrix. TPA mRNA increased as much as twofold when HUVECs were treated with NR1, whereas expression of PAI-Ϊ -specific mRNA was not significantly affected by NR1. Since most studies on P. notoginseng have involved its mixture with other herbs, the researchers noted that it was difficult to assess how their results relate to the situation in vivo when is used therapeutically in humans (Id ., at 1045, second column, first paragraph). In addition, since the researchers only studied one riiajor componentOf the herb, it is not possible to ascertain the molecular effect of the whole herb or the interactions among components ofthe herb from this study. Dobashi et al. (1995.; euroscience Letters 197:235-238) studied the effect of two of the main components of saiko agents, a Chinese herbal drag used to treat nephrotic syndrome, bronchial asthma and chronic rheumatoid arthritis. Administration of SS-d increased plasma adrenocorticotropin (ACTH) levels; proopiomelanocortin mRNA levels in the anterior pituitary and the CRF mRNA level in the rat hypothalamus in a dose dependent manner. In contrast, treatment with SS-a tailed to affect the levels of these molecular markers. While this study indicates that administration of SS-d may have an important role in saiko agents-induced CRF release and CRF gene, expression in rat hypothalamus, it fails to address the molecular effect of the herbal medication' as a whole. Kojima et al ( 998. __ioi. Pharm. Bull. 4:426-428) describe the utilization of differential display of mRNA to isolate and identify genes transcriptionally regulated in mouse liver by sho-saiko-to, an herbal medicine used for treating various inflammatory diseases in Japan. These researchers limited ' their study to the use of mRNA differential display techniques in investigating the molecular mechanisms of herbal medicine. It also failed to address effects in multiple organs Of treated animals and did not provide any guidance for quality control, new use, and standardization of effects. In addition, the study failed to analyze the individual components of the herb and compare the individual results with the results obtained using the whole herbal mixture. Ma Ji et al. (1998: Chinese Medical Journal 111(1): 17-23) investigated the therapeutic effect ofthe herb Astragali membranaceus on sodium and water retention in rats experiencing aortocaval fistula-caused experimental congestive heart failure. Chronic heart failure rats with and without Astraglia treatment were compared for changes in various morphological characteristics (e.g., body weight, serum sodium concentration); physiological characteristics (e.g., mean arterial pressure, heart rate, hematoprit and plasma osmolality); mRNA expression levels (e.g., hypothalamic arginine vasopressin (AVP), AVP V^ receptor, renal AVP V 2 receptor, aquaporin-2 (AXP2)) and protein excretion (e.g., plasma atrial monophosphate peptide (ANP) and urinary cyclic guanidino monophosphate (cGMP)). The researchers found that treatment with Astraglia improved cardiac and renal functions, partially corrected abnormal mRNA expressions of the AVP system and AQP2, and improved the renal reaction to ANP. This study did not address using the collected data to guide the development of new formulations or for elucidating the synergistic or other interactions among various herbs in a formula, or validate the differential power ofthe effects for quality control purposes. As shown by the abσve'review of relevant scientific articles, molecular-based technology has not been used to explore and validate cellular and molecular responses in biological systems that are treated' or challenged with multiple chemicals at the same time, such as herbal medicines arid TCM. Furthermore, these recent advances have not been integrated with other technologies and methods to produce a process for the systematic exploration of biological effe'cts of herbal medicines and TCM. SUMMARY OF THE INVENTION This invention provides' the tools and methodologies for creating, maintaining, improving and utilizing Herbal BioResponse Arrays (HBR Arrays), wherein the HBR Arrays constitμte data sets associated with particular herbal compositions. The HBR Arrays ofthe present invention may include^information on'the plant-related parameters ofthe herbal constituents, marker information collected following the exposure of a biosystem to the herbal composition, and biological response information collected following the exposure of a biosystem to the herbal cόrnposition.. The present invention provides the tools and methodologies necessary for establishing standardized HBR Arrays for particular herbal compositions, wherein the standardized HBR Arrays are used as benchmarks by which to evaluate batches of similar or different herbal compositions. The present invention further provides the tools and methodologies necessary to update and maintain the standardized HBR Arrays. Particular embodiments ofthe present invention involve iterative processes whereby data for additional batches ofthe herbal composition, additional plant-related data, additional marker information, and/or additional BioResponse information is periodically added to the standardized HBR Arrays. Thus, the present invention provides the tools and methodologies for creating, maintaining, updating and using HBR Arrays on an oμgoing basis. The present invention provides the tools and methodologies necessary to guide the standardization of herbal compositions; to determine which specific components of herbal compositions are responsible or particular biological activities; to predict the biological activities of herbal compositions; for the development of improved herbal therapeutics; for adjusting or modifying an herbal corriposition; for measuring the relatedness of different herbal compositions; for identifying specific molecules in the batch herbal composition which retain the desired biological activity; for determining which herbal components of a known herbal composition can be eliriiinated from the known herbal composition while maintaining or improving the desired biological activity ofthe known herbal composition; for identifying new uses and previously unknown biological activities for the batch herbal composition; and for using the predicted biological activity ofthe batch herbal composition to aid in the design of therapeutics which include herbal components and synthetic chemical drags, including the design of therapeutics using ;the methods of combinatorial chemistry. More specifically, the present mvention provides methods of establishing standardized Herbal BioResponse Arrays (JΪBR Arrays) for herbal compositions, wherein the methods comprise: ■ ■ ' " ' ' '• . • ' a) selecting a characterize herbal composition; b) exposing a biosysteni td a batch ofthe characterized herbal composition and collecting data on two or more markers, wherein one ofthe markers is a change in gene expression determined through the use' of a nucleic acid microarray, produced by the steps comprising: i) producing a cell bankirig system; ii) profiling the- geήe'exfiression pattern of cells from the cell banking system before and after exposure to the herbal composition; iii) selecting as' markers those genes whose expression levels are changed by exposure to the herbal composition; c) storing the marker data of step b) as a standardized HBR array. The present invention further provides such methods which further comprise exposing a biosystem to one ' or morenatches ofthe herbal composition, collecting the data on one or more BioResponses, and adding the collected BioResponse data io the standardized HBR Array for that herbal composition. . The present invention provides methods of evaluating herbal compositions, wherein the methods comprise exposing a biosystem to a batch ofthe herbal composition and collecting data on two or more markers; and comparing the collected marker data with a standardized ( HBR Array for the same or a substantially same herbal composition as that ofthe batch herbal compositions. The present invention provides a system for predicting the biological activity of an herbal composition comprising: . 1). a biosystem comprising one or more different types of cells, tissues, organs or in vitro assays; . . . ' ; , '- , 2). a batch h rbdl composition; 3). two or more molecular markers; 4). a means for/exposing the biosystem to the batch herbal composition and measuring the differential responses of the molecular markers; 5). a computerprpcessor, including memory, for analyzing and storing the differential response measμrements ofthe molecular markers so as to create an Herbal BioResponse Array (HBR Array) data set for the batch herbal composition; 6). a computed processor, including memory, for comparing the HBR Array ofthe batch herbal composition to one or more previously-stored HBR Arrays so as to predict the biological activity ofthe batch herbal composition, wherein the biological activities ofthe herbal compositions used tp generate thb one or more previously-stored HBR Arrays are known. ' " ' , " . ', .■ ' BRIEF DESCRIPTION OF THE DRAWINGS Figure 1. Figure 1 provides a schematic ofthe basic method steps for constructing a Standardized Herbal BioResponse Array (HBR Array) for any selected herbal composition. The figure is shown in its mos basic form for ease of understanding. As discussed herein, each ofthe pathways of the schernatic can be done iteratively. Furthermore, any information contained in one box can be used to guide decisions regarding gathering information for any other box. In this way, numerous feedback loops are also possible throughout the scheme. Figure 2. Figure 2 provides la schematic ofthe basic method steps for constructing a an Herbal BioResponse Array (HBR Array) for any batch herbal composition and for comparing this batch HRB Array to a selected subset of information from the Standardized HBR Array. The figure is shown in its most basic form for ease of understanding. As discussed herein, each ofthe pathways of the schematic can be done iteratively. Furthermore, any information contained in one box can be.μsed.to guide decisions regarding gathering information for any other box. In this way, numerous feedback loops are possible throughout the scheme. Figure 3. Figure 3 provides a schematic ofthe basic method steps for establishing and using a major data set. The figure is shown in its most basic form for ease of understanding. As discussed herein, each ofthe pathways ofthe schematic can be done iteratively. Furthermore, any information contained in one box can be used to guide decisions regarding gathering information for any other box. In this way, numerous feedback loops are possible throughout the scheme. Figure 4. Western blot for arious herbal compositions. A. No herbal composition. . B. Huang Qing Tang A, (HQT A) (0.2 mg/ml). C. HQT A (4 mg/ml). D. HQT B (0.2 mg/ml). E. HQT B (4 mg/ml). . F. Scute (0.2 mg/ml). ' G. Scute (4 mg/mlj... .. ' _ Figure 5. HPLC for Paeonie lactiflorapallus. Figure 6. HPLC for ZizipHifructus. Figure 7. Figure 7 provides a schematic for establishing a bio-response data set for an herbal composition. The data set is based on differentially expressed gene induced by the herbal medicine for more than three different concentrations in a mammalian cell culture. Figure 8. Figure 8 provides a schematic for establishing a characteristic expression profile database or HBR Array for an herbal medicine or a complex herbal preparation. Figure 9. Figure . provides a schematic for identifying an unknown herbal composition. The expression profiles induced by the unknown herbal medicine are aligned with the expression profile database and statistical method is employed to score the possible identities of herbal medicines archived in the database. Figure 10. Figure 10 provides a schematic for extracting signature genes for an herbal composition or a complex herbal preparation. Figure 11. Figure 11 " rovides a schematic for extracting signature genes for individual chemical constituents in an herbal medicine or a complex herbal preparation. Figure 12. Clustered display of gene expression data from cells treated with three types of single-element herbal extracts (Cordyceps Sinensis Mycelium(CSM), ST024, ST117) with high and low concentrations (indicated with H and L, respectively). (A) Cluster analysis was performed by the program "Cluster" (Eisen et al, 1999) with 492 selected genes (see text). (B) Enlarged image of genes up-regulated by ST 117 treatment but down-regulated by other herbal extract treatments. The clone ID and pμtative gene name are indicated. (C) The clustering algorithm separated CSM, ST024 and ST117 into 3 distinct clusters. The distance between each cluster as displayed by the hierarchical dendrogram can be viewed as the difference between'-tjie expression profiles ofthe three herbal extracts treated cells. . • ■ ■ •, . . ' , ... . Figure 13. ' ■ ' :. ' ;" Λ - , ' : '- - (A) Pseudo-color encoded display of clustering results as calculated based on the selected 492 genes. The boxes in (A) indicate the positions ofthe three clusters of genes described above. ' ■ ". . ' . (B) Enlarged image of genes down-regulated by the CSM but up-regulated by the others. " (C) Genes up-regύlated by all kinds of herbal treatments. (D) Genes dόwn-regύlated by CSM and up-regulated by the others. The IMAGE clone ID and putative- gene name are indicated. Figure 14. Clustered display of expression data from 2 batches of multi-element herbal preparations of the Huang Chin Tang (PHY906-303503 (#11) and PHY906-284003 (#12)) treated cells w ttt mgn and low concentrations (indicated with H and L, respectively). The data were averaged based on three' repeated experiments on three different dates. Cluster analysis was performed based on the selected 500 genes (see text). (B) The clustering algorithm separated #11-L, #11H and (#12-H and #12-L) into 3 distinct clusters. Distance between clusters or resemblance coefficient is indicated by the hierarchical clustering dendrogram. Figure 15. Enlarged image of (A) averaged and (B) individual gene expression levels measured by three independent experiments. Boxl encloses genes that were down regulated in #11-L treated cells but up regulated in others, Box2 encloses the genes that were up regulated by all the herbal treatments. Box3 enclosed the genes that showed no response by #11 -L treatment but were down regulated by the others. Box4 encloses the genes highly down regulated by low concentration herbal treatments but show mild response at high concentration herbal treatrnents. The clone ID arid putative gene name are indicated beside each gene. Figure 16. Classification of gene expression profiles in the cells treated by herbal medicines. Hierarchical clμstering of (A) the data sets normalized with the expression data of the untreated control cells and (B) data sets standardized to have zero-mean and unit- variance. (C) The result of a non-hierarchical flustering by the self-organizing maps algorithm. Figure 17. Candidate class predictors for the classification of herbal medicines based on the gene expression profiles induced by the medicines. 50 class predictors with their expression profiles for discriminating #11 and #12 herbal preparations are shown in this figure. The IMAGE clone ID and pμtative gene name are indicated beside each gene. Figure 18. The gene expression profiles induced by a batch of a complex herbal preparation of five different concentrations. A 6x4 clustering of expression profiles is shown in (A), and the details of thέ gene expression profiles for the selected clusters are shown in (B). Figure 19. Figure 19 illustrates how the expression profiles in Figure 18 are categorized into three different .groups for subsequent hamming distance calculation. Figure 20. Figure 20 shows the analysis results of gene expression profiles induced by five batches of a complex herbal preparation. The numbers in the table are hamming distance. The smaller the distance, the more similar are the expression profiles. Figure 21. Shown in . (A) is a table of integrated peak intensities of 4 chemical constituents in HPLC analyses, of five batches of a complex herbal preparation. Two additional parameters, BG+B and BG/B are introduced td the table and a 6 parameter radial plot is shown in (B) to illustrate that one .batch is more similar to a second batch #18 than to the other batches by the HPLC analysis. Figure 22. A display ofthe signature genes induced by a complex herbal preparation, the Huang Chin Tang, in Jurkat T cells. Figure 23. Figure 23 illustrates the principle of identifying signature genes induced by individual chemical constituents in a mix of herbal rriedicines. The signature genes are those whose expression levels correlate with the amount of chemical constituents in the herbal medicine and that the correlation. cpefficient is larger than 0.99 or smaller than -0.99. (A) shows that the R value between the, gene and Glycyrrhizin was 0.998, and (B) shows that the gene whose, expression levels increase with the decrease of Wogonin has an R value of -0.997. ' ! „ . Figure 24. The signature genes induced by the chemical constituent Albiflorin in a complex herbal preparation, Huang Chin Tang, in Jurkat T cell. (A) show the genes that were positively correlated with Aibiflorin, and (B) shows the genes that were negatively correlated with Albiflorin. .. ' Figure 25. Correlation of gene expression profiles to a control group. (A) is the gene expression profile of a control group, and (B) is the gene expression profile of a sample group. (C) shows the number of genes with a differential expression ratio having greater than 2-fold increase with concentration of herbal treatment. Figure 26. Clusters of expression profiles clustered by a non-hierarchical analysis program, wherein the program is based on a self-organizing map (SOM) principle. The X-axis represents the herbal concentration from low to high and the Y-axis is the gene-expression ratio. Figure 27. Figure 27 shows the induced and repressed genes commonly found in two batches of Huang Chin Tang, Figure 28. SOM clustering results for two batches of Huang Chin Tang. (A) shows the SOM clustering results for the expression profiles of two batches of Huang Chin Tang. (B) shows that ten genes have similarly responded to the two batches, and (C) shows how the weighing factor decreases as cluster I and cluster j become more different. Figure 29. Calculation of S score between pairs of herbal preparations in cluster analysis. (A) is a tabulation of the scores, and (B) is demonstrates how 5 batches of similar herbal preparations are related. DETAILED DESCRIPTION OF THE INVENTION Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although any methods arid materials similar or equivalent to those described herein can be used in the practice or testing ofthe present invention, the preferred methods and materials are described. . , - , Overview of the Invention As set forth above, the present invention is directed to tools and methods useful for predicting the biological response f an herbal composition. More particularly, this invention provides methods of creating Herbal BioResponse Array (HBR Array) databases as well as methods for using such databases to improve the design of effective herbal-based therapeutics. The goal ofthe present invention is the overall design, creation, improvement and use of HBR Arrays for the preparation,.testing and administration of herbal compositions, and guide development of new herbal cornpositipns and novel uses of existing herbal compositions. Phytomics. As used herein, depending on the context in which it is used, "phytomics" refers to using bioinformatics and statistical approaches to address the qualitative and quantitative aspects ofthe components of herbal compositions or to the actual data bases which are developed for addressing such aspects. Herbal BioResponse ' Array. As used herein, an HBR Array constitutes a data set of two or more observations or measurements associated with an herbal composition. The HBR Array may include qualitative arid quantitative data on the plants in the composition (plant- related data), marker information obtained after exposure of a biosystem to the herbal composition including a dose dependent study, and a database of BioResponse data obtained after exposμre of a biosystem to the herbal cόrriposition. The data in any particular HBR Array can be statistically analyzed ri ' either.2- or 3-dimensional space. HBR Arrays may be designated as batch HBR Arrays and standardized HBR Arrays. Batch HBR Arrays are arrays of data associated with specific batches of an herbal composition. Standardized HBR Arrays are arrays of data associated with a standardized herbal composition. ■ • ' ' , Major Data Set. As used herein, the term "major data set" refers to the data set which acts as the baseline set of data by which various other sets of data are compared or otherwise analyzed for the same or different herbal compositions. Generally, the major data set is created using biotechnological techniques to ascertain some genetic or protein aspect ofthe herbal compositions. Thus, the major data set will usually, but not always, be based on a genomic or proteomic set of data. For example, nucleic acid microarray results could be the major data set which is used to compare to other, dependent or minor data sets. Minor or Dependent Data Set. As used herein, the "minor data set" or "dependent data set" refers to one or more data sets which are used for comparing to the major data set. Generally, but hot always, the minor data set will consist of information on an herbal composition which are collected by more traditional methods. For example, the minor, or dependent, data set may consist of a collection of plant-related data obtained by more conventional means. Exaniples of plant-related data include, but are not limited to, the genus/species ofthe herb(s) in the herbal composition, the particular plant parts ofthe herb(s) in the composition and the geographic location where the herb(s) were located. Another example of a minor data set might consist of a set of biological responses of a cell, tissue, organ or organism after treatment with, one or more different amounts ofthe herbal composition. Exaniples of such biological data or a whole organism may include, but are not limited to, cell toxicity studies; enzyme treatment studies, growth rates, weight gain or loss, changes in motor skills and changes in mental abilities. Herb. Technicall speaking an herb is a small, non-woody (i.e., fleshy stemmed), annual or perennial seed-bearing plant in which all the aerial parts die back at the end of each growing season. Herbs are valued for their medicinal, savory or aromatic qualities. As the word is more generally used and as the word is used herein, an "herb" refers to any plant or plant part which has a food supplement, medicinal, drag, therapeutic or life-enhancing use. Thus, as used herein, an herb is not limited to the botanical definition of an herb but rather to any botanical, plant or plant part used for such purposes, including any plant or plant part of any plant species or subspecies ofthe Metaphyta kingdom, including herbs, shrubs, subshrabs, and trees. Plant parts used iri herbal compositions include, but are not limited to, seeds, leaves, stems, twigs, branches, buds, flowers, bulbs, corms, tubers, rhizomes, runners, roots, fruits, cones, berries, cambium and bark. Herbal Composition. As used herein, an "herbal composition" refers to any composition which includes herbs, herbal plants or herbal plant parts. Thus, as used herein, an herbal composition is any herbal preparation, including herbal food supplements, herbal medicines, herbal drugs and medical foods. Examples of herbal compositions include, but are not limited to, the following cpmppnents: a whple plant pr a plant part pf a single plant species; whple plants pr plant parts of multiple plant species; multiple components derived from a single plant species; multiple components derived from multiple plant species; or any combination of these various, components. For a thorough review of various herbal compositions, see, for exahi le, Kee, Chang Huang, The Pharmacology of Chinese Herbs. CRC Press (1993), herein irico orated.in its. entirety. Representative examples of various herbal compositions are provided iri 'the following paragraphs. Herbal compositions \yhich mclude the bark ofthe willow tree have been used to treat fever since the mid-eighteenth century in England. The active ingredient in willow bark is a bitter glycoside called saliciri, which on hydrolysis yields glucose and salicylic alcohol. Aspirin (acetylsalicylic acid)' . and aspirin-like drags (e.g., ibuprofen), all of which are often called nonsteroidal antiinflammatory drugs (NSAIDs), are frequently used to treat pain, fever, and inflammation. Meadowsweet is another herb that contains salicylates. Treatment of arthritic and arthritic-like syriiptoms with willow bark or meadowsweet requires the consumption of large quantities of herbal teas made from these plants. The entire Populus species (i.e., poplar trees and; shrubs) also contains salicylate precursors and poplar-buds have been used in antiinflammatoiy, antipyretic and analgesic medications. U.S. Patents have been.issϊied for herbal compositions used for the treatment of various diseases and other healthτrelated problems afflicting humans and animals. For example, U.S. Patent No. 5,417,979 discloses a composition comprising a mixture of herbs, including species of Stephania and Glycyrrhiza, as well as their extracts, which is used as an appetite stimulant and for the treatment of pain. Herbal compositions which include Glycyrrhiza uralensis have been found useful for treating eczema; psoriasis, pruritis and inflammatory reactions ofthe skin (U.S. Patent No. 5,466,452). U.S. Patent No " . 5,595,743 discloses various herbal compositions which include licorice extract (Glycyrrhiza) and siegesbeckia, sophora, stemona and tetrandra herbs used for the treatment of various mammalian diseases, including inflammation and rheumatoid arthritis. Ocular inflammation can be treated with a pharmaceutical composition containing the plant alkaloid tetrandrine (U.S. Patent No. 5,627,195). U.S. Patent No. 5,683,697 discloses a pharmaceutical composition having anti- inflammatory, anti-fever, expectorant or anti-tussive action, wherein the composition includes plant parts from the, species Meli , Angepica, Dendrobium, Impatiens, Citrus, Loranthus, Celosia, Cynanchum and Glehnia. An herbal composition which includes extracts ofthe roots, rhizomes, and/or vegetation iAlphinia, Smilax, Tinospora, Tribulus, Withania and Zingiber has been found to reduce or alleviate the symptoms associated with rheumatoid arthritis, osteoarthritis, reactive arthritis and for reducing the production of proinflammatory cytokines (U.S. Patent No. 5,683,698). , Herbal compositipns are available in many forms, including capsules, tablets, or coated tablets; pellets; extracts or tinctures? powders; fresh or dried plants or plant parts; prepared teas; juices; creams and ointments; essential oils;. or, as combinations of any of these forms. Herbal medicines are administered by any one of various methods, including orally, rectally, parenterally, enterally, transdermally, intravenously, via feeding tubes, and topically. Herbal compositions encompassed by the present invention include herbal compositions which also contain non-herbal components. Examples of such non-herbal components include, but arenot limited to, whole insects and insect parts, worms, animal or insect feces, natural or petroleu 'oils, carbonate of ammonia, salt of tartar, liquor, water, glycerin, steroids, pharmaceuticals, vitamins, nutrient extracts, whey, salts, and gelatin. For oral administration, the herbal compositions disclosed may take the form of, for example, tablets or capsules prepared by conventional means in admixture with generally acceptable excipients such as binding agents (e.g., pregelatinised maize starch, polyvinylpyrrolidone or hydrbxypropyl methy cellulose); fillers (e.g., lactose, microcrystalline cellulose or calcium phosphate); lubricants (e.g., magnesium stearate, talc or silica); disintegrants (e.g., potato starch or sodium starch glycolate); or wetting agents (e.g., sodium lauryl sulphate); glidanfs; artificial and natural flavors and sweeteners; artificial or natural colors and dyes; and sόlubilizers. The herbal comppsitions may be additionally formulated to release the active agents in a time-release manner as is known in the art and as discussed in U.S. Patent Nos. 4,690,825 and ~ 5,055,300. The tablets may be coated by methods well known in the art. ■ '• ' " ' Liquid preparations for oral administration may take the form of, for example, solutions, syrups, suspensions', of slurries (such, as the liquid nutritional supplements described in Mulchandani et al, 1992 U.S. Patent No. 5,108,767), or they may be presented as a dry product for reconstitution with water or other suitable vehicles before use. Liquid preparations of folic acid, and other vitamins and minerals may come in the form of a liquid nutritional supplement specifically designed for ESRD patients. Such liquid preparations may be prepared by conventional means with pharmaceutically acceptable additives such as suspending agents (e.g., sorbitol syrup, methyl cellulose or hydrogenated edible fats); emulsifying agents (e.g., .lecithin or acacia); non-aqueous vehicles (e.g., almond oil, oily esters or ethyl alcohol); preservatives (e.g., methyl or propyl p-hydroxybenzoates or sorbic acid); and artificial or natural colors and or sweeteners. For topical administration, herbal components may be combined in admixture with at least one other ingredient constituting an acceptable carrier, diluent or excipient in order to provide a composition, such as a creafn, gel, solid, paste, salve, powder, lotion, liquid, aerosol treatment, or the like, which is most suitable for topical application. Sterile distilled water alone and simple cream, ointment and gel bases may be employed as carriers ofthe herbal components. Preservatives and buffers may also be added. The formulation may be applied to a sterile dressing, biodegradable, absorbable patches or dressings for topical application, or to slow release implant systems with a high initial release decaying to slow release. For a more complete overview and discussion of herbal-based compositions see Earl Mindell, Earl Mindell's Herb Bible. Simon & Schuster (1992); Culpeper's Complete Herbal. W. Foulsham & Co., Ltd. (Originally published in the mid 1600's); and, Rodale's Illustrated Encyclopedia of Herbs. Rbdal'e Press (1987). Standardized Herbail Composition. As used herein, a "standardized herbal composition" or a "characterized herbal composition" refers to a particular herbal composition wliich is chosen as the standardrierbal composition for evaluating batch herbal compositions which have the same, similar p different components as the components ofthe standardized herbal composition. Sometimes herein also. referred to as the "master herbal composition." Standardized herbal compositions: are generally herbal compositions which have been well characterized and which demonstrate the desired biological responses in a particular biosystem. Standardized herbal compositions are' usually standardized by chemical tests well known to one skilled in the a t and are properly stored for long term usage and reference. The standardized herbal composition is used to establish a standardized HBR Array based on observations and me'asurerhenis for the plants (i:e., plant-related data), markers and BioResponses so as to characterize the herbal composition. Batch Herbal Composition. As used herein, a "batch herbal composition" refers to any test herbal composition which is used to establish a HBR Array based on observations and measurements for the plants and markers so as to characterize the herbal composition. Sometimes herein also referred to as a "test" or "batch" herbal composition. Observations and measurements of BioResponses may or may not be included. The herbal compositions used to establish the standardized herbal composition may also be referred to as "batch herbal compositions" until designated as "standardized herbal compositions." Batch. As used herein, a "batch" refers to a particular quantity of an herbal composition which can be identified as to some particular attribute so as to distinguish it from any other particular quantity of that same herbal composition. For example, one batch of an herbal cpmppsitipn may differ from, another batch of that same herbal composition in that one ofthe batches was harvested, at. a different time, or in a different geographical location than the other batch. Other differences that distinguish particular batches may include, but are not limited to, the following;' 1-) the particular plant part used (e.g., the root of an herb was used in one batch while the leaves of that sanie herb were used in a different batch); 2) the post-harvest treatment ofthe iridividμal herbs or herbal composition (e.g., one batch may be processed with distilled water while a different batch may be processed with Hydrogen Chloride to simulate the acidity ofthe human stomach); and, 3) the relative proportions ofthe individual herbs in an herbal composition (e.g., one batch may have equal parts by weight or volume of three different herbs while another batch has proportionally more of one herb than the other two). Biosystem. As used herein, a "biosystem" refers to any biological entity for which biological responses may be observed of measured. Thus, a biosystem includes, but is not limited to, any cell, tissue^ organ, whole organism or in vitro assay. Biological Activity; As used herein, the "biological activity" of an herb refers to the specific biological effect peculiar' to an herbal composition on a given biosystem. Plant-Related Data! As used herein, "Plant-related data" refers to the data collected on the herbal composition, including, but not limited to, data about the plants, their growing conditions and the handling of the plants during and after harvesting. The plant-related data also includes the relative propprtioris ofthe components in an herbal compositions, wherein the components may be different plant parts, different plant species, other non-plant ingredients (e.g., insect parts, chemical' drags) or any combinations of these variables. Plant-related data wfricfr may be gathered for an herbal composition includes, but is not limited to, the following:' 1) the plant species (and, if available, the specific plant variety, cultivar, clone, line, etc.) and specific plant parts used in the composition; 2) the geographic origin ofthe herbs, including the longitude/latitude and elevation; 3) the growth conditions of the herbs, including fertilizer types and amounts, amounts and times of rainfall and irrigation, average microEinsteins received per day, pesticide usage, including herbicides, insecticides, miticides and fungicides, and tillage methods; 4) methods and conditions used for processing the herbs, including age/maturity ofthe herbs, soaking times, drying times, extraction methods and grinding methods; and 5) storing methods and conditions for the herbal components and the final herbal composition. Additionally, the standardized herbal composition may be analyzed chemically. Chemical characterization may be accomplished by any chemical analysis method generally known by one skilled in the art. Examples of applicable chemical analyses include, but are not limited to, HPLC, TLC, chQiriical fingerprinting, mass spectrophotometer analyses and gas chromatography. Cell Banking Systerii. .As used herein, a "cell banking system" includes a Master Cell Bank (MCB) and a Working'.Cell Bank (WCB) of cells. The use of a cell banking system minimizes cell variability for herbal medicine testing, and is used for all types of cells in nucleic acid microarray studies. . Bioinformatics. As used herein, "bioinformatics" refers to the use and organization of information of biological interest. Bioinformatics covers, among other things, the following: (1) data acquisition and analysis; (2) database development; (3) integration and links; and (4) further analysis of the .resulting database. Nearly all bioinformatics resources were developed as public domain freeware until the .early 1990s, and much is still available free over the Internet. Some companies^ have developed proprietary databases or analytical software. Genomic or Genόmics. As used herein, the term "genomics" refers to the study of genes and their function. Gehomics emphasizes the integration of basic and applied research in comparative gene mappmg^molecular cloning, large-scale restriction mapping, and DNA sequencing and computational analysis. Genetic information is extracted using fundamental techniques, such as DNA sequencing, protein sequencing and PCR. Gene function is deterrriiried (1) by analyzing the effects of DNA mutations in genes on normal development and health ofthe cell, tissue, organ or organism; (2) by analyzing a variety of signals encoded ' in trie DNA sequence; and (3) by studying the proteins produced by a gene or system of related.genes. Proteoniic or Proteomics. As used herein, the term "proteomics", also called "proteome fesearch" or "pheriόrhe", refers to the quantitative protein expression pattern of a genome under defined conditions. As used generally, proteomics refers to methods of high throughput, automated analysis using protein biochemistry. Conducting proteome fesearch in addition to genome research is necessary for a number of reasons. First, the level of gene expression does not necessarily represent the amount of active protein in a cell. Also, the gene sequence does not describe post-tranlsational modifications which are essential for the function and activity of a protein. In addition, the genome itself does not describe the dynamic cell processes which alter the protein level either up or down. Proteome programs see to characterize all the proteins in a cell, identifying at least part of their amino acid sequence of an isolated protein. In general, the proteins are first separated using 2D gels or HPLC and then the peptides or proteins are sequenced using high throughput mass spectrometry. Using a computer, the output ofthe mass spectrometry can be analyzed so as to link a gerie.arid the particμlar protein for which it codes. This overall process is sometimes referred to as""functiόnal genomics". A number of commercial ventures now offer proteomic services (e.g., /Pharmaceutical Proteomics™, The ProteinChip™ System from Ciphergen Biosystem; PefSepfive Biosystems). For general information on proteome research, see, for example, J.S. Fraton, 1999, Proteins. Enzymes. Genes:. The Interplay of Chemistry and Biology. Yale Univ. Pr.; Wilkins et al., 1997, Proteome Research: New Frontiers in Functional Genomics (Principles and Practice). Springer Verlag; A.J. Link, 1999, 2-D Proteome Analysis Protocals (Methods in Molecular Biology. 112. Humana Pr.; ιKariιp et all, 1999, Proteome and Protein Analysis. Springer Verlag. . ' . ' ' ' ' ' Signal Transduction^ ;As used herein, "signal transduction", also known as cellular sigrial transduction, refers tb the pathways through which cells receive external signals and transmit, amplify and direct .therir internally. Signaling pathways require intercommunicating chains of proteins that, transmit the signal in a stepwise fashion. Protein kinases often participate in this cascade of reactions, since many signal transductions involve receiving an extracellular chemical signal, which triggers the phosphorylation of cytoplasmic proteins to amplify the signal. ..-. • . .. Post-translational -Modification. As used herein, "post-translational modification" is a blanket term used to cover the alterations thathappen to a protein after it has been synthesized as a primary polypeptide. Such post-translational modifications include, but are not limited to, glycosylation,. removal ofthe N-terminal methionine (or N-formyl methionine), signal peptide removal, acetylation, formylation, amino acid modifications, internal cleavage of peptide chains to release smaller proteins or peptides, phosphorylation, and modification of methionine. Array or Microarray.. As used herein, an "array" or "microarray" refers to a grid system which has each position or probe cell occupied by a defined nucleic acid fragment. The arrays themselves are sometimes referred to as "chips", "biochips", "DNA chips" or "gene chips". High-density nucleic acid microarrays often have thousands of probe cells in a variety of grid styles. Once the array is fabricated, ' DNA or protein molecules derived from a biosystem are added and some form.of chemistry occurs between the DNA or protein molecules and the array to give some recognition pattern that is particular to that array and biosystem. Autoradiography of radiolabeled batches is a traditional detection strategy, but other options are available, including fluorescence, colorimetry, and electronic signal transduction. Markers. As used herein, the term "markers" refers to any biological-based measurement or observation fpr a particular herbal composition that is characteristic of a particular biosystem which is being exposed to a particular batch of an herbal composition. The term "marker" encompasses both qualitative and qualitative measurements and observations of a biosystefn.' The marker database constitutes a data set which characterizes gene expression patterns in response to herbal therapies, wherein the patterns show which genes are turned on, off, up, or down in response to specific herbal compositions. Thus, "markers" refers to any bfologically^based measurement or observation whose up- and down- or temporal regulation ' s; of qualitative or quantitative changes of expression levels in a biosystem are used to characterize differential biological responses of a biosystem to an herbal composition. The particular batch of an herbal composition to which the biosystem is exposed may be an unknown herbal composition, a known herbal composition, or a standardized herbal composition. Examples of markers useful in accomplishing the present invention include, but are not limited to, molecular markers, cytogenetic markers, biochemical markers or macromolecular markers. Macrofnόlecular markers include, but are not limited to, enzymes, polypeptides, peptides, sugars-, antibodies, DNA, RNA, proteins (both translational proteins and post-translational proteins),:nucleic : acids, polysaccharides. Any marker that satisfies the definition of "marker" herein is appropriate for conducting the present invention. The term "markers" includes related, alternative terms, such as "biomarker" or "genetic marker" or "gene marker." There may be one or more primary markers along with secondary markers, or a hierarchy of markers for achieving the purposes of increasing the discriminating power of a HBR array. Thus, selected molecular markers may be combined with various other rnolecular, cytogenetic, biochemical or macromolecular markers to enable an even more accurate, extended HBR Array. A molecular marker comprises one or more microscopic molecules from one or more classes of molecular compounds, such as DNA, RNA, cDNA, nucleic acid fragments, proteins, protein fragments, lipids, fatty acids, carbohydrates, and glycoproteins. The establishment;,, generation and use of applicable molecular markers are well known to one skilled in the art. Examples of particularly useful technologies for the characterization of molecular markers include differential display, reverse transcriptase polymerase chain reactions (RT-PCR), large-scale sequencing of expressed sequence tags (ESTs), serial analysis of gene expression (SAGE), Western immunoblot or 2D, 3D study of proteins, and microarray technology. One skilled in the art ..of molecular marker technology is familiar with the methods and uses of such technology; (_ee, e.g. Bernard R. Glick and Jack J. Pasternak, Molecular Biotechnology. Principles and Applications of Recombinant DNA. Second Edition, ASM Press (1998); Mathew R. Walkef arid Ralph Rapley, Route Maps in Gene Technology. Blackwell Science (1997); Roe et al., DNA Isolation and Sequencing. John Wiley & Sons (1996)James D. Watson et al, Recoήibinant DNA.' Second Edition. Scientific American Books (1992)). DNA, RNA and protein isolation and sequencing methods are well known to those skilled in the art. Examples- of such well known techniques can be found in Molecular Cloning: A Laboratory Manual 2nd Edition. Sambrook et al., Cold Spring Harbor, N.Y. (1989); Hanspeter Saluz ari J. P. Jcist, A Laboratory Guide to Genomic Sequencing: The Direct Sequencing of Nativb Uricloned DNA (Biomethods Vol 1). Birkhauser (1988); and B. Roe et al., DNA Isolation and Sequencing. Wiley (1996). Examples of conventional molecular biology techniques include, but are not limited to, in vitro ligation, restriction endonuclease digestion, PCR, cellular trarisformation, hybridization, electrophoresis, DNA sequencing, cell culture, and the like. Specific kits and tools available commercially for use in the present invention include, but are not limited to, those useful for RNA isolation, PCR cDNA library construction, retroviral expresSiori libraries, vectors, gene expression analyses, protein antibody purification, " cytotoxicity assays, protein expression and purification, and high- throughput plasmid purification (see, e.g., CLONTECHniques prpduct catalog, XIϋ(3), 1-32 (1998) or www.clontech.com; Atlas™ cDNA Expression Assays product catalog (1998); SIGMA® product catalog (1997)). ' For discussions, methodologies and applications of oligonucleotide arrays, microarrays, DNA chips or biocriips, see, for example, U.S. Patent Numbers .5,445,934, 5,605,662, 5,631,134, 5,736,257, 5,74i,644, 5,744,305, 5,795,714; Schena et al., Parallel human genome analysis: Microarray-based expression monitoring of 1000 genes, Proc. Natl. Acad. Sci. USA 93, 10614-10619 (1996); DeRisi et al., Exploring the Metabolic and Genetic Control of Gene Expression on a Genomic Scale. Science 278, 680-686 (1997); Wodicka, et al., Genome-wide Expression Monitoring in Saccharomyces cerevisiae, Nature Biotechnology 15, 1359-1367 (1997); Pardee, Complete Genome Expression Monitoring: The Human Race, Nature Biotechnology 15, 1343-1344 (1997); Schafer et al., DNA Variation and the Future of Human Genetics, Nature Biotechnology 16,.3.3-39 (1998); DeRisi et al., Use of a cDNA Microarray to Analyze Gene Expression Patterns in Human Cancer, Nature Genetics 14, 457-460 (1996); Heller et al., Discovery arid Analysis of Inflammatory Disease-Related Genes Using cDNA Microarrays, Prόc. Natl. Acad. ' Sci. USA 94, 2150-2155 (1997); Marshall et al, DNA Chips: An Array of Possibilities. Nature Biotechnology 16, 27-31 (1998); Schena et al., Microarrays: Biotechnology's Discovery Platform for Functional Genomics, Tibtech 16, 301-306 (1998); Ramsay, DNA Chips: State- f-the-art, Nature Biotechnology 16, 40-44 (1998); Chee et al., Accessing Genetic Information with, High-Density DNA Arrays, Science 274, 610-614 (1996); and Chen et al., Profiling Expression Patterns and Isolating Differentially Expressed Genes by cDNA Microarray System with Colori netry Detection, Genomics 50, 1-12 (1998); P. Andrew Outinen et al., Characterization ofthe stress-inducing effects of homocysteine, Biochem. J. 332, 213-221 (1998); and Gilbert et al., Will genetics really revolutionize the drag discovery process, Curr Opin Biotechrioϊ 8(6), 669-674 (1997). Other, more specific, references applicable to the instant invention include, but are not limited to, those addressing the expression technologies, such as ESTs (see, e.g., Michael R. Fannon, Gene expression ϊri normal and disease states - identification of therapeutic targets, TIBTECH 14, 294-298 (1996)); the generation of protein profiles (see, e.g., Robinson et al., A Tyrosine Kinase Profile of Prostate Carcinoma, Proc. Natl. Acad. Sci. USA 93, 5958-5962 (1996)); chemical and ' spectrbscόpic methods for identifying components of herbal compositions (Kojima et al., Saponins from Gliricidia sepium, Phvtochemistry 48(5), 885-888 (1998)); the determination of furictibrial antigens (see, e.g., Aris Persidis, Functional antigenics, Nature Biotechnology 16, 305-307 (1998)); HPLCs (see, e.g., Milton T. W. Hearn (Editor), HPLC of Proteins. Pepties. and Polynucleotides: Contemporary Topics and Applications (Analytical Techniques in Clinical Chemistry and Laboratory Manual). VCH Pub. (1991); electrophoresis (see, e.g., Westermeier et al, Electrophoresis in Practice: A Guide to Methods and Applications of DNA and Protein Separations. John Wiley & Sons (1997)); and cross-reactivity marker assays (see, e.g., Irving Millman et al, Woodchuck Hepatitis Virus: Experimental Infection and Natural Occurrence, Hepatology 4(5):817-823 (1984)). The use of structural genomics for solving the structures of all the proteins encoded for in completed genomes, wherein the.methodology includes high-throughput direct structure determinations and computational methods, is discussed by Terry Gaasterland, Structural genomics: Bioinformatics in the driver's seat, Nature Biotechnology 16, 625-627. For bioinformatics methodologies, see, for example, Andreas Baxevanis (Editor), Bioinformatics: A Practical Guide to the Analysis of Genes and Proteins. John Wiley & Sons (1998) and Luke Alphey, DNA Sequencing: From Experimental Methods to Bioinformatics (Introduction to Biotechniques Series). Springer Verlag (1997). Cytogenetic parameters inclu le, but are not limited to, karyotype analyses (e.g., relative chromosome lengths,, centromere positions, presence or absence of secondary constrictions), ideograms; (i.e. a diagrammatic representation ofthe karyotype of an organism), the behavior of chromosomes during mitosis and meiosis, chromosome staining and banding patterns, DNA-protein interactions (also known as nuclease protection assays), neutron scattering studies, rolling circles (A.M. Diegelriian and E.T. Kool, Nucleic Acids Res 26(13):3235-3241 (1998); Backert et al:, Mol. Cell. Biol. 16(ll):6285-6294 (1996); Skaliter et al, J. Viol. 70(2):1132-1136 (1996); A. Fire and S.Q. Xu, Proc. Natl. Acad. Sci. USA 92(10):4641 -4645 (1995)), and aUtoradiography of whole nuclei following incubation with radiolabelled ribonucleotides. . Biochemical parameters include, but are not limited to, specific pathway analyses, such as signal transduction, protein 'synthesis and transport, RNA transcription, cholesterol synthesis and degradation, glucogenesis arid glycolysis. Fingerprinting. As used herein, the term "finge rinting" as used herein refers to the means of making a characteristic, profile of a substance, particularly an herb, in order to identify it. The term "fingerprint" as used herein refers to the display ofthe result ofthe particular means employed for the fingerprinting. Examples ofthe various types of fingerprinting means include, but are not limited to, DNA fingerprinting, protein fingerprinting, chemical fingerprinting and footprinting. DNA fingerprinting, of profiling, refers to a way of making a unique pattern from the DNA of particular biological source (e.g., a particular plant, plant species, genus of plant, plant part or plant tissue). The DNA fingefprint, or profile, can be used to distinguish that particular biological source from a different biological source. The patterns obtained by analyzing a batch using microarrays, όligόήμcletide arrays, DNA chips or biochips are also referred to as "fingerprints". Protein fingerprinting refers to generating a pattern of proteins in a cell, tissue, organ or organism, such as a plant, which provides a cpmpletely characteristic "fingerprint" of that cell, tissue, organ or organism at that time. Chemical fingerprinting' refers to the analysis ofthe low molecular weight chemicals in a cell and the resulting pattern' used to identify a cell, tissue, organ or organism, such as a plant. The analysis is usually done μsing Gas Chromatography (GC), HPLC or mass spectrometry. Footprinting refers to a r efhod of finding how two molecules stick together. In the case of DNA, a protein is bbμnd to a labeled piece of DNA, and then the DNA is broken down, by enzymes or by chemical attack. This process produces a "ladder" of fragments of all sizes. Where the DNA is protected y the bound protein it is degraded less, and so the "ladder" appears fainter. Footprinting. is a common technique for homing in on where the proteins that regulate gene activity actually bind to the DNA. The means, or methods, used to accomplish each type of fingerprinting are described in detail elsewhere herein. ; . , • BioResponses. As used herein, a "BioResponse" refers to any observation or measurement of a biological response of a biosystem following exposure to an herbal composition. Sometimes herein also referred to as a "biological effect." A BioResponse is a qualitative or quantitative data ppint for the biological activity of a particular herbal composition! BioResponse data includes both dosage and temporal information, wherein such information is well known' to one skilled in the art of measuring responses of biosystems to various treatments. Thus, BipResppnse data includes information on the specific biological response of a specific biosystem to a specific dosage of herbal composition administered in a particular manner for a specific period of time. BioResponses include, but are not limited to, physiological responses, morphological responses, cognitive responses, motivational responses, autonomic responses and post- translational modifications, sμch as signal transduction measurements. Many herbal compositions demonstrate more than one BioResponse (see, e.g., Kee Chang Huang, The Pharmacology of Chinese Herbs, CRC Press (1993)). Some particular BioResponses may be included in more than " pne pf the delineated groups or have aspects or components ofthe response that encompass ore than one group. BioResponses applicable to the instant invention are well known to one skilled in the art. The following references are representative ofthe state of art in the field: Kee Chang Huang, The Pharmacology of Chinese Herbs. CRC Press (1993); Earl Mindell. E rl Mindell's Herb Bible. Simon & Schuster (1992); Goodman & Gilman's The Pharmacological Basis of Therapeutics. Ninth Edition. Joel G. Hardman, et. al. (eds.), McGraw Hill, Health Professions Division (1996); P. J. Bentley, Elements of pharmacology. A primer on drug action. Cambridge University Press (1981); P. T. Marshall and G. M. Hughes, Physiology pf matnmals and other vertebrates. Second Edition. Cambridge University Press (1980); Report ofthe Committee on Infectious Diseases. American Academy of Pediatrics (1991); Knut Schmidt-Nielsen, Animal Physiology: Adaptation and Environment. 5th Edition, Cambridge University ' Press (1997); Elain N. Marieb, Human Anatomy & Physiology. Addison-Wessϊey Pub. 1 Co. (1997); William F. Ganong, Review of Medical Physiology (18th EdV Appleton & Lange (1997); Arthur C. Guyton and John E. Hall, Textbook of Medical Physiology. W. B. Saunders Co. (1995). A "physiological response" refers to any characteristic related to the physiology, or functioning, of a biosystem; Physiological responses on a cellular, tissue or organ level include, but are not limited to, temperature, blood flow rate, pulse rate, oxygen concentration, bioelectric potential^ pH value, cholesterol levels, infection state (e.g., viral, bacterial) and ion flux. Physiological responses; on a whole organism basis include gastrointestinal functioning (e.g., ulcers, upset stomach, ' indigestion, heartburn), reproductive tract functioning (e.g., physiologically-based iriφotence μterine cramping, menstraal cramps), excretory functions (e.g., urinary tract problems, kidney ailments, diarrhea, constipation), blood circulation (e.g., hypertension, heart disorders), oxygen consumption, skeletal health (e.g., osteoporosis), condition ofthe cartilage and connective' tissues (e.g., joint pain and inflammation), locomotion, eyesight (e.g.j ' nϊyppia, blindness), muscle tone (e.g. wasting syndrome, muscle strains), presence or absence of pain, epidermal and dennal health (e.g., skin irritation, itching, skin wounds), functioning ofthe endocrine system, cardiac functioning, nervous coordination, head-related health (e.g., headaches, dizziness), age (e.g., life span, longevity) and respiration (e.g., congestion, respiratory ailments). A "morphological response" refers to any characteristic related to the morphology, or the form and structure, of a biosystem following exposure to an herbal composition. Morphological responses, regardless ofthe type of biosystem, include, but are not limited to, size, weight, height, width, color, degree of inflammation, general appearance (e.g., opaqueness, transparency, -paleness), degree of wetoess or dryness, presence or absence of cancerous growths, and the presence or lack of parasites or pests (e.g., mice, lice, fleas). Morphological responses ' on a ' whole organism basis include, but are not limited to, the amount and location of hair growth (e.g.,'hirsutism, baldness), presence or absence of wrinkles, type and degree of nail and skiri growth, degree of blot clotting, presence or absence of sores or wounds, and presence or absence of hemorrhoids. A "cognitive response" refers to any characteristic related to the cognitive, or mental state, of a biosystem following exposure to an herbal composition. Cognitive responses include, but are not limited tό, ! perceiving, recognizing, conceiving, judging, memory, reasoning and imagining: ' A "motivational response" refers to any' characteristic related to the motivation, or induces action, of a biosysterii following exposure to an herbal composition. Motivational responses include, but are iiot limited, to, emotion (e.g., cheerfulness), desire, learned drive, particular physiological needs, (e.g., appetite, sexual drive) or similar impulses that act as incitements to action (e.g., stamina, sex drive). An "autonomic response" refers to any characteristic related to autonomic responses of a biosystem following exposure to an herbal composition. Autonomic responses are related to the autonomic nervous system ofthe bipsystem. Examples of autonomic responses include, but art not limited to, involuntary functioning (e.g., nervousness, panic attacks), or physiological needs (e.g., respiration, cardiac rhythm, hormone release, immune responses, insomnia, narcolepsy). ; BioResponses of cells, tissues, organs and whole organisms treated with various herbal compositions or herbal components are well known in the herbal arts. For example, the herbal compositions Sairei-tό (TJ-l i,4), alismatis rhizpma (Japanese name 'Takusha') and hoelen (Japanese name 'Bukuryόu') were, each found to inhibit the synthesis and expression of endothelin-1 in rats (Hattori et al, Sairei-to may inhibit the synthesis of endothelin-1 in nephritic glomeruli. Nippon Jinzo Gakkai Shi 39(2), 121-128 (1997)). Interieukin (IL)-1 alpha production was significantly promoted by treatment of cultured human epidermal keratinocytes with the herbal medicine Shp-saiko-to (Matsumoto et al, Enhancement of interieukin- 1 alpha mediated autocrine growth of cultured human keratinocytes by sho-saiko-to, Jpn J. Pharmacol 73(4), 333-336 (1997). Adding Shb-saiko-to to a culture of peripheral blood mononuclear cells obtained from healthy volunteers resulted in a dose-dependent increase in the production of granulocyte colony-stiriiulating factor (G-CSF) (Yamashiki et al, Herbal medicine "sho- saiko-to" induces in vitro granulocyte colony-stimulating factor production on peripheral blood mononuclear cells, J Clin Lab Immunol 37(2), 83-90 (1992)). These researchers concluded that the administration of Shorsaikό-to may be useful for the treatment of chronic liver disease, malignant diseases and acute infectious diseases where G-CSF is efficacious. Plasminogen activator inhibitor type 1 (PAI-1 )-specific nι[_NA expression decreased and tissue-type plasminogen activator (t-PA)-specific mRNA increased after treatment of human umbilical vein endothelial cells (HUVECs) with the saponin astragaloside IV (AS-IV) purified from the Chinese herb Astragalus meriibranaceus' (Zhang et al, Regulation ofthe fibrinolytic potential of cultured human umbilical vein- endothelial cells: astragaloside IV down regulates plasminogen activator iήhϊbitor-l and up regulates tissue-type plasminogen activator expression, J Vase Res 34(4), 273-280 (1997)). One component out of four components isolated from the roots o Panax ginseng was found to be a potent inducer of IL-8 production by human monocytes and THP-1 cells, and this induction was accompanied by increased IL-8 mRNA expression (Sonoda et al; Stimulation of interleukin-8 production by acidic polysaccharides from the root of panax ginseng, Immunopharmacology 38(3), 287-294 (1998)). By flow cytometric analysis, the expression of Fc gamma 11/111 receptors and complement receptor 3 (CR3) oil macrophages were found to be increased by treatment with the Kampo-herbal medicine Toki-shakuyakusan (TSS) (Cyong, New BRM from kampo-herbal medicine, Nippon Yakurigaku " Zasshi 110 Suppl 1, 87P-92P (1997)). Using computer image analysis, Chen et al (Image analysis for intercellular adhesion molecule- 1 expression in MRI/lpr mice: effects of Chinese herb medicine, Chung Hua I Hsueh Tsa Chih 75(4), 204-206 (1995)) found that the distribution intensity of intercellular adhesion molecule-1 (ICAM-1), immunoglobulins and C3.were significantly decreased in MRL/lpr mice after treatment with the Chinese herb stragaliri. Western blot analysis showed that tetradrine, isolated from a natural Chinese herbal medicine, inhibited signal-induced NF-kappa B activation in rat alveolar macrophages (Cheά έt al , Tetrandrine inhibits signal-induced NF-kappa B activation in rat alveolar macrophages, Biocherή Biophys Res Commun 231(1), 99-102 (1997)). Algorithm:. As used herein, an "algorithm" refers to a step-by-step problem-solving procedure, especially an established, recursive computational procedure with a finite number of steps. Appropriate algorithms for two- and three-dimensional analyses ofthe plant-related, marker and BioResponse data sets are well known to one skilled in the computational arts. Such algorithms are useful in constructing the Herbal BioResponse Arrays ofthe present invention. For general information on algorithms, see, for example, Jerrod H. Zar, Biostatistical Analysis, second edition. Prentice Hall (1984); Robert A. Schowengerdt, Techniques for image processing and classification in remote sensing. Academic Press (1983); Steven Gold et al., New Algorithms fpr 2D and 3D Point Matching: Pose Estimation and Correspondence. Pattern Recognition. 31(8):1019-1031 (1998); Berc Rustem, Algorithms for Nonlinear Programming and Multiple-Objective Decisions. Wiley-Interscience Series in Systems and Optimization, John Wiley & Sons (1998); Jeffrey H. Kingston, Algorithms and Data Structures: Design. Correctness. Analysis. International Computer Science Series. Addison- Wesley Pub. Co. (1997); Steven S. Skiena, The Algorithm Design Manual. Springer Verlag (1997); and Marcel F.'Neuts, Algorithm Probability: A Collection of Problems (Stochastic Modeling). Chapman & Hall (1995). For information more specific to the application of algorithms to genetic-based data, see, for example, Dan Gusfield, Algoritlims on Strings. Trees, and Sequences: Computer Science and Computational Biology. Cambridge University Press (1997); Melanie Mitchell, An Introduction to Genetic Algorithms (Complex Adaptive Systems). MIT Press (1996); David E. Goldberg, Genetic Algorithms in Search. Optimization and Machine Learning. Addison- Wessley Pub. Co. (1989); Zbigniew Michalewicz, Genetic Algorithms + Data Stractures = Evolution Programs. Springer Verlag (1996); Andre G. Uitterlinderi'and Jan Vijg. Two-Dimensional DNA Typing: A Parallel Approach to Genome Analysis. Ellis Horwood Series in Molecular Biology, Ellis Horwood Ltd. (1994); and Pierre Baldi and So en Brunak, Bioinformatics: The Machine Learning Approach (Adaptive Computation arid Machine Learning). MIT Press (1998). Combinatorial Chemistry. As used herein, "combinatorial chemistry" refers to the numerous technologies μsed to' create hundreds or thousands of chemical compounds, wherein each ofthe chemical compounds differ for one or more features, such as their shape, charge, and/or hydrophobic characteristics Combinatorial chemistry can be utilized to generate compounds which are chemical variations of herbs or herbal components. Such compounds can be evaluated using the methods of the present invention. Basic combinatorial chemistry concepts are well known to one of ordinary skill in the chemical arts and can also be found in Nichplas K. Terrett, Combinatorial Chemistry (Oxford Chemistry. Masters). Oxford Univ. Press (1998); Anthony W. Czarnik and Sheila Hobbs Dewitt (Editors), A Practical Guide to Combinatorial Chemistry. Amer. Chemical Society (1997); Stephen R. Wilson (Editor) and Anthony W. Czarnik (Contributor), Combinatorial Chemistry: Synthesis and Application. John Wiley & Sons (1997); Eric M. Gordon and James F. Kerwin (Editors), Combinatorial Chemistry and Molecular Diversity in Drug Discovery. Wiley-Liss (1998); Shπiuel Cabilly (Editor), Combinatorial Peptide Library Protocols (Methods in Molecular Biology). Human Press (1997); John P. Devlin, High Throughput Screening. Marcel Dekker (1998); Larry Gold and Joseph Alper, Keeping pace with genomics through combinatorial chemistry; Nature Biotechnology 15, 297 (1997); Aris Persidis, Combinatorial chemistry. Nature Biotechnology 16, 691-693 (1998). EXAMPLES Example 1. Establishing a Standardized HBR Array for Selected Herbal Compositions. ■ ■■■ : .;; ' , . ■ •: '• ': ■ '• ' , The basic scheme for establishing a Standardized HBR Array is provided in Figure 1. Definitions of each component of the schematic are provided above. Following selection ofan herbal composition of interest, data is collected for various traits associated with the herbal composition, including, but not limited to plant-related characteristics and marker and BioResponse information. Plant-related data includes, but is not limited to, the plant species, specific plant parts, geographic origin of the plants in the herbal composition, the growth conditions ofthe plants, the processing methods used to'prepare the herbal components, storage methods and conditions, and various chemical analyses ofthe herbal composition. Marker information includes qualitative and qμarititative data for markers collected after exposure of a biosystem to the herbal cprnppst. Applicable rnaikers include, but are npt limited to, molecular markers, cytogenetic markers, biochemical '.markers and macromolecular markers. BioResponse information includes qualitatiy'e' and quantitative data for biological responses collected after exposure of a biosystem to; the herbal comppsition. Each type of data (e.g. 1 , chemical, marker, BioResponse) can be obtained using one or more assays " on the same,'si_riilar, substantially similar, or different batches ofthe herbal composition of intefest. Such diffefent assays can be conducted at the same or different times. In addition, data can be collected for the same or different markers at the same or different times. Similarly, BioResponse' data can be collected for the same or different biological responses at the same or different times/ Thus, collection ofthe data for the HBR Array is either collected at one time'or ' collected on an on-going basis. Where a biosystem is exposed to an herbal composition so as to collect data, information is recorded on the administered dosages ofthe herbal composition as well as treatment times. BioResponse data may also consist of post-translational modifications, such as measurements of signal transduction. After collection of two or more types of data (e.g., data for two or more markers and a BioResponse; data for plarit-related traits and data for a BioResponse), the data is analyzed using algorithms so as to create 2- and/or 3 -dimensional Herbal BioResponse Arrays. Various statistical parafneters may be calculated for the HBR Array and may become part ofthe HBR Array data set. These statistical parameters may include, but are not limited to, means, standard deviations, correlation or match (or mismatch) matrices, ratios, regression coefficients, and transformed values (e.g., arcsin percentage transformations ofthe raw data). Thus, the HBR Array may consist pf the raw data as well as certain calculations, distributions, graphical presentations and other data manipulations associated with the raw data. Particular examples of such information include, but are not limited to, digital images, scatter graphs, cluster analyses and large scale gene expression profiles for marker data. The total accumulated data and resultant analyses constitute a standardized HBR Array for the particular herbal cόrnposition used to establish the HBR Array data set. Due to the iterative nature of the process rised to establish and maintain an HBR Array for an herbal composition, such arrays can be viewed as either static at any one point in time or dynamic over time. ■ The resulting analyses can identify subsets ofthe standardized HBR Arrays which are correlated (positively or negatively) or associated (i.e., showing a general trend) with one or more specific biological activities cf ariy particular herbal compositipn. Example 2. Establishing a Batch HBR Array for Batch Herbal Compositions. The basic scheme for establishing a HBR Array for a batch of an herbal composition is provided in Figure 2. Definitions of each component ofthe schematic are provided above. The procedure for establishing such an array is the same as that set forth immediately above for the standardized HBR Array; Generally, the amount, of data collected for a batch HBR Array will be less than that collected to establish a standardized HBR Array. However, data collected for a batch herbal composition may be added to an established HBR Array or used to establish a new standardized HBR Array. . "' - Generally, the only data collected for a batch herbal composition is that data which has been found to be highly correlated or associated with the desired biological activities ofthe herbal composition beingiested. Tfdr example, if it has been determined that a particular subset of plant-related and marker data is highly correlated to a desired biological activity of a particular herbal composition (based on the standardized HBR Array data and analyses discussed above), it is only necessary to test the batch herbal composition for that subset of traits in order to determine whether or not the batch has the desired biological activity. By comparing the data obtained-for that subset of traits obtained from the batch (i.e., the batch HBR Array) with the standardized HBR Array. for that particular herbal composition, one skilled in the art can determine whether or not that particular batch has the desired biological activity. • ' ; • Example s. Establishing and Using a Major Data Set. The basic scheme for establishing and using a major data set for an herbal composition is provided in Figure 3. Definitions, of each component ofthe schematic are provided above. The first step, is the establishment of a major data set for a selected herbal composition or batch herbal composition; This is accomplished by exposing a biosystem to the herbal composition and collecting the resultant marker information which will constitute the major data set. In most, but not instances, the major data set will consist of genomics and/or proteomics data in the form of an array, such as an array obtained with a DNA biocbip. Next, the majoi data set is analyzed to see if differential expression/results have been obtained for the tested herbal composition. Differential expression/results are necessary in order to generate meaningful algorithms in the next step. Examples of such differential expression/results include, but are riot limited to, indications that certain genes are up- or down-regulated in respohsέ to exposure to the herbal composition or that the levels of certain proteins have been increased of decreased in response to the exposure. If no meaningful of useful differential expression/results are obtained, then it is necessary to repeat the exposure and marker collection step. If it is believed that experimental error lead to the lack of a adequate result the first time then the exposure/data collection step can be repeated with all ofthe variables the same as the first time (e.g., same biosystem, same marker set, same experimental protocol, etc.). However, it may be necessary to vary the biosystem sampling (e.g., type of cell utilized, stage of cell growth), use a different marker set and/or change the experimental protocol in order to get differential expression/result.. Example 4. Using HBR Array Information. The HBR Array information discussed herein can be used for many different purposes including, but not limited to, the following: 1) evaluating the components of an herbal composition; 2) predicting the BioResponse of an herbal composition; 3) determining which marker information is most highly correlated with a particular BioResponse of an herbal composition; 3) determining what data set of information (i.e., plant-related data, marker data, and BioResponse data) is most correlated with a particular BioResponse of an herbal compost; 4) determining which type of biosystem is best for evaluating the biological activity of an herbal composition; 5) adjusting or changing the components of a herbal composition so that the HBR Array of that herbal composition corresponds to a standardized HBR Array for the same or substantially the same herbal composition; 6) adjusting or changing the components of an herbal composition so that the herbal composition will have the desired biological activity; 7) measuring the relateddess of different herbal compositions; 8) creating and updating standardized HBR Arrays; 9) identifying specific components (e.g., plant parts, proteins, molecules) which retain the desired biological activity of an herbal composition; 10) determining which components of ah herbal composition can be eliminated while maintaining or improving the desired biological activity ofthe herbal composition; 11) identifying one or more previously unknown biological activities for an herbal composition; 12) aiding in the design of therapeutics which include herbal and non-herbal components, such as chemically- synthesized drugs or pharmaceuticals and 13) utilizing the HBR Array information to complement combinatorial chehiistfy methods of designing therapeutics. Each of these embodiments ofthe present invention can be accomplished by one skilled in the applicable art using the methods and tools provided herein. Example 5. Quality Control. The HBR Array technology ofthe present invention is used to correlate or to determine a substantial equivalence of a specific batch of an herbal composition (single herb or multiple herbs of a formula) to a standardized, or master, batch of a same or substantial similar herbal composition. The HBR Ar ays utilized in this process include the acceptable range of quantitative variation for each ofthe biological effects (i.e., BioResponse), and possibly a global score composed of weighted values assigned to each ofthe biological effects, which may consist of markers from multiple biochemical pathways of a biosystem. "Data mining" refers to a'process used to determine or select which subset of biological effects is the minimum riuriiber of biological effects required in any specific HBR Array. The information for data mining results from exposing a biosystem (e.g., a cell line) in a dose dependent manner to a standardized herbal composition to establish a standardized HBR Array. This standardized HBR Array can then be compared to various HBR Arrays established for test herbal compositions. These test herbal compositions include, but are not limited to, different batches prepared at different dates; different batches prepared from raw herbs collected at different times; and different batches prepared from raw herbs collected at different locations. . ' • ' ' Example 6. Improving an Herbal Composition or Identifying New Uses for an Herbal Composition. HBR Arrays are generated by exposing biosystems to either extracts from individual herbs of a formula, or to extracts from the whole formula, and examining the biological effects ofthe extracts. The observed, biological effects can be from multiple biochemical pathways of a biosystem and/or from mμltiple tissues of an animal, wherein various markers are evaluated for their corresponding qualitative and/or quantitative changes. The resulting HBR Arrays can be compared to novel HBR. Arrays or to similar HBR Arrays from different herbal compositions or herbal corripositioris. prepared by different processes. This procedure is useful for selecting a given set of biological effects and the minimum number of markers required to predict that a given batch herbal composition has the given set of biological effects. In order to construct !HBR Arrays, one skilled in the art utilizes various data mining tools including, but are not limited to, statistical analyses, artificial intelligence, and database research on neural work. The statistical methods of choice include, but are not limited to, basic exploratory data analysis (EDA), graphic EDA (such as bushing) and multivariate exploratory techniques (e.g., cluster analysis, discriminating factor analyses, stepwise linear on non-linear regression, classification tree) (see, e.g., STATISTICA™, software packages from StatSoft, Tulsa, OK 74104; Tel: 918-749-1119; Fax: 918-749-2217; www.statsoft.com). Data mining tools are used to explore laf ge amounts of HBR Array data in search of constructing an HBR Array and consistent pattern within, between or among various HBR Arrays. The procedure consists of exploration, construction of an HBR array, and validation. This procedure is typically .repeated iteratively until a robust HBR Array, or standardized HBR Array, is identified. ' . ' Example 7. Establishing a Standardized HBR Array for Ginseng Recipes. For the purposes of this example, standard ginseng is chosen to be Panax Ginseng CA. Meyer Gl 15 grown either in Manchuria or in Korea. The climate for growth is between -10 to +10°C with an annual rainfall of 50-100 cm ' (see Huang in The Pharmacology of Chinese Herbs. (1993) pp 21-45, CRC Press, Boca Raton, FL, fully incoφorated by reference). Ginseng batches will first be characterized by geographic origin, species, plant part (e.g., rhizome, root, leaf skin, seed, bud and flower); growth conditions, processing methods and storage conditions both before and after processing. Verification of chemical content for these batches will be performed by qualitative HPLC analysis for determination of ginsenoside saponins (e.g., Ro, Ral, Ra2, Rbl, Rb2, Rb3, Re, Rgl, Rg2, Rd, Re, Rf, Rhl, Rh2, NG-R2 and Z-Rl), including TLC qualitative analysis for lipophilic constituents (see, Elkin et al, Chumg Kuo Yao Li Hsueh Pao (1993) 14: 97-100 and Yoshikawa et al, Yakugaku Zasshi (1993) 113: 460-467). The saponin content of different herbs should be between 2.1 and 20.6% (by weight) depending on the species (see Table 1). These data will then be stored, preferably in the memory pf a cpmputer processpr, fbr further manipulatipn. Table 1. Saponin Content of Different Ginseng Herbs. Species Total saponins (% by weight) Panax ginseng CA. Meyer 2.1-4.4% Panax quiquefolius , ' ' 4.9% Panax notoginseng and Panax japonica 13.6-20.6% Panax japonica γax. ma]ox ' ■ ' ' 9.34% *from Huang in The Pharmacology of Chinese Herbs. (1993) page 29, CRC Press, Boca Raton, FL. Expression biomarkers for standard ginseng (i.e., G115) include the following: IL-8, IL-2, GM-CSF, NfkB, ICAM-1, interferon gamma, choline acetyl transferase, trk A, nerve growth factor (Kim et α/...PlantaMed (1998) 64: 110-115; Sonoda et αl, Immunopharmacology (1998) 38: 287-294;.Baum et αl, Eur J Appl Physiol (1997) 76: 165- 169; Iwangawa et αl., Free. Radio Biol Med (1998) 24: 1256-1268; Rhind et αl., Eur J Appl Phvsiol (1996) 74: 348-360). Alternatively, for a broader batch size, the 400,000 oligonucleotide group/1.6 cm 2 chip of Affymetrix can be used (U.S. Pat. No.5,556,752). The expression biomarkers fof standard ginseng will be prepared by nucleic acid microarray technology using either phptplithpgraphy, mechanical microspotting or ink jet application (see Schena et αl, TIBTECH (1998) 16: 301-306). .Selected sets of cells will be contacted with standard ginseng for varying periods of time, under varying conditions to generate multiple microarray sets. The microarray set's will then be analyzed by hybridization-based expression monitoring of biochemical extracts via deduction of steady state mRNA levels from fluorescence intensity at each position on the microarrays (Schena et αl., Science (1995) 270: 467-470; Schena et άl, Proc Natl Acad Sci USA (1996) 93: 10614-10619; Lockhart et αl. Nat Biotechnol (1996) 14: 1675-1680; DeRisi et al, Nat Genet (1996) 14: 457-460; Heller et al, Proc Natl Acad Sci USA (1997) . 94: 2150-2155). The array data sets are then input into algorithms to generate statistical expression biomarker values for standard ginseng. Biochemical biomarkers for . standard ginseng include quantitative analysis for increases in cycloheximide sensitive [ H]-leucine incorporation proportional to protein synthesis and [ H]- thymidine incorporation .reflective of mitosis, (see Yamamoto et al, Arzneimittelforschung (1977) 27: 1169-1173). For biochemical biomarkers, bone marrow cells will be contacted with standard ginseng for varying time periods under varying conditions in the presence of [ H]- thymidine (for DNA synthesis) or in the presence and absence of cycloheximide and [ 3 H]- leucine (for protein synthesis) to perform multiple quantitative analysis of biochemical biomarkers (i.e., BBM sets). The BBM sets are then input into algorithms to generate statistical biochemical biomarker Values, fpf standard ginseng. Statistical data will then be stored, preferably in the memory of a computer processor, for further manipulation. Biological response, pf a biosystem (i.e., BioResponses) will be determined using cells and whole animals. For cells, ginseng batches will be exposed to specific cell types, including, but not limited to, fibrόblasts, macrophages, monocytes, PMNL, LAK cells, B16-F10 melanoma cells, THP-1 cells- and hippocampal neurons at a concentration of 0.5 mg/ml to 100 mg/ml. For animal treatments, ' 0.5-100 mg/kg of ginseng herbal extract will be administered orally, by intraperitoneal injection or subcutaneous injection. To determine a biological response of a biosystem to standardized ginseng, human ovarian cancer cells will be inoculated into nude mice, which results in the formation of palpable tumors. After tumor' formation the mice will be treated by co-administration of cis- diamminecichloropiatinum arid standard ginseng. Mice will be examined for tumor growth inhibition, increase in survival time and lowered adverse side-effects on hematocrit values and body weight (Nakata et al, ' Jpri J Cancer Res (1998) 89:733-740). The assay will be repeated using various concentrations of standard ginseng to generate measures of central tendency, dispersion and variability for each variable. The data collected 'will then be subjected to multidimensional analysis to generate multivarianf normal distribution sets as a means of determining a baseline correlation between biological activity and standard ginseng (see Zar, J. H., in Biostatistical Analysis. 2nd ed. (1984), pp 328-360, Prentice-Hall, Englewood Cliffs, NJ). A second independent determination of a biological response of a biosystem to standard ginseng will be the effect of standard ginseng on physical perforrriahce during exercise. Rats will be treated for 4 day ' s with standard ginseng at various " concentrations (between 0.5-100 mg/kg/day) and animals will be tested for increased plasma free fatty acid level and maintenance of glucose level during exercise at approximately 70% VO2max (see Wang et al, Planta Med (1998) 64130-133). The data generated will be collected and then subjected to multidimensional analysis to generate multivariant normal distribμtion sets as a means of determining a baseline correlation between biological activity and standard ginseng (see Zar, J. H., in Biostatistical Analysis. 2nd ed. (1984), pp 328-360, Prentice Hall, Englewood Cliffs, NJ, Herein, fully incorporated by reference). The distribution sets for each BioResponse are then put into algorithms to generate statistical values for standard ginseng. Statistical data will then be stored, preferably in a memory of a computer processor, for further manipulation. Each of these steps (i.e., chemical analysis, generation of biomarker information and determination of responsesjof a biosystem) is reiterated to generate a large database of statistical values. These values are compiled and input into an algorithm to generate 2- and 3- dimensional Herbal Response' rrays (HBR Array) for standardized ginseng. Through reiteration, the resulting arrays ' (t.e. ' ,' Standardized Arrays) display the highest correlation between composition (including growth conditions), biomarker information and biological response for standardized ginseng.' By determining two or more known associated variables for composition and biomarker information values via display on an HBR Array for a test batch, the values for biological response variables can be predicted for the test batch by comparing test values against Standardized HBR Array values for standardized ginseng. The resulting prediction will be used to evaluate the quality of a given ginseng batch without necessitating the use of an observed biological response of a biosystem (see Example 2). Example 8. Evialuation of a Selected Herbal Composition of Ginseng Using a Subset of Variables Correlated with a Specific Biological Response. To evaluate the quality of 'a test batch herbal composition, data is first collected concerning the plarit-related.parametefs for the herbs in the selected herbal composition (e.g., plant species, plant parts, ge ' pgfaphic origin, growth conditions, processing methods and storage conditions). The selected herbal composition is then manipulated such that chemical analysis can be performed to ' determine the chemical content ofthe herb (see Elkin et al, Chumg Kuo Yao Li Hsueh Pab (1993)- 14: 97-100 and Yoshikawa et al, Yakugaku Zasshi (1993) 113: 460-467). Previously obtained ginseng data has demonstrated a strong correlation between oxygen consumption during aerobic exercise performance and the presence of a subset of saponin components, especially Rgl and Rbl (Wang et al, Planta Med (1998) 64: 130-133). '" _ The test batch is then exposed to test cells including, but not limited to, fibroblasts, macrophages, monocytes, PMNL, LAK cells, B16-F10 melanoma cells, THP-1 cells and hippocampal neurons at a concentration of 0.5 mg/ml to 100 mg/ml to determine expression biomarker values. - mRNA is isolated from exposed cells which is subsequently manipulated to serve as a substrate for hybridization-based expression monitoring of biochemical extracts using microarrays comprising IL-8, IL-2 and Interferon gamma cDNA (Schena et al, Science (1995) 270: 467-470; Schena et al', Proc Natl Acad Sci USA (1996) 93: 10614-10619; Lockhart et al, Nat Biotechnol (1996) 14: 1675-1680; DeRisi et al, Nat Genet (1996) 14: 457- 460; Heller et al, Proc Natl Acad Sci USA (1997) 94: 2150-2155). Previously obtained ginseng data has demonstrated a strong correlation between oxygen consumption during aerobic exercise performance.and the induction ofthe expression biomarkers IL-8, IL-2 and Interferon gamma in test cells (Venkatraman et al, Med Sci Sports Exerc (1997) 29: 333-344 and Wang et al, Planta Med (1998) 64: 130-133). For biochemical biomarkers, rat bone marrow cells will then be exposed to the test batch and assayed for [ H]-thymιdine incorporation reflective of mitosis. Previously obtained ginseng data has demonstrated that Rbl and Rgl show a strong correlation with DNA synthesis in rat bone marrow cells (Yamamoto et al, Arzneimittelforschung (1978) 28: 2238-2241). After reiterative analysis, data from each assay will be input into an algorithm to generate a test HBR areay- or the selected herbal composition based on the enumerated plant- related data, including chemical analyses, and data concerning the subset of biomarkers. The quality of a test batch will be determined by comparing test HBR and standard ginseng Standardized HBR Array variables directed toward arialysis ofthe above observations and subsets, wherein the demonstration of the induction of IL-2, IL-8 and INF gamma mRNA in vitro and an increase in [ 3 H],-thyriιidine incorporation in rat bone marrow cells (including data collected on growth conditions, origin, and verification ofthe saponins Rgl and Rbl) is predictive of an equivalent BioResponse effect ofthe test batch on oxygen consumption as that exhibited by standard ginseng. Based on this procedure it can be determined whether or not the test batch is of a similar of- different quality than that of the standard for the given biological response or biological fesponse of interest. Exampϊe 9. Establishing a Standardized HBR Array for Huang Ling (HL) Recipes. For the purposes of this example, standard huang ling (HL) is chosen to be Coptis chinesis France, from soμthwest Asia, wherein growth conditions are well known to one skilled in the art (see Hμang in The Pharmacology of Chinese Herbs. (1993), pp 69 and 287- 288, CRC Press, Boca Raton, FL).. Dried rhizomes of Coptis chinesis France will be verified for chemical content by quantitative chemical analysis for determination arsenic, berberine, caeraleic acid, columbamine, cppsine, cpptine, coptiside-I, coptiside-II, coptisine, coreximine, epiberberine, ferulic acid, greerilandicine, isocoptisine, lumicaerulic acid, magnoflorine, oxybererine, thalifendine, umbellatine, urbenine, worenine, palmatine, jatrorrhizine and colubamine (see also Zhu M., Chung Yao Tung Pao (1984) 9: 63-64). Content ofthe alkaloid berberine of different herbs ' should be between 7-9% (by weight). These data will be stored, preferably in the mempry όf a computer processor, for further manipulation. Expression biomarkers -for standard HL include the following: Nf B; bcl-2 analog, Al; zinc finger protein, A20; ILr2 receptor; cell cycle probes; c-Ki-ras2; growth regulators probes and glucocorticoid receptor dependent apoptosis probes (see Chi et al, Life Sci (1994) 54: 2099-2107; Yang et al. Nauhyn Schmiedebergs Arch Pharmacol (1996) 354: 102-108; Miura et al.. Biochem Pharmacol (l ' 997) 53; Chang K.S., J Formos Med Assoc (1991) 90: 10-14). Alternatively, for a broader batch size, the 400,000 oligonucleotide group/1.6 cm 2 chip of Affymetrix can be used (U.S. Pat. No.5,556,752). The expression biomarkers for standard HL will be prepared by microarray technology as described in Example 1, including analysis and statistical data genefation. Biochemical biomarkers for standard HL include increase in glucocorticoid receptor and inhibitiph of alpha-fetoprotein secretion in HL exposed HepG2 cells (see Chi et al, Life Sci (1994^ 54: 2099-2107). BBM sets are generated and analyzed as described in Example 1. Statistical data will then be stored, preferably in the memory of a computer processor, for further manipulation. Biological response of a biosystem will be determined using cells and whole animals. Batches ofthe selected herbal composition will be exposed to specific cell types, including but not limited to, human HepG2 hepatoma cells, human embryonal carcinoma cells and thymocytes at concentrations ! from 0.1-100mg/ml. For animal treatments 0.1mg-2g/kg of Coptic herbal composition (i.e., HL) will be administered orally, by intraperitoneal injection or subcutaneous injection. To determine a biological response of a biosystem to standardized HL, human embryonal carcinoma _;lόne, NT2/D1 is exposed to various concentrations of standard HL and cells will be examined for differentiation into cells with neuronal-like cell morphology (Chang K.S„ J Formos Med Assoc (1991) 90: 10-14). The assay will be repeated to generate measures and analysis will be performed as described for ginseng in Example 1. A second independent determination Of a biological response of a biosystem to standard HL will be the effect of standard HL on diarrhea due to enterόfoxigenic Escherichia coli (ETEC). Patients with active diarrhea due tp ETEC will be treated with varipus cpncentratipns of HL (e.g., 2g/kg) and stool volumes will be determined (see, e.g., Rabbani G.H., Dan Med Bull (1996) 43: 173-185). The assay will be repeated to generate measures and analysis will be performed as described for ginseng in Example 1. The distribution sets for each biological system are then put into algorithms to generate statistical values for standard HL. Statistical data will then be stored, preferably in the memory of a computer processor, for further manipulation. Lastly, as in Example.!, the steps are reiterated to generate HBR arrays for standardized HL, wherein the resulting HBR arrays will then be used to predict biological activity and evaluate batch quality. Using this method, a Standardized HBR Array can be generated and updated pefipdically. Example 10. Evaluation of a Selected Herbal Composition of Huang Ling Using a Subset of Variables Correlated with a Specific Biological Response. To evaluate the quality of a selected test batch of an herbal composition of Huang Ling, data is first collected con erriing the plant-related characteristics (e.g., plant species, plant parts, geographic origin, growth conditions, processing methods and storage conditions). The herbal composition is then manipμlated such that chemical analysis can be performed to determine the chemical content ofthe composition (see also Zhu M.. Chung Yao Tung Pao (1984) 9: 63-64). ' ' , ' / \ Previously obtained HL data has demonstrated terminal differentiation of human embryonal carcinoma clpries intp neuronal-like cells is strongly correlated with the presence of berberine (see Chang K.S., J Formos Med Assoc (1991) 90: 10-14). The test batch is then exposed to test cells Including liuman 'embryonal carcinoma clone, NT2/D1 at a concentration starting at a non-toxic cOiicentf ation (determination of which is within the skill ofthe ordinary artisan). mRNA is isolated from exposed cells which is subsequently manipulated to serve as substrate for hybridization based expression monitoring of biochemical extracts using microarrays comprising IL-2 receptor and NfkB; (see Chi et al, Life Sci (1994) 54: 2099- 2107; Yang et al, Naunyn Schmie ' debergs Arch Pharmacol (1996) 354: 102-108; Miura et al., Biochem Pharmacol (1997) 53; Chang K.S., J Formos Med Assoc (1991) 90: 10-14; U.S. Pat No.5,556,752), and which can be used to determine down regulation of c-Ki-ras2 gene expression in said cells. Previously obtained HL data has demonstrated terminal differentiation of human embryonal carcinoma clones into neuronal-like cells is strongly correlated with induction of m ' itόgen probes and down regulation of c-Ki-ras2 gene expression (see Chang K.S., J Formos Med Assόc (1991) 90: 10-14). For biochemical markers, HepG2 cells are exposed to the test composition and cells are assayed for increase in glucocorticoid receptor and inhibition of alpha-fetoprotein secretion (see Chi et al, Life Sci (1994) 54: 2099-2107). Previously obtained HL data has demonstrated that inhibition of glucocorticoid induced apoptosis is strongly correlated with berberine-type alkaloids (see Miura et al.;. Biochem Pharmacol (1997) 53 : 1315-1322). After reiterative analysis, data from each assay will be input into an algorithm to generate a test HBR array based on the enumerated observational data, chemical data and data concerning the subset of biomarkers. \ The quality of a. test batch will be determined by comparing test HBR and standard HL HBR Array variables directed toward analysis ofthe above observations and subsets, wherein the demonstration of the induction of IL-2 receptor and NfkB, the down regulation of c-Ki- ras2 gene expression, an increase in glucocorticoid receptor and inhibition of alpha-fetoprotein secretion for HepG2 cells (tp including data cpllected pn growth conditions, origin, and verification of berberine alkaloid) is predictive of an equivalent BioResponse effect ofthe test batch on terminal diffefentiatiori of human embryonal carcinoma clones into neuronal-like cells and inhibition of dexamethasone induced apoptosis as that exhibited by standard HL. Based on this procedure it can be determined whether or not the test batch is of a similar or different quality than that of the. standard. Example 11. Evaluation of Xiao Chai Hu Tang (sho-saiko-to) Using Two Bioassays. To evaluate the . quality of three sources of Xiao Chai Hu Tang, two bioassays were used: 1) cell growth inhibitiori'arid 2) hepatitis B virus secretion from infected cells. The Xiao Chai Hu Tang composition is made f ofn a mixture of 6-7 herbal plants (Radix Bupeuri, Rhizoma Pinelliae, Rhizoma ' Zingiberis, Radix Scutellariae, Fructus Ziziphi, Codonopsis Pilosula, Radix Ginseng and Radix Glycyrrhizae, see Table 2 for relative amounts, by weight). The three "recipes" 'originate in either Singapore, Korea or Taiwan. Batches were evaluated for toxicity and : fόr the ability to inhibit hepatitis B virus as detected by DNA quantitation or detection of hepatitis B surface antigen (HbsAg) (see Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499\ - Briefly, pne gram of preparation was added with 10 ml of water. The mixture was boiled for 30 minutes. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter, wo cell types were used: a) 2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005- 1009) and b) HepG2 cells (ATCC cat # HB-8065). One to fifty dilutions were used for each assay. The cell growth inhibition assay was performed for 72 hours. All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499. The results ofthe assays using the three batches is displayed in Table 3. Based on these data, the Taiwan source would be selected as a standard herbal composition because of its low toxicity combined with its effectiveness iri reducing secretion HbsAG (which is proportional to viral release) by more thari half. Table 3. Bioassay of Xiao Chai Hu Tang (Sho-saiko-to). The data presented in Tables 2 and 3 for the Taiwan herbal composition constitute the initial data for the standardized HBR Array for this herbal composition. Therefore, this data set would initially include the source ofthe herbal composition, the plant species and relative amounts of each the herbal, composition, and two BioResponses (i.e., cell growth inhibition and hepatitis B virus secretion froin infected cells). Using the procedures set forth in the schematic of Figure 1 and in Examples 1 and 3, additional data can be collected pn plant-related data, markers and BioResponses for the standard herbal composition. This additional data is added to the initial standardized HBR Array to generate an expanded standardized HBR Array. Appropriate analyses ofthe resulting database can be conducted as set forth in the detailed description and the examples in order to ascertain the subset of variables which is most highly correlated or associated with the BioResponse of interest. Batch HBR Arrays may be determined using the methods depicted in Figure 2 and in the procedures of Examples 2 and 4. The resultant batch HBR Array can be compared to the standardized HBR Array so as to predict the BioResponse of the'batch herbal compositions. Example 12. Herbal Preparation The standardized protocol Tor the herbal extract preparatipn was as follows: The ingredients of herbal raw materials, with proper ratios were placed in a jacketed reactor and extracted with water at an elevated constant temperature with mixing. The solid was separated from the liquid with a„120-mesh screen. The resultant filtrate was collected and then concentrated by evaporating the water under reduced pressure. The concentrated liquor was spray dried at elevated temperature to yield granulated powder. This bulk substance was then formulated into the desired dosage form. Example 13. Evaluation of Huang Qing Tang Huang Qing Tang (HQT) is an ancient Chinese botanical formula composed of four distinct herbs: Scutellariae (scute), Glycyrrhizae (licorice), Paeonie lactiflora pallus (whitepeony root), andFructus ziziphό (date): (Table 4). This herbal formula has been long used in Asia to treat a variety of gastrointestinal ailments since 300 AD. Table 4. Herbal Ingredients of TCM Formula HQT Biological and Enzyme Assays Table 5. Batch Properties (HQT) Briefly, one gram of each batch of Huang Qing Tang (HQT) was added with 10 ml of water (1 mg/ml). The jnixture was treated as oμtlined in Table 5. The supernatant was collected after cenfrifugation and filtered through a 0.22 μm filter. Two cell types were used to test for biological effects of each batch of HQT: a) Jurkat T cells (ATCC cat #ΗB-152) and b) HepG2 cells (ATCC cat # HB-8065). One to fifty dilutions were used for each assay. Frozen cells (10 7 /ml) were quickly thawed in a water bath at 37 °C. The cells were then diluted in 10 ml of pre-warmed media (see Life Technologies, Inc., Catalogue and Reference Guide, 1998- 1999, Cell Culture section) followed by cenfrifugation at 1500 rpm for 5 min. The supernatant was then discarded and the cells were cultured in 100 ml media at 37 °C, 5% CO 2 . After 2 days, the cells were counted (approximately 8 x 10 5 /ml, total 100 ml). Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by DNA quantitation (see Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, one gram of preparation was added with 10 ml of water. The mixture was boiled for 30 minutes. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter. HepG2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. One to fifty dilutions were used for each assay. The cell growth inhibition assay was performed for 72 hours. All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA ( 99n 88: 8495-8499. β-glucuronidase was assayed as HQT is known for its anti-diarrhea properties. Different HQT extracts were added to triplicate wells of a 96-well plate which contained O.lmM phenolphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of dialyzed β- glucuronidase (from E. Coli, purchased from Sigma™) to a final volume of 80 μl. After 2 hr incubation, at 37°C, the reactions ere terminated with 200 μl of stopping solution which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was monitored with a kinetic microplate reader at 540 nm. , The results ofthe assays using the three batches are displayed in Table 6. Based on these data, HQT sources A and B have relatively low toxicities cpmbined with higher inhibitory activity relative to batch HQT C (i.e., approximately 5 fold greater toxicity toward HepG2 cells and 3.3 fold less inhibitory activity against β-glucuronidase than either HQT A or B, see Table 6). T able 6. Biolo2ical Assay of Three Preparations of HOT* ' E. Coli HepG2 Jurkat HBVJ β-Glucuronidase DNA HQT A ' ■ : r 0.6 . 1.50 0.76 None HQT B • ' ,ι ' °- 7 - 1.6 0.81 ND HQT C ." ■ ' 1.1 ■ - ' . 0.32 ND ND *Values t, % of 'Control represent IC 50 values. ND, not determined. Evaluation of HOT Effects on Protein E> roression HeρG2 cells (1 x 10 6 ) were seeded in 25 cm 2 flasks in 3.0 ml of RPMI-1640 medium (see Life Technologies,, Inc. ' ,; Catalpgue and Reference Guide, 1998-1999, Cell Culture section) 24 hr before the drug addition.. The cells were treated with or without herbal medicine, where the former is added at two final concentrations of 0.2 mg/ml or 4 mg/ml, respectively, and incubated at 37°C for 24 hours ' : The medium was removed and the cells were washed twice with cold PBS. The cells were harvested into 1 ml of PBS and centrifuged at 10,000 rpm for 2 minutes, extracted on ice with a buffer containing 50 mM Tris-Cl (pH 7.5), 0.2 mM PMSF and 10% glycerol, followed by three freeze-thaw cycles. Potassium chloride was added to the cell lysate at a final concentration of 0.15 M prior to centrifugation. The protein concentration was determined and the cell extract was electrophor.esed according to the method of Laemmli (Nature (1970) 227:680-685). Western blots were performed by standard techniques known in the art, see for example Sariibrook, et al (1989). The antibodies used were directed to the following proteins: Topo ϊ; Stέf (20707); Cyclin Bl; MAPK (Ab2) and Nm 23 HI. Figure 4 demonstiatέs that the higher concentrations of HQT A or HQT B differentially effects the expression of cyclin B 1 , polypeptide. HPLC Analysis ; " ' "'' ' ' ', "■ '" " The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere™ column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin Elmer). The wavelengths fof UV detection were monitored at 280 nm and 340 nm. The mobile phase was pumped _ιf 1 ml/min and consisted of Solvent A: H O and Solvent B : 20% MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 minutes. This was followeσb ' y the addition of 100 % solvent A for 5 minutes. The HPLC markers are baicalin aridbaicaleih. Mass Spectrometry The herbal extract was analyzed by Mariner™ ESI-TOF Mass Spectrometry (MS) from PE Biosystems. Control tracings were generated using baicalein and baicalin, two known active ingredients in HQT. ; ' HQT samples in water and acid treated batches were been analyzed by HPLC and Mass Spectrometry. While water treated HQT batches A and B had distinct HPLC and MS tracings, acid treated batches gave almost ideritical patterns (data not shown). Algorithm " • ' * • The data collected 'form part ofthe multidimensional analysis used to generate multivariant normal distribution sets as a means pf determining a baseline correlation between biological activity and standard HQT chemical (HPLC and Mass Spec), arid origin/growth characteristics. ' : Example 14. Individual components A. Licorice. Evaluation of Glycyrrhizae Radix (licorice) Licorice is useful for moistening the lungs and reducing coughs, helps to relax spasm and pain. The properties of the licorice batches used in this example are presented in Table 7. Table 7. Batch Properties (licorice) Biological and Enzyme Assays To assay the quality, of herbal sources, each herbal extract supernatant was assayed and the analysis was repeated three times. For a given sample to be assayed, 1 gram of herbal powder was dissolved in 10 ml of 80° C deionized water (neutral pH) in a polypropylene tube. The tube was then incubated as outlined in Table 7, then centrifuged to obtain the supernatant. Batches of licorice were tested against either HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. Cells were cultured for 24 hours as described above. Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by DNA quantitation (see Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, one gram of preparation was added with 10 ml of water. The mixture was boiled for 30 minutes. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter.2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. One to fifty dilutions were used for each' assay. The cell growth inhibition assay was performed for 72 hours. All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-84991: ; '■ ' ] ■ Again, β-glucuronidase was assayed. Different licorice extracts were added to triplicate wells of a 96-well plate which contained O.lmM phenolphthalein glucuronidate, 70 mM Tris- HCl (pH 6.8) and 0.8 ήg of dialyzed beta-glucuronidase (from E.Coli, purchased from Sigma) to a final volume of 80 μl and ' assayed as above. The results ofthe assays using the two batches is displayed in Table 8. Based on these data, licorice batch A was much more toxic to Jurkat cells than batches B (approximately 9 fold) and a more effective inhibitor of β-glucuronidase (see Table 8). Table 8. Biological Assay of Four Preparations of Licorice* • ' . " E. Coli ' ' .. HepG2 Jurkat HBVJ β-Glucuronidase DNA Licorice A • ,. ' . ' ' 1.07 0.41 None Licorice B '• : ' \; ND. ' ND 3.6 ND Licorice C . . . .. -; . 2.1 ND ND ND Licorice D • ' ' ' ' . ' •■ND ,. " ND • >2.0 53.8 *Values represent IC 50 . %, % of Control values. ND, not determined. Expression Assay In order to assay gene expression, Jurkat T cells were treated with herbal extract as follows: Jurkat cells (107ml) were quickly thawed in a water bath at 37 °C. The cells were then diluted in 10 ml of pre-warmed media (see Life Technologies, Inc., Catalogue and Reference Guide, 1998-1999, Cell Culture section) followed by centrifugation at 1500 rpm for 5 min. The supernatant was. then discarded and the cells were cultured in 100 ml media at 37 °C, 5% CO 2 . After 2 days,;the cells were counted (approximately 8 x 10 5 /ml, total 100 ml). The herbal extract solution was prepared as outlined above (e.g., 2 g of an herbal powder to obtain 20 ml of sterile solution (0.1 g/ml). The cells were divided into 3 flasks at a density of 2.5 x 10 5 /ml, 100 ml/each flask. Assays were carried out with control (no extract), and 10 ml of extract at 10 mg/ml, and 1 mg/ml. Again, toxicity results were used to determine the "high" and "low" concentrations for any given extract. After extract addition, cell cultures were incubated for 24 hours uridef conditions as outlined above. The cells were counted and subsequently collected in 50 ml centrifuge tubes. The resulting cell pellet was treated with an RNA isolation means to extract mRNA (see, for example, Sambrook et al, 1989 at pages 7.3- 7.39). ' \ " '' ; " - :J Λ ■' ■ •' . " • • • . Microarray ' , . " ' ' ' • Microarray printing was carried out as follows: Human gene clones ' were obtained from the IMAGE Consortium libraries through its distributors and comprise genes from, various tissues. Most clones have been partially sequenced and the sequences wefe ; available as expressed sequence tags in the dbEST database of GenBank. Clones were cultured and amplified using commercially available primers prior to application on nylon meriibraries : (Chen et al, Genomics (1998) 51:313-324). Approximately 10 ng of each amplified target was applied on a positively charged nylon membrane using a PC (personal computer) contf oiled arraying system. The arraying system allows high density spotting and is capable of depositing 31,000 spots on a piece of nylon membrane measuring 18 by 27 mm using a 24-pin arraying tool. cDNA probe and Membrane Hybridization Two microgram of each mRNA sample (mRNA was isolated as outlined above) was labeled with biotin and/of digόxigemn using random primed reverse transcription. The labeled samples were treated with alkali and the resulting labeled nucleic acids were precipitated prior to use m hybridizatiori. Membrane hybridization and washing were carried out using the labeled probes as disclosed in .Chen et al (1998). To detect the spots on the membrane in dual color mode (i.e., both biotin and digoxigenin), β-galactosidase-conjugated streptavidin (Strept- Gal) and alkaline phosphatase-conjugated digoxigenin antibody (anti-Dig- AP) were employed. After color development, iriiage digitization using an imaging means was employed (e.g., a flatbed scanner or digital .cariiera). Quantitative measurements were determined by computer analysis which uses a program that measures the integrated density of the primary color components of each spot, performs regression analysis ofthe integrated density data and locates statistical outliers as differentially expressed genes. Gene expression data for samples 1. 2 and licorice (ST117) Extract 1, 2 and 6 corresponding to extract of Cordyceps sinensis, Poria cocos (ST 027) and licorice, respectively, were assayed by the following method: Batches were evaluated for toxicity using Jurkat T cells. The extracts were prepared as putlined in Example 6. The cells were divided into 24 well culture plates by adding 1- ml of Jurkat cells at a density of 5 x 10 5 /ml. Assays were carried out with control (no extract), and 5 concentrations of extracts as described (see Table 9). The high and low concentrations for the cell culture assays were varied between 10 mg/ml and 0.05 mg/ml (i.e., mg dry weight of herbal extract per ml) depending on the toxicity ofthe extract to cells. For certain sarήples the toxicities at 10 mg/ml were such that "high" and "low" concentrations were adjusted downward, nevertheless, at least pne prder pf magnitude between extremes was maintained. - For example, for licorice (ST117) the "high" was 0.5 mg/ml and the "low" was 0.05 mg/ml (see Table 9). After extract addition, cell cultures were incubated for 24 hours under conditions as. outlined in Example 6. The cells were counted and the resulting data tabulated to demonstrate extract toxicity. The resulting data is shown in Table 9. Table 9. Survival Cell Number at Different Concentration of Herbal Extract Solution experiment concentration no.xl0 5 ml ■ 10 mg ml 5 mg/ml 2 mg ml 1 mg/ml 0.5 mg/ml No drug High cone Low cone, (mg/ml) (mg/ml) 1 Cordyceps sinensis 8. 4 ■ 11.9 , 11.5 9.2 9.0 12.4 10 1 mycelium 2 ST027 4.2 8 . 10 7.5 10 10 10 1 3 STO-44 ' - . •' -5.9 8.4 9.9 9.4 5 0.5 4 ST051 - .-. - • 1.7 5.4 0.5 0.05 5 ST093 - 1.9 3.8 4.4 0.5 0.05 6 ST117 i ' .β ' 3.6 4.6 5.8 0.5 0.05 7 ST123 3.4 '" 6 - 4 8 9.3 7.8 5 0.5 8 ST128 3.5 • ■ .7.7 ' • : 7.9 7.7 8.3 5 0.5 9 ST13 4 2.9' ' '■ ■ ' .6.1 , 11.2 9.6 9.8 5 0.5 10 ST237 2.5 6.6 8.7 1 0.1 Note: original c ell numb er is 5x 10 ml and the number t 0 lOxlC )7ml after 24h incubatic describes all dead cells; Protocol: ..- ,. ' • ' 1. Add 1 ml of 5xl0 5 /ml Jurkat cells into 24 well culture plates. 2. Prepare 12 kinds of herbal extract solutions and sterilize. 3. Test 5 concentratipn per sample. 10 mg/ml, 5 mg/ml, 2 mg/ml, 1 mg/ml, 0.5 mg/ml 4. Culture the cells fpr 24 h in 37C with 5% CO incubator. In each analysis, 144X96 genes (i.e., 13,824 genes) were analyzed (data not shown) and about 100 genes showed significant differences in comparison with that of control (Table 10). Some of the. genes were' up-regulated and others were down-regulated. The magnitude of the difference with the control sometirnes varied depending on the relative amount ofthe herbal composition to which the particular cells were exposed. Numbers under Cl (control treatment) and H or L (herbs) represent intensities of mRNA expressed after subtraction of background (Table 10). The gene designation is encoded in Array AD, which can be traced to a specific GenBank clone.' The level of expression was determined by H or L divided by C. Only a fraction of 13,824 genes in each herb treated samples showed significant changes, namely, up, down or unchanged (see" Table 10). rem The value show in colume K,L,M,N,O,P represent the ratio of mRNA express lev Title Cl IH ArrayAD TRANSDUCIN-LIKE ENHANCER PROTEIN 1 Sg-Bk : Sg-Bk <94,030> Adenylosuccinate lyase 55.29 8.55 <91,097> Neurotrophic tyrosine kinase, receptor, type 1 89.55 98.29 <89,083> MHC class I protein HLA-A (HLA-A28,-B40, -Cw3) 80.68 61.91 <87,083> 85.33 50.96 <86, 131> ESTs, Highly similar to HELIX-LOOP-HELIX PROTEIN 99.48 109.46 <85, 129> INTEGRAL MEMBRANE PROTEIN El 6 90.64 72.12 <85,1 12> 35.23 60.82 <84,075> 89.56 77.54 c <83,129> 118.83 108.15 m <82,082> Homo sapiens androgen receptor associated protein 24 (Aϊ ' 60.23 47.87 <82,027> 65.11 90.01 <80,129> 96.27 99.96 <80,035> Human small GTP binding protein Rab7 mRNA, complete 44.91 46.23 <74,108> 65.2 77.25 <74,012> 34.68 60.92 <73, 132> Human mRNA for KIAA0078 gene, complete eds 109.86 103.38 <72,101> Human trans-Golgi p230 mRNA, complete eds 46.96 25.76 <72,014> Chlordecone reductase 4.37 15.08 <70,101> 60.86 13.55 <69,107> 37.92 31.41 <68,064> Human mRNA for KIAA0034 gene, complete eds 2.24 1.27 90 © r- © ( H, L for high ccoonnee.. A And low cone) compared with untreated cell (Cl) © CΛ 6H 6L IH IL 2h 2L 6H 6L P Sg-Bk Sg-Bk COMPARE TO CONTROL ClonelD H U a- 29.21 30.82 0.154639 0.129318 0.845361 0.688551 0.528305 0.557424 504540 36.02 89.63 1.097599 0.78459 0.947292 1.047795 0.402233 1.000893 1013392 46.16 32.35 0.767353 0.571765 0.741943 0.367377 0.572137 0.400967 510048 37.5 49.82 0.597211 0.4424 0.624048 0.299309 0.43947 0.583851 1188706 35.65 75.78 1.100322 1.061721 0.569763 0.67893 0.358363 0.761761 72215 48.47 59.36 0.795675 1.023279 0.41891 0.773169 0.534753 0.654898 118548 28.25 59.33 1.72637 1.38717 1.250071 0.843315 0.801873 1.684076 115134 43.37 42.97 0.865788 0.869808 0.738053 0.663689 0.484256 0.47979 172767 55.8 66.27 0.910124 0.993436 0.638223 0.539342 0.469578 0.557687 1184183 44.8 46.5 0.794787 0.580442 0.621119 0.163374 0.743815 0.772041 118235 oo yr, 39.3 78.37 1.38243 1.025035 1.018277 0.98971 0.603594 1.203655 171864 48.32 53.87 1.03833 1.135348 0.850525 0.910564 0.501922 0.559572 1172135 37.48 42.49 1.029392 0.5609 1.741483 0.90737 0.834558 0.946114 1172266 51.52 70.19 1.184816 1.044325 0.923313 1.015491 0.790184 1.076534 37502 42.34 65.02 1.756632 1.664648 0.709631 1.4109 1.220877 1.874856 45672 61.91 32.44 0.941016 0.74959 0.782541 0.629437 0.563535 0.295285 79342 1.27 17.51 0.548552 0,390758 0.224446 0.380537 0.027044 0.372871 563992 20.64 30.11 3.450801 3.045767 9.810069 5.208238 4.723112 6.89016 21759 16.36 18.7 0.222642 0.57887 0.548636 0.667926 0.268814 0.307263 564469 36.95 50.55 0.828323 0.688819 2.184072 0.755274 0.97442 1.33307 29009 © 4.2 41.64 0.566964 2.370536 0.09375 1.116071 1.875 18.58929 51927 90 <62,117> 06.71 94.67 83.98 81.39 61.45 <62,084> NA 1.13 0 0 13.55 0.58 <62,001 > Heterogeneous nuclear ribonucleoprotein A 1 97.11 71.92 73.16 42.01 78.94 <61,118> NA 89.52 59.97 54.11 66.86 43.05 <59,142> ATPase, Na+/K+ transporting, alpha 1 polypeptide 13.03 15.48 47.75 27.71 27.86 <59,135> 60.85 90.23 39.83 93.01 38.17 <58,087> ESTs, Weakly similar to KIAA0062 [H.sapiens] 77.72 46.86 76.33 64.43 22.52 <57,016> Neuroblastoma RAS viral (v-ras) oncogene homolog 71.71 16.47 44.61 39.51 66.89 <56,088> 3-HYDROXY-3-METHYLGLUTARYL-COENZYME A 56.77 52.66 59.53 59.6 16.89 <56,060> Translation elongation factor 1 -alpha- 1 81.26 92.28 92.08 36.86 79.09 <55,040> 0 43.58 39.47 53.67 46.42 <54, 134> CARBONIC ANHYDRASE III 0 0 4.38 43.47 8.23 in <53,048> H.sapiens mRNA for interferon regulatory factor 3 0.25 0 0 33.47 0 <53,023> Human mRNA for proteasome subunit HsC7-I, complete c 67.76 14.98 42.33 38.48 66.52 <52,135> 49.26 46.81 44.4 72.95 55.5 <51,131> ATP citrate Iyase 2.02 0 52.06 39.33 6.03 <51,101> POLYADENYLATE-BINDING PROTEIN 8.65 0.17 0.45 38.99 14.31 <50,035> SERUM ALBUMIN PRECURSOR 0.88 0 10.21 44.14 0 <47,055> ESTs 38.53 31.69 19.45 43.8 56.1 <46,097> 34.26 35.95 32.49 39.81 63.19 <45 , 100> Ribosomal protein L5 26.15 33.6 26.3 55.86 48.59 <44,103> NA 72.01 54.36 70.27 75.43 91.53 <43,039> 55.86 18.02 20.83 11.68 17.8 <43,035> 14-3-3 PROTEIN TAU 86.86 23.18 36.23 18.14 12.76 <43,019> 72.48 10.4 23.38 27.59 12.83 90 © yo 77 41.8 0.887171 0.786993 0.762721 0.57586 0.721582 0.391716 241351 © © 32.81 17.71 0 0 11.99115 0.513274 29.0354 15.67257 28012 CΛ 83.94 96.87 0.740603 0.753372 0.432602 0.812893 0.864381 0.997529 236388 H U 65.65 34.92 0.669906 0.604446 0.746872 0.480898 0.733356 0.39008 240018 α. 23.22 16.58 1.188028 3.66462 2.126631 2.138143 1.782041 1.272448 121270 86.39 31.2 1.482827 0.65456 1.528513 0.62728 1.419721 0.512736 510863 56.78 63.4 0.602934 0.982115 0.829002 0.289758 0.730571 0.815749 645239 50.43 35.52 0.229675 0.622089 0.550969 0.932785 0.703249 0.495328 47526 72.87 60.28 0.927603 1.048617 1.04985 0.297516 1.2836 1.061828 509520 71.74 70.47 1.135614 1.133153 0.453606 0.973296 0.882845 0.867216 31027 40 65.93 328351 32.61 9.74 287006 0 0.53 0 0 133.88 0 0 2.12 116915 27.24 37.24 0.221074 0.624705 0.567887 0.9817 0.402007 0.549587 221285 63.21 30.35 0.950264 0.90134 1.480918 1.126675 1.283191 0.616119 286222 22.36 5.6 0 25.77228 19.4703 2.985149 11.06931 2.772277 624420 38.25 11.64 0.019653 0.052023 4.507514 1.654335 4.421965 1.345665 529138 0.79 0.01 0 11.60227 50.15909 0 0.897727 0.011364 510245 59.36 53.2 0.822476 0.504801 1.136777 1.456008 1.540618 1.380742 26801 30.86 51.76 1.049329 0.948336 1.161996 1.844425 0.900759 1.5108 511591 31.26 43.29 1.284895 1.005736 2.136138 1.858126 1.195411 1.655449 511410 70.5 80.86 0.754895 0.975837 1.047493 1.271073 0.979031 1.1229 26099 11.82 41.06 0.322592 0.372897 0.209094 0.318654 0.2116 0.735052 592919 12.99 42.28 0.266866 0.417108 0.208842 0.146903 0.149551 0.48676 120011 3.47 15.99 0.143488 0.322572 0.380657 0.177014 0.047875 0.220613 488154 30.25 33.73 36.83 61 51.56 ESTs 0 0 0.2 38.37 12.78 Thioredoxin 4.26 0.34 2.83 37.28 8.03 Pre-alpha (globulin) inhibitor, H3 polypeptide 9.46 24.06 15.82 26.08 25.48 <38,097> NA 11.81 6.94 23.82 80.14 16.65 <38,081> Human GP36b glycoprotein mRNA, complete eds 101.51 82.61 89.86 22.88 71.4 <35,096> 47.35 24.52 0 23.2 32.88 <34,143> 61.23 2.99 25.05 30.61 16.93 <32,118> ESTs, Weakly similar to CASEIN KINASE I HOMOLOG 21.98 8.91 14.11 53.54 24.01 <32,034> Pyruvate kinase, muscle 78.41 54.15 79.36 29.56 40.5 <30,102> RETINOBLASTOMA BINDING PROTEIN P48 27.3 23.1 16.06 62.35 34.39 <30,033> 70.49 30.32 65.13 17.18 33.65 <29,058> NA 25.77 16.66 0 59.85 33.24 <29,035> 60.19 57.12 86.66 46.47 56.59 <28,104> Ribosomal protein S13 0 0 0 37.19 1.37 <28,035> 71.46 35.93 51.86 11.27 30.75 <27,119> 0 0 0 48.27 1.66 <27,097> ESTs, Moderately similar to Etr-3 [H.sapiens] 16.06 24.53 9.26 76.01 24.34 <25,126> ESTs 42.93 22.03 29.88 67.76 40.39 <25,106> 61.9 48.38 43.41 91.26 58.89 <24,098> Human splicesomal protein (SAP 61) mRNA, complete cd 24.26 25.43 21.8 56.18 37.28 <24,071> HEAT SHOCK 70 KD PROTEIN 1 63.92 58.19 61.1 50.86 66.27 <23,126> 5.18 1.1 10.08 54.71 12.11 © 90 <23,099> yo Ribosomal protein S3A 5.13 1.93 1.25 38.99 14.35 yo <22,013> EUKARYOTIC INITIATION FACTOR 4B 61.8 49.95 41.52 22.99 49.11 <22,006> 3-HYDROXY-3-METHYLGLUTARYL-COENZYME A 83.17 57.84 39.8 40.04 44.33 <20,126> MITOCHONDRIAL STRESS-70 PROTEIN PRECURSO: 27.02 22.19 33.86 56.87 32.97 <20, 104> Human eukaryotic translation initiation factor (eIF3) mRN 22.09 15.17 17.33 74.98 19.17 <20, 103> Laminin receptor (2H5 epitope) 1.45 1.26 0 39.66 8.67 <19,098> Human semaphorin (CD 100) mRNA, complete eds 17.79 14.38 10.79 52.94 28.97 <18, 103> Ribosomal protein L5 31.47 21.13 30.07 78.6 27.18 <18,084> ESTs 54.61 75.17 46.16 79.5 60.01 <17,103> NA 0 0 0 33.96 2.87 <17,097> NA 1.52 2.1 1.21 37.77 11.59 <17,062> 92.72 73.27 103.64 69.79 69.35 <16,103> Integrin, beta 1 (fibronectin receptor, beta polypeptide, an 0 0.92 0 41.3 1.84 < 16, 100> Homo sapiens E 1 B 19K/Bcl-2-binding protein Nip3 mRN, 0.79 0 0.7 38.06 8.28 <15,1 12> 16.47 11.65 11.87 46.52 22.98 £ <14,103> NA 27.61 34.26 30.18 55.67 38.45 <14,011> NA 47.42 34.87 68.28 38.71 53.6 <13, 1 14> Homo sapiens clone 24689 mRNA sequence 20.64 13.87 19.2 51.26 32.34 <12, 144> Homo sapiens NADH:ubiquinone oxidoreductase NDUFS 76.78 67.2 76.11 64.5 50.04 <12,111> 40.04 55.04 46.34 77.78 46.07 <11,135> 90.34 114.29 87.88 80.49 87.22 < 11 , 105> Ribosomal protein L3 8.27 14.41 15.59 36.82 25.87 <1 1 ,053> Calpain, large polypeptide L2 10.48 23.58 14.8 15.69 8.87 <10, 136> Malic enzyme 2, mitochondrial 74.34 81.63 90.94 83.18 74.85 <10,101> ESTs 13.65 16.34 13.97 47.57 26.68 2 <09,136> NA 59.82 95 87.16 77.34 83.51 I <09,110> ER LUMEN PROTEIN RETAINING RECEPTOR 1 17.3 41.73 48.92 55.32 35.56 -59/5- O © b U ) to — — © © _- o o © © © © © © 00 4s. © NJ 4s. bo o ^_ SO SO 4s. o o bo l—L ^ bo SO ft ^ to 4s. to Ul to to 00 UJ so Ul SO UJ UJ so to 00 SO 4s. Ul © 00 Ul -O -J UJ t Ul ■ o -J © so UJ © ft ■ o ft Ul Ui ft 4s. UJ 00 -J © 4s. t to to UJ to to 00 © -J o -o © to Ui Ul U Ul 4s. UJ 4s. so to -J 4s. -o UJ so 00 © -J -J Ul ft NJ 4s. UJ ft 00 UJ 00 4s. -o -J UJ 4s. UJ SO to 4s. ft 4s. UJ ft UJ © 00 ft 00 UJ > _ UJ _ 4s. © © t © t to 4s. © to NJ NJ NJ UJ NJ © to 4s- i_t 4s. 4s. bo SO bo 4s. 00 © bo 00 ^ 4s> 4s. 4s. SO -4 UJ 4s. so so 00 »__ SO Ui so 4s. 4s. 00 t >—_ Ul bo Ul o -o UJ so © 00 •o to 4s. 00 •o t © to © J ft ft 4s. ~J to 4s. Ul -J Ul Ul 4s. 4s. (—- ft 00 SO so to o Ui o Ui UJ to Ul -J ft 00 l ft 00 >--. NJ -J 4s. 00 -J 00 UJ UJ ft Ul ft to to so to t so ft ■ o »—- NJ -o SO UJ NJ 00 so IO U ) ^J ^J • ^1 »o UJ •o to 00 O to ft 00 ^ so NJ ft ■ o 4s. to © UJ © H—. © >—* —. _ © _- © © © UJ SO 1—4 bo ^_ SO ft Ϊ -* UJ UJ © l © © ft Ul bo to Ui ui so Ul © 4s. to ft Ul Ul Ui UJ so SO 4s. 4s. NJ so ft NJ UJ Ul t 4s. © ft 00 Ul © >—- ft o t Ul 00 -o -J SO bo f oo ft 00 ■ o -o © UJ 4s. © u< ft UJ »—-. 4s. Ul ~o ft UJ ft to SO ft 00 UJ 4s. SO 00 NJ © yo t - oo -o •o ft SO UJ 00 to ft © Ul K) 00 ft 4s. UJ © © *- SO ft 4s. 4s. 4s. SO NJ ft Ul → 4s. " Ul UJ 00 UJ 1-1 UJ -o Ul © © >- * © Ul UJ © Ul © ~o o Ul UJ Ul 00 -J 4s. ~J >—. © 4s. UJ 4s. Ul 00 ft NJ NJ UJ t © to Ul © © ,— >—» to © so ft H-* © 4s. so bo 1 ft SO UJ © UJ ^_> bo !—- bo 4s. © UJ ^J to ft to 4s. ft SO 4s. ~j Ul y o ft SO SO 4s. © to UJ © so 4s. UJ -0 Ul NJ NJ Ul NJ SO 00 © so 00 Ul ft © © to ft ft Ul ft © UJ Ul ft NJ UJ ft SO to ft 00 Ul to OO >— * © ft Ul »— * © ft -o UJ Ul ft Ul 4s. 1—* -O 4s. ft Ul -o 00 00 to U to Ul 4s. 00 ft -J 00 NJ so © 00 Ul 4s. to *— * so ft 00 4s. 00 so 00 4s. Ul to -o UJ 00 so -J ft Ul ft 4s. ft ft Ui Ui UJ Ui ft ft OO ft NJ NJ NJ UJ SO 3 -O 00 NJ SO ft ft 00 NJ ft ft o <09,085> Human phospholipid transfer protein mRNA, complete eds 80.14 88.83 96.22 80.4 77.44 I <08,138> 41.57 60.4 69.53 73.27 81.26 <08,132> 67.22 102.84 103.65 90.56 88.43 H <08, 122> Homo sapiens mRNA for DRAK 1 , complete eds 32.01 23;21 47.74 63.02 35.64 * <08,105> 69.58 58.7 84.49 90.5 68.91 <08,055> 19.5 20.01 22.83 27.87 34.66 <07,1 11> 30.8 44.67 48.75 73.84 51.31 <07,107> ESTs 4.57 3.23 9.01 41.32 8.83 <06,135> 71.42 112.68 128.81 102.8 101.48 <06,102> ESTs 39.82 36.03 32.57 67.74 46 _..0_4. <06,101> Intercellular adhesion molecule 2 4 4..4433 2 2..9955 4 4..2211 3 388..9955 1 166..66 <06,034> Ribosomal protein L5 0 0 0 0 33.93 <06,006> Human fragile X mental retardation syndrome related protø 0 31.71 0 0 0 <03,099> ESTs 22.97 18.57 12.1 52.83 23.48 <02,111> ESTs, Highly similar to HYPOTHETICAL 64.5 KD PRO 34.89 23.63 25 77.41 28.84 <01 , 136> Homo sapiens mRNA for 5-aminoimidazole-4-carboxamic 21.36 34.84 23.95 53.14 45.66 <01 ,058> Acyl-Coenzyme A dehydrogenase, C-4 to C-12 straight chi 60.68 53.35 73.26 62.56 56.58 Total 116 90 © yo 76.05 82.81 1.108435 1.200649 1.003244 0.966309 0.948964 1.033317 613511 © 59.49 66.73 1.452971 1.6726 1.762569 1.954775 1.43108 1.605244 632480 © CΛ 76.85 81.15 1.529902 1.541952 1.347218 1.315531 1.143261 1.20723 66427 H 32.48 46.85 0.725086 1.491409 1.96876 1.113402 1.014683 1.463605 308580 Uα. 89.98 94.96 0.843633 1.214286 1.300661 0.990371 1.293188 1.36476 30607 48.86 50.43 1.026154 1.170769 1.429231 1.777436 2.505641 2.586154 611114 59.76 68.29 1.450325 1.582792 2.397403 1.665909 1.94026 2.217208 46880 22.9 23.56 0.706783 1.971554 9.041575 1.932166 5.010941 5.155361 29970 93.72 94.2 1.577709 1.803556 1.439373 1.420891 1.312237 1.318958 286790 59.07 57.13 0.904822 0.817931 1.701155 1.156203 1.483425 1.434706 530700 15.85 32.59 0.665914 0.950339 8.792325 3.747178 3.577878 7.356659 530665 0 0 666094 0 0 530672 64.83 64.43 1.399436 1.19421 1.839611 1.507046 1.66103 1.650781 285979 56.63 55.12 1.512706 1.393901 2.969504 1.393901 3.59784 3.501906 31309 ft 74.42 81.16 1.575071 1.092903 1.660128 1.320257 1.327033 1.447218 624595 53.26 41.55 0.808446 0.526774 2.299956 1.022203 2.318677 1.808881 42281 56.6 60.36 0.677271 0.716538 2.218687 0.826598 1.622241 1.730009 45327 40.67 51.53 1.631086 1.121255 2.487828 2.13764 1.904026 2.412453 567287 72.81 91.71 0.879202 1.207317 1.030982 0.932432 1.199901 1.511371 30921 -60- In this manner, we are able to correlate specific gene expression with the exposure of a cell to no, low (L) or high (H) amounts of an herbal composition. Many ofthe genes identified in this way code for proteins important in known metabolic or biochemical pathways. Many of these proteins have direct and indirect effects on certain physiological, morphological and psychological parameters. Thus, this method permits the association of a particular genetic fingerprint of an herbal composition with its array biological effects. Such associations can be used to profile or characterize an herbal composition for the purposes of Quality Control and Quality Assurance and evaluating pharmacological or toxicological properties. The role of primary and secondary herbs in an herbal formula can also be assessed by this approach. HPLC Analysis The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere™ column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The mobile phase was pumped at 1 ml/min and consisted of Solvent A: H 2 O and Solvent B: 20% MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 minutes. This was followed by the addition of 100 % solvent A for 5 minutes. The HPLC marker is glycyrrhizin. ' ' ' ' ' " Algorithm The data collected form part ofthe multidimensional analysis used to generate multivariant normal distribution sets as a means of determining a baseline correlation between biological activity and standard licorice molecular, chemical (HPLC and Mass Spec), and origin/growth characteristics. B. Scute Evaluation of Radix Scutellariae (Scute) Scute has been found to be useful in reducing capillary permeability and inflammation. It can also be used treat enteritis and dysentery, increases the secretion of bile to treat jaundice; to relieve muscle spasms; to treat coughing and to expel parasites. The properties ofthe scute batches used in this example are presented in Table 11. -61- Table 11. Batch Properties (Scute) Biological and Enzyme Assays Briefly, one gram of each preparation of scute extract was added with 10 ml of water (1 mg/ml). The mixture was treated as Outlined in Table 11. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter. Batches of scute were tested against either HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. One to fifty dilutions were used for each assay. Cells were cultured for 24 hours as described above. Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by DNA quantitation (see Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, one gram of preparation was added with 10 ml of water. The mixture was treated as outlined in Table 11. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter.2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad S ' ci USA(1987) 84: 1005-1009) were used in this assay. One to fifty dilutions were used for each assay. The cell growth inhibition assay was performed for 72 hours. All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499. : ' For β-glucuronidase, different scute extracts were added to triplicate wells of a 96-well plate which contained O.imM phenόlphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of dialyzed β-glucuroriidase (from E. Coli, purchased from Sigma) to a final volume of 80 μl. After 2 hr incubation at 37°C, the reactions were terminated with 200 μl of stopping solution which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was monitored with a kinetic microplate reader at 540 nm. -62- The results ofthe assays using the three batches is displayed in Table 12. Table 12. Biological Assay of Four Preparations of Scute* Ε. Coli HepG2 Jurkat HBVJ DNA β-Glμcuronidase Scute A , ...I.5. - 0.33 0.45 None Scute B ' ■* . ■ ' _ ■"• _ : - ' ' : ' ' ND ND ND Scute C ' ' , ;' "θ.3 ' ' ND ND ND Scute D ND 0.65 ND 27.5 *Values represent %, % of Control IC 5 o values. ND, not determined. Evaluation of Scute Effects oh Protein Exi pression HepG2 cells (1 10?) were seeded in 25 cm 2 flasks in 3.0 ml of RPMI-1640 medium (see Life Technologies, Inc.j Catalogue and Reference Guide, 1998-1999, Cell Culture section) 24 hr before the extract addition. The cells were treated with or without herbal medicine, where the former is added at two final concentrations of 0.2 mg/ml or 4 mg/ml, respectively, and incubated at 37°C for 24 hours. The medium was removed and the cells were washed twice with cold PBS. The cells were harvested into 1 ml of PBS and centrifuged at 10,000 rpm for 2 minutes, extracted on ice with a buffer containing 50 mM Tris-Cl (pH 7.5), 0.2 mM PMSF and 10% glycerol, followed by three freeze-thaw cycles. Potassium chloride was added to the cell lysate at a final concentration of 0.15 M prior to centrifugation. The protein concentration was determined arid the cell extract was electrophoresed according to the method of Laemmli U.K. (Nature (1970) 227:680-685). Western blots were performed by standard techniques known in the art, : seie for example Sambrook, et al (1989). The antibodies used were directed to the following proteins: Topo I; Stat (20707); Cyclin Bl; MAPK (Ab2) and Nm 23 Hl. Figure 4 demonstrates that scute batches A and B do not differentially affect the expression ofthe polypeptides resolved on Western blots. HPLC Analysis ■ ; ■ • • ' " The herbal batches wer analyzed by HPLC with a Beckman ODS Ultrasphere™ column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The mobile phase was pumped atT πil/min and consisted of Solvent A: H 2 O and Solvent B: 20% -63- MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; and 3) the solvent was changed to 10% solvent A / 90% solvent B for the next 40 minutes. This was followed by the addition of 100 % solvent A for 5 minutes. The HPLC markers are baicalin and baicalein. . Scute batches in wafer and acid treated samples were analyzed by HPLC. Water and acid treated batches were virtually indistinguishable. Algorithm The data collected form part ofthe multidimensional analysis used to generate multivariant normal distribution sets as a means of determining a baseline correlation between biological activity and standard scute chemical (HPLC), and origin/growth characteristics. C. White Peony Root " . Evaluation of Paeonie lactiflora pallus radix (Peony) Peony is used to suppress and soothe pain. It is also known to soothe ligaments and purify the blood. The properties ofthe peony batches used in this example are presented in Table 13. Table 13. Batch Properties (Peony) Biological and Enzyme Assays Briefly, one gram of each preparation of scute extract was added with 10 ml of water (1 mg/ml). The mixture was treated as outlined. in Table 13. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter. Batches of peony were tested against either HepG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #ΗB-152) or both. One to fifty dilutions were used for each assay. Cells were cultured for 24 hours as described above. Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by DNA quantitation (see Don et al:, Proc Natl Acad Sci USA (1991) 88: 8495-8499). -64- Briefly, one gram of prepafation was added with 10 ml of water. The mixture was treated as outlined in Table 13. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter ' . . 2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. One to fifty dilutions were used for each assay. The cell growth inhibition assay was performed for 72 hours: All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA ( J991 . 88: 8495-8499. Different peony extracts were added to triplicate wells of a 96-well plate which contained O.lmM phenolphthalein glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of dialyzed beta-glucuronidase (from E.coli, purchased from Sigma) to a final volume of 80 μl. After 2 hr incubation at 37°C,-the reactions were terminated with 200 μl of stopping solution which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was monitored with a kinetic microplate reader at: 540 nm. Results are shown in Table 14. Table 14. Biological Assay of Two Preparations of Peony* " ; : E. .Cpli , HepG2 Jurkat HBVJ β-Glucuronidase Peony A ' . _ ' ■ 2.8 >1.5 1.1 None Peony B >2.5 ND ND ND *Values represent IC 50 J, % of Control values. ND, not determined. HPLC Analysis • • ■ •■ . The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin Elmer). The wavelengths for UV detection were monitored at 280 nm and 340 nm. The mobile phase was pumped at 1 ml/riiin arid consisted of Solvent A: H_O and Solvent B: 20% MeOH with the following gradient: l).the solvent was 100% solvent A for the first 5 minutes; 2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; and 3) the solvent was changed to 1 % solvent A / 90% solvent B for the next 40 minutes. This was followed by the addition of 100 % solvent A for 5 minutes. HPLC marker is paeoniflorin. Peony batches were analyzed * by HPLC as shown in Figure 5. Algorithm -65- The data collected form part of the multidimensional analysis used to generate multivariant normal distribution sets as a means of determimng a baseline correlation between biological activity and standard peony chemical (HPLC), and origin/growth characteristics. P. Date Evaluation of Ziziphi Fructus (Date) Date has been used for diuretic properties and strengthening effects. The properties of the date batches used in this example are presented in Table 15. Table 15. Batch Properties (Date) Biological and Enzyme Assays Briefly, one gram of each batch of scute extract was added with 10 ml of water (1 mg/ml). The mixture was treated as. outlined in Table 15. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter. Batches of date were tested against either HeρG2 cells (ATCC cat # HB-8065) or Jurkat T cells (ATCC cat #TIB-152) or both. One to fifty dilutions were used for each, assay. Cells were cultured for 24 hours as described above. Batches were also evaluated for the ability to inhibit hepatitis B virus as detected by DNA quantitation (see Dong etal, Proc Natl Acad Sci USA (1991) 88: 8495-8499). Briefly, one gram of preparation was added with lO ml of water. The mixture was treated as outlined in Table 15. The supernatant was collected after centrifugation and filtered through a 0.22 μm filter. HepG2.2.15 cells which secrete hepatitis B virons (kindly provided by Professor G. Ace; see Ace et al. Proc Natl Acad Sci USA (1987) 84: 1005-1009) were used in this assay. One to fifty dilutions were usect for each assay. The cell growth inhibition assay was -66- performed for 72 hours.' All other procedures were performed as described by Dong et al, Proc Natl Acad Sci USA (1991) 88: 8495-8499. Different peony extracts were added to triplicate wells of a 96-well plate which contained O.lmM phenolphthalein.glucuronidate, 70 mM Tris-HCl (pH 6.8) and 0.8 ng of dialyzed beta-glucuronidase (from E. Coli, purchased from Sigma) to a final volume of 80 μl. After 2 hr incubation at 37°G, the reactions were terminated with 200 μl of stopping solution which contained 0.2 M Glycine and 0.2 M NaCl (pH 10.4), and the OD was momtored with a kinetic microplate reader at 540 hm. Results are shown in Table 16. Table 16. Biological Assay of Three Preparations of Date* Εi Coli HepG2 Jurkat HBVJ DNA β-Glucuronidase •' Date A 1,2 ' . ' 1.5 -, 5.1 None Date B • ND >2.0 ND 52.3 Date C ;; 2.5 ■ ND ND ND *Values represent %, % of Control . IC 5 o values. ND, not determined. HPLC Anal vsis ' ' • • " - ■ .. - ' The herbal batches were analyzed by HPLC with a Beckman ODS Ultrasphere column (5 micron particles, 4.6 mm X 25 cm) and detected with an UV spectrophotometer (Perkin Elmer). The wavelengths fo UV detection were monitored at 280 nm and 340 nm. The mobile phase was pumped at 1 ml/ min and consist of Solvent A: H 2 O and Solvent B: 20% MeOH with the following gradient: 1) the solvent was 100% solvent A for the first 5 minutes; 2) the solvent composition was changed to 10% solvent A / 90% solvent B for the next 10 minutes; and 3) the solvent was changed to 10% solvent A / 90%) solvent B for the next 40 minutes. This was followed by the addition of 100 % solvent A for 5 minutes. HPLC markers for date are chelidonic acid and cAMP. Date batches samples were analyzed by HPLC as shown in Figure 6. Algorithm The data collected form part ofthe multidimensional analysis used to generate multivariant normal distribution sets as a means of determining a baseline correlation between biological activity and standard peony chemical (HPLC), and origin/growth characteristics. -67- Example 15. Characterizing Herbal Medicines by Nucleic Acid Microarray Analysis. Introduction. The rapid developrrient of nucleic acid microarray technology has led to an explosion of gene expression data (Lander, 1999, Duggan et al, 1999). Four characteristics ofthe gene expression account for the great value of using nucleic acid microarrays to study the gene expression profiles, (i) Nucleic acid microarray makes it easier to measure the transcripts of thousands of genes at once, (ii) Close association between the function of a gene product and its expression pattern makes gene function predictable, (iii) Cells respond to the micro- environmental changes by changing the expression level of specific genes, (iv) The sets of genes expressed in a cell determine, what the cell is derived of, what biochemical and regulatory systems are involved, and so on (Brown and Botstein, 1999). By using a microarray system, the above features cah be studied in an ensemble manner. The expression of any. desired number of genes can be detected using the nucleic acid microarray technology. For; example, up to about 20,000 genes may be placed on a single array. We have developed a nucleic acid microarray with colorimetric detection system (Microarray/CD) (Chen et al, 1998). Gene expression profiles of different cell lines were studied using microarray filter membranes (2.7 cm x 1.8 cm) with about 10,000 cDNA representing approximately 10,000 distinct human transcripts. The sensitivity and detection limits ofthe microarray/CD. system , have been characterized and are comparable to the system with radioactive detection or ' file system with laser induced fluorescence detection (Bertucci et al, 1999). « ' , • ' •: '■ ■• As previously described, cellular gene expression profiles portray the origin, the present differentiation ofthe cell, and the cellular responses to external stimulants. In other words, the gene expression profiles reveal the. state ofthe cell and microarray is a perfect tool for the rendering purpose. In the pfeseήt studies, we apply the microarray/CD system to characterize cellular responses to external stimulants, in this case, the Chinese herbal medicines. Conversely, we also based on the stimulated gene expression profiles to classify different herbal medicines: " "' Figure 7 is a flowchart depicting a general method that may be used for establishing an expression response data set for cells treated with an herbal composition. The method comprises the steps of: -68- (a) Determine the IC 50 concentration ofian herbal composition by incubating various concentrations of the herbal medicine in mammalian cell cultures and identify the concentration that leaves 50% of survival cells after a predetermined time. (b) Incubate the mariimalian cell cultures with herbal extracts of various fractions of IC 50 concentrations. (c) Harvest and count the cultured cells after a predetermined culture time. (d) Immediately lyse the cells after they are removed from the incubator and extract mRNA from cell lysate. (e) Label the mRNA by reverse transcription reaction to turn mRNA into labeled cDNA. (f) Mix the labeled cDNA with control cDNA of plant origin and perform hybridization to a microarray of mammalian gene probes, (g) Measure expression level of genes by analyzing digitized images ofthe microarray hybridization results, (h) Perform data pre-processing to select data for statistical analysis. (i) Acquired expression data generated by microarray experiments of an herbal composition with various concentrations, (j) Data pre-processing to select the genes with statistical significance in cells treated with different concentrations ofthe herbal medicine, (k) Categorize expression profiles into clusters by statistical methods such as the self- organizing-map algorithm. (1) Deduce the characteristic expression profiles for the herbal medicine based on the expression profile clusters. Figure 8 is a flowchart defnoristratifig how data sets of expression data for various batches of the herbal composition are integrated to make an expression profile database for the particular herbal composition. The expression profile database then becomes part of the HBR Array. HBR Arrays confaihiϋg expfession profiles may also be used to identify an unknown herbal composition. Figure 9 is a flowchart depicting a general method for identifying an unknown herbal composition, the method comprising the steps of: (a) Construct an HBR Array containing characteristic expression profiles for an herbal medicine or a collection of expression profiles of various herbal medicines by the aforementioned steps. (b) Obtain the characteristic expression profile data set ofthe unknown herbal composition. -69- (c) Compare the HB__ Array containing the characteristic expression profile induced by the said unknown herbal composition with a standardized HBR Array containing expression data by algorithms such as the Hamming distance algorithm. (d) Score possible alignments to identify the most probable herbal composition whose characteristic expression profiles are archived in the said HBR Array. Scoring possible alignments of HBR Arrays containing expression profiles may be performed using hierarchical cluster analysis ofthe Hamming distance matrix. Use of hierarchical cluster analysis for the Hamming distance matrix is well known in the art. The gene expression profiles may also be incorporated into the standardized HBR Array. As has been already discussed, the standardized HBR Array containing such gene expression profiles induced by an herbal composition can be used for studying the pharmacological mechanisms ofthe herbal composition, for discovering new application ofthe herbal composition, and for designing optimized formulation of a complex herbal preparation. As can be seen from the flowchart of Figure 10, the method may be generally outlined as comprising the steps of: (a) Construct a data set containing the characteristic gene expression profiles for an herbal composition. • '' (b) Score each gene by the consistency of its expression profiles in the data set using known statistical parameters, such as the coefficient of variation. (c) Based on the statistical scoring, gene expression profiles for an herbal medicine are selected to be incorporated into the standardized HBR Array. HBR Arrays containing gene expression profiles may also be used to identify signature gene expression profiles induce by. " individual chemical constituents in an herbal composition consisting of complex chemical constituents, as outlined in the flowchart of Figure 11. The method comprises the steps of: (a) Construct a HBR Array- containing characteristic gene expression profiles for an herbal composition by the aforementioned steps. (b) Determine the composition of chemical constituents in an herbal medicine by high performance liquid chromatography (HPLC) or liquid chromatography mass spectrometry (LC- ASS). (c) Repeat the step (b) for various batches of herbal medicine preparations. (d) Score the correlation ' coefficients between the expression levels of each gene with the amount of individual chemical constituent in an herbal preparation. -70- (e) The signature gene expression profiles for individual chemical constituent are selected with a Pearson correlation coefficient exceeding 0.99 or smaller than -0.99. Any herbal composition can then be characterized through the use of gene expression profiles generated through the use of nucleic acid microarrays. Moreover, one can choose any number of genes that are differentially expressed to be included in the data set represented the gene expression profiles. For example, one may choose about 10 genes, about 100 genes, about 500 genes, about 1000 genes, about 1500 genes, abut 2000 genes, about 2500 genes or more, or any number in between. The prescription ofthe Chinese herbal medicine Scute and Licorice combination (Huang Chin Tang) stops diarrhea, relieves spasms and clears fever. The ingredients of Huang Chin Tang are Scute, Peony, Licorice and Jujube. This recipe has been used for more than 1000 years but the chemical and biomedical studies on the prescription have not been carried out until recent decades. In this study we used the nucleic acid microarray technology to study the gene expression profiles oϊherbal medicines treated cells. Our aims are to demonstrate the feasibility of using the microarray/CD system for classification of different herbal compositions or different'preparatiόns and to find the predictor genes (marker genes) for the Huang Chin Tang prescription, the long-terfn goals are to find the correlation ofthe biochemical ingredients in each herbal composition with the gene expression profiles of various treated cells and to decipher the molecular pharmacological mechanisms ofthe Chinese herbal medicines in a rational fashion. Materials and Methods. 1. Development of. a "cell banking system. Purpose: Microarray system is a sensitive detection method to monitor gene expression patterris of .cells: If is necessary to build a Cell Banking System with a Master Cell Bank " (MCB) -ιnd a Working Cell Bank (WCB) to minimize cell variability for herbal medicine testing. Scope: The Cell Barik System is used for all types of cells in microarray studies. Apparatus: CO 2 Air- Jacketed Incubator (NUAIRE™ DH autoflow) Centrifuge (KUBOTA 2100) Freezing vial (Cofning Costar, Cat. #430659) Tissue culture flask 750 ml (Falcon, Cat. #3045) Tissue culture dish 150x25 mm (Falcon, Cat. #3025) -71- Cell: Jurkat T cell from Dr. Alexandra Ho Reagents: RPMI Medium 1640 (GIBCO BRL, Cat. #31800-014) Dimethyl Sulphoxide (DMSO) (Sigma, Cat. #D-2650) Fetal Bovine Serum (HyClone, Cat. #SH30070.03, Lot#AGL7258) 2-mercaptoethanoi (GIBCO BRL, Cat. #21985-023, 5xl0 "2 M) Media: I. Culture medium: 90% RPMI + 10% fetal bovine seru + 2- mercaptoethanol (5x10 "5 M) II. Freezing medium: 90% RPMI 1640 + 10% DMSO Procedure: A. Master Cell Bank . . . 1. Follow the standard sterile procedure of cell culture. 2. Seed Jurkat T cells in culture medium in flask at 37°C with 5% CO 2 incubator. 3. After incubatiori for two days, count the cell number and spread the cells to two flasks. Note. Cell density is kept about 5x10 4 - 2x10 6 per ml. 4. Culture and count the cell hurnber until the number reaches 2x10 . 5. Collect the cells in 50 ml centrifuge tubes and spin at 1300 rpm (300xg) for 5 min. 6. Discard the supernatant and re-suspend the cell pellet with chilled freezing media. The cell number in each vial is about lxl 0 6 per ml. 7. Slowly freeze the cells by the following temperature profile: -20°C for 2 hr, -80°C for 24 hr and then place the cells in liquid N storage. Store a total of 20 frozen vials in MCB. ' . ■ ■. . ' B. Working Cell Bank . 1. Retrieve one vial o cell from MCB in liquid N 2 tank and quickly thaw at 37°C water bath. 2. Transfer the cells into, 10 ml of warm culture medium. 3. Spin down the cells at 1300 rpm (300xg) for 5 min. Discard the supernatant. Culture the cells with 20 ml medium in a flask. 4. Sub-culture the cells to 2 flasks. • . 5. Seed 5xl0 7 cells with 500 ml culture medium in each flask with stirring for a total of 2 flasks. Culture the celϊs'for 2 days. 6. Culture until the cell density reaches lxl0 6 /ml and a total volume of 1 L. 7. Prepare freezing media by adding 100 ml of fetal bovine serum and 10 ml of DMSO. -72- 8. Centrifuge, discard the supernatant and re-suspend the cells to 110 ml of freezing medium. 9. Dispense 1 rill to every freezing vial (10 million cells per vial) for a total of 100 vials. Slowly freeze the cells by the temperature profile described above. 2. Determination of growth inhibition concentration of herbal extract in cell cultures. Purpose: Most drugs are toxic to cells. This experiment is designed to examine the toxicity of herbal extracts in Jurkat T cells and to determine the growth inhibition concentration of herbal extracts that keeps the cells alive. Scope: This assay can be used in all kinds of herbal extracts to examine the toxicity. Apparatus: CO 2 Air- Jacketed Incubator (NUAIRE™ DH autoflow) Counting chamber (Hemacytometer, Reichert, USA) Microscope (Zeiss, Axiovert 100) Cell: Jurkat T cell _ ■ \ / • Reagents: RPMI Medium 1640 (GIBCO BRL, Cat. #31800-014) Fetal Bovine S<=ram (HyClone, Cat. #SH30070.03, Lot #AGL7258) 2-mercaptoethanol (GIBCO BRL, Cat. #21985-023, 5xl0 "2 M) Culture media: 90% RPMI + 10% fetal bovine serum + 2-mercaptoethanol (5xlO "5 M) ; ' ; ' Disposable sterile syringe filters (0.2 m, Corning, Cat. #21052-25) Herbal extracts: 1 Cordyceps Sinensis Myceliurii 2 ST 024 3 ST 044 4 ST 051 5 ST 093 6 ST 117 7 ST 123 8 ST 128 9 ST 134 1 0. ST 237 1 1. PHY90 6-303503: Complex mix composed of 4, 6, 7, 10 -73- 12. PHY906-284003 : Complex mix composed of 4, 6, 7, 10 Procedure: A . Herbal Extract Preparation 1. Dissolve 1 gram of herbal powder in 10 ml of 80 °C deionized water (neutral pH) in a polypropylene tube. 2. Incubate the tube at 80 °C water bath for 30 minutes with gentle shaking then centrifuge at 4000 rpm (1500xg) for 5 min to obtain the supernatant. 3. Centrifuge at 11000 rpm ( 1 . 4000 xg) for 10 min to collect the supernatant. 4. Using the disposable sterile syringe filter to filter the supernatant. B. Cell Survival Test 1. Culture Jurkat T cells as described above. 2. Dispense 1 ml of5xi0 5 /ml cells per well to 24-well culture plates. 3. Prepare 12 kinds of herbal extract solutions. The extract solutions must be freshly prepared and used inimediately. 4. Add 100, 50, 20, 10, 5 μl of each herbal extract solution into the 24 well culture plates to get the five different concentrations: 10, 5, 2, 1, 0.5 mg/ml. 5. Culture the cells for 24 h at 37°C in an ' incubator filled with 5% CO 2 . 6. Count the number of cells per well. Mix 10 μl of cell solution with 10 μl of Trypan blue dye and load into ceil counting chamber. 7. Count the four majof squafe areas to calculate the cell number. (number of cells in 4 areas)/4 x IO 4 x dilution factor = number of cells per ml 3. Profiling gene expression patterns of Jurkat T cells treated with herbal extracts. Purpose: Profile the gene expression patterns of Jurkat T cells treated with herbal extracts. A high-density nucleic acid microarray with colorimetric detection system is used. Apparatus: Heat block (Boekel, Model 110002) Spectrophotometer (Beckman, DV640) Centrifuge (KUBOTA 1910) Water bath (SLM. AMINCO, Model 800) Hybridization incubator (YIH DER OH-800) Heat sealer (TISH-300, TEW) Reagents: RNAzσl™ B (Tel-Test, Cat. #CS-104) -74- Oligotex mRNA Midi Kit (Qiagen, Hilden, Germany) Hybridizatiόri Bags (GLBCO BRL, Cat. #18278-010) EasiSeal (Hybaid, Cat. #HBOSSSEZlE) Glass slidόs (Matsunami, S2214, Japan) Aerosol Resistant Tips (ART tip) (molecular BlO-Products, Cat. #2139) Random hexariier primer (GIBCO BRL, Cat. #48190-011) Reverse transcriptase and 5x buffer (GIBCO BRL, Cat. #18064-014) RNase inhibitor (GIBCO BRL, Cat. #10777-019) Biotin- 16-dUTP (Boehringer Mannheim, Cat. #1093070) Dig-11-dUTP (Boehringer Mannheim, Cat. #1558706,) Blocking powder for hybridization (Boehringer Mannheim, Cat. #1096176) Bovine serum albumin (Sigma, Cat. #A2153) 20X SSC (Amresco, Cat.#0918S-2-20XPTM5L) SDS (Merck, Caf l 13760) Dextran Sulfate (Sigma, Cat. #D6001) Streptavidin-βτgalactosidase (GIBCO BRL, Cat. #19536-010) Anti-digόxigeriin-AP Fab fragments (Boehringer Mannheim, Cat. #1093274,) X-gal (GIBCO BRL, Cat: #15520-018) Maleic acid (Sigma, Cat. #M1125) N-lauroylsarcosin (Sigma, Cat. #L5777) Fast red TR/AS-MX substrate kit (PIERCE, Cat. #34034) Polyethylene glycol (Sigma, Cat.#P2139) Reagent Preparation: lx hybridization buffer (4X SSC, 0.1% N-lauroylsarcosine, 0.02% SDS, 1% BM blocking reagent) > „ . ' ., • > . 20x SSC , ' ■' .. ' , : • • ' 16 m 1% N-lauroylsarcosine 8 ml 10% SDS .. " • . 160 μl BM blocking powder .' ' . ' 0.8 g H 2 O 51 ' ml total 8,0 ml -75- Heat to 65°C to dissolve the powder then store at -20°C. 50% PEG-8000 (Polyethylene Glycol) PEG-8000 10 g. H 2 O up to 20 ml Heat to 65°C to dissolve then autoclave.. Aliquot and store at -20°C. 10x TBS (100 mM Tris, 1Ϊ5M NaCl, pH 7.4) NaCl 87.6 g H 2 O uρ to 1000 mi - 120 mM X-gal X-gal lOO mg * ' DMF 2 mi Store at -20°C. , X-gal Substrate Buffer (ImM MgCl 2 , 3 mM K 3 Fe(CN) 6 , 3 mM K 4 Fe(CN) 6 in IX TBS buffer) ' . ' ' " . 500 ml 10X TBS/pH7.4 ' : 50 ml Potassium Ferrocyanide 633.5 mg Potassium Ferricyariide 493.9 mg MgCl 2 101.6 mg Filter and store at -20°C. BM Blocking Dilution Buffer/pH7.5 (0.1 M maleic acid, 0.15 M NaCl) 1M Maleic Acid ' 100 ml 5M NaCl ' - • ' • • ■ 30 ml Solid NaOH 7.5 g H 2 O up to 1000 ml 10% Blocking Reagent 100 ml Blocking Powder 10 g -76- Blocking Dilution Buffer (no tween 20) 100 ml Heat to 70°C then autoclave: Store at 4°C. 20% Dextran Sulfate Dextran Sulfate ' 2 g H 2 O up to 8 ml Autoclave then store at -20°C. DEPC-treated water 800 ml Diethyl Pyrocarbonate (store at 4°C) 400 μl H 2 O 800 ml Put in 37°C shaking water bath for 4 hrs and then in 37°C warm room for overnight. Autoclave the solution for 45 minutes of 25 minutes each for two times. Procedure: A. Preparation of Jurkat T Cells 1. Thaw one vial of Jurkat T cells. Transfer to 10 ml growth medium. The number of vials to thaw depends ori the number of test to perform. In general, one vial of cells is needed for 2 herbal extract tests. 2. Re-suspend the cells in 50 ml of culture medium. 3. Incubate the cells for one day. Add 150 ml of culture medium and divide into two flasks with 100 ml each. 4. Culture the cells for 3 days. 5. Change the medium arid distribute 4 flasks with 100 ml of medium each. 6. Culture the cells for 2 days. 7. Count the cell number. Collect the cells and spin down. Re-suspend the cell pellet with culture medium to 5x1,0^ cells per ml and 100 ml per flask. 8. Culture the cells for 3 hrs before adding the herbal extract. B. Herbal extract treatment' Prepare the herbal extract (see herbal extract preparation). According to the growth inhibition concentration of herbal extract determined by the cell survival experiments, calculate the ' 50% growth inhibition concentration of each herbal extract. ■ : -77- 3. Define the 50% growth inhibition concentration as H. Treat the cell with serial dilution of herbal extracts with the following concentrations: H, H/2.5, H/5, H/10, H/20. 4. Culture the cell for 24. hrs. 5. Collect cells and count the number of cells. Centrifuge to obtain cell pellet. 6. Wash the cell pellet with 1 xPBS once. 7. Discard the supernatant. The cell pellet is ready for total RNA isolation. C. Isolation of total RNA . , ' . 1. Add 1 ml of RNAzol^ B per IO 7 cells. Homogenize the cell pellet but do not vortex. 2. Add 0.1 ml chloroform per ml of homogenate, cover the samples tightly, shake vigorously for 1 min (do not vortex). Place on ice for 15 min. 3. Centrifuge at 12000 . rpm (1350Qxg) at 4°C for 15 min. 4. After centrifugation, the hornogenate develops two phases: a lower blue phenol- chloroform phase and a colorless upper aqueous phase. DNA and proteins are in the interphase and the organic phase. Transfer the aqueous phase to a new tube, add an equal volume of isόpropanol and store the samples at -80°C. Note. The range of isopropanol addition, i_ rom 0.7 to 1 volume ofthe aqueous phase solution. 5. Keep the samples at -80°C until rise. Let the sample completely thaw before centrifugiήg and mix 2 to 3 times by inverting the tube. Centrifuge samples for 15 min at 13000 rpm (15000 ' x^). 6. Remove the supernatan and wash the RNA pellet once with 1 ml of 75% ethanol. Centrifuge for 3 ' niiri at; 13000 rpm (15000 χ g) and at 4°C. 7. Discard the supernatant. Dry the pellet under vacuum for 1 min. Note. Do not let the RNA pellet dry completely. It will greatly decrease its solubility. 8. Dissolve the RNA pellet in' 50- 100 μl of diethylpyrocarbonate (DEPC) - treated water by pipetting. Note. If the pellet is hard to dissolve, incubating the pellet for 10 - 15 min at 60°C may heϊp ' ' •" > '■ '• '" " " " 9. Measure absorbance af 26θ nm (A 260 ) and 280nm (A 280 ) with a spectrophotometer. Concentration analysis: OD 260 ' x 40 ' ng/μl dilution factor = total RNA (ng/μl). D. Isolation of oly-A+ mRNA from total RNA 1. Determine the amount of starting RNA and the appropriate volume of Buffer OBB and Oligotex Suspensiori solution to be added in the RNA solution according to the Table 17. -78- Table 17. Buffer amounts for Oligotex mRNA Spin-Column Protocol. The following procedures are based on using 500 μg total RNA as an example. 2. Add 500 μl of 2x Binding Buffer and 30 μl Oligotex Suspension to the total RNA sample. Mix the contents thoroughly by flipping the tube. 3. Incubate the sample for 10 min at 70°C. 4. Incubate for 20 min at room temperature. 5. Centrifuge for 2 miή at maximum speed (14000 to 18000 χ g) and aspirate the supernatant. . ' ' ■ ' - 6. Re-suspend the pellet in 400 μl of Wash Buffer OW2 and transfer onto spin column and centrifuge the spin column for 1 min. 7. Wash with 400 μl Of OW2 and centrifuge as above. 8. Add 20 μl of preheated (70°C) Elution Buffer onto the column and re-suspend the resin. Close the microcentrifuge tube. 9. Put the spin column with 1.5 ml microcentrifuge tube at 70°C for 3 min. 10. Centrifuge the column at maximum speed for 2 min at room temperature. 11. Elute again, (repeat step 8 to 10 to get better yield) . E. cDNA Labeling '. •' ' 1. Mix 2μg mRNA, 1 μl of control plants mRNA for single color label (Hat22: lxlO 9 , Rbcl: 5xl0 8 , Ga4: MO 8 , Rca: 5xl0 7 , Asal: lxlO 7 , Atps: 5xl0 6 molecule/μl), 6 μl of 50 mM random hexamer and DEPC-H 2 O to 28.88 μl final volume. For dual-color mode, use 2 μg of mRNA each in Biotin or Dig labeling and individual addition of control plants mRNA: 1. biotin labeling: Hat22: lxlO 8 , Rbcl: 5 χ l0 7 , Ga4: 2 χ l0 7 , Rca: lxlO 7 , Asal: lxlO 7 , Atps:- 1- χ lO 7 , Hat4: lxl0 7 /μl. 2. Dig labeling: Hat22: lxlO 7 , Rbcl: lxlO 7 , Ga4: lxlO 7 , Rca: lx.10 7 , Asal: MO 8 , Atps: 5xl0 7 , Hat4: 2xl0 7 /μl. 2. Denature for 10 min at 70°C, then chill quickly in ice for 5 min. 3. Add 10 μl of 5x first strand buffer, 5 μl of 0.1 M DTT, 1 μl of 25 mM dATP, dCTP, dGTP mixture, 1 μl of 2 mM dTTP , 2 μl of 1 mM biotin- 16-dUTP, or Dig-11-dUTP (1 -79- mM), 0.63 μl of 40 U/μl RNAsin and 1.5 μl of Superscript II (reverse transcriptase, GIBCO BRL)(200 U/μl). 4. Mix well and incubate for 10 min at 25°C, then for 90 min at 42°C. 5. Stop the reaction for 5 min at 94°C. 6. Add 5.5 μl of 3 M NaOH for 30 min at 50°C. 7. Add 5.5 μl of 3 M CH 3 COOH for 30 min at 50°C. 8. Precipitate the labeled cDNA by adding 34 μl of water, 50 μl of 7.5 M ammonia acetate, 10 μg of linear polyacrylamide as carrier and 380 μl of absolute alcohol. 9. Incubate the sample for 30 min at -80°C. Centrifuge at 13000 rpm for 15 min. 10. Wash the pellet with 1 ml of 70% ethanol and centrifuge at 13000rpm for 5 min. 11. Dissolve the pellet in 36 μl of autoclaved H 2 O. For dual color, combine two labeled cDNA together. F. Array Hybridization 1. The filter membrane carrying the 9600 EST PCR products is pre-hybridized in 5 ml of lx hybridization buffer (4X SSC, 0.1% N-lauroylsarcosine, 0.02% SDS, 1% BM blocking reagent (Boehringer Mannheim)), and 50 μg/ml salmon sperm DNA (GIBCO BRL) at 63°C for 1.5 hours. Note. You can prepare 80 ml of lx hybridization buffer and store it at -20 o C. thaw the buffer at 60°C before use. 2. Stick one side of adhesive EasiSeal ® square to a clean glass slide and place the pre- hybridized membrane in the center ofthe square with the spots facing up. 3. Mix the probe with 2 μl of poly-d(A)10 (10 μg/μl) and 2 μl of human Cot-1 DNA (10 μg/μl) (GIBCO BRL) and 40 μl of 2x hybridization buffer to 80 μl final volume. 4. Denature the probe mixture at 95°C for 5 min and then cool on ice. 5. Seal the filter membrane with the probe solution in the hybridization bag. 6. Incubate at 95°C for 5 fniή and then at 63°C for 12-16 h (overnight). 7. Wash the filter rriembrane twice with 5 ml of 2x SSC, 0.1 % SDS for 5 min at room temperature. - -, ' 8. Wash three times for 15 min each with 5 ml of O.lx SSC, 0.1% SDS at 63°C. 9. Block the filter membrane with 5 ml of 1% BM blocking reagent containing 2% dextran sulfate at room temperature for 1 h. 10. Incubate with 5 ml , mixture containing 700x diluted Streptavidin-β-galactosidase (1.38U/ml, enzyme activity)(GIBCO BRL), lOOOOx diluted anti-Digoxigenin-alkaline -80- phosphatase (0.075U/ml, .enzyme activity)(Boehringer Mannheim), 4% polyethylene glycol 8000 (Sigma), and 0:3% BSA in lx TBS buffer at room temperature for 2 hours. Note. This formula is for dual-color mode. For single color mode, anti-Dig- AP is not needed and the incubation time can be reduced to 1 hour. 11. Wash with lx TBS buffer three times for 5 minutes each. 12. Freshly prepare X-gal . substrate solution (1.2 mM X-galj ImM MgCl 2 , 3 mM K 3 Fe(CN) 6 , 3 mM __ 4 Fe(CN) 6 in lx TBS buffer) by mixing 50 μl of 120 mM X-Gal and 5 ml of X-Gal substrate buffer. Immerse the filter membrane in the X-gal substrate solution for 45 min at 37°C with gentle shaking. 13. Wash with lx TBS. 14. Dual color development: stain the membrane with 5 ml of Fast red TR/naphthol AS- MX substfate (Pierce- Rockford, IL) at room temperature for 30 minutes with gentle shaking. i . 15. Wash with deionized watef. Stop the reaction with lx PBS containing 20 mM EDTA for 20 min. ^ ' ' 16. Air dry the filter membrane. Results. 1. Determine the growth inhibition concentrations of herbal extracts in cell cultures. Each herbal extract has different cellular toxicity, thus it is necessary to determine the growth inhibition concentration of every herbal medicine before treating the cells. Five serial dilutions of herbal extract ' (10,' '5, 2, 1, 0.5 mg/ml) were added to 5xl0 5 /ml cultured cells and incubated for 24 hours in an incubator at 37°C with 5% CO . The numbers of survival cells at different concentrations of herbal extracts are shown in Table. 18. -81- Note: The original cell number was 5x10 ml. The number increased to 10x10 /ml after 24h incubation.."-" indicates all. dead cells. . The cell number withpμt herbal extract addition doubled after 24 hours incubation. On the other hand, the number, of survival cells varies with different herbal medicine treatments. We chose the 50% growth inhibition concentration (IC 50 ) as the high concentration and one- tenth of it as the low concentration. In order to maintain consistence ofthe Jurkat T cell line, a cell banking system was established. In the cell bank, a total of 100 vials of cells (10 million cells per vial) were frozen in a.-150°C freezer. 2. Molecular classification of herbal medicines by nucleic acid microarray analysis. Analysis of 3 single-element herbal medicines. Three single-element herbal medicines, Cordyceps Sinensis Mycelium (CSM), ST024, and ST117 were used to treat the cell cultures as described in the methods section. Gene expression measurements were performed by using microarrays of 13824 cDNA fragments each representing a distinct human transcript. For the data analyses, gene spots of high data quality were selected. The selection was based on signal to background ratio greater than -2.5 or the coefficient of variation (CV) of spot area smaller than 10%. All the data sets were normalized with the control cells, which received no herbal treatments. The spot intensity was rounded up to 10 for those intensities that were less than 10. Based on the selection criteria, a total of 492 genes with differential expression ratio greater -82- , than 1.5 were selected for cluster analysis. These data pre-processing procedures were performed by the program "DataExfr act" and "Ratio2" developed in-house. These 492 genes ' were cluster analyzed by the average-linkage method. The distance between genes is used as the linear correlation coefficient or resemblance coefficient. The cluster analysis programs, Cluster and TreeView, were based on hierarchical clustering method and were written by Dr. Michael Eisen of Stanford University (Eisen, 1999, Eisen et al., 1998). The results are shown in Figure 12. From Figure 12C, one can clearly identify that three different herbal medicines of high and low concentrations are each clustered together. For instance, CSM-L is closer to ;CSM-H and less similar to ST024 or ST117 in the clustering tree. A different clustering algorithm, the self-organizing map algorithm, which is based on non- hierarchical method, yields the same results (data not shown). From the clustering results shown in Figures 12A & 13 A several features are noted. (1) 4 genes were up regulated by ST117 treatment but down regulated by other herbal treatments (Figure 12B). (2) 34 genes were down regulated by CSM.treatment but up regulated by others (Figure 13B). (3) 2 genes were up regulated by all the three herbal treatments, one is Malic enzyme 2 and the other one is an anonymous gerie (clone ID: 328351) (Figure 13C). (4) 12 genes were highly induced by the high concentration treatment and less induced by the low concentration treatment in all the three herbal medicines (Fig. 13D). Analysis of 2 preparations of multi-eleriient herbal medicines. Two batches of Huang Chin Tang, PHY906-303503 (#11) & PHY906-284003 (#12), each with low and high concentrations were used to treat the cell cultures in three independent experiments. The gene expression profiles were acquired with microarrays of 9600 non-redundant cDNA elements. After the data pre-processing procedures as described above about 5000 genes were selected for the subsequent data analysis ' /There were 3 repeats for each herbal treatment. For data analysis, we use a modified method . based on the one reported by Slonim et al. (Slonim, 1999). The following algorithm is designed to search for the candidate marker genes that have high differential expression ratios, but low deviation in the three repeats. We designate a P(z ' ) value to account for the gene i with the aforementioned features. P(z ' ) = square root of (Σ(μ m μ c ) )/ (σ c +∑σ m ) μ : Mean expression levels in three repeat experiments for herbal treated cells (μ m ) or untreated control cells (μ c ).., -83- σ : Standard deviation " of the expression levels in three repeat experiments for herbal treated cells (σ m ) or untreated control cells (σ c ). We calculated the P(z ' ) value for each gene and selected 500 genes with the highest scores as candidate genes ; for cluster analysis (Figure 14). The values of each gene were averaged over the 3 repeats: As shown in Figure 14B, two different concentrations of #12 are clustered together (12-H & 12-L). The higher concentration ofthe #11 preparation is closer to the #12 preparation cluster than the lower concentration of #11 preparation. However, all these clusters have similar resemblance coefficient (distance between clusters) compared with the tree shown in Figure 12. These results suggest that the gene expression profiles of #11 and #12 preparations of Huang Chin'. Tang are similar. The results are justified based on the fact that these two preparations are based on the same herbal medicine mix. Several features are noted in the expression profiles illustrated in Figures 14A & 15. The averaged gene expressiori levels are shown in Figure 15A The Boxl encloses genes that were down regulated in #11-L treated cells but up-regulated in others. These genes include 2 tRNA synthetase (isoleucine ahcl methion), RNA polymerase II polypeptide B (Clone ID 42020), KIAA0212 gene (Clone tt) 310497, containing ATP/GTP-binding site motif A), and KIAA0577 (Clone lD 29263,- ATP-dependent RNA helicase). It is interesting to note that 3 out ofthe 6 genes were involved ih the RNA replication. Box2 encloses the genes that were up regulated by all the #11 arid #12 treatments. Box3 encloses the genes which showed no response by #11-L treatment but were down regulated by the others. Box4 encloses the genes that were highly repressed ' by low concentration herbal treatment but were less repressed by high concentration herbal treatment. Finally, in Boxl and Box3, the expression profiles of #11 treated cells are different from the profiles generated by the other 3 treatments. This result is consistent with the finding depicted in Figure 14B. Combining the data sets ofthe gene expression profiles ofthe 3 single-elements and the 2 preparations of multi-elenierit hefbal medicines together, a couple features are noted as an illustration. The KiAA0212 gerie (Clone ID 310497, containing ATP/GTP-binding site motif A) was highly induced by all the 'high concentration herbal treatments except that it was only mildly induced by the #11-ϊi and the CSM treatments. Two genes, an anonymous gene, (Clone ID 510908) and Proteasorrie chain 7 precursor (Clone ID 70088) were highly up-regulated on all the treatments except down-regulated by the CSM treatment. -84- We next worked on the crux to cluster analyze the gene expression profiles ofthe 5 different types of herbal medicine treatments. The data pre-processing procedures were performed as aforementioned and 500 genes were selected for cluster analysis. A hierarchical clustering was performed by the program "Cluster" described above. The hierarchical tree was cut at the position where the range of distance between clusters is the greatest (Romesburg, 1989) and the result is shown in Figure 16A. The 3 single-element herbal medicines, CSM, ST024, and STl 17 are clustered together. Out ofthe 2 different batches ofthe multi-element herbal medicine, #12-H, #12-L and #1 IH are clustered together and the #11-L stands by itself. The result suggests that higher similarity exists between the #11 and #12 as compared with that of CSM, ST024, and STl 17. In order to better classify the different herbal medicines, the data analysis algorithm was improved- by standardizing all the data sets so that the expression level of each gene across the different data sets has zero-mean and unit- variance (Tavazoie et al., 1999; Chen et al., 1999). This yields the transformed variables: χi = (xι-μ x )/ σ x μ x : mean expression levels in the' data set σ x : standard, deviation ofthe expression levels in the data set Xi : un-transformed gene expression level χi : transformed gene expression level After standardizing the data set,, #11 and #12 are clustered together as shown in Figure 16B. The ST024 and ST 117 are clustered together and the CSM is in an independent cluster. Furthermore, the clustering also suggests that CSM is more similar to #11 and #12 than to the ST024 and STl 17. Another clustering algorithm, self-organizing maps, was performed with the same standardized data sets and yields the same result as the hierarchical clustering (Figure 16C). Class predictors for discriminating #11 and #12 herbal treated expression profiles. The above cluster analyses for .the #11 and #12 show that they are similar and further classification is difficult by the hierarchical clustering or self-organizing maps methods with the data set containing the 500 genes ofthe highest P(z) values. We then modified the algorithm to select genes with larger expression ratio difference between #11 and #12 herbal treated cells, but -85- smaller variation in the two herbal treated cells. The T(z ' ) value is defined to score this feature as following: ' • • ' .' T(i) = log(μ 11 ) - log(μ 12 ) / (σ π +σ 12 ) μ: Mean expression. ratio? in three time experiments for #11 treated cells (μπ) or #12 treated cells (μ 12 ) : • * • σ : Standard deviation of the expression ratios for #11 treated cells (σπ) or #12 treated cells (σπ) We calculated the T.(t) . value for each gene and selected 50 genes with the highest scores as class predictors.(Figure 17). 18 genes were up-regulated by the #11 treatment and down-regulated by. the #12,treatment. The rest ofthe genes were up-regulated by the #12 treatment and down-regulated by the #11 treatment. We then used these class predictors to classify two test herbal prepafations based on a modified method described by Golub et al, 1999. ■ Two different batches ofHuang Chin Tang preparations, PHY010401 (#16) & PHY010402 (#17) were obtained from Sun Ten Pharmaceutical Co. and were used for the class prediction test. The gene expression profiles of #16 and #17 preparations were normalized with the expression profile ofthe untreated control cells and standardized with the class predictors. Each predictor g; votes for either #11 or #12 herbal preparation depending on whether its expression level Xj is closer to #11 or #12. The vote for each gene is given by Vi = |xι - (μ m +μ c )/2|, where ; . ■ ' •" • ' ' •" μ : Mean expression ratio, in three repeat experiments for #11 (μπ) or #12 herbal treated cells (μ 12 ) : l . The average votes V ' π and V 12 were collected from the predictor genes correlated with the predictor on #11 and #12, respectively. The prediction strength (PS) reflects the margin of victory and it was defined as PS = (Vπ-V 12 )/ (Vπ+V 12 ). If the PS was greater than 0, it indicated that the herbal prepa ation was more similar to #11 and less similar to #12. The results obtained from the analyses on,#16 and #17 indicated that #16-H was similar to #11 -86- (PS=0.1) and #16-L, #17-H, and #17-L were similar to #12 (PS=-0.29, -0.21 and -0.2, respectively). Based on the information of #16 and #17 preparations, this test failed to correctly identify #16-L as more similar to #11. Discussion. Characteristic gene expression profiles for herbal medicines treated cells. The predictor genes are selected from the differentially expressed genes in two herbal treated cells. These genes represent the cellular responses to the herbal medicine treatments. In this study, we have identified some interesting genes based on their responses across different herbal treatments (Figures 13, 15, and 17). These genes are valuable assets in studying the signaling pathways of cells in response to the herbal stimulation and in deciphering the molecular pharmacological mechanism of herbal medicines. Classification of herbal medicines by nucleic acid microarray analysis. In summary, a two-step classification procedure is proposed. An initial classification procedure based on the standardized data sets and the clustering algorithms is performed and followed by a final classification procedure with' the class predictors. All the procedures can be integrated in a computer program. In these preliminary studies, all the genes have the same contribution for classification. When the data sets are large enough, the weight for each gene (or predictor) can be acquired from the linear correlation coefficient (Golub et al., 1999, Chen et al., 1999). The #11-L failed to be clustered with #11-H by significant association. A similar preparation of #11 , the #16 preparation yielded the same results. It was interesting to discover that no matter what clustering algorithms were applied, the #11 and the #16 preparations did not yield the expected result ' s. Even with independent experiments, the results remained the same. The reason behind the failure will be investigated with detailed information of #11 and #16 preparations. Quality control and arialysis in microarray system. The quality of acquired microarray data and the choice ofthe' statistical analysis methods are both important factors for achieving meaningful results. We have recognized that the variation among arrays contributes to errors in measuring gene expression levels. Based on the data in this report, we have found that for every herbal preparation, the high concentration treated expression profiles always cluster with its lower concentration counterpart (Figures 16B and 16C) and we could classify the STl 17, ST024, CSM and Huang Chin'Tang with two different clustering methods. In the past three months, the array quality has-been improved to have less than 7% CV. We have also set up a standard procedure for assessing' the quality of every batch of arrays fabricated in the lab. All . -87- these experimental findirigs and improvements in the microarray technologies suggest that classification and characterization ofthe Chinese herbal medicines by microarray system are feasible. References for Example 15.' Bertucci-F; Bernard-K; Loriόd-B; Chang- YC; Granjeaud-S; Birnbaum-D; Nguyen-C; Peck-K; Jordan-BR (1999) Sensitivity issues in DNA array-based expression measurements and performance of nylon microarrays' for small samples Human Mol. Genetics 8(9): 1715-1722. Bittner-L, Trent- J, Meltzer-P (1999) Data analysis and integration: of steps and arrows. Nature Genet. 22, 213-215 Brown-PO; Botstein-D (1999) Exploring the new world ofthe genome with DNA microarrays. Nature genetics 21 (1) suppleriient, 33-37. Chen-JJ; Wu-R; Yang-PC; Huang-JY; Sher-YP; Han-MH; Kao-WC; Lee-PJ; Chiu-TF; Chang- F; Chu-YW; Wu-CW; Peck-K (i 998) expression patterns and isolating differentially expressed genes by cDNA microarrSy systeήi 'with colorimetry detection. Genomics. 51: 313-24. Chen-Y, Bittner-M, Dougherty-ER (1999) Issues associated with microarray data analysis and integration, Informatiori supplementary to article by Michael Bittner, Jeffrey Trent and Paul ι Meltzer (Nature Genet. 22, 213-215). Duggan-DJ; Bittner-M; Cheri-Y; Meltzer-P Trent -JM (1999) Expression profiling using cDNA microarrays. Naturfe genetics 21 (1) supplement, 10-14. : ' . ' . . . Eisen-M (1999) Cluster and ; Tfeeview manual, (rana.stanford.edu/software) ft ,' •'.! •- . . Eisen-M, Spellman-PT, BrOwn-PO, Botstein-D. (1998) Cluster analysis and display of gene- wide expression patterns. Proc. Natl. Acad. Sci. USA 99:14863-14868 ' ". "' Golub-TR, Slonim-DK, Tamayo-P, Huard-C, Gaasenbeek-M, Mesirov-JP, Coller-H, Loh-ML, Downing- JR, CaligiurirMA, Bloomfield-CD, and Lander-ES (1999) Molecular classification -88- of cancer: class discovery and class prediction by gene expression monitoring. Science Oct 15: 531-537 " '• : ' .. .. " > '■ ■ - . • , Lander-ES (1999) Arr y of hope. Nature genetics 21 (1) supplement, 3-4. Romesburg-HC (1989) Cluster analysis for researchers. Chapter 16: How to make classifications. P203-216, Kriegef Publishing Co. Malabar, Florida, USA Slonim-DK, Tamayo-P, Mesirov-JP, Golub-TR, Lander-ES (1999) Class prediction and discovery using gene expression data. (www.genome.wi.mit.edu/MPR) Tavazole-S, Hughes- JD, Campbέll-MJ, Cho-RJ, Church-GM. (1999) Systemic determination of genetic network architecture. Nature genetic 22: 281-285 Example 16. Identify characteristic gene expression profiles induced by an herbal medicine As stated in Example 15, the prescription ofthe Chinese herbal medicine Scute and Licorice combination (Huang Chin Tang) stops diarrhea, relieves spasms and clears fever. The ingredients of Huang Chin Tarig are Scute, Peony, Licorice and Jujube. In this study we used the nucleic acid microarray technology to study the gene expression profiles induced by herbal medicines in mammalian bells. To investigate the characteristic expression profiles induced by the Huang Chin Tang, Jurkat' T cells were treated with 5 batches of Huang Chin Tang (PHY01040; #16, PHY010402; #17, PHY03061; #18, PHY03062; #19 and PHY02231; #20 obtained from Sun Ten Pharmaceutical Co.) by 5 concentrations (1/2, 1/2.5, 1/5, 1/10, and 1/20 ofIC 50 ). • • '• ' -. . • ■ ' Nucleic Acid microarray with a two color detection method was employed to measure the expression profiles The mRNA extracted from herbal treated cells was labeled with digoxigenin and the mRNA: extracted from untreated cells was labeled with biotin. Arrays of 9600 features were employed and the procedures described by Chen et al. (Genomics, 51, 313- 324, 1998) were adopted in the 'experiments. For data pre-processing, only array spots of high data quality were selected: The selection was based on signal to background ratio greater than 2.5 and 1.5X differential expression ratio. By these criteria, 1081 genes were selected for further statistical analysis. Non-hiefarchical cluster analysis programs such as the Genecluster program developed in Massachusetts Institute of Technology (Tamayo et al, 1999) was -89- employed to categorize the expression profiles. The Genecluster program is based on self organizing map (SOM) principle. A 6X4 SOM clustering of expression profiles are shown in Figure 18A. The details of gene expression profiles for the selected clusters are shown in Figure 18B. In these clusters, clusters 3 and 20 (labeled c3 and c20) were selected for that the gene expression levels increase with higher herbal concentrations. Similarly, clusters 5 and 9 for that the gene expression leyels decrease with higher herbal concentrations. Cluster 23 collects genes whose expression levels were up regulated compared with that in untreated cells and cluster 0 for down regulated genes. These expression profile clusters are further condensed into two major groups, A & B. Group A collects genes up regulated by herbal treatment and Group B collects. genes down regulated by herbal treatment. The expression profiles in Group A & B form the basis of a characteristic expression data set for an herbal preparation. The same procedures were repeated for 5 different batches of Huang Chin Tang, and 952 genes were selected to establish trie characteristic expression profile database ofthe Huang Chin Tang. As shown in Figure l9, ' by the : aforementioned procedures, a gene can be categorized as Group A, B or none (non-A arid ήon-B) and its expression profile can be represented by 1, -1, and 0 respectively. The number of different gene expression profiles between batch #1 and batch #2 are 3 in Group A (Gene 6,'7, and 8) and 2 in Group B (Gene 15 &16). By the same principle, the number of different expression profiles between batch #1 and #3 are 10 in Group A and B and the number is '.ii' between batch #2 and batch #3. These numbers indicate that batch #1 and #2 are more similar than batch #3. This principle was applied to classify 5 different batches of herbal preparations. The following algorithm is designed to calculate the distance between a pair of herbal preparation batches, i and j. dj j = Σ δ(Xι, X j ) . . (Hamming distance) The gene X in i batch of preparation is assigned to Group A, B or none if X i o X j , δ(X i ;X j j^ l.': ' if X i = X j , δ(X i ,X j j - 0. ' ' We calculated all the dy value between pairs of herbal preparations for cluster analysis. The analysis programs, Kitscri'Ciuster was based on hierarchical clustering principle and was written by Dr. Joseph Felsenstein of Washington University (http://evolution.genetics.washirigton.edu/phylip.html). From the Hamming distance table (Figure 20), one can clearly identify that the shortest distance lie between batch #17 and batch #18 and that batch #17 is similar to #18. Batch #16 also similar with batch #17 and #18 but -90- batch #19 is dissimilar to the st ofbatches. The results were confirmed by HPLC analyses as described below. By HPLC, the chemical composition ofthe 5 batches of herbal preparations was analyzed. Four major peaks (BG, B, Gly, and Pf) in the chromatograms were selected for statistical analyses. Two additional parafrieters, BG+B and BG/B, were included in plotting the 6- coordinate radar graph as shown in Figure 21. The distance on each coordinate is the integrated intensity of that particular chemical constituent in the chromatogram. In general, #16, #17 and #18 were similar in their constituent content (baicalin and baicalein) of Scutellariae Radix, which are within 33.55-36:08; while the amount ofthe same constituents are higher in #19 and #20 (42.49 and 44.96, respectively). The resemblance of #16, #17, and #18 can be seen from the coincident radar plots in Figure 21B. In order to identify, the unknown herbal medicine based on the characteristic expression profile database established as described as above, Jurkat T cells were treated with a tester sample #17 in 5 concentrations to set up the characteristic expression data set for the tester. The Hamming distances between the tester and each ofthe data sets (#16, #17, #18, #19 and #20) in the characteristic expression database were calculated and the scores are: #16: 502, #17: 405, #18: 402, #19: 699, and #20: 531. These data show that the tester is most similar to #17 having a lowest Hamming distance score of 405. The example demonstrates that this invention teaches a method ' to "identify unknown herbal medicine based on the gene expression profiles induced by the herbal-medicine in mammalian cells. The identity ofthe unknown herbal medicine can be infefred by aligning the characteristic expression profiles with a collection of characteristic expression profiles of herbal medicines in an HBR Array. Based on the characteristic expression database, marker genes and signature expression profiles can be deduced for ari herbal medicine for studying its pharmacological mechanisms and for optimizing the formulation of a complex herbal preparation. For this example, 5 different batches of Huang Chin Tang' preparations (#16, #17, #18, #19 and #20) were obtained from Sun Ten Pharmaceutical Co. and a characteristic expression profile database was constructed based on aforementioned procedures. For each gene, the consistency of expression profiles in the database was scored by the coefficient of variation (CV value): CV = σ / (Σ μi /n) ' " : - μ;: Mean expression ratios for.#i treated cells. n : Number ofthe data set, n = 5 in this case. σ : Standard deviation of trie expression ratios for #16, #17, #18, #19 and #20. -91- Since the CV reflects the variation of data, the marker genes for an herbal medicine were selected based on the CV score. The top 50 genes with the minimum CV scores were selected. Figure 22 shows 25 marker genes with up regulated signature profiles and 25 marker genes with down regulated signature profiles for Huang Chin Tang. The characteristic expression profile database can be used to infer the expression profiles of individual chemical constituents in a mixture as complex as an herbal medicine if the amount ofthe chemical constituents can be semi-quantitatively determined. In this example, the chemical composition of an herbal medicine is determined by high performance liquid chromatography. The integrated intensities of 4 chemical constituents in five batches of Huang Chin Tang preparation were quantified by HPLC analysis. The gene expression ratios for each batch of herbal preparation were calculated by taking the median ofthe expression ratios induced by 5 concentrations of herbal preparation. The correlation between a constituent and a gene expression profile was quantified by the Pearson correlation coefficient. The Pearson correlation coefficient for gerie x and the constituent y is: R = (1/n) Σ(Xi-μ x )( yf-μy) " / σ x c_ y , i= 1 to n n : Number of the herbal preparatiori,. n = 5 in this case. μ x : Mean expression ratios in five herbal preparation for gene x. μ y : Mean integrated intensity in five batches of herbal preparation for constituent y. x;: Gene expression ratios in #i herbal preparation for gene x. yi.- Integrated intensity in #i herbal preparation for constituent y. σ: Standard deviation of the expression ratios (σ x ) or integrated intensities (σ y ) for five herbal preparations. • ' • v' ;- Several genes whose expression levels highly correlated (with |R|>0.99) with the amount of chemical constituents in Huang Chin Tang were identified for each constituent. For example, the R value between the gene (clonelD: 67185) and Glycyrrhizin was 0.998 (Figure 23A). One the other hand, the gene (clorielD: 344720) whose expression levels increase with the decrease of Wogonin(WG) has an R value of -0.997 (Figure 23B). In addition to the above two examples, 191 and 170 genes were highly correlated with individual constituents with R value > 0.9 and R value < -0.9, respectively. For instance, 17 and 18 genes were positively and negatively, respectively, correlated with Albiflorin (Af) (Figure 24). This example teaches a -92- method to profile gene expression for individual constituents in a mixture without isolating them to perform the expression analyses one constituent by another. References for Example 16. US Patent documents . Stoughton-Roland and Friend-SH, USP# 5965352: Methods for identifying pathways of drug action . , ' ■ ... ' ■ • • ; Brown-PO and Shalon-TD USP#5807522: Methods for fabricating microarrays of biological samples ■ •■ _ • ; "■ . : , Lockhart-DJ, Brown-EL, Won r .GG, Chee-MS, and Gingeras-TR. USP#6040138: Expression monitoring by hybridization:, to high density oligonucleotide arrays Ladunga-I. USP#5987390: Methods and systems for identification of protein classes • •• ' ■' • • . ' -. . ' ,. : ' Mahant-S, Shivaling-S, Vivek-G. USP#5951711: Method and device for determining hamming distance between two multi-bit digital words Foreign . Patent documents Brown-PO and Shalon-TD ,EP#913485A1 : Method and apparatus for fabricating microarrays of biological samples Other Publications Bertucci-F; Bernard-K; Loriod-B; Chang- YC; Granjeaud-S; Birnbaum-D; Nguyen-C; Peck-K; Jordan-BR (1999) Sensitivity issues in DNA array-based expression measurements and performance of nylon microaprays for small samples Human Mol. Genet. 8(9): 1715-1722. Bittner-L, Trent- , Meltzer-P (1999) Data analysis and integration: of steps and arrows. Nature Genet. 22, 213-215 ■ • Brown-PO; Botsteih-D (1999) Exploring the new world ofthe genome with DNA microarrays. Nature genet. 21 (1) supplement, 33-37. . -93- Chen-JJ; Wu-R; Yang-PC; Huang-JY; Sher-YP; Han-MH; Kao-WC; Lee-PJ; Chiu-TF; Chang- F; Chu-YW; Wu-CW; Peck-K (1998) expression patterns and isolating differentially expressed genes by cDNA microarray, system with colorimetry detection. Genomics. 51: 313-24. Chen-Y, Bittner-M, Dougherfy-ER (1999) Issues associated with microarray data analysis and integration. Information supplementary to article by Michael Bittner, Jeffrey Trent and Paul Meltzer, Nature Genet. 22, 213-215. Duggan-DJ; Bittner-M; Cheh-Y; Meltzer-P Trent -JM (1999) Expression profiling using cDNA microarrays. Nature genet. 21 (1) supplement, 10-14. Eisen-M (1999) Cluster and Treeview manual, (ftp://rana.stanford.edu/software) Eisen-M, Spellman-PT, Brow^i-PO, Botsteiri-D. (1998) Cluster analysis and display of gene- wide expression patterns. Proc. Natl. Acad. Sci. USA 99:14863-14868 Golub-TR, Slonim-DK, Tairiayb-P, Huard-C, Gaasenbeek-M, Mesirov-JP, Coller-H, Loh-ML, Downing- JR, Caligiuri-MA, : Bloomfield-CD, and Lander-ES (1999) Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science Oct 15: 531-537 ' ■ '■■ : ' Lander-ES (1999) Array of hope. Nature genet 21 (1) supplement, 3-4. Tamayo-P, Slonim-D, Merirov-J, Zhu-Q, Kitareewan-S, Dmitrovsky-E, Lander-ES, and Golub-TR (1999) Interpreting patterns of gene expression with self-organizing maps: Methods and application to hematopoietic differentiation. Proc. Natl. Acad. Sci. USA 99:2907-2912 Romesburg-HC (1989) Cluster analysis for researchers. Chapter 16: How to make classifications. P203-216, Krieger Publishing Co. Malabar, Florida, USA Scherf-U, Ross-DT, Waltham- , Smith-LH, Lee-JK, Tanabe-L, Kohn-KW, Reinhold-WC, Myers-TG, Andrews-DT, Scudiero-DA, Eisen-MB, Sausville-EA, Pommier-Y, Botstein-D, -94- Brown-PO, Weinstein-JN (2000) A gene expression database for the molecular pharmacology of cancer. Nature genet. 24: 236-44; Slonim-DK, Tamayo-P, Mesirov-JP, Golub-TR, Lander-ES (1999) Class prediction and discovery using gene expression data. (ht ://www.genome.wi.mit.edu/MPR) Tavazole-S, Hughes- JD, Cariipbell-MJ, Cho-RJ, Church-GM. (1999) Systemic determination of genetic network architecture. Nature genet. 22: 281-285 Example 17. Identification of the bioresponses and the signature genes of an herbal composition. To further investigate the expression profiles induced by the. Huang Chin Tang, Jurkat T cells were treated with Huang Chih Tang (PHY906#2 obtained from Sun Ten Pharmaceutical Co. Taiwan) by 5 concentrations (1/20, 1/10, 1/5, 1/2.5, and 1 of IC 50 ). Nucleic acid microarray with dual-color detection was employed to measure the expression profiles. The mRNA extracted from herbal treated and untreated cells were labeled with Biotin- 16-dUTP and Dig-11-dUTP, respectively. A control group was established by using mRNA extracted from untreated cells and labeled the mRNA with biotin and dig respectively in equal proportions. Five sample groups with different concentrations of herbal treatment and one control group were employed for the study by the procedures described by Chen et al. (1998) with minor modifications. For data pre-processing, only array spots of high data quality were selected. The selection criteria were spots with signal to background ratio greater than 2.5-fold and more than 1.5-fold in differential expression ratio. By these criteria, 1044 genes were selected for further statistical analysis. The gene expression profiles ofthe control group were highly correlated with only 48 genes listed as statistical outliers that lie beyond the 2-fold differential expression range (Figure 25A). For the sample group, many differentially expressed genes are evident as illustrated in Figure 25B. The number of genes whose differential expression ratio greater than 2-fold increases with concentration of herbal treatment as shown in Figure 25C. These results prove that the identified differential expressed genes are truly inducted by Huang Chin Tang treatments. To identify the geries that were specifically induced by PHY906 (signature genes for PHY906), the expression profiles were clustered by a non-hierarchical cluster analysis programs "GeneCluster" developed ' by Tamayo et al., 1999. The computer program is based on self-organizing map (SOM) principle and the clusters of expression profiles are shown in -95- Figure 26. The X-axis represents the herbal concentration from low to high and the Y-axis is the gene-expression ratio.. The' signature genes were selected from the expression profiles which exhibit dosage response to the PHY906#2. The induced and repressed genes were selected from cluster 3 & 4 and cluster 18 & 19, respectively. In order to identify signature genes for PHY906, another batch of Huang Chin Tang, PHY906#3, with the same formula and manufacturing process were performed as described for PHY906#2. The induced and repressed genes commonly found/in both batches are shown in Figure 27. Score similarity of bioresponses by self-organizing map (SOM) To differentiate.herbal medicines of similar compositions, a scoring method is developed and the score S represents the difference in bioresponses of a biosystem to two different herbal compositions. S = Σ PøWø, ' - " ' " ■ . ■ " ■ Where P, is the number' of the common genes induced both by herbal prep. A and herbal prep. B in cluster i and in Clustery. For example, the SOM clustering results for the expression profiles of both batches of PHY906 are shown in Figure 28A. In cluster C13 and C14, 17 and 25 genes share the same expression, profiles for both batches of PHY906, respectively. In addition, 10 genes whose expression profiles induced by PHY906#2 are clustered in C13 but are clustered in C14 for PHY906#3. Therefore, P 7 ___=17, P W7 = 25, and P J J J ^=10. A weighing factor, W tj , describes the distance between the cluster i andj to indicate the similarity ofthe two expression profile clusters. In the case of C13 and C14, these 10 genes have similarly response to PHY906#2 and PHY906#3 (Figure 28B). The weighing factor is defined as: W y = 1 - E y /Max(E, ), where E, is the Euclidean distance between the cluster i andj and the value is normalized by Ey /Max(E y ). When i —j, W ; is 1. The number decreases as cluster i and clustery become more different (Figure 28C). Classification of 5 batches of herbal medicines To test how well the above method performs in classifying 5 batches of similar herbal preparations, Jurkat T cells were treated with 5 batches of Huang Chin Tang (PHY01040; #16, PHY010402; #17, PHY03061;',#Ϊ8, PHY03062; #19 and PHY02231; #20 obtained from Sun Ten Pharmaceutical Co.) by 5 concentrations (1, 1/2.5, 1/5, 1/10, and 1/20 of IC 50 ). The S ; - scores were calculated between pairs of herbal preparations in cluster analysis (Figure 29). The analysis programs, Kitsch Cluster was based on hierarchical clustering principle and was written by Dr. Joseph Fel≤en ' stein of Washington University (http://evolution.genetics.washington.edu/phylip.html). The S scores (distance) are tabulated -96- (Figure 29A), one can clearly identify that the shortest distance lie between batch #17 and batch #18 and that batch #17 is. similar to #18. Batch #16 also similar with batch #17 and #18 but batch #19 is dissimilar to the rest of batches. The results were confirmed by HPLC analyses. • ' ",-, ' Characterize an unknown herbal medicine based on the expression profiles To identify an unknown herbal medicine based on the characteristic expression profile database established as described as above, Jurkat T cells were treated with a tester sample #17 in 5 concentrations to set up the characteristic expression data set for the tester. The S score between the tester and each of the data sets (#16, #17, #18, #19 and #20) in the characteristic expression database were calculated and the S scores are: #16: 0.78, #17: 0.85, #18: 0.84, #19: 0.77, and #20: 0.79. These data show that the tester is most similar to #17 having a higher S score of 0.84. The exafnpl ' e deriioristrates that one can apply the method to identify an unknown herbal medicine based on the gene expression profiles induced by the herbal medicine in mammalian Cells. The identity of the. unknown herbal medicine can be inferred by aligning the characteristic expression profiles with a collection of characteristic expression profiles of herbal medicines in an HBR Array. The property of an herb can be described by four natures and five flavors (in Chinese Herbal Phramaceuticals. Ed. Zheng Hua Yen, People's Health publications, Beijing, China, 1997; Book of Ben Cao Gan IVlu by Shi Zeng Li, Ming Dynasty, China). Each ofthe four herbs in PHY906 may relate ;to anothef set of herbs with similar property (see Table 19). Or herbs with similar property rriay exhibit similar bioresponse. HBR Arrays may be used to determine or measure the relatedness hi terms ofthe property of herbs. Such information may be useful in creating a new herbal formulation. -97- Table 19: Herbs with PHY906 Properties Properties: Four natures-cold, hot;, warm, cool. Five flavors-acrid, bitter,, sweet, sour, bland. Example 18. Evaluation of an herbal medicines by HBR Array. As stated in Example! 6!, the component herbs of Huang Chin Tang are Scute, Peony, Licorice and Jujube. The gerie expression profiles induced by five batches of Huang Chin Tang in mammalian cells were characterized. A standard formula for Huang Chin Tang can be defined and characterized with animal studies or with clinical studies. For example, the #17 was used as the standard formula fof Huang Chin Tang based on the quality control and other standards set up by Sun Ten Pharmaceutical Co. The bioresponses of #17 were used to build the HBR Array for Huang Chin Tang. The marker genes in the HBR Array were selected to evaluate other preparations of Huang Chin Tang composition. A tester Huang Chin Tang may contain the same herb compositions but the component herbs may be grown under various environmental characteristics.. Comparing the bioresponses ofthe tester with the marker genes of standardized HBR Array,, the biological activities ofthe tester were evaluated. Furthermore, the marker 'genes whose expression levels are highly correlated (with |R|>0.99) with the dosage of component herbs in Huang Chin Tang (as stated in Example 16 and Figure 25) are selected for evaluation purpose. The tester Huang Chin Tang can be evaluated by comparing the specific bioresponses or expression levels ofthe selected set of marker genes with the HBR Array. If the expression levels or bioresponses ofthe selected marker genes are beyond ' the acceptable variation region, the amount or characteristics ofthe -98- component herbs are adjusted or modified to meet the acceptable variation. The process is repeated until the bioresponses induced by the revised herbal composition are within the acceptable variation range by comparing with the standard HBR Array. Example 19. Predicting biological activity and therapeutic applications of an herbal composition. According to the identified marker genes for PHY 906 (Figure 27), these genes can be used to predict the biological activities ofthe herbal composition. For example, the following underlined marker genes of PHY906 have been reported to involve in the following biological activities and therapeutic effects. The only effective drug against ALL is to inhibit the asparagine synthetase due to increased cellular apoptosis (Nandy et al, 1998). Long-acting drug somatostatin analogs are applied in the treatment of neurofibroma for their tumor growth inhibitory effect because they induce antiproliferative action mediated by the inhibition of G6PD, transketolase. or both (Boros et al., 1998). Ephrin-Al is a new melanoma growth factor and is highly expressed during melanoma progression (Easty et al., 1999). Mitogen- activated protein kinase (MAPK) family members have been recently reported to have opposing effects on apoptosis (Dabrowski et al., 2000). The expressions of asparagine synthetase, transketolase, ejphriri-Ai- and MAPK are repressed with the higher concentration of PHY906 treatments. The down-regulation of these genes are involved in cell apoptosis. The expression ofthe enzyme argininosuccinate synthetase, cathepsin G and chemokine RANTES are highly induced in inflammatory mechanism. By the PHY906 treatment the inflammatory involved genes are suppressed: . These literature reports provide a basis for predicting the biological activities or therapeutic effects of an herbal composition. References for Example 1.9. Nandy-P; Periclou-AP; Avramis-VI (1998) The synergism of 6-mercapopurine plus cytosine arabinoside followed by PEG-asparaginase in human leukemia cell lines (CCRF/CEM/0 and CCRF/CfiM/ara-C/7A) is due to increased cellular apoptosis. Anticancer Research 18: 727-737 Boros-LG; Brandes-JL; " Yusuf-FI; Cascante-M; Williams-RD; Schirmer-WJ (1998) Inhibition ofthe oxidative arid nonoxidative pentose phosphate pathways by somatostatin: a possible mechanism of antiturήor actiori. Medical Hypotheses 50: 501-506 Easty-DJ; Hill-SP; Hsu-MY; Fallowfield-ME; Florenes-VA; Herlyn-M; Benett-DC -99- (1999) Up-regulation of ephrin-Al during melanoma progression. Int. J. Cancer 84:494-501 Dabrowski-A; Tribillo-I; Dabrowska MI; Wereszczynska-SU; Gabryelewicz-A (2000) Activation of mitogen-activated protein kinases in different models of pancreatic acinar cell damage. Z-Gastroenterol. 38 : 469-481 The foregoing detaile description has been given for clearness of understanding only and no unnecessary limitations should be understood therefrom as modifications will be obvious to those skilled in the art. : ; • • ; While the invention has been described in connection with specific embodiments thereof, it will be understood that it: is capable of further modifications and this application is intended to cover any variations,' uses, or ' adaptations ofthe invention following, in general, the principles ofthe invention and including such departures from the present disclosure as come within nown or custori ary'pfactice Within the art to which the invention pertains and as may be applied to the 'essential features hereinbefore set forth and as follows in the scope ofthe appended claims.

Claims: -1,00- WHAT IS CLAIMED IS: ' . 1. A method of establishing a standardized Herbal BioResponse Array (HBR Array) for an herbal composition comprising: d) selecting a characterized herbal composition; e) exposing a biosystem to. a batch of the characterized herbal composition and collecting data on two or more markers, wherein one ofthe markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising: iv) producing a cell banking system; v) profiling the gene expression pattern of cells from the cell banking system before and after exposure to the herbal composition; vi) selecting as iriarkers those genes whose expression levels are changed by exposure to the herbal composition; " • ■" : • ">" . ' f) storing the marker data of step b) as a standardized HBR array. 2. The method of claim 1 , further comprising: g) repeating steps b) and c) for one more batches ofthe herbal composition using two or more ofthe same or differerit markers than used in step b); h) combining the HBR Arrays obtained in steps c) and d); and i) analyzing the combined HBR Array of step e) to generate a standardized HBR Array for the characterized herbal composition. -101- 3. The method of claims 1 or 2, wherein the characterized herbal composition has at least one known BioResponse. 4. The method of claims 1 or 2, wherein one or more ofthe following is known for the characterized herbal composition: chemical testing, the part ofthe plant used, the growing conditions of one or more ofthe individual herbs in the characterized herbal composition, the pre-harvest treatment of one or more ofthe individual herbs in the characterized herbal composition, the post-harvest treatment of one or more ofthe individual herbs in the characterized herbal composition, the post-harvest treatment ofthe characterized herbal composition, and the relative proportions ofthe individual herbs in the herbal composition. 5. The method of claims 1. or 2, wherein the cell banking system comprises a master cell bank and a working cell bank. 6. The method of claim 5, wherein the cells ofthe working cell bank are obtained from the master cell bank. ' 7. The method of claim 5, wherein the step of profiling the gene expression pattern of cells from the cell banking system before and after exposure to the herbal composition is performed using cells from the working cell bank. 8. The method of claims 1 or . 2, wherein the change in gene expression is determined using a nucleic acid microarray. : , ' 9. The method of claim 8, wherein the said genes whose expression levels are changed by exposure to the herbal coriiposition are selected based on the criteria of having a signal to noise ratio of about 2.5 or greater in the nucleic acid microarray and having an about 1.5 or greater change in the differential expression ratio. 10. The method of claim 8, • wherein data regarding between about 10 and about 20,000 genes whose expression levels are changed is stored as part ofthe HBR Array. -102- 11. The method of claim 1.0, wherein data between about 10 and about 1,500 genes whose expression levels are changed is stored as part ofthe HBR Array. 12. A method of evaluating an herbal composition comprising: a) exposing a biosystem to a batch of the herbal composition and collecting data on two or more markers, wherein one ofthe markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising: i) producing a cell banking system; ii) profiling the gene expression pattern of cells from the cell banking system before and aftef exposure to the herbal composition; iii) selecting as ' markers those genes whose expression levels are changed by exposure to. the herbal composition; b) comparing the collected'friarker data with a standardized HBR Array for the same or a substantially same herbal composition as that ofthe batch herbal composition, wherein the standardized HBR Array contains one ofthe markers data on gene expression. " ' ' • 13. A method of determining if an herbal composition meets a standard specification comprising: a) exposing a biosystem to a batch ofthe herbal composition and collecting data on two or more markers, wherein one of the markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising: i) producing a' cell banking system; ii) profiling the gerie expression pattern of cells the cell banking system before and after exposure to the herbal composition; -103- iii) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition; b) comparing the collected marker data with a standardized HBR Array for the same or a substantially same herbal composition as that ofthe batch herbal composition, wherein the standardized HBR Array contains as one of the markers data on gene expression; and c) determining which herbal compositions have marker data that is similar to that ofthe standardized HBR Array within an acceptable level. 14. The method of claim 113, wherein said determining which herbal compositions have marker data which is similar to that ofthe standardized HBR Array within an acceptable level is determined quantitatively or qualitatively. 15. The method of claims 13 or 14, wherein the standardized HBR Array includes an acceptable range of variation for each marker. 16. A method of adjusting the components of a herbal composition so that it meets a standard specification for the same or substantially the same herbal composition, comprising: , ' " " . ' ' ■' a) exposing a biosystem to a batch ofthe herbal composition and collecting data on two or more markers, wherein one ofthe markers is a change in gene expression determined through the use of a nucleic acid : miCroarray, produced by the steps comprising: i) producing a cell banking system; ii) profiling the ' gene expression pattern of cells the cell banking system before and after exposure to the hefbal composition; iii) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition; -104- b) comparing the collected marker data with a standardized HBR Array for the same or a substantially same herbal composition as that ofthe batch herbal composition, wherein the standardized HBR Array contains as one of markers data on gene expression, and wherein the standardized HBR Array also includes an acceptable range of variation for each marker; c) determining whether the herbal composition has marker data that is within the acceptable level of variation for the standardized HBR Array; and d) if the marker data is not within the acceptable level of variation for the standardized HBR Array, adjusting the components ofthe herbal composition. 17. The method of claim.16, whefein steps (a) through (d) are repeated until the marker data of the herbal composition is within the acceptable level of variation ofthe standardized HBR Array. ' . ; 18. A method of changing the components of a herbal composition so that it meets a standard specification of another herbal composition, comprising: a) exposing a biosystem to a batch ofthe herbal composition and collecting data on two or more markers, wherein one' of the markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising: i) producing a cei bankhig system; ii) profiling.the gene expression pattern of cells the cell banking system before and after exposure to the herbal composition; iii) selecting as riiarkers those genes whose expression levels are changed by exposure to the herbal composition; b) comparing the Collected, marker data with a standardized HBR Array for the other herbal compositions as that of the. batch herbal composition, wherein the standardized HBR Array -105- contains as one of markefs data on gene expression, and wherein the standardized HBR Array also includes an acceptable range of variation for each marker; c) determining whether the herbal composition has marker data that is within the acceptable level of variation for the standardized HBR Array; and d) if the marker data is not within the acceptable level of variation for the standardized HBR Array, changing the components ofthe herbal composition. 19. The method of claim 18, wherein steps a) through d) are repeated until the marker data of the herbal composition is within the acceptable level of variation ofthe standardized HBR Array. ' ." . : - • " ' • ' 20. A method for predicting the biological activity of an herbal composition comprising: a) exposing a biosystem to a batch ofthe herbal composition and measuring the differential responses of two or more markers, wherein one ofthe markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising; • . ' • i) producing a cell banking system; ii) profiling the gene expression pattern of cells the cell banking system before and after exposure to the herbal composition; iii) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition; wherein the set of differential response measurements constitute an Herbal BioResponse Array (HBR Array) data set; -106- b) comparing the HBR Array of the batch herbal composition to at least one previously- obtained HBR Array of a characterized herbal composition, wherein the previously- obtained HBR Array contains as one ofthe markers data on gene expression; and c) predicting the biological activity of the batch herbal composition based on the HBR Array comparison made in step b). 21. A method of measuring the relatedness of a herbal composition to a characterized herbal composition comprising: a) exposing a biosystem to a. batch ofthe herbal composition and measuring the differential responses of two or mofe markers, whefein one ofthe markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising; i) producing a cell banking system; ii) profiling the gene expression pattern of cells the cell banking system before and after exposure to the herbal composition; i iii) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition; wherein the set of differential response measurements constitute an Herbal BioResponse Array (HBR Array) data set; b) comparing the HBR Array ofthe batch herbal composition to at least one previously- obtained HBR Array of a characterized herbal composition, wherein the previously- obtained HBR Array contains as one ofthe markers data on gene expression; and c) determining the relatedness ofthe herbal composition to the characterized herbal composition based on the HBR Array comparison made in step b). -107- 22. A method for predicting new therapeutic applications of an herbal comprising: a) exposing a biosystem to a batch ofthe herbal composition and measuring the differential responses of two or more markers, wherein one ofthe markers is a change in gene expression determined through the use of a nucleic acid microarray, produced by the steps comprising; i) producing a cell banking system; ii) profiling the gene expression pattern of cells the cell banking system before and after exposure to the herbal composition; iii) selecting as markers those genes whose expression levels are changed by exposure to tile herbal composition; wherein the set of diffeferitial response measurements constitute an Herbal BioResponse Array (HBR Array) data set; : •■■: b) predicting the new therapeutic applications based on the predicted biological activity ofthe markers in the HBR Array. 23. A method for determining the gene expression profile induced by individual chemical entities in an herbal comppsition comprising: a) producing a cell banking system; b) profiling the gene expression pattern of cells from the cell banking system before and after exposure to the herbal composition; -108- c) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition and placing into an HBR Array; d) comparing the HBR Array generated in step (c) with a standardized HBR Array for a similar or modified herbal composition; e) determining the relative amourits ofthe individual chemical entities ofthe herbal composition; and ; V. ' . '■ • • f) comparing the ariiόunt of the individual chemical entities to the result of step (b) to identify those genes whose expression levels change as the amount ofthe individual chemical entity in the herbal composition changes. 24. A method for determining the gene expression profile induced by individual chemical entities in a complex mixture without extracting the chemicals from the complex mixture such as an herbal composition, comprising a) producing a cell banking system; b) profiling the gene expfessiόn pattern of cells from the cell banking system before and after exposure to the herbal composition; c) selecting as markers those genes whose expression levels are changed by exposure to the herbal composition; d) comparing the collected marker data with a standardized HBR Array for a substantially same or modified herbal composition; -109- e) characterizing the chemical components ofthe said herbal compositions; f) comparing the identified chemical compositions to identify the differential levels of individual chemical components in herbal compositions; g) correlating the differential chemical cohiponent amounts with the differential bioresponses in the HBR Array to identify the characteristic bioresponses for each chemical entity.

Inventors: CHENG YUNG-CHI, KUNG PATRICK C, LEE YUN-SHIEN, PECK KONAN, SHER YUH-PYNG, CHENG, YUNG-CHI, KUNG, PATRICK, C., LEE, YUN-SHIEN, PECK, KONAN, SHER, YUH-PYNG
Applicants: CHENG YUNG CHI, KUNG PATRICK C, LEE YUN SHIEN, PECK KONAN, SHER YUH PYNG, UNIV YALE, CHENG, YUNG-CHI, KUNG, PATRICK, C., LEE, YUN-SHIEN, PECK, KONAN, SHER, YUH-PYNG, YALE UNIVERSITY
Application year: 2001
Main classification: C12Q1
IPC codes: C12N15, C12Q1, G01N33, G01N37, G06F17