End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design

Oleanolic acid purified from the stem bark of Olax subscorpioidea Oliv. inhibits the function and catalysis of human 17β-hydroxysteroid dehydrogenase 1

Cancer is a leading cause of global death. Medicinal plants have gained increasing attention in cancer drug discovery. In this study, the stem bark extract of Olax subscorpioidea, which is used in ethnomedicine to treat cancer, was subjected to phytochemical investigation leading to the isolation of oleanolic acid (OA). The structure was elucidated by 1-dimensional and 2-dimensional nuclear magnetic resonance spectroscopic (NMR) data, and by comparing its data with previously reported data. Molecular docking was used to investigate the interactions of OA with nine selected cancer-related protein targets. OA docked well with human 17β-hydroxysteroid dehydrogenase type-1 (17βHSD1), caspase-3, and epidermal growth factor receptor tyrosine kinase (binding affinities: -9.8, -9.3, and -9.1 kcal/mol, respectively). OA is a triterpenoid compound with structural similarity to steroids. This similarity with the substrates of 17βHSD1 gives the inhibitor candidate an excellent opportunity to bind to 17βHSD1. The structural and functional dynamics of OA-17βHSD1 were investigated by molecular dynamics simulations at 240 ns. Molecular mechanics/Poisson-Boltzmann surface area (MMPBSA) studies showed that OA had a binding free energy that is comparable with that of vincristine (-52.76, and -63.56 kcal/mol, respectively). The average C-α root mean square of deviation (RMSD) value of OA (1.69 Å) compared with the unbound protein (2.01 Å) indicated its high stability at the protein's active site. The binding energy and stability at the active site of 17βHSD1 recorded in this study indicate that OA exhibited profound inhibitory potential. OA could be a good scaffold for developing new anti-breast cancer drugs.

Molecularly imprinted polymer-based sensors for identification volatile compounds in pharmaceutical products: in silico rational design

The present study aimed to strategically design a Molecularly Imprinted Polymer (MIP) with selective extraction capabilities for volatile compounds found in pork. These specific volatile compounds, such as 3-methyl-1-butanol, 1-nonanal, octanal, hexanal, 2-pentyl-furan, 1-penten-3-one, N-morpholinomethyl-isopropyl-sulfide, methyl butyrate, and (E,E)-2,4-decadienal, are primarily responsible for the distinctive aroma and flavor characteristics associated with pork. Molecular dynamics simulations were employed to investigate the stability of the pre-polymerization system, simulating the interactions between the volatile compounds as templates, 4-hydroxyethyl methacrylate (HEMA) as monomers, and ethylene glycol dimethacrylate (EGDMA) as crosslinkers. Computational simulations revealed that the optimal mole ratio of 1:4:20 for templates, monomers, and crosslinkers resulted in the most favorable functional radial distribution and exhibited the strongest interactions. To validate the computational findings, additional analyses were performed utilizing Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA), radial distribution function (RDF), and hydrogen bond (HBond) occupancy. The calculated binding free energy demonstrated that all template molecules were capable to bind with both the monomers and crosslinkers, including 1-penten-3-one and N-morpholinomethyl-isopropyl-sulfide displaying the strongest interactions, with values of -12,674 kJ/mol and -11,646 kJ/mol, respectively. The congruence between the results obtained from the molecular simulation analyses highlights the crucial role of molecular dynamics simulations in the study and development of MIP for the analysis of marker compounds present in pork.Communicated by Ramaswamy H. Sarma.

Naturally Occurring Plant-Based Anticancerous Candidates as Potential ERK2 Inhibitors: In-Silico Database Mining and Molecular Dynamics Simulations

The evolutionarily conserved extracellular signal-regulated kinase 2 (ERK2) is involved in regulating cellular signaling in both normal and pathological conditions. ERK2 expression is critical for human development, while hyperactivation is a major factor in tumor progression. Up to now, there have been no approved inhibitors that target ERK2, and as such, here we report on screening of a naturally occurring plant-based anticancerous compound-activity-target (NPACT) database for prospective ERK2 inhibitors. More than 1,500 phytochemicals were screened using in-silico molecular docking and molecular dynamics (MD) approaches. NPACT compounds with a docking score lower than a co-crystallized LHZ inhibitor (calc. -10.5 kcal/mol) were subjected to MD simulations. Binding energies (ΔGbinding) of inhibitor-ERK2 complexes over the MD course were estimated using an MM-GBSA approach. Based on MM-GBSA//100 ns MD simulations, the steroid zhankuic acid C (NPACT01034) demonstrated greater binding affinity against ERK2 protein than LHZ, with ΔGbinding values of -50.0 and -47.7 kcal/mol, respectively. Structural and energetical analyses throughout the MD course demonstrated stabilization of zhankuic acid C complexed with ERK2 protein. The anticipated ADMET properties of zhankuic acid C indicated minimal toxicity. Moreover, in-silico evaluation of fourteen ERK2 inhibitors in clinical trials demonstrated the higher binding affinity of zhankuic acid C towards ERK2 protein.

Computational Prediction of 3,5-Diaryl-1H-Pyrazole and spiropyrazolines derivatives as potential acetylcholinesterase inhibitors for alzheimer disease treatment by 3D-QSAR, molecular docking, molecular dynamics simulation, and ADME-Tox

The efficacy of 40 synthesized variants of 3,5-diaryl-1H-pyrazole and spiropyrazoline' derivatives as acetylcholinesterase inhibitors is verified using a quantitative three-dimensional structure-activity relationship (3D-QSAR) by comparative molecular field analysis (CoMFA) and molecular similarity index analysis (CoMSIA) models. In this research, different field models proved that CoMSIA/SE model is the best model with high predictive power compared to several models (Qved2 = O.65; R2 = 0.980; R2test = 0.727). Also, contour maps produced by CoMSIA/SE model have been employed to prove the key structural needs of the activity. Consequently, six new compounds have been generated. Among these compounds, M4 and M5 were the most active but remained toxic and had poor absorption capacities. While the M1, M2, M3 and M6 remained highly active while respecting ADMET's characteristics. Molecular docking results showed compound M2 better with acetylcholinesterase than compound 22. The interactions are classical hydrogen bonding with residues TYR:124, TYR:72, and SER:293, which play a critical role in the biological activity as AChE inhibitors. MD results confirmed the docking results and showed that compound M2 had satisfactory stability with (ΔGbinding = -151.225 KJ/mol) in the active site of AChE receptor compared with compound 22 (ΔGbinding = -133.375 KJ/mol). In addition, both compounds had good stability regarding RMSD, Rg, and RMSF. The previous results show that the newly designed compound M2 is more active in the active site of AChE receptor than compound 22.Communicated by Ramaswamy H. Sarma.

Microsecond Molecular Dynamics Simulation to Gain Insight Into the Binding of MRTX1133 and Trametinib With KRASG12D Mutant Protein for Drug Repurposing

The Kirsten Rat Sarcoma (KRAS) G12D mutant protein is a primary driver of pancreatic ductal adenocarcinoma, necessitating the identification of targeted drug molecules. Repurposing of drugs quickly finds new uses, speeding treatment development. This study employs microsecond molecular dynamics simulations to unveil the binding mechanisms of the FDA-approved MEK inhibitor trametinib with KRASG12D, providing insights for potential drug repurposing. The binding of trametinib was compared with clinical trial drug MRTX1133, which demonstrates exceptional activity against KRASG12D, for better understanding of interaction mechanism of trametinib with KRASG12D. The resulting stable MRTX1133-KRASG12D complex reduces root mean square deviation (RMSD) values, in Switch I and II domains, highlighting its potential for inhibiting KRASG12D. MRTX1133's robust interaction with Tyr64 and disruption of Tyr96-Tyr71-Arg68 network showcase its ability to mitigate the effects of the G12D mutation. In contrast, trametinib employs a distinctive binding mechanism involving P-loop, Switch I and II residues. Extended simulations to 1 μs reveal sustained network interactions with Tyr32, Thr58, and GDP, suggesting a role of trametinib in maintaining KRASG12D in an inactive state and impede the further cell signaling. The decomposition binding free energy values illustrate amino acids' contributions to binding energy, elucidating ligand-protein interactions and molecular stability. The machine learning approach reveals that van der Waals interactions among the residues play vital role in complex stability and the potential amino acids involved in drug-receptor interactions of each complex. These details provide a molecular-level understanding of drug binding mechanisms, offering essential knowledge for further drug repurposing and potential drug discovery.

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

The human transmembrane protein Transferrin Receptor-1 is regarded as a promising target for the systemic delivery of therapeutic agents, particularly of nucleic acid therapeutics, such as short double stranded RNAs. This ubiquitous receptor is involved in cellular iron uptake, keeping intracellular homeostasis. It is overexpressed in multiple cancer cell types and is internalized via clathrin-mediated endocytosis. In previous studies, a human transferrin receptor-1 RNA aptamer, identified as TR14 ST1-3, was shown to be capable of effectively internalizing into cells in culture and to deliver small, double stranded RNAs in vitro and in vivo, via systemic administration. To understand, at the molecular level, the aptamer binding to the receptor and the impact of conjugation with the therapeutic RNA, a multi-level in silico protocol was employed, including protein-aptamer docking, molecular dynamics simulations and free energy calculations. The competition for the binding pocket, between the aptamer and the natural ligand human Transferrin, was also evaluated. The results show that the aptamer binds to the same region as Transferrin, with residues from the helical domain showing a critical role. Moreover, the conjugation to the therapeutic RNA, was shown not to affect aptamer binding. Overall, this study provides an atomic-level understanding of aptamer association to human Transferrin Receptor-1 and of its conjugation with a short model-therapeutic RNA, providing also important clues for futures studies aiming to deliver other oligonucleotide-based therapeutics via Transferrin Receptor.

