Dynamical Cross-correlation Matrix Research Articles

4D-QSAR, ADMET properties, and molecular dynamics simulations for designing N-substituted urea/thioureas as human glutaminyl cyclase inhibitors

Human glutaminyl cyclase (hQC) inhibitors have great potential to be used as anti- Alzheimer's disease (AD) agents by reducing the toxic pyroform of β-amyloid in the brains of AD patients. The four-dimensional quantitative structure activity relationship (4D-QSAR) model of N-substituted urea/thioureas was established with satisfying predictive ability and statistical reliability (Q2 = 0.521, R2 = 0.933, R2prep = 0.619). By utilizing the developed 4D-QSAR model, a set of new N-substituted urea/thioureas was designed and evaluated for their Absorption Distribution Metabolism Excretion and Toxicity (ADMET) properties. The results of molecular dynamics (MD) simulations, Principal component analysis (PCA), free energy landscape (FEL), dynamic cross-correlation matrix (DCCM) and molecular mechanics generalized Born Poisson-Boltzmann surface area (MM-PBSA) free energy calculations, revealed that the designed compounds were remained stable in protein binding pocket and compounds b ∼ f (−35.1 to −44.55 kcal/mol) showed higher binding free energy than that of compound 14 (−33.51 kcal/mol). The findings of this work will be a theoretical foundation for further research and experimental validation of urea/thiourea derivatives as hQC inhibitors.

Unveiling the potential of Traditional Chinese Medicines in combating NorA-mediated S. aureus drug resistance. A molecular dynamic study

Drug resistance Staphylococcus aureus (S. aureus) is emerging as the most prevalent pathogen in surgical wounds and orthopedic surgeries. NorA is a major drug efflux involved in the drug resistance, reducing intracellular drug concentrations. The current study aimed to find some potential inhibitors derived from Traditional Chinese Medicines (TCM) against NorA using a 500 ns molecular dynamics (MD) simulation. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), free energy, dynamic cross correlation matrix (DCCM) was computed. Three TCM compound; piperine, demethoxycurcumin, and tectorigenin exhibited many hydrogen bonds and hydrophobic contacts in the binding pockets residues, especially with ionizable residues (Arg 98, Glu 222, Asp 307, and Arg 310), forming two different charged regions. RMSD of NorA-piperine complex shows some initial instability at 5.2 Å, followed by an increase to 5.8 Å at 222 ns, and a further rise to 6.1 Å by the end of the 500 ns simulation period. Similarly, the NorA-demethoxycurcumin and NorA-tectorigenin complexes display elevated terminal RMSD of 7.9 Å each. RMSF, Rg, free energy, and binding pockets interactions confirmed the potential activity of piperine, demethoxycurcumin, and tectorigenin against NorA inhibition. In apo state NorA exhibited dispersion within the ranges of −100 to 150 (PC1: 39.97 %) while in complexed with ligands it illustrated significant differences in the motion patterns on PC1 (46.04 % and 27.39 %). The NorA-ligands DCCM displayed predominantly negatively correlated motions, particularly at the terminal resides from 200 to 370 where the pocket binding residues are located. These compounds also demonstrated drug-likeness and fitting best in the binding pocket residues (Asn > 137, Phe > 140, Phe > 303 Asp > 307, Arg > 310, Glu > 222, and Thr > 223). In conclusion, the compounds piperine, demethoxycurcumin, and tectorigenin interact with NorA through hydrogen bonds and hydrophobic contacts, particularly with key ionizable residues, Arg 98, Glu 222, Asp 307, and Arg 310, indicating potential inhibitory activity.

Impact of Phosphorylation at Various Sites on the Active Pocket of Human Ferrochelatase: Insights from Molecular Dynamics Simulations.

