Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Distribution of events of positive selection and population differentiation in a metabolic pathway: the case of asparagine N-glycosylation

Background: Asparagine N-Glycosylation is one of the most important forms of protein post-translational modification in eukaryotes. This metabolic pathway can be subdivided into two parts: an upstream sub-pathway required for achieving proper folding for most of the proteins synthesized in the secretory pathway, and a downstream sub-pathway required to give variability to trans-membrane proteins, and involved in adaptation to the environment andinnate immunity. Here we analyze the nucleotide variability of the genes of this pathway in human populations, identifying which genes show greater population differentiation and which genes show signatures of recent positive selection. We also compare how these signals are distributed between the upstream and the downstream parts of the pathway, with the aim of exploring how forces of population differentiation and positive selection vary among genes involved in the same metabolic pathway but subject to different functional constraints. Results:Our results show that genes in the downstream part of the pathway are more likely to show a signature of population differentiation, while events of positive selection are equally distributed among the two parts of the pathway. Moreover, events of positive selection arefrequent on genes that are known to be at bifurcation points, and that are identified as beingin key position by a network-level analysis such as MGAT3 and GCS1.Conclusions: These findings indicate that the upstream part of the Asparagine N-Glycosylation pathway has lower diversity among populations, while the downstream part is freer to tolerate diversity among populations. Moreover, the distribution of signatures of population differentiation and positive selection can change between parts of a pathway, especially between parts that are exposed to different functional constraints. Our results support the hypothesis that genes involved in constitutive processes can be expected to show lower population differentiation, while genes involved in traits related to the environment should show higher variability. Taken together, this work broadens our knowledge on how events of population differentiation and of positive selection are distributed among different parts of a metabolic pathway.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

Evaluation of Different Methods to Calculate Heavy-Truck VMT, Final Report, December 2004

Reliable estimates of heavy-truck volumes are important in a number of transportation applications. Estimates of truck volumes are necessary for pavement design and pavement management. Truck volumes are important in traffic safety. The number of trucks on the road also influences roadway capacity and traffic operations. Additionally, heavy vehicles pollute at higher rates than passenger vehicles. Consequently, reliable estimates of heavy-truck vehicle miles traveled (VMT) are important in creating accurate inventories of on-road emissions. This research evaluated three different methods to calculate heavy-truck annual average daily traffic (AADT) which can subsequently be used to estimate vehicle miles traveled (VMT). Traffic data from continuous count stations provided by the Iowa DOT were used to estimate AADT for two different truck groups (single-unit and multi-unit) using the three methods. The first method developed monthly and daily expansion factors for each truck group. The second and third methods created general expansion factors for all vehicles. Accuracy of the three methods was compared using n-fold cross-validation. In n-fold cross-validation, data are split into n partitions, and data from the nth partition are used to validate the remaining data. A comparison of the accuracy of the three methods was made using the estimates of prediction error obtained from cross-validation. The prediction error was determined by averaging the squared error between the estimated AADT and the actual AADT. Overall, the prediction error was the lowest for the method that developed expansion factors separately for the different truck groups for both single- and multi-unit trucks. This indicates that use of expansion factors specific to heavy trucks results in better estimates of AADT, and, subsequently, VMT, than using aggregate expansion factors and applying a percentage of trucks. Monthly, daily, and weekly traffic patterns were also evaluated. Significant variation exists in the temporal and seasonal patterns of heavy trucks as compared to passenger vehicles. This suggests that the use of aggregate expansion factors fails to adequately describe truck travel patterns.

7

