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

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.

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

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.

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.

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.

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

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.

Markov chain Monte Carlo sampling propagators, including numerical integrators for stochastic dynamics, are central to the calculation of thermodynamic quantities and determination of structure for molecular systems. Efficiency is paramount, and to a great extent, this is determined by the integrated autocorrelation time (IAcT). This quantity varies depending on the observable that is being estimated. It is suggested that it is the maximum of the IAcT over all observables that is the relevant metric. Reviewed here is a method for estimating this quantity. For reversible propagators (which are those that satisfy detailed balance), the maximum IAcT is determined by the spectral gap in the forward transfer operator, but for irreversible propagators, the maximum IAcT can be far less than or greater than what might be inferred from the spectral gap. This is consistent with recent theoretical results (not to mention past practical experience) suggesting that irreversible propagators generally perform better if not much better than reversible ones. Typical irreversible propagators involve a parameter controlling the mix of ballistic and diffusive movement. To gain insight into the effect of the damping parameter for Langevin dynamics, its optimal value is obtained here for a multidimensional quadratic potential energy function.

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.

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.

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.

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.

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.

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

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

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.

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

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

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

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

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

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

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

The mechanical properties of α-chitin and chitosan biocomposite are investigated due to their broad range of potential applications in bioengineering, materials science, and environmental technology. Employing molecular dynamics (MD), and stress-strain analyses, this study comprehensively investigates biopolymers mechanical response and anisotropic mechanical behavior under uniaxial tensile loading conditions in an aqueous environment. Our computational models were validated against existing experimental data, showing high degrees of correspondence. Uniaxial tensile tests revealed that α-chitin exhibits a remarkable ultimate tensile strength (UTS) of 10.07 GPa in the crystalline chain direction (y-axis) at a high strain rate of 0.636, compared to a significantly lower UTS of 2.78 GPa in the perpendicular direction (x-axis) at a lower strain rate of 0.066. The biocomposite nanostructure, encompassing randomly distributed chitosan, reduced stiffness, maintaining remarkable flexibility as α-chitin, with a UTS of 5.03 GPa in the y-axis and 2.34 GPa in the x-axis. Further, the directional elastic modules for α-chitin were calculated as 51.76 GPa and 39.76 GPa in the x and y directions, respectively. In contrast, these values for the α-chitin-chitosan biocomposite were estimated as 31.66 GPa and 26.00 GPa, respectively. The findings provide valuable insights into the distinct mechanical properties of α-chitin and chitosan biocomposite, making a substantial contribution to optimizing these materials for specialized applications.

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

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

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

Spatially resolved, line-of-sight measurements of aluminum monoxide emission spectra in laser ablation plasma are used with Abel inversion techniques to extract radial plasma temperatures. Contour mapping of the radially deconvolved signal intensity shows a ring of AlO formation near the plasma boundary with the ambient atmosphere. Simulations of the molecular spectra were coupled with the line profile fitting routines. Temperature results are presented with simultaneous inferences from lateral, asymmetric radial, and symmetric radial AlO spectral intensity profiles. This analysis indicates that we measured shockwave phenomena in the radial profiles, including a temperature drop behind the blast wave created during plasma initiation.

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.

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.

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.

Lipid mesophases are being intensively studied as potential candidates for drug-delivery purposes. Extensive experimental characterization has unveiled a wide palette of release features depending on the nature of the host lipids and of the guest molecule, as well as on the environmental conditions. However, only few simulation works have addressed the matter, which hampers a solid rationalization of the richness of outcomes observed in experiments. Particularly, to date there are no theoretical works addressing the impact of hydropathy on the transport of a molecule within lipid mesophases, despite the significant fraction of hydrophobic molecules among currently-available drugs. Similarly, the high heterogeneity of water mobility in the nanoscopic channels within lipid mesophases has also been neglected. To fill this gap, we introduce here a minimal model to account for these features in a lamellar geometry, and systematically study the role played by hydropathy and water-mobility heterogeneity by Brownian-dynamics simulations. We unveil a fine interplay between the presence of free-energy barriers, the affinity of the drug for the lipids, and the reduced mobility of water in determining the net molecular transport. More in general, our work is an instance of how multiscale simulations can be fruitfully employed to assist experiments in release systems based on lipid mesophases.

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.

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.

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.

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.