A comprehensive study of experimental and theoretical characterization and in silico toxicity analysis of new molecules

In this study, for the first time in the literature, a 2-(3-methoxyphenylamino)-2-oxoethyl acrylate (3MPAEA) molecule was synthesized in two steps, and a 2-chloro-N-(3-methoxyphenyl)acetamide (m-acetamide) was obtained in the first step. Experimental results were obtained using FTIR, 1H, and 13C NMR spectroscopy methods for m-acetamide and 3MPAEA compounds created in the laboratory environment and compared with theoretical results. Band gap (BG) energy, chemical hardness, electronegativity, chemical potential, and electrophilicity index were calculated. With vibration spectroscopic analysis, atom-molecule vibrations of the theoretical and experimental peaks of the spectrum were observed. The locations of C and H atoms were determined by nuclear magnetic resonance spectroscopy. The green, blue, and red regions of the potential energy map (MEP) map were examined. Some observed that the energy thermal, heat capacity, and entropy graphs increased in direct proportion to increasing the temperature in Kelvin, which is known as thermochemistry. The changes in the rotation, translation, and vibration of the molecule as its temperature increased were examined. When the thermochemistry surface map was examined, some observed that the temperature was high in the middle binding site of the molecules. Covalent interactions were graphed using the non-covalent interactions (NCIs) calculation method. In silico toxicity studies were carried out for m-acetamide and 3MPAEA molecules: fathead minnow LC50 (96 h), Daphnia magna LC50 (48 h), Tetrahymena pyriformis IGC50 (48 h), oral rat LD50, water solubility, bioconcentration factor, developmental toxicity, mutation, normal boiling point, flash point, melting point, density, thermal conductivity, viscosity, vapor pressure, etc. parameters were investigated.

Multiple Functional Protein-Protein Interaction Interfaces Allosterically Regulate ATP-Binding in Cyclin-Dependent Kinase-1

The ATP-dependent phosphorylation activity of cyclin-dependent kinase 1 (CDK1), an essential enzyme for cell cycle progression, is regulated by interactions with Cyclin-B, substrate, and Cks proteins. We have recently shown that active site acetylation in CDK1 abrogated binding to Cyclin-B which posits an intriguing long-range communication between the catalytic site and the protein-protein interaction (PPI) interface. Now, we demonstrate a general allosteric link between the CDK1 active site and all three of its PPI interfaces through atomistic molecular dynamics (MD) simulations. Specifically, we examined ATP binding free energies to CDK1 in native nonacetylated (K33wt) and acetylated (K33Ac) forms as well as the acetyl-mimic K33Q and the acetyl-null K33R mutant forms, which are accessible in vitro. In agreement with experiments, ATP binding is stronger in K33wt relative to the other three perturbed states. Free energy decomposition reveals, in addition to expected local changes, significant and selective nonlocal entropic responses to ATP binding/perturbation of K33 from theαC -helix, activation loop (A-loop), andαG -α H segments in CDK1 which interface with Cyclin-B, substrate, and Cks proteins, respectively. Statistical analysis reveals that while entropic responses of protein segments to active site perturbations are on average correlated with their dynamical changes, such correlations are lost in about 9%-48% of the dataset depending on the segment. Besides proving the bi-directional communication between the active site and the CDK1:Cyclin-B interface, our study uncovers a hitherto unknown mode of ATP binding regulation by multiple PPI interfaces in CDK1.

New oxomethacrylate and acetamide: synthesis, characterization, and their computational approaches: molecular docking, molecular dynamics, and ADME analyses

The compounds 2-chloro-N-(3-methoxyphenyl)acetamide (m-acetamide) and 2-(3-methoxyphenylamino)-2-oxoethyl methacrylate (3MPAEMA) were synthesized in this study for the first time in the literature. FTIR, 1H, and 13C NMR spectroscopic techniques were used to characterize it. Subsequently, computational techniques were used to assess various ADME factors, such as drug-likeness properties, bioavailability score, and adherence to Lipinski's rule. Finally, molecular docking experiments were conducted with the human topoisomerase α2 (TOP2A) protein to verify and validate the reliability and stability of the docking procedure. The results of the docking scores, which quantify binding affinity, indicated that these derivatives exhibited a stronger affinity for TOP2A.

Novel benzothiazole derivatives target the Gac/Rsm two-component system as antibacterial synergists against Pseudomonas aeruginosa infections

The management of antibiotic-resistant, bacterial biofilm infections in skin wounds poses an increasingly challenging clinical scenario. Pseudomonas aeruginosa infection is difficult to eradicate because of biofilm formation and antibiotic resistance. In this study, we identified a new benzothiazole derivative compound, SN12 (IC50 = 43.3 nmol/L), demonstrating remarkable biofilm inhibition at nanomolar concentrations in vitro. In further activity assays and mechanistic studies, we formulated an unconventional strategy for combating P. aeruginosa-derived infections by targeting the two-component (Gac/Rsm) system. Furthermore, SN12 slowed the development of ciprofloxacin and tobramycin resistance. By using murine skin wound infection models, we observed that SN12 significantly augmented the antibacterial effects of three widely used antibiotics-tobramycin (100-fold), vancomycin (200-fold), and ciprofloxacin (1000-fold)-compared with single-dose antibiotic treatments for P. aeruginosa infection in vivo. The findings of this study suggest the potential of SN12 as a promising antibacterial synergist, highlighting the effectiveness of targeting the two-component system in treating challenging bacterial biofilm infections in humans.

A head-to-head comparison of MM/PBSA and MM/GBSA in predicting binding affinities for the CB1 cannabinoid ligands

CONTEXT

The substantial increase in the number of active and inactive-state CB1 receptor experimental structures has provided opportunities for CB1 drug discovery using various structure-based drug design methods, including the popular end-point methods for predicting binding free energies-Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA). In this study, we have therefore evaluated the performance of MM/PBSA and MM/GBSA in calculating binding free energies for CB1 receptor. Additionally, with both MM/PBSA and MM/GBSA being known for their highly individualized performance, we have evaluated the effects of various simulation parameters including the use of energy minimized structures, choice of solute dielectric constant, inclusion of entropy, and the effects of the five GB models. Generally, MM/GBSA provided higher correlations than MM/PBSA (rMM/GBSA = 0.433 - 0.652 vs. rMM/PBSA = 0.100 - 0.486) regardless of the simulation parameters, while also offering faster calculations. Improved correlations were observed with the use of molecular dynamics ensembles compared with energy minimized structures and larger solute dielectric constants. Incorporation of entropic terms led to unfavorable results for both MM/PBSA and MM/GBSA for a majority of the dataset, while the evaluation of the various GB models exerted a varying effect on both the datasets. The findings obtained in this study demonstrate the utility of MM/PBSA and MM/GBSA in predicting binding free energies for the CB1 receptor, hence providing a useful benchmark for their applicability in the endocannabinoid system as well as other G protein-coupled receptors.

METHODS

The study utilized the docked dataset (Induced Fit Docking with Glide XP scoring function) from Loo et al., consisting of 46 ligands-23 agonists and 23 antagonists. The equilibrated structures from Loo et al. were subjected to 30 ns production simulations using GROMACS 2018 at 300 K and 1 atm with the velocity rescaling thermostat and the Parinello-Rahman barostat. AMBER ff99SB*-ILDN was used for the proteins, General Amber Force Field (GAFF) was used for the ligands, and Slipids parameters were used for lipids. MM/PBSA and MM/GBSA binding free energies were then calculated using gmx_MMPBSA. The solute dielectric constant was varied between 1, 2, and 4 to study the effect of different solute dielectric constants on the performance of MM/PB(GB)SA. The effect of entropy on MM/PB(GB)SA binding free energies was evaluated using the interaction entropy module implemented in gmx_MMPBSA. Five GB models, GBHCT, GBOBC1, GBOBC2, GBNeck, and GBNeck2, were evaluated to study the effect of the choice of GB models in the performance of MM/GBSA. Pearson correlation coefficients were used to measure the correlation between experimental and predicted binding free energies.

Discovery of Cyclic Peptide Inhibitors Targeted on TNFα-TNFR1 from Computational Design and Bioactivity Verification

Activating tumor necrosis factor receptor 1 (TNFR1) with tumor necrosis factor alpha (TNFα) is one of the key pathological mechanisms resulting in the exacerbation of rheumatoid arthritis (RA) immune response. Despite various types of drugs being available for the treatment of RA, a series of shortcomings still limits their application. Therefore, developing novel peptide drugs that target TNFα-TNFR1 interaction is expected to expand therapeutic drug options. In this study, the detailed interaction mechanism between TNFα and TNFR1 was elucidated, based on which, a series of linear peptides were initially designed. To overcome its large conformational flexibility, two different head-to-tail cyclization strategies were adopted by adding a proline-glycine (GP) or cysteine-cysteine (CC) to form an amide or disulfide bond between the N-C terminal. The results indicate that two cyclic peptides, R1_CC4 and α_CC8, exhibit the strongest binding free energies. α_CC8 was selected for further optimization using virtual mutations through in vitro activity and toxicity experiments due to its optimal biological activity. The L16R mutant was screened, and its binding affinity to TNFR1 was validated using ELISA assays. This study designed a novel cyclic peptide structure with potential anti-inflammatory properties, possibly bringing an additional choice for the treatment of RA in the future.