Climate change data and predictions for the Himalayas are very sparse and uncertain, characterized by a ?Himalayan data gap? and difficulties in predicting changes due to topographic complexity. A few reliable studies and climate change models for Nepal predict considerable changes: shorter monsoon seasons, more intensive rainfall patterns, higher temperatures, and drought. These predictions are confirmed by farmers who claim that temperatures have been increasing for the past decade and wonder why the rains have ?gone mad.? The number of hazard events, notably droughts, floods, and landslides are increasing and now account for approximately 100 deaths in Nepal annually. Other effects are drinking water shortages and shifting agricultural patterns, with many communities struggling to meet basic food security before climatic conditions started changing. The aim of this paper is to examine existing gaps between current climate models and the realities of local development planning through a case study on flood risk and drinking water management for the Municipality of Dharan in Eastern Nepal. This example highlights current challenges facing local-level governments, namely, flood and landslide mitigation, providing basic amenities ? especially an urgent lack of drinking water during the dry season ? poor local planning capacities, and limited resources. In this context, the challenge for Nepal will be to simultaneously address increasing risks caused by hazard events alongside the omnipresent food security and drinking water issues in both urban and rural areas. Local planning is needed that integrates rural development and disaster risk reduction (DRR) with knowledge about climate change considerations. The paper concludes with a critical analysis of climate change modeling and the gap between scientific data and low-tech and low capacities of local planners to access or implement adequate adaptation measures. Recommendations include the need to bridge gaps between scientific models, the local political reality and local information needs.

Bridging the Faraoni and Selli oceanic anoxic events: late Hauterivian to early Aptian dysaerobic to anaerobic phases in the Tethys

A detailed geochemical analysis was performed on the upper part of the Maiolica Formation in the Breggia (southern Switzerland) and Capriolo sections (northern Italy). The analysed sediments consist of well-bedded, partly siliceous, pelagic carbonate, which lodges numerous thin, dark and organic-rich layers. Stable-isotope, phosphorus, organic-carbon and a suite of redox-sensitive trace-element contents (RSTE: Mo, U, Co, V and As) were measured. The RSTE pattern and C-org:P-tot ratios indicate that most organic-rich layers were deposited under dysaerobic rather than anaerobic conditions and that latter conditions were likely restricted to short intervals in the latest Hauterivian, the early Barremian and the pre-Selli early Aptian. Correlations are both possible with organic-rich intervals in central Italy (the Gorgo a Cerbara section) and the Boreal Lower Saxony Basin, as well as with the facies and drowning pattern in the Helvetic segment of the northern Tethyan carbonate platform. Our data and correlations suggest that the latest Hauterivian witnessed the progressive installation of dysaerobic conditions in the Tethys, which went along with the onset in sediment condensation, phosphogenesis and platform drowning on the northern Tethyan margin, and which culminated in the Faraoni anoxic episode. This episode is followed by further episodes of dysaerobic conditions in the Tethys and the Lower Saxony Basin, which became more frequent and progressively stronger in the late early Barremian. Platform drowning persisted and did not halt before the latest early Barremian. The late Barremian witnessed diminishing frequencies and intensities in dysaerobic conditions, which went along with the progressive installation of the Urgonian carbonate platform. Near the Barremian-Aptian boundary, the increasing density in dysaerobic episodes in the Tethyan and Lower Saxony Basins is paralleled by a change towards heterozoan carbonate production on the northern Tethyan shelf. The following return to more oxygenated conditions is correlated with the second phase of Urgonian platform growth and the period immediately preceding and corresponding to the Selli anoxic episode is characterised by renewed platform drowning and the change to heterozoan carbonate production. Changes towards more humid climate conditions were the likely cause for the repetitive installation of dys- to anaerobic conditions in the Tethyan and Boreal basins and the accompanying changes in the evolution of the carbonate platform towards heterozoan carbonate-producing ecosystems and platform drowning.

The Toarcian in the Subbetic basin (southern Spain): Bio-events (ammonite and calcareous nannofossils) and carbon-isotope stratigraphy