Explosive boiling is a fast phase transition from an ultra-thin liquid film to vapor under an extremely high heat flux, which typically has been studied using molecular dynamics simulation (MDS) method. The present MDS study investigated the explosive boiling of a liquid argon nanofilm over different solid copper surfaces with different nanowall patterns including parallel and cross nanowalls. For each surface, atomic motion trajectories, number of liquid and vapor argon atoms, heat flux, and, mainly, onset time of explosive boiling were investigated. The simulation results indicated that explosive boiling occurs earlier on parallel and cross nanowall surfaces than on an ideally smooth surface, regardless of the topology and configuration of nanowalls. Moreover, the results revealed that by using the cross nanowall surfaces, the onset time of explosive boiling decreased by 0.7–4% compared to the parallel nanowall surfaces. In addition, it was found that the onset time of explosive boiling strongly depends on the potential energy barrier and the movement space between nanowalls for both parallel and cross nanowall surfaces. Furthermore, simulation findings showed that even though increasing the height of cross nanowalls increases heat flux and temperature of the fluid argon domain, it does not necessarily result in a shorter onset time of explosive boiling. These findings demonstrate the capability of cross nanowall surfaces for explosive boiling, thereby could be utilized in future surface design for thermal management applications.

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

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

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

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.

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.

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

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

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

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

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

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.

This study presents a detailed computational investigation aimed at identifying and evaluating potential small organic molecular scaffolds as inhibitors targeting the Pks13 protein, a critical enzyme in tuberculosis (TB) drug discovery. Our approach integrates multiple computational techniques, including quantitative structure-activity relationship (QSAR) modelling, compound library generation, molecular docking, biological activity prediction, re-evaluation of screened compounds, molecular dynamics (MD) analysis, and binding free energy estimation using the MM-GBSA approach. We developed a robust QSAR model using multiple linear regression (MLR) techniques, specifically tailored to predict the minimum inhibitory concentration (MIC) values of potential anti-TB agents. This model demonstrated excellent predictive performance, validated through rigorous internal and external validation procedures. Through molecular docking simulations, we identified five promising inhibitors Ethyl (R,Z)-2-(2-((2-chlorobenzyl) oxy)benzylidene)-5-(4-chlorophenyl)-7-methyl-3-oxo-2,3-dihydro-5H-thiazolo [3,2-a] pyrimidine-6-carboxylate (PKD1), (7R,9aR)-7-(2-hydroxyphenyl)-4-(2-methoxyphenyl)-3-methyl-7,8,9,9a-tetrahydroisoxazolo[5,4-b]quinolin-5(6H)-one (PKD2), (7-(4-methoxy phenyl)-2-methyl-3-phenylpyrazolo[1,5-a]pyrimidin-5-yl)(4-(3-(trifluoromethyl)phenyl) piperazin-1-yl)methanone (PKD3), (8S,10S)-10-(((2R,4S,5S,6S)-4-amino-5-hydroxy-6-methyl tetrahydro-2H-pyran-2-yl)oxy)-6,8,11-trihydroxy-8-(2-hydroxyacetyl)-1-methoxy-7,8,9,10-tetra-hydrotetracene-5,12-dione (PKD4), 3-(3,4-dimethylphenyl)-5-((2-(4-fluorophenyl)-2-oxoethyl) thio)-6-phenyl-2-thioxo-2,3,5,7a-tetrahydrothiazolo[4,5-d]pyrimidin-7(6H)-one (PKD5)with high binding affinities, prioritizing them based on their binding energies. These selected compounds were further assessed for their biological activity predictions, guiding the selection of candidates for experimental validation. Subsequent re-evaluation using absolute binding free energy estimation and MD simulations provided insights into the stability and dynamics of the protein-ligand complexes over time. Our findings highlighted the potential of identified compounds as stable binders within the Pks13-TE domain binding pocket, underscoring their viability as drug candidates. Additionally, binding free energy analysis using MM-GBSA reaffirmed the strong affinity of the proposed compounds towards Pks13. Our integrated computational approach offers a systematic and detailed framework for the identification and characterization of potential Pks13 inhibitors with specific relevance to TB drug discovery.

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

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

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

