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.

Rice blast and brown planthopper (BPH) are two of the most destructive diseases and pests in rice (Oryza sativa L.) production in China. The development of multi-resistant varieties is widely recognized as the most efficient and environmentally friendly approach to controlling crop diseases and pests. Functional molecular markers (FMMs) have been developed from functional variants in the genic region associated with trait variation, greatly enhancing the efficiency of identifying and pyramiding valuable genes in crop breeding. In this study, two FMMs and a multiplex PCR amplification system were developed for two major broad-spectrum BPH resistance genes, Bph6 and Bph9. By assistance of FMMs in multi-resistant rice restorer line development pipeline, two lines (Huahui7713 and Huahui3006) with blast and BPH resistance were developed by pyramiding three resistance genes Pigm, Bph6 and Bph9. Huahui7713 and Huahui3006 as male lines, two new hybrid rice varieties (WeiLiangYou7713 and Xuanliangyou3006) were certificated and commercial released in China. WeiLiangYou7713 and Xuanliangyou3006 showed improved blast and BPH resistance, and no penalty on yield and grain-quality. Annual promotion area of WeiLiangYou7713 has exceeded 53.3 thousand hectares, indicating its widespread adoption and success in the market. Huahui7713 and Huahui3006 could have great potential for utilization in blast and BPH resistance rice breeding.

Familial Alzheimer’s disease (FAD) can be caused by mutations in PSEN1 that encodes presenilin-1, a component of the gamma-secretase complex that cleaves amyloid precursor protein. Alterations in calcium (Ca2+) homeostasis and glutamate signaling are implicated in the pathogenesis of FAD, however it has been difficult to assess in humans whether or not these phenotypes are the result of amyloid or tau pathology. This study aimed to assess early calcium and glutamate phenotypes of FAD by measuring the Ca2+ response of induced pluripotent stem cell (iPSC)-derived neurons bearing PSEN1 mutations to glutamate and the ionotropic glutamate receptor agonists NMDA, AMPA and kainate compared to isogenic control and healthy lines. The data shows that in early neurons, even in the absence of amyloid and tau phenotypes, FAD neurons exhibit increased Ca2+ responses to glutamate and AMPA, but not NMDA or kainate. Together this suggests that PSEN1 mutations alter Ca2+ and glutamate signaling as an early phenotype of FAD.

This work presents a novel approach for the study of the movement of droplets on inclined surfaces under the influence of gravity and chemical heterogeneities. The developed numerical methodology uses data-driven modeling to extend the applicability limits of an analytically derived reduced-order model for the contact line velocity. More specifically, while the reduced-order model is able to capture the effects of the chemical heterogeneities to a satisfactory degree, it does not account for gravity. To alleviate this shortcoming, datasets generated from direct numerical simulations are used to train a data-driven model for the contact line velocity, which is based on the Fourier neural operator and corrects the reduced-order model predictions to match the reference solutions. This hybrid surrogate model, which comprises of both analytical and data-driven components, is then integrated in time to simulate the droplet movement, offering a speedup of five orders of magnitude compared to direct numerical simulations. The performance of this hybrid model is quantified and assessed in different wetting scenarios, by considering various inclination angles and values for the Bond number, demonstrating the accuracy of the predictions as long as the adopted parameters lie within the ranges considered in the training dataset.

As the most active plateau on the Earth, the Qinghai-Tibet Plateau has a complex crust-mantle structure. Knowledge of the distribution of such a structure provides information for understanding the underlying geodynamic processes. We obtains a three-dimensional density model of crustal and upper mantle beneath Qinghai-Tibet plateau and its surrounding areas from the residual geoid anomalies using the Earth Gravitational Model (EGM) 2008. We estimate a refined density model by iterations, using an initial density contrast model. We confirm that the EGM2008 mission products can be used to constrain the crust-mantle density structures. Our major findings are: (1). At 300-400 km depth, high-D anomalies terminate around Jinsha River Suture (JRS) in the central TP, suggesting that the Indian plate has been reached over the Bangong Nujiang Suture (BNS) and almost reach to the JRS. (2). On the eastern TP, low-D anomalies at the depth of 0-300 km together with high-D anomalies at 400-670 km further verified the current eastward subduction of Indian plate. The ongoing subduction provides forces to the occurrences of frequent earthquakes and volcano. (3). At 600 km depth, low-D anomalies inside the TP illustrate the existence of hot weak material beneath there, contributing to the external material inward-thrusting.

Music has the innate potential to reach all parts of the brain, stimulates certain brain areas which are not achievable through other modalities. Music Therapy (MT) is being used for more than a century to treat individuals who needs personalized care. MT optimizes motor, speech and language responsibilities of the brain and improves cognitive performance. Pervasive developmentdisorder (PDD) is a multifaceted, neuro developmental disorder and autism spectrum disorder (ASD) comes under PDD, which is defined by deficiencies in three principal spheres: social connection with others, communicative and normal movement skills. The conventional imaging studies illustrate reduced brain area connectivity in people with ASD, involving selected parts of the brain cortex. People with ASD express much interest in musical activities which engages the brain network areas and improves communication and social skills.The main objective of this review is to analyze the potential role of MT in treating the neurological conditions, particularly ASD. Evidence based studies have reported the extensive therapeutic application of music on various part of the brain in a nonverbal child with autism through hearing or making music.Hence we hypothesized that MT intervention can improve the communication capacity in people with ASD, than customary neurorestoration therapy alone.

Fibroin-chitosan composites, especially those containing nanohydroxyapatite, show potential for bone tissue regeneration. The physicochemical properties of these biocomposites depend on the compatibility between their components. In this study, the intermolecular interactions of fibroin and chitosan were analyzed using a molecular dynamics approach. Two types of systems were investigated: one containing acetic acid and the other containing calcium (Ca²⁺) and hydrogen phosphate (HPO₄²⁻) ions mimicking hydroxyapatite conditions. After obtaining the optimal equilibrium structures, the distributions of several types of interactions including hydrogen bonds, ionic contacts, and hydrophobic contacts along with the structural and energetical features were examined. The calculated binding energy values for the fibroin-chitosan complexes confirm their remarkable stability. The high affinity of fibroin for chitosan can be explained by the formation of a dense network of interactions between considered biopolymers. These interactions were found to be primarily hydrogen bonds and ionic contacts involving ALA, ARG, ASN, ASP, GLN, GLU, GLY, LEU, PRO, SER, THR, TYR, and VAL residues. As it was established, the complexation of fibroin with chitosan maintains the β-sheet conformation of the peptide. β-Sheet fragments in fibroin are involved in the formation of significant amount of hydrogen bonds and ionic contacts with chitosan.

In this paper, the boiling flow inside a nanochannel with 700000 argon particle has been simulated by molecular dynamics (MD) simulation. This approach has been employed to analysis the superheated flow and its heat transfer pattern as well. For all simulations an external thrust force varying from 1 PN to 12 PN is exerted on inlet nanoparticles along the channel to have the forced annular boiling flow. Computations reveal that saturation condition and superheat degree have significant impacts on the liquid-vapor interface. Furthermore, because of the major influence of surface tension throughout a nanochannel, the x-velocity of liquid film and vapor core has not considerable fluctuations and stay smooth. All provided results show the behaviors completely similar to the available outcomes in the literature.

Circulating tumor cells (CTCs) are cells released from the primary tumor and intravasate into the blood or lymphatic vessels, where they are transported to distant sites and act as seeds that initiate cancer metastases. Recent advances in CTC research have shown their relevance as prognostic markers for early and metastatic disease detection, predictive biomarkers for relapse, and response to medical intervention or therapy. The rapidly evolving landscape of CTC biology has opened new avenues for understanding cancer progression, metastasis, and treatment response. Additionally, translating these findings into clinical applications holds promise for improving cancer diagnostics, prognosis, and personalized therapeutic strategies. This review discusses the significance of CTCs in cancer research and their associated challenges. We explore recent developments in the detection and characterization of CTCs and their implications in cancer research and clinical practice.

The BET family of proteins is characterized by the presence of two tandem bromodomains and an extra-terminal domain. The mammalian BET family of proteins comprises Brd2, Brd3, Brd4, and Brdt, which are encoded by paralogous genes that may have been generated by repeated duplication of an ancestral gene during evolution. Bromodomains that can specifically bind acetylated lysine residues in histones serve as chromatin-targeting modules that decipher the histone acetylation code. BET proteins play a crucial role in regulating gene transcription through epigenetic interactions between bromodomains and acetylated histones during cellular proliferation and differentiation processes. On the other hand, BET proteins have been reported to mediate latent viral infection in host cells and be involved in oncogenesis. Human BRD4 is involved in multiple processes of the DNA virus life cycle, including viral replication, genome maintenance, and gene transcription through interaction with viral proteins. Aberrant BRD4 expression contributes to carcinogenesis by mediating hyperacetylation of the chromatin containing the cell proliferation-promoting genes. BET bromodomain blockade using small-molecule inhibitors gives rise to selective repression of the transcriptional network driven by c-Myc. These inhibitors are expected to be potential therapeutic drugs for a wide range of cancers. This review presents an overview of the basic roles of BET proteins and highlights the pathological functions of BET and the recent developments in cancer therapy targeting BET proteins in animal models.

The information of seismic damage of buildings in SAR images of different time phase, especially in SAR images after earthquake, is easily disturbed by other factors, which affects the accuracy of information discrimination. In order to identify and evaluate the distribution information of the seismic damage accurately and make full use of the abundant texture features in the SAR image. The conventional method of change detection based on texture features usually takes the pixel as the calculating unit. In this paper, a method of texture feature change detection of SAR images based on watershed segmentation algorithm is proposed. Based on the optimization of texture feature parameters, the feature parameters are segmented by the watershed segmentation algorithm, and the feature object image is obtained. This method introduces the idea of object oriented, and carries out the calculation of the difference map at the object level, Finally, the classification threshold value of different types of seismic damage types is selected, and the recognition of building damage is achieved. Taking the ALOS data before and after the earthquake in Yushu as an example to verify the effectiveness of the method, the overall accuracy of the building extraction is 88.9%, Compared with pixel-based methods, it is proved that the proposed method is effective.

Reconstructing the Tree of Life remains a central goal in biology. Early methods, relying on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the Tree of Life, two recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in two or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. We identify and discuss specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the Tree of Life.

The p53 tumor suppressor protein shows antiviral activity. The replication of viruses is also antagonized by interferons (IFNs) – cytokines, which control immunity genes by the STAT transcription factors. The best studied interferons belong to type I (e.g. IFNα) and type II (IFNγ) groups. IFNα and IFNγ both induce activating phosphorylation of STAT1 on Tyr701. Previously, we noticed that p53 activates SOCS1, which is a negative regulator of STAT1 phosphorylation. Thus, we hypothesized that p53 by activating SOCS1 expression reduces phosphorylation of STAT1 and attenuates activation of genes stimulated either by IFNα or IFNγ. To test this hypothesis, we exposed p53-proficient and p53-deficient cells to strong p53 activators, and either to IFNα or IFNγ. Subsequently, we estimated the phosphorylation status of STAT1 and the expression of genes regulated by either of these cytokines. Strong activation of p53 reduced the phosphorylation of STAT1 on Tyr701, however surprisingly it did not lead to reduced expression of most tested interferon-stimulated genes. On contrary, p53 and IFNγ synergized in the activation of CASP1, IFIT1 and IFIT3 genes. We conclude that the interactions between p53 and interferon-activated pathways are more complicated than initially expected, but the apparent cooperation between p53 and IFNγ deserves further studies.

