Variation in per capita interaction strength: Thresholds due to nonlinear dynamics and nonequilibrium conditions

I measured the strength of interaction between a marine herbivore and its growing resource over a realistic range of absolute and relative abundances. The herbivores (hermit crabs: Pagurus spp.) have slow and/or weak functional and numerical responses to epiphytic diatoms (Isthmia nervosa), which show logistic growth in the absence of consumers. By isolating this interaction in containers in the field, I mimicked many of the physical and biological variables characteristic of the intertidal while controlling the densities of focal species. The per capita effects of consumers on the population dynamics of their resource (i.e., interaction strength) were defined by using the relationship between hermit crab density and proportional change in the resource. When this relationship is fit by a Weibull function, a single parameter distinguishes constant interaction strength from one that varies as a function of density. Constant interaction strength causes the proportion of diatoms to fall linearly or proportionally as hermit crab density increases whereas per capita effects that increase with density cause an accelerating decline. Although many mathematical models of species interactions assume linear dynamics and invariant parameters, at least near equilibrium, the per capita effects of hermit crabs on diatoms varied substantially, apparently crossing a threshold from weak to strong when consumption exceeded resource production. This threshold separates a domain of coexistence from one of local extinction of the resource. Such thresholds may help explain trophic cascades, resource compensation, and context-dependent interaction strengths, while indicating a way to predict trophic effects, despite nonlinearities, as a function of vital rates.

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.

Dynamics and elasticity of the fibronectin matrix in living cell culture visualized by fibronectin–green fluorescent protein

Fibronectin (FN) forms the primitive fibrillar matrix in both embryos and healing wounds. To study the matrix in living cell cultures, we have constructed a cell line that secretes FN molecules chimeric with green fluorescent protein. These FN–green fluorescent protein molecules were assembled into a typical matrix that was easily visualized by fluorescence over periods of several hours. FN fibrils remained mostly straight, and they were seen to extend and contract to accommodate movements of the cells, indicating that they are elastic. When fibrils were broken or detached from cells, they contracted to less than one-fourth of their extended length, demonstrating that they are highly stretched in the living culture. Previous work from other laboratories has suggested that cryptic sites for FN assembly may be exposed by tension on FN. Our results show directly that FN matrix fibrils are not only under tension but are also highly stretched. This stretched state of FN is an obvious candidate for exposing the cryptic assembly sites.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

Dynamics and orientation of N+(CD3)3-bromoacetylcholine bound to its binding site on the nicotinic acetylcholine receptor

Dynamic and structural information has been obtained for an analogue of acetylcholine while bound to the agonist binding site on the nicotinic acetylcholine receptor (nAcChoR), using wide-line deuterium solid-state NMR. Analysis of the deuterium lineshape obtained at various temperatures from unoriented nAcChoR membranes labeled with deuterated bromoacetylcholine (BAC) showed that the quaternary ammonium group of the ligand is well constrained within the agonist binding site when compared with the dynamics observed in the crystalline solids. This motional restriction would suggest that a high degree of complementarity exists between the quaternary ammonium group of the ligand and the protein within the agonist binding site. nAcChoR membranes were uniaxially oriented by isopotential centrifugation as determined by phosphorous NMR of the membrane phospholipids. Analysis of the deuterium NMR lineshape of these oriented membranes enriched with the nAcChoR labeled with N+(CD3)3-BAC has enabled us to determine that the angle formed between the quaternary ammonium group of the BAC and the membrane normal is 42° in the desensitized form of the receptor. This measurement allows us to orient in part the bound ligand within the proposed receptor binding site.

The role of structure, energy landscape, dynamics, and allostery in the enzymatic function of myoglobin

The grail of protein science is the connection between structure and function. For myoglobin (Mb) this goal is close. Described as only a passive dioxygen storage protein in texts, we argue here that Mb is actually an allosteric enzyme that can catalyze reactions among small molecules. Studies of the structural, spectroscopic, and kinetic properties of Mb lead to a model that relates structure, energy landscape, dynamics, and function. Mb functions as a miniature chemical reactor, concentrating and orienting diatomic molecules such as NO, CO, O2, and H2O2 in highly conserved internal cavities. Reactions can be controlled because Mb exists in distinct taxonomic substates with different catalytic properties and connectivities of internal cavities.

Constraints on nebular dynamics and chemistry based on observations of annealed magnesium silicate grains in comets and in disks surrounding Herbig Ae/Be stars