Ferrochelatase (FECH) is the terminal enzyme in human heme biosynthesis, catalyzing the insertion of ferrous iron into protoporphyrin IX (PPIX) to form protoheme IX (Heme). Phosphorylation increases the activity of FECH, and it has been confirmed that the activity of FECH phosphorylated at T116 increases. However, it remains unclear whether the T116 site and other potential phosphorylation modification sites collaboratively regulate the activity of FECH. In this study, we identified a new phosphorylation site, T218, and explored the allosteric effects of unphosphorylated (UP), PT116, PT218, and PT116 + PT218 states on FECH in the presence and absence of substrates (PPIX and Heme) using molecular dynamics (MD) simulations. Binding free energies were evaluated with the MM/PBSA method. Our findings indicate that the PT116 + PT218 state exhibits the lowest binding free energy with PPIX, suggesting the strongest binding affinity. Additionally, this state showed a higher binding free energy with Heme compared to UP, which facilitates Heme release. Moreover, employing multiple analysis methods, including free energy landscape (FEL), principal component analysis (PCA), dynamic cross-correlation matrix (DCCM), and hydrogen bond interaction analysis, we demonstrated that phosphorylation significantly affects the dynamic behavior and binding patterns of substrates to FECH. Insights from this study provide valuable theoretical guidance for treating conditions related to disrupted heme metabolism, such as various porphyrias and iron-related disorders.

Quercetin derivatives as potential inhibitors of Nipah virus phosphoprotein through in silico drug design approaches

The Nipah virus (NiV) is a pathogenic infection presenting a substantial risk to both human and animal communities. Despite its virulence, there are currently no available drugs or vaccines for NiV. To fight against this challenge, we applied a multi-step in silico drug design strategy, including PASS prediction, molecular docking, absorption, distribution, metabolism and excretion (ADMET) analysis, molecular dynamics simulations, frontier molecular orbitals calculations, dynamic cross-correlation matrix (DCCM), and principal component analysis (PCA). To find out the potential drug candidate against NiV target protein, more than 100 derivatives were taken from the PubChem database. After that, the Pa (probability “to be active"), and Pi (probability “to be inactive") calculation was conducted. Based on the maximum probability “to be active score (Pa), the top nine compounds were studied against NiV. Our findings suggest that the quercetin derivatives exhibit promising binding affinities with the Nipah virus phosphoprotein, particularly compounds 03 (−7.8 kcal/mol), 22 (−6.9 kcal/mol), and 89 (−7.3 kcal/mol). The molecular dynamics simulations over the 100 had confirmed excellent stability and ADMET profiles has reported that all the mentioned ligands are free from Hepatotoxicity, and AMES toxicity, suggesting them as promising candidates for anti-NiV drugs. Based on the results of molecular dynamics simulations, it was observed that the reported drug candidates 03, 22, and 89 exhibited a high degree of stability in forming protein-ligand complexes. These complexes showed only minor fluctuations in RMSF over the 100 ns simulation period. The utilization of frontier molecular orbitals aided in the understanding of the electronic structure behavior, providing strong evidence for their suitability. The study indicates that certain quercetin derivatives, such as 03, 22, and 89, show promise as potentially effective antiviral drug candidates against NiV. Future research should focus on experimental testing of these compounds to develop new antiviral drugs against NiV.

Investigating natural plant products as potential inhibitors to disrupt NS1 β-roll domain polymerisation in DENV2: a detailed computational chemistry approach

ABSTRACT The aim of this investigation is to design a novel drug, screened through various computational methods, that harnesses the historical evidence of plant products inhibiting viral infections. The beta-roll domain of Flavivirus NS1, crucial for its stability, membrane interaction and viral RNA replication, underscores its significance in disease pathogenesis, being the sole non-structural protein detected in the sera of DENV-infected patients. Our goal is to develop optimised compounds that disrupt NS1 β-roll domain polymerisation, mitigating hydrophobic core formation, to effectively inhibit DENV2 cytopathic infection by leveraging natural interference with viral replication. We created 10 plant-based ligand libraries and applied scaffold hopping and bioisosteric replacement to improve them. Computational methodologies, including virtual screening, molecular dynamics simulation, contact matrix and dynamic cross correlation matrix, were employed for analysing the dynamics of the protein with the association of ligand. ADMET analysis of ligand 126 revealed a favourable ADME profile with no observed toxicity properties. Calculations of pKd, pKi and pIC50 values indicated a heightened affinity of the ligand for the receptor protein, suggesting enhanced inhibitory action and a reduced dissociation rate. These encouraging findings suggest a greater potency of the drug, warranting further in vitro testing for validation.

Identification of Helicobacter pylori-carcinogenic TNF-alpha-inducing protein inhibitors via daidzein derivatives through computational approaches.

