Perturbation Biology: Inferring Signaling Networks in Cellular Systems

We present a powerful experimental-computational technology for inferring network models that predict the response of cells to perturbations, and that may be useful in the design of combinatorial therapy against cancer. The experiments are systematic series of perturbations of cancer cell lines by targeted drugs, singly or in combination. The response to perturbation is quantified in terms of relative changes in the measured levels of proteins, phospho-proteins and cellular phenotypes such as viability. Computational network models are derived de novo, i.e., without prior knowledge of signaling pathways, and are based on simple non-linear differential equations. The prohibitively large solution space of all possible network models is explored efficiently using a probabilistic algorithm, Belief Propagation (BP), which is three orders of magnitude faster than standard Monte Carlo methods. Explicit executable models are derived for a set of perturbation experiments in SKMEL-133 melanoma cell lines, which are resistant to the therapeutically important inhibitor of RAF kinase. The resulting network models reproduce and extend known pathway biology. They empower potential discoveries of new molecular interactions and predict efficacious novel drug perturbations, such as the inhibition of PLK1, which is verified experimentally. This technology is suitable for application to larger systems in diverse areas of molecular biology.     

Prediction of Drug-Target Interactions for Drug Repositioning Only Based on Genomic Expression Similarity

Small drug molecules usually bind to multiple protein targets or even unintended off-targets. Such drug promiscuity has often led to unwanted or unexplained drug reactions, resulting in side effects or drug repositioning opportunities. So it is always an important issue in pharmacology to identify potential drug-target interactions (DTI). However, DTI discovery by experiment remains a challenging task, due to high expense of time and resources. Many computational methods are therefore developed to predict DTI with high throughput biological and clinical data. Here, we initiatively demonstrate that the on-target and off-target effects could be characterized by drug-induced in vitro genomic expression changes, e.g. the data in Connectivity Map (CMap). Thus, unknown ligands of a certain target can be found from the compounds showing high gene-expression similarity to the known ligands. Then to clarify the detailed practice of CMap based DTI prediction, we objectively evaluate how well each target is characterized by CMap. The results suggest that (1) some targets are better characterized than others, so the prediction models specific to these well characterized targets would be more accurate and reliable; (2) in some cases, a family of ligands for the same target tend to interact with common off-targets, which may help increase the efficiency of DTI discovery and explain the mechanisms of complicated drug actions. In the present study, CMap expression similarity is proposed as a novel indicator of drug-target interactions. The detailed strategies of improving data quality by decreasing the batch effect and building prediction models are also effectively established. We believe the success in CMap can be further translated into other public and commercial data of genomic expression, thus increasing research productivity towards valid drug repositioning and minimal side effects.     

Structure of a Bimodular Botulinum Neurotoxin Complex Provides Insights into Its Oral Toxicity

Botulinum neurotoxins (BoNTs) are produced by Clostridium botulinum and cause the fatal disease botulism, a flaccid paralysis of the muscle. BoNTs are released together with several auxiliary proteins as progenitor toxin complexes (PTCs) to become highly potent oral poisons. Here, we report the structure of a ∼760 kDa 14-subunit large PTC of serotype A (L-PTC/A) and reveal insight into its absorption mechanism. Using a combination of X-ray crystallography, electron microscopy, and functional studies, we found that L-PTC/A consists of two structurally and functionally independent sub-complexes. A hetero-dimeric 290 kDa complex protects BoNT, while a hetero-dodecameric 470 kDa complex facilitates its absorption in the harsh environment of the gastrointestinal tract. BoNT absorption is mediated by nine glycan-binding sites on the dodecameric sub-complex that forms multivalent interactions with carbohydrate receptors on intestinal epithelial cells. We identified monosaccharides that blocked oral BoNT intoxication in mice, which suggests a new strategy for the development of preventive countermeasures for BoNTs based on carbohydrate receptor mimicry.     

Pathways-Driven Sparse Regression Identifies Pathways and Genes Associated with High-Density Lipoprotein Cholesterol in Two Asian Cohorts

Standard approaches to data analysis in genome-wide association studies (GWAS) ignore any potential functional relationships between gene variants. In contrast gene pathways analysis uses prior information on functional structure within the genome to identify pathways associated with a trait of interest. In a second step, important single nucleotide polymorphisms (SNPs) or genes may be identified within associated pathways. The pathways approach is motivated by the fact that genes do not act alone, but instead have effects that are likely to be mediated through their interaction in gene pathways. Where this is the case, pathways approaches may reveal aspects of a trait''s genetic architecture that would otherwise be missed when considering SNPs in isolation. Most pathways methods begin by testing SNPs one at a time, and so fail to capitalise on the potential advantages inherent in a multi-SNP, joint modelling approach. Here, we describe a dual-level, sparse regression model for the simultaneous identification of pathways and genes associated with a quantitative trait. Our method takes account of various factors specific to the joint modelling of pathways with genome-wide data, including widespread correlation between genetic predictors, and the fact that variants may overlap multiple pathways. We use a resampling strategy that exploits finite sample variability to provide robust rankings for pathways and genes. We test our method through simulation, and use it to perform pathways-driven gene selection in a search for pathways and genes associated with variation in serum high-density lipoprotein cholesterol levels in two separate GWAS cohorts of Asian adults. By comparing results from both cohorts we identify a number of candidate pathways including those associated with cardiomyopathy, and T cell receptor and PPAR signalling. Highlighted genes include those associated with the L-type calcium channel, adenylate cyclase, integrin, laminin, MAPK signalling and immune function.     

