Cancer has substantial economic ramifications for healthcare systems. PIM kinases, specifically PIM-1, are commonly upregulated in different types of cancers, thereby promoting cancer development. PIM-1 inhibitors have garnered interest for their potential efficacy in cancer therapy. This study used computational methods to screen a library of 7,600 natural compounds targeting the PIM-1 active site. Five top candidates-ZINC00388658, ZINC00316459, ZINC00197401, ZINC00001673, and ZINC00316479-were selected for subsequent interaction studies, which involved molecular dynamics simulations (MDS) and free energy calculation using the MMPBSA method. These compounds interacted with key PIM-1 residues and had multiple common binding site interactions with the co-crystallized ligand 6YN, which was used as a control. Furthermore, the selected compounds exhibited favorable drug-like properties and stable docked complexes during a 200-ns molecular dynamics simulation,followed by MMPBSA analysis. Among the candidates, ZINC00388658 had the most favorable binding energy profile, indicating exceptional stability andintense interaction with PIM 1. This makes ZINC00388658 the most promising candidate for further development as a PIM-1 inhibitor. These findings suggest that ZINC00388658 and other promising compounds holdsignificant potential for developingnew cancer therapies that target PIM-1.

Hypoxia-inducible factors (HIFs) are transcription factors that regulate erythropoietin (EPO) synthesis and red blood cell (RBC) production. Prolyl-4-hydroxylase domain (PHD) enzymes are key regulators of HIF’s stability and activity. Inhibiting PHD enzymes can enhance HIF-mediated responses and have therapeutic potential for diseases such as anemia, cancer, stroke, ischemia, neurodegeneration, and inflammation. In this study, we searched for novel PHD inhibitors from four databases of natural products and synthetic compounds: AfroDb Natural Products, AnalytiCon Discovery Natural Product (NP), HIM-Herbal Ingredients In-Vivo Metabolism, and Herbal Ingredients’ Targets, with a total number of 13,597 compounds. We screened the candidate compounds by molecular docking and validated them by molecular dynamics simulations and free energy calculations. We identified four target hits (ZINC36378940, ZINC2005305, ZINC31164438, and ZINC67910437) that showed stronger binding affinity to PHD2 compared to the positive control, Vadadustat (AKB-6548), with docking scores of − 13.34 kcal/mol, − 12.76 kcal/mol, − 11.96 kcal/mol, − 11.41 kcal/mol, and − 9.04 kcal/mol, respectively. The target ligands chelated the active site iron and interacted with key residues (Arg 383, Tyr329, Tyr303) of PHD2, in a similar manner as Vadadustat. Moreover, the dynamic stability-based assessment revealed that they also exhibited stable dynamics and compact trajectories. Then the total binding free energy was calculated for each complex which revealed that the control has a TBE of − 31.26 ± 0.30 kcal/mol, ZINC36378940 reported a TBE of − 38.65 ± 0.51 kcal/mol, for the ZINC31164438 the TBE was − 26.16 ± 0.30 kcal/mol while the ZINC2005305 complex reported electrostatic energy of − 32.75 ± 0.58 kcal/mol. This shows that ZINC36378940 is the best hit than the other and therefore further investigation should be performed for the clinical usage. Our results suggest that these target hits are promising candidates that reserve further in vitro and in vivo validations as potential PHD inhibitors for the treatment of renal anemia, cancer, stroke, ischemia, neurodegeneration, and inflammation.

Exploring phytoconstituent for confronting the symptoms of polycystic ovarian syndrome: molecular dynamics simulation, quantum studies, free energy calculations and network analysis approaches

