Optimal prediction of underresolved dynamics

A method is presented for computing the average solution of problems that are too complicated for adequate resolution, but where information about the statistics of the solution is available. The method involves computing average derivatives by interpolation based on linear regression, and an updating of a measure constrained by the available crude information. Examples are given.

Optimal selectin-mediated rolling of leukocytes during inflammation in vivo requires intercellular adhesion molecule-1 expression

Leukocyte interactions with vascular endothelium during inflammation occur through discrete steps involving selectin-mediated leukocyte rolling and subsequent firm adhesion mediated by members of the integrin and Ig families of adhesion molecules. To identify functional synergy between selectin and Ig family members, mice deficient in both L-selectin and intercellular adhesion molecule 1 (ICAM-1) were generated. Leukocyte rolling velocities in cremaster muscle venules were increased significantly in ICAM-1-deficient mice during both trauma- and tumor necrosis factor α-induced inflammation, but rolling leukocyte flux was not reduced. Elimination of ICAM-1 expression in L-selectin-deficient mice resulted in a sharp reduction in the flux of rolling leukocytes during tumor necrosis factor α-induced inflammation. The observed differences in leukocyte rolling behavior demonstrated that ICAM-1 expression was required for optimal P- and L-selectin-mediated rolling. Increased leukocyte rolling velocities presumably translated into decreased tissue emigration because circulating neutrophil, monocyte, and lymphocyte numbers were increased markedly in L-selectin/ICAM-1-deficient mice. Furthermore, neutrophil emigration during acute peritonitis was reduced by 80% in the double-deficient mice compared with either L-selectin or ICAM-1-deficient mice. Thus, members of the selectin and Ig families function synergistically to mediate optimal leukocyte rolling in vivo, which is essential for the generation of effective inflammatory responses.

An optimal Mg2+ concentration for kinetic folding of the Tetrahymena ribozyme

Divalent metal ions, such as Mg2+, are generally required for tertiary structure formation in RNA. Although the role of Mg2+ binding in RNA-folding equilibria has been studied extensively, little is known about the role of Mg2+ in RNA-folding kinetics. In this paper, we explore the effect of Mg2+ on the rate-limiting step in the kinetic folding pathway of the Tetrahymena ribozyme. Analysis of these data reveals the presence of a Mg2+-stabilized kinetic trap that slows folding at higher Mg2+ concentrations. Thus, the Tetrahymena ribozyme folds with an optimal rate at 2 mM Mg2+, just above the concentration required for stable structure formation. These results suggest that thermodynamic and kinetic folding of RNA are cooptimized at a Mg2+ concentration that is sufficient to stabilize the folded form but low enough to avoid kinetic traps and misfolding.

Direct physical interaction between DnaG primase and DnaB helicase of Escherichia coli is necessary for optimal synthesis of primer RNA

The primase DnaG of Escherichia coli requires the participation of the replicative helicase DnaB for optimal synthesis of primer RNA for lagging strand replication. However, previous studies had not determined whether the activation of the primase or its loading on the template was accomplished by a helicase-mediated structural alteration of the single-stranded DNA or by a direct physical interaction between the DnaB and the DnaG proteins. In this paper we present evidence supporting direct interaction between the two proteins. We have mapped the surfaces of interaction on both DnaG and DnaB and show further that mutations that reduce the physical interaction also cause a significant reduction in primer synthesis. Thus, the physical interaction reported here appears to be physiologically significant.

Combined transgenic expression of α-galactosidase and α1,2-fucosyltransferase leads to optimal reduction in the major xenoepitope Galα(1,3)Gal

