Biohydrogen is a versatile energy carrier for the generation of electric energy from renewable sources. Hydrogenases can be used in enzymatic fuel cells to oxidize dihydrogen. The rate of electron transfer (ET) at the anodic side between the [NiFe]-hydrogenase enzyme distal iron-sulfur cluster and the electrode surface can be described by the Marcus equation. All parameters for the Marcus equation are accessible from DFT calculations. The distal cubane FeS-cluster has a three cysteine and one histidine coordination [Fe₄S₄](His)(Cys)₃ first ligation sphere. The reorganization energy (inner- and outer-sphere) is almost unchanged upon a histidine-to-cysteine substitution. Differences in rates of electron transfer between the wild-type enzyme and the all-cysteine mutant can be rationalized by a diminished electronic coupling between the donor and acceptor molecules in the [Fe₄S₄](Cys)₄ case. The fast and efficient electron transfer from the distal iron-sulfur cluster is realized by a fine-tuned protein environment which facilitates the flow of electrons. This study enables the design and control of electron transfer rates and pathways by protein engineering.

Purification of biologically-derived therapeutics is a major cost contributor to the production of this rapidly growing class of pharmaceuticals. Monoclonal antibodies comprise a large percentage of these products therefore new antibody purification tools are needed. Small peptides, as opposed to traditional antibody affinity ligands such as Protein A, may have advantages in stability and production costs. Multiple heptapeptides that demonstrate Fc binding behavior that have been identified from a combinatorial peptide library using M13 Phage Display are presented herein. Seven unique peptide sequences of diverse hydrophobicity and charge were identified. All seven peptides showed strong binding to the four major human IgG isotypes, human IgM, as well as binding to canine, rat, and mouse IgG. These seven peptides were also shown to bind human IgG4 from DMEM cell culture media with 5% FCS and 5 g/L ovalbumin present. These peptides may be useful as surface ligands for antibody detection and purification purposes. Molecular docking and classical molecular dynamics (MD) simulations were conducted to elucidate the mechanisms and energetics for the binding of these peptides to the Fc region. The binding site was found to be located between the two glycan chains inside the Fc fragment. Both hydrogen bonding and hydrophobic interactions were found to be crucial for the binding interactions. Excellent agreement for the binding strength was obtained between experimental results and simulations.

The attitude dynamics and orbit dynamics of large rigid space structures are coupled with each other under gravity gradient, which may affect the stability of the large rigid space structures. In this paper, the gravity gradient stability of a large rigid space structure is studied under both small and large disturbances. Based on the rigid body dynamics and orbit dynamics, an accurate dynamic model without any linearization of the large rigid space structure is established via the natural coordinate formulation (NCF), which is able to describe the large overall motions of the structures. By using the generalized-α algorithm, the gravity gradient stability of the large rigid space structure is simulated and analyzed via various examples, including the influence of large disturbance angles, the positions at the stabilization and unstabilization regions. Finally, the relationship between spinning stability and gravity gradient stability is also investigated via a large spinning space structure.

A U.S. government health agency conducted epidemiological studies to evaluate whether exposures to drinking water contaminated with volatile organic compounds at U.S. Marine Corps Base Camp Lejeune, North Carolina, were associated with increased health risks to children and adults. These health studies required knowledge of contaminant concentrations in drinking water—at monthly intervals—delivered to family housing, barracks, and other facilities within the study area. Because concentration data were limited or unavailable during much of the period of contamination (1950s–1985), the historical reconstruction process was used to quantify estimates of monthly mean contaminant-specific concentrations. This paper integrates many efforts, reports, and papers into a synthesis of the overall approach to, and results from, a drinking-water historical reconstruction study. Results show that at the Tarawa Terrace water treatment plant (WTP) reconstructed (simulated) tetrachloroethylene (PCE) concentrations reached a maximum monthly average value of 183 micrograms per liter (ug/L) compared to a one-time maximum measured value of 215 ug/L and exceeded the U.S. Environmental Protection Agency’s current maximum contaminant level (MCL) of 5 ug/L during the period November 1957–February 1987. At the Hadnot Point WTP, reconstructed trichloroethylene (TCE) concentrations reached a maximum monthly average value of 783 ug/L compared to a one-time maximum measured value of 1,400 ug/L during the period August 1953–December 1984. The Hadnot Point WTP also provided contaminated drinking water to the Holcomb Boulevard housing area continuously prior to June 1972, when the Holcomb Boulevard WTP came on line (maximum reconstructed TCE concentration of 32 ug/L) and intermittently during the period June 1972–February 1985 (maximum reconstructed TCE concentration of 66 ug/L). Applying the historical reconstruction process to quantify contaminant-specific monthly drinking-water concentrations is advantageous for epidemiological studies when compared to using the classical exposed versus unexposed approach.

Nuclear export of viral mRNAs, is an essential step in the HIV replication cycle. This role is played by a small regulatory protein of HIV-1 called Rev.The N-terminal region of Rev contains an arginine-rich sequence. The arginine-rich motif (ARM) is located between amino acids 38-50 and forms an alpha-helical secondary structure. Expression of the structural proteins of human immunodeficiency virus type 1 requires the direct interaction of multiple copies of the viral Rev protein with its highly structured RNA target sequence, the Rev Response element (RRE). The major viral proteins are not produced if this transport of RNA is stopped. Therefore, knowledge of Rev structure is essential for understanding of its cooperative binding to the RRE, for understanding the mechanism of HIV infection and for the development of antiviral drugs that interfere with Rev’s essential functions and for acknowledgment of good candidate drugs for treatment of AIDS. To understand how REV interact with RRE element of HIV-RNA and its formation of oligomeric complex it is better to characterize the domain wise structure of REV with regard in function of each domain. Due to lack of structural data on Rev no single compound is reported as inhibitor of REV expect antiviral drugs. Identification of a high-affinity RNA-binding site for the human immunodeficiency virus type 1 Rev Protein is much more important. The ARM is a highly specific sequence which allows for the multimerization of Rev and also binding of REV with RNA. Here we are first time exploring the structural characteristics of REV protein both in free form and in complex with RNA at domain function level especially explore the role of ARM motif in REV HIV-1 protein as RNA binding sites by molecular dynamics (MD) simulation and homology modeling studies. Results indicate that the arginine-rich motif (ARM) is crucial in stability of this complex. The residues ARG38, 39, 41, 43, 44, 48, 50, and ASN40 are most interacting with nucleobases of RRE in Crystal structure of Rev and Rev-response-element RNA complex. Our study plays a major role in elaboration of binding of RNA with REV and pave the way for further investigation for therapeutically agent for HIV.

The β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively probe identification mechanism of BACE1 on three inhibitors 60W, 954 and 60X. The analyses on MD trajectories indicate that the presence of three inhibitors influences structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward treatment of AD.

Boltzmann codes are essential tools for studying cosmology. Most of the codes are based on a seminal work by Ma and Bertschinger. We found that the formalism employed in those codes has at least three possible problems which need to be understood. i) The equation of motion for baryons are gauge incompatible, ii) they break the Bianchi identity, and iii) it is not clear from the equations of motion which physical system is considered. In this work we revisit the baryon physics and address all the above mentioned issues, based solely on conservation of stress-energy tensor, resulting in taking into account a correction that is usually neglected and that is numerically small. We also study the tight coupling approximation up to the second order without choosing any gauge. We implement the improved baryon equations in a Boltzmann code and investigate the change in the estimate of cosmological parameters by performing an MCMC analysis. While in this paper we study the Lambda-CDM model only, our baryon equations can be easily implemented in other models and various modified gravity theories.

The discrete element method (DEM) coupled with computational dynamics (CFD) has been considered one of the most sensitive ways of studying fluidized beds. This paper presents gas-solid fluidized bed simulations by use of CFD-DEM. The two-phase interaction is calculated in a way of decoupling. Within the effective neighborhood of the particles, the kernel approximation method is used to calculate the local porosity, thereby calculating the drag force on the particles. At the grid scale, the momentum exchange coefficient is calculated using a two fluid model based method, thereby calculating the source term of the phase interaction in the fluid control equation. The simulated form, size and motion process of the big bubble are all consistent with experimental observations. The simulated solid volume fraction and relative pressure averaged in a horizontal plane at 45 mm above the bottom are in reasonable agreement with the experimental data. The fluctuation of the simulated bed height over time is similar to the experimental measurement results. The reproduction of these qualitative and quantitative results indicates that the proposed method effectively improves simulation performance and accuracy.

Tight junctions (TJ) are cell-cell adhesive structures that define the permeability of barrier-forming epithelia and endothelia. In contrast to this seemingly static function, TJs display a surprisingly high molecular complexity and an unexpected dynamic regulation, which allows the TJs to maintain a barrier in the presence of physiological forces and in response to perturbations. Cell-cell adhesion receptors play key roles during the dynamic regulation of TJs. They connect individual cells within cellular sheets and link sites of cell-cell contacts to the underlying actin cytoskeleton. Recent findings support roles of adhesion receptors in transmitting mechanical forces and promoting phase separation. In this review, we discuss newly discovered functions of cell adhesion receptors localized at the TJs and their role in the regulation of the barrier function.

In the course of SARS-CoV-2 infection, the 3CL or nsp5 protease plays a pivotal role as the most important viral protease required for the maturation of viral proteins during host infection. Herein, we simulated for 500 ns 3CLproWT, 3CLproH41A, , 3CLproBeta, and 3CLproOmicron, in complex with the substrates nsp 4|5 and nsp 5|6. Our results show that mutations in the 3CLpro present in the SARS-CoV-2 variant of concerns (VOCs) did not lead to significant conformational changes or changes in substrate binding affinities. However, significantly high cleavage rates for the boundary between nsp4 and nsp5 were obtained for 3CLproBeta and 3CLproOmicron and may play a key role in the viral replication and virus fitness gain. Our molecular dynamics data suggest that the cleavage rate of nsp4|5 may be related to the increased amount of viral load observed for these VOCs, releasing more nsp4 than other non-structural proteins. Based on our hydrogen bonding analyses, we also discovered that Gly143 and Glu166 are key residues in substrate recognition, suggesting that these residues may be incorporated as pharmacophoric centers for Beta and Omicron variants in drug design. Our results suggest that Gly143 and Glu166 are essential residues to interact with Gln6 of the different substrates and, therefore, are potential broad-spectrum pharmacophoric centers of SARS-CoV-2 3CLpro.

Hydraulic fracturing is the most effective method of stimulation for hydrocarbon reservoirs. However the use of water-based fracturing fluids, can be a problem in water-sensitive formations due to the permeability damage hazard caused by clay minerals swelling. For this reason, the foamed fracturing fluids with addition of natural, fast hydrating guar gum were examined. The rheology and filtration coefficients of foamed fracturing fluids were examined and compared to the properties of conventional water-based fracturing fluid. Laboratory results provided the input for numerical simulation of the fractures geometry for water-based fracturing fluids and 50% N₂ foamed fluids. The results show, that the foamed fluids were able to create shorter and thinner fractures compared to the fractures induced by the non-foamed fluid. The simulation proved that the concentration of proppant in the fracture and its conductivity are similar or slightly higher when using the foamed fluid. Moreover such fluids are able to significantly reduce the amount of water necessary for fracturing treatments, limiting clay minerals swelling, and reducing the reservoir permeability damage. The foamed fluids, when injected to the reservoir, provide additional energy, that allow for more effective flowback, and maintain the proper fracture geometry and proppant placing. The results of laboratory work in combination with the 3D simulation showed, that the foamed fluids have suitable viscosity which allows opening the fracture, and transport the proppant into the fracture, providing successful fracturing operation.

The tricellular tight junction protein angulin-1/lipolysis-stimulated lipoprotein receptor (LSR) is associated with numerous signal transduction pathways that regulate gene expression, epithelial cell function, and morphogenesis. The impacts of particulate matter on LSR and asthma are not known. Therefore, the aim of the present study was to assess the impacts of particulate matter as one of air pollution on LSR expression in asthma. To this end, we investigated the involvement of LSR in the pathogenesis of bronchial asthma in an ovalbumin (OVA) mouse model of asthma exposed to titanium dioxide (TiO2), and examined plasma LSR levels and their relationship with clinical variables in asthma patients and healthy controls. Plasma LSR levels were lower in asthma patients than in healthy controls. Although plasma levels of receptor for advanced glycation end-products (RAGE) were lower in asthma patients than in healthy controls, high mobility group box protein 1 (HMGB1) levels did not differ between these groups. In asthma patients, LSR levels were correlated with the forced vital capacity percentage, ratio of forced expiration (first second)-to-forced vital capacity, white blood cell count, and blood lymphocyte proportion. In mice, LSR expression was significantly decreased in OVA mice compared with control mice, and changed in OVA mice exposed to TiO2. Transforming growth factor beta levels were increased in OVA mice compared with control mice, and changed in OVA mice exposed to TiO2. Overall, plasma LSR levels were decreased in asthma patients, and were influenced by exposure to air pollution. These findings suggest the involvement of LSR in the pathogenesis of asthma.

