Approaching the CAPRI challenge with an efficient geometry-based docking

The last 3 rounds (3-5) of CAPRI included a wide range of docking targets. Several targets were especially challenging, since they involved large-scale movements and symmetric rearrangement, while others were based on homology models. We have approached the targets with a variety of geometry-based docking algorithms that include rigid docking, symmetric docking, and flexible docking with symmetry constraints. For all but 1 docking target, we were able to submit at least 1 acceptable quality prediction. Here, we detail for each target the prediction methods used and the specific biological data employed, and supply a retrospective analysis of the results. We highlight the advantages of our techniques, which efficiently exploit the geometric shape complementarity properties of the interaction. These enable them to run only few minutes on a standard PC even for flexible docking, thus proving their scalability toward computational genomic scale experiments. We also outline the major required enhancements, such as the introduction of side-chain position refinement and the introduction of flexibility for both docking partners.     

Intracellular delivery of phosphatidylinositol (3,4,5)-trisphosphate causes incorporation of glucose transporter 4 into the plasma membrane of muscle and fat cells without increasing glucose uptake

Insulin stimulates glucose uptake into muscle and fat cells by translocating glucose transporter 4 (GLUT4) to the cell surface, with input from phosphatidylinositol (PI) 3-kinase and its downstream effector Akt/protein kinase B. Whether PI 3,4,5-trisphosphate (PI(3,4,5)P(3)) suffices to produce GLUT4 translocation is unknown. We used two strategies to deliver PI(3,4,5)P(3) intracellularly and two insulin-sensitive cell lines to examine Akt activation and GLUT4 translocation. In 3T3-L1 adipocytes, the acetoxymethyl ester of PI(3,4,5)P(3) caused GLUT4 migration to the cell periphery and increased the amount of plasma membrane-associated phospho-Akt and GLUT4. Intracellular delivery of PI(3,4,5)P(3) using polyamine carriers also induced translocation of myc-tagged GLUT4 to the surface of intact L6 myoblasts, demonstrating membrane insertion of the transporter. GLUT4 translocation caused by carrier-delivered PI(3,4,5)P(3) was not reproduced by carrier-PI 4,5-bisphosphate or carrier alone. Like insulin, carrier-mediated delivery of PI(3,4,5)P(3) elicited redistribution of perinuclear GLUT4 and Akt phosphorylation at the cell periphery. In contrast to its effect on GLUT4 mobilization, delivered PI(3,4,5)P(3) did not increase 2-deoxyglucose uptake in either L6GLUT4myc myoblasts or 3T3-L1 adipocytes. The ability of exogenously delivered PI(3,4,5)P(3) to augment plasma membrane GLUT4 content without increasing glucose uptake suggests that input at the level of PI 3-kinase suffices for GLUT4 translocation but is insufficient to stimulate glucose transport.     

QSOX contains a pseudo-dimer of functional and degenerate sulfhydryl oxidase domains

Quiescin sulfhydryl oxidase (QSOX) catalyzes formation of disulfide bonds between cysteine residues in substrate proteins. Human QSOX1 is a multi-domain, monomeric enzyme containing a module related to the single-domain sulfhydryl oxidases of the Erv family. A partial QSOX1 crystal structure reveals a single-chain pseudo-dimer mimicking the quaternary structure of Erv enzymes. However, one pseudo-dimer “subunit” has lost its cofactor and catalytic activity. In QSOX evolution, a further concatenation to a member of the protein disulfide isomerase family resulted in an enzyme capable of both disulfide formation and efficient transfer to substrate proteins.     

Shape similarity measures for the design of small RNA switches

Conformational switching in the secondary structure of RNAs has recently attracted considerable attention, fostered by the discovery of 'riboswitches' in living organisms. These are genetic control elements that were found in bacteria and offer a unique regulation mechanism based on switching between two highly stable states, separated by an energy barrier between them. In riboswitches, the energy barrier is crossed by direct metabolite binding, which facilitates regulation by allosteric means. However, other event triggers can cause switching to occur, such as single-point mutations and slight variations in temperature. Examples of switches with these event triggers have already been reported experimentally in the past. Here, the goal is to computationally design small RNA switches that rely on these triggers. Towards this end, our computer simulations utilize a variety of different similarity measures to assess the distances between an initial state and triggered states, based on the topology of the secondary structure itself. We describe these combined similarity measures that rely on both coarse-grained and fine-grained graph representations of the RNA secondary structure. As a result of our simulations, we provide some candidate sequences of approximately 30-50 nt, along with the exact triggers that drive the switching. The event triggers under consideration can be modelled by Zuker's mfold or the Vienna package. The proposed methodology that rely on shape measures can further be used to computationally generate more candidates by simulating various event triggers and calculating their effect on the shape.     

