Additive functions in boolean models of gene regulatory network modules.

Gene-on-gene regulations are key components of every living organism. Dynamical abstract models of genetic regulatory networks help explain the genome's evolvability and robustness. These properties can be attributed to the structural topology of the graph formed by genes, as vertices, and regulatory interactions, as edges. Moreover, the actual gene interaction of each gene is believed to play a key role in the stability of the structure. With advances in biology, some effort was deployed to develop update functions in Boolean models that include recent knowledge. We combine real-life gene interaction networks with novel update functions in a Boolean model. We use two sub-networks of biological organisms, the yeast cell-cycle and the mouse embryonic stem cell, as topological support for our system. On these structures, we substitute the original random update functions by a novel threshold-based dynamic function in which the promoting and repressing effect of each interaction is considered. We use a third real-life regulatory network, along with its inferred Boolean update functions to validate the proposed update function. Results of this validation hint to increased biological plausibility of the threshold-based function. To investigate the dynamical behavior of this new model, we visualized the phase transition between order and chaos into the critical regime using Derrida plots. We complement the qualitative nature of Derrida plots with an alternative measure, the criticality distance, that also allows to discriminate between regimes in a quantitative way. Simulation on both real-life genetic regulatory networks show that there exists a set of parameters that allows the systems to operate in the critical region. This new model includes experimentally derived biological information and recent discoveries, which makes it potentially useful to guide experimental research. The update function confers additional realism to the model, while reducing the complexity and solution space, thus making it easier to investigate.

White matter maturation reshapes structural connectivity in the late developing human brain.

From toddler to late teenager, the macroscopic pattern of axonal projections in the human brain remains largely unchanged while undergoing dramatic functional modifications that lead to network refinement. These functional modifications are mediated by increasing myelination and changes in axonal diameter and synaptic density, as well as changes in neurochemical mediators. Here we explore the contribution of white matter maturation to the development of connectivity between ages 2 and 18 y using high b-value diffusion MRI tractography and connectivity analysis. We measured changes in connection efficacy as the inverse of the average diffusivity along a fiber tract. We observed significant refinement in specific metrics of network topology, including a significant increase in node strength and efficiency along with a decrease in clustering. Major structural modules and hubs were in place by 2 y of age, and they continued to strengthen their profile during subsequent development. Recording resting-state functional MRI from a subset of subjects, we confirmed a positive correlation between structural and functional connectivity, and in addition observed that this relationship strengthened with age. Continuously increasing integration and decreasing segregation of structural connectivity with age suggests that network refinement mediated by white matter maturation promotes increased global efficiency. In addition, the strengthening of the correlation between structural and functional connectivity with age suggests that white matter connectivity in combination with other factors, such as differential modulation of axonal diameter and myelin thickness, that are partially captured by inverse average diffusivity, play an increasingly important role in creating brain-wide coherence and synchrony.

SwissSidechain: a molecular and structural database of non-natural sidechains.

