Disentangling the Complexity of HGF Signaling by Combining Qualitative and Quantitative Modeling.

Signaling pathways are characterized by crosstalk, feedback and feedforward mechanisms giving rise to highly complex and cell-context specific signaling networks. Dissecting the underlying relations is crucial to predict the impact of targeted perturbations. However, a major challenge in identifying cell-context specific signaling networks is the enormous number of potentially possible interactions. Here, we report a novel hybrid mathematical modeling strategy to systematically unravel hepatocyte growth factor (HGF) stimulated phosphoinositide-3-kinase (PI3K) and mitogen activated protein kinase (MAPK) signaling, which critically contribute to liver regeneration. By combining time-resolved quantitative experimental data generated in primary mouse hepatocytes with interaction graph and ordinary differential equation modeling, we identify and experimentally validate a network structure that represents the experimental data best and indicates specific crosstalk mechanisms. Whereas the identified network is robust against single perturbations, combinatorial inhibition strategies are predicted that result in strong reduction of Akt and ERK activation. Thus, by capitalizing on the advantages of the two modeling approaches, we reduce the high combinatorial complexity and identify cell-context specific signaling networks.

Whole exome sequencing for a patient with Rubinstein-Taybi syndrome reveals de novo variants besides an overt CREBBP mutation.

Rubinstein-Taybi syndrome (RSTS) is a rare condition with a prevalence of 1 in 125,000-720,000 births and characterized by clinical features that include facial, dental, and limb dysmorphology and growth retardation. Most cases of RSTS occur sporadically and are caused by de novo mutations. Cytogenetic or molecular abnormalities are detected in only 55% of RSTS cases. Previous genetic studies have yielded inconsistent results due to the variety of methods used for genetic analysis. The purpose of this study was to use whole exome sequencing (WES) to evaluate the genetic causes of RSTS in a young girl presenting with an Autism phenotype. We used the Autism diagnostic observation schedule (ADOS) and Autism diagnostic interview revised (ADI-R) to confirm her diagnosis of Autism. In addition, various questionnaires were used to evaluate other psychiatric features. We used WES to analyze the DNA sequences of the patient and her parents and to search for de novo variants. The patient showed all the typical features of Autism, WES revealed a de novo frameshift mutation in CREBBP and de novo sequence variants in TNC and IGFALS genes. Mutations in the CREBBP gene have been extensively reported in RSTS patients, while potential missense mutations in TNC and IGFALS genes have not previously been associated with RSTS. The TNC and IGFALS genes are involved in central nervous system development and growth. It is possible for patients with RSTS to have additional de novo variants that could account for previously unexplained phenotypes.

Transmission dynamics of oral polio vaccine viruses and vaccine-derived polioviruses on networks.

One drawback of oral polio vaccine (OPV) is the potential reversion to more transmissible, virulent circulating vaccine-derived polioviruses (cVDPVs), which may cause outbreaks of paralytic poliomyelitis. Previous modeling studies of the transmission of cVDPVs assume an unrealistic homogeneous mixing of the population and/or ignore that OPV viruses and cVDPVs compete for susceptibles, which we show is a key to understanding the dynamics of the transmission of cVDPVs. We examined the transmission of OPV viruses and cVDPVs on heterogeneous, dynamic contact networks using differential equation-based and individual-based models. Despite the lower transmissibility, OPV viruses may outcompete more transmissible cVDPVs in the short run by spreading extensively before cVDPVs emerge. If viruses become endemic, however, cVDPVs eventually dominate and force OPV viruses to extinction. This study improves our understanding of the emergence of cVDPVs and helps develop more detailed models to plan a policy to control paralytic polio associated with the continued use of OPV in many countries.

An engineered viral protease exhibiting substrate specificity for a polyglutamine stretch prevents polyglutamine-induced neuronal cell death.

