On a resampling approach for tests on the number of clusters with mixture model-based clustering of tissue samples

We consider the problem of assessing the number of clusters in a limited number of tissue samples containing gene expressions for possibly several thousands of genes. It is proposed to use a normal mixture model-based approach to the clustering of the tissue samples. One advantage of this approach is that the question on the number of clusters in the data can be formulated in terms of a test on the smallest number of components in the mixture model compatible with the data. This test can be carried out on the basis of the likelihood ratio test statistic, using resampling to assess its null distribution. The effectiveness of this approach is demonstrated on simulated data and on some microarray datasets, as considered previously in the bioinformatics literature. (C) 2004 Elsevier Inc. All rights reserved.

Tractable molecular theory of transport of Lennard-Jones fluids in nanopores

We present here a tractable theory of transport of simple fluids in cylindrical nanopores, which is applicable over a wide range of densities and pore sizes. In the Henry law low-density region the theory considers the trajectories of molecules oscillating between diffuse wall collisions, while at higher densities beyond this region the contribution from viscous flow becomes significant and is included through our recent approach utilizing a local average density model. The model is validated by means of equilibrium as well nonequilibrium molecular dynamics simulations of supercritical methane transport in cylindrical silica pores over a wide range of temperature, density, and pore size. The model for the Henry law region is exact and found to yield an excellent match with simulations at all conditions, including the single-file region of very small pore size where it is shown to provide the density-independent collective transport coefficient. It is also shown that in the absence of dispersive interactions the model reduces to the classical Knudsen result, but in the presence of such interactions the latter model drastically overpredicts the transport coefficient. For larger micropores beyond the single-file region the transport coefficient is reduced at high density because of intermolecular interactions and hindrance to particle crossings leading to a large decrease in surface slip that is not well represented by the model. However, for mesopores the transport coefficient increases monotonically with density, over the range studied, and is very well predicted by the theory, though at very high density the contribution from surface slip is slightly overpredicted. It is also seen that the concept of activated diffusion, commonly associated with diffusion in small pores, is fundamentally invalid for smooth pores, and the apparent activation energy is not simply related to the minimum pore potential or the adsorption energy as generally assumed. (C) 2004 American Institute of Physics.

Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation. (C) 2004 American Institute of Physics.

Coherent molecular bound states of bosons and fermions near a Feshback resonance

We analyze molecular bound states of atomic quantum gases near a Feshbach resonance. A simple, renormalizable field theoretic model is shown to have exact solutions in the two-body sector, whose binding energy agrees well with observed experimental results in both Bosonic and Fermionic cases. These solutions, which interpolate between BEC and BCS theories, also provide a more general variational ansatz for resonant superfluidity and related problems.

Enzymatic characterization and homology model of a catalytically active recombinant West Nile virus NS3 protease

West Nile Virus (WNV) is a mosquito-borne flavivirus with a rapidly expanding global distribution. Infection causes severe neurological disease and fatalities in both human and animal hosts. The West Nile viral protease (NS2B-NS3) is essential for post-translational processing in host-infected cells of a viral polypeptide precursor into structural and functional viral proteins, and its inhibition could represent a potential treatment for viral infections. This article describes the design, expression, and enzymatic characterization of a catalytically active recombinant WNV protease, consisting of a 40-residue component of cofactor NS2B tethered via a noncleavable nonapeptide (G(4)SG(4)) to the N-terminal 184 residues of NS3. A chromogenic assay using synthetic para-nitroanilide (pNA) hexapeptide substrates was used to identify optimal enzyme-processing conditions (pH 9.5, I < 0.1 M, 30% glycerol, 1 mM CHAPS), preferred substrate cleavage sites, and the first competitive inhibitor (Ac-FASGKR- H, IC50 &SIM; 1 μM). A putative three-dimensional structure of WNV protease, created through homology modeling based on the crystal structures of Dengue-2 and Hepatitis C NS3 viral proteases, provides some valuable insights for structure-based design of potent and selective inhibitors of WNV protease.

Use of multiple biomarkers for a molecular diagnosis of prostate cancer

