Inactivation of the putative suppressor gene DOK1 by promoter hypermethylation in primary human cancers

The DOK1 gene is a putative tumour suppressor gene located on the human chromosome 2p13 which is frequently rearranged in leukaemia and other human tumours. We previously reported that the DOK1 gene can be mutated and its expression down-regulated in human malignancies. However, the mechanism underlying DOK1 silencing remains largely unknown. We show here that unscheduled silencing of DOK1 expression through aberrant hypermethylation is a frequent event in a variety of human malignancies. DOK1 was found to be silenced in nine head and neck cancer (HNC) cell lines studied and DOK1 CpG hypermethylation correlated with loss of gene expression in these cells. DOK1 expression could be restored via demethylating treatment using 5-aza-2'deoxycytidine. In addition, transduction of cancer cell lines with DOK1 impaired their proliferation, consistent with the critical role of epigenetic silencing of DOK1 in the development and maintenance of malignant cells. We further observed that DOK1 hypermethylation occurs frequently in a variety of primary human neoplasm including solid tumours (93% in HNC, 81% in lung cancer) and haematopoietic malignancy (64% in Burkitt's lymphoma). Control blood samples and exfoliated mouth epithelial cells from healthy individuals showed a low level of DOK1 methylation, suggesting that DOK1 hypermethylation is a tumour specific event. Finally, an inverse correlation was observed between the level of DOK1 gene methylation and its expression in tumour and adjacent non tumour tissues. Thus, hypermethylation of DOK1 is a potentially critical event in human carcinogenesis, and may be a potential cancer biomarker and an attractive target for epigenetic-based therapy.

Repression of osteocyte Wnt/beta-catenin signaling is an early event in the progression of renal osteodystrophy

Chronic kidney diseasemineral bone disorder (CKD-MBD) is defined by abnormalities in mineral and hormone metabolism, bone histomorphometric changes, and/or the presence of soft-tissue calcification. Emerging evidence suggests that features of CKD-MBD may occur early in disease progression and are associated with changes in osteocyte function. To identify early changes in bone, we utilized the jck mouse, a genetic model of polycystic kidney disease that exhibits progressive renal disease. At 6 weeks of age, jck mice have normal renal function and no evidence of bone disease but exhibit continual decline in renal function and death by 20 weeks of age, when approximately 40% to 60% of them have vascular calcification. Temporal changes in serum parameters were identified in jck relative to wild-type mice from 6 through 18 weeks of age and were subsequently shown to largely mirror serum changes commonly associated with clinical CKD-MBD. Bone histomorphometry revealed progressive changes associated with increased osteoclast activity and elevated bone formation relative to wild-type mice. To capture the early molecular and cellular events in the progression of CKD-MBD we examined cell-specific pathways associated with bone remodeling at the protein and/or gene expression level. Importantly, a steady increase in the number of cells expressing phosphor-Ser33/37-beta-catenin was observed both in mouse and human bones. Overall repression of Wnt/beta-catenin signaling within osteocytes occurred in conjunction with increased expression of Wnt antagonists (SOST and sFRP4) and genes associated with osteoclast activity, including receptor activator of NF-?B ligand (RANKL). The resulting increase in the RANKL/osteoprotegerin (OPG) ratio correlated with increased osteoclast activity. In late-stage disease, an apparent repression of genes associated with osteoblast function was observed. These data confirm that jck mice develop progressive biochemical changes in CKD-MBD and suggest that repression of the Wnt/beta-catenin pathway is involved in the pathogenesis of renal osteodystrophy. (C) 2012 American Society for Bone and Mineral Research.

Multiple Quaternary Refugia in the Eastern Guiana Shield Revealed by Comparative Phylogeography of 12 Frog Species

