PCA-correlated SNPs for structure identification in worldwide human populations

Existing methods to ascertain small sets of markers for the identification of human population structure require prior knowledge of individual ancestry. Based on Principal Components Analysis (PCA), and recent results in theoretical computer science, we present a novel algorithm that, applied on genomewide data, selects small subsets of SNPs (PCA-correlated SNPs) to reproduce the structure found by PCA on the complete dataset, without use of ancestry information. Evaluating our method on a previously described dataset (10,805 SNPs, 11 populations), we demonstrate that a very small set of PCA-correlated SNPs can be effectively employed to assign individuals to particular continents or populations, using a simple clustering algorithm. We validate our methods on the HapMap populations and achieve perfect intercontinental differentiation with 14 PCA-correlated SNPs. The Chinese and Japanese populations can be easily differentiated using less than 100 PCA-correlated SNPs ascertained after evaluating 1.7 million SNPs from HapMap. We show that, in general, structure informative SNPs are not portable across geographic regions. However, we manage to identify a general set of 50 PCA-correlated SNPs that effectively assigns individuals to one of nine different populations. Compared to analysis with the measure of informativeness, our methods, although unsupervised, achieved similar results. We proceed to demonstrate that our algorithm can be effectively used for the analysis of admixed populations without having to trace the origin of individuals. Analyzing a Puerto Rican dataset (192 individuals, 7,257 SNPs), we show that PCA-correlated SNPs can be used to successfully predict structure and ancestry proportions. We subsequently validate these SNPs for structure identification in an independent Puerto Rican dataset. The algorithm that we introduce runs in seconds and can be easily applied on large genome-wide datasets, facilitating the identification of population substructure, stratification assessment in multi-stage whole-genome association studies, and the study of demographic history in human populations.     

Spatiotemporal distribution of Fras1/Frem proteins during mouse embryonic development

The Fras1/Frem gene family encodes for structurally similar, developmentally regulated extracellular matrix proteins. Mutations in Fras1, Frem1 and Frem2 have been identified in different classes of mouse bleb mutants, while defects in the human orthologs FRAS1 and FREM2 are causative for Fraser syndrome. The hallmark phenotypic feature of bleb mice is embryonic skin blistering due to dermal-epidermal detachment. The similarity of the phenotypic characteristics among the bleb mouse mutants, together with the fact that Fras1/Frem proteins are co-localized in embryonic epithelial basement membranes, suggest that they operate in a common pathway. Here, we report for the first time the immunofluorescence pattern of Frem3 and provide a comparative analysis of the spatiotemporal localization of all Fras1/Frem proteins during mouse embryonic development. We demonstrate their overall co-localization in embryonic epithelial basement membranes, with emphasis on areas of phenotypic interest such as eyelids, limbs, kidneys, lungs and organs of the gastrointestinal tract and the central nervous system. We further studied collagen VII, impairment of which produces dystrophic epidermolysis bullosa, a postnatal skin blistering disorder. We show that basement membrane levels of collagen VII rise at late embryonic life, concomitant with descending Fras1/Frem immunolabeling.     

Optic chiasm presentation of Semaphorin6D in the context of Plexin-A1 and Nr-CAM promotes retinal axon midline crossing

At the optic chiasm, retinal ganglion cells (RGCs) project ipsi- or contralaterally to establish the circuitry for binocular vision. Ipsilateral guidance programs have been characterized, but contralateral guidance programs are not well understood. Here, we identify a tripartite molecular system for contralateral RGC projections: Semaphorin6D (Sema6D) and Nr-CAM are expressed on midline radial glia and Plexin-A1 on chiasm neurons, and Plexin-A1 and Nr-CAM are also expressed on contralateral RGCs. Sema6D is repulsive to contralateral RGCs, but Sema6D in combination with Nr-CAM and Plexin-A1 converts repulsion to growth promotion. Nr-CAM functions as a receptor for Sema6D. Sema6D, Plexin-A1, and Nr-CAM are all required for efficient RGC decussation at the optic chiasm. These findings suggest a mechanism by which a complex of Sema6D, Nr-CAM, and Plexin-A1 at the chiasm midline alters the sign of Sema6D and signals Nr-CAM/Plexin-A1 receptors on RGCs to implement the contralateral RGC projection.     