We characterize the self-assembly and phase behavior of Janus rods over a broad range of temperatures and volume fractions, using Langevin dynamics simulations and free energy calculations. The Janus rods consist of a line of fused overlapping spheres that interact via a soft-core repulsive potential, with the addition of an attractive pseudo-square-well tail to a fraction of the spheres (the coverage) ranging from 5% to 100% of sites. Competition between the stability of liquid crystal phases originating from shape anisotropy and assembly driven by directional interactions gives rise to a rich polymorphism that depends on the coverage. At low densities near the Boyle temperature, we observe the formation of spherical and tubular micelles at low coverages, while at higher coverages, randomly oriented monolayers form as the attractive parts of the rods overlap. At higher densities, bilayer structures appear and merge to form smectic and crystalline lamellar phases. All these structures gradually become unstable as the temperature is increased until eventually regular nematic and smectic phases appear, consistent with the hard rod limit. Our results indicate that the intermediate regime where shape-entropic effects compete with anisotropic attractions provided by site specificity is rich in structural possibilities and should help guide the design of rod-like colloids for specific applications.

Using a dual immunoinformatics and bioinformatics approach to design a novel and effective multi-epitope vaccine against human torovirus disease

ABSTRACT There is a consensus that epigenomic changes play a significant role in carcinogenesis. The effect of DNA methylation and its increased modulations on gene expression during carcinogenesis is significant for diagnostic and therapeutic purposes. Potential phytochemicals as anticancer approach modulators of epigenetic pathways are the focus of this investigation. Molecular docking studies were conducted to foretell how phytochemicals will interact with DNMT3A. As part of our effort to identify novel anti-cancer treatments, we examined a wide variety of phytochemicals from many chemical classes, including cannabinoids, carotenoids, chalcones, fatty acids, lignans, polysaccharides, saponins, steroids, and tannins. The docking scores for Dihydromyricetin are -10.207 kcal/mol, -9.313 kcal/mol for S-Adenosyl-L-methionine (SAM), and -8.757 kcal/mol for Zeylenol were determined to be superior than to phytochemical inhibitors against DNMT3A in the virtual screening method. Docking scores, hydrogen bond interactions, ADME characteristics, and DFT. Molecular Dynamics Simulation and free energy calculations demonstrated that these compounds bind stably to DNMT3A, potentially inhibiting its activity. Network Pharmacology suggested Dihydromyricetin's specific anticancer potential against cervical cancer. These findings provide insights into protein-Ligand interactions with Dnmt3A and highlights the need for In-vitro validation.

In silico analysis of aptamer-RNA conjugate interactions with human transferrin receptor

The dynamic restructuring of Cu has been observed under electrochemical conditions, and it has been hypothesized to underlie the unique reactivity of Cu toward CO2 electroreduction. Roughening is one of the key surface phenomena for Cu activation, whereby numerous atomic vacancies and adatoms form. However, the atomic structure of such surface motifs in the presence of relevant adsorbates has remained elusive. Here, we explore the chemical space of Cu surface restructuring under coverage of CO and H in realistic electroreduction conditions, by combining grand canonical DFT and global optimization techniques, from which we construct a potential-dependent grand canonical ensemble representation. The regime of intermediate and mixed CO and H coverage─where structures exhibit some elevated surface Cu─is thermodynamically unfavorable yet kinetically inevitable. Therefore, we develop a quasi-kinetic Monte Carlo simulation to track the system's evolution during a simulated cathodic scan. We reveal the evolution path of the system across coverage space and identify the accessible metastable structures formed along the way. Chemical bonding analysis is performed on the metastable structures with elevated Cu*CO species to understand their formation mechanism. By molecular dynamics simulations and free energy calculations, the surface chemistry of the Cu*CO species is explored, and we identify plausible mechanisms via which the Cu*CO species may diffuse or dimerize. This work provides rich atomistic insights into the phenomenon of surface roughening and the structure of involved species. It also features generalizable methods to explore the chemical space of restructuring surfaces with mixed adsorbates and their nonequilibrium evolution.