Structural and biochemical studies have recently revealed a range of rationally engineered nanobodies with efficient neutralizing capacity against SARS-CoV-2 virus and resilience against mutational escape. In this work, we combined atomistic simulations and conformational dynamics analysis with the ensemble-based mutational profiling of binding interactions for a diverse panel of SARS-CoV-2 spike complexes with nanobodies. Using this computational toolkit, we identified dynamic signatures and binding affinity fingerprints for the SARS-CoV-2 spike protein complexes with nanobodies Nb6 and Nb20, VHH E, a pair combination VHH E+U, a biparatopic nanobody VHH VE, and a combination of CC12.3 antibody and VHH V/W nanobodies. Through ensemble-based deep mutational profiling of stability and binding affinities, we identify critical hotspots and characterize molecular mechanisms of SARS-CoV-2 spike protein binding with single ultra-potent nanobodies, nanobody cocktails and biparatopic nanobodies. By quantifying dynamic and energetic determinants of the SARS-CoV-2 S binding with nanobodies, we also examine the effects of circulating variants and escaping mutations. We found that mutational escape mechanisms may be controlled through structurally and energetically adaptable binding hotspots located in the host receptor-accessible binding epitope that are dynamically coupled to the stability centers in the distant epitope targeted by VHH U/V/W nanobodies. The results of this study suggested a mechanism in which through cooperative dynamic changes, nanobody combinations and biparatopic nanobody can modulate the global protein response and induce the increased resilience to common escape mutants.

The study of fractures in the subsurface is very important in unconventional reservoirs since they are the main conduits for hydrocarbon flow. For this reason, a variety of equivalent medium theories have been proposed for the estimation of fracture and fluid properties within reservoir rocks. We experimentally investigated the feasibility of applying the Galvin model in fractured tight stones. For this proposal, three artificial fractured tight sandstone samples with the same background porosity (11.7% ± 1.2%) but different fracture densities of 0.00, 0.0312, and 0.0624 were manufactured. The fracture thickness was 0.06 mm and the fracture diameter was 3 mm in all the fractured samples. Ultrasonic P- and S-wave velocities were measured at 0.5 MHz in a laboratory setting in dry and water-saturated conditions in directions at 0°, 45°, and 90° to the fracture normal. The results were compared with theoretical predictions of the Galvin model. The comparison showed that model predictions significantly underestimated P- and S- wave velocities as well as P-wave anisotropy in water-saturated conditions, but overestimated P-wave anisotropy in dry conditions. By analyzing the differences between the measured results and theoretical predictions, we modified the Galvin model by adding the squirt flow mechanism to it and used the Thomsen model to obtain the elastic moduli in high- and low-frequency limits. The modified model predictions showed good fits with the measured results.

With the popularization of natural gas and the requirements for environmental protection, the development and utilization of natural gas is particularly important. The status of natural gas in China's oil and gas exploration and development is constantly improving, and the country is paying more and more attention to the exploitation and utilization of natural gas. The Upper Paleozoic tight sandstone in the Ordos Basin is characterized by low porosity, low permeability and large area of concealed gas reservoirs. By injecting CO2 into the formation, the recovery of natural gas can be improved, and at the same time, the stable storage of CO2 can be achieved to achieve a win-win situation of CO2 emission reduction and utilization. Injecting greenhouse gas CO2 into gas reservoirs for storage and improving recovery has also become a hot research issue. In order to improve the recovery efficiency of tight sandstone gas reservoir, this paper takes the complex tight sandstone of Upper Paleozoic in Ordos Basin as the research object, through indoor physical simulation experiments, carried out the influence of displacement rate, fracture dip angle, core permeability, core dryness and wetness on CO2 gas displacement efficiency and storage efficiency, and analyzed the influence of different factors on CO2 gas displacement efficiency and storage efficiency to improve the recovery and storage efficiency. The research results show that under different conditions, when the CO2 injection pore volume is less than 1PV, the relationship between the CH4 recovery rate and the CO2 injection pore volume is linear, and the tilt angle is 45 °. When the CO2 injection pore volume exceeds 1PV, the CH4 recovery rate increases slightly with the increase of displacement speed, the recovery rate of CO2 displacement CH4 is between 87% - 97%, and the CO2 breakthrough time is 0.7PV-0.9PV. In low-permeability and low-speed displacement cores, the diffusion of CO2 molecules is more significant. The lower the displacement speed is, the earlier the breakthrough time of CO2 is, and the final recovery of CH4 slightly decreases. Gravity has a great impact on CO2 storage and enhanced recovery. The breakthrough of high injection and low recovery of CO2 is earlier, and the recovery of CH4 is about 3.3% lower than that of low injection and high recovery. The bound water makes the displacement phase CO2 partially dissolved in the formation water, and the CO2 breakthrough lags about 0.1PV. Ultimately, CH4 recovery factor and CO2 storage rate are higher than those of dry core displacement. The research results provide theoretical data support for CO2 injection to improve recovery and storage efficiency in complex tight sandstone gas reservoirs.

In the intricate dance between viruses and host cells, RNA-binding proteins (RBPs) serve as crucial orchestrators of gene expression and cellular processes. We will delve into the riveting realm of viral hijacking, where viruses deftly exploit host RBPs to manipulate cellular machinery for their replication and pathogenesis. Through a masterstroke of molecular subterfuge, viruses co-opt RBPs involved in various facets of RNA metabolism - from transcription to degradation - to promote viral gene expression and evade host immune defenses. This manipulation leads to global alterations in cellular RNA metabolism, impacting essential host genes vital for immune responses and homeostasis.Unveiling this clandestine alliance between viruses and RBPs is not just a scientific pursuit, but a critical imperative for devising innovative antiviral strategies to disrupt these interactions. By unraveling the intricate interplay between viral proteins and host RBPs throughout the viral life cycle - from entry to assembly - researchers aim to identify vulnerabilities that can be targeted for therapeutic intervention. Strategies such as disrupting viral protein-RBP interactions hold promise for inhibiting viral replication and curbing pathogenesis, offering a beacon of hope in the battle against viral infections.Our manuscript elucidates the indispensable roles played by RBPs in viral replication, translation, and pathogenesis, shedding light on the molecular mechanisms driving viral success. Delving deeper, it explores how viruses intricately entwine with host RBPs, manipulating cellular signaling pathways to create a hospitable environment for viral spread. By dissecting these manipulative tactics, researchers uncover new targets for antiviral therapy, envisioning a future where tailored interventions disrupt viral-host RBP interactions with precision and efficacy.As the narrative unfolds, the therapeutic implications of targeting RBPs or their interactions with viral proteins emerge as a promising frontier in the fight against viral infections. From small molecule inhibitors to RNA-based therapeutics, innovative approaches are on the horizon, offering new avenues for combatting viral diseases. We set the stage for future research, beckoning researchers to delve deeper into the molecular intricacies of viral hijacking of RBPs and charting a course towards novel therapeutic interventions that promise to reshape the landscape of antiviral therapy.\Ultimately, it beckons the scientific community to embark on a voyage of discovery, unraveling the secrets of viral hijacking of RBPs and paving the way for transformative advances in antiviral therapeutics. The stage is set for a new chapter in the battle against viral infections, where knowledge becomes the sword and innovation the shield against the pernicious machinations of viral pathogens.

Scholars of the English Literature unanimously say that J.R.R. Tolkien and C.S. Lewis influenced each other writings. For the first time we have investigated this issue mathematically by using an original multi–dimensional analysis of linguistic parameters, based on surface deep–language variables and linguistic channels. To set our investigation in the framework of the English Liter-ature, we have considered some novels written by earlier authors, such as C. Dickens, G. Mac-Donald and others. The deep–language variables and the linguistic channels, discussed in the paper, are likely due to writers’ unconscious design and have revealed connections between texts far beyond writers’ awareness. In summary, the capacity of the extended short–term memory required to readers, the universal readability index of texts, the geometrical representation of texts and the fine tuning of linguistic channels within texts – all tools largely discussed in the paper – have revealed strong connections between The Lord of the Rings (Tolkien), The Chronicles of Narnia, The Space Trilogy (Lewis) and novels by MacDonald, therefore in agreement with what scholars of the English Literature say.

In 2001, the first author of this paper reported that schizophrenia is associated with an increased frequency of the haptoglobin (Hp)-2 gene. The precursor of Hp-2 is zonulin, a molecule that affects intercellular tight junction integrity. Recently, we reported increased plasma IgA/IgM responses to Gram-negative bacteria in deficit schizophrenia indicating leaky gut and gut dysbiosis. The current study was performed to examine the integrity of the paracellular (tight and adherens junctions) and transcellular (cytoskeletal proteins) pathways in deficit versus non-deficit schizophrenia. We measured IgM responses to zonulin, occludin, E-cadherin, talin, actin and vinculin in association with IgA responses to Gram-negative bacteria, CCL-11, IgA responses to tryptophan catabolites (TRYCATs), immune activation and IgM to malondialdehyde (MDA) and NO-cysteinyl in 78 schizophrenia patients and 40 controls. We found that the ratio of IgM to zonulin + occudin / talin + actin + viculin (PARA/TRANS) was significantly greater in deficit than in non-deficit schizophrenia and higher in schizophrenia than controls and was significantly associated with increased IgA responses to Gram-negative bacteria. IgM responses to zonulin were positively associated with schizophrenia (versus controls), while IgM to occludin was significantly associated with deficit schizophrenia (versus non-deficit schizophrenia and controls). A large part of the variance (90.8%) in negative and PHEM (psychosis, hostility, excitation and mannerism) symptoms was explained by PARA/TRANS ratio, IgA to Gram-negative bacteria, IgM to E-cadherin and malondialdehyde (MDA) and memory dysfunctions, while 53.3% of the variance in the latter was explained by PARA/TRANS ratio, IgA to Gram-negative bacteria, CCL-11, TRYCATs and immune activation. The results show an upregulated paracellular pathway with breakdown of the tight and ahherens junctions and increased bacterial translocation in deficit schizophrenia. These dysfunctions in the intestinal paracellular route together with lowered natural IgM, immune activation and production of CCL-11 and TRYCATs contribute to the phenomenology of deficit schizophrenia.

The temporo-spatial theory of consciousness postulates that the brain implements its own inner time and space for conscious processing of outside world. In line with this hypothesis, I present a mechanistic model of mutually connected processes that encode the inner space and time. The model is used to elaborate the binding mechanism between two sets of processes representing the inner space and time, respectively. Further, a stochastic version of the model is developed to investigate the interplay between the binding strength and noise. The spectral entropy is used to characterize the noise effects on the systems of processes when the binding strength is varied. The stochastic modeling results reveal that the spectral entropy values for strongly bound systems have similar entropy values for weakly bound or even decoupled systems. Thus, the analysis performed in this study allows to conclude that the binding mechanism is noise-resilient.

The Source Region of the Yellow River (SRYR), renowned as the "Water Tower of the Yellow River", serves as an important water conservation domain in the upper reaches of the Yellow River, significantly influencing water resources within the basin. Based on the Weather Research and Forecasting (WRF) Model Hydrological modeling system (WRF-Hydro) model, the key variables of the atmosphere-land-hydrology coupling processes over the SRYR during the 2013 rainy season are analyzed. The investigation involves a comparative analysis between the coupled WRF-Hydro and the standalone WRF simulations, focusing on the hydrological response to the atmosphere. The results reveal that the WRF-Hydro model's proficiency in depicting streamflow variations over the SRYR, yielding a Nash Efficiency Coefficient (NSE) of 0.44 and of  0.61 during the calibration and validation periods, respectively. Compared to the standalone WRF simulations, the coupled WRF-Hydro model demonstrates enhanced performance in soil heat fluxes simulations, reducing the mean Root Mean Square Error (RMSE) of surface soil temperature by 0.96 K and soil moisture by 0.01 m\textsuperscript{\textsuperscript{3}}/m\textsuperscript{\textsuperscript{3}}. Furthermore, the coupled model adeptly captures the streamflow variation characteristics with the NSE of 0.33. This underscores the significant potential of the coupled WRF-Hydro model for delineating coupled atmosphere-land-hydrology processes in regions characterized by cold climates and intricate topography.

The deep Cretaceous tight sandstone in Kuqa Depression of Tarim foreland basin is an ultra-low porosity and ultra-deep gas-bearing reservoir, which is characterized by small pores, fine throats, and poor connectivity. The wireline logging responses are so complex, and especially, it is difficult to identify fluid types from resistivity logs. Based on acoustic, density, and neutron logs response differences in gas and water layers, effective fluid sensitivity factors are constructed for gas layer identification. From conventional logs, acoustic-neutron porosity difference, density-neutron porosity difference, and triple-porosity ratio are all sensitive parameters to the gas layer. From the NMR logging response mechanism, the density and NMR porosity difference, and T2 geometric mean of the movable fluid are also two sensitive parameters to the gas layer. Based on these parameters, a series of fluid typing charts are constructed and their adaptabilities are analyzed and compared. By contrast, NMR log interpretation is better, and triple-porosity comprehensive method from conventional logs is also effective when NMR logging is not available. Finally, the comprehensive fluid typing strategy by combining some methods for ultra-low porosity tight sandstone is summarized and optimized. This study is another alternative for fluid identification using non-electrical logs.