Simulations of (an)isotropic diffusion on curved biological surfaces

We present a computational particle method for the simulation of isotropic and anisotropic diffusion on curved biological surfaces that have been reconstructed from image data. The method is capable of handling surfaces of high curvature and complex shape, which are often encountered in biology. The method is validated on simple benchmark problems and is shown to be second-order accurate in space and time and of high parallel efficiency. It is applied to simulations of diffusion on the membrane of endoplasmic reticula (ER) in live cells. Diffusion simulations are conducted on geometries reconstructed from real ER samples and are compared to fluorescence recovery after photobleaching experiments in the same ER samples using the transmembrane protein tsO45-VSV-G, C-terminally tagged with green fluorescent protein. Such comparisons allow derivation of geometry-corrected molecular diffusion constants for membrane components from fluorescence recovery after photobleaching data. The results of the simulations indicate that the diffusion behavior of molecules in the ER membrane differs significantly from the volumetric diffusion of soluble molecules in the lumen of the same ER. The apparent speed of recovery differs by a factor of approximately 4, even when the molecular diffusion constants of the two molecules are identical. In addition, the specific shape of the membrane affects the recovery half-time, which is found to vary by a factor of approximately 2 in different ER samples.     

Insight into the protein and solvent contributions to the reduction potentials of [4Fe–4S]2+/+ clusters: crystal structures of the Allochromatium vinosum ferredoxin variants C57A and V13G and the homologous Escherichia coli ferredoxin

The crystal structures of the C57A and V13G molecular variants of Allochromatium vinosum 2[4Fe–4S] ferredoxin (AlvinFd) and that of the homologous ferredoxin from Escherichia coli (EcFd) have been determined at 1.05-, 1.48-, and 1.65-Å resolution, respectively. The present structures combined with cyclic voltammetry studies establish clear effects of the degree of exposure of the cluster with the lowest reduction potential (cluster I) towards less negative reduction potentials (E°). This is better illustrated by V13G AlvinFd (high exposure, E° = ?594 mV) and EcFd (low exposure, E° = ?675 mV). In C57A AlvinFd, the movement of the protein backbone, as a result of replacing the noncoordinating Cys57 by Ala, leads to a +50-mV upshift of the potential of the nearby cluster I, by removal of polar interactions involving the thiolate group and adjustment of the hydrogen-bond network involving the cluster atoms. In addition, the present structures and other previously reported accurate structures of this family of ferredoxins indicate that polar interactions of side chains and water molecules with cluster II sulfur atoms, which are absent in the environment of cluster I, are correlated to the approximately 180–250 mV difference between the reduction potentials of clusters I and II. These findings provide insight into the significant effects of subtle structural differences of the protein and solvent environment around the clusters of [4Fe–4S] ferredoxins on their electrochemical properties.     

Development of a CP P NMR Broadline Simulation Methodology for Studying the Interactions of Antihypertensive AT1 Antagonist Losartan with Phospholipid Bilayers

A cross-polarization (CP) ³¹P NMR broadline simulation methodology was developed for studying the effects of drugs in phospholipids bilayers. Based on seven-parameter fittings, this methodology provided information concerning the conformational changes and dynamics effects of losartan in the polar region of the dipalmitoylphosphatidylcholine bilayers. The test molecule for this study was losartan, an antihypertensive drug known to exert its effect on AT₁ transmembrane receptors. The results were complemented and compared with those of differential scanning calorimetry, solid-state ¹³C NMR spectroscopy, Raman spectroscopy, and electron spin resonance. More specifically, these physical chemical methodologies indicated that the amphipathic losartan molecule interacts with the hydrophilic-head zone of the lipid bilayers. The CP ³¹P NMR broadline simulations showed that the lipid molecules in the bilayers containing losartan displayed greater collective tilt compared to the tilt displayed by the load-free bilayers, indicating improved packing. The Raman results displayed a decrease in the trans/gauche ratio and increased intermolecular interactions of the acyl chains in the liquid crystalline phase. Additional evidence, suggesting that losartan possibly anchors in the realm of the headgroup, was derived from upfield shift of the average chemical shift σ^iso of the ³¹P signal in the presence of losartan and from shift of the observed peak at 715 cm⁻¹ attributed to C-N stretching in the Raman spectra.     

