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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The permeation of dioxin-like pollutants, namely, chlorinated dibenzodioxins and dibenzofurans, through lipid membranes has been simulated using classic molecular dynamics (CMD) combined with the umbrella sampling approach. The most toxic forms of chlorinated dibenzodioxin and dibenzofuran, 2,3,7,8-tetrachloro-p-dibenzodioxin (TCDD) and 2,3,7,8-tetrachlorodibenzofuran (TCDF), and a dioleyl-phosphatidylcholine (DOPC) lipid membrane of 50 Å wide have been chosen for our study. The free energy profile shows the penetration process is largely favoured thermodynamically (DG≈-12 kcal/mol), with a progressively decrease of the free energy until reaching the energy minima at distances of 8Å and 9.5Å from the centre of the membrane for, respectively, TCDD and TCDF. At the centre of the membrane, both molecules display subtle local maxima with free energy differences of 0.5 and 1 kcal/mol with respect to the energy minima for TCDD and TCDF, respectively. Furthermore, the intermolecular interactions between the molecules and the lipid membrane have been characterized at the minima and the local maxima using hybrid quantum mechanics/molecular mechanics energy decomposition analysis (QM/MM-EDA). Total interaction energies of -17.5 and -16.5 kcal/mol have been found at the energy minima for TCDD and TCDF, respectively. In both cases, the dispersion forces govern the molecule-membrane interactions, no significant changes have been found at the local maxima, in agreement with the classical free energy profile. The small differences found in the results obtained for TCDD and TCDF point out the adsorption and diffusion processes through the cell membrane are not related to the different toxicity shown by these pollutants.

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.

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

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.

CO2 geological sequestration in coal seams can achieve the dual objectives of CO2 emission reduction and enhanced coalbed methane production, making it a highly promising carbon capture and storage technology. However, the injection of CO2 into coal reservoirs in the form of supercritical fluid (ScCO2) leads to complex physicochemical reactions with the coal seam, altering the properties of the coal reservoir and impacting the effectiveness of CO2 sequestration and methane production enhancement.In this paper, theoretical calculations based on ReaxFF-MD were conducted to study the interaction mechanism between ScCO2 and the macromolecular structure of both low-rank and high-rank coal, to address the limitations of experimental methods. The reaction of ScCO2 with low-rank coal and high-rank coal exhibits significant differences. At swelling stage, the low-rank coal experiences a decrease in aromatic structure and aliphatic structure, high-rank coal shows an increase in aromatic structure and a decrease in aliphatic structure, while swelling phenomenon is more pronounced in high-rank coal. At dissolution stage, low-rank coal is initially decomposed into two secondary molecular fragments, and then recombine to form new molecular structure, the aromatic structure increased and the aliphatic structure decreased. In contrast, high-rank coal occurs stretches-breakage-movement-reconnection, a reduction in aromatic structure and an increase in aliphatic structure.The primary reasons for these variations lie in distinct molecular structure compositions and the properties of ScCO2, leading to different reaction pathways of functional group and aromatic structure. The reaction pathways of functional groups and aromatic structures in coal can be summarized as follows: the breakage of the O-H bond in hydroxyl groups, the breakage of the C-OH bond in carboxyl groups, transformation of aliphatic structures into smaller hydrocarbon compounds or formation of long-chain alkenes, and various pathways involving broken, rearrangement, recombination of aromatic structures. In low-rank coal, there is a higher abundance of oxygen-containing functional groups and aliphatic structures. The breakage of O-H and C-OH chemical bonds results in the formation of free radical ions, while some aliphatic structures detach to produce hydrocarbons. Additionally, some of these aliphatic structures combine with carbonyl groups and free radical ions to generate new aromatic structures. Conversely, in high-rank coal, the lower content of oxygen-containing functional groups and aliphatic structures, along with stronger intramolecular forces, results in fewer chemical bond breakages and makes it less conducive to the formation of new aromatic structures.These results elucidate the specific deformations of different chemical groups, offering a molecular-level understanding of the interaction between CO2 and coal.

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

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.

