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.

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

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

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

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.

Natural products from plants such as, chemopreventive agents attract huge attention because of their low toxicity and high specificity. The rational drug design in combination with structure based modeling and rapid screening methods offer significant potential for identifying and developing lead anticancer molecules. Thus, the molecular docking method plays an important role in screening a large set of molecules based on their free binding energies and proposes structural hypotheses of how the molecules can inhibit the target. Several peptide based therapeutics have been developed to combat several health disorders including cancers, metabolic disorders, heart-related, and infectious diseases. Despite the discovery of hundreds of such therapeutic peptides however, only few peptide-based drugs have made it to the market. Moreover, until date the activities of cyclic peptides towards molecular targets such as protein kinases, proteases, and apoptosis related proteins have never been explored. In this study we explore the in silico kinase and protease inhibitor potentials of cyclosaplin as well as study the interactions of cyclosaplin with other cancer-related proteins. Previously, the structure of cyclosaplin was elucidated by molecular modeling associated with dynamics that was used in the current study. Docking studies showed strong affinity of cyclosaplin towards cancer-related proteins. The binding affinity closer to 10 indicated efficient binding. Cyclosaplin showed strong binding affinities towards protein kinases such as EGFR, VEGFR2, PKB and p38 indicating its potential role in protein kinase inhibition. Moreover, it displayed strong binding affinity to apoptosis related proteins and revealed the possible role of cyclosaplin in apoptotic cell death. The protein-ligand interactions using LigPlot displayed some similar interactions between cyclosaplin and peptide-based ligands especially in case of protein kinases and a few apoptosis related proteins. Thus, the in silico analyses gave an insight of cyclosaplin as a potential apoptosis inducer and protein kinase inhibitor.

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

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

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.

Some α-Pyrone derivatives isolated from Alternaria phragmospora fungus showed promising anti leukemic activities, while others were inactive. CRM1/XPO1 (chromosome region maintenance 1 protein, also called exportin1 or PO1 in humans) has been chosen as a target for antileukemic molecular docking study for those compounds to understand their modes of interaction and structure activity relationships. The results showed that two (2 and 4), out of six, natural α-Pyrone derivatives exhibited well-established interactions with the amino acids of the receptor, which was in agreement with the experimental anti-leukemic results of these compounds. Moreover, twenty hypothetical chemically modified α-Pyrone derivatives (7-27) have been designed. Compounds 7, 8, 22 and 24 showed more efficient docking properties than the previously considered natural compounds.

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.

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.

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

We have synthesized quinoxaline analogs (1-25), characterized by 1HNMR and HREI-MS and evaluated for thymidine phosphorylase inhibition. Among the series, nineteen analogs showed better inhibition when compared with the standard inhibitor 7-Deazaxanthine (IC50 = 38.68 ± 4.42 µM). The most potent analog among the series is analog 25 with IC50 value 3.20 ± 0.10 µM. Sixteen analogs 1, 2, 3, 4, 5, 6, 7, 12, 13, 14, 15, 16, 17, 18, 21and 24 showed outstanding inhibition which is many folds better than the standard 7-Deazaxanthine. Two analogs 8 and 9 showed moderate inhibition. A structure- activity relationship has been established mainly based upon the substitution pattern on the phenyl ring. The binding interactions of the active compounds were confirmed through molecular docking studies.

A sedentary lifestyle associated with unregulated diets rich in high-calorie foods have contributed to the great prevalence of metabolic dysfunction-associated steatotic liver disease (MASLD) latterly, with up to 60% in high-risk population and 25% in the general population. The absence of specific pharmacological strategies for this syndrome represents one of the major problems in the management of MASLD patients. Lifestyle interventions and adherence to a healthy diet are the main cornerstones of current therapies. The identification of nutraceuticals useful in the treatment of MASLD appears to be one of the most promising strategies for the development of new effective and safe treatments for this disease. The onion, one of the most widely studied foods in the field of nutraceuticals, serves as an inexhaustible reservoir of potent compounds with various beneficial effects. The following preliminary study analyses, mediating in silico studies, the iteration of a library of typical onion compounds with 3-hydroxy-3-methylglutaryl-coenzyme A reductase, Liver receptors X α and β as well as Peroxisome proliferator-activated receptors α and γ. In this study fort the first time promising smart molecules from the onion that could have a beneficial action in MASLD patients were identified.

In the search of potent α-amylase inhibitors, we have synthesized seventeen derivatives of 2-mercaptobenzimidazole bearing sulfonamide (1-17) and evaluated for their α-amylase inhibitory potential. All compounds display a variable degree of α-amylase activity having IC50 values ranging between 0.90 ± 0.05 to 11.20 ± 0.30 µM when compared with the standard drug acarbose having IC50 value 1.70 ± 0.10 µM. Compound 1, 2, 11, 12 and 14 having IC50 values 1.40 ± 0.10, 1.30 ± 0.05, 0.90 ± 0.05, 1.60 ± 0.05 and 1.60 ± 0.10 µM respectively were found many folds better than the standard drug acarbose. The remaining analogs showed good inhibitory potentials. All the synthesized compounds were characterized by HREI-MS, 1H and 13C-NMR. Structure activity relationship (SAR) has been recognized for all newly synthesized analogs. Through molecular docking study, binding mode of active analogs with α-amylase enzyme was confirmed.

New class of triazinoindole bearing thiosemicarbazide (1-25) was synthesized and evaluated for α-glucosidase inhibitory potential. All synthesized analogues exhibited excellent inhibitory potential having IC50 values ranging from 1.30 ± 0.01 to 35.80 ± 0.80 µM when compared with the standard acarbose having IC50 value 38.60 ± 0.20 µM. Among series the analogues 1 and 23 was found the most potent having IC50 values 1.30 ± 0.05 and 1.30 ± 0.01 µM respectively. Structure activity relationship (SAR) was mainly based upon by bring about difference of substituents on phenyl rings. To confirm the binding interactions, molecular docking study was performed. Synthesized analogues were characterized through HREI-MS, 1H and 13C-NMR analysis.

