Data visualization during the early stages of drug discovery

Multidimensional compound optimization is a new paradigm in the drug discovery process, yielding efficiencies during early stages and reducing attrition in the later stages of drug development. The success of this strategy relies heavily on understanding this multidimensional data and extracting useful information from it. This paper demonstrates how principled visualization algorithms can be used to understand and explore a large data set created in the early stages of drug discovery. The experiments presented are performed on a real-world data set comprising biological activity data and some whole-molecular physicochemical properties. Data visualization is a popular way of presenting complex data in a simpler form. We have applied powerful principled visualization methods, such as generative topographic mapping (GTM) and hierarchical GTM (HGTM), to help the domain experts (screening scientists, chemists, biologists, etc.) understand and draw meaningful decisions. We also benchmark these principled methods against relatively better known visualization approaches, principal component analysis (PCA), Sammon's mapping, and self-organizing maps (SOMs), to demonstrate their enhanced power to help the user visualize the large multidimensional data sets one has to deal with during the early stages of the drug discovery process. The results reported clearly show that the GTM and HGTM algorithms allow the user to cluster active compounds for different targets and understand them better than the benchmarks. An interactive software tool supporting these visualization algorithms was provided to the domain experts. The tool facilitates the domain experts by exploration of the projection obtained from the visualization algorithms providing facilities such as parallel coordinate plots, magnification factors, directional curvatures, and integration with industry standard software. © 2006 American Chemical Society.

Interocular masking and summation indicate two stages of divisive contrast gain control

Our understanding of early spatial vision owes much to contrast masking and summation paradigms. In particular, the deep region of facilitation at low mask contrasts is thought to indicate a rapidly accelerating contrast transducer (eg a square-law or greater). In experiment 1, we tapped an early stage of this process by measuring monocular and binocular thresholds for patches of 1 cycle deg-1 sine-wave grating. Threshold ratios were around 1.7, implying a nearly linear transducer with an exponent around 1.3. With this form of transducer, two previous models (Legge, 1984 Vision Research 24 385 - 394; Meese et al, 2004 Perception 33 Supplement, 41) failed to fit the monocular, binocular, and dichoptic masking functions measured in experiment 2. However, a new model with two-stages of divisive gain control fits the data very well. Stage 1 incorporates nearly linear monocular transducers (to account for the high level of binocular summation and slight dichoptic facilitation), and monocular and interocular suppression (to fit the profound 42 Oral presentations: Spatial vision Thursday dichoptic masking). Stage 2 incorporates steeply accelerating transduction (to fit the deep regions of monocular and binocular facilitation), and binocular summation and suppression (to fit the monocular and binocular masking). With all model parameters fixed from the discrimination thresholds, we examined the slopes of the psychometric functions. The monocular and binocular slopes were steep (Weibull ߘ3-4) at very low mask contrasts and shallow (ߘ1.2) at all higher contrasts, as predicted by all three models. The dichoptic slopes were steep (ߘ3-4) at very low contrasts, and very steep (ß>5.5) at high contrasts (confirming Meese et al, loco cit.). A crucial new result was that intermediate dichoptic mask contrasts produced shallow slopes (ߘ2). Only the two-stage model predicted the observed pattern of slope variation, so providing good empirical support for a two-stage process of binocular contrast transduction. [Supported by EPSRC GR/S74515/01]

β-amyloid plaques:Stages in life history or independent origin?

Several types of discrete β-amyloid (Aβ) deposit or senile plaque have been identified in the brains of individuals with Alzheimer's disease and Down's syndrome. The majority of these plaques can be classified into four morphological types: diffuse, primitive, classic and compact. Two hypotheses have been proposed to account for these plaques. Firstly, that the diffuse, primitive, classic and compact plaques develop in sequence and represent stages in the life history of a single plaque type. Secondly, that the different Aβ plaques develop independently and therefore, unique factors are involved in the formation of each type. To attempt to distinguish between these hypotheses, the morphology, ultrastructure, composition, and spatial distribution in the brain of the four types of plaque were compared. Although some primitive plaques may develop from diffuse plaques, the evidence suggests that a unique combination of factors is involved in the pathogenesis of each plaque type and, therefore, supports the hypothesis that the major types of Aβ plaque develop independently.

Diagenesis and aquifer characteristics of the Sherwood Sandstone Group in North Nottinghamshire

The Sherwood Sandstone Group forms an important aquifer in Eastern England, which in North Nottinghamshire comprises the Nottingham Castle and Lenton Sandstone Formations. The aquifer is formed by an alluvial red-bed sequence dominated by medium-coarse grained sandstones which are texturally immature to submature and have only been subjected to shallow burial diagenesis. These sandstones reached the mature stage of the meso diagenetic regime, and four stages are recognized in their diagenetic history depending upon the physical/chemical processes prevailing and the subsequent effect on porosity and permeability. Stage "One" represents changes including dissolution of unstable silicates, clay replacement, red colouration and precipitation of authigenic minerals (quartz, feldspar, illite, l/S, kaolinite, dolomite, ferroan calcite, calcite). The net result of these changes was porosity reduction. Stage "Two" included changes due to mechanical compaction which resulted in minor porosity reduction. Stage "Three" was the main phase of secondary porosity enhancement. Stage "Four" represents changes taking place in the present groundwater where porosity and permeability may have been increased by dissolution and partly reduced by kaolinite precipitation. Porosity measured by water-resaturation and Hg-injection gave average values of 25.63% and 24.85% respectively. The results are comparable and showed marked correlation especially in highly porous/permeable rocks. Porosity measurements from photomicrographs were markedly offset from laboratory results. Horizontal Kw ranged between 1.43 x 10-5 and 1.13 x 10-1 mm/sec, with an average of  1.68 x 10-2 mm/sec. The estimated KHg ranged between 7.29 x 10-6 and 6.99 x 10-2 mm/sec with an average of 1.47 x 10-2 mm/sec. Both results are significantly correlated for highly porous/permeable rocks. The hydraulic properties are highly dependent upon the diagenetic properties (as most of the pores present are of secondary origin) as well as the pore size distribution. The chemistry of these groundwaters indicates that they are under-saturated with respect to dolomite, calcite, K-feldspar, l/S clay, and montmorillonite. The precipitation of kaolinite,and to a lesser extent illite, is favoured in the present groundwater regime.

