Activation of a nuclear factor kappaB/interleukin-1 positive feedback loop by amphiregulin in human breast cancer cells.

We have recently shown that an amphiregulin-mediated autocrine loop is responsible for growth factor-independent proliferation, motility, and invasive capacity of some aggressive breast cancer cells, such as the SUM149 breast cancer cell line. In the present study, we investigated the mechanisms by which amphiregulin activation of the epidermal growth factor receptor (EGFR) regulates these altered phenotypes. Bioinformatic analysis of gene expression networks regulated by amphiregulin implicated interleukin-1alpha (IL-1alpha) and IL-1beta as key mediators of amphiregulin's biological effects. The bioinformatic data were validated in experiments which showed that amphiregulin, but not epidermal growth factor, results in transcriptional up-regulation of IL-1alpha and IL-1beta. Both IL-1alpha and IL-1beta are synthesized and secreted by SUM149 breast cancer cells, as well as MCF10A cells engineered to express amphiregulin or MCF10A cells cultured in the presence of amphiregulin. Furthermore, EGFR, activated by amphiregulin but not epidermal growth factor, results in the prompt activation of the transcription factor nuclear factor-kappaB (NF-kappaB), which is required for transcriptional activation of IL-1. Once synthesized and secreted from the cells, IL-1 further activates NF-kappaB, and inhibition of IL-1 with the IL-1 receptor antagonist results in loss of NF-kappaB DNA binding activity and inhibition of cell proliferation. However, SUM149 cells can proliferate in the presence of IL-1 when EGFR activity is inhibited. Thus, in aggressive breast cancer cells, such as the SUM149 cells, or in normal human mammary epithelial cells growing in the presence of amphiregulin, EGFR signaling is integrated with NF-kappaB activation and IL-1 synthesis, which cooperate to regulate the growth and invasive capacity of the cells.

Structure alignment based on coding of local geometric measures.

BACKGROUND: A structure alignment method based on a local geometric property is presented and its performance is tested in pairwise and multiple structure alignments. In this approach, the writhing number, a quantity originating from integral formulas of Vassiliev knot invariants, is used as a local geometric measure. This measure is used in a sliding window to calculate the local writhe down the length of the protein chain. By encoding the distribution of writhing numbers across all the structures in the protein databank (PDB), protein geometries are represented in a 20-letter alphabet. This encoding transforms the structure alignment problem into a sequence alignment problem and allows the well-established algorithms of sequence alignment to be employed. Such geometric alignments offer distinct advantages over structural alignments in Cartesian coordinates as it better handles structural subtleties associated with slight twists and bends that distort one structure relative to another. RESULTS: The performance of programs for pairwise local alignment (TLOCAL) and multiple alignment (TCLUSTALW) are readily adapted from existing code for Smith-Waterman pairwise alignment and for multiple sequence alignment using CLUSTALW. The alignment algorithms employed a blocked scoring matrix (TBLOSUM) generated using the frequency of changes in the geometric alphabet of a block of protein structures. TLOCAL was tested on a set of 10 difficult proteins and found to give high quality alignments that compare favorably to those generated by existing pairwise alignment programs. A set of protein comparison involving hinged structures was also analyzed and TLOCAL was seen to compare favorably to other alignment methods. TCLUSTALW was tested on a family of protein kinases and reveal conserved regions similar to those previously identified by a hand alignment. CONCLUSION: These results show that the encoding of the writhing number as a geometric measure allow high quality structure alignments to be generated using standard algorithms of sequence alignment. This approach provides computationally efficient algorithms that allow fast database searching and multiple structure alignment. Because the geometric measure can employ different window sizes, the method allows the exploration of alignments on different, well-defined length scales.

From microarray to biological networks: Analysis of gene expression profiles.

