A method and program for mass quantum chemical calculations of protein—ligand docking complexes

A new fast computational method for mass calculations of docking complexes by the AM1/PM3 semiempirical methods is proposed. The computation time is shortened by at least an order of magnitude compared to alternative schemes of quantum chemical calculations. The root-mean-square deviation of the AM1 calculated energies of formation of complexes from the results obtained by conventional diagonalization procedure is at most 0.4 kcal mol⁻¹. Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 418–420, February, 2008.     

Molecular modeling of cytochrome P450 3A4

The three-dimensional structure of human cytochrome P450 3A4 was modeled based on crystallographic coordinates of four bacterial P450s: P450 BM-3, P450cam, P450terp, and P450eryF. The P450 3A4 sequence was aligned to those of the known proteins using a structure-based alignment of P450 BM-3, P450cam, P450terp, and P450eryF. The coordinates of the model were then calculated using a consensus strategy, and the final structure was optimized in the presence of water. The P450 3A4 model resembles P450 BM-3 the most, but the B helix is similar to that of P450eryF, which leads to an enlarged active site when compared with P450 BM-3, P450cam, and P450terp. The 3A4 residues equivalent to known substrate contact residues of the bacterial proteins and key residues of rat P450 2B1 are located in the active site or the substrate access channel. Docking of progesterone into the P450 3A4 model demonstrated that the substrate bound in a 6-orientation can interact with a number of active site residues, such as 114, 119, 301, 304, 305, 309, 370, 373, and 479, through hydrophobic interactions. The active site of the enzyme can also accommodate erythromycin, which, in addition to the residues listed for progesterone, also contacts residues 101, 104, 105, 214, 215, 217, 218, 374, and 478. The majority of 3A4 residues which interact with progesterone and/or erythromycin possess their equivalents in key residues of P450 2B enzymes, except for residues 297, 480 and 482, which do not contact either substrate in P450 3A4. The results from docking of progesterone and erythromycin into the enzyme model make it possible to pinpoint residues which may be important for 3A4 function and to target them for site-directed mutagenesis.     

Virtual Screening and Hit Selection of Natural Compounds as Acetylcholinesterase Inhibitors

Acetylcholinesterase (AChE) is one of the classical targets in the treatment of Alzheimer’s disease (AD). Inhibition of AChE slows down the hydrolysis of acetycholine and increases choline levels, improving the cognitive function. The achieved success of plant-based natural drugs acting as AChE inhibitors, such as galantamine (GAL) from Galanthus genus and huperzine A from Huperzia serrate (approved drug in China), in the treatment of AD, and the fact that natural compounds (NCs) are considered as safer and less toxic compared to synthetic drugs, led us to screen the available NCs (almost 150,000) in the ZINC12 database for AChE inhibitory activity. The compounds were screened virtually by molecular docking, filtered for suitable ADME properties, and 32 ligands from 23 structural groups were selected. The stability of the complexes was estimated via 1 μs molecular dynamics simulation. Ten compounds formed stable complexes with the enzyme and had a vendor and a reasonable price per mg. They were tested for AChE inhibitory and antioxidant activity. Five compounds showed weak AChE inhibition and three of them exhibited high antioxidant activity.     

Chemical and Pharmacological Profiling of Wrightia coccinea (Roxb. Ex Hornem.) Sims Focusing Antioxidant,Cytotoxic, Antidiarrheal,Hypoglycemic, and Analgesic Properties

The aim of the study was to conduct phytochemical and pharmacological investigations of Wrightia coccinea (Roxb. ex Hornem.) Sims via several in vitro, in vivo, and in silico models. A total of four compounds were identified and isolated from the methanol extract of the bark and the methanol extract of the seed pulp of W. coccinea through successive chromatographic techniques and were characterized as 3β-acetyloxy-olean-12-en-28-ol (1), wrightiadione (2), 22β-hydroxylupeol (3), and β-sitosterol (4) by spectroscopic analysis. The aqueous fraction of the bark and chloroform fraction of the fruits provided the most potent antioxidant capacity (IC₅₀ = 7.22 and 4.5 µg/mL, respectively) in DPPH free radical scavenging assay compared with the standard ascorbic acid (IC₅₀ = 17.45 µg/mL). The methanol bark extract and the methanol fruit coat extract exerted anti-diarrheal activity by inhibiting 74.55 ± 0.67% and 77.78 ± 1.5% (mean ± SEM) of the diarrheal episode in mice, respectively, after four hours of loading the samples. In the hypoglycemic test, the methanol bark extract and the methanol fruit coat extract (400 mg/kg) produced a significant (p < 0.05) reduction in the blood glucose level in mice. Both doses of the plant extracts (200 mg/kg and 400 mg/kg) used in the study induced a significant (p < 0.05) increase in pain reaction time. The in vitro and in vivo findings were supported by the computational studies. The isolated compounds exhibited higher binding affinity compared with the standard drugs towards the active binding sites of glutathione reductase, epidermal growth factor receptor (EGFR), kappa opioid receptor, glucose transporter 3 (GLUT 3), Mu opioid receptor, and cyclooxygenase 2 (COX-2) proteins due to their potent antioxidant, cytotoxic, anti-diarrheal, hypoglycemic, and central and peripheral analgesic properties, respectively. The current findings concluded that W. coccinea might be a potential natural source for managing oxidative stress, diarrhea, hyperglycemia, and pain. Further studies are warranted for extensively phytochemical screening and establishing exact mechanisms of action.     