Alpha-1 antitrypsin inhibits pertussis toxin

Pertussis (whooping cough) is a vaccine-preventable but re-emerging, highly infectious respiratory disease caused by Bordetella pertussis. There are currently no effective treatments for pertussis, complicating care for nonvaccinated individuals, especially newborns. Disease manifestations are predominantly caused by pertussis toxin (PT), a pivotal virulence factor classified as an ADP-ribosylating AB-type protein toxin. In this work, an unbiased approach using peptide libraries, bioassay-guided fractionation and mass spectrometry revealed α1-antitrypsin (α1AT) as a potent PT inhibitor. Biochemistry-, cell culture-, and molecular modeling-based in vitro experimentation demonstrated that the α1AT mode of action is based on blocking PT-binding to the host target cell surface. In the infant mouse model of severe pertussis, α1AT expression was reduced upon infection. Further, systemic administration of α1AT significantly reduced B. pertussis-induced leukocytosis, which is a hallmark of infant infection and major risk factor for fatal pertussis. Taken together our data demonstrates that α1AT is a novel PT inhibitor and that further evaluation and development of α1AT as a therapeutic agent for pertussis is warranted. Importantly, purified α1AT is already in use clinically as an intravenous augmentation therapy for those with genetic α1AT deficiency and could be repurposed to clinical management of pertussis.

Drug repurposing: An antidiabetic drug Ipragliflozin as Mycobacterium tuberculosis sirtuin-like protein inhibitor that synergizes with anti-tuberculosis drug isoniazid

The surge of drug-resistant Mycobacterium tuberculosis (DR-TB) impedes the World Health Organization's efforts in ending TB and calls for new therapeutic formulations. M. tuberculosis sirtuin-like protein Rv1151c is a bifunctional enzyme with both deacetylation and desuccinylation activities, which plays an important role in M. tuberculosis drug resistance and stress responses. Thus, it appears to be a promising target for the development of new TB therapeutics. In this study, we screened 31,057 ligand compounds from seven compound libraries in silico to identify inhibitors of Rv1151c. Ipragliflozin can bind to Rv1151c and interact stably. Ipragliflozin can change the acylation level of M. tuberculosis by inhibiting Rv1151c and effectively inhibit the growth of M. tuberculosis H37Rv and M. smegmatis. It can potentiate the first-front anti-TB drug isoniazid. As an antidiabetic drug, Ipragliflozin can be potentially included in the regimen to treat diabetes-tuberculosis comorbidity.

Applying Spatiotemporal Modeling of Cell Dynamics to Accelerate Drug Development

Cells act as physical computational programs that utilize input signals to orchestrate molecule-level protein-protein interactions (PPIs), generating and responding to forces, ultimately shaping all of the physiological and pathophysiological behaviors. Genome editing and molecule drugs targeting PPIs hold great promise for the treatments of diseases. Linking genes and molecular drugs with protein-performed cellular behaviors is a key yet challenging issue due to the wide range of spatial and temporal scales involved. Building predictive spatiotemporal modeling systems that can describe the dynamic behaviors of cells intervened by genome editing and molecular drugs at the intersection of biology, chemistry, physics, and computer science will greatly accelerate pharmaceutical advances. Here, we review the mechanical roles of cytoskeletal proteins in orchestrating cellular behaviors alongside significant advancements in biophysical modeling while also addressing the limitations in these models. Then, by integrating generative artificial intelligence (AI) with spatiotemporal multiscale biophysical modeling, we propose a computational pipeline for developing virtual cells, which can simulate and evaluate the therapeutic effects of drugs and genome editing technologies on various cell dynamic behaviors and could have broad biomedical applications. Such virtual cell modeling systems might revolutionize modern biomedical engineering by moving most of the painstaking wet-laboratory effort to computer simulations, substantially saving time and alleviating the financial burden for pharmaceutical industries.

Molecular Mechanism-Driven Discovery of Novel Small Molecule Inhibitors against Drug-Resistant SARS-CoV-2 Mpro Variants

Under the selective pressure of nirmatrelvir, a peptidomimetic covalent drug targeting SARS-CoV-2 Mpro, various drug-resistant mutations on Mpro have been acquired in vitro. Among the mutations, L50F and E166V, along with the combination of L50F and E166V, are particularly representative and pose considerable obstacles to the effective treatment of COVID-19. Our previous study identified NMI-001 and NMI-002 as novel nonpeptide inhibitors that target SARS-CoV-2 Mpro, possessing unique scaffolds and binding modes different from those of nirmatrelvir. In view of these findings, we proposed a drug design strategy aimed at rapidly identifying inhibitors that can combat mutation-induced drug resistance. Initially, molecular dynamics (MD) simulation was employed to investigate the binding mechanisms of NMI-001 and NMI-002 against the three drug-resistant mutants (Mpro_L50F, Mpro_E166V, and Mpro_L50F+E166V). Then, we conducted two phases of high-throughput virtual screening. In the first phase, NMI-001 served as a template to perform scaffold hopping-based similarity search in a library of 15,742,661 compounds. In the second phase, 968 compounds exhibiting similarity to NMI-001 were evaluated via molecular docking and MD simulations. Six compounds that may be effective against at least one mutant were identified, and five compounds were procured for conducting in vitro assays. Finally, the compound Z1557501297 (NMI-003) exhibiting inhibitory effects against the E166V (IC50 = 27.81 ± 2.65 μM) and L50F+E166V (IC50 = 8.78 ± 0.74 μM) mutants was discovered. The binding modes referring to NMI-003-Mpro_E166V and NMI-003-Mpro_L50F+E166V were further elucidated at the atomic level. In summary, NMI-003 reported herein is the first compound with activity against E166V and L50F+E166V, which provides a good starting point to design novel antiviral drugs for the treatment of drug-resistant SARS-CoV-2.

Dynamic interactions and inhibitory mechanisms of Artemisia annua terpenoids with carbonic anhydrase IX

This study evaluates the inhibitory potential of terpenoids isolated from Artemisia annua against carbonic anhydrase IX (CAIX), a crucial enzyme overexpressed in hypoxic tumor environments. Employing a multidisciplinary approach, we utilized in vitro assays, enzyme kinetics, molecular docking, and molecular dynamics (MD) simulations to comprehensively assess the efficacy of these compounds. Among the terpenoids tested, manool emerged as the most potent inhibitor, exhibiting the lowest IC50 value of 160.2 ± 15.2 nM. This was followed by labda-8(17),12-diene-15,16-dial with an IC50 of 297.9 ± 8.84 nM. Enzyme kinetics revealed a mixed inhibition mode for both compounds. Molecular docking aligned well with in vitro data, showing extensive polar and hydrophobic interactions within the CAIX binding site. Further insights were gained through 300 ns MD simulations, which highlighted the dynamic interactions and stability of these complexes. Manool demonstrated the most significant stabilization of CAIX, as evidenced by favorable RMSD, Rg, SASA profiles, and the strongest hydrogen bonding interactions. Additionally, MM/PBSA calculations confirmed manool's superior binding affinity. These findings underscore the therapeutic potential of manool as a potent CAIX inhibitor, providing a foundation for the development of effective anticancer agents targeting hypoxic tumor environments. ADMET analysis revealed favorable pharmacokinetic profiles for the terpenoids, with manool demonstrating high lipophilicity and BBB permeability, though potential CYP-mediated interactions were noted.

Structural basis for hepatitis B virus restriction by a viral receptor homologue

Macaque restricts hepatitis B virus (HBV) infection because its receptor homologue, NTCP (mNTCP), cannot bind preS1 on viral surface. To reveal how mNTCP loses the viral receptor function, we here solve the cryo-electron microscopy structure of mNTCP. Superposing on the human NTCP (hNTCP)-preS1 complex structure shows that Arg158 of mNTCP causes steric clash to prevent preS1 from embedding onto the bile acid tunnel of NTCP. Cell-based mutation analysis confirms that only Gly158 permitted preS1 binding, in contrast to robust bile acid transport among mutations. As the second determinant, Asn86 on the extracellular surface of mNTCP shows less capacity to restrain preS1 from dynamic fluctuation than Lys86 of hNTCP, resulting in unstable preS1 binding. Additionally, presence of long-chain conjugated-bile acids in the tunnel induces steric hindrance with preS1 through their tailed-chain. This study presents structural basis in which multiple sites in mNTCP constitute a molecular barrier to strictly restrict HBV.

Design, Synthesis, and Pharmacological Evaluation of Dual FXR-LIFR Modulators for the Treatment of Liver Fibrosis