The Guiana Shield (GS) is one of the most pristine regions of Amazonia and biologically one of the richest areas on Earth. How and when this massive diversity arose remains the subject of considerable debate. The prevailing hypothesis of Quaternary glacial refugia suggests that a part of the eastern GS, among other areas in Amazonia, served as stable forested refugia during periods of aridity. However, the recently proposed disturbance-vicariance hypothesis proposes that fluctuations in temperature on orbital timescales, with some associated aridity, have driven Neotropical diversification. The expectations of the temporal and spatial organization of biodiversity differ between these two hypotheses. Here, we compare the genetic structure of 12 leaf-litter inhabiting frog species from the GS lowlands using a combination of mitochondrial and nuclear sequences in an integrative analytical approach that includes phylogenetic reconstructions, molecular dating, and Geographic Information System methods. This comparative and integrated approach overcomes the well-known limitations of phylogeographic inference based on single species and single loci. All of the focal species exhibit distinct phylogeographic patterns highlighting taxon-specific historical distributions, ecological tolerances to climatic disturbance, and dispersal abilities. Nevertheless, all but one species exhibit a history of fragmentation/isolation within the eastern GS during the Quaternary with spatial and temporal concordance among species. The signature of isolation in northern French Guiana (FG) during the early Pleistocene is particularly clear. Approximate Bayesian Computation supports the synchrony of the divergence between northern FG and other GS lineages. Substructure observed throughout the GS suggests further Quaternary fragmentation and a role for rivers. Our findings support fragmentation of moist tropical forest in the eastern GS during this period when the refuge hypothesis would have the region serving as a contiguous wet-forest refuge.

Iodide treatment acutely increases pendrin (SLC26A4) mRNA expression in the rat thyroid and the PCCl3 thyroid cell line by transcriptional mechanisms

Iodine is a critical element involved in thyroid hormone synthesis. Its efflux into the follicular lumen is thought to occur, in part, through pendrin at the apical membrane of thyrocytes. This study attempted to investigate whether iodide administration affects SLC26A4 mRNA expression in rat thyroid and in PCCl3 cells. Rats and cells were treated or not with Nal from 30 min up to 48 h. One group was concomitantly treated with sodium perchlorate. SLC26A4 mRNA expression was also investigated in PCCl3 cells treated with actinomycin D prior to Nal treatment. Iodide administration significantly increased SLC26A4 mRNA content in both models. The simultaneous administration of Nal and perchlorate, as well as the treatment of PCCl3 cells with actinomycin D prevented this effect, indicating that intracellular iodide is essential for this event, which appears to be triggered by transcriptional mechanisms. These data show that intracellular iodide rapidly upregulates SLC26A4 mRNA expression. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Inclusive charged hadron elliptic flow in Au+Au collisions at root s(NN)=7.7-39 GeV

A systematic study is presented for centrality, transverse momentum (p(T)), and pseudorapidity (eta) dependence of the inclusive charged hadron elliptic flow (v(2)) at midrapidity (vertical bar eta vertical bar < 1.0) in Au + Au collisions at root s(NN) = 7.7, 11.5, 19.6, 27, and 39 GeV. The results obtained with different methods, including correlations with the event plane reconstructed in a region separated by a large pseudorapidity gap and four-particle cumulants (v(2){4}), are presented to investigate nonflow correlations and v(2) fluctuations. We observe that the difference between v(2){2} and v(2){4} is smaller at the lower collision energies. Values of v(2), scaled by the initial coordinate space eccentricity, v(2)/epsilon, as a function of p(T) are larger in more central collisions, suggesting stronger collective flow develops in more central collisions, similar to the results at higher collision energies. These results are compared to measurements at higher energies at the Relativistic Heavy Ion Collider (root s(NN) = 62.4 and 200 GeV) and at the Large Hadron Collider (Pb + Pb collisions at root s(NN) = 2.76 TeV). The v(2)(pT) values for fixed pT rise with increasing collision energy within the pT range studied (<2 GeV/c). A comparison to viscous hydrodynamic simulations is made to potentially help understand the energy dependence of v(2)(pT). We also compare the v(2) results to UrQMD and AMPT transport model calculations, and physics implications on the dominance of partonic versus hadronic phases in the system created at beam energy scan energies are discussed.

Energy and system-size dependence of two- and four-particle nu(2) measurements in heavy-ion collisions at root S-NN=62.4 and 200 GeV and their implications on flow fluctuations and nonflow

We present STAR measurements of azimuthal anisotropy by means of the two- and four-particle cumulants nu(2) (nu(2){2} and nu(2){4}) for Au + Au and Cu + Cu collisions at center-of-mass energies root S-NN = 62.4 and 200 GeV. The difference between nu(2){2}(2) and nu(2){4}(2) is related to nu(2) fluctuations (sigma(nu 2)) and nonflow (delta(2)). We present an upper limit to sigma(nu 2)/nu 2. Following the assumption that eccentricity fluctuations sigma(epsilon) dominate nu(2) fluctuations nu(2)/sigma nu(2) approximate to epsilon/sigma epsilon we deduce the nonflow implied for several models of eccentricity fluctuations that would be required for consistency with nu(2){2} and nu(2){4}. We also present results on the ratio of nu(2) to eccentricity.