Powerful new methods, such as expression profiles using cDNA arrays, have been used to monitor changes in gene expression levels as a result of a variety of metabolic, xenobiotic, or pathogenic challenges. This potentially vast quantity of data enables, in principle, the dissection of the complex genetic networks that control the patterns and rhythms of gene expression in the cell. Here we present a general approach to developing dynamic models for analyzing time series of whole-genome expression. The parameters in the model show the influence of one gene expression level on another and are calculated using singular value decomposition as a means of inverting noisy and near-singular matrices. Correlative networks can then be generated based on these parameters with a simple threshold approach. We also demonstrate how dynamic models can be used in conjunction with cluster analysis to analyze microarray time series. Using the parameters from the dynamic model as a metric, two-way hierarchical clustering could be performed to visualize how influencing genes affect the expression levels of responding genes. Application of these approaches is demonstrated using gene expression data in yeast cell cycle.

Evaluation of soluble CD30 as an immunologic marker in heart transplant recipients.

CD30 is an immunologic molecule that belongs to the TNF-R superfamily. CD30 serves as a T-cell signal transducing molecule that is expressed by a subset of activated T lymphocytes, CD45RO+ memory T cells. Augmentation of soluble CD30 during kidney transplant rejection has been reported. Our study sought to determine whether the level of sCD30 prior to heart transplant could categorize patients into high versus low immunologic risk for a poor outcome. A significant correlation was observed between high levels of soluble CD30 and a reduced incidence of infection. None of the 35 patients with high pretransplant levels of sCD30 level (>90 U/mL) developed infections posttransplantation. However, 9 of 65 patients who had low levels of sCD30 (<90 U/mL) developed infections posttransplantation (P < .02). No remarkable differences were noted among the other clinical parameters. The results also showed that the high-definition flow-bead (HDB) assay detected both weak and strong class I and class II HLA antibodies, some of which (weak class II HLA Abs) were undetectable by the anti-human globulin cytotoxicity method. In addition, more antibody specificities were detected by HDB. In conclusion, we have observed that high levels of sCD30 prior to heart transplant may be associated with greater immunologic ability and therefore produce a protective effect on the development of infection post heart transplant. We have also shown that the HDB assay is superior to the visual cytotoxicity method to detect HLA antibodies, especially those to class II HLA antigens.

Left mini-thoracotomy for beating heart bypass grafting: a safe alternative to high-risk intervention for selected grafting of the circumflex artery distribution.

BACKGROUND: Progression of disease and bypass graft attrition results in a population of patients who require repeated coronary interventions. Frequently, these patients have patent internal mammary artery grafts and require isolated intervention to the circumflex distribution. As an alternative to high-risk repeated sternotomy and conventional bypass surgery or catheter-based intervention, the circumflex marginal vessels may be approached by thoracotomy. We reviewed our experience in revascularizing the circumflex distribution with off-pump techniques via left mini-thoracotomy. METHODS AND RESULTS: Thirty-two patients underwent off-pump bypass grafting of the circumflex vessels via thoracotomy from December 1995 to April 2000. Twenty-seven patients presented with circumflex disease after having previous bypass grafting. Five patients, who presented with circumflex disease and either nondiseased or ungraftable disease in their other arteries, were revascularized as a primary procedure. There was no observed mortality. Seven patients (22%) required inotropes on leaving the operating room, and 3 patients (9.4%) received transfusion of packed red blood cells. There was 1 reoperation for bleeding and 1 patient with a postoperative neurological deficit. There were no perioperative myocardial infarctions. The average length of stay was 4.8 days from time of surgery to discharge. CONCLUSIONS: Off-pump grafting via thoracotomy provides a safe and effective alternative approach for patients requiring limited revascularization. Potential cardiac injury and danger to viable grafts from repeated sternotomy is minimized, and manipulation of the diseased ascending aorta is avoided. Morbidity, hospital length of stay, and cost are less than for conventional repeated coronary bypass surgery.

Improved survival rates support left ventricular assist device implantation early after myocardial infarction.