The Wilms Tumor Gene,Wt1, Is Critical for Mouse Spermatogenesis via Regulation of Sertoli Cell Polarity and Is Associated with Non-Obstructive Azoospermia in Humans

Azoospermia is one of the major reproductive disorders which cause male infertility in humans; however, the etiology of this disease is largely unknown. In the present study, six missense mutations of WT1 gene were detected in 529 human patients with non-obstructive azoospermia (NOA), indicating a strong association between WT1 mutation and NOA. The Wilms tumor gene, Wt1, is specifically expressed in Sertoli cells (SCs) which support spermatogenesis. To examine the functions of this gene in spermatogenesis, Wt1 was deleted in adult testis using Wt1^flox and Cre-ER^TM mice strains. We found that inactivation of Wt1 resulted in massive germ cell death and only SCs were present in most of the seminiferous tubules which was very similar to NOA in humans. In investigating the potential mechanism for this, histological studies revealed that the blood–testis barrier (BTB) was disrupted in Wt1 deficient testes. In vitro studies demonstrated that Wt1 was essential for cell polarity maintenance in SCs. Further studies found that the expression of cell polarity associated genes (Par6b and E-cadherin) and Wnt signaling genes (Wnt4, Wnt11) were downregulated in Wt1 deficient SCs, and that the expression of Par6b and E-cadherin was regulated by Wnt4. Our findings suggest that Wt1 is important in spermatogenesis by regulating the polarity of SCs via Wnt signaling pathway and that WT1 mutation is one of the genetic causes of NOA in humans.     

Hepatitis C virus inhibits AKT-tuberous sclerosis complex (TSC), the mechanistic target of rapamycin (MTOR) pathway,through endoplasmic reticulum stress to induce autophagy

Hepatitis C virus (HCV) is able to induce autophagy via endoplasmic reticulum (ER) stress, but the exact molecular signaling pathway is not well understood. We found that the activity of the mechanistic target of rapamycin complex 1 (MTORC1) was inhibited in Huh7 cells either harboring HCV-N (genotype 1b) full-genomic replicon or infected with JFH1 (genotype 2a) virus, which led to the activation of UNC-51-like kinase 1 (ULK1) and thus to autophagy. We then analyzed activity upstream of MTORC1, and found that both protein kinase, AMP-activated, α (PRKAA, including PRKAA1 and PRKAA2, also known as AMP-activated protein kinase, AMPKα) and AKT (refers to pan AKT, including three isoforms of AKT1-3, also known as protein kinase B, PKB) were inhibited by HCV infection. The inhibition of the AKT-TSC-MTORC1 pathway contributed to upregulating autophagy, but inhibition of PRKAA downregulated autophagy. The net effect on autophagy was from AKT, which overrode the inhibition effect from PRKAA. It was further found that HCV-induced ER stress was responsible for the inhibition of the AKT pathway. Metformin, a PRKAA agonist, inhibited HCV replication not only by activating PRKAA as previously reported, but also by activating AKT independently of the autophagy pathway. Taken together, our data suggested HCV inhibited the AKT-TSC-MTORC1 pathway via ER stress, resulting in autophagy, which may contribute to the establishment of the HCV-induced autophagy.     

The molecular basis for the inhibition of phosphodiesterase-4D by three natural resveratrol analogs. Isolation,molecular docking,molecular dynamics simulations,binding free energy,and bioassay

The phosphodiesterase-4 (PDE4) enzyme is a promising therapeutic target for several diseases. Our previous studies found resveratrol and moracin M to be natural PDE4 inhibitors. In the present study, three natural resveratrol analogs [pterostilbene, (E)-2′,3,5′,5-tetrahydroxystilbene (THSB), and oxyresveratrol] are structurally related to resveratrol and moracin M, but their inhibition and mechanism against PDE4 are still unclear. A combined method consisting of molecular docking, molecular dynamics (MD) simulations, binding free energy, and bioassay was performed to better understand their inhibitory mechanism. The binding pattern of pterostilbene demonstrates that it involves hydrophobic/aromatic interactions with Phe340 and Phe372, and forms hydrogen bond(s) with His160 and Gln369 in the active site pocket. The present work also reveals that oxyresveratrol and THSB can bind to PDE4D and exhibits less negative predicted binding free energies than pterostilbene, which was qualitatively validated by bioassay (IC₅₀ = 96.6, 36.1, and 27.0 μM, respectively). Additionally, a linear correlation (R² = 0.953) is achieved for five PDE4D/ligand complexes between the predicted binding free energies and the experimental counterparts approximately estimated from their IC₅₀ values (≈RT ln IC₅₀). Our results imply that hydrophobic/aromatic forces are the primary factors in explaining the mechanism of inhibition by the three products. Results of the study help to understand the inhibitory mechanism of the three natural products, and thus help the discovery of novel PDE4 inhibitors from resveratrol, moracin M, and other natural products.     

