Linear Aggregation In The Social Accounting Matrix Framework

In economic literature, information deficiencies and computational complexities have traditionally been solved through the aggregation of agents and institutions. In inputoutput modelling, researchers have been interested in the aggregation problem since the beginning of 1950s. Extending the conventional input-output aggregation approach to the social accounting matrix (SAM) models may help to identify the effects caused by the information problems and data deficiencies that usually appear in the SAM framework. This paper develops the theory of aggregation and applies it to the social accounting matrix model of multipliers. First, we define the concept of linear aggregation in a SAM database context. Second, we define the aggregated partitioned matrices of multipliers which are characteristic of the SAM approach. Third, we extend the analysis to other related concepts, such as aggregation bias and consistency in aggregation. Finally, we provide an illustrative example that shows the effects of aggregating a social accounting matrix model.

Complementarities between mathematical and computational modeling: the case of the repeated Prisoners' Dilemma

We study the properties of the well known Replicator Dynamics when applied to a finitely repeated version of the Prisoners' Dilemma game. We characterize the behavior of such dynamics under strongly simplifying assumptions (i.e. only 3 strategies are available) and show that the basin of attraction of defection shrinks as the number of repetitions increases. After discussing the difficulties involved in trying to relax the 'strongly simplifying assumptions' above, we approach the same model by means of simulations based on genetic algorithms. The resulting simulations describe a behavior of the system very close to the one predicted by the replicator dynamics without imposing any of the assumptions of the mathematical model. Our main conclusion is that mathematical and computational models are good complements for research in social sciences. Indeed, while computational models are extremely useful to extend the scope of the analysis to complex scenarios hard to analyze mathematically, formal models can be useful to verify and to explain the outcomes of computational models.

Line search multilevel optimization as computational methods for dense optical flow

We evaluate the performance of different optimization techniques developed in the context of optical flowcomputation with different variational models. In particular, based on truncated Newton methods (TN) that have been an effective approach for large-scale unconstrained optimization, we develop the use of efficient multilevel schemes for computing the optical flow. More precisely, we evaluate the performance of a standard unidirectional multilevel algorithm - called multiresolution optimization (MR/OPT), to a bidrectional multilevel algorithm - called full multigrid optimization (FMG/OPT). The FMG/OPT algorithm treats the coarse grid correction as an optimization search direction and eventually scales it using a line search. Experimental results on different image sequences using four models of optical flow computation show that the FMG/OPT algorithm outperforms both the TN and MR/OPT algorithms in terms of the computational work and the quality of the optical flow estimation.

Mutations at a single codon in Mad homology 2 domain of SMAD4 cause Myhre syndrome.

Myhre syndrome (MIM 139210) is a developmental disorder characterized by short stature, short hands and feet, facial dysmorphism, muscular hypertrophy, deafness and cognitive delay. Using exome sequencing of individuals with Myhre syndrome, we identified SMAD4 as a candidate gene that contributes to this syndrome on the basis of its pivotal role in the bone morphogenetic pathway (BMP) and transforming growth factor (TGF)-β signaling. We identified three distinct heterozygous missense SMAD4 mutations affecting the codon for Ile500 in 11 individuals with Myhre syndrome. All three mutations are located in the region of SMAD4 encoding the Mad homology 2 (MH2) domain near the site of monoubiquitination at Lys519, and we found a defect in SMAD4 ubiquitination in fibroblasts from affected individuals. We also observed decreased expression of downstream TGF-β target genes, supporting the idea of impaired TGF-β-mediated transcriptional control in individuals with Myhre syndrome.

Different approaches to modelling the cost-effectiveness of schistosomiasis control