The outbreak of novel coronavirus (COVID-19) infections in 2019 is in dire need of finding potential therapeutic agents. In this study, we used molecular docking to repurpose HIV protease inhibitors and nucleoside analogues for COVID-19, with evaluations based on docking scores calculated by AutoDock Vina and RosettaCommons. Our results suggest that Indinavir and Remdesivir possess the best docking scores, and comparison of the docking sites of the two drugs reveal a near perfect dock in the overlapping region of the protein pockets. After further investigation of the functional regions inferred from the proteins of SARS coronavirus, we discovered that Indinavir does not dock on any active sites of the protease, which may give rise to concern in regards to the efficacy of Indinavir. On the other hand, the docking site of Remdesivir is not compatible with any known functional regions, including template binding motifs, polymerization motifs and nucleoside triphosphate (NTP) binding motifs. However, when we tested the active form (CHEMBL2016761) of Remdesivir, the docking site revealed a perfect dock in the overlapping region of the NTP binding motif. This result suggests that Remdesivir could be a potential therapeutic agent. Clinical trials still must be done in order to confirm the curative effect of these drugs.

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.

Background: COVID-19, a member of corona virus family is spreading its tentacles across the world due to lack of drugs at present. Associated with its infection are cough, fever and respiratory problems causes more than 15% mortality worldwide. It is caused by a positive, single stranded RNA virus from the enveloped coronaviruse family. However, the main viral proteinase (Mpro/3CLpro) has recently been regarded as a suitable target for drug design against SARS infection due to its vital role in polyproteins processing necessary for coronavirus reproduction.Objectives: The present in silico study was designed to evaluate the effect of Eucalyptol (1,8 cineole), a essential oil component from eucalyptus oil, on Mpro by docking study.Methods: In the present study, molecular docking studies were conducted by using 1-click dock and swiss dock tools. Protein interaction mode was calculated by Protein Interactions Calculator.Results: The calculated parameters such as RMSD, binding energy, and binding site similarity indicated effective binding of eucalyptol to COVID-19 proteinase. Active site prediction further validated the role of active site residues in ligand binding. PIC results indicated that, Mpro/eucalyptol complexes forms hydrophobic interactions, hydrogen bond interactions and strong ionic interactions.Conclusions: Therefore, eucalyptol may represent potential treatment potential to act as COVID-19 Mpro inhibitor. However, further research is necessary to investigate their potential medicinal use.

Curcumin, a major constituent of turmeric, has beneficial effects against several diseases. In cystic fibrosis (CF), this compound potentiates the activity of a number of cystic fibrosis transmembrane conductance regulator (CFTR) mutants. Despite holding promise in the treatment of CF, the curcumin binding site in CFTR and the molecular mechanism of activation of this channel are still unknown. The results of the study, based on docking and molecular dynamics (MD) simulations, allow to propose that curcumin binds the closed ATP-free CFTR near the NBD1/ ICl1/ICl4 interface. The bound ligand, once approached by NBD2 during transient channel opening, lays at a multiple interdomain cross point. Thereafter, curcumin can bridge NBD1 and NBD2, and also ICL1/ICL4 and ICL2/ICL3 finally tightening the same interdomain interactions that normally uphold the active conformation in the wild type ATP-bound CFTR. The proposed binding site is compatible with biochemical observations made in previous CFTR-curcumin interaction studies. These findings provide the framework for the design of novel drugs that activate CFTR mutants characterized by defects in ATP binding and/or NBD dimerization, or even lacking NBD2.

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

Leucyl-aminopeptidase (LAP), an important metallopeptidase, hydrolyses amino acid residues from the N-terminal of polypeptides and proteins, acting preferentially on the peptide bond formed by N- terminus leucine. A new leucyl-aminopeptidase was found in Bacillus cereus CZ. Its gene (bclap) contained a 1485 bp ORF encoding 494 amino acids with a molecular weight of 54 kDa. The bcLAP protein was successfully expressed in E. coli BL21(DE3). Optimal activity is obtained at pH 9.0 and 58°C. The bcLAP displays a moderate thermostability and an alkaline pH adaptation range. Enzymatic activity is dramatically enhanced by Ni2+. EDTA significantly inhibits the enzymatic activity, and bestatin and SDS also show strong inhibition. The three-dimensional model of bcLAP monomer and homohexamer is simulated by Phyre2 server and SWISS-MODEL server. The docking of bestatin, Leu-Trp, Asp-Trp and Ala-Ala-Gly to bcLAP is performed using AutoDock4.2.5 , respectively. Molecular docking results show that the residues Lys260, Asp265, Lys272, Asp283, Asp342, Glu344, Arg346, Gly372, His437 are involved in the hydrogen bonding with the ligands and zinc ions. There may be two nucleophilic catalytic mechanisms in bcLAP, one involving His 437 or Arg346 and the other involving His437 and Arg346. The bcLAP can hydrolyse the peptide bonds in Leu-Trp, Leu-Trp and Ala-Ala-Gly.

Introduction: Oral candidiasis is an opportunistic fungal infection, often caused by Candida albicans. One of the virulence factors of C. albicans is sterol 14a-demethylase or CYP51. Compounds of anthocyanins (ternatin, cyanin) and flavonoids (rutin, isoquercetin, kaempferol, myricetin, limocitrin) from the extract of the Clitoria ternatea (butterfly pea) plant are believed to have potential to inhibit the CYP51 virulence factor of C. albicans to suppress oral candidiasis. The purpose of this research was to investigate the interaction between the ligands of anthocyanins and flavonoids from C. ternatea (butterfly pea) with CYP51. Methods: The in silico study applied molecular docking to examine the molecular interactions (binding affinity) between the anthocyanin and flavonoid ligands of interest from C. ternatea and the CYP51 receptor of C. albicans. Results: Predicted interactions exist between the ligands of anthocyanins (ternatin, cyanin) and flavonoids (rutin, isoquercetin, kaempferol, myricetin, limocitrin) from C. ternatea and the CYP51 receptor of C. albicans. Conclusion: C. ternatea has the potential to provide anti-fungal agents that can inhibit or slow down the occurrence of oral candidiasis, as indicated by in silico studies. However, it is important to conduct additional research regarding the properties of the plant to determine the unknown aspects of its effects to the human body.

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.

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

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