Investigation Driven by Network Pharmacology on Potential Components and Mechanism of DGS,a Natural Vasoprotective Combination,for the Phytotherapy of Coronary Artery Disease

Phytotherapy offers obvious advantages in the intervention of Coronary Artery Disease (CAD), but it is difficult to clarify the working mechanisms of the medicinal materials it uses. DGS is a natural vasoprotective combination that was screened out in our previous research, yet its potential components and mechanisms are unknown. Therefore, in this study, HPLC-MS and network pharmacology were employed to identify the active components and key signaling pathways of DGS. Transgenic zebrafish and HUVECs cell assays were used to evaluate the effectiveness of DGS. A total of 37 potentially active compounds were identified that interacted with 112 potential targets of CAD. Furthermore, PI3K-Akt, MAPK, relaxin, VEGF, and other signal pathways were determined to be the most promising DGS-mediated pathways. NO kit, ELISA, and Western blot results showed that DGS significantly promoted NO and VEGFA secretion via the upregulation of VEGFR2 expression and the phosphorylation of Akt, Erk1/2, and eNOS to cause angiogenesis and vasodilation. The result of dynamics molecular docking indicated that Salvianolic acid C may be a key active component of DGS in the treatment of CAD. In conclusion, this study has shed light on the network molecular mechanism of DGS for the intervention of CAD using a network pharmacology-driven strategy for the first time to aid in the intervention of CAD.     

Three New Acrylic Acid Derivatives from Achillea mellifolium as Potential Inhibitors of Urease from Jack Bean and α-Glucosidase from Saccharomyces cerevisiae

(1) Background: Achillea mellifolium belongs to a highly reputed family of medicinal plants, with plant extract being used as medicine in indigenous system. However, limited data is available regarding the exploitation of the medicinal potential of isolated pure compounds from this family; (2) Methods: A whole plant extract was partitioned into fractions and on the basis of biological activity, an ethyl acetate fraction was selected for isolation of pure compounds. Isolated compounds were characterized using different spectroscopic techniques. The compounds isolated from this study were tested for their medicinal potential using in-vitro enzyme assay, coupled with in-silico studies; (3) Results: Three new acrylic acid derivatives (1–3) have been isolated from the ethyl acetate fraction of Achillea mellifolium. The characterization of these compounds (1–3) was carried out using UV/Vis, FT-IR, 1D and 2D-NMR spectroscopy (¹H-NMR, ¹³C-NMR, HMBC, NOESY) and mass spectrometry. These acrylic acid derivatives were further evaluated for their enzyme inhibition potential against urease from jack bean and α glucosidase from Saccharomyces cerevisiae, using both in-silico and in-vitro approaches. In-vitro studies showed that compound 3 has the highest inhibition against urease enzyme (IC₅₀ =10.46 ± 0.03 μΜ), followed by compound 1 and compound 2 with percent inhibition and IC₅₀ value of 16.87 ± 0.02 c and 13.71 ± 0.07 μΜ, respectively, compared to the standard (thiourea-IC₅₀ = 21.5 ± 0.01 μΜ). The investigated IC₅₀ value of compound 3 against the urease enzyme is two times lower compared to thiourea, suggesting that this compound is twice as active compared to the standard drug. On the other hand, all three compounds (1–3) revealed mild inhibition potential against α-glucosidase. In-silico molecular docking studies, in combination with MD simulations and free energy, calculations were also performed to rationalize their time evolved mode of interaction inside the active pocket. Binding energies were computed using a MMPBSA approach, and the role of individual residues to overall binding of the inhibitors inside the active pockets were also computed; (4) Conclusions: Together, these studies confirm the inhibitory potential of isolated acrylic acid derivatives against both urease and α-glucosidase enzymes; however, their inhibition potential is better for urease enzyme even when compared to the standard.     

Sarcorucinine-D Inhibits Cholinesterases and Calcium Channels: Molecular Dynamics Simulation and In Vitro Mechanistic Investigations

Acetylcholinesterase (AChE) inhibitors and calcium channel blockers are considered effective therapies for Alzheimer’s disease. AChE plays an essential role in the nervous system by catalyzing the hydrolysis of the neurotransmitter acetylcholine. In this study, the inhibition of the enzyme AChE by Sarcorucinine-D, a pregnane type steroidal alkaloid, was investigated with experimental enzyme kinetics and molecular dynamics (MD) simulation techniques. Kinetics studies showed that Sarcorucinine-D inhibits two cholinesterases—AChE and butyrylcholinesterase (BChE)—noncompetitively, with K_i values of 103.3 and 4.66 µM, respectively. In silico ligand-protein docking and MD simulation studies conducted on AChE predicted that Sarcorucinine-D interacted via hydrophobic interactions and hydrogen bonds with the residues of the active-site gorge of AChE. Sarcorucinine-D was able to relax contractility concentration-dependently in the intestinal smooth muscles of jejunum obtained from rabbits. Not only was the spontaneous spasmogenicity inhibited, but it also suppressed K⁺-mediated spasmogenicity, indicating an effect via the inhibition of voltage-dependent Ca²⁺ channels. Sarcorucinine-D could be considered a potential lead molecule based on its properties as a noncompetitive AChE inhibitor and a Ca²⁺ channel blocker.     