Cholinesterase (ChE) and secretase (BACE) inhibitors, and fibril- and β-amyloid-suppressing medicines are used to treat Alzheimer's disease (AD) symptomatically. The prevalence and complex nature of AD have increased the urgent need for multi-targeted directed ligands (MTDLs). This is because MTDLs can prevent potential drug-drug interactions during poly-therapy and have a better therapeutic profile than single targeted agents. Using piperazine linker or spacer and two primary scaffolds, benzofuran and 8-hydroxyquinoline, a unique class of multi-targeted medicines was recently reported by our group. These compounds showed re-markable effectiveness in inhibiting Aβ1-42 aggregation and acting as iron chelators.The results prompted us to synthesize an additional series of compounds as multimodal anti-AD agents and investigate its inhibitor activities to ChE (AChE/BChE) and BACE1 enzymes. The resulting in-hibitory effects suggested that the compounds under study might be used to improve cognitive function. The docking analysis' findings revealed that the compounds bind to AChE and BChE by forming H-bond interactions with amino acid residues at binding sites and µ-stacking inter-actions with aromatic residues, whereas the binding to BACE1 only revealed H-bond interac-tions with amino acid residues at binding sites.

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

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

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

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

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.

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

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

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

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

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.

We investigated the prevalence and genetic diversity of Wolbachia pipientis strains in Aedes albopictus populations in Greece. Using a combination of PCR and Sanger sequencing techniques, we genotyped Wolbachia strains in 105 mosquitoes collected across eight different prefectures in 2021. We found a high prevalence of Wolbachia in both male (90%) and female (97%) mosquitoes. Among the infected samples, 84% had double infections with both wAlbA and wAlbB strains, while 16% had infections with only wAlbB. Our comparison of the Multi-Locus Sequence Typing (MLST) profile, employing gatB-coxA-hcpA-ftsZ-fbpA genotyping, revealed a single MLST profile for each wAlbA and wAlbB strain in Greek populations. The same MLST profiles were also reported in populations from China, Russia, and Argentina, suggesting low levels of global diversity in wAlbA and wAlbB strains. Furthermore, our results indicated a significant association between temperature and the prevalence of single infections (P = 6.498E-07), with higher temperatures correlating with an increased likelihood of single infections. Although male-bias showed a tendency towards single infections, the effect was marginally non-significant (P = 0.053). These results were confirmed using a bootstrap-with-replacement analysis approach. Overall, our findings offer novel insights into the distribution and species diversity of Wolbachia strains in Greek Ae. albopictus populations, emphasizing the importance of understanding the short-term plastic and adaptive responses of these organisms to environmental stressors and rapid climate change.

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

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

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

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

In this work we discuss the insights from activity landscape, docking and molecular dynamics towards the understanding of the structure-activity relationships of dual inhibitors of major epigenetic targets: lysine metiltransferase (G9a) and DNA metiltranferase 1 (DNMT1). The study was based on a novel data set of 50 published compounds with reported experimental activity for both targets. The activity landscape analysis revealed the presence of activity cliffs, e.g., pairs of compounds with high structure similarity but large activity difference. Activity cliffs were further rationalized at the molecular level by means of molecular docking and dynamics simulations that led to the identification of interactions with key residues involved in the dual activity or selectivity with the epigenetic targets.

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

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

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

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

In this review we summarize the recent development in understanding the role of PIP2 in cellular function and signaling. We first discuss the effect of PIP2 on actin binding proteins addressing the mechanism of the actin cytoskeletal dynamics such as polymerization or depolymerization of the filamentous network or the coupling to membrane to generate forces. Next, we outline the role of PIP2 in membrane dynamics. We summarized how the membrane organization depends upon PIP2 in the presence of ions or transmembrane proteins that are sensitive to membrane curvature. We discuss how clathrin coated pits interact with adaptor proteins during the endocytosis process, which is facilitated by PIP2. Finally, we discuss the role of PIP2 in cell signaling and diseases.

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.

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

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

Newton's second law, sometimes called the fundamental principle of dynamics, although considered as an irreducible axiom of mechanics based on and validated by experience, might paradoxically depend on the energetic content of our cosmos, through the existence of a dynamical form of aether appearing, in the ancient theory of Fatio de Duillier and Lesage, as an explanation of gravitational forces. In such a framework, Newton's theory is an approximation valid for small velocities as it was the case in special relativity theory.

Most works related to the identification of mathematical nonlinear systems make one suppose that such approaches can always be directly applied to any nonlinear system. This misconception produces a great deal of discouragement when the obtained results are not the expected ones. Thus, the current work hypothesizes that the most information one has about the mathematical structure of the model, the most precise the identification result should be. Therefore, a variant of the Sparse Identification of Nonlinear Dynamics (SINDY) approach is presented to obtain the full mathematical nonlinear model of a high-order system with coupled dynamics, namely, a commercial quadcopter. Furthermore, due to its high sensitivity to inputs, a control system is devised using the identified model to stabilize the quadcopter. This illustrates the effectiveness of the proposed identification method.

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.