Although multiple approaches have been suggested, treating mild-to-severe fibrosis in the context of metabolic dysfunction associated with liver disease (MASLD) remains a challenging area in drug discovery. Pathogenesis of liver fibrosis is multifactorial, and pathogenic mechanisms are deeply intertwined; thus, it is well accepted that future treatment requires the development of multitarget modulators. Harnessing the 3,4,5-trisubstituted isoxazole scaffold, previously described as a key moiety in Farnesoid X receptor (FXR) agonism, herein we report the discovery of a novel class of hybrid molecules endowed with dual activity toward FXR and the leukemia inhibitory factor receptor (LIFR). Up to 27 new derivatives were designed and synthesized. The pharmacological characterization of this series resulted in the identification of 3a as a potent FXR agonist and LIFR antagonist with excellent ADME properties. In vitro and in vivo characterization identified compound 3a as the first-in-class hybrid LIFR inhibitor and FXR agonist that protects against the development of acute liver fibrosis and inflammation.

In Silico Modeling of Myelin Oligodendrocyte Glycoprotein Disulfide Bond Reduction by Phosphine-Borane Complexes

BACKGROUND

Neurodegenerative diseases can cause vision loss by damaging retinal ganglion cells in the optic nerve. Novel phosphine-borane compounds (PBs) can protect these cells from oxidative stress via the reduction of disulfide bonds. However, the specific targets of these compounds are unknown. Proteomic evidence suggests that myelin oligodendrocyte glycoprotein (MOG) is a potential target. MOG is of significant interest due to its role in anti-MOG optic neuritis syndrome.

METHODS

We used in silico modeling to explore the structural consequences of cleaving the extracellular domain MOG disulfide bond, both in isolation and in complex with anti-MOG antibodies. The potential binding of PBs to this bond was examined using molecular docking.

RESULTS

Cleaving the disulfide bond of MOG altered the structure of MOG dimers and reduced their energetic favorability by 46.13 kcal/mol. The energy profiles of anti-MOG antibody complexes were less favorable when the disulfide bond of MOG was reduced in the monomeric state by 55.21 kcal/mol, but the reverse was true in the dimeric state. PBs exhibited reducing capabilities with the MOG extracellular disulfide bond, with this best-scoring compound binding with an energy of -28.54 kcal/mol to the MOG monomer and -24.97 kcal/mol to the MOG dimer.

CONCLUSIONS

These findings suggest that PBs can affect the structure of MOG dimers and the formation of antibody complexes by reducing the MOG disulfide bond. Structural changes in MOG could have implications for neurodegenerative diseases and anti-MOG syndrome.

Advances, opportunities, and challenges in methods for interrogating the structure activity relationships of natural products

Time span in literature: 1985-early 2024Natural products play a key role in drug discovery, both as a direct source of drugs and as a starting point for the development of synthetic compounds. Most natural products are not suitable to be used as drugs without further modification due to insufficient activity or poor pharmacokinetic properties. Choosing what modifications to make requires an understanding of the compound's structure-activity relationships. Use of structure-activity relationships is commonplace and essential in medicinal chemistry campaigns applied to human-designed synthetic compounds. Structure-activity relationships have also been used to improve the properties of natural products, but several challenges still limit these efforts. Here, we review methods for studying the structure-activity relationships of natural products and their limitations. Specifically, we will discuss how synthesis, including total synthesis, late-stage derivatization, chemoenzymatic synthetic pathways, and engineering and genome mining of biosynthetic pathways can be used to produce natural product analogs and discuss the challenges of each of these approaches. Finally, we will discuss computational methods including machine learning methods for analyzing the relationship between biosynthetic genes and product activity, computer aided drug design techniques, and interpretable artificial intelligence approaches towards elucidating structure-activity relationships from models trained to predict bioactivity from chemical structure. Our focus will be on these latter topics as their applications for natural products have not been extensively reviewed. We suggest that these methods are all complementary to each other, and that only collaborative efforts using a combination of these techniques will result in a full understanding of the structure-activity relationships of natural products.

Mechanistic insights into carbonic anhydrase IX inhibition by coumarins from Calendula officinalis: in vitro and in silico approaches

Given the critical role of carbonic anhydrase IX (CA IX) in various pathological conditions, there is a significant demand for new inhibitors to enhance patient outcomes and clinical management. In this study, we examined the inhibitory effectiveness of five coumarins derived from Calendula officinalis against CA IX using in vitro assays and computational modeling. Among the coumarins tested, xeroboside and isobaisseoside were identified as the most potent inhibitors. Kinetic studies indicated that xeroboside and isobaisseoside exhibit a mixed inhibition mode. Molecular docking analysis showed that the tested coumarins exhibit binding affinities and extensive polar interactions with CA IX. These coumarins demonstrated significant hydrophobic interactions and occupied the same binding site as acetazolamide (AAZ). Molecular dynamics (MD) indicated that xeroboside and isobaisseoside exhibited consistent trajectories and notable energy stabilization during their interaction with CA IX. MM/PBSA calculations showed that xeroboside displayed the lowest binding free energy (-27.26 ± 2.48 kJ mol-1). Potential Energy Landscape (PEL) analysis revealed distinct and stable conformational states for the CA IX-ligand complexes, with xeroboside exhibiting the most stable and lowest energy configuration. These computational findings are consistent with the experimental results, highlighting the potential efficacy of xeroboside and isobaisseoside as CA IX inhibitors. In conclusion, Calendula officinalis-derived coumarins are promising candidates as effective CA IX inhibitors.

Identification of Novel PPARγ Partial Agonists Based on Virtual Screening Strategy: In Silico and In Vitro Experimental Validation

Thiazolidinediones (TZDs) including rosiglitazone and pioglitazone function as peroxisome proliferator-activated receptor gamma (PPARγ) full agonists, which have been known as a class to be among the most effective drugs for the treatment of type 2 diabetes mellitus (T2DM). However, side effects of TZDs such as fluid retention and weight gain are associated with their full agonistic activities toward PPARγ induced by the AF-2 helix-involved "locked" mechanism. Thereby, this study aimed to obtain novel PPARγ partial agonists without direct interaction with the AF-2 helix. Through performing virtual screening of the Targetmol L6000 Natural Product Library and utilizing molecular dynamics (MD) simulation, as well as molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis, four compounds including tubuloside b, podophyllotoxone, endomorphin 1 and paliperidone were identified as potential PPARγ partial agonists. An in vitro TR-FRET competitive binding assay showed podophyllotoxone displayed the optimal binding affinity toward PPARγ among the screened compounds, exhibiting IC50 and ki values of 27.43 µM and 9.86 µM, respectively. Further cell-based transcription assays were conducted and demonstrated podophyllotoxone's weak agonistic activity against PPARγ compared to that of the PPARγ full agonist rosiglitazone. These results collectively demonstrated that podophyllotoxone could serve as a PPARγ partial agonist and might provide a novel candidate for the treatment of various diseases such as T2DM.

Driving forces and large scale affinity calculations for piperine/γ-cyclodxetrin complexes: Mechanistic insights from umbrella sampling simulation and DFT calculations

Piperine (PiP), a bioactive molecule, exhibits numerous health benefits and is frequently employed as a co-delivery agent with various phytomedicines (e.g., curcumin) to enhance their bioavailability. This is attributed to PiP's inhibitory activity against drug-metabolizing proteins, notably CYP3A4. Nevertheless, PiP encounters solubility challenges addressed in this study using cyclodextrins (CDs). Specifically, γ-CD and its derivatives, Hydroxypropyl-γ-CD (HP-γ-CD), and Octakis (6-O-sulfo)-γ-CD (Octakis-S-γ-CD), were employed to form supramolecular complexes with PiP. The conformational space of the complexes was assessed through 1 μs molecular dynamics simulations and umbrella sampling. Additionally, quantum mechanical calculations using wB97X-D dispersion-corrected DFT functional and 6-311 + G(d,p) basis set were conducted on the complexes to examine the thermodynamics and kinetic stability. Results indicated that Octakis-S-γ-CD exhibits superior host capabilities for PiP, with the most favorable complexation energy (-457.05 kJ/mol), followed by HP-γ-CD (-249.16 kJ/mol). Furthermore, two conformations of the Octakis-S-γ-CD/PiP complex were explored to elucidate the optimal binding orientation of PiP within the binding pocket of Octakis-S-γ-CD. Supramolecular chemistry relies significantly on non-covalent interactions. Therefore, our investigation extensively explores the critical atoms involved in these interactions, elucidating the influence of substituted groups on the stability of inclusion complexes. This comprehensive analysis contributes to emphasizing the γ-CD derivatives with improved host capacity.

Characterizing the Cooperative Effect of PROTAC Systems with End-Point Binding Free Energy Calculation

Proteolytic targeting chimeras (PROTACs), as an emerging type of drug, function by proximity-based modalities that narrow the distance between a target protein and the E3 ubiquitin ligase to facilitate the ubiquitination labeling of the target protein for degradation. Although it is evidenced that the cooperativity of the PROTAC ternary interaction is one of the key factors affecting the degradation rate of a target protein, PROTAC design utilizing this indicator is still challenging because of the complicated/flexible interactions in a target-PROTAC-E3 ternary system. Therefore, developing reliable and practicable computational methods is of great interest for PROTAC design. Hence, in this study, we investigate the feasibility of using the end-point binding free energy calculation method, represented by molecular mechanics/Poisson-Boltzmann (generalized-Born) surface area (MM/PB(GB)SA), for characterizing cooperativity (including the stabilization and hook effects) of the PROTAC systems. The result shows that MM/GBSA is a good predictor in characterizing these effects under a relatively long molecular dynamics adjustment (50-100 ns) and low dielectric constant (εin = 1), with the Pearson correlation coefficient (rp) > 0.5 and 0.6 for the stabilization and hook effect, respectively. This study provides a feasible strategy for characterizing the cooperativity of the PROTAC systems, facilitating the rational design of PROTAC molecules.

