Search

RecA family recombinases are the core enzymes in the process of homologous recombination, and their normal operation ensures the stability of the genome and the healthy development of organisms. The UvsX protein from bacteriophage T4 is a member of the RecA family recombinases and plays a central role in T4 phage DNA repair and replication, which provides an important model for the biochemistry and genetics of DNA metabolism. UvsX shares a high degree of structural similarity and function with RecA, which is the most deeply studied member of the RecA family. However, the detailed molecular mechanism of UvsX has not been resolved. In this study, a comprehensive all-atom molecular dynamics simulation of the UvsX protein dimer complex was carried out in order to investigate the conformational and binding properties of UvsX in combination with ATP and DNA, and the simulation of RecA was synchronized with the property comparison learning for UvsX. This study confirmed the highly conserved molecular structure characteristics and catalytic centers of RecA and UvsX, and also discovered differences in regional conformation, volatility and the ability to bind DNA between the two proteins at different temperatures, which would be helpful for the subsequent understanding and application of related recombinases.

Sialidases are increasingly used in the production of sialyloligosaccharides, a significant component of human milk oligosaccharides. Elucidating the catalytic mechanism of sialidases is critical for the rational design of better biocatalysts, thereby facilitating the industrial production of sialyloligosaccharides. Through comparative all-atom molecular dynamics simulations, we investigated the structural dynamics of sialidases in Glycoside Hydrolase family 33 (GH33). Interestingly, several sialidases displayed significant conformational transition and formed a new cleft in the simulations. The new cleft was adjacent to the innate active site of the enzyme, which serves to accommodate the glycosyl acceptor. Furthermore, the residues involved in the specific interactions with the substrate were evolutionarily conserved in the whole GH33 family, highlighting their key roles in the catalysis of GH33 sialidases. Our results enriched the catalytic mechanism of GH33 sialidases, with potential implications in the rational design of sialidases.

The Ebola virus (EBOV) is still highly infectious and causes severe hemorrhagic fevers in primates. However, there are no regulatorily approved drugs against the Ebola virus disease (EVD). The highly virulent and lethal nature of EVD highlights the need to develop therapeutic agents. Viral protein 40 kDa (VP40), the most abundantly expressed protein during infection, coordinates the assembly, budding, and release of viral particles into the host cell. It also regulates viral transcription and RNA replication. This study sought to identify small molecules that could potentially inhibit the VP40 protein by targeting the N-terminal domain using an in silico approach. The statistical quality of AutoDock Vina’s capacity to discriminate between inhibitors and decoys was determined, and an area under the curve of the receiver operating characteristic (AUC-ROC) curve of 0.791 was obtained. A total of 29,519 natural-product-derived compounds from Chinese and African sources as well as 2738 approved drugs were successfully screened against VP40. Using a threshold of −8 kcal/mol, a total of 7, 11, 163, and 30 compounds from the AfroDb, Northern African Natural Products Database (NANPDB), traditional Chinese medicine (TCM), and approved drugs libraries, respectively, were obtained after molecular docking. A biological activity prediction of the lead compounds suggested their potential antiviral properties. In addition, random-forest- and support-vector-machine-based algorithms predicted the compounds to be anti-Ebola with IC50 values in the micromolar range (less than 25 μM). A total of 42 natural-product-derived compounds were identified as potential EBOV inhibitors with desirable ADMET profiles, comprising 1, 2, and 39 compounds from NANPDB (2-hydroxyseneganolide), AfroDb (ZINC000034518176 and ZINC000095485942), and TCM, respectively. A total of 23 approved drugs, including doramectin, glecaprevir, velpatasvir, ledipasvir, avermectin B1, nafarelin acetate, danoprevir, eltrombopag, lanatoside C, and glycyrrhizin, among others, were also predicted to have potential anti-EBOV activity and can be further explored so that they may be repurposed for EVD treatment. Molecular dynamics simulations coupled with molecular mechanics Poisson–Boltzmann surface area calculations corroborated the stability and good binding affinities of the complexes (−46.97 to −118.9 kJ/mol). The potential lead compounds may have the potential to be developed as anti-EBOV drugs after experimental testing.

Fibroblast growth factor receptors (FGFRs) play critical roles in the regulation of cell growth, differentiation, and proliferation. Specifically, FGFR2 gene amplification has been implicated in gastric and breast cancer. Pan-FGFR inhibitors often cause large toxic side effects, and the highly conserved ATP-binding pocket in the FGFR1/2/3 isoforms poses an immense challenge in designing selective FGFR2 inhibitors. Recently, an indazole-based inhibitor has been discovered that can selectively target FGFR2. However, the detailed mechanism involved in selective inhibition remains to be clarified. To this end, we performed extensive molecular dynamics simulations of the apo and inhibitor-bound systems along with multiple analyses, including Markov state models, principal component analysis, a cross-correlation matrix, binding free energy calculation, and community network analysis. Our results indicated that inhibitor binding induced the phosphate-binding loop (P-loop) of FGFR2 to switch from the open to the closed conformation. This effect enhanced extensive hydrophobic FGFR2-inhibitor contacts, contributing to inhibitor selectivity. Moreover, the key conformational intermediate states, dynamics, and driving forces of this transformation were uncovered. Overall, these findings not only provided a structural basis for understanding the closed P-loop conformation for therapeutic potential but also shed light on the design of selective inhibitors for treating specific types of cancer.

The single-point mutation D26E in human β-tubulin is associated with drug resistance seen with two anti-mitotic taxanes (paclitaxel and docetaxel) when used to treat cancers. The molecular mechanism of this resistance remains elusive. However, docetaxel and a third-generation taxane, cabazitaxel, are thought to overcome this resistance. Here, structural models of both the wildtype (WT) and D26E mutant (MT) human β-tubulin were constructed based on the crystal structure of pig β-tubulin in complex with docetaxel (PDB ID: 1TUB). The three taxanes were docked into the WT and MT β-tubulin, and the resulting complexes were submitted to three independent runs of 200 ns molecular dynamic simulations, which were then averaged. MM/GBSA calculations revealed the binding energy of paclitaxel with WT and MT β-Tubulin to be −101.5 ± 8.4 and −90.4 ± 8.9 kcal/mol, respectively. The binding energy of docetaxel was estimated to be −104.7 ± 7.0 kcal/mol with the WT and −103.8 ± 5.5 kcal/mol with the MT β-tubulin. Interestingly, cabazitaxel was found to have a binding energy of −122.8 ± 10.8 kcal/mol against the WT and −106.2 ± 7.0 kcal/mol against the MT β-tubulin. These results show that paclitaxel and docetaxel bound to the MT less strongly than the WT, suggesting possible drug resistance. Similarly, cabazitaxel displayed a greater binding propensity against WT and MT β-tubulin than the other two taxanes. Furthermore, the dynamic cross-correlation matrices (DCCM) analysis suggests that the single-point mutation D26E induces a subtle dynamical difference in the ligand-binding domain. Overall, the present study revealed how the single-point mutation D26E may reduce the binding affinity of the taxanes, however, the effect of the mutation does not significantly affect the binding of cabazitaxel. © 2023 Elsevier Inc.