Amino acids form the building blocks of all proteins. Naturally occurring amino acids are restricted to a few tens of sidechains, even when considering post-translational modifications and rare amino acids such as selenocysteine and pyrrolysine. However, the potential chemical diversity of amino acid sidechains is nearly infinite. Exploiting this diversity by using non-natural sidechains to expand the building blocks of proteins and peptides has recently found widespread applications in biochemistry, protein engineering and drug design. Despite these applications, there is currently no unified online bioinformatics resource for non-natural sidechains. With the SwissSidechain database (http://www.swisssidechain.ch), we offer a central and curated platform about non-natural sidechains for researchers in biochemistry, medicinal chemistry, protein engineering and molecular modeling. SwissSidechain provides biophysical, structural and molecular data for hundreds of commercially available non-natural amino acid sidechains, both in l- and d-configurations. The database can be easily browsed by sidechain names, families or physico-chemical properties. We also provide plugins to seamlessly insert non-natural sidechains into peptides and proteins using molecular visualization software, as well as topologies and parameters compatible with molecular mechanics software.

Enzymic, phylogenetic, and structural characterization of the unusual papain-like protease domain of Plasmodium falciparum SERA5.

Serine repeat antigen 5 (SERA5) is an abundant antigen of the human malaria parasite Plasmodium falciparum and is the most strongly expressed member of the nine-gene SERA family. It appears to be essential for the maintenance of the erythrocytic cycle, unlike a number of other members of this family, and has been implicated in parasite egress and/or erythrocyte invasion. All SERA proteins possess a central domain that has homology to papain except in the case of SERA5 (and some other SERAs), where the active site cysteine has been replaced with a serine. To investigate if this domain retains catalytic activity, we expressed, purified, and refolded a recombinant form of the SERA5 enzyme domain. This protein possessed chymotrypsin-like proteolytic activity as it processed substrates downstream of aromatic residues, and its activity was reversed by the serine protease inhibitor 3,4-diisocoumarin. Although all Plasmodium SERA enzyme domain sequences share considerable homology, phylogenetic studies revealed two distinct clusters across the genus, separated according to whether they possess an active site serine or cysteine. All Plasmodia appear to have at least one member of each group. Consistent with separate biological roles for members of these two clusters, molecular modeling studies revealed that SERA5 and SERA6 enzyme domains have dramatically different surface properties, although both have a characteristic papain-like fold, catalytic cleft, and an appropriately positioned catalytic triad. This study provides impetus for the examination of SERA5 as a target for antimalarial drug design.

Structural and functional studies of nonstructural protein 2 of the hepatitis C virus reveal its key role as organizer of virion assembly.

Non-structural protein 2 (NS2) plays an important role in hepatitis C virus (HCV) assembly, but neither the exact contribution of this protein to the assembly process nor its complete structure are known. In this study we used a combination of genetic, biochemical and structural methods to decipher the role of NS2 in infectious virus particle formation. A large panel of NS2 mutations targeting the N-terminal membrane binding region was generated. They were selected based on a membrane topology model that we established by determining the NMR structures of N-terminal NS2 transmembrane segments. Mutants affected in virion assembly, but not RNA replication, were selected for pseudoreversion in cell culture. Rescue mutations restoring virus assembly to various degrees emerged in E2, p7, NS3 and NS2 itself arguing for an interaction between these proteins. To confirm this assumption we developed a fully functional JFH1 genome expressing an N-terminally tagged NS2 demonstrating efficient pull-down of NS2 with p7, E2 and NS3 and, to a lower extent, NS5A. Several of the mutations blocking virus assembly disrupted some of these interactions that were restored to various degrees by those pseudoreversions that also restored assembly. Immunofluorescence analyses revealed a time-dependent NS2 colocalization with E2 at sites close to lipid droplets (LDs) together with NS3 and NS5A. Importantly, NS2 of a mutant defective in assembly abrogates NS2 colocalization around LDs with E2 and NS3, which is restored by a pseudoreversion in p7, whereas NS5A is recruited to LDs in an NS2-independent manner. In conclusion, our results suggest that NS2 orchestrates HCV particle formation by participation in multiple protein-protein interactions required for their recruitment to assembly sites in close proximity of LDs.

Structural determinants for membrane association and dynamic organization of the hepatitis C virus NS3-4A complex.

Hepatitis C virus (HCV) NS3-4A is a membrane-associated multifunctional protein harboring serine protease and RNA helicase activities. It is an essential component of the HCV replication complex and a prime target for antiviral intervention. Here, we show that membrane association and structural organization of HCV NS3-4A are ensured in a cooperative manner by two membrane-binding determinants. We demonstrate that the N-terminal 21 amino acids of NS4A form a transmembrane alpha-helix that may be involved in intramembrane protein-protein interactions important for the assembly of a functional replication complex. In addition, we demonstrate that amphipathic helix alpha(0), formed by NS3 residues 12-23, serves as a second essential determinant for membrane association of NS3-4A, allowing proper positioning of the serine protease active site on the membrane. These results allowed us to propose a dynamic model for the membrane association, processing, and structural organization of NS3-4A on the membrane. This model has implications for the functional architecture of the HCV replication complex, proteolytic targeting of host factors, and drug design.

X-ray structures of Sap1 and Sap5: structural comparison of the secreted aspartic proteinases from Candida albicans.

Proteolytic activity is an important virulence factor for Candida albicans (C. albicans). It is attributed to the family of the secreted aspartic proteinases (Saps) from C. albicans with a minimum of 10 members. Saps show controlled expression and regulation for the individual stages of the infection process. Distinct isoenzymes can be responsible for adherence and tissue damage of local infections, while others cause systemic diseases. Earlier, only the structures of Sap2 and Sap3 were known. In our research, we have now succeeded in solving the X-ray crystal structures of the apoenzyme of Sap1 and Sap5 in complex with pepstatin A at 2.05 and 2.5 A resolution, respectively. With the structure of Sap1, we have completed the set of structures of isoenzyme subgroup Sap1-3. Of subgroup Sap4-6, the structure of the enzyme Sap5 is the first structure that has been described up to now. This facilitates comparison of structural details as well as inhibitor binding modes among the different subgroup members. Structural analysis reveals a highly conserved overall secondary structure of Sap1-3 and Sap5. However, Sap5 clearly differs from Sap1-3 by its electrostatic overall charge as well as through structural conformation of its entrance to the active site cleft. Design of inhibitors specific for Sap5 should concentrate on the S4 and S3 pockets, which significantly differ from Sap1-3 in size and electrostatic charge. Both Sap1 and Sap5 seem to play a major part in superficial Candida infections. Determination of the isoenzymes' structures can contribute to the development of new Sap-specific inhibitors for the treatment of superficial infections with a structure-based drug design program.

Malaria vaccine candidate: design of a multivalent subunit α-helical coiled coil poly-epitope.

A new strategy for the rapid identification of new malaria antigens based on protein structural motifs was previously described. We identified and evaluated the malaria vaccine potential of fragments of several malaria antigens containing α-helical coiled coil protein motifs. By taking advantage of the relatively short size of these structural fragments, we constructed different poly-epitopes in which 3 or 4 of these segments were joined together via a non-immunogenic linker. Only peptides that are targets of human antibodies with anti-parasite in vitro biological activities were incorporated. One of the constructs, P181, was well recognized by sera and peripheral blood mononuclear cells (PBMC) of adults living in malaria-endemic areas. Affinity purified antigen-specific human antibodies and sera from P181-immunized mice recognised native proteins on malaria-infected erythrocytes in both immunofluorescence and western blot assays. In addition, specific antibodies inhibited parasite development in an antibody dependent cellular inhibition (ADCI) assay. Naturally induced antigen-specific human antibodies were at high titers and associated with clinical protection from malaria in longitudinal follow-up studies in Senegal.

Expanding molecular modeling and design tools to non-natural sidechains.

Protein-protein interactions encode the wiring diagram of cellular signaling pathways and their deregulations underlie a variety of diseases, such as cancer. Inhibiting protein-protein interactions with peptide derivatives is a promising way to develop new biological and therapeutic tools. Here, we develop a general framework to computationally handle hundreds of non-natural amino acid sidechains and predict the effect of inserting them into peptides or proteins. We first generate all structural files (pdb and mol2), as well as parameters and topologies for standard molecular mechanics software (CHARMM and Gromacs). Accurate predictions of rotamer probabilities are provided using a novel combined knowledge and physics based strategy. Non-natural sidechains are useful to increase peptide ligand binding affinity. Our results obtained on non-natural mutants of a BCL9 peptide targeting beta-catenin show very good correlation between predicted and experimental binding free-energies, indicating that such predictions can be used to design new inhibitors. Data generated in this work, as well as PyMOL and UCSF Chimera plug-ins for user-friendly visualization of non-natural sidechains, are all available at http://www.swisssidechain.ch. Our results enable researchers to rapidly and efficiently work with hundreds of non-natural sidechains.

Extreme prematurity and intra uterine growth restriction effects in brain network topology at school age

Higher risk for long-term behavioral and emotional sequelae, with attentional problems (with or without hyperactivity) is now becoming one of the hallmarks of extreme premature (EP) birth and birth after pregancy conditions leading to poor intra uterine growth restriction (IUGR) [1,2]. However, little is know so far about the neurostructural basis of these complexe brain functional abnormalities that seem to have their origins in early critical periods of brain development. The development of cortical axonal pathways happens in a series of sequential events. The preterm phase (24-36 post conecptional weeks PCW) is known for being crucial for growth of the thalamocortical fiber bundles as well as for the development of long projectional, commisural and projectional fibers [3]. Is it logical to expect, thus, that being exposed to altered intrauterine environment (altered nutrition) or to extrauterine environment earlier that expected, lead to alterations in the structural organization and, consequently, alter the underlying white matter (WM) structure. Understanding rate and variability of normal brain development, and detect differences from typical development may offer insight into the neurodevelopmental anomalies that can be imaged at later stages. Due to its unique ability to non-invasively visualize and quantify in vivo white matter tracts in the brain, in this study we used diffusion MRI (dMRI) tractography to derive brain graphs [4,5,6]. This relatively simple way of modeling the brain enable us to use graph theory to study topological properties of brain graphs in order to study the effects of EP and IUGR on childrens brain connectivity at age 6 years old.

Structural organization of alpha-subunit from purified and microsomal toad kidney (Na+ + K+)-ATPase as assessed by controlled trypsinolysis

The membrane organization of the alpha-subunit of purified (Na+ + K+)-ATPase ((Na+ + K+)-dependent adenosine triphosphate phosphorylase, EC 3.6.1.3) and of the microsomal enzyme of the kidney of the toad Bufo marinus was compared by using controlled trypsinolysis. With both enzyme preparations, digestions performed in the presence of Na+ yielded a 73 kDa fragment and in the presence of K+ a 56 kDa, a 40 kDa and small amounts of a 83 kDa fragment from the 96 kDa alpha-subunit. In contrast to mammalian preparations (Jørgensen, P.L. (1975) Biochim. Biophys. Acta 401, 399-415), trypsinolysis of the purified amphibian enzyme led to a biphasic loss of (Na+ + K+)-ATPase activity in the presence of both Na+ and K+. These data could be correlated with an early rapid cleavage of 3 kDa from the alpha-subunit in both ionic conditions and a slower degradation of the remaining 93 kDa polypeptide. On the other hand, in the microsomal enzyme, a 3 kDa shift of the alpha-subunit could only be produced in the presence of Na+. Our data indicate that (1) purification of the amphibian enzyme with detergent does not influence the overall topology of the alpha-subunit but produces a distinct structural alteration of its N-terminus and (2) the amphibian kidney enzyme responds to cations with similar conformational transitions as the mammalian kidney enzyme. In addition, anti alpha-serum used on digested enzyme samples revealed on immunoblots that the 40 kDa fragment was better recognized than the 56 kDa fragment. It is concluded that the NH2-terminal of the alpha-subunit contains more antigenic sites than the COOH-terminal domain in agreement with the results of Farley et al. (Farley, R.A., Ochoa, G.T. and Kudrow, A. (1986) Am. J. Physiol. 250, C896-C906).

Structural assessment of single amino acid mutations: application to TP53 function.

Single amino acid substitution is the type of protein alteration most related to human diseases. Current studies seek primarily to distinguish neutral mutations from harmful ones. Very few methods offer an explanation of the final prediction result in terms of the probable structural or functional effect on the protein. In this study, we describe the use of three novel parameters to identify experimentally-verified critical residues of the TP53 protein (p53). The first two parameters make use of a surface clustering method to calculate the protein surface area of highly conserved regions or regions with high nonlocal atomic interaction energy (ANOLEA) score. These parameters help identify important functional regions on the surface of a protein. The last parameter involves the use of a new method for pseudobinding free-energy estimation to specifically probe the importance of residue side-chains to the stability of protein fold. A decision tree was designed to optimally combine these three parameters. The result was compared to the functional data stored in the International Agency for Research on Cancer (IARC) TP53 mutation database. The final prediction achieved a prediction accuracy of 70% and a Matthews correlation coefficient of 0.45. It also showed a high specificity of 91.8%. Mutations in the 85 correctly identified important residues represented 81.7% of the total mutations recorded in the database. In addition, the method was able to correctly assign a probable functional or structural role to the residues. Such information could be critical for the interpretation and prediction of the effect of missense mutations, as it not only provided the fundamental explanation of the observed effect, but also helped design the most appropriate laboratory experiment to verify the prediction results.

Structural abnormalities in the thalamus of migraineurs with aura: a multiparametric study at 3 T.

BACKGROUND AND OBJECTIVES: The thalamus exerts a pivotal role in pain processing and cortical excitability control, and migraine is characterized by repeated pain attacks and abnormal cortical habituation to excitatory stimuli. This work aimed at studying the microstructure of the thalamus in migraine patients using an innovative multiparametric approach at high-field magnetic resonance imaging (MRI). DESIGN: We examined 37 migraineurs (22 without aura, MWoA, and 15 with aura, MWA) as well as 20 healthy controls (HC) in a 3-T MRI equipped with a 32-channel coil. We acquired whole-brain T1 relaxation maps and computed magnetization transfer ratio (MTR), generalized fractional anisotropy, and T2* maps to probe microstructural and connectivity integrity and to assess iron deposition. We also correlated the obtained parametric values with the average monthly frequency of migraine attacks and disease duration. RESULTS: T1 relaxation time was significantly shorter in the thalamus of MWA patients compared with MWoA (P < 0.001) and HC (P ≤ 0.01); in addition, MTR was higher and T2* relaxation time was shorter in MWA than in MWoA patients (P < 0.05, respectively). These data reveal broad microstructural alterations in the thalamus of MWA patients compared with MWoA and HC, suggesting increased iron deposition and myelin content/cellularity. However, MWA and MWoA patients did not show any differences in the thalamic nucleus involved in pain processing in migraine. CONCLUSIONS: There are broad microstructural alterations in the thalamus of MWA patients that may underlie abnormal cortical excitability control leading to cortical spreading depression and visual aura.

Survey of prenatal screening policies in Europe for structural malformations and chromosome anomalies, and their impact on detection and termination rates for neural tube defects and Down's syndrome.

OBJECTIVE: To 'map' the current (2004) state of prenatal screening in Europe. DESIGN: (i) Survey of country policies and (ii) analysis of data from EUROCAT (European Surveillance of Congenital Anomalies) population-based congenital anomaly registers. SETTING: Europe. POPULATION: Survey of prenatal screening policies in 18 countries and 1.13 million births in 12 countries in 2002-04. METHODS: (i) Questionnaire on national screening policies and termination of pregnancy for fetal anomaly (TOPFA) laws in 2004. (ii) Analysis of data on prenatal detection and termination for Down's syndrome and neural tube defects (NTDs) using the EUROCAT database. MAIN OUTCOME MEASURES: Existence of national prenatal screening policies, legal gestation limit for TOPFA, prenatal detection and termination rates for Down's syndrome and NTD. RESULTS: Ten of the 18 countries had a national country-wide policy for Down's syndrome screening and 14/18 for structural anomaly scanning. Sixty-eight percent of Down's syndrome cases (range 0-95%) were detected prenatally, of which 88% resulted in termination of pregnancy. Eighty-eight percent (range 25-94%) of cases of NTD were prenatally detected, of which 88% resulted in termination. Countries with a first-trimester screening policy had the highest proportion of prenatally diagnosed Down's syndrome cases. Countries with no official national Down's syndrome screening or structural anomaly scan policy had the lowest proportion of prenatally diagnosed Down's syndrome and NTD cases. Six of the 18 countries had a legal gestational age limit for TOPFA, and in two countries, termination of pregnancy was illegal at any gestation. CONCLUSIONS: There are large differences in screening policies between countries in Europe. These, as well as organisational and cultural factors, are associated with wide country variation in prenatal detection rates for Down's syndrome and NTD.

Structural connectomics in brain diseases.

Imaging the connectome in vivo has become feasible through the integration of several rapidly developing fields of science and engineering, namely magnetic resonance imaging and in particular diffusion MRI on one side, image processing and network theory on the other side. This framework brings in vivo brain imaging closer to the real topology of the brain, contributing to narrow the existing gap between our understanding of brain structural organization on one side and of human behavior and cognition on the other side. Given the seminal technical progresses achieved in the last few years, it may be ready to tackle even greater challenges, namely exploring disease mechanisms. In this review we analyze the current situation from the technical and biological perspectives. First, we critically review the technical solutions proposed in the literature to perform clinical studies. We analyze for each step (i.e. MRI acquisition, network building and network statistical analysis) the advantages and potential limitations. In the second part we review the current literature available on a selected subset of diseases, namely, dementia, schizophrenia, multiple sclerosis and others, and try to extract for each disease the common findings and main differences between reports.