OBJECTIVES: Implantation of left ventricular assist devices (LVADs) early after acute myocardial infarction (MI) has traditionally been thought to be associated with high mortality rates due to technical limitations and severe end-organ dysfunction. At some experienced centers, doctors have refrained from earlier operation after MI to allow for a period of hemodynamic and end-organ stabilization. METHODS: We retrospectively investigated the effect of preoperative MI on the survival rates of 25 patients who received a Thermocardiosystems Incorporated LVAD either <2 weeks (Early) (n = 15) or >2 weeks (Late) (n = 10) after MI. Outcome variables included perioperative right ventricular assistance (and right-sided circulatory failure), hemodynamic indexes, percent transplanted or explanted, and mortality. RESULTS: No statistically significant differences were demonstrated between demographic, perioperative or hemodynamic variables between the Early and Late groups. Patients in the Early group demonstrated a lower rate of perioperative mechanical right ventricular assistance, but had a higher rate of perioperative inhaled nitric oxide use. In addition, 67% of patients in the Early group survived to transplantation and 7% to explantation, findings comparable to those in the Late group (60% and 0% respectively). CONCLUSIONS: This clinical experience suggests that patients may have comparable outcomes whether implanted early or late after acute MI. These data therefore support the early identification and timely application of this modality in post-MI LVAD candidates, as this strategy may also reveal a subgroup of patients for whom post-MI temporary LVAD insertion may allow for full ventricular recovery.

A random polymer model of the statistical significance of structure alignment.

A theory for assessing the statistical significance of structure alignment is developed using a random or Gaussian chain model. In this model, we consider the statistical distribution of the root mean square distance (rmsd) of the alignment between two random chains of equal length and common center of mass (referred to as Case 1). We demonstrate that the rmsd2 is distributed as a sum of independent Gamma variables. Analytic results on the mean and variance of the rmsd2 are presented. Since rmsd is strongly dependent on the length, we define the dimensionless quantity, reduced rmsd, as the rmsd divided by the radius of gyration. We find that the reduced rmsd can be accurately approximated by an extreme value distribution (EVD) that is independent of chain length and of bond length. The parameters of the EVD can be calculated from the mean and the variance of the rmsd2. We also consider the case of two chains with a common center of mass that are then rotated to minimize the rmsd (Case 2). In this case, the distribution of reduced rmsd can again be accurately approximated by an EVD, which is independent of the chain length and expected bond length. This distribution is used to calculate the p-value for a given reduced rmsd. Performing an analogous comparison for proteins, we find that approximately M(nu) and nu = 0.28 and 0.32 for Case 1 and Case 2, respectively, where M is the chain length. This result for Case 2 exactly matches with previous scaling results and suggests that rmsd/M(nu)is an appropriate metric for protein structure alignment and will be independent of chain length. We also find that the new score roughly follows the EVD.

From microarrays to networks: mining expression time series.

Over the past few years, powerful new methods have been devised that enable researchers to study the expression dynamics of many genes simultaneously (e.g. gene expression profiles using cDNA microarrays). In principle, this potentially vast quantity of data enables the dissection of the complex genetic networks that control the patterns and rhythms of gene expression in the cell. Finding the patterns in those data represents the next major phase in our understanding of the programming and functioning of the living cell. Simple dynamic models can be used to generate gene expression networks. These networks reveal the phenomenological link between the expression of different genes. This review discuss how these networks are generated and outlines several data-mining techniques for extracting relationships and hypotheses in gene expression. These emerging methods can be applied to a range of biological problems.

A sequence alignment algorithm with an arbitrary gap penalty function.

An algorithm for aligning biological sequences is presented that is an adaptation of the sequence generating function approach used in the statistical mechanics of biopolymers. This algorithm uses recursion relationships developed from a partition function formalism of alignment probabilities. It is implemented within a dynamic programming format that closely resembles the forward algorithm used in hidden Markov models (HMM). The algorithm aligns sequences or structures according to the statistically dominant alignment path and will be referred to as the SDP algorithm. An advantage of this method over previous ones is that it allows more complicated and physically realistic gap penalty functions to be incorporated into the algorithm in a facile manner. The performance of this algorithm in a case study of aligning the heavy and light chain from the variable region of an immunoglobulin is investigated.