Geochemical evidence of the 8.2 ka event and other Holocene environmental changes recorded in paleolagoon sediments, southeastern Brazil

The paleoclimatic record of Jureia Paleolagoon, coastal southeastem Brazil, includes cyclic and gradual changes with different intensities and frequencies through geological time, and it is controlled by astronomical, geophysical, and geological phenomena. These variations are not due to one single cause, but they result from the interaction of several factors, which act at different temporal and spatial scales. Here, we describe paleoenvironmental evidence regarding climatic and sea level changes from the last 9400 cal yr BP at the Jureia Paleolagoon - one of the main groups of protected South Atlantic ecosystems. Geochemical evidences were used to identify anomalies from multi-proxy analyses of a paleolagoon sediment core. The anomalies of centennial scale were correlated to climate and transgression-regression cycles from the Holocene period. Decadal scale anomalous oscillations in the Quaternary paleolagoon sediments occur between 9400 and 7500 cal yr BP, correlated with long- and short-term natural events, which generated high sedimentation rates, mainly between 8385 and 8375 cal yr BP (10 cm/yr). Our results suggest that a modem-day short-duration North Atlantic climatic event, such as the 82 ka event, could affect the environmental equilibrium in South America and intensify the South American Summer Monsoon. (C) 2011 University of Washington. Published by Elsevier Inc. All rights reserved.

A Robust Structural PGN Model for Control of Cell-Cycle Progression Stabilized by Negative Feedbacks

The cell division cycle comprises a sequence of phenomena controlled by a stable and robust genetic network. We applied a probabilistic genetic network (PGN) to construct a hypothetical model with a dynamical behavior displaying the degree of robustness typical of the biological cell cycle. The structure of our PGN model was inspired in well-established biological facts such as the existence of integrator subsystems, negative and positive feedback loops, and redundant signaling pathways. Our model represents genes interactions as stochastic processes and presents strong robustness in the presence of moderate noise and parameters fluctuations. A recently published deterministic yeast cell-cycle model does not perform as well as our PGN model, even upon moderate noise conditions. In addition, self stimulatory mechanisms can give our PGN model the possibility of having a pacemaker activity similar to the observed in the oscillatory embryonic cell cycle.

Functional clustering of time series gene expression data by Granger causality

Background: A common approach for time series gene expression data analysis includes the clustering of genes with similar expression patterns throughout time. Clustered gene expression profiles point to the joint contribution of groups of genes to a particular cellular process. However, since genes belong to intricate networks, other features, besides comparable expression patterns, should provide additional information for the identification of functionally similar genes. Results: In this study we perform gene clustering through the identification of Granger causality between and within sets of time series gene expression data. Granger causality is based on the idea that the cause of an event cannot come after its consequence. Conclusions: This kind of analysis can be used as a complementary approach for functional clustering, wherein genes would be clustered not solely based on their expression similarity but on their topological proximity built according to the intensity of Granger causality among them.

Fe valence fluctuations and magnetoelastic coupling in Pb-based multiferroics perovskites

Lead-based multiferroics perovskites with nominal compositions Pb(Fe1/2Nb1/2)O3 and Pb(Fe2/3W1/3)O3 were synthesized following a two-stage method. Magnetic proprieties were investigated and correlated to anelastic proprieties, measured by the conventional pulse-echo method. The discussions are focused in the region around 250 K, where magnetoelectroelastic instabilities have been observed. X-ray absorption nearedge structure (XANES) study further indicates that the edge position varies with temperature revealing a fluctuation on the valence of iron ions with the temperature, which can be related to a variation in anelastic and magnetic proprieties.

Expression of the repressor element-1 silencing transcription factor (REST) is regulated by IGF-I and PKC in human neuroblastoma cells