Alchemical Transformations and Beyond: Recent Advances and Real-World Applications of Free Energy Calculations in Drug Discovery

Computational methods constitute efficient strategies for screening and optimizing potential drug molecules. A critical factor in this process is the binding affinity between candidate molecules and targets, quantified as binding free energy. Among various estimation methods, alchemical transformation methods stand out for their theoretical rigor. Despite challenges in force field accuracy and sampling efficiency, advancements in algorithms, software, and hardware have increased the application of free energy perturbation (FEP) calculations in the pharmaceutical industry. Here, we review the practical applications of FEP in drug discovery projects since 2018, covering both ligand-centric and residue-centric transformations. We show that relative binding free energy calculations have steadily achieved chemical accuracy in real-world applications. In addition, we discuss alternative physics-based simulation methods and the incorporation of deep learning into free energy calculations.

Based on Molecular Docking, Molecular Dynamics Simulation and MM/PB(GB)SA to Study Potential Inhibitors of PRRSV-Nsp4

Porcine reproductive and respiratory syndrome (PRRS) is one of the most serious infectious immunosuppressive diseases in the world. The nonstructural protein Nsp4 can be used as an ideal target for anti-PRRSV replication inhibitors. However, little is known about potential inhibitors that target Nsp4 to affect PRRSV replication. The purpose of this study was to screen potential natural inhibitors that affect PRRSV replication by inhibiting Nsp4. Five compounds with strong binding affinity to Nsp4 were selected by structure-based molecular docking method. The complexes of naringin dihydrochalcone (NDC), agathisflavone (AGT), and amentoflavone (AMF) with Nsp4 were stable throughout the molecular dynamics simulation. According to MM/PBSA analysis, the free energies of binding of NDC, AGT, and AMF to Nsp4 were less than-30 Kcal/mol. In conclusion, these three compounds are worthy of further investigation as novel inhibitors of PRRSV. This study provides a theoretical basis for the development of anti-PRRSV natural drugs.

A Method for Treating Significant Conformational Changes in Alchemical Free Energy Simulations of Protein-Ligand Binding

Relative binding free energy (RBFE) simulation is a rigorous approach to the calculation of quantitatively accurate binding free energy values for protein-ligand binding in which a reference binder is gradually converted to a target binder through alchemical transformation during the simulation. The success of such simulations relies on being able to accurately sample the correct conformational phase space for each alchemical state; however, this becomes a challenge when a significant conformation change occurs between the reference and target binder-receptor complexes. Increasing the simulation time and using enhanced sampling methods can be helpful, but effects can be limited, especially when the free energy barrier between conformations is high or when the correct target complex conformation is difficult to find and maintain. Current RBFE protocols seed the reference complex structure into every alchemical window of the simulation. In our study, we describe an improved protocol in which the reference structure is seeded into the first half of the alchemical windows, and the target structure is seeded into the second half of the alchemical windows. By applying information about the relevant correct end point conformations to different simulation windows from the beginning, the need for large barrier crossings or simulation prediction of the correct structures during an alchemical simulation is in many cases obviated. In the diverse cases we examine below, the simulations yielded free energy predictions that are satisfactory as compared to experiment and superior to running the simulations utilizing the conventional protocol. The method is straightforward to implement for publicly available FEP workflows.

Computational exploration of acefylline derivatives as MAO-B inhibitors for Parkinson's disease: insights from molecular docking, DFT, ADMET, and molecular dynamics approaches

Monoamine oxidase B (MAO-B) plays a pivotal role in the deamination process of monoamines, encompassing crucial neurotransmitters like dopamine and norepinephrine. The heightened interest in MAO-B inhibitors emerged after the revelation that this enzyme could potentially catalyze the formation of neurotoxic compounds from endogenous and exogenous sources. Computational screening methodologies serve as valuable tools in the quest for novel inhibitors, enhancing the efficiency of this pursuit. In this study, 43 acefylline derivatives were docked against the MAO-B enzyme for their chemotherapeutic potential and binding affinities that yielded GOLD fitness scores ranging from 33.21 to 75.22. Among them, five acefylline derivatives, namely, MAO-B14, MAO-B15, MAO-B16, MAO-B20, and MAO-B21, displayed binding affinities comparable to the both standards istradefylline and safinamide. These derivatives exhibited hydrogen-bonding interactions with key amino acids Phe167 and Ile197/198, suggesting their strong potential as MAO-B inhibitors. Finally, molecular dynamics (MD) simulations were conducted to evaluate the stability of the examined acefylline derivatives over time. The simulations demonstrated that among the examined acefylline derivatives and standards, MAO-B21 stands out as the most stable candidate. Density functional theory (DFT) studies were also performed to optimize the geometries of the ligands, and molecular docking was conducted to predict the orientations of the ligands within the binding cavity of the protein and evaluate their molecular interactions. These results were also validated by simulation-based binding free energies via the molecular mechanics energies combined with generalized Born and surface area solvation (MM-GBSA) method. However, it is necessary to conduct in vitro and in vivo experiments to confirm and validate these findings in future studies.

Activation of polycystin-1 signaling by binding of stalk-derived peptide agonists

Polycystin-1 (PC1) is the protein product of the PKD1 gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.

A Small Molecule Impedes the Aβ1-42 Tetramer Neurotoxicity by Preserving Membrane Integrity: Microsecond Multiscale Simulations

Amyloid-β (Aβ1-42) peptides aggregated into plaques deposited in the brain are the main hallmark of Alzheimer's disease (AD), a social and economic burden worldwide. In this context, insoluble Aβ1-42 fibrils are the main components of plaques. The recent trials that used approved AD drugs show that they can remove the fibrils from AD patients' brains, but they did not halt the course of the disease. Mounting evidence envisages that the soluble Aβ1-42 oligomers' interactions with the neuronal membrane trigger higher cell death than Aβ1-42 fibril interactions. Developing a compound that can alleviate the oligomer's toxicity is one of the most demanding tasks for curing the disease. We performed two molecular dynamics (MD) simulations in an explicit solvent model. In the first case, 55-μs of multiscale all-atom (AA)/coarse-grained (CG) MD simulations were carried out to decipher the impact of a previously described small anti-Aβ molecule, termed M30 (2-octahydroisoquinolin-2(1H)-ylethanamine), on an Aβ1-42 tetramer structure in close contact with a DMPC bilayer. In the second case, 15-μs AA/CG MD simulations were performed to rationalize the dynamics between Aβ1-42 and Aβ1-42-M30 tetramer complexes embedded in DMPC. On the membrane bilayer, we found that the Aβ1-42 tetramer penetrates the bilayer surface due to unrestricted conformational flexibility and many contacts with the membrane phosphate groups. In contrast, no Aβ1-42-M30 tetramer penetration was observed during the entire course of the simulation. In the case of the membrane-embedded Aβ1-42 tetramer, the integrity of the bottom bilayer leaflet was severely affected by the interactions between the negatively charged phosphate groups and the positively charged residues of the Aβ1-42 tetramer, resulting in a deep tetramer penetration into the bilayer hydrophobic region. These contacts were not observed in the case of the membrane-embedded Aβ1-42-M30 tetramer. It was noted that M30 molecules bind to Aβ1-42 tetramer through hydrogen bonds, resulting in a conformational stable Aβ1-42-M30 complex. The associated complex has reduced conformational changes and an enhanced rigidity that prevents the tetramer dissociation by interfering with the tetramer-membrane contacts. Our findings suggest that the M30 molecules could bind to Aβ1-42 tetramer resulting in a rigid structure, and that such complexes do not significantly perturb the membrane bilayer organization. These observations support the in vitro and in vivo experimental evidence that the M30 molecules prevent synaptotocity, improving AD-affected mice memory.

A Novel spice-antioxidant-based nano-vehicle as a putative green alternative of synthetic AChE inhibitor drugs

The present treatment for Alzheimer's disease (AD) involves well known synthetic acetylcholine esterase (AChE) inhibitor drugs which besides having short duration of action also have deleterious impact on human health. Therefore, there is a need for natural plant-based biomolecule(s) with potential AChE inhibition activity (ies). The aim of the work is to design a spice-based nano-vehicle as a novel green alternative of synthetic AD drugs by nanoencapsulating a solvent-less supercritical CO2 extract of small cardamom seeds (SCE) having a synergistic consortium of five antioxidant molecules, using polyethylene glycol and emulsifiers, selected based on Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) analyses. Ellman's assay and enzyme inhibition kinetics of the antioxidant molecules as well as the extract and its nanoliposomal formulation (SCE-NL) were performed, followed by rigorous molecular docking and dynamics studies using MM-PBSA and umbrella sampling. The antioxidants exhibited significant AChE inhibition in vitro, individually with 1, 8-cineole having the least IC50 value of 65.53 ± 0.05 µg/mL. . Although SCE-NL had higher IC50 value (575.67 ± 0.5 µg/mL) vis-à-vis that of rivastigmine (67.52 ± 0.02 µg/mL), it is safer for usage being 'green'.The Lineweaver-Burk plots (Vmax ∼1.04 mM/min) revealed competitive mode(s) of inhibition of AChE with each of these antioxidants. Binding energy analyses suggested very good binding free energies and stable docking/binding complexes (between the antioxidants and AChE). This study has delivered a nanoliposomal vehicle of food antioxidants as a putative 'green' alternative of synthetic AChE inhibitor drugs.Communicated by Ramaswamy H. Sarma.