This paper reviews three different approaches to modelling the cost-effectiveness of schistosomiasis control. Although these approaches vary in their assessment of costs, the major focus of the paper is on the evaluation of effectiveness. The first model presented is a static economic model which assesses effectiveness in terms of the proportion of cases cured. This model is important in highlighting that the optimal choice of chemotherapy regime depends critically on the level of budget constraint, the unit costs of screening and treatment, the rates of compliance with screening and chemotherapy and the prevalence of infection. The limitations of this approach is that it models the cost-effectiveness of only one cycle of treatment, and effectiveness reflects only the immediate impact of treatment. The second model presented is a prevalence-based dynamic model which links prevalence rates from one year to the next, and assesses effectiveness as the proportion of cases prevented. This model was important as it introduced the concept of measuring the long-term impact of control by using a transmission model which can assess reduction in infection through time, but is limited to assessing the impact only on the prevalence of infection. The third approach presented is a theoretical framework which describes the dynamic relationships between infection and morbidity, and which assesses effectiveness in terms of case-years prevented of infection and morbidity. The use of this model in assessing the cost-effectiveness of age-targeted treatment in controlling Schistosoma mansoni is explored in detail, with respect to varying frequencies of treatment and the interaction between drug price and drug efficacy.

Protective immunity induced in mice by F8.1 and F8.2 antigens purified from Schistosoma mansoni eggs

Schistosoma mansoni soluble egg antigens (SEA) were fractionated by isoelectric focusing, resulting in 20 components, characterized by pH, absorbance and protein concentration. The higher absorbance fractions were submitted to electrophoresis, and fraction 8 (F8) presented a specific pattern of bands on its isoelectric point. Protein 3 was observed only on F8, and so, it was utilized to rabbit immunization, in order to evaluate its capacity of inducing protective immunity. IgG antibodies from rabbit anti-F8 serum were coupled to Sepharose, and used to obtain the specific antigen by affinity chromatography. This antigen, submitted to electrophoresis, presented two proteic bands (F8.1 and F8.2), which were transferred to nitrocellulose membrane (PVDF) and sequenciated. The homology of F8.2 to known proteins was determined using the Basic Local Alignment Search Tool program (BLASTp). Significant homologies were obtained for the rabbit cytosolic Ca2+ uptake inhibitor, and for the bird a1-proteinase inhibitor. Immunization of mice with F8.1 and F8.2, in the presence of Corynebacterium parvum and Al(OH)3 as adjuvant, induced a significant protection degree against challenge infection, as observed by the decrease on worm burden recovered from portal system.

Population pharmacokinetic modelling and evaluation of different dosage regimens for darunavir and ritonavir in HIV-infected individuals.

OBJECTIVES: Darunavir is a protease inhibitor that is administered with low-dose ritonavir to enhance its bioavailability. It is prescribed at standard dosage regimens of 600/100 mg twice daily in treatment-experienced patients and 800/100 mg once daily in naive patients. A population pharmacokinetic approach was used to characterize the pharmacokinetics of both drugs and their interaction in a cohort of unselected patients and to compare darunavir exposure expected under alternative dosage regimens. METHODS: The study population included 105 HIV-infected individuals who provided darunavir and ritonavir plasma concentrations. Firstly, a population pharmacokinetic analysis for darunavir and ritonavir was conducted, with inclusion of patients' demographic, clinical and genetic characteristics as potential covariates (NONMEM(®)). Then, the interaction between darunavir and ritonavir was studied while incorporating levels of both drugs into different inhibitory models. Finally, model-based simulations were performed to compare trough concentrations (Cmin) between the recommended dosage regimen and alternative combinations of darunavir and ritonavir. RESULTS: A one-compartment model with first-order absorption adequately characterized darunavir and ritonavir pharmacokinetics. The between-subject variability in both compounds was important [coefficient of variation (CV%) 34% and 47% for darunavir and ritonavir clearance, respectively]. Lopinavir and ritonavir exposure (AUC) affected darunavir clearance, while body weight and darunavir AUC influenced ritonavir elimination. None of the tested genetic variants showed any influence on darunavir or ritonavir pharmacokinetics. The simulations predicted darunavir Cmin much higher than the IC50 thresholds for wild-type and protease inhibitor-resistant HIV-1 strains (55 and 550 ng/mL, respectively) under standard dosing in >98% of experienced and naive patients. Alternative regimens of darunavir/ritonavir 1200/100 or 1200/200 mg once daily also had predicted adequate Cmin (>550 ng/mL) in 84% and 93% of patients, respectively. Reduction of darunavir/ritonavir dosage to 600/50 mg twice daily led to a 23% reduction in average Cmin, still with only 3.8% of patients having concentrations below the IC50 for resistant strains. CONCLUSIONS: The important variability in darunavir and ritonavir pharmacokinetics is poorly explained by clinical covariates and genetic influences. In experienced patients, treatment simplification strategies guided by drug level measurements and adherence monitoring could be proposed.