Hyperacute rejection of pig organs by humans involves the interaction of Galα(1,3)Gal with antibodies and complement. Strategies to reduce the amount of xenoantigen Galα(1,3)Gal were investigated by overexpression of human lysosomal α-galactosidase in cultured porcine cells and transgenic mice. The overexpression of human α-galactosidase in cultured porcine endothelial cells and COS cells resulted in a 30-fold reduction of cell surface Galα(1,3)Gal and a 10-fold reduction in cell reactivity with natural human antibodies. Splenocytes from transgenic mice overexpressing human α-galactosidase showed only a 15–25% reduction in binding to natural human anti-Galα(1,3)Gal antibodies; however, this decrease was functionally significant as demonstrated by reduced susceptibility to human antibody-mediated lysis. However, because there is residual Galα(1,3)Gal and degalactosylation results in the exposure of N-acetyllactosamine residues and potential new xenoepitopes, using α-galactosidase alone is unlikely to overcome hyperacute rejection. We previously reported that mice overexpressing human α1,2-fucosyltransferase as a transgene had ≈90% reduced Galα(1,3)Gal levels due to masking of the xenoantigen by fucosylation; we evaluated the effect of overexpressing α-galactosidase and α1,2-fucosyltransferase on Galα(1,3)Gal levels. Galα(1,3)Gal-positive COS cells expressing α1,3-galactosyltransferase, α1,2-fucosyltransferase, and α-galactosidase showed negligible cell surface staining and were not susceptible to lysis by human serum containing antibody and complement. Thus, α1,2-fucosyltransferase and α-galactosidase effectively reduced the expression of Galα(1,3)Gal on the cell surface and could be used to produce transgenic pigs with negligible levels of cell surface Galα(1,3)Gal, thereby having no reactivity with human serum and improving graft survival.

c-Abl is required for development and optimal cell proliferation in the context of p53 deficiency

The c-Abl tyrosine kinase and the p53 tumor suppressor protein interact functionally and biochemically in cellular genotoxic stress response pathways and are implicated as downstream mediators of ATM (ataxia-telangiectasia mutated). This fact led us to study genetic interactions in vivo between c-Abl and p53 by examining the phenotype of mice and cells deficient in both proteins. c-Abl-null mice show high neonatal mortality and decreased B lymphocytes, whereas p53-null mice are prone to tumor development. Surprisingly, mice doubly deficient in both c-Abl and p53 are not viable, suggesting that c-Abl and p53 together contribute to an essential function required for normal development. Fibroblasts lacking both c-Abl and p53 were similar to fibroblasts deficient in p53 alone, showing loss of the G1/S cell-cycle checkpoint and similar clonogenic survival after ionizing radiation. Fibroblasts deficient in both c-Abl and p53 show reduced growth in culture, as manifested by reduction in the rate of proliferation, saturation density, and colony formation, compared with fibroblasts lacking p53 alone. This defect could be restored by reconstitution of c-Abl expression. Taken together, these results indicate that the ATM phenotype cannot be explained solely by loss of c-Abl and p53 and that c-Abl contributes to enhanced proliferation of p53-deficient cells. Inhibition of c-Abl function may be a therapeutic strategy to target p53-deficient cells selectively.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

A resource range invariance rule for optimal offspring size predicts patterns of variability in parental phenotypes

Previous analysis of the rules regarding how much more a female should invest in a litter of size C rather than producing a litter with one more offspring revealed an invariance relationship between litter size and the range of resources per offspring in any litter size. The rule is that the range of resources per offspring should be inversely proportional to litter size. Here we present a modification of this rule that relates litter size to the total resources devoted to reproduction at that litter size. The result is that the range of resources devoted to reproduction should be the same for all litter sizes. When parental phenotypes covary linearly with resources devoted to reproduction, then those traits should also show equal ranges within each litter size category (except for litters of one). We tested this prediction by examining the range in body size (=total length) of female mosquito fish (Gambusia hubbsi) at different litter sizes. Because resources devoted to reproduction may take many forms (e.g., nest defense), this prediction may have broad applicability.

Minimal and optimal mechanisms for GroE-mediated protein folding