Computational exploration of natural compounds targeting Staphylococcus aureus: inhibiting AgrA promoter binding for antimicrobial intervention

Staphylococcus aureus is a highly virulent nosocomial pathogen that poses a significant threat to individuals exposed to healthcare settings. Due to its sophisticated machinery for producing virulence factors, S. aureus can cause severe and potentially fatal infections in humans. This study focuses on the response regulator AgrA, which plays a crucial role in regulating the production of virulence factors in S. aureus. The objective is to identify natural compounds that can inhibit the binding of AgrA to its promoter site, thus inhibiting the expression of virulence genes. To achieve this, a pharmacophore model was generated using known drugs and applied to screen the ZINC natural product database. The resulting compounds were subjected to molecular docking-based virtual screening against the C-terminal DNA binding domain of AgrA. Three compounds, namely ZINC000077269178, ZINC000051012304, and ZINC000004266026, were shortlisted based on their strong affinity for key residues involved in DNA binding and transcription initiation. Subsequently, the unbound and ligand-bound complexes were subjected to a 200 ns molecular dynamics simulation to assess their conformational stability. Various analyses, including RMSD, RMSF, Rg, SASA, Principal Component Analysis, and Gibbs free energy landscape, were conducted on the simulation trajectory. The RMSD profile indicated similar fluctuations in both bound and unbound structures, while the Rg profile demonstrated the compactness of the protein without any unfolding during the simulation. Furthermore, Principal component analysis revealed that ligand binding reduced the overall atomic motion of the protein whereas free energy landscape suggested the energy variations obtained in complexes.Communicated by Ramaswamy H. Sarma.

Lipogenic stearoyl-CoA desaturase-1 (SCD1) targeted virtual screening for chemical inhibitors: molecular docking / dynamics simulation and in vitro assessment of anti-NAFLD efficacy

Amidst rising global prevalence of metabolic syndrome, the associated risk of non-alcoholic fatty liver disease (NAFLD) is also rapidly increasing. The pathogenesis of NAFLD starts with fat accumulation and progresses through inflammation and fibrotic sequel, often involving complex molecular mechanisms involving de novo lipogenesis. Stearoyl-CoA desaturase 1 (SCD1) enzyme, expressed in liver and adipose tissue, converts saturated fatty acids to monounsaturated fatty acids (MUFAs), contributing to triglyceride and cholesterol ester formation. In this study, potential SCD1 inhibitors were screened using the ZINC database of curated medically-approved drugs by virtual screening, molecular docking, and molecular dynamics simulations. The top-scoring five ligands with strong binding affinity against SCD1 were ZINC000003831151 > ZINC000001540998 > ZINC000003830713 > ZINC000000897251 > ZINC000002005305, which showed stable protein-ligand complexation and favorable pharmacokinetic attributes. The top ligand, Montelukast, was experimentally validated for its pharmacological efficacy in an in vitro cell culture model of steatosis (NAFLD). Montelukast showed a dose-dependent decrease in hepatic fat accumulation, reduced levels of free radicals, and lowered oxidative stress (P < 0.05). These outcomes suggest Montelukast to be a potential SCD1 inhibitor, with anti-NAFLD efficacy. These findings open new avenues for therapeutic development of the top 5 ligands in metabolic disorders involving SCD1.

Identification of efflux pump inhibitors for Pseudomonas aeruginosa MexAB-OprM via ligand-based pharmacophores, 2D-QSAR, molecular docking, and molecular dynamics approaches

Efflux pumps have been reported as one of the significant mechanisms by which bacteria evade the effects of multiple antibiotics. The tripartite efflux pump MexAB-OprM in Pseudomonas aeruginosa is one of the most significant multidrug efflux systems due to its broad resistance to antibiotics such as chloramphenicol, fluoroquinolones, lipophilic β-lactam antibiotics, nalidixic acid, novobiocin, rifampicin, and tetracycline. A promising strategy to overcome this resistance mechanism is to combine antibiotics with efflux pump inhibitors (EPIs), which can increase their intracellular concentration to enhance their biological activities. Based on 143 EPIs with chemically diverse skeletons, the 3D pharmacophore and 2D-QSAR modelings were developed and used for the virtual screening on 9.2 million compounds including ZINC15, DrugBank, and Traditional Chinese Medicine databases to identify new EPIs. The molecular docking was also performed to evaluate the binding affinity of potential EPIs to the distal-binding pocket of MexB and resulted in 611 potential EPIs. The structure-activity relationship analyses suggested that nitrogen heterocyclic compounds, piperazine and pyridine scaffolds, and amide derivatives are the most favorable chemically features for MexAB inhibitory activities. The results from molecular dynamics analysis in 100 ns indicated that ZINC009296881 and ZINC009200074 were the most potential MexB inhibitors with strong binding affinity to the distal pocket and MM/GBSA ∆Gbind values of - 38.97 and - 30.19 kcal mol-1, respectively. The predicted pharmacokinetic properties and toxicity of these compounds indicated their potential oral drugs.

Computational Analysis of the Accumulation of Mutations in Therapeutically Important RNA Viral Proteins During Pandemics with Special Emphasis on SARS-CoV-2

Single stranded RNA viruses are primary causative agents for pandemics, causing extensive morbidity and mortality worldwide. A pivotal question in pandemic preparedness and therapeutic intervention is what are the specific mutations which are more likely to emerge during such global health crises? This study aims to identify markers for mutations with the highest probability of emergence in these pandemics, focusing on the SARS-CoV-2 spike protein, an essential and therapeutically significant viral protein, starting from sequence information from the onset of the pandemic until July 2022. Quite consistently, we observed that emerged mutations tended to demonstrate a high genetic score, which reflects high similarity of the type of codon required for translation between an amino acid and to the mutated one. Further, this pattern is also observed in therapeutically significant proteins of other ssRNA pandemic viruses, including influenza (HA, NA), spike proteins of Ebola, envelope of Dengue and Chikungunya. We propose that the genetic score serves as an initial indicator, preceding the actual impact of the mutation on viral fitness. Finally, we developed a comprehensive computational pipeline to further explore and predict the subsequent effects of mutations on viral fitness. We believe that our pipeline can narrow down and predict future mutations in therapeutically important viral proteins during a pandemic.

Molecular docking and molecular dynamics studies of natural products unravel potential inhibitors against OmpA of Acinetobacter baumannii

Emerging antimicrobial resistance has highlighted the need to design more effective antibiotics to treat deadly bacterial infections. Acinetobacter baumannii's outer membrane protein A (OmpA) is a critical virulence component involved in biofilm formation, immunomodulation, and antibiotic resistance, which characterizes it as a potential therapeutic target. The present study aimed to screen the natural product database (>1,00,000) to identify the potential inhibitor against OmpA. Molecular docking studies revealed that 10 compounds had good docking scores (≤ -7 kcal/mol) compared to the reported inhibitor epiestriol (-3.079). Further, these 10 compounds were subjected to ADME analysis and MMGBSA analysis. Based on MMGBSA results, we selected 5 compounds [NP-1 (MolPort-039-337-117), NP-5(MolPort-019-932-973), NP-6 (MolPort-005-948-336), NP-8(MolPort-042-673-978) and NP-9(MolPort-042-673-766)] with high binding affinity. Molecular dynamics simulation found that NP-5, NP-8, and NP-9 were stable after analysing their RMSD, RMSF, the radius of gyration, and hydrogen interactions of complexes. Our study revealed that NP-5, NP-8, and NP-9 bind perfectly with OmpA and can act as its potential inhibitors. The results of this study imply that the identified inhibitors have the potential for further investigation.Communicated by Ramaswamy H. Sarma.

Potential suppression of multidrug-resistance-associated protein 1 by coumarin derivatives: an insight from molecular docking and MD simulation studies

Human MRP1 protein plays a vital role in cancer multidrug resistance. Coumarins show promising pharmacological properties. Virtual screening, ADMET, molecular docking and molecular dynamics (MD) simulations were utilized as pharmacoinformatic tools to identify potential MRP1 inhibitors among coumarin derivatives. Using in silico ADMET, 50 hits were further investigated for their selectivity toward the nucleotide-binding domains (NBDs) of MRP1 using molecular docking. Accordingly, coumarin, its symmetrical ketone derivative Lig. No. 4, and Reversan were candidates for focused docking study with the NBDs domains compared with ATP. The result indicates that Lig. No. 4, with the best binding score, interacts with NBDs via hydrogen bonds with residues: GLN713, LYS684, GLY683, CYS682 in NBD1, and GLY1432, GLY771, SER769 and GLN1374 in NBD2, which mostly overlap with ATP binding residues. Moreover, doxorubicin (Doxo) was docked to the transmembrane domains (TMDs) active site of MRP1. Doxo interaction with TMDs was subjected to MD simulation in the NBDs free and occupied with Lig. No. 4 states. The results showed that Doxo interacts more strongly with TMD residues in inward facing feature of TMDs helices. However, when Lig. No. 4 exists in NBDs, Doxo interactions are different, and TMD helices show more outward-facing conformation. This result may suggest a partial competitive inhibition mechanism for the Lig. No. 4 on MRP1 compared with ATP. So, it may inhibit active complex formation by interfering with ATP entrance to NBDs and locking MRP1 conformation in outward-facing mode. This study suggests a valuable coumarin derivative that can be further investigated for potent MRP1 inhibitors.Communicated by Ramaswamy H. Sarma.