BACKGROUND: Polyglutamine (polyQ)-induced protein aggregation is the hallmark of a group of neurodegenerative diseases, including Huntington's disease. We hypothesized that a protease that could cleave polyQ stretches would intervene in the initial events leading to pathogenesis in these diseases. To prove this concept, we aimed to generate a protease possessing substrate specificity for polyQ stretches. METHODOLOGY/PRINCIPAL FINDINGS: Hepatitis A virus (HAV) 3C protease (3CP) was subjected to engineering using a yeast-based method known as the Genetic Assay for Site-specific Proteolysis (GASP). Analysis of the substrate specificity revealed that 3CP can cleave substrates containing glutamine at positions P5, P4, P3, P1, P2', or P3', but not substrates containing glutamine at the P2 or P1' positions. To accommodate glutamine at P2 and P1', key residues comprising the active sites of the S2 or S1' pockets were separately randomized and screened. The resulting sets of variants were combined by shuffling and further subjected to two rounds of randomization and screening using a substrate containing glutamines from positions P5 through P3'. One of the selected variants (Var26) reduced the expression level and aggregation of a huntingtin exon1-GFP fusion protein containing a pathogenic polyQ stretch (HttEx1(97Q)-GFP) in the neuroblastoma cell line SH-SY5Y. Var26 also prevented cell death and caspase 3 activation induced by HttEx1(97Q)-GFP. These protective effects of Var26 were proteolytic activity-dependent. CONCLUSIONS/SIGNIFICANCE: These data provide a proof-of-concept that proteolytic cleavage of polyQ stretches could be an effective modality for the treatment of polyQ diseases.

Inhibition of melanogenesis by 5,7-dihydroxyflavone (chrysin) via blocking adenylyl cyclase activity.

Due to its multiple biological activities, 5,7-dihydroxyflavone (chrysin) in propolis has gained attention as potentially useful therapeutics for various diseases. However, the efficacy of chrysin for the use of dermatological health has not been fully explored. To clarify the action mechanism of the skin protecting property of chrysin, we firstly investigated the molecular docking property of chrysin on the mammalian adenylyl cyclase, which is the key enzyme of cAMP-induced melanogenesis. We also examined the involvement of chrysin in alpha-MSH and forskolin-induced cAMP signaling within a cell based assay. In addition, we inquired into the inhibitory effect of chrysin on melanogenesis and found that the pretreatment with chrysin inhibited the forskolin-induced melanin contents significantly without annihilating the cell viability. These results strongly suggest that chrysin directly inhibits the activity of adenylyl cyclase, downregulates forskolin-induced cAMP-production pathway, consequently inhibiting melanogenesis. Thus, chrysin may also be used as an effective inhibitor of hyperpigmentation.

Engineering of protease variants exhibiting altered substrate specificity.

By using an improved genetic screening system, variants of the HAV 3CP protease which exhibit altered P2 specificity were obtained. We randomly mutated the His145, Lys146, Lys147, and Leu155 residues that constitute the S2 pocket of 3CP and then isolated variants that preferred substrates with Gln over the original Thr at the P2 position using a yeast-based screening method. One of the isolated variants cleaved the Gln-containing peptide substrate more efficiently in vitro, proving the efficiency of our method in isolating engineered proteases with desired substrate selectivity.

Identification of mouse heart transcriptomic network sensitive to various heart diseases.

Exploring biological systems from highly complex datasets is an important task for systems biology. The present study examined co-expression dynamics of mouse heart transcriptome by spectral graph clustering (SGC) to identify a heart transcriptomic network. SGC of microarray data produced 17 classified biological conditions (called condition spectrum, CS) and co-expression patterns by generating bi-clusters. The results showed dynamic co-expression patterns with a modular structure enriched in heart-related CS (CS-1 and -13) containing abundant heart-related microarray data. Consequently, a mouse heart transcriptomic network was constructed by clique analysis from the gene clusters exclusively present in the heart-related CS; 31 cliques were used for constructing the network. The participating genes in the network were closely associated with important cardiac functions (e. g., development, lipid and glycogen metabolisms). Online Mendelian Inheritance in Man (OMIM) database indicates that mutations of the genes in the network induced serious heart diseases. Many of the tested genes in the network showed significantly altered gene expression in an animal model of hypertrophy. The results suggest that the present approach is critical for constructing a heart-related transcriptomic network and for deducing important genes involved in the pathogenesis of various heart diseases.

An overview of cardiac systems biology.

The cardiac system has been a major target for intensive studies in the multi-scale modeling field for many years. Reproduction of the action potential and the ionic currents of single cardiomyocytes, as well as the construction of a whole organ model is well established. Still, there are major hurdles to overcome in creating a realistic and predictive functional cardiac model due to the lack of a profound understanding of the complex molecular interactions and their outcomes controlling both normal and pathological cardiophysiology. The recent advent of systems biology offers the conceptual and practical frameworks to tackle such biological complexities. This review provides an overview of major themes in the developing field of cardiac systems biology, summarizing some of the high-throughput experiments and strategies used to integrate the datasets, and various types of computational approaches used for developing useful quantitative models capable of predicting complex biological behavior.