Oligopeptide transporters (OPTs) are a group of membrane localized proteins that have a broad range of substrate transport capabilities and contribute to numerous biological processes. However, limited information has been reported on OPTs in higher plants. In this study, a comprehensive analysis of the OPT gene family in Medicago truncatula was performed. A total of 26 OPT genes (MtOPT01-MtOPT26) have been identified in the Medicago truncatula genome. Phylogenetic analyses indicated that MtOPTs consisted of two distinct subgroups, 12 MtOPTs belonged to the peptide transport subgroup (PT-OPT) based on their predicted amino acid sequences containing the two highly conserved motifs (NPG and KIPPR) and 14 MtOPTs belonged to yellow stripe subgroup (YS-OPT). The MtOPTs distributed on each of 8 chromosomes in Medicago truncatula. Sequence analysis verified that MtOPTs significant similar to those in other plants. The copy number of MtOPTs was low and the multiply of MtOPTs was simple relatively. Gene structure analysis showed that most of the MtOPTs have various numbers of introns. The multiple of MtOPTs might play different biological roles which were supported by the fact that MtOPTs have a distinct tissue-specific expression pattern. The data obtained in this study will help to better understand the complexity of the MtOPTs gene family and provide new evidence for the function and evolution of the OPT gene family in higher plants.

Non-coding RNAs have been described as crucial regulators of gene expression and guards for cellular homeostasis. Some recent papers are focused on vault RNAs, one of the classes of non-coding RNA, and their role in cell proliferation, tumorigenesis, apoptosis, cancer response to therapy, and autophagy, which made them potential therapy targets in oncology. In the human genome, four vault RNA paralogues can be distinguished. They are associated with vault complexes, considered the largest ribonucleoprotein complex. The protein part of these complexes consists of a major vault protein (MVP) and two minor vault proteins (vPARP and TEP1). The name of the complex, as well as vault RNA, comes from the hollow barrel-shaped structure that reminds a vault. Their sequence and structure are highly evolutionary conserved and show many similarities in comparison with different species, but vault RNAs show various roles. Vaults were discovered in 1986, and their functions remained unclear for many years. Although not much is known about their contribution to cell metabolism, it has become clear that vault RNAs are involved in various processes and pathways associated with cancer progression and modulating cell functioning in normal and pathological stages. In this review, we discuss known functions of the human vault RNAs in the context of cellular metabolism, emphasizing processes related to cancer and cancer therapy efficacy.

From their discovery in 2002, BRAF mutations have been identified as clear drivers of oncogenesis in several cancer types. Currently, their incidence rate is nearly 7% of all solid tumors with BRAF V600E constituting approximately 90% of these diagnoses [1]. In melanoma, thyroid cancer, and histiocytic neoplasms, BRAF hotspot mutations are found at a rate of about 50%, while in lung and colorectal cancers they range from 3 to 10% of reported cases [2-5]. Though present in other malignancies such as breast and ovarian cancers, they constitute a small portion of diagnoses (<1%). Given their frequency along with advancements in screening technologies, various methods are used for the detection of BRAF-mutant cancers. Among these are targeted next-generation sequencing (NGS) on tumor tissue or circulating tumor DNA (ctDNA) and immunohistochemistry (IHC)-based assays. With advancements in detection technologies, several approaches have been taken for the treatment of BRAF mutant cancers. In this review we retrace the milestones that led to the clinical development of targeted therapies currently available for these tumors.

Huperzine A is an herbal reversible inhibitor of Acetylcholinesterase (AChE). A molecular docking analysis on Huperzine A molecule has been carried out to understand its structure, conformational flexibility, intermolecular interaction and the binding affinity in the active site of AChE enzyme. Further, the charge density distribution of huperzine A molecule (lifted from the active site of AChE) was determined from the high level quantum chemical calculations coupled with charge density analysis. The binding affinity of Huperzine A towards AChE was calculated from the molecular docking; the lowest docked energy is -8.46 kcal/mol. In the active site, huperzine A molecule interacts with acyl binding pocket-Phe330 of AChE, that is, the bicyclo ring group of huperzine A forms an intermolecular interaction with the oxygen atom of main chain of the amino acid residue Phe330 at the distances 3.02 and 3.25 Å respectively. On the other hand, a gas phase study on huperzine A molecule also performed using HF and DFT (B3LYP) methods with the basis set 6-311G**. The molecular structure, conformation, and the charge density distribution of huperzine A molecule in the gas phase have determined using quantum chemical calculations and the charge density analysis. The comparative studies between the gas phase and the active site forms of huperzine A molecule, explicitly reveals the degree of conformational modification and the charge density redistribution of huperzine A when present in the active site. The dipole moment of the molecule in the active site is 6.85 D, which is slightly higher than its gas phase value (5.91 D). The electrostatic potential (ESP) surface of active site molecule clearly shows the strong electronegative and positive ESP regions of the molecule, which are the expected strong reactive locations of the molecule.

Food systems will need to undergo considerable transformation. To be better prepared for and resilient to uncertainty and disturbances in the future, resource users and managers need to further develop knowledge about the food and farming system, with its dominating feedback structures and complexities, and to test robust and integrated system-based solutions. This paper investigates how participatory system dynamics modeling can be adapted to groups at the community level with low or no formal educational background. The paper also analyzes the refinement of workshop participants’ mental models as a consequence of a participatory system dynamics intervention. For this purpose, we ran two workshops with small-scale farmers in Zambia. Analysis of workshop data and post-workshop interviews shows that participatory system dynamics is well adaptable to support an audience-specific learning-by-doing approach. The use of pictures, objects and water glasses in combination with the basic aspects of causal loop diagramming makes for a well-balanced toolbox. Participants acquire understanding that is also relevant beyond systems thinking in that is offers a range of practical insights such as a critical evaluation of common food security strategies.

More than 20 years after the discovery of the first anelloviruses (AVs), they still remain a mysterious group of viruses. To date, there is no clear understanding of their impact on the host organism, although information is emerging about their participation in various pathologies. However, it is well known that AVs are an important component of the human virome. In the current study we assessed the abundance of AVs in viromes from four nasal and pharyngeal samples from children with respiratory symptoms using a metagenomic approach. It is worth noting that three patients had known respiratory viruses and one did not have any. Regarding AVs, they were detected in all samples, namely Torque teno mini virus (TTMV). The greatest abundance and diversity of AVs (Torque teno virus (TTV), TTMV and Torque teno midi virus (TTMDV)) were shown in a diseased patient with an absence of respiratory viruses, which may indicate the potential involvement of AVs in the maintenance of respiratory symptoms. An activated immune system to combat the main pathogens, against the background of which the number of AVs decreases, is discussed as a possible reason for the reduced diversity of anellome in other patients. We also compared the results of the metagenomic approach and qPCR for major respiratory viruses, which were in agreement. In addition, metagenomic sequencing made it possible to obtain additional data on viral genomes in order to establish their taxonomic identification and characterize individual viral profiles in patients from the Siberian part of Russia.

Enveloped viruses represent a significant category of pathogens that cause serious diseases in animals. These viruses express envelope glycoproteins that are singularly important during infection of host cells by mediating fusion between the viral envelope and host cell membranes. Despite low homology at protein levels, three classes of viral fusion proteins have, as of yet, been identified based on structural similarities. Their incorporation into viral particles is dependent upon their proper sub-cellular localization after being expressed and folded properly in the endoplasmic reticulum (ER). However, viral protein expression can cause stress in the ER, and host cells respond to alleviate the ER stress in the form of the unfolded protein response (UPR); the effects of which have been observed potentiating or inhibiting viral infection. One important arm of UPR is to elevate the capacity of the ER-associated protein degradation (ERAD) pathway, which is comprised of host quality control machinery that ensures proper protein folding. In this review, we provide relevant details regarding viral envelope glycoproteins, UPR, ERAD, and their interactions in host cells.

Silencing of the fragile X mental retardation 1 (FMR1) gene and consequently lack of synthesis of FMR protein (FMRP) are associated with fragile X syndrome, which is one of the most prevalent inherited intellectual disabilities. FMRP is a multifunctional protein involved in many cellular functions in nearly all subcellular compartments under normal conditions and under conditions of cellular stress in both neuronal and non-neuronal cell types. This is achieved through its trafficking signals, nuclear localization signal (NLS), nuclear export signal (NES), and nucleolar localization signal (NoLS), as well as its RNA and protein binding domains, and is modulated by various post-translational modifications such as phosphorylation, ubiquitination, sumoylation, and methylation. This review summarizes recent advances in understanding the interaction networks of FMRP with a special focus on FMRP stress-related functions, including stress granule formation, mitochondrion and endoplasmic reticulum plasticity, ribosome biogenesis, cell cycle control, and DNA damage response.

A large diversity of epigenetic factors, such as microRNAs and histones modifications, are known to be capable of regulating gene expression without altering DNA sequence itself. In particular, miR-1 is considered the first essential microRNA in cardiac development. In this study, miR-1 potential role in early cardiac chamber differentiation is analyzed through specific signaling pathways. For this, we performed in chick embryos functional experiments by means of miR-1 microinjections into the posterior cardiac precursors -of both primitive endocardial tubes- committed to sinoatrial region fates. Subsequently, embryos were subjected to whole mount in situ hybridization, immunohistochemistry and RT-qPCR analysis. As a relevant nov-elty, our results revealed that miR-1 increases Tbx5, Gata4 and AMHC1, while this microRNA diminishes Mef2c expression, during early differentiation of cardiac sinoatrial region. Further-more, we observed in this developmental context that miR-1 upregulates CRABPII and RARß, and downregulates CRABPI, which are three crucial factors in retinoic acid signaling pathway. Interestingly, we also noticed that miR-1 directly interacts with HDAC4 and Calmodulin, as well as with Erk2/MAPK1, which are three key factors actively involved in Mef2c regulation. Our study shows, for the first time, a key role of miR-1 as an epigenetic regulator in early differenti-ation of cardiac sinoatrial region through orchestrating opposite actions between retinoic acid and Mef2c, fundamental to properly assign cardiac cells to their respective heart chambers. A better understanding of those molecular mechanisms modulated by miR-1 will definitely help in fields applied to therapy and cardiac regeneration and repair.

Some peptide-based drugs including oxytocin, vasopressin, ziconotide, pramlintide, nesiritide, and octreotide, contain one intramolecular disulfide bonds. A novel and reusable monodispersed silica nanosphere-supported Pt(IV) complex (SiO2@TPEA@Pt(IV)) was synthesized via a four-step procedure and used for formation of intramolecular disulfide bonds in peptides. TEM and chemical mapping results for the Pt(II) intermediates and for SiO2@TPEA@Pt(IV) show that the silica nanospheres possess a monodisperse spherical structure and contain uniformly distributed Si, O, C, N, Cl, and Pt. The valence state of Pt on the silica nanospheres was characterized by XPS. The Pt(IV) loaded on SiO2@TPEA@Pt(IV) was 0.15 mmol/g, as determined by UV-vis spectrometry. Formation of intramolecular disulfides in six dithiol-containing peptides of variable lengths by use of SiO2@TPEA@Pt(IV) was investigated and the relative oxidation yields were determined by HPLC. In addition, peptide 1 (Ac-CPFC-NH2) was utilized to study the reusability of SiO2@TPEA@Pt(IV). No significant decrease in the relative oxidation yield was observed after ten reaction cycles. Moreover, the structure of SiO2@TPEA@Pt(IV) after being used ten cycles was determined to be similar to its initial one, demonstrating the cycling stability of the complex.

Many fundamental molecular techniques (PCR, Microarray, Southern and northern hybridization, siRNA, CRISPR/Cas9 etc.) developed so far shows errors. I wish to highlight these molecular techniques are developed on basis of Watson-Crick DNA model, ignoring the concept of parallel stranded DNA. Through this opinion article, I wish to highlight specificity and accuracy of these molecular techniques can be enhanced by considering both parallel and anti parallel hybridization of DNA. Hopefully my views will also solve issue of irreproducibility in life science research.

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.

Systems Chemistry has been a key component of origin of life research, invoking models of life’s inception based on evolving molecular networks. One such model is the Graded Autocatalysis Replication Domain (GARD) formalism embodied in a Lipid World scenario, which offers rigorous computer simulation based on defined chemical kinetics equations. GARD suggests that the first pre-RNA life-like entities could have been homeostatically-growing assemblies of amphiphiles, undergoing compositional replication and mutations, as well as rudimentary selection and evolution. Recent progress in Molecular Dynamics has provided an experimental tool to study complex biological phenomena such as protein folding, ligand-receptor interactions and micellar formation, growth and fission. The detailed molecular definition of GARD and its inter-molecular catalytic interactions make it highly compatible with Molecular Dynamics analyses. We present a roadmap for simulating GARD’s kinetic and thermodynamic behavior using various Molecular Dynamics methodologies. We review different approaches for testing the validity of the GARD model, by following micellar accretion and fission events and examining compositional changes over time. Near future computational advances could provide empirical delineation for further system complexification, from simple compositional non-covalent assemblies towards more life-like protocellular entities with covalent chemistry that underlies metabolism and genetic encoding.

