We present a critical analysis of the mechanism of reversible hydrogen evolution reaction based on thermodynamics of hydrogen processes considering atomic and ionic species as intermediates. Clear distinction between molecular hydrogen evolution/oxidation (H$_2$ER and H$_2$OR) and atomic hydrogen evolution/oxidation (HER and HOR) reactions is made. It is suggested that the main reaction describing reversible H$_2$ER and H$_2$OR in acidic and basic solutions is: $\ce{H3O^+} + 2e^- \xrightleftharpoons{\mathrm{({H}_{2}^{+}})_{ad}} \ce{H2} + \ce{OH^-}$, which proceeds via the formation of hydrogen molecular ion as intermediate and its potential is $E^{0}$ = 0.414~V (at the standard Hydrogen electrode SHE scale). We analyse experimentally reported data with models which provide \emph{quantitative} match (R.J.Kriek et al., Electrochem. Sci. Adv. e2100041 (2021)). The experimentally observed slope of $\lg(Current)$ vs. $Potential$ equal to $\approx 30$~mV/dec corroborates the presented model of H$_2$ER/H$_2$ER based on \ce{H2$^+$} as intermediate.

This work communicates laser-plasma experiments in a gaseous mixture of hydrogen and nitrogen. Time-resolved spectroscopy measures the first four Balmer series hydrogen lines together with selected neutral and ionized nitrogen lines. Optical breakdown plasma is generated in a 1:1 hydrogen:nitrogen mixture at ambient temperature and 0.27-atm pressure. Time-resolved spectroscopy records emitted radiation with spatial resolution along the slit height for the H_α, H_β, H_γ, and H_δ lines. For 13 selected time delays from 0.25 µs to 3.25 µs and 0.025 µ gate-widths, micro-plasma diagnostics is evaluated. Of interest are the peak-separation and width of H_δ and width of H_γ for electron densities in the range of 0.1 to 1.0×10¹⁷ cm^-3, and comparisons with H_β and H_α diagnostics. Integral inversions interrogate spatial distributions of the plasma expansion. Applications include laboratory and stellar astrophysics plasma diagnosis.

This work examines atomic and molecular signatures in laser-induced plasma in standard ambient temperature and pressure environments, including background contributions to the spectra that depend on the laser pulse-width. Investigations include solids, gases, and nano-particles. Abel inversions of measured line-of-sight data reveal insight into the radial plasma distribution. For nominal 6 nanosecond laser pulses and for pulse-energies in the range of 100 to 800 mJ, expansion dynamics and turbulence due to shock phenomena are elucidated to address local equilibrium details that are frequently assumed in spatially averaged emission spectroscopy. Chemical equilibrium computations reveal temperature dependence of selected plasma species. Specific interests include atomic hydrogen (H) and cyanide (CN). Atomic hydrogen spectra indicate axisymmetric shell structures and isentropic expansion of the plasma kernel over and above the usual shockwave. The recombination radiation of CN emanates within the first 100 nanoseconds for laser-induced breakdown in a 1:1 CO2:N2 gas mixture when using nanosecond laser pulses to create the micro-plasma. The micro-plasma is generated using 1064 nm, 150 mJ, 6 ns Q-switched Nd:YAG laser radiation. Measurements of the optical emission spectra utilize a 0.64 m Czerny-Turner type spectrometer and an intensified charge-coupled device.

The removal of H2S from natural gas and other gases such as biogas, refinery gases, and coal gas is required because H2S is toxic and corrosive, even in traces. Zeolites are widely used in the removal of H2S from the above-mentioned gases. In this work, we prepared an Ag-exchanged 13X molecular sieve by using different concentrations of AgNO3 to increase its adsorption properties. XRD, SEM, and BET techniques were used to characterize samples. When compared to unmodified 13X, the adsorption capacity of AgII-13X increased by about 50 times. The results of this study suggest that the silver-modified 13X molecular sieve is highly effective at extracting H2S from natural gas.

Phenol is added to acetophenone (methyl phenyl ketone) and to six of its halogenated derivatives in a supersonic jet expansion to determine the hydrogen bonding preference of the cold and isolated 1:1 complexes by linear infrared spectroscopy. Halogenation is found to have a pronounced effect on the docking site in this intermolecular ketone balance experiment. The spectra unambiguously decide between competing variants of phenyl group stacking due to their differences in hydrogen bond strength. Structures where the phenyl group interaction strongly distorts the hydrogen bond are more difficult to quantify in the experiment. For unsubstituted acetophenone, phenol clearly prefers the methyl side despite a predicted sub-kJ/mol advantage which is nearly independent of zero point vibrational energy, turning this complex into a challenging benchmark system for electronic structure methods which include long range dispersion interactions in some way.

An H-bond involves the sharing of a hydrogen atom between an electronegative atom to which it is covalently bound (the donor) and another electronegative atom serving as an electron acceptor. Such bonds represent a critically important geometrical force in biological macromolecules and, as such, have been characterized extensively. H-bond formation invariably leads to a weakening within the acceptor moiety due to the pulling exerted by the donor hydrogen. This phenomenon can be compared to a spring connecting two masses; pulling one mass stretches the spring, similarly affecting the bond between the two masses. Herein, we describe the opposite phenomenon when investigating the energetics of the C-H...O=C bond. This bond underpins the most prevalent protein transmembrane dimerization motif (GxxxG) in which a glycine C-H on one helix forms a hydrogen bond with a carbonyl in a nearby helix. We use isotope-edited FT-IR spectroscopy and corroborating computational approaches to demonstrate a surprising strengthening of the acceptor C=O bond upon binding with the glycine C-H. We show that electronic factors associated with the C-H bond strengthen the C=O oscillator by increasing the s-character of the -bond, lowering the hyperconjugative disruption of the -bond. In addition, a reduction of the acceptor C=O bond’s polarity is observed upon the formation of the C-H...O=C bond. Our findings challenge the conventional understanding of H-bond dynamics and provide new insights into the structural stability of inter-helical protein interactions.

Knowledge of the biological effects of molecular hydrogen (H2), hydrogen gas, is constantly expanding giving a reason for optimism for several healthcare practitioners regarding management of multiple diseases, including socially significant ones. However, mechanisms underlying the biological effects of H2 are still being actively argued. In this review, we focus on mast cells as a potential target for H2 at the level of a specific tissue microenvironment. H2 regulates processing of pro-inflammatory components of the mast cell secretome and their entry into the extracellular matrix; this can significantly affect the state of the integrated-buffer metabolic environment and the structure of the immune landscape of the local tissue microenvironment. The analysis performed highlights several potential mechanisms for developing the biological effects of H2 and offers great opportunities for translating the obtained findings into clinical practice.

In regenerative medicine, mesenchymal stem cells (MSCs) have shown their importance and potential in tissue reconstruction and immune system modification. However, such cells’ potential is often diminished by factors such as oxidative stress, immune rejection, and inadequate engraftment. This review highlights the role of molecular hydrogen (H₂) and cold atmospheric plasma (CAP) as adjunct therapies to improve the effectiveness of MSC therapy has a strong antioxidative and anti-inflammatory action as it quenches reactive oxygen species and positively stimulates the Nrf2 pathway that promotes MSC survival and life. CAP, being a modulated source of ROS and RNS, also assists MSCs by altering the cellular redox balance, thus facilitating cellular adaptation, migration and differentiation. H₂ and CAP in conjunction with each other assist in establishing an ambience favorable for promoting MSC survival and growth ability and reducing the healing time in various pathways such as wound, neuroprotection, and ischemia. Besides these concerns, this review also covers the best administration routes and doses of H₂ and CAP together with MSCs in therapy. This study informs on a novel dual method aimed at improving the outcome of MSC therapy while adding several molecular targets and relevant clinical use concerning these therapies. Research of the future has to deal with bettering these protocols so that the therapeutic benefits can be maximized without long-term implications for clinical applications.In regenerative medicine, mesenchymal stem cells (MSCs) have shown their importance and potential in tissue reconstruction and immune system modification. However, such cells’ potential is often diminished by factors such as oxidative stress, immune rejection, and inadequate engraftment. This review highlights the role of molecular hydrogen (H₂) and cold atmospheric plasma (CAP) as adjunct therapies to improve the effectiveness of MSC therapy has a strong antioxidative and anti-inflammatory action as it quenches reactive oxygen species and positively stimulates the Nrf2 pathway that promotes MSC survival and life. CAP, being a modulated source of ROS and RNS, also assists MSCs by altering the cellular redox balance, thus facilitating cellular adaptation, migration and differentiation. H₂ and CAP in conjunction with each other assist in establishing an ambience favorable for promoting MSC survival and growth ability and reducing the healing time in various pathways such as wound, neuroprotection, and ischemia. Besides these concerns, this review also covers the best administration routes and doses of H₂ and CAP together with MSCs in therapy. This study informs on a novel dual method aimed at improving the outcome of MSC therapy while adding several molecular targets and relevant clinical use concerning these therapies. Research of the future has to deal with bettering these protocols so that the therapeutic benefits can be maximized without long-term implications for clinical applications.

With the rapid increase in diabetes worldwide, the number of patients with diabetic kidney disease (DKD), a complication of diabetes mellitus, is also on the rise. DKD is a major cause of chronic kidney disease progressing to end-stage renal failure; however, current medical treatments for DKD do not achieve satisfactory outcomes. Molecular hydrogen (H2) is an antioxidant that selectively reduces hydroxyl radicals, a reactive oxygen species with a very potent oxidative capacity. H2 was recently shown to exert not only antioxidant, but also anti-inflammatory, cell lethality-regulating, and signal transduction-regulating effects, and is now being applied clinically. Many factors contribute to the onset and progression of DKD, with mitochondrial dysfunction, oxidative stress, and inflammation being strongly implicated. Recent preclinical and clinical trials reported that substances with antioxidant properties may slow the progression of DKD. Therefore, we conducted a literature review on animal models and human clinical trials in which H2 showed efficacy against various renal diseases. This literature review and our previous findings collectively suggest that H2 exerts therapeutic effects in patients with DKD by improving mitochondrial function. Future large-scale clinical studies are needed to confirm these findings.

Prereactive complexes in noncovalent organocatalysis are sensitive to the relative chirality of the binding partners and to hydrogen bond isomerism. Both effects are present when a transiently chiral alcohol docks on a chiral α-hydroxy ester, turning such 1:1 complexes into elementary, non-reactive model systems for chirality induction in the gas phase. With the help of linear infrared and Raman spectroscopy in supersonic jet expansions, conformational preferences are investigated for benzyl alcohol in combination with methyl lactate, also exploring p-chlorination of the alcohol and the achiral homolog methyl glycolate to identify potential London dispersion and chirality effects on the energy sequence. Three of the four combinations prefer barrierless complexation via the hydroxy group of the ester (association). In contrast, the lightest complex shows predominantly insertion into the intramolecular hydrogen bond, like the analogous lactate and glycolate complexes of methanol. The experimental findings are rationalized with computations and a uniform helicality induction in the alcohol by the lactate is predicted, independent on insertion into or association with the internal lactate hydrogen bond. p-Chlorination of benzyl alcohol has a stabilizing effect on association, because the insertion motif prevents a close contact between the chlorine and the hydroxy ester. After simple anharmonicity and substitution corrections, the B3LYP-D3 approach offers a fairly systematic description of the known spectroscopic data on alcohol complexes with α-hydroxy esters.

The human body contains 60-70% of water, depending on age. As a body fluid, it is not only a medium in which physical and chemical processes take place, but it is also one of the active mediators. Water is the richest substance with non-covalent hydrogen bonds. Water molecules, by themselves (in vacuum), are diamagnetic, but when organized into clusters, become diamagnetic or paramagnetic. Also, biomolecules (DNA, collagen, clathrin, and other proteins) have non-covalent hydrogen bonds in their structure. The interaction, as well as signal transmission, between water and biomolecules is achieved through the vibrations of covalent and non-covalent hydrogen bonds, which determine the state and dynamics of conformational changes of biomolecules. Disruptive conformational changes of biomolecules, cells, and tissues lead to their dysfunctionality, so they are a frequent cause of many disorders and diseases. For example, the rearrangement of hydrogen bonding due to mitochondrial disease mutation in cytochrome bc1 disturbs heme bH redox potential and spin state. In order to prevent and repair the dysfunctional conformational changes, a liquid substance was developed based on the second derivative of the C₆₀ molecule (SD-C₆₀), which has classical and quantum properties. Characterization of SD-C₆₀ by UV-VIS-NIR, FTIR, TEM, and AFM/MFM was done and it is shown that SD-C₆₀ water layers generate vibrations with near-zero phase dispersion which are transmitted through the Fibonacci’s water chains to biomolecules. In comparison to previously published SD-C₆₀ derivate (3HFWC, size until 10 nm, and 1-5 water layers), improved formulation (3HFWC-W, size 10-25 nm, and 6-9 water layers) showed multiplied cytotoxic activity against melanoma cell lines of different aggressiveness. Apart from this, the mode of action was preserved and based on induction of senescence rather than cell death. Importantly, high selectivity toward malignant phenotype was detected. Observed effects can be ascribed to a machinery of hydrogen bonds, which are generated in SD-C₆₀ and transmitted through water to biomolecules. This approach may open a new field in science and healthcare - a “water-based nanomedicine”.

The aim of this study was to investigate the relationship between postprandial glycemic levels based on flash-monitoring and the production of intestinal hydrogen (H2) and methane (CH4) gases based on measurement of the amount of these gases in exhaled air. We studied 14 subjects with type 2 diabetes mellitus (T2DM) and 10 individuals without diabetes (control) with 3 food-load tests with simultaneous determination of gases in exhaled air and assessment of glucose levels. In 15.4% of patients with T2DM, there was a non-significant increase in exhaled H2. How-ever, all subjects in the control group had significant increase in exhaled H2. Moreover, the level of H2 in exhaled breathafter food load in patients with T2DM was lower than in normoglycemic subjects. Also, in normoglycemic subjects there was an inverse correlation between maximum glucose rise and maximum H2 in exhaled air after food load. At the same time in patients with T2DM the level of CH4 in exhaled air directly correlated with the parameters of postprandial glycemia in the test with lactulose. Confirmation of a causal relationship between decreased H2 production, increased intestinal CH4 production, and more severe postprandial glycemia may identify new therapeutic targets in the correction of postprandial glycemia in patients with T2DM.

