Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Chemistry and accumulation rates of Jurassic to early Cretaceous sediments from the central equatorial Pacific

Sedimentation in the central Pacific during the Jurassic and Early Cretaceous was dominated by abundant biogenic silica. A synthesis of the stratigraphy, lithology, petrology, and geochemistry of the radiolarites in Sites 801 and 800 documents the sedimentation processes and trends in the equatorial central Pacific from the Middle Jurassic through the Early Cretaceous. Paleolatitude and paleodepth reconstructions enable comparisons with previous DSDP sites and identification of the general patterns of sedimentation over a wide region of the Pacific. Clayey radiolarites dominated sedimentation on Pacific oceanic crust within tropical paleolatitudes from at least the latest Bathonian through Tithonian. Radiolarian productivity rose to a peak within 5° of the paleoequator, where accumulation rates of biogenic silica exceeded 1000 g/cm**2/m.y. Wavy-bedded radiolarian cherts developed in the upper Tithonian at Site 801 coinciding with the proximity of this site to the paleoequator. Ribbon-bedding of some radiolarian cherts exposed on Pacific margins may have formed from silicification of radiolarite deposited near the equatorial high-productivity zone where radiolarian/clay ratios were high. Silicification processes in sediments extensively mixed by bioturbation or enriched in clay or carbonate generally resulted in discontinuous bands or nodules of porcellanite or chert, e.g., a "knobby" radiolarite. Ribbon-bedded cherts require primary alternations of radiolarian-rich and clay-rich layers as an initial structural template, coupled with abundant biogenic silica in both layers. During diagenesis, migration of silica from clay-rich layers leaves radiolarian "ghosts" or voids, and the precipitation in adjacent radiolarite layers results in silicification of the inter-radiolarian matrix and infilling of radiolarian tests. Alternations of claystone and clay-rich radiolarian grainstone were deposited during the Callovian at Site 801 and during the Berriasian-Valanginian at Site 800, but did not silicify to form bedded chert. Carbonate was not preserved on the Pacific oceanic floor or spreading ridges during the Jurassic, perhaps due to an elevated level of dissolved carbon dioxide. During the Berriasian through Hauterivian, the carbonate compensation depth (CCD) descended to approximately 3500 m, permitting the accumulation of siliceous limestones at near-ridge sites. Carbonate accumulation rates exceeded 1500 g/cm**2/m.y. at sites above the CCD, yet there is no evidence of an equatorial carbonate bulge during the Early Cretaceous. In the Barremian and Aptian, the CCD rose, coincident with the onset of mid-plate volcanic activity. Abundance of Fe and Mn and the associated formation of authigenic Fe-smectite clays was a function of proximity to the spreading ridges, with secondary enrichments occurring during episodes of spreading-center reorganizations. Callovian radiolarite at Site 801 is anomalously depleted in Mn, which resulted either from inhibited precipitation of Mn-oxides by lower pH of interstitial waters induced by high dissolved oceanic CO2 levels or from diagenetic mobilization of Mn. Influx of terrigenous (eolian) clay apparently changed with paleolatitude and geological age. Cyclic variations in productivity of radiolarians and of nannofossils and in the influx of terrigenous clay are attributed to Milankovitch climatic cycles of precession (20,000 yr) and eccentricity (100,000 yr). Diagenetic redistribution of biogenic silica and carbonate enhanced the expression of this cyclic sedimentation. Jurassic and Lower Cretaceous sediments were deposited under oxygenated bottom-water conditions at all depths, accompanied by bioturbation and pervasive oxidation of organic carbon and metals. Despite the more "equable" climate conditions of the Mesozoic, the super-ocean of the Pacific experienced adequate deep-water circulation to prevent stagnation. Efficient nutrient recycling may have been a factor in the abundance of radiolarians in this ocean basin.

Core alignment and composite depth scale of DSDP Site 86-577

Deep Sea Drilling Project Site 577 on Shatsky Rise (North Pacific Ocean) recovered a series of cores at three holes that contain calcareous nannofossil ooze of latest Cretaceous (late Maastrichtian) through early Eocene age. Several important records have been generated using samples from these cores, but the stratigraphy has remained outdated and confusing. Here we revise the stratigraphy at Site 577. This includes refining several age datums, realigning cores in the depth domain, and placing all stratigraphic markers on a current time scale. The work provides a template for appropriately bringing latest Cretaceous and Paleogene data sets at old drill sites into current paleoceanographic literature for this time interval. While the Paleocene Eocene Thermal Maximum (PETM) lies within core gaps at Holes 577* and 577A, the sedimentary record at the site holds other important events and remains crucially relevant to understanding changes in oceanographic conditions from the latest Cretaceous through early Paleogene.