We have analyzed the effects of different components of the GroE chaperonin system on protein folding by using a nonpermissive substrate (i.e., one that has very low spontaneous refolding yield) for which rate data can be acquired. In the absence of GroES and nucleotides, the rate of GroEL-mediated refolding of heat- and DTT-denatured mitochondrial malate dehydrogenase was extremely low, but some three times higher than the spontaneous rate. This GroEL-mediated rate was increased 17-fold by saturating concentrations of ATP, 11-fold by ADP and GroES, and 465-fold by ATP and GroES. Optimal refolding activity was observed when the dissociation of GroES from the chaperonin complex was dramatically reduced. Although GroEL minichaperones were able to bind denatured mitochondrial malate dehydrogenase, they were ineffective in enhancing the refolding rate. The spectrum of mechanisms for GroE-mediated protein folding depends on the nature of the substrate. The minimal mechanism for permissive substrates (i.e., having significant yields of spontaneous refolding), requires only binding to the apical domain of GroEL. Slow folding rates of nonpermissive substrates are limited by the transitions between high- and low-affinity states of GroEL alone. The optimal mechanism, which requires holoGroEL, physiological amounts of GroES, and ATP hydrolysis, is necessary for the chaperonin-mediated folding of nonpermissive substrates at physiologically relevant rates under conditions in which retention of bound GroES prevents the premature release of aggregation-prone folding intermediates from the chaperonin complex. The different mechanisms are described in terms of the structural features of mini- and holo-chaperones.

The optimal measure of allelic association

Allelic association between pairs of loci is derived in terms of the association probability ρ as a function of recombination θ, effective population size N, linear systematic pressure v, and time t, predicting both ρrt, the decrease of association from founders and ρct, the increase by genetic drift, with ρt = ρrt + ρct. These results conform to the Malecot equation, with time replaced by distance on the genetic map, or on the physical map if recombination in the region is uniform. Earlier evidence suggested that ρ is less sensitive to variations in marker allele frequencies than alternative metrics for which there is no probability theory. This robustness is confirmed for six alternatives in eight samples. In none of these 48 tests was the residual variance as small as for ρ. Overall, efficiency was less than 80% for all alternatives, and less than 30% for two of them. Efficiency of alternatives did not increase when information was estimated simultaneously. The swept radius within which substantial values of ρ are conserved lies between 385 and 893 kb, but deviation of parameters between measures is enormously significant. The large effort now being devoted to allelic association has little value unless the ρ metric with the strongest theoretical basis and least sensitivity to marker allele frequencies is used for mapping of marker association and localization of disease loci.

Vaccinia virus DNA replication: two hundred base pairs of telomeric sequence confer optimal replication efficiency on minichromosome templates.

Vaccinia virus is a complex DNA virus that exhibits significant genetic and physical autonomy from the host cell. Most if not all of the functions involved in replication and transcription of the 192-kb genome are virally encoded. Although significant progress has been made in identifying trans-acting factors involved in DNA synthesis, the mechanism of genome replication has remained poorly understood. The genome is a linear duplex with covalently closed hairpin termini, and it has been presumed that sequences and/or structures within these termini are important for the initiation of genome replication. In this report we describe the construction of minichromosomes containing a central plasmid insert flanked by hairpin termini derived from the viral genome and their use as replication templates. When replication of these minichromosomes was compared with a control substrate containing synthetic hairpin termini, specificity for viral telomeres was apparent. Inclusion of > or = 200 bp from the viral telomere was sufficient to confer optimal replication efficiency, whereas 65-bp telomeres were not effective. Chimeric 200-bp telomeres containing the 65-bp terminal element and 135 bp of ectopic sequence also failed to confer efficient replication, providing additional evidence that telomere function is sequence-specific. Replication of these exogenous templates was dependent upon the viral replication machinery, was temporally coincident with viral replication, and generated covalently closed minichromosome products. These data provide compelling evidence for specificity in template recognition and utilization in vaccinia virus-infected cells.