Facial patterns in Cercopithecoidea and Hominoidea: a geometric approach

The maxillofacial and orbital compartments of the primate skull contribute to the ontogenetic and phylogenetic variability of the viscerocranium and are of crucial evolutionary relevance. As the form of organisms changes depending on endo- and exogenous factors, metrical evaluation of specific adaptations and incorporation of the results into a biological framework could be helpful in identifying valid characters for separation of taxa (e.g. family, genus, and species) and in understanding divergence and convergence. During the last two decades a morphometric "revolution" heralded by Rohlf & Marcus (1993), Adams et al. (2004) and Oxnard (2004) brought about a synthesis of traditional quantitative-morphometrical with modern methods. This approach is called "Geometric Morphometrics (GM)" and constitutes the coremethod applied here. Based on standardized photographs (in Norma frontalis), landmarks (LM) were set and two-dimensional coordinates (X, Y) recorded for the facial cranium in selected representatives of the superfamilies Cercopithecoidea and Hominoidea. The comparison of two datasets by means of factor analysis and distance computation for the complete maxillofacial complex on the one hand, and circumorbital and orbital features on the other, indicate that morphological differences between super-families and genera are valid for separating them even in a heterogeneous sample like the one presented here. Including more landmarks and therewith capturing the morph in a more complex way optimizes separation within the sample.     

Scoring protein interaction decoys using exposed residues (SPIDER): a novel multibody interaction scoring function based on frequent geometric patterns of interfacial residues

Accurate prediction of the structure of protein-protein complexes in computational docking experiments remains a formidable challenge. It has been recognized that identifying native or native-like poses among multiple decoys is the major bottleneck of the current scoring functions used in docking. We have developed a novel multibody pose-scoring function that has no theoretical limit on the number of residues contributing to the individual interaction terms. We use a coarse-grain representation of a protein-protein complex where each residue is represented by its side chain centroid. We apply a computational geometry approach called Almost-Delaunay tessellation that transforms protein-protein complexes into a residue contact network, or an undirectional graph where vertex-residues are nodes connected by edges. This treatment forms a family of interfacial graphs representing a dataset of protein-protein complexes. We then employ frequent subgraph mining approach to identify common interfacial residue patterns that appear in at least a subset of native protein-protein interfaces. The geometrical parameters and frequency of occurrence of each "native" pattern in the training set are used to develop the new SPIDER scoring function. SPIDER was validated using standard "ZDOCK" benchmark dataset that was not used in the development of SPIDER. We demonstrate that SPIDER scoring function ranks native and native-like poses above geometrical decoys and that it exceeds in performance a popular ZRANK scoring function. SPIDER was ranked among the top scoring functions in a recent round of CAPRI (Critical Assessment of PRedicted Interactions) blind test of protein-protein docking methods.     

Matching Pursuit with Asymmetric Functions for Signal Decomposition and Parameterization

The method of adaptive approximations by Matching Pursuit makes it possible to decompose signals into basic components (called atoms). The approach relies on fitting, in an iterative way, functions from a large predefined set (called dictionary) to an analyzed signal. Usually, symmetric functions coming from the Gabor family (sine modulated Gaussian) are used. However Gabor functions may not be optimal in describing waveforms present in physiological and medical signals. Many biomedical signals contain asymmetric components, usually with a steep rise and slower decay. For the decomposition of this kind of signal we introduce a dictionary of functions of various degrees of asymmetry – from symmetric Gabor atoms to highly asymmetric waveforms. The application of this enriched dictionary to Otoacoustic Emissions and Steady-State Visually Evoked Potentials demonstrated the advantages of the proposed method. The approach provides more sparse representation, allows for correct determination of the latencies of the components and removes the "energy leakage" effect generated by symmetric waveforms that do not sufficiently match the structures of the analyzed signal. Additionally, we introduced a time-frequency-amplitude distribution that is more adequate for representation of asymmetric atoms than the conventional time-frequency-energy distribution.     

Towards a calculus of biomolecular complexes at equilibrium