Age models and stable isotope analysis on sediment core Site 199-1218 from the equatorial Pacific

A 13-million-year continuous record of Oligocene climate from the equatorial Pacific reveals a pronounced "heartbeat" in the global carbon cycle and periodicity of glaciations. This heartbeat consists of 405,000-, 127,000-, and 96,000-year eccentricity cycles and 1.2-million-year obliquity cycles in periodically recurring glacial and carbon cycle events. That climate system response to intricate orbital variations suggests a fundamental interaction of the carbon cycle, solar forcing, and glacial events. Box modeling shows that the interaction of the carbon cycle and solar forcing modulates deep ocean acidity as well as the production and burial of global biomass. The pronounced 405,000-year eccentricity cycle is amplified by the long residence time of carbon in the oceans.

Miocene microflora of Belorado (Spain)

Neogene deposits from Belorado (Province de Burgos, Spain) in the NE border of the Duero Basin have been analyzed. The palynologic analysis of the samples suggests the existence of arboreal landscapes associated with the herbaceous cover creating belts of vegetation around restricted aquatic areas (marshy and lacustrine environments). The characteristics of the palynomorph assemblages allow to define a warm-template climate with strong seasons. The age of these deposits could be included in the middle Miocene (Aragonian).

Seawater carbonate chemistry and calcification rates of Porites rus and Hydrolithon onkodes in experiments of Moorea

Central to evaluating the effects of ocean acidification (OA) on coral reefs is understanding how calcification is affected by the dissolution of CO2 in sea water, which causes declines in carbonate ion concentration [CO3]2- and increases in bicarbonate ion concentration [HCO3]-. To address this topic, we manipulated [CO3]2- and [HCO3]- to test the effects on calcification of the coral Porites rus and the alga Hydrolithon onkodes, measured from the start to the end of a 15-day incubation, as well as in the day and night. [CO3]2- played a significant role in light and dark calcification of P. rus, whereas [HCO3]- mainly affected calcification in the light. Both [CO3]2- and [HCO3]- had a significant effect on the calcification of H. onkodes, but the strongest relationship was found with [CO3]2-. Our results show that the negative effect of declining [CO3]2- on the calcification of corals and algae can be partly mitigated by the use of [HCO3]- for calcification and perhaps photosynthesis. These results add empirical support to two conceptual models that can form a template for further research to account for the calcification response of corals and crustose coralline algae to OA.

(Table 1) Stratigraphic datums from ODP Site 130-806

A 2 m.y. oxygen isotope record of Globigerinoides sacculifer from the Ontong Java Plateau, based on cores from Ocean Drilling Program Leg 130, is dated by matching variations to an orbital template. The procedure allows us to present the most complete Quaternary record available for the western equatorial Pacific. The template-generating algorithm describes a balance between growth and melting of ice. Following basic Milankovitch theory, ice growth is taken as constant, while melting is taken to depend on summer insolation, current ice mass, and average past ice mass. Template settings must be changed once, between 1 and 1.2 Ma, to reflect a major shift in climate. Template fits are strikingly good over much of the record and can be used to detect and fill gaps from core breaks and other disturbances. One result of template dating is an exact age for the Brunhes-Matuyama magnetic reversal boundary, at 790+/-5 ka, as well as several other precise dates (900 ka for the middle Pleistocene climate shift; 1070, 1240, and 1450 ka for isotope stages 31, 37, and 47, respectively). Sedimentation rates fluctuate between 18 and 28 m/m.y., a ca. 400 ka cycle being the most prominent. Major anomalies arise within the transitional regime (1.2 to 1 Ma). The origin of the cycles is unknown; we propose productivity variations in the western equatorial Pacific.

Video analysis based vehicle detection and tracking using an MCMC sampling framework