The repressor element 1-silencing transcription factor (REST) was first identified as a protein that binds to a 21-bp DNA sequence element (known as repressor element 1 (RE1)) resulting in transcriptional repression of the neural-specific genes [Chong et al., 1995; Schoenherr and Anderson, 1995]. The original proposed role for REST was that of a factor responsible for restricting neuronal gene expression to the nervous system by silencing expression of these genes in non-neuronal cells. Although it was initially thought to repress neuronal genes in non-neuronal cells, the role of REST is complex and tissue dependent. In this study I investigated any role played by REST in the induction and patterning of differentiation of SH-SY5Y human neuroblastoma cells exposed to IGF-I. and phorbol 12- myristate 13-acetate (PMA) To down-regulate REST expression we developed an antisense (AS) strategy based on the use of phosphorothioate oligonucleotides (ODNs). In order to evaluate REST mRNA levels, we developed a real-time PCR technique and REST protein levels were evaluated by western blotting. Results showed that nuclear REST is increased in SH-SY5Y neuroblastoma cells cultured in SFM and exposed to IGF-I for 2-days and it then declines in 5-day-treated cells concomitant with a progressive neurite extension. Also the phorbol ester PMA was able to increase nuclear REST levels after 3-days treatment concomitant to neuronal differentiation of neuroblastoma cells, whereas, at later stages, it is down-regulated. Supporting these data, the exposure to PKC inhibitors (GF10923X and Gö6976) and PMA (16nM) reverted the effects observed with PMA alone. REST levels were related to morphological differentiation, expression of growth coneassociated protein 43 (GAP-43; a gene not regulated by REST) and of synapsin I and βIII tubulin (genes regulated by REST), proteins involved in the early stage of neuronal development. We observed that differentiation of SH-SY5Y cells by IGF-I and PMA was accompanied by a significant increase of these neuronal markers, an effect that was concomitant with REST decrease. In order to relate the decreased REST expression with a progressive neurite extension, I investigated any possible involvement of the ubiquitin–proteasome system (UPS), a multienzymatic pathway which degrades polyubiquinated soluble cytoplasmic proteins [Pickart and Cohen, 2004]. For this purpose, SH-SY5Y cells are concomitantly exposed to PMA and the proteasome inhibitor MG132. In SH-SY5Y exposed to PMA and MG 132, we observed an inverse pattern of expression of synapsin I and β- tubulin III, two neuronal differentiation markers regulated by REST. Their cytoplasmic levels are reduced when compared to cells exposed to PMA alone, as a consequence of the increase of REST expression by proteasome inhibitor. The majority of proteasome substrates identified to date are marked for degradation by polyubiquitinylation; however, exceptions to this principle, are well documented [Hoyt and Coffino, 2004]. Interestingly, REST degradation seems to be completely ubiquitin-independent. The expression pattern of REST could be consistent with the theory that, during early neuronal differentiation induced by IGF-I and PKC, it may help to repress the expression of several genes not yet required by the differentiation program and then it declines later. Interestingly, the observation that REST expression is progressively reduced in parallel with cell proliferation seems to indicate that the role of this transcription factor could also be related to cell survival or to counteract apotosis events [Lawinger et al., 2000] although, as shown by AS-ODN experiments, it does not seem to be directly involved in cell proliferation. Therefore, the decline of REST expression is a comparatively later event during maturation of neuroroblasts in vitro. Thus, we propose that REST is regulated by growth factors, like IGF-I, and PKC activators in a time-dependent manner: it is elevated during early steps of neural induction and could contribute to down-regulate genes not yet required by the differentiation program while it declines later for the acquisition of neural phenotypes, concomitantly with a progressive neurite extension. This later decline is regulated by the proteasome system activation in an ubiquitin-indipendent way and adds more evidences to the hypothesis that REST down-regulation contributes to differentiation and arrest of proliferation of neuroblastoma cells. Finally, the glycosylation pattern of the REST protein was analysed, moving from the observation that the molecular weight calculated on REST sequence is about 116 kDa but using western blotting this transcription factor appears to have distinct apparent molecular weight (see Table 1.1): this difference could be explained by post-translational modifications of the proteins, like glycosylation. In fact recently, several studies underlined the importance of O-glycosylation in modulating transcriptional silencing, protein phosphorylation, protein degradation by proteasome and protein–protein interactions [Julenius et al., 2005; Zachara and Hart, 2006]. Deglycosilating analysis showed that REST protein in SH-SY5Y and HEK293 cells is Oglycosylated and not N-glycosylated. Moreover, using several combination of deglycosilating enzymes it is possible to hypothesize the presence of Gal-β(1-3)-GalNAc residues on the endogenous REST, while β(1-4)-linked galactose residues may be present on recombinant REST protein expressed in HEK293 cells. However, the O-glycosylation process produces an immense multiplicity of chemical structures and monosaccharides must be sequentially hydrolyzed by a series of exoglycosidase. Further experiments are needed to characterize all the post-translational modification of the transcription factor REST.