Minocycline, a repurposed approved medication, shows promise in treating neurodegeneration. However, the specific pathways targeted by minocycline remain unclear despite the identification of molecular targets. This study explores minocycline's potential protective effects against TNF-α-mediated neuronal death in PC12 cells, with a focus on unraveling its interactions with key molecular targets. The study begins by exploring minocycline's protective role against TNF-α-mediated neuronal death in PC12 cells, showcasing a substantial reduction in cleaved caspase-3 expression, DNA fragmentation, and intracellular ROS levels following minocycline pretreatment. Subsequently, a comprehensive analysis utilizing pull-down assays, computational docking, mutation analysis, molecular dynamics simulations, and free energy calculations is conducted to elucidate the direct interaction between minocycline and p47phox-the organizer subunit of NADPH oxidase-2 (NOX2) complex. Computational insights, including a literature survey and analysis of key amino acid residues, reveal a potential binding site for minocycline around Trp193 and Cys196. In silico substitutions of Trp193 and Cys196 further confirm their importance in binding with minocycline. These integrated findings underscore minocycline's protective mechanisms, linking its direct interaction with p47phox to the modulation of NOX2 activity and attenuation of NOX-derived ROS generation. Minocycline demonstrates protective effects against TNF-α-induced PC12 cell death, potentially linked to its direct interaction with p47phox. This interaction leads to a reduction in NOX2 complex assembly, ultimately attenuating NOX-derived ROS generation. These findings hold significance for researchers exploring neuroprotection and the development of p47phox inhibitors.

Human β-cardiac myosin plays a critical role in generating the mechanical forces necessary for cardiac muscle contraction. This process relies on a delicate dynamic equilibrium between the disordered relaxed state (DRX) and the super-relaxed state (SRX) of myosin. Disruptions in this equilibrium due to mutations can lead to heart diseases. However, the structural characteristics of SRX and the molecular mechanisms underlying pathogenic mutations have remained elusive. To bridge this gap, we conducted molecular dynamics simulations and free energy calculations to explore the conformational changes in myosin. Our findings indicate that the size of the phosphate-binding pocket can serve as a valuable metric for characterizing the transition from the DRX to SRX state. Importantly, we established a global dynamic coupling network within the myosin motor head at the residue level, elucidating how the pathogenic mutation E483K impacts the equilibrium between SRX and DRX through allosteric effects. Our work illuminates molecular details of SRX and offers valuable insights into disease treatment through the regulation of SRX.

Objective: Estrogen receptor (ER) inhibitors have significant therapeutic potential for hormone-dependent cancers and related disorders. Tamoxifen, a well-known selective estrogen receptor modulator, has been widely used as adjuvant therapy for estrogen receptor-positive breast cancer. However, tamoxifen may exhibit a tendency to develop resistance with prolonged usage and particularly elevate the risk of uterine cancer. Therefore, there is a need for the discovery and development of new ER modulators or inhibitors. In this study, we identified potential estrogen receptor inhibitors through computational drug repositioning. Methods: A set of 2048 compounds, encompassing FDA-approved drugs and active metabolites, were subjected to molecular docking, molecular dynamics simulations, and free energy calculations to evaluate their interaction with estrogen receptor α (ERα). Results: Among the compounds evaluated, conivaptan, atogepant, and lomitapide exhibited the highest affinities for ERα. Lomitapide displayed a superior docking score (-12 kcal/mol) compared to the established ER inhibitor, tamoxifen (-10 kcal/mol). Further investigation using molecular dynamics simulations and free energy calculations disclosed lomitapide's heightened binding affinity of -380.727 kJ/mol, surpassing tamoxifen's binding affinity of -352.029 kJ/mol. Conclusion: This comprehensive computational exploration underscores lomitapide's potential as a compelling candidate with an envisaged stronger estrogen receptor affinity than the acknowledged standard, tamoxifen. To validate lomitapide's promise as a novel ER inhibitor, essential in vitro and in vivo studies are suggested. These investigations will provide essential insights into lomitapide's reposition in addressing the challenges tied to hormone-dependent cancers and associated maladies.

