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

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

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

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

Ti-Al based alloys, particularly TiAl and Ti3Al two-phase alloys, have garnered significant attention as potential replacements for various high-temperature structural materials due to their exceptional properties, including low density, oxidation resistance, and high strength at elevated temperatures. Despite these advantages, experimental studies on the microstructure evolution of Ti-Al based alloys under complex conditions remain challenging to observe and characterize. This review article examines the current research on molecular dynamics (MD) simulations of Ti-Al based alloys, focusing on two-phase Ti-Al alloys, Ti-Al amorphous alloys, Ti-Al composite materials, and the welding and multi-layer/film applications of Ti-Al alloys. The review highlights the unique capabilities of MD simulations in predicting the behavior of Ti-Al based alloys and addresses existing scientific challenges. Furthermore, the article discusses future research directions and development prospects in this field.

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

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.

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.

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

We performed for the first time to our knowledge fully ab initio molecular dynamics simulations of additive tribochemistry in boundary lubrication conditions. We consider an organophosphourus additive that has been experimentally shown to reduce friction in steel-on-steel sliding contacts thanks to the tribologically-induced formation of an iron phosphide tribofilm. The simulations allow us to observe in real time the molecular dissociation at the sliding iron interface under pressure and to understand the mechanism of iron phosphide formation. We discuss the role played by the mechanical stress by comparing the activation times for molecular dissociation observed in the tribological simulations at different applied loads with that expected on the basis of the dissociation barrier.

Mutation of a single amino acid residue may significantly affect the structure and function of an entire protein. The effect of single amino acid substitutions can be assessed by examining the physicochemical environment surrounding the amino acid of interest, an emerging form of quantification of which is multidimensional tensors. However, the effect with respect to a protein variant’s inherent dynamics in tensor space is rarely assessed despite the potential importance of this form of analysis in revealing local physicochemical properties of the protein and response to mutation. Using the wild-type and 936 mutant structures of the protein domain 1pga, the present research evaluated the effects of local protein context and single amino acid substitutions on molecular dynamics simulation-derived structural distributions via the use of tensors capturing a range of biochemical properties. It was observed that the extent of simulated physicochemical variation local to a substituted amino acid is positively associated with local mechanical stiffness, loss of protein thermostability and decreased local hydrophobicity. In addition, it was observed that the largest tensor variation occurs in densely-packed, hydrophobic core-associated regions of protein structures. In summary, the pattern of tensor change aligns with prior knowledge about protein stability and physicochemical properties.

The human serotonin transporter (hSERT) terminates neurotransmission by removing serotonin from the synaptic cleft, which is an essential process plays an important role in depression. In addition to substrate serotonin, hSERT is also the target of drugs of abuse like cocaine and clinically used antidepressants such as escitalopram and paroxetine. To date, few studies attempt to investigate the unbinding mechanism underlying the orthosteric and allosteric modulation of hSERT. The high-resolution X-ray structure of hSERT resolved recently enables us to theoretically study the unbinding of the above four ligands against the S1 or S2 site of hSERT, by means of molecular docking, molecular dynamics (MD) and potential of mean force (PMF) simulations. We proved that for either the S1 or S2 site, the other three ligands (cocaine, escitalopram and paroxetine) are much more favorable than the original substrate serotonin, whether in kinetics along the unbinding pathways or in thermodynamics at the equilibrium states. Furthermore, the S1 site is much more favorable than the S2 site, for each ligand. Interestingly, inspection revealed that there are ~ 3Å lengths between the allosteric site of serotonin and cocaine, and an unseen un-binding pathway for escitalopram at the S1 site except for verification of the broadest trail.

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

The adsorption of proteins on polymer is widely used in biosensors. Here, molecular dynamics (MD) simulation was used to study the immobilization of angiotensin converting enzyme II (ACE2) with six initial orientations proposed on polystyrene (PS) at the ambient conditions of pH (4.5, 6, 7, 8, 9.5) and NaCl (0.01, 0.05, 0.1, 0.15, 0.2, 0.25 M). ACE2 immobilization under favorable ambient conditions was characterized by minimum distance (short), settlement time (fast), interaction energy (substantial) and protein configuration (stable). ACE2 orientations proposed in 0.15M NaCl were respectively preferable to (90, 0, 0), (0, 0, 0) and (0, 270, 0), (180, 0, 0) and (0, 90, 0), (90, 0, 0) at pH 4.5, 6, 7, 9.5. ACE2 immobilization was further evaluated at pH 7 by optimizing NaCl concentration. Its proposed orientations of (i) (0, 270, 0), (ii) (0, 0, 0) and (90, 0, 0), and (iii) (0, 90, 0) and (90, 0, 0) were preferable in 0.05, 0.1 and 0.2 M NaCl, respectively. The great significance of cubic-orientation settlement mode provides tangible improvement for the microfabrication of biochips for rapid diagnosis of severe acute respiratory syndrome coronavirus II (SARS-CoV-2).

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.

Alzheimer's disease (AD) is a progressive neurodegenerative brain disorder. One of the important therapeutic approaches of AD is the inhibition of β‐site APP cleaving enzyme‐1 (BACE1). This enzyme plays a central role in the synthesis of the pathogenic β-amyloid peptides (Aβ) in Alzheimer's disease. A group of potent BACE1 inhibitors with known x-ray structures (PDB ID 5i3X, 5i3Y, 5iE1, 5i3V, 5i3W, 4LC7, 3TPP) were studied by molecular dynamics simulation and binding energy calculation employing MM_GB(PB)SA. The calculated binding energies gave Kd values 0.139 µM, 1.39 nM, 4.39 mM, 24.3 nM, 1.39 mM, 29.13 mM and 193.07 nM, respectively. These inhibitors showed potent inhibitory activities in enzymatic and cell assays. The Kd values were compared with experimental values, the structures were discussed in view of the energy contributions to binding. Drug likeness of these inhibitors is also discussed. Accommodation of ligands in the catalytic site of BACE1 is discussed depending on the type of fragment involved in each structure. Molecular dynamics (MD) simulations and energy studies were used to explore the recognition of the selected BACE1 inhibitors by Asp 32, Asp228 and the hydrophobic flap. The results show that selective BACE1 inhibition may be due to the formation of strong electrostatic interactions with Asp32 and Asp228 and a large number of hydrogen bonds, π-π and Van der Waals interactions with the amino acid residues located inside the catalytic cavity. Interactions with the ligands show a similar binding mode with BACE1. These results help to rationalize the design of selective BACE1 inhibitors.

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

Molecular dynamics(MD) simulations are carried out to characterize the mechanicalproperties and failure behavior of van der Waals heterostructures composed ofgraphene and hexagonal boron nitride(hBN) single layer sheets. The graphene–hBNand hBN–graphene–hBN heterostructures simulations are carried out under tensileand shear deformation. Accordingly, stress versus strain curves of each systemare plotted and various properties of heterostructures, namely elastic modulusand shear modulus as well as failure stress and failure strain, are evaluatedand compared with one another as well as with the pristine graphene and hBNsheets. Subsequently, the failure mechanism/characterization of each sheet in theheterostructures is described. Alternatively, the elastic and shear modulus of eachheterostructure are computed by means of rule of mixture(ROM) which are in goodagreement with results that are obtained from MD simulations.