New Policies,New Technologies: Modelling the Potential for Improved Smear Microscopy Services in Malawi

Background

To quantify the likely impact of recent WHO policy recommendations regarding smear microscopy and the introduction of appropriate low-cost fluorescence microscopy on a) case detection and b) laboratory workload.

Methodology/Principal Findings

An audit of the laboratory register in an urban hospital, Lilongwe, Malawi, and the application of a simple modelling framework. The adoption of the new definition of a smear-positive case could directly increase case detection by up to 28%. Examining Ziehl-Neelsen (ZN) sputum smears for up to 10 minutes before declaring them negative has previously been shown to increase case detection (over and above that gained by the adoption of the new case definition) by 70% compared with examination times in routine practice. Three times the number of staff would be required to adequately examine the current workload of smears using ZN microscopy. Through implementing new policy recommendations and LED-based fluorescence microscopy the current laboratory staff complement could investigate the same number of patients, examining auramine-stained smears to an extent that is equivalent to a 10 minutes ZN smear examination.

Conclusions/Significance

Combined implementation of the new WHO recommendations on smear microscopy and LED-based fluorescence microscopy could result in substantial increases in smear positive case-detection using existing human resources and minimal additional equipment.     

Purification and biophysical characterization of the core protease domain of anthrax lethal factor

Anthrax lethal toxin (LeTx) stands for the major virulence factor of the anthrax disease. It comprises a 90 kDa highly specific metalloprotease, the anthrax lethal factor (LF). LF possesses a catalytic Zn²⁺ binding site and is highly specific against MAPK kinases, thus representing the most potent native biomolecule to alter and inactivate MKK [MAPK (mitogen-activated protein kinase) kinases] signalling pathways. Given the importance of the interaction between LF and substrate for the development of anti-anthrax agents as well as the potential treatment of nascent tumours, the analysis of the structure and dynamic properties of the LF catalytic site are essential to elucidate its enzymatic properties. Here we report the recombinant expression and purification of a C-terminal part of LF (LF_672-776) that harbours the enzyme’s core protease domain. The biophysical characterization and backbone assignments (¹H, ¹³C, ¹⁵N) of the polypeptide revealed a stable, well folded structure even in the absence of Zn²⁺, suitable for high resolution structural analysis by NMR.     

Conformational dynamics of the anthrax lethal factor catalytic center

Anthrax lethal factor (LF) is a zinc-metalloprotease that together with the protective antigen constitutes anthrax lethal toxin, which is the most prominent virulence factor of the anthrax disease. The solution nuclear magnetic resonance and in silico conformational dynamics of the 105 C-terminal residues of the LF catalytic core domain in its apo form are described here. The polypeptide adopts a compact structure even in the absence of the Zn(2+) cofactor, while the 40 N-terminal residues comprising the metal ligands and residues that participate in substrate and inhibitor recognition exhibit more flexibility than the C-terminal region.     

Discovery of a potent and selective Bcl-2 inhibitor using SAR by NMR

The Bcl-2 family of proteins plays a major role in the regulation of apoptosis, or programmed cell death. Overexpression of the anti-apoptotic members of this family (Bcl-2, Bcl-x_L, and Mcl-1) can render cancer cells resistant to chemotherapeutic agents and therefore these proteins are important targets for the development of new anti-cancer agents. Here we describe the discovery of a potent, highly selective, Bcl-2 inhibitor using SAR by NMR and structure-based drug design which could serve as a starting point for the development of a Bcl-2 selective anti-cancer agent. Such an agent would potentially overcome the Bcl-x_L mediated thrombocytopenia observed with ABT-263.