Fractional flow reserve (FFR) is well recognized as a gold standard measure for the estimation of functional coronary stenosis. Technological progressions in image processing have empowered the reconstruction of three-dimensional models of the coronary arteries via both non-invasive imaging modalities, such as Coronary Computed Tomography Angiography (CCTA), and invasive coronary angiography or intravascular imaging. The application of computational fluid dynamics (CFD) techniques to coronary 3D anatomical models allows the virtual evaluation of the hemodynamic significance of a coronary lesion with high diagnostic accuracy in comparison to the invasive FFR using pressure wires. These opportunities providing hemodynamic information based on anatomy have given rise to a new era of functional angiography and coronary imaging. However, further validations are needed to overcome several scientific and computational challenges before these methods are applied in everyday clinical practice. This review aims to present all the data on the studies and methodologies that have been presented for the virtual estimation of FFR, and the diagnostic execution of the virtual tools for the hemodynamic assessment of coronary lesions compared to the invasive FFR.

Candidiasis is considered an emerging public health concern because of the occurrence of drug-resistant Candida strains and the lack of available structurally diverse antifungal drug armamentariums. The indole alkaloid globospiramine from the anticandidal Philippine medicinal plant Voacanga globosa exhibits a variety of biological activities; however, its antifungal properties remain to be explored. In this study, we report the in vitro antifungal activities of globospiramine against Candida albicans and Candida tropicalis and explore its possible target proteins using in silico methods. The colony-forming unit (CFU) viability assay revealed time- and concentration-dependent anticandidal effects of the alkaloid by decreasing almost 50% of the viable CFUs 60 min after treatment. Results of the MIC and MFC assays indicated inhibitory and fungicidal effects of globospiramine against C. albicans (MIC = 8 µg/mL; MFC = 4 µg/mL) and potential fungistatic effects against C. tropicalis at lower concentrations (MIC = 4 µg/mL; MFC > 64 µg/mL). The FAM-FLICA poly-caspase assay showed metacaspase activation in C. albicans cells at concentrations of 16 and 8 µg/mL, which agreed with the MIC and MFC values. Molecular docking and molecular dynamics simulation experiments indicated stable, strong binding of globospiramine with the target proteins 1,3-β-glucan synthase and Als3 adhesin enzymes, which are indirectly involved in apoptotic-driven candidal inhibition.

This study presents the complexity and sensitivity of chaotic system dynamics in the case of the double pendulum. It applied detailed numerical analyses of the double pendulum in multiple computing platforms in order to demonstrate the complexity in behavior of the system of double pendulums. The equations of motion were derived from the Euler-Lagrange formalism, in order to capture the system's dynamics, which is coupled nonlinearly. These were solved numerically using the efficient Runge-Kutta-Fehlberg method, implemented in Python, R, GNU Octave, and Julia, while runtimes and memory usage were extensively benchmarked across these environments. Time series analyses, including the calculation of Shannon entropy and the Kolmogorov-Smirnov test, quantified the system's unpredictability and sensitivity to infinitesimal perturbations of the initial conditions. Phase space diagrams illustrated the intricate trajectories and strange attractors, as further confirmation of the chaotic nature of the double pendulum. All the findings have a clear indication of the importance of accurate measurements of the initial condition in a chaotic system, contributing to an increased understanding of nonlinear dynamics. Future research directions are faster simulations using Numba and GPU computing, stochastic effects, chaotic synchronization, and applications in climate modeling. This work will be useful for understanding chaos theory and efficient computational approaches in complex systems of dynamical nature.

State transition in the multiple-levels system has the great potential applications in the quantum technology. In this article we employ a deterministic approach in complex space to analyze the dynamics of the 1s-2p electron transition in the hydrogen atom. The electron’s spin motion is embodied in the framework of quantum Hamilton mechanics that allows us to examine the transition dynamics more precisely. The transition is driven by an oscillating electric field in the z-direction. The electron’s transition process can be visualized by monitoring its motion in the complex space. The quantum potential and the total energy proposed in this paper provide new indices to observe the dynamic changes of electrons in the transition process.

Complex human behavior, including interlimb and interpersonal coordination, has been studied from a dynamical system perspective. We review the applications of a dynamical system approach to a sporting activity, which includes continuous, discrete, and switching dynamics. Continuous dynamics identified switching between in- and anti-phase synchronization, controlled by an interpersonal distance of 0.1 m during expert kendo matches, using a relative phase analysis. As discrete dynamics, return map analysis was applied to the time series of movements during kendo matches, and six coordination patterns were classified. Furthermore, state transition probabilities were calculated based on the two states, which clarified the coordination patterns and switching behavior. We introduced switching dynamics with temporal inputs to clarify the simple rules underlying the complex behavior corresponding to switching inputs in a striking action as a non-autonomous system. As a result, we determined that the time evolution of the striking action was characterized as fractal-like movement patterns generated by a simple Cantor set rule with rotation. Finally, we propose a switching hybrid dynamics to understand both court-net sports, as strongly coupled interpersonal competition, and weakly coupled sports, such as martial arts.