The analysis of high-temperature physicochemical properties of theK₃AlF₆-Na₃AlF₆-AlF₃ molten salt system is crucial for practical production of electrolytic aluminum. In this paper, the physicochemical properties of the K₃AlF₆-Na₃AlF₆-AlF₃ molten salt electrolyte system at 1173 K and standard atmospheric pressure were simulated using molecular dynamics calculations. The effects of Zr、B、C、W、Ti、Hf、Nb、and Ta on the radial distribution function, coordination number, viscosity and conductivity of the electrolyte system are discussed in detail. The simulation results show that the coordination number between Al-F and the conductivity is decreasing as the content of B²- increases. And the coordination number between Al-F and the viscosity is decreasing as the content of Zr⁴⁺ increases. With the addition of impurity elements such as C, W, Ti, Hf, Nb, and Ta, the data show that C has no significant effect on this electrolyte system, while the presence of W will decrease the coordination number between Al-F in the system. The presence of Ti, Hf, Nb, and Ta may intensify the decomposition reaction between one of the Al-F ligands, which increases the coordination number between Al-F and the number of ions in the electrolyte. This work provides theoretical support for the subsequent study of inert anodes containing Zr, B, Ti, Hf, Nb, and Ta.

Nitrile butadiene rubber (NBR) and its various composites are extensively studied as matrix materials for oil-free lubricated mechanical friction surfaces. While numerous macroscopic friction experiments can compare the frictional performance of different NBR composites, the complex material preparation and low experimental efficiency pose challenges. Nanoscale molecular dynamics (MD) simulations can help develop materials with improved friction properties, but they cannot directly provide the friction coefficients of actual macroscopic friction surfaces. Therefore, it is important to integrate the advantages of macroscopic and nanoscale friction studies and perform a synergistic analysis to modify the friction properties of composite materials. In this study, a multi-scale approach is proposed to simulate the frictional characteristics of carbon nanotube (CNT)-reinforced NBR by combining the MD friction simulation， micromechanical bridging model, and macroscopic finite element analysis (FEA) methods. The results of the multi-scale friction simulation of copper (Cu)-CNT/NBR composites show that the addition of CNTs significantly improves the frictional properties of NBR, and the extent of the improvement is related to the mass fraction of CNTs. As the mass fraction of CNTs increases (0%, 1.25%, 2.5%, 5%), the coefficient of friction decreases from 0.50 to 0.38. The results are validated by ring-block friction experiments at the same scale as the simulations.

Small molecule modulators of neurotensin receptor 1 (NTSR1), a class A G protein-coupled receptor (GPCR), has been emerging as promising therapeutics for psychiatric disorders and cancer. Interestingly, a chemical group substitution of NTSR1 modulators can launch different downstream regulation, highlighting the significance of deciphering the internal fine-tuning mechanism. Here, we conducted synergistic application of Gaussian accelerated molecular dynamics simulation, conventional molecular dynamics simulation, and Markov state models (MSM) to investigate the underpinning mechanism of ‘driver chemical groups’ of modulators triggering inverse signaling. The result indicated that the flexibility of leucine moiety in NTSR1 agonists contributes to the inward displacement of TM7 through a loosely coupled allosteric pathway, while the rigidity of adamantane moiety in NTSR1 antagonists leads to unfavorable downward transduction of agonistic signaling. Furthermore, we found that R3226.54, Y3196.51, F3537.42, R1483.32, S3567.45 and S3577.46 may play a key role in inducing the activation of NTSR1. Together, our findings not only highlight the ingenious signal transduction within class A GPCRs, but also lay a foundation for the development of targeted drugs harboring different regulatory function of NTSR1.

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.

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 main regulatory factors during the adaptation of cancer cells to hypoxic stress are the hypoxia-inducible factors (HIFs), which are being increasingly recognized as an interesting and challenging target for the design of new chemotherapeutic molecules. HIF2A was found to have an large internal hydrophilic cavity within its PAS-B domain, unique to this sub-unit and is suggested to be a possible ligand-binding site. Regulation of HIF2A by cellular molecules is still greatly unknown. In This paper we have employed in-silico techniques, such as molecular docking and dynamic simulation, to design new direct inhibitors against HIF-2A subunit via targeting one of its critical domains and the final top screened molecules have been tested on hypoxic cancer cells for further validation of their inhibitory potential. we targeted the hydrophilic cavity inside the PAS-B domain of the HIF2A to identify novel molecules with a high binding capacity. Virtual Screening methodology was used for molecular docking of NSC library against the target domain inside the HIF2A PAS-B domain with the top 5% compounds with significant MolDock and Re-rank scores were selected for further analysis. The NSC 106416, NSC 217021, NSC 217026, and NSC 215639 compounds were selected based on their docking scores. NSC 215639 had the minimum polar solvation energy and also had a relative strong binding energy. NSC 217026 had the strongest binding energy among other compounds.

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.

This work is devoted to consideration and analysis of the application of molecular dynamics simulation (MDS) methods to the study of nanosized polymer polyvinylidene fluoride (PVDF) thin ferroelectric films (two-dimensional ferroelectrics) and their composites with graphene layers: 1) to study and calculations of the polarization switching time depending on the electric field and PVDF film thickness; 2) to study and calculations of the polarization switching time depending changes of the PVDF on PVDF-TrFE film; 3) to study the polarization switching time in PVDF under influence of graphene layers. All calculations at each MDS step were carried out using quantum semi-empirical methods PM3. A comparison and analysis of the results of these calculations and the kinetics of polarization switching within the framework of the Landau-Ginzburg-Devonshire theory for homogeneous switching in ferroelectric polymer films is carried out. The study of the composite heterostructures of the “graphene-PVDF” type and calculations of their polarization switching times were presented too. It is shown that replacing PVDF with PVDF-TrFE significantly changes the polarization switching times in these thin polymer films, and that introducing various graphene layers into the PVDF layered structure leads to both an increase and a decrease in the polarization switching time. Here everything also depends on the position and displacement of the ferroelectric coercive field depending on the system damping parameters.

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

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

Molecular dynamics simulations were carried out to investigate the atomic structure of a model Cu64Zr36 bulk metallic glass (BMG). It is found that the amount of icosahedral content of the system is significantly increased in a well relaxed structure. While we considered four connection types of vertex-, edge-, face-, and volume-sharing, the huge cluster in the relaxed samples mainly involve volume-type connection and exhibits a remarkable athermal plasticity that great stiffness and great yield strength compared to the as-quenched samples. In addition, the bond-angle distribution of annealed sample shows sharp peaks at specific bond angles which is an evidence of crystallized Laves-phase formed by icosahedral atoms, however the peaks are to be broaden after loading, which indicates decreasing amount of icosahedral content and their shape distortion. These results suggest that icosahedral content in a bulk metallic glasses plays a key role to determine the mechanical properties such as rigidity and maximum stress carrying capacity.

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.

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.

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.

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.

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.

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.

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

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

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

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.

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

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

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

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.

