A modified gambler's ruin model of polyethylene chains in the amorphous region.

Polyethylene chains in the amorphous region between two crystalline lamellae M unit apart are modeled as random walks with one-step memory on a cubic lattice between two absorbing boundaries. These walks avoid the two preceding steps, though they are not true self-avoiding walks. Systems of difference equations are introduced to calculate the statistics of the restricted random walks. They yield that the fraction of loops is (2M - 2)/(2M + 1), the fraction of ties 3/(2M + 1), the average length of loops 2M - 0.5, the average length of ties 2/3M2 + 2/3M - 4/3, the average length of walks equals 3M - 3, the variance of the loop length 16/15M3 + O(M2), the variance of the tie length 28/45M4 + O(M3), and the variance of the walk length 2M3 + O(M2).

Fos and Jun bend the AP-1 site: effects of probe geometry on the detection of protein-induced DNA bending.

The effect of Fos and Jun binding on the structure of the AP-1 recognition site is controversial. Results from phasing analysis and phase-sensitive detection studies of DNA bending by Fos and Jun have led to opposite conclusions. The differences between these assays, the length of the spacer between two bends and the length of the sequences flanking the bends, are investigated here using intrinsic DNA bend standards. Both an increase in the spacer length as well as a decrease in the length of flanking sequences resulted in a reduction in the phase-dependent variation in electrophoretic mobilities. Probes with a wide separation between the bends and short flanking sequences, such as those used in the phase-sensitive detection studies, displayed no phase-dependent mobility variation. This shape-dependent variation in electrophoretic mobilities was reproduced by complexes formed by truncated Fos and Jun. Results from ligase-catalyzed cyclization experiments have been interpreted to indicate the absence of DNA bending in the Fos-Jun-AP-1 complex. However, truncated Fos and Jun can alter the relative rates of inter- and intramolecular ligation through mechanisms unrelated to DNA bending, confounding the interpretation of cyclization data. The analogous phase- and shape-dependence of the electrophoretic mobilities of the Fos-Jun-AP-1 complex and an intrinsic DNA bend confirm that Fos and Jun bend DNA, which may contribute to their functions in transcription regulation.

Inhibition of neoplastic development in the liver by hepatocyte growth factor in a transgenic mouse model.

Overexpression of the c-myc oncogene is associated with a variety of both human and experimental tumors, and cooperation of other oncogenes and growth factors with the myc family are critical in the evolution of the malignant phenotype. The interaction of hepatocyte growth factor (HGF) with c-myc during hepatocarcinogenesis in a transgenic mouse model has been analyzed. While sustained overexpression of c-myc in the liver leads to cancer, coexpression of HGF and c-myc in the liver delayed the appearance of preneoplastic lesions and prevented malignant conversion. Furthermore, tumor promotion by phenobarbital was completely inhibited in the c-myc/HGF double transgenic mice, whereas phenobarbital was an effective tumor promoter in the c-myc single transgenic mice. The results indicate that HGF may function as a tumor suppressor during early stages of liver carcinogenesis, and suggest the possibility of therapeutic application for this cytokine.

Long-lasting reduction of inhibitory function and gamma-aminobutyric acid type A receptor subunit mRNA expression in a model of temporal lobe epilepsy.

This study evaluated hippocampal inhibitory function and the level of expression of gamma-aminobutyric acid type A (GABAA) receptor mRNA in an in vivo model of epilepsy. Chronic recurrent limbic seizures were induced in rats using injections of pilocarpine. Electrophysiological studies performed on hippocampal slices prepared from control and epileptic animals 1 to 2 months after pilocarpine injections demonstrated a significant hyperexcitability in the epileptic animals. Reduced levels of mRNA expression for the alpha 2 and alpha 5 subunits of the GABAA receptors were evident in the CA1, CA2, and CA3 regions of the hippocampus of epileptic animals. No decrease in mRNA encoding alpha 1, beta 2, or gamma 2 GABAA receptor subunits was observed. In addition, no change in the mRNA levels of alpha CaM kinase II was seen. Selective decreases in mRNA expression did not correlate with neuronal cell loss. The results indicate that selective, long-lasting reduction of GABAA subunit mRNA expression and increased excitability, possibly reflecting loss of GABAergic inhibition, occur in an in vivo model of partial complex epilepsy.

A model of self-thinning through local competition.