Molecular hydrogen (H2) has been shown to have antioxidant and anti-inflammatory activities that may reduce the development and progression of many diseases. In this study, Hydrogen-rich water (HRW) was obtained by reacting hybrid magnesium-carbon hydrogen storage materials with water. Then the effects of intake of HRW on the activities of xenobiotic-metabolizing enzymes, membrane transporters, and oxidative stress in rats were investigated. Rats were given HRW ad libitum for 4 weeks. Results showed that intake of HRW had no significant effect on the activities of various cytochrome P450 (CYP) enzymes (CYP1A1, 1A2, 2B, 2C, 2D, 2E1, 3A, 4A), glutathione-S-transferase and UDP-glucuronosyltransferase. Except for a slight lower plasma glucose concentration, intake of HRW had no effect on other plasma biochemical parameters in rats. P-Glycoprotein and multidrug resistance-associated protein (Mrp)2 protein expressions in liver were elevated after intake of HRW. However, HRW had no significant effects on glutathione, glutathione peroxidase, or lipid peroxidation in liver. Results from this study suggest that consumption of HRW may not affect xenobiotic metabolism or oxidative stress in liver. However, intake of HRW may increase the efflux of xenobiotics or toxic substances from the liver into bile by enhancing the expression of p-glycoprotein and Mrp2 protein.

Urinary tract infections (UTIs) affect nearly 50% of women in their lifetime. Uropathogenic Escherichia coli (UPEC) expresses F9/Fml pili tipped with protein FmlH that specifically bind to terminal galactoside and galactosaminoside units in glycoproteins in kidney and bladder cells and colonize host tissues. Traditional UTI treatment using excessive antibiotics has led to the rise of various UPEC antibiotic-resistant strains. An alternative therapeutic approach prevents the initial bacterial attachment on the host cells using competitive FmlH-binding inhibitors. In this study, we have used computer-aided drug design techniques to identify novel glycomimetics that are predicted to bind strongly to UPEC FmlH and inhibit binding to host glycoproteins. We have performed in-silico receptor-based and ligand-based scaffold hopping, and molecular docking to predict novel FmlH-binding glycomimetics with high chemical synthesizability. We have replaced the two major scaffolds of the most potent known FmlH-binding ligand to obtain novel compounds. Additionally, we have applied global machine-learning models to predict ADMET properties of the molecules. Compounds with calculated low ADMET risks were subjected to molecular dynamics simulations and detailed investigation of the FmlH–glycomimetic interactions was performed. We have prepared and supplied a library of 58 novel glycomimetics that can be subjected to further biological activity studies.

Mink astrovirus infection remains a poorly understood disease entity, and the aetiological agent itself causes disease with a heterogeneous course, including gastrointestinal and neurological symptoms. This paper presents cases of astrovirus infection in mink from continental Europe. RNA was isolated from the brains and intestines of animals showing symptoms typical of shaking mink syndrome (n = 6). RT-PCR was used to detect astrovirus genetic material, and the reaction products were separated on a 1% agarose gel. The specificity of the reaction was confirmed by sequencing all samples. The presence of astrovirus RNA was detected in each of the samples tested. Sequencing and bioinformatic analysis indicated the presence of the same variant of the virus in all samples. Comparison of the variant with the sequences available in bioinformatics databases confirmed that the Polish isolates form a separate clade, closely related to Danish isolates. The similarity of the Polish variant to those isolated in other countries ranged from 2.4% (in relation to Danish isolates) to 7.1% (in relation to Canadian isolates). Phylogenetic relationships between variants appear to be associated with the geographic distances between them. To our knowledge, this work describes the first results on the molecular epidemiology of MAstV in continental Europe. The detection of MAstV in Central Europe indicates the need for further research to broaden our understanding of the molecular epidemiology of MAstV in Europe.

Herein, we are proposing two chalcone molecules, (E)-1-(4-methoxyphenyl)-3-(p-tolyl) prop-2-en-1-one and (E)-3-(4-hydroxyphenyl)-1-(2,4,6-trihydroxyphenyl) prop-2-en-1-one, based on the anticancer bioactive molecule Xanthohumol, which are suitable for further in vitro and in vivo studies. Their ability to create stable complexes with the antiapoptotic X-linked IAP (XIAP) protein makes them promising anticancer agents. The calculations were based on ligand-based and structure-based virtual screening combined for the pharmacophore built. Additionally, the structures passed Lipinski's rule for drug use, and their reactivity was confirmed using density functional theory studies. The candidates were chosen between 10639400 compounds, and the docking protocols were evaluated using molecular dynamics simulations.

Infection-related cardiovascular diseases (CVDs) represent a significant health challenge, necessitating novel therapeutic strategies targeting key receptors involved in inflammation and infection. Antimicrobial peptides (AMPs) offer a promising approach, potentially disrupting pathogenic processes. This study aimed to investigate the efficacy of AMPs as therapeutic agents by examining their interactions with critical CVD-related receptors. A comprehensive computational approach was utilized to assess the interactions between AMPs and receptors associated with CVDs. Molecular docking studies were conducted to evaluate AMP binding to target receptors: ACE2, CRP, MMP9, NLRP3, and TLR4. The top-performing AMPs were further analyzed through 100 ns molecular dynamics (MD) simulations, and their binding affinities were quantified using MM/PBSA calculations. The analysis revealed that Tachystatin, Pleurocidin, and Subtilisin A exhibit strong binding affinities to key CVD-related receptors, including ACE2, CRP, and MMP9. These AMPs demonstrated the potential for disrupting receptor-peptide interactions critical to infection and inflammation. MD simulations confirmed the stability of AMP-receptor complexes, with MM/PBSA calculations showing significant binding energies. Future directions include conducting in vitro and in vivo studies to validate the therapeutic efficacy and safety of these AMPs in clinical settings.

The inhibition of kynurenine aminotransferase II (KATII) is a target for treating several diseases with an excess of kynurenic acid (KYNA). KATII inhibitors exist, but they have adverse side effects, given their irreversible inhibition mechanism. This study aimed to discover potent and safe inhibitors of KATII by using computational and in vitro methods. We identified herbacetin and (-)-Epicatechin as reversible KATII inhibitors from a small library of natural compounds, whose Glide docking scores (-8.66 kcal mol-1 and-8.16 kcal mol-1, respectively) and MM/GBSA binding energies (-50.30 kcal mol-1 and -51.35 kcal mol-1, respectively), which are higher scores than the standard inhibitor, PF-04859989 (-7.12 kcal mol-1 and -38.41 kcal mol-1, respectively). The ADMET analysis showed that herbacetin and (-)-Epicatechin have suitable pharmacokinetic parameters, are moderately bioavailable, and have a safe toxicity profile. Their KATII inhibitory activity is in the range of IC50 values of 5.98± 0.18 µM and 8.76 ± 0.76 µM, respectively. Furthermore, the MTT assay for cell toxicity showed that the two compounds do not alter HepG2 cell viability at their required concentration for KATII inhibition. These results show that herbacetin and(-)-Epicatechin are suitable for KATII inhibition and appropriate candidates for further development of KATII inhibitors.

Cellular aging is considered as one of the main factors leading to female infertility. The process called “inflammaging”, which was proposed by C. Franceschi in 2014 is characterized by chronic asymptomatic inflammation state resulting to whole body aging and age-associated diseases. In this work the inflammaging model was developed; markers expression associated with SASP (senescence associated secretory phenotype) and cell aging were verified in endometrial cells; the correlation analysis with the known inflammaging markers was done; as well as the age dynamic of markers expression in the groups of young reproductive age and late reproductive age was evaluated. For the first time, the comparative study of the expression of the key signalling molecules involved in the mechanisms of the aging was carried out: IL-6, IL-8, IL-1a, MMP3, SIRT1,6, TERF-1, CALR in endometrial cells in the age aspect and under the influence of the genotoxic stress. The SASP phenotype characterizing the phenomenon of inflammaging can be significantly expanded by the including of the studied molecules - SIRT-1, SIRT-6, TERF, CALR, that we demonstrated for the first time. These signalling molecules may be new therapeutic targets for the age-associated female reproductive system diseases treatment.

Leishmaniasis is a group of neglected infectious, non-contagious diseases caused by parasites of the Leishmania genus which affects millions of people worldwide. The drugs available for its treatment, especially in cases of visceral leishmaniasis, are old, outdated, have serious side effects and in many cases are ineffective due to the emergence of resistant strains. In this context, the search for new candidates for more effective and safer drugs which can become new therapeutic alternatives is a constant need. In this work, 10 new chalcones were synthesized in high yields through Claisen-Schmidt condensation and were evaluated in vitro against Leishmania infantum promastigotes and amastigotes. CP04 and CP06 compounds were the most promising, respectively showing IC50 values = 13.64 ± 0.25 and 11.19 ± 0.22 against promastigotes, and IC50 = 18.92 ± 0.05 and 22.42 ± 0.05 against amastigotes. The compounds did not show cytotoxicity in erythrocytes, showing selectivity indexes greater than 74 and 90. Molecular docking studies carried out under the sterol 14-alpha demethylase (CYP-51) (PDB: 3L4D) and trypanothione reductase (PDB: 5EBK) enzymes from L. infantum evidenced the great affinity of the CP04 molecule for these targets, showing Moldock score values of -94.0758 and -50.5692 KJ/mol-1, with these values being much lower than the co-crystallized ligand energies. Molecular dynamics simulations demonstrated great stability of CP04 in its targets, thus indicating a real interaction possibility. Our findings confirm the potential that chalcones have as candidates for anti-leishmania drugs, and suggest that the CP04 compound may become a promising drug candidate to design and develop novel chalcone-based analogs for leishmaniosis therapy.

Oleanolic acid (OA) being the ubiquitous compound in the plant kingdom is studied for both the neuroprotective and neurotoxic properties. The mechanism of acetylcholinesterase (AChE) inhibitory potential of OA is investigated using the molecular dynamic simulations (MD) and docking as well as membrane-like and in vitro tests. There are also applied SH-SY5Y human neuroblastoma cells as well as in vivo on the zebrafish tests. The inhibitory potential towards the AChE enzyme is examined using the TLC-bioautography assay (the IC50 value is 9.22 μM). The CH- interactions between the central fragment of the ligand molecule and the aromatic cluster created by the His440, Phe288, Phe290, Phe330, Phe331, Tyr121, Tyr334, Trp84 and Trp279 side chains are observed. The in vitro tests carried out on the SH-SY5Y cells indicates that the viability rate is reduced to 71.5%, 61%, and 43% at the concentrations of 100 µg/mL, 300 µg/mL, and 1000 µg/mL, respectively after 48 hours of incubation whereas cytotoxicity against the tested cell line with the IC50 value is 714.32 ± 32.40 µg/mL. The in vivo tests on the zebrafish prove that there is no difference between the control and experimental groups as regards the mortality rate and morphology (P>0.05).

Molecular classifications of colorectal cancer (CRC) are benefitting cancer research by providing insights into subtype-specific disease prognosis and better therapeutic intervention. So far different conventional DNA markers such as microsatellite instability (MSI), CpG island methylator phenotype (CIMP), chromosomal instability (CIN), and BRAF and KRAS mutations have been used to classify CRC patients but have not shown promising prognostic values. Here, for the first time, we show classification of CRC tumors from Saudi Arabian patients based on gene expression profile (GEP). An existing method of CRC subtyping has been applied to the GEP of tumors from Saudi CRC patients. Survival analysis was carried out on predicted CRC subtypes. In-silico functional analyses were conducted on the gene signature used for subtype prediction. The predicted subtypes showed distinct but statistically insignificant overall survival distribution (log-rank test, p = 0.069). Comparison of predicted subtypes in Saudi CRC patients with that of the French one showed significant dissimilarity in the two populations (Chi-square test, p = 0.0091). Functional analyses of the gene signature used for subtyping suggest their association with “cancer” and “gastrointestinal diseases”. Most of the signature genes were found differentially expressed in CRC tumors compared to adjacent normal tissues. Such a classification framework might help improve the treatment of colorectal cancer patients.

Amorphous solid dispersions are considered a promising formulation strategy for the oral delivery of poorly soluble drugs. The limiting factor for the applicability of this approach is the physical (in)stability of the amorphous phase in solid samples. Minimizing the risk of reduced shelf life for a new drug by establishing a suitable excipient/polymer-type from first principles would be desirable to accelerate formulation development. Here we perform Molecular Dynamics simulations to determine properties of blends of eight different polymer-small molecule drug combinations for which stability data is available from a consistent set of literature data. We calculate thermodynamic factors (mixing energies) as well as mobilities (diffusion rates and roto-vibrational fluctuations). We find that either of the two factors, mobility and energetics, can determine the relative stability of the amorphous form for a given drug. Which factor is rate limiting depends on physico-chemical properties of the drug and the excipients/polymers. The methods outlined here can be readily employed for an in-silico pre-screening of different excipients for a given drug to establish a qualitative ranking of the expected relative stabilities, thereby accelerating and streamlining formulation development.

This study employs a comprehensive computational approach to identify and evaluate potential inhibitors of gamma D crystallin (CRYGD), a key protein implicated in cataract formation. A library of compounds, including natural products and FDA-approved drugs, was screened using molecular docking simulations and binding energy calculations. The top-scoring compounds were further assessed for drug-likeness properties and potential toxicity. Molecular dynamics simulations were conducted on the most promising candidates to investigate their interactions with CRYGD at the atomic level. The study identified CNP0185903, CID 3741, and CNP0151361 as the most promising compounds, with CNP0151361 demonstrating the most favorable combination of binding affinity, structural stabilization, and conformational modulation of CRYGD. These findings provide a foundation for the development of novel pharmacological interventions for cataract treatment, potentially offering alternatives or complements to current surgical approaches.

Background: Selectively targeting dopamine receptors has been a persistent challenge in the last years for the development of new treatments to combat the large variety of diseases evolving these receptors. Although, several drugs have been successfully brought to market, the subtype-specific binding mode on a molecular basis has not been fully elucidated. Methods: Homology modeling and molecular dynamics were applied to construct robust conformational models of all dopamine receptor subtypes (D1-like and D2-like receptors). Fifteen structurally diverse ligands were docked to these models. Contacts at the binding pocket were fully described in order to reveal new structural findings responsible for DR sub-type specificity. Results: We showed that the number of conformations for a receptor:ligand complex was associated to unspecific interactions > 2.5 Å and hydrophobic contacts, while the decoys binding energy was influenced by specific electrostatic interactions. Known residues such as 3.32Asp, the serine microdomain and the aromatic microdomain were found interacting in a variety of modes (HB, SB, π-stacking). Purposed TM2-TM3-TM7 microdomain was found to form a hydrophobic network involving Orthosteric Binding Pocket (OBP) and Secondary Binding Pocket (SBP). T-stacking interactions revealed as especially relevant for some large ligands such as apomorphine, risperidone or aripiprazole. Conclusions: This in silico approach was successful in showing known receptor-ligand interactions as well as in determining unique combinations of interactions, key for the design of more specific ligands.