Gastric cancer is considered a class 1 carcinogen that is closely linked to infection with Helicobacter pylori (H. pylori), which affects over 1 million people each year. However, the major challenge to fight against H. pylori and its associated gastric cancer due to drug resistance. This research gap had led our research team to investigate a potential drug candidate targeting the Helicobacter pylori-carcinogenic TNF-alpha-inducing protein. In this study, a total of 45 daidzein derivatives were investigated and the best 10 molecules were comprehensively investigated using in silico approaches for drug development, namely pass prediction, quantum calculations, molecular docking, molecular dynamics simulations, Lipinski rule evaluation, and prediction of pharmacokinetics. The molecular docking study was performed to evaluate the binding affinity between the target protein and the ligands. In addition, the stability of ligand-protein complexes was investigated by molecular dynamics simulations. Various parameters were analysed, including root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), radius of gyration (Rg), hydrogen bond analysis, principal component analysis (PCA) and dynamic cross-correlation matrix (DCCM). The results has confirmed that the ligand-protein complex CID: 129661094 (07) and 129664277 (08) formed stable interactions with the target protein. It was also found that CID: 129661094 (07) has greater hydrogen bond occupancy and stability, while the ligand-protein complex CID 129664277 (08) has greater conformational flexibility. Principal component analysis revealed that the ligand-protein complex CID: 129661094 (07) is more compact and stable. Hydrogen bond analysis revealed favourable interactions with the reported amino acid residues. Overall, this study suggests that daidzein derivatives in particular show promise as potential inhibitors of H. pylori.

Insights into novel inhibitors intending HCMV protease a computational molecular modelling investigation for antiviral drug repurposing

Human cytomegalovirus (HCMV) is an important opportunistic pathogen that is the most significant viral cause of congenital birth abnormalities and is responsible for a high morbidity and death rate in immunocompromised patients. HCMV's rising severity and the limitations of current vaccines in preventing infection highlight the urgent need for effective antiviral drugs. This study aimed to identify small compounds using extensive & appropriate in-silico drug design techniques, including molecular docking, Post-docking MM-GBSA, ADMET, MD simulation, PCA, and DCCM were employed in this study capable of binding to the viral protease and inhibiting its activity, thereby preventing proteolytic processing during capsid maturation. 3516 compounds from life chemical were used in molecular docking, and the top four compounds having high binding affinity and promising ADMET properties with life chemicals ID: F3407-0101 (CID: 49667672), F6559-3323 (CID: 121022124), F6456-1266 (CID: 71810903), and F3411-7969 (CID: 50785034). MD simulation was also used to assess the stability of the protein-ligand complex structure. Finally, after MD simulation, principal component analysis (PCA), and dynamic cross-correlation matrix (DCCM) analysis were performed using trajectories, and we can suggest the best drug candidate which is CID: 50785034 (N-(3-fluoro-4-methylphenyl)-2-({7-oxo-8-[(oxolan-2-yl) methyl]-5-thia-1,8,10,11-tetraazatricyclo [7.3.0.0^ {2,6}] dodeca-2(6),3,9,11-tetraen-12-yl}sulfanyl)acetamide), another two compounds CID: 121022124, and CID: 71810903 which comes after CID: 50785034. All three compounds may have the potential to be developed as a therapy option for HCMV infection.

In silico computational prediction of Saussurea pulchella compounds with inhibitory effects on plasmepsin X in Plasmodiumfalciparum

Malaria poses serious problems in countries located in tropical and subtropical regions worldwide. Plasmodium falciparum is responsible for approximately 90 % of malaria-related deaths worldwide. Therefore, alternative treatments are required to effectively reduce the number of malaria cases. Therefore, the objective of this study was to explore alternative malaria treatments by identifying compounds targeting P. falciparum plasmepsin X and characterizing the interactions under various temperatures. Molecular docking was used to evaluate interactions between polyphenolic compounds from Saussurea pulchella and plasmepsin X, followed by stability analyses using molecular dynamics (MD) simulations. Selected compounds were analyzed for biological activity, reactivity, and toxicity. RMSD (root mean square deviation), RMSF (root mean square fluctuation), Rg (radius of gyration), SASA (solvent-accessible surface area), and H-bonds were evaluated. Furthermore, MM-PBSA (molecular mechanics Poisson–Boltzmann surface area), PCA (principal component analysis), and DCCM (dynamic cross-correlation matrix) were use to evaluate effects of various temperatures. Pinoresinol 4-glucoside had a binding energy of 10.006 kcal/mol and was predicted to be nontoxic. MD simulations of pinoresinol 4-glucoside were conducted at different temperatures (298 K, 300 K, 305 K, 310 K, and 320 K) over a trajectory of 150 ns, revealing a strong binding energy using the MM-PBSA method. Our findings suggest that pinoresinol 4-glucoside is an alternative treatment for malaria and should be further evaluated at the laboratory and clinical levels.