This article presents a probabilistic method for vehicle detection and tracking through the analysis of monocular images obtained from a vehicle-mounted camera. The method is designed to address the main shortcomings of traditional particle filtering approaches, namely Bayesian methods based on importance sampling, for use in traffic environments. These methods do not scale well when the dimensionality of the feature space grows, which creates significant limitations when tracking multiple objects. Alternatively, the proposed method is based on a Markov chain Monte Carlo (MCMC) approach, which allows efficient sampling of the feature space. The method involves important contributions in both the motion and the observation models of the tracker. Indeed, as opposed to particle filter-based tracking methods in the literature, which typically resort to observation models based on appearance or template matching, in this study a likelihood model that combines appearance analysis with information from motion parallax is introduced. Regarding the motion model, a new interaction treatment is defined based on Markov random fields (MRF) that allows for the handling of possible inter-dependencies in vehicle trajectories. As for vehicle detection, the method relies on a supervised classification stage using support vector machines (SVM). The contribution in this field is twofold. First, a new descriptor based on the analysis of gradient orientations in concentric rectangles is dened. This descriptor involves a much smaller feature space compared to traditional descriptors, which are too costly for real-time applications. Second, a new vehicle image database is generated to train the SVM and made public. The proposed vehicle detection and tracking method is proven to outperform existing methods and to successfully handle challenging situations in the test sequences.

InN/InGaN multiple quantum wells emitting at 1.5 mu m grown by molecular beam epitaxy

This work reports on the growth by molecular beam epitaxy and characterization of InN/InGaN multiple quantum wells (MQWs) emitting at 1.5 μm. X-ray diffraction (XRD) spectra show satellite peaks up to the second order. Estimated values of well (3 nm) and barrier (9 nm) thicknesses were derived from transmission electron microscopy and the fit between experimental data and simulated XRD spectra. Transmission electron microscopy and XRD simulations also confirmed that the InGaN barriers are relaxed with respect to the GaN template, while the InN MQWs grew under biaxial compression on the InGaN barriers. Low temperature (14 K) photoluminescence measurements reveal an emission from the InN MQWs at 1.5 μm. Measurements as a function of temperature indicate the existence of localized states, probably due to InN quantum wells’ thickness fluctuations as observed by transmission electron microscopy.

InN/GaN heterojunction electrical behavior

Indium nitride (InN) has been the subject of intense research in recent years. Some of its most attractive features are its excellent transport properties such as its small band edge electron effective mass, high electron mobilities and peak drift velocities, and high frequency transient drift velocity oscillations [1]. These suggest enormous potential applications for InN in high frequency electronic devices. But to date the high unintentional bulk electron concentration (n~1018 cm-3) of undoped InN samples and the surface electron accumulation layer make it a hard task to create a reliable metalsemiconductor Schottky barrier. Some attempts have been made to overcome this problem by means of material oxidation [2] or deposition of insulators [3]. In this work we present a way to obtain an electrical rectification behaviour by means of heterojunction growth. Due to the big band gap differences among nitride semiconductors, it’s possible to create a structure with high band offsets. In InN/GaN heterojunctions, depending on the GaN doping, the magnitude of conduction and valence band offset are critical parameters which allow distinguishing among different electrical behaviours. The earliest estimate of the valence band offset at an InN–GaN heterojunction in a wurtzite structure was measured to be ~0.85 eV [4], while the Schottky barrier heights were determined to be ~ 1,4 eV [5].We grew In-face InN layer with varying thickness (between 150 nm and 1 mm) by plasma assisted molecular beam epitaxy (PA-MBE) on GaNntemplates (GaN/Al2O3), with temperatures ranging between 300°C and 450°C. The different doping in GaN template (Si doping, Fe doping and Mg doping) results in differences in band alignments of the two semiconductors changing electrical barriers for carriers and consequently electrical conduction behaviour. The processing of the devices includes metallization of the ohmic contacts on InN and GaN, for which we used Ti/Al/Ni/Au. Whereas an ohmic contact on InN is straightforward, the main issue was the fabrication of the contact on GaN due to the very low decomposition temperature of InN. A standard ohmic contact on GaN is generally obtained by high temperature rapid thermal annealing (RTA), typically done between 500ºC and 900ºC[6]. In this case, the limitation due to the presence of In-face InN imposes an upper limit on the temperature for the thermal annealing process and ohmic contact formation of about 450°C. We will present results on the morphology of the InN layers by X-Ray diffraction and SEM, and electrical measurements, in particular current-voltage and capacitance-voltage characteristics.