Nanoporous membranes promise energy-efficient water desalination. Hexagonal boron nitride (h-BN), like graphene, exhibits outstanding physical and chemical properties, making it a promising candidate for water treatment. We employed Car–Parrinello molecular dynamics simulations to establish accurate modeling of Na+ and Cl- permeation through hydrogen passivated nanopores in graphene and h-BN membranes. In these complex systems the classical potentials cannot properly account for the physics at play. We demonstrate that ion separation works well for the h-BN by imposing a barrier of 0.13 eV and 0.24 eV for Na+ and Cl- permeation respectively, compared to the graphene where the Cl- ion faces a negative barrier of 0.68 eV and is prone toward blockade and Na+ permeation is associated with a slightly negative barrier of 0.03 eV. From the change in the solvation shell, we get the following trend: Na/h-BN > Na/G > Cl/G > Cl/h-BN. We argue that the trend in the permeation barrier changes to Cl-h-BN > Na-h-BN > Na-G > Cl-G due to a combination of several interactions, including the distortion of the water network induced by the ions, the ion-water interaction, and the ion-nanopore interaction. Overall, the desalination performance of h-BN surpasses that of their graphene counterparts.

The p53 tumor suppressor protein is an activator of transcription. Diverse stress factors lead to various sets of posttranslational modifications of p53 what results in different sets of upregulat-ed genes. We noticed that actinomycin D and nutlin-3a (A+N) synergize in inducing activating phosphorylations of p53 and upregulation of selected p53-target genes. Here we found that one of these genes is DUSP13, which codes for poorly-studied, dual-specificity phosphatase having at least two isoforms, one expressed in testis and the other in skeletal muscles. In cancer cells ex-posed to A+N, DUSP13 is expressed from an alternative promoter in intron, what results in ex-pression of isoform named TMDP-L1. The luciferase reporter tests demonstrated that this pro-moter is activated by both endogenous and ectopically expressed p53. We showed for the first time that mRNA expressed from this promoter actually produces the protein, which can be de-tected by Western blotting in all examined cancer cell lines with wild-type p53 exposed to A+N. In some cell lines it is also induced by clinically relevant camptothecin or by nutlin-3a acting alone. This isoform, fused with green fluorescent protein localizes in perinuclear region of cells. Thus, TMDP-L1 isoform may be an important element of p53-regulated stress response system.

This study focuses on the potential medicinal properties of Campsis radicans L. extracts and its aim is to investigate the antioxidant and antiproliferative effects of these extracts on the human colon adenocarcinoma cell line (HT-29). Cytotoxic activity in the HT-29 cell line was evaluated, and the 3-dimensional structures of phenolic compounds were obtained using the GC-MS method. Docking simulations were performed with AutoDock Vina, targeting proteins that are crucial in the apoptotic pathway (Bax, Bcl-2, and Caspase-3). Drug similarities and ADMET estimates of the phenolic compounds were obtained from databases. The highest antioxidant activity of Campsis radicans L. plant was observed in ethanol and water extracts, with notable inhibition capacities. While the total phenolic content was highest in ethanol, water, and ether extracts respectively the total flavonoid content was highest in ethanol, then ether, and water extracts. Certain elements were undetectable in mineral determination, and others exceeded the toxic doses for plants. The water extract exhibited dose-dependent cytotoxic activity in the HT-29 cell line, with an IC50 value of 4.47 mg/mL. Molecular docking results revealed the binding of phenolic molecules to the active sites of apoptotic pathway proteins. The findings suggest that Campsis radicans L. plant possesses cytotoxic properties, significant phenolic and flavonoid content, and a high antioxidant capacity. These properties make it a potential candidate for treatment use, emphasizing its potential therapeutic value in the HT-29 cell line. Further research and exploration on the identified phenolic compounds could contribute to the development of novel therapeutic agents.

Transition-metal-free synthesis of 4-pyrones via TfOH-promoted nucleophilic addition/cyclization of diynones and water has been developed. This transformation is simple, atom economical and environmentally benign, providing rapid and efficient access to substituted 4-pyrones.

Aurora kinase A (AURKA) is a normal cell proliferation-inducing enzyme encoded by AURKA gene, with over-expression observed in different types of malignancies. Hence, the goal is to find potential inhibitors against AURKA. In this study, molecular docking, Standard Precision and Extra Precision methods were employed. After the docking study, the ligands showed an extremely low binding score which suggested very high binding affinity of the ligands. Furthermore, Quantum polarized ligand docking (QPLD) was performed to predict the binding status of the molecules. Based on the binding affinity, the top four compounds were chosen for further analysis. The docked complexes were further analyzed in explicit water conditions using 100 ns molecular dynamics simulations and binding free energy calculation. Then, density functional theory (DFT) calculation was used to calculate the molecular properties of the molecules. Finally, systems biology experiments validated the molecular docking and molecular dynamics simulation studies and indicated that quercetin, kaempferol, luteolin and rutin could inhibit the AURKA. The results show that, these four molecules have high binding affinity to the AURKA and significant interactions (LEU139, GLU211and ALA213) were also identified with the hinge region of Aurora kinase A. Thus, LEU139, GLU211, and ALA213 were identified as the crucial protein mechanisms.

Recurrent pregnancy loss (RPL) affects around 2 % of women of reproductive age. Primary RPL is defined by ≥2 pregnancy losses and no normal birth delivery. In secondary RPL, the losses are after a normal pregnancy and delivery. Most cases have no clear aetiology, although primary cases are the most complex. Several gene polymorphisms have been associated with RPL. Some polymorphisms are increased in women of Asian or Caucasian race. In admixed RPL populations, the information on possible genetic links is scarce and contradictory. The study aimed to assess the frequency of two SNP in different enzymes involved in medical conditions observed during pregnancy. Microsomal epoxy hydrolase (mEPH) is involved in detoxifying xenobiotics. It is present in ovaries, and it is hormonally regulated. The endothelial nitric oxide synthase forms nitric oxide from arginine and is involved in vascular tone. Two SNPs, rs1051740 and rs1799983, were assessed. The study included 26 controls and 63 primary RPL patients. Significant differences were observed in the frequencies of both SNP (p<0.05). The frequency of mutated alleles was higher in patients. Double mutated homozygotes were encountered only in RLP patients (p<0.05). Genetic polymorphisms rs1051740 and rs1799983 may be involved in primary RPL in the Venezuelan admix population. Genetic studies could provide crucial information on the aetiology of primary RPL.

The African continent has the highest prevalence of hypertension globally, with South Africa reporting the highest prevalence in Southern Africa. While the influence of genetic variability in the pathogenesis of hypertension is well described internationally, limited reports are available for African populations. This study aimed to assess the association of genetic variants and essential hypertension in a cohort of two ethnic South African population groups. Two hundred and seventy-seven hypertensive and one hundred and seventy-six normotensive individuals were genotyped for 79 variants. Genotyping was performed using the Illumina GoldenGate Assay and allele-specific polymerase chain reaction. The association of variants was assessed using the Fisher Exact test under the additive and allelic genetic models, while multivariate logistic regression was used to predict the development of hypertension. Five variants (CYP11B2 rs179998, AGT rs5051 and rs699, AGTR1 rs5186 and ACE rs4646994) were significantly associated with essential hypertension in the cohort under study. Furthermore, AGTR1 rs5186 and AGT rs699 were identified as risk factors for the development of hypertension in both ethnic groups. In two ethnic South African populations an association was observed between Renin-Angiotensin-Aldosterone System (RAAS) related genes and the development of hypertension.

This article develops a multi-perspective view on motivations and methods for tobamovirus purification through the ages, and presents a novel efficient, easy-to-use approach that can be well-adapted to different species of native and functionalized virions. We survey the various driving forces prompting researchers to enrich tobamoviruses: from the search for the causative agents of mosaic diseases in plants to their increasing recognition as versatile nanocarriers in biomedical and engineering applications. Best practice and also rarely applied options for the serial processing steps required for successful isolation of tobamoviruses are then reviewed. Adaptations for distinct particle species, pitfalls and 'forgotten' or underrepresented technologies are considered as well. The article is topped off with our own development of a method for virion preparation, rooted in historical protocols. It combines selective re-solubilization of polyethylene glycol (PEG) virion raw precipitates with density step gradient centrifugation in biocompatible iodixanol formulations, yielding ready-to-use particle suspensions. This newly established protocol and some considerations for perhaps worthwile further developments could serve as putative stepping stones towards preparation procedures appropriate for routine practical uses of the multivalent soft-matter nanorods.

Chronic kidney disease patients treated by hemodialysis present a high cardio-vascular morbidity and mortality. There is an imperative need for novel biomarkers in identifying and to offer possible therapeutically intervention in these patients. We performed a prospective observational cohort study on 77 patients in the period October 2021 – October 2023. We measured serum plasma levels of Interleukin 1-Beta, Galectin 3, human suppression of tumorigenicity factor 2, bone morphogenetic protein 2 and fibroblastic growth factor 23 at the inclusion. We evaluated the correlations of these biomarkers with cardiac function and structure evaluated by echocardiography. Mean age was 61.02 (±11.81) years, 45 (56.2%) males and with a dialysis vintage of 4.95 (2.4-7.8) years. Median ejection fraction was 51 (43-54 %) and more than two thirds of the patients presented valvular calcifications. Overall mortality was 22%. Interleukin 1-Beta correlated positively with ejection fraction and global longitudinal strain and negatively with left atrium diameter and left ventricle tele systolic diameter. Galectin 3 values negatively correlated with aortic valve fibrosis and mitral valve calcifications and human suppression tumorigenity factor 2 negatively correlated with mitral valve calcifications. Some of these novel biomarkers could be used to better asses the cardio-vascular disease in patients on maintenance hemodialysis.

Iron is an essential nutrient and necessary for biological functions from DNA replication and repair to transcriptional regulation, mitochondrial respiration, electron transfer, oxygen transport, photosynthesis, enzymatic catalysis, and nitrogen fixation. However, due to iron’s propensity to generate toxic radicals which can cause damage to DNA, proteins, and lipids, multiple processes regulate uptake and distribution of iron in living systems. Understanding how intracellular iron metabolism is optimized and how iron is utilized to regulate other intracellular processes is important to our overall understanding of a multitude of biological processes. One of the tools that the cell utilizes to regulate a multitude of functions is ligation of the iron-sulfur (Fe-S) cluster cofactor. Fe-S clusters comprised of iron and inorganic sulfur are ancient components of living matter on earth that are integral for physiological function in all domains of life. Here we will explore the ways in which the cell utilizes Fe-S clusters to sense the intracellular environment and respond to maintain equilibrium.

Shigella and Salmonella have been implicated in neonatal sepsis, food poisoning, and bacteremia. Their increasing resistance to antibiotics has undermined treatment options, constituting a public health crisis. Also, the scarcity of new antibiotics in the development pipeline compounds multidrug-resistance (MDR). Hence, there is a need to potentiate antibiotics to mitigate resistance. This study leveraged riboflavin to re-sensitize fluoroquinolones and carbapenems against MDR Shigella and Salmonella spp. Disk diffusion, microbroth-dilution, chequerboard, and time-kill assays were used to determine resistance profiles and riboflavin-antibiotic activity. Biofilm-formation and bacterial viability (after riboflavin-antibiotic challenge) was determined with crystal-violet and Macrophage infection assays, respectively. Strains exhibited &gt;80% resistance to fluoroquinolones and carbapenems, with 160–320 µg MIC ranges and 0.5–2 FIC, indicating no synergy. At 40–80 µg riboflavin, there was two-fold reduction in bacterial survival, with meropenem-riboflavin exhibiting significant reduction and FIC of 0.075 indicating high synergistic activity. Biofilm-formation was altered (from strong to weak phenotype) with 40 µg riboflavin-antibiotic treatment. There was a 40–70% reduction in bacterial survival in THP-1 cells treated with meropenem-riboflavin after 3 hours, and complete clearance between 24-48 hpi. Riboflavin potentiated fluoroquinolones and meropenem, leading to modulation of the host immune response for the clearance of the strains.

