Semiparametric estimation of treatment effect in a pretest-posttest study

Inference on treatment effects in a pretest-posttest study is a routine objective in medicine, public health, and other fields. A number of approaches have been advocated. We take a semiparametric perspective, making no assumptions about the distributions of baseline and posttest responses. By representing the situation in terms of counterfactual random variables, we exploit recent developments in the literature on missing data and causal inference, to derive the class of all consistent treatment effect estimators, identify the most efficient such estimator, and outline strategies for implementation of estimators that may improve on popular methods. We demonstrate the methods and their properties via simulation and by application to a data set from an HIV clinical trial.     

The molecular basis of the temperature- and pH-induced conformational transitions in elastin-based peptides

Elastin undergoes an inverse temperature transition and collapses at high temperatures in both simulation and experiment. We investigated a pH-dependent modification of this transition by simulating a glutamic acid (Glu)-substituted elastin at varying pHs and temperatures. The Glu-substituted peptide collapsed at higher temperature than the unsubstituted elastin when Glu was charged. The charge effects could be reversed by neutralization of the Glu carboxyl groups at low pH, and in that case the peptide collapsed at a lower temperature. The collapse was accompanied by the formation of beta-turns and short distorted beta-sheets. Formation of contacts between hydrophobic side chains drives the collapse at high temperature, but interactions between water and polar groups (Glu and main chain) can attenuate this effect at high pH. The overall competition and balance of the polar and nonpolar groups determined the conformational states of the peptide. Water hydration contributed to the conformational transition, and the peptide and its hydration shell must be considered. Structurally, waters near polar residues mainly formed hydrogen bonds with the protein atoms, while waters around the hydrophobic side chains tended to be parallel to the peptide groups to maximize water-water interactions.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Detecting differential gene expression with a semiparametric hierarchical mixture method

Mixture modeling provides an effective approach to the differential expression problem in microarray data analysis. Methods based on fully parametric mixture models are available, but lack of fit in some examples indicates that more flexible models may be beneficial. Existing, more flexible, mixture models work at the level of one-dimensional gene-specific summary statistics, and so when there are relatively few measurements per gene these methods may not provide sensitive detectors of differential expression. We propose a hierarchical mixture model to provide methodology that is both sensitive in detecting differential expression and sufficiently flexible to account for the complex variability of normalized microarray data. EM-based algorithms are used to fit both parametric and semiparametric versions of the model. We restrict attention to the two-sample comparison problem; an experiment involving Affymetrix microarrays and yeast translation provides the motivating case study. Gene-specific posterior probabilities of differential expression form the basis of statistical inference; they define short gene lists and false discovery rates. Compared to several competing methodologies, the proposed methodology exhibits good operating characteristics in a simulation study, on the analysis of spike-in data, and in a cross-validation calculation.     

Prediction and classification of protein subcellular location-sequence-order effect and pseudo amino acid composition

Given a protein sequence, how to identify its subcellular location? With the rapid increase in newly found protein sequences entering into databanks, the problem has become more and more important because the function of a protein is closely correlated with its localization. To practically deal with the challenge, a dataset has been established that allows the identification performed among the following 14 subcellular locations: (1) cell wall, (2) centriole, (3) chloroplast, (4) cytoplasm, (5) cytoskeleton, (6) endoplasmic reticulum, (7) extracellular, (8) Golgi apparatus, (9) lysosome, (10) mitochondria, (11) nucleus, (12) peroxisome, (13) plasma membrane, and (14) vacuole. Compared with the datasets constructed by the previous investigators, the current one represents the largest in the scope of localizations covered, and hence many proteins which were totally out of picture in the previous treatments, can now be investigated. Meanwhile, to enhance the potential and flexibility in taking into account the sequence‐order effect, the series‐mode pseudo‐amino‐acid‐composition has been introduced as a representation for a protein. High success rates are obtained by the re‐substitution test, jackknife test, and independent dataset test, respectively. It is anticipated that the current automated method can be developed to a high throughput tool for practical usage in both basic research and pharmaceutical industry. © 2003 Wiley‐Liss, Inc.     

Comparison of Bayesian and maximum likelihood bootstrap measures of phylogenetic reliability

Owing to the exponential growth of genome databases, phylogenetic trees are now widely used to test a variety of evolutionary hypotheses. Nevertheless, computation time burden limits the application of methods such as maximum likelihood nonparametric bootstrap to assess reliability of evolutionary trees. As an alternative, the much faster Bayesian inference of phylogeny, which expresses branch support as posterior probabilities, has been introduced. However, marked discrepancies exist between nonparametric bootstrap proportions and Bayesian posterior probabilities, leading to difficulties in the interpretation of sometimes strongly conflicting results. As an attempt to reconcile these two indices of node reliability, we apply the nonparametric bootstrap resampling procedure to the Bayesian approach. The correlation between posterior probabilities, bootstrap maximum likelihood percentages, and bootstrapped posterior probabilities was studied for eight highly diverse empirical data sets and were also investigated using experimental simulation. Our results show that the relation between posterior probabilities and bootstrapped maximum likelihood percentages is highly variable but that very strong correlations always exist when Bayesian node support is estimated on bootstrapped character matrices. Moreover, simulations corroborate empirical observations in suggesting that, being more conservative, the bootstrap approach might be less prone to strongly supporting a false phylogenetic hypothesis. Thus, apparent conflicts in topology recovered by the Bayesian approach were reduced after bootstrapping. Both posterior probabilities and bootstrap supports are of great interest to phylogeny as potential upper and lower bounds of node reliability, but they are surely not interchangeable and cannot be directly compared.     

Pedigree analysis of eight Spanish beef cattle breeds

The genetic structure of eight Spanish autochthonous populations (breeds) of beef cattle were studied from pedigree records. The populations studied were: Alistana and Sayaguesa (minority breeds), Avileña – Negra Ibérica and Morucha ("dehesa" breeds, with a scarce incidence of artificial insemination), and mountain breeds, including Asturiana de los Valles, Asturiana de la Montaña and Pirenaica, with extensive use of AI. The Bruna dels Pirineus breed possesses characteristics which make its classification into one of the former groups difficult. There was a large variation between breeds both in the census and the number of herds. Generation intervals ranged from 3.7 to 5.5 years, tending to be longer as the population size was larger. The effective numbers of herds suggest that a small number of herds behaves as a selection nucleus for the rest of the breed. The complete generation equivalent has also been greatly variable, although in general scarce, with the exception of the Pirenaica breed, with a mean of 3.8. Inbreeding effective population sizes were actually small (21 to 127), especially in the mountain-type breeds. However, the average relatedness computed for these breeds suggests that a slight exchange of animals between herds will lead to a much more favourable evolution of inbreeding. The effective number of founders and ancestors were also variable among breeds, although in general the breeds behaved as if they were founded by a small number of animals (25 to 163).