Dynamical principles in biological processes: A model of charge migration in proteins and DNA

The generalized master equations (GMEs) that contain multiple time scales have been derived quantum mechanically. The GME method has then been applied to a model of charge migration in proteins that invokes the hole hopping between local amino acid sites driven by the torsional motions of the floppy backbones. This model is then applied to analyze the experimental results for sequence-dependent long-range hole transport in DNA reported by Meggers et al. [Meggers, E., Michel-Beyerle, M. E., & Giese, B. (1998) J. Am. Chem. Soc. 120, 12950–12955]. The model has also been applied to analyze the experimental results of femtosecond dynamics of DNA-mediated electron transfer reported by Zewail and co-workers [Wan, C., Fiebig, T., Kelley, S. O., Treadway, C. R., Barton, J. K. & Zewail, A. H. (1999) Proc. Natl. Acad. Sci. USA 96, 6014–6019]. The initial events in the dynamics of protein folding have begun to attract attention. The GME obtained in this paper will be applicable to this problem.

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.

A knock-out model of paroxysmal nocturnal hemoglobinuria: Pig-a(-) hematopoiesis is reconstituted following intercellular transfer of GPI-anchored proteins.

We created a "knockout" embryonic stem cell via targeted disruption of the phosphatidylinositol glycan class A (Pig-a) gene, resulting in loss of expression of cell surface glycosyl phosphatidylinositol-anchored proteins and reproducing the mutant phenotype of the human disease paroxysmal nocturnal hemoglobinuria. Morphogenesis of Pig-a- embryoid bodies (EB) in vitro was grossly aberrant and, unlike EB derived from normal embryonic stem cells, Pig-A EB produced no secondary hematopoietic colonies. Chimeric EB composed of control plus Pig-A- cells, however, appeared normal, and hematopoiesis from knock-out cells was reconstituted. Transfer in situ of glycosyl phosphatidylinositol-anchored proteins from normal to knock-out cells was demonstrated by two-color fluorescent analysis, suggesting a possible mechanism for these functional effects. Hematopoietic cells with mutated PIG-A genes in humans with paroxysmal nocturnal hemoglobinuria may be subject to comparable pathophysiologic processes and amenable to similar therapeutic protein transfer.

A general model of invariant chain association with class II major histocompatibility complex proteins.

The binding of invariant chain to major histocompatibility complex (MHC) proteins is an important step in processing of MHC class II proteins and in antigen presentation. The question of how invariant chain can bind to all MHC class II proteins is central to understanding these processes. We have employed molecular modeling to predict the structure of class II-associated invariant chain peptide (CLIP)-MHC protein complexes and to ask whether the predicted mode of association could be general across all MHC class II proteins. CLIP fits identically into the MHC class II alleles HLA-DR3, I-Ak, I-Au, and I-Ad, with a consistent pattern of hydrogen bonds, contacts, and hydrophobic burial and without bad contacts. Our model predicts the burial of CLIP residues Met-91 and Met-99 in the deep P1 and P9 anchor pockets and other detailed interactions, which we have compared with available data. The predicted pattern of I-A allele-specific effects on CLIP binding is very similar to that observed experimentally by alanine-scanning mutations of CLIP. Together, these results indicate that CLIP may bind in a single, general way across products of MHC class II alleles.

A model for the structure of the P domains in the subtilisin-like prohormone convertases

The proprotein convertases are a family of at least seven calcium-dependent endoproteases that process a wide variety of precursor proteins in the secretory pathway. All members of this family possess an N-terminal proregion, a subtilisin-like catalytic module, and an additional downstream well-conserved region of ≈150 amino acid residues, the P domain, which is not found in any other subtilase. The pro and catalytic domains cannot be expressed in the absence of the P domains; their thermodynamic instability may be attributable to the presence of large numbers of negatively charged Glu and Asp side chains in the substrate binding region for recognition of multibasic residue cleavage sites. Based on secondary structure predictions, we here propose that the P domains consist of 8-stranded β-barrels with well-organized inner hydrophobic cores, and therefore are independently folded components of the proprotein convertases. We hypothesize further that the P domains are integrated through strong hydrophobic interactions with the catalytic domains, conferring structural stability and regulating the properties and activity of the convertases. A molecular model of these interdomain interactions is proposed in this report.

Interaction of SP100 with HP1 proteins: A link between the promyelocytic leukemia-associated nuclear bodies and the chromatin compartment