Carotid atherosclerotic lesions are a major cause of cerebrovascular disease (CVD). Identification and quantification of carotid plaques are important for categorizing the vulnerability of plaques for rupture and assessing the impact of treatments. The irregularity of plaque surface is associated with previous plaque rupture and plays an important role in the risk of stroke. Thus, the aim of this study is to develop and validate novel vulnerability biomarkers from three-dimensional ultrasound (3DUS) images by analyzing the surface morphological characterization of carotid plaque using fractal geometry features. 3D box-counting and 3D blanket are the two types of 3D fractal dimension that were employed to describe the smoothness of plaques. This fractal dimension analysis tool was used to evaluate the effect of atorvastatin using 3DUS carotid images, which were acquired from 6 patients treated with atorvastatin with 80 mg daily and 5 patients with placebo. The Student's T Test results showed that those two fractal features were effective for detecting the statin-related changes in carotid atherosclerosis with p<0.0068 and p<0.015 respectively, suggesting that 3D fractal dimension measurements can be used effectively to analyze the surface characteristics of carotid plaques, especially for evaluating the impact of the treatment.

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.

Studies performed at single-cell resolution have demonstrated the physiological significance of cell-to-cell variability. Various types of mathematical models and systems analyses of biological networks have further been used to gain a better understanding of the sources and regulatory mechanisms of such variability. In this work, we present a novel sensitivity analysis method, called molecular density function perturbation (MDFP), for the dynamical analysis of cellular heterogeneity. The proposed analysis is based on introducing perturbations to the density or distribution function of the cellular state variables at specific time points, and quantifying how such perturbations affect the state distribution at later time points. We applied the MDFP analysis to a model of signal transduction pathway involved in TRAIL (tumor necrosis factor-related apoptosis-inducing ligand)-induced apoptosis in HeLa cells. The MDFP analysis showed that caspase-8 activation regulates the timing of the switch-like increase of cPARP (cleaved poly(ADP-ribose) polymerase), an indicator of apoptosis. Meanwhile, the cell-to-cell variability in the commitment to apoptosis depended on mitochondrial outer membrane permeabilization (MOMP) and events following MOMP, including the release of Smac (second mitochondria-derived activator of caspases) and cytochrome-C from mitochondria, the inhibition of XIAP (X-linked inhibitor of apoptosis) by Smac and the formation of apoptosome.

Glioblastoma multiforme (GBM) is a highly heterogeneous brain tumor with limited treatment options and a poor prognosis. Cancer stem cells (CSCs) have emerged as a critical factor in GBM resistance and management, contributing to tumor growth, heterogeneity, and immunosuppression. The transcription factor FOXM1 has been identified as a key player in the progression, spread, and therapy resistance of various cancers, including GBM. In this study, researchers conducted structure-based in silico screening to identify natural compounds that could target the DNA-binding domain (DBD) of the FOXM1 protein. Through molecular docking analyses, identified Silybin B as a potential inhibitor of FOXM1, exhibiting strong interaction with the protein. MD simulations were performed to validate the binding stability of the FOXM1-Silybin B complex. The study provides valuable insights into the potential of Silybin B as a FOXM1 inhibitor and its ability to induce senescence in GBM stem cells. These findings contribute to the development of structure-based design strategies for FOXM1 inhibitors and innovative therapeutic approaches for the treatment of Glioblastoma.

In order to overcome the limitation of the traditional adaptive Unscented Kalman Filtering (UKF) algorithm in noise covariance estimation for statement and measurement, we propose a hybrid adaptive UKF algorithm based on combining Maximum a posteriori (MAP) criterion and Maximum likelihood (ML) criterion, in this paper. First, to prevent the actual noise covariance deviating from the true value which can lead to the state estimation error and arouse the filtering divergence, a real-time covariance matrices estimation algorithm based on hybrid MAP and ML is proposed for obtaining the statement and measurement noises covariance, respectively; and then, a balance equation the two kinds of covariance matrix is structured in this proposed to minimize the statement estimation error. Compared with the UFK based MAP and based ML, the proposed algorithm provides better convergence and stability.

Raman spectra of [TeX6](Y2) (X, Y - Cl, Br, I) were recorded in the temperature range 5-300 K. The vibration frequency of molecular halogens as function on temperature and the type of halogen X- in the [TeX6]2- complex was obtained. The vibration frequency of the molecular halogen Y2 depends on the type of halogen X- complex [TeX6]2- The temperature dependence of the Y2 frequency is anomalous compared to what was expected for a conventional anharmonic oscillator. It is shown that the interpretation of the experimental results based on the concept of traditional chemical bonding is quite satisfactory in all cases. The study was motivated by the need to consider the mechanism of interaction in the X-Y2 contacts.

During the last years there has been a high interest in the development of new purely-organic single-component conductors. Very recently, we have reported a new neutral radical conductor based on the perchlorotriphenylmethyl (PTM) radical moiety linked to a monopyrrolo-tetrathiafulvalene (MPTTF) unit by a -conjugated bridge (1). Interestingly, this system behaves as a semiconductor with high conductivity and small energy gap under high pressure. With the aim of developing a new material with improved conducting properties, we have designed and synthesized the radical dyad 2 which was functionalized with an ethylenedithio (EDT) group in order to improve the intermolecular interactions of the TTF subunits. The physical properties of the new radical dyad 2 were studied in detail in solution to further analyze its electronic structure as well as its potential use as radical conductor in the solid state.

The details of antigen-antibody interactions and the identification of epitopes are critical for the development of monoclonal antibody drugs. Ab42 is a native human-derived anti-CFH monoclonal antibody. In this study, the interaction between antigen pCFH and antibody (Ab42) was theoretically demonstrated by molecular docking and MD simulation, combined with free energy calculation and computational alanine scanning (CAS), and key amino acids and epitopes were identified. Experimental alanine scanning (EAS) was then carried out to verify the results of the calculation, and our results indicated that Ab42 antibody forms hydrogen bonds and interacts hydrophobically with pCFH through the Tyr315, Ser100, Gly33, and Tyr53 residues on its CDR, while the main pCFH epitopes are located at the six sites of Pro441, Ile442, Asp443, Asn444, Ile447, and Thr448. In conclusion, this study has explored the mechanism of antigen-antibody interaction from both theoretical and experimental aspects, and our results have important theoretical significance for the design and development of relevant antibody drugs.

The effect of graphene (GR) on Ni surface relaxation and reconstruction in three different substrate orientations: {111}, {001}, {011}, at two different temperatures: 300K and 400K was studied using Molecular Dynamics simulation. The change of the interplanar distances of the substrate and redistribution of Ni and C atoms in direction perpendicular to the surface were compared with the equilibrium state of GR and bulk Ni, in the absence of the counterpart. The surface reconstruction for GR/Ni system was analyzed base on the calculated radial pair distribution function of Ni and C atoms. The surface roughness was visualized using 2D atomic distribution maps. For the studied substrate orientations and temperatures, it was found that the most densely packed orientation of the Ni base {111} provides minimal changes of the structural parameters of both counterparts at 400 K.

Zirconium nitride (Zr-N) is a candidate material for inert matric fuel elements in breeder reactors. They have shown exceptional resistance to radiative corrosive environments. They have shown excellent resistance to radiation stability at operational conditions in current generation reactors. But being a candidate nuclear fuel material in generation reactors, their radiation stability at higher temperature conditions observed has not been reported. We have employed molecular dynamics simulations to study the damage cascades and radiation stability of ZrN at 600K. The process of primary damage evolution by considering varying overlapping cascades with 10, 20 keV incident PKA’s reported here. The defects interface interaction and diffusion mechanism at the interface are reported.

In this paper, we introduce a simple yet powerful and working version of the molecular dynamics code using the Python 3.9 language. The code contents are published in the link given in the appendix 1. The structure and components of the program is given in detail using flowcharts and code snippets. The program consists of major features like velocity verlet integrator, thermostats, COM removal, input and output modules, virial, pressure, and other thermodynamic quantities estimation etc. The author believes that this program will be helpful to graduate students who perform research in molecular dynamics simulations who intend to write their own code instead of the sophisticated open source packages.

Hybrid lead halide perovskites have been revolutionary in the photovoltaic research field, reaching efficiencies comparable with the most established photovoltaic technologies, although they do not yet reach their competitor stability. The search for a stable configuration required the engineering of the charge extraction layers; in this work, molecular doping is used as an efficient method for small molecule and polymer, employed as hole transport materials in planar heterojunction configuration on compact-TiO2. We proved the viability of this approach, obtaining significantly increased performances and reduced hysteresis on compact titania-based devices. We investigated the photovoltaic performance, correlating to the hole transport material structure. We have demonstrated that the molecular doping mechanism is more reliable than the oxidative doping, and verified that molecular doping in polymeric hole transport materials leads to highly efficient perovskite solar cell, with long-term stability.

Cyclooxygenases (COX), or prostaglandin endoperoxide synthases (PTGS), are key enzimes in the synthesis of prostaglandins, which are chemical species critical in mediating inflammatory processes. There are two highly homologous COX isoforms: COX-1 and COX-2. COX-1 is involved in the production of prostaglandins, chemical compuounds that take part in physiological processes such as: protection of the gastric epithelium, maintenance of renal flow, platelet aggregation, neutrophil migration and, also, are expressed in the vascular endothelium. Meanwhile, COX-2 is induced by proinflammatory stimuli. It is very frequent the use of nonsteroidal antiinflammatory drugs (NSAIDs) to counteract the symptoms of inflammatory processes. These drugs, in addition to its benefits, can cause side effects on people’s health, such as cardiovascular and respiratory problems, among others. In the past years, it has been recognized the potential of plants secondary metabolites as pharmacological agents, prompting the need for investigations that shed light into its mechanism of action. In this work we have applied computational techniques, based on quantum chemistry and mechanical statistics, to study the protein-ligand interaction involving COX’s and secondary metabolites from natural sources. Our aim is to determine the structure activity interplay in processes involving the participation of secondary plant metabolites such as luteolin, galangin, kaempferol, apigenin, morine and quercetin on the inactivation of COX’s. From molecular docking analysis, we have extracted the energetics of the COX-(1,2)/metabolite coupling. By defining energy based factors, we have determined a procedure that predicts the chemical species with highest stability and selectivity towards inactivation of COX-2 over COX-1. The results are discussed with regard to conformational features of the selected ligands and its intermolecular strong/weak interactions inside the active-sites of the COX’s hosts.

The novel coronavirus disease 19 (COVID-19) has resulted in an estimated 20 million excess deaths and the recent resurgence of COVID-19 in China is predicted to result in up to 1 million deaths over the next few months. With vaccines unable to halt transmission it is important to continue our quest for safe, effective, affordable drugs that will be available to all countries. Drug repurposing is one of the strategies being explored in this context. Recently, out of 7,817 approved drugs, 214 candidates were systematically down-selected using a combination of 11 filters including approval status, assay data against SARS-CoV-2, pharmacokinetic, pharmacodynamic and toxicity profiles. These drugs were subjected in this study to virtual screening against various targets of SARS-CoV-2 followed by molecular dynamic studies of the best scoring ligands against each target. The chosen molecular targets were Spike receptor binding domain, Nucleocapsid protein RNA binding domain, and key non-structural proteins 3, 5, 12, 13 and 14. Four drugs approved for other indications — alendronate, cromolyn, natamycin and treprostinil — look sufficiently promising from our in silicostudies to warrant further in vitro and in vivo investigations as appropriate to ascertain their extent of anti-viral activities.