The spread of the global COVID-19 pandemic, the lack of specific treatment and the urgent situation requires use of all resources to remedy this scourge. In the present study, using molecular docking, we identify new probable inhibitors of COVID-19 by molecules from Nigella sativa L, which is highly reputed healing herb in North African societies and both Islamic and Christian traditions. The discovery of the M^pro protease structure in COVID-19 provides a great opportunity to identify potential drug candidates for treatment. Focusing on the main proteases in CoVs (3CL^pro/M^pro) (PDB ID 6LU7 and 2GTB); docking of compounds from Nigella Sativa and drugs under clinical test was performed using Molecular Operating Environment software (MOE). Nigelledine docked into 6LU7 active site gives energy complex about -6.29734373 Kcal/mol which is close to the energy score given by chloroquine (-6.2930522 Kcal/mol) and better than energy score given by hydroxychloroquine (-5.57386112 Kcal/mol) and favipiravir (-4.23310471 kcal/mol). Docking into 2GTB active site showed that α- Hederin gives energy score about-6.50204802 kcal/mol whcih is better energy score given by chloroquine (-6.20844936 kcal/mol), hydroxychloroquine (-5.51465893 kcal/mol)) and favipiravir (-4.12183571kcal/mol). Nigellidine and α- Hederin appeared to have the best potential to act as COVID-19 treatment. Further, researches are necessary to testify medicinal use of identified and to encourage preventive use of Nigella Sativa against coronavirus infection.

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

The outbreak of Coronavirus, as well as its emerging potential consequences, has become a global challenge worldwide, demanding effective and controlled therapeutic strategies. Potential drug candidates could achieve that with minimal toxicity. The current investigation selected 2-Deoxy-D-glucose prescribed by Defense Research and Development Organization (DRDO), India. The derivative and modified form was tested through in silico analysis against COVID-19 main protease complex with N3 inhibitor. The derived form of 2-Deoxy-D-glucose was generated by replacing the hydroxy group with a hydrogen atom, and Cypate 2-Deoxy-D-glucose was chosen as a derivative against the COVID-19 main protease complex. A molecular docking approach was adopted to identify the stable and competent form among the modified and derivative forms of 2-Deoxy-D-glucose based on binding energy. By further promoting the stabilized complex, these compounds' toxicity was also scrutinized through ADMET analysis to predict the potential candidate. The current investigation suggested that the modified version of 2-Deoxy-D-glucose was more stable with minimal toxicity against COVID-19 main protease.

Corona Virus Infectious Disease-2019 (COVID-19) outbreak originated recently at Wuhan, China in December 2019. It has already spread rapidly to more than 200 countries and has been declared a pandemic by WHO. It is caused by a beta-coronavirus named as SARS-CoV-2. There is no definitive cure, either drug or vaccine, to treat or prevent this viral disease. Recently, the crystal structure of the main protease M^pro has been determined. M^pro is responsible for the proteolytic maturation of the polyprotein essential for the viral replication and transcription, which makes it an important drug target. The discovery of new drug molecules may take years before getting to the clinics. So, considering urgency we performed molecular docking studies using FDA approved drugs to identify molecules that could potentially bind to the substrate-binding site and inhibit SARS-CoV-2 main protease (M^pro). We used the Glide module in Schrodinger software suite to perform molecular docking studies followed by MM-GBSA based energy calculations to score the hit molecules. Molecular docking and manual analysis suggest that several drugs may bind and potentially inhibit M^pro. We also performed molecular simulations studies for selected compounds to evaluate protein-drug interactions. Interestingly, we observed only one antiviral compound, Adefovir, in the top50 list of compounds. Considering bioavailability, lesser toxicity, route of administration some of the top-ranked drugs including lumefantrine (antimalarial), dipyridamole (coronary vasodilator), dihydroergotamine (used for treating migraine), hexoprenaline (anti- asthmatic), riboflavin (vitamin B2) and pantethine (vitamin B5) may be taken forward for further in vitro and in vivo experiments to investigate their therapeutic potential.

The article describes facile one-pot, hi-yielding reactions to synthesize substituted 3,4-dimethyl-1H-pyrrole-2-carboxamide (3a–m) and carbohydrazide analogues (5a–l) as potential antifungal and antimicrobial agents. The structural integrity and purity of the synthesized compounds were assigned based on appropriate spectroscopic techniques. Synthesized compounds were assessed in vitro for antifungal and antibacterial activity. The compound 5h, 5i and 5j were found to be the most potent against A. fumigatus, with MIC value of 0.031 mg/mL. The compound 5f bearing a 2,6-dichloro group on the phenyl ring was found to be the most active broad spectrum antibacterial agent with MIC value of 0.039 mg/mL. The mode of action of the most promising antifungal compounds (one representative from each series; 3j and 5h) was established by their molecular docking to the active site of sterol 14α-demethylase. Molecular docking studies revealed a highly spontaneous binding ability of the tested compounds in the access channel away from catalytic heme iron of the enzyme, which suggested that the tested compounds inhibit this enzyme and would avoid heme iron related deleterious side effects observed with existing antifungal compounds.