Understanding dynamic conditions in the Solar Nebula is the key to prediction of the material to be found in comets. We suggest that a dynamic, large-scale circulation pattern brings processed dust and gas from the inner nebula back out into the region of cometesimal formation—extending possibly hundreds of astronomical units (AU) from the sun—and that the composition of comets is determined by a chemical reaction network closely coupled to the dynamic transport of dust and gas in the system. This scenario is supported by laboratory studies of Mg silicates and the astronomical data for comets and for protoplanetary disks associated with young stars, which demonstrate that annealing of nebular silicates must occur in conjunction with a large-scale circulation. Mass recycling of dust should have a significant effect on the chemical kinetics of the outer nebula by introducing reduced, gas-phase species produced in the higher temperature and pressure environment of the inner nebula, along with freshly processed grains with “clean” catalytic surfaces to the region of cometesimal formation. Because comets probably form throughout the lifetime of the Solar Nebula and processed (crystalline) grains are not immediately available for incorporation into the first generation of comets, an increasing fraction of dust incorporated into a growing comet should be crystalline olivine and this fraction can serve as a crude chronometer of the relative ages of comets. The formation and evolution of key organic and biogenic molecules in comets are potentially of great consequence to astrobiology.

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.

Spindle Dynamics and the Role of γ-Tubulin in Early Caenorhabditis elegans EmbryosV⃞

γ-Tubulin is a ubiquitous and highly conserved component of centrosomes in eukaryotic cells. Genetic and biochemical studies have demonstrated that γ-tubulin functions as part of a complex to nucleate microtubule polymerization from centrosomes. We show that, as in other organisms, Caenorhabditis elegans γ-tubulin is concentrated in centrosomes. To study centrosome dynamics in embryos, we generated transgenic worms that express GFP::γ-tubulin or GFP::β-tubulin in the maternal germ line and early embryos. Multiphoton microscopy of embryos produced by these worms revealed the time course of daughter centrosome appearance and growth and the differential behavior of centrosomes destined for germ line and somatic blastomeres. To study the role of γ-tubulin in nucleation and organization of spindle microtubules, we used RNA interference (RNAi) to deplete C. elegans embryos of γ-tubulin. γ-Tubulin (RNAi) embryos failed in chromosome segregation, but surprisingly, they contained extensive microtubule arrays. Moderately affected embryos contained bipolar spindles with dense and long astral microtubule arrays but with poorly organized kinetochore and interpolar microtubules. Severely affected embryos contained collapsed spindles with numerous long astral microtubules. Our results suggest that γ-tubulin is not absolutely required for microtubule nucleation in C. elegans but is required for the normal organization and function of kinetochore and interpolar microtubules.

Toxoplasma invasion: the parasitophorous vacuole is formed from host cell plasma membrane and pinches off via a fission pore.

Most intracellular pathogens avoid lysing their host cells during invasion by wrapping themselves in a vacuolar membrane. This parasitophorous vacuole membrane (PVM) is often retained, serving as a critical transport interface between the parasite and the host cell cytoplasm. To test whether the PVM formed by the parasite Toxoplasma gondii is derived from host cell membrane or from lipids secreted by the parasite, we used time-resolved capacitance measurements and video microscopy to assay host cell surface area during invasion. We observed no significant change in host cell surface area during PVM formation, demonstrating that the PVM consists primarily of invaginated host cell membrane. Pinching off of the PVM from the host cell membrane occurred after an unexpected delay (34-305 sec) and was seen as a 0.219 +/- 0.006 pF drop in capacitance, which corresponds well to the predicted surface area of the entire PVM (30-33 microns2). The formation and closure of a fission pore connecting the extracellular medium and the vacuolar space was detected as the PVM pinched off. This final stage of parasite entry was accomplished without any breach in cell membrane integrity.

Host-parasite dynamics and outgrowth of virus containing a single K70R amino acid change in reverse transcriptase are responsible for the loss of human immunodeficiency virus type 1 RNA load suppression by zidovudine.

The association between human immunodeficiency virus type I (HIV-1) RNA load changes and the emergence of resistant virus variants was investigated in 24 HIV-1-infected asymptomatic persons during 2 years of treatment with zidovudine by sequentially measuring serum HIV-1 RNA load and the relative amounts of HIV-1 RNA containing mutations at reverse transcriptase (RT) codons 70 (K-->R), 41 (M-->L), and 215 (T-->Y/F). A mean maximum decline in RNA load occurred during the first month, followed by a resurgence between 1 and 3 months, which appeared independent of drug-resistance. Mathematical modeling suggests that this resurgence is caused by host-parasite dynamics, and thus reflects infection of the transiently increased numbers of CD4+ lymphocytes. Between 3 and 6 months of treatment, the RNA load returned to baseline values, which was associated with the emergence of virus containing a single lysine to arginine amino acid change at RT codon 70, only conferring an 8-fold reduction in susceptibility. Despite the relative loss of RNA load suppression, selection toward mutations at RT codons 215 and 41 continued. Identical patterns were observed in the mathematical model. While host-parasite dynamics and outgrowth of low-level resistant virus thus appear responsible for the loss of HIV-1 RNA load suppression, zidovudine continues to select for alternative mutations, conferring increasing levels of resistance.