Comprehending conformational changes in EmrE, multidrug transporter at different pH: insights from molecular dynamics simulations

EmrE is a small multidrug resistance (SMR) pump of antiparallel topology that confers resistance to a broad range of polyaromatic cations in Escherichia coli. Atomic-level understanding of conformational changes for the selectivity of substrate and transport of a diverse array of drugs through the smallest known efflux pumps is crucial to multi-drug resistance. Therefore, the present study aims to provide insights into conformational changes during the transport through EmrE transporter at different pH. Molecular dynamics simulations have been carried out on the complete structure of EmrE in the absence of substrate. Computational analyses such as secondary structure, principal component, dynamic cross-correlation matrix, and hydrogen bond calculations have been performed. Analysis of MD trajectories in this study revealed pH-dependent interactions that influenced the structural dynamics of EmrE. Notably, at high pH, Glu14 and Tyr60 in both monomers formed electrostatic interactions, while these interactions decreased significantly at a low pH. Interestingly, a kink at helix 3 (H3) and dual open conformation of EmrE at low pH were also observed in contrast to a closed state discerned towards the periplasmic side at high pH. Significant interactions between C-terminal residues and residues at the edge of H1 & Loop1 and H3 & Loop3 were identified, suggesting their role in stabilizing the closed conformation of EmrE at the periplasmic end under high pH conditions. The present study enhances our understanding of EmrE’s conformational changes, shedding light on the pH-dependent mechanisms that are likely to impact its function in multi-drug resistance. Communicated by Ramaswamy H. Sarma

Novel computational and drug design strategies for the inhibition of human T-cell leukemia virus 1-associated lymphoma by Astilbin derivatives

Human T-cell leukemia virus 1 (HTLV-1) associated lymphoma is a devastating malignancy triggered by HTLV-1 infections. We employeda comprehensive drug design and computational strategy in this work to explore the inhibitory activitiesof Astilbin derivatives against HTLV-1-associated lymphoma. We evaluated the stability, binding affinities, and various computational analysis of Astilbin derivatives against target proteins, such as HTLV-1 main protease and HTLV-1 capsid protein. The root mean square deviation (RMSD), root mean square fluctuation, radius of gyration, hydrogen bond analysis, principal component analysis (PCA) and dynamic cross-correlation matrix (DCCM) were applied to characterize these protein–ligand interactions further. Ligand-03 and ligand-04 exhibited notable binding affinity to HTLV-1 capsid protein, while ligand-05 displayed high binding affinity to HTLV-1 protease. MD simulation analysis revealed that ligand-03, bound to HTLV-1 capsid protein, demonstrated enhanced stability with lower RMSD values and fewer conformational changes, suggesting a promising binding orientation. Ligand-04, despite stable binding, exhibited increased structural deviations, making it less suitable. Ligand-05 demonstrated stable binding to HTLV-1 protease throughout the simulation period at 100 nanoseconds. Hydrogen bond analysis indicated that ligand-05 formed persistent hydrogen bonds with significantresidues, contributing to its stability. PCA highlighted ligand-03's more remarkable conformational changes, while DCCM showed ligand-05's distinct dynamics, indicating its different behavior in the complex. Furthermore, binding free energy calculations supported the favorable interactions of ligand-03 and ligand-04 with HTLV-1 capsid protein, while ligand-05 showed weaker interactions with HTLV-1 protease. Molecular electrostatic potential and frontier molecular orbital analyses provided insights into these compounds’ charge distribution and stability. In conclusion, this research found Astilbin derivatives as potential inhibitors of HTLV-1-associated lymphoma. Future attempts at drug development will benefit from the steady interaction landscape provided by Ligand-03, Ligand-04 and Ligand-05, which showed the most attractive binding profile with the target protein. These results open up new opportunities for innovative drug development, and more experimental testing should be done between Astilbin derivatives and HTLV-1-associated lymphoma. Communicated by Ramaswamy H. Sarma

Structure-oriented engineering of nitrile hydratase: Reshaping of substrate access tunnel and binding pocket for efficient synthesis of cinnamamide