We present a generic and open-source framework for the numerical modeling of the expected transport and storage mechanisms in unconventional gas reservoirs. These unconventional reservoirs typically contain natural fractures at multiple scales. Considering the importance of these fractures in shale gas production, we perform a rigorous study on the accuracy of different fracture models. The framework is validated against an industrial simulator and is used to perform a history-matching study on the Barnett shale. This work presents an open-source code that leverages cutting-edge numerical modeling capabilities like automatic differentiation, stochastic fracture modeling, multi-continuum modeling and other explicit and discrete fracture models. We modified the conventional mass balance equation to account for the physical mechanisms that are unique to organic-rich source rocks. Some of these include the use of an adsorption isotherm, a dynamic permeability-correction function, and an embedded discrete fracture model (EDFM) with fracture-well connectivity. We explore the accuracy of the EDFM for modeling hydraulically-fractured shale-gas wells, which could be connected to natural fractures of finite or infinite conductivity, and could deform during production. Simulation results indicates that although the EDFM provides a computationally efficient model for describing flow in natural and hydraulic fractures, it could be inaccurate under these three conditions: 1. when the fracture conductivity is very low. 2. when the fractures are not orthogonal to the underlying Cartesian grid blocks, and 3. when sharp pressure drops occur in large grid blocks with insufficient mesh refinement. Each of these results are very significant considering that most of the fluids in these ultra-low matrix permeability reservoirs get produced through the interconnected natural fractures, which are expected to have very low fracture conductivities. We also expect sharp pressure drops near the fractures in these shale gas reservoirs, and it is very unrealistic to expect the hydraulic fractures or complex fracture networks to be orthogonal to any structured grid. In conclusion, this paper presents an open-source numerical framework to facilitate the modeling of the expected physical mechanisms in shale-gas reservoirs. The code was validated against published results and a commercial simulator. We also performed a history-matching study on a naturally-fractured Barnett shale-gas well considering adsorption, gas slippage & diffusion and fracture closure as well as proppant embedment, using the framework presented. This work provides the first open-source code that can be used to facilitate the modeling and optimization of fractured shale-gas reservoirs. To provide the numerical flexibility to accurately model stochastic natural fractures that are connected to hydraulically-fractured wells, it is built atop other related open-source codes. We also present the first rigorous study on the accuracy of using EDFM to model both hydraulic fractures and natural fractures that may or may not be interconnected.

This study investigates the stochastic dynamics of hepatitis B virus (HBV) infection using a newly proposed stochastic model. Contrary to deterministic models that fail to encapsulate the inherent randomness and fluctuations in biological processes, our stochastic model provides a more realistic representation of HBV infection dynamics. It incorporates random variability, thereby acknowledging the changes in viral and cellular populations and uncertainties in parameters such as infection rates and immune responses. We examine the solution's existence, uniqueness, and positivity for the proposed model, followed by a comprehensive stability analysis. We provide the necessary and sufficient conditions for local and global stability, offering deep insight into the infection dynamics. Furthermore, we utilize numerical simulations to corroborate our theoretical results. This research provides a robust tool for understanding the complex behavior of HBV dynamics, contributing significantly to the ongoing quest for more effective HBV control and prevention strategies.

Into the frame of the French TANDEM project (Tsunamis in the Atlantic and the English ChaNnel: Definition of the Effects through numerical Modelling) Principia has been working on the development and qualification of two in-house CFD software’s: the 2D EOLE-SV (Saint Venant) model for simulation of large scale tsunami propagation from the source up to coastal scale and the 3D EOLE-NS (Navier-Stokes) model dedicated to tsunami coastal impact modelling. This paper presents a large range of test cases carried out into the frame of the project and dedicated to the validation of numerical codes in various tsunami wave conditions. The main aspects of phenomena such as wave generation, propagation and coastal impact are investigated on academic situations. A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.

Efficiency improvement is the new challenge in all fields of design. In this scenario the reduction of power losses is becoming more and more a main concern also in the design of power transmissions. Appropriate models to predict power losses are therefore from the earliest stages of the design phase. The aim of the project is to carry on lubrication simulations of several variants of a cylindrical-roller-bearing to understand the lubricant distribution and the related churning power losses. Several strategies to reduce the computational effort have been used. Among them the sectorial symmetry and three innovative meshing strategies (purely analytical with and without interfaces and analytical/subtractive) that have been implemented in the OpenFOAM® environment. The results of the different approaches were compared among them and with experimental observations showing good agreement and reasonable savings in terms of computational effort.

Pore structure characterization and fractal analysis have great significance for understanding and evaluating tight limestone reservoirs. In this work, the pore structure of tight limestone is characterized combined with low-temperature nitrogen adsorption (LTNA) and low-field nuclear magnetic resonance (NMR), and the fractal dimension of tight limestone pore structure is discussed based on LTNA and NMR data. The results indicate that the pores of tight limestone have H3 and H4 types, the pore size distribution (PSD) of H3 type is a wave distribution ranging from 2 to 10 nm, and the PSD of H4 type is a unimodal distribution ranging from 2 to 10 nm. The transverse relaxation time (T2) spectrum of tight limestone shows single peak (DF), double peak (SF) and three peak (TF), and the T2 spectrums for micropores, mesopores and macropores range from 0.1 to 10 ms, from 10 to 100 ms, and greater than 100 ms, respectively. The LTNA fractal dimension of tight limestone (DL) ranges between 2.4446 and 2.7688, with an average of 2.5729, and the NMR fractal dimensions of micropores (DNMR1), mesopores (DNMR2) and macropores (DNMR3) are distributed between 0.3744 to 1.1293, 2.4263 to 2.9395 and 2.6582 to 2.9989, respectively. Moreover, there is no correlation between fractal dimension and permeability and porosity, a negative correlation between DL and average pore radius, a positive correlation between DL and specific surface area, and a positive correlation between DNMR2 and DNMR3 and micropore content, while DNMR2 and DNMR3 are negatively correlated with the content of mesopores and macropores.

Many studies discuss carbon-based materials because of the versatility of its element. They include different opinions for scientific problems and discuss fairly at convincing and compelling levels within the scope and application. A gas-state carbon atom converts into various states depending on its conditions of processing. The electron transfer mechanism in the gas-state carbon atom is responsible to convert it into various states, namely, graphite, nanotube, fullerene, diamond, lonsdaleite and graphene. The shape of ‘energy trajectory’ enables transferring electrons from the left- and right-sides of an atom is like a parabola. That ‘energy trajectory’ is linked to states (filled state and suitable nearby unfilled state) where force-exertion along the poles of transferring electrons is remained balance. So, the mechanism of originating different states of a gas-state carbon atom is under the involvement of energy first. This is not the case for atoms executing confined inter-state electron-dynamics as the force is involved first. Graphite-, nanotube- and fullerene-state atoms ‘partially evolve partially develop’ (form) their structures. These possess one-dimensional, two-dimensional and four-dimensional ordering of atoms, respectively. Their structural formation also comprises ‘energy curve’ having a shape-like parabola. Transferring suitable filled state electron to suitable nearby unfilled state is under a balance force exerting along the poles. The graphite structure under only attained-dynamics of atoms can also be formed but in two-dimension. Here, binding energy between graphite-state carbon atoms is for a small difference of exerting forces along their opposite poles. Structural formation in diamond, lonsdaleite and graphene atoms involve energy to gain required infinitesimal displacements of electrons through which they maintain orientationally-controlled exerting forces along dedicated poles. In this study, the growth of diamond is found to be south to east-west (ground) where atoms bound ground to south. Thus, diamond atoms merge for a tetra-electron ground to south topological structure. Lonsdaleite atoms merge for a bi-electron ground to just-south topological structure. The growth of graphene is found to be north to ground where atoms bound ground to north. Thus, graphene atoms merge for a tetra-electron ground to north topological structure. Glassy carbon exhibits layered-topological structure where, tri-layers of gas-, graphite- and lonsdaleite-state atoms successively bind in repetitive order. Nanoscale hardness is also sketched based on different force-energy behaviors of different state carbon atoms. Here, structure evolution in each carbon state atom explores its own science.

Microemulsions with excellent properties have shown significant potential for enhancing oil recovery from tight formations through spontaneous imbibition. This work present a high-temperature-tolerant lower-phase microemulsion (HTLP-ME) using a microemulsion dilution method. The properties and morphological characteristics of HTLP-ME were evaluated and proposed a enhanced spontaneous imbibition oil recovery mechanism using imbibition tests and CT scanning technology. The results show that the optimum concentration of HTLP-ME is 0.2wt% has a superior thermal stability, small particle size, lower IFT, good wettability alteration ability, and minimum adsorption loss. Imbibition and CT experiments showed that the longer the imbibition distance, the better is the effect on spontaneous imbibition oil recovery. Furthermore, the solubilization coefficient and self-driving force were defined and calculated to quantitatively analyze the imbibition mechanisms. The imbibition mechanisms demonstrate that the reduction of oil/solid adhesion is due to the synergistic effect of IFT reduction and wettability alteration, and the ability to increase the imbibition distance through a larger self-driving force.

Comparison with conventional reservoirs, the prediction of pressure in ultra-deep tight sandstone formations is difficult. The prediction of seismic pressure is more challenging than well logging pressure prediction. The main methods for seismic pressure prediction include equivalent depth method, Eaton method, Fillippone formula, and its modified versions. Among them, the Eaton method is widely used and has good effectiveness. However, this method relies on difficult-to-obtain normal compaction trend line, which leads to low prediction accuracy in space. To address this issue, a method combining the Eaton method and Collocated cokriging is proposed. The Eaton formula is used to predict formation pressure at the well, with compressional wave velocity as the covariate for predicting the main variable, formation pressure. By simulating the shear wave velocity based on seismic compressional wave velocity, the influence of various parameters on the prediction results is analyzed, and the accuracy of this method is verified by comparing it with other methods. The proposed method is then applied to predict formation pressure in the ultra-deep formations of the Junggar Basin. The simulation results show that the Collocated cokriging method achieves higher planar accuracy and better matches the experimental expectations in terms of prediction results. The application results also demonstrate the scientific effectiveness of the combined method, which has achieved good results in practical production applications.

Intermediate-luminosity-optical-transients (ILOTs) are stellar outbursts with luminosity between those of classical novae and supernovae. They are divided into a number of sub-groups depending on the erupting progenitor and the properties of the eruption. Many of the ILOTs sit on the slanted Optical Transient Stripe (OTS) in the Energy-Time Diagram (ETD) that shows their total energy vs. duration of their eruption. We describe the different kinds of ILOTs that populate the OTS and other parts of the ETD. We also stand on similarities between Planetary Nebulae (PN) to ILOTs, and suggest that some PNe were formed in an ILOT event. The high energy part of the OTS is reserved to the supernova impostors -- giant eruption of very massive stars. We show results of 3D hydrodynamical simulations of supernova impostors that expose the mechanism behind these giant eruptions, and present new models for recent ILOTs. We stand on the connection between different kinds of ILOTs, and suggest that they are powered by a similar source of energy -- gravitational energy released by mass transfer.

Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can provide predictions, reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.

In this paper we show how to efficiently implement parallel discrete simulations on Multi-Core and GPU architectures through a real example of application: a cellular automata model of laser dynamics. We describe the techniques employed to build and optimize the implementations using OpenMP and CUDA frameworks. We have evaluated the performance on two different hardware platforms that represent different target market segments: high-end platforms for scientific computing, using an Intel Xeon Platinum 8259CL server with48cores and also an NVIDIA Tesla6V100 GPU, both running on Amazon Web Server (AWS) Cloud, and on a consumer-oriented platform, using an Intel Core i9 9900k CPU and an NVIDIA GeForce GTX 1050 TI GPU. Performance results are compared and analysed in detail. We show that excellent performance and scalability can be obtained in both platforms, and we extract some important issues that imply a performance degradation for them. We also found that current Multi-Core CPUs with large core numbers can bring a performance very near to that of GPUs, even similar in some cases.

