Complex dynamics of multilocus systems subjected to cyclical selection.

Earlier we have shown that oscillations with a long period ("supercycles") may arise in two-locus systems experiencing cyclical selection with a short period. However, this mode of complex limiting behavior appeared to be possible for narrow ranges of parameters. Here we demonstrate that a multilocus system subjected to stabilizing selection with cyclically moving optimum can generate ubiquitous complex limiting behavior including supercycles, T-cycles, and chaotic-like phenomena. This mode of multilocus dynamics far exceeds the potential attainable under ordinary selection models resulting in simple behavior. It may represent a novel evolutionary mechanism increasing genetic diversity over long-term time periods.

Mantle dynamics and seismic tomography

Three-dimensional imaging of the Earth's interior, called seismic tomography, has achieved breakthrough advances in the last two decades, revealing fundamental geodynamical processes throughout the Earth's mantle and core. Convective circulation of the entire mantle is taking place, with subducted oceanic lithosphere sinking into the lower mantle, overcoming the resistance to penetration provided by the phase boundary near 650-km depth that separates the upper and lower mantle. The boundary layer at the base of the mantle has been revealed to have complex structure, involving local stratification, extensive structural anisotropy, and massive regions of partial melt. The Earth's high Rayleigh number convective regime now is recognized to be much more interesting and complex than suggested by textbook cartoons, and continued advances in seismic tomography, geodynamical modeling, and high-pressure–high-temperature mineral physics will be needed to fully quantify the complex dynamics of our planet's interior.

Reorganization of an arid ecosystem in response to recent climate change

Natural ecosystems contain many individuals and species interacting with each other and with their abiotic environment. Such systems can be expected to exhibit complex dynamics in which small perturbations can be amplified to cause large changes. Here, we document the reorganization of an arid ecosystem that has occurred since the late 1970s. The density of woody shrubs increased 3-fold. Several previously common animal species went locally extinct, while other previously rare species increased. While these changes are symptomatic of desertification, they were not caused by livestock grazing or drought, the principal causes of historical desertification. The changes apparently were caused by a shift in regional climate: since 1977 winter precipitation throughout the region was substantially higher than average for this century. These changes illustrate the kinds of large, unexpected responses of complex natural ecosystems that can occur in response to both natural perturbations and human activities.

Drought-induced shift of a forest–woodland ecotone: Rapid landscape response to climate variation

In coming decades, global climate changes are expected to produce large shifts in vegetation distributions at unprecedented rates. These shifts are expected to be most rapid and extreme at ecotones, the boundaries between ecosystems, particularly those in semiarid landscapes. However, current models do not adequately provide for such rapid effects—particularly those caused by mortality—largely because of the lack of data from field studies. Here we report the most rapid landscape-scale shift of a woody ecotone ever documented: in northern New Mexico in the 1950s, the ecotone between semiarid ponderosa pine forest and piñon–juniper woodland shifted extensively (2 km or more) and rapidly (<5 years) through mortality of ponderosa pines in response to a severe drought. This shift has persisted for 40 years. Forest patches within the shift zone became much more fragmented, and soil erosion greatly accelerated. The rapidity and the complex dynamics of the persistent shift point to the need to represent more accurately these dynamics, especially the mortality factor, in assessments of the effects of climate change.

Crystal structure of Taq DNA polymerase in complex with an inhibitory Fab: The Fab is directed against an intermediate in the helix-coil dynamics of the enzyme