The PML/SP100 nuclear bodies (NBs) were first described as discrete subnuclear structures containing the SP100 protein. Subsequently, they were shown to contain the PML protein which is part of the oncogenic PML-RARα hybrid produced by the t(15;17) chromosomal translocation characteristic of acute promyelocytic leukemia. Yet, the physiological role of these nuclear bodies remains unknown. Here, we show that SP100 binds to members of the heterochromatin protein 1 (HP1) families of non-histone chromosomal proteins. Further, we demonstrate that a naturally occurring splice variant of SP100, here called SP100-HMG, is a member of the high mobility group-1 (HMG-1) protein family and may thus possess DNA-binding potential. Both HP1 and SP100-HMG concentrate in the PML/SP100 NBs, and overexpression of SP100 leads to enhanced accumulation of endogenous HP1 in these structures. When bound to a promoter, SP100, SP100-HMG and HP1 behave as transcriptional repressors in transfected mammalian cells. These observations present molecular evidence for an association between the PML/SP100 NBs and the chromatin nuclear compartment. They support a model in which the NBs may play a role in certain aspects of chromatin dynamics.

The Escherichia coli SOS mutagenesis proteins UmuD and UmuD′ interact physically with the replicative DNA polymerase

The Escherichia coli umuDC operon is induced in response to replication-blocking DNA lesions as part of the SOS response. UmuD protein then undergoes an RecA-facilitated self-cleavage reaction that removes its N-terminal 24 residues to yield UmuD′. UmuD′, UmuC, RecA, and some form of the E. coli replicative DNA polymerase, DNA polymerase III holoenzyme, function in translesion synthesis, the potentially mutagenic process of replication over otherwise blocking lesions. Furthermore, it has been proposed that, before cleavage, UmuD together with UmuC acts as a DNA damage checkpoint system that regulates the rate of DNA synthesis in response to DNA damage, thereby allowing time for accurate repair to take place. Here we provide direct evidence that both uncleaved UmuD and UmuD′ interact physically with the catalytic, proofreading, and processivity subunits of the E. coli replicative polymerase. Consistent with our model proposing that uncleaved UmuD and UmuD′ promote different events, UmuD and UmuD′ interact differently with DNA polymerase III: whereas uncleaved UmuD interacts more strongly with β than it does with α, UmuD′ interacts more strongly with α than it does with β. We propose that the protein–protein interactions we have characterized are part of a higher-order regulatory system of replication fork management that controls when the umuDC gene products can gain access to the replication fork.

Rapid treadmilling of brain microtubules free of microtubule-associated proteins in vitro and its suppression by tau

We have determined the treadmilling rate of brain microtubules (MTs) free of MT-associated proteins (MAPs) at polymer mass steady state in vitro by using [3H]GTP-exchange. We developed buffer conditions that suppressed dynamic instability behavior by ≈10-fold to minimize the contribution of dynamic instability to total tubulin-GTP exchange. The MTs treadmilled rapidly under the suppressed dynamic instability conditions, at a minimum rate of 0.2 μm/min. Thus, rapid treadmilling is an intrinsic property of MAP-free MTs. Further, we show that tau, an axonal stabilizing MAP involved in Alzheimer’s disease, strongly suppresses the treadmilling rate. These results indicate that tau’s function in axons might involve suppression of axonal MT treadmilling. We describe mathematically how treadmilling and dynamic instability are mechanistically distinct MT behaviors. Finally, we present a model that explains how small changes in the critical tubulin subunit concentration at MT minus ends, caused by intrinsic differences in rate constants or regulatory proteins, could produce large changes in the treadmilling rate.

A family of fibrinogen-related proteins that precipitates parasite-derived molecules is produced by an invertebrate after infection

After infection with the digenetic trematode Echinostoma paraensei, hemolymph of the snail Biomphalaria glabrata contains lectins comprised of 65-kDa subunits that precipitate polypeptides secreted by E. paraensei intramolluscan larvae. Comparable activity is lacking in hemolymph of uninfected snails. Three different cDNAs with sequence similarities to peptides derived from the 65-kDa lectins were obtained and unexpectedly found to encode fibrinogen-related proteins (FREPs). These FREPs also contained regions with sequence similarity to Ig superfamily members. B. glabrata has at least five FREP genes, three of which are expressed at increased levels after infection. Elucidation of components of the defense system of B. glabrata is relevant because this snail is an intermediate host for Schistosoma mansoni, the most widely distributed causative agent of human schistosomiasis. These results are novel in suggesting a role for invertebrate FREPs in recognition of parasite-derived molecules and also provide a model for investigating the diversity of molecules functioning in nonself-recognition in an invertebrate.