Based on screening in computational biology and biological in vitro assays, five natural products isolated from extracts of the herbal medicine toad skin, such as cinobufagin (CBFi), bufalin (BFi), arenobufagin (ABFi), telocinobufagin (TBFi), bufotalin (BFTi), were subjected to molecular docking calculations with the use of SARS-CoV-2 main protease (PDB 6LU7 and 7BTF) and top-scoring ligand-receptor complexes were obtained. The results showed that the binding energy of ABFi to the 3CL protein was -17.044kcal/mol, which was higher than CBFi and TBFi. However, the binding energy of ABFi to the RdRp protease was -23.250 kcal/mol, which was much lower than that of CBFi and TBFi, EVEN lower than that of ABFi to the 3CL protein. ABFi also has polar interactions with amino acids such as Glu811, Ser814, Ser681 and Thr680 of RdRp enzyme. The results revealed that ABFi had a moderate inhibitory effect on the cell proliferation of SARS-CoV-2 in vitro, with an inhibition rate of 61.12%, even weaker than Remdesivir. This new discovery provides us with new ideas for in-depth studies on the development of natural products with this class of structural generalizations as inhibitors of SARS-CoV-2, and provides an experimental basis for the next step of mechanistic studies.

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.

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.

Marburg virus disease (MVD) is a hemorrhagic fever and cause death up to 88% of people and is the same as Ebola virus, It is transmitted through skin contact or mucous membrane in the eyes, nose or mouth with blood or body fluids like urine, saliva sweat, feces and by object contaminated with body fluids. Despite these facts, yet there is no approved vaccine have been developed for the eradication of Marburg virus infections. Therefore, this study described a multi epitope-based peptide vaccine against Marburg virus viral protein 35, using several immunoinformatics tools combined with molecular docking studies. Utilizing Vaxijen 2.0 server, the V35 protein revealed to be antigenic with a score of 0.4316. Prediction of the T-cell and B-cell epitopes was then conducted. RTFDAFLGV epitope was found to be the most promising one with binding affinity to the highest numbers of MHC I alleles, a positive score in the Class I immunogenicity study and non-allergen. These results were further confirmed by the good interaction of RTFDAFLGV to the groove of HLA-A*02:01 with binding energy of -8.1 kcal/mole. Finally, the vaccine was cloned in silico to ensure its validity and efficiency of expression. To ensure its safety and immunogenic profile, in-vitro and in-vivo bioassays are recommended to confirm these findings.

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

This study employs molecular docking experiments to investigate the binding interactions of Hypericin with key proteins implicated in Alzheimer's disease (AD). Hypericin demonstrated a significant binding energy of -10.1 kcal/mol with ABAD (3-hydroxyacyl-CoA dehydrogenase type II). The analysis revealed multiple binding interactions, including hydrogen bonds with arginine 6, serine 7, and glutamic acid 140, as well as hydrophobic bonds involving leucine 11, aspartic acid 142, and aspartic acid 86. These findings illuminate the intricate molecular interactions between Hypericin and ABAD, suggesting its potential as a lead compound for further therapeutic development in the context of Alzheimer's disease. Additionally, a parallel docking investigation was conducted for Beta-secretase 1 (BACE1), a protein pivotal in the generation of beta-amyloid peptides—a hallmark of Alzheimer's disease. Inhibiting BACE1 activity is considered a promising therapeutic strategy for AD. Computational analyses suggest that Hypericin may serve as a potential inhibitor of BACE1. Molecular Docking results reveal a robust binding affinity between Hypericin and BACE1, supported by a significant binding energy score of -10.3 kcal/mol. Notably, both ABAD and BACE1 exhibit high binding affinity with Hypericin. This theoretical exploration provides an initial understanding of the intricate molecular pathways involved in Alzheimer's disease and proposes Hypericin as a potential therapeutic compound for mitigating its progression.

Investigations are currently in progress to assess the potential of elevated FABP3 levels as biomarkers for assessing cardiac and cerebral injury post-stroke. Ongoing research endeavors are concentrated on elucidating the specific mechanisms by which FABPs, particularly FABP3, contribute to the complex processes involved in ischemic stroke. This study employs molecular docking techniques within a Virtual Screening framework to examine the interaction between various natural substances and the active site of human heart fatty acid-binding protein (PDB Code 4TKH). Utilizing Pyrx and Autodock Vina, the research aims to pinpoint compounds demonstrating robust binding interactions, characterized by high binding energy scores (kcal/mol) and the formation of multiple chemical bonds. The anticipated outcomes have the potential to yield valuable insights into the molecular interactions influencing human heart fatty acid-binding protein, suggesting possible implications for therapeutic interventions using these natural substances.From the docking and predicted toxicity results, Resistomycin has demonstrated superior characteristics compared to Morusin.

This dissertation presents a method for the synthesis of substituted indoles bearing heterocyclic substituents onto the 3- position. The route for the synthesis of the heterocyclic substituents indoles starts from 3-acetyl indole nucleus depend on two parts. One part: Reaction of 3-acetylindole 1 with hydrazide compounds such as phenyl hydrazine, hydrazine hydrate, and thio-semicarbazide, to yield 3-(1-(2-phenylhydrazono) ethyl)-1H-indole 2, 3-(1-hydrazonoethyl)-1H-indole 6, and thiosemicarbazone 10 respectively. 3-(1-Hydrazonoethyl)-1H-indole 6 reacted with thiophene-2-carboxaldeyde, isatin and 3-acetyl indole. In the same way, 3-(1-(2-phenylhydrazono) ethyl)-1H-indole 2 reacted with thioglycollic acid, glycine and benzaldehyde. 3-Acetylindole 1 thiosemicarbazone reacted with acetic anhydride, piperidine, concentration of hydrochloric acid and thiophene-2-carboxaldeyde. Second part. Reaction of 3-acetyl indole 1 with amines compounds such as p-nitroaniline to formed Schiff base 15 which it reacted with thioglycollic acid to give compound 16. 3-Acetyl indole 1 reacted with ethylene diamine to afford bis imine indole to afford compound 17. The reports of this docking study revealed that the new compounds exhibit good antibacterial activity. The synthesized compounds screened in vitro for their antibacterial activity revealed remarkable inhibitory effects against the selected microorganisms. And also an anti-oxidant activity studying for some synthesized compounds. Structures of the newly synthesized compounds examination by spectral data (IR, 1H- NMR and 13C -NMR).

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.