E-pharmacophore and deep learning based high throughput virtual screening for identification of CDPK1 inhibitors of Cryptosporidium parvum

Cryptosporidiosis, a prevalent gastrointestinal illness worldwide, is caused by the protozoan parasite Cryptosporidium parvum. Calcium-dependent protein kinase 1 (CpCDPK1), crucial for the parasite's life cycle, serves as a promising drug target due to its role in regulating invasion and egress from host cells. While potent Pyrazolopyrimidine analogs have been identified as candidate hit molecules, they exhibit limitations in inhibiting Cryptosporidium growth in cell culture, prompting exploration of alternative scaffolds. Leveraging the most potent compound, RM-1-95, co-crystallized with CpCDPK1, an E-pharmacophore model was generated and validated alongside a deep learning model trained on known CpCDPK1 compounds. These models facilitated screening Enamine's 2 million HTS compound library for novel CpCDPK1 inhibitors. Subsequent hierarchical docking prioritized hits, with final selections subjected to Quantum polarized docking for accurate ranking. Results from docking studies and MD simulations highlighted similarities in interactions between the cocrystallized ligand RM-1-95 and identified hit molecules, indicating comparable inhibitory potential against CpCDPK1. Furthermore, assessing metabolic stability through Cytochrome 450 site of metabolism prediction offered crucial insights for drug design, optimization, and regulatory approval processes.

Computational simulations uncover enantioselective metabolism of chiral triazole fungicides by human CYP450 enzymes: A case study of tebuconazole

Tebuconazole (TEB), a prominent chiral triazole fungicide, has been extensively utilized for plant pathogen control globally. Despite experimental evidence of TEB metabolism in mammals, the enantioselectivity in the biotransformation of R- and S-TEB enantiomers by specific CYP450s remains elusive. In this work, integrated in silico simulations were employed to unveil the binding interactions and enantioselective metabolic fate of TEB enantiomers within human CYP1A2, 2B6, 2E1, and 3A4. Molecular dynamics (MD) simulations clearly delineated the binding specificity of R- and S-TEB to the four CYP450s, crucially determining their differences in metabolic activity and enantioselectivity. The primary driving force for robust ligand binding was identified as van der Waals interactions with CYP450s, particularly involving the hydrophobic residues. Mechanistic insights derived from quantum mechanics/molecular mechanics (QM/MM) calculations established C2-methyl hydroxylation as the predominant route of R-/S-TEB metabolism, while C6-hydroxylation and triazol epoxidation were deemed kinetically infeasible pathways. Specifically, the resulting hydroxy-R-TEB metabolite primarily originates from R-TEB biotransformation by 1A2, 2E1 and 3A4, whereas hydroxy-S-TEB is preferentially produced by 2B6. These findings significantly contribute to our comprehension of the binding specificity and enantioselective metabolic fate of chiral TEB by CYP450s, potentially informing further research on human health risk assessment associated with TEB exposure.

Efficiency enhancement in Aspergillus niger α-L-rhamnosidase reverse hydrolysis by using a tunnel site rational design strategy

There has been ongoing interest in improving the efficiency of glycoside hydrolase for synthesizing glycoside compounds through protein engineering, given the potential applications of glycoside compounds. In this study, a strategy of modifying the substrate access tunnel was proposed to enhance the efficiency of reverse hydrolysis catalyzed by Aspergillus niger α-L-rhamnosidase. Analysis of the tunnel dynamics identified Tyr299 as a key modifiable residue in the substrate access tunnel. The location of Tyr299 was near the enzyme surface and at the outermost end of the substrate access tunnel, suggested its role in substrate recognition and throughput. Based on the properties of side chains, six mutants were designed and expressed by Pichia pastoris. Compared to WT, the reverse hydrolysis efficiencies of mutants Y299P and Y299W were increased by 21.3 % and 11.1 %, respectively. The calculation results of binding free energy showed that the binding free energy was inversely proportional to the reverse hydrolysis efficiency. Further, when binding free energy levels were comparable, the mutants with shorter side chains displayed a higher reverse hydrolysis efficiency. These results proved that substrate access tunnel modification was an effective method to improve the reverse hydrolysis efficacy of α-L-rhamnosidase and also provided new insights for modifying other glycoside hydrolases.

Exploiting the bile acid binding protein as transporter of a Cholic Acid/Mirin bioconjugate for potential applications in liver cancer therapy

Bioconjugation is one of the most promising strategies to improve drug delivery, especially in cancer therapy. Biomolecules such as bile acids (BAs) have been intensively explored as carriers, due to their peculiar physicochemical properties and biocompatibility. BAs trafficking is regulated by intracellular lipid-binding proteins and their transport in the liver can be studied using chicken liver Bile Acid-Binding Proteins (cL-BABPs) as a reference model. Therefore, we conceived the idea of developing a BA-conjugate with Mirin, an exonuclease inhibitor of Mre11 endowed with different anticancer activities, to direct its transport to the liver. Following computational analysis of various BAs in complex with cL-BABP, we identified cholic acid (CA) as the most promising candidate as carrier, leading to the synthesis of a novel bioconjugate named CA-M11. As predicted by computational data and confirmed by X-ray crystallographic studies, CA-M11 was able to accommodate into the binding pocket of BABP. Hence, it can enter BAs trafficking in the hepatic compartment and here release Mirin. The effect of CA-M11, evaluated in combination with varying concentrations of Doxorubicin on HepG2 cell line, demonstrated a significant increase in cell mortality compared to the use of the cytotoxic drug or Mirin alone, thus highlighting chemo-sensitizing properties. The promising results regarding plasma stability for CA-M11 validate its potential as a valuable agent or adjuvant for hepatic cancer therapy.

Targeting the phosphatidylglycerol lipid: An amphiphilic dendrimer as a promising antibacterial candidate

The rapid emergence and spread of multidrug-resistant bacterial pathogens require the development of antibacterial agents that are robustly effective while inducing no toxicity or resistance development. In this context, we designed and synthesized amphiphilic dendrimers as antibacterial candidates. We report the promising potent antibacterial activity shown by the amphiphilic dendrimer AD1b, composed of a long hydrophobic alkyl chain and a tertiary amine-terminated poly(amidoamine) dendron, against a panel of Gram-negative bacteria, including multidrug-resistant Escherichia coli and Acinetobacter baumannii. AD1b exhibited effective activity against drug-resistant bacterial infections in vivo. Mechanistic studies revealed that AD1b targeted the membrane phospholipids phosphatidylglycerol (PG) and cardiolipin (CL), leading to the disruption of the bacterial membrane and proton motive force, metabolic disturbance, leakage of cellular components, and, ultimately, cell death. Together, AD1b that specifically interacts with PG/CL in bacterial membranes supports the use of small amphiphilic dendrimers as a promising strategy to target drug-resistant bacterial pathogens and addresses the global antibiotic crisis.

A Study on the Binding Mechanism and the Impact of Key Residue Mutations between SND1 and MTDH Peptide through Molecular Dynamics Simulations

Metastasis of breast cancer is the main cause of death for patients with breast cancer. The interaction between metadherin (MTDH) and staphylococcal nuclease domain 1 (SND1) plays a pivotal role in promoting breast cancer development. However, the binding details between MTDH and SND1 remain unclear. In this study, we employed all-atom molecular dynamics simulations (MDs) and conducted binding energy calculations to investigate the binding details and the impact of key residue mutations on binding. The mutations in key residues have not significantly affected the overall stability of the structure and the fluctuation of residues near the binding site; they have exerted a substantial impact on the binding of SND1 and MTDH peptide. The electrostatic interactions and van der Waals interactions play an important role in the binding of SND1 and the MTDH peptide. The mutations in the key residues have a significant impact on electrostatic and van der Waals interactions, resulting in weakened binding. The energy contributions of key residues mainly come from the electrostatic energy and van der Waals interactions of the side chain. In addition, the key residues form an intricate and stable network of hydrogen bonds and salt-bridge interactions with the MTDH peptide. The mutations in key residues have directly disrupt the interactions formed between SND1 and MTDH peptide, consequently leading to changes in the binding mode of the MTDH peptide. These analyses unveil the detailed atomic-level interaction mechanism between SND1 and the MTDH peptide, providing a molecular foundation for the development of antibreast cancer drugs.

Insight from atomistic molecular dynamics simulations into the supramolecular assembly of the aldo-keto reductase from Trypanosoma cruzi

CONTEXT