Modulation of the conductance-voltage relationship of the BK Ca channel by mutations at the putative flexible interface between two RCK domains.

Calcium-dependent gating of the large-conductance Ca(2+)-activated K(+) (BK(Ca)) channel is conferred by the large cytosolic carboxyl terminus containing two domains of the regulator of K(+) conductance (RCK) and the high-affinity Ca(2+)-binding site (the Ca(2+)-bowl). In our previous study, we located the putative second RCK domain (RCK2) and demonstrated that it interacts directly with RCK1 via a hydrophobic "assembly interface". In this study, we tested the structural model of the other interface, the "flexible interface", by strategically positioning charge pairs across the putative interface. Several charge mutations on RCK2 affected the voltage-dependent activation of the channel. In particular, the Gly-to-Asp substitution at position 803 profoundly affected channel activation by stabilizing the open conformation of the channel with minimal effects on its Ca(2+) affinity and voltage sensitivity. Various mutations at Gly-803 shifted the channel's conductance-voltage curve either left or right over a 145-mV range. Since this residue is predicted to be in the first loop of RCK2 these results strongly suggest that this loop plays a critical role in determining the intrinsic equilibrium constant for channel opening, and they support the hypothesis that this loop is part of an interface that mediates conformational coupling between RCK1 and RCK2.

Overexpression and purification of the RyR1 pore-forming region.

Ryanodine receptor 1 (RyR1) is a large homotetrameric calcium channel that plays a pivotal role in skeletal muscle contraction. Sequence comparison and mutagenesis studies indicate that the pore architecture of RyR1, including the last two transmembrane helices and the luminal loop linking them, is similar to that of the bacterial KcsA K(+) channel. Here, we describe the overexpression and purification of the C-terminal polyhistidine-tagged RyR1 pore-forming region. The nonionic detergent lauryldimethylamine oxide (LDAO) was selected for solubilization of the protein based on its ability to extract the protein from the membrane and to maintain it in a monodisperse state. The protein was then purified using nickel-affinity chromatography and gel filtration. Gel filtration analysis confirmed that the RyR1 fragment containing the pore-forming region (amino acids 4829-5037) is sufficient to form a tetramer.

Development of potent inhibitors of the coxsackievirus 3C protease.

Coxsackievirus B3 (CVB3) 3C protease (3CP) plays essential roles in the viral replication cycle, and therefore, provides an attractive therapeutic target for treatment of human diseases caused by CVB3 infection. CVB3 3CP and human rhinovirus (HRV) 3CP have a high degree of amino acid sequence similarity. Comparative modeling of these two 3CPs revealed one prominent distinction; an Asn residue delineating the S2' pocket in HRV 3CP is replaced by a Tyr residue in CVB3 3CP. AG7088, a potent inhibitor of HRV 3CP, was modified by substitution of the ethyl group at the P2' position with various hydrophobic aromatic rings that are predicted to interact preferentially with the Tyr residue in the S2' pocket of CVB3 3CP. The resulting derivatives showed dramatically increased inhibitory activities against CVB3 3CP. In addition, one of the derivatives effectively inhibited the CVB3 proliferation in vitro.

Functional importance of polymerization and localization of calsequestrin in C. elegans.

Dual roles of calsequestrin (CSQ-1) being the Ca2+ donor and Ca2+ acceptor make it an excellent Ca2+-buffering protein within the sarcoplasmic reticulum (SR). We have isolated and characterized a calsequestrin (csq-1)-null mutant in Caenorhabditis elegans. To our surprise, this mutant csq-1(jh109) showed no gross defects in muscle development or function but, however, is highly sensitive to perturbation of Ca2+ homeostasis. By taking advantage of the viable null mutant, we investigated the domains of CSQ-1 that are important for polymerization and cellular localization, and required for its correct buffering functions. In transgenic animals rescued with various CSQ-1 constructs, the in vivo patterns of polymerization and localization of several mutated calsequestrins were observed to correlate with the structure-function relationship. Our results suggest that polymerization of CSQ-1 is essential but not sufficient for correct cellular localization and function of CSQ-1. In addition, direct interaction between CSQ-1 and the ryanodine receptor (RyR) was found for the first time, suggesting that the cellular localization of CSQ-1 in C. elegans is indeed modulated by RyR through a physical interaction.

Hydrophobic interface between two regulators of K+ conductance domains critical for calcium-dependent activation of large conductance Ca2+-activated K+ channels.