An overview is presented of the construction and use of algebraic partition functions to represent the equilibrium statistical mechanics of multimolecular complexes and their action within a larger regulatory network. Unlike many applications of equilibrium statistical mechanics, multimolecular complexes may operate with various subsets of their components present and connected to the others, the rest remaining in solution. Thus they are variable-structure systems. This aspect of their behavior may be accounted for by the use of 'fugacity' variables as a representation within the partition functions. Four principles are proposed by which the combinatorics of molecular complex construction can be reflected in the construction of their partition functions. The corresponding algebraic operations on partition functions are multiplication, addition, function composition and a less commonly used operation called contraction. Each has a natural interpretation in terms of probability distributions on multimolecular structures. Possible generalizations to nonequilibrium statistical mechanics are briefly discussed.     

SELDI-TOF-MS analysis of chloroquine resistant and sensitive Plasmodium falciparum strains

The resistance of the malarial parasite Plasmodium falciparum to chloroquine represents an emerging problem since neither mode of drug action nor mechanisms of resistance are fully elucidated. We describe a protein expression profiling approach by SELDI-TOF-MS as a useful tool for studying the proteome of malarial parasites. Reproducible and complex protein profiles of the P. falciparum strains K1, Dd2, HB3 and 3D7 were measured on four array types. Hierarchical clustering led to a clear separation of the two major subgroups "resistant" and "sensitive" as well as of the four parasite strains. Our study delivers sets of regulated proteins derived from extensive comparative analyses of 64 P. falciparum protein profiles. A group of 12 peaks reflecting proteome changes under chloroquine treatment and a set of 10 potential chloroquine resistance markers were defined. Three of these regulated peaks were preparatively enriched, purified and identified. They were shown to represent the plasmodial EXP-1 protein, also called circumsporozoite-related antigen, as well as the alpha- and beta- (delta-) chains of human hemoglobin.     

A new model of developmental constraints as applied to the Drosophila system

Von Baer's laws of development observe that an embryo, in the course of its ontogeny, progresses through a series of forms which diverge increasingly from the embryonic forms of related species, and in an evolutionary interpretation, from those of its phylogenetic ancestors. This observation on the relation of phylogeny to ontogeny is explained by Wimsatt's (1986) "Developmental Lock" model of complex generative systems, which proposes that evolution is constrained to alter developmental programs in a manner that usually modifies or adds new complexity to pre-existent developmental functions at positions relatively "downstream" in the causal structure. If the Developmental Lock model is correct, (1) evolution should have resulted in hierarchically ordered developmental programs, and (2) the most important developmental functions in the hierarchy should be ancient. Wimsatt also suggests that developmental functions be analyzed according to a degree property called "generative entrenchment", which replaces the temporal analysis in the traditional formulation of von Baer's laws. Herein, a substantial body of data on Drosophila ontogeny is analyzed according to generative entrenchment, in order to try the effectiveness of this form of analysis, and also to empirically test these two main predictions of the Developmental Lock model. The novel analytic approach proves to be fruitful, both in generating experimental hypotheses and in ordering existing data. Moreover, data concerning the developmental functions discussed here indicate that the order of the Drosophila developmental program conforms to the predictions of Wimsatt's model with few deviations. Explanations of the anomalies are offered, along with proposals for experiments to test some of those explanations.     

Foam separation of DNA and proteins from solutions

Foam separation on BSA-DNA (bovine serum albumine/deoxyribonucleic acid) and Lysozyme-DNA systems is performed. The separation of the total protein from DNA is evaluated for dissociated chromatin solution. Foam separation for the same systems is done also by a new method of creating a pressure gradient in the Plateau-Gibbs borders in the foam and obtaining a "dry" foam. It is shown that the effectiveness of the foam separation can be improved significantly by the application of the latter method. Some factors (pH, initial concentration of the solution, expansion factor of the foam) influencing the separation of proteins from DNA in the foam and in the residual solution are studied as well.     

Cleavage of cohesin by the CD clan protease separin triggers anaphase in yeast