Graphical abstract The interleukin-1 receptor like ST2 has emerged as a potential drug discovery target since it was identified as the receptor of the novel cytokine IL-33, which is involved in many inflammatory and autoimmune diseases. For the treatment of such IL-33-related disorders, efforts have been made to discover molecules that can inhibit the protein–protein interactions (PPIs) between IL-33 and ST2, but to date no drug has been approved. Although several anti-ST2 antibodies have entered clinical trials, the exploration of small molecular inhibitors is highly sought-after because of its advantages in terms of oral bioavailability and manufacturing cost. The aim of this study was to discover ST2 receptor inhibitors based on its PPIs with IL-33 in crystal structure (PDB ID: 4KC3) using virtual screening tools with pharmacophore modeling and molecular docking. From an enormous chemical space ZINC, a potential series of compounds has been discovered with stronger binding affinities than the control compound from a previous study. Among them, four compounds strongly interacted with the key residues of the receptor and had a binding free energy

The VUV photoionisation and photofragmentation of cyclo-alanine-alanine (cAA) has been studied in a joint experimental and theoretical work. The photoelectron spectrum and the photoelectronphotoion coincidence (PEPICO) measurements, which enable control of the energy being deposited, combined with quantum chemistry calculations, provide direct insight into the cAA molecular stability after photoionisation. The analysis of the ion-neutral coincidence experiments with the molecular dynamics simulations and the exploration of the potential energy surface allows a complete identification of the fragmentation pathways. It has been found that the fragmentation always start with the ring opening through the C–C bond cleavage, followed by release of neutral moieties CO or HNCO.

Tridax procumbens is a flowering plant of the Asteraceae family with a wide range of medicinal uses like anti-inflammatory, anti-diabetic, anti-microbial, immunomodulatory, etc. This study aimed to investigate the anti-cancerous activity of human lung cancer for targeting luteolin, a phytochemical of Tridax procumbens. The computational study has been done for studying the structural properties of luteolin. The drug-likeness of the molecule has been predicted by virtual screening of ADMET properties. The molecular docking technique of the in-silico method is performed to check the complex formation between protein and ligand. The reactivity and stability of the molecule are investigated with the help of molecular dynamics (MD) simulations. In the present work, we have tried to establish a strong candidature of any of the phytochemical of Tridax Procumbens as an inhibitor against human lung cancer. Supplementary information The online version contains supplementary material available at 10.1007/s11224-022-01882-7.

The transmission and infectivity of COVID-19 have caused a pandemic that has lasted for several years. This is due to the constantly changing variants and subvariants that have evolved rapidly from SARS-CoV-2. To discover drugs with therapeutic potential for COVID-19, we focused on the 3CL protease (3CLpro) of SARS-CoV-2, which has been proven to be an important target for COVID-19 infection. Computational prediction techniques are quick and accurate enough to facilitate the discovery of drugs against the 3CLpro of SARS-CoV-2. In this paper, we used both ligand-based virtual screening and structure-based virtual screening to screen the traditional Chinese medicine small molecules that have the potential to target the 3CLpro of SARS-CoV-2. MD simulations were used to confirm these results for future in vitro testing. MCCS was then used to calculate the normalized free energy of each ligand and the residue energy contribution. As a result, we found ZINC15676170, ZINC09033700, and ZINC12530139 to be the most promising antiviral therapies against the 3CLpro of SARS-CoV-2.

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

HIV-1 maturation is the final step in the retroviral lifecycle that is regulated by the proteolytic cleavage of the Gag precursor protein. As a first-in-class HIV-1 maturation inhibitor (MI), bevirimat blocks virion maturation by disrupting capsid-spacer peptide 1 (CA-SP1) cleavage, which acts as the target of MIs. Previous alterations of beesioside I (1) produced (20S,24S)-15ꞵ,16ꞵ-diacetoxy-18,24; 20,24-diepoxy-9,19-cyclolanostane-3ꞵ,25-diol 3-O-3′,3′-dimethylsuccinate (3, DSC), showing similar anti-HIV potency compared to bevirimat. To ascertain the binding modes of this derivative, further modification of compound 1 was conducted. Three-dimensional quantitative structure–activity relationship (3D-QSAR) analysis combined with docking simulations and molecular dynamics (MD) were conducted. Five new derivatives were synthesized, among which compound 3b showed significant activity against HIV-1NL4-3 with an EC50 value of 0.28 µM. The developed 3D-QSAR model resulted in great predictive ability with training set (r2 = 0.99, q2 = 0.55). Molecular docking studies were complementary to the 3D-QSAR analysis, showing that DSC was differently bound to CA-SP1 with higher affinity than that of bevirimat. MD studies revealed that the complex of the ligand and the protein was stable, with root mean square deviation (RMSD) values

The androgen receptor (AR) is an important drug target in prostate cancer and a driver of castration-resistant prostate cancer (CRPC). A significant challenge in designing effective drugs lies in targeting constitutively active AR variants and, most importantly, nearly all AR variants lacking the ligand-binding domain (LBD). Recent findings show that an AR’s constitutive activity may occur in the presence of somatic DNA mutations within non-coding regions, but the role of these mutations remains elusive. The discovery of new drugs targeting CRPC is hampered by the limited molecular understanding of how AR binds mutated DNA sequences, frequently observed in prostate cancer, and how mutations within the protein and DNA regulate AR-DNA interactions. Using atomistic molecular dynamics (MD) simulations and quantum mechanical calculations, we focused our efforts on (i) rationalising the role of several activating DBD mutations linked to prostate cancer, and (ii) DBD interactions in the presence of abasic DNA lesions, which frequently occur in CRPC. Our results elucidate the role of mutations within DBD through their modulation of the intrinsic dynamics of the DBD-DNA ternary complex. Furthermore, our results indicate that the DNA apurinic lesions occurring in the androgen-responsive element (ARE) enhance direct AR-DNA interactions and stabilise the DBD homodimerisation interface. Moreover, our results strongly suggest that those abasic lesions may form reversible covalent crosslinks between DNA and lysine residues of an AR via a Schiff base. In addition to providing an atomistic model explaining how protein mutations within the AR DNA-binding domain affect AR dimerisation and AR-DNA interactions, our findings provide insight into how somatic mutations occurring in DNA non-coding regions may activate ARs. These mutations are frequently observed in prostate cancer and may contribute to disease progression by enhancing direct AR-DNA interactions.