There are currently worldwide efforts to reduce sugar intake due to the various adverse health effects linked with the overconsumption of sugars. Artificial sweeteners have been used as an alternative to nutritive sugars in numerous applications; however, their long-term effects on human health remain controversial. This led to a shift in consumer preference towards non-caloric sweeteners from natural sources. Thaumatins are a class of intensely sweet proteins found in arils of the fruits of the West-African plant Thaumatococcus danielli. Thaumatins’ current production method through aqueous extraction from this plant and uncertainty of the harvest from tropical rainforests limits its supply while the demand is increasing. Despite successful recombinant expression of the protein in several organisms, no large-scale bioproduction facilities exist. We present preliminary process design, process simulation, and economic analysis for a large-scale (50 metric tons/year) production of thaumatin II variant by several different molecular farming platforms.

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

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

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

Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To answer this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In the present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32 % improvement over the docked structures in terms of the change of root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements was also discussed to help their implementation in future methods and applications.

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

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

Rotavirus (RV) replicates within viroplasms, membrane-less electron-dense globular cytosolic inclusions with liquid-liquid phase properties. In these structures occur the virus transcription, replication, and packaging of the virus genome in newly assembled double-layered particles. The viroplasms are composed of virus proteins (NSP2, NSP5, NSP4, VP1, VP2, VP3, and VP6), single and double-stranded virus RNAs, and host components such as microtubules, perilipin-1 and chaperonins. The formation, coalescence, maintenance, and perinuclear localization of viroplasms rely on their association with the cytoskeleton. A stabilized microtubule network involving microtubules and kinesin Eg5 and dynein molecular motors is associated with NSP5, NSP2, and VP2, facilitating dynamic processes such as viroplasm coalescence and perinuclear localization. Key post-translation modifications, particularly phosphorylation events of RV proteins NSP5 and NSP2, play pivotal roles in orchestrating these interactions. Actin filaments also contribute, triggering the formation of the viroplasms through the association of soluble cytosolic VP4 with actin and the molecular motor myosin. This review explores the evolving understanding of RV replication, emphasizing the host requirements essential for viroplasm formation and highlighting their dynamic interplay within the host cell.

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

Complex colloidal fluids, depending on particulates’ shapes and packing fractions, may have a wide range of shear thinning and thickening behaviors. A particular interesting way to transition between different types of such behavior is by infusing functional complex particles that can be manufactured using modern techicques such as 3D printing. In this paper, we display 2D molecular dynamics simulations of such fluids with infused star-shaped functional particles, with variable leg length and number of legs, as they are infused in a non-interacting, coarse-grained fluid. We vary the packing fraction (ϕ) of the system, and for each different system we apply shear with various strain rate that turns the fluid into a jammed state and rise the apparent viscosity of fluid. We demonstrate the dependence of viscosity with the particles’ packing fraction. We show the role of shape and design dependence of the functional particles towards the transition to a shear thickening fluid .

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

The 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.

The aim of this study was the reconstruction of SARS-CoV-2 evolutionary dynamics in time and space in Italy and Europe between February and June 2020. The cluster analysis showed that pure Italian clusters were observed mainly after the lockdown and distancing measures were adopted. Lineage B and B.1 spread between late January and early February 2020, from China to Veneto and Lombardy, respectively. Lineage B.1.1 most probably evolved within Italy and spread from central to south Italian regions, and to European countries. The lineage B.1.1.1 entered Italy only in the second half of March and remained localized in Piedmont until June 2020. In conclusion, the reconstructed ancestral scenario suggests a central role of China and Italy in the widespread diffusion of the D614G variant in Europe in the early phase of the pandemic and more dispersed exchanges involving several European countries from the second half of March 2020.

In this study, modified red mud after phosphorus adsorption was used as the adsorbent, and hydrochloric acid and deionized water were used as desorbents to desorb phosphorus. The components in the adsorbent were optimized based on density functional theory, and adsorbent and desorbent models were established. Molecular dynamics simulation was performed to determine the phosphorus concentration before and after desorption, interaction energies, radial distribution function(RDF), mean-square displacement(MSD), and diffusion coefficient. The Monte Carlo method was used to simulate the desorption isotherm, desorption site, heat of desorption, and desorption energy. Simulation results showed that deionized water could only desorb phosphorus on the adsorbent surface, and the stability of the system deteriorated upon adding hydrochloric acid. Hydrochloric acid destroyed the ionic and hydrogen bonding between the O atoms in H2PO4− and reactive metal and oxygen atoms in the activated red mud particles. Moreover, the van der Waals force decreased considerably. The ionic and hydrogen bonds between H2PO4− and the surface of activated red mud particles were broken by hydrochloric acid, which accelerated the desorption of phosphorus from the adsorbent surface. The interaction between hydrochloric acid and phosphorus accelerated the diffusion, which decreased the adsorption capacity. Moreover, the desorption capacity increased with increasing temperature.

The chemical feedback between building blocks in templated polymerization of diblock copolymers and their consecutive micellization was studied for the first time by means of coarse-grained molecular dynamics simulations. Using a stochastic polymerization model, we were able to reproduce the experimental findings on the effect of chemical feedback on the polymerization rates at low and high solution concentrations. The size and shape of micelles were computed using a newly development software in python conjugated with graph theory. In full agreement with the experiments, our simulations revealed that micelles formed by the templated micellization are more spherical and have lower radius of gyration than those formed by the traditional two-step micellization method. Understanding the underlying mechanisms in templated reaction/assembly of polymers will help for rational design of new synthetic supramolecular materials.

The diffusion process of water molecules within a polyetherimide (PEI) glassy matrix has been analyzed by combining the experimental analysis of water sorption kinetics performed by FTIR spectroscopy with theoretical information gathered from Molecular Dynamics simulations and with the expression of water chemical potential provided by a non-equilibrium lattice fluid model able to describe the thermodynamics of glassy polymers. This approach allowed to construct a convincing description of the diffusion mechanism of water in PEI providing molecular details of the process related to the effects of the cross- and self-hydrogen bondings established in the system on the dynamics of water mass transport.

The conformational state of the activation loop (A-loop) is pivotal for the activation of most protein kinases. Hence, the characterization of the conformational dynamics of the A-loop is important to increase our understanding of the molecular processes related to diseases and to support the discovery of small molecule kinase inhibitors. Here, we carry out a combination of Molecular Dynamics (MD) and Essential Dynamics (ED) analyses to fully map the effects of phosphorylation, ADP, and conformation disrupting (CD) inhibitors (i.e., CD532 and MLN8054) on the dynamics of the A-loop of Aurora-A. MD revealed that the stability of the A-loop in an open conformation is enhanced by single phospho-Thr-288, while paradoxically, the presence of second phosphorylation at Thr-287 decreases such stability and renders the A-loop more fluctuant in time and space. Moreover, we found that such post-translational modification has a significant effect on the direction of A-loop motions. ED analysis suggests that the presence of the phosphate moiety induces the dynamics of Aurora-A to sample two distinct energy minima, instead of a single large minimum, as in unphosphorylated Aurora-A states. This suggests that the conformational distributions of Aurora-A with both single and double phospho-threonine modifications are remarkably different from the unphosphorylated state. In the closed states, binding of CD532 and MLN8054 inhibitors has the effect of increasing the distance of the N- and C-lobes of the kinase domain of Aurora-A, and the angle analysis between those two lobes during MD simulations showed that the N- and C-lobes are kept more open in presence of CD532, compared to MLN8054. As the A-loop is a common feature of Aurora protein kinases, our studies provide a general description of the conformational dynamics of this structure upon phosphorylation and different ligands binding.

