The interspecific mass–density relationship and plant geometry

We present an a priori theoretical framework for the interspecific allometric relationship between stand mass and plant population density. Our model predicts a slope of −\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\frac{1}{3}\end{equation*}\end{document} between the logarithm of stand mass and the logarithm of stand density, thus conflicting with a previously assumed slope of −½. Our model rests on a heuristic separation of resource-limited living mass and structural mass in the plant body. We point out that because of similar resource requirements among plants of different sizes, a nonzero plant mass–density slope is primarily defined by structural mass. Specifically, the slope is a result of (i) the physical size-dependent relationship between stem width and height, (ii) foliage-dependent demands of conductance, and (iii) the cumulative nature of structural mass. The data support our model, both when the potential sampling bias of taxonomic relatedness is accounted for and when it is not. Independent contrasts analyses show that observed relationships among variables are not significantly different from the assumptions made to build the model or from its a priori predictions. We note that the dependence of the plant mass–density slope on the functions of structural mass provides a cause for the difference from the zero slope found in the animal population mass–density relationship; for the most part, animals do not have a comparable cumulative tissue type.

Localization of the C terminus of the assembly domain of hepatitis B virus capsid protein: Implications for morphogenesis and organization of encapsidated RNA

The capsid protein of hepatitis B virus, consisting of an “assembly” domain (residues 1–149) and an RNA-binding “protamine” domain (residues 150–183), assembles from dimers into icosahedral capsids of two different sizes. The C terminus of the assembly domain (residues 140–149) functions as a morphogenetic switch, longer C termini favoring a higher proportion of the larger capsids, it also connects the protamine domain to the capsid shell. We now have defined the location of this peptide in capsids assembled in vitro by engineering a mutant assembly domain with a single cysteine at its C terminus (residue 150), labeling it with a gold cluster and visualizing the cluster by cryo-electron microscopy. The labeled protein is unimpaired in its ability to form capsids. Our density map reveals a single undecagold cluster under each fivefold and quasi-sixfold vertex, connected to sites at either end of the undersides of the dimers. Considering the geometry of the vertices, the C termini must be more crowded at the fivefolds. Thus, a bulky C terminus would be expected to favor formation of the larger (T = 4) capsids, which have a greater proportion of quasi-sixfolds. Capsids assembled by expressing the full-length protein in Escherichia coli package bacterial RNAs in amounts equivalent to the viral pregenome. Our density map of these capsids reveals a distinct inner shell of density—the RNA. The RNA is connected to the protein shell via the C-terminal linkers and also makes contact around the dimer axes.

Quantification of the hydrophobic interaction by simulations of the aggregation of small hydrophobic solutes in water

The hydrophobic interaction, the tendency for nonpolar molecules to aggregate in solution, is a major driving force in biology. In a direct approach to the physical basis of the hydrophobic effect, nanosecond molecular dynamics simulations were performed on increasing numbers of hydrocarbon solute molecules in water-filled boxes of different sizes. The intermittent formation of solute clusters gives a free energy that is proportional to the loss in exposed molecular surface area with a constant of proportionality of 45 ± 6 cal/mol⋅Å2. The molecular surface area is the envelope of the solute cluster that is impenetrable by solvent and is somewhat smaller than the more traditional solvent-accessible surface area, which is the area transcribed by the radius of a solvent molecule rolled over the surface of the cluster. When we apply a factor relating molecular surface area to solvent-accessible surface area, we obtain 24 cal/mol⋅Å2. Ours is the first direct calculation, to our knowledge, of the hydrophobic interaction from molecular dynamics simulations; the excellent qualitative and quantitative agreement with experiment proves that simple van der Waals interactions and atomic point-charge electrostatics account for the most important driving force in biology.

Expression of the Yeast FRE Genes in Transgenic Tobacco