Lung cancer is one of the most important health risks worldwide for human. Non-small cell lung cancer (NSCLC) is the most common cause of premature death from malignant disease. This study provides in-depth insights from systems biology analyses to identify molecular to inform systemic drug targeting in NSCLC. Gene expression profiles from non small cell lung cancer were analyzed with genome-scale biomolecular networks (I,e., protein-protein interaction, transcriptional and post transcriptional regulatory networks). The aim of the study was to determine the pathways and expression profile of the genes to discover molecular signature at RNA and protein levels which could serve as potential drug targets for therapeutics innovation and the identification of novel targets. Eight proteins, six TFs and seven miRNAs came into prominence as potential drug targets. The differential expression profiles of these reporter biomolecules were cross-validated by independent RNA-Seq and miRNA-Seq. Risk discrimination performance of the reporter biomolecules NPR3, JUN, PPARG, TP53, CKMT1A, SP3 and TFAP2A were also evaluated. Total 213 drugs and 7 proteins were found for non small cell lung cancer through dgidb. Among these identified drugs seven drugs such as- Gemcitabine, Carboplatin, paclitaxel, Docetaxel, Crizotinib, Bevacizumab and Gemcitabine is used for NSCLC which is approved by National Cancer Institute. The molecular signatures and repurposed drugs presented here permit further attention for experimental studies which are offer significant potential as biomarkers and candidate therapeutics for precision medicine approaches to clinical management of NSCLC.

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

Background: Colorectal cancer is considered the third most fetal among all type of cancer. Spleen tyrosine kinase (SYK) is a non-receptor type tyrosine-protein that plays crucial role in signaling mediated via immune receptor. We adopted an onco-informatics analysis to evaluate the SYK expression and prognostic value of SYK in colorectal cancer, and identification of potential phytochemicals which may inhibit overexpression of SYK protein as well as minimized colorectal cancer. Materials & Methods: Differential expression of SYK gene was analyzed using the several transcriptomic databases including Oncomine, UALCAN, GENT2 and GEPIA2. The server, cBioPortal was used to analyze mutation and copy number alterations whereas GENT2, GEPIA, OncoLnc and PrognoScan were employed to examine the survival rate. A protein-protein interaction network of SYK and co-expressed genes of SYK was conducted via GeneMANIA. Considering SYK gene encoding protein as drug target, selected phytochemicals were assessed by molecular docking using PyRx 0.8 packages. YASARA molecular dynamics simulators were applied for the post validation of the molecular docking data. Results: We have observed significant overexpression of mRNA expression levels of SYK gene colorectal adenocarcinoma (COAD) samples compared with normal tissues. Significant methylation level and various genetic alterations are assembled in SYK gene which can lead to the development of colorectal cancer. As a result, lower level of SYK expression was related to the more chances of patients’ survival by which all the outcomes from the multiple bioinformatics platforms and web resources have demonstrated the significant evidences that the SYK kinsase can possess as a potential biomarker for the treatment of colorectal cancer. Here, aromatic phytochemicals namely, Kaempferol and Glabridin targeting SYK showed more stability compared to controls and may be useful for the treatment of colorectal cancer. Conclusion: Our study showed dysregulated expression of SYK in colorectal cancer and potentiality to act as a biomarker for the prognosis of CRC. Moreover, we have shown phytochemicals (Kaempferol and Glabridin) target SYK as potential treatment strategies and drug repositioning potentiality in colorectal cancer.

Tuberculosis (TB) is a contagious disease, caused by Mycobacterium tuberculosis (MTB) that has infected and killed a lot of people in the past. At present treatments against TB are available at a very low cost. Since these chemical drugs have many adverse effects on health, more attention is now given on the plant-derived phytochemicals as potential agents to fight against TB. In this study, 5 phytochemicals, 4-hydroxybenzaldehyde, benzoic acid, bergapten, psoralen, and p-hydroxybenzoic acid, are selected to test their potentiality, safety, and efficacy against two potential targets, the MTB RNA polymerase and enoyl-acyl carrier protein (ACP) reductase, the InhA protein, using various tools of in silico biology. The molecular docking experiment, drug-likeness property test, ADME/T-test, P450 SOM prediction, pharmacophore mapping, and modeling, solubility testing, DFT calculations, and PASS prediction study had confirmed that all the molecules had the good potentiality to inhibit the two targets. However, two agents, 4-hydroxybenzaldehyde and bergapten were considered as the best agents among the five selected agents and they also showed far better results than the two currently used drugs, that function in these pathways, rifampicin (MTB RNA polymerase) and isoniazid (InhA protein). These two agents can be used effectively to treat tuberculosis.

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.