Comparative Analysis by Polymerase Chain Reaction Amplified Minicircles of Kinetoplast DNA of a Stable Strain of Trypanosoma cruzi from São Felipe, Bahia, its Clones and Subclones: Possibility of Predominance of a Principal Clone in this Area

Molecular characterization of one stable strain of Trypanosoma cruzi, the 21 SF, representative of the pattern of strains isolated from the endemic area of São Felipe, State of Bahia, Brazil, maintained for 15 years in laboratory by serial passages in mice and classified as biodeme Type II and zymodeme 2 has been investigated. The kinetoplast DNA (kDNA) of parental strain, 5 clones and 14 subclones were analyzed. Schizodeme was established by comparative study of the fragments obtained from digestion of the 330-bp fragments amplified by polymerase chain reaction (PCR) from the variable regions of the minicicles, and digested by restriction endonucleases Rsa I and Hinf I. Our results show a high percentual of similarity between the restriction fragment lenght polymorphism (RFLP) for the parental strain and its clones and among these individual clones and their subclones at a level of 80 to 100%.This homology indicates a predominance of the same "principal clone" in the 21SF strain and confirms the homogeneity previously observed at biological and isozymic analysis. These results suggest the possibility that the T. cruzi strains with similar biological and isoenzymic patterns, circulating in this endemic area, are representative of one dominant clone. The presence of "principal clones" could be responsible for a predominant tropism of the parasites for specific organs and tissues and this could contribute to the pattern of clinico-pathological manifestations of Chagas's disease in one geographical area.

Three-dimensional models of 14α-sterol demethylase (Cyp51A) from Aspergillus lentulus and Aspergillus fumigatus: an insight into differences in voriconazole interaction.

Aspergillus lentulus, an Aspergillus fumigatus sibling species, is increasingly reported in corticosteroid-treated patients. Its clinical significance is unknown, but the fact that A. lentulus shows reduced antifungal susceptibility, mainly to voriconazole, is of serious concern. Heterologous expression of cyp51A from A. fumigatus and A. lentulus was performed in Saccharomyces cerevisiae to assess differences in the interaction of Cyp51A with the azole drugs. The absence of endogenous ERG11 was efficiently complemented in S. cerevisiae by the expression of either Aspergillus cyp51A allele. There was a marked difference between azole minimum inhibitory concentration (MIC) values of the clones expressing each Aspergillus spp. cyp51A. Saccharomyces cerevisiae clones expressing A. lentulus alleles showed higher MICs to all of the azoles tested, supporting the hypothesis that the intrinsic azole resistance of A. lentulus could be associated with Cyp51A. Homology models of A. fumigatus and A. lentulus Cyp51A protein based on the crystal structure of Cyp51p from Mycobacterium tuberculosis in complex with fluconazole were almost identical owing to their mutual high sequence identity. Molecular dynamics (MD) was applied to both three-dimensional protein models to refine the homology modelling and to explore possible differences in the Cyp51A-voriconazole interaction. After 20ns of MD modelling, some critical differences were observed in the putative closed form adopted by the protein upon voriconazole binding. A closer study of the A. fumigatus and A. lentulus voriconazole putative binding site in Cyp51A suggested that some major differences in the protein's BC loop could differentially affect the lock-up of voriconazole, which in turn could correlate with their different azole susceptibility profiles.

Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Transcriptional Activation and Receptor Dimerization is Affected by Mutations in the NR2E3 Ligand-Binding Domain