One of the strategies to reduce the contents of the low density lipoproteins (LDLs) in a blood is a hemoperfusion, in which case they are selectively retracted from plasma by an adsorber located outside the patient body. Recently, a photo-controllable smart surface was developed for this purpose characterised by high selectivity and reusability [C.Guo et al., ACS Applied Materials & Interfaces 2022, DOI:10.1021/acsami.2c07193]. We present a mesoscopic model for such a setup involving the azobenzene-containing polymer brush and the model LDL particles. The latter comprise an uniform spherical core covered by a shell of elongated particles representing phospholipids. The system is studied using the coarse-grained molecular dynamics simulation. We examined the dependences of the binding energy on both the length of polymer chains and the grafting density of a brush, and established optimal conditions for the adsorption. These are explained by a competition between the concentration of azobenzenes and phospholipids in the same spatial region, flexibility of polymer chains, and excluded volume effects.

The selectivity in the simultaneous detection of ascorbic acid (AA), dopamine (DA), and uric acid (UA) has been an open problem in the biosensing field. Many surface modification methods were carried out for glassy carbon electrodes (GCE), including the use of graphene oxide and amino acids as a selective layer. In this work, molecular dynamics (MD) simulations were performed to investigate the role of serine oligomers on the selectivity of the AA, DA, UA analytes. Our models consisted of a graphene oxide (GO) sheet under a solvent environment. Serine tetramers were added into the simulation box and were adsorbed on the GO surface. Then, the adsorption of each analyte on the mixed surface was monitored from MD trajectories. It was found that the adsorption of AA was preferred by serine oligomers due to the largest number of hydrogen-bond forming functional groups of AA, while UA was the least preferred due to its highest aromaticity. Finally, the role of hydrogen bonds on the electron transfer selectivity of biosensors was discussed with some previous studies.

A derivative series of 3,4-dimethoxy-β-nitrostyrene were synthesized and identified including new compound 6. The effect of antimicrobial activity of 3,4-alkyloxy modification of β-nitrostyrene was investigated. A molecular docking was also performed to obtain information about their interactions with Protein Tyrosine Phosphatase 1B (PTP1B). PTP1B containing cysteine 215 and arginine 221 as essential active residues plays a key role in signaling pathways that regulate various cell functions of microorganisms, which also act as negative regulator in signaling pathways of insulin that are involved in type 2 diabetes and other metabolic diseases. Compound 5 and 6 were the most potent as fragment of PTP1B inhibitor based on molecular docking, but compound 5 was more effective against Candida albicans. These compounds interact with serine 216 and arginine 221 residues. However, further research is needed to investigate their potential medicinal use.

The World Health Organization (WHO) documents malaria as one of the leading causes of high morbidity and mortality worldwide. The disease affects millions and kills thousands of people annually. Efforts to reduce the global burden of malaria have prompted WHO to recommend prevention strategies like using anti-malarial drugs. However, these strategies have been ineffective because of anti-malarial drug resistance. The only efficacious malaria treatment is Artemisinin-based Combination Therapy (ACT). However, the extended ACTs clearance times, linked to the emergence of artemisinin monotherapy resistance recorded most recently in Africa and the Great Mekong region, pose a danger to its efficacy. Therefore, better efficacious malaria drugs are required. Since P. falciparum heat shock protein 90 (PfHsp90) is a well-characterized malaria drugs target, this study uses it to discover more efficacious malaria drugs. An in-silico approach was used to discover PfHsp90 inhibitors with pharmacological properties against Plasmodium malaria using molecular dynamics simulation (MDS) and hierarchical virtual screening. Geldanamycin (GDM), a well-known anti-PfHsp90 compound, was used to identify PfHsp90 inhibitors with pharmacological properties against Plasmodium malaria by screening it against the ZINC database via the ZINCPHARMER web server. This virtual screening process resulted in 17 hits. These ZINCPHARMER hits were subjected to drug-likeness and pharmacokinetics properties analysis in the SwissADME web server, and 9 of them satisfied the requirements. The 9 ZINC compounds were docked with PfHsp90 using the PyRx software to understand their interactions. From the molecular docking results, ZINC09060002 (-8.2 kcal/mol), ZINC72133064 (-7.8 kcal/mol), ZINC72163401 (-7.7 kcal/mol), ZINC72358537 (-8.1 kcal/mol), and ZINC72358557 (-7.6 kcal/mol) had better binding affinities to PfHsp90 than GDM (-7.5 kcal/mol). The stability of these molecularly docked protein-inhibitor complexes was assessed through MDS using GROMACS 2022. ZINC72163401, ZINC72358537, and ZINC72358557 formed stable complexes with PfHsp90. The lead compounds were subjected to in vitro validation for their inhibitory capability. They showed promising inhibition of parasite growth with IC50 values ranging between 200 – 400 ng/ml. In this regard, the three PfHsp90 inhibitors can be used as anti-malarial drugs. However, further structural optimization studies and clinical (in vivo) tests are necessary to ascertain the antimalarial activity of these compounds in humans.

Cathepsin D (CatD; EC 3.4.23.5) family peptidases of parasitic organisms are regarded as potential drug targets as they play critical roles in the physiology and pathobiology of parasites. Previously, we characterized the biochemical features of cathepsin D isozyme 2 (CatD2) in the carcinogenic liver fluke Clonorchis sinensis (CsCatD2). In this study, we performed all-atomic molecular dynamics simulations by applying different systems for the ligand-free/bound forms under neutral and acidic conditions to investigate the pH-dependent structural alterations and associated functional changes in CsCatD2. CsCatD2 showed several distinctive characteristics as follows: 1) CsCatD2-inhibitor complex formed more hydrogen bonds; 2) acidic pH caused major conformational transitions from open to closed state in this enzyme; 3) neutral pH induced displacement of the N-terminal part to hinder the accessibility of the active site and open allosteric site of this enzyme; and 4) the flap dynamics metrics, including distance (d1), TriCα angles (θ1 and θ2), and dihedral angle (ϕ), account for the asymmetrical twisting motion of the active site of this enzyme. These findings provide an in-depth understanding of the pH-dependent structural dynamics of CsCatD2 and basic information for the rational design of an inhibitor as a drug targeting CsCatD2.

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

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

This study provides a comprehensive analysis of the conformational dynamics of pseudo-endocannabinoid receptors (p-CBRs) CB₁ and CB₂ in response to ligand binding. The ligands studied include N-arachidonoyl dopamine (NADA), Cannabidiol (CBD), known agonists (FUB-MDMB and E3R), and inverse agonists (Taranabant and AM630). Using homology models based on crystal structures (PDB ID: 5U09 for CB₁ and 5ZTY for CB₂), molecular docking, and molecular dynamics (MD) simulations, the research validates the structural integrity and dynamic behaviours of these pseudo-receptors. Induced-fit docking revealed NADA as having the highest binding affinity for both receptors, with docking scores of −13.92 kcal/mol for CB₁ and −13.32 kcal/mol for CB₂, facilitated by key hydrogen bond interactions. MD simulations over 250 ns at ~298 K showed successful equilibration, indicating the most stable interaction for NADA, thereby further validating the models. CB₁ and CB₂ exhibited distinct conformational responses to ligands, with Cα RMSD and Cα RMSF analyses suggesting a more rigid binding pocket for CB₁ and a flexible binding pocket for CB₂. These findings highlight the significance of receptor behaviour and dynamics in the development of cannabinoid-based therapeutics, supporting the potential of CBD and NADA as modulators for CB₁ and CB₂, respectively.