Duplication models for biological networks.

Are biological networks different from other large complex networks? Both large biological and nonbiological networks exhibit power-law graphs (number of nodes with degree k, N(k) approximately k(-beta)), yet the exponents, beta, fall into different ranges. This may be because duplication of the information in the genome is a dominant evolutionary force in shaping biological networks (like gene regulatory networks and protein-protein interaction networks) and is fundamentally different from the mechanisms thought to dominate the growth of most nonbiological networks (such as the Internet). The preferential choice models used for nonbiological networks like web graphs can only produce power-law graphs with exponents greater than 2. We use combinatorial probabilistic methods to examine the evolution of graphs by node duplication processes and derive exact analytical relationships between the exponent of the power law and the parameters of the model. Both full duplication of nodes (with all their connections) as well as partial duplication (with only some connections) are analyzed. We demonstrate that partial duplication can produce power-law graphs with exponents less than 2, consistent with current data on biological networks. The power-law exponent for large graphs depends only on the growth process, not on the starting graph.

A new scoring function and associated statistical significance for structure alignment by CE.

A new scoring function for assessing the statistical significance of protein structure alignment has been developed. The new scores were tested empirically using the combinatorial extension (CE) algorithm. The significance of a given score was given a p-value by curve-fitting the distribution of the scores generated by a random comparison of proteins taken from the PDB_SELECT database and the structural classification of proteins (SCOP) database. Although the scoring function was developed based on the CE algorithm, it is portable to any other protein structure alignment algorithm. The new scoring function is examined by sensitivity, specificity, and ROC curves.

Cluster analysis of dynamic parameters of gene expression.

Cluster analysis has proven to be a valuable statistical method for analyzing whole genome expression data. Although clustering methods have great utility, they do represent a lower level statistical analysis that is not directly tied to a specific model. To extend such methods and to allow for more sophisticated lines of inference, we use cluster analysis in conjunction with a specific model of gene expression dynamics. This model provides phenomenological dynamic parameters on both linear and non-linear responses of the system. This analysis determines the parameters of two different transition matrices (linear and nonlinear) that describe the influence of one gene expression level on another. Using yeast cell cycle microarray data as test set, we calculated the transition matrices and used these dynamic parameters as a metric for cluster analysis. Hierarchical cluster analysis of this transition matrix reveals how a set of genes influence the expression of other genes activated during different cell cycle phases. Most strikingly, genes in different stages of cell cycle preferentially activate or inactivate genes in other stages of cell cycle, and this relationship can be readily visualized in a two-way clustering image. The observation is prior to any knowledge of the chronological characteristics of the cell cycle process. This method shows the utility of using model parameters as a metric in cluster analysis.

Alternative technique of right-sided outflow cannula insertion for right ventricular support.

Right ventricular assist devices are an important part of the armamentarium of cardiac surgeons for the treatment of right-sided circulatory failure after cardiac transplantation or insertion of a left ventricular assist device. However, right ventricular assist device insertion can be technically challenging in the setting of pulmonary hypertension because of a number of concomitant anatomic and physiologic phenomena. We present a technique for the insertion of the right ventricular assist device outflow cannula that is easier and faster to insert, and safer to explant, especially if cardiopulmonary bypass is to be avoided.

Off-pump bypass grafting is safe in patients with left main coronary disease.

BACKGROUND: Because of a concern about the ability to tolerate beating heart grafting, patients with left main coronary artery stenosis have been excluded from off-pump bypass. We reviewed our experience with off-pump coronary artery bypass grafting for patients with left main coronary artery disease. METHODS: Eight hundred twenty-three patients underwent bypass grafting for left main coronary artery disease from January 1998 to October 1999. One hundred patients were revascularized without the use of cardiopulmonary bypass and compared with a contemporaneous cohort of 723 patients who underwent grafting with the aid of cardiopulmonary bypass. All patients had multivessel grafting performed through a sternotomy. RESULTS: There was one death (1%) in the group undergoing off-pump grafting as compared with a 30-day mortality of 4.7% (p = 0.059) in the on-pump group. Univariate analysis established that patients revascularized without cardiopulmonary bypass were significantly less likely to require postoperative inotropic support (23% versus 62%, p < 0.001) and transfusion (35% versus 67%, p < 0.001). Logistic regression analysis revealed that cardiopulmonary bypass was an independent risk factor for mortality (odds ratio, 7.3; 95% confidence interval, 1.34 to 138.4). CONCLUSIONS: Coronary artery bypass grafting using off-pump techniques are safe and effective in left main coronary artery disease.