A developmentally regulated lipocalin-like gene is overexpressed in Tomato yellow leaf curl virus-resistant tomato plants upon virus inoculation, and its silencing abolishes resistance

To discover genes involved in tomato resistance to Tomato yellow leaf curl virus (TYLCV), we previously compared cDNA libraries from susceptible (S) and resistant (R) tomato lines. Among the genes preferentially expressed in R plants and upregulated by TYLCV infection was a gene encoding a lipocalin-like protein. This gene was termed Solanum lycopersicum virus resistant/susceptible lipocalin (SlVRSLip). The SlVRSLip structural gene sequence of R and S plants was identical. SlVRSLip was expressed in leaves during a 15-day window starting about 40?days after sowing (20?days after planting). SlVRSLip was upregulated by Bemisia tabaci (the TYLCV vector) feeding on R plant leaves, and even more strongly upregulated following whitefly-mediated TYLCV inoculation. Silencing of SlVRSLip in R plants led to the collapse of resistance upon TYLCV inoculation and to a necrotic response along the stem and petioles accompanied by ROS production. Contrary to previously identified tomato lipocalin gene DQ222981, SlVRSLip was not regulated by cold, nor was it regulated by heat or salt. The expression of SlVRSLip was inhibited in R plants in which the hexose transporter gene LeHT1 was silenced. In contrast, the expression of LeHT1 was not inhibited in SlVRSLip-silenced R plants. Hence, in the hierarchy of the gene network conferring TYLCV resistance, SlVRSLip is downstream of LeHT1. Silencing of another gene involved in resistance, a Permease-I like protein, did not affect the expression of SlVRSLip and LeHT1; expression of the Permease was not affected by silencing SlVRSLip or LeHT1, suggesting that it does not belong to the same network. The triple co-silencing of SlVRSLip, LeHT1 and Permease provoked an immediate cessation of growth of R plants upon infection and the accumulation of large amounts of virus. SlVRSLip is the first lipocalin-like gene shown to be involved in resistance to a plant virus.     

Hypoxia associated NMDA receptor 2 subunit composition: developmental comparison between the hypoxia-tolerant subterranean mole-rat, Spalax, and the hypoxia-sensitive rat

Vertebrate brains are sensitive to oxygen depletion, which may lead to cell death. Hypoxia sensitivity originates from the high intrinsic rate of ATP consumption of brain tissue, accompanied by the release of glutamate, leading to the opening of ionotropic glutamate receptors, such as N-methyl-D-aspartate (NMDA) receptors (NMDARs). The relative expression levels of the four NMDAR-2 (NR2) subunits change during mammalian development with higher levels of units NR2B and NR2D observed during early development and correlated with hypoxic tolerance during embryonic and neonatal stages of development. Higher levels of NR2D are also abundant in brains of hypoxia tolerant species such as the crucian carp. The subterranean mole-rat, Spalax spends its life underground in sealed burrows and has developed a wide range of adaptations to this special niche including hypoxia-tolerance. In this study, we compared the in vivo mRNA expression of NR2 subunits in the brains of embryonic, neonatal and adult Spalax and rat. Our results demonstrate that under normoxic conditions, mRNA levels of NR2D are higher in Spalax than in rat at all developmental stages studied and are similar to levels in neonatal rat and in other hypoxia/anoxia tolerant species. Furthermore, under hypoxia Spalax NR2D mRNA levels increase while no response was observed in rat. Similarly, hypoxia induces an increase in mRNA levels of Spalax NR2A, claimed to promote neuronal survival. We suggest that indeed the proportional combinations of NMDAR-2 subunits contribute to the ability of the Spalax brain to cope with hypoxic environments.     

Can we use diffusion MRI as a bio-marker of neurodegenerative processes?

Magnetic resonance imaging (MRI) is an imaging technique with a rapidly expanding application range. This methodology, which relies on quantum physics and substance magnetic properties, is now being routinely used in the clinics and medical research. With the advent of measuring functional brain activity with MRI (functional MRI), this methodology has reached a larger section of the neuroscience community (e.g. psychologists, neurobiologists). In the past, the use of MRI as a biomarker or as an assay to probe tissue pathophysiological condition was limited. However, with the new applications of MRI: molecular imaging, contrast-enhanced imaging and diffusion imaging, MRI is turning into a powerful tool for in vivo characterization of tissue pathophysiology. This review focuses on the diffusion MRI. Although it only measures the averaged Brownian translational motion of water molecules, using different analysis schemes, one can extract a wide range of quantitative indices that represent tissue morphology and compartmentalization. Statistical and visualization routines help to relate these indices to biologically relevant measures such as cell density, water content and size distribution. The aim of this review is to shed light on the potential of this methodology to be used in biological research. To that end, this review is intended for the non-MRI specialists who wish to pursue biological research with this methodology. We will overview the current applications of diffusion MRI and its relation to cellular biology of brain tissue.