It has been suggested that the large conductance Ca(2)+-activated K(+) channel contains one or more domains known as regulators of K(+) conductance (RCK) in its cytosolic C terminus. Here, we show that the second RCK domain (RCK2) is functionally important and that it forms a heterodimer with RCK1 via a hydrophobic interface. Mutant channels lacking RCK2 are nonfunctional despite their tetramerization and surface expression. The hydrophobic residues that are expected to form an interface between RCK1 and RCK2, based on the crystal structure of the bacterial MthK channel, are well conserved, and the interactions of these residues were confirmed by mutant cycle analysis. The hydrophobic interaction appears to be critical for the Ca(2+)-dependent gating of the large conductance Ca(2+)-activated K(+) channel.

HCNet: a database of heart and calcium functional network.

SUMMARY: The Heart and Calcium functional Network (HCNet) database is a collection of functional gene modules calculated from the microarray data compendium available from the GEO database. It is a specialized database designed to assist experimentalists for cardiac calcium signaling research by providing the pre-calculated gene clusters and their potential correlation network in heart. In the current release of HCNet, 57 functional modules from 786 target genes obtained by a bi-clustering analysis of 381 microarray datasets are available. Detailed information of the clusters such as expression profiles, network diagrams is provided in two categories, heart-specific genes and heart-specific genes along with calcium toolkit genes. Overrepresented gene ontological categories and transcription factors in each cluster are also provided to infer the biological implications of the detected functional modules. AVAILABILITY: HCNet is available at http://sbrg2.gist.ac.kr/hcnet.

Crystal structure and functional studies reveal that PAS factor from Vibrio vulnificus is a novel member of the saposin-fold family.

PAS factor is a novel putative bacterial secretion factor thought to induce secretion of periplasmic proteins. We solved the crystal structure of PAS factor from Vibrio vulnificus at 1.8A resolution and found it to be comprised of five alpha helices that form an antiparallel bundle with an up-and-down topology, and to adopt the saposin-fold characteristic of a family of proteins that bind to membranes and lipids. PAS factor lacks the disulfide bridge characteristic of mammalian saposin-fold proteins; in fact, it shows no sequence homology with mammalian proteins. Nevertheless, the molecular architectures are similar, and the shared propensity for membrane interaction suggests strongly that PAS factor is another member of the saposin-fold family. Analysis of the CD spectra showed that PAS factor binds to membranes directly, while measurement of calcein dye leakage showed that PAS factor interacts strongly with liposomes composed of anionic phospholipids, making them leaky, but binds very weakly with liposomes composed of zwitterionic phospholipids. Moreover, by analyzing tryptophan fluorescence emission from four single-tryptophan mutants (V10W, T22W, F35W, and L70W), we identified the putative phospholipid-binding site of PAS factor. The resultant membrane destabilization likely mediates secretion of periplasmic proteins required for the in vivo survival and pathogenesis of V.vulnificus.

Engineering of Kex2 variants exhibiting altered substrate specificity.

Engineering of secreted protease variants exhibiting altered substrate specificity is a challenging task because effective screening methods for the desired property are not available yet. In this study, we sought to obtain variants of Kex2, a yeast Golgi protease, which exhibit altered P2 specificity. We first randomly mutated three Asp residues (D176, D210, and D211) that constitute the S2 pocket of Kex2 and then isolated from the resulting library Kex2 variants that preferred substrates with Met (poorly preferred by wild type Kex2) at the P2 position using a yeast-based screening method. The Kex2 variants isolated from this initial screening were further tested against various substrate sequences. Four out of the 16 isolated Kex2 variants showed greater preference for Met than for Lys (preferred by the wild-type Kex2) at the P2 position. We therefore suggest that our method might serve as an efficient tool for engineering and directing the evolution of secreted proteases.

Correction of X-ray intensities from an HslV-HslU co-crystal containing lattice-translocation defects.

Because of lattice-translocation defects, two identical but translated lattices can coexist as a single coherent mosaic block in a crystal. The observed structure in such cases is a weighted sum of two identical but translated structures, one from each lattice; the observed structure factors are a weighted vector sum of the structure factors with identical unit amplitudes but shifted phases. The correction of X-ray intensities from a single crystal containing these defects of the hybrid HslV-HslU complex, which consists of Escherichia coli HslU and Bacillus subtilis HslV (also known as CodW), is reported. When intensities are not corrected, a biologically irrelevant complex (with CodW from one lattice and HslU from another) is implied to exist. Only upon correction does a biologically functional CodW-HslU complex structure emerge.