Estimating biomarker-based HIV incidence using prevalence data in high risk groups with missing outcomes

The novel two-step serologic sensitive/less sensitive testing algorithm for detecting recent HIV seroconversion (STARHS) provides a simple and practical method to estimate HIV-1 incidence using cross-sectional HIV seroprevalence data. STARHS has been used increasingly in epidemiologic studies. However, the uncertainty of incidence estimates using this algorithm has not been well described, especially for high risk groups or when missing data is present because a fraction of sensitive enzyme immunoassay (EIA) positive specimens are not tested by the less sensitive EIA. Ad hoc methods used in practice provide incorrect confidence limits and thus may jeopardize statistical inference. In this report, we propose maximum likelihood and Bayesian methods for correctly estimating the uncertainty in incidence estimates obtained using prevalence data with a fraction missing, and extend the methods to regression settings. Using a study of injection drug users participating in a drug detoxification program in New York city as an example, we demonstrated the impact of underestimating the uncertainty in incidence estimates using ad hoc methods. Our methods can be applied to estimate the incidence of other diseases from prevalence data using similar testing algorithms when missing data is present.     

Investigation of the Role of a Conserved Glycine Motif in the Saccharomyces cerevisiae Xylose Reductase

All yeast xylose reductases, with the exception of that from Schizosaccharomyces pombe, possess the catalytic and coenzyme-binding elements from both the aldo–keto reductase and short-chain dehydrogenase–reductase (SDR) enzyme families in their primary sequences. In the Saccharomyces cerevisiae xylose reductase (XR), the SDR-like coenzyme-binding GXXXGXG motif (Gly motif) is located between residues 128 and 134, with the third Gly residue being replaced by an Asp. We used site-directed mutagenesis to study the role of this SDR-like Gly motif in the S. cerevisiae xylose reductase. Site-directed mutagenesis of the individual conserved Gly residue positions (G128A, G132A, D134G, and D134A) did not significantly affect the specific activity, kinetic constants (K_m, K_cat, and K_cat/K_m), or dissociation constants (K_d) in any of the variants compared with the wild type. Deletion of the entire Gly motif produced an unstable protein that could not be purified. These results indicate that the SDR-like Gly motif likely provides support to the overall structure of the enzyme, but it does not contribute directly to coenzyme binding in this XR.     

Expression and Aggregation of Recombinant Human Consensus Interferon-α Mutant by Pichia pastoris

A recombinant human consensus interferon-α mutant (cIFN) was expressed in Pichia pastoris. The maximum dry cell weight, cIFN concentration and antiviral activity were 160 g l⁻¹, 1.24 g l⁻¹ and 4.1 × 10⁷ IU ml⁻¹, respec tively. The cIFN secreted into the medium was in the form of aggregates dominantly by non-covalent interaction and partially by disulphide bond. When the fermentation supernatant was disaggregated with 6 M guanidine hydrochloride, the antiviral activity of cIFN achieved 2.2 × 10⁸ IU ml⁻¹.     

Sequence analysis and functional study of thymidylate synthase from zebrafish, Danio rerio

The thymidylate synthase (TS), an important target for many anticancer drugs, has been cloned from different species. But the cDNA property and function of TS in zebrafish are not well documented. In order to use zebrafish as an animal model for screening novel anticancer agents, we isolated TS cDNA from zebrafish and compared its sequence with those from other species. The open reading frame (ORF) of zebrafish TS cDNA sequence was 954 nucleotides, encoding a 318-amino acid protein with a calculated molecular mass of 36.15 kDa. The deduced amino acid sequence of zebrafish TS was similar to those from other organisms, including rat, mouse and humans. The zebrafish TS protein was expressed in Escherichia coli and purified to homogeneity. The purified zebrafish TS showed maximal activity at 28 degrees C with similar K(m) value to human TS. Western immunoblot assay confirmed that TS was expressed in all the developmental stages of zebrafish with a high level of expression at the 1-4 cell stages. To study the function of TS in zebrafish embryo development, a short hairpin RNA (shRNA) expression vector, pSilencer 4.1-CMV/TS, was constructed which targeted the protein-coding region of zebrafish TS mRNA. Significant change in the development of tail and epiboly was found in zebrafish embryos microinjected pSilencer4.1-CMV/TS siRNA expression vector.     

Case-control and case-only designs with genotype and family history data: estimating relative risk, residual familial aggregation, and cumulative risk

In case-control studies of inherited diseases, participating subjects (probands) are often interviewed to collect detailed data about disease history and age-at-onset information in their family members. Genotype data are typically collected from the probands, but not from their relatives. In this article, we introduce an approach that combines case-control analysis of data on the probands with kin-cohort analysis of disease history data on relatives. Assuming a marginally specified multivariate survival model for joint risk of disease among family members, we describe methods for estimating relative risk, cumulative risk, and residual familial aggregation. We also describe a variation of the methodology that can be used for kin-cohort analysis of the family history data from a sample of genotyped cases only. We perform simulation studies to assess performance of the proposed methodologies with correct and mis-specified models for familial aggregation. We illustrate the proposed methodologies by estimating the risk of breast cancer from BRCA1/2 mutations using data from the Washington Ashkenazi Study.     

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.     

Role of nuclear chromogranin B in inositol 1,4,5-trisphosphate-mediated nuclear Ca2+ mobilization

Recently, secretory granule Ca(2+) storage protein chromogranin B (CGB) was shown to be present in the nucleoplasm proper in a complex structure that consists of the inositol 1,4,5-trisphosphate receptor (IP(3)R)/Ca(2+) channels and the phospholipids. Further, the amounts of IP(3)Rs present in the nucleus of bovine chromaffin cells were shown to be comparable to that of the endoplasmic reticulum. Therefore, we investigated here the potential contribution of nuclear CGB on the IP(3)-dependent Ca(2+) mobilization in the nucleus, using both neuroendocrine PC12 and nonneuroendocrine NIH3T3 cells. Chromogranin A (CGA) expression in the NIH3T3 cells, which do not contain intrinsic chromogranins, increased the IP(3)-induced Ca(2+) releases in the nucleus by 45%, while CGB expression in the same cells increased the IP(3)-induced Ca(2+) releases in the nucleus by 80%. Microinjection of IP(3) into the nucleus of CGB-expressing NIH3T3 cells increased the IP(3)-dependent nuclear Ca(2+) mobilization approximately 3-fold, whereas in CGA-expressing cells it remained the same as that of control cells. In contrast, inhibition of CGA expression in PC12 cells by siRNA treatment decreased the IP(3)-induced Ca(2+) releases in the nucleus by 17%, while inhibition of CGB expression decreased the IP(3)-induced Ca(2+) releases in the nucleus by 55%. Microinjection of IP(3) into the nucleus of siCGB-treated PC12 cells decreased the IP(3)-dependent nuclear Ca(2+) mobilization by approximately 75%, whereas in siCGA-treated cells it remained the same as that of control cells. Given the presence of CGB in the nucleus, these results further highlight the critical contribution of nuclear CGB in the IP(3)-induced Ca(2+) release in the nucleus.