The microscopic pore structure of tight sandstone reservoir significantly influences the characteris-tics of CO2 flooding after water flooding. The research was conducted using various techniques such as casting thin sections, high pressure mercury injection, scanning electron microscopy, Nu-clear magnetic resonance (NMR) testing, and a self-designed high-temperatures and high-pressures microscopic visualization displacement system. Three type cores with different pore structures were selected for the flooding experiments and the microscopic visualization displace-ment experiments, including CO2 immiscible flooding, near-miscible flooding, and miscible flood-ing after the conventional water flooding. The characteristics of CO2 flooding and the residual oil distribution after water flooding were quantitatively analyzed and evaluated. The results show that: (1) During the water flooding process, the oil produced from type I and type III samples mainly comes from large and some medium pores. Oil utilization of all pores is significant for type II samples. The physical properties and pore types have a greater impact on water flooding. Type I and II samples are more suitable for near-miscible flooding after water flooding. Type III samples are more suitable for miscible flooding after water flooding. (2) In the CO2 flooding, oil recovery increases gradually with increasing pressure for all three types of samples. Type II core samples have the highest recovery. Before miscibility, the oil recovered from type I and type II samples was primarily from large pores, however, the oil recovery mainly comes from medium pores when reaching miscibility. As for the type III samples, the oil produced in the immiscible state mainly comes from large and medium pores, and the enhanced oil recovery mainly comes from medium and small pores after reaching the near-miscible phase. (3) It can be seen from the microscopic re-sidual oil distribution that oil recovery will increase as the petrophysical properties of the rock model improve. The oil recovery of near-miscible flooding after water flooding regarding the type II model is up to 68.11%. The oil recovery of miscible flooding after water flooding about type III model is the highest at 74.57%. With increasing pressure, the proportion of flake residual oil grad-ually decreases, while the proportion of droplet-like and film-like residual oil gradually increases. Type II samples have a relatively large percentage of reticulated residual oil in the near-miscible stage.

Polymers' controlled pyrolysis is an economical and environmentally friendly solution to prepare activated carbon. However, due to the experimental difficulty in measuring the dependence between tissues and pyrolysis parameters at high temperatures, the unknown pyrolysis mechanism hinders access to the target products with desirable morphologies and performances. In this study, we investigate the pyrolysis process of polystyrene under different heating rates and temperatures employing reactive molecular dynamics (ReaxFF-MD) simulations. A clear profile of the generation of pyrolysis products determined by the temperature and heating rate is constructed. It is found that the heating rate affects the type and amount of pyrolysis intermediates and their timing, and low-rate heating helps yield more diverse pyrolysis intermediates. While the temperature affects the pyrolytic tissues of the final equilibrium products, either too low or too high a target temperature is detrimental to generating large areas of graphitized tissue. The established theoretical evolution process matches experiments well, thus contributing to preparing target activated carbons by referring to the regulatory mechanism of pyrolytic tissues.

Eukaryotic cells contain membranes with various curvatures, from the near-plane plasma membrane to the highly curved membranes of organelles, vesicles, and membrane protrusions. These curvatures are generated and sustained by curvature-inducing proteins, peptides, and lipids, and describing these mechanisms is an important scientific challenge. In addition to that, some molecules can sense membrane curvature and a thereby be trafficked to specific locations. The description of curvature-sensing is another fundamental challenge. Curved lipid membranes and their interplay with membrane-associated proteins can be investigated with molecular dynamics (MD) simulations, given the right tools. Various methods for simulating curved membranes with MD are discussed here, including tools for setting up simulation of vesicles, and methods for sustaining membrane curvature. The latter are divided into methods that exploit scaffolding virtual beads, methods that use curvature-inducing molecules, and methods applying virtual forces. The variety of simulation tools allow the researcher to closely match the conditions of experimental studies of membrane curvatures.

Abstract: Plasma hemopexin (HPX) is the key antioxidant protein of the endogenous clearance pathway that limits the deleterious effects of heme released from hemoglobin and myoglobin. During intra-vascular hemolysis, heme partitioning to protein and lipid increases as the plasma concentration of HPX declines. Therefore, the development of HPX as a replacement therapy during high heme stress could be a relevant intervention for hemolytic disorders. A logical approach to enhance HPX yield involves recombinant production strategies from human cell lines. The present study focuses on a biophysical assessment of heme binding to recombinant human HPX (rhHPX) produced in the Expi293FTM (HEK293) cell system. In this report, we examine rhHPX in comparison with plasma HPX using a systematic analysis of protein structural and functional characteristics related to heme binding. Analysis of rhHPX by UV/Vis absorption spectroscopy, circular dichroism (CD), SEC-HPLC and catalase-like activity demonstrated a similarity to HPX fractionated from plasma. In particular, the titration of HPX apo-protein(s) with heme was performed for the first time using a wide range of heme concentrations to model HPX-heme interactions to approximate physiological conditions (from extremely low to more than 2-fold heme excess). The CD titration data showed an induced bisignate CD Soret band pattern typical for plasma and rhHPX versions at low heme-to-protein molar ratios and demonstrated that further titration is dependent on the amount of protein-bound heme to the extent that the arising opposite CD couplet results in a complete inversion of the observed CD pattern. The data generated in this study suggests more than one binding site in both plasma and rhHPX. Further, our study provides a useful analytical platform for detailed characterization of HPX-heme interactions and potentially, novel HPX fusion constructs.

Numerical models widely used for hydrocarbon phase behavior and compositional flow simulations are based on assumption of thermodynamic equilibrium. However, it is not uncommon for oil and gas-condensate reservoirs to exhibit essentially non-equilibrium phase behavior, e.g., in the processes of secondary recovery after pressure depletion below saturation pressure, or during gas injection, or for condensate evaporation at low pressures. In many cases the ability to match field data with equilibrium model depends on simulation scale. The only method to account for non-equilibrium phase behavior adopted by the majority of flow simulators is the option of limited rate of gas dissolution (condensate evaporation) in black oil models. For compositional simulations no practical yet thermodynamically consistent method has been presented so far except for some upscaling techniques in gas injection problems. Previously reported academic non-equilibrium formulations have a common drawback of doubling the number of flow equations and unknowns compared to the equilibrium formulation. In the paper a unified thermodynamically-consistent formulation for compositional flow simulations with non-equilibrium phase behavior model is presented. Same formulation and a special scale-up technique can be used for upscaling of an equilibrium or non-equilibrium model to a coarse-scale non-equilibrium model. A number of test cases for real oil and gas-condensate mixtures are given. Model implementation specifics in a flow simulator are discussed and illustrated with test simulations. A non-equilibrium constant volume depletion algorithm is presented to simulate condensate recovery at low pressures in gas-condensate reservoirs. Results of satisfactory model matching to field data are reported and discussed.

Computational modeling (CM) is a versatile scientific methodology used to examine the properties and behavior of complex systems, such as polymeric materials for biomedical bioengineering. CM has emerged as a primary tool for predicting, setting up, and interpreting experimental results. Integrating in silico and in vitro experiments accelerates scientific advancements, yielding quicker results at a reduced cost. While CM is a mature discipline, its recent prominence in Bio-medical Engineering for biopolymer materials has recently gained prominence. In biopolymer biomedical engineering, CM focuses on three key research areas: A) Computer-aided design (CAD/CAM), B) Finite Element Analysis, and C) Molecular Dynamic Simulations. This review centers on Molecular Dynamics (MD) simulations in biopolymers, analyzing structural, functional, and evolutionary aspects of biomolecular systems over time. MD simulations solve Newton's equations of motion for all-atom systems, generating spatial trajectories for each atom, providing insights into properties like water absorption on biopolymer surfaces and interactions with solid surfaces, crucial for biomaterial assessment. This review offers readers a comprehensive overview of the diverse applications of Molecular Dynamics (MD) simulations on biopolymers. Furthermore, it emphasizes the flexibility, robustness, and synergistic relationship between in silico and experimental techniques.

X-rays production through betatron radiation emission from electron bunches is a valuable resource for several research fields. The EuAPS (EuPRAXIA Advanced Photon Sources) project, within the framework of the EuPRAXIA project, aims to provide 1-10 keV photons (soft X-rays) developing a compact plasma-based system designed to exploit self-injection processes that occur in highly nonlinear laser-plasma interaction (LWFA) to drive electron betatron oscillations. Since the emitted radiation spectrum, intensity, angular divergence, and possible coherence strongly depend on the properties of the self-injected beam, accurate preliminary simulations of the process are necessary to evaluate the optimal diagnostic device specifications. Dedicated tools for these tasks are under development: electron trajectories from particle-in-cell (PIC) simulations are currently undergoing numerical analysis through the calculation of retarded fields and spectra for various plasma and laser parameter combinations. At the same time, an analytical approach is being examined to obtain an expression for the emitted spectrum, assuming linear or quadratic electron energy variation. Code structure and some trajectories analysis results are presented.

Modeling the flow of polymer solutions requires knowledge at various length and time scales. The macroscopic behavior is described by the overall velocity, pressure and stress. The polymeric contribution to the stress requires knowledge of the evolution of polymer chains. In this work, we use a microscopic model, the finitely extensible nonlinear elastic (FENE) model, to capture the polymer behavior. The benefit of microscopic models is they remain faithful to the polymer dynamics without information loss via averaging. The downside is the computational cost needed to solve the thousands to millions of differential equations describing the microstructure. We thus describe a multiscale flow solver that utilizes GPUs for massively parallel, efficient simulations. We compare and contrast the microscopic model with its macroscopic counterpart under various flow conditions. In particular, significant differences are observed under nonlinear flow conditions, where the polymers become highly stretched and oriented.

Tsunami modeling and simulation has changed in the past few years more than it had in decades, especially so with respect to coastal inundation. Among other things, this change is supported by the approaching era of exa-scale computing, whether via GPU or more likely forms of hybrid computing whose presence is growing across the geosciences. For reasons identified across this review, exa-scale computing efforts will impact the on-shore, highly turbulent r\'egime to higher degree than the 2D shallow water equations used to model tsunami propagation in the open ocean. This short review describes the different approaches to tsunami modeling from generation to impact and underlines the limits of each model based on the flow r\'egime. Moreover, from the perspective of a future comprehensive multi-scale modeling infrastructure to simulate a full tsunami, we underline the current challenges associated with this approach and review the few efforts that are currently underway to achieve this goal. A table of existing tsunami software packages is provided along with an open Github repository to allow developers and model users to update the table with additional models as they are published and help with model discoverability.

Reservoir in the Central Structural Zone of the Xihu Sag is the Huagang Formation, dominated by natural gas reservoirs, and the reservoir in the Western Slope Zone is the Pinghu Formation, dominated by light oil and wet gas. In this paper, the Gas Drive Water Displacement- Magnetic Resonance Imaging (GWD-MRI) experiments are used to simulate the charging characteristics of the sandstone migration layer. The centrifugal- Magnetic Resonance (Cen-NMR) experiments simulate the short-term fast trap charging process, and the Semi-Permeable Baffle (SPB) charging experiment simulates the slow trap accumulation process. Studies have shown that there is a start-up pressure for the migration layer charging, and the start-up pressure of the core with a permeability of 0.3mD is about 0.6MPa. Experimental simulations have confirmed that there is a front zone with changed water saturation, with a width of about 1-1.5cm. The migration layer charging mainly has two actions, drive effect and carry effect. The drive effect reduces the water saturation to 70%-80%, and the carry effect can still reduce the water saturation by 5%-10%. The water saturation of rapid charging is mainly affected by the petro-physical. The porosity is high, the water saturation is low. The water saturation decreased significantly with the increase of centrifugal force when the centrifugal force is smaller. If the centrifugal force is greater than 0.8 MPa, the water saturation decreases slowly. The water saturation of the trap slowly charging is basically maintained at 40%-50%, which matches the water saturation of the airtight coring from 40%-55%, the petro-physical do not affect the final water saturation.

The selectivity in selective autophagy pathways is achieved via the selective autophagy receptors (SARs) – proteins that bind a ligand on the substrate to be degraded and the Atg8-family protein on the growing autophagic membrane, phagophore, effectively bridging them. In mammals, the most common ligand of SARs is ubiquitin, a small protein modifier that tags substrates for their preferential degradation by autophagy. Consequently, the most common SARs are the ubiquitin-binding SARs, such as SQSTM1/p62 (sequestosome 1). Surprisingly, there is only one SAR of this type in yeast – Cue5, which acts as a receptor for aggrephagy and proteaphagy – pathways that remove the ubiquitinated protein aggregates and proteasomes, respectively. However, recent studies described the ubiquitin-dependent autophagic pathways that do not require Cue5, e.g. stationary phase lipophagy for intracellular lipid droplets and nitrogen starvation-induced mitophagy for mitochondria. What is the role of ubiquitin in these pathways? Here, we propose that the ubiquitinated lipid droplets and mitochondria are recognized by the alternative ubiquitin-binding SARs. Our analysis identifies the proteins that could potentially fulfill this role in yeast. We think that matching of the ubiquitin-dependent (but Cue5-independent) autophagic pathways with the ubiquitin-and-Atg8-binding proteins enlisted here might uncover the novel ubiquitin-binding SARs in yeast.