p300/CBP tyrosine phosphorylation in response to dna damage activated by c-Abl tyrosine kinase

The nuclear signaling that is triggered in response to DNA damage entails the recruitment and assembly of repair proteins and the induction of genes involved in the activation of cell cycle checkpoint, apoptosis or senescence. The extensive changes in chromatin structure underlying these processes suggest that chromatin-modifying enzymes could be relevant targets of DNA damage-activated signaling. The acetyltransferases p300 and CBP participate in DNA damage-activated responses, including local histone hyperacetylation, cell cycle regulation, and co-activation of DNA damage activated proteins, such as p53, p73 and BRCA1. However, the link between DNA damage and p300/CBP activation has not been identified.We have detected p300 tyrosine phosphorylation in response to DNA damage. We show that the DNA damage-activated cAbl tyrosine kinase enters the nuclei of cells exposed to genotoxic agents and phosphorylates p300 on a tyrosine residue within the bromodomain that is conserved in p300, CBP and many other bromodomain-containing proteins. Antibodies against tyrosine phosphorylated p300/CBP show a DNA damage-inducible nuclear staining, suggesting that p300 tyrosine phosphorylation is an event linking DNA damage and chromatin modifications.

Slope stability analysis by multi-temporal DEMs and 3D modelling: The 2002 and 2007 Stromboli landslide events

Natural hazard related to the volcanic activity represents a potential risk factor, particularly in the vicinity of human settlements. Besides to the risk related to the explosive and effusive activity, the instability of volcanic edifices may develop into large landslides often catastrophically destructive, as shown by the collapse of the northern flank of Mount St. Helens in 1980. A combined approach was applied to analyse slope failures that occurred at Stromboli volcano. SdF slope stability was evaluated by using high-resolution multi-temporal DTMMs and performing limit equilibrium stability analyses. High-resolution topographical data collected with remote sensing techniques and three-dimensional slope stability analysis play a key role in understanding instability mechanism and the related risks. Analyses carried out on the 2002–2003 and 2007 Stromboli eruptions, starting from high-resolution data acquired through airborne remote sensing surveys, permitted the estimation of the lava volumes emplaced on the SdF slope and contributed to the investigation of the link between magma emission and slope instabilities. Limit Equilibrium analyses were performed on the 2001 and 2007 3D models, in order to simulate the slope behavior before 2002-2003 landslide event and after the 2007 eruption. Stability analyses were conducted to understand the mechanisms that controlled the slope deformations which occurred shortly after the 2007 eruption onset, involving the upper part of slope. Limit equilibrium analyses applied to both cases yielded results which are congruent with observations and monitoring data. The results presented in this work undoubtedly indicate that hazard assessment for the island of Stromboli should take into account the fact that a new magma intrusion could lead to further destabilisation of the slope, which may be more significant than the one recently observed because it will affect an already disarranged deposit and fractured and loosened crater area. The two-pronged approach based on the analysis of 3D multi-temporal mapping datasets and on the application of LE methods contributed to better understanding volcano flank behaviour and to be prepared to undertake actions aimed at risk mitigation.

Mechanics and dynamics of liposomes and cells studied by QCM and ECIS

