A mixture model for occupational exposure mean testing with a limit of detection

Information from detectable exposure measurements randomly sampled from a left-truncated log-normal distribution may be used to evaluate the distribution of nondetectable values that fall below an analytic limit of detection. If the proportion of nondetects is larger than expected under log normality, alternative models to account for these unobserved data should be considered. We discuss one such model that incorporates a mixture of true zero exposures and a log-normal distribution with possible left censoring, previously considered in a different context by Moulton and Halsey (1995, Biometrics 51, 1570-1578). A particular relationship is demonstrated between maximum likelihood parameter estimates based on this mixture model and those assuming either left-truncated or left-censored data. These results emphasize the need for caution when choosing a model to fit data involving nondetectable values. A one-sided likelihood ratio test for comparing mean exposure under the mixture model to an occupational exposure limit is then developed and evaluated via simulations. An example demonstrates the potential impact of specifying an incorrect model for the nondetectable values.     

Statistical estimators of frontal sinus cross section ontogeny from very noisy data

Cross-sectional areas of human frontal sinuses in the occipitofrontal projection are a good surrogate for frontal sinus volumes. This study looks at these areas in a dataset of some 200 children and 100 adults of both sexes. As measured by planimetry of roentgenograms, the areas are extremely variable ("noisy" in a statistical sense). In fact, they appear to be distributed log-normally with quite high variance. The mean of the distribution is evidently a function of age and the variances differ by sex. After logarithmic transformation, the data are adequately fitted by one sigmoid curve for each sex. Our discussion highlights implications of this finding for the biological aspects of frontal sinuses and methodological issues in ontogenetic analysis of data so noisy.     

Solvation effects and driving forces for protein thermodynamic and kinetic cooperativity: how adequate is native-centric topological modeling?

What energetic and solvation effects underlie the remarkable two-state thermodynamics and folding/unfolding kinetics of small single-domain proteins? To address this question, we investigate the folding and unfolding of a hierarchy of continuum Langevin dynamics models of chymotrypsin inhibitor 2. We find that residue-based additive Gō-like contact energies, although native-centric, are by themselves insufficient for protein-like calorimetric two-state cooperativity. Further native biases by local conformational preferences are necessary for protein-like thermodynamics. Kinetically, however, even models with both contact and local native-centric energies do not produce simple two-state chevron plots. Thus a model protein's thermodynamic cooperativity is not sufficient for simple two-state kinetics. The models tested appear to have increasing internal friction with increasing native stability, leading to chevron rollovers that typify kinetics that are commonly referred to as non-two-state. The free energy profiles of these models are found to be sensitive to the choice of native contacts and the presumed spatial ranges of the contact interactions. Motivated by explicit-water considerations, we explore recent treatments of solvent granularity that incorporate desolvation free energy barriers into effective implicit-solvent intraprotein interactions. This additional feature reduces both folding and unfolding rates vis-à-vis that of the corresponding models without desolvation barriers, but the kinetics remain non-two-state. Taken together, our observations suggest that interaction mechanisms more intricate than simple Gō-like constructs and pairwise additive solvation-like contributions are needed to rationalize some of the most basic generic protein properties. Therefore, as experimental constraints on protein chain models, requiring a consistent account of protein-like thermodynamic and kinetic cooperativity can be more stringent and productive for some applications than simply requiring a model heteropolymer to fold to a target structure.     

Investigations into sequence and conformational dependence of backbone entropy,inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.     

间种蔬菜对甘蔗地生态环境和甘蔗生长的影响