Effective scaffolding of immunogens is crucial for generating conformationally selective antibodies through active immunization, particularly in the treatment of protein misfolding diseases such as Alzheimer's and Parkinson's disease. Previous computational work has revealed that a disorder-prone region of the tau protein, when in a stacked form, is predicted to structurally resemble a small, soluble protofibril, having conformational properties similar to those of experimental in vitro tau oligomers. Such an oligomeric structural mimic has the potential to serve as a vaccine immunogen design for Alzheimer's disease. In this study, we developed a cyclization scaffolding method in Rosetta, in which multiple cyclic peptides are stacked into a protofibril. Cyclization results in significant stabilization of protofibril-like structures by constraining the conformational space. Applying this method to the disorder-prone region of the tau fibril, we evaluated the metastability of the cyclized tau immunogen using molecular dynamics simulations, and we identified sequences of two cyclic constructs having high metastability in the protofibril. We then assessed their thermodynamic stability by computing the free energy required to separate a distal chain from the rest of the stacked structure. Our computational results, based on molecular dynamics simulations and free energy calculations, demonstrate that two cyclized constructs, cyclo-(VKSEKLDFKDRVQSKIFyN) and cyclo-(VKSEKLDFKDRVQSKIYvG) (lowercase letters indicate d-form amino acids), possess significantly increased thermodynamic stability in the protofibril over an uncyclized linear construct VKSEKLDFKDRVQSKI. The cyclization scaffolding approach proposed here holds promise as a means to effectively design immunogens for protein misfolding diseases, particularly those involving liposome-conjugated peptide constructs.

SARS-CoV-2 (Severe Acute Respiratory Syndrome Coronavirus 2) is the etiological agent responsible for the global outbreak of COVID-19 (Coronavirus Disease 2019). The main protease of SARS-CoV-2, Mpro, is a key enzyme that plays a vital role in mediating viral replication and transcription. In this study, a comprehensive computational approach was employed to investigate the binding affinity, selectivity, and stability of natural product candidates as potential new antivirals acting on the viral polyprotein processing mediated by SARS-CoV-2 Mpro. A library of 288 flavonoids extracted from Brazilian biodiversity was screened to select potential Mpro inhibitors. An initial filter based on Lipinski's rule of five was applied, and 204 compounds that did not violate any of the Lipinski rules were selected. The compounds were then docked into the active site of Mpro using the GOLD program, and the poses were subsequently re-scored using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) binding free energy calculations performed by AmberTools23. The top five flavonoids with the best MM-GBSA binding free energy values were selected for analysis of their interactions with the active site residues of the protein. Next, we conducted a toxicity and drug-likeness analysis, and non-toxic compounds were subjected to molecular dynamics simulation and free energy calculation using the MM-PBSA (Molecular Mechanics Poisson-Boltzmann Surface Area) method. It was observed that the five selected flavonoids had lower MM-GBSA binding free energy with Mpro than the co-crystal ligand. Furthermore, these compounds also formed hydrogen bonds with two important residues, Cys145 and Glu166, in the active site of Mpro. Two compounds that passed the drug-likeness filter showed stable conformations during the molecular dynamics simulations. Among these, NuBBE_867 exhibited the best MM-PBSA binding free energy value compared to the crystallographic inhibitor. Therefore, this study suggests that NuBBE_867 could be a potential inhibitor against the main protease of SARS-CoV-2 and may be further examined to confirm our results.

The major histocompatibility complex (MHC) can recognize and bind to external peptides to generate effective immune responses by presenting the peptides to T cells. Therefore, understanding the binding modes of peptide-MHC complexes (pMHC) and predicting the binding affinity of pMHCs play a crucial role in the rational design of peptide vaccines. In this study, we employed molecular dynamics (MD) simulations and free energy calculations with an Alanine Scanning with Generalized Born and Interaction Entropy (ASGBIE) method to investigate the protein-peptide interaction between HLA-A*02:01 and the G9209 peptide derived from the melanoma antigen gp100. The energy contribution of individual residue was calculated using alanine scanning, and hotspots on both the MHC and the peptides were identified. Our study shows that the pMHC binding is dominated by the van der Waals interactions. Furthermore, we optimized the ASGBIE method, achieving a Pearson correlation coefficient of 0.91 between predicted and experimental binding affinity for mutated antigens. This represents a significant improvement over the conventional MM/GBSA method, which yields a Pearson correlation coefficient of 0.22. The computational protocol developed in this study can be applied to the computational screening of antigens for the MHC1 as well as other protein-peptide binding systems.