Li₂MnO₃ is critical component in the well-studied Li-excess cathode materials xLi₂MnO₃•(1- x)LiMO₂ for achieving high lithium storage capacity. In this article, the diffusion of Li ions in Li₂MnO₃ is studied using molecular dynamics (MD) simulation with well-behaved empirical force fields obtained by fitting against the crystal structure from experiment and phonons calculated using Density Function Theory (DFT). We have found two possible tetrahedral hopping channels, 0-TM and 1-TM(Mn⁴⁺) channel, which are differentiated by the face sharing octahedral cations. Simulation results show that the 0-TM channel is active for Li hopping, while 1-TM(Mn4+) channel is inactive. During the delithiation process, the Li ions in the TM layered are firstly removed, then those in the Li layer. However, the Li ions will be trapped in the tetrahedral 0-TM channels as long as the four face sharing octahedral sites are cleared. Up to x=1.0 for Li₂-xMnO₃, almost all the Li ions are located at the tetrahedral sites, forming a regular array along a axis. The de-intercalation of tetrahedral Li ions requires a high voltage (>5.2 V vs. Li/Li+), limiting the practical capacities measured in lab. The diffusion of Mn ions into the Li layers is observed in a deeper delithiated structure (x=1.2 for Li₂-xMnO₃), indicating an initial phase transformation to a spinel-like structure. However, the Mn ions are mainly trapped in the tetrahedral sites in the Li layer, instead of the octahedral sites fin spinel-like structure. A few of Mn ions diffusing into the octahedral sites in Li layers have no face sharing tetrahedral Li ions, revealing a further Li de-intercalation is imperative for the complete phase transformation. Our model is not stable for x≥1.4 in Li₂-xMnO₃. Other charge compensation mechanism should be considered in this high delithiation stage, eg. oxygen release.

Stimulation of the alpha 7 nicotinic acetylcholine receptor (α7nAChR) has shown beneficial effects in several acute inflammatory disease models. This study aims to examine whether treatment with the selective α7nAChR agonist PHA 568487 can dampen inflammation and thereby improve cardiac function after myocardial infarction in mice. The possible anti-inflammatory properties of α7nAChR agonist PHA 568487 were tested in vivo using the air pouch model and in a permanent occlusion model of acute myocardial infarction in mice. Hematology parameters and cytokine levels were determined. Infarct size and cardiac function were assessed via echocardiography 24 h and one week after the infarction. Treatment with α7nAChR agonist PHA 568487 decreased 12 (CCL27, CXCL5, IL6, CXCL10, CXCL11, CXCL1, CCL2, MIP1a, MIP2, CXCL16, CXCL12 and CCL25) out of 33 cytokines in the air pouch model of acute inflammation. However, α7nAChR agonist PHA 568487 did not alter infarct size, ejection fraction, cardiac output or stroke volume at 24 h nor at 7 days after the myocardial infarction compared with control mice. In conclusion, despite promising immunomodulatory effects in the acute inflammatory air pouch model, α7nAChR agonist PHA 568487 did not affect infarct size or cardiac function after a permanent occlusion model of acute myocardial infarction in mice. Consequently, this study does not strengthen the hypothesis that stimulation of the α7nAChR is a future treatment strategy for acute myocardial infarction when reperfusion is lacking. However, whether other agonists of the α7nAChR can have different effects remains to be investigated.

Histone Lysine Specific Demethylase 1 (LSD1) is overexpressed in many cancers and become a new target for anticancer drugs. In recent years, the small molecule inhibitors with various structures targeting LSD1 have been reported. Here we report the binding interaction modes of a series of thieno[3,2-b]pyrrole-5-carboxamides LSD1 inhibitors using molecular docking, three dimensional quantitative structure-activity relationship (3D-QSAR). Comparative molecular field analysis (CoMFA q²=0.783, r²=0.944, r²_pred=0.851) and Comparative molecular similarity indices analysis (CoMSIA q²=0.728, r²=0.982, r²_pred=0.814) were used to establish 3D-QSAR models, which had good verification and prediction capabilities. Based on the contour maps and the information of molecular docking, 8 novel small molecules were designed in silico, among which compounds D4, D5 and D8 with high predictive activity were subjected to further molecular dynamics simulations (MD), and their possible binding modes were explored. It was found that Asn535 plays a crucial role in stabilizing the inhibitors. Furthermore, the ADME and bioavailability prediction for D4, D5 and D8 were carried out. The results would provide valuable guidance for designing new reversible LSD1 inhibitors in the future.

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.

Plant transcriptomes are an extremely dynamic entities shaped spatially and temporally by many intracellular and environmental cues. In this review, we first summarize the complexity and diversity of plant genomes and transcriptomes as a start point for the multitude of transcriptomic responses. Numerous alterations within various tissue and organ-specific transcriptomes as well as the most relevant transcriptomic responses associated with plant acclimation to selected abiotic and biotic stress conditions, from the current studies employing high-throughput transcriptomic analysis are widely discussed. Understanding changes within plant transcriptomes, revealed by in silico functional analysis, allows for the characterization of stress-affected genes and stress acclimatory mechanisms, as well as allows to perform plant metabolic engineering. The latter allow cultivars to produce more secondary metabolites in the future, which are often desirable substances in the biomedical industry. Accordingly, in this review special attention was also paid to characterize the potential of transcriptomic analyses of medicinal species, particularly to search for new cultivars. Extensive characterization of transcriptomic responses in stress would also result in the development of new cultivars that display physiological and molecular mechanisms that allow them to cope with adverse environmental conditions more adequately.

A major public health emergency has been created by the COVID-19 pandemic, which is brought on by the SARS-CoV-2 virus. Due to its crucial function in viral replication, the primary protease (Mpro) of the virus is a prime target for therapeutic research. In this study, we used molecular modeling and molecular dynamics simulations to examine the potential therapeutic uses of Mpro inhibitors for the treatment of COVID-19. Using induced fit docking and molecular dynamics simulations, we confirmed the top compounds after screening a library of compounds for their ability to bind to Mpro. Simulation interaction diagrams were used to investigate protein-ligand interactions, and MM-GBSA was used to determine binding energies. The Swiss ADME server was used to predict ADME properties. According to our findings, numerous substances are strong COVID-19 medication candidates since they have excellent ADME features and high binding affinities. This work serves as a foundation for additional experimental research and drug development initiatives aimed at Mpro.

Ab initio density functional theory molecular dynamics extensive calculations were used for the first time to evaluate stability conditions for relevant phases of InN. In particular, thermal decomposition of InN with formation of N2 molecules and pressure induced wurtzite – rocksalt solid-solid phase transition were established. The results consistent with available experimental data allowed a critical evaluation of potential and limitations of the proposed simulation method. The ab initio molecular dynamics as an efficient tool for simulations of phase transformations of InN including solid-solid structural transition and thermal decomposition with formation of N2 molecules has been demonstrated. It is of high interest because InN being important component of epitaxial quantum structures has still not been obtained as a bulk single crystal making difficult a determination of its basic physical properties and designing new applications.

This study comprehensively investigates Al2O3's mechanical properties, focusing on fracture toughness, surface energy, Young's modulus, and crack propagation. The density functional theory (DFT) is employed to model the vacancies in Al2O3, providing essential insights into this material’s structural stability and defect formation. The DFT simulations reveal a deep understanding of vacancy-related properties and their impact on mechanical behavior. In conjunction with molecular dynamics (MD) simulations, the fracture toughness and crack propagation in Al2O3 are explored, offering valuable information on material strength and durability. The surface energy of Al2O3 is also assessed using DFT, shedding light on its interactions with the surrounding environment. The combination of DFT and MD simulations provides a robust framework for a comprehensive understanding of Al2O3's mechanical properties, with implications for material science and engineering applications.

Nuclear magnetic resonance (NMR) techniques are widely used to identify pure substances and probe protein dynamics. Edible oil is a complex mixture composed of hydrocarbons, which have a wide range of molecular size distribution. In this research, low field NMR (LF-NMR) relaxation characteristic data of various sample oils were analyzed. We also suggest a new method for prediction size of edible oil molecules using LF-NMR relaxation time distributions. According to the relative molecular mass, carbon chain length, transverse relaxation time of different sample oils, combined with oil viscosity and other factors, the relationship between carbon chain length and transverse relaxation time rate was analyzed. Various oils and fats in the mixed fluid are displayed, reflecting the composition information of different oils. To study the correlation between the rotation correlation time and molecular information of oil molecules. The molecular composition of the resulting fluid determines its properties, such as viscosity and phase behavior. The results show that low-field NMR can obtain information on the composition, macromolecular aggregation, and molecular dynamics of complex fluids. Measurements of grease in the free fluid state show that the relaxation time can reflect the intrinsic properties of the fluid. It is shown that the composition characteristics and states of complex fluids can be measured by low-field nuclear magnetic resonance.

Our work uses Iterative Boltzmann Inversion (IBI) to study the coarse-grained interaction between 20 amino acids and the representative carbon nanotube CNT55L3. IBI is a multi-scale simulation method that has attracted the attention of many researchers in recent years. It can effectively modify the coarse-grained model derived from the Potential of Mean Force (PMF). IBI is based on the distribution result obtained by All-Atom molecular dynamics simulation, that is, the target distribution function, the PMF potential energy is extracted, and then the initial potential energy extracted by the PMF is used to perform simulation iterations using IBI. Our research results have gone through more than 100 iterations, and finally, the distribution obtained by coarse-grained molecular simulation (CGMD) can effectively overlap with the results of all-atom molecular dynamics simulation (AAMD). In addition, our work lays the foundation for the study of force fields for the simulation of the coarse-graining of super-large proteins and other important nanoparticles.

This study delves into the dynamic field of medical education research, emphasizing the integration of research training within medical curricula. It seeks to understand the impact of such integration on the competencies of future medical professionals. The primary aim is to systematically categorize and analyze the current trends and future directions in medical education research. This involves assessing the influence of research training on medical students' skills and the methodologies used in such research. Methods: The research employs an extensive bibliographic literature review across multiple databases. It classifies studies into categories like experiential or case studies, editorials, and original research articles. This classification is based on criteria such as geographical location, research objectives, theoretical frameworks, and methodologies. Results: Findings reveal a diverse landscape in medical education research, with a significant emphasis on research training. The research showcases varying methodologies and approaches used globally, highlighting the thematic focus and geographical distribution of these studies. Conclusion: Medical education research is a globally expansive and evolving field. It underscores the importance of continuous investigation, particularly focusing on integrating research elements at curricular levels and exploring innovative educational strategies. The study also points out potential gaps in research, especially in underrepresented regions, indicating directions for future research efforts.

Species of the genus Drosophila have served as favorite models in speciation studies; however, genetic factors of interspecific reproductive incompatibility are underinvestigated. Here we performed the analysis of hybrid female sterility by crossing Drosophila melanogaster females and Drosophila simulans males. Using transcriptomic data analysis, molecular, cellular, and genetic approaches, we analyzed differential gene expression, transposable element (TE) activity, piRNA biogenesis, and functional defects of oogenesis in hybrids. Premature germline stem cell loss was the most prominent defect of oogenesis in hybrid ovaries. Due to the differential expression of genes encoding piRNA pathway components, rhino and deadlock, the functional RDCmel complex in hybrid ovaries was not assembled. However, the activity of the RDCsim complex was maintained in hybrids, independent of the genomic origin of piRNA clusters. Despite the identification of a cohort of overexpressed TEs in hybrid ovaries, we found no evidence that their activity can be considered the main cause of hybrid sterility. We revealed a complicated pattern of Vasa protein expression in the hybrid germline, including partial AT-chX piRNA targeting of the vasasim allele and a significant zygotic delay in vasamel expression. We arrived at the conclusion that the hybrid sterility phenotype was caused by intricate multi-locus differences between the species.