Cinnamamide and its derivatives are the most common and important building blocks widely present in natural products. Currently, nitrile hydratase (NHase, EC 4.2.1.84) has been widely used in large-scale industrial production of nicotinamide and acrylamide, while its catalytic activity is extremely low or inactive for bulky nitrile substrates such as cinnamonitrile. Therefore, beneficial variant βF37P/L48P/F51N were obtained from PtNHase of Pseudonocardia thermophila JCM3095 by reshaping of substrate access tunnel and binding pocket, which exhibited 14.88-fold improved catalytic efficiency compared to the wild-type PtNHase. Structure analysis, molecular dynamics simulations and dynamical cross-correlation matrix (DCCM) analysis revealed that the introduced mutations enlarged the substrate access tunnel and binding pocket, enhanced overall anti-correlated movements of enzymes, which would promote product release during the dynamic process of catalysis. In a hydration process, the complete conversion of 5 mM cinnamonitrile was achieved by βF37P/L48P/F51N in a 50 mL reaction, with cinnamamide yield of almost 100 % and productivity of 0.736 g L−1 h−1. The study demonstrates the co-evolution of substrate access tunnel and binding pocket is an effective strategy, and provides a valuable reference for future research. Furthermore, NHases have huge potential for catalyzing bulky nitriles to form corresponding amides in large-scale industrial production.

Computational screening of potential bromodomain-containing protein 2 inhibitors for blocking SARS-CoV-2 infection through pharmacophore modeling, molecular docking and molecular dynamics simulation

Recently, human bromodomain-containing protein 2 (BRD2) has been reported as a potential drug target in host cells to combat SARS-CoV-2 as well as to treat COVID-19. In this study, we aimed to screen potential BRD2 inhibitors as anti- SARS-CoV-2 agents through computational approach. A pharmacophore model was generated with crystal structure of BRD2 in complex with RVX208 using ZINCPharmer server, and 10,842 compounds were screened from ZINC15 database against the pharmacophore model. Next, 102 compounds exhibited more affinity to BRD2 than control compound (RVX208) were screened using molecular docking. The top five hits (ZINC20417563, ZINC12322175, ZINC12391479, ZINC13098192, and ZINC20407881), which showed docking affinity to BRD2 with −10, −9.9, −9.8, −9.7 and −9.7 kcal/mol, respectively, were subjected to 100 ns molecular dynamics simulation to evaluate the stability of docked complexes and their interaction mechanisms. The appropriate stability of these complexes was displayed, and these compounds were found to be closely adhered to the binding site of the BRD2 protein based on the global dynamics analysis with RMSD, RMSF, Rg, SASA, DSSP and H-bond. The essential dynamics analysis with principal component analysis, dynamic cross-correlation matrix, and free energy landscape showed BRD2 binding of the screened compounds remained the stable conformation throughout molecular dynamics simulation. In addition, the pharmacokinetic and ADMET properties of the five compounds showed these compounds could be potential drug candidates. The findings of this study may contribute to further examinational validation and rational design of novel BRD2 inhibitors as anti-SARS-CoV-2 agents.

Polyphenolic flavonoid compounds act as the inhibitory potential of aggregation process: Implications for the prevention and therapeutics against FALS-associated D101G SOD1 mutant

Aggregation of proteins is a biological phenomenon caused by misfolded proteins. Human superoxide dismutase (hSOD1) misfolding and aggregation underlie the neurological illness amyotrophic lateral sclerosis (ALS). The most significant contributing factor to ALS is genetic point mutations in SOD1. particularly, D101G mutant is the most harmful because it significantly reduces the life expectancy of patients. Subsequently, the use of natural polyphenolic flavonoids is strongly recommended to reduce the amyloidogenic behavior of protopathic proteins. In this study, using computational parameters such as protein-ligand interaction and molecular dynamics (MD) simulation analyses, we are trying to identify a pharmacodynamically promising flavonoid compound that can effectively inhibit the pathogenic behavior of the D101G mutant. Epigallocatechin-gallate (EGCG), Hesperidin, Isorhamnetin, and Diosmetin were identified as potential leads in a preliminary screening of flavonoids to anti-amyloid action. The results of MD showed that the binding of flavonoids to D101G mutant caused changes in stability, hydrophobicity of protein, and flexibility, as well as significantly led to the restoration of lost hydrogen bonds. Secondary structure analysis showed that protein destabilization and the increased propensity of β-sheet caused by the mutation were restored to the wild-type state upon binding of flavonoids. Besides, to differentiate aggregation, we elucidated alterations in the free energy landscape (FEL) and dynamic cross-correlation matrix (DCCM) of WT-SOD1 and mutant (unbound /bound) states. Among flavonoids, Epigallocatechin-gallate and Hesperidin had the most therapeutic efficacy against the D101G mutant. Therefore, Epigallocatechin-gallate and Hesperidin promise considerable therapeutic potential to develop highly effective inhibitors in reducing fatal and irreversible ALS.