甘蔗 (ＳａｃｃｈａｒｕｍｏｆｆｉｃｉｎａｒｕｍＬ .)是广西重要的经济作物 ,种植面积达 4 5 33× 10 6ｈａ[1] ,甘蔗的前期生长比较缓慢 ,行间有比较大的可利用空间 ,到6月下旬～ 7月上旬才进入迅速生长期并逐渐封行。有研究表明 ,在甘蔗园间种玉米对甘蔗的生长发育和产量有一定的影响[5] ,而合理间种豆科作物 (黄豆和花生 )不仅不影响甘蔗的正常生长发育 ,还可以增加甘蔗的产量[4 ,6] 。但间种作物对甘蔗园生态环境的影响未见报道 ,我们于 1997- 1999年在桂西北开展了甘蔗园间种经济效益比较高的蔬菜作物的试验 ,并对间种甘蔗园生态环…     

Determination of the torsion angles of alanine and glycine residues of model compounds of spider silk (AGG)(10) using solid-state NMR methods

Spiders synthesize several kinds of silk fibers. In the primary structure of spider silk, one of the major ampullate (dragline, frame) silks, spidroin 1, and flagelliform silk (core fibers of adhesive spiral), there are common repeated X-Gly-Gly (X = Ala, Leu, Pro, Tyr, Glu, and Arg) sequences, which are considered to be related to the elastic character of these fibers. In this paper, two dimensional spin diffusion solid-state NMR under off magic angle spinning (OMAS), ¹³C chemical shift contour plots, and Rotational Echo DOuble Resonance (REDOR) were applied to determine the torsion angles of one Ala and two kinds of Gly residues in the Ala-Gly-Gly sequence of ¹³C=O isotope-labeled (Ala-Gly-Gly)₁₀. The torsion angles were determined to be (, ) = (–90°, 150° ) within an experimental error of ±10° for each residue. This conformation is characterized as 3₁ helix which is in agreement with the structure proposed from the X-ray powder diffraction pattern of poly(Ala-Gly-Gly). The 3₁ helix of (Ala-Gly-Gly)₁₀ does not change by formic acid treatment although (Ala-Gly)₁₅ easily changes from the silk I conformation (the structure of Bombyx mori silk fibroin before spinning in the solid state) to silk II conformation (the structure of the silk fiber after spinning) by such treatment. Thus, the 3₁ helix conformation of (Ala-Gly-Gly)₁₀ is considered very stable. Furthermore, the torsion angles of the 16th Leu residue of (Leu-Gly-Gly)₁₀ were also determined as (, ) = (–90°, 150° ) and this peptide is also considered to take 3₁ helix conformation.     

Thermodynamic and structural characterization of Asn and Ala residues in the disallowed II' region of the Ramachandran plot

Residue Asn47 at position L1 of a type II' beta-turn of the alpha-spectrin SH3 domain is located in a disallowed region of the Ramachandran plot (phi = 56 +/- 12, psi = -118 +/- 17). Therefore, it is expected that replacement of Asn47 by Gly should result in a considerable stabilization of the protein. Thermodynamic analysis of the N47G and N47A mutants shows that the change in free energy is small (approximately 0.7 kcal/mol; approximately 3 kJ/mol) and comparable to that found when mutating a Gly to Ala in a alpha-helix or beta-sheet. X-ray structural analysis of these mutants shows that the conformation of the beta-turn does not change upon mutation and, therefore, that there is no relaxation of the structure, nor is there any gain or loss of interactions that could explain the small energy change. Our results indicate that the energetic definition of II' region of the Ramachandran plot (phi = 60 +/- 30, psi = -115 +/- 15) should be revised for at least Ala and Asn in structure validation and protein design.     

Applications of likelihood asymptotics for nonlinear regression in herbicide bioassays

Dose-response models are intensively used in herbicide bioassays. Despite recent advancements in the development of new herbicides, statistical analyses are commonly based on asymptotic approximations that are sometimes poor. This paper presents the use of recent results in higher order asymptotics for likelihood-based inference in nonlinear regression. The methods presented provide accurate approximation for the distribution of test statistics and for prediction limits. Analyses of the fit and measures of detection limits of the bioassays are considered, and the potential of the methods is illustrated by examples with real data.