We develop an action principle to construct the field equations for dissipative/resistive general relativistic two-temperature plasmas, including a neutrally charged component. The total action is a combination of four pieces: an action for a multi-fluid/plasma system with dissipation/resistivity and entrainment; the Maxwell action for the electromagnetic field; the Coulomb action with a minimal coupling of the four-potential to the charged fluxes; and the Einstein-Hilbert action for gravity (with the metric being minimally coupled to the other action pieces). We use a pull-back formalism from spacetime to abstract matter spaces to build unconstrained variations for the neutral, positively, and negatively charged fluid species and for three associated entropy flows. The full suite of field equations is recast in the so-called ``$3 + 1$'' form (suitable for numerical simulations), where spacetime is broken up into a foliation of spacelike hypersurfaces and a prescribed ``flow-of-time''. A previously constructed phenomenological model for the resistivity is updated to include the modified heat flow and the presence of a neutrally charged species. We impose baryon number and charge conservation as well as the Second Law of Thermodynamics in order to constrain the number of free parameters in the resistivity. Finally, we take the Newtonian limit of the ``$3 + 1$'' form of the field equations which can be compared to existing non-relativistic formulations. Applications include main sequence stars, neutron star interiors, accretion disks, and the early universe.

Drawing inspiration from a recent construction of a polyhedral structure associated with an icosahedrally-symmetric map on the Riemann sphere, the article shows how to build such "dynamical polyhedra" for other icosahedral maps. First, icosahedral algebra is used to determine a special family of maps with 60 periodic critical points. The topological behavior of each map is worked out and results in a geometric algorithm that constructs a system of edges---the dynamical polyhedron---in natural correspondence to a map's topology. It turns out that the maps' descriptions fall into classes the presentation of which concludes the paper.

Evolution of cooperation has traditionally been studied by assuming that individuals adopt either of two pure strategies, to cooperate or defect. Recent work have considered continuous cooperative investments, turning full cooperation and full defection into two opposing ends of a spectrum and sometimes allowing for the emergence of the traditionally-studied pure strategies through evolutionary diversification. These studies have typically assumed a well-mixed population in which individuals are encountered with equal probability, Here, we allow for the possibility of assortative interactions by assuming that, with specified probabilities, an individual interacts with one or more other individuals of the same strategy. A closely related assumption has previously been made in evolutionary game theory and has been interpreted in terms of relatedness. We systematically study the effect of relatedness and find, among other conclusions, that the scope for evolutionary branching is reduced by either higher average degree of, or higher uncertainty in, relatedness with interaction partners. We also determine how different types of non-linear dependencies of benefits and costs constrain the types of evolutionary outcomes that can occur. While our results overall corroborate the conclusions of earlier studies, that higher relatedness promotes the evolution of cooperation, our investigation gives a comprehensive picture of how relatedness affects the evolution of cooperation with continuous investments.

In this paper we introduce an ODE system to model the interaction between nature resources and human exploitation in a rich consumeristic society. In this model the rate of change of population depends on wealth per capita, and the rate of consumption has a quadratic growth with respect to population and wealth. As in our previous work \cite{BaCra_1}, we distinguish between renewable and non renewable resources, and we introduce a replenishment term for non renewable resources. We first get some information on the asymptotic behaviour of wealth and population, then we compute all the equilibria of the system and we study their stability when the parameter of exploitation of resources is small. Some numerical simulations confirm the theoretical results, and give suggestions for future research.

Reservoir computing (RC) systems can efficiently forecast chaotic time series using nonlinear dynamical properties of an artificial neural network of random connections. The versatility of RC systems has motivated further research on both hardware counterparts of traditional RC algorithms and more efficient RC-like schemes. Inspired by the nonlinear processes in a living biological brain and using solitary waves excited on the surface of a flowing liquid film, in this paper we experimentally validate a physical RC system that substitutes the effect of randomness for a nonlinear transformation of input data. Carrying out all operations using a microcontroller with a minimal computational power, we demonstrate that the so-designed RC system serves as a technically simple hardware counterpart to the `next-generation’ improvement of the traditional RC algorithm.