Janus kinases (JAKs) are involved in numerous cellular signaling processes related to immune cell functions. JAK2 and JAK3 are associated with the pathogenesis of leukemia and common lymphoid-derived illnesses. JAK2/3 inhibitors could reduce the risk of various diseases by targeting this pathway. Herein, the naphthoquinones were experimentally and theoretically investigated to identify novel JAK2/3 inhibitors. Napabucasin and 2′-methyl napabucasin exhibited potent cell growth inhibition in TF1 (IC50 = 9.57 and 18.10 μM) and HEL (IC50 = 3.31 and 6.65 μM) erythroleukemia cell lines, and they significantly inhibited JAK2/3 kinase activity (in a nanomolar range) better than the known JAK inhibitor, tofacitinib. Flow cytometric analysis revealed that these two compounds induced apoptosis in TF1 cells in a time and dose-dependent manner. From the molecular dynamics study, both compounds formed hydrogen bonds with Y931 and L932 residues and hydrophobically contacted with the conserved hinge region, G loop, and catalytic loop of the JAK2. Our obtained results suggested that napabucasin and its methylated analog were potential candidates for further development of novel anticancer drug targeting JAKs.

Lower activity of the histaminergic system is associated with neurological disorders, including Alzheimer’s disease (AD). Thus, the enhancement of histaminergic neurotransmission by inhibition of histamine N-methyl transferase (HNMT), which degrades histamine, appears as an important approach. For this purpose, rigid and flexible molecular docking studies of 185 FDA-approved drugs with the HNMT enzyme were carried out to select two compounds to perform molecular dynamics (MD) simulations to evaluate the binding free energies and stability of the enzyme–drug complexes. Finally, an HNMT inhibition assay was performed to corroborate their effect towards HNMT. Molecular docking studies with HNMT allowed the selection of dihydroergotamine and vilazodone since these molecules showed the lowest Gibbs free energy values. Analysis of the binding mode of vilazodone showed interactions with the binding pocket of HNMT with Glu28, Gln143, and Asn283. In contrast, dihydroergotamine binds to the HNMT active site in a different location, apparently because it is overall the more rigid ligand compared to flexible vilazodone. HNMT inhibitory activity for dihydroergotamine and vilazodone was corroborated (IC50 = 72.89 μM and 45.01 μM, respectively) by in vitro assays. Drug repurposing of HNMT was achieved by employing computational studies.

Ewing sarcoma (ES) is a highly malignant carcinoma prevalent in children and most frequent in the second decade of life. It mostly occurs due to t(11;22) (q24;q12) translocation. This translocation encodes the oncogenic fusion protein EWS/FLI (Friend leukemia integration 1 transcription factor), which acts as an aberrant transcription factor to deregulate target genes essential for cancer. Traditionally, flavonoids from plants have been investigated against viral and cancerous diseases and have shown some promising results to combat these disorders. In the current study, representative flavonoid compounds from various subclasses are selected and used to disrupt the RNA-binding motif of EWS, which is required for EWS/FLI fusion. By blocking the RNA-binding motif of EWS, it might be possible to combat ES. Therefore, molecular docking experiments validated the binding interaction patterns and structural behaviors of screened flavonoid compounds within the active region of the Ewing sarcoma protein (EWS). Furthermore, pharmacogenomics analysis was used to investigate potential drug interactions with Ewing sarcoma-associated genes. Finally, molecular dynamics simulations were used to investigate the stability of the best selected docked complexes. Taken together, daidzein, kaempferol, and genistein exhibited a result comparable to ifosfamide in the proposed in silico study and can be further analyzed as possible candidate compounds in biological in vitro studies against ES.

Plastic pollution has become a global environmental threat, which leads to an increasing concern over the consequences of plastic exposition on global health. Plastic nanoparticles have been shown to influence the folding of proteins and influence the formation of aberrant amyloid proteins, therefore potentially triggering the development of systemic and local amyloidosis. This work aims to study the interaction between nanoplastics and β-amyloid fibrils to better understand the potential role of nanoplastics in the outbreak of neurodegenerative disorders. Using microsecond-long coarse-grained molecular dynamics simulations, we investigated the interactions between neutral and charged nanoparticles made of the most common plastic materials (i.e., polyethylene, polypropylene, and polystyrene) and β-amyloid fibrils. We observe that the occurrence of contacts, region of amyloid fibril involved, and specific amino acids mediating the interaction depend on the type and charge of the nanoparticles.

The polyextremophilic β-galactosidase enzyme of the haloarchaeon Halorubrum lacusprofundi functions in extremely cold and hypersaline conditions. To better understand the basis of polyextremophilic activity, the enzyme was studied using steady-state kinetics and molecular dynamics at temperatures ranging from 10 °C to 50 °C and salt concentrations from 1 M to 4 M KCl. Kinetic analysis showed that while catalytic efficiency (kcat/Km) improves with increasing temperature and salinity, Km is reduced with decreasing temperatures and increasing salinity, consistent with improved substrate binding at low temperatures. In contrast, kcat was similar from 2–4 M KCl across the temperature range, with the calculated enthalpic and entropic components indicating a threshold of 2 M KCl to lower the activation barrier for catalysis. With molecular dynamics simulations, the increase in per-residue root-mean-square fluctuation (RMSF) was observed with higher temperature and salinity, with trends like those seen with the catalytic efficiency, consistent with the enzyme’s function being related to its flexibility. Domain A had the smallest change in flexibility across the conditions tested, suggesting the adaptation to extreme conditions occurs via regions distant to the active site and surface accessible residues. Increased flexibility was most apparent in the distal active sites, indicating their importance in conferring salinity and temperature-dependent effects.

Fatty acid amide hydrolase (FAAH) plays a key role in the control of cannabinoid signaling and it represents a promising therapeutic strategy for the treatment of a wide range of diseases, including neuropathic pain and chronic inflammation. Starting from kinetics experiments carried out in our previous work for the most potent inhibitor 2-amino-3-chloropyridine amide (TPA14), we have investigated its non-competitive mechanism of action using molecular dynamics, thermodynamic integration and QM-MM/GBSA calculations. The computational studies highlighted the impact of mutations on the receptor binding pockets and elucidated the molecular basis of the non-competitive inhibition mechanism of TPA14, which prevents the endocannabinoid anandamide (AEA) from reaching its pro-active conformation. Our study provides a rationale for the design of non-competitive potent FAAH inhibitors for the treatment of neuropathic pain and chronic inflammation.

The inelastic scattering of N2 molecules from silica surfaces, taken at 100 K, has been investigated by adopting a semiclassical collision model in conjunction with the appropriate treatment of the long-range interaction forces. Such forces promote the formation of the precursor state that controls all basic elementary processes occurring at the gas–surface interphase. The probabilities for the different elementary surface processes triggered by quartz are determined and compared with those recently obtained for another silica polymorph (cristobalite). In addition, the final roto-vibrational distributions of N2 molecules undergoing inelastic scattering have been characterized. N2 molecules, impinging on both considered surfaces in low-medium vibrational states, preserve the initial vibrational state, while those inelastically scattered are rotationally excited and translationally colder. The surface temperature effect, investigated by raising the temperature itself from 100 K up to 1000 K, emerges more sharply for the cristobalite polymorph, mainly for the molecules impinging in the ground roto-vibrational state and with low collision energies.