Background: Pancreatic ductal adenocarcinoma (PDAC) is the third leading cause of death by cancer worldwide. Mostly diagnosed with locally advanced or metastatic disease, patients lack treatment options. Gene alterations (GAs) are frequently observed in PDAC; some of which are considered for molecular targeted therapies (MTTs), with potential clinical benefits and im-proved outcomes. However, the applicability of molecular profiling (MP) for precision medicine in PDAC remains to be demonstrated. Methods: We conducted a retrospective analysis of all patients, aged ≥18 years with histologi-cally confirmed PDAC, who underwent tumor MP between 2010 and 2020 in our institution as part of personalized medicine trials. The primary study endpoint was overall survival (OS), and (minimal follow-up was 6 months after MP. Results: Of 115 eligible patients, MP was successful in 102 patients (89%). KRAS mutations were the most frequent GAs, mostly G12D. Based on ESCAT classification, actionable GAs were found in 29 patients (28%), involving mainly BRCA1 or BRCA2 (5 [18%]), HER2 (5 [18%]), MTAP (5 [18%]), and FGFR (3 [11%]). Only 12 of these 29 patients (41%, or 10% of the whole population) received MTTs, with a median progression-free survival of 1.6 months. Median OS was 19 months in patients with actionable GAs treated with MTTs (n=12 [11.8%]), 14 months in patients with actionable GAs treated with standard therapies (n=17 [16.7%]), and 17 months in patients without actionable GAs treated with standard therapies (n=73 [71.5%]; p=0.26). The absence of liver metastases was associated with better OS (HR=0.471, p=0.01). The highest OS following MTT was observed in patients with BRCA mutations treated with olaparib. Interpretation: Actionable GAs were found in more than a quarter of patients with advanced PDAC. Overall, targeting actionable GAs with MTTs was not associated with improved OS in this retrospective study with limited patient numbers. However, selected GA/MTT combinations (e.g., BRCA mutations/olaparib) were associated with a better outcome.

A series of 2-phenylamino-3-acyl-1,4-naphtoquinones were evaluated regarding their in vitro antiproliferative activities using DU-145, MCF-7 and T24 cancer cells. Such activities were discussed in terms of molecular descriptors like half-wave potentials, hydrophobicity and molar refractivity. Compounds 4 and 11 display the highest antiproliferative activity against the three cancer cells, therefore, they were subject to further studies. The in silico prediction of drug likeness, using pkCSM and SwissADME explorer online, shows that compound 11 is a suitable lead molecule to be developed. Furthermore, the expression of some key genes was studied in DU-145 cancer cells. They include genes involved in apoptosis (Bcl-2), tumor metabolism regulation (mTOR), redox homeostasis (GSR), cell cycle regulation (CDC25A), cell cycle progression (TP53), epigenetic (HDAC4), cell-cell communication (CCN2) and inflammatory pathways (TNF). Compound 11 displays an interesting profile because among these genes, mTOR was significantly less expressed as compared to control conditions. Molecular docking show that compound 11 has good affinity with mTOR, unraveling a potential inhibitory effect on this protein. Due to the key role of mTOR on tumor metabolism, we suggest that impaired DU-145 cells proliferation by compound 11 is caused by a reduced mTOR expression (less mTOR protein) and inhibitory activity on mTOR protein.

Melatonin is a hormone used for treating several disorders. However, its oral administration is problematic due to the variable absorption and extensive first-pass metabolism. The topical application of this drug does not present these disadvantages and can be used for therapeutic and skin protection purposes. The adsorption and release of melatonin in the montmorillonite clay was investigated experimentally and theoretically. The drug-clay interaction products were prepared and characterized, showing a modified and improved diffusion and release of melatonin. The crystal structure and spectroscopic properties of melatonin and melatonin-water-mineral materials were studied by molecular modeling, finding that the adsorption of this drug is energetically favourable, and the results were consistent with the experimental data. Ab initio molecular dynamics simulations with the metadynamics methodology showed the release of melatonin from the confined interlayer nanospace of montmorillonite to the liquid bulk water with a 20 kcal/mol free energy barrier.

Salmonella, especially invasive non-typhoidal Salmonella (iNTS) are responsible for developing various invasive diseases, and possess higher mortality rate, due to their higher antibiotic resistance profile than the other bacteria. Therefore, the present study was concerned to develop epitope based peptide vaccine against iNTS species as a successive and substitute protective measures. The study considered comprehensive Immunoinformatic approaches, followed by molecular docking and molecular dynamics simulation to predict the efficient vaccine candidate T cell and B cell epitopes, based on the outer membrane proteins. The study identified two best epitopes YGIFAITAL and KVLYGIFAI from total iNTS outer membrane proteins, which showed higher immunity, non-allergenicity, non-toxicity and also showed higher conservancy and population coverage values. Both epitopes showed higher binding affinity and stability towards HLA-C* 03:03. The MM-PBSA binding free energy showed the YGIFAITAL epitope binds more tightly with both MHC-I and MHC-II molecules. The total contact, H-bond analysis and RMSF results also validate the efficiency of these epitopes as vaccine candidate. The projected B cell epitopes AAPVQVGEAAGS, TGGGDGSNT and TGGGDGSNTGTTTT showed higher antigenicity. Overall, the study concluded that these epitopes can be considered as the potential vaccine candidate to make a successive vaccine against iNTS species. However, this result further needs to be validate by wet lab research to make successive vaccine with these projected epitopes.

Ubiquitin-specific protease 14 (USP14) is a member of the Deubiquitinating enzymes (DUBs) involved in disrupting the regulation of the ubiquitin-proteasome system, responsible for the degradation of impaired and misfolded proteins which is an essential mechanism in eukaryotic cells. The involvement of USP14 in cancer progression and neurodegenerative disorders has been reported. Thereof USP14 is a prime therapeutic target; hence, designing efficacious inhibitors of USP14 is central in curbing these conditions. Herein, we relied on structural bioinformatics methods incorporating molecular docking, Molecular Mechanics Generalized Born Surface Area (MM-GBSA), Molecular dynamics simulation (MD simulation) and ADME to identify potential allosteric USP14 inhibitors. A library of over 733 compounds from the PubChem repository with >90% match to the IU1 chemical structure was screened in a multi-step framework to attain prospective drug-like inhibitors. Two potential lead compounds (CID 43013232 and CID 112370349) were shown to record better binding affinity compared to IU1, but with subtle difference to IU1-47, a 10-fold potent compound when compared to IU1. The stability of the lead molecules complexed with USP14 was studied via MD simulations. The molecules were found to be stable within the binding site throughout the 50ns simulation time. Moreover, the protein-ligand interactions across the simulation run time suggest Phe331, Tyr476, and Gln197 as crucial residues for USP14 inhibition. Furthermore, in silico pharmacological evaluation revealed the lead compounds as pharmacological sound molecules. Overall, the methods deployed in this study revealed two novel candidates that may show selective inhibitory activity against USP14, which could be exploited to produce potent and harmless USP14 inhibitors.

Molecular dynamics(MD) simulations are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals the solvent structure and preferential interaction sites (hot spots) on the protein surface. This information is largely transferable across all classes of protein ligands (from water to drugs) and can be used very effectively to understand ligand recognition and improve the predictive capability of well-established methods such as molecular docking. MD simulations for protein and drug or drug-like compounds are now being used but are still computationally expensive and can only be applied to specific cases. On the other hand, MDmix simulations can now be used in SBDD and we will describe the latest developments and implementations. We expect to see an increase in the application of these techniques to a plethora of protein targets to identify new drug candidates with the advent of new tools and faster computers.

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.

The growing interest in Kv7.2/7.3 agonists originates from the involvement of these channels in several brain hyperexcitability disorders. In particular, Kv7.2/7.3 mutants have been clearly associated with epileptic encephalopathies (DEE) as well as with a spectrum of focal epilepsy disorders, often associated with developmental plateauing or regression. Nevertheless, there is a lack of available therapeutic options, considering that retigabine, the only molecule used in clinic as broad-spectrum Kv7 agonist, has been withdrawn from the market in the late 2016. This is why several efforts have been made both by both academia and industry in the search for suitable chemotypes acting as Kv7.2/7.3 agonists. In this context, in silico methods have played a major role, since the precise structures of different Kv7 homotetramers have been only recently disclosed. In the present review the computational methods used for the design of Kv.7.2/7.3 small molecule agonists and the underlying medicinal chemistry are discussed in the context of their biological and structure-function properties.

Innovations in drug delivery systems are crucial for enhancing therapeutic efficiency. Our research presents a novel approach based on Electro-fluid dynamic atomization (EFDA) to fabricate core-shell monophasic particles (CSMp) from sodium alginate blends of varying molecular weights. This study explores the morphological characteristics of these particles in relation to material properties and process conditions, highlighting their potential in drug delivery applications. A key aspect of our work is the development of a mathematical model that simulates the release kinetics of small molecules, specifically sodium diclofenac. By assessing the diffusion properties of different molecules and gel formulations through transport and rheological models, we have created a predictive tool for evaluating the efficiency of these particles in drug delivery. Our findings underscore two critical independent parameters for optimizing drug release: the external shell thickness and the diffusivity ratios within the dual layers. This allows for precise control over the timing and intensity of the release profile. The study advances our understanding of EFDA in the fabrication of CSMp and offers promising avenues for enhancing drug delivery systems by tailoring release profiles through particle characteristic manipulation.

The type III receptor tyrosine kinase FLT3 is a pivotal kinase for hematopoietic progenitor cell regulation, with significant implications in acute myeloid leukemia (AML) through mutations like internal tandem duplication (ITD). This study delves into the structural intricacies of FLT3, the roles of activation loop mutants, and their interaction with tyrosine kinase inhibitors. Coupled with this, the research leverages molecular contrastive learning and protein language modeling to examine interactions between small molecule inhibitors and FLT3 activation loop mutants. Utilizing the ConPLex platform, over 5.7 million unique FLT3 activation loop mutants-small molecule pairs were analyzed. The binding free energies of three inhibitors were assessed, and cellular apoptotic responses were evaluated under drug treatments. Notably, the introduction of the Xepto50 scoring system provides a nuanced metric for drug efficacy. The findings underscore the modulation of molecular interactions and cellular responses by Y842 mutations in FLT3-KD, highlighting the need for tailored therapeutic approaches in FLT3-ITD-related malignancies.

The process of viral integration into the host genome is an essential step of the HIV-1 life cycle. The viral Integrase (IN) enzyme catalyses integration. IN is an ideal therapeutic enzyme targeted by several drugs; raltegravir (RAL), elvitegravir (EVG), dolutegravir (DTG) and bictegravir (BIC) having been approved by the USA Food and Drug Administration (FDA). Due to high HIV-1 diversity, it is not well understood how specific naturally occurring polymorphisms (NOPs) in IN may affect the structure/function and binding affinity of Integrase Strand Transfer Inhibitors (INSTIs). In this study, we applied computational methods of molecular modelling and docking to analyse the effect of NOPs on the full-length IN structure and INSTI binding. We identified 16 NOPs within the Cameroonian derived CRF02_AG IN sequences and further identified 17 NOPs within HIV-1C South African sequences. The NOPs in the IN structures did not show any effect on INSTI binding. INSTIs displayed similar binding affinities to each IN structure. All INSTIs are clinically effective against diverse HIV-1 strains from INSTI treatment naïve populations. This study supports the use of second-generation INSTI DTG as part of first-line combination antiretroviral therapy (cART) regimens, due to DTG possessing a stronger genetic barrier to the emergence of drug resistance.

The increasingly common occurrence of antibiotic-resistant bacteria has become an urgent public health issue. There are currently some infections without any effective treatment, which require new therapeutic strategies. An attractive alternative is the design of compounds capable of disrupting bacterial communication known as quorum sensing (QS). In gram-negative bacteria, such communication is regulated by acyl-homoserine lactones (AHLs). QS allows bacteria to proliferate before expressing virulence factors. Our group previously reported that hexyloxy phenyl imidazoline (9) demonstrated 71% inhibitory activity of QS at 100 µM (IC50=90.9 µM) in Chomobacterium violaceum, a gram-negative bacterium. The aim of the present study was to take 9 as a lead compound to design and synthesize three 2-imidazolines (13-15) and three 2-oxazolines (16-18), to be evaluated as quorum sensing inhibitors on C. violaceum CV026. We were looking for compounds with a higher affinity towards the Cvi receptor of this bacterium and the ability to inhibit QS. The binding mode of the test compounds on the Cvi receptor was explored with docking studies and molecular dynamics. It was found that 8-pentyloxyphenyl-2-imidazoline 13 reduced the production of violacein (IC50=56.38 µM) without affecting bacterial growth, suggesting inhibition of quorum sensing. Indeed, compound 13 is apparently one of the best QS inhibitors known to date. Molecular docking revealed the affinity of compound 13 for the orthosteric site of the N-hexanoyl homoserine lactone (C₆-AHL) on the CviR protein. Ten aminoacid residues in the active site of C₆-AHL interacted with 13, and 7 of these are the same as those interacting with AHL. Contrarily, 8-octyloxyphenyl-2-imidazoline 14, 8-decyloxyphenyl-2-imidazoline 15 and 9-decyloxyphenyl-2-oxazoline 18 bound only to an allosteric site and thus did not compete with C₆-AHL for the orthosteric site.

Flavonoids are widely recognized as natural polydrugs, given their anti-inflammatory, antioxidant, sedative and antineoplastic activities. Recently, different studies have shown that flavonoid have the potential to inhibit BET bromodomains. Previous reports suggest that flavonoids are putative inhibitors of the ZA channel due to their orientation and interactions with P86, V87, L92, L94 and N140. Herein, a comprehensive characterization of the binding mode of the biflavonoid amentoflavone and fisetin is discussed. To this end, both compounds were docked with BRD4 using four docking programs. Results were post-processed with protein-ligand interaction fingerprints. To gain further insights into the binding mode of the two natural products, docking results were further analyzed with molecular dynamics. Results showed that amentoflavone makes numerous contacts in the ZA channel, as previously described for flavonoids and kinase inhibitors. It was also found that amentoflavone can potentially make contacts with non-canonical residues for BET inhibition. Most of these contacts were not observed with fisetin. Based on these results, amentoflavone was tested for BRD4 inhibition, showing activity in the micromolar range. This work may serve as basis for scaffold optimization and further characterization of flavonoids as BET inhibitors.

Caveolin-1 (Cav-1) is 22 kDa caveolae protein, acts as a scaffold within caveolar membranes. It interacts with alpha subunits of G-protein and thereby regulates their activity. Earlier studies reported elevated or up-regulated levels of caveolin-1 in the serum of prostate cancer patients. Secreted Cav-1 promotes angiogenesis, cell proliferation and anti-apoptotic activities in prostate cancer patients. Cav-1 upregulation is mainly related to prostate cancer metastasis. Keeping above facts in view, the present study was designed to explore Cav-1 as a target for prostate cancer therapy using computational approach. Molecular docking, structural base molecular modelling and molecular dynamics simulations were performed to investigate Cav-1 inhibitors. A predictive model was generated and validated to establish a stable structure. ZINC database of biogenic compounds was used for induced fit docking (IFD) and high throughput virtual screening. The H-bond interactions of the compounds with active site residues of Cav-1 were estimated by IFD and 100 ns long molecular dynamic simulations. The reported compounds showed significant binding and thus can be considered as potent therapeutic inhibitors of Cav-1. This study provides a valuable insight into biochemical interactions of Cav-1 for therapeutic applications and warrants for experimental validation of the predicted ‘active(s)’.

The bipartite entanglement is comprehensively investigated in the mononuclear molecular complex (Et3NH)[Ni(hfac)2L], where HL denotes 2-(2-hydroxy-3-methoxy-5-nitrophenyl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazol-3-oxide-1-oxyl and hfacH stands for hexafluoroacetylacetone. From the magnetic point of view, the molecular compound (Et3NH)[Ni(hfac)2L] consists of an exchange-coupled spin-1 Ni2+ magnetic ion and a spin-12 nitronyl-nitroxide radical substituted nitrophenol. The nickel-radical molecular complex affords an experimental realization of a mixed spin-(12, 1) Heisenberg dimer with a strong antiferromagnetic exchange coupling J/kB= 505 K and two distinct g-factors gRad=2.005 and gNi=2.275. By adopting this set of magnetic parameters we demonstrate that the Zeeman splitting of a quantum ferrimagnetic ground-state doublet due to a weak magnetic field may substantially reinforce the strength of bipartite entanglement at low temperatures. The molecular compound (Et3NH)[Ni(hfac)2L] maintains sufficiently strong thermal entanglement even at room temperature, vanishing only above 546 K. Specifically, the thermal entanglement in the nickel-radical molecular complex retains approximately 40% of the maximum value corresponding to perfectly entangled Bell states at room temperature, which implies that this magnetic compound provides suitable platform of a molecular qubit with potential implications for room-temperature quantum computation and quantum information processing.

Drug nanosuspensions offer a promising approach to improve bioavailability for poorly soluble drug candidates. Such formulations often necessitate the inclusion of an excipient to stabilize the drug nanoparticles. However, a rationale to choose the correct excipient for a given drug candidate remains unclear. To gain molecular insight into formulation design, this work first utilizes molecular dynamics simulation to computationally investigate drug-excipient interactions for a number of combinations that have been previously studied experimentally. We find that hydrophobic interactions drive excipient adsorption to drug nanoparticles, and that the fraction of polar surface area serves as a predictor for experimental measurements of nanosuspension stability. To test these ideas prospectively, we applied our model to an uncharacterized drug compound, GDC-0810. Our simulations predicted that salt forms of GDC-0810 would lead to more stable nanosuspensions than the neutral form; therefore, we tested the stability of salt GDC-0810 nanosuspensions and found that the salt form readily formed nanosuspensions even without excipient. To avoid computationally expensive simulations in the future, we extended our model by showing that simple, two-dimensional properties of the single drug molecule can be used to rationalize nanosuspension design without simulation. In all, our work demonstrates how computational tools can provide molecular insight into drug-excipient interactions, and aid in rational formulation design.

We extend a previously developed model for the Stark resonances of the water molecule. The method employs a partial-wave expansion of the single-particle orbitals using spherical harmonics. To find the resonance positions and decay rates, we use the exterior complex scaling approach which involves the analytic continuation of the radial variable into the complex plane and yields a non-hermitian Hamiltonian matrix. The real part of the eigenvalues provides the resonance positions (and thus the Stark shifts), while the imaginary parts $-\Gamma/2$ are related to the decay rates $\Gamma$, i.e., the full-widths at half-maximum of the Breit-Wigner resonances. We focus on the three outermost (valence) orbitals, as they are dominating the ionization process. We find that for forces directed along the three Cartesian co-ordinates, the fastest ionizing orbital always displays a non-monotonic Stark shift. For the case of fields along the molecular axis we show results as a function of the number of spherical harmonics included ($\ell_{\max}=3,4$). Comparison is made with total molecule resonance parameters from the literature obtained with Hartree-Fock and coupled cluster methods.

Biological membranes are fascinating. Santiago Ramón y Cajal, who received the Nobel prize in 1906 together with Camillo Golgi for their work on the nervous system, wrote “[..]in the study of this membrane[..] I felt more profoundly than in any other subject of study the shuddering sensation of the unfathomable mystery of life”[1]. The visualization and conceptualization of these biological objects have profoundly shaped many aspects of modern biology, drawing inspiration from experiments, computer simulations, as well as from the imagination of scientists and artists. The aim of this review is to provide a fresh look on current ideas of biological membrane organization and dynamics by discussing selected examples across fields [1] The full quotation is “I must not conceal the fact that in the study of this membrane I for the first time felt my faith in Darwinism (hypothesis of natural selection) weakened, being amazed and confounded by the supreme constructive ingenuity revealed not only in the retina and in the dioptric apparatus of the vertebrates but even in the meanest insect eye. There, in fine, I felt more profoundly than in any other subject of study the shuddering sensation of the unfathomable mystery of life.” from the autobiography Recollections of My Life.

The formation of atroposelective biaryl compounds in plants and fungi is well understood; however, the polyketide aglycone synthesis and dimerization in bacteria remain unclear. Thus, the biosynthetic gene cluster (BGC) responsible for antibacterial setomimycin production from Streptomyces nojiriensis JCM3382 was examined in comparison with the BGCs of spectomycin, julichromes, lincolnenins, and huanglongmycin. Setomimycin BGC includes post-PKS assembly/cycling enzymes StmD (C-9 ketoreductase), StmE (aromatase), and StmF (thioesterase) as key components. The heterodimeric TcmI-like cyclases StmH and StmK are proposed to aid in forming the setomimycin monomer. In addition, StmI (P-450) is predicted to catalyze biaryl coupling of two monomeric setomimycin units, with StmM (ferredoxin) specific to setomimycin BGC. The roles of StmL and StmN, part of the nuclear transport factor 2 (NTF-2)-like protein family and unique to setomimycin BGCs, could particularly interest biochemists and combinatorial biologists. α-Glucosidase, a key enzyme in type 2 diabetes, hydrolyzes carbohydrates into glucose, thereby elevating blood glucose levels. This study aimed to assess the α-glucosidase inhibitory activity of EtOAc extracts of JCM 3382 and setomimycin. The JCM 3382 EtOAc extract and setomimycin exhibited greater potency than the standard inhibitor, acarbose, with IC50 values of 285.14 ± 2.04 μg/mL and 231.26 ± 0.41 μM, respectively. Kinetic analysis revealed that setomimycin acts as a competitive inhibitor of α-glucosidase. Molecular docking demonstrated two hydrogen bonds with maltase-glucoamylase chain A residues Thr205 and Lys480 (binding energy = −6.8 kcal·mol⁻¹), two π-π interactions with Trp406 and Phe450, and one π-cation interaction with Asp542. Residue-energy analysis highlighted Trp406 and Phe450 as key in setomimycin's binding to maltase-glucoamylase. These findings suggest that setomimycin is a promising candidate for further enzymological research and potential antidiabetic therapy.

It’s been reported that Mori Folium (MF) and Eucommiae Cortex (EC) exhibit pharmacological effects on the treatment of immunosuppression. However, the mechanism of MF and EC against immunosuppression remains unclear. This study aims to explore the mechanism of action of MF and EC for the treatment of immunosuppression through network pharmacology, molecular docking, molecular dynamics simulations and animal experiments. As a result, 11 critical components, 9 hub targets, and related-signaling pathways of treatment of immunosuppression were obtained based on network pharmacology. The molecular docking suggested that 11 critical components exhibited great binding affinity to 9 hub targets of immunosuppression. The molecular dynamics simulations results showed that (-)-Tabernemontanine-AR, beta-sitosterol-AR and Dehydrodieugenol-HSP90AA1 complexes are stably bound. Additionally, in the animal experiments, the treated group results compared to the control group suggest that MF and EC have a significant effect on the treatment of immunosuppression. Therefore, the MF and EC treatment for immunosuppression may take effects in a multi-component, multi-target, and multi-pathway manner. The results would provide novel insights into the treatment of immunosuppression in humans.

Organic-organic interactions play important roles in the secondary organic aerosol formation, but the interactions are complex and poorly understood. Here we use environmental molecular beam experiments combined with molecular dynamics simulations to investigate the interactions between methanol and nopinone, as atmospheric organic proxies. In the experiments, methanol monomers and clusters are sent to collide with three types of surfaces, i.e., graphite, thin nopinone coating on graphite and multilayer surfaces, at temperatures between 140 K and 230 K. Methanol monomers are efficiently scattered from the graphite surface, whereas the scattering is substantially suppressed from nopinone surfaces. The desorption from the three surfaces is similar, suggesting that all the surfaces have weak or similar influences on the methanol desorption. The molecular dynamics results show that upon collisions the methanol clusters shatter, and the shattered fragments quickly diffuse and recombine to clusters. The desorption involves a series of processes, including detaching from clusters and desorbing as monomers. The experimental results also show that all trapped methanol molecules completely desorb within a short experimental time scale at temperatures of 180 K and above. At lower temperatures, the desorption rate decreases, and a long experimental time scale is used to resolve the desorption, where three desorption components are identified. The fast component is beyond the experimental detection limit. The intermediate component exhibits multi-step desorption character and has an activation energy of E_a = 0.18 ± 0.03 eV, in good agreement with simulation results. The slow desorption component is related to diffusion processes due to the weak temperature dependence.

Kidneys of mice, rats and humans possess progenitors that maintain daily homeostasis and take part in endogenous regenerative processes following injury, owing to their capacity to proliferate and differentiate. In the glomerular and tubular compartments of the nephron, consistent studies demonstrated that well-characterized, distinct populations of progenitor cells, localized in the parietal epithelium of Bowman capsule and scattered in the proximal and distal tubules, could generate segment-specific cells in physiological conditions and following tissue injury. However, defective or abnormal regenerative responses of these progenitors can contribute to pathologic conditions. The molecular characteristics of renal progenitors have been extensively studied, revealing that numerous classical and evolutionarily conserved pathways, such as Notch or Wnt/β-catenin, play a major role in cell regulation. Others, such as retinoic acid, renin-angiotensin-aldosterone system, TLR2 (Toll-Like Receptor 2) and leptin, are also important in this process. In this review, we summarize the plethora of molecular mechanisms directing renal progenitor responses during homeostasis and following kidney injury. Finally, we will explore how single cell RNA sequencing could bring the characterization of renal progenitors to the next level, while knowing their molecular signature is gaining relevance in the clinic.

Aquaporin-3 (AQP3) is one of the aquaglyceroporins, which is expressed in the basolateral layer of the skin membrane. Studies have reported that human skin squamous cell carcinoma overexpresses AQP3 and inhibition of its function may alleviate skin tumorigenesis. In the present study, we have applied a virtual screening method that encompasses filters for physicochemical properties and molecular docking to select potential hit compounds that bind to the Aquaporin-3 protein. Based on molecular docking results, the top 20 hit compounds were analyzed for stability in the binding pocket using unconstrained molecular dynamics simulations and further evaluated for binding free energy. Furthermore, examined the ligand-unbinding pathway of the inhibitor from its bound form to explore possible routes for inhibitor approach to the ligand-binding site. With a good docking score, stability in the binding pocket, and free energy of binding, these hit compounds can be developed as Aquaporin-3 inhibitors in the near future.

For condensed systems, the incorporation of quantum chemical solvent effects into molecular dynamics simulations has been a major concern. To this end, quantum mechanical/molecular mechanical (QM/MM) techniques are popular and powerful options to treat gigantic systems. However, they cannot be directly applied because of temporal and spatial discontinuity problems. To overcome these problems, in a previous study, we proposed a corrective QM/MM method, size-consistent multipartitioning (SCMP) QM/MM, and successfully demonstrated that, using SCMP, it is possible to perform stable molecular dynamics simulations by effectively taking into account solvent quantum chemical effects. The SCMP method is characterized by two original features: size-consistency of a QM region among all QM/MM partitioning and partitioning update. However, in our previous study, the performance was not fully elicited compared to the theoretical upper bound, and the optimal partitioning update protocol and parameters were not fully verified. To elicit the potential performance, in the present study, we simplified the theoretical framework and modified the partitioning protocol.

(1) Background: Human papillomaviruses (HPV) are primarily related to cervical cancer around the world. This study presented an in silico design of a chimeric L1/L2 mRNA-based vaccine candidate against a diverse spectrum of human papillomaviruses; (2) Methods: We employed a computational strategy to predict and screen for antigenicity, allergenicity, toxicity, and conservation of immunodominant cytotoxic T cells, helper T cells, and linear B-cell epitopes from the most conserved HPV L1 and L2 proteins. The selected epitope peptides were coupled together with appropriate linkers to create a multi-epitope chimeric peptides. The mRNA-based vaccine was developed by combining an optimized open reading frame with co-translation regulation and co-transcription components; (3) Results: Twenty five (25) non-toxic, non-allergenic, highly antigenic, and conserved epitopes were identified. The final vaccine had a molecular weight of 60161.29 kDa, a theoretical pI of 9.44, and a solubility of 0.451. The mRNA-based vaccine had a stable mRNA structure with a minimum free energy of -731.10 kcal/mol and an estimated molecular weight of 644.98 kDa. The suggested vaccine had no cross-affinity for the human genome and achieved a global population coverage rate of 86.24%. The vaccine demonstrated a stable docking complex and robust interactions with Toll-like receptor 4 (TLR4) with binding affinities of -1377.5 kcal/mol and was capable of stimulating immunoglobulin IgG, IgM, T and B cells lymphocytes, and cytokines by in silico simulation; (4) Conclusion: The vaccine is considered a preventive and therapeutic candidate; however, further in vitro and in vivo tests are needed to prove its safety and efficacy.

The elderly population is particularly susceptible to myocardial damage caused by heat stress (HS). With the continuous increase in global temperatures, this poses an undeniable challenge to human health. This study aimed to investigate the therapeutic effect of Shikimic Acid (SA) in treating HS-induced myocardial damage through network pharmacology, molecular docking, molecular dynamics (MD) simulations, and in vitro experiments. Network pharmacology analysis revealed that SA alleviates the inflammatory response in the heart muscle associated with HS by targeting 60 disease-related targets, such as TP53, TNF, and IL-6. The tumor proteoglycan, the PI3K-Akt signaling pathway, and hepatitis B were identified as the primary therapeutic pathways. Molecular docking and MD simulation results confirmed that SA can stably bind to both TNF (-6.642 kcal/mol) and IL-6 (-7.261 kcal/mol), and the complexes of Shikimic Acid-IL-6 and Shikimic Acid-TNF did not undergo significant changes within 100 ns. In addition, the in vitro experiment showed that 250 μM SA significantly increased the survival rate of HL-1 cells exposed to 43°C for 2 h by 71.1% and promoted the proliferation of HL-1 cells by 30.7%. Overall, these findings provide substantial evidence for the therapeutic potential of SA in HS-induced myocardial damage.

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.

The monkeypox virus (MPXV) has become a major threat due to the increasing global caseload and the ongoing multi-country outbreak in non-endemic territories. Due to limited research in this avenue and the lack of intervention strategies, the present study was aimed to virtually screen bioactive phytochemicals against envelope proteins of MPXV via rigorous computational approaches. Molecular docking and molecular dynamic (MD) simulations were used to investigate the binding affinity of 12 phytochemicals against three envelope proteins of MPXV, viz., D13, A26, and H3. Silibinin, oleanolic acid, and ursolic acid were computationally identified as potential phytochemicals that showed strong binding affinity towards all the tested structural proteins of MPXV through molecular docking. The stability of the docked complexes was also confirmed by MD simulations. ADME analysis also computationally confirmed the drug-like properties of the phytochemicals, thereby asserting their suitability for consumption. Hence, this study envisions the candidature of bioactive phytochemicals as promising inhibitors against the MPXV, serving as template molecules that could further be experimentally evaluated for their efficacy against monkeypox.

Zika virus (ZIKV) has been characterized as one of the potential pathogens and placed under future epidemic outbreaks by the WHO. However, lack of potential therapeutics can result in uncontrolled pandemic like other human pandemic viruses; therefore, prioritized effective therapeutics development has been recommended against ZIKV. In this context, the present study adopted the strategy to explore the lead compounds from Azadirachta indica against ZIKV via concurrent inhibition of the ZIKVpro and ZIKVRdRp proteins using molecular simulations. Initially, structure-based virtual screening of 44 bioflavonoids reported in Azadirachta indica against the crystal structures of targeted ZIKV proteins resulted in the identification of top four common bioflavonoids, viz. rutin, nicotiflorin, isoquercitrin, and hyperoside. These compounds showed substantial docking energy (-7.9 to -11.01 kcal/mol) and intermolecular interactions with essential residues of ZIKVpro (His51, Asp74, and Ser135) and ZIKVRdRp (Asp540, Ile799, and Asp665) by comparison to the reference compounds, O7N inhibitor (ZIKVpro) and Sofosbuvir inhibitor (ZIKVRdRp). Long interval molecular dynamics simulation (500 ns) on the selected docked poses reveals the stability of respective docked complexes contributed by intermolecular hydrogen bonds and hydrophobic interactions. The predicted complex stability was further supported by calculated end-point binding free energy using molecular mechanics generalized born surface area (MM/GBSA) method. Consequently, the identified common bioflavonoids are recommended as promising therapeutic inhibitors of ZIKVpro and ZIKVRdRp for further experimental assessment.

An important stage in SARS-CoV-2 life cycle is the fusion of spike(S) protein with the ACE2 host-cell receptor. Therefore, to explore conserved features in S protein dynamics and to identify potentially novel regions for drugging, we measured variability derived from 791 viral genomes and studied its properties by MD simulation. The findings indicated that S2 subunit (HR1, CH, and CD domains) showed low variability, low fluctuations in MD, and displayed a trimer cavity. By contrast, the RBD domain, which is typically targeted in drug discovery programmes, exhibits more sequence variability and flexibility. Interpretations from MD suggest that the monomer is in constant motion showing transitions up-to-down state, and the trimer cavity may function as a 'bouncing spring' that may facilitates S protein interactions with ACE2. Feasibility of trimer cavity for potential drug target was examined by SBVS screening. Several hits that have already been validated or suggested to inhibit the SARS-CoV-2 virus in cell systems were identified; in particular, the data suggest an action mechanism for such molecules including Chitosan and macrolide types. These findings identify a novel binding-site formed by the S protein, that might assist in future drug discovery programmes aimed at targeting the CoV family of viruses.

The global pandemic caused by the new SARS-COV-2 coronavirus makes it necessary to search for drugs for its control. Within of this research it has been known that the ivermectin drug, a FDA-approved drugs which is formulated as an 80:20 mixture of ivermectin B1a and B1b and used commonly for parasitic infections, has an inhibitory effect on viruses, includes SARS-COV-2 at in vitro level. In the particular case of SARS-COV-2 its mechanism of action remains elusive and controversial. Interestingly, the energy of interaction of ivermectin with any of the proteins the SARS-CoV-2 and the possible structural alterations at the protein level that this drug can cause have not been reported. In this sense, we carried out a bioinformatics study with docking strategies and molecular dynamics to predict the binding and disturbance induced by ivermectin in proteins associated with SARS-CoV-2. We use DockThor and Molegro docking scores. The DockThor server and myPresto software were used to build complexes and dynamics studies, respectively. The results obtained suggested that ivermectin is capable of docking with the 3CL protease and the HR2 domain, and may promote structural changes in these proteins by inducing unfolding/folding. Specifically, ivermectin brings protease to a significantly more deployed conformational state and the HR2 domain to a more compact state compared to the native state. Finally, it is shown that B1a and B1b macrocyclic lactones have a behavior different from to each target protein. These results suggest a possible inhibitory effect against SARS-CoV-2 due to a synergistic role of this drug to spontaneously bind to two important proteins involve in the proliferation of this virus. However, more studies are required on this possible mechanism of action.

The SARS-CoV-2 was confirmed to cause the regional outbreak of coronavirus disease 2019 (COVID-19) in Wuhan, China. The 3C-like protease (3CL^pro), an essential enzyme for viral replication, is a valid target to compacts SARS-CoV and MERS-CoV. In this research, an integrated library consisting of 1000 compounds from Asinex Focused Covalent (AFCL) library and 16 FDA-approved protease inhibitors were screened against SARS-CoV-2 3CL^pro. Top compounds with significant docking scores and making stable interactions with catalytic dyad residues were obtained. The screening results in identification of compound 621 from AFCL library as well as Paritaprevir and Simeprevir from FDA-approved protease inhibitors as potential inhibitors of SARS-CoV-2 3CL^pro. The mechanism and dynamic stability of binding between the identified compounds and SARS-CoV-2 3CL^pro were characterized using 50 nanoseconds (ns) molecular dynamic (MD) simulation approach. The identified compounds are potential inhibitors worthy of further development as SARS-CoV-2 3CLpro inhibitors/drugs. Importantly, the identified FDA-approved therapeutics could be ready for clinical trials to treat infected patients and help to curb the COVID-19.

The antioxidant activity of molecules constitutes an important factor for the regulation of redox homeostasis and reduction of oxidative stress. Cells affected by oxidative stress can undergo genetic alteration, causing structural changes and promoting the onset of chronic diseases, such as cancer. The in silico study performed here was developed to evaluate the antioxidant potential of two molecules, ZINC08706191 (Z91) and ZINC08992920 (Z20), with recognized epithelial anticancer potential. Molecular docking, chemical-quantum calculations and Pearson's correlation were performed. The Z91 and Z20 molecules showed lower binding free energy (ΔG) values for the receptor-ligand interaction than the reference molecules (caffeine – CAF and ascorbic acid – AA), and better results for values of molecular descriptors correlated with ΔG, resulting in a decrease of ΔG. Strong correlations were observed between ΔG values for the five receptors evaluated and ΔG values of the potential epithelial anticancer activity available in literature. These results attest to the significant antioxidant potential of the Z91 and Z20 molecules and their strong relation with the potential epithelial anticancer activity and may be indicated for further analysis in relation to the control of oxidative stress and epithelial anticancer activity.

Borneol and menthol were two terpenes wildly used as penetrate enhancer in transdermal drug delivery. To explore their penetration enhancement effect towards hydrophilic drug, 5-FU was selected as model drug. A method combined vitro permeation studies and coarse grain molecular dynamics was used to investigate their penetration enhancement effect towards 5-FU. As a result, although both borneol and menthol showed a penetration enhancement effect towards 5-FU, they differed a lot in the penetration enhancement mechanism, which was also thought to account for their different penetration enhancement effect. As for menthol, SC bilayer disrupting effect seemed to be its main mechanism. While for borneol, its mechanism seemed to be more complicated. Except for disrupting the SC bilayer, it could also increase the permeation of 5-FU by enhancing the diffusion rate of 5-FU or inducing the formation of transient pore. All of this enable us a molecular understanding of borneol and menthol’s penetration enhancement effect towards hydrophilic drug, which might provide some guidance in the latter research and application.

Any homeostatic protometabolism would have required orchestration of disparate biochemical pathways into integrated circuits. Extraordinarily specific moieties were also required. Assembly Theory conflated with its cousins—Complexity Theory, Chaos theory, Quantum Mechanics, Irreversible Nonequilibrium Thermodynamics and Molecular Evolution theory— collectively have great naturalistic appeal in hopes of their providing the needed exquisite steering and controls. They collectively offer the best hope of circumventing the need for active selection required to formally orchestrate bona fide formal organization (as opposed to the mere self-ordering of chaos theory) [1]. This paper focuses specifically on Assembly Theory’s contribution to life-origin models.

Here, we present a general framework for computing the infrared anharmonic vibrational spectra of polyatomic molecules using the Born-Oppenheimer molecular dynamics (BOMD) with the software PyRAMD. To account for the nuclear quantum effects, we suggest using the simplified Wigner sampling (SWS) approach [\textit{Phys. Chem. Chem. Phys.}, 2023, \textbf{25}, 18406-18423] simultaneously coupled with the Andersen and Berendsen thermostats. We propose a new criterion for selecting the parameter of the SWS based on the molecules' harmonic vibrational frequencies and usage of the large-time-step blue shift correction, allowing for a decrease in computational expenses. For the Fourier transform of the dipole moment autocorrelation function, we propose to use the regularized least-squares analysis, which allows us to obtain higher frequency resolution than the direct application of fast Fourier transform. Finally, we suggest the usage of the pre-parameterized scaling factors for the IR spectra from BOMD, also providing the scaling factors for the spectra at the BLYP-D3(BJ)/6-31G, PBE-D3(BJ)/6-31G, and PBEh-3c levels of theory.

Benzazoles, such as benzoxazoles and benzothiazoles, are compounds with important biological and pharmacological activities and important intermediaries in synthesis. This report presents the synthesis of a butanamide derived from linking 5-chloro-2-aminobenzoxazole and 2-aminobenzothiazole via 4-chlorobutanoyl chloride. The corresponding compound N-(benzothiazol-2-yl)-4-((5-chlorobenzoxazol-2-yl)aminobutanamide was obtained 76 % yield global using accessible starting materials and methodology in two reaction steps. Furthermore, we made docking studies of this compound on 3-TOP protein to explore its potential as an antidiabetic agent.

Interferons were the original prototype cytokine system discovered in 20th-century research. As the name implies, they were originally thought to be synthesised and secreted between cells. Thanks to technological advances, the processes involved in protein secretion can be explained comparatively more clearly at both the genetic and biochemical levels. The discovery of interferon (IFN) occurred when genetic research was still in its infancy. Franklin and Wilkins discovered the structure and function of deoxyribonucleic acid (DNA) at the same time as Crick and Watson; however, Isaacs and Lindemann, two scientists, described the first IFN in 1957. Mutations can be caused by inherent genetic protein synthesis and during infection as well as within IFN regulation pathways affecting cell proliferation. This remains central to host cell IFN synthesis and effects through IFN protein receptor subunits defined by 6 protein domains. Type II IFN is key to immune cell function secreted by a variety of immune cells, mainly natural killer (NK) as well as T cells. Single–stranded and/or double–stranded RNA/DNA viruses, as well as bacterial infections (e.g., Escherichia coli) and fungal infections (e.g., Aspergillus), also affect IFN regulation. Pathogenic proteins utilise intra/extracellular proteins that sense foreign antigens like Toll–like Receptors (TLRs), affected by mutations within the human cellular IFN transduction pathways. Since the discovery of the third IFN type in 2003, when immune cell phenotypes were further characterised, questions remain about the immunological mechanisms contributing to the regulation of the innate and adaptive host immune system. Alterations in the synthesis of type I/II/III host IFNs can differentially and beneficially alter homeostatic cellular pathways in pathological disease, with type I IFN being synthesised in cancer as well as by homeostatic cells. Therefore, considered here are the overall IFN molecular, cell regulatory mechanisms in the context of immune cell research developments.

A self-consistent charge density-functional tight-binding method combined with molecular dynamics simulations is employed to reveal the effect of polymorph on the thermal decomposition stability of 1,1-Diamino-2,2-Dinitroethylene (FOX-7). Two types of heating, constant temperature heating and temperature-programmed heating, are adopted. Potential evolution indicates that γ-FOX-7 possesses the lowest thermal stability, as it is closer to the decomposition state. Crystal form has an important influence on the thermal decomposition of FOX-7, resulting in different decomposition rates and initial reactions. In general, β- and γ-FOX-7 always decompose more completely than α-FOX-7. This work emphases the importance of polymorph dependent initial decay of an energetic polymorphic compound once heated in a volume constrained condition.

CoFe2O4 is a promissing catalytic material for many chemical reactions. We have used ab initio molecular dynamic simulations to study the structure and reactivity of the A- and B-terminations of the low index CoFe2O4 (001) surfaces to water adsorption at room temperature. Upon adsorption, water partly dissociates on both termination with a higher dissociation degree on the A- termination (30% versus 19%). The 2-fold coordinated Fe3+(tet) in the tetrahedral voids and the 5-fold coordinated Fe3+(oct) in the octahedral voids are the main active sites for water dissociation on the A- and B-termination, respectively. Molecular water, hydroxydes and surface OH resulting from proton transfer to surface oxygens are present on the surfaces. Both water free sur- face terminations undergo reconstruction. The outermost Fe3+(tet) on the A-termination and B-termination move towards the nearby unoccupied octahedral voids. In the presence of a thin film of 32 water molecules, the reconstructions are partially and completely lifted on the A- and B-termination, respectively.

We research the interaction between six representative carbon-based nanoparticles (CBNs) and 20 standard amino acids through all-atom molecular dynamics simulations. The six carbon-based nanoparticles are fullerene(C60), CNT55L3, CNT1010L3, CNT1515L3, CNT2020L3, and two-dimensional-graphene(Graphene33). Their curvatures decrease sequentially, and all of CNT are single-walled carbon nanotubes. We have observed that as the curvature of CBNs decreases, the adsorption effect of 20 amino acids with them has an increasing trend. In addition, we also used multi-dimensional clustering to analyze the adsorption effects of 20 amino acids on six carbon-based nanoparticles. We observed that the π-π interaction still plays an extremely important role in the adsorption of amino acids on carbon-based nanoparticles. Individual long-chain amino acids and “Benzene-like” Pro also have a strong adsorption effect with carbon-based nanoparticles.

Drug discovery has been initially attributed to coincidence or trial and error where the traditional approach was complex, lengthy, and expensive. Conventional drug discovery methods require the costly random screening of synthesized compounds or natural products. Another downside for this approach is the wide dependency on the experimental use of animals for in vi-vo testing. Currently, in silico modeling has become a vital tool for drug discovery and repurposing, and molecular docking is being used to find the best matching between a ligand and a molecule. Practical application of in silico docking will predict the biomolecular interactions between the drug and the target host. Beauvericin (BEA) is an emerging mycotoxin produced by the entomopathogenic fungus Beauveria bassiana. Originally investigated for its pesticide capability, BEA is now considered as a molecule of interest for its potentially diverse biotechnological applications in the pharmacological industry and the field of medicine. In this manuscript, we will provide an overview of the repurposing of BEA into a potentially superior therapeutic molecule in a broad range of diseases. Furthermore, considerable attention has been given to the fundamental role of in silico techniques to i) further investigate the spectrum of this secondary metabolite and ii) elucidate the pathways of BEA for its promising therapeutic action

Bioluminescence – i.e., the emission of visible light by living organisms - is defined as a biochemical reaction involving, at least, a luciferin substrate, an oxygen derivative, and a specialised luciferase enzyme. In some cases, the enzyme and the substrate are durably associated and form a photoprotein. While this terminology is educatively useful to explain bioluminescence, it gives a false idea that all luminous organisms are using identical or homologous molecular tools to achieve light emission. As usually observed in biology, reality is more complex. To date, 11 different luciferins have indeed been discovered, and several non-homologous luciferases lato sensu have been identified which, all together, confirms that bioluminescence emerged independently multiple times during the evolution of living organisms. While some phylogenetically related organisms may use non-homologous luciferases (e.g., at least four convergent luciferases are found in Pancrustacea), it has also been observed that phylogenetically distant organisms may use homologous luciferases (e.g., parallel evolution observed in some cnidarians, tunicates and echinoderms that are sharing a homologous luciferase-based system). The evolution of luciferases then appears puzzling. The present review takes stock of the diversity of known “bioluminescent proteins”, their evolution and potential evolutionary origins. A total of 134 luciferase and photoprotein sequences have been investigated (from 75 species and 11 phyla), and our analyses identified 12 distinct types – defined as a group of homologous bioluminescent proteins. The literature review indicated that genes coding for luciferases and photoproteins have potentially emerged as new genes or have been co-opted from ancestral non-luciferase/photoprotein genes. In this latter case, the homologous gene’s co-options may occur independently in phylogenetically distant organisms.

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.

An important parameter in the development of a new drug is the drug's affinity to the identified target (protein/enzyme). Predicting the ligand binding to the target (protein/enzyme) by molecular simulation would allow the synthesis to be restricted to the most promising compounds.A restricted hybrid HF-DFT calculation was performed in order to obtain the most stable conformer of each ligand and a series of DFT calculations using the B3LYP levels with 6-31G* basis set has been conducted. The docking studies of the quinolone compounds will be performed with the CLC Drug Discovery Workbench to identify and visualize the ligand-receptor interaction mode.

Introduction of exotic species in new areas through anthropic action is one of the major problems that can affect biodiversity. Branchiomma luctuosum is known for its highly invasive potential and the actual occurrence of species commonly associated with port activity areas is an extra evidence that this anthropogenic activity should not be underestimated. In order to develop suitable molecular markers for future studies on colonization routes and population dynamics of the invading individuals of B. luctuosum, nine highly polymorphic microsatellite loci were isolated and their polymorphism levels were evaluated. These loci showed a range of number of alleles per locus from five to ten and all loci had a high level of genetic diversity, and exhibited significant heterozygote deficiencies probably due to the presence of null alleles. Significant deviations from the Hardy-Weinberg equilibrium were detected at several loci and most of them were related to a heterozygous deficit. Heterozygous deficiency can be expected in this case due to the biology and history of this invasive species, in relation to its recent introduction in Brazilian coast and possible action of multiple introductory events.

Stabilization of a G-quadruplex (G4) in the promotor of the c-MYC proto-oncogene leads to inhibition of gene expression, and it thus represents a potentially attractive new strategy for cancer treatment. However, most G4 stabilizers show little selectivity among the many G4s present in the cellular complement of DNA and RNA. Intriguingly, a crescent-shaped cell-penetrating thiazole peptide, TH3, preferentially stabilizes the c-MYC G4 over other promotor G4s, but the mechanisms leading to this selective binding remain obscure. To investigate these mechanisms at the atomic level, we performed an in silico comparative investigation of the binding of TH3 and its analogue TH1 to the G4s from the promotors of c-MYC, c-KIT1, c-KIT2, and BCL2. Molecular docking and molecular dynamics simulations, combined with in-depth analyses of non-covalent interactions and bulk and per-nucleotide binding free energies, revealed that both TH3 and TH1 can induce the formation of a sandwich-like framework through stacking with both the top and bottom G-tetrads of the c-MYC G4 and the adjacent terminal capping nucleotides. This framework produces enhanced binding affinities for the c-MYC G4 relative to other promotor G4s, with TH3 exhibiting an outstanding binding priority. Van der Waals interactions were identified to be the key factor in complex formation in all cases. Collectively, our findings fully agree with available experimental data. Therefore, the identified mechanisms leading to specific binding of TH3 towards the c-MYC G4 provide valuable information to guide the development of new selective G4 stabilizers.

Tuberculosis (TB) continues to pose a global health challenge, exacerbated by the rise of drug-resistant strains. The development of new TB therapies is an arduous and time-consuming process. To expedite the discovery of effective treatments, computational structure-based drug re-purposing has emerged as a promising strategy. From this perspective, conditionally essential targets present a valuable opportunity, and the mycobactin biosynthesis pathway stands out as a prime example highlighting the intricate response of Mycobacterium tuberculosis (Mtb) to changes in iron availability. This study focuses on the re-purposing and revival of FDA-approved drugs (library) as potential inhibitors of MbtA, a crucial enzyme in mycobactin biosynthesis in Mtb conserved among all species of mycobacteria. Literature suggests this pathway to be associated with drug efflux pumps, which potentially contribute to drug resistance. This makes it a potential target for antitubercular drug discovery. Herein we utilized cheminformatics and structure-based drug repurposing approaches, viz., molecular docking, dynamics, and PCA analysis, to decode the intermolecular interactions and binding affinity of the FDA-reported molecules against MbtA. The virtual screening revealed ten molecules with significant binding affinities and interactions with MbtA. These drugs, originally designed for different therapeutic indications (4: antiviral, 3: anticancer, 1: CYP450 inhibitor, 1: ACE inhibitor, and 1: leukotriene antagonist), are repurposed as potential MbtA inhibitors. Furthermore, our study explores the binding modes and interactions between these drugs and MbtA, shedding light on the structural basis of their inhibitory potential. Principal component analysis highlighted significant motions in MbtA-bound ligands, emphasizing the stability of the top Protein-Ligand Complexes (PLCs). This computational approach provides a swift and cost-effective method for identifying new MbtA inhibitors, which can subsequently undergo validation through experimental assays. This streamlined process is facilitated by the fact that these compounds are already FDA-approved and have established safety and efficacy profiles. This study has the potential to lay the groundwork for addressing the urgent global health challenge at hand, specifically in the context of combating Antimicrobial Resistance (AMR) and Tuberculosis (TB).

Polo-like kinase 1 (PLK1) plays a pivotal role in cell division regulation and emerges as a promising therapeutic target for cancer treatment. Consequently, the development of small-molecule inhibitors targeting PLK1 has become a focal point in contemporary research. The adenosine triphosphate (ATP)-binding site and the polo-box domain in PLK1 present crucial interaction sites for these inhibitors, aiming to disrupt the protein's function. However, designing potent and selective small-molecule inhibitors can be challenging, requiring a deep understanding of protein-ligand interaction mechanisms at these binding sites. In this context, our study leverages the fragment molecular orbital (FMO) method to explore these site-specific interactions in-depth. Through the FMO approach, we used the FMO method to elucidate the molecular mechanisms of small-molecule drug binding to these sites to design PLK1 inhibitors that are both potent and selective. Our investigation further entailed a comparative analysis of various PLK1 inhibitors, each characterized by distinct structural attributes. This comparison enhanced our understanding of the structure-activity relationships of these inhibitors, underscoring the efficacy of the FMO method in identifying critical binding features and predicting binding modes for small-molecule ligands. Furthermore, Our research spotlighted "hot spot" residues instrumental for selective and robust PLK1 binding. Our findings offer profound insights, priming the rational design of innovative potential PLK1 inhibitors with significant implications for developing anticancer therapeutics.

In this study, we investigate the structural changes, electronic properties, and charge redistribution within azo-bithiophene (Azo-BT) chemisorbed monolayers under different light stimuli using density functional theory and molecular dynamics simulations. We consider two types of switches, Azo-BT and BT-Azo, with different arrangements of the Azo and BT blocks counting from the anchor thiol group. The chemisorbed monolayers of pure cis- and trans-isomers with a surface concentration of approximately 2.7 molecules per nm2 are modeled on a gold surface using classical all-atom molecular dynamics. Our results reveal significant shrinkage of the BT-Azo layer under UV illumination, whereas the thicknesses of the Azo-BT layer remain comparable for both isomers. This difference in behavior is attributed to the ordering of the trans-molecules in the layers, which is more pronounced for Azo-BT, leading to a narrow distribution of the inclination angle to the gold surface. Conversely, both layers consisting of cis-switches exhibit disorder, resulting in similar brush heights. To study charge transfer within the immobilized layers, we analyze each snapshot of the layer and calculate mean charge transfer integrals using Nelsen’s algorithm for a number of interacting neighboring molecules. Combining these integrals with reorganization energies defined for the isolated molecules, we evaluate the charge transfer rates and mobilities for electron and hole hopping within the layers at room temperature based on Marcus’ theory. This research offers new perspectives for the innovative design of electrode surface modifications and provides insights into controlling charge transfer within immobilized layers using light triggers. Additionally, we identify molecular properties that are enhanced through specific molecular design, which contributes to the development of more efficient molecular switches for various electronic applications.

Infection of hosts by morbilliviruses is facilitated by the interaction between viral hemagglutinin (H-protein) and the signaling lymphocytic activation molecule (SLAM). Recently, the functional importance of the N-terminal region of human SLAM as a measles virus receptor was demonstrated. However, the functional roles of this region in the infection process by other morbilliviruses and host range determination remain unknown partly because this region is highly flexible, which has hampered accurate structure determination of this region by X-ray crystallography. In this study, we analyzed the interaction between the H-protein from canine distemper virus (CDV-H) and SLAMs by a computational chemistry approach. Molecular dynamics simulations and fragment molecular orbital analysis demonstrated that the unique His28 in the N-terminal region of SLAM from Macaca is a key determinant that enables formation of a stable interaction with CDV-H, providing a basis for CDV infection in Macaca. The computational chemistry approach presented should enable determination of molecular interactions involving regions of proteins that are difficult to predict from crystal structures because of their high flexibility.

Three-dimensional structure of six closely related hydrogenases from purple bacteria has been modeled by combining template-based and ab initio modeling approach. The results lead to conclusion that there should be 4Fe3S-cluster in the structure of these enzymes. Thus, these hydrogenases could drive attention for exploring their oxygen tolerance and practical applicability in hydrogen fuel cells. Analysis of 4Fe3S-cluster's microenvironment showed intragroup heterogeneity. Possible function of the C-terminal part of the small subunit in membrane binding has been discussed. Comparison of the built models with existing hydrogenases of the same subgroup (membrane-bound oxygen-tolerant hydrogenases) has been carried out. Analysis of intramolecular interactions in the large subunits showed statistically reliable differences in number of hydrophobic interactions and number of ionic interactions. Molecular tunnels were mapped in the models and compared with structures from PDB. Protein-protein docking showed that these enzymes could exchange electrons in oligomeric state, which is important for oxygen-tolerant hydrogenases. Molecular docking with model electrode compounds showed mostly the same results as with hydrogenases from E.coli, H. marinus, R. eutropha, S. enterica; some interesting results were shown in case of HupSL from Rba. sphaeroides and Rvx. gelatinosus.

Prototype of a family of at least nine members, c-src tyrosine kinase is a therapeutically interesting target, because its inhibition might be of interest not only in a number of malignancies, but also in a diverse array of conditions, from neurodegenerative pathologies to certain viral infections. Computational methods in drug discovery are considerably cheaper than conventional methods and offer opportunities of screening very large numbers of compounds in conditions that would be simply impossible within the wet lab experimental settings. We have explored the use of global QSAR models and molecular ligand docking in the discovery of new c-src tyrosine kinase inhibitors. Using a data set of 1038 compounds from ChEMBL and 19 blocks of molecular descriptors, we have developed over 200 QSAR classification models, based on six machine learning algorithms and 17 feature selection methods. We have selected 49 with reasonably good performance (positive predictive value and balanced accuracy higher than 70% in nested cross validation) and the models were assembled by stacking with a simple majority vote and used for the virtual screening of over the “named” ZINC data set (over 100,000 compounds). 744 compounds were predicted by at least 50% of the QSAR models as active, 147 compounds were within the applicability domain and predicted by at least 75% of the models to be active. The latter 147 compounds were submitted to molecular ligand docking using Vina and Ledock, and a number of 90 were predicted to be active based on the binding energy. External data from CHEMBL and PUBCHEM confirmed that at least 7.83% (in the case of QSAR) or 6.67% (in the case of integrated QSAR and molecular docking) of the compounds are active on the c-src target.

Precise classification of sarcomas is crucial to drive optimal clinical management. In this prospective, multicenter, observational study within the Hellenic Group of Sarcoma and Rare Cancers (HGSRC) we assessed the effect of expert pathology review, coupled with the application of molecular diagnostics, on diagnosis and management of sarcoma patients. Newly diagnosed sarcoma patients were addressed by their physicians to one of the two sarcoma pathologists of HGSRC for histopathological diagnostic assessment. RNA next-generation sequencing was per-formed to all samples, using a platform targeting 86 sarcoma gene fusions. Additional molecular methods were performed on the opinion of the expert pathologist. Therefore, the expert pathologist provided a final diagnosis based on the histopathological findings and, when necessary on molecular tests. In total, 128 specimens from 122 patients were assessed. Among the 119 cases in which there was a preliminary diagnosis by a non-sarcoma pathologist, there were 37 modifications in diagnosis (31.1%) by the sarcoma pathologist, resulting in 17 (14.2%) modifications in management. Among the 110 cases in which molecular tests were performed, there were 29 modifications in diagnosis (26.4%) through the genomic results, resulting in 12 (10.9%) modifications in management. Our study confirms that expert pathology review is of utmost importance for the optimal sarcoma diagnosis and management and should be assisted by molecular methods in selected cases.

Exploring therapeutic options is crucial in the ongoing COVID-19 pandemic caused by SARS-CoV-2. Nirmatrelvir, a potent inhibitor targeting the SARS-CoV-2 Mpro, shows promise as an antiviral treatment. Additionally, Ivermectin, a broad-spectrum antiparasitic drug, has demonstrated effectiveness against the virus in laboratory settings. However, its clinical implications are still debated. Using computational methods such as molecular docking and 100 ns molecular dynamics simulations, we investigated how Nirmatrelvir and Ivermectin interact with SARS-CoV-2 Mpro(A). Calculations using density functional theory have been instrumental in elucidating the behavior of isolated molecules, primarily by analyzing the frontier molecular orbitals. Our analysis revealed distinct binding patterns: Nirmatrelvir formed strong interactions with amino acids like MET49, MET165, HIS41, HIS163, HIS164, PHE140, CYS145, GLU166, and ASN142, showing stable binding with a root mean square deviation (RMSD) of around 2.0 Å. On the other hand, Ivermectin interacted with THR237, THR239, LEU271, LEU272, and LEU287, displaying an RMSD of 1.87 Å, indicating enduring interactions. Both ligands stabilized Mpro(A), with Ivermectin showing stability and persistent interactions despite forming fewer hydrogen bonds. These findings offer detailed insights into how Nirmatrelvir and Ivermectin bind to the SARS-CoV-2 main protease, providing valuable information for potential therapeutic strategies against COVID-19.

The novel coronavirus disease (COVID-19) pandemic has resulted in over 720 million confirmed cases and 7 million deaths worldwide, with insufficient treatment options. Innumerable efforts are being made around the world for faster identification of therapeutic agents to treat the deadly disease. Postacute sequelae of SARS-CoV-2 infection or COVID-19 (PASC), also called Long COVID, is still being understood and lacks treatment options as well. A growing list of drugs are being suggested by various in silico, in vitro and ex vivo models, however currently only two treatment options are widely used: the RNA-dependent RNA polymerase (RdRp) inhibitor remdesivir, and the main protease inhibitor nirmatrelvir in combination with ritonavir. Computational drug development tools and in silico studies involving molecular docking, molecular dynamics, entropy calculations and pharmacokinetics can be useful to identify new targets to treat COVID-19 and PASC, as shown in this work and our recent paper that identified alendronate as a promising candidate. We have now investigated all bisphosphonates which can bind competitively to nidovirus RdRp-associated nucleotidyl (NiRAN) transferase domain, and systematically down selected seven candidates (CHEMBL608526, CHEMBL196676, CHEMBL164344, CHEMBL4291724, CHEMBL4569308, CHEMBL387132, CHEMBL98211), two of which closely resemble the approved drugs minodronate and zoledronate. This work and our recent paper together provide an in silico mechanistic explanation for alendronate and zoledronate users having dramatically reduced odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations, and indicate that similar observational studies with minodronate could be valuable.

Dysregulation of the arsenal of cells, growth factors, cytokines, proteases, and proteins involved in wound healing (WH), could delay the healing process. Plants and plant-derived molecules (PDMs) have been used as medicines for skin wounds. This study aimed to select and evaluate the pro-healing capacity of PMDs. Thirty-three PDMs and 10 therapeutic targets (TTs) were chosen for molecular docking and PDMs: aristolochic acid (AA), α-copaene, selinenes, β-caryophyllene (BC), and BCoxide, and TTs: metalloproteinase (MMP) 3, MMP13, and tumor necrosis factor (TNF α) for molecular dynamics. The bests inhibitor was AA, but because of its toxicity, BC and BCoxide were chosen for evaluation in HaCat cells and BALB/c mice excisional model. Both compounds were not toxic, and higher percentages of lesion area closure (LAC) were induced by BC treatment (p<0.05). All mice healed; however, an initial delay of LAC was induced by 0.05, 0.1% BCoxide, 0.1% BC, and allantoin. At the end of treatment, lower numbers of mast cells, level of cellular infiltrates (except 1% BCoxide and 0.05, 0.5%BC), epidermal thickness tendency (except 0.1, 0.5% BCoxide), and higher values of SC-thickness, and collagen fibers only by 0.5% BC treatment compared with untreated control (p<0.05) were observed. Our results show that BC is a promising multitarget candidate for wound healing treatment and optimization campaigns.

Most existing cancer therapies negatively affect normal tissue as well as cancerous tissue. A potentially effective strategy for treating cancer that precludes off-target damage and could be an option for most patients would involve targeting one or more mutations that are ubiquitous in the given patient’s tumor(s). To effect this strategy, one would employ multi-region sequencing of a patient’s primary tumor and metastases to seek out mutations that are shared between all or at least most regions. Once the target or targets are known, one would ideally rapidly generate a molecular switch for at least one of said ubiquitous mutations that can distinguish the mutated DNA, RNA, or protein from the wild-type version and subsequently trigger a therapeutic response. I propose that the therapeutic response involve the replication of an oncolytic virus or intracellular bacterium, as any mutation can theoretically be detected by a vector that enters the cell - and automatic propagation could be very helpful. Moreover, the mutation “signal” can be easily enhanced through transcriptional and translational (if the target is an intracellular protein) enhancement. Importantly, RNA may make the best target for the molecular switches in terms of amplification of the signal and ease of targeting.

Thrombosis has threatened human health in past decades. Bacillus nattokinase is a potential thrombolytic drug without side-effect and low cost and has been introduced into the consumer market as a functional food or dietary supplement. This review firstly summarizes the biodiversity sources and fermentation process of nattokinase, and systematically expounds the structure, cata-lytic mechanism and enzymatic properties of nattokinase; In view of the problems of low fermen-tation yield, insufficient activity and stability of nattokinase, this review discusses the heterologous expression of nattokinase in different microbial hosts, and summarizes the protein and genetic en-gineering progress of nattokinase-producing strains; Finally, this review summarizes the clinical application of nattokinase.

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.

This research aims to find out whether the synthetic 1, 2, 4-triazine and its derivatives have antifungal effects and can protect humans from infection with Candida albicans. Molecular docking and molecular dynamic simulation are widely used in modern drug design to target a We are interested in using molecular docking and molecular dynamics modelling to investigate the interaction between the derivatives of 1, 2, 4-triazine and the resulting lanosterol 14 - demethylase (CYP51) of Candida albicans The inhibition of Candida albicans CYP51 is the main goal of our research. The 1, 2, 4-triazine and its derivatives have been docked to the CYP51 enzyme, which is involved in Candida albicans Multidrug Drug Resistance (MDR). Autodock tools were used to identifying the binding affinities of molecules against the target proteins. Compared to conventional fluconazole, the molecular docking results indicated that each drug has a high binding affinity for CYP51 proteins and forms unbound interactions and hydrogen bonds with their active residues and surrounding allosteric residues. The docking contacts were made using a 10 ns MD simulation with nine molecules. RMSD, RMSF, hydrogen bonds, and the Rg all confirm these conclusions. In addition, these compounds were expected to have a favorable pharmacological profile and low toxicity. The compounds are being offered as scaffolds for the development of new antifungal drugs and as candidates for future in vitro testing.

In this study, 2 chalcone analogues were synthesized and assessed through in silico and experimental methods for their potential to inhibit the lipoxygenase enzyme, which plays a role in the inflammation pathway. The structure elucidation of each chalcone was conducted through a combination of Nuclear Magnetic Resonance (NMR) and Density Functional Theory (DFT). Α “LOX-chalcone” complex, predicted by docking studies, was further examined using Molecular Dynamics (MD) simulations to evaluate the stability of the complex. After fully characterizing the “LOX-chalcone” complexes in silico, the atomic details of each chalcone’s interaction with LOX-1 and 5-LOX were revealed through Saturation Transfer Difference (STD) NMR (Nuclear Magnetic Resonance). Finally, their selectivity profile was investigated for human 15-LOX-1 and general Lipoxidase activity. The in silico methods suggest that chalcones could be promising lead compounds for drug design targeting the LOX enzyme.

Microorganisms play a pivotal role in ecology, medicine, agriculture, and biotechnology. To comprehend their diversity, physiology, and functions, isolation and characterization are imperative. This review provides a concise overview of prevalent molecular biotechniques for microbial isolation. It describes the advantages, constraints, and practical applications of both culture-dependent and culture-independent approaches. Discussed methods encompass selective/differential media, 16S rRNA gene sequencing, metagenomics, fluorescence in situ hybridization (FISH), polymerase chain reaction (PCR)-based techniques, and single-cell genomics. Emerging high-throughput sequencing and single-cell isolation methods exemplify cutting-edge procedures explored herein. Integration of these molecular tools offers comprehensive insights into microbial-diversity, community structure, metabolic capabilities, and functional interplays.

Quinoa (Chenopodium quinoa, Amaranthaceae) is a pseudocereal native to the Andes of South America that contains high protein content and adequate nutrient levels. Peru possesses an abundant morphological diversity of quinoas and is among the top producers and exporters worldwide of this precious crop. However, knowledge about the genetic and population components of quinoa from the Peruvian Andes is still limited. We here used 13 microsatellite markers to determine the genetic diversity and population structure of 105 landraces of quinoa cultivated in 11 provinces of Ayacucho, the southern Peruvian Andes. A total of 285 bands were manually scored, generating a 105 x 285 presence/absence data set. Principal coordinate analysis, and similar to a dendrogram using UPGMA clustering algorithm, showed that quinoa from Ayacucho is grouped into three clusters without a clear geographic component. Estimation of the genetic diversity indices was conducted considering the three populations (C1: south 1, C2: south 2, C3: north) determined by STRUCTURE analysis, showing mean expected heterozygosity was 0.08. The highest population divergence (FST) was exhibited for C2 and C3 (0.03), whereas the lowest for C1 and C3 (0.02). Analysis of molecular variance revealed the greatest variation within populations (80.07%) and indicated that variability between populations is 19.93%. Microsatellite markers were effective; however, more studies of the genetic components of quinoa from other Peruvian Andean localities are still needed. We expect this work helps pave the way towards the development of modern breeding programs of quinoa in Peru with accurate strategies of conservation of this nutritious crop.

Cartilage microfracturation is a surgical technique specifically designed to address chondral defects, which are injuries to the cartilage that covers the ends of bones in joints. These defects can result from traumatic injuries, degenerative conditions such as osteoarthritis, or congenital abnormalities. The primary objective of microfracture surgery is to promote the regeneration of functional cartilage tissue, thereby restoring joint function, alleviating pain, and enhancing mobility. The procedure involves creating small, controlled perforations, or microfractures, in the subchondral bone plate beneath the damaged cartilage. This process, performed with precision to minimize damage to surrounding healthy tissue, penetrates the subchondral bone to reach the bone marrow, which is rich in mesenchymal stem cells (MSCs). MSCs are multipotent cells capable of differentiating into various cell types, including chondrocytes, which are essential for cartilage production. The microfractures provide an access pathway for MSCs, allowing them to migrate from the bone marrow to the defect site. Once at the site, these stem cells differentiate into chondrocytes and initiate the formation of new cartilage tissue. This newly formed tissue helps fill the defect, aiming to restore the smooth surface of the cartilage and improve the joint's structural integrity. Understanding the molecular processes involved in cartilage repair is crucial for optimizing rehabilitation strategies and improving clinical outcomes. These processes include the regulation of various signaling pathways that control stem cell migration, differentiation, and extracellular matrix (ECM) synthesis. Growth factors such as transforming growth factor-beta (TGF-β), bone morphogenetic proteins (BMPs), and fibroblast growth factors (FGFs) play pivotal roles in these pathways, guiding MSCs in their transformation into functional cartilage cells. Additionally, the role of the ECM, composed primarily of collagen type II and proteoglycans like aggrecan, is paramount in providing structural support and biochemical signals. The ECM supports the newly differentiated chondrocytes and helps maintain the cartilage's mechanical properties. Matrix metalloproteinases (MMPs) and their inhibitors, tissue inhibitors of metalloproteinases (TIMPs), are crucial in ECM remodeling, ensuring the balance between synthesis and degradation of the matrix to support tissue repair and integrity. Angiogenesis and vascularization are also critical aspects of the molecular repair process. The initial phase of cartilage repair often includes transient vascularization facilitated by vascular endothelial growth factor (VEGF), ensuring that the newly formed tissue receives adequate nutrients and oxygen. However, as the repair tissue matures, it transitions back to an avascular state characteristic of healthy articular cartilage, which is essential for maintaining its mechanical properties and longevity. Rehabilitation strategies following microfracture surgery are meticulously crafted to facilitate recovery and functional restoration. These strategies include early mobilization techniques, controlled mechanical loading, biological augmentation, nutritional support, and advanced therapeutic modalities. Early mobilization with gentle passive range-of-motion exercises promotes synovial fluid circulation within the joint, preventing adhesions and joint stiffness, and delivering essential nutrients and growth factors. Controlled mechanical loading through gradual weight-bearing activities helps condition the repaired tissue, stimulating chondrocyte proliferation and ECM production while ensuring proper integration and function. Biological augmentation, including platelet-rich plasma (PRP) and hyaluronic acid (HA), supports cartilage repair by enhancing cellular activities and ECM formation. PRP, rich in growth factors, promotes cell proliferation and matrix synthesis, while HA injections improve joint lubrication, reduce pain, and provide a scaffold for new tissue growth. Nutritional support, including amino acids, vitamins, minerals, glucosamine, chondroitin sulfate, and antioxidants, plays a vital role in collagen synthesis and overall cartilage maintenance, protecting chondrocytes from oxidative stress and modulating inflammation. Advanced therapeutic modalities such as low-intensity pulsed ultrasound (LIPUS), pulsed electromagnetic fields (PEMF), photobiomodulation therapy, and cryotherapy offer additional benefits by enhancing chondrocyte activity, reducing inflammation, and promoting ECM production. These modalities work through various molecular mechanisms, such as activating signaling pathways and modulating gene expression, to improve the quality and durability of the repaired cartilage. By comprehensively understanding the molecular and practical aspects of cartilage repair, clinicians can develop more effective treatment plans that optimize cartilage healing and improve patient outcomes. The integration of molecular insights with advanced rehabilitation techniques holds the promise of revolutionizing cartilage repair, offering hope for those suffering from debilitating joint conditions. This review aims to bridge the gap between molecular biology and clinical practice, providing a roadmap for optimizing cartilage repair strategies and ensuring long-term success for patients.

Phasins are proteins, which are located on the surface of Polyhydroxybutyrate (P3HB) granules and that carry out activities that affect the polymer metabolism. This study aimed to evaluate the effect of inactivating phasins PhbP2 or PhbP3 on the growth rate, production and molecular mass of P3HB in cultures under low or high oxygen transfer rates (OTR). The results reveal that under high OTRₘₐₓ conditions (between 8.1 and 8.9 mmol L⁻¹ h⁻¹ depending on the strain evaluated), the absence of phasins PhbP2 and PhbP3 resulted in a strong negative effect on the growth rate. In contrast, the rates of specific oxygen consumption increased in the case of cultures with strains where phasins were absent. This behavior was not observed under conditions at a low oxygen transfer rate (3.9 ± 0.71 mol L⁻¹ h⁻¹); where both growth and oxygen consumption were the same for the different strains evaluated. It was observed, at high OTRₘₐₓ that the absence of phasin PhbP3 harms the P3HB synthesis, finding a decrease of 30 % in the polymer production at the end of cultivation; however, the molecular weight of the polymer remained constant during the course of cultivation. In summary, it can be concluded that both growth and polymer synthesis were affected when phasin PhbP3 was absent from the P3HB granule, especially at high OTRₘₐₓ.

The anterior cruciate ligament (ACL) is a critical structure within the knee joint, primarily responsible for maintaining stability and enabling the knee to perform a wide range of movements essential for both everyday activities and athletic performance. It functions as a key stabilizer by preventing excessive forward movement of the tibia relative to the femur and by controlling rotational forces that occur during activities such as pivoting, cutting, and jumping. The biomechanical integrity provided by the ACL is crucial, especially during dynamic activities that place significant strain on the knee joint, making it indispensable for athletes and active individuals. However, due to its pivotal role and the considerable mechanical demands placed upon it, the ACL is particularly vulnerable to injury. ACL injuries are prevalent, especially among athletes involved in high-risk sports such as soccer, basketball, football, and skiing, where rapid changes in direction, sudden stops, and direct impacts are common. The incidence of ACL injuries is a significant concern in sports medicine, not only because of the immediate impact on an athlete's career but also due to the long-term health implications associated with such injuries. These injuries often result in complex, multifaceted complications, including persistent joint instability that can compromise the structural and functional integrity of the knee, leading to recurrent injuries and reduced performance levels. One of the more insidious consequences of ACL injuries is the loss of proprioception, the body's ability to sense the position, movement, and action of the joints. Proprioception is vital for coordinated movement and balance, and its impairment can lead to a diminished capacity to control knee movements, further increasing the risk of re-injury. Additionally, the long-term sequelae of ACL injuries often include the development of post-traumatic osteoarthritis, a degenerative condition characterized by the gradual breakdown of joint cartilage and underlying bone. This condition not only causes chronic pain and stiffness but can also significantly impair an individual's quality of life, limiting their ability to engage in physical activities and, in severe cases, leading to disability. Despite the remarkable progress made in the fields of orthopedic surgery and sports rehabilitation, the recovery process following an ACL injury remains challenging. Surgical techniques, such as ACL reconstruction, have evolved to offer more reliable outcomes, and rehabilitation protocols have been refined to promote quicker and more effective recovery. However, the reality is that many patients do not achieve the level of recovery needed to return to their pre-injury level of activity or to maintain long-term knee health. This suboptimal recovery is often attributed to the complex interplay of neurophysiological and molecular factors that are not yet fully understood. Neurophysiologically, the recovery from an ACL injury involves not only the physical repair of the ligament but also the re-establishment of neuromuscular control and the adaptation of the central nervous system (CNS) to the altered mechanics of the knee. The CNS plays a crucial role in modulating the body's response to injury, influencing the way muscles activate and coordinate to protect the injured joint. However, disruptions in neuromuscular control post-injury can lead to compensatory movement patterns that may predispose individuals to further injury or hinder their recovery. Understanding these neurophysiological changes is essential for developing rehabilitation strategies that address not just the mechanical aspects of recovery but also the re-training of the neuromuscular system. On the molecular level, the healing process of the ACL involves a cascade of cellular and biochemical events that are critical for tissue repair and regeneration. These processes include the inflammatory response, cell proliferation, and the remodeling of the extracellular matrix, all of which are regulated by various molecular signals and pathways. Advances in molecular biology have begun to unravel the complexities of these processes, offering insights into how they can be manipulated to enhance healing and reduce the risk of complications such as fibrosis or improper tissue regeneration. By integrating molecular biology with clinical practice, there is potential to develop targeted therapies that can optimize the healing environment, improve surgical outcomes, and accelerate the recovery process. This review aims to provide a comprehensive examination of the neurophysiological aspects of ACL injury and the rehabilitation process, while also exploring the potential contributions of molecular biology to improving treatment outcomes. By delving into the intricate mechanisms that underlie ACL injury and recovery, the review seeks to highlight the challenges that remain in the field and to propose new directions for research that could lead to more effective and individualized treatment approaches. The ultimate goal is to enhance our understanding of ACL injuries and to develop strategies that not only restore knee function but also preserve long-term joint health, allowing individuals to return to their desired level of activity with confidence.

Indole derivatives are key components of natural products and possesses a wide range of biological and pharmaceutical applications. Here we present the synthesis of a new indole derivative namely 2-(1-ethyl-5-nitro-1H-indole-7-carbonyl) butyl 4-methylbenzenesulfonate (5). The structural elucidation of 5 was accomplished through comprehensive spectroscopic analysis, including Fourier-transform infrared spectroscopy (FTIR), nuclear magnetic resonance (NMR), and high-resolution mass spectrometry (HRMS). Our molecular docking study revealed that compound 5 exhibits strong affinity towards Tyrosinase, making it a promising candidate as an antioxidant agent. Keywords: Indole derivative; Tosyl group; Molecular docking; Antioxidant agent.

Mutations affecting codon 172 of the isocitrate dehydrogenase 2 (IDH2) gene define a subgroup of sinonasal undifferentiated carcinomas (SNUCs) with a relatively favorable prognosis and a globally hypermethylated phenotype. They are also recurrent (along with IDH1 mutations) in gliomas, acute myeloid leukemia, and intrahepatic cholangiocarcinoma. Commonly reported mutations, all associated with aberrant IDH2 enzymatic activity, include R172S, R172T, R172G, and R172M. We present a case of SNUC with a never-before-described IDH2 mutation, R172A. Our report compares the methylation pattern of our sample to other cases from the GEO database. Hierarchical clustering suggests a strong association between our sample and other IDH-mutant SNUCs and a clear distinction between sino-nasal normal tissues and tumors. Principal component analysis (PCA), using 100 principal components explaining 94.5% of the variance, shows the position of our sample within 1.02 standard deviation of the other IDH-mutant SNUCs. Molecular modeling analysis of the IDH2 R172A versus other R172 variants provides a structural explanation to how they affect the protein active site. Our findings thus suggest that the R172A mutation in IDH2 confers a gain of function similar to other R172 mutations in IDH2, resulting in a similar hypermethylated profile.

Auxins are chemical compounds of wide interest, mostly due to their role in plant metabolism and development. Since the latter half of the XX century, synthetic auxins have been used as herbicides, presenting notable advantages, such as low toxicity in humans. As of today, extensive descriptions of syntetic and natural auxins have been made; yet studies centered on the molecular recognition process of auxins are scarce. Herein, we present a comprehensive in silico characterization of auxin receptor TIR1, using structure based methods to aid de novo design of putative ligands. In addition, we conducted binding site characterization using cosolvent molecular dynamics and enhanced sampling simulations with reference ligands. We also present a database of putative auxins and decoy compounds which can aid further studies focusing on synthetic auxin design. To the best of our knowledge, this study is the first of its kind focusing on TIR1.

Abstract: The Willows, belonging to the genus Salix, encompass a vast array of species and hybrids, totaling between 330 to 500 species along with approximately 200 hybrids. These plants are found across continents such as Africa, North America, Europe, and Asia, exhibiting diverse forms ranging from trees to shrubs and prostrate plants. Traditionally, Willows have been utilized in folk medicine due to their rich reservoir of biologically active compounds, with salicin being a notable example, serving as a prodrug for salicylic acid.The chemical composition of Willows is highly diverse, with 322 secondary metabolites identified thus far. These include various classes such as flavonoids (94 types encompassing flavonols, flavones, flavanones, isoflavones, flavan-3-ols like catechins and procyanidins, chalcones, dihydrochalcone, anthocyanins, and dihydroflavonols), phenolic glycosides (76), organic acids (28), non-phenolic glycosides (17), sterols, and terpenes (17), simple phenolics (13), lignans (7), as well as volatiles and fatty acids (69). This rich chemical di-versity underpins the various pharmacological activities exhibited by Willows. The pharmacological properties attributed to Willows are extensive and include analgesic, an-ti-inflammatory, antioxidant, anticancer, cytotoxic, antidiabetic, antimicrobial, antiobesity, neuroprotective, and hepatoprotective effects. These properties make Willows a subject of con-siderable interest in the fields of medicine and pharmacology. In summary, Willows represent a valuable natural resource with a diverse chemical composition and a wide range of pharmacological activities. This study aimed to enhance our comprehension of the medicinal characteristics of willow bark extract, focusing on its antimicrobial properties. This investigation holds significance in validating the traditional usage of willow bark and investi-gating its potential utility in contemporary medicine.

Calicioids form a research field that has encompassed ascomycetous fungi with stalked ascomata similar to those of the lichen genus Calicium. Early generic circumscriptions in calicioid lichens and fungi were mainly based on morphological and secondary chemistry information. After the introduction of molecular data, taxonomy in the group has been reconsidered. Here, based on a broad geographical sampling, Coniocybe Ach. was revised using molecular and morphological features. Three loci (ITS, LSU and rpb1) were compared to infer its phylogenetic position, and a total of 52 new sequences (14 ITS, 24 LSU and 14 rpb1) were produced. Coniocybe was revised and emended to include, apart for its type C. furfuracea, also C. brachypoda and C. confusa. In addition, a new species Coniocybe eufuracea, was described, and a key to the species of Coniocybe provided.

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

Diffuse intrinsic pontine glioma (DIPG), affecting children aged 4-7 years, is a rare, aggressive tumor that originates in the pons and then spreads to nearby tissue. DIPG is the leading cause of death for pediatric brain tumors due to its infiltrative nature and inoperability. Radiotherapy has only a palliative effect in stabilizing symptoms. In silico and preclinical studies identified ONC201, as a cytotoxic agent against some human cancer cell lines, including DIPG ones. Single crystal X-ray analysis of the complex of the human mitochondrial caseinolytic serine protease type C (hClpP) and ONC201 (PDB ID: 6DL7) allowed to identify hClpP as its main target. Hyperactivation of hClpP causes damage to mitochondrial oxidative phosphorylation and cell death. In some DIPG patients, receiving ONC201, an acquired resistance was observed. In this context, a wide program was initiated to discover original scaffolds for new hClpP activators to treat ONC201-non-responding patients. Harmaline, a small molecule belonging to the chemical class of β -carboline, was identified through Fingerprints for Ligands and Proteins (FLAP), a structure-based Virtual Screening approach. Molecular dynamics simulations and a deep in vitro investigation showed interesting information on the interaction and activation of hClpP by harmaline.

NSD3 is a member of the NSD histone methyltransferase family of proteins. In recent years, it has been identified as a potential oncogene in certain types of cancer. NSD3 gene encodes three isoforms, the full length, a short (NSD3S) and WHISTLE isoforms. Importantly, NSD3S isoform corresponds to the N-terminal of the full-length protein, lacking the methyltransferase domain. The chromosomal location of NSD3 is frequently amplified across cancer types, such as breast, lung, colon, among others. Recently, this amplification has been correlated to a chromothripsis event, that could explain the different NSD3 alterations found in cancer. Fusion proteins containing NSD3 have also been reported, such as leukemia NSD3-NUP98, and in NUT midline carcinoma (NMC), NSD3-NUT fusion. Its role as an oncogene has been described by modulating different cancer pathways through its methyltransferase activity, or the short isoform of the protein, through protein interactions. Specifically, in this review we will focus on stating the functions that have been characterized as methyltransferase dependent, and those that have been correlated with the expression of the NSD3S isoform. Altogether, there is evidence that both isoforms of NSD3 are relevant for cancer progression, establishing NSD3 as a therapeutic target. However, further functional studies are needed to differentiate NSD3 oncogenic activity as dependent or independent of the catalytic domain of the protein, as well as the contribution of each isoform and its clinical significance in cancer progression.