New class of triazinoindole bearing thiosemicarbazide (1-25) was synthesized and evaluated for α-glucosidase inhibitory potential. All synthesized analogues exhibited excellent inhibitory potential having IC50 values ranging from 1.30 ± 0.01 to 35.80 ± 0.80 µM when compared with the standard acarbose having IC50 value 38.60 ± 0.20 µM. Among series the analogues 1 and 23 was found the most potent having IC50 values 1.30 ± 0.05 and 1.30 ± 0.01 µM respectively. Structure activity relationship (SAR) was mainly based upon by bring about difference of substituents on phenyl rings. To confirm the binding interactions, molecular docking study was performed. Synthesized analogues were characterized through HREI-MS, 1H and 13C-NMR analysis.

Stueckelberg-Horwitz-Piron (SHP) electrodynamics formalizes the distinction between coordinate time (measured by laboratory clocks) and chronology (temporal ordering) by defining 4D spacetime events x^μ as functions of an external evolution parameter τ. Classical spacetime events x^μ (τ) evolve as τ grows monotonically, tracing out particle worldlines dynamically and inducing the five U(1) gauge potentials through which events interact. Since Lorentz invariance imposes time reversal symmetry on x⁰ but not τ, the formalism resolves grandfather paradoxes and related problems of irreversibility. The action involves standard first order field derivatives but includes a higher order τ derivative that while preserving gauge and Lorentz invariance removes certain singularities and makes the related QFT super-renormalizable. The resulting field equations are Maxwell-like but τ-dependent and sourced by a current that represents a statistical ensemble of instantaneous events distributed along the worldline. The width λ of this distribution defines a correlation time for the interactions and a mass spectrum for the photons that carry the interaction. As λ becomes very large, the photon mass goes to zero and the field equations become τ-independent Maxwell’s equations. Maxwell theory thus emerges as an equilibrium limit of SHP, in which λ is larger than any other relevant time scale. Particles and fields are not constrained to mass shells in SHP theory, and by exchanging mass may produce pair creation/annihilation processes at the classical level. On-shell evolution with fixed particle masses is restored through a self-interaction associated with the 5D wave equation.

In patients with heart failure, the evaluation of left ventricular (LV) diastolic function is vital, offering crucial insights into hemodynamic impact and prognostic accuracy. Echocardiography stands as the primary imaging modality for diastolic function assessment, and using it effectively requires a profound understanding of the underlying pathology. This review covers four main topics: first, the fundamental driving forces behind each phase of normal diastolic dynamics, along with the physiological basis of two widely used echocardiographic assessment parameters, E/e' and mitral annulus early diastolic velocity (e'); second, the intricate functional relationship between the left atrium and LV in patients with varying degrees of LV diastolic dysfunction (LVDD); third, the role of stress echocardiography in diagnosing LVDD and the significance of parameter changes in this context; and fourth, the clinical utility of evaluating diastolic function from echocardiography images across diverse cardiovascular care areas.

Attribute reduction is a core technique in the rough set domain and an important step in data preprocessing. Researchers have proposed numerous innovative methods to enhance the capability of attribute reduction, such as the emergence of multi-granularity rough set models, which can effectively process distributed and multi-granularity data. However, these innovative methods still have many shortcomings, such as how to deal with complex constraints and how to perform multi-angle effectiveness evaluations. Based on the multi-granularity model, this study proposes a new method of attribute reduction, namely using multi-granularity neighborhood information gain ratio as the measurement criterion. This method combines both supervised and unsupervised perspectives, and by integrating multi-granularity technology with neighborhood rough set theory, constructs a model that can adapt to multi-level data features. This model can select the optimal granularity level or attribute set according to the requirements. Finally, the proposed method is experimented with on 15 UCI datasets against six other attribute reduction algorithms. The experimental outcomes, upon analysis and comparative review, substantiate the reliability and consistency of the attribute reduction approach introduced in this study.

A facilitative transport system exists on the blood brain barrier (BBB) that has been tacitly assumed to be a path for glutamate entry to brain. But glutamate is a non-essential amino acid whose brain content is much greater than plasma, and studies in vivo show that glutamate does not enter brain in material quantities except in those small regions with fenestrated capillaries (circumventricular organs). The situation became understandable when luminal (blood facing) and abluminal (brain facing) membranes were isolated and studied separately. Facilitative transport of glutamate and glutamine exist only on the luminal membranes whereas Na+-dependent transport systems for glutamate, glutamine and some other amino acids are present only on the abluminal membrane. The Na+-dependent cotransporters of the abluminal membrane are in a position to actively transport amino acids from the extracellular fluid (ECF) into the endothelial cells of the BBB. These powerful secondary active transporters couple the energy of the Na+-gradient to move glutamate and glutamine into the ECF whereupon glutamate can exit to blood on the luminal facilitative glutamate transporter. Glutamine may also exit brain on a separate facilitative transport system that exists on the luminal membranes or glutamine can be hydrolyzed to glutamate within the BBB thereby releasing ammonia that is freely diffusible. The γ-glutamyl participate cycle participates indirectly by producing oxoproline (pyroglutamate) that stimulates almost all secondary active transporters yet discovered in the abluminal membranes of the BBB.

The River Habitat Survey (RHS) is a method for assessing the conditions of a river channel and valley, taking into account hydromorphological factors, and assessing the degree of naturalness or the degree of anthropogenic alteration within the channel and valley based on physical and landscape data. The hypothesis tested in this study was that the quality of river habitats demonstrated using the RHS should translate into the composition of their fauna. We investigated aquatic bugs (Heteroptera Aquatica), which are a group with considerable dispersal abilities and, at the same time, composed of stenotopic species associated with a particular environment as well as of eurytopic elements inhabiting a wide range of environments (REF). The research was carried out in a small lowland river, the River Krąpiel (north-western Poland) at 6 sites with 43 sub-sites included. Our analyses suggest the RHS method can predict the composition and quality of the aquatic bug fauna inhabiting a given river. Indicators of river habitat quality (HQA and RHQ) will be associated with different habitat types depending on the geographic location and type of river, once indicating habitats associated with the mainstream river and mineral substrate (mountain and small lowland rivers) and at other times habitats covered by vegetation with relatively slow water flow (large and medium lowland rivers). Therefore, RHS may be a good method for estimating the water bug fauna, linking it to the size and character of the river.

The use of carbon nanotubes to improve the mechanical properties of polymers is one of the promising directions in materials science. The addition of single-walled carbon nanotubes (SWCNT) to a polymer results in significant improvement of mechanical, electrical, optical, and structural properties. However, the addition of SWCNTs does not always improve the polymer properties. Also, when a certain content of SWCNTs is exceeded, the mechanical properties of the nanocomposite are getting worse. This article reports the results of computer simulations for predicting the mechanical properties of polymer/single-walled carbon nanotube nanocomposites. The efficiency of reinforcing polymer composites is considered depending on the concentration of carbon nanotubes in the polymer matrix, their size and structure. The elastic moduli of the nanocomposites were predicted using computer simulations for unit cell tension (0.1%). General trends in the mechanical properties of composites with polypropylene, poly (ethyl methacrylate), polystyrene matrices and SWCNTs are shown.

The construction of molecular robotic-like objects that imitate living things is an important challenge for current chemists. Such molecular devices are expected to perform their duties robustly to carry out mechanical motion, process information, and make independent decisions. Dissipative self-organization plays an essential role in meeting these purposes. To produce a micro-robot that can perform the above tasks autonomously as a single entity, a function generator is required. Although many elegant review articles featuring chemical devices that mimic biological mechanical functions have been published recently, the dissipative structure, which is the minimum requirement, has not been sufficiently discussed. This article aims to show clearly that dissipative self-organization is a phenomenon involving autonomy, robustness, mechanical functions, and energy transformation. Moreover, the author details the recent experimental results of an autonomous light-driven molecular device that achieves all of these features. In addition, a chemical model of cell-amplification is also discussed to focus on the generation of hierarchical movement by dissipative self-organization. By reviewing this research, it may be perceived that mainstream approaches to synthetic chemistry have not always been appropriate. In summary, the author proposes that the integration of catalytic functions is a key issue for the creation of autonomous microarchitecture.

2-Hydroxypropyl-β-cyclodextrin (HPβCD) has unique properties to enhance the stability and the solubility of low water-soluble compounds by inclusion complexation. Understanding of the structural properties of HPβCD and its derivatives based on the number of 2-hydroxypropyl (HP) substituents at the a-D-glucopyranose subunits is rather important. In this work, replica exchange molecular dynamics simulations were performed to investigate the conformational changes of single- and double-sided HP-substitution called as 6-HPβCDs and 2,6-HPβCDs, respectively. The results show that glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have lower chance to flip than in βCD. Also, HP groups are occasionally blocking the hydrophobic cavity of HPβCDs, thus hindering the drug inclusion. We found that HPβCDs with high number of HP-substitutions are more likely to be blocked, while HPβCDs with double-sided HP-substitution are even more probable to be blocked. Overall, 6-HPβCDs with three and four HP-substitutions are highlighted as the most suitable structures for guest encapsulation based on our conformational analyses such as structural distortion, radius of gyration, circularity and cavity self-closure of the HPβCDs.

Micro-pollutants such as 17β-Estradiol (E2) have been detected in different water resources and their negative effects on the environment and organisms have been observed. Aptamers are established as a possible detection tool, but the underlying ligand binding is largely unexplored. In this study, a previously described 35-mer E2-specific aptamer was used to analyse the binding characteristics between E2 and the aptamer with a MD simulation in an aqueous medium. Because there is no 3D structure information available for this aptamer, it was modeled using coarse-grained modeling method. The E2 ligand was positioned inside a potential binding area of the predicted aptamer structure, the complex was used for an 25 ns MD simulation, and the interactions were examined for each time step. We identified E2-specific bases within the interior loop of the aptamer and also demonstrated the influence of frequently underestimated water-mediated hydrogen bonds. The study contributes to the understanding of the behavior of ligands binding with aptamer structure in an aqueous solution. The developed workflow allows generating and examining further appealing ligand-aptamer complexes.

In this paper we introduce a computational model for the simulation of hypervelocity impact (HVI) phenomena which is based on the Discrete Element Method (DEM). Our paper constitutes the first application of DEM to the modeling and simulating of impact events for velocities beyond 5 kms⁻¹. We present here the results of a systematic numerical study on HVI of solids. For modeling the solids, we use discrete spherical particles that interact with each other via potentials. In our numerical investigations we are particularly interested in the dynamics of material fragmentation upon impact. We model a typical HVI experiment configuration where a sphere strikes a thin plate and investigate the properties of the resulting debris cloud. We provide a quantitative computational analysis of the resulting debris cloud caused by impact and a comprehensive parameter study by varying key parameters of our model. We compare our findings from the simulations with recent HVI experiments performed at our institute. Our findings are that the DEM method leads to very stable, energy–conserving simulations of HVI scenarios that map the experimental setup where a sphere strikes a thin plate at hypervelocity speed. Our chosen interaction model works particularly well in the velocity range where the local stresses caused by impact shock waves markedly exceed the ultimate material strength.

Currently, obesity is a critical global public health burden. Numerous studies have demonstrated the regulation of the pathogenesis of obesity and metabolic abnormalities by the gut microbiota and microbial factors; however, their involvement in the various degrees of obesity is not yet understood well. Previously, obesity has been shown to be associated with decreased levels of vitamin B12. Considering exclusive microbial production of vitamin B12, we hypothesised that a decrease in cobalamin levels in obese individuals may be at least partially caused by its depleted production in the intestinal tract by the commensal microbiota. In the present study, our aim was to estimate the abundance of enzymes and metabolic pathways for vitamin B12 synthesis in the gut microbiota of mouse models of alimentary and genetically determined obesity, to evaluate the contribution of the obesogenic microbiome to vitamin B12 synthesis in the gut. We have defined a significantly lower representation of enzymes and metabolic pathways for vitamin B12 biosynthesis in obese mice compared to non-obese mice, wherein enzyme depletion was more pronounced in lepr(-/-) (db/db) mice, which developed severe obesity. The representation of enzymes involved in cobalamin synthesis is strongly correlated with representation of several microbes in high-fat diet fed mice, whereas there were almost no correlations in db/db mice. Therefore, the degree of obesity and the composition of the correspondent microbiota are the main contributors to the representation of genes and pathways for cobalamin biosynthesis in the mouse gut.