Cervical cancer is a significant global health concern, with human papillomavirus (HPV) infection being a crucial risk factor. Photodynamic therapy (PDT) is an attractive, minimally invasive treatment for HPV-related cervical lesions, which uses photosensitizers and light to selectively de-stroy abnormal cells. We aim to provide an overview of the various sorts of dyes that are utilized in PDT for decreasing the incidence and mortality of cervical cancer. Besides, the article discusses ongoing clinical trials for PDT in low-grade squamous intraepithelial neoplasia (LSIL) and high-grade squamous intraepithelial lesions (HSIL), as well as preclinical approaches for PDT in cervical cancer using various dyes. Moreover, it highlights potential dyes for PDT, examine them in in-silico condition, their pros and cons and also the solutions for enhancing their anticancer com-patibility. We also display that PDT is a promising therapeutic strategy for the diagnosis and treatment of HPV-associated cervical lesions. Moreover, it shows that using different classes of dyes improves the anticancer effects of PDT. Finally, among all dyes which are used in PDT, Fullerene demonstrated high tendency to over expressed receptors in cervical cancer cells and seems to be a proper candidate for be used in PDT more than before, but further research is necessary to assess its long-term efficacy and safety.

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.

We report a detailed investigation of vibrational modes, structure, and dynamics of glutathione (GSH) solutions using ultrasonic relaxation spectroscopy, FT-IR vibrational spectroscopy and electronic absorption measurements. The experimental data were analyzed in view of density functional theory (DFT) and molecular docking calculations. Three distinct Debye-type relaxation processes are resolved in the acoustic spectra, which are assigned to conformational changes between GSH conformers, self-association of GSH, and protonation processes. The standard volume changes for each process were estimated both experimentally and theoretically revealing a close resemblance among them. The higher the effect of the relaxation process in the structure, the greater the induced volume changes. From the temperature dependence of specific acoustic parameters, the thermodynamic characteristics of each process were determined. The experimental FT-IR spectra were compared with the corresponding theoretically predicted vibrational spectra revealing that the GSH dimers and extended conformers dominate the structure of GSH solutions in the high concentration region. The absorption spectra in the ultraviolet region confirmed the gradual aggregation mechanism that takes place in the aqueous GSH solutions. The results of the present study were discussed and analyzed in the framework of the current phenomenological status of the field.

[Co(AMPY)(DAPY)Cl2(H2O)].H2O (1) [Cu(AMPY)(DAPY)Cl2(H2O)].H2O (2) [Zn(AMPY)(DAPY)Cl2(H2O)] (3) were prepared from the ligands; AMPY = 2-amino-4-methylpyrimidine (L1), DAPY = 2,3-diaminopyridine (L2) and CoCl2.6H2O, CuCl2.2H2O and ZnCl2 in water/ethanol solutions and the three products characterized by elemental analysis, ultraviolet-visible spectroscopy (UV–Vis), Fourier-transform infrared spectroscopy (FT-IR), magnetic susceptibility, molar conductivity methods, and TGA analysis. The X-ray powder diffraction of the Co(II), Cu(II), and Zn(II) compounds showed that the geometry of monoclinic and SEM analysis revealed their morphology with a smooth surface. Molecular modeling was performed for all compounds using the density functional method DFT/B3LYP to study the structures and the frontier molecular orbitals (HOMO and LUMO). We have used Gaussian09 software for the calculations. In this study, different complexes were tested against Gram negative and Gram positive bacterial species to give insight into their broad-spectrum effects. The used pathogenic strains were two Gram positive species "Staphylococcus aureus and Micrococcus luteus" and two Gram negative species "Salmonella thyphimurium and Escherichia coli. The antifungal activity was evaluated against a pathogenic reference strain of the yeast Candida albicans. The antimicrobial and antioxidant assay results demonstrate that the tested compounds are effective against Gram positive and negative bacteria. Additionally, the compounds have an antifungal effect against Candida albicans with a maximum inhibitory zone of 2.5cm. The results demonstrated high antioxidant potential for the Zn(II) complex with a DPPH scavenging of 91.5%, however, the Cu(II) complex was low (16.5%). The data of docking with tyrosyl-tRNA synthetase presented that all compounds fit very well in the catalytic pockets of the proteins of the receptor.