A conserved residue in the tip proteins of CS1 and CFA/I pili of enterotoxigenic Escherichia coli that is essential for adherence

Enterotoxigenic Escherichia coli associated with human diarrheal disease utilize any of a limited group of serologically distinguishable pili for attachment to intestinal cells. These include CS1 and CFA/I pili. We show here that chemical modification of arginyl residues in CS1 pili abolishes CS1-mediated agglutination of bovine erythrocytes, which serves as a model system for attachment. Alanine substitution of the single arginyl residue in CooA, the major pilin, had no effect on the assembly of pili or on hemagglutination. In contrast, substitution of alanine for R181 in CooD, the minor pilin associated with the pilus tip, abolished hemagglutination, and substitution of R20 reduced hemagglutination. Neither of these substitutions affected CS1 pilus assembly. This shows that CooD is essential for CS1-mediated attachment and identifies specific residues that are involved in receptor binding but not in pilus assembly. In addition to mediating agglutination of bovine erythrocytes, CFA/I also mediates agglutination of human erythrocytes. Substitution of R181 by alanine in the CooD homolog, CfaE, abolished both of these reactions. We conclude that the same region of the pilus tip protein is involved in adherence of CS1 and CFA/I pili, although their receptor specificities differ. This suggests that the region of the pilus tip adhesin protein that includes R181 might be an appropriate target for therapeutic intervention or for a vaccine to protect against human diarrhea caused by enterotoxigenic E. coli strains that have serologically different pili.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Evolution of the folding ability of proteins through functional selection

An evolutionary process is simulated with a simple spin-glass-like model of proteins to examine the origin of folding ability. At each generation, sequences are randomly mutated and subjected to a simulation of the folding process based on the model. According to the frequency of local configurations at the active sites, sequences are selected and passed to the next generation. After a few hundred generations, a sequence capable of folding globally into a native conformation emerges. Moreover, the selected sequence has a distinct energy minimum and an anisotropic funnel on the energy surface, which are the imperative features for fast folding of proteins. The proposed model reveals that the functional selection on the local configurations leads a sequence to fold globally into a conformation at a faster rate.

Identification of the prooxidant site of human ceruloplasmin: A model for oxidative damage by copper bound to protein surfaces

Free transition metal ions oxidize lipids and lipoproteins in vitro; however, recent evidence suggests that free metal ion-independent mechanisms are more likely in vivo. We have shown previously that human ceruloplasmin (Cp), a serum protein containing seven Cu atoms, induces low density lipoprotein oxidation in vitro and that the activity depends on the presence of a single, chelatable Cu atom. We here use biochemical and molecular approaches to determine the site responsible for Cp prooxidant activity. Experiments with the His-specific reagent diethylpyrocarbonate (DEPC) showed that one or more His residues was specifically required. Quantitative [14C]DEPC binding studies indicated the importance of a single His residue because only one was exposed upon removal of the prooxidant Cu. Plasmin digestion of [14C]DEPC-treated Cp (and N-terminal sequence analysis of the fragments) showed that the critical His was in a 17-kDa region containing four His residues in the second major sequence homology domain of Cp. A full length human Cp cDNA was modified by site-directed mutagenesis to give His-to-Ala substitutions at each of the four positions and was transfected into COS-7 cells, and low density lipoprotein oxidation was measured. The prooxidant site was localized to a region containing His426 because CpH426A almost completely lacked prooxidant activity whereas the other mutants expressed normal activity. These observations support the hypothesis that Cu bound at specific sites on protein surfaces can cause oxidative damage to macromolecules in their environment. Cp may serve as a model protein for understanding mechanisms of oxidant damage by copper-containing (or -binding) proteins such as Cu, Zn superoxide dismutase, and amyloid precursor protein.

Folding and aggregation of designed proteins

Protein aggregation is studied by following the simultaneous folding of two designed identical 20-letter amino acid chains within the framework of a lattice model and using Monte Carlo simulations. It is found that protein aggregation is determined by elementary structures (partially folded intermediates) controlled by local contacts among some of the most strongly interacting amino acids and formed at an early stage in the folding process.