We report the crystal structure of Thermus aquaticus DNA polymerase I in complex with an inhibitory Fab, TP7, directed against the native enzyme. Some of the residues present in a helical conformation in the native enzyme have adopted a γ turn conformation in the complex. Taken together, structural information that describes alteration of helical structure and solution studies that demonstrate the ability of TP7 to inhibit 100% of the polymerase activity of the enzyme suggest that the change in conformation is probably caused by trapping of an intermediate in the helix-coil dynamics of this helix by the Fab. Antibodies directed against modified helices in proteins have long been anticipated. The present structure provides direct crystallographic evidence. The Fab binds within the DNA binding cleft of the polymerase domain, interacting with several residues that are used by the enzyme in binding the primer:template complex. This result unequivocally corroborates inferences drawn from binding experiments and modeling calculations that the inhibitory activity of this Fab is directly attributable to its interference with DNA binding by the polymerase domain of the enzyme. The combination of interactions made by the Fab residues in both the polymerase and the vestigial editing nuclease domain of the enzyme reveal the structural basis of its preference for binding to DNA polymerases of the Thermus species. The orientation of the structure-specific nuclease domain with respect to the polymerase domain is significantly different from that seen in other structures of this polymerase. This reorientation does not appear to be antibody-induced and implies remarkably high relative mobility between these two domains.

Arp2/3 complex and actin dynamics are required for actin-based mitochondrial motility in yeast

The Arp2/3 complex is implicated in actin polymerization-driven movement of Listeria monocytogenes. Here, we find that Arp2p and Arc15p, two subunits of this complex, show tight, actin-independent association with isolated yeast mitochondria. Arp2p colocalizes with mitochondria. Consistent with this result, we detect Arp2p-dependent formation of actin clouds around mitochondria in intact yeast. Cells bearing mutations in ARP2 or ARC15 genes show decreased velocities of mitochondrial movement, loss of all directed movement and defects in mitochondrial morphology. Finally, we observe a decrease in the velocity and extent of mitochondrial movement in yeast in which actin dynamics are reduced but actin cytoskeletal structure is intact. These results support the idea that the movement of mitochondria in yeast is actin polymerization driven and that this movement requires Arp2/3 complex.

Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.

Cloning of inversion breakpoints in the Anopheles gambiae complex traces a transposable element at the inversion junction

Anopheles arabiensis, one of the two most potent malaria vectors of the gambiae complex, is characterized by the presence of chromosomal paracentric inversions. Elucidation of the nature and the dynamics of these inversions is of paramount importance for the understanding of the population genetics and evolutionary biology of this mosquito and of the impact on malaria epidemiology. We report here the cloning of the breakpoints of the naturally occurring polymorphic inversion 2Rd′ of A. arabiensis. A cDNA clone that cytologically mapped on the proximal breakpoint was the starting material for the isolation of a cosmid clone that spanned the breakpoint. Analysis of the surrounding sequences demonstrated that adjacent to the distal breakpoint lies a repetitive element that exhibits distinct distribution in different A. arabiensis strains. Sequencing analysis of that area revealed elements characteristic of transposable element terminal repeats. We called this presumed transposable element Odysseus. The presence of Odysseus at the junction of the naturally occuring inversion 2Rd′ suggests that the inversion may be the result of the transposable element’s activity. Characteristics of Odysseus’ terminal region as well as its cytological distribution in different strains may indicate a relatively recent activity of Odysseus.

Identification of two distinct human SMC protein complexes involved in mitotic chromosome dynamics

The structural maintenance of chromosomes (SMC) family member proteins previously were shown to play a critical role in mitotic chromosome condensation and segregation in yeast and Xenopus. Other family members were demonstrated to be required for DNA repair in yeast and mammals. Although several different SMC proteins were identified in different organisms, little is known about the SMC proteins in humans. Here, we report the identification of four human SMC proteins that form two distinct heterodimeric complexes in the cell, the human chromosome-associated protein (hCAP)-C and hCAP-E protein complex (hCAP-C/hCAP-E), and the human SMC1 (hSMC1) and hSMC3 protein complex (hSMC1/hSMC3). The hCAP-C/hCAP-E complex is the human ortholog of the Xenopus chromosome-associated protein (XCAP)-C/XCAP-E complex required for mitotic chromosome condensation. We found that a second complex, hSMC1/hSMC3, is required for metaphase progression in mitotic cells. Punctate vs. diffuse distribution patterns of the hCAP-C/hCAP-E and hSMC1/hSMC3 complexes in the interphase nucleus indicate independent behaviors of the two complexes during the cell cycle. These results suggest that two distinct classes of SMC protein complexes are involved in different aspects of mitotic chromosome organization in human cells.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