In eukaryotic cells, replicated DNA strands remain physically connected until their segregation to opposite poles of the cell during anaphase. This "sister chromatid cohesion" is essential for the alignment of chromosomes on the mitotic spindle during metaphase. Cohesion depends on the multisubunit cohesin complex, which possibly forms the physical bridges connecting sisters. Proteolytic cleavage of cohesin's Sccl subunit at the metaphase to anaphase transition is essential for sister chromatid separation and depends on a conserved protein called separin. We show here that separin is a cysteine protease related to caspases that alone can cleave Sccl in vitro. Cleavage of Sccl in metaphase arrested cells is sufficient to trigger the separation of sister chromatids and their segregation to opposite cell poles.     

Determination of constant rates of adsorption of ligand on DNA: analysis of correlation functions

Correlation functions and spectral density of the number of molecules of ligand bound to DNA are calculated theoretically. Kinetics of rates of formation and decomposition of the complex are determined by calculating the dependence of correlation function on concentration of ligand in solution. The analysis of spectral density allows to distinguish "fast" and "slow" adsorption of ligands on macromolecule.     

Graphical approach to model reduction for nonlinear biochemical networks

Model reduction is a central challenge to the development and analysis of multiscale physiology models. Advances in model reduction are needed not only for computational feasibility but also for obtaining conceptual insights from complex systems. Here, we introduce an intuitive graphical approach to model reduction based on phase plane analysis. Timescale separation is identified by the degree of hysteresis observed in phase-loops, which guides a "concentration-clamp" procedure for estimating explicit algebraic relationships between species equilibrating on fast timescales. The primary advantages of this approach over Jacobian-based timescale decomposition are that: 1) it incorporates nonlinear system dynamics, and 2) it can be easily visualized, even directly from experimental data. We tested this graphical model reduction approach using a 25-variable model of cardiac β(1)-adrenergic signaling, obtaining 6- and 4-variable reduced models that retain good predictive capabilities even in response to new perturbations. These 6 signaling species appear to be optimal "kinetic biomarkers" of the overall β(1)-adrenergic pathway. The 6-variable reduced model is well suited for integration into multiscale models of heart function, and more generally, this graphical model reduction approach is readily applicable to a variety of other complex biological systems.     

Functional morphology and evolutionary origin of the three-part pharynx in nematodes

The distinct morphological regions of the typical tripartite pharynx found in the nematode taxon Secernentea have distinctive functions. Besides the basic functions of sucking and pumping food against the pressure in the body cavity, the pharynx of Secernentea such as rhabditids serves two additional functions restricted to two pharyngeal subunits. The corpus traps bacteria behind the stoma and at its posterior end. The newly discovered pharyngeal pocket valve helps to trap particles behind the corpus in the rhabditid Poikilolaimus oxycercus and the cephalobid Acrobeles ciliatus (both Secernentea). The grinder of the terminal bulb serves for chewing trapped bacteria. The separated sites of trapping and chewing are connected by the isthmus that transports bacteria towards the grinder. It is likely that this complex feeding structure originated step by step from a two-part pharynx comprising a propharynx and the terminal bulb as in "Plectidae" (that probably include the closest relatives of the Secernentea within the "Adenophorea"). Analysis of video sequences of feeding rhabditids and plectids provided new data to reconstruct this transformation. Within the "Plectidae" two types of grinders occur. The first type or "parietinus type" has triangular chewing plates that can bulge medially and crush food particles. When they retract, new ingested particles are drawn into the grinder. The second type with more solid chewing plates called "butterfly valves" occurs in Ceratoplectus, Plectus parvus, and Wilsonema and can be homologized with the grinder in Secernentea ("Plectidae" is a paraphyletic taxon). Because butterfly valves cannot be retracted, the evolution of such valves required the evolution of an alternative mechanism to fill the grinder with bacteria. The differentiated closing pattern of the dilated pharynx lumen in Ceratoplectus, Plectus parvus, and Wilsonema can be interpreted as the first step in the development of a functional separation of trapping bacteria and of transporting them towards the grinder, which led to the morphologically discernible units of corpus and isthmus found in the Secernentea.     