The natural product Salvinorin A (SalA) was the first nitrogen-lacking agonist discovered for the opioid receptors and exhibits high selectivity for the kappa opioid receptor (KOR) turning SalA into a promising analgesic to overcome the current opioid crisis. Since SalA’s suffers from poor pharmacokinetic properties, particularly the absence of gastrointestinal bioavailability, fast metabolic inactivation and subsequent short duration of action, the rational design of new tailored analogs with improved clinical usability is highly desired. Despite being known for decades, the binding mode of SalA within the KOR remains elusive as several conflicting binding modes of SalA were proposed hindering the rational design of new analgesics. In this study, we rationally determined the binding mode of SalA to the active state KOR by in silico experiments (docking, molecular dynamics simulations, dynophores) in the context of all available mutagenesis studies and structure-activity relationship (SAR) data. To the best of our knowledge this is the first comprehensive evaluation of SalA’s binding mode since the determination of the active state KOR crystal structure. SalA binds above the morphinan binding site with its furan pointing towards the intracellular core while the C2-acetoxy group is oriented towards the extracellular loop 2 (ECL2). SalA is solely stabilized within the binding pocket by hydrogen bonds (C210ECL2, Y3127.35, Y3137.36) and hydrophobic contacts (V1182.63, I1393.33, I2946.55, I3167.39). With the disruption of this interaction pattern or the establishment of additional interactions within the binding site we were able to rationalize the experimental data for selected analogs. We surmise the C2-substituent interactions as important for SalA and its analogs to be experimentally active, albeit moderate frequency within MD simulations of SalA. We further identified the non-conserved residues 2.63, 7.35, and 7.36 responsible for the KOR subtype selectivity of SalA. We are confident that the elucidation of the SalA binding mode will promote the understanding of KOR activation and the facilitate the development of novel analgesics that are urgently needed.

Previous research has indicated that various metabolites belonging to phenolic acids (PAs), produced by gut microflora through the breakdown of polyphenols, help in promoting bone development and protecting bone from degeneration. Results have also suggested that G-protein-coupled receptor 109A (GPR109A) functions as a receptor for those specific PAs such as hippuric acid (HA) and 3-(3-hydroxyphenyl) propionic acid (3-3-PPA). Indeed, HA has a molecular structural similarity with nicotinic acid (niacin) which has been shown previously to bind to GPR109A receptor and to mediate antilipolytic effects; however, the binding pocket and the structural nature of the interaction remain to be recognized. In the present study, we employed a computational strategy to elucidate the molecular structural determinants of HA binding to GPR109A and GPR109B homology models in understanding the regulation of osteoclastogenesis. Based on the docking and molecular dynamics simulation studies, HA binds to GPR109A similarly to niacin. Specifically, the transmembrane helices 3, 4 and 6 (TMH3, TMH4 and TMH6) and Extracellular loop 1 and 2 (ECL1 and ECL2) residues of GRP109A; R111 (TMH3), K166 (TMH4), ECL2 residues; S178 and S179, and R251 (TMH6), and residues of GPR109B; Y87, Y86, S91 (ECL1) and C177 (ECL2) contribute for HA binding. Simulations and Molecular Mechanics Poisson-Boltzmann solvent accessible area (MM-PBSA) calculations reveal that HA has higher affinity for GPR109A than for GPR109B. Additionally, in silico mutation analysis of key residues have disrupted the binding and HA exited out from the GPR109A protein. Furthermore, measurements of time-resolved circular dichroism spectra revealed that there are no major conformational changes in the protein secondary structure on HA binding. Taken together, our findings suggest a mechanism of interaction of HA with both GPR109A and GPR109B receptors.

Lysine-specific demethylase 1 (LSD1) is a histone-modifying enzyme, which is a significant target for anticancer drug research. In this work, 40 reported tetrahydroquinoline-derivative inhibitors targeting LSD1 were studied to establish the three-dimensional quantitative structure–activity relationship (3D-QSAR). The established models CoMFA (Comparative Molecular Field Analysis (q2 = 0.778, Rpred2 = 0.709)) and CoMSIA (Comparative Molecular Similarity Index Analysis (q2 = 0.764, Rpred2 = 0.713)) yielded good statistical and predictive properties. Based on the corresponding contour maps, seven novel tetrahydroquinoline derivatives were designed. For more information, three of the compounds (D1, D4, and Z17) and the template molecule 18x were explored with molecular dynamics simulations, binding free energy calculations by MM/PBSA method as well as the ADME (absorption, distribution, metabolism, and excretion) prediction. The results suggested that D1, D4, and Z17 performed better than template molecule 18x due to the introduction of the amino and hydrophobic groups, especially for the D1 and D4, which will provide guidance for the design of LSD1 inhibitors.

The interactions of lipid molecules with various solvent molecules is of utmost importance in the formulation of various drug delivery and personal care formulations. In this manuscript, a series of all-atom molecular dynamics simulations were used to investigate how the structural and interfacial properties of a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) monolayer change when interacting with a range of diols that have varying carbon chain lengths and patterns of hydroxylation. In comparison to water, we find that all of the diols studied result in a more disordered and thinner monolayer. Additionally, we find that the shorter diols with the hydroxyl groups on neighbouring carbons (1,2-ethanediol and 1,2-propanediol) are able to penetrate deeper into the head group region of the lipid monolayers and as a result significantly disorder and thin the monolayers. Like water, we find that the diols also form hydrogen-bonded networks that connect the DMPC head groups in neighbouring molecules. Interestingly, we find that the number of butanediol molecules that form these solvent-mediated interactions between the DMPC head groups is directly affected by the distribution of the hydroxyl groups within the diol molecules. The results presented here provide a mechanistic description of how the chemistry of diol solvent molecules will affect the structural and interfacial properties of lipid structures in solution.

COVID-19 can cause different neurological symptoms in some people, including smell, inability to taste, dizziness, confusion, delirium, seizures, stroke, etc. Owing to the issue of vaccine effectiveness, update and coverage, we still need one or more diversified strategies as the backstop to manage illness. Characterizing the structural basis of ligand recognition in the main protease (Mpro) of SARS-CoV-2 will facilitate its rational design and development of potential drug candidates with high affinity and selectivity against COVID-19. Up to date, covalent-, non-covalent inhibitors and allosteric modulators have been reported to bind to different active sites of Mpro. In the present work, we applied the molecular dynamics (MD) simulations to systematically characterize the potential binding features of catalytic active site and allosteric binding sites in Mpro using a dataset of 163 3D structures of Mpro-inhibitor complexes, in which our results are consistent with the current studies. In addition, umbrella sampling (US) simulations were used to explore the dissociation processes of substrate pathway and allosteric pathway. All the information provided new insights into the protein features of Mpro and will facilitate its rational drug design for COVID-19.