Medicinal plants possess a surplus of novel and biologically active secondary metabolites that are responsible for counteracting diseases. Traditionally, Gomphandra tetrandra (Wall.) Sleumer is used to treat mental disorders. The present research was designed to explore phytochemicals from the ethanol leaf extract of Gomphandra tetrandra (Wall.) Sleumer to identify the potential pharmacophore(s) in the treatment of neurological disorders. The chemical compounds of the experimental plant were identified through GC-MS analysis. Besides, in-vitro antioxidant activity was assessed using different methods. Furthermore, in-vivo neurological activity was assessed in Swiss-albino mice. Computer aided analysis was appraised to determine the best-fit phytoconstituent of a total of fifteen identified compounds in the experimental plant extract against beta-amyloid precursor protein. The experimental extract revealed fifteen compounds in GC-MS analysis and the highest content was 9, 12, 15-octadecatrienoic acid (Z, Z, Z). Also, the extract showed potent anti-oxidant activity in in-vitro assays. Furthermore, in in-vivo neurological assays, the extract disclosed significant (p<0.05) neurological activity. The most favorable phytochemicals as neurological agents were selected via ADMET profiling and molecular docking was studied with beta-amyloid precursor protein. Moreover, in the computer aided study, 1, 5-Diphenyl-2H-1, 2, 4-triazoline-3-thione (Pub Chem CID: 2802516) was more active than other identified compounds with strong binding affinity to beta-amyloid precursor protein. The present in vivo and in silico studies revealed neuropharmacological features of G. tetrandra leaves extract as a natural agent against neurological disorders, especially Alzheimer’s disease.

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.

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

Breast cancer (BC) is the most common type of cancer in both men and women alike, but it is more prevalent in women. Natural compounds that can modulate the oncogenic process can be considered a significant anti-cancer agent for treating BC. These natural compounds are more effective than synthetic drugs, which have profound side effects on the normal cell and resistance to cancer cells. Genistein is a type of dietary phytoestrogen included in the flavonoid group with a similar structure of estrogen that might provide a strong alternative and complementary medicine to existing chemotherapeutics drugs. Several research studies demonstrated that it can target the estrogen receptor (ER), Human epidermal growth factor receptor-2 (HER2), and Breast cancer gene-1 (BRCA-1) in multiple BC cell lines, as well as sensitize cancer cell lines to this compound when used at an optimal inhibitory concentration. Genistein effectively showed anti-cancer activities through apoptosis induction, arresting cell cycle, inhibiting angiogenesis with metastasis, reducing inflammation, mammosphere formation, tumor growth, up-regulating tumor suppressor gene, and downregulating oncogene in suppressing cancer progression in vitro and animal model study. In addition, research studies have also suggested that these phytochemicals synergistically reverse the resistance mechanism of chemotherapeutic drugs, increasing the efficacy of some chemoinformatics drugs. Our review article aims to unbox and validate the molecular pharmacology in breast tissue, cell-specific anti-cancer mechanism with synergistic activity, and possible pharmacokinetic parameters of Genistein as a potential alternative therapeutic option for the treatment of BC.