Developmental expression of P2X5 receptors in the mouse prenatal central and peripheral nervous systems

The functions of P2X purinoceptors (P2X1-7) in the nervous system of adults have been widely studied. However, little is known about their roles during embryonic development. Our previous work has reported an extensive expression of P2X5 receptors in the adult mouse central nervous system. In the present study, we have examined the expression pattern of P2X5 receptor mRNA and protein during prenatal development of the mouse nervous system (from embryonic day E8 to E17). P2X5 receptors appeared in the neural tube as early as E8 and were gradually confined to new-born neurons in the cortical plate and ventral horn of the spinal cord. Heavy signals for P2X5 receptors were also found in dorsal root ganglia (DRG), retina, olfactory epithelium, and nerve fibers in skeletal muscles. In conclusion, P2X5 receptors were strongly represented in the developing mouse nervous system. The transient high expression pattern of P2X5 receptors in epithelium-like structures suggests a role during early neurogenesis.     

Block of P2X7 receptors could partly reverse the delayed neuronal death in area CA1 of the hippocampus after transient global cerebral ischemia

Transient global ischemia (which closely resembles clinical situations such as cardiac arrest, near drowning or severe systemic hypotension during surgical procedures), often induces delayed neuronal death in the brain, especially in the hippocampal CA1 region. The mechanism of ischemia/reperfusion (I/R) injury is not fully understood. In this study, we have shown that the P2X7 receptor antagonist, BBG, reduced delayed neuronal death in the hippocampal CA1 region after I/R injury; P2X7 receptor expression levels increased before delayed neuronal death after I/R injury; inhibition of the P2X7 receptor reduced I/R-induced microglial microvesicle-like components, IL-1β expression, P38 phosphorylation, and glial activation in hippocampal CA1 region after I/R injury. These results indicate that antagonism of the P2X7 receptor and signaling pathways of microglial MV shedding, such as src-protein tyrosine kinase, P38 MAP kinase and A-SMase, might be a promising therapeutic strategy for clinical treatment of transient global cerebral I/R injury.     

Long-term tillage effects on soil organic carbon and dissolved organic carbon in a purple paddy soil of Southwest China

The impact of conservation tillage practices on soil carbon has been of great interest in recent years. Conservation tillage might have the potential to enhance soil carbon accumulation and alter the depth distribution of soil carbon compared to conventional tillage based systems. Changes in the soil organic carbon (SOC) as influenced by tillage, are more noticeable under long-term rather than short-term tillage practices. The objective of this study was to determine the impacts of long-term tillage on SOC and dissolved organic carbon (DOC) status after 19 years of four tillage treatments in a Hydragric Anthrosol. In this experiment four tillage systems included conventional tillage with rotation of rice and winter fallow system (CTF), conventional tillage with rotation of rice and rape system (CTR), no-till and ridge culture with rotation of rice and rape system (NT) and tillage and ridge culture with rotation of rice and rape system (TR). Soils were sampled in the spring of 2009 and sectioned into 0–10, 10–20, 20–30, 30–40, 40–50 and 50–60 cm depth, respectively.Tillage effect on SOC was observed, and SOC concentrations were much larger under NT than the other three tillage methods in all soil depths from 0 to 60 cm. The mean SOC concentration at 0–60 cm soil depth followed the sequence: NT (22.74 g kg^?1) > CTF (14.57 g kg^?1) > TR (13.10 g kg^?1) > CTR (11.92 g kg^?1). SOC concentrations under NT were significantly higher than TR and CTR (P < 0.01), and higher than CTF treatment (P < 0.05). The SOC storage was calculated on equivalent soil mass basis. Results showed that the highest SOC storage at 0–60 cm depth presented in NT, which was 158.52 Mg C ha^?1, followed by CTF (106.74 Mg C ha^?1), TR (93.11 Mg C ha^?1) and CTR (88.60 Mg C ha^?1). Compared with conventional tillage (CTF), the total SOC storage in NT increased by 48.51%, but decreased by 16.99% and 12.77% under CTR and TR treatments, respectively. The effect of tillage on DOC was significant at 0–10 cm soil layer, and DOC concentration was much higher under CTF than the other three treatments (P < 0.01). Throughout 0–60 cm soil depth, DOC concentrations were 32.92, 32.63, 26.79 and 22.10 mg kg^?1 under NT, CTF, CTR and TR, and the differences among the four treatments were not significant (P > 0.05). In conclusion, NT increased SOC concentration and storage compared to conventional tillage operation but not for DOC.