Explanations of self-thinning in plant populations have focused on plant shape and packing. A dynamic model based on the structure of local interactions successfully reproduces the pattern and can be approximated to identify key parameters and relationships. The approach generates testable new explanations for differences between species and populations, unifies self-thinning with other patterns in plant population dynamics, and indicates why organisms other than plants can follow the law.

Zinc- and iron-rubredoxins from Clostridium pasteurianum at atomic resolution: a high-precision model of a ZnS4 coordination unit in a protein.

The Zn(Scys)4 unit is present in numerous proteins, where it assumes structural, regulatory, or catalytic roles. The same coordination is found naturally around iron in rubredoxins, several structures of which have been refined at resolutions of, or near to, 1 A. The fold of the small protein rubredoxin around its metal ion is an excellent model for many zinc finger proteins. Zn-substituted rubredoxin and its Fe-containing counterpart were both obtained as the products of the expression in Escherichia coli of the rubredoxin-encoding gene from Clostridium pasteurianum. The structures of both proteins have been refined with an anisotropic model at atomic resolution (1.1 A, R = 8.3% for Fe-rubredoxin, and 1.2 A, R = 9.6% for Zn-rubredoxin) and are very similar. The most significant differences are increased lengths of the M-S bonds in Zn-rubredoxin (average length, 2.345 A) as compared with Fe-rubredoxin (average length, 2.262 A). An increase of the CA-CB-SG-M dihedral angles involving Cys-6 and Cys-39, the first cysteines of each of the Cys-Xaa-Xaa-Cys metal binding motifs, has been observed. Another consequence of the replacement of iron by zinc is that the region around residues 36-46 undergoes larger displacements than the remainder of the polypeptide chain. Despite these changes, the main features of the FeS4 site, namely a local 2-fold symmetry and the characteristic network of N-H...S hydrogen bonds, are conserved in the ZnS4 site. The Zn-substituted rubredoxin provides the first precise structure of a Zn(Scys)4 unit in a protein. The nearly identical fold of rubredoxin around iron or zinc suggests that at least in some of the sites where the metal has mainly a structural role-e.g., zinc fingers-the choice of the relevant metal may be directed by its cellular availability and mobilization processes rather than by its chemical nature.

An information theory model of hydrophobic interactions.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

Analysis of genetic instability during mammary tumor progression using a novel selection-based assay for in vivo mutations in a bacteriophage lambda transgene target.

Genetic instability is thought to be responsible for the numerous genotypic changes that occur during neoplastic transformation and metastatic progression. To explore the role of genetic instability at the level of point mutations during mammary tumor development and malignant progression, we combined transgenic mouse models of mutagenesis detection and oncogenesis. Bitransgenic mice were generated that carried both a bacteriophage lambda transgene to assay mutagenesis and a polyomavirus middle T oncogene, mammary gland-targeted expression of which led to metastatic mammary adenocarcinomas. We developed a novel assay for the detection of mutations in the lambda transgene that selects for phage containing forward mutations only in the lambda cII gene, using an hfl- bacterial host. In addition to the relative ease of direct selection, the sensitivity of this assay for both spontaneous and chemically induced mutations was comparable to the widely used mutational target gene, lambda lacI, making the cII assay an attractive alternative for mutant phage recovery for any lambda-based mouse mutagenesis assay system. The frequencies of lambda cII- mutants were not significantly different in normal mammary epithelium, primary mammary adenocarcinomas, and pulmonary metastases. The cII mutational spectra in these tissues consisted mostly of G/C-->A/T transitions, a large fraction of which occurred at CpG dinucleotides. These data suggest that, in this middle T oncogene model of mammary tumor progression, a significant increase in mutagenesis is not required for tumor development or for metastatic progression.

Iron chelates bind nitric oxide and decrease mortality in an experimental model of septic shock.