As part of valorisation of agricultural waste into bioactive compounds, a series of structurally novel oleanolic acid ((3β-hydroxyolean-12-en-28-oic acid, OA-1)-phtalimidines (isoindolinones) conjugates 18a-v bearing 1,2,3-triazole moieties were designed and synthesized by treating a previously azide 4 prepared from OA-1 isolated from olive pomace (Olea europaea L.) with a wide range of propargylated phtalimidines using the Cu(I)-catalyzed click chemistry approach. OA-1 and its newly prepared analogues 18a-v were screened in vitro for their antibacterial activity against two Gram-positive bacteria, Staphylococcus aureus and Listeria monocytogenes; and two Gram-negative bacteria, Salmonella thyphimurium and Pseudomonas aeruginosa. Attractive results were obtained notably against L. monocytogenes. Compounds 18d, 18g and 18h exhibited the highest antibacterial activity when compared with OA-1 and other compounds in the series against tested pathogenic bacterial strains. Molecular docking study was performed to explore the binding mode of the most active derivatives against the target protein from L. monocytogenes. Results showed the importance of both hydrogen bonding and hydrophobic interactions with the target protein and are in favor to the experimental data.

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.

Chemotherapy represents the most applied approach to cancer treatment. Owing to the frequent onset of chemoresistance and tumor relapses, there is an urgent need to discover novel and more effective anticancer drugs. In search for therapeutic alternatives to treat the cancer disease, a se-ries of hybrid pyrazolo[3,4-d]pyrimidin-4(5H)-ones tethered with hydrazide-hydrazones, 5a-h, was synthesized from condensation reaction of pyrazolopyrimidinone-hydrazide 4 with a series of arylaldehydes in EtOH, in acid catalysis. In vitro assessment of antiproliferative effects against MCF-7 breast cancer cells, unveiled that 5a, 5e, 5g, and 5h were the most effective compounds of the series and exerted their cytotoxic activity through apoptosis induction and G0/G1 phase cell-cycle arrest. To explore their mechanism at a molecular level, 5a, 5e, 5g, and 5h were evalu-ated for their binding interactions with two well-known anticancer targets, namely the epider-mal growth factor receptor (EGFR) and the G-quadruplex DNA structures. Molecular docking simulations highlighted high binding affinity of 5a, 5e, 5g, and 5h towards EGFR. CD experi-ments suggested 5a as a stabilizer agent of the G-quadruplex from the KRAS oncogene promot-er. In the light of these findings, we propose the pyrazolo-pyrimidinone scaffold bearing a hy-drazide-hydrazone moiety as a lead skeleton for designing novel anticancer compounds.

Natural and synthetic small molecules from the NCI Developmental Therapeutics Program (DTP) were employed in molecular dynamics-based docking with DNA repair proteins whose RNA-Seq based expression was associated with overall cancer survival (OS) after adjustment for the PCNA metagene. The compounds employed were required to elicit a sensitive response (vs. resistance) in more than half of the cell lines tested for each cancer. Methodological approaches included peptide sequence alignments and homology modeling for 3D protein structure determination, ligand preparation, docking, toxicity and ADME prediction. Docking was performed for unique lists of DNA repair proteins which predict OS for AML, cancers of the breast, lung, colon, and ovaries, GBM, melanoma, and renal papillary cancer. Results indicate hundreds of drug-like and lead-like ligands with best-pose binding energies less than -6 kcal/mol. Ligand solubility for the top 20 drug-like hits approached lower bounds, while lipophilicity was acceptable. Most ligands were also blood-brain barrier permeable with high intestinal absorption rates. While the majority of ligands lacked positive prediction for Herg channel blockage and Ames carcinogenicity, there was considerable variation for predicted fathead minnow, honey bee, and Tetrahymena pyriformis toxicity. The computational results suggest the potential for new targets and mechanisms of repair inhibition and can be directly employed for in vitro and in vivo confirmatory laboratory experiments to identify new targets of therapy for cancer survival.

The aim of the study was to determine the possibility of suppression of the prostaglandin synthesis by new derivatives of N-phenylanthranilic acids; they inhibit the activity of the microsomal prostaglandin E synthase-1 (mPGES-1) enzyme using the method of a flexible molecular docking. For the docking studies the crystallographic structural models with high resolution from Protein Data Bank were used: mPGES-1 in the complex with glutathione (pdb code 4AL0). A flexible molecular docking was carried out using the Molecular Operating Environment (MOE) software package. According to the results of the docking studies four scoring functions were calculated (Affinity dG Scoring, Alpha HB Scoring, London dG Scoring, GBVI/WSA dG Scoring). The values of the scoring functions calculated indicate the thermodynamic probability and energy favorability of forming complexes between molecules of the substances under research and the specified receptor, in which arrangement of ligands in the active site of the receptor and residues of amino acids of side chains are of similar geometry and types of binding of the known inhibitors of mPGES-1 determined on the basis of the crystallographic studies.

According to our previous investigation the total methanol extract from Graptopetalum paraguayense E. Walther demonstrates a significant inhibitory effect on HSV-1. To clarify what causes this inhibitory activity on HSV-1, a metabolic profile of the plant was performed. Three main fractions: non-polar substances, polar metabolites and phenolic compounds were obtained and GC-MS analysis was carried out. Since it is well known that phenolic compounds show a significant anti-herpes effect and that viral DNA polymerase (DNApol) appears to play a key role in HSV virus replication, we present a docking and quantum-chemical analysis of the binding of these compounds to viral DNApol amino acids. Fourteen different phenolic acids found by GS/MS analyses, were used in molecular docking simulations. According to the interaction energies of all fourteen ligands in the DNApol pockets based on docking results, DFT calculations were performed on the five optimally interacting with the receptor acids. It was found that hydroxybenzoic acids from phenolic fraction of Graptopetalum paraguayense E. Walther show a good binding affinity to the amino acids from the active site of the HSV DNApol, but significantly lower than that of acyclovir. The mode of action on virus replication of acyclovir (by DNApol) is different from that of the plant phenolic acids one, probably.

Aims: Angiotensin-converting enzyme 2 (ACE2) plays an important role in the entry of coronaviruses into host cells. This paper described how carnosine, a naturally occurring supplement, can be an effective drug candidate for coronavirus disease (COVID-19) on the basis of molecular docking and modeling to host ACE2 co-crystallized with COVID-19 spike protein. Methods: First, the starting point was ACE2 inhibitors and their structure-activity relationship (SAR). Next, chemical similarity (or diversity) and PubMed searches made it possible to repurpose and assess approved or experimental drugs for COVID-19. In parallel, at all stages, authors performed bioactivity scoring to assess potential repurposed inhibitors at ACE2. Finally, investigators performed molecular docking and modeling of the identified drug candidate to host ACE2 co-crystallized with COVID-19 spike protein. Results: Carnosine emerged as the best known drug candidate to match ACE2 inhibitor structure. Preliminary docking was more optimal to ACE2 than the known typical angiotensin-converting enzyme 1 (ACE1) inhibitor (enalapril) and quite comparable to known or presumed ACE2 inhibitors. Viral spike protein elements binding to ACE2 were retained in the best carnosine pose in SwissDock at 1.75 Angstroms. Out of the three main areas of attachment expected to the co-crystallized protein structure, carnosine bind with higher affinity to two compared to the known ACE2 active site. LibDock score was 92.40 for site 3, 90.88 for site 1, and inside the active site 85.49. Conclusion: Carnosine has promising inhibitory interactions with host ACE2 co-crystallized with COVID-19 spike protein and hence could offer potential mitigating effect against current COVID-19 pandemic.

The Novel Coronavirus (COVID-19) is a positive-sense single-stranded RNA ((+)ssRNA) virus. The COVID-19 Main Proteases play very important role in the propagation of the Novel Coronavirus (COVID-19). It has already killed more than 8000 people around the world and thousands of people are getting infected every day. Therefore, it is very important to identify a potential inhibitor against COVID-19 Main Proteases to inhibit the propagation of the Novel Coronavirus (COVID-19). We have applied a drug repurposing approach of computational methodology, depending on the synergy of molecular docking and virtual screening techniques, aimed to identify possible potent inhibitors against Novel Coronavirus (COVID-19) from FDA approved antiviral compounds and from the library of active phytochemicals. On the basis of recently resolved COVID-19 Main Protease crystal structure (PDB:6LU7), the library of 100 FDA approved antiviral compounds and 1000 active components of Indian Medicinal Plants extracted for screening against COVID-19 Main Protease. The compounds were further screened using Pyrex virtual screening tool and then best inhibitors, top 19 compounds optimally docked to the COVID-19 Main Protease structure to understand the participation of specific amino acids with inhibitors at active sites. Total 19 best compounds were identified after screening based on their highest binding affinity with respect to the other screened compounds. Out of 19, 6 best compounds were further screened based on their binding affinity and best ADME properties. Nelfinavir exhibited highest binding energy -8.4 kcal/mol and strong stability with the TRP207, ILE281, LEU282, PHE3, PHE291, GLN127, ARG4, GLY283, GLU288, LYS5, LYS137, TYR126, GLY138, TYR126, SER139 and VAL135 amino acid residues of COVID-19 Main Protease participating in the interaction at the binding pocket. In addition to Nelfinavir (-8.4), Rhein (-8.1), Withanolide D (-7.8), Withaferin A (-7.7), Enoxacin (-7.4), and Aloe-emodin (-7.4) also showed good binding affinity and best ADME properties. Our findings suggest that these compounds can be used as potential inhibitors against COVID-19 Main Protease, which could be helpful in inhibiting the propagation of the Novel Coronavirus (COVID-19). Moreover, further in vitro and in vivo validation of these findings would be very helpful to bring these inhibitors to next level study.

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

Structurally related 7-acetyl-2-aryl-5-bromoindoles 2a–d and the 7-acetamido-2-aryl-5-bromoindoles 4a–d as well as their corresponding 3-trifluoroacetyl–substituted derivatives 5a–d and 5e–h were evaluated for potential antigrowth effect in vitro against the human lung cancer (A549) and cervical cancer (HeLa) cells. All of the 3-trifluoroacetyl–substituted 7-acetamido-2-aryl-5-bromoindoles 5e–h were found to be more active against both cell lines when compared to the chemotherapeutic drug, Melphalan. The most active compound 5g induced apoptosis in a caspase dependent manner for both cell lines. Compounds 5e–h were found to significantly inhibit tubulin polymerization. Molecular docking of 5g into the colchicine-binding site suggests that the compounds bind to tubulin by different type of interactions including pi-alkyl, amide-pi stacked and alkyl interactions as well as hydrogen bonding with the protein residues to elicit anticancer activity.

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

The 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-(Adamantan-1-yl)-2H-isoindole-1-carbonitrile (1) has been identified as a neurobiological fluorescent ligand that may be used to develop receptor and enzyme binding affinity assays. Compound 1 was synthesised using an optimised microwave irradiation reaction and crystallised from ethanol. Crystallization occurred in the orthorhombic space group P212121 with unit cell parameters: a = 6.4487(12) Å, b = 13.648(3) Å, c = 16.571(3) Å, V = 1458(5) Å3, Z = 4. Density functional theory (DFT) (B3LYP/6-311++G (d,p)) calculations of 1 were carried out. Results showed that the optimised geometry is similar to the crystal structure parameters with a root-mean-squared deviation of 0.143 Å. Frontier molecular orbital energies and net atomic charges are discussed with a focus on potential biological interactions. Docking experiments within the active site of the neuronal nitric oxide synthase (nNOS) protein crystal structure were carried out and analysed. Important binding interactions between the DFT optimised structure and amino acids within the nNOS active site were identified that explain the strong NOS binding affinity reported. Fluorescent properties of 1 were studied using aprotic solvents of different polarities. Compound 1 showed the highest fluorescence intensity in polar solvents with excitation and emission values of 336 nm and 380 nm, respectively.

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

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.

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

The current study describes the preparation of chitosan nanoparticles (CNPs) using hydroxychloroquine (HCQ), widely used in traditional medicine due to its diverse phar-macological and medicinal uses. This work aims to combine the HCQ drug with CS NPs to generate a novel nanocomposite with improved characteristics and bioavailability. HCQ@CS NPs is roughly shaped like roadways and has a smooth surface with an average size of 159.3±7.1 nm, a PdI of 0.224±0.101, and a zeta potential of +46.6±0.8 mV. To aid in the development of pharmaceutical systems for use in cancer therapy, the binding mech-anism and affinity of the interaction between HCQ and HCQ@CS NPs and BSA were ex-amined using stopped-flow, other spectroscopic approaches, supplemented by molecular docking analysis. HCQ and HCQ@CS NPs binding with BSA is driven by a ground-state complex formation that may be accompanied by a non-radiative energy transfer process, and binding constants indicated that HCQ@CS NPs-BSA was more stable than HCQ-BSA. The stopped-flow analysis demonstrated that, in addition to increasing BSA affinity, the nano formulation HCQ@CS NPS changes the binding process and may open up new routes for interaction. Docking experiments verified the development of the HCQ-BSA complex, with HCQ binding to the site I on the BSA structure, primarily with the amino acids Thr 578, Gln 579, Gln 525, Tyr 400, and Asn 404. Furthermore, the nano-formulation HCQ@CS NPS not only increased cytotoxicity against the A549 lung cancer cell line (IC50 = 28.57±1.72 g/ml) compared to HCQ (102.21±0.67) g/ml), but also exhibited higher anti-bacterial activity against both Gram-positive and Gram-negative bacteria when compared to HCQ and chloramphenicol which in agreement with the binding constants. The nano formulation developed in this study may offer a viable therapy option for A549 lung cancer.

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

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.

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.