Novel computational and drug design strategies for inhibition of human papillomavirus-associated cervical cancer and DNA polymerase theta receptor by Apigenin derivatives

The present study deals with the advanced in-silico analyses of several Apigenin derivatives to explore human papillomavirus-associated cervical cancer and DNA polymerase theta inhibitor properties by molecular docking, molecular dynamics, QSAR, drug-likeness, PCA, a dynamic cross-correlation matrix and quantum calculation properties. The initial literature study revealed the potent antimicrobial and anticancer properties of Apigenin, prompting the selection of its potential derivatives to investigate their abilities as inhibitors of human papillomavirus-associated cervical cancer and DNA polymerase theta. In silico molecular docking was employed to streamline the findings, revealing promising energy-binding interactions between all Apigenin derivatives and the targeted proteins. Notably, Apigenin 4′-O-Rhamnoside and Apigenin-4′-Alpha-l-Rhamnoside demonstrated higher potency against the HPV45 oncoprotein E7 (PDB ID 2EWL), while Apigenin and Apigenin 5-O-Beta-d-Glucopyranoside exhibited significant binding energy against the L1 protein in humans. Similarly, a binding affinity range of − 7.5 kcal/mol to − 8.8 kcal/mol was achieved against DNA polymerase theta, indicating the potential of Apigenin derivatives to inhibit this enzyme (PDB ID 8E23). This finding was further validated through molecular dynamic simulation for 100 ns, analyzing parameters such as RMSD, RMSF, SASA, H-bond, and RoG profiles. The results demonstrated the stability of the selected compounds during the simulation. After passing the stability testing, the compounds underwent screening for ADMET, pharmacokinetics, and drug-likeness properties, fulfilling all the necessary criteria. QSAR, PCA, dynamic cross-correlation matrix, and quantum calculations were conducted, yielding satisfactory outcomes. Since this study utilized in silico computational approaches and obtained outstanding results, further validation is crucial. Therefore, additional wet-lab experiments should be conducted under in vivo and in vitro conditions to confirm the findings.

Identification of PPA1 inhibitor candidates for potential repurposing in cancer medicine.

Inorganic pyrophosphatase 1 (PPA1) is pivotal to cellular metabolism as it facilitates the hydrolysis of PPi-a by-product of various metabolic processes that influence cell growth and differentiation. Overexpression of PPA1 enzyme has been linked to diminished patient survival and was shown to influence tumor cell dynamics, thereby positioning it as a potential therapy target for a variety of cancers including colorectal cancer, diffuse large B-cell lymphoma, and lung adenocarcinoma. Despite this therapeutic promise, there are no known inhibitors of PPA1 as of today. In this study, we searched for potential PPA1 inhibitors using a molecular docking screen of 30 470 compounds with a history of clinical trials and/or US Food and Drug Administration approval. We specifically targeted the active pocket that coincides with the established catalytic domain. Our screen identified promising hits, which we further subjected to ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering. Subsequent molecular dynamics (MD) analyses were conducted on devazepide, quinotolast, and tarazepide-the three substances that successfully navigated all filters. MD analyses reinforced the stability of the protein-ligand complexes and confirmed ligand binding, as substantiated by our root mean square deviation, radius of gyration and secondary structures of proteins analyses. Furthermore, Molecular Mechanics Poisson-Boltzmann Surface Area calculations post-MD identified devazepide and quinotolast as showing higher binding affinities; being supported by principal component analysis, free energy landscape, and dynamic cross-correlation matrix results. Overall, our study reveals devazepide and quinotolast as potential candidates for PPA1 inhibition which could be considered for repurposing studies that need further experimental validation. These results not only reveal a potential for clinical repurposing for PPA1 inhibition but they also offer valuable insights into the development of future compounds for targeting the crucial PPA1 enzyme.

In vitro and in silico analysis of Solanum torvum fruit and methyl caffeate interaction with cholinesterases