Introduction: Papatasi fever, or three-day fever, is a viral disease transmitted by a species of female sandfly called Phlebotomus papatasi that has symptoms similar to influenza. Robat Karim County, due to the high abundance of Sandfly and climate change in recent years, provides a suitable environment for the growth and reproduction of Sandfly, especially Phlebotomus papatasi Therefore, the aim of this study was molecular monitoring of Papatasi fever virus in Sandfly of Robat Karim County, Tehran province. Materials and methods: In 7 selected areas in Robat Karim County, 679 Sandfly using sampling methods including traps and aspirators were caught and 500 samples were randomly selected according to the fishing location and transferred to the laboratory of the Medical Entomology Department. After fixing Sandfly on the slide surface, identification and percentage of frequency of sandfly species were determined and then by Real-time PCR method, molecular examination of Papatasi virus infection was performed at Tehran Pasteur Institute. Findings: Out of 679 samples taken, 310 (45.65%) belonged to males and 369 (54.35%) belonged to females. 7 species including 3 species of Sandfly of the genus Phlebotomus and 4 species of the genus Sergentomyia were identified among the samples, the highest percentage of catch belonged to Phlebotomus sergenti (28.57%). Of the 500 samples sent to the laboratory for testing for phleboviruses by Real-time PCR, no positive cases of Papatasi virus infection were observed in them. Conclusion: Robat Karim is one of the endemic areas of infectious diseases such as leishmaniasis. Due to climate change in this region that conditions for the growth and reproduction of Sandfly, especially Phlebotomus papatasi as a Vector of sandfly fever and since the sandfly fever disease is not vaccinated therefore, preventive measures to control and reduce diseases transmitted by Sandfly are important.

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

This work aims to investigate the effect of Pt concentration on crystal growth mechanism of Platinum-Palladium (Pt-Pd) binary alloy system during the annealing process starting from amorphous phase until some definite temperatures. The calculations have been performed by using molecular dynamic (MD) simulations. Interatomic interactions are described by on Sutton-Chen type Embedded Atom Potential Energy function. In order to understand the main structural properties at the stable and unstable phases, changes in RDF curves versus time have been analysed for different annealing temperatures. Crystalline type bonded pairs have been determined using MD calculations which is required for the computation of Avrami coefficients and for understanding crystal growth mechanism. The results demonstrate that the increase in concentration of Pt during annealing leads to migration of atoms in the crystal lattice points, elimination of dislocations and formation of perfect crystal structure.

A better understanding of high temperature processes in slags contributes to facilitate knowledge based design of the solidified product. Here, a slag analogue with a nominal composition of 17 wt% LiMnO2 and 83 wt% Ca2SiO4 was synthesized encountering fairly high cooling rates. The Mn species from 1223 K to 1773 K was simulated using a thermodynamic model assuming a homogeneous melt. The micro-composition including the Mn species of the solidified slag was determined experimentally and was used as basis for molecular dynamics (MD) simulation. The MD simulation provides information on structure and viscosity at high temperatures, otherwise difficult to access. These parameters significantly influence oxidation state of redox-active elements and the solidified product. The micro-composition analyzed by electron probe micro analysis (EPMA) and synchrotron based micro-X-ray fluorescence (micro-XRF) showed that Mn-rich and Ca-Si-rich phases are separated. While the Mn-O phases did not contain noticeable Ca, the Ca2SiO4 phase had incorporated 0.6 wt% of Mn. The slag solidified into round shaped and droplet shaped grains of a Li-Mn-oxide, some Mn3O4 and Ca2SiO4. The powder X-ray diffraction (PXRD) confirmed the formation of larnite, the identity of the Li-Mn-oxide however remained inconclusive. The Mn oxidation state (OS) was identified using synchrotron based micro-X-ray absorption near edge spectroscopy (micro-XANES). The Mn-O grains, matched well with Li-Mn-oxides and a Mn OS: +3 e.g. LiMn3+O2. Small areas matching Hausmannite (Mn2+Mn23+O4) were also identified. The OS of Mn in the silicate phase could not be identified. For comparison a slowly cooled slag analogue with similar composition however higher Si content, was also subjected to micro-XANES. The slowly cooled slag formed long Mn-rich needles in a matrix of large calcium silicate crystals. The Mn-rich crystals matched well with the XANES spectrum of a Mn3+ Li-oxide like LiMn3+O2. At the rim of the needles the Mn-spectra matched well the Hausmannite (Mn2+Mn23+O4) reference. In the silicate phases Mn had a OS: +2, unambiguously. The melt structure at different temperatures of two compositions i.e. LiMn3+O2 and Ca2SiO4 was simulated using molecular dynamics (MD). They serve as model compositions assuming a heterogeneous melt. The results show significant different degrees of polymerization and viscosity. Information from MD simulations can support the identification of potentially different oxygen permeability and with that prediction of oxidation states. The bulk composition was identified by inductively coupled plasma optical emission spectrometry (ICP-OES), bulk structure by PXRD and bulk species by lab-XANES. The synchrotron micro analysis including micro-XRD were performed at the microfocus beamline I18 at the Diamond Light Source. Pure reference compounds were prepared and characterized with the same multi-modal approach.

Curcumin is one of the bioactive metabolites of turmeric (Curcuma longa), known for its pleiotropic pharmacological actions, including antioxidant and anti-inflammation, anticholinesterase, immunomodulation, and neuroprotection. Substantial evidence suggests the therapeutic benefits of curcumin against neurodegenerative disorders, including Alzheimer’s disease (AD), acting on a diverse array of brain targets that make the molecular mechanisms complicated. System biology level-investigation could potentially present a comprehensive molecular mechanism to delineate the neuropharmacological action of curcumin. In this study, we used integrated system pharmacology and molecular simulation analysis to gain insights into the underlying mechanism of curcumin action against AD. Network pharmacology study identified curcumin-targeted potential cellular pathways such as phosphatidylinositol 3-kinase/protein kinase B (PI3K/Akt) signaling, neurotrophin signaling, toll-like receptor (TLR) signaling, and autophagy, and proteins such as tropomyosin receptor kinase B (TrkB), liver X-receptor-beta (LXR-β), estrogen receptor-β (ER-β), mammalian target of rapamycin (mTOR), TLR-2, N-methyl-D-acetate receptor subunit 2B (GluN2B), β-secretase and glycogen synthase kinase-3β (GSK-3β), which are intimately associated with neuronal growth and survival, immune response, and inflammation. Moreover, the molecular modeling further verified that curcumin showed a significant binding affinity to mTOR, TrkB, LXR-β, TLR-2, ER-β, GluN2B, β-secretase, and GSK-3β, which are the crucial regulators of molecular and cellular processes associated with AD. Together, the present system pharmacology and in silico findings demonstrate that curcumin might play a significant role in modulating AD-pathobiology, supporting its therapeutic application for the prevention and treatment of AD.