Periodontitis is a disease that originates from gingival tissue, and if left untreated, can cause inflammation to spread to deeper tissues, altering bone homeostasis, and even leading to tooth loss. When the bacteria in plaque reach sufficient numbers, the immune system responds by activating immune cells that will produce cytokines, including interleukin-1 (IL-1), interleukin-6 (IL-6), and tumor necrosis factor-alpha (TNF-α). These immune signaling molecules will coordinate the inflammatory responses. IL-1 binds to the IL-1 receptor (IL-1R) on the surface of target cells, initiating an inflammatory response by activating intracellular signalling pathways. Inhibition of IL-1R can prevent the interaction of IL-1 with this receptor, thereby halting the activation of the inflammatory pathway. By substituting this connection with active anti-inflammatory substances, it is possible to prevent the contact between IL-1 and IL-1R. The licorice plant, Glycyrrhiza glabra, includes eight types of flavonoids with significant anti-inflammatory potential, namely isoliquiritigenin, glyzaglabrin, prunetin, shinpterocarpin, licochalcone A, glabridin, glisoflavone, and isoangustone A. This study aims to assess the potential of such flavonoids as candidates for periodontitis treatment through molecular docking, using AutoDock Tools v.1.5.6. All tested compounds showed highly stable bonds to IL-1R, with a binding affinity values below -7 kcal/mol. Van der Waals forces, hydrogen bonds, and hydrophobic bonds all contribute to the stability of these interactions. Among the eight tested compounds, glisoflavone, which has a Ki value of 90.89 nM and a ΔG value of -9.61 kcal/mol, has the highest binding stability to IL-1R.