Background: Currently, limited number of therapeutic options are available to treat Diabetes mellitus, and to find a potential candidate, laboratory work takes time and also needs animal studies. So, the aim of this study was to determine the efficacy of some natural phytocompounds with the aid of molecular docking, bioinformatics and in silico drug design approaches. Method: Two proteins (Human CYP3A4 linked to metformin and Human dipeptidyl peptidase-IV) were selected and molecular docking studies were conducted using Pymol, AutoDock Vina, PyRx, and Discovery Studio. Different important pharmacokinetic parameters like ADME and toxicity data were obtained from online databases SwissADME and pkCSM program. Results: It was found that human dipeptidyl peptidase-IV (PDB ID: 4A5S) has exhibited a maximal affinity of -9.7 Kcal/mol for bryonolic acid and hesperidin, but only -6.7 kcal/mol for metformin hydrochloride. Similarly, the highest affinity of hesperidin for human CYP3A4 bound to metformin (PDB ID: 5G5J) is -10.7 Kcal/mol compared to metformin hydrochlor (-6.3 Kcal/mol). Besides, all the compounds have been documented outstanding ADMET profile, and accepted by drug-likeness or Lipinski rule. Conclusions: The present study suggested that these compounds can be further investigated in vitro and in vivo to establish them as lead compounds against Diabetes mellitus..

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

The approval of sucrose fatty acid esters (SFAEs) as food additives/preservatives has triggered enormous interest in discovering new applications for these materials. Accordingly, many researchers reported that SFAEs consist of various sugar moieties, and hydrophobic side chains are highly active against certain fungal species. The combination of chain length and site of acylation is crucial in endowing the SFAE with high antimicrobial potency against Aspergillus species. Following several important studies, we herein present the synthesis and an assessment of the effects of acylation site and chain length (i.e., C-6 vs. C-2, C-3, C-4, and long-chain vs. short-chain) on the antimicrobial activity of mannopyranoside fatty acid esters. In vitro tests revealed that the fatty acid chain length in mannopyranoside esters significantly affects the antifungal activity where C12 chains are more potent against Aspergillus species. In terms of acylation site, mannopyranoside esters with a C8 chain substituted at the C-6 position are more active in antifungal inhibition. Molecular docking also revealed that these mannopyranoside esters had comparatively better stable binding energy, and hence better inhibition, with the fungal enzymes lanosterol 14-alpha-demethylase (3LD6), urate oxidase (1R51), and glucoamylase (1KUL) than the standard antifungal drug fluconazole. Additionally, the thermodynamic, orbital, drug-likeness, and safety profiles of these mannopyranoside esters were calculated and discussed, along with the structure-activity relationships (SAR). This study thus highlights the importance of the acylation site and lipid-like fatty acid chain length that govern the antimicrobial activity of mannopyranoside-based SFAE.

The illicit use of fentanyl has led to hundreds of thousands of opioid-related deaths worldwide. Therefore, the detection of fentanyl by law enforcement and recreational users is of utmost importance. However, current detection methods are expensive, time-consuming, require special storage conditions and necessitate complex instrumentation that is generally unportable and requires skilled personnel to operate. An alternative approach would be using molecularly imprinted polymers (MIPs) as the recognition component of a handheld sensor, testing strip, or color-based assay. In this work, a molecularly imprinted polymer was made by using methacrylic acid (MAA) and ethyleneglycol dimethacrylate (EGDMA) as the functional monomer and cross-linking monomer respectively, and benzylfentanyl (Bfen) as the template. The use of benzylfentanyl is advantageous because it closely resembles the structure of fentanyl but does not have any physiological effects as a narcotic. Important studies herein determined the optimum ratio of the template:functional monomer, with subsequent evaluations of selectivity versus the commonly encountered narcotics such as heroin, cocaine, and methamphetamine. The data obtained from the HPLC analysis showed that the Bfen-MIP was successful in selectively binding the template and actual fentanyl, better than other common narcotics.