Cardiac dysfunction and cell damage across clinical stages of severity in growth-restricted fetuses

Objective - The purpose of this study was to assess cardiac function and cell damage in intrauterine growth-restricted (IUGR) fetuses across clinical Doppler stages of deterioration. Study Design - One hundred twenty appropriate-for-gestational-age and 81 IUGR fetuses were classified in stages 1/2/3 according umbilical artery present/absent/reversed end-diastolic blood flow, respectively. Cardiac function was assessed by modified-myocardial performance index, early-to-late diastolic filling ratios, cardiac output, and cord blood B-type natriuretic peptide; myocardial cell damage was assessed by heart fatty acid–binding protein, troponin-I, and high-sensitivity C-reactive protein. Results - Modified-myocardial performance index, blood B-type natriuretic peptide, and early-to-late diastolic filling ratios were increased in a stage-dependent manner in IUGR fetuses, compared with appropriate-for-gestational-age fetuses. Heart fatty acid–binding protein levels were higher in IUGR fetuses at stage 3, compared with control fetuses. Cardiac output, troponin-I, and high-sensitivity C-reactive protein did not increase in IUGR fetuses at any stage. Conclusion - IUGR fetuses showed signs of cardiac dysfunction from early stages. Cardiac dysfunction deteriorates further with the progression of fetal compromise, together with the appearance of biochemical signs of cell damage.

Regional competitiveness, economic growth and stages of development

This paper positions the concept of regional competitiveness within theories concerning regional economic growth and stages of economic development. It examines the sources of regional competitiveness encompassing an analysis based on the particular stage of economic development that the nations within which regions are situated have reached. As a means to achieve this, the paper undertakes an empirical analysis of data stemming from the World Competitiveness Index of Regions, and identifies regional competitiveness as a dual concept that explains relative differences in rates of economic development across regions, as well as an understanding of the future economic growth trajectories of regions at a similar stage of economic development. As with endogenous growth and development theory, the notion of regional competitiveness presented here places knowledge, innovation and entrepreneurship at the forefront of conceptualisations of regional economic differentiation.

Immunological synapse: a mathematical model of the bond formation process:mathematical modelling of the immature synapse

The cell:cell bond between an immune cell and an antigen presenting cell is a necessary event in the activation of the adaptive immune response. At the juncture between the cells, cell surface molecules on the opposing cells form non-covalent bonds and a distinct patterning is observed that is termed the immunological synapse. An important binding molecule in the synapse is the T-cell receptor (TCR), that is responsible for antigen recognition through its binding with a major-histocompatibility complex with bound peptide (pMHC). This bond leads to intracellular signalling events that culminate in the activation of the T-cell, and ultimately leads to the expression of the immune eector function. The temporal analysis of the TCR bonds during the formation of the immunological synapse presents a problem to biologists, due to the spatio-temporal scales (nanometers and picoseconds) that compare with experimental uncertainty limits. In this study, a linear stochastic model, derived from a nonlinear model of the synapse, is used to analyse the temporal dynamics of the bond attachments for the TCR. Mathematical analysis and numerical methods are employed to analyse the qualitative dynamics of the nonequilibrium membrane dynamics, with the specic aim of calculating the average persistence time for the TCR:pMHC bond. A single-threshold method, that has been previously used to successfully calculate the TCR:pMHC contact path sizes in the synapse, is applied to produce results for the average contact times of the TCR:pMHC bonds. This method is extended through the development of a two-threshold method, that produces results suggesting the average time persistence for the TCR:pMHC bond is in the order of 2-4 seconds, values that agree with experimental evidence for TCR signalling. The study reveals two distinct scaling regimes in the time persistent survival probability density prole of these bonds, one dominated by thermal uctuations and the other associated with the TCR signalling. Analysis of the thermal fluctuation regime reveals a minimal contribution to the average time persistence calculation, that has an important biological implication when comparing the probabilistic models to experimental evidence. In cases where only a few statistics can be gathered from experimental conditions, the results are unlikely to match the probabilistic predictions. The results also identify a rescaling relationship between the thermal noise and the bond length, suggesting a recalibration of the experimental conditions, to adhere to this scaling relationship, will enable biologists to identify the start of the signalling regime for previously unobserved receptor:ligand bonds. Also, the regime associated with TCR signalling exhibits a universal decay rate for the persistence probability, that is independent of the bond length.