Currently, Chagas disease represents an important public health problem affecting more than 8 million people worldwide. The vector of this disease is the Trypanosoma cruzi (Tc) parasite. Our research specifically focuses on the structure and aggregation states of the enzyme aldo-keto reductase of Tc (TcAKR) reported in this parasite. TcAKR belongs to the aldo-keto reductase (AKR) superfamily, enzymes that catalyze redox reactions involved in crucial biological processes. While most AKRs are found in monomeric forms, some have been reported to form dimeric and tetrameric structures. This is the case for some TcAKR. To better understand how TcAKR multimers form and remain stable, we conducted a comprehensive computational analysis using molecular dynamics (MD) simulations. Our approach to elucidating the aggregation states of TcAKR involved two strategies. Initially, we explored the dynamic behaviour of pre-assembled TcAKR dimers. Subsequently, we examined the self-aggregation of eight monomers. This investigation led to the identification of crucial residues that contribute to the stabilization of protein-protein interactions. It was also found that TcAKRs can form stable supramolecular assemblies, with each monomer typically surrounded by three first neighbours. These findings align with experimental reports of tetrameric or more complex supramolecular structures. Our computational studies could guide further experimental investigations aiming at drug development and assist in designing strategies to modulate aggregation.

METHOD

Atomistic molecular dynamics simulations were carried out. The TcAKR 3D model structure was obtained by homology modelling using the Swiss Model for the TcAKR sequence (GenBank accession no. EU558869). Further, we checked the model with Alphafold2 and found a high degree of similarity between models. Several tools were used to build the dimers including CLUSPRO, GRAMM-Docking, Hdock, and Py-dock. Protein superstructures were built using the PACKMOL package. CHARMM-GUI was used to set up the simulation systems. GROMACS version 2020.5 was used to perform the simulations with the CHARMM36 force field for the protein and ions and the TIP3P model for water. Further analyses were performed using VMD, GROMACS, AMBER tools, MDLovoFit, bio3d, and in-house programs.

Analysis of factors affecting microbial degradation of cyanobacterial toxins based on theoretical calculations

Cyanobacterial toxins are the most common algal toxins, which are highly toxic and can persist in the aquatic environment without easy degradation, posing risks to the ecosystem and human health that cannot be ignored. Although microbiological methods for the removal of cyanobacterial toxins from aqueous environments are highly efficient, their degradation efficiency is susceptible to many abiotic environmental factors. In this paper, Microcystin-LR (MC-LR) and its microbial degrading enzymes were selected to study the effects of common environmental factors (temperature (T), NO3-, NH4+, Cu2+, Zn2+) and their levels during microbial degradation of cyanobacterial toxins in aqueous environments by using molecular docking, molecular dynamics simulation, analytical factor design, and the combined toxicokinetics of TOPKAT (toxicity prediction). It was found that the addition of T, NO3- and Cu2+ to the aqueous environment promoted the microbial degradation of MC-LR, while the addition of NH4+ and Zn2+ inhibited the degradation; The level effect study showed that the microbial degradation of MC-LR was promoted by increasing levels of added T and NO3- in the aqueous environment, whereas it was inhibited and then promoted by increasing levels of NH4+, Cu2+ and Zn2+. In addition, the predicted toxicity of common Microcystins (MCs) showed that MC-LR, Microcystin-RR (MC-RR) and Microcystin-YR (MC-YR) were not carcinogenic, developmentally toxic, mutagenic or ocular irritants in humans. MC-LR and MC-RR are mild skin irritants and MC-YR is not a skin irritant. MC-YR has a higher chronic and acute toxicity in humans than MC-LR and MC-RR. Acute/chronic toxicity intensity for aquatic animals: MC-YR > MC-LR > MC-RR and for aquatic plants: MC-LR > MC-YR > MC-RR. This suggests that MC-YR also has a high environmental health risk. This paper provides theoretical support for optimizing the environmental conditions for microbial degradation of cyanobacterial toxins by studying the effects of common environmental factors and their level effects in the aquatic environment.

Molecular Photoswitching Unlocks Glucose Oxidase for Synergistically Reinforcing Fenton Reactions for Antitumor Chemodynamic Therapy

We have developed a new type of nanoparticles with potent antitumor activity photoactivatable via the combination of molecular photoswitching of spiropyran (SP) and enzymatic reaction of glucose oxidase (GOx). As two key processes involved therein, Fe(III)-to-Fe(II) photoreduction in Fe(III) metal-organic frameworks (MOFs) brings about the release of free Fe2+/Fe3+ while the photoswitching of SP to merocyanine (MC) unlocks the enzymatic activity of GOx that was pre-passivated by SP. The release of free Fe3+ boosts its hydrolysis and therefore enables the acidification of microenvironment, which is further reinforced by one of the products of the GOx-mediated glucose oxidation reaction, gluconic acid (GlcA). Based on the generation of Fe2+ and acidic milieu together with another product of the oxidation reaction, hydrogen peroxide (H2O2), these two processes jointly present triple enabling factors for generating lethal hydroxyl radicals (⋅OH) species via Fenton reactions and therefore oxidative stress capable of inhibiting tumor. The antitumor potency of such nanoparticle is verified in tumor-bearing model mice in vivo, proclaiming its potential as a potent and safe agent based on the unique mechanism of optically manipulating enzyme activity for synergistic antitumor therapeutics with high spatial precision, enhanced efficacy and minimized side effects.

Stereo-selectivity of enantiomeric inhibitors to ubiquitin-specific protease 7 (USP7) dissected by molecular docking, molecular dynamics simulations, and binding free energy calculations

The ubiquitin-specific protease 7 (USP7), as a member of deubiquitination enzymes, represents an attractive therapeutic target for various cancers, including prostate cancer and liver cancer. The change of the inhibitor stereocenter from the S to R stereochemistry (S-ALM → R-ALM34) markedly improved USP7 inhibitory activity. However, the molecular mechanism for the stereo-selectivity of enantiomeric inhibitors to USP7 is still unclear. In this work, molecular docking, molecular dynamics (MD) simulations, molecular mechanics/Generalized-Born surface area (MM/GBSA) calculations, and free energy landscapes were performed to address this mystery. MD simulations revealed that S-ALM34 showed a high degree of conformational flexibility compared to the R-ALM34 counterpart, and S-ALM34 binding led to the enhanced intradomain motions of USP7, especially the BL1 and BL2 loops and the two helices α4 and α5. MM/GBSA calculations showed that the binding strength of R-ALM34 to USP7 was stronger than that of S-ALM34 by - 4.99 kcal/mol, a similar trend observed by experimental data. MM/GBSA free energy decomposition was further performed to differentiate the ligand-residue spectrum. These analyses not only identified the hotspot residues interacting with R-ALM34, but also revealed that the hydrophobic interactions from F409, K420, H456, and Y514 play the major determinants in the binding of R-ALM34 to USP7. This result is anticipated to shed light on energetic basis and conformational dynamics information to aid in the design of more potent and selective inhibitors targeting USP7.

The potential of Mitragyna speciosa leaves as a natural source of antioxidants for disease prevention

Mitragyna speciosa is famous for its addictive effect. On the other hand, this plant has good potential as an antioxidant agent, and so far, it was not explicitly explained what the most contributing compound in the leaves to that activity is. This study has been conducted using several computational methods to determine which compounds are the most active in interacting with cytochrome P450, myeloperoxidase, and NADPH oxidase proteins. First, virtual screening was carried out based on molecular docking, followed by profiling the properties of adsorption, distribution, metabolism, excretion, and toxicity (ADMET); the second one is the molecular dynamics (MD) simulations for 100 ns. The virtual screening results showed that three compounds acted as inhibitors for each protein: (-)-epicatechin, sitogluside, and corynoxeine. The ADMET profiles of the three compounds exhibit good drug ability and toxicity. The trajectories study from MD simulations predicts that the complexes of these three compounds with their respective target proteins are stable. Furthermore, these compounds identified in this computational study can be a potential guide for future experiments aimed at assessing the antioxidant properties through in vitro testing.

Repurposing of USFDA-approved drugs to identify leads for inhibition of acetylcholinesterase enzyme: a plausible utility as an anti-Alzheimer agent

In the quest to identify new anti-Alzheimer agents, we employed drug repositioning or drug repositioning techniques on approved USFDA small molecules. Herein, we report the structure-based virtual screening (SBVS) of 1880 USFDA-approved drugs. The in silico-based identification was followed by calculating Prime MMGB-SA binding energy and molecular dynamics simulation studies. The cumulative analysis led to identifying domperidone as an identified hit. Domperidone was further corroborated in vitro using anticholinesterase-based assessment, keeping donepezil as a positive control. The analysis revealed that the identified lead (domperidone) could induce an inhibitory effect on AChE in a dose-dependent manner with an IC50 of 3.67 μM as compared to donepezil, which exhibited an IC50 of 1.37 μM. However, as domperidone is known to have poor BBB permeability, we rationally proposed new analogues utilizing the principles of bioisosterism. The bioisostere-clubbed analogues were found to have better BBB permeability, affinity, and stability within the catalytic domain of AChE via molecular docking and dynamics studies. The proposed bioisosteres may be synthesized in the future. They may plausibly be explored for their implication in the developmental progress of new anti-Alzheimer agent achieved via repurposing techniques in future.