Quercetin (QUA), a flavonoid compound, is ubiquitously found in plants and has demonstrated a diverse range of biological activities. The primary objective of the present study is to assess the potential antiemetic properties of QUA using an in vivo and in silico approach. In this experiment, emesis was induced in 4-day-old chicks by orally administering copper sulfate pentahydrate (CuSO4.5H2O) at a dose of 50 mg/kg (orally). Domperidone (6 mg/kg), ondansetron (5 mg/kg), and hyoscine (21 mg/kg) were treated as positive controls (PCs), and distilled water was employed as a negative control (NC). QUA was administered orally at two distinct doses (25 and 50 mg/kg). Additionally, QUA (50 mg/kg) and PCs were administered separately or in combination to assess their antagonistic or synergistic effects on the chicks. The binding affinity of QUA and referral ligands towards 5HT3, D2, D3, and M1–M5 receptors was estimated, and ligand-receptor interactions were visualized through various computational tools. In vivo findings indicate that QUA (25 and 50 mg/kg) has a significant effect on reducing the number of retches (16.50 ± 4.65 and 10.00 ± 4.19 times) and increasing the chick latency period (59.25 ± 4.75 and 94.25 ± 4.01 seconds), respectively. Additionally, QUA (50 mg/kg) in combination with domperidone and ondansetron exhibited superior antiemetic effects, reducing the number of retches and increasing the onset of emesis-inducing time. Furthermore, it is worth noting that QUA demonstrated the strongest binding affinity towards the D2 receptor with a value of -9.7 kcal/mol through the formation of hydrogen and hydrophobic bonds. In summary, the study found that QUA exhibited antiemetic activity in chicks, potentially by interacting with the D2 receptor pathway.

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

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.

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

Copepods are the most abundant metazoans on Earth, exhibit remarkable species diversity and global distribution. Accurate species identification of copepods presents a challenge due to their cryptic nature, necessitating meticulous examination and differentiation. To address this challenge, innovative approaches utilizing molecular markers have emerged, enabling more reliable and precise species identification. Mitochondrial genes, which are characterized by high variability, have proven invaluable in uncovering taxonomically significant variations among copepod species. In this study, we present the first complete mitochondrial genome analysis of the marine calanoid copepod Bestiolina similis (Sewell, 1914). Recent investigations employing whole mitochondrial genome-based phylogenetic analyses have generated considerable interest in exploring the deep branching patterns of arthropods. The complete mitochondrial genome of B. similis spans 23,704 base pairs and includes 13 protein-coding genes and 20 transfer RNA genes. The nucleotide composition of the mitochondrial genome is characterized by 43.42% adenine, 43.22% thymine, 6.48% guanine, and 6.88% cytosine. Maximum likelihood-based phylogenetic analysis of the mitochondrial genome were constructed of the available copepod species. Mitochondrial ge-nome-based phylogenies often reveal intriguing and distinctive topologies compared to traditional phylogenetic approaches. However, the sequencing of copepod mitochondrial genomes, particularly within the calanoid copepods remains limited. Therefore, this study aims to construct a phylogenetic tree based on the extracted mitochondrial genome, incorporating available data, to enhance our understanding of the evolutionary relationships among calanoid copepods.

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

Neuromorphic computing, reconfigurable optical metamaterials operational over a wide spectral range, holographic and nonvolatile displays of extremely high resolution, integrated smart photonics and many other applications need phase-change materials (PCMs) of the next generation with better energy efficiency, wider temperature and spectral range for reliable operation compared to current flagship PCMs as Ge2Sb2Te5 or doped Sb2Te. Gallium tellurides are promising candidates to achieve the necessary requirements because of higher melting and crystallization temperatures, combined with a low switching power and fast switching rate. At the same time, Ga2Te3 and non-stoichiometric alloys appear to be atypical PCMs, characterized by regular tetrahedral structure and the absence of metavalent bonding. The sp3 gallium hybridization in cubic and amorphous Ga2Te3 is also different from conventional p-bonding in flagship PCMs, raising a question of the phase-change mechanism. Besides, gallium tellurides exhibit a number of unexpected and highly unusual phenomena as nanotectonic compression or viscosity anomaly just above melting. Using high-energy X-ray diffraction, supported by first-principles simulations, we will unravel the atomic structure of amorphous Ga2Te5 PLD films, compare it with the crystal structure of tetragonal gallium pentatelluride, and investigate electrical, optical and thermal properties of these two materials to estimate their potential for memory applications as well as for other fields.