The carboxyl terminus of the bacteriophage T4 DNA polymerase is required for holoenzyme complex formation

To further elucidate the mechanism and dynamics of bacteriophage T4 holoenzyme formation, a mutant polymerase in which the last six carboxyl-terminal amino acids are deleted, was constructed, overexpressed, and purified to homogeneity. The mutant polymerase, designated ΔC6 exo−, is identical to wild-type exo− polymerase with respect to kcat, kpol, and dissociation constants for nucleotide and DNA substrate. However, unlike wild-type exo− polymerase, the ΔC6 exo− polymerase is unable to interact with the 45 protein to form the stable holoenzyme. A synthetic polypeptide corresponding to the carboxyl terminus of the wild-type exo− polymerase was tested as an in vitro inhibitor of bacteriophage T4 DNA replication. Surprisingly, the peptide does not directly inhibit holoenzyme complex formation by disrupting the interaction of the polymerase with the 45 protein. On the contrary, the peptide appears to disrupt the interaction of the 44/62 protein with the 45 protein, suggesting that the 44/62 protein and the polymerase use the same site on the 45 protein for functional interactions. Data presented are discussed in terms of a model correlating the functionality of the carboxyl terminus of the polymerase for productive interactions with the 45 protein as well as in terms of the 45 protein concomitantly interacting with the 44/62 protein and polymerase.

A theory for controlling cell cycle dynamics using a reversibly binding inhibitor

We demonstrate, by using mathematical modeling of cell division cycle (CDC) dynamics, a potential mechanism for precisely controlling the frequency of cell division and regulating the size of a dividing cell. Control of the cell cycle is achieved by artificially expressing a protein that reversibly binds and inactivates any one of the CDC proteins. In the simplest case, such as the checkpoint-free situation encountered in early amphibian embryos, the frequency of CDC oscillations can be increased or decreased by regulating the rate of synthesis, the binding rate, or the equilibrium constant of the binding protein. In a more complex model of cell division, where size-control checkpoints are included, we show that the same reversible binding reaction can alter the mean cell mass in a continuously dividing cell. Because this control scheme is general and requires only the expression of a single protein, it provides a practical means for tuning the characteristics of the cell cycle in vivo.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

The Yeast Dynactin Complex Is Involved in Partitioning the Mitotic Spindle between Mother and Daughter Cells during Anaphase B

Although vertebrate cytoplasmic dynein can move to the minus ends of microtubules in vitro, its ability to translocate purified vesicles on microtubules depends on the presence of an accessory complex known as dynactin. We have cloned and characterized a novel gene, NIP100, which encodes the yeast homologue of the vertebrate dynactin complex protein p150glued. Like strains lacking the cytoplasmic dynein heavy chain Dyn1p or the centractin homologue Act5p, nip100Δ strains are viable but undergo a significant number of failed mitoses in which the mitotic spindle does not properly partition into the daughter cell. Analysis of spindle dynamics by time-lapse digital microscopy indicates that the precise role of Nip100p during anaphase is to promote the translocation of the partially elongated mitotic spindle through the bud neck. Consistent with the presence of a true dynactin complex in yeast, Nip100p exists in a stable complex with Act5p as well as Jnm1p, another protein required for proper spindle partitioning during anaphase. Moreover, genetic depletion experiments indicate that the binding of Nip100p to Act5p is dependent on the presence of Jnm1p. Finally, we find that a fusion of Nip100p to the green fluorescent protein localizes to the spindle poles throughout the cell cycle. Taken together, these results suggest that the yeast dynactin complex and cytoplasmic dynein together define a physiological pathway that is responsible for spindle translocation late in anaphase.