Xyloglucan endotransglycosylases (XETs) play key roles in the remodelling and reconstruction of plant cell walls. These enzymes catalyse homo-transglycosylation reactions with xyloglucan-derived donor and acceptor substrates and hetero-transglycosylation reactions with a variety of structurally diverse polysaccharides. In this work, we describe the basis of acceptor substrate binding specificity in non-specific Tropaeolum majus (TmXET6.3) and specific Populus tremula x tremuloides (PttXET16A) XETs, using molecular docking and molecular dynamics (MD) simulations combined with binding free energy calculations. The data indicate that the enzyme-donor (xyloglucan heptaoligosaccharide or XG-OS7)/acceptor complexes with the linear acceptors, where a backbone consisted of glucose (Glc) moieties linked via (1,4)- or (1,3)-β-glycosidic linkages, were bound stably in the active sites of TmXET6.3 and PttXET16A. Conversely, the acceptors with the (1,6)-β-linked Glc moieties were bound stably in TmXET6.3 but not in PttXET16A. When in the (1,4)-β-linked Glc containing acceptors, the saccharide moieties were replaced with mannose or xylose, they bound stably in TmXET6.3 but lacked stability in PttXET16A. MD simulations of the XET-donor/acceptor complexes with acceptors derived from (1,4;1,3)-β-glucans highlighted the importance of (1,3)-β-glycosidic linkages and side chain positions in the acceptor substrates. Our findings explain the differences in acceptor binding specificity between non-specific and specific XETs and associate theoretical to experimental data.

In the present study, the anti-diabetic potential of Ocimum tenuiflorum was investigated using computational techniques for α-glucosidase, α-amylase, aldose reductase, and glycation at multiple stages. It aimed to elucidate the mechanism by which phytocompounds of O. tenuiflorum treat diabetes mellitus using concepts of druglikeness and pharmacokinetics, molecular docking simulations, molecular dynamics simulations, and binding free energy studies. Isoeugenol is a phenylpropene, propenyl-substituted guaiacol found in the essential oils of plants. During molecular docking modelling, isoeugenol was found to inhibit all the target enzymes, with a higher binding efficiency than standard drugs. Furthermore, molecular dynamic experiments revealed that isoeugenol was more stable in the binding pockets than the standard drugs used. Since our aim was to discover a single lead molecule with a higher binding efficiency and stability, isoeugenol was selected. In this context, our study stands in contrast to other computational studies that report on more than one compound, making it difficult to offer further analyses. To summarize, we recommend isoeugenol as a potential widely employed lead inhibitor of α-glucosidase, α-amylase, aldose reductase, and glycation based on the results of our in silico studies, therefore revealing a novel phytocompound for the effective treatment of hyperglycemia and diabetes mellitus.

The recent successes of AlphaFold and RoseTTAFold have demonstrated the value of AI methods in highly accurate protein structure prediction. Despite these advances, the role of these methods in the context of small-molecule drug discovery still needs to be thoroughly explored. In this study, we evaluated whether the AI-based models can reliably reproduce the three-dimensional structures of protein–ligand complexes. The structure we chose was NLRP3, a challenging protein target in terms of obtaining a three-dimensional model both experimentally and computationally. The conformation of the binding pockets generated by the AI models was carefully characterized and compared with experimental structures. Further molecular docking results indicated that AI-predicted protein structures combined with molecular dynamics simulations offers a promising approach in small-molecule drug discovery.

Eight gallic acid alkyl esters (1–8) were synthesized via Fischer esterification and evaluated for their trypanocidal and leishmanicidal activity using bloodstream forms of Trypanosoma brucei and promastigotes of Leishmania major. The general cytotoxicity of the esters was evaluated with human HL-60 cells. The compounds displayed moderate to good trypanocidal but zero to low leishmanicidal activity. Gallic acid esters with alkyl chains of three or four carbon atoms in linear arrangement (propyl (4), butyl (5), and isopentyl (6)) were found to be the most trypanocidal compounds with 50% growth inhibition values of ~3 μM. On the other hand, HL-60 cells were less susceptible to the compounds, thus, resulting in moderate selectivity indices (ratio of cytotoxic to trypanocidal activity) of >20 for the esters 4–6. Modeling studies combining molecular docking and molecular dynamics simulations suggest that the trypanocidal mechanism of action of gallic acid alkyl esters could be related to the inhibition of the T. brucei alternative oxidase. This suggestion is supported by the observation that trypanosomes became immobile within minutes when incubated with the esters in the presence of glycerol as the sole substrate. These results indicate that gallic acid alkyl esters are interesting compounds to be considered for further antitrypanosomal drug development.

In this paper, Amomi fructus (Latin) was used to explore the mechanism of ionic liquids (ILs) in improving the extraction efficiency of essential oils. Microwave assisted ionic liquid treatment followed by a hydro-distillation (MILT-HD) process for isolating Amomi fructus essential oil was optimized by multi-objective optimization. Under optimum operating conditions, the IL-assisted extraction method not only enhances extraction efficiency but also reduces energy demands and CO2 emissions. Since the hydrogen bond structure network of cellulose in the cell wall is an important reason for hindering diffusion of essential oils, the mechanism of ILs was explored by density functional theoretical (DFT) and molecular dynamics (MD) simulations. According to DFT calculations, ILs can facilitate the cleavage of cellulose chains and have strong non-covalent interactions with cellulose. Based on the MD simulations, the degree of destruction of the cellulose hydrogen bond structure was explored. According to the DFT and MD simulations, the ILs can significantly destroy cellulose structure, thereby promoting essential oil release from the plant. These results were confirmed by scanning electron microscopy (SEM) and Fourier transform infrared spectroscopy (FTIR). This work is conducive to better understand the MILT-HD process for isolating essential oil and comprehensively understand the mechanism of ILs in the extraction process.

Human interferon-gamma (hIFNγ) is a crucial signaling molecule with an important role in the initialization and development of the immune response of the host. However, its aberrant activity is also associated with the progression of a multitude of autoimmune and other diseases, which determines the need for effective inhibitors of its activity. The development of such treatments requires proper understanding of the interaction of hIFNγ to its cell-surface receptor hIFNGR1. Currently, there is no comprehensive model of the mechanism of this binding process. Here, we employ molecular dynamics simulations to study on a microscopic level the process of hIFNγ–hIFNGR1 complex formation in different scenarios. We find that the two molecules alone fail to form a stable complex, but the presence of heparan-sulfate-like oligosaccharides largely facilitates the process by both demobilizing the highly flexible C-termini of the cytokine and assisting in the proper positioning of its globule between the receptor subunits. An antiproliferative-activity assay on cells depleted from cell-surface heparan sulfate (HS) sulfation together with the phosphorylation levels of the signal transducer and activator of transcription STAT1 confirms qualitatively the simulation-based multistage complex-formation model. Our results reveal the key role of HS and its proteoglycans in all processes involving hIFNγ signalling.

The interaction between coiled carbon nanotubes (CCNT) and the polymer matrix is important in the mechanical, thermal, and electrical properties of the CCNT reinforced nanocomposite. In this study, molecular dynamics (MD) simulations were performed to study the interfacial characteristics of polymer nanocomposites (PNCs). Furthermore, the influence of the geometries of the CCNTs on the load transfer mechanism is evaluated. Pullout simulations considering different geometries of CCNTs are carried out to examine the tensile force and the interfacial shear stress (ISS). The results reveal that the maximal tensile force is reduced by increasing CCNT inner diameters, increasing the helix angles, and decreasing nanotube diameters. The distance between CCNTs and the polymer matrix is varied, and the interfacial distance favors greater ISS. Decreasing the inner diameter of the CCNT, the helix angle, and the tube diameter increases the ISS. The enhancement mechanism of CCNT/polymer composites has also been illustrated. Due to a lack of experimental results, only numerical results are given. The present study helps to understand the interfacial adhesion behavior between the polymer matrix and CCNTs and is expected to contribute to the development of CCNT reinforced polymer composites.