Single-stranded breaks (SSBs) are the most frequent DNA lesions threatening genomic integrity-understanding how DNA sensor proteins recognize certain SSB types is crucial for studies of the DNA repair pathways. During repair of damaged DNA the final SSB that is to be ligated contains a 5'-phosphorylated end. The present work employed molecular simulation (MD) of DNA with a phosphorylated break in solution to address multiple questions regarding the dynamics of the break site. How does the 5'-phosphate group behave before it initiates a connection with other biomolecules? What is the conformation of the SSB site when it is likely to be recognized by DNA repair factors once the DNA repair response is triggered? And how is the structure and dynamics of DNA affected by the presence of a break? For this purpose, a series of MD simulations of 20 base pair DNAs, each with either a pyrimidine-based or purine-based break, were completed at a combined length of over 20,000 ns simulation time and compared with intact DNA of the same sequence. An analysis of the DNA forms, translational and orientational helical parameters, local break site stiffness, bending angles, 5’-phosphate group orientation dynamics, and the effects of the protonation state of the break site phosphate group provides insights into the mechanism for the break site recognition.

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

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

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.

Multivariate statistical approaches (MSA), such as principal components analysis and multidimensional scaling, seek to uncover meaningful patterns within datasets by considering multiple response variables in a concerted fashion. Although these techniques are readily used by ecologists to visualize and explain differences between study sites, they could theoretically be employed to differentiate organisms within an experimental framework while simultaneously identifying response variables that drive documented experimental differences. Therefore, MSA were used herein to attempt to understand the response of the common, Indo-Pacific reef coral Seriatopora hystrix to temperature changes using data from laboratory-based temperature challenge studies performed in Southern Taiwan. Gene expression and physiological data partitioned experimental specimens by time of sampling, treatment temperature, and site of origin upon employing MSA, signifying that S. hystrix and its dinoflagellate endosymbionts display physiological and molecular signatures that are characteristic of sampling time, site of colony origin, and/or temperature regime. These findings promote the utility of MSA for documenting biologically meaningful shifts in the physiological and/or sub-cellular response of marine invertebrates exposed to environmental change.

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.

The drug discovery process is a complex and time-consuming endeavor that can greatly benefit from advancements in computational biology. This study explores the implementation of GPU-accelerated computational methods to enhance molecular docking simulations and virtual screening techniques, two crucial components of drug discovery. By leveraging the parallel processing capabilities of Graphics Processing Units (GPUs), we significantly accelerate the analysis of protein-ligand interactions and the identification of potential lead compounds. Our approach utilizes optimized algorithms and software frameworks to achieve substantial speedups over traditional CPU-based methods. Results demonstrate that GPU-accelerated molecular docking simulations and virtual screening can reduce computation times by up to 10-fold, enabling rapid and accurate identification of promising drug candidates. This accelerated workflow has the potential to revolutionize the drug discovery pipeline, facilitating the development of novel therapeutics and reducing the time and cost associated with bringing new treatments to market.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Plants have a natural protective process of altering their genetic molecules in response to changing environments. To uncover the genetic potential of plants, it is crucial to understand how they adapt to adverse conditions by analyzing their genetic molecules. In the study, we focused on understanding the responsive genes of tomatoes under drought conditions. We analyzed RNASeq data from different Tomato genotypes, tissue types, and different drought durations. We used a time series scale to identify early and late drought-responsive gene modules and applied a machine learning method to identify the best responsive genes. We found six candidate genes of Tomato (ASCT, FLA2, BAG5, DCL2b, NFP7.3, and ADC1) that were responsive to drought. We further constructed their protein-protein interaction network to identify their potential interactors and found them drought responsive proteins. The candidate genes can help to explore the adaptation of tomato plants under drought conditions. The identification of these candidate genes and modules can have far-reaching implications for molecular breeding and genome editing in Tomato, providing insights into the molecular mechanisms that underlie drought adaptation. This research underscores the importance of the genetic basis of plant adaptation, particularly in changing climates and growing populations.

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.

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.

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.

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.

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.

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 .

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

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

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.