The hydroxamic acid siderophore ferrioxamine B [FeIII(HDFB)+] and the iron complex of diethylenetri-aminepentaacetic acid [FeIII(DTPA)2-] protected mice against death by septic shock induced by Corynebacterium parvum + lipopolysaccharide. Although FeIII(DTPA)2- was somewhat more effective than FeIII(HDFB)+, the iron-free ligand H4DFB+ was significantly more effective than DTPA. The hydroxamic acid chelator has a much higher iron affinity than the amine carboxylate, allowing for more efficient formation of the FeIII(HDFB)+ complex upon administration of the iron-free ligand. Electrochemical studies show that FeIII(DTPA)2- binds NO stoichiometrically upon reduction to iron(II) at biologically relevant potentials to form a stable NO adduct. In contrast, FeIII(HDFB)+ is a stable and efficient electrocatalyst for the reduction of NO to N2O at biologically relevant potentials. These results suggest that the mechanism of protection against death by septic shock involves NO scavenging and that particularly effective drugs that operate a low dosages may be designed based on the principle of redox catalysis. These complexes constitute a new family of drugs that rely on the special ability of transition metals to activate small molecules. In addition, the wealth of information available on siderophore chemistry and biology provides an intellectual platform for further development.

Trials of the beta model for complex inheritance.

Theoretical advantages of nonparametric logarithm of odds to map polygenic diseases are supported by tests of the beta model that depends on a single logistic parameter and is the only model under which paternal and maternal transmissions to sibs of specified phenotypes are independent. Although it does not precisely describe recurrence risks in monozygous twins, the beta model has greater power to detect family resemblance or linkage than the more general delta model which describes the probability of 0, 1, or 2 alleles identical by descent (ibd) with two parameters. Available data on ibd in sibs are consistent with the beta model, but not with the equally parsimonious but less powerful gamma model that assumes a fixed probability of 1/2 for 1 allele ibd. Additivity of loci on the liability scale is not disproven. A simple equivalence extends the beta model to multipoint analysis.

Ras-catalyzed hydrolysis of GTP: a new perspective from model studies.

Despite the biological and medical importance of signal transduction via Ras proteins and despite considerable kinetic and structural studies of wild-type and mutant Ras proteins, the mechanism of Ras-catalyzed GTP hydrolysis remains controversial. We take a different approach to this problem: the uncatalyzed hydrolysis of GTP is analyzed, and the understanding derived is applied to the Ras-catalyzed reaction. Evaluation of previous mechanistic proposals from this chemical perspective suggests that proton abstraction from the attacking water by a general base and stabilization of charge development on the gamma-phosphoryl oxygen atoms would not be catalytic. Rather, this analysis focuses attention on the GDP leaving group, including the beta-gamma bridge oxygen of GTP, the atom that undergoes the largest change in charge in going from the ground state to the transition state. This leads to a new catalytic proposal in which a hydrogen bond from the backbone amide of Gly-13 to this bridge oxygen is strengthened in the transition state relative to the ground state, within an active site that provides a template complementary to the transition state. Strengthened transition state interactions of the active site lysine, Lys-16, with the beta-nonbridging phosphoryl oxygens and a network of interactions that positions the nucleophilic water molecule and gamma-phosphoryl group with respect to one another may also contribute to catalysis. It is speculated that a significant fraction of the GAP-activated GTPase activity of Ras arises from an additional interaction of the beta-gamma bridge oxygen with an Arg side chain that is provided in trans by GAP. The conclusions for Ras and related G proteins are expected to apply more widely to other enzymes that catalyze phosphoryl (-PO(3)2-) transfer, including kinases and phosphatases.

A joint hazard and time scaling model to compare survival curves.

To provide a more general method for comparing survival experience, we propose a model that independently scales both hazard and time dimensions. To test the curve shape similarity of two time-dependent hazards, h1(t) and h2(t), we apply the proposed hazard relationship, h12(tKt)/ h1(t) = Kh, to h1. This relationship doubly scales h1 by the constant hazard and time scale factors, Kh and Kt, producing a transformed hazard, h12, with the same underlying curve shape as h1. We optimize the match of h12 to h2 by adjusting Kh and Kt. The corresponding survival relationship S12(tKt) = [S1(t)]KtKh transforms S1 into a new curve S12 of the same underlying shape that can be matched to the original S2. We apply this model to the curves for regional and local breast cancer contained in the National Cancer Institute's End Results Registry (1950-1973). Scaling the original regional curves, h1 and S1 with Kt = 1.769 and Kh = 0.263 produces transformed curves h12 and S12 that display congruence with the respective local curves, h2 and S2. This similarity of curve shapes suggests the application of the more complete curve shapes for regional disease as templates to predict the long-term survival pattern for local disease. By extension, this similarity raises the possibility of scaling early data for clinical trial curves according to templates of registry or previous trial curves, projecting long-term outcomes and reducing costs. The proposed model includes as special cases the widely used proportional hazards (Kt = 1) and accelerated life (KtKh = 1) models.