Influence of diabetes on mortality and morbidity: off-pump coronary artery bypass grafting versus coronary artery bypass grafting with cardiopulmonary bypass.

BACKGROUND: Myocardial revascularization in diabetic patients is challenging with no established optimum treatment strategy. We reviewed our coronary artery bypass grafting experience to determine the impact of eliminating cardiopulmonary bypass on outcomes in diabetic patients relative to nondiabetic patients. METHODS: From January 1995 through December 1999, 9,965 patients, of whom 2,891 (29%) had diabetes, underwent isolated coronary artery bypass grafting. Diabetic and nondiabetic patients were further divided into groups on the basis of cardiopulmonary bypass use. Twelve percent (346 of 2,891) of diabetic patients and 12% (829 of 7,074) of nondiabetic patients underwent coronary artery bypass grafting without cardiopulmonary bypass; the remainder had coronary artery bypass grafting with cardiopulmonary bypass. Nineteen preoperative variables were compared among treatment groups by univariate analysis. RESULTS: Patients undergoing coronary artery bypass grafting without cardiopulmonary bypass compared with those having coronary artery bypass grafting with cardiopulmonary bypass had higher mean predicted mortalities (diabetic, 3.96% versus 3.72%, p = 0.83; nondiabetic, 3.03% versus 2.86%, p = 0.79). In nondiabetic patients, coronary artery bypass grafting without cardiopulmonary bypass provides an actual and risk-adjusted survival advantage over coronary artery bypass grafting with cardiopulmonary bypass (1.81% versus 3.44%, p = 0.0127; risk-adjusted mortality, 1.79% versus 3.61%, p = 0.007). This survival benefit of coronary artery bypass grafting without cardiopulmonary bypass was not seen in diabetic patients (2.89% versus 3.69%, p = 0.452; risk-adjusted mortality, 2.19% versus 2.98%, p = 0.42). Diabetic patients undergoing coronary artery bypass grafting without cardiopulmonary bypass had fewer complications, including decreased blood product use (34.39% versus 58.4%, p = 0.001), and reduced incidence of prolonged ventilation (6.94% versus 12.10%, p = 0.005), atrial fibrillation (15.90% versus 23.26%, p = 0.002), and renal failure requiring dialysis (0.87% versus 2.75%, p = 0.036). CONCLUSIONS: The survival advantage in nondiabetic patients treated with coronary artery bypass grafting without cardiopulmonary bypass is not apparent in diabetic patients. Coronary artery bypass grafting without cardiopulmonary bypass in diabetic patients is nevertheless associated with a significant reduction in morbidity.

Beating heart surgery reduces mortality in the reoperative bypass patient.

BACKGROUND: Reoperative coronary artery bypass grafting (redo-CABG) has an increased operative morbidity and mortality compared to patients undergoing primary revascularization. In an effort to reduce the hazards of reoperative CABG, we commenced revascularizing selected patients without cardiopulmonary bypass (CPB) as an alternative to conventional approaches. METHODS: From January 1998 to Dec. 2000, 432 patients underwent reoperative CABG, 153 patients (35%) without the aid of CPB. Treatment groups were compared by means of univariate analysis for preoperative risk factors and postoperative complications. Predicted risk and risk-adjusted mortality were determined by the Society of Thoracic Surgeons risk algorithm. RESULTS: There was a significant difference in the preoperative predicted risk scores between the two treatment groups (off pump 6.5% vs. on pump 5.4%, p=0.0343). There was a significant difference in the off pump observed mortality (2.61%) versus the on pump group (9.68%, p=0.0065). Decreased morbidity in the off pump group was evidenced by a reduced need for blood products (25% vs. 67%, p<0.0001), and the incidence of prolonged ventilation (4% vs. 14%, p=0.0032). The off pump group also had shorter hospital stays (6.2 +/- 5.96 days vs. 8.0 +/- 7.82, p=0.0091). No significant differences between the two groups were seen in the prevalence of perioperative myocardial infarction, stroke, renal failure, or reoperation for bleeding. CONCLUSION: Bypass grafting without CPB significantly decreases mortality and morbidity in selected reoperative patients, and should be considered a viable alternative to conventional approaches.