There is considerable interest in cellulose solubilisation in deep eutectic solvents (DESs) as a new emerging class of green solvents. However, the experimental data have indicated that cellulose is almost insoluble in most of DESs. Carboxylic acid-based DESs are a category of DESs that exhibit a modest ability to solubilise cellulose. However, the molecular-level reasons why carboxylic acid-based DESs can dissolve cellulose are unknown. Such information is crucial in developing DESs with enhanced cellulose solubilising ability. Here, we use molecular dynamics simulations to explore the molecular mechanism of cellulose solubility in oxaline (choline chloride/oxalic acid) as a typical carboxylic acid-based DES and cellulose insolubility in reline (choline chloride/urea) as the most famous DES. It was revealed that oxalic acid possesses a greater ability to weaken the interaction between the chloride anion and the choline cation when compared to other urea in reline. Self-interactions of hydrogen bond donors are also highly attenuated in oxaline than in reline. These results explain the observed more significant melting point depression of oxaline in comparison to reline. We depict a molecular mechanism indicating how oxalic acid synergically works with chloride ion to dissolve cellulose and why urea cannot play such a role in reline.

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.

Electrostatic interactions have a determining role in conformational and dynamic behavior of polyelectrolyte molecules [1]. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in water solution with the most commonly used K+ or Na+ counterions were investigated using atomistic molecular dynamics (MD) simulations. Seven common force fields, AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and with the non-bonded fix (NBFIX) and electronic continuum corrections (ECC) to were studied. These corrections have bene introduced to correct for the problem of overbinding of ions to the charged groups of polyelectrolytes. Physical properties, such as molecular sizes, local structure and dynamics, were studied using two types of common counterions, potassium and sodium. The results show that in some cases, the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of ions and charged side chain groups. The local structures and dynamics are more sensitive on dihedral angle parameterization resulting in a preference for defined monomer conformations amd the type of correction used.

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 &amp; 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.

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.

Lumpy skin disease (LSD), caused by a virus within the Poxviridae family and Capripoxvirus genus, induces nodular skin lesions in cattle, spreading through direct contact and insect vectors, significantly affecting global cattle farming. Despite available vaccines, their effectiveness is limited by poor bioavailability and adverse effects. Our study aimed to identify potential inhibitors targeting LSD-associated DNA polymerase protein, selecting LSDV039 for further investigation through comprehensive analysis and computational methods. Virtual screening revealed rhein and taxifolin as potent binders among 380 phytocompounds, with respective affinities of -8.965 and -7.195 kcal/mol. Canagliflozin and tepotinib exhibited strong affinities (-9.858 and -8.856 kcal/mol) among 718 FDA-approved antiviral drugs. Molecular dynamics simulations of canagliflozin, tepotinib, rhein, and taxifolin highlighted taxifolin's superior stability and binding energy. Rhein displayed compactness in RMSD and RMSF but fluctuations in Rg and SASA, while canagliflozin demonstrated stability compared to tepotinib. The study highlights the promising potential of repurposed drugs and phytocompounds as LSD therapeutics. However, extensive validation through in-vitro, in-vivo, and clinical trials is crucial for their practical application.

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.

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.

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.

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

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.

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

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.

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.

The dispersive and polar properties of materials, and especially of polymers and copolymers play an important role in several engineering applications implying their surfaces and interfaces. The surface energetic properties of poly(styrene-co-butadiene) have never been studied. We proposed in this study an accurate determination of such properties by using the inverse gas chromatography (IGC) at infinite dilution. Background: The IGC surface technique led to the dispersive and polar properties of poly(styrene-co-butadiene) rubber (SBR) by adsorption of organic solvents at various temperatures. Methods: Our new methodology based on the thermal Hamieh model and the London dispersion interaction energy was used to determine the London dispersion surface energy, the polar acid-base surface energy, and the Lewis acid-base properties of the copolymer. Results: The different surface energy parameters of the SBR were obtained as a function of temperature from the chromatographic measurements. Conclusions: The dispersive and polar free energies of adsorption of the various n-alkanes and polar molecules on poly(styrene-co-butadiene) were determined at different temperatures. A decrease of the London dispersive surface energy and the polar Lewis acid-base surface energies of SBR was highlighted when the temperature increased. It was showed a Lewis amphoteric character of poly(styrene-co-butadiene) with a highest basic constant 10 times larger than its acidic constant.

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.

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.