Oxidative stress along with dysfunction in cholinergic neurotransmission primarily underlies cognitive impairment. A significant approach to mitigate cognitive dysfunction involves the inhibition of cholinesterases, namely acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Exploring the potential antioxidant and anticholinesterase effects of edible plants holds promise for their utilization as botanicals to enhance cognition. Solanum torvum fruit with vast biological properties are used as food. In the present study, butanolic extract of S. torvum fruits (BESTF) was prepared. Additionally, the study investigated into the properties of methyl caffeate (MC), a compound present in S. torvum, obtained in its pure form. In vitro antioxidant and anticholinesterases activity of BESTF and MC were determined. BESTF and MC showed potent antioxidant property. BESTF and MC dose-dependently inhibited AChE (IC50 values: 166.6 µg/ml and 680.6 µM, respectively) and BChE (IC50 values: 161.55 µg/ml and 413 µM, respectively). BESTF and MC inhibited AChE and BChE in competitive mode. Active site gorge of AChE/BChE was occupied by MC which formed interaction with amino acids present in catalytic site and PAS in in silico. Further, molecular dynamics simulations followed by free energy calculation, principal component analysis and dynamic cross-correlation matrix provided the compelling evidence that that MC maintained stable interactions during MD simulation with AChE and BChE. Collectively, the results from the present study underlines the cognitive-enhancing effect of BESTF and MC.

Antimalarial drug discovery against malaria parasites through haplopine modification: An advanced computational approach.

The widespread emergence of antimalarial drug resistance has created a major threat to public health. Malaria is a life-threatening infectious disease caused by Plasmodium spp., which includes Apicoplast DNA polymerase and Plasmodium falciparum cysteine protease falcipain-2. These components play a critical role in their life cycle and metabolic pathway, and are involved in the breakdown of erythrocyte hemoglobin in the host, making them promising targets for anti-malarial drug design. Our current study has been designed to explore the potential inhibitors from haplopinederivatives against these two targets using an in silico approach. A total of nine haplopine derivatives were used to perform molecular docking, and the results revealed that Ligands 03 and 05 showed strong binding affinity compared to the control compound atovaquone. Furthermore, these ligand-protein complexes underwent molecular dynamics simulations, and the results demonstrated that the complexes maintained strong stability in terms of RMSD (root mean square deviation), RMSF (root mean square fluctuation), and Rg (radius of gyration) over a 100 ns simulation period. Additionally, PCA (principal component analysis) analysis and the dynamic cross-correlation matrix showed positive outcomes for the protein-ligand complexes. Moreover, the compounds exhibited no violations of the Lipinski rule, and ADMET (absorption, distribution, metabolism, excretion, and toxicity) predictions yielded positive results without indicating any toxicity. Finally, density functional theory (DFT) and molecular electrostatic potential calculations were conducted, revealing that the mentioned derivatives exhibited better stability and outstanding performance. Overall, this computational approach suggests that these haplopine derivatives could serve as a potential source for developing new, effective antimalarial drugs to combat malaria. However, further in vitro or in vivo studies might be conducted to determine their actual effectiveness.

An in-silico receptor-pharmacophore based multistep molecular docking and simulation study to evaluate the inhibitory potentials against NS1 of DENV-2

DENV-2 strain is the most fatal and infectious of the five dengue virus serotypes. The non-structural protein NS1 encoded by its genome is the most significant protein required for viral pathogenesis and replication inside the host body. Thus, targeting the NS1 protein and designing an inhibitor to limit its stability and secretion is a propitious attempt in our fight against dengue. Four novel inhibitors are designed to target the conserved cysteine residues (C55, C313, C316, and C329) and glycosylation sites (N130 and N207) of the NS1 protein in an attempt to halt the spread of the dengue infection in the host body altogether. Numerous computer-aided drug designing techniques including molecular docking, molecular dynamics simulation, virtual screening, principal component analysis, and dynamic cross-correlation matrix were employed to determine the structural and functional activity of the NS1-inhibitor complexes. From our analysis, it was evident that the extent of structural and atomic level fluctuations of the ligand-bound protein exhibited a declining trend in contrast to unbound protein which was prominently noticeable through the RMSD, RMSF, Rg , and SASA graphs. The ADMET analysis of the four ligands revealed a promising pharmacokinetics and pharmacodynamic profile, along with good bioavailability and toxicity properties. The proposed drugs when bound to the targeted cavities resulted in stable conformations in comparison to their unbound state, implying they have good affinity promising effective drug action. Thus, they can be tested in vitro and used as potential anti-dengue drugs. Communicated by Ramaswamy H. Sarma