The oncogenic potential of high-risk HPVs is focused on producing the E6 and E7 oncoproteins responsible for disrupting the control of the cell cycle. Epidemiological studies propose the presence of the N29S and H51N variants of the HPV16 E7 protein as a significant association with cervical cancer. It has been suggested that changes in the amino acid sequence of E7 variants may affect the oncoprotein 3D structure; however, this remains unknown. Analysis of the structural differences of the HPV16 E7 protein and its variants (N29S and H51N) was performed through homology modeling and structural refinement by molecular dynamics simulation. We propose for the first time a 3D structure of the E7 reference protein and two of its variants (N29S and H51N) and conclude that the mutations induced by the variants in N29S and H51N have a significant influence on the 3D structure of the E7 protein of HPV16, which could be related to the oncogenic capacity of this protein.

Meta-aminobenzoic acid, an important model system in the study of polymorphism and crystallization of active pharmaceutical ingredients, exist in water in both the nonionic (mABA) and zwitterionic (mABA^±) forms. However, the constituent molecules of the polymorph that crystallizes from aqueous solutions are zwitterionic. This study reports atomistic simulations of the events surrounding the early stage of crystal nucleation of meta-aminobenzoic acid from aqueous solutions. Ab initio molecular dynamics was used to simulate the hydration of mABA± and mABA, and to quantify the interaction of these molecules with the surrounding water molecules. Density functional theory calculations were conducted to determine the low-lying energy conformers of meta-aminobenzoic acid dimers and compute the Gibbs free energies in water of nonionic, (mABA)₂, zwitterionic, (mABA^±)₂ and nonionic-zwitterionic, (mABA)(mABA^±), species. Classical molecular dynamics simulations of mixed mABA–mABA^± aqueous solutions were carried out to examine the aggregation of meta-aminobenzoic acid. According to these simulations the selective crystallization of the polymorph which constituent molecules are zwitterionic is driven by the formation of zwitterionic dimers in solution, which are thermodynamically more stable than (mABA)₂ and (mABA)(mABA^±) pairs. This work represents a paradigm of the role of molecular processes during the early stages of crystal nucleation in affecting polymorph selection during crystallization from solution.

Reactive oxygen species (ROS) play a pivotal role in cellular function, and their overproduction contributes to oxidative stress, fatigue, aging, and diseases. To counteract these effects, incorpo-rating antioxidant-rich nutrients is recommended. Functional waters, including electrolyzed hydrogen water (EHW), are gaining attention for potential antioxidant effects through daily in-take. This study employs aquaphotomics and near-infrared spectroscopy (NIRS) to elucidate the molecular mechanisms behind the antioxidant activity of EHW. Previous investigation reported a 100% loss in antioxidant functionality through autoclaving, and a 40% decrease with hydrogen gas removal, prompting the current investigation into the mechanisms of antioxidant effects and especially the role of water molecular structure in it. Aquaphotomics analysis of near infrared (NIR) spectral data revealed that the primary characteristic of EHW associated with anti-radical function is presence of water solvation shells. Autoclaving and degassing led to loss of these water molecular species and a significant increase in free water molecules. The study suggests that water shells (solvation shells) in EHW play a crucial role in neutralizing ROS, although the specific ions or nanoparticles responsible for creating these shells remain unidentified. This concept aligns well with existing proposals of the crucial role of water solvation shells in ROS neutralization.

Efficient thermal management of modern electronics requires the use of thin films with highly anisotropic thermal conductivity. Such films enable the effective dissipation of excess heat along one direction while simultaneously providing thermal insulation along the perpendicular direction. In this study, we investigate the thermal conductivity of bilayer graphene (BLG) sheets using non-equilibrium molecular dynamics, examining both in-plane and cross-plane thermal conductivities. The in-plane thermal conductivity of 10 nm × 10 nm BLG with zigzag and armchair edges at room temperature is found to be around 204 W/m.K and 124 W/m.K, respectively. The in-plane thermal conductivity of BLG increases with sheet length. For all sizes, BLG with a zigzag edge exhibits higher thermal conductivity than that of the armchair edge, with this difference becoming more pronounced as the lengths of the sheets increase. In addition, increasing temperature from 300 K to 600 K decreases the in-plane thermal conductivity of a 10 nm× 10 nm zigzag BLG by about 34%. Similarly, the application of a 12.5% tensile strain induces a 51% reduction in its thermal conductivity compared to the strain-free values. Armchair configurations exhibit similar responses to variations in temperature and strain, but with less sensitivity. Furthermore, the cross-plane thermal conductivity of BLG at 300 K is estimated to be 0.05 W/m·K, significantly lower than the in-plane results by approximately four orders of magnitude. The cross-plane thermal conductance of BLG decreases with increasing temperatures, specifically, at 600 K, its value is almost 16% of that observed at 300 K.

The past two decades have witnessed the introduction of and then a steady increase in the use of computational techniques in the study and development of molecularly imprinted polymers (MIPs). Molecular dynamics (MD) based studies have had a significant role in this development as they can provide insights concerning the mechanisms governing the molecular level events underlying MIP synthesis and MIP-ligand interactions and can be used for the identification of preferred monomer compositions and for the prediction of MIP properties. We here review the role that MD has played in the development of molecular imprinting and examine the different types of MD strategies that have been used, including their advantages and challenges. Recent trends in the application of MD to the study of MIPs are presented, along with a perspective on the future importance of MD-based studies for the development of molecular imprinting science and technology.

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 alignment of nanostructures in materials such as lyotropic liquid crystals (LLCs) templated materials has the potential to signicantly improve their performances. However, accurately characterising and quantifying the alignement of such fine structures remains very challenging. In situ small angle X-ray scattering (SAXS) and molecular dynamics were employed for the first time to understand the hexagonal LLC alignment process with magnetic nanoparticles under a magnetic field. The enhanced alignment has been illustrated from the distribution of azimuthal intensity in the samples exposed to magnetic field. Molecular dynamics simulations reveal the relationship between the imposed force of the magenetic nanoparticles under magnetic field and the force transferred to the LLC cylinders which leads to the LLC alignment. The combinational study with experimental measurement and computational simulation will enable the development and control of nanostructures in novel materials for various applications.

COX II containing a dual core CuA active site is one of the three core subunits of mitochondrial Cco, which plays a significant role in the physiological process. In this report, the full-length cDNA of COXⅡ gene was cloned from Sitophilus zeamais, which had an ORF of 684 bp encoding 227 amino acids residues. The predicted COXⅡ protein had a molecular mass of 26.2 kDa with pI value of 6.37, and multiple sequence alignment and phylogenetic analysis indicated that Sitophilus zeamais COXⅡ had high sequence identity 78.51% with the COXⅡ of other insect species, especially similarity to sitophilus oryzae. This gene was subcloned into the prokaryotic expression vector pET-32a, and induced by IPTG in E.coli Transetta (DE3) expression system. Finally the COXⅡ with 6-His tag was purified using affinity chromatography with Ni²⁺-NTA agarose. WB showed the recombinant COXⅡ was about 44 kD, and the concentration of fusion protein was 50μg/mL. UV-spectrophotometer and infrared spectrometer analysis showed that recombinant COXⅡ could catalyze the oxidation of substrate Cytc, and influenced by AITC. It was found that AITC could form a hydrophobic region with COXⅡ protein via molecular docking, besides, a sulfur atom of AITC structure could form a length of 2.9 Å hydrogen bond with Leu-31. These results will provide valuable information for elucidating the role of COXⅡ in Sitophilus zeamais responses to AITC, meanwhile, it will helpful to carry out a point mutation in AITC binding sites for the future research.

Members of the genus Armillaria are widespread forest pathogens against which effective protection has not yet been developed. Due to their longevity and the creation of large-scale cloning of Armillaria individuals, the use of mycoviruses as biocontrol agents (BCAs) against these pathogens could be an effective alternative. This work describes the screening of double-stranded RNA banding patterns as a first approach for potential virus discovery and the use of stranded total RNA-seq for viral detection and characterization in Armillaria spp. collected in the Czech Republic. A total of five single-stranded RNA viruses were detected in Armillaria ostoyae and A. cepistipes, including the first virus of the family Tymoviridae found infecting a member the genus Armillaria and four viruses belonging to the recently described virus group "ambiviruses" with circular ambisense genome organization. Both hammerhead (HHRz) and hairpin (HpRz) ribozymes were detected in all the ambiviricot sequences. Armillaria viruses were compared by phylogenetic analysis and confirmed their specific host by direct RT-PCR, which showed the occurrence of interspecies transmission.

The concept of vulnerable carotid plaques is pivotal in understanding the pathophysiology of is-chemic stroke secondary to large-artery atherosclerosis. At the macroscopic evaluation, vulnera-ble plaques are characterized by one or more of the following features: microcalcification, neo-vascularization, lipid-rich necrotic core (LRNC), intraplaque hemorrhage (IPH), thin fibrous caps, plaque surface ulceration, huge dimension determining stenosis, and plaque rupture. Recognizing these macroscopic characteristics is crucial to estimating the risk of cerebrovascular events, also in case of non-significant (less than 50%) stenosis. Inflammatory biomarkers, such as cytokines and adhesion molecules, lipid-related markers like oxidized low-density lipoprotein (LDL), and proteolytic enzymes capable of degrading extracellular matrix components are among the key molecules that will be scrutinized for their associative roles in plaque instability. Through their quantification and evaluation, these biomarkers reveal intricate molecular cross-talking govern-ing plaque inflammation, rupture potential, and thrombogenicity. Current evidence demonstrates that plaque vulnerability phenotypes are multiple and heterogeneous and are associated with as many highly complex molecular pathways that determine the activation of an immune-mediated cascade that culminates in thromboinflammation. This narrative review provides a comprehen-sive analysis of current knowledge on molecular biomarkers expressed by symptomatic carotid plaques. It explores the association of these biomarkers with the structural and compositional at-tributes that characterize vulnerable plaques

Knowledge is not merely the characteristic of certain pieces of information being deemed true and justifiable, but rather a process. Within an axiomatic system, the set of information that reaches the status of knowledge forms a universe that provides an interpretation of the reality one aims to describe. Language exchanges information through communication and the process of justification (exchange of information), if it consistently justifies a particular piece of information, it elevates it to knowledge. In this paper, we provide a framework for knowledge within an axiomatic model where we define the fundamental elements of communication and grammar. As a result, the structure dictates that knowledge, to be considered as such, must be investigateable by anyone who utilizes it in its entirety, both in form and content. Ultimately, the development of knowledge processing, referred to as \emph{knowmatics}, represents a clear evolution of the computer science that processes information (subset of knowledge).

A single water pollution by dyes is a huge environmental problem; there is a necessity to produce new decolorization methods that are effective, cost-attractive, and suitable in industrial use. Magnetic chitosan polymers offer the advantage of easy separation from the dye solution. In this work, a Chitosan-Fe polymer was synthesized, characterized, and tested for removal of the azo dye Direct Red 83:1 from water. The polymer was characterized in terms of surface morphology (SEM), elemental analysis (EA), differential scanning calorimetry (DSC), thermal gravimetric analysis (TGA), infrared spectrophotometry (IR), X-ray powder diffraction (XRD) and density functional theory (DFT). The reported results hint that 1 g and pH 7.0 were the best conditions to carry out both kinetic and isotherm models. To understand the adsorption behavior of Direct Red 83:1 onto the polymers, adsorption rate and maximum adsorption capacity were calculated using kinetic tests and isotherm curves respectively. The kinetic data and mechanism of the adsorption process were analyzed by three models and the equilibrium data by three adsorption isotherms. The results indicated that the adsorption process follows a pseudo-second-order kinetics (R2= 0.999) and fits the Temkin isotherm (R2= 0.946), indicating that the adsorption occurred on heterogeneous surfaces. The newly synthesized Chitosan-Fe polymer exhibited good adsorption properties and easy separation of cleaned water.