The cytochrome bc1 complex in the mitochondrial respiratory chain functions according to the Q cycle hypothesis of Mitchell: the proof using a stochastic approach?

The bc1 complex is a central complex in the mitochondrial respiratory chain. It links the electrons transfer from ubiquinol (or coenzyme Q) to cytochrome c and proton translocation across the inner mitochondrial membrane. It is widely agreed that the "Q-cycle mechanism" proposed by Mitchell correctly describes the bc1 complex working. It is based on an unexpected separation of the two electrons coming from the coenzyme Q bound at the Q0 site of the bc1 complex. Using the stochastic approach of Gillespie and the known spatial structure of bc1 complexes with the kinetic parameters described by Moser and Dutton we demonstrated the natural emergence of the Q-cycle mechanism and the quasi absence of short-circuits in the functional dimer of bc1 complex without the necessity to invoke any additional mechanism. This approach gives a framework which is well adapted to the modelling of all oxido-reduction reactions of the respiratory chain complexes, normal or mutant.     

Capillary electrophoresis of peptides and proteins with plug of Pluronic gel

Electromigration capillary methods are promising techniques in proteomics and they are still under research. We used a partial filling approach, i.e. a combination of gel and non-gel separation mechanisms in a single dimension. We tried using an interesting gel, Pluronic F 127, which can be considered as a surfactant capable of self-association both with isotropic and anisotropic gels. The Pluronic was inserted inside the capillary as a plug at the start of the capillary, and it provided separation at the first time. Separation by this gel was achieved according to molecular weight and/or hydrophobicity. The applicability of this method was demonstrated in the separation of real samples-peptides arising from collagen after CNBr or collagenase cleavage and albumin after trypsin cleavage (peptide mapping). Some peptides and proteins were selectively retained by the Pluronic gel. These interactions with the gel did not depended on their molecular weight alone, but they probably depend on a combination of both principles. It was confirmed that capillary electrophoresis with Pluronic plug can give us another new separation option, complementary to free solution capillary electrophoresis. The CE method presented here, consisting of a partial filling approach with combine gel and non-gel separation mechanisms seemed to be a promising method for the separation of complex mixtures of peptides.     

Glycosylation of hepatitis C virus envelope proteins

Enveloped viruses are surrounded by a membrane derived from the host-cell that contains proteins called "envelope proteins". These proteins play a major role in virus assembly and entry. In most of the enveloped viruses, they are modified by N-linked glycosylation which is supposed to play a role in their stability, antigenicity and biological functions. Glycosylation is also known to play a major role in the biogenesis of proteins by being directly and/or indirectly involved in protein folding. Recent studies on hepatitis C virus (HCV) envelope proteins have revealed a complex interplay between cleavage by signal peptidase, folding and glycosylation. The knowledge that has been accumulated on the early steps of glycosylation of these proteins is presented in this review.     

Functional characterization of individual genomic condensin-binding sites in <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis> using minichromosome mitotic segregation stability assay

Proper chromatin compaction in mitosis (condensation) is required for equal chromosome distribution and the precise inheritance of genetic information. A protein complex called condensin is responsible for mitotic chromosome condensation, chromosome individualization, the timely separation of sister chromatids in mitosis, and proper tension in the mitotic spindle. The mitotic function of condensin depends on the recognition of specific binding sites in the chromosome. The mechanism for binding condensin to individual sites of mitotic chromosomes, as well as the molecular anatomy of these sites, remains to be elucidated. Even less is known about the process that translates condensin binding to individual sites into the segregation of chromosomes in anaphase. In the present work, by using minichromosome assay, we analyze seven individual condensin-binding sites in S. cerevisiae identified in the whole-genome ChIP-on-chip screening. This approach allowed us to estimate the individual contribution of condensin-binding sites to the segregation fidelity of minichromosomes.     

Two-dimensional separation of human plasma proteins using iterative free-flow electrophoresis