Two yeast genes, FRE1 and FRE2 (encoding Fe(III) reductases) were placed under the control of the cauliflower mosaic virus 35S promoter and introduced into tobacco (Nicotiana tabacum L.) via Agrobacterium tumefaciens-mediated transformation. Homozygous lines containing FRE1, FRE2, or FRE1 plus FRE2 were generated. Northern-blot analyses revealed mRNA of two different sizes in FRE1 lines, whereas all FRE2 lines had mRNA only of the expected length. Fe(III) reduction, chlorophyll contents, and Fe levels were determined in transgenic and control plants under Fe-sufficient and Fe-deficient conditions. In a normal growth environment, the highest root Fe(III) reduction, 4-fold higher than in controls, occurred in the double transformant (FRE1 + FRE2). Elevated Fe(III) reduction was also observed in all FRE2 and some FRE1 lines. The increased Fe(III) reduction occurred along the entire length of the roots and on shoot sections. FRE2 and double transformants were more tolerant to Fe deficiency in hydroponic culture, as shown by higher chlorophyll and Fe concentrations in younger leaves, whereas FRE1 transformants did not differ from the controls. Overall, the beneficial effects of FRE2 were consistent, suggesting that FRE2 may be used to improve Fe efficiency in crop plants.

Dynamic friction and the origin of the complexity of earthquake sources.

We study a simple antiplane fault of finite length embedded in a homogeneous isotropic elastic solid to understand the origin of seismic source heterogeneity in the presence of nonlinear rate- and state-dependent friction. All the mechanical properties of the medium and friction are assumed homogeneous. Friction includes a characteristic length that is longer than the grid size so that our models have a well-defined continuum limit. Starting from a heterogeneous initial stress distribution, we apply a slowly increasing uniform stress load far from the fault and we simulate the seismicity for a few 1000 events. The style of seismicity produced by this model is determined by a control parameter associated with the degree of rate dependence of friction. For classical friction models with rate-independent friction, no complexity appears and seismicity is perfectly periodic. For weakly rate-dependent friction, large ruptures are still periodic, but small seismicity becomes increasingly nonstationary. When friction is highly rate-dependent, seismicity becomes nonperiodic and ruptures of all sizes occur inside the fault. Highly rate-dependent friction destabilizes the healing process producing premature healing of slip and partial stress drop. Partial stress drop produces large variations in the state of stress that in turn produce earthquakes of different sizes. Similar results have been found by other authors using the Burridge and Knopoff model. We conjecture that all models in which static stress drop is only a fraction of the dynamic stress drop produce stress heterogeneity.

Three distinct structural environments of a transmembrane domain in the inwardly rectifying potassium channel ROMK1 defined by perturbation.

To probe the protein environment of an ion channel, we have perturbed the structure of a transmembrane domain by substituting side chains with those of two different sizes by using site-specific mutagenesis. We have used Trp and Ala as a high- and a low-impact perturbation probe, respectively, to replace each of 18 consecutive residues within the putative second transmembrane segment, M2, of an inwardly rectifying potassium channel, ROMK1. Our rationale is that a change in the channel function as a consequence of these mutations at a particular position will reflect the structural environment of the altered side chain. Each position can then be assigned to one of three classes of environments, as grated by different levels of perturbation: very tolerant (channel functions with both Trp and Ala substitutions), tolerant (function preserved with Ala but not with Trp substitution), and intolerant (either Ala or Trp substitution destroys function). We identify the very tolerant environment as being lipid-facing, tolerant as protein-interior-facing, and intolerant as pore-facing. We observe a strikingly ordered pattern of perturbation of all three environmental classes. This result indicates that M2 is a straight alpha-helix.

Ectopic expression of the minimal whiE polyketide synthase generates a library of aromatic polyketides of diverse sizes and shapes

The single recombinant expressing the Streptomyces coelicolor minimal whiE (spore pigment) polyketide synthase (PKS) is uniquely capable of generating a large array of well more than 30 polyketides, many of which, so far, are novel to this recombinant. The characterized polyketides represent a diverse set of molecules that differ in size (chain length) and shape (cyclization pattern). This combinatorial biosynthetic library is, by far, the largest and most complex of its kind described to date and indicates that the minimal whiE PKS does not independently control polyketide chain length nor dictate the first cyclization event. Rather, the minimal PKS enzyme complex must rely on the stabilizing effects of additional subunits (i.e., the cyclase whiE-ORFVI) to ensure that the chain reaches the full 24 carbons and cyclizes correctly. This dramatic loss of control implies that the growing polyketide chain does not remain enzyme bound, resulting in the spontaneous cyclization of the methyl terminus. Among the six characterized dodecaketides, four different first-ring cyclization regiochemistries are represented, including C7/C12, C8/C13, C10/C15, and C13/C15. The dodecaketide TW93h possesses a unique 2,4-dioxaadamantane ring system and represents a new structural class of polyketides with no related structures isolated from natural or engineered organisms, thus supporting the claim that engineered biosynthesis is capable of producing novel chemotypes.

Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.