The folding mechanism of larger model proteins: role of native structure.

The folding mechanism of a 125-bead heteropolymer model for proteins is investigated with Monte Carlo simulations on a cubic lattice. Sequences that do and do not fold in a reasonable time are compared. The overall folding behavior is found to be more complex than that of models for smaller proteins. Folding begins with a rapid collapse followed by a slow search through the semi-compact globule for a sequence-dependent stable core with about 30 out of 176 native contacts which serves as the transition state for folding to a near-native structure. Efficient search for the core is dependent on structural features of the native state. Sequences that fold have large amounts of stable, cooperative structure that is accessible through short-range initiation sites, such as those in anti-parallel sheets connected by turns. Before folding is completed, the system can encounter a second bottleneck, involving the condensation and rearrangement of surface residues. Overly stable local structure of the surface residues slows this stage of the folding process. The relation of the results from the 125-mer model studies to the folding of real proteins is discussed.

Perturbations in maturation of secretory proteins and their association with endoplasmic reticulum chaperones in a cell culture model for epithelial ischemia.

The effects of ischemia on the maturation of secretory proteins are not well understood. Among several events that occur during ischemia-reperfusion are a rapid and extensive decrease in ATP levels and an alteration of cellular oxidative state. Since the normal folding and assembly of secretory proteins are mediated by endoplasmic reticulum (ER) molecular chaperones, the function of which depends on ATP and maintenance of an appropriate redox environment, ischemia might be expected to perturb folding of secretory proteins. In this study, whole animal and cultured cell models for the epithelial ischemic state were used to examine this possibility. After acute kidney ischemia, marked increases in the mRNA levels of the ER chaperones glucose-regulated protein (grp)78/immunoglobulin-binding protein (BiP), grp94, and ER protein (ERp)72 were noted. Likewise, when cellular ATP was depleted to less than 10% of control with antimycin A, mRNA levels of BiP, ERp72, and grp94 were increased in kidney and thyroid epithelial cell culture models. Since the signal for the up-regulation of these stress proteins is believed to be the accumulation of misfolded/misassembled secretory proteins in the ER, their induction after ischemia in vivo and antimycin treatment of cultured cells suggests that maturation of secretory proteins in the ER lumen might indeed be perturbed. To analyze the effects of antimycin A on the maturation of secretory proteins, we studied the fate of thyroglobulin (Tg), a large oligomeric secretory glycoprotein, the folding and assembly of which seems to require a variety of ER chaperones. Treatment of cultured thyroid epithelial cells with antimycin A greatly inhibited ( > 90%) the secretion of Tg. Sucrose density gradient analysis revealed that in antimycin A-treated cells Tg associates into large macromolecular complexes which, by immunofluorescence, appeared to localize to the ER. Furthermore, coimmunoprecipitation studies after antimycin A treatment demonstrated that Tg stably associates with BiP, grp94, and ERp72. Together, our results suggest that a key cellular lesion in ischemia is the misfolding of secretory proteins as they transit the ER, and this leads not only to increased expression of ER chaperones but also to their stable association with and the subsequent retention of at least some misfolded secretory proteins.

Evolutionarily conserved Galphabetagamma binding surfaces support a model of the G protein-receptor complex.

The pivotal role of G proteins in sensory, hormonal, inflammatory, and proliferative responses has provoked intense interest in understanding how they interact with their receptors and effectors. Nonetheless, the locations of the receptors and effector binding sites remain poorly characterized, although nearly complete structures of the alphabetagamma heterotrimeric complex are available. Here we apply evolutionary trace (ET) analysis [Lichtarge, O., Bourne, H. R. & Cohen, F. E. (1996) J. Mol. Biol. 257, 342-358] to propose plausible locations for these sites. On each subunit, ET identifies evolutionarily selected surfaces composed of residues that do not vary within functional subgroups and that form spatial clusters. Four clusters correctly identify subunit interfaces, and additional clusters on Galpha point to likely receptor or effector binding sites. Our results implicate the conformationally variable region of Galpha in an effector binding role. Furthermore the range of predicted interactions between the receptor and Galphabetagamma, is sufficiently limited that we can build a low resolution and testable model of the receptor-G protein complex.