Structural Investigation of the Interaction Mechanism between Chlorogenic Acid and AMPA Receptor via In Silico Approaches

Chlorogenic acid (CGA), an important metabolite in natural plant medicines such as honeysuckle and eucommia, has been shown to have potent antinociceptive effects. Nevertheless, the mechanism by which CGA relieves chronic pain remains unclear. α-amino-3-hydroxy-5-methyl-4-isooxazolpropionic acid receptor (AMPAR) is a major ionotropic glutamate receptor that mediates rapid excitatory synaptic transmission and its glutamate ionotropic receptor AMPA type subunit 1 (GluA1) plays a key role in nociceptive transmission. In this study, we used Western blot, surface plasmon resonance (SPR) assay, and the molecular simulation technologies to investigate the mechanism of interaction between CGA and AMPAR to relieve chronic pain. Our results indicate that the protein expression level of GluA1 showed a dependent decrease as the concentration of CGA increased (0, 50, 100, and 200 μM). The SPR assay demonstrates that CGA can directly bind to GluA1 (K_D = 496 μM). Furthermore, CGA forms a stable binding interaction with GluA1, which is validated by molecular dynamics (MD) simulation. The binding free energy between CGA and GluA1 is −39.803 ± 14.772 kJ/mol, where van der Waals interaction and electrostatic interaction are the major contributors to the GluA1–CGA binding, and the key residues are identified (Val-32, Glu-33, Ala-36, Glu-37, Leu-48), which play a crucial role in the binding interaction. This study first reveals the structural basis of the stable interaction between CGA and GluA1 to form a binding complex for the relief of chronic pain. The research provides the structural basis to understand the treatment of chronic pain and is valuable to the design of novel drug molecules in the future.     

DNA Interaction,DNA Photocleavage,Photocytotoxicity In Vitro,and Molecular Docking of Naphthyl-Appended Ruthenium Complexes

With the development of metal-based drugs, Ru(II) compounds present potential applications of PDT (photodynamic therapy) and anticancer reagents. We herein synthesized two naphthyl-appended ruthenium complexes by the combination of the ligand with naphthyl and bipyridyl. The DNA affinities, photocleavage abilities, and photocytotoxicity were studied by various spectral methods, viscosity measurement, theoretical computation method, gel electrophoresis, and MTT method. Two complexes exhibited strong interaction with calf thymus DNA by intercalation. Production of singlet oxygen (¹O₂) led to obvious DNA photocleavage activities of two complexes under 365 nm light. Furthermore, two complexes displayed obvious photocytotoxicity and low dark cytotoxicity towards Hela, A549, and A375 cells.     

An Integrated Mass Spectrometry-Based Glycomics-Driven Glycoproteomics Analytical Platform to Functionally Characterize Glycosylation Inhibitors

Cancer progression is linked to aberrant protein glycosylation due to the overexpression of several glycosylation enzymes. These enzymes are underexploited as potential anticancer drug targets and the development of rapid-screening methods and identification of glycosylation inhibitors are highly sought. An integrated bioinformatics and mass spectrometry-based glycomics-driven glycoproteomics analysis pipeline was performed to identify an N-glycan inhibitor against lung cancer cells. Combined network pharmacology and in silico screening approaches were used to identify a potential inhibitor, pictilisib, against several glycosylation-related proteins, such as Alpha1-6FucT, GlcNAcT-V, and Alpha2,6-ST-I. A glycomics assay of lung cancer cells treated with pictilisib showed a significant reduction in the fucosylation and sialylation of N-glycans, with an increase in high mannose-type glycans. Proteomics analysis and in vitro assays also showed significant upregulation of the proteins involved in apoptosis and cell adhesion, and the downregulation of proteins involved in cell cycle regulation, mRNA processing, and protein translation. Site-specific glycoproteomics analysis further showed that glycoproteins with reduced fucosylation and sialylation were involved in apoptosis, cell adhesion, DNA damage repair, and chemical response processes. To determine how the alterations in N-glycosylation impact glycoprotein dynamics, modeling of changes in glycan interactions of the ITGA5–ITGB1 (Integrin alpha 5-Integrin beta-1) complex revealed specific glycosites at the interface of these proteins that, when highly fucosylated and sialylated, such as in untreated A549 cells, form greater hydrogen bonding interactions compared to the high mannose-types in pictilisib-treated A549 cells. This study highlights the use of mass spectrometry to identify a potential glycosylation inhibitor and assessment of its impact on cell surface glycoprotein abundance and protein–protein interaction.