Purpose: Mutations in the ligand-binding domain (LBD) of NR2E3 cause recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS), Goldmann-Favre syndrome (GFS) and clumped pigmentary retinal degeneration (CPRD). In addition to ligand binding, the LBD contains also essential amino acid sequences for the oligomerization of nuclear receptors. The aim of our studies is to characterize the impact of mutations in the LBD on receptor oligomerization and transcriptional activity of NR2E3. Methods: The different NR2E3 mutants were generated by QuickChange mutagenesis and analyzed in 293T-based transactivation studies and BRET2 (bioluminescence resonance electron transfer) assays. In silico homology modeling of mutant proteins was also performed using available crystallographic data of related nuclear receptors. Results: The mutants p.W234S, p.A256V, p.A256E, p.L263P, p.R309G, p.R311Q, p.R334G, p.L336P, p.L353V, p.R385P and p.M407K, all located in the LBD, showed impaired receptor dimerization at various degrees. Impaired repressor dimerization as assessed by BRET2 assays did not always correlate with impaired repressor function of NR2E3 as assessed by cell-based reporter assays. There were minor differences of transcriptional activity of mutant proteins on mouse S-opsin (opn1sw), mouse cone arrestin (arr3) and human cone arrestin, suggesting that the effect of LBD mutations was independent of the promoter context. Conclusions: Mutational analysis and homology modeling allowed the characterization of potential oligomerization interfaces of the NR2E3 LBD. Additionally, mutations in NR2E3 LBD may cause recessive retinal degenerations by different molecular mechanisms.

Characterization of human fibroleukin, a fibrinogen-like protein secreted by T lymphocytes.

We have recently cloned the human homologue of the murine pT49 cDNA (hpT49h), a transcript encoding a protein homologous to the beta- and gamma-chains of fibrinogen. Here, we report the identification of the hpT49h gene product using mAbs generated against a peptide corresponding to the carboxyl-terminal end of the deduced protein and a recombinant protein fragment expressed in Escherichia coli. mAbs 23A6, 7B12, and 3F4 specifically recognized a protein of 70 kDa in reducing SDS-PAGE in the culture supernatant of 293T cells transiently transfected with the full length hpT49h cDNA and freshly isolated PBMC. Under nonreducing conditions, the material migrated with a molecular mass of 250 to 300 kDa, indicating that the 70-kDa protein forms a disulfide bonded complex. Because of its homology with fibrinogen, we have termed this protein fibroleukin. Fibroleukin is spontaneously secreted in vitro by freshly isolated CD4+ and CD8+ T lymphocytes. RT-PCR analysis revealed preferential expression of fibroleukin mRNA in memory T lymphocytes (CD3+/CD45R0+) compared with naive T lymphocytes (CD3+/CD45RA+). Fibroleukin production by PBMC was rapidly lost in culture. Production could be partially maintained in the presence of IFN-gamma, while T lymphocyte activation had no effect. To demonstrate fibroleukin production in vivo, we analyzed colon mucosa by immunohistology. Fibroleukin staining was detected in the extracellular matrix of the T lymphocyte-rich upper portion of the lamina propria mucosa. While the exact function of fibroleukin remains to be defined, these data suggest that fibroleukin may play a role in physiologic lymphocyte functions at mucosal sites.

Mexican Trypanosoma cruzi stocks: analysis of minicircle kDNA homologies by cross-hybridization