The present thesis introduces a novel sensitive technique based on TSM resonators that provides quantitative information about the dynamic properties of biological cells and artificial lipid systems. In order to support and complement results obtained by this method supplementary measurements based on ECIS technique were carried out. The first part (chapters 3 and 4) deals with artificial lipid systems. In chapter 3 ECIS measurements were used to monitor the adsorption of giant unilamellar vesicles as well as their thermal fluctuations. From dynamic Monte Carlo Simulations the rate constant of vesicle adsorption was determined. Furthermore, analysis of fluctuation measurements reveals Brownian motion reflecting membrane undulations of the adherent liposomes. In chapter 4 QCM-based fluctuation measurements were applied to quantify nanoscopically small deformations of giant unilamellar vesicles with an external electrical field applied simultaneously. The response of liposomes to an external voltage with shape changes was monitored as a function of cholesterol content and adhesion force. In the second part (chapters 5 - 8) attention was given to cell motility. It was shown for the first time, that QCM can be applied to monitor the dynamics of living adherent cells in real time. QCM turned out to be a highly sensitive tool to detect the vertical motility of adherent cells with a time resolution in the millisecond regime. The response of cells to environmental changes such as temperature or osmotic stress could be quantified. Furthermore, the impact of cytochalasin D (inhibits actin polymerization) and taxol (facilitate polymerization of microtubules) as well as nocodazole (depolymerizes microtubules) on the dynamic properties of cells was scrutinized. Each drug provoked a significant reduction of the monitored cell shape fluctuations as expected from their biochemical potential. However, not only the abolition of fluctuations was observed but also an increase of motility due to integrin-induced transmembrane signals. These signals were activated by peptides containing the RGD sequence, which is known to be an integrin recognition motif. Ultimately, two pancreatic carcinoma cell lines, derived from the same original tumor, but known to possess different metastatic potential were studied. Different dynamic behavior of the two cell lines was observed which was attributed to cell-cell as well as cell-substrate interactions rather than motility. Thus one may envision that it might be possible to characterize the motility of different cell types as a function of many variables by this new highly sensitive technique based on TSM resonators. Finally the origin of the broad cell resonance was investigated. Improvement of the time resolution reveals the "real" frequency of cell shape fluctuations. Several broad resonances around 3-5 Hz, 15-17 Hz and 25-29 Hz were observed and that could unequivocally be assigned to biological activity of living cells. However, the kind of biological process that provokes this synchronized collective and periodic behavior of the cells remains to be elucidated.

Thermal fluctuations and boundary conditions in the lattice Boltzmann method

The lattice Boltzmann method is a popular approach for simulating hydrodynamic interactions in soft matter and complex fluids. The solvent is represented on a discrete lattice whose nodes are populated by particle distributions that propagate on the discrete links between the nodes and undergo local collisions. On large length and time scales, the microdynamics leads to a hydrodynamic flow field that satisfies the Navier-Stokes equation. In this thesis, several extensions to the lattice Boltzmann method are developed. In complex fluids, for example suspensions, Brownian motion of the solutes is of paramount importance. However, it can not be simulated with the original lattice Boltzmann method because the dynamics is completely deterministic. It is possible, though, to introduce thermal fluctuations in order to reproduce the equations of fluctuating hydrodynamics. In this work, a generalized lattice gas model is used to systematically derive the fluctuating lattice Boltzmann equation from statistical mechanics principles. The stochastic part of the dynamics is interpreted as a Monte Carlo process, which is then required to satisfy the condition of detailed balance. This leads to an expression for the thermal fluctuations which implies that it is essential to thermalize all degrees of freedom of the system, including the kinetic modes. The new formalism guarantees that the fluctuating lattice Boltzmann equation is simultaneously consistent with both fluctuating hydrodynamics and statistical mechanics. This establishes a foundation for future extensions, such as the treatment of multi-phase and thermal flows. An important range of applications for the lattice Boltzmann method is formed by microfluidics. Fostered by the "lab-on-a-chip" paradigm, there is an increasing need for computer simulations which are able to complement the achievements of theory and experiment. Microfluidic systems are characterized by a large surface-to-volume ratio and, therefore, boundary conditions are of special relevance. On the microscale, the standard no-slip boundary condition used in hydrodynamics has to be replaced by a slip boundary condition. In this work, a boundary condition for lattice Boltzmann is constructed that allows the slip length to be tuned by a single model parameter. Furthermore, a conceptually new approach for constructing boundary conditions is explored, where the reduced symmetry at the boundary is explicitly incorporated into the lattice model. The lattice Boltzmann method is systematically extended to the reduced symmetry model. In the case of a Poiseuille flow in a plane channel, it is shown that a special choice of the collision operator is required to reproduce the correct flow profile. This systematic approach sheds light on the consequences of the reduced symmetry at the boundary and leads to a deeper understanding of boundary conditions in the lattice Boltzmann method. This can help to develop improved boundary conditions that lead to more accurate simulation results.