Objective: Pharmacoinformatics is an innovative approach rapidly evolving in pharmaceutical research and drug development. This study focuses on analysing Morus macroura, a plant species with untapped pharmacological potential. This investigation aims to leverage pharmacoinformatics techniques to unveil the hidden potential of Morus macroura in drug discovery and development. Methods: The study includes analyses of protein-protein interactions, deep learning docking, adsorption tests, distribution, metabolism, excretion, molecular dynamics simulations and free energy calculation using Molecular Mechanics Generalized Born Surface Area (MMGBSA). Results: Nine active compounds were identified in Morus macroura, namely Andalasin A, Guangsangon K, Guangsangon L, Guangsangon M, Guangsangon N, Macrourone C, Mulberrofuran G, Mulberrofuran K, and Mulberroside C. These compounds exhibit protein-protein interaction activities against a cytochrome P450 monooxygenase that catalyses the conversion of C19 androgens. These plant compounds influence aromatase excess syndrome, deficiency, and ovarian dysgenesis. Regarding drug-likeness, Mulberroside C and Macrourone C demonstrated good absorption potential by adhering to Lipinski's rule of five. Deep learning docking simulations yielded affinity results of-9.62 kcal/mol for Guangsangon M,-10.44 kcal/mol for Macrourone C, and-10.99 kcal/mol for Guangsangon L. Subsequent molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during a 100 ns simulation. Conclusion: Morus macroura interacts with important proteins, particularly CYP19A1, which might influence health conditions like aromatase excess syndrome and ovarian dysgenesis. These findings provide potential paths for addressing specific health issues and advancing drug development. Molecular dynamics simulations indicated that Guangsangon L and Macrourone C remained stable during simulation.

The allosteric modulation of the homodimeric H10-03-6 protein to glycan ligands L1 and L2, and the STAB19 protein to glycan ligands L3 and L4, respectively, has been studied by molecular dynamics simulations and free energy calculations. The results revealed that the STAB19 protein has a significantly higher affinity for L3 (-11.38 ± 2.32 kcal/mol) than that for L4 (-5.51 ± 1.92 kcal/mol). However, the combination of the H10-03-6 protein with glycan L2 (1.23 ± 6.19 kcal/mol) is energetically unfavorable compared with that of L1 (-13.96 ± 0.35 kcal/mol). Further, the binding of glycan ligands L3 and L4 to STAB19 would result in the significant closure of the two CH2 domains of the STAB19 conformation with the decrease of the centroid distances between the two CH2 domains compared with the H10-03-6/L1/L2 complex. The CH2 domain closure of STAB19 relates directly to the formation of new hydrogen bonds and hydrophobic interactions between the residues Ser239, Val240, Asp265, Glu293, Asn297, Thr299, Ser337, Asp376, Thr393, Pro395, and Pro396 in STAB19 and glycan ligands L3 and L4, which suggests that these key residues would contribute to the specific regulation of STAB19 to L3 and L4. In addition, the distance analysis revealed that the EF loop in the H10-03-6/L1/L2 model presents a high flexibility and partial disorder compared with the stabilized STAB19/L3/L4 complex. These results will be helpful in understanding the specific regulation through the asymmetric structural characteristics in the CH2 and CH3 domains of the H10-03-6 and STAB19 proteins.