A detailed carbon-isotope stratigraphic study for the uppermost Pliensbachian lowermost Aalenian interval in the Median Subbetic palaeogeographic domain (External zones of the Betic Cordillera, southern Spain) has been carried out. During the Early Jurassic, the Median Subbetic, which represents a typical basin of the Hispanic Corridor connecting the Tethys and the Eastern Pacific, was located in the westernmost Tethys. The analyzed sections encompass the entire Toarcian stage as represented in the southern Iberian palaeomargin. Rocks are mainly rhythmic sequences of grey marls and marly limestones containing a rich ammonite fauna, nannofossils, and benthic foraminifers-all these provide an accurate biostratigraphic control. The lower and upper Toarcian boundaries are well represented in some of these sections and therefore represent optimal sites to link the carbon-isotope curves to ammonite zones, and to nannofossil events. delta C-13 values of bulk carbonates from the different localities of the Subbetic basin have similar variations from the uppermost Pliensbachian to the lowermost Aalenian, suggesting changes in the original DIC carbon isotope composition along the Hispanic corridor. The transition from Pliensbachian to Toarcian is marked by increasing delta C-13 values from similar to 12 to 2.0 parts per thousand, interrupted in the Serpentinum Zone by a negative shift concomitant with the Toarcian oceanic anoxic event (T-OAE), with the major ammonite extinction event of the Toarcian, and an important turnover of calcareous nannoplankton. The negative shift observed in the Serpentinum Zone confirms the global perturbation of the carbon cycling documented along the Tethys and the palaeo-Pacific in organic material and in marine carbonates. However, the amplitude of the negative excursion (similar to - 1.5 parts per thousand) is not compatible with an isotopic homogeneous seawater DIC and/or CO2 atmospheric reservoirs. The interval from the middle to the top of the Toarcian delta C-13 shows relatively constant values, minor ammonite turnovers, and is associated with increasing diversity of calcareous nannoplankton. (c) 2012 Elsevier B.V. All rights reserved.

The risk of subsequent cancer and arterial cardiovascular events in patients with superficial vein thrombosis in the legs.

Although it has been clearly demonstrated that venous thromboembolism is associated with an increased risk of subsequent overt cancer and arterial cardiovascular events in comparison with control populations, whether this association also applies to patients with isolated (ie, without concomitant involvement of the deep vein system) superficial vein thrombosis (SVT) in the legs is unknown. In 737 consecutive patients with isolated SVT not involving the sapheno-femoral junction, we conducted a retrospective investigation to assess the rate of cancer and that of arterial cardiovascular events occurring during follow-up. The event rates were compared with those occurring in 1438 controls having comparable characteristics. Both cases and controls were followed-up for an average period of 26 ± 8 months (range, 3-45). Malignancy was diagnosed in 26 cases (3.5%) and 56 controls (3.9%), leading to a hazard ratio of 0.86 (95% confidence interval, 0.55%-1.35%). Arterial cardiovascular events occurred in 32 cases (4.3%) and 63 controls (4.4%), leading to a hazard ratio of 0.97 (95% confidence interval, 0.63%-1.50%). We conclude that the occurrence of isolated SVT in the legs does not place patients at an increased risk of malignancies or arterial cardiovascular events. Whether this conclusion also applies to patients whose thrombosis involves the sapheno-femoral junction remains to be demonstrated.

Causes for the underreporting of adverse drug events by health professionals: a systematic review

Objective: Identifying the main causes for underreporting of Adverse Drug Reaction (ADR) by health professionals. Method: A systematic review carried out in the following databases: LILACS, PAHO, SciELO, EMBASE and PubMed in the period between 1992 and 2012. Descriptors were used in the search for articles, and the identified causes of underreporting were analyzed according to the classification of Inman. Results: In total, were identified 149 articles, among which 29 were selected. Most studies were carried out in hospitals (24/29) for physicians (22/29), and pharmacists (10/29). The main causes related to underreporting were ignorance (24/29), insecurity (24/29) and indifference (23/29). Conclusion: The data show the eighth sin in underreporting, which is the lack of training in pharmacovigilance. Therefore, continuing education can increase adherence of professionals to the service and improve knowledge and communication of risks due to drug use.


Fábrica de "Metalheads" : arquetipo de la representación del aficionado al Heavy Metal en España

El presente estudio expone los resultados del tema desarrollado en la línea investigativa de culturas urbanas contemporáneas y procesos comunicacionales. Trabajo que describe (construye) el arquetipo del aficionado al heavy metal en España.

Subsampling the mean of heavy-tailed dependent observations

We establish the validity of subsampling confidence intervals for themean of a dependent series with heavy-tailed marginal distributions.Using point process theory, we study both linear and nonlinear GARCH-liketime series models. We propose a data-dependent method for the optimalblock size selection and investigate its performance by means of asimulation study.

Snapshot Events, 2005

This monthly report from the Iowa Department of Transportation is about the water quality management of Iowa's rivers, streams and lakes.