Influenza virus infections continue to be a significant and recurrent public health problem. Although vaccine efficacy varies, regular immunisation is the most effective method for suppressing the influenza virus. Antiviral drugs are available for influenza, although two of the four FDA-approved antiviral treatments have resulted in significant drug resistance. Therefore, new treatments are being sought to reduce the burden of flu-related illness. The time-consuming development of treatments for new and re-emerging diseases such as influenza and the high failure rate are increasing concerns. In this context, we used an in silico-based drug repurposing method to repurpose FDA-approved drugs as potential therapies against the H7N9 virus. To find potential inhibitors, a total of 2568 drugs were screened. Promacta, tucatinib, and lurasidone were identified as promising hits in the DrugBank database. According to the calculations of MM-GBSA, tucatinib (−54.11 kcal/mol) and Promacta (−56.20 kcal/mol) occupied the active site of neuraminidase with a higher binding affinity than the standard drug peramivir (−49.09 kcal/mol). Molecular dynamics (MD) simulation studies showed that the C-α atom backbones of the complexes of tucatinib and Promacta neuraminidase were stable throughout the simulation period. According to ADME analysis, the hit compounds have a high gastrointestinal absorption (GI) and do not exhibit properties that allow them to cross the blood–brain barrier (BBB). According to the in silico toxicity prediction, Promacta is not cardiotoxic, while lurasidone and tucatinib show only weak inhibition. Therefore, we propose to test these compounds experimentally against the influenza H7N9 virus. The investigation and validation of these potential H7N9 inhibitors would be beneficial in order to bring these compounds into clinical settings.

© 2022 Elsevier Inc.The connection of Epstein Barr virus (EBV) with diseases such as Burkitt Lymphoma, Hodgkin disease, multiple sclerosis, systemic lupus erythematosus and various B-cell lymphomas made EBV glycoproteins one of the most popular vaccine immunogens. As a protein being encoded by EBV, the viral membrane envelope protein gp350 is studied extensively due to its abundancy on the surface and its interaction with complementary receptor, CR2. The binding of CR2 and gp350 not only leads to the entrance of the virus to the B-cells, but also prevents CR2 and C3d protein interactions that are required for immune response. Thus, understanding the inhibition of gp350 activity is crucial for vaccine development. Although, the active residues on gp350 structure were determined by several mutational studies, the exact mechanism of CR2 binding is still not clear. To this end, we have performed molecular docking followed by molecular dynamics simulations and MM-PBSA on wildtype and several mutated gp350 and CR2 structures. Apart from identifying crucial amino acids, the results of per-residue decomposition energy analysis clarified the individual energy contributions of amino acids and were also found to be accurate in differentiating the active site residues in CR2 binding. Here, we highlight the role of binding region residues (linker-1) but more interestingly, the dynamic relation between the distant sites of gp350 (linker-2 and D3 residues) and CR2. These findings can lead further vaccine development strategies by pointing to the importance of computationally found novel regions that can be potentially used to modulate gp350 activity.

Chia seed peptides (CSP) can be a source of multifunctional biopeptides to treat non-communicable diseases. However, interactions and binding affinity involved in targeting specific receptors remains unexplored. In this study, molecular simulation techniques were used as virtual screening of CSP to determine drug-like candidates using a multi-target-directed ligand approach. CSP fraction with the best bioactivities in vitro was sequenced. Then, a prediction model was built using physicochemical descriptors (hydrophobicity, hydrophilicity, intestinal stability, antiangiogenic, antihypertensive, and anti-inflammatory) to calculate potential scores and rank possible biopeptides. Furthermore, molecular dynamics simulations (MDS) and ensemble molecular docking analysis were carried out using four human protein targets (ACE, angiotensin converting enzyme; VEGF, vascular endothelial growth factor; GLUC, glucocorticoid and MINC, mineralocorticoid receptors). Five known-sequence peptides (NNVFYPF, FNIVFPG, SRPWPIDY, QLQRWFR, GSRFDWTR) and five de novo peptides (DFKF, DLRF, FKAF, FRSF, QFRF) had the lowest energy score and higher affinity for ACE and VEGF. The therapeutic effects of these selected peptides can be related to the inhibition of the enzymes involved in angiogenesis and hypertension, due to formation of stable complexes with VEGF and ACE binding sites, respectively. The application of MDS is a good resource for identifying bioactive peptides for future experimental validation.

International audience; Properties of polyhydroxycompounds, such as sugars and polyols, which are used to stabilize various pharmaceuticals and biological objects during freezing, desiccation, and storage, depend on water content. In this investigation, glycerol-water mixtures containing 40 to 0 wt% water are studied by FTIR spectroscopy and molecular dynamics simulation. The HOH bending (approx. 1640 cm−1) band in the FTIR spectra is essentially independent of water content in the range of 0 to approx. 12 wt% water, whereas a red shift to lower wavenumbers is observed as water content increases beyond 12 wt%. Water clustering patterns, which obtained from MD simulations, show a gradual decrease of unclustered water molecules at a higher water concentration, with the corresponding increase in small (3–4 water molecules) water clusters. Large water clusters (>30 molecules) appear at 9 to 12 wt% water; at 25 wt% water, smaller water clusters disappear, with essentially all water molecules belonging to continuous water domain. The water clustering and hydrogen bonding patterns are discussed in view of fundamental physical chemical processes including proton transfer, water catalysis of chemical reactions, and freezing behavior.

