Excited-state dynamics in photosystem II: Insights from the x-ray crystal structure

The heart of oxygenic photosynthesis is photosystem II (PSII), a multisubunit protein complex that uses solar energy to drive the splitting of water and production of molecular oxygen. The effectiveness of the photochemical reaction center of PSII depends on the efficient transfer of excitation energy from the surrounding antenna chlorophylls. A kinetic model for PSII, based on the x-ray crystal structure coordinates of 37 antenna and reaction center pigment molecules, allows us to map the major energy transfer routes from the antenna chlorophylls to the reaction center chromophores. The model shows that energy transfer to the reaction center is slow compared with the rate of primary electron transport and depends on a few bridging chlorophyll molecules. This unexpected energetic isolation of the reaction center in PSII is similar to that found in the bacterial photosystem, conflicts with the established view of the photophysics of PSII, and may be a functional requirement for primary photochemistry in photosynthesis. In addition, the model predicts a value for the intrinsic photochemical rate constant that is 4 times that found in bacterial reaction centers.

Intrinsic compressibility and volume compression in solvated proteins by molecular dynamics simulation at high pressure.

Constant pressure and temperature molecular dynamics techniques have been employed to investigate the changes in structure and volumes of two globular proteins, superoxide dismutase and lysozyme, under pressure. Compression (the relative changes in the proteins' volumes), computed with the Voronoi technique, is closely related with the so-called protein intrinsic compressibility, estimated by sound velocity measurements. In particular, compression computed with Voronoi volumes predicts, in agreement with experimental estimates, a negative bound water contribution to the apparent protein compression. While the use of van der Waals and molecular volumes underestimates the intrinsic compressibilities of proteins, Voronoi volumes produce results closer to experimental estimates. Remarkably, for two globular proteins of very different secondary structures, we compute identical (within statistical error) protein intrinsic compressions, as predicted by recent experimental studies. Changes in the protein interatomic distances under compression are also investigated. It is found that, on average, short distances compress less than longer ones. This nonuniform contraction underlines the peculiar nature of the structural changes due to pressure in contrast with temperature effects, which instead produce spatially uniform changes in proteins. The structural effects observed in the simulations at high pressure can explain protein compressibility measurements carried out by fluorimetric and hole burning techniques. Finally, the calculation of the proteins static structure factor shows significant shifts in the peaks at short wavenumber as pressure changes. These effects might provide an alternative way to obtain information concerning compressibilities of selected protein regions.

Protonation dynamics of the extracellular and cytoplasmic surface of bacteriorhodopsin in the purple membrane.

The dynamics of proton binding to the extracellular and the cytoplasmic surfaces of the purple membrane were measured by laser-induced proton pulses. Purple membranes, selectively labeled by fluorescein at Lys-129 of bacteriorhodopsin, were pulsed by protons released in the aqueous bulk from excited pyranine (8-hydroxy-1,3,6-pyrenetrisulfonate) and the reaction of protons with the indicators was measured. Kinetic analysis of the data imply that the two faces of the membrane differ in their buffer capacities and in their rates of interaction with bulk protons. The extracellular surface of the purple membrane contains one anionic proton binding site per protein molecule with pK = 5.1. This site is within a Coulomb cage radius (approximately 15 A) from Lys-129. The cytoplasmic surface of the purple membrane bears 4-5 protonable moieties (pK = 5.1) that, due to close proximity, function as a common proton binding site. The reaction of the proton with this cluster is at a very fast rate (3.10(10) M-1.s-1). The proximity between the elements is sufficiently high that even in 100 mM NaCl they still function as a cluster. Extraction of the chromophore retinal from the protein has a marked effect on the carboxylates of the cytoplasmic surface, and two to three of them assume positions that almost bar their reaction with bulk protons. The protonation dynamics determined at the surface of the purple membrane is of relevance both for the vectorial proton transport mechanism of bacteriorhodopsin and for energy coupling, not only in halobacteria, but also in complex chemiosmotic systems such as mitochondrial and thylakoid membranes.

Ultrasensitivity in the mitogen-activated protein kinase cascade.

The mitogen-activated protein kinase (MAPK) cascade is a highly conserved series of three protein kinases implicated in diverse biological processes. Here we demonstrate that the cascade arrangement has unexpected consequences for the dynamics of MAPK signaling. We solved the rate equations for the cascade numerically and found that MAPK is predicted to behave like a highly cooperative enzyme, even though it was not assumed that any of the enzymes in the cascade were regulated cooperatively. Measurements of MAPK activation in Xenopus oocyte extracts confirmed this prediction. The stimulus/response curve of the MAPK was found to be as steep as that of a cooperative enzyme with a Hill coefficient of 4-5, well in excess of that of the classical allosteric protein hemoglobin. The shape of the MAPK stimulus/ response curve may make the cascade particularly appropriate for mediating processes like mitogenesis, cell fate induction, and oocyte maturation, where a cell switches from one discrete state to another.

Ultraviolet-induced movement of the human DNA repair protein, Xeroderma pigmentosum type G, in the nucleus.

Xeroderma pigmentosum type G (XPG) is a human genetic disease exhibiting extreme sensitivity to sunlight. XPG patients are defective XPG endonuclease, which is an enzyme essential for DNA repair of the major kinds of solar ultraviolet (UV)-induced DNA damages. Here we describe a novel dynamics of this protein within the cell nucleus after UV irradiation of human cells. Using confocal microscopy, we have localized the immunofluorescent, antigenic signal of XPG protein to foci throughout the cell nucleus. Our biochemical studies also established that XPG protein forms a tight association with nuclear structure(s). In human skin fibroblast cells, the number of XPG foci decreased within 2 h after UV irradiation, whereas total nuclear XPG fluorescence intensity remained constant, suggesting redistribution of XPG from a limited number of nuclear foci to the nucleus overall. Within 8 h after UV, most XPG antigenic signal was found as foci. Using beta-galactosidase-XPG fusion constructs (beta-gal-XPG) transfected into HeLa cells, we have identified a single region of XPG that is evidently responsible both for foci formation and for the UV dynamic response. The fusion protein carrying the C terminus of XPG (amino acids 1146-1185) localized beta-gal specific antigenic signal to foci and to the nucleolus regions. After UV irradiation, antigenic beta-gal translocated reversibly from the subnuclear structures to the whole nucleus with kinetics very similar to the movements of XPG protein. These findings lead us to propose a model in which distribution of XPG protein may regulate the rate of DNA repair within transcriptionally active and inactive compartments of the cell nucleus.

Interaction of the histone (H3-H4)2 tetramer of the nucleosome with positively supercoiled DNA minicircles: Potential flipping of the protein from a left- to a right-handed superhelical form.

We have studied the ability of the histone (H3-H4)2 tetramer, the central part of the nucleosome of eukaryotic chromatin, to form particles on DNA minicircles of negative and positive superhelicities, and the effect of relaxing these particles with topoisomerase I. The results show that even modest positive torsional stress from the DNA, and in particular that generated by DNA thermal fluctuations, can trigger a major, reversible change in the conformation of the particle. Neither a large excess of naked DNA, nor a crosslink between the two H3s prevented the transition from one form to the other. This suggested that during the transition, the histones neither dissociated from the DNA nor were even significantly reshuffled. Moreover, the particles reconstituted on negatively and positively supercoiled minicircles look similar under electron microscopy. These data agree best with a transition involving a switch of the wrapped DNA from a left- to a right-handed superhelix. It is further proposed, based on the left-handed overall superhelical conformation of the tetramer within the octamer [Arents, G., Burlingame, R. W., Wang, B. C., Love, W. E. & Moudrianakis, E. N. (1991) Proc. Natl.Acad. Sci. USA 88, 10148-10152] that this change in DNA topology is mediated by a similar change in the topology of the tetramer itself, which may occur through a rotation (or a localized deformation) of the two H3-H4 dimers about their H3-H3 interface. Potential implications of this model for nucleosome dynamics in vivo are discussed.

Direct observation of fast protein folding: the initial collapse of apomyoglobin.

The rapid refolding dynamics of apomyoglobin are followed by a new temperature-jump fluorescence technique on a 15-ns to 0.5-ms time scale in vitro. The apparatus measures the protein-folding history in a single sweep in standard aqueous buffers. The earliest steps during folding to a compact state are observed and are complete in under 20 micros. Experiments on mutants and consideration of steady-state CD and fluorescence spectra indicate that the observed microsecond phase monitors assembly of an A x (H x G) helix subunit. Measurements at different viscosities indicate diffusive behavior even at low viscosities, in agreement with motions of a solvent-exposed protein during the initial collapse.

N2 fixation in marine heterotrophic bacteria: dynamics of environmental and molecular regulation.

Molecular and immunological techniques were used to examine N2 fixation in a ubiquitous heterotrophic marine bacterium, the facultative anaerobic Vibrio natriegens. When batch cultures were shifted from aerobic N-replete to anaerobic N-deplete conditions, transcriptional and post-translational regulation of N2 fixation was observed. Levels of nifHDK mRNA encoding the nitrogenase enzyme were highest at 140 min postshift and undetectable between 6 and 9 h later. Immunologically determined levels of nitrogenase enzyme (Fe protein) were highest between 6 and 15 h postshift, and nitrogenase activity peaked between 6 and 9 h postshift, declining by a factor of 2 after 12-15 h. Unlike their regulation in cyanobacteria, Fe protein and nitrogenase activity were present when nifHDK mRNA was absent in V. natriegens, indicating that nitrogenase is stored and stable under anaerobic conditions. Both nifHDK mRNA and Fe protein disappeared within 40 min after cultures were shifted from N2-fixing conditions (anaerobic, N-deplete) to non- N2-fixing conditions (aerobic, N-enriched) but reappeared when shifted to conditions favoring N2 fixation. Thus, unlike other N2-fixing heterotrophic bacteria, nitrogenase must be resynthesized after aerobic exposure in V. natriegens. Immunological detection based on immunoblot (Western) analysis and immunogold labeling correlated positively with nitrogenase activity; no localization of nitrogenase was observed. Because V. natriegens continues to fix N2 for many hours after anaerobic induction, this species may play an important role in providing "new" nitrogen in marine ecosystems.

Cloning, expression, and properties of the microtubule-stabilizing protein STOP.

Nerve cells contain abundant subpopulations of cold-stable microtubules. We have previously isolated a calmodulin-regulated brain protein, STOP (stable tubule-only polypeptide), which reconstitutes microtubule cold stability when added to cold-labile microtubules in vitro. We have now cloned cDNA encoding STOP. We find that STOP is a 100.5-kDa protein with no homology to known proteins. The primary structure of STOP includes two distinct domains of repeated motifs. The central region of STOP contains 5 tandem repeats of 46 amino acids, 4 with 98% homology to the consensus sequence. The STOP C terminus contains 28 imperfect repeats of an 11-amino acid motif. STOP also contains a putative SH3-binding motif close to its N terminus. In vitro translated STOP binds to both microtubules and Ca2+-calmodulin. When STOP cDNA is expressed in cells that lack cold-stable microtubules, STOP associates with microtubules at 37 degrees C, and stabilizes microtubule networks, inducing cold stability, nocodazole resistance, and tubulin detyrosination on microtubules in transfected cells. We conclude that STOP must play an important role in the generation of microtubule cold stability and in the control of microtubule dynamics in brain.

Myosin dynamics in live Dictyostelium cells.

Conventional myosin plays a key role in the cytoskeletal reorganization necessary for cytokinesis, migration, and morphological changes associated with development in nonmuscle cells. We have made a fusion between the green fluorescent protein (GFP) and the Dictyostelium discoideum myosin heavy chain (GFP-myosin). The unique Dictyostelium system allows us to test the GFP-tagged myosin for activity both in vivo and in vitro. Expression of GFP-myosin rescues all myosin null cell defects. Additionally, GFP-myosin purified from these cells exhibits the same ATPase activities and in vitro motility as wild-type myosin. GFP-myosin is concentrated in the cleavage furrow during cytokinesis and in the posterior cortex of migrating cells. Surprisingly, GFP-myosin concentration increases transiently in the tips of retracting pseudopods. Contrary to previous thinking, this suggests that conventional myosin may play an important role in the dynamics of pseudopods as well as filopodia, lamellipodia, and other cellular protrusions.

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

Sequential cytokine dynamics in chronic rejection of rat renal allografts: roles for cytokines RANTES and MCP-1.

Chronic rejection, the most important cause of long-term graft failure, is thought to result from both alloantigen-dependent and -independent factors. To examine these influences, cytokine dynamics were assessed by semiquantitative competitive reverse transcriptase-PCR and by immunohistology in an established rat model of chronic rejection lf renal allografts. Isograft controls develop morphologic and immunohistologic changes that are similar to renal allograft changes, although quantitatively less intense and at a delayed speed; these are thought to occur secondary to antigen-independent events. Sequential cytokine expression was determined throughout the process. During an early reversible allograft rejection episode, both T-cell associated [interleukin (IL) 2, IL-2 receptor, IL-4, and interferon gamma] and macrophage (IL-1 alpha, tumor necrosis factor alpha, and IL-6) products were up-regulated despite transient immunosuppression. RANTES (regulated upon activation, normal T-cell expressed and secreted) peaked at 2 weeks; intercellular adhesion molecule (ICAM-1) was maximally expressed at 6 weeks. Macrophage products such as monocyte chemoattractant protein (MCP-1) increased dramatically (to 10 times), presaging intense peak macrophage infiltration at 16 weeks. In contrast, in isografts, ICAM-1 peaked at 24 weeks. MCP-1 was maximally expressed at 52 weeks, commensurate with a progressive increase in infiltrating macrophages. Cytokine expression in the spleen of allograft and isograft recipients was insignificant. We conclude that chronic rejection of kidney allografts in rats is predominantly a local macrophage-dependent event with intense up-regulation of macrophage products such as MCP-1, IL-6, and inducible nitric oxide synthase. The cytokine expression in isografts emphasizes the contribution of antigen-independent events. The dynamics of RANTES expression between early and late phases of chronic rejection suggest a key role in mediating the events of the chronic process.

Shigella flexneri surface protein IcsA is sufficient to direct actin-based motility.

Shigella flexneri is a Gram-negative bacterial pathogen that can grow directly in the cytoplasm of infected host cells and uses a form of actin-based motility for intra- and intercellular spread. Moving intracellular bacteria are associated with a polarized "comet tail" composed of actin filaments. IcsA, a 120-kDa outer membrane protein necessary for actin-based motility, is located at a single pole on the surface of the organism, at the junction with the actin tail. Here, we demonstrate that stable expression of IcsA on the surface of Escherichia coli is sufficient to allow actin-dependent movement of E. coli in cytoplasmic extracts, at rates comparable to the movement of S. flexneri in infected cells. Thus, IcsA is the sole Shigella-specific factor required for actin-based motility. Continuous protein synthesis and polarized distribution of the protein are not necessary for actin tail formation or movement. Listeria monocytogenes is an unrelated bacterial pathogen that exhibits similar actin-based intracytoplasmic motility. Actin filament dynamics in the comet tails associated with the two different organisms are essentially identical, which indicates that they have independently evolved mechanisms to interact with the same components of the host cytoskeleton.

A sampling problem in molecular dynamics simulations of macromolecules.

Correlations in low-frequency atomic displacements predicted by molecular dynamics simulations on the order of 1 ns are undersampled for the time scales currently accessible by the technique. This is shown with three different representations of the fluctuations in a macromolecule: the reciprocal space of crystallography using diffuse x-ray scattering data, real three-dimensional Cartesian space using covariance matrices of the atomic displacements, and the 3N-dimensional configuration space of the protein using dimensionally reduced projections to visualize the extent to which phase space is sampled.

Toward an outline of the topography of a realistic protein-folding funnel.

Experimental information on the structure and dynamics of molten globules gives estimates for the energy landscape's characteristics for folding highly helical proteins, when supplemented by a theory of the helix-coil transition in collapsed heteropolymers. A law of corresponding states relating simulations on small lattice models to real proteins possessing many more degrees of freedom results. This correspondence reveals parallels between "minimalist" lattice results and recent experimental results for the degree of native character of the folding transition state and molten globule and also pinpoints the needs of further experiments.