Beta-endorphins and mood states during resistance exercise.

Research examining the relationship between beta-endorphin immunoreactivity and affect has resulted in equivocal findings. To examine this relationship further, 10 male and 10 female college-age students participated in a controlled free-weight exercise session. Blood plasma samples were procured both prior to and following the exercise bout and beta-endorphin levels were determined by immunoassay. Analysis indicated that plasma beta-endorphin levels following exercise were significantly decreased from pre-exercise levels. Pearson correlations showed no significant relationship between pre- or postexercise plasma beta-endorphins and either total mood disturbance or Profile of Mood States subscores. The data support previous failure of resistance exercise to produce an increase in beta-endorphin immunoreactivity as well as the lack of a significant relationship between affect and the beta-endorphin response to exercise.

Oncogene activation induces metabolic transformation resulting in insulin-independence in human breast cancer cells.

Normal breast epithelial cells require insulin and EGF for growth in serum-free media. We previously demonstrated that over expression of breast cancer oncogenes transforms MCF10A cells to an insulin-independent phenotype. Additionally, most breast cancer cell lines are insulin-independent for growth. In this study, we investigated the mechanism by which oncogene over expression transforms MCF10A cells to an insulin-independent phenotype. Analysis of the effects of various concentrations of insulin and/or IGF-I on proliferation of MCF10A cells demonstrated that some of the effects of insulin were independent from those of IGF-I, suggesting that oncogene over expression drives a true insulin-independent proliferative phenotype. To test this hypothesis, we examined metabolic functions of insulin signaling in insulin-dependent and insulin-independent cells. HER2 over expression in MCF10A cells resulted in glucose uptake in the absence of insulin at a rate equal to insulin-induced glucose uptake in non-transduced cells. We found that a diverse set of oncogenes induced the same result. To gain insight into how HER2 oncogene signaling affected increased insulin-independent glucose uptake we compared HER2-regulated gene expression signatures in MCF10A and HER2 over expressing MCF10A cells by differential analysis of time series gene expression data from cells treated with a HER2 inhibitor. This analysis identified genes specifically regulated by the HER2 oncogene, including VAMP8 and PHGDH, which have known functions in glucose uptake and processing of glycolytic intermediates, respectively. Moreover, these genes specifically implicated in HER2 oncogene-driven transformation are commonly altered in human breast cancer cells. These results highlight the diversity of oncogene effects on cell regulatory pathways and the importance of oncogene-driven metabolic transformation in breast cancer.

Toward a classification of isodynamic feed-forward motifs.

A preceding study analysed how the topology of network motifs affects the overall rate of the underlying biochemical processes. Surprisingly, it was shown that topologically non-isomorphic motifs can still be isodynamic in the sense that they exhibit the exact same performance rate. Because of the high prevalence of feed-forward functional modules in biological networks, one may hypothesize that evolution tends to favour motifs with faster dynamics. As a step towards ranking the efficiency of feed-forward network motifs, we use a linear flow model to prove theorems establishing that certain classes of motifs are isodynamic. In partitioning the class of all motifs on n nodes into equivalence classes based upon their dynamics, we establish a basis for comparing the efficiency/performance rates of different motifs. The potential biological importance of the theorems is briefly discussed and is the subject of an ongoing large-scale project.