In this study, we explored the structural and chemical modifications of cellulose fibres subjected to chemical and mechanical treatments through an innovative analytical approach. We employed photoacoustic spectroscopy (PAS) and reversed double-beam photoacoustic spectroscopy (RDB-PAS) to examine the morphological changes and the chemical integrity of the treated fibres. The methodology provided enhanced sensitivity and specificity in detecting subtle alterations in the treated cellulose structure. Additionally, we applied Coifman wavelet transformation to the PAS signals, which facilitated a refined analysis of the spectral features indicative of chemical and mechanical modifications at a molecular level. This advanced signal processing technique allowed for a detailed decomposition of the PAS signals, revealing hidden characteristics that are typically overshadowed in raw data analyses. Further, we utilized the concept of energy trap distribution to interpret the wavelet-transformed data, providing insights into the distribution and density of energy states within the fibres. Our results indicated significant differences in the energy trap spectra between untreated and treated fibres, reflecting the impact of chemical and mechanical treatments on the fibre’s physical properties. The combination of these sophisticated analytical techniques elucidated the complex interplay between mechanical and chemical treatments and their effects on the structural integrity and chemical composition of cellulose fibres.

Asymmetric lithium salts, such as Lithium (difluoromethanesulfonyl)(trifluoromethanesulfonyl)imide (LiDFTFSI), have been demonstrated to surpass traditional symmetric lithium salts with improved Li+ conductivity and the capacity to generate stable solid electrolyte interphase (SEI) while maintaining compatibility with aluminum (Al0) current collector. However, the intrinsic reductive mechanism through which LiDFTFSI influences battery performance remains unclear and under debate. Herein, the detailed SEI reactions of LiDFTFSI-based electrolytes were investigated by combining density functional theory and molecular dynamics, aiming to clarify the formation process and atomic structure of SEI. Our results show that asymmetric DFTFSI– weakens the interaction between carbonate solvents and Li+, and substantially alter the solvation structure, exhibiting a well-balanced coordination capacity to bis(trifluoromethanesulfonyl)imide (TFSI–). Nanoseconds hybrid molecular dynamics simulation further reveals that the preferential decomposition of LiDFTFSI produces sufficient LiF and Li2O to facilitate a robust SEI. Moreover, the abundant F– generated from LiDFTFSI decomposition accumulates on the Al surface and subsequently combines with Al3+ from the current collector to form AlF3, potentially inhibiting corrosion of the current collector. Overall, these findings elucidate how LiDFTFSI regulates the solvation sheath and SEI structure, advancing the development of high-performance electrolytes compatible with current collector.

This study investigates the innovative use of natural porous clays from the Bejaad region in Morocco as a support matrix for the encapsulation and controlled release of lemon essential oil (Citrus limonum, EOCL), a natural compound with well-documented insecticidal properties. The research aims to address the inherent challenges of essential oils, particularly their high volatility and rapid degradation, by improving their stability and insecticidal efficiency against the grain pest Sitophilus granarius. By anchoring EOCL onto clay matrices, the study seeks to achieve a sustained and controlled release of the active components, thereby enhancing their practical application as biopesticides. The clays were comprehensively characterized using advanced analytical techniques, including X-ray diffraction (XRD), X-ray fluorescence (XRF), Fourier-transform infrared spectroscopy (FT-IR), scanning electron microscopy with energy-dispersive X-ray analysis (SEM-EDX), and thermogravimetric analysis (TGA). These techniques revealed the mineralogical composition, thermal properties, and morphology of the clays, demonstrating their suitability for effectively adsorbing and retaining EOCL. The insecticidal performance of the clay/EOCL composites was rigorously tested under controlled conditions, revealing a marked improvement in efficacy, with significantly lower lethal doses required to achieve high mortality rates in Sitophilus granarius. The diffusion of EOCL through the clay matrix was modeled using Fick law of diffusion, and the results were further refined through statistical optimization to identify key parameters that influence the release and effectiveness of the active compounds. Complementing the experimental approach, a bioinformatics analysis was conducted to explore the molecular interactions between limonene, the primary active component of EOCL, and target proteins in insects. This theoretical investigation provided insights into the potential mechanisms of action, reinforcing the empirical findings. The study concludes that encapsulating EOCL within porous clay matrices not only enhances the stability and controlled release of the oil but also significantly boosts its insecticidal effectiveness. This approach presents a promising, environmentally sustainable strategy for crop protection, integrating material science, theoretical modeling, and bioinformatics to develop more efficient and durable biopesticides.

SOCS2, an E3 ubiquitin ligase, regulates the JAK/STAT signaling pathway, essential for cytokine signaling and immune responses. Its dysregulation contributes to cardiovascular diseases (CVDs) by promoting abnormal cell growth, inflammation, and resistance to cell death. This study aimed to elucidate the molecular mechanisms underlying the interactions between Lumbricus-derived proteins and peptides and SOCS2, with a focus on identifying potential therapeutic candidates for CVDs. Utilizing a multifaceted approach, advanced computational methodologies, including 3D structure modeling, protein-protein docking, 100 ns molecular dynamics (MD) simulations, and MM/PBSA calculations, were employed to assess the binding affinities and functional implications of Lumbricus-derived proteins on SOCS2 activity. The findings revealed that certain proteins, such as Lumbricin, Chemoattractive glycoprotein ES20, and Lumbrokinase-7T1, exhibited similar activities to standard antagonists in modulating SOCS2 activity. Furthermore, MM/PBSA calculations were employed to assess the binding free energies of these proteins with SOCS2. Specifically, Lumbricin exhibited an average ΔGbinding of -59.25 kcal/mol, Chemoattractive glycoprotein ES20 showed -55.02 kcal/mol, and Lumbrokinase-7T1 displayed -69.28 kcal/mol. These values suggest strong binding affinities between these proteins and SOCS2, reinforcing their potential therapeutic efficacy in cardiovascular diseases. Further in vitro and animal studies are recommended to validate these findings and explore broader applications of Lumbricus-derived proteins.

Shut-in after fracturing is generally adopted for wells in shale oil reservoirs, and imbibition occurring in matrix nanopores has been proven as an effective way to improve recovery. In this research, the molecular dynamics (MD) simulation was used to investigate the effects of wettability and pressure on nanopores imbibition during shut-in for a typical shale reservoir Jimsar. The results indicate that the microscopic advancement mechanism of the imbibition front is the competitive adsorption between “interfacial water molecules” at the imbibition front and “adsorbed oil molecules” on the pore wall. The essence of spontaneous imbibition involves the adsorption and aggregation of water molecules onto the hydroxyl groups on the pore wall. The flow characteristics of shale oil suggest that the overall push of the injected water to the oil phase is the main reason for the displacement of adsorbed oil molecules. Thus, shale oil especially the heavy hydrocarbon component in the adsorbed layer tends to slip on the walls. However, the weak slip ability of heavy components on the wall surface is an important reason that restricts the displacement efficiency of shale oil during spontaneous imbibition. The effectiveness of spontaneous imbibition is strongly dependent on the hydrophilicity of the matrix pore's wall. The better hydrophilicity of the matrix pore's wall facilitates higher levels of adsorption and accumulation of water molecules on the pore wall and requires less time for “interfacial water molecules” to compete with adsorbed oil molecules. During the forced imbibition process, the pressure difference acts on both the bulk oil and the boundary adsorption oil, but mainly on the bulk oil, which leads to the occurrence of wetting hysteresis. Meanwhile, shale oil still existing in the pore always maintains a good stratified adsorption structure. Because of the wetting hysteresis phenomenon, as the pressure difference increases, the imbibition effect gradually increases, but the actual capillary pressure gradually decreases and there is a loss in the imbibition velocity relative to the theoretical value. Simultaneously, the decline in hydrophilicity further weakens the synergistic effect of pressure difference because of the more pronounced wetting hysteresis. Thus, selecting an appropriate well pressure enables cost savings and maximizes the utilization of the formation's natural power for EOR.

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