Target specific inhibition of West Nile virus envelope glycoprotein and methyltransferase using phytocompounds: an in silico strategy leveraging molecular docking and dynamics simulation.

Mosquitoes are the primary vector for West Nile virus, a flavivirus. The virus's ability to infiltrate and establish itself in increasing numbers of nations has made it a persistent threat to public health worldwide. Despite the widespread occurrence of this potentially fatal disease, no effective treatment options are currently on the market. As a result, there is an immediate need for the research and development of novel pharmaceuticals. To begin, molecular docking was performed on two possible West Nile virus target proteins using a panel of twelve natural chemicals, including Apigenin, Resveratrol, Hesperetin, Fungisterol, Lucidone, Ganoderic acid, Curcumin, Kaempferol, Cholic acid, Chlorogenic acid, Pinocembrin, and Sanguinarine. West Nile virus methyltransferase (PDB ID: 2OY0) binding affinities varied from -7.4 to -8.3 kcal/mol, whereas West Nile virus envelope glycoprotein affinities ranged from -6.2 to -8.1 kcal/mol (PDB ID: 2I69). Second, substances with larger molecular weights are less likely to be unhappy with the Lipinski rule. Hence, additional research was carried out without regard to molecular weight. In addition, compounds 01, 02, 03, 05, 06, 07, 08, 09, 10 and 11 are more soluble in water than compound 04 is. Besides, based on maximum binding affinity, best three compounds (Apigenin, Curcumin, and Ganoderic Acid) has been carried out molecular dynamic simulation (MDs) at 100 ns to determine their stability. The MDs data is also reported that these mentioned molecules are highly stable. Finally, advanced principal component analysis (PCA), dynamics cross-correlation matrices (DCCM) analysis, binding free energy and dynamic cross correlation matrix (DCCM) theoretical study is also included to established mentioned phytochemical as a potential drug candidate. Research has indicated that the aforementioned natural substances may be an effective tool in the battle against the dangerous West Nile virus. This study aims to locate a bioactive natural component that might be used as a pharmaceutical.

Binding Mechanism of CD47 with SIRPα Variants and Its Antibody: Elucidated by Molecular Dynamics Simulations.

The intricate complex system of the differentiation 47 (CD47) and the signal-regulatory protein alpha (SIRPα) cluster is a crucial target for cancer immunotherapy. Although the conformational state of the CD47-SIRPα complex has been revealed through crystallographic studies, further characterization is needed to fully understand the binding mechanism and to identify the hot spot residues involved. In this study, molecular dynamics (MD) simulations were carried out for the complexes of CD47 with two SIRPα variants (SIRPαv1, SIRPαv2) and the commercially available anti-CD47 monoclonal antibody (B6H12.2). The calculated binding free energy of CD47-B6H12.2 is lower than that of CD47-SIRPαv1 and CD47-SIRPαv2 in all the three simulations, indicating that CD47-B6H12.2 has a higher binding affinity than the other two complexes. Moreover, the dynamical cross-correlation matrix reveals that the CD47 protein shows more correlated motions when it binds to B6H12.2. Significant effects were observed in the energy and structural analyses of the residues (Glu35, Tyr37, Leu101, Thr102, Arg103) in the C strand and FG region of CD47 when it binds to the SIRPα variants. The critical residues (Leu30, Val33, Gln52, Lys53, Thr67, Arg69, Arg95, and Lys96) were identified in SIRPαv1 and SIRPαv2, which surround the distinctive groove regions formed by the B2C, C'D, DE, and FG loops. Moreover, the crucial groove structures of the SIRPα variants shape into obvious druggable sites. The C'D loops on the binding interfaces undergo notable dynamical changes throughout the simulation. For B6H12.2, the residues Tyr32LC, His92LC, Arg96LC, Tyr32HC, Thr52HC, Ser53HC, Ala101HC, and Gly102HC in its initial half of the light and heavy chains exhibit obvious energetic and structural impacts upon binding with CD47. The elucidation of the binding mechanism of SIRPαv1, SIRPαv2, and B6H12.2 with CD47 could provide novel perspectives for the development of inhibitors targeting CD47-SIRPα.