Blood plasma is the most complex human-derived proteome, containing other tissue proteomes as subsets. This proteome has only been partially characterized due to the extremely wide dynamic range of the plasma proteins of more than ten orders of magnitude. Thus, the reduction in sample complexity prior to mass spectrometric analysis is particularly important and alternative separation methodologies are required to more effectively mine the lower abundant plasma proteins. Here, we demonstrated a novel separation approach using 2-D free-flow electrophoresis (FFE) separating proteins and peptides in solution according to their pI prior to LC-MS/MS. We used the combination of sequential protein and peptide separation by first separating the plasma proteins into specific FFE fractions. Tryptic digests of the separated proteins were generated and subsequently separated using FFE. The protein separation medium was optimized to segregate albumin into specific fractions containing only few other proteins. An optimization of throughput for the protein separation reduced the separation time of 1 mL of plasma to approximately 3 h providing sufficient material for digestion and the subsequent peptide separation. Our approach revealed low-abundant proteins (e.g., L-selectin at 17 ng/mL and vascular endothelial-cadherin precursor at 30 ng/mL) and several tissue leakage products, thus providing a powerful orthogonal separation step in the proteomics workflow.     

GCTA: a tool for genome-wide complex trait analysis

For most human complex diseases and traits, SNPs identified by genome-wide association studies (GWAS) explain only a small fraction of the heritability. Here we report a user-friendly software tool called genome-wide complex trait analysis (GCTA), which was developed based on a method we recently developed to address the "missing heritability" problem. GCTA estimates the variance explained by all the SNPs on a chromosome or on the whole genome for a complex trait rather than testing the association of any particular SNP to the trait. We introduce GCTA's five main functions: data management, estimation of the genetic relationships from SNPs, mixed linear model analysis of variance explained by the SNPs, estimation of the linkage disequilibrium structure, and GWAS simulation. We focus on the function of estimating the variance explained by all the SNPs on the X chromosome and testing the hypotheses of dosage compensation. The GCTA software is a versatile tool to estimate and partition complex trait variation with large GWAS data sets.     

A new method for the extraction of fetal ECG from the dependent abdominal signals using blind source separation and adaptive noise cancellation techniques

Background

The electrocardiogram (ECG) is a diagnostic tool that records the electrical activity of the heart, and depicts it as a series of graph-like tracings, or waves. Being able to interpret these details allows diagnosis of a wide range of heart problems. Fetal electrocardiogram (FECG) extraction has an important impact in medical diagnostics during the mother pregnancy period. Since the observed FECG signals are often mixed with the maternal ECG (MECG) and the noise induced by the movement of electrodes or by mother motion, the separation process of the ECG signal sources from the observed data becomes quite complicated. One of its complexity is when the ECG sources are dependent, thus, in this paper we introduce a new approach of blind source separation (BSS) in the noisy context for both independent and dependent ECG signal source. This approach consist in denoising the observed ECG signals using a bilateral total variation (BTV) filter; then minimizing the Kullbak-Leibler divergence between copula densities to separate the FECG signal from the MECG one.

Results

We present simulation results illustrating the performance of our proposed method. We will consider many examples of independent/dependent source component signals. The results will be compared with those of the classical method called independent component analysis (ICA) under the same conditions. The accuracy of source estimation is evaluated through a criterion, called again the signal-to-noise-ratio (SNR). The first experiment shows that our proposed method gives accurate estimation of sources in the standard case of independent components, with performance around 27 dB in term of SNR. In the second experiment, we show the capability of the proposed algorithm to successfully separate two noisy mixtures of dependent source components - with classical criterion devoted to the independent case - fails, and that our method is able to deal with the dependent case with good performance.

Conclusions

In this work, we focus specifically on the separation of the ECG signal sources taken from skin two electrodes located on a pregnant woman’s body. The ECG separation is interpreted as a noisy linear BSS problem with instantaneous mixtures. Firstly, a denoising step is required to reduce the noise due to motion artifacts using a BTV filter as a very effective one-pass filter for denoising. Then, we use the Kullbak-Leibler divergence between copula densities to separate the fetal heart rate from the mother one, for both independent and dependent cases.