Computerized Three-Dimmensional Craniofacial Reconstruction from Skulls Based on Landmarks.

Human identification from a skull is a critical process in legal and forensic medicine, specially when no other means are available. Traditional clay-based methods attempt to generate the human face, in order to identify the corresponding person. However, these reconstructions lack of objectivity and consistence, since they depend on the practitioner. Current computerized techniques are based on facial models, which introduce undesired facial features when the final reconstruction is built. This paper presents an objective 3D craniofacial reconstruction technique, implemented in a graphic application, without using any facial template. The only information required by the software tool is the 3D image of the target skull and three parameters: age, gender and Body Mass Index (BMI) of the individual. Complexity is minimized, since the application database only consists of the anthropological information provided by soft tissue depth values in a set of points of the skull.

Study of the effect of blurriness in image acquisition for Hand Biometrics in Mobile Devices

This paper presents a study on the effect of blurred images in hand biometrics. Blurred images simulates out-of-focus effects in hand image acquisition, a common consequence of unconstrained, contact-less and platform-free hand biometrics in mobile devices. The proposed biometric system presents a hand image segmentation based on multiscale aggregation, a segmentation method invariant to different changes like noise or blurriness, together with an innovative feature extraction and a template creation, oriented to obtain an invariant performance against blurring effects. The results highlight that the proposed system is invariant to some low degrees of blurriness, requiring an image quality control to detect and correct those images with a high degree of blurriness. The evaluation has considered a synthetic database created based on a publicly available database with 120 individuals. In addition, several biometric techniques could benefit from the approach proposed in this paper, since blurriness is a very common effect in biometric techniques involving image acquisition.

Secure access control by means of human stress detection

This paper proposes a stress detection system based on fuzzy logic and the physiological signals heart rate and galvanic skin response. The main contribution of this method relies on the creation of a stress template, collecting the behaviour of previous signals under situations with a different level of stress in each individual. The creation of this template provides an accuracy of 99.5% in stress detection, improving the results obtained by current pattern recognition techniques like GMM, k-NN, SVM or Fisher Linear Discriminant. In addition, this system can be embedded in security systems to detect critical situations in accesses as cross-border control. Furthermore, its applications can be extended to other fields as vehicle driver state-of-mind management, medicine or sport training.

A comparative study on unconstrained hand biometrics

Biometrics applied to mobile devices are of great interest for security applications. Daily scenarios can benefit of a combination of both the most secure systems and most simple and extended devices. This document presents a hand biometric system oriented to mobile devices, proposing a non-intrusive, contact-less acquisition process where final users should take a picture of their hand in free-space with a mobile device without removals of rings, bracelets or watches. The main contribution of this paper is threefold: firstly, a feature extraction method is proposed, providing invariant hand measurements to previous changes; second contribution consists of providing a template creation based on hand geometric distances, requiring information from only one individual, without considering data from the rest of individuals within the database; finally, a proposal for template matching is proposed, minimizing the intra-class similarity and maximizing the inter-class likeliness. The proposed method is evaluated using three publicly available contact-less, platform-free databases. In addition, the results obtained with these databases will be compared to the results provided by two competitive pattern recognition techniques, namely Support Vector Machines (SVM) and k-Nearest Neighbour, often employed within the literature. Therefore, this approach provides an appropriate solution to adapt hand biometrics to mobile devices, with an accurate results and a non-intrusive acquisition procedure which increases the overall acceptance from the final user.

Analysis of pattern recognition techniques for in-air signature biometrics

As a result of advances in mobile technology, new services which benefit from the ubiquity of these devices are appearing. Some of these services require the identification of the subject since they may access private user information. In this paper, we propose to identify each user by drawing his/her handwritten signature in the air (in-airsignature). In order to assess the feasibility of an in-airsignature as a biometric feature, we have analysed the performance of several well-known patternrecognitiontechniques—Hidden Markov Models, Bayes classifiers and dynamic time warping—to cope with this problem. Each technique has been tested in the identification of the signatures of 96 individuals. Furthermore, the robustness of each method against spoofing attacks has also been analysed using six impostors who attempted to emulate every signature. The best results in both experiments have been reached by using a technique based on dynamic time warping which carries out the recognition by calculating distances to an average template extracted from several training instances. Finally, a permanence analysis has been carried out in order to assess the stability of in-airsignature over time.