The behavior against temperature and thermal stability of enzymes is a topic of importance for industrial biocatalysis. This study focuses on the kinetics and thermodynamics of the thermal inactivation of Lipase PS from B. cepacia and Palatase from R. miehei. Thermal inactivation was investigated using eight inactivation models at a temperature range of 40–70 °C. Kinetic modeling showed that the first-order model and Weibull distribution were the best equations to describe the residual activity of Lipase PS and Palatase, respectively. The results obtained from the kinetic parameters, decimal reduction time (D and tR), and temperature required (z and z’) indicated a higher thermal stability of Lipase PS compared to Palatase. The activation energy values (Ea) also indicated that higher energy was required to denature bacterial (34.8 kJ mol−1) than fungal (23.3 kJ mol−1) lipase. The thermodynamic inactivation parameters, Gibbs free energy (ΔG#), entropy (ΔS#), and enthalpy (ΔH#) were also determined. The results showed a ΔG# for Palatase (86.0–92.1 kJ mol−1) lower than for Lipase PS (98.6–104.9 kJ mol−1), and a negative entropic and positive enthalpic contribution for both lipases. A comparative molecular dynamics simulation and structural analysis at 40 °C and 70 °C were also performed.

Uperin 3.5 is a remarkable natural peptide obtained from the skin of toadlets comprised of 17 amino acids which exhibits both antimicrobial and amyloidogenic properties. Molecular dynamics simulations were performed to study the β-aggregation process of uperin 3.5 as well as two of its mutants, in which the positively charged residues Arg7 and Lys8 have been replaced by alanine. All three peptides rapidly underwent spontaneous aggregation and conformational transition from random coils to beta-rich structures. The simulations reveal that the initial and essential step of the aggregation process involves peptide dimerization and the formation of small beta-sheets. A decrease in positive charge and an increase in the number of hydrophobic residues in the mutant peptides lead to an increase in the rate of their aggregation.

Classical molecular dynamics simulations of aqueous methanesulfonic acid solutions have been con- ducted at room temperature in the entire composition range. The dissociation of the acid has been con- sidered according to the available experimental data. Then, the systems are constituted by the following molecular species: water and methanesulfonic acid molecules, hydronium cations, and mesylate anions. The simulations have been carried out employing a reliable force field, which provides density values that show an excellent agreement with the available experimental data. The shear viscosity, the diffusion coefficients of the molecular species, and the radial distribution functions, which involve water mole- cules, have also been computed. We have observed that water molecules diffuse faster than the other species at all concentrations. Moreover, the shear viscosity and the diffusion coefficients exhibit a notice- able unimodal concentration dependence with extrema located at 90 wt% (0.628 mol fraction). A detailed hydrogen bond analysis, concerning water molecules, has also been made. At low dilutions, water exhi- bits its well-known hydrogen bond tetrahedral structure, which vanishes as the acid concentration increases. The most labile water molecules are those bonded to the anions. At large concentrations, the lability of the water molecules, bonded to other water molecules, increases, and that of the water molecules, bonded to the cations, decreases. The interrupted lifetimes also show a unimodal dependence with maxima located at 90 wt%. This behavior could be related to the appearance of the methanesulfonic acid hydrogen bond network, which emerges at large weight percentage values because the acid almost completely dissociates up to concentrations of 70 wt% (0.304 mole fraction).

Molecular dynamics simulations employing the all-atom optimized potential for liquid simulations (OPLS-AA) force field were performed for determining self-diffusion coefficients (D11) of ethanol and tracer diffusion coefficients (D12) of solutes in ethanol at several temperature and pressure conditions. For simulations employing the original OPLS-AA diameter of ethanol’s oxygen atom (σOH), calculated and experimental diffusivities of protic solutes differed by more than 25%. To correct this behavior, the σOH was reoptimized using the experimental D12 of quercetin and of gallic acid in liquid ethanol as benchmarks. A substantial improvement of the calculated diffusivities was found by changing σOH from its original value (0.312 nm) to 0.306 nm, with average absolute relative deviations (AARD) of 3.71% and 4.59% for quercetin and gallic acid, respectively. The new σOH value was further tested by computing D12 of ibuprofen and butan-1-ol in liquid ethanol with AARDs of 1.55% and 4.81%, respectively. A significant improvement was also obtained for the D11 of ethanol with AARD = 3.51%. It was also demonstrated that in the case of diffusion coefficients of non-polar solutes in ethanol, the original σOH=0.312 nm should be used for better agreement with experiment. If equilibrium properties such as enthalpy of vaporization and density are estimated, the original diameter should be once again adopted.

The Rhodococcus erythropolis (strain IGTS8) bacterium has a tremendous industrial interest as it can remove sulfur from crude oil through its four-enzyme (DszA-D) 4S metabolic pathway. DszC is one of the rate-limiting enzymes of the pathway and the one that most suffers from feedback inhibition. We have combined molecular docking and molecular dynamics simulations to identify binding sites through which two products of the 4S pathway, 2-hydroxybiphenyl and 2′-hydroxybiphenyl-2-sulfinate, induce DszC feedback inhibition. We have identified four potential binding sites: two adjacent binding sites close to the 280–295 lid loop proposed to contribute to DszC oligomerization and proper binding of the flavin mononucleotide cofactor, and two other close to the active site of DszC and the substrate binding site. By considering (i) the occupancy of the binding sites and (ii) the similar inhibitor poses, we propose that the mechanism of feedback inhibition of DszC occurs through disturbance of the DszC oligomerization and consequent binding of the flavin mononucleotide due to the weakening of the interactions between the 280–295 lid loop, and both the 131–142 loop and the C-terminal tail. Nevertheless, inhibitor binding close to the active site or the substrate binding sites also compromises critical interactions within the active site of DszC. The disclosed molecular details provide valuable insight for future rational enzyme engineering protocols to develop DszC mutants more resistant against the observed feedback inhibition mechanism.

Protein–protein interactions (PPI) represent attractive targets for drug design. Thus, aiming at a deeper insight into the HSV-1 envelope glycoprotein D (gD), protein–protein docking and dynamic simulations of gD-HVEM and gD-Nectin-1 complexes were performed. The most stable complexes and the pivotal key residues useful for gD to anchor human receptors were identified and used as starting points for a structure-based virtual screening on a library of both synthetic and designed 1,2,3-triazole-based compounds. Their binding properties versus gD interface with HVEM and Nectin-1 along with their structure-activity relationships (SARs) were evaluated. Four [1,2,3]triazolo[4,5-b]pyridines were identified as potential HSV-1 gD inhibitors, for their good theoretical affinity towards all conformations of HSV-1 gD. Overall, this study suggests promising basis for the design of new antiviral agents targeting gD as a valuable strategy to prevent viral attachment and penetration into the host cell.

The global pandemic of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become an urgent public health problem. Spike (S) protein mediates the fusion between the virus and the host cell membranes, consequently emerging as an important target of drug design. The lack of comparisons of in situ full-length S homotrimer structures in different states hinders understanding the structures and revealing the function, thereby limiting the discovery and development of therapeutic agents. Here, the steady-state structures of the in situ full-length S trimer in closed and open states (Sclosed and Sopen) were modeled with the constraints of density maps, associated with the analysis of the dynamic structural differences. Subsequently, we identified various regions with structure and property differences as potential binding pockets for ligands that promote the formation of inactive trimeric protein complexes. By using virtual screening strategy and a newly defined druggable cavity, five ligands were screened with potential bioactivities. Then molecular dynamic (MD) simulations were performed on apo protein structures and ligand bound complexes to reveal the conformational changes upon ligand binding. Our simulation results revealed that sulforaphane (SFN), which has the best binding affinity, could inhibit the conformational changes of S homotrimer that would occur during the viral membrane fusion. Our results could aid in the understanding of the regulation mechanism of S trimer aggregation and the structure-activity relationship, facilitating the development of potential antiviral agents.

Human lipoxygenase 12 (hALOX12) catalyzes the conversion of docosahexaenoic acid (DHA) into mainly 14S-hydroperoxy-4Z,7Z,10Z,12E,16Z,19Z-docosahexaenoic acid (14S-H(p)DHA). This hydroperoxidation reaction is followed by an epoxidation and hydrolysis process that finally leads to maresin 1 (MaR1), a potent bioactive specialized pro-resolving mediator (SPM) in chronic inflammation resolution. By combining docking, molecular dynamics simulations, and quantum mechanics/molecular mechanics calculations, we have computed the potential energy profile of DHA hydroperoxidation in the active site of hALOX12. Our results describe the structural evolution of the molecular system at each step of this catalytic reaction pathway. Noteworthy, the required stereospecificity of the reaction leading to MaR1 is explained by the configurations adopted by DHA bound to hALOX12, along with the stereochemistry of the pentadienyl radical formed after the first step of the mechanism. In pig lipoxygenase 15 (pigALOX15-mini-LOX), our calculations suggest that 14S-H(p)DHA can be formed, but with a stereochemistry that is inadequate for MaR1 biosynthesis.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for COVID-19, which was declared a global pandemic in March 2020 by the World Health Organization (WHO). Since SARS-CoV-2 main protease plays an essential role in the virus’s life cycle, the design of small drug molecules with lower molecular weight has been a promising development targeting its inhibition. Herein, we evaluated the novel peptidomimetic azatripeptide and azatetrapeptide nitriles against SARS-CoV-2 main protease. We employed molecular dynamics (MD) simulations to elucidate the selected compounds’ binding free energy profiles against SARS-CoV-2 and further unveil the residues responsible for the drug-binding properties. Compound 8 exhibited the highest binding free energy of −49.37 ± 0.15 kcal/mol, followed by compound 7 (−39.83 ± 0.19 kcal/mol), while compound 17 showed the lowest binding free energy (−23.54 ± 0.19 kcal/mol). In addition, the absorption, distribution, metabolism, and excretion (ADME) assessment was performed and revealed that only compound 17 met the drug-likeness parameters and exhibited high pharmacokinetics to inhibit CYP1A2, CYP2C19, and CYP2C9 with better absorption potential and blood-brain barrier permeability (BBB) index. The additional intermolecular evaluations suggested compound 8 as a promising drug candidate for inhibiting SARS-CoV-2 Mpro. The substitution of isopropane in compound 7 with an aromatic benzene ring in compound 8 significantly enhanced the drug’s ability to bind better at the active site of the SARS-CoV-2 Mpro.

© 2022 Elsevier Inc.Type IV (T4) pilus is among the virulence factors with a key role in serious bacterial diseases. Specifically, in Neisseria meningitidis and Pseudomonas aeruginosa, it determines pathogenicity and causes infection. Here, a computational approach has been pursued to find piperidine-based inhibitor molecules against the elongation ATPase of T4 pili in these two selected pathogens. Using the modeled structures of the PilF and PilB ATPases of N. meningitidis and P. aeruginosa, virtual library screening via molecular docking has returned inhibitor molecule candidates. The dynamics of the best three binders have further been investigated in detail via molecular dynamic simulations. Among these, ligands with COCONUT IDs CNP0030078 and CNP0051517 were found to have higher potential in the inhibition of ATPases based on molecular dynamic simulation analysis and biological activity information. The obtained results will guide future efforts in antivirulence drug development against T4 pili of N. meningitidis and P. aeruginosa.

Guanine quadruplex nucleic acids (G4s) are involved in key biological processes such as replication or transcription. Beyond their biological relevance, G4s find applications as biotechnological tools since they readily bind hemin and enhance its peroxidase activity, creating a G4-DNAzyme. The biocatalytic properties of G4-DNAzymes have been thoroughly studied and used for biosensing purposes. Despite hundreds of applications and massive experimental efforts, the atomistic details of the reaction mechanism remain unclear. To help select between the different hypotheses currently under investigation, we use extended explicit-solvent molecular dynamics (MD) simulations to scrutinize the G4/hemin interaction. We find that besides the dominant conformation in which hemin is stacked atop the external G-quartets, hemin can also transiently bind to the loops and be brought to the external G-quartets through diverse delivery mechanisms. The simulations do not support the catalytic mechanism relying on a wobbling guanine. Similarly, catalytic role of the iron-bound water molecule is not in line with our results, however, given the simulation limitations, this observation should be considered with some caution. The simulations rather suggest tentative mechanisms in which the external G-quartet itself could be responsible for the unique H2O2-promoted biocatalytic properties of the G4/hemin complexes. Once stacked atop a terminal G-quartet, hemin rotates about its vertical axis while readily sampling shifted geometries where the iron transiently contacts oxygen atoms of the adjacent G-quartet. This dynamics is not apparent from the ensemble-averaged structure. We also visualize transient interactions between the stacked hemin and the G4 loops. Finally, we investigated interactions between hemin and on-pathway folding intermediates of the parallel-stranded G4 fold. The simulations suggest that hemin drives the folding of parallel-stranded G4s from slip-stranded intermediates, acting as a G4 chaperone. Limitations of the MD technique are briefly discussed.For Table of Contents Only

Z-DNA and Z-RNA have long appeared as oddities to nucleic acid scientists. However, their Z-step constituents are recurrently observed in all types of nucleic acid systems including ribosomes. Z-steps are NpN steps that are isostructural to Z-DNA CpG steps. Among their structural features, Z-steps are characterized by the presence of a lone pair…π contact that involves the stacking of the ribose O4′ atom of the first nucleotide with the 3′-face of the second nucleotide. Recently, it has been documented that the CpG step of the ubiquitous r(UNCG) tetraloops is a Z-step. Accordingly, such r(UNCG) conformations were called Z-turns. It has also been recognized that an r(GAAA) tetraloop in appropriate conditions can shapeshift to an unusual Z-turn conformation embedding an ApA Z-step. In this report, we explore the multiplicity of RNA motifs based on Z-steps by using the WebFR3D tool to which we added functionalities to be able to retrieve motifs containing lone pair…π contacts. Many examples that underscore the diversity and universality of these motifs are provided as well as tutorial guidance on using WebFR3D. In addition, this study provides an extensive survey of crystallographic, cryo-EM, NMR, and molecular dynamics studies on r(UNCG) tetraloops with a critical view on how to conduct database searches and exploit their results.

This study was aimed to perform the mechanistic investigations of chalcone scaffold as inhibitors of acetylcholinesterase (AChE) enzyme using molecular docking and molecular dynamics simulation tools. Basic chalcones (C1–C5) were synthesized and their in vitro AChE inhibition was tested. Binding interactions were studied using AutoDock and Surflex-Dock programs, whereas the molecular dynamics simulation studies were performed to check the stability of the ligand–protein complex. Good AChE inhibition (IC50 = 22 ± 2.8 to 37.6 ± 0.75 μM) in correlation with the in silico results (binding energies = −8.55 to −8.14 Kcal/mol) were obtained. The mechanistic studies showed that all of the functionalities present in the chalcone scaffold were involved in binding with the amino acid residues at the binding site through hydrogen bonding, π–π, π–cation, π–sigma, and hydrophobic interactions. Molecular dynamics simulation studies showed the formation of stable complex between the AChE enzyme and C4 ligand.