Principal component analysis (PCA) is used to reduce the dimensionalities of high dimensional datasets in a variety of research areas. For example, biological macromolecules, such as proteins, exhibit many degrees of freedom, allowing them to adopt intricate structures and exhibit complex functions by undergoing large conformational changes. Therefore, molecular simulations of and experiments on proteins generate a large number of structure variations in high dimensional space. PCA and many PCA-related methods have been developed to extract key features from such structural data, and these approaches have been widely applied for over 30 years to elucidate macromolecular dynamics. This review mainly focuses on the methodological aspects of PCA and related methods, and their applications for investigating protein dynamics.

Protein dynamics play a crucial role in biological function, encompassing motions ranging from atomic vibrations to large-scale conformational changes. Recent advancements in experimental techniques, computational methods, and artificial intelligence have revolutionized our understanding of protein dynamics. Nuclear magnetic resonance spectroscopy provides atomic-resolution insights, while molecular dynamics simulations offer detailed trajectories of protein motions. The integration of machine learning, exemplified by AlphaFold2, has accelerated structure prediction and dynamics analysis. These approaches have revealed the importance of protein dynamics in allosteric regulation, enzyme catalysis, and intrinsically disordered proteins. The shift towards ensemble representations of protein structures and the application of single-molecule techniques have further enhanced our ability to capture the dynamic nature of proteins. Understanding protein dynamics is essential for elucidating biological mechanisms, designing drugs, and developing novel biocatalysts, marking a significant paradigm shift in structural biology and drug discovery.

Molecular dynamics simulation is currently the theoretical technique eligible to simulate a wide range of systems from soft condensed matter to biological systems. However, of the excellent results that the technique has arrogated, this approach remains computationally expensive, but with the emergence of the new supercomputing technologies bases on graphics processing units graphical processing units-based systems GPUs, the perspective has changed. The GPUs allow performing large and complex simulations at a significantly reduced time. In this work, we present recent innovations in the acceleration of molecular dynamics in GPUs to simulate non-Hamiltonian systems. In particular, we show the performance of measure-preserving geometric integrator in the canonical ensemble, that is, at constant temperature. We provide a validation and performance evaluation of the code by calculating the thermodynamic properties of a Lennard-Jones fluid. Our results are in excellent agreement with reported data reported from literature, which were calculated with CPUs. The scope and limitations for performing simulations of high-throughput MD under rigorous statistical thermodynamics in the canonical ensemble are discussed and analyzed.

Dysregulation of the Cyclin-dependent kinase 8 (CDK8) activity has been associated to many diseases including the colorectal and breast cancer. As usual in the CDK family, the activity of CDK8 is controlled by a regulatory protein called cyclin C (CycC). But while human CDK family members are generally activated in two steps that is, the binding of the cyclin to CDK and the phosphorylation of a residue in the CDK activation loop, CDK8 does not require the phosphorylation step to be active. Another peculiarity of CDK8 is its ability to be associated to CycC while adopted an inactive form. These specificities raise the question of the role of CycC in the complex CDK8-CycC, which appears to be more complex than for the other members of CDK family. Through MD simulations and binding free energy calculations, we investigated the effect of CycC on the structure and dynamics of CDK8 on the one hand, and the structural molecular basis of the protein-protein interaction between the two partners on the other hand. We found that CycC has a stabilizing effect on CDK8 and identified specific interaction hotspots within its interaction surface compared to other human CDK/Cyc pairs. Targeting these specific interaction hotspots could be a promising approach in terms of specificity, to effectively disrupt the interaction between CDK8. The simulation of the conformational transition from the inactive to the active form of CDK8 suggests that the residue Glu99 of the CycC may assume the role of the missing phosphorylation step in the activation mechanism of CDK8. In a more general view, these results point the importance of keeping the CycC in computational studies when studying the human CDK8 protein in both the active and the inactive form.

The molecular dynamics simulation has been used to investigate the structural and transport properties of (Ba_0.5-xSr_x)La_0.5InO_3-δ (x=0, 0.1, 0.2) oxygen-ion conductor. The previous studies reported that the ionic conductivity of Ba-doped LaInO3 decreases because Ba dopant forms narrow oxygen path in the lattice, which could hinder the diffusion of oxygen ion. In this study, we reveal the mechanism to improve the ionic conductivity by Ba and Sr co-doping on La site in LaInO3 perovskite oxide. The results show that the ionic conductivity of (Ba_0.5-xSr_x)La_0.5InO_3-δ increases with increasing numbers of Sr ions, which oxygen diffusion paths including Sr ion have larger critical radius than Ba ions. The RDF calculations showed the heights of peak in composition including Sr ions is lower and broaden, so oxygen ions moved easily into other oxygen sites.

In this paper, the nanoscale dissipative mechanisms of a Cu pad in a Ball Grid Array (BGA) packaging structure during isothermal ageing and uniaxial tension were investigated by the molecular dynamics (MD) method and experiments. From the result of the isothermal ageing test, a nonuniform consumption of Cu and large amount of Kirkendall voids were observed at the interface of Cu and Cu₃Sn. To study the effect of pressure and orientation on this phenomenon, MD simulations were conducted on four types of Cu-Cu₃Sn interface structures with different orientations of Cu. By comparing the diffusion coefficients of atoms in those cases, it was found that the tensile stress would inhibit the diffusion of atoms, whereas compressive stress would accelerate it, and this would be more significant under a larger magnitude of stress and temperature. Note that, in the model with the (101) surface Cu at the interface, both Cu and Cu₃Sn have a higher diffusion coefficient compared with the model with (001) surface Cu. Thus, the orientation of Cu will also contribute to the uniform consumption of the pad. Uniaxial tension simulation combined with DXA and CSP analyses on those models also shows the model with (001) surface Cu has a greater mechanical reliability in our simulations and related experiments.

The current outbreak of coronavirus disease (COVID-19) caused by SARS-CoV-2 in Wuhan, China has killed more than 2600 people since December 2019. Currently there is no effective treatment for this epidemic. Drug for anti SARS-CoV-2 are urgently needed. In this study we evaluated two compound libraries containing launched drugs and compounds from 300 kinds of Traditional Chinese Medicine in order to find anti SARS-CoV-2 drugs. Docking and then calculating binding free energy were performed as workflow against four key anti-SARS-CoV-2 drug targets, 3CLpro, PLpro and RdRp from SARSCoV-2, and AAK1 from human as well. As a result, drugs launched with potential for antiviral usage were selected in the hope of providing some knowledge for future drug discovery.

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

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

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

