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

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.

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

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

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

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.

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.

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.

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

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

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

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

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

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

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

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

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

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.

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.

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

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.

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.

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 .

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.

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

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

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.

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

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

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

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

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

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.

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

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.

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.

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.

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.

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

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.

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

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.

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.

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

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

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.

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

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

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

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

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.

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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