Abstract: Background: Hemopexin (Hx) is a plasma glycoprotein that scavenges heme (Fe(III) protoporphyrin IX), Hx have important implication in hemolytic disorders and hemorrhagic condition because the release of hemoglobin increase labile heme, which is potentially toxic producing oxidative stress. Hx has been considered for therapeutic use and diagnostics. In this work, we analyzed and mapped interaction sequences of Hx with hemin and hemoglobin (2) Methods: Spot-synthesis technique was used to map human hemopexin (P02790) binding to hemin and human hemoglobin, a library of 15 amino acid peptides with a 10-amino acid overlap was designed to represent the entire coding region (aa 1-462) of hemopexin and synthesized onto cellulose membranes. In silico approach was performed to analyze amino acid frequency in identified interaction regions, and molecular docking was applied for protein-protein interaction (3) Results: Seven linear peptide sequences in Hx were identified to bind hemin (H1-H7), and five were described for Hb (Hb1-Hb5) interaction, with just two sequences shared between hemin and Hb. Amino acid composition of identified sequences demonstrated that Histidine residues are relevant for heme binding, H105, H293, H373, H400, H429, and H462 was distributed in H1-H7 peptide sequences, but other residues may also play an important role. Molecular docking analysis demonstrated Hx association with the β-chain of Hb, with several hot spot amino acids that coordinated interaction. (4) Conclusions: This study highlights new insights on Hx-hemin binding motifs and protein-protein interactions with Hb. Binding sequences and identified specific peptides can be used for therapeutic purposes and diagnostics, as hemopexin is under investigation to treat different diseases, and there is an urgent need for diagnostics of labile heme for monitoring hemolysis.

Dictyostelium discoideum is gaining increasing attention as a model organism for the study of calcium binding and calmodulin function in basic biological events as well as human diseases. After a short overview of calcium-binding proteins, the structure of Dictyostelium calmodulin and the conformational changes effected by calcium ion binding to its four EF hands is compared to its human counterpart, emphasizing the highly conserved nature of this central regulatory protein. The calcium-dependent and -independent motifs involved in calmodulin binding to target proteins are discussed with examples of the diversity of calmodulin binding proteins that have been studied in this amoebozoan. The methods used to identify and characterize calmodulin binding proteins is covered followed by the ways Dictyostelium is currently being used as a system to study several neurodegenerative diseases and how it could serve as a model for studying calmodulinopathies such as those associated with specific types of heart arrythmia. Because of its rapid developmental cycles, its genetic tractability, and a richly endowed stock center, Dictyostelium is in a position to become a leader in the field of calmodulin research.

Faced with difficulties stemming from the complex interactions between tight gas sand bodies and fractures when describing and identifying reservoirs, a composite reservoir model was established. By setting the supply boundary to characterize the superposition characteristics of sand bodies, a mathematical model of unstable seepage in fractured vertical wells in tight sandstone gas reservoirs was constructed considering the factors such as stress sensitivity, fracture density and fracture symmetry. The seepage law and pressure response characteristics of gas well in tight sandstone discontinuous reservoir with stress sensitivity, semi-permeable supply boundary and complex fracture topology were determined, and the reliability of the model was verified. The research results more accurately depict the pressure characteristic curve of the superposed sand body complex fracture vertical well, and provide a more comprehensive model for tight gas production dynamic analysis and well test data analysis, which can more accurately guide the dynamic inversion of reservoir and fracture parameters.

In this study, a new structure is proposed based on the body-centered cubic (BCC) lattice structure: adding a cubic truss in the center of the BCC structure, which is denoted by TLC (Truss-lattice-cube). The different dimensions of the central cube can affect the mechanical properties of the lattice structure notably. With a fixed length (15mm) of a unit cell, the optimal size for the central cube is determined to be 5 mm. Quasi-static compressive tests are performed on specimens made of Polylactic Acid (PLA) using additive manufacturing technology. The deformation characteristics of the new structure is analyzed in detail by experiments and numerical simulations. Compared to the BCC structure, the mechanical properties of the TLC structure are significantly improved. The initial flow stress of the TLC has been increased by 122% at a strain of 0.1; the specific strength has been enhanced by 293% at a strain of 0.5; and the specific energy absorption has been improved by 312% at the strain of 0.6. Printing defects in the lattice structure could deteriorate the mechanical property remarkably. In this work, incorporation of microcracks into the finite element model makes the simulation capture the influence of printing defects and improve the predictive accuracy of simulation significantly.

Turbulent dissipation is a central issue in the star and galaxy formation process. Its fundamental property of space-time intermittency, well characterised in incompressible laboratory experiments, remains elusive in cosmic turbulence. Progress requires the combination of state-of-the-art modelling, numerical simulations and observations. The power of such a combination is illustrated here where the statistical method intended to locate extrema of velocity shears in a turbulent field is applied to numerical simulations of compressible magneto-hydrodynamical (MHD) turbulence dedicated to dissipation scales and to a nearby turbulent diffuse molecular cloud. We demonstrate that short-spacing increments of observables can detect strongly dissipative structures. In our simulations, we compute structure functions of various synthetic observables and show that they verify Extended Self-Similarity. This allows to compute their intermittency exponents and we show how they could help constraining some properties of the underlying three-dimensional turbulence. In observations of a turbulent cloud close to the Sun in our Galaxy, a remarkable coherent structure of velocity shear extremum is disclosed . At the location of the largest velocity shear, this coherent structure is found to foster 10 × more carbon monoxide molecules than standard diffuse molecular gas, an enrichment supported by models of non-equilibrium warm chemistry triggered by turbulent dissipation. The power of the combination between modelling and observations is also illustrated by observations of the CH+ cation that provide unique quantitative informations on the kinetic energy trail in the massive, multi-phase and turbulent circum-galactic medium of a galaxy group at redshift z = 2.8.

Shale gas and coalbed methane are energy sources that partly or even mainly consist of methane stored in an adsorbed state in the pores of the organic-rich rock and coal seams. In this study, the graphene nanoslit model is employed to model the nanometer slit pores in shale and coal. The Grand Canonical Monte Carlo (GCMC) and Molecular Dynamics (MD) modelling methods are used to investigate the mechanisms of adsorption and displacement of methane in the slit pore. It was found that the CVFF forcefield resulted in the largest adsorption amount, while the PCFF forcefield resulted in the least. The COMPASS and COMPASS II force fields led to similar results. As the width of the slit pore increases, the adsorption amount of gas molecules increases, and the number density profile of adsorbed methane molecules alters from a single adsorption layer to multi-adsorption layers. The minimum slit pore width at which methane molecules can penetrate the slit pore was found to be 0.7 nm. Moreover, it is demonstrated that by lowering the temperature, the adsorption rate of the methane increases since the adsorption is an exothermic process. Enhancing methane recovery was investigated by the injection of gases such as CO2 and N2 to displace the adsorbed methane. The comparison of adsorption isotherms of gas molecules provided the following order in terms of the amount of adsorption: CO2>CH4>N2.

This study explores wind energy resources in different locations through the Gulf of Oman and also their future variability due climate change impacts. In this regard, EC-EARTH near-surface wind outputs obtained from CORDEX-MENA simulations are used for historical and future projection of the energy. The ERA5 wind data are employed to assess the suitability of the climate model. Moreover, the ERA5 wave data over the study area are applied to compute sea surface roughness as an important variable for converting near-surface wind speeds to those of wind speed at turbine hub height. Considering the power distribution, bathymetry and distance from the coats, some spots as tentative energy hotspots to provide a detailed assessment of directional and temporal variability and also to investigate climate change impact studies. RCP8.5 is a common climatic scenario is used to project and extract future variation of the energy in the selected sites. The results of this study demonstrate that the selected locations have a suitable potential for wind power turbine plans and constructions.

This study has developed a Matlab application for simulating statistical models project (SMp) probabilistic distributions that are similar to binomial and Poisson, which were created by mathematical procedures. The simulated distributions are graphically compared with these popular distributions. The application allows to obtain many probabilistic distributions, and shows the trend (τ ) for n trials with success probability p, i.e. the maximum probability as τ=np. While the Poisson distribution PD(x;µ) is a unique probabilistic distribution, where PD=0 in x=+∞, the application simulates many SMp(x;µ,Xmax) distributions, where µ is the Poisson parameter and value of x with generally the maximum probability, and Xmax is the upper limit of x with SMp(x;µ,Xmax) ≥ 0 and limit of the stochastic region of a random discrete variable. It is shown that by simulation via, one can get many and better probabilistic distributions than by mathematical one.

This study has developed a Matlab application for simulating statistical models project (SMp) probabilistic distributions that are similar to binomial and Poisson, which were created by mathematical procedures. The simulated distributions are graphically compared with these legendary distributions. The application allows to obtain many probabilistic distributions, and shows the trend (τ ) for n trials with success probability p, i.e. the maximum probability as τ=np. While the Poisson distribution PD(x;µ) is a unique probabilistic distribution, where PD=0 in x=+∞, the application simulates many SMp(x;µ,Xmax) distributions, where µ is the Poisson parameter and value of x with generally the maximum probability, and Xmax is upper limit of x with SMp(x;µ,Xmax) ≥ 0 and limit of the stochastic region of the random discrete variable X. It is shown that by simulation via, one can get many and better probabilistic distributions than by mathematical one.

Transient analysis of gas flow in pipeline was studied. Finite Element Method based on ANSYS FLOTRAN was used to account for changes in pressure, temperature and flow rate. Compressibility factor function of temperature and pressure was considered. For non-isothermal transient results, the pressure and the flowrate gave wave propagation as a result of slow transients created by demand condition at the outlet end of custody transfer. Results obtained were in agreement with the demand restrictions at the outlet end of custody transfer indicating that the predictions are accurate and reliable. The results demonstrated that the Finite Element Method gave accurate prediction of pressure, temperature and flowrate in transient gas studies. For steady state non-isothermal model, results showed that the magnitude of the average pressure drop was higher when pressure was predicted with constant compressibility factor, but the same average pressure drop was reduced when the pressure was predicted with variable compressibility factor, z. Since compressibility factor is a function of temperature and pressure, the above findings signifies that in the case when gas temperature does not stabilize, the prediction of pressure with isothermal model and constant compressibility factor will lead to significant errors.

It is always a difficult task to assign the peaks recorded from vibrational spectrum. Herein we explored a new pathway of density functional theory (DFT) simulation to present three kinds of spectra of ice XIV and can be referenced as inelastic neutron scattering (INS), IR, and Raman experimental spectrum. The INS spectrum is proportional to phonon density of states (PDOS) while the photon scattering signals reflect the normal vibration frequencies near the Brillouin zone (BZ) center. Based on good agreements with experimental data, we identified the relative frequency and made scientific assignments by normal vibration modes analysis. The debating two hydrogen bond (H-bond) peaks among ice phases from INS have been discussed and the dynamic process of H-bond vibrations were found to be classified into two basic modes. we deduce that two H-bond modes are a general rule among ice family and more works are ongoing to investigate this subject.

Biomass gasification acquired great interest over the past decades as an effective and trustable technology to produce energy and fuels with net-zero carbon emissions. Moreover, using biomass waste as feedstock enables to recycle organic wastes and to fit the circular economy goals thus reducing the environmental impacts of waste management. Even if many studies have been already carried on, this kind of process must still be investigated and optimized with the final aim to develop industrial plants for different applications, from the hydrogen production to net-negative emission strategies. Modeling and developing of process simulations became an important tool to investigate the chemical and physical behavior of the plant, allowing to make a first raw optimization of the process and to define heat and material balances of the plant as well as to define optimal geometrical parameters with cost and time effective approaches. The present review paper focuses on the main literature models developed until now to describe biomass gasification process, and in particular on kinetic models, thermodynamic models, and computational fluid dynamic models. The aim of this study is to point out the strengths and the weakness of those models, comparing them and indicating in which situation is better to use an approach instead of another. Moreover, theoretical shortcut models and software simulations, not explicitly addressed by prior reviews, are taken into account. For researchers and designers this review provides a detailed methodology characterization as a guide to develop innovative study or project.

Background and motivation: Combining Deep Reinforcement Learning (Deep RL) and Health Systems Simulations has significant potential, for both research into improving Deep RL performance and safety, and in operational practice. While individual toolkits exist for Deep RL and Health Systems Simulations, no framework to integrate the two has been established. Aim: Provide a framework for integrating Deep RL Networks with Health System Simulations, and to ensure this framework is compatible with Deep RL agents that have been developed and tested using OpenAI Gym. Methods: We developed our framework based on the OpenAI Gym framework, and demonstrate its use on a simple hospital bed capacity model. We built the Deep RL agents using PyTorch, and the Hospital Simulation using SimPy. Results: We demonstrate example models using a Double Deep Q Network or a Duelling Double Deep Q Network as the Deep RL agent. Conclusion: SimPy may be used to create Health System Simulations that are compatible with agents developed and tested on OpenAI Gym environments. GitHub repository of code: https://github.com/MichaelAllen1966/learninghospital