Chlorotoxin (CTX) is a 36–amino acid peptide with 8 Cys residues that forms four Cys-Cys bonds. It has high affinity for the glioma-specific chloride channel and matrix metalloprotease-2. Structural and binding properties of CTX analogs with various Cys residue substitutions with L-a-aminobutyric acid (Abu) have been previously reported. Using 4.2 ìs molecular dynamics, we compared the conformational and essential space sampling of CTX and analogs [Abu/Ser^2,19]CTX, [Abu/Ser^5,28]CTX, [Abu/Ser^16,33]CTX, [Abu/Ser^20,35]CTX, and [Abu/Ser^{2,5,16,19,20,28,33,35}]CTX. The native and substituted peptides maintained a high degree of a-helix propensity from residues 8 through 21, with the exception of [Ser^5,28]CTX and [Abu^16,33]CTX. In agreement with previous circular dichroism spectropolarimetry results, the C-terminal b-sheet content varied less from residues 25 through 29 and 32 through 36 and was well conserved in all analogs, except: [Abu^16,33]CTX, [Ser^20,35]CTX, [Abu^{2,5,16,19,20,28,33,35}]CTX, and [Ser^{2,5,16,19,20,28,33,35}]CTX. The Cys¹⁶-Cys³³ and Cys²⁰-Cys³⁵ Cys-Cys bonds appear to be required to maintain the ab-motif of CTX. Their selective substitution with Ser^16,33 however, may mitigate the destabilizing effect of Cys¹⁶-Cys³³ substitution by the formation of an inter residue H-bond from OgH of Ser¹⁶ to OgH of Ser³³ bridged by a water molecule. All peptides shared considerable sampled conformational space, which explains the retained receptor binding of the nonnative analogs.

ND1 subunit possesses the majority of the inhibitor binding domain of the human MRC-I. This is an attractive target for the search for new inhibitors that seek mitochondrial dysfunction. It is known, from in vitro experiments, some metabolites from Annona muricata called acetogenins have important biological activities such as anticancer, antiparasitic, and insecticide. Previous studies propose an inhibitory activity of bovine MRC-I by bis-THF acetogenins such as annocatacin B, however, there are few studies on its inhibitory effect on human MRC-I. In this work, we evaluate the molecular and energetic affinity of the annocatacin B molecule with the human ND1 subunit in order to elucidate its potential capacity to be a good inhibitor of this subunit. For this purpose, QM optimizations, MD simulations and MM/PBSA analysis were performed. As a control to compare our outcomes, the molecule rotenone, which is a known MRC-I inhibitor, was chosen. Our results show that annocatacin B has a greater affinity for the ND1 structure, its size and folding were probably the main characteristics that contributed to stabilize the molecular complex. Furthermore, the MM/PBSA calculations showed a 35% stronger BFE compared to the rotenone complex. Detailed analysis of the BFE shows that the aliphatic chains of annocatacin B play a key role in molecular coupling by distributing favorable interactions throughout the major part of the ND1 structure. These results are consistent with experimental studies that mention that acetogenins may be good inhibitors of MRC-I.

In its prefusion state, the SARS-CoV-2 Spike protein (similarly to other class I viral fusion proteins) is metastable, which is considered to be an important feature for optimizing or regulating their functions. After the binding process of its S1 subunit (S1) with the ACE2, the Spike protein (S) suffers a dramatic conformational change where S1 splits from the S2 subunit, which then penetrates the membrane of the host cell, promoting the fusion of the viral and cell membranes. This results in the infection of the host cell. In a previous work, we showed -using large scale molecular dynamics simulations- that the application of external electric fields (EF) induces drastic changes and damage in the receptor-binding domain (RBD) of the wild type Spike protein, as well of the Alpha, Beta and Gamma variants, leaving a structure which cannot be recognized any more by ACE2. In this work, we first extend the study to the Delta and Omicron variants and confirm the high sensitivity and extreme vulnerability of S to moderate EF (as weak as 104 V m-1), but, more importantly, we also show that, in contrast, the postfusion state of the spike protein does not suffer structural damage even if electric field intensities four orders of magnitude higher are applied. These results provide a solid scientific basis to confirm the prefusion state metastability roots of the SARS-CoV-2 Spike protein susceptibility to be damaged by EF. After the virus docks to the ACE2 receptor, the stable and robust postfusion conformation develops, which exhibits a similar resistance to EF (damage threshold higher than 108 V m-1) like most globular proteins.

A series of novel substituted isoxazolidine derivatives was synthesized through 1,3-dipolar cycloaddition reactions of α-aryl-N-methyl nitrones with various ethylene moieties, resulting in the formation of isoxazolidine derivatives. The structures of the synthesized compounds were characterized using spectroscopic methods. The [3+2] cycloaddition (32CA) reactions of α-aryl-N-methyl nitrones with diverse ethylene substituents were computed employing Density Functional Theory (DFT) at the B3LYP/6–31G(d) level, upgrading bond parameters. An investigation of the electron localization function (ELF) unveiled the zwitterionic nature of α-aryl-N-methyl nitrones, devoid of pseudoradical or carbenoid centers. Utilization of CDFT indices enabled the prediction of global electronic flux from the strongly nucleophilic α-aryl-N-methyl nitrones to the varied electrophilic ethylene functionalities. Molecular docking and molecular dynamics simulations were conducted on the isoxazolidine derivatives against the EGFR receptor (PDB: 4ZAU). In silico assessment of the synthesized isoxazolidines involved molecular docking, wherein the compound with the highest docking score was subjected to molecular dynamics simulations. Compound 11a, out of the thirty isoxazolidine derivatives, displayed the most substantial binding energy, yielding a g score of -6.57 kcal/mol. These observations highlighted the functional role of amide NH2 and NO2 in the binding process with the target protein, influencing its activity. Furthermore, prediction of ADME drug-likeness properties was performed to identify a suitable candidate.

Atomic force microscopes are used, beside their principal function as surface imaging tools, in surface manipulation and measurement of interfacial properties. In particular, they can be modified to measure lateral friction forces that occur during the sliding of the tip against the underlying substrate. However, the shape, size, and deformation of the tips profoundly affect the measurements in a manner that is difficult to predict. In this work, we investigate the contribution of these effect to the magnitude of the lateral forces during sliding. The surface substrate is chosen to be few-layer AB-stacked graphene surface whereas the tip is initially constructed from face-centered cubic gold. In order to separate the effect of deformation from shape, rigid tips of three different shapes were considered first, namely a cone, a pyramid and a hemisphere. The shape was seen to dictate all aspects of the interface during sliding, from temperature dependence to stick-slip behavior. Deformation was investigated next by comparing a rigid hemispherical tip to one of identical shape and size but with all but the top three layers of atoms free to move. The deformation, as also verified by an indentation analysis, occurs by means of the lower layers collapsing on the upper ones, thereby increasing the contact area. This collapse mitigates the friction force and decreases it with respect to the rigid tip for the same vertical distance. Finally, the size effect is studied by means of calculating the friction forces for a much larger hemispherical tip whose atoms are free to move. In this case, the deformation is found to be much smaller but the stick-slip behavior is much more clearly seen.

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

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.

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.