Homologies of minicircle kDNA of 27 Mexican stocks were studied by cross-hybridization with four kDNA probes derived from three reference stocks belonging to groups Trypanosoma cruzi I (SO34 cl4 and Silvio) and T. cruzi II (MN) and one Mexican stock. High homologies were only observed with Silvio (six stocks) and Mexican probes (11 stocks). After 30 min exposure (low homology) additional stocks were recognized with SO34 cl4 (three stocks) and Silvio (six stocks) probes; with the Mexican probe only five stocks remained non-reactive. All the stocks were typed by isoenzyme (16 loci) and Mexican stocks belonged to T. cruzi I. Hybridization patterns were not strictly correlated with the observed clustering and cross-hybridization of kDNA minicircles is not available to distinct Mexican stocks.

LEDGF/p75 TATA-less promoter is driven by the transcription factor Sp1.

PSIP1 (PC4 and SFRS1 interacting protein 1) encodes two splice variants: lens epithelium-derived growth factor or p75 (LEDGF/p75) and p52. PSIP1 gene products were shown to be involved in transcriptional regulation, affecting a plethora of cellular processes, including cell proliferation, cell survival, and stress response. Furthermore, LEDGF/p75 has implications for various diseases and infections, including autoimmunity, leukemia, embryo development, psoriasis, and human immunodeficiency virus integration. Here, we reported the first characterization of the PSIP1 promoter. Using 5' RNA ligase-mediated rapid amplification of cDNA ends, we identified novel transcription start sites in different cell types. Using a luciferase reporter system, we identified regulatory elements controlling the expression of LEDGF/p75 and p52. These include (i) minimal promoters (-112/+59 and +609/+781) that drive the basal expression of LEDGF/p75 and of the shorter splice variant p52, respectively; (ii) a sequence (+319/+397) that may control the ratio of LEDGF/p75 expression to p52 expression; and (iii) a strong enhancer (-320/-207) implicated in the modulation of LEDGF/p75 transcriptional activity. Computational, biochemical, and genetic approaches enabled us to identify the transcription factor Sp1 as a key modulator of the PSIP1 promoter, controlling LEDGF/p75 transcription through two binding sites at -72/-64 and -46/-36. Overall, our results provide initial data concerning LEDGF/p75 promoter regulation, giving new insights to further understand its biological function and opening the door for new therapeutic strategies in which LEDGF/p75 is involved.

Desenvolupament de codis numèrics per a la resolució de fluxos turbulents en ordinadors paral·lels emprant discretitzacions tipus symmetry-preserving

Projecte de recerca elaborat a partir d’una estada a la University of Groningen, Holanda, entre 2007 i 2009. La simulació directa de la turbulència (DNS) és una eina clau dins de la mecànica de fluids computacional. Per una banda permet conèixer millor la física de la turbulència i per l'altra els resultats obtinguts són claus per el desenvolupament dels models de turbulència. No obstant, el DNS no és una tècnica vàlida per a la gran majoria d'aplicacions industrials degut al elevats costos computacionals. Per tant, és necessari cert grau de modelització de la turbulència. En aquest context, s'han introduïts importants millores basades en la modelització del terme convectiu (no lineal) emprant symmetry-preserving regularizations. En tracta de modificar adequadament el terme convectiu a fi de reduir la producció d'escales més i més petites (vortex-stretching) tot mantenint tots els invariants de les equacions originals. Fins ara, aquest models s'han emprat amb èxit per nombres de Rayleigh (Ra) relativament elevats. En aquest punt, disposar de resultats DNS per a configuracions més complexes i nombres de Ra més elevats és clau. En aquest contexte, s'han dut a terme simulacions DNS en el supercomputador MareNostrum d'una Differentially Heated Cavity amb Ra=1e11 i Pr=0.71 durant el primer any dels dos que consta el projecte. A més a més, s'ha adaptat el codi a fi de poder simular el fluxe al voltant d'un cub sobre una pared amb Re=10000. Aquestes simulacions DNS són les més grans fetes fins ara per aquestes configuracions i la seva correcta modelització és un gran repte degut la complexitat dels fluxes. Aquestes noves simulacions DNS estan aportant nous coneixements a la física de la turbulència i aportant resultats indispensables per al progrés de les modelitzacións tipus symmetry-preserving regularization.