Global prevalence of neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease is increasing gradually, whereas approvals of successful therapeutics for central nervous system disorders are inadequate. Accumulating evidence suggests pivotal roles of the receptor-interacting serine/threonine-protein kinase 1 (RIPK1) in modulating neuroinflammation and necroptosis. Discoveries of potent small molecule inhibitors for RIPK1 with favorable pharmacokinetic properties could thus address the unmet medical needs in treating neurodegeneration. In a structure-based virtual screening, we performed site-specific molecular docking of 4,858 flavonoids against the kinase domain of RIPK1 using AutoDock Vina. We predicted physicochemical descriptors of the top ligands using the SwissADME webserver. Binding interactions of the best ligands and the reference ligand L8D were validated using replicated 500-ns Gromacs molecular dynamics simulations and free energy calculations. From Vina docking, we shortlisted the top 20 flavonoids with the highest binding affinities, ranging from -11.7 to -10.6 kcal/mol. Pharmacokinetic profiling narrowed down the list to three orally bioavailable and blood-brain-barrier penetrant flavonoids: Nitiducarpin, Pinocembrin 7-O-benzoate, and Paratocarpin J. Next, trajectories of molecular dynamics simulations of the top protein-ligand complexes were analyzed for binding interactions. The root-mean-square deviation (RMSD) was 1.191 Å (±0.498 Å), 1.725 Å (±0.828 Å), 1.923 Å (±0.942 Å), 0.972 Å (±0.155 Å) for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D, respectively. The radius of gyration (Rg) was 2.034 nm (±0.015 nm), 2.0.39 nm (± 0.025 nm), 2.053 nm (±0.021 nm), 2.037 nm (±0.016 nm) for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D, respectively. The solvent accessible surface area (SASA) was 159.477 nm2 (±3.021 nm2), 159.661 nm2 (± 3.707 nm2), 160.755 nm2 (±4.252 nm2), 156.630 nm2 (±3.521 nm2), for Nitiducarpin, Pinocembrin 7-O-benzoate, Paratocarpin J, and L8D complexes, respectively. Therefore, lower RMSD, Rg, and SASA values demonstrated that Nitiducarpin formed the most stable complex with the target protein among the best three ligands. Finally, 2D protein-ligand interaction analysis revealed persistent hydrophobic interactions of Nitiducarpin with the critical residues of RIPK1, including the catalytic triads and the activation loop residues, implicated in the kinase activity and ligand binding. Our target-based virtual screening identified three flavonoids as strong RIPK1 inhibitors, with Nitiducarpin exhibiting the most potent inhibitory potential. Future in vitro and in vivo studies with these ligands could offer new hope for developing effective therapeutics and improving the quality of life for individuals affected by neurodegeneration.

(De)carboxylation mechanisms of heteroaromatic substrates catalyzed by prenylated FMN-dependent UbiD decarboxylases: An in-silico study

The human Na+/H+ exchanger (NHE1) plays a crucial role in maintaining intracellular pH by regulating the electroneutral exchange of a single intracellular H+ for one extracellular Na+ across the plasma membrane. Understanding the molecular mechanisms governing ion transport and the binding of inhibitors is of importance in the development of anticancer therapeutics targeting NHE1. In this context, we performed molecular dynamics (MD) simulations based on the recent cryo-electron microscopy (cryo-EM) structures of outward- and inward-facing conformations of NHE1. These simulations allowed us to explore the dynamics of the protein, examine the ion-translocation pore, and confirm that Asp267 is the ion-binding residue. Our free energy calculations did not show a significant difference between Na+ and K+ binding at the ion-binding site. Consequently, Na+ over K+ selectivity cannot be solely explained by differences in ion binding. Our MD simulations involving NHE1 inhibitors (cariporide and amiloride analogues) maintained stable interactions with Asp267 and Glu346. Our study highlights the importance of the salt bridge between the positively charged acylguanidine moiety and Asp267, which appears to play a role in the competitive inhibitory mechanism for this class of inhibitors. Our computational study provides a detailed mechanistic interpretation of experimental data and serves the basis of future structure-based inhibitor design.

Nanoparticles' (NPs) permeation through cell membranes, whether it happens via passive or active transport, is an essential initial step for their cellular internalization. The NPs' surface coating impacts the way they translocate through the lipid bilayer and the spontaneity of the process. Understanding the molecular details of NPs' interaction with cell membranes allows the design of nanosystems with optimal characteristics for crossing the lipid bilayer: computer simulations are a powerful tool for this purpose. In this work, we have performed coarse-grained molecular dynamics simulations and free energy calculations on spherical titanium dioxide NPs conjugated with polymer chains of different chemical compositions. We have demonstrated that the hydrophobic/hydrophilic character of the chains, more than the nature of their terminal group, plays a crucial role in determining the NPs' interaction with the lipid bilayer and the thermodynamic spontaneity of NPs' translocation from water to the membrane. We envision that this computational work will be helpful to the experimental community in terms of the rational design of NPs for efficient cell membrane permeation.