As the natural resources are getting exhausted, the concept of sustainable development of region has received increasing attentions, especially for resource-dependent cities. In this paper an innovative method that emergy analysis and IPAT (Human Impact Population Affluence Technology) model were combined in order to analyze the quantitative relationship of economic growth and energy consumption and further to evaluate its overall sustainability level. Taiyuan, a traditional, resource-dependent city in China, is selected as the case study region. The main results show that total emergy of Taiyuan increased from 9.023× 1023 sej in 2007 to 9.116× 1023sej in 2014, with 38% reduction on non-renewable emergy and 125% growth on imported emergy. Regional emergy money ratio (RMB) was reduced by 48% from 5.31× 1013sej/$ in 2007 to 2.74× 1013sej/$ in 2014, indicating that the increasing speed of consuming resources and energy was faster than the increase of GDP, and that Taiyuan’s money purchasing power declined. The lower emergy sustainability index (ESI) indicates that Taiyuan was explored and produced large quantities of mineral resources, which puts more stress on the environment as a consequence, and that this is not sustainable in the long run. The IPAT analysis demonstrates that Taiyuan sticks to the efforts of energy conservation and environmental protection, in order to promote regional sustainable development, it is necessary to take an integrated effort. Policy insights suggest that resourceful regions should improve include energy and resource efficiency, optimizing energy and resourceful structure.

Extracellular vesicle (EV) release, which occurs in most eukaryotic cells, has recently been associated with peptidylarginine deiminase (PAD)-driven protein deimination. Evidence points to the involvement of deiminated cytoskeletal proteins and changes in histone deimination. Both PADs and EVs are associated with various pathologies including cancers, autoimmune and neurodegenerative diseases. The elevated PAD expression observed in cancers may contribute to increase in EV shedding observed from cancer cells, contributing to cancer progression. Similarly, elevated PAD expression observed in neurodegenerative diseases may cause increased EV shedding and spread of neurodegenerative EV cargo, contributing to disease progression and pathologies. Pharmacological inhibition of PAD-mediated deimination using pan-PAD inhibitor Cl-amidine, reduced cellular EV release in prostate cancer cells, rendering them significantly more susceptible to chemotherapeutic drugs. Studies on models of central nervous system damage have demonstrated critical functional roles for PADs and neuroprotective effects using PAD inhibitors in vivo, while human neurodegenerative iPSC in vitro models showed evidence of increased protein deimination. Besides using refined PAD inhibitors to selectively manipulate EV biogenesis for novel combination therapies in cancer treatment, we also speculate how EV biogenesis could be targeted via the newly identified PAD-pathway to ameliorate neurodegenerative disease progression.

Brain aging causes a wide variety of changes at molecular and cellular levels, leading to the decline of cognitive functions and increased vulnerability to neurodegenerative disorders. The research aimed at understanding the aging of the brain has made much progress in recent decades. Technological innovations such as single-cell RNA-sequencing (scRNA-seq), proteomic analyses, and spatial transcriptomic analyses have facilitated the research on the dynamic changes occurring within neurons, glia, and other cells along with their impacts on intercellular communication during aging. In this review, we introduce recent trends of how neurons and glia change during aging and discuss the impact of the brain microenvironment such as the blood-brain barrier (BBB).

Background: Acute kidney disease (AKD) is a known risk factor for increased mortality and evolution towards chronic kidney disease (CKD) in adults. The data regarding children are scarce. The purpose of our study is to explore the risk factor for AKD and impact of AKD following an acute kidney injury (AKI) episode in children on mortality and new-onset CKD. Methods: We retrospectively evaluated 736 children with ages between 2 and 18 years old with AKI admitted in a tertiary care hospital from the west Romania over a 9 years period. Results: AKD incidence following AKI was 17%. Patients that developed AKD were older, with higher baseline serum creatinine values, urea, C reactive protein and lower proteins, haemoglobin and sodium. They were more often exposed to sepsis, critical illness and nephrotoxic medications and presented certain susceptibilities like CKD, neoplasia, anemia, heart failure, arterial hypertension, stem cell transplant and female gender. AKI causes that increased the risk of AKD were intrinsic ones. AKD increased the overall mortality by three times. The factors that independently increased the risk of evolution towards CKD were AKD, acute tubular necrosis and higher baseline serum creatinine values. Conclusions: AKD is an independent risk factor for mortality and new-onset CKD in children.

In this paper, a novel permanent magnet-assisted synchronous reluctance machine (PMASynRM) with rare-earth PMs and ferrite magnets is proposed. The performance of PMASynRM is discussed with respected to the different magnet ratio of rare-earth PMs and ferrite magnets. Some characteristics including the flux density, output torque, cogging torque, output power, power factor, torque ripple, loss, efficiency, and demagnetization are calculated by 2-D finite element analysis (FEA). The analysis results show that the excellent performance can be obtained by using hybrid magnet of rare-earth PMs and ferrite magnets with the suitable magnet ratio, and provide some desirable cost-performance trade-off.

Sustainable development of our society demands recycling of metallic wastes which may contain a variety of metallic impurities. Thus, it is desirable to develop a ‘universal’ structure, such as an amorphous alloy tolerating high concentrations of metallic pollutants. We performed ab initio molecular dynamics simulations to manufacture a novel multicomponent amorphous alloy of twenty elements (20e-HEAA). Analysis reveals both topological and chemical compositional short-range ordering in the obtained 20e-HEAA and formation of B-B clusters. Electronic structure calculations reveal the metallic nature of the amorphous alloy with its Fermi level falling in a valley in the density of states. Moreover, the Zn 3d and Cd 4d electrons exhibit a localized nature and thus, should be excluded from the valence electron concentration (VEC) analysis in the study of high-entropy alloys. The information obtained here is helpful to get insight into multicomponent high entropy materials, especially multicomponent high-entropy amorphous alloys as potential circular metals.

A large amount of steel slag (SS) and CO2 are generated in the steelmaking process. The indirect CO2 capture using SS is a promising way for co-treatment of the wastes. Ammonium salt solution is widely used to extract Ca2+ from the SS since it is recyclable. Several works have focus on improving the carbonation rate by adjusting various parameters (e.g., temperature and pH). However, there is little detail information about the associating behaviors and interaction strength between the various ions in the ammonia solution during the carbonation process. In this work, the Ca2+–CO32-–NH4+–Cl-–H2O system was established by using Material studio software. The effects of temperature and concentration of CO32- on CaCO3 growth were explored at the atomic scale by calculating the binding energy, mean square displacement, and diffusion coefficient between particles. Furthermore, the microstructure, bonding characteristics, and occurrence behavior of each particle were studied though molecular dynamics simulation methods. The results showed that with the increase of temperature (20–80 ℃), the binding ability and diffusion coefficients of Ca2+ and CO32- increase in the system, which is beneficial to the formation of CaCO3 clusters. With the increase of the concentration of CO32- (15–25 vol.%), the binding ability and diffusion coefficient of Ca2+ and CO32- in the system are enhanced, which is beneficial to the formation of CaCO3 clusters.

When liquids are confined into nanometer-scale slit, the induced layering-like film structure allows the liquid to sustain non-isotropic stresses and thus being load-bearing. Such anisotropic characteristics of liquid under confinement arise naturally from the liquids’ wave number dependent compressibility that does not need solidification to take place as a prerequisite. In other words, liquids under confinement can still remain fluidity with molecules being (sub-)diffusive. However, the extensively prolonged structural relaxation time can cause hysteresis of stress relaxation of confined molecules in response to the motions of confining walls and thereby yield the quasi-static stress tensor history-dependent. In this work, by means of molecule dynamics, the discrepancy of stress tensor of a highly confined key base-oil component, i.e. 1-decene trimer, is captured after its relaxation from being compressed and decompressed. The results indicate that among the effects (e.g. confinement, molecular structure, and film density) that can potentially affect confined stress tensor, the ordering status of the confined molecules plays a predominant role.

The thermal lubrication of an entangled polymeric liquid in wall-driven shear flows between parallel plates is investigated by using a multiscale hybrid method coupling molecular dynamics and the hydrodynamics (i.e., the synchronized molecular dynamics method). The temperature of the polymeric liquid rapidly increases due to viscous heating once the drive force exceeds a certain threshold value. The rheological properties of the polymeric liquid drastically change at around the critical drive force. In the weak viscous-heating regime, the conformation of polymer chains is dominated by the local shear flow so that the anisotropy of the bond orientation tensor grows as the drive force increases. However, in the large viscous-heating regime, the conformation dynamics is dominated by the thermal agitation of polymer chains so that the bond orientation tensor recovers more uniform and random structures as the drive force increases, even though the local shear flows are further enhanced. Remarkably, these counter-intuitive transitional behaviors give an interesting re-entrant transition in the stress--optical relation, where a linear formalism in the stress--optical relation approximately holds even though each of the macroscopic quantities behaves nonlinearly. The robustness of the linear stress--optical relation is also confirmed in the spatiotemporal evolution at the hydrodynamic level.

By manipulating the various physico-chemical properties of amino acids, design of peptides with specific self-assembling properties has been emerging since more than a decade. In this context, short peptides possessing detergent properties (so-called “peptergents”) have been developed to self-assemble into well-ordered nanostructures that can stabilize membrane proteins for crystallization. In this study, the peptide with “peptergency” properties, called ADA8 extensively described by Zhang et al., is studied by molecular dynamics for its self-assembling properties in different conditions. In water, it spontaneously forms beta sheets with a β barrel-like structure. We next simulated the interaction of this peptide with a membrane protein, the bacteriorhodopsin, in the presence or absence of a micelle of dodecylphosphocholine. According to the literature, the peptergent ADA8 is thought to generate a belt of β structures around the hydrophobic helical domain that could help stabilize purified membrane proteins. Molecular dynamics is here used to challenge this view and to provide further molecular details for the replacement of detergent molecules around the protein. To our best knowledge, this is the first molecular mechanism proposed for ''peptergency''. In addition, our calculation approach should serve as a predicting tool for the design of beta peptergent with diverse amphipathic properties.

Graphene hybrid-filler polymer composites have emerged as prominent materials that revolutionize heavy industries. This study encapsulates the in-depth analysis of different influential factors such as filler/graphene type, aspect ratios, dispersion methods, filler-matrix compatibility, fiber orientation, synergistic effects, different processing techniques, and post-curing conditions that affect the processing, properties and determine the resultant applications of graphene hybrid polymer composites. Additionally, it discusses the substantial role of graphene collaboration with other fillers such as carbon nanotubes, silica, nano-clays and metal oxides to produce functionalized hybrid polymer composites with synergistically enhanced tailored properties to offer solutions for heavy industries like aerospace, automotive, electronics, and energy harvesting. The review concludes with some suggestions and outlook on the future of these composite materials by emphasizing the need for continued research to fully optimize their potential.

The apicomplexan parasite, Sarcocystis neurona causes the degenerative neurological equine protozoal myeloencephalitis (EPM) disease of horses. Due to its host range expansion, S. neurona is an emerging threat that requires close monitoring. In apicomplexans, protein kinases (PKs) have been implicated in a myriad of critical functions such as host cell invasion, cell cycle progression and host immune responses evasion. Here, we used various bioinformatics methods to define the kinome of S. neurona and phylogenetic relatedness of its PKs to other apicomplexans. Further, three-dimensional (3D) homology models for selected S. neurona putative PKs were constructed and evaluated for inhibitor docking. We identified 92 putative PKs clustering within the AGC, CAMK, CK1, CMGC, STE, TKL, aPK and OPK groups. Although containing the universally conserved PKA (AGC group), S. neurona kinome was devoid of PKB and PKC, but contained the six apicomplexan conserved CDPKs (CAMK group). The OPK group was represented by ROPKs 19A, 27, 30, 33, 35 and 37, but was devoid of the virulence-associated ROPKs 5, 6, 18 and 38. Two out of the three S. neurona CK1 enzymes had high sequence similarities to T. gondii TgCK1-α and TgCK1-β and the Plasmodium PfCK1. Docking of four inhibitors onto homology models of putative ROP27 and PKA indicated that inhibition of S. neurona PKs is feasible, but needs to be experimentally tested. The essentiality of apicomplexan PKs makes the elucidation of S. neurona kinome a key milestone for development of novel therapeutics for EPM.