Blazars are highly variable and display complex characteristics. A key characteristic is the flux distribution or flux PDF whose shape depends upon the form of the underlying physical process driving variability. The BL Lacertae, Mrk 421 is one of the brightest and most variable blazars across the electromagnetic spectrum. It has been reported to show hints of lognormality across the spectrum from observed histograms of radio to gamma-ray fluxes. This would imply that the underlying mechanisms may not conform to the "standard" additive, multi-zone picture, but could potentially have multiplicative processes. This is investigated by testing the observed lightcurves at different wavelengths with time-series simulations. We find that the simulations reveal a more complex scenario, than a single lognormal distribution explaining the multiwavelength lightcurves of Mrk 421.

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. There is evidence that some ligands that bind to the GPCRs activate different downstream signaling pathways. G protein activation or -arrestin biased signaling involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. Molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions between receptors and ligands. Here, we highlight recent simulation studies and methods used to study biased G protein-coupled receptor signaling and their conformational dynamics. We also highlight applications of MD simulations to drug discovery.

The log interpretation and evaluation of Chang 8 tight sandstone reservoir in Longdong West area, Ordos Basin, is facing great challenges due to the existence of low resistivity contrast oil pays. To better guide the exploration of oil resources in this area, the characteristics and control mechanism of low resistivity contrast oil pays were studied in this research. Firstly, according to the relative value of apparent resistivity increase rate of the target formation, the reservoir was divided into two types: low resistivity contrast oil pay (LRCP) and high resistivity oil pay (HRP). Then, the reservoir characteristics were studied by comparing and analyzing the experimental data, formation water data and logging data collected from the two reservoir types. On this basis, the control mechanism of LRCP was studied and summarized into reservoir micro factors and regional macro factors, respectively. It is found that the reservoir rock composition between LRCP and HRP is basically the same. Compared with HRP reservoir, the average porosity and permeability of LRCP reservoir are relatively better, but the pore structure is relatively poorer because of the high content of micro pores. When the filling pressure of crude oil is sufficient enough, the high oil saturation can be formed in LRCP reservoir. The high irreducible water saturation and high formation water salinity are the main micro factors reduced the reservoir resistivity. Besides, the difference of hydrocarbon expulsion capacity of source rock and the regional difference of formation water salinity are the macro factors influenced the distribution of LRCP and HRP from vertical and horizontal of the region. The formation of LRCP is controlled by the comprehensive influence of reservoir micro factors and regional macro factors. And the comprehensive consideration of the influence of micro factors and macro factors on LRCP is suggested during the log interpretation and evaluation.

In an attempt to correct misunderstandings this article brings together the observations on Calcium, Myosin Binding Protein-C and Hypertrophic Cardiomyopathy in the basic function of cardiac muscle. A finding of many years ago is reiterated in a novel enzyme kinetic format with defined rate limiting step which makes CaATP the apparent substrate of the actomyosin cross-bridge. The relationship of these kinetics to recent observations on disruption of myosin binding protein-C is described along with how this bears on the understanding of the related cardiomyopathies.

At nuclear scale, we present three heuristic relations pertaining to strong and electroweak coupling constants. With these relations, close to beta stability line, it is possible to study nuclear binding energy with a single energy coefficient of magnitude $\left(\frac{1}{{\alpha }_s}\right)\left[\frac{e^2}{4\pi {\epsilon }_0{R}_0}\right]\approx 10.0\text{MeV}\text{.}$ With reference to up and down quark masses, it is also possible to interpret that, nuclear binding energy is proportional to the mean mass of $\left[\left(2m_u+m_d\right)\text{and}\left(m_u+2m_d\right)\right]\approx 10.0\text{MeV}.$

Computational simulations of three-dimensional flow around a NACA 0018 wing with an aspect ratio of AR=5 are carried out using the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations with the Shear-Stress Transport turbulence model closure. Simulations were performed to capture the aerodynamic stall hysteresis using a developed Pseudo-Transient Continuation (PTC) method based on a dual time step approach in the CFD code OpenFOAM. The flow was characterized by an incompressible Mach number of M=0.12 and a moderate Reynolds number Re=0.67×10^6. The results obtained indicate the presence of a noticeable aerodynamic hysteresis in the static dependencies of the force and moment coefficients, as well as the manifestation of bi-stable flow separation patterns, accompanied by the development of asymmetry in the stall zone. The URANS simulation results are in good agreement with the experimental data obtained for the NACA 0018 finite aspect ratio wing from the low-speed wind tunnel under the same test conditions. A new phenomenological bifurcation model of aerodynamic stall hysteresis under static and dynamic conditions is formulated, which is proven to be able to closely match the experimental data.

Fires are complex multi-physics problems that span wide spatial scale ranges. Capturing this complexity in computationally affordable numerical simulations for process studies and “outer-loop” techniques (e.g., optimization and uncertainty quantification) is a fundamental challenge in reacting flow research. Further complications arise for propagating fires where a priori knowledge of the fire spread rate and direction is typically not available. In such cases, static mesh refinement at all possible fire locations is a computationally inefficient approach to bridging the wide range of spatial scales relevant to fire behavior. In the present study, we address this challenge by incorporating adaptive mesh refinement (AMR) in fireFoam, an OpenFOAM solver for simulations of complex fire phenomena involving pyrolyzing solid surfaces. The AMR functionality in the extended solver, called fireDyMFoam, is load balanced, models gas, solid, and liquid phases, and allows us to dynamically track regions of interest, thus avoiding inefficient over-resolution of areas far from a propagating flame. We demonstrate the AMR capability and computational efficiency for fire spread on vertical panels, showing that the AMR solver reproduces results obtained using much larger statically refined meshes, but at a substantially reduced computational cost. We then leverage the computational efficiency of the AMR solver to demonstrate an optimization framework for fire suppression based on the open-source Dakota toolkit.

Introduction: This study aims to elucidate the clinical diagnostic value of S100 Calcium Binding Protein beta (S100B) in epilepsy, a topic of ongoing debate among European and Asian epilepsy research centres.Methods and analysis: A structured long-term research plan was developed, comprising three phases. Phase one involves observational studies investigating the correlation between S100B levels in serum or cerebrospinal fluid and epilepsy. Phase two integrates observational and randomised controlled data to establish baseline S100B levels across diverse epilepsy demographics, underpinning the predictive value of S100B for epilepsy. Phase three focuses on randomised controlled trials to evaluate the diagnostic accuracy of S100B as an epilepsy biomarker.Ethics and dissemination: All enrolled patients will provide written consent. The ethics committee of the Second People’s Hospital of Gansu Province, Affiliated Hospital of Northwest Minzu University, China, approved the study (GSSEY2022-KY009-01).Registration: The study is registered with PROSPERO (CRD42023425431).

Intermolecular interactions are the fabrics underlying almost all processes in living organisms, where two cornerstone concepts, intermolecular binding affinity (K d ) and binding energy (ΔG), have long been established to physically describe the strengths of biomolecular interactions, e.g., drug-target K d and ΔG to describe the strength of drug-target interaction. The past two-three years saw a big step forward in the use of artificial intelligence (AI) in structural biology (e.g., AlphaFold for protein structure prediction) and drug discovery & design. In light of the roles of K d and ΔG in drug discovery & design, the speed of this AI progress raises a question of what’s next for its practical application in the pharmaceutical industry, in addition to a system-wide account of biomolecular structures and motions. Last August, the concept of a general intermolecular binding affinity calculator (GIBAC) was for the first time coined and proposed in an MDPI-published preprint. Here, this article puts forward an updated conceptual and practical framework of GIBAC, including its inception, definition, construction, practical applications, technical challenges and limitations, and future directions. Moreover, this article argues that the time is now ripe for the construction of such an accurate, precise and efficient GIBAC to be on the agenda of the entire drug discovery & design community, to ensure its applicability & reliability, and to enhance its value in drug R&D in future.

The brain’s ability to create a unified conscious representation of an object by integrating information from multiple perception pathways is called perceptual binding. Binding is crucial for normal cognitive function. Some perceptual binding errors and disorders have been linked to certain neurological conditions, brain lesions, and conditions that give rise to illusory conjunctions. However, the mechanism of perceptual binding remains elusive. Here, we present a computational model of binding using two sets of coupled oscillatory processes that are assumed to occur in response to two different percepts. We use the model to study the dynamic behavior of coupled processes to characterize how these processes can modulate each other and reach a temporal synchrony. We identify different oscillatory dynamic regimes that depend on coupling mechanisms and parameter values. The model can also discriminate different combinations of initial inputs that are set by initial states of coupled processes. Decoding brain signals that are formed through perceptual binding is a challenging task, but our modeling results demonstrate how crosstalk between two systems of processes can possibly modulate their outputs. Therefore, our mechanistic model can help one gain a better understanding of how crosstalk between perception pathways can affect the dynamic behavior of the systems that involve perceptual binding.

Drilling operation is the major cost factor for the industry. Appropriately designed operations are essential for successful drilling. Optimized drilling operations also allow enhanced drilling performance and reduce the drilling cost. This is achieved by increasing the bit life (minimizing premature bit wear), drilling faster that allows reducing drilling time, and also reducing tripping operations. The paper presents two parts. The first part compares the parametric physic-based swab and surge simulation results obtained from Bingham Plastic, Power Law, and Robertson and Stiff models. The aim is to show how the model predictions deviate from each other. Two 80-20 Oil Water ratios (OWR) and two 90-10 OWR oil-based drilling fluids with 1.96 sg and 2.0 sg were considered in vertical- and deviated wells. Simulation results analysis revealed that the deviations depend on the drilling fluid’s physical and rheological parameters as well as the well trajectory. Moreover, the model’s predictions were inconsistent. Data-driven machine learning (ML) modeling makes up the second section. Data-driven modeling was done using both software-generated datasets and field datasets. Results show that the Random Forest Regressor (RF), Artificial Neural Network (ANN), Long Short-Term Memory (LSTM), LightGBM, XGboost, and Multivariate regression models predict the training and test datasets with higher R-squared and minimum root means square values. Deploying the ML model in real-time applications and the planning phase would have the potential application of artificial intelligence for well planning and optimization processes.

Large office buildings are responsible for a substantial portion of energy consumption in urban districts. However, thorough assessments regarding the Nordic countries are still lacking. In this paper we analyse the largest dataset to date for a Nordic office building, by considering a case study located in Stockholm, Sweden, that is occupied by nearly a thousand employees.Distinguishing the lighting and occupants’ appliances energy use from heating and cooling, we can estimate the impact of occupancy without any schedule data. A standard frequentist analysis is compared with Bayesian inference, and the according regression formulas are listed in tables that are easy to implement into building performance simulations (BPS). Monthly as well as seasonal correlations are addressed, showing the critical importance of occupancy.A simple method, grounded on the power drain measurements aimed at generating boundary conditions for the BPS, is also introduced; it shows how, for this type of data and number of occupants, no more complexities are needed in order to obtain reliable predictions. For an average year, we overestimate the measured cumulative consumption by only 4.7%. The model can be easily generalised to a variety of datasets.

Deficit schizophrenia is characterized by leaky tight and adherens junctions and bacterial translocation. Here we examine whether (deficit) schizophrenia is accompanied by leaky paracellular, transcellular and vascular barriers in the gut and blood brain barriers. We measured IgA responses to occludin, claudin-5, E-cadherin and β-catenin (paracellular pathway, PARA), talin, actin, vinculin and epithelial intermediate filament (transcellular pathway, TRANS) and plasmalemma vesicle-associated protein (PLVAP, vascular pathway) in 78 schizophrenia patients and 40 controls. IgA responses to claudin-5, E-cadherin and β-catenin, the sum of the four PARA proteins and the ratio PARA/TRANS were significantly higher in deficit schizophrenia than in non-deficit schizophrenia and controls. A large part of the variance in PHEMN (psychosis, hostility, excitation, mannerism and negative) symptoms, psychomotor retardation, formal thought disorders, verbal fluency, word list memory, word list recall and executive functions was explained by the PARA/TRANS ratio coupled with plasma IgA responses to Gram-negative bacteria, IgM to malondialdehyde, CCL-11 (eotaxin), IgA levels of the ratio of noxious to more protective tryptophan catabolites (NOX/PRO TRYCATs) and a plasma immune activation index. Moreover, IgA levels to Gram-negative bacteria were significantly associated with IgA to E-cadherin, β-catenin and PLVAP, while IgA levels to claudin-5 were significantly predicted by IgA to E-cadherin, NOX/PRO TRYCAT ratio, Gram-negative bacteria and CCL11. The phenomenology of the deficit syndrome is to a large extent explained by the cumulative effects of lowered natural IgM, breakdown of the paracellular and vascular pathways, increased bacterial translocation, peripheral immune-inflammatory responses and indices of BBB breakdown.