New Schiff base ofN'-((E)-3,3-diphenyl-2,3-dihydro-1H-inden-1-ylidene)-2-((Z)-2-oxoindolin-3-ylidene)hydrazine-1-carbothiohydrazide (HL) and its Co(II),Ni(II), Fe(III), Cu(II), Cd(II), Zn(II) and Mn(II) complexes were created, isolated, and characterized by analytical analysis in good yields of 60 to 90%, UV–Vis, magnetic measurements, molar conductivity, and FT-IR. Additionally, mass spectrometry was used to examine the new ligand complexes' structural characteristics.1H NMR spectroscopy and elemental analysis. Depending on the electronic spectra and magnetic susceptibility measurements, these complexes are predicted to have octahedral molecular geometry. As an example, the quantum chemical properties, bond lengths, and bond angles were discussed. The structural formula for the examined ligand was improved using the Gaussian09 tool. The DFT/B3LYP method was used to calculate the energy gaps and other crucial theoretical features. The prepared hydrazine ligand complexes Co (II), Ni (II), Fe (III), Cu (II), Cd (II), Zn (II) and Mn (II) complexes were evaluated for its antimicrobial properties. As antimicrobial agents, it was found that the organic ligand was less biologically active than the targeted complexes. Gram positive bacteria (3ty7-Staphylococcus), Gram negative bacteria (3t88-Escherichia coli), and fungi (5k04-Candida albicans) all have crystal structures that interact well with the chemicals in the title, according to molecular docking.

The development of saltiness enhancement peptides is important in reducing salt intake and promoting dietary health. Our study aimed to analyze and identify straw mushroom enzymatic peptides with salty taste enhancing effects. The taste-related peptides within the straw mush-room extract were isolated via ultrafiltration and identified using UPLC-Q-TOF-MS/MS. Sub-sequent, E-tongue and sensory analyses unveiled that the ultrafiltration fraction (500-2000 Da) of straw mushroom peptides exhibited a saltiness enhancement effect. The UPLC-Q-TOF-MS/MS analysis identified 220 peptides within the ultrafiltration fractions (500-2000 Da), which were then assessed for their interaction with the T1R1/T1R3 receptor. The investigation underscored the significant involvement of Asp223, Gln243, Leu232, Asp251, and Pro254 in the binding of peptides from triple enzymatically hydrolyzed straw mushrooms to T1R1/T1R3. Based on the binding energy and active site analysis, three peptides were selected for synthesis: DFNALPFK (-9.2 kcal/mol), YNEDNGIVK (-8.8 kcal/mol), and VPGGQEIKDR (-8.9 kcal/mol). Notably, 3.2 mmol of VPGGQEIKDR enhanced a 0.05% NaCl solution to the saltiness level of a 0.15% NaCl solution. Also, 0.8 mmol of YNEDNGIVK elevated a 0.05% NaCl solution to the saltiness of a 0.1% NaCl solution.

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.

New [1,2]dithiolo[3,4-b]pyridine-5-carboxamides were synthesized through the reaction of dithiomalondianilide (N,N'-diphenyldithiomalondiamide) with 3-aryl-2-cyanoacrylamides or by a three-component reaction involving aromatic aldehydes, cyanoacetamide and dithiomalon- dianilide in the presence of morpholine. The structure of 6-amino-4-(2,4-dichlorophenyl)- 7-phenyl-3-(phenylimino)-4,7-dihydro-3H-[1,2]dithiolo[3,4-b]pyridine-5-carboxamide was con-firmed using X-ray crystallography. To understand the reaction mechanism in detail, density functional theory (DFT) calculations were performed with a Grimme B97-3c composite computa-tional scheme. The results revealed that the rate limiting step is a cyclization process leading to the closure of 1,4-dihydropyridine ring with an activation barrier of 28.8 kcal/mol. Some of dithiolo[3,4-b]pyridines exhibited moderate herbicide safening effects against 2,4-D . Additional-ly, ADMET parameters were calculated and molecular docking studies were performed to identify potential protein targets.

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.

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

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.

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

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

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

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

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.

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

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

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

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

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

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

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.

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.