The identification of biomarkers capable of providing a reliable molecular diagnostic test for prostate cancer (PCa) is highly desirabie clinically. We describe here 4 biomarkers, UDP-N-Acetyl-alpha-D-galactosamine transferase (GalNAc-T3; not previously associated with PCa), PSMA, Hepsin and DD3/PCA3, which, in combination, distinguish prostate cancer from benign prostate hyperplasia (BPH). GalNAc-T3 was identified as overexpressed in PCa tissues by microarray analysis, confirmed by quantitative real-time PCR and shown immunohistochemically to be localised to prostate epithelial cells with higher expression in malignant cells. Real-time quantitative PCR analysis across 21 PCa and 34 BPH tissues showed 4.6-fold overexpression of GalNAc-T3 (p = 0.005). The noncoding mRNA (DD3/PCA3) was overexpressed 140-fold (p = 0.007) in the cancer samples compared to BPH tissues. Hepsin was overexpressed 21-fold (p = 0.049, whereas the overexpression for PSMA was 66-fold (p = 0.047). When the gene expression data for these 4 biomarkers was combined in a logistic regression model, a predictive index was obtained that distinguished 100% of the PCa samples from all of the BPH samples. Therefore, combining these genes in a real-time PCR assay represents a powerful new approach to diagnosing PCa by molecular profiling. (c) 2005 Wiley-Liss, Inc.

Model for growth hormone receptor activation based on subunit rotation within a receptor dimer

Growth hormone is believed to activate the growth hormone receptor (GHR) by dimerizing two identical receptor subunits, leading to activation of JAK2 kinase associated with the cytoplasmic domain. However, we have reported previously that dimerization alone is insufficient to activate full-length GHR. By comparing the crystal structure of the liganded and unliganded human GHR extracellular domain, we show here that there is no substantial change in its conformation on ligand binding. However, the receptor can be activated by rotation without ligand by inserting a defined number of alanine residues within the transmembrane domain. Fluorescence resonance energy transfer ( FRET), bioluminescence resonance energy transfer (BRET) and coimmunoprecipitation studies suggest that receptor subunits undergo specific transmembrane interactions independent of hormone binding. We propose an activation mechanism involving a relative rotation of subunits within a dimeric receptor as a result of asymmetric placement of the receptor-binding sites on the ligand.

Design strategies for controlling the molecular weight and rate using reversible addition-fragmentation chain transfer mediated living radical polymerization

Living radical polymerization has allowed complex polymer architectures to be synthesized in bulk, solution, and water. The most versatile of these techniques is reversible addition-fragmentation chain transfer (RAFT), which allows a wide range of functional and nonfunctional polymers to be made with predictable molecular weight distributions (MWDs), ranging from very narrow to quite broad. The great complexity of the RAFT mechanism and how the kinetic parameters affect the rate of polymerization and MWD are not obvious. Therefore, the aim of this article is to provide useful insights into the important kinetic parameters that control the rate of polymerization and the evolution of the MWD with conversion. We discuss how a change in the chain-transfer constant can affect the evolution of the MWD. It is shown how we can, in principle, use only one RAFT agent to obtain a poly-mer with any MWD. Retardation and inhibition are discussed in terms of (1) the leaving R group reactivity and (2) the intermediate radical termination model versus the slow fragmentation model. (c) 2005 Wiley Periodicals, Inc.

Studies on the membrane interactions of the cyclotides kalata B1 and kalata B6 on model membrane systems by surface plasmon resonance

In this study we have demonstrated the interactions of kalata B1 and its naturally occurring analogue kalata B6 with five model lipid membranes and have analyzed the binding kinetics using surface plasmon resonance. Two kalata peptides showed a higher affinity for the phosphatidylethanolamine-containing membranes, indicating that the peptides would bind selectively to bacterial membranes. Also we have optimized the procedure for the immobilization of five liposome mixtures and have shown that the procedure provides reproducible levels of immobilized liposomes and could be used to screen the selective binding of putative antimicrobial peptides to model mammalian or microbial phospholipid membranes. (C) 2004 Elsevier Inc. All rights reserved.

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.

The pathophysiology of 'brain shrinkage' in alcoholics - Structural and molecular changes and clinical implications

This article represents a symposium of the 2004 ISBRA Congress held in Mannheim. The presentations were: Review of the neuropathological and neurochemical changes seen in alcohol-related ' brain shrinkage ' by Clive Harper; In Vivo Detection of Macrostructural and Microstructural Markers of Brain Integrity in Human Alcoholism and a Rodent Model of Alcoholism by Adolf Pfefferbaum, Elfar Adalsteinsson and Edith Sullivan; Gene and Protein Changes in the Brains of Alcoholics with ' Brain Shrinkage ' by Joanne Lewohl and Peter Dodd; Cross sectional and longitudinal MR spectroscopy studies of chronic adult alcoholics by Michael Taylor; Brain Atrophy Associated with Impairment on a Simulated Gambling Task in Long-Term Abstinent Alcoholics by George Fein and Bennett Landman.

Adaptation is not required to explain the long-term response of axons to molecular gradients

It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.

Zebrafish as a model for caveolin-associated muscle disease; caveolin-3 is required for myofibril organization and muscle cell patterning

Caveolae are an abundant feature of many animal cells. However, the exact function of caveolae remains unclear. We have used the zebrafish, Danio rerio, as a system to understand caveolae function focusing on the muscle-specific caveolar protein, caveolin-3 (Cav3). We have identified caveolin-1 (alpha and beta), caveolin-2 and Cav3 in the zebrafish. Zebrafish Cav3 has 72% identity to human CAV3, and the amino acids altered in human muscle diseases are conserved in the zebrafish protein. During embryonic development, cav3 expression is apparent by early segmentation stages in the first differentiating muscle precursors, the adaxial cells and slightly later in the notochord. cav3 expression appears in the somites during mid-segmentation stages and then later in the pectoral fins and facial muscles. Cav3 and caveolae are located along the entire sarcolemma of late stage embryonic muscle fibers, whereas beta-dystroglycan is restricted to the muscle fiber ends. Down-regulation of Cav3 expression causes gross muscle abnormalities and uncoordinated movement. Ultrastructural analysis of isolated muscle fibers reveals defects in myoblast fusion and disorganized myofibril and membrane systems. Expression of the zebrafish equivalent to a human muscular dystrophy mutant, CAV3P104L, causes severe disruption of muscle differentiation. In addition, knockdown of Cav3 resulted in a dramatic up-regulation of eng1a expression resulting in an increase in the number of muscle pioneer-like cells adjacent to the notochord. These studies provide new insights into the role of Cav3 in muscle development and demonstrate its requirement for correct intracellular organization and myoblast fusion.