Ubiquitin binding domains (UBDs) are modular elements that bind non-covalently to ubiquitin and act as downstream effectors and amplifiers of the ubiquitination signal. With few exceptions, UBDs recognize the hydrophobic path centered on Ile44 (Leu-8, Ile-44, Val-70). Nevertheless, a variety of different orientations, which can be attributed to specific contacts between each UBD and surface residues surrounding the hydrophobic patch, specify how each class of UBD recognizes ubiquitin. Here, we describe the structure of a novel ubiquitin-binding domain that we identified in NEDD4 binding protein 1 (N4BP1). By performing protein sequence analysis, mutagenesis and NMR spectroscopy of the ¹⁵N isotopically labelled protein, we demonstrate that a Phe-Pro motif in N4BP1 recognizes the canonical hydrophobic patch of ubiquitin. This recognition mode resembles the molecular mechanism evolved in the CUE (Coupling of ubiquitin conjugation to ER degradation) domain family, where an invariant proline, usually following a phenylalanine, is required for binding to ubiquitin. Interestingly, the UBD of N4BP1 is evolutionary related to CUBAN (Cullin binding domain associating with NEDD8) (40% identity and 47% similarity), a protein module that also recognizes the ubiquitin-like NEDD8, which is the closest relative of ubiquitin (58% identity and 80% similarity). By performing circular dichroism and ¹⁵N NMR chemical shift perturbation of N4BP1 in complex with ubiquitin, we demonstrate that the UBD of N4BP1 lacks the NEDD8 binding properties observed in CUBAN and it recognizes the Ile44-centered patch of ubiquitin through a dedicated binding site, which share some of the features observed in the CUE domain family. Moreover, we show that, in addition to mediating the interaction with ubiquitin and ubiquitinated substrates, both the CUBAN and the UBD of N4BP1 are poly-ubiquitinated in cells. This modification is dependent on the presence of a functional domain. We believe that the structural and functional characterization of this novel UBD will allow a deeper understanding of the molecular mechanisms governing N4BP1 function, while at the same time providing a valuable tool for clarifying how the discrimination between ubiquitin and the highly related NEDD8 is achieved.

Initially, retinol-binding protein (RBP), was thought to be a biomarker for proximal convoluted tubule dysfunction could be important in chronic kidney diseases (CKD). Membranous glomerulonephritis (MGN) is an important cause of CKD and end-stage renal disease (ESRD). Therefore, monitoring MGN patients using urinary RBP is important in effective treatment decision making and prognostication of MGN patients. Enzyme-linked immunosorbent assay (ELISA) technique was used to detect the RBP in the urine samples of 69 MGN patients comprising 47 primary and 22 secondary MGN, at the end of the follow-up period. The test for the urinary biomarker gave the following results: urinary RBP was detected in 27 (39.1%) and 6 (8.7%) of the primary and secondary MGN patients, respectively. The correlation analysis demonstrated a significant relationship between urinary RBP and renal function test parameters, in addition to a logistic regression analysis that proved urinary RBP as a prognostic non-invasive biomarker for primary MGN. Therefore, urinary RBP could be employed to monitor and provide effective prognosis and early treatment decisions in primary MGN.

Physics behind collimated highly directional nature of lasers, and factors that keep the seven coloured waves that form white light together during their journey from Sun to Earth, in the face of the natural disruptive forces, is not fully understood. Energy levels were measured, in terms of alterations in induced current and voltage, in beams from a red laser, white LED light and the Sunlight before and during their disruption by diffusers (frosted glass for the lasers) and diffractors (diffraction grating for the white light) using a photovoltaic solar cell panel attached to a digital multimeter. Results show that disruption of the beams results in release of extra energy named as ‘Latent Light Binding’ Energy’. It is hypothesized that the ‘binding’ energy keeps laser waves firmly bound together both end-on and side-on enabling laser beams to travel long distances in collimated manner. Likewise, the 7 coloured waves that constitute white light are kept together, probably side-on, in their journey from the Sun to the Earth. The observation that diffraction of sunbeam is associated with increased power generation provides a new lead to improve harnessing of solar energy, where, currently, the focus is mainly on improving efficiency of photovoltaic cell.

A pivotal strategy in decarbonizing the transport sector is by transitioning from internal combustion engine vehicles (ICEVs) to electric vehicles (EVs), especially if they are charged with renewable energy. To provide an incentive to the common populace to switch to EVs, cities should set up - or enable or even encourage the private sector to do so - electric vehicle charging stations (EVCS) at multiple and suitably accessible locations. This paper aims to present a digitaltwin-enabled methodology study for selecting the most suitable possible locations to establish EVCS with consideration of the low voltage (LV) network operation in terms of power losses, active and reactive power in EV charging demand and the effect of the EVCS on the transformers, based on examples in Newcastle upon Tyne (England, UK). It also evaluates the most suitable type of charger that could be installed at these locations by studying the estimated traffic, cost and impact of the EVCS on the power grid, both at peak and off hours. This paper demonstrates how digitaltwin technologies, such as the web-based tool AutoDesign by Northern Powergrid and EA Technology and open source Pandapower can be used for identifying the most viable and cost-effective locations to connect EVCSs to a Distribution Network Operators (DNO)’s LV network. The results of the analysis indicate that with the installation of the most suitable type of EVCS in a particular location, the local power grid would not be significantly affected at any hour of the day.

3D bioprinting, like 3D printing, is a process that builds structures by depositing material layer by layer. In the specific case of bioprinting, these layers are composed of a biocompatible material together with living cells, allowing the creation of three-dimensional structures that recreate functional tissues and organs. Bioprinting cell-laden structures is complicated by the high rate of cell damage and stress that occurs during the process, caused by the high pressures and stresses to which they are subjected. To minimise damage during the process, it is important to study and optimise certain bioprinting conditions beforehand and to analyse how they affect the cells by means of computational simulations. In this review, both quantitative and qualitative data are collected to improve the nozzle geometry through computational simulation studies. Optimal ranges for nozzle diameter (0.2-1mm) and length (8-10mm; 300-900μm) have been defined, and recommendations for improving nozzle performance during bioprinting, such as the provision of an inner angle of 20-30°, an inner EDTA coating and a shear stress of less than 10kPa, have been gathered. Finally, based on the data collected, a design of experiments is proposed to obtain an optimal bioprinting configuration for a particular bioink.

This study investigates the magnetization of Cobalt (Co) nanoparticles with two different crystal structures: a centered cubic (fcc) and an hexagonal close packed (hcp), under the magnetic field of the Earth. For that reason, the hysteresis loops of the two material samples are exported and compared with the corresponding value by experimental procedure in the laboratory under the Earth’s magnetization. For both samples, the nanoparticle number variation and the nearest neighbor distance effect, in the Co hysteresis loop are investigated, structural information that is useful for the sample preparation. For the simulations the Finite Element Method (FEM) was used, while experiments are conducted using a measurement set-up based on fluxgate magnetometers. From the simulation, it is proven that, for both samples, the nanoparticle number effect does not change the hysteresis loop. In contrast, when the distance between the nanoparticles increases, the magnetization obtains smaller value. This study aims to build a novel, simple and economical measurement arrangement enabling in situ magnetic measurements for ferromagnetic nano-material characterization. As the need for continuous knowledge of the condition and the phase of ferromagnetic nanoparticles that is crucial in many applications, like catalytic processes, this work aims to introduce an alternative approach that will add another dimension to conventional chemical reaction engineering this of the magnetic state of the material resolved in time and space and to contribute in the field of the monitoring and the early fault diagnosis of nanoparticle materials.

Air pollution is responsible for any adverse effects on human beings. Thermal discomfort, on the other hand, is able to overload the human body and eventually provoke health implications due to the heat imbalance. Methods: The aim of the present work is to study the behavior of two bio-climatic indexes and statistical characteristics of the air quality index for Sofia city - the capital of Bulgaria for the period 2008 - 2014. The study is based on WRF-CMAQ model system simula-tions with a spatial resolution of 1 kilometer. The air quality is estimated by the air quality index, taking into account the influence of different pollutants and the thermal conditions by two indi-ces, respectively, for hot and cold weather. Results: It was found that half of the heat and cold index categories are present in the simulations. Their distribution has some spatial features. All air quality categories are present in the domain, with dominance only of the O3 and PM10. Conclu-sions: It was found that Sofia is not so hot and air polluted place, but in some situations, people have to have some concerns when intend to be outdoors for a prolonged time.

The availability of high quality surface observations of precipitation and volume observations by polarimetric operational radars make it possible to constrain, evaluate and validate numerical models with a wide variety of microphysical schemes. In this article, a novel particle-based Monte-Carlo microphysical model (called “McSnow”) is used to simulate the outer rain bands of Hurricane Dorian which traversed the densely instrumented precipitation research facility operated by NASA at Wallops Island, Virginia. The rain bands showed steady stratiform vertical profiles with radar signature of dendritic growth layers near −15 °C and peak reflectivity in the bright band of 55 dBZ along with polarimetric signatures of wet snow with sizes inferred to exceed 15 mm. A 2D-video disdrometer measured frequent occurrences of large drops >5 mm and combined with an optical array probe the drop size distribution was well-documented in spite of uncertainty for drops <0.5 mm due to high wind gusts and turbulence. The 1D McSnow control run and four numerical “experiments” were conducted and compared with observations. One of the main findings is that even at the moderate rain rate of 10 mm/h collisional breakup is essential for the shape of the drop size distribution.

In this manuscript, an automated framework dedicated to design space exploration and design optimization studies is presented. The framework integrates a set of numerical simulation, computer-aided design, numerical optimization, and data analytics tools using scripting capabilities. The tools used are open-source and freeware, and can be deployed on any platform. The main feature of the proposed methodology is the use of a cloud-based parametrical computer-aided design application, which allows the user to change any parametric variable defined in the solid model. We demonstrate the capabilities and flexibility of the framework using computational fluid dynamics applications; however, the same workflow can be used with any numerical simulation tool (e.g., a structural solver or a spread-sheet) that is able to interact via a command line interface or using scripting languages. We conduct design space exploration and design optimization studies using quantitative and qualitative metrics, and to reduce the high computing times and computational resources intrinsic to these kinds of studies, concurrent simulations and surrogate-based optimization are used.

The present studies deliver the computational investigations of a 10 MW turbine with a diameter of 205.8 m developed within the framework of the AVATAR (Advanced Aerodynamic Tools for Large Rotors) project. The simulations were carried out using two methods with different fidelity levels, namely the computational fluid dynamics (CFD) and blade element and momentum (BEM) approaches. For this purpose, a new BEM code namely B-GO was developed employing several correction terms and three different polar and spatial interpolation options. Several flow conditions were considered in the simulations, ranging from the design condition to the off-design condition where massive flow separation takes place, challenging the validity of the BEM approach. An excellent agreement is obtained between the BEM computations and the 3D CFD results for all blade regions, even when massive flow separation occurs on the blade inboard area. The results demonstrate that the selection of the polar data can influence the accuracy of the BEM results significantly, where the 3D polar datasets extracted from the CFD simulations are considered the best. The BEM prediction depends on the interpolation order and the blade segment discretization.

The total enthalpies of the 16 different spin configurations that can be realized in the unit cell of the archetype spin crossover complex [Fe(phen)₂(NCS)₂] (phen=1,2-phenanthroline) were calculated applying periodic density functional theory combined with the Hubbard model and the Grimme-D2 dispersion correction (DFT+U/D2). The obtained enthalpy differences between the individual spin configurations were used to determine spin couplings of an Ising-like model and subsequent Monte Carlo simulations for this model allowed to estimate the phenomenological interaction parameter Г of the Slichter-Drickamer model, which is commonly used to describe the cooperativity of the spin transition. The calculational procedure described here, which led to an estimate of about 3 kJ/mol⁻¹ for Г, in good agreement with experiment, may be used to predict from first principles how modifications of spin crossover complexes can change the character of the spin transition from gradual to abrupt and vice versa.

This study explores the impact of multi-user interaction, group work, and simulations within virtual worlds on language learning and the development of intergroup empathy. Two distinct research projects were conducted, involving a total of 241 participants aged 19-45. The language learning study engaged 116 participants in diverse interactive experiences, while the intercultural study had 125 participants collaborating in multicultural groups and participating in perspective-taking simulations. Both studies employed qualitative data collection methods, including surveys, interviews, and observations. The findings suggest that the combination of networking strategies, collaborative learning, and simulations within virtual worlds contributes to improvements in learners' language proficiency, confidence, and empathy towards diverse social groups. Participants reported increased motivation and engagement, attributed to the immersive and interactive nature of the virtual environments. The studies highlight the importance of collaboration and reflection in facilitating language acquisition and intercultural understanding. Technical challenges were identified as potential barriers to implementation. The results demonstrate the potential of virtual worlds to enhance language education and foster empathy in diverse societies, offering valuable insights for educators and researchers. However, the findings may be limited by the specific contexts and sample sizes of the studies, warranting further research to explore the generalizability and long-term impact of virtual world interventions. not exaggerate the main conclusions.

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this mini-review, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions and, on other hand, highly simplified structure-based Elastic Network Models of protein flexibility. In contrast to classical Molecular Dynamics the modeling strategies discussed here allow quite accurate modeling of much larger systems and longer time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Point defects induced by doping rare earth elements RE (Nd, Er) into the magnesium oxide host were investigated by classical atomistic simulations utilising the General Utility Lattice Program (GULP). Formation and association energies were calculated for the potential defect structures. Both isolated defects and defect complexes were considered. The most energetically favourable structures of defect complexes were found for rare earth doped and Li co-doped systems. The correlation between the association energy and the structure of the defect complex was investigated. The influences of Li were revealed with respect to energy and structure. The simulation results contribute to the understanding of the point defects of doped MgO and how Li influences the doping of rare earth elements in the MgO host.

In a previous companion study, satellite data were used to describe peculiar characteristics of ocean surface wave fields, generated by two extra-tropical cyclones (ETCs) rapidly propagating in the North Atlantic. Based on a 2D parametric wave model, further details are now provided to analyse and interpret the spatio-temporal evolution of ETC very intense generated waves. Significant wave height and wavelength values are shown to reach extreme values, 18 m and 500 m, respectively. Resulting energetic swell systems waves then radiate in the whole eastern part of North Atlantic, and more particularly in the Norwegian sea region. Moving to higher latitudes, wind forcing characteristics of ETCs evolve, with the shape of the wind field changing from quasi-cyclonic to "air jets/Icelandic lows". In this paper, the resulting swell generation and propagation, after the deformation of an individual ETC, are studied, as well. Confirmed with comparisons with multi-satellite observations, the application of the parametric-2D wave-ray model is demonstrated to provide robust and highly detailed information on wave generation under very complex wind regime changes.

Formation in soccer is among the most important tactical choices for a successful match.Herein, the simulations of 420000 match-plays have been performed varying the formation, the number of opponents that are actively pressing the team, the speed of the opponents in attempting a pass interception. Dribbling has been neglected. The match-play ends either with a successful series of passes from a central back to the line of the strikers or with the opponents that steal the ball. In this work, I demonstrate that 3-4-3 formation, which is among the most employed formations, relates to the highest probability of success.

In this work, a tool for estimating the contact angle from the molecular dynamics simulations is developed and presented. The tool (Achilles) can detect water droplet on hydrophobic and hydrophilic surfaces. The tool can reconstruct the droplets broken across the periodic boundaries. Further a neighbor density based accurate filter is used to find the droplet liquid vapor interface and a circle is fitted using it after removing the dense layers of water next to solid surface. This fitted circle is solved for contact angle and results are outputted in the form of graphical images and text. The entire content of the internal computations of the tool is broken down into 4 phases and users can monitor the outcomes at every phase through output images. The tool is tested using sample molecular dynamics results of water droplet on hydrophobic and hydrophilic surfaces. We believe this tool can be a good addition to the molecular dynamics simulation community who work on the interfacial physics, droplet evaporation, super hydrophobic surfaces, and wettability etc.

A common reactor type in the chemical and process industry is the fixed-bed reactor. Accurate modeling can be achieved with particle-resolved Computational Fluid Dynamic (CFD) simulations. However, the underlying bed morphology plays a paramount role. Synthetic bed-generation methods are much more flexible and faster than image-based approaches. In this study, we look critically at the two different bed generation methods: Discrete Element Method (DEM) (in the commercial software STAR-CCM+) and the rigid-body model (in the open-source software Blender). The two approaches are compared in terms of synthetically generated beds with experimental data of overall and radial porosity, particle orientation, as well as radial velocities. Both models show accurate agreement for the porosity. However, only Blender shows similar particle orientation than the experimental results. The main drawback of the DEM is the long calculation time and the shape approximation with composite particles.

Computational Fluid Dynamics (CFD) simulations, based on Reynolds Averaged Navier Stokes (RANS) models, are a useful tool for a wide range of coastal and offshore applications, providing a high fidelity representation of the underlying hydrodynamic processes. Generating input waves in the CFD simulation is performed by a numerical wavemaker (NWM), with a variety of different NWM methods existing for this task. While NWMs, based on impulse source methods, have been widely applied for wave generation in depth averaged, shallow water models, they have not seen the same level of adoption in the more general RANS based CFD simulations, due to difficulties in relating the required impulse source function to the resulting free surface elevation for non-shallow water cases. This paper presents an implementation of an impulse source wavemaker, which is able to self-calibrate the impulse source function to produce a desired wave series in deep or shallow water at a specific point in time and space. Example applications are presented, for a numerical wave tank (NWT), based on the opensource CFD software OpenFOAM, for wave packets in deep and shallow water, highlighting the correct calibration of phase and amplitude. Also, the suitability for cases requiring very low reflection from NWT boundaries is demonstrated. Possible issues in the use of the method are discussed and guidance for good application is given.

β-dicalcium silicate (β-Ca₂SiO₄, or β-C₂S in cement chemistry notation) is one of the most important minerals in cement. An improvement of its hydration rate would be the key point for developing environmentally friendly cements with lower energy consumption and CO₂ emissions. However, there is a lack of fundamental understanding on the water/β-C₂S surface interactions. In this work we aim to evaluate the water adsorption on three β-C₂S surfaces at the atomic scale using density functional theory (DFT) calculations. Our results indicate that thermodynamically favorable water adsorption takes place in several surface sites, with a broad range of adsorption energies (−0.78 to −1.48 eV), depending on the particular mineral surface and adsorption site. To clarify the key factor governing the adsorption, the electronic properties of water at the surface were analyzed. The partial density of states (DOS), charge analysis, and electron density difference analyses suggest a dual interaction of water with β-C₂S (100) surface: a nucleophilic interaction of the water oxygen lone pair with surface calcium atoms, and an electrophilic interaction (hydrogen bond) of one water hydrogen with surface oxygen atoms. Despite the elucidation of the adsorption mechanism, no correlation was found between the electronic structure and the adsorption energies.

The catalytic parameters describing enzyme catalysis are discussed in the context of the design of de novo enzymes. It is shown that high values of the catalytic efficiency may in some cases reflect tight substrate binding and lead to an immediate levelling-off of rates to reach a low plateau, an undesirable scenario for many industrial applications. It is argued that increasing the catalytic rate constant should be a main target in de novo enzyme engineering.

Zinc plays a crucial role in the process of virion maturation inside the host cell. The accessory Cys-rich proteins expressed in SARS-CoV-2 by genes ORF7a and ORF8 are likely involved in zinc binding and in interactions with cellular antigens activated by extensive disulfide bonds. In this report we provide a proof of concept for the feasibility of a structural study of orf7a and orf8 proteins. A conceivable hypothesis is that lack of cellular zinc, or substitution thereof, might lead to a significant slowing down of viral maturation.

Affordances are an important concept in cognition, which can be applied to robots in order to perform a successful human-robot interaction (HRI). In this paper we explore and discuss the idea of emotional affordances and propose a viable model for implementation into HRI. We consider “2-ways” affordances: perceived object triggering an emotion, and perceived human emotion expression triggering an action. In order to make the implementation generic, the proposed model includes a library that can be customised depending on the specific robot and application’s scenario. We present the AAA (Affordance-Appraisal-Arousal) model, which incorporates Plutchik’s Wheel of Emotions, and show some examples of simulation and possible scenarios.

Ice-binding proteins (IBPs) have several functions that permit their hosts to thrive in the presence of ice. The ability of IBPs to control ice growth makes them potential additives in various industries ranging from food storage and cryopreservation to anti-icing systems. For IBPs to be used in commercial applications, however, methods are needed to produce sufficient quantities of high-quality proteins. Here, we describe a new method for IBP purification, termed falling water ice purification (FWIP). The method is based on the affinity of IBPs for ice. A crude IBP solution is allowed to flow continuously over the large chilled vertical surface of a commercial ice machine. The temperature of the surface is lowered gradually until ice crystals are produced, to which the IBPs bind but other solutes do not. As in other ice affinity methods, FWIP does not require molecular tags and is suitable for purifying recombinant IBPs as well as IBPs from natural sources. The advantage of FWIP over other ice affinity methods is that it exploits an ice machine designed to produce large volumes of clear ice daily. This system can be easily scaled up and suits the purification of industrial quantities of IBPs. The FWIP method significantly advances the use of IBPs in research and industry.

Ferric uptake regulation protein is a repressor protein which binds an AT rich region of DNA (the iron box). Fur binds as a dimer in a helix turn helix mode and it is activated by iron(II) and other divalent transition metal ions at elevated concentrations in a process to regulate the ion uptake. Each transition metal ion induces certain conformational changes to aid the Fur binding, both the N-terminal and C-terminal domains take part in binding to DNA in addition to His 88 and His 86 residues. The process is discussed in view of experimental reports. Fe(II), Mn(II) and Co(II) activate Fur to bind DNA experimentally but Zinc plays a structural role and does not activate Fur to bind DNA.

A new approach was presented in this study to simulate a stable atom along with subatomic particles in it. Subatomic particles in an atom are treated as conventional particles, exhibiting attraction and repulsion between them, giving a stable atom simulation. This new approach leads to developing new potentials for subatomic particles to simulate a stable atom. Subatomic particle dynamics (SPD) simulations were developed for within-atom simulations that give a stable atom structure for infinite time and provide many insights into electron trajectories and atom properties like energies and radius of the atom. The effect of protons and neutrons movement in the nucleus on electron trajectories was captured in SPD. Hydrogen and Carbon atom SPD were reported in this study. The Pilot-wave theory was implemented in this study to simulate the wave nature of subatomic particles in an atom. Electrons motion was guided by the interference pattern produced by the electron’s aether medium waves in SPD. SPD was implemented on molecules, resulting in a stable molecular structure giving insights into electron trajectories shared by multiple atoms.

In this study, we introduce the dynamics of a Hepatitis B virus (HBV) model with the class of asymptomatic carriers and conduct a comprehensive analysis to explore its theoretical aspects and examine the crossover effect within the HBV model. To investigate the crossover behavior of the operators, we divide the study interval into two subintervals. In the first interval, the classical derivative is employed to study the qualitative properties of the proposed system, while in the second interval, we utilize the ABC fractional differential operator. Consequently, the study is initiated using the piecewise Atangana-Baleanu derivative framework for the systems. The HBV model is then analyzed to determine the existence, Hyers-Ulam (HU) stability, and disease-free equilibrium point of the model. Moreover, we showcase the application of an Adams-type predictor-corrector (PC) technique for Atangana-Baleanu derivatives and an extended Adams-Bashforth-Moulton (ABM) method for Caputo derivatives through numerical results. Subsequently, we employ computational methods to numerically solve the models and visually present the obtained outcomes using different fractional-order values. This network is designed to provide more precise information for disease modeling, considering that communities often interact with one another, and the rate of disease spread is influenced by this factor.

The successive steps of the oxidation mechanism of crocin, major compound of saffron, by the free OH• radical are investigated by the pulse radiolysis, steady-state (gamma) radiolysis methods and molecular simulations. The optical absorption properties of the transient species and their reaction rate constants are determined. The absorption spectrum of the oxidized radical of crocin resulting from the H-abstraction presents a maximum at 678 nm and a band at 441 nm, almost as intense as that of crocin. The spectrum of the covalent dimer of this radical contains an intense band at 441 nm and a weaker one at 330 nm. The final oxidized crocin, issued from radical disproportionation, absorbs weaker with a maximum at 330 nm. The molecular simulation results suggest that the OH• radical is electrostatically attracted by the terminal sugar and is scavenged predominantly by the neighbor methyl site of the polyene chain as in a sugar-driven mechanism. Based on detailed experimental and theoretical investigations, the antioxidant properties of crocin are highlighted.

Molecular Dynamics (MD) simulations model motion of molecules in atomistic detail and aid in drug design. While simulations on large systems may require several days to complete, analysis of terabytes of data generated in the process could also be time consuming. Recent studies captured exciting and dramatic drug-receptor interactions under cell-like complex conditions. Such advances make simulations of biomolecular interactions more realistic, insightful, and informative and have potential to make drug design more realistic. However, currently available resources and techniques do not provide, in reasonable time, a comprehensive understanding of events seen in simulations. We demonstrate that big data approach results in significant speedups, and provides rapid insights into simulations performed. Advancing this improvement, we propose a scalable, self-tuning, and responsive framework based on Cloud-infrastructure to accomplish the best possible MD studies with given priorities and within available resources.

The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling disordered protein interactions and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein-peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein-peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution and studies of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.

The combined efect of shape and electrostatic complementarities (Sc, EC) at the interface of the interacting protein partners (PPI) serves as the physical basis for such associations and is a strong detterminant of their binding energetics. EnCPdock (https://www.scinetmol.in/EnCPdock/) presents a comprehensive web-platform for the direct conjoint comparative analyses of complementarity and binding energetics in PPIs. It elegantly interlinks the dual nature of local (Sc) and non-local complementarity (EC) in PPIs using the Complementarity Plot. It further derives an AI-based ∆G_binding with a prediction accuracy comparable to the state-of-the-art. This book-chapter presents a practical manual to conceptualize and implement EnCPdock with its various features and functionalities, collectively having the potential to serve as a valuable protein engineering tool in the design of novel protein interfaces.