1
|
Yıldırım Akdeniz G, Timuçin AC. Structure based computational RNA design towards MafA transcriptional repressor implicated in multiple myeloma. J Mol Graph Model 2024; 132:108839. [PMID: 39096645 DOI: 10.1016/j.jmgm.2024.108839] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2024] [Revised: 07/24/2024] [Accepted: 07/29/2024] [Indexed: 08/05/2024]
Abstract
Multiple myeloma is recognized as the second most common hematological cancer. MafA transcriptional repressor is an established mediator of myelomagenesis. While there are multitude of drugs available for targeting various effectors in multiple myeloma, current literature lacks a candidate RNA based MafA modulator. Thus, using the structure of MafA homodimer-consensus target DNA, a computational effort was implemented to design a novel RNA based chemical modulator against MafA. First, available MafA-consensus DNA structure was employed to generate an RNA library. This library was further subjected to global docking to select the most plausible RNA candidates, preferring to bind DNA binding region of MafA. Following global docking, MD-ready complexes that were prepared via local docking program, were subjected to 500 ns of MD simulations. First, each of these MD simulations were analyzed for relative binding free energy through MM-PBSA method, which pointed towards a strong RNA based MafA binder, RNA1. Second, through a detailed MD analysis, RNA1 was shown to prefer binding to a single monomer of the dimeric DNA binding domain of MafA using higher number of hydrophobic interactions compared with positive control MafA-DNA complex. At the final phase, a principal component analyses was conducted, which led us to identify the actual interaction region of RNA1 and MafA monomer. Overall, to our knowledge, this is the first computational study that presents an RNA molecule capable of potentially targeting MafA protein. Furthermore, limitations of our study together with possible future implications of RNA1 in multiple myeloma were also discussed.
Collapse
Affiliation(s)
- Güneş Yıldırım Akdeniz
- Department of Molecular Biology, Genetics and Bioengineering, Faculty of Engineering and Natural Sciences, Sabancı University, 34956, Tuzla, İstanbul, Turkey.
| | - Ahmet Can Timuçin
- Department of Molecular Biology and Genetics, Faculty of Engineering and Natural Sciences, Acıbadem Mehmet Ali Aydınlar University, 34752, Ataşehir, İstanbul, Turkey.
| |
Collapse
|
2
|
Ganapathy Vilasam Sreekala A, Gupta KK, Nathan VK. Identification of coastal pesticide pollutants as potent inhibitors of Bacillus pasteurii urease mediated calcium carbonate precipitation: a computational approach. J Biomol Struct Dyn 2024; 42:9628-9638. [PMID: 37691444 DOI: 10.1080/07391102.2023.2252089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2023] [Accepted: 08/21/2023] [Indexed: 09/12/2023]
Abstract
Microbially induced calcite precipitation (MICP) through urease enzyme has attained a lot of recognition in various fields of civil engineering and geotechnology for stabilizing the strength of soil and various concrete materials. The activity of urease has been found to be affected by various factors like temperature, substrate concentrations, pH of the medium, presence of inhibitors, etc. Through this study, the outcome of the interaction of pesticides (commonly found in Indian coastal regions) on Bacillus pasteurii urease, a major organism reported for MICP studies has been investigated in silico. The results from the study revealed that the enzyme has higher interactions of -4.1, -3.2, and -3.4 kJ/mol with common pesticides like dichloro diphenyl dichloro ethane(DDD), dichloro diphenyl trichloroe thane (DDT), and methyl parathion of organochlorides and organophosphates class. From the molecular dynamics simulation analysis, complex 1 (DDD -receptor) has been found to have the highest and more compact structure followed by methyl parathion -receptor. Prime MM-GBSA analysis also revealed the highest binding energy of -27.8 kcal/mol with the protein and DDD. Thus, it can be inferred from the current study that pesticides, particularly, DDD, DDT, and methyl parathion present in the coastal areas may have an impact on urease. This interaction can result in the inhibition of the urease activity of B. pasteurii, thus preventing the biomineralization process. This study would be the first report on the computational approach to understanding the interaction of prominent pesticides on the coastal region and B. pasteurii urease.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
| | - Krishna Kant Gupta
- School of Chemical and Biotechnology, SASTRA Deemed to be University, Thirumalasamudram, India
| | - Vinod Kumar Nathan
- School of Chemical and Biotechnology, SASTRA Deemed to be University, Thirumalasamudram, India
| |
Collapse
|
3
|
Dowerah D, V N Uppuladinne M, Paul S, Das D, Gour NK, Biswakarma N, Sarma PJ, Sonavane UB, Joshi RR, Ray SK, Deka RC. A Study Modeling Bridged Nucleic Acid-Based ASOs and Their Impact on the Structure and Stability of ASO/RNA Duplexes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:21407-21426. [PMID: 39370641 DOI: 10.1021/acs.langmuir.4c02171] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/08/2024]
Abstract
Antisense medications treat diseases that cannot be treated using traditional pharmacological technologies. Nucleotide monomers of bare and phosphorothioate (PS)-modified LNA, N-MeO-amino-BNA, 2',4'-BNANC[NH], 2',4'-BNANC[NMe], and N-Me-aminooxy-BNA antisense modifications were considered for a detailed DFT-based quantum chemical study to estimate their molecular-level structural and electronic properties. Oligomer hybrid duplex stability is described by performing an elaborate MD simulation study by incorporating the PS-LNA and PS-BNA antisense modifications onto 14-mer ASO/RNA hybrid gapmer type duplexes targeting protein PTEN mRNA nucleic acid sequence (5'-CTTAGCACTGGCCT-3'/3'-GAAUCGUGACCGGA-5'). Replica sets of MD simulations were performed accounting to two data sets, each set simulated for 1 μs simulation time. Bulk properties of oligomers are regulated by the chemical properties of their monomers. As such, the primary goal of this work focused on establishing an organized connection between the monomeric BNA nucleotide's electronic effects observed in DFT studies and the macroscopic behavior of the BNA antisense oligomers, as observed in MD simulations. The results from this study predicted that spatial orientation of MO-isosurfaces of the BNA nucleotides are concentrated in the nucleobase region. These BNA nucleotides may become less accessible for various electronic interactions when coupled as ASOs forming duplexes with target RNAs and when the ASO/RNA duplexes further bind with the RNase H. Understanding such electronic interactions is crucial to design superior antisense modifications with specific electronic properties. Also, for the particular nucleic acid sequence solvation of the duplexes although were higher compared to the natural oligonucleotides, their binding energies being relatively lower may lead to decreased antisense activity compared to existing analogs such as the LNAs and MOEs. Fine tuning these BNAs to obtain superior binding affinity is thus a necessity.
Collapse
Affiliation(s)
- Dikshita Dowerah
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Mallikarjunachari V N Uppuladinne
- HPC - Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Subrata Paul
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Department of Chemistry, Assam University, Silchar, Assam 788011, India
| | - Dharitri Das
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Nand K Gour
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Nishant Biswakarma
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Plaban J Sarma
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Department of Chemistry, Gargaon College, Simaluguri, Sivasagar, Assam 785686, India
| | - Uddhavesh B Sonavane
- HPC - Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Rajendra R Joshi
- HPC - Medical & Bioinformatics Applications Group, Centre for Development of Advanced Computing (C-DAC), Panchavati, Pashan, Pune 411008, India
| | - Suvendra K Ray
- Department of Molecular Biology and Biotechnology, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Center for Multidisciplinary Research, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| | - Ramesh Ch Deka
- CMML-Catalysis and Molecular Modelling Lab, Department of Chemical Sciences, Tezpur University, Napaam, Sonitpur, Assam 784028, India
- Center for Multidisciplinary Research, Tezpur University, Napaam, Sonitpur, Assam 784028, India
| |
Collapse
|
4
|
Xu K, Wang Z, Xiang S, Tang R, Deng Q, Ge J, Jiang Z, Yang K, Hou T, Sun H. Characterizing the Cooperative Effect of PROTAC Systems with End-Point Binding Free Energy Calculation. J Chem Inf Model 2024; 64:7666-7678. [PMID: 39361611 DOI: 10.1021/acs.jcim.4c01227] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2024]
Abstract
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.
Collapse
Affiliation(s)
- Kexin Xu
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Zhe Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058 Zhejiang, P. R. China
- School of Pharmacy, Hangzhou Normal University, Hangzhou 311121, P. R. China
| | - Sutong Xiang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Rongfan Tang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Qirui Deng
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Jingxuan Ge
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058 Zhejiang, P. R. China
| | - Zhiliang Jiang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Kaimo Yang
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058 Zhejiang, P. R. China
| | - Huiyong Sun
- Department of Medicinal Chemistry, China Pharmaceutical University, Nanjing 210009 Jiangsu, P. R. China
| |
Collapse
|
5
|
Pant P. Design and Characterization of Neutral Linker-Based Bis-Intercalator via Computer Simulations: Balancing DNA Binding and Cellular Uptake. Chem Biodivers 2024; 21:e202400768. [PMID: 38980964 DOI: 10.1002/cbdv.202400768] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Revised: 06/23/2024] [Accepted: 07/09/2024] [Indexed: 07/11/2024]
Abstract
Bis-intercalators refer to a class of chemical compounds known for their unique ability to simultaneously intercalate, or insert, into DNA at two distinct sites. These molecules typically feature two intercalating moieties connected by a linker, allowing them to engage with DNA base pairs at multiple locations. The bis-intercalation phenomenon plays a significant role in altering the DNA structure, affecting its stability, and potentially influencing various cellular processes. These compounds have gained considerable attention in medicinal chemistry and biochemistry due to their potential applications in cancer therapy, where they may interfere with DNA replication and transcription, leading to anticancer effects. Traditionally, these molecules often possess a high positive charge to enhance their affinity for the negatively charged DNA. However, due to a high positive charge, their cellular uptake is compromised, along with their enhanced potential off-target effects. In this study, we utilized bis-intercalator TOTO and replaced the charged linker segment (propane-1,3-diammonium) with a neutral peroxodisulphuric acid linker. Using molecular modeling and computer simulations (500 ns, 3 replicas), we investigated the potential of the designed molecule as a bis-intercalator and compared the properties with the control bis-intercalator bound to DNA. We observed that the designed bis-intercalator exhibited improved DNA binding (as assessed through MM-PBSA and Delphi methods) and membrane translocation permeability. With an overall reduced charge, significantly less off-target binding of the designed molecule is also anticipated. Consequently, bis-intercalators based on peroxodisulphuric linkers can potentially target DNA effectively, and their role in the future design of bis-intercalators is foreseen.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Biotechnology, School of Engineering and Applied Sciences, Bennett University, Greater Noida, UP, India
| |
Collapse
|
6
|
Parvin R, Habib Ullah Masum M, Ferdous J, Mahdeen AA, Shafiqul Islam Khan M. Designing of a chimeric multiepitope vaccine against bancroftian lymphatic filariasis through immunoinformatics approaches. PLoS One 2024; 19:e0310398. [PMID: 39298468 DOI: 10.1371/journal.pone.0310398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2024] [Accepted: 09/01/2024] [Indexed: 09/21/2024] Open
Abstract
The filarial worms of Wuchereria bancrofti are the primary cause of lymphatic filariasis (LF), a mosquito-borne disease among the neglected tropical parasitic diseases. Considering the global endemic consequences of the disease, there is a need to develop a successful vaccine candidate against LF. Using advanced immunoinformatics approaches, we designed two multiepitope vaccines targeting W. bancrofti's glutathione S-transferase and thioredoxin. Therefore, we predicted several MHC-1, MHC-2, and B-cell epitopes from these proteins and mapped two vaccine candidates (V1 and V2). The vaccines were subsequently employed for physicochemical analysis, structural prediction and validation, docking and normal mode analysis, codon optimization, and immune simulation. The selected MHC-1, MHC-2, and B-cell epitopes were antigenic without allergenicity or toxicity. The designed vaccines were expected to be soluble, stable proteins under physiological conditions. Compared to V2, V1's secondary and tertiary structures were simultaneously favorable, with Ramachandran plot analysis revealing 95.6% residues in favored areas. Subsequently, the molecular docking analysis indicated that the V1 had a high binding affinity for the TLR-2, TLR-4 and TLR-5, as suggested by the docking scores of -1248.7, -1038.5 and -1562.8, respectively. The NMA of these complexes further indicated their structural flexibility. Molecular dynamics simulations of V1-TLR complexes revealed V1-TLR-4 as the most stable, with the lowest free energy and minimal fluctuations, indicating the strongest binding affinity. The results of the codon optimization showed high levels of expression, with a favorable CAI score (<1.0). A three-dose vaccination analysis showed significant and persistent immunological responses, including adaptive and innate immune responses. The findings emphasize the potential of the V1 against W. bancrofti, but further validation is required through in vitro, in vivo, and clinical trials.
Collapse
Affiliation(s)
- Rehana Parvin
- Department of Pathology and Parasitology, Faculty of Veterinary Medicine, Chattogram Veterinary and Animal Sciences University (CVASU), Chattogram, Bangladesh
| | - Md Habib Ullah Masum
- Department of Genomics and Bioinformatics, Faculty of Biotechnology and Genetic Engineering, Chattogram Veterinary and Animal Sciences University (CVASU), Chattogram, Bangladesh
| | - Jannatul Ferdous
- Department of Obstetrics and Gynecology, Chittagong Medical College, Chittagong, Bangladesh
| | - Ahmad Abdullah Mahdeen
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | - Md Shafiqul Islam Khan
- Department of Cellular and Molecular Biology, Faculty of Biotechnology and Genetic Engineering, Chattogram Veterinary and Animal Sciences University (CVASU), Chattogram, Bangladesh
| |
Collapse
|
7
|
Wang P, Song X, Liang Q. Study on the Inhibitory Effect of Bioactive Peptides Derived from Yak Milk Cheese on Cholesterol Esterase. Foods 2024; 13:2970. [PMID: 39335898 PMCID: PMC11431439 DOI: 10.3390/foods13182970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2024] [Revised: 09/14/2024] [Accepted: 09/18/2024] [Indexed: 09/30/2024] Open
Abstract
The bioactive peptides derived from yak milk cheese exhibited cholesterol-lowering properties. However, there was limited research on their inhibitory effects on cholesterol esterase (CE) and elucidation of their potential inhibitory mechanisms. In this study, we identified CE-inhibiting peptides through virtual screening and in vitro assays. Additionally, molecular docking and molecular dynamics studies were conducted to explore the mechanisms. The results indicated that peptides RK7 (RPKHPIK), KQ7 (KVLPVPQ), QP13 (QEPVLGPVRGPFP), TL9 (TPVVVPPFL), VN10 (VYPFPGPIPN), LQ10 (LPPTVMFPPQ), and SN12 (SLVYPFPGPIPN) possessed molecular weights of less than 1.5 kDa and a high proportion of hydrophobic amino acids, demonstrating notable inhibitory effects on CE. Molecular docking and dynamics revealed that peptides RK7, KQ7, QP13, and VN10 bound to key amino acid residues Arg423, His435, and Ser422 of CE through hydrogen bonds, hydrophobic interactions, salt bridges, and π-π stacking, occupying the substrate-binding site and exerting inhibitory effects on CE. The four peptides were further synthesized to verify their CE-inhibitory effects in vitro. RK7, KQ7, QP13, and VN10 exhibited inhibitory activity on CE with IC50 values of 8.16 × 10-7 mol/L, 8.10 × 10-7 mol/L, 4.63 × 10-7 mol/L, and 7.97 × 10-7 mol/L; RK7, KQ7, QP13, and VN10 were effective in inhibiting CE after simulated gastrointestinal digestion, especially with a significant increase in the inhibitory activity of KQ7 and RK7, respectively. Our findings suggested that bioactive peptides from yak milk cheese represented a novel class of potential CE inhibitors.
Collapse
Affiliation(s)
- Peng Wang
- Functional Dairy Products Engineering Laboratory of Gansu Province, College of Food Science and Engineering, Gansu Agricultural University, Lanzhou 730070, China
| | - Xuemei Song
- Functional Dairy Products Engineering Laboratory of Gansu Province, College of Food Science and Engineering, Gansu Agricultural University, Lanzhou 730070, China
| | - Qi Liang
- Functional Dairy Products Engineering Laboratory of Gansu Province, College of Food Science and Engineering, Gansu Agricultural University, Lanzhou 730070, China
| |
Collapse
|
8
|
Pant P. Flexible RNA aptamers as inhibitors of Bacillus anthracis ribosomal protein S8: Insights from molecular dynamics simulations. Biophys Chem 2024; 312:107273. [PMID: 38850843 DOI: 10.1016/j.bpc.2024.107273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Accepted: 06/02/2024] [Indexed: 06/10/2024]
Abstract
Bacillus anthracis, the causative agent of anthrax, poses a substantial threat to public health and national security, and is recognized as a potential bioweapon due to its capacity to form resilient spores with enduring viability. Inhalation or ingestion of even minute quantities of aerosolized spores can lead to widespread illness and fatalities, underscoring the formidable lethality of the bacterium. With an untreated mortality rate of 100%, Bacillus anthracis is a disconcerting candidate for bioterrorism. In response to this critical scenario, we employed state-of-the-art computational tools to conceive and characterize flexible RNA aptamer therapeutics tailored for anthrax. The foundational structure of the flexible RNA aptamers was designed by removing the C2'-C3' in each nucleotide unit. Leveraging the crystal structure of Bacillus anthracis ribosomal protein S8 complexed with an RNA aptamer, we explored the structural, dynamic, and energetic aspects of the modified RNA aptamer - S8 protein complexes through extensive all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicas each), followed by drawing comparisons to the control system. Our findings demonstrate the enhanced binding competencies of the flexible RNA aptamers to the S8 protein via better shape complementarity and improved H-bond network compared to the control RNA aptamer. This research offers valuable insights into the development of RNA aptamer therapeutics targeting Bacillus anthracis, paving the way for innovative strategies to mitigate the impact of this formidable pathogen.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Biotechnology, School of Engineering and Applied Sciences, Bennett University, Greater Noida, U.P., India.
| |
Collapse
|
9
|
da Fonseca AM, Caluaco BJ, Madureira JMC, Cabongo SQ, Gaieta EM, Djata F, Colares RP, Neto MM, Fernandes CFC, Marinho GS, Dos Santos HS, Marinho ES. Screening of Potential Inhibitors Targeting the Main Protease Structure of SARS-CoV-2 via Molecular Docking, and Approach with Molecular Dynamics, RMSD, RMSF, H-Bond, SASA and MMGBSA. Mol Biotechnol 2024; 66:1919-1933. [PMID: 37490200 DOI: 10.1007/s12033-023-00831-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Accepted: 07/17/2023] [Indexed: 07/26/2023]
Abstract
Severe Acute Respiratory Syndrome caused by a coronavirus is a recent viral infection. There is no scientific evidence or clinical trials to indicate that possible therapies have demonstrated results in suspected or confirmed patients. This work aims to perform a virtual screening of 1430 ligands through molecular docking and to evaluate the possible inhibitory capacity of these drugs about the Mpro protease of Covid-19. The selected drugs were registered with the FDA and available in the virtual drug library, widely used by the population. The simulation was performed using the MolAiCalD algorithm, with a Lamarckian genetic model (GA) combined with energy estimation based on rigid and flexible conformation grids. In addition, molecular dynamics studies were also performed to verify the stability of the receptor-ligand complexes formed through analyses of RMSD, RMSF, H-Bond, SASA, and MMGBSA. Compared to the binding energy of the synthetic redocking coupling (-6.8 kcal/mol/RMSD of 1.34 Å), which was considerably higher, it was then decided to analyze the parameters of only three ligands: ergotamine (-9.9 kcal/mol/RMSD of 2.0 Å), dihydroergotamine (-9.8 kcal/mol/RMSD of 1.46 Å) and olysio (-9.5 kcal/mol/RMSD of 1.5 Å). It can be stated that ergotamine showed the best interactions with the Mpro protease of Covid-19 in the in silico study, showing itself as a promising candidate for treating Covid-19.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Mestrado Acadêmico em Sociobiodiversidades e Tecnologias Sustentáveis - MASTS, Instituto de Engenharias e Desenvolvimento Sustentável, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Bernardino Joaquim Caluaco
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | | | - Sadrack Queque Cabongo
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Eduardo Menezes Gaieta
- Fundação Oswaldo Cruz - Fiocruz, R. São José, S/N - Precabura, Eusébio, Ceará, 61773-270, Brazil
| | - Faustino Djata
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Regilany Paulo Colares
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Moises Maia Neto
- Curso de Graduação Em Farmácia, Centro Universitário Fametro, Fortaleza, CE, Brazil
| | | | - Gabrielle Silva Marinho
- Faculdade de Filosofia, Dom Aureliano Matos - FAFIDAM, Universidade Estadual Do Ceará, Centro, Limoeiro Do Norte, CE, Brazil
| | | | - Emmanuel Silva Marinho
- Faculdade de Filosofia, Dom Aureliano Matos - FAFIDAM, Universidade Estadual Do Ceará, Centro, Limoeiro Do Norte, CE, Brazil
| |
Collapse
|
10
|
Masum MHU, Wajed S, Hossain MI, Moumi NR, Talukder A, Rahman MM. An mRNA vaccine for pancreatic cancer designed by applying in silico immunoinformatics and reverse vaccinology approaches. PLoS One 2024; 19:e0305413. [PMID: 38976715 PMCID: PMC11230540 DOI: 10.1371/journal.pone.0305413] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2024] [Accepted: 05/30/2024] [Indexed: 07/10/2024] Open
Abstract
Pancreatic ductal adenocarcinoma is the most prevalent pancreatic cancer, which is considered a significant global health concern. Chemotherapy and surgery are the mainstays of current pancreatic cancer treatments; however, a few cases are suitable for surgery, and most of the cases will experience recurrent episodes. Compared to DNA or peptide vaccines, mRNA vaccines for pancreatic cancer have more promise because of their delivery, enhanced immune responses, and lower proneness to mutation. We constructed an mRNA vaccine by analyzing S100 family proteins, which are all major activators of receptors for advanced glycation end products. We applied immunoinformatic approaches, including physicochemical properties analysis, structural prediction and validation, molecular docking study, in silico cloning, and immune simulations. The designed mRNA vaccine was estimated to have a molecular weight of 165023.50 Da and was highly soluble (grand average of hydropathicity of -0.440). In the structural assessment, the vaccine seemed to be a well-stable and functioning protein (Z score of -8.94). Also, the docking analysis suggested that the vaccine had a high affinity for TLR-2 and TLR-4 receptors. Additionally, the molecular mechanics with generalized Born and surface area solvation analysis of the "Vaccine-TLR-2" (-141.07 kcal/mol) and "Vaccine-TLR-4" (-271.72 kcal/mol) complexes also suggests a strong binding affinity for the receptors. Codon optimization also provided a high expression level with a GC content of 47.04% and a codon adaptation index score 1.0. The appearance of memory B-cells and T-cells was also observed over a while, with an increased level of helper T-cells and immunoglobulins (IgM and IgG). Moreover, the minimum free energy of the mRNA vaccine was predicted at -1760.00 kcal/mol, indicating the stability of the vaccine following its entry, transcription, and expression. This hypothetical vaccine offers a groundbreaking tool for future research and therapeutic development of pancreatic cancer.
Collapse
Affiliation(s)
- Md Habib Ullah Masum
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
- Microbiology, Cancer and Bioinformatics Research Group, Noakhali Science and Technology University, Noakhali, Bangladesh
| | - Shah Wajed
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
- Microbiology, Cancer and Bioinformatics Research Group, Noakhali Science and Technology University, Noakhali, Bangladesh
- Infectiology: Biology of Infectious Diseases, Universite Paris-Saclay, Gif-sur-Yvette, France
| | - Md Imam Hossain
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
| | - Nusrat Rahman Moumi
- Medical Sciences, University of Central Lancashire, Preston, Lancashire, United Kingdom
| | - Asma Talukder
- Microbiology, Cancer and Bioinformatics Research Group, Noakhali Science and Technology University, Noakhali, Bangladesh
- Department of Biotechnology and Genetic Engineering, Noakhali Science and Technology University, Noakhali, Bangladesh
- School of Pharmacy and Medical Sciences, and Menzies Health Institute Queensland, Griffith University, Brisbane, Queensland, Australia
| | - Md Mijanur Rahman
- Department of Microbiology, Noakhali Science and Technology University, Noakhali, Bangladesh
- Microbiology, Cancer and Bioinformatics Research Group, Noakhali Science and Technology University, Noakhali, Bangladesh
- School of Pharmacy and Medical Sciences, and Menzies Health Institute Queensland, Griffith University, Brisbane, Queensland, Australia
| |
Collapse
|
11
|
Silva de Freitas Cesário HP, das Chagas Lima Pinto F, Marques Canuto K, Rocha Silveira E, Veras Wilke D, Gois Ferreira E, Marques da Fonseca A, Alves de Vasconcelos M, Teixeira EH, Deusdênia Loiola Pessoa O. Further Polycyclic Quinones of Micromonospora sp. Chem Biodivers 2024; 21:e202301771. [PMID: 38628065 DOI: 10.1002/cbdv.202301771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Accepted: 03/31/2024] [Indexed: 06/01/2024]
Abstract
The crude acetone extract of a marine Micromonospora sp. strain associated with Eudistoma vannnamei was fractioned with hexane and ethyl acetate. The crude extract and both soluble fractions were assayed against several bacteria strains. The new polycyclic quinones 12-hydroxy-9-propyltetracene-6,1-dione (1), 5,12-dihydroxy-4-methoxy-9-propyltetracene-5,12-dione (2), and 4,6-dihydroxy-3-methoxycarbonyl- methyl-6a-(oxobutyl)-5,12-anthraquinone (3), along with the known 4,6-dihydroxy-3-methoxycarbonyl-methyl-6a-(oxo-3-methyl-butyl)-5,12-anthraquinone (4) and 4,6-dihydroxy-3-methoxycarbonyl-methyl-6a-(oxopentyl)-5,12-anthraquinone (5) were isolated from the hexane-soluble fraction, while from the active ethyl acetate fraction were isolated the known 4,6,11-trihydroxy-9-propyltetracene-5,12-dione (6), 4-methoxy-9-propyltetracene-6,11-dione (7), 7,8,9,10-tetrahydro-9-hydroxy-4-methoxy-9-propyltetracene-6,11-dione (8), and 10β-carbomethoxy-7,8,9,10-tetrahydro-4,6,7α,9α,11-pentahydroxy-9-propyltetracene-5,12-dione (9). The structures of the new compounds were established by interpretation of HRMS and NMR techniques. A study of molecular docking was performed with the compounds from the active ethyl acetate fraction to correlate tentatively with the antimicrobial activity. Molecular docking, RMSD, RMSF, and MM-GBSA evaluations were performed to investigate the inhibitory activity of 6-8 against the protein PDB-codex 1MWT, being considered a promising target for studying drug development responsible for inhibiting replication of Staphylococcus aureus. Penicillin G was used as the standard inhibitory. Anthracyclinones 6-8 were the best hydrolase inhibitor with affinity energy -8.1 to -7.9 kcal/mol compared to penicillin G, which presented -6.9 kcal/mol. Both 8 and 7 present potent inhibitory effects against hydrolase through molecular dynamics simulation and exhibit favorable drug-like properties, promising new hydrolase blockers to fight bacterial infections from Staphylococcus aureus.
Collapse
Affiliation(s)
| | - Francisco das Chagas Lima Pinto
- Institute of Exact and Natural Sciences, University of International Integration of Afro-Brazilian Lusofonia, 62785-000, Acarape, CE, Brazil
| | | | - Ediberto Rocha Silveira
- Department of Organic and Inorganic Chemistry, Science Center, Federal University of Ceará, Fortaleza, CE, 60455-760, Brazil
| | - Diego Veras Wilke
- Department of Physiology and Pharmacology, Federal University of Ceará, 60165-085, Fortaleza, Ceará, Brazil
| | - Elthon Gois Ferreira
- Department of Physiology and Pharmacology, Federal University of Ceará, 60165-085, Fortaleza, Ceará, Brazil
| | - Aluísio Marques da Fonseca
- Institute of Exact and Natural Sciences, University of International Integration of Afro-Brazilian Lusofonia, 62785-000, Acarape, CE, Brazil
| | - Mayron Alves de Vasconcelos
- Integrated Laboratory of Biomolecules (LIBS), Department of Pathology and Legal Medicine, Federal University of Ceará, Fortaleza, CE, 62042-280, Brazil
- Faculty of Education of Itapipoca (FACEDI), State University of Ceará, Itapipoca, CE, 62500-000, Brazil
| | - Edson Holanda Teixeira
- Faculty of Education of Itapipoca (FACEDI), State University of Ceará, Itapipoca, CE, 62500-000, Brazil
| | - Otilia Deusdênia Loiola Pessoa
- Department of Organic and Inorganic Chemistry, Science Center, Federal University of Ceará, Fortaleza, CE, 60455-760, Brazil
| |
Collapse
|
12
|
Quitté L, Leclercq M, Prunier J, Scott-Boyer MP, Moroy G, Droit A. A Machine Learning Approach to Identify Key Residues Involved in Protein-Protein Interactions Exemplified with SARS-CoV-2 Variants. Int J Mol Sci 2024; 25:6535. [PMID: 38928241 PMCID: PMC11204244 DOI: 10.3390/ijms25126535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2024] [Revised: 05/24/2024] [Accepted: 06/03/2024] [Indexed: 06/28/2024] Open
Abstract
Human infection with the coronavirus disease 2019 (COVID-19) is mediated by the binding of the spike protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to the human angiotensin-converting enzyme 2 (ACE2). The frequent mutations in the receptor-binding domain (RBD) of the spike protein induced the emergence of variants with increased contagion and can hinder vaccine efficiency. Hence, it is crucial to better understand the binding mechanisms of variant RBDs to human ACE2 and develop efficient methods to characterize this interaction. In this work, we present an approach that uses machine learning to analyze the molecular dynamics simulations of RBD variant trajectories bound to ACE2. Along with the binding free energy calculation, this method was used to characterize the major differences in ACE2-binding capacity of three SARS-CoV-2 RBD variants-namely the original Wuhan strain, Omicron BA.1, and the more recent Omicron BA.5 sublineages. Our analyses assessed the differences in binding free energy and shed light on how it affects the infectious rates of different variants. Furthermore, this approach successfully characterized key binding interactions and could be deployed as an efficient tool to predict different binding inhibitors to pave the way for new preventive and therapeutic strategies.
Collapse
Affiliation(s)
- Léopold Quitté
- Centre de Recherche du CHU de Québec, Université Laval, Québec, QC G1V 0A6, Canada; (L.Q.); (M.L.); (J.P.); (M.-P.S.-B.)
| | - Mickael Leclercq
- Centre de Recherche du CHU de Québec, Université Laval, Québec, QC G1V 0A6, Canada; (L.Q.); (M.L.); (J.P.); (M.-P.S.-B.)
| | - Julien Prunier
- Centre de Recherche du CHU de Québec, Université Laval, Québec, QC G1V 0A6, Canada; (L.Q.); (M.L.); (J.P.); (M.-P.S.-B.)
| | - Marie-Pier Scott-Boyer
- Centre de Recherche du CHU de Québec, Université Laval, Québec, QC G1V 0A6, Canada; (L.Q.); (M.L.); (J.P.); (M.-P.S.-B.)
| | - Gautier Moroy
- Université Paris Cité, CNRS, INSERM, Unité de Biologie Fonctionnelle et Adaptative, F-75013 Paris, France
| | - Arnaud Droit
- Centre de Recherche du CHU de Québec, Université Laval, Québec, QC G1V 0A6, Canada; (L.Q.); (M.L.); (J.P.); (M.-P.S.-B.)
- Département de Médecine Moléculaire, Université Laval, Québec, QC G1V 0A6, Canada
| |
Collapse
|
13
|
Moltrasio C, Silva CA, Tricarico PM, Marzano AV, Sueleman M, Crovella S. Biosensing circulating MicroRNAs in autoinflammatory skin diseases: Focus on Hidradenitis suppurativa. Front Genet 2024; 15:1383452. [PMID: 38655054 PMCID: PMC11035790 DOI: 10.3389/fgene.2024.1383452] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2024] [Accepted: 03/25/2024] [Indexed: 04/26/2024] Open
Abstract
MicroRNAs (miRNAs) play a crucial role in the early diagnosis of autoinflammatory diseases, with Hidradenitis Suppurativa (HS) being a notable example. HS, an autoinflammatory skin disease affecting the pilosebaceous unit, profoundly impacts patients' quality of life. Its hidden nature, with insidious initial symptoms and patient reluctance to seek medical consultation, often leads to a diagnostic delay of up to 7 years. Recognizing the urgency for early diagnostic tools, recent research identified significant differences in circulating miRNA expression, including miR-24-1-5p, miR-146a-5p, miR26a-5p, miR-206, miR338-3p, and miR-338-5p, between HS patients and healthy controls. These miRNAs serve as potential biomarkers for earlier disease detection. Traditional molecular biology techniques, like reverse transcription quantitative-polymerase chain reaction (RT-qPCR), are employed for their detection using specific primers and probes. Alternatively, short peptides offer a versatile and effective means for capturing miRNAs, providing specificity, ease of synthesis, stability, and multiplexing potential. In this context, we present a computational simulation pipeline designed for crafting peptide sequences that can capture circulating miRNAs in the blood of patients with autoinflammatory skin diseases, including HS. This innovative approach aims to expedite early diagnosis and enhance therapeutic follow-up, addressing the critical need for timely intervention in HS and similar conditions.
Collapse
Affiliation(s)
- Chiara Moltrasio
- Dermatology Unit, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico, Milan, Italy
| | | | - Paola Maura Tricarico
- Department of Advanced Diagnostics, Institute for Maternal and Child Health-IRCCS Burlo Garofolo, Trieste, Italy
| | - Angelo Valerio Marzano
- Dermatology Unit, Fondazione IRCCS Ca’ Granda Ospedale Maggiore Policlinico, Milan, Italy
- Department of Pathophysiology and Transplantation, Università degli Studi di Milano, Milan, Italy
| | | | - Sergio Crovella
- Laboratory of Animal Research (LARC), Qatar University, Doha, Qatar
| |
Collapse
|
14
|
Qiu X, Zhang Q, Li Z, Zhang J, Liu H. Revealing the Interaction Mechanism between Mycobacterium tuberculosis GyrB and Novobiocin, SPR719 through Binding Thermodynamics and Dissociation Kinetics Analysis. Int J Mol Sci 2024; 25:3764. [PMID: 38612573 PMCID: PMC11011931 DOI: 10.3390/ijms25073764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 03/22/2024] [Accepted: 03/24/2024] [Indexed: 04/14/2024] Open
Abstract
With the rapid emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), various levels of resistance against existing anti-tuberculosis (TB) drugs have developed. Consequently, the identification of new anti-TB targets and drugs is critically urgent. DNA gyrase subunit B (GyrB) has been identified as a potential anti-TB target, with novobiocin and SPR719 proposed as inhibitors targeting GyrB. Therefore, elucidating the molecular interactions between GyrB and its inhibitors is crucial for the discovery and design of efficient GyrB inhibitors for combating multidrug-resistant TB. In this study, we revealed the detailed binding mechanisms and dissociation processes of the representative inhibitors, novobiocin and SPR719, with GyrB using classical molecular dynamics (MD) simulations, tau-random acceleration molecular dynamics (τ-RAMD) simulations, and steered molecular dynamics (SMD) simulations. Our simulation results demonstrate that both electrostatic and van der Waals interactions contribute favorably to the inhibitors' binding to GyrB, with Asn52, Asp79, Arg82, Lys108, Tyr114, and Arg141 being key residues for the inhibitors' attachment to GyrB. The τ-RAMD simulations indicate that the inhibitors primarily dissociate from the ATP channel. The SMD simulation results reveal that both inhibitors follow a similar dissociation mechanism, requiring the overcoming of hydrophobic interactions and hydrogen bonding interactions formed with the ATP active site. The binding and dissociation mechanisms of GyrB with inhibitors novobiocin and SPR719 obtained in our work will provide new insights for the development of promising GyrB inhibitors.
Collapse
Affiliation(s)
- Xiaofei Qiu
- School of Pharmacy, Lanzhou University, Lanzhou 730000, China; (X.Q.); (Z.L.); (J.Z.)
| | - Qianqian Zhang
- Faculty of Applied Science, Macao Polytechnic University, Macao SAR, China;
| | - Zhaoguo Li
- School of Pharmacy, Lanzhou University, Lanzhou 730000, China; (X.Q.); (Z.L.); (J.Z.)
| | - Juan Zhang
- School of Pharmacy, Lanzhou University, Lanzhou 730000, China; (X.Q.); (Z.L.); (J.Z.)
| | - Huanxiang Liu
- Faculty of Applied Science, Macao Polytechnic University, Macao SAR, China;
| |
Collapse
|
15
|
Jiang D, Du H, Zhao H, Deng Y, Wu Z, Wang J, Zeng Y, Zhang H, Wang X, Wang E, Hou T, Hsieh CY. Assessing the performance of MM/PBSA and MM/GBSA methods. 10. Prediction reliability of binding affinities and binding poses for RNA-ligand complexes. Phys Chem Chem Phys 2024; 26:10323-10335. [PMID: 38501198 DOI: 10.1039/d3cp04366e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
Ribonucleic acid (RNA)-ligand interactions play a pivotal role in a wide spectrum of biological processes, ranging from protein biosynthesis to cellular reproduction. This recognition has prompted the broader acceptance of RNA as a viable candidate for drug targets. Delving into the atomic-scale understanding of RNA-ligand interactions holds paramount importance in unraveling intricate molecular mechanisms and further contributing to RNA-based drug discovery. Computational approaches, particularly molecular docking, offer an efficient way of predicting the interactions between RNA and small molecules. However, the accuracy and reliability of these predictions heavily depend on the performance of scoring functions (SFs). In contrast to the majority of SFs used in RNA-ligand docking, the end-point binding free energy calculation methods, such as molecular mechanics/generalized Born surface area (MM/GBSA) and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA), stand as theoretically more rigorous approaches. Yet, the evaluation of their effectiveness in predicting both binding affinities and binding poses within RNA-ligand systems remains unexplored. This study first reported the performance of MM/PBSA and MM/GBSA with diverse solvation models, interior dielectric constants (εin) and force fields in the context of binding affinity prediction for 29 RNA-ligand complexes. MM/GBSA is based on short (5 ns) molecular dynamics (MD) simulations in an explicit solvent with the YIL force field; the GBGBn2 model with higher interior dielectric constant (εin = 12, 16 or 20) yields the best correlation (Rp = -0.513), which outperforms the best correlation (Rp = -0.317, rDock) offered by various docking programs. Then, the efficacy of MM/GBSA in identifying the near-native binding poses from the decoys was assessed based on 56 RNA-ligand complexes. However, it is evident that MM/GBSA has limitations in accurately predicting binding poses for RNA-ligand systems, particularly compared with notably proficient docking programs like rDock and PLANTS. The best top-1 success rate achieved by MM/GBSA rescoring is 39.3%, which falls below the best results given by docking programs (50%, PLNATS). This study represents the first evaluation of MM/PBSA and MM/GBSA for RNA-ligand systems and is expected to provide valuable insights into their successful application to RNA targets.
Collapse
Affiliation(s)
- Dejun Jiang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Hongyan Du
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Huifeng Zhao
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Yafeng Deng
- Hangzhou Carbonsilicon AI Technology Co., Ltd, Hangzhou, Zhejiang 310018, China
| | - Zhenxing Wu
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Jike Wang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Yundian Zeng
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Haotian Zhang
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Xiaorui Wang
- China State Key Laboratory of Quality Research in Chinese Medicines, Macau University of Science and Technology, Macau 999078, China
| | - Ercheng Wang
- Zhejiang Laboratory, Hangzhou, Zhejiang 311100, China.
| | - Tingjun Hou
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| | - Chang-Yu Hsieh
- Innovation Institute for Artificial Intelligence in Medicine of Zhejiang University, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China.
| |
Collapse
|
16
|
Meng T, Wang Z, Zhang H, Zhao Z, Huang W, Xu L, Liu M, Li J, Yan H. In Silico Investigations on the Synergistic Binding Mechanism of Functional Compounds with Beta-Lactoglobulin. Molecules 2024; 29:956. [PMID: 38474468 DOI: 10.3390/molecules29050956] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2023] [Revised: 02/03/2024] [Accepted: 02/20/2024] [Indexed: 03/14/2024] Open
Abstract
Piceatannol (PIC) and epigallocatechin gallate (EGCG) are polyphenolic compounds with applications in the treatment of various diseases such as cancer, but their stability is poor. β-lactoglobulin (β-LG) is a natural carrier that provides a protective effect to small molecule compounds and thus improves their stability. To elucidate the mechanism of action of EGCG, PIC, and palmitate (PLM) in binding to β-LG individually and jointly, this study applied molecular docking and molecular dynamics simulations combined with in-depth analyses including noncovalent interaction (NCI) and binding free energy to investigate the binding characteristics between β-LG and compounds of PIC, EGCG, and PLM. Simulations on the binary complexes of β-LG + PIC, β-LG + EGCG, and β-LG + PLM and ternary complexes of (β-LG + PLM) + PIC, (β-LG + PLM) + EGCG, β-LG + PIC) + EGCG, and (β-LG + EGCG) + PIC were performed for comparison and characterizing the interactions between binding compounds. The results demonstrated that the co-bound PIC and EGCG showed non-beneficial effects on each other. However, the centrally located PLM was revealed to be able to adjust the binding conformation of PIC, which led to the increase in binding affinity with β-LG, thus showing a synergistic effect on the co-bound PIC. The current study of β-LG co-encapsulated PLM and PIC provides a theoretical basis and research suggestions for improving the stability of polyphenols.
Collapse
Affiliation(s)
- Tong Meng
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng 252059, China
| | - Zhiguo Wang
- Institute of Ageing Research, School of Basic Medical Sciences, Hangzhou Normal University, Hangzhou 311121, China
| | - Hao Zhang
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng 252059, China
| | - Zhen Zhao
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng 252059, China
| | - Wanlin Huang
- School of Basic Medical Sciences, Hangzhou Normal University, Hangzhou 311121, China
| | - Liucheng Xu
- School of Basic Medical Sciences, Hangzhou Normal University, Hangzhou 311121, China
| | - Min Liu
- School of Chemistry and Chemical Engineering, Liaocheng University, Liaocheng 252059, China
| | - Jun Li
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng 252059, China
| | - Hui Yan
- School of Pharmaceutical Sciences, Liaocheng University, Liaocheng 252059, China
| |
Collapse
|
17
|
Hu Z, Zeng Y, Zhang Y, Zhang Q, Xu J, Liu L. Discovery of small molecule c-Maf inhibitors using molecular docking-based virtual screening, molecular dynamics simulation, and biological evaluation. Chem Biol Drug Des 2024; 103:e14403. [PMID: 37984986 DOI: 10.1111/cbdd.14403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 10/18/2023] [Accepted: 11/08/2023] [Indexed: 11/22/2023]
Abstract
Multiple myeloma (MM) is a prevalent plasma cell malignancy in the blood system that remains incurable. Given the abnormally high expression of c-Maf in most MM patients, targeting c-Maf presents an attractive therapeutic approach for treating MM malignancies. In this study, we employed a combined strategy involving molecular docking-based virtual screening, molecular dynamics (MD) simulation, and molecular mechanics/generalized Born surface area (MM/GBSA) free energy calculation on existing FDA-approved drugs. Six compounds were selected for further experimental assay: vemurafenib, sorafenib, sildenafil, fluvastatin, erlotinib, and glimepiride. Among these compounds, sorafenib and glimepiride exhibited significant inhibition of myeloma cell proliferation in the RPMI-8226 cell line. Moreover, both compounds simultaneously downregulated c-Maf protein expression to induce G1 phase arrest and apoptosis in myeloma cells. Collectively, sorafenib and glimepiride may be considered promising candidates for developing more potent c-Maf inhibitors in the future.
Collapse
Affiliation(s)
- Zhiwei Hu
- College of Medical Imaging, Xuzhou Medical University, Xuzhou, Jiangsu, China
| | - Yindi Zeng
- Blood Diseases Institute, Xuzhou Medical University, Xuzhou, Jiangsu, China
| | - Yaxin Zhang
- Blood Diseases Institute, Xuzhou Medical University, Xuzhou, Jiangsu, China
| | - Qiurong Zhang
- The Second Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, China
| | - Jinge Xu
- The Second Affiliated Hospital of Xuzhou Medical University, Xuzhou, Jiangsu, China
| | - Linlin Liu
- College of Medical Imaging, Xuzhou Medical University, Xuzhou, Jiangsu, China
| |
Collapse
|
18
|
Pandya K, Singh N. In silico study reveals unconventional interactions between MDC1 of DDR and Beclin-1 of autophagy. Mol Divers 2023; 27:2789-2802. [PMID: 36482226 DOI: 10.1007/s11030-022-10579-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2022] [Accepted: 12/02/2022] [Indexed: 12/14/2022]
Abstract
DNA damage response (DDR) and autophagy are concerned with maintaining cellular homeostasis and dysregulation of these two pathways lead to pathologic conditions including tumorigenesis. Autophagy is activated as a protective mechanism during DDR which is indicative of their functional cooperativity but the molecular mechanism leading to the convergence of these two pathways during genotoxic stress remains elusive. In this study, through in silico analysis, we have shown an interaction between the Mediator of DNA damage checkpoint 1 (MDC1), an important DDR-associated protein, and Beclin-1, an autophagy inducer. MDC1 is an adaptor or scaffold protein known to regulate DDR, apoptosis, and cell cycle progression. While, Beclin-1 is involved in autophagosome nucleation and exhibits affinity for binding to Fork-head-associated domain (FHA) containing proteins. The FHA domain is commonly conserved in DDR-related proteins including MDC1. Through molecular docking, we have predicted the modeled complex between the MDC1 FHA domain and the Beclin-1 Coiled coil domain (CCD). The docking complex was modeled using ClusPro2.0, based on the crystal structure for the dimerized MDC1 FHA domain and Beclin-1 CCD. The complex stability and binding affinities were assessed using a Ramachandran plot, MD simulation, MM/GBSA, and PRODIGY webserver. Finally, the hot-spot residues at the interface were determined using computational alanine scanning by the DrugScorePPI webserver. Our analysis unveils significant interaction between MDC1 and Beclin-1, involving hydrogen bonds, non-bonded contacts, and salt bridges and indicates MDC1 possibly recruits Beclin-1 to the DSBs, as a consequence of which Beclin-1 is able to modulate DDR.
Collapse
Affiliation(s)
- Kavya Pandya
- Department of Biotechnology and Bioengineering, Indian Institute of Advanced Research, Gandhinagar, India
| | - Neeru Singh
- Department of Biotechnology and Bioengineering, Indian Institute of Advanced Research, Gandhinagar, India.
| |
Collapse
|
19
|
da Fonseca AM, Cabongo SQ, Caluaco BJ, Colares RP, Fernandes CFC, Dos Santos HS, de Lima-Neto P, Marinho ES. The search for new efficient inhibitors of SARS-COV-2 through the De novo drug design developed by artificial intelligence. J Biomol Struct Dyn 2023; 41:9890-9906. [PMID: 36420665 DOI: 10.1080/07391102.2022.2148128] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Accepted: 11/10/2022] [Indexed: 11/25/2022]
Abstract
The pandemic caused by Sars-CoV-2 is a viral infection that has generated one of the most significant health problems worldwide. Previous studies report the main protease (Mpro) as a potential target for this virus, as it is considered a crucial enzyme in mediating replication and viral transcription. This work presented the construction of new bioactive compounds for possible inhibition. The De novo molecular design of drugs method in the incremental construction of a ligant model within a receptor model was used, producing new structures with the help of artificial intelligence. The research algorithm and the scoring function responsible for predicting orientation and affinity in the molecular target at the time of coupling showed, as a result of the simulation, the compound with the highest bioaffinity value, Hit 998, with the energy of -17.62 kcal/mol, and synthetic viability close to 50%. While hit 1103 presented better synthetic viability (80%), its affinity energy of -10.28 kcal/mol. Both were compared with the reference linker N3, with a binding affinity of -7.5 kcal/mol. ADMET tests demonstrated that simulated compounds have a low risk of metabolic activation and do not exert effective distribution in the CNS, suggesting a pharmacokinetic mechanism based on local action, even with high topological polarity, which resulted in low oral bioavailability. In conclusion, MMGBSA, H-bonds, RMSD, SASA, and RMSF values were also obtained through molecular dynamics to verify the stability of the receptor-ligant complex within the active protein site to seek new therapeutic propositions in the fight against the pandemic.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Mestrado Acadêmico em Sociobiodiversidades e Tecnologias Sustentáveis - MASTS, Instituto de Engenharias e Desenvolvimento Sustentável, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Sadrack Queque Cabongo
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Bernardino Joaquim Caluaco
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Regilany Paulo Colares
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | | | | | - Pedro de Lima-Neto
- Department of Analytical Chemistry and Physical Chemistry, Science Center, Federal University of Ceara, Fortaleza, CE, Brazil
| | - Emmanuel Silva Marinho
- Grupo de química Teorica e Eletroquimica-GQTE, Universidade Estadual do Ceará, Limoeiro do Norte, CE, Brazil
| |
Collapse
|
20
|
da Fonseca AM, Luthierre Gama Cavalcante A, Mendes AMDS, da Silva FDFC, Ferreira DCL, Ribeiro PRV, Dos Santos JCS, Dos Santos HS, Gaieta EM, Marinho GS, Colares RP, Marinho ES. Phytochemical study of Lantana camara flowers, ecotoxicity, antioxidant, in vitro and in silico acetylcholinesterase: molecular docking, MD, and MM/GBSA calculations. J Biomol Struct Dyn 2023; 41:9282-9296. [PMID: 36326114 DOI: 10.1080/07391102.2022.2141883] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Accepted: 10/25/2022] [Indexed: 11/05/2022]
Abstract
Lantana camara L. (Verbenaceae), commonly called lead cambará, has often been used in folk medicine as antiseptic, antispasmodic, against hemorrhages, flu, colds, and diarrheic. This plant is considered a weed and an ornamental and medicinal plant and is an essential source of natural organic compounds, mainly flavonoids. This work aims to investigate the chemical composition and evaluate the biological properties such as antioxidant and acetylcholinesterase of the constituents from L. camara flowers. In addition, the computational simulation was carried out with the constituents identified. The results showed that methanolic extract of the flowers of L. camara presents toxicity, antioxidant activity with 97.8% inhibition percentage in the concentration of 0.25 mg mL-1 against the DPPH radical, and acetylcholinesterase activity. The phytochemical study of extract from L. camara flowers resulted in LC-MS identification of 18 polyphenolic compounds, such as phenolic acid derivatives, phenylethanoid glycosides, and flavonoids. In the in silico study, flavonoid isoverbascoside showed affinity energy of -9.9 kcal.mol-1 with the AChE enzyme. Their phytochemical content, mainly the presence of flavonoids and phenolic compounds in L. camara extracts, may be related to the antioxidant and anticholinesterase potential observed.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusophony, Redenção, Brazil
| | | | - Antônia Mayara Dos Santos Mendes
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusophony, Redenção, Brazil
| | | | - Débora Cristina Lima Ferreira
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusofonia, Redenção, Brazil
| | | | - José Cleiton Sousa Dos Santos
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusofonia, Redenção, Brazil
| | | | - Eduardo Menezes Gaieta
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, Brazil
| | - Gabrielle Silva Marinho
- Theoretical Chemistry and Electrochemistry Group, State University of Ceará, Limoeiro do Norte, Brazil
| | - Regilany Paulo Colares
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, Brazil
| | - Emmanuel Silva Marinho
- Theoretical Chemistry and Electrochemistry Group, State University of Ceará, Limoeiro do Norte, Brazil
| |
Collapse
|
21
|
Pant P, Leese F. Probing the Nucleic Acid Flexibility to Disarm the Viral Counter-Defense Machinery: Design and Characterization of Potent p19 Inhibitors. J Phys Chem B 2023; 127:8842-8851. [PMID: 37797202 DOI: 10.1021/acs.jpcb.3c04788] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
Plant viruses are highly destructive and significant contributors to several global pandemics and epidemics in plants. A viral disease outbreak in plants can cause a scarcity of food supply and is a severe concern to humanity. The siRNA (small interfering RNA)-mediated RNA-induced silencing complex (RISC) formation is a primary defense mechanism in plants against viruses, where the RISC binds and degrades viral mRNAs. As a counter-defense, many viruses encode RNA-silencing suppressor proteins (e.g., the p19 protein from the Tombusviridae family) for viral proliferation in plants. The functional form of p19 (homodimer) binds to plant siRNA with high affinities, thereby interrupting the RISC formation and thus preventing the viral mRNA silencing in plants. By altering the RISC formation, the p19 protein helps the virus invasion in the plant and ultimately stunts host growth. In this study, we designed several modified siRNA-based molecules for p19 inhibition. The viral p19 protein is known to interact predominantly through H-bonds with 2'-OH and phosphates of the plant siRNA. We utilized this information and in silico-designed flexible substituents of siRNA, where we removed the C2'-C3' bond in each nucleotide unit. We performed all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicates each) for control/modified siRNA─p19 complexes (8 in total) followed by energetic estimations. Strikingly, in a few modified complexes, the siRNA not only retained the double-helical structural integrity but also displayed remarkably enhanced p19 binding compared to the control siRNA; hence, we consider it important to perform biological and chemical in vitro and in vivo studies on proposed flexible nucleic acids as p19 inhibitors for crop protection.
Collapse
Affiliation(s)
- Pradeep Pant
- Faculty of Biology, University of Duisburg Essen, Essen 45141, Germany
| | - Florian Leese
- Faculty of Biology, University of Duisburg Essen, Essen 45141, Germany
| |
Collapse
|
22
|
Rocha Aguiar G, Leda Gomes de Lemos T, Braz-Filho R, Marques da Fonseca A, Silva Marinho E, Vasconcelos Ribeiro PR, Marques Canuto K, Queiroz Monte FJ. Synthesis and in silico study of chenodeoxycholic acid and its analogues as an alternative inhibitor of spike glycoprotein of SARS-CoV-2. J Biomol Struct Dyn 2023; 41:8334-8348. [PMID: 36218138 DOI: 10.1080/07391102.2022.2133010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Accepted: 09/30/2022] [Indexed: 10/17/2022]
Abstract
COVID-19, caused by SARS-CoV-2, is a viral infection that has generated one of the most significant health problems in the world. Spike glycoprotein is a crucial enzyme in viral replication and transcription mediation. There are reports in the literature on using bile acid in the fight against this virus through in vitro tests. This work presents the synthesis of nine chenodeoxycholic acid derivatives (1-9), which were prepared by oxidation, acetylation, formylation, and esterification reactions, and the analogs 6-9 have not yet been reported in the literature and the possibility of conducting an in silico study of bile acid derivatives as a therapeutic alternative to combat the virus using glycoprotein as a macromolecular target. As a result, five compounds (1, 6-9) possessed favorable competitive interactions with the lowest energies compared to the native ligand (BLA), and the highlighted compound 9 got the best scores. At the same time, analog 1 presented the best ADME filter result. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. Physical and spectroscopic data have fully characterized all the compounds.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Gisele Rocha Aguiar
- Departamento de Química Orgânica, Universidade Federal do Ceará, Fortaleza-CE, Brazil
| | | | - Raimundo Braz-Filho
- Laboratório de Ciências Químicas, Universidade Estadual do Norte Fluminense Darcy Ribeiro, Rio de Janeiro-RJ, Brazil
| | - Aluísio Marques da Fonseca
- Instituto de Ciências Exatas e Naturais, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Redenção-CE, Brazil
| | - Emmanuel Silva Marinho
- Faculdade de Filosofia Dom Aureliano Matos, Universidade Estadual do Ceará, Limoeiro do Norte-CE, Brazil
| | | | - Kirley Marques Canuto
- Laboratório multiusuário de Química de Produtos Naturais, Embrapa Agroindústria Tropical, Fortaleza-CE, Brazil
| | | |
Collapse
|
23
|
Dean E, Dominique A, Palillero A, Tran A, Paradis N, Wu C. Probing the Activation Mechanisms of Agonist DPI-287 to Delta-Opioid Receptor and Novel Agonists Using Ensemble-Based Virtual Screening with Molecular Dynamics Simulations. ACS OMEGA 2023; 8:32404-32423. [PMID: 37720760 PMCID: PMC10500586 DOI: 10.1021/acsomega.3c01918] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/21/2023] [Accepted: 08/16/2023] [Indexed: 09/19/2023]
Abstract
Pain drugs targeting mu-opioid receptors face major addiction problems that have caused an epidemic. The delta-opioid receptor (DOR) has shown to not cause addictive effects when bound to an agonist. While the active conformation of the DOR in complex with agonist DPI-287 has been recently solved, there are still no FDA-approved agonists targeting it, providing the opportunity for structure-based virtual screening. In this study, the conformational plasticity of the DOR was probed using molecular dynamics (MD) simulations, identifying two representative conformations from clustering analysis. The two MD conformations as well as the crystal conformation of DOR were used to screen novel compounds from the ZINC database (17 million compounds), in which 69 drugs were picked as potential compounds based on their docking scores. Notably, 37 out of the 69 compounds were obtained from the simulated conformations. The binding stability of the 69 compounds was further investigated using MD simulations. Based on the MM-GBSA binding energy and the predicted drug properties, eight compounds were chosen as the most favorable, six of which were from the simulated conformations. Using a dynamic network model, the communication between the crystal agonist and the top eight molecules with the receptor was analyzed to confirm if these novel compounds share a similar activation mechanism to the crystal ligand. Encouragingly, docking of these eight compounds to the other two opioid receptors (kappa and mu) suggests their good selectivity toward DOR.
Collapse
Affiliation(s)
- Emily Dean
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| | - AnneMarie Dominique
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| | - Americus Palillero
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| | - Annie Tran
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| | - Nicholas Paradis
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| | - Chun Wu
- Department of Molecular &
Cellular Biosciences, College of Science and Mathematics, Rowan University, Glassboro, New Jersey 08028, United States
| |
Collapse
|
24
|
Williams AH, Zhan CG. Staying Ahead of the Game: How SARS-CoV-2 has Accelerated the Application of Machine Learning in Pandemic Management. BioDrugs 2023; 37:649-674. [PMID: 37464099 DOI: 10.1007/s40259-023-00611-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/28/2023] [Indexed: 07/20/2023]
Abstract
In recent years, machine learning (ML) techniques have garnered considerable interest for their potential use in accelerating the rate of drug discovery. With the emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, the utilization of ML has become even more crucial in the search for effective antiviral medications. The pandemic has presented the scientific community with a unique challenge, and the rapid identification of potential treatments has become an urgent priority. Researchers have been able to accelerate the process of identifying drug candidates, repurposing existing drugs, and designing new compounds with desirable properties using machine learning in drug discovery. To train predictive models, ML techniques in drug discovery rely on the analysis of large datasets, including both experimental and clinical data. These models can be used to predict the biological activities, potential side effects, and interactions with specific target proteins of drug candidates. This strategy has proven to be an effective method for identifying potential coronavirus disease 2019 (COVID-19) and other disease treatments. This paper offers a thorough analysis of the various ML techniques implemented to combat COVID-19, including supervised and unsupervised learning, deep learning, and natural language processing. The paper discusses the impact of these techniques on pandemic drug development, including the identification of potential treatments, the understanding of the disease mechanism, and the creation of effective and safe therapeutics. The lessons learned can be applied to future outbreaks and drug discovery initiatives.
Collapse
Affiliation(s)
- Alexander H Williams
- Molecular Modeling and Biopharmaceutical Center, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA
- GSK Upper Providence, 1250 S. Collegeville Road, Collegeville, PA, 19426, USA
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY, 40536, USA.
| |
Collapse
|
25
|
Lima JPO, da Fonseca AM, Marinho GS, da Rocha MN, Marinho EM, dos Santos HS, Freire RM, Marinho ES, de Lima-Neto P, Fechine PBA. De novo design of bioactive phenol and chromone derivatives for inhibitors of Spike glycoprotein of SARS-CoV-2 in silico. 3 Biotech 2023; 13:301. [PMID: 37588795 PMCID: PMC10425314 DOI: 10.1007/s13205-023-03695-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2023] [Accepted: 06/29/2023] [Indexed: 08/18/2023] Open
Abstract
This work presents the synthesis of 12 phenol and chromone derivatives, prepared by the analogs, and the possibility of conducting an in silico study of its derivatives as a therapeutic alternative to combat the SARS-CoV-2, pathogen responsible for COVID-19 pandemic, using its S-glycoprotein as a macromolecular target. After the initial screening for the ranking of the products, it was chosen which structure presented the best energy bond with the target. As a result, derivative 4 was submitted to a molecular growth study using artificial intelligence, where 8436 initial structures were obtained that passed through the interaction filters and similarity to the active glycoprotein pocket through the MolAICal computational package. Thus, 557 Hits with active configuration were generated, which is very promising compared to the BLA reference link for inhibiting the biological target. Molecular dynamics also simulated these compounds to verify their stability within the active protein site to seek new therapeutic propositions to fight against the pandemic. The Hit 48 and 250 are the most active compounds against SARS-CoV-2. In summary, the results show that the Hit 250 would be more active than the natural compound, which could be further developed for further testing against SARS-CoV-2. The study employs the de novo approach to design new drugs, combining artificial intelligence and molecular dynamics simulations to create efficient molecular structures. This research aims to contribute to the development of effective therapeutic strategies against the pandemic.
Collapse
Affiliation(s)
- Joan Petrus Oliveira Lima
- Advanced Materials Chemistry Group (GQMat)-Department of Analytical Chemistry and Physical Chemistry, Federal University of Ceará, Campus Pici, Fortaleza, Ceará 60455-970 Brazil
| | - Aluísio Marques da Fonseca
- Mestrado Acadêmico em Sociobiodiversidades e Tecnologias Sustentáveis-MASTS, Instituto de Engenharias e Desenvolvimento Sustentável, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE 62785-000 Brazil
| | - Gabrielle Silva Marinho
- Faculdade de Filosofia Dom Aureliano Matos-FAFIDAM, Universidade Estadual do Ceará, Centro, Limoeiro do Norte, CE 62930-000 Brazil
| | - Matheus Nunes da Rocha
- Faculdade de Filosofia Dom Aureliano Matos-FAFIDAM, Universidade Estadual do Ceará, Centro, Limoeiro do Norte, CE 62930-000 Brazil
| | - Emanuelle Machado Marinho
- Advanced Materials Chemistry Group (GQMat)-Department of Analytical Chemistry and Physical Chemistry, Federal University of Ceará, Campus Pici, Fortaleza, Ceará 60455-970 Brazil
| | | | | | - Emmanuel Silva Marinho
- Faculdade de Filosofia Dom Aureliano Matos-FAFIDAM, Universidade Estadual do Ceará, Centro, Limoeiro do Norte, CE 62930-000 Brazil
| | - Pedro de Lima-Neto
- Advanced Materials Chemistry Group (GQMat)-Department of Analytical Chemistry and Physical Chemistry, Federal University of Ceará, Campus Pici, Fortaleza, Ceará 60455-970 Brazil
| | - Pierre Basílio Almeida Fechine
- Advanced Materials Chemistry Group (GQMat)-Department of Analytical Chemistry and Physical Chemistry, Federal University of Ceará, Campus Pici, Fortaleza, Ceará 60455-970 Brazil
| |
Collapse
|
26
|
Hong X, Tong X, Xie J, Liu P, Liu X, Song Q, Liu S, Liu S. An updated dataset and a structure-based prediction model for protein-RNA binding affinity. Proteins 2023; 91:1245-1253. [PMID: 37186412 DOI: 10.1002/prot.26503] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Revised: 03/08/2023] [Accepted: 04/12/2023] [Indexed: 05/17/2023]
Abstract
Understanding the process of protein-RNA interaction is essential for structural biology. The thermodynamic process is an important part to uncover the protein-RNA interaction mechanism. The regulatory networks between protein and RNA in organisms are dominated by the binding or dissociation in the cells. Therefore, determining the binding affinity for protein-RNA complexes can help us to understand the regulation mechanism of protein-RNA interaction. Since it is time-consuming and labor-intensive to determine the binding affinity for protein-RNA complexes by experimental methods, it is necessary and urgent to develop computational methods to predict that. To develop a binding affinity prediction model, first we update the dataset of protein-RNA binding affinity benchmark (PRBAB), which includes 145 complexes now. Second, we extract the structural features based on complex structure, and then we analyze and select the representative structural features to train the regression model. Third, we random select the subset from the PRBAB2.0 to fit the protein-RNA binding affinity determined by experiment. In the end, we tested our model on the nonredundant PDBbind dataset, and the results showed that Pearson correlation coefficient r = .57 and RMSE = 2.51 kcal/mol. The Pearson correlation coefficient achieves 0.7 while removing 5 complex structures with modified residues/nucleotides and metal ions. While testing on ProNAB, the results showed that 71.60% of the prediction achieves Pearson correlation coefficient r = .61 and RMSE = 1.56 kcal/mol with experiment values.
Collapse
Affiliation(s)
- Xu Hong
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Xiaoxue Tong
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Juan Xie
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Pinyu Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Xudong Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Qi Song
- Key Laboratory of Fermentation Engineering (Ministry of Education), Hubei University of Technology, Wuhan, China
| | - Sen Liu
- Key Laboratory of Fermentation Engineering (Ministry of Education), Hubei University of Technology, Wuhan, China
| | - Shiyong Liu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
27
|
Tang R, Wang Z, Xiang S, Wang L, Yu Y, Wang Q, Deng Q, Hou T, Sun H. Uncovering the Kinetic Characteristics and Degradation Preference of PROTAC Systems with Advanced Theoretical Analyses. JACS AU 2023; 3:1775-1789. [PMID: 37388700 PMCID: PMC10301679 DOI: 10.1021/jacsau.3c00195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 05/16/2023] [Accepted: 05/16/2023] [Indexed: 07/01/2023]
Abstract
Proteolysis-targeting chimeras (PROTACs), which can selectively induce the degradation of target proteins, represent an attractive technology in drug discovery. A large number of PROTACs have been reported, but due to the complicated structural and kinetic characteristics of the target-PROTAC-E3 ligase ternary interaction process, the rational design of PROTACs is still quite challenging. Here, we characterized and analyzed the kinetic mechanism of MZ1, a PROTAC that targets the bromodomain (BD) of the bromodomain and extra terminal (BET) protein (Brd2, Brd3, or Brd4) and von Hippel-Lindau E3 ligase (VHL), from the kinetic and thermodynamic perspectives of view by using enhanced sampling simulations and free energy calculations. The simulations yielded satisfactory predictions on the relative residence time and standard binding free energy (rp > 0.9) for MZ1 in different BrdBD-MZ1-VHL ternary complexes. Interestingly, the simulation of the PROTAC ternary complex disintegration illustrates that MZ1 tends to remain on the surface of VHL with the BD proteins dissociating alone without a specific dissociation direction, indicating that the PROTAC prefers more to bind with E3 ligase at the first step in the formation of the target-PROTAC-E3 ligase ternary complex. Further exploration of the degradation difference of MZ1 in different Brd systems shows that the PROTAC with higher degradation efficiency tends to leave more lysine exposed on the target protein, which is guaranteed by the stability (binding affinity) and durability (residence time) of the target-PROTAC-E3 ligase ternary complex. It is quite possible that the underlying binding characteristics of the BrdBD-MZ1-VHL systems revealed by this study may be shared by different PROTAC systems as a general rule, which may accelerate rational PROTAC design with higher degradation efficiency.
Collapse
Affiliation(s)
- Rongfan Tang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Zhe Wang
- Innovation
Institute for Artificial Intelligence in Medicine of Zhejiang University,
College of Pharmaceutical Sciences, Zhejiang
University, Hangzhou 310058, Zhejiang, P. R. China
| | - Sutong Xiang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Lingling Wang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Yang Yu
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Qinghua Wang
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Qirui Deng
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| | - Tingjun Hou
- Innovation
Institute for Artificial Intelligence in Medicine of Zhejiang University,
College of Pharmaceutical Sciences, Zhejiang
University, Hangzhou 310058, Zhejiang, P. R. China
| | - Huiyong Sun
- Department
of Medicinal Chemistry, China Pharmaceutical
University, Nanjing 210009, Jiangsu, P. R. China
| |
Collapse
|
28
|
Akkus E, Tayfuroglu O, Yildiz M, Kocak A. Revisiting MMPBSA by Adoption of MC-Based Surface Area/Volume, ANI-ML Potentials, and Two-Valued Interior Dielectric Constant. J Phys Chem B 2023; 127:4415-4429. [PMID: 37171911 DOI: 10.1021/acs.jpcb.3c00834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Here, we report the accuracy improvements of molecular mechanics Poisson-Boltzmann surface area (MMPBSA) calculations by adoption of ANI-ML potentials in replacement of MM terms, the use of solvent-accessible surface area (SASA) and volume (SAV) values from the Monte Carlo sampling of the probe, and introducing two different interior dielectric constants for electrostatic interactions of protein-ligand (P-L) and polar solvation term in the MMPBSA calculations. Our results show that the Pearson correlation coefficients of MMPBSA-calculated values with respect to experimental binding free energies can be drastically improved from 0.48 to 0.90 by adoption of ANI-ML potentials in replacement of MM energy terms in the equation, referred to as ANI-PBSA. Moreover, we show that the SASA/SAV-combined equation in the scaled particle theory (SPT) can be a better choice to model nonpolar solvation term, reaching nearly the same accuracy by ANI-PBSA calculations. Finally, we introduce two different values of interior dielectric constants, which could be an alternative strategy between the single and variable constant definitions.
Collapse
Affiliation(s)
- Ebru Akkus
- Department of Chemistry, Gebze Technical University, 41400 Kocaeli, Turkey
| | - Omer Tayfuroglu
- Department of Chemistry, Gebze Technical University, 41400 Kocaeli, Turkey
| | - Muslum Yildiz
- Department of Molecular Biology and Genetics, Gebze Technical University, 41400 Kocaeli, Turkey
| | - Abdulkadir Kocak
- Department of Chemistry, Gebze Technical University, 41400 Kocaeli, Turkey
| |
Collapse
|
29
|
Marques da Fonseca A, Freire da Silva A, Barbosa da Silva FL, Caluaco BJ, Gaieta EM, Nunes da Rocha M, Colares RP, Sobczak JF, Marinho GS, Dos Santos HS, Marinho ES. Isolation, characterization and in silico study of propenamide alkaloids from Hymenoepmecis bicolor poison against active μ-opioid receptor. J Biomol Struct Dyn 2023; 41:14621-14637. [PMID: 36815273 DOI: 10.1080/07391102.2023.2183043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 02/16/2023] [Indexed: 02/24/2023]
Abstract
Some insects produce venoms to defend against predators and directly interact with opioid receptors. In the present study, it was identified two alkaloids in the wasp venom species Hymenoepimecis bicolor. It was demonstrated that these could act as potential inhibitors of opioid receptors through their robust affinity to the receptors. The interaction profile was given to opioid receptors (μOR), with 60% of targets similar to alkaloid 1, with 0.25 probability, and 46.7% of targets similar to alkaloid 2, with a probability 0.17 of affinity as a target, which is considered signaling macromolecules and can mediate the most potent analgesic and addictive properties of opiate alkaloids. Notably, both alkaloids showed -7.6 kcal/mol affinity to the morphine agonies through six residues, Gly124, Asp147, Trp293, Ile296, Ile322, and Tyr326. These observations suggest further research on opioid receptors using in vitro studies of possible therapeutic applications.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, CE, Brazil
| | - Ananias Freire da Silva
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusofonia, Redenção, CE, Brazil
| | - Francisco Lennon Barbosa da Silva
- Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusofonia, Redenção, CE, Brazil
| | - Bernardino Joaquim Caluaco
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, CE, Brazil
| | - Eduardo Menezes Gaieta
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, CE, Brazil
| | - Matheus Nunes da Rocha
- Faculty of Philosophy Dom Aureliano Matos, State University of Ceará, Limoeiro do Norte, CE, Brazil
| | - Regilany Paulo Colares
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, CE, Brazil
| | - Jober Fernando Sobczak
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Redenção, CE, Brazil
| | - Gabrielle Silva Marinho
- Faculty of Philosophy Dom Aureliano Matos, State University of Ceará, Limoeiro do Norte, CE, Brazil
| | | | - Emmanuel Silva Marinho
- Faculty of Philosophy Dom Aureliano Matos, State University of Ceará, Limoeiro do Norte, CE, Brazil
| |
Collapse
|
30
|
Sun Y, He X, Hou T, Cai L, Man VH, Wang J. Development and test of highly accurate endpoint free energy methods. 1: Evaluation of ABCG2 charge model on solvation free energy prediction and optimization of atom radii suitable for more accurate solvation free energy prediction by the PBSA method. J Comput Chem 2023; 44:1334-1346. [PMID: 36807356 DOI: 10.1002/jcc.27089] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2022] [Revised: 12/22/2022] [Accepted: 01/29/2023] [Indexed: 02/23/2023]
Abstract
Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson-Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1-BCC-GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1-BCC (Austin Model 1-bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types, on, oi, hn1, hn2, hn3, were introduced to further improve the fitting performance. Then, we tuned the parameters in the nonpolar SFE model using the experimental SFE data and the PB calculation results. By adopting the new radius parameters and new nonpolar SFE model, the root mean square error (RMSE) of the SFE calculation for the 544 molecules decreased from 2.38 to 1.05 kcal/mol. Finally, the new radius parameters were applied in the prediction of protein-ligand binding free energies using the MM-PBSA method. For the eight systems tested, we could observe higher correlation between the experiment data and calculation results and smaller prediction errors for the absolute binding free energies, demonstrating that our new radius parameters can improve the free energy calculation using the MM-PBSA method.
Collapse
Affiliation(s)
- Yuchen Sun
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Tingjun Hou
- College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, China
| | - Lianjin Cai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Viet Hoag Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| |
Collapse
|
31
|
da Fonseca AM, Soares NB, Colares RP, Macedo de Oliveira M, Santos Oliveira L, Marinho GS, Raya Paula de Lima M, da Rocha MN, Dos Santos HS, Marinho ES. Naphthoquinones biflorin and bis-biflorin ( Capraria biflora) as possible inhibitors of the fungus Candida auris polymerase: molecular docking, molecular dynamics, MM/GBSA calculations and in silico drug-likeness study. J Biomol Struct Dyn 2023; 41:11564-11577. [PMID: 36597918 DOI: 10.1080/07391102.2022.2163702] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Accepted: 12/24/2022] [Indexed: 01/05/2023]
Abstract
A new worldwide concern has emerged with the recent emergence of infections caused by Candida auris. This reflects its comparative ease of transmission, substantial mortality, and the increasing level of resistance seen in the three major classes of antifungal drugs. Efforts to create a better design for structure-based drugs that described numerous modifications and the search for secondary metabolic structures derived from plant species are likely to reduce the virulence of several fungal pathogens. In this context, the present work aimed to evaluate in silico two naphthoquinones isolated from the roots of Capraria biflora, biflorin, and its dimmer, bis-biflorin, as potential inhibitors of Candida auris polymerase. Based on the simulation performed with the two naphthoquinones, biflorin and bis-biflorin, it can be stated that bis-biflorin showed the best interactions with Candida auris polymerase. Still, biflorin also demonstrated favorable coupling energy. Predictive pharmacokinetic assays suggest that biflorin has high oral bioavailability and more excellent metabolic stability compared to the bis-biflorin analogue. constituting a promising pharmacological tool.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Mestrado Acadêmico em Sociobiodiversidades e Tecnologias Sustentáveis - MASTS, Instituto de Engenharias e Desenvolvimento Sustentável, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Neidelenio Baltazar Soares
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | - Regilany Paulo Colares
- Instituto de Ciências Exatas e da Natureza, Universidade da Integração Internacional da Lusofonia Afro-Brasileira, Acarape, CE, Brazil
| | | | | | - Gabrielle Silva Marinho
- Grupo de Química Teórica e Eletroquímica - GQTE, Universidade Estadual de Ceará, Limoiro do Norte, CE, Brazil
| | - Mira Raya Paula de Lima
- Instituto Federal de Educação Ciência e Tecnologia do Ceará - Campus Juazeiro do Norte, Juazeiro do Norte, CE, Brazil
| | - Matheus Nunes da Rocha
- Grupo de Química Teórica e Eletroquímica - GQTE, Universidade Estadual de Ceará, Limoiro do Norte, CE, Brazil
| | | | - Emmanuel Silva Marinho
- Grupo de Química Teórica e Eletroquímica - GQTE, Universidade Estadual de Ceará, Limoiro do Norte, CE, Brazil
| |
Collapse
|
32
|
Durgam L, Guruprasad L. Molecular mechanism of ATP and RNA binding to Zika virus NS3 helicase and identification of repurposed drugs using molecular dynamics simulations. J Biomol Struct Dyn 2022; 40:12642-12659. [PMID: 34516356 DOI: 10.1080/07391102.2021.1973909] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Congenital Zika virus syndrome has caused a public health emergency of international concern. So far, there are no drugs available to prevent or treat the infection caused by Zika virus. The Zika virus NS3 helicase is a potential protein target for drug discovery due to its vital role in viral genome replication. NS3 helicase unwinds the viral RNA to enable the reproduction of the viral genome by the NS5 protein. NS3 helicase has two crucial binding sites; the ATP binding site and the RNA binding site. We used molecular docking and molecular dynamics (MD) simulations to study the structural behavior of Zika virus NS3 helicase in its apo form and in the presence of ATP, single-stranded RNA, and both ATP-RNA to understand their potential implications in NS3 helicase activity. Further, we have carried out virtual screening of FDA approved drugs, followed by molecular docking to identify the ATP-competitive hit molecules as probable Zika virus NS3 helicase inhibitors. The MD simulations trajectories were analyzed using normal mode analysis and principal component analysis that reveals fluctuations in the R-loop. These findings aid in understanding the molecular mechanisms of the simultaneous binding of ATP and RNA, and guide the design and discovery of new inhibitors of the Zika virus NS3 helicase as a promising drug target to treat the Zika virus infection. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Laxman Durgam
- School of Chemistry, University of Hyderabad, Hyderabad, India
| | | |
Collapse
|
33
|
Pavan M, Bassani D, Sturlese M, Moro S. Investigating RNA-protein recognition mechanisms through supervised molecular dynamics (SuMD) simulations. NAR Genom Bioinform 2022; 4:lqac088. [PMID: 36458023 PMCID: PMC9706429 DOI: 10.1093/nargab/lqac088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Revised: 10/20/2022] [Accepted: 11/09/2022] [Indexed: 12/03/2022] Open
Abstract
Ribonucleic acid (RNA) plays a key regulatory role within the cell, cooperating with proteins to control the genome expression and several biological processes. Due to its characteristic structural features, this polymer can mold itself into different three-dimensional structures able to recognize target biomolecules with high affinity and specificity, thereby attracting the interest of drug developers and medicinal chemists. One successful example of the exploitation of RNA's structural and functional peculiarities is represented by aptamers, a class of therapeutic and diagnostic tools that can recognize and tightly bind several pharmaceutically relevant targets, ranging from small molecules to proteins, making use of the available structural and conformational freedom to maximize the complementarity with their interacting counterparts. In this scientific work, we present the first application of Supervised Molecular Dynamics (SuMD), an enhanced sampling Molecular Dynamics-based method for the study of receptor-ligand association processes in the nanoseconds timescale, to the study of recognition pathways between RNA aptamers and proteins, elucidating the main advantages and limitations of the technique while discussing its possible role in the rational design of RNA-based therapeutics.
Collapse
Affiliation(s)
- Matteo Pavan
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences University of Padova, via Marzolo 5, 35131 Padova, Italy
| | - Davide Bassani
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences University of Padova, via Marzolo 5, 35131 Padova, Italy
| | - Mattia Sturlese
- Molecular Modeling Section (MMS), Department of Pharmaceutical and Pharmacological Sciences University of Padova, via Marzolo 5, 35131 Padova, Italy
| | - Stefano Moro
- To whom correspondence should be addressed. Tel: +39 0498275704; Fax: +39 0498275366;
| |
Collapse
|
34
|
Frota LS, da Rocha MN, Bezerra LL, da Fonseca AM, Marinho ES, de Morais SM. HIF1 inhibition of the biflavonoids against pancreas cancer: drug-likeness, bioavailability, ADMET, PASS, molecular docking, molecular dynamics, and MM/GBSA calculations. J Biomol Struct Dyn 2022:1-12. [PMID: 36002285 DOI: 10.1080/07391102.2022.2112619] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Abstract
Pancreatic cancer is an aggressive disease with a high death rate and is difficult to treat. This disease, in the most cases, is asymptomatic until it progresses to an advanced stage. Therefore, the search for bioactive molecules is urgent to combat pancreatic cancer. Then, this work analyzed the interaction potential of agathisflavone and amentoflavone molecules against the HIF1 target using the ADMET, molecular docking, and molecular dynamics simulations. More recent drug-likeness filters that combine physicochemical and physiological parameters have shown that high polar surface area (TPSA > 75 Å2) drives biflavonoids out of the toxic drug space of Pfizer dataset. Regarding the pharmacokinetic descriptors, it was possible to notice that Amentoflavone showed a better order of passive cell permeability (Papp = 8 × 10-6 cm/s) and better metabolic stability, biotransformed by aromatic hydroxylation reactions by the CYP3A4 isoenzyme on the human liver, that favor its hepatic clearance. The molecular docking and molecular dynamics simulations indicated the high interaction potential and stability between the ligands analyzed (highlighted the amentoflavone molecule), respectively. The MM/GBSA calculations showed that the amentoflavone ligand registered the highest ΔG binding value of -32.6957 kcal/mol with the HIF1 target. Then, this molecule may be used as a potential inhibitor of pancreatic cancer. In this perspective, the present work represents an initial step in the virtual bioprospecting a pharmacological tool for treating of pancreatic cancer.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Lucas Soares Frota
- Graduate Program in Biotechnology, Northeast Biotechnology Network, Faculty of Veterinary Medicine, State University of Ceará, Campus do Itaperi, Fortaleza , Ceara, Brazil
| | - Matheus Nunes da Rocha
- Group of Theoretical Chemistry and Electrochemistry, Faculty of Philosophy Dom Aureliano Matos, State University of Ceara, Limoeiro do Norte, Ceará, Brazil
| | - Lucas Lima Bezerra
- Group of Theoretical Chemistry and Electrochemistry, Faculty of Philosophy Dom Aureliano Matos, State University of Ceara, Limoeiro do Norte, Ceará, Brazil
| | - Aluísio Marques da Fonseca
- Academic Master in Sociobiodiversity and Sustainable Technologies, Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusophony, Auroras Campus, Redenção, CE, Brazil
| | - Emmanuel Silva Marinho
- Group of Theoretical Chemistry and Electrochemistry, Faculty of Philosophy Dom Aureliano Matos, State University of Ceara, Limoeiro do Norte, Ceará, Brazil
| | - Selene Maia de Morais
- Chemistry Course, Science and Technology Center, State University of Ceará, Itaperi Campus, 60714-903, Fortaleza, Ceará, Brazil
| |
Collapse
|
35
|
Yu Y, Wang Z, Wang L, Tian S, Hou T, Sun H. Predicting the mutation effects of protein–ligand interactions via end-point binding free energy calculations: strategies and analyses. J Cheminform 2022; 14:56. [PMID: 35987841 PMCID: PMC9392442 DOI: 10.1186/s13321-022-00639-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2022] [Accepted: 08/08/2022] [Indexed: 12/04/2022] Open
Abstract
Protein mutations occur frequently in biological systems, which may impact, for example, the binding of drugs to their targets through impairing the critical H-bonds, changing the hydrophobic interactions, etc. Thus, accurately predicting the effects of mutations on biological systems is of great interests to various fields. Unfortunately, it is still unavailable to conduct large-scale wet-lab mutation experiments because of the unaffordable experimental time and financial costs. Alternatively, in silico computation can serve as a pioneer to guide the experiments. In fact, numerous pioneering works have been conducted from computationally cheaper machine-learning (ML) methods to the more expensive alchemical methods with the purpose to accurately predict the mutation effects. However, these methods usually either cannot result in a physically understandable model (ML-based methods) or work with huge computational resources (alchemical methods). Thus, compromised methods with good physical characteristics and high computational efficiency are expected. Therefore, here, we conducted a comprehensive investigation on the mutation issues of biological systems with the famous end-point binding free energy calculation methods represented by MM/GBSA and MM/PBSA. Different computational strategies considering different length of MD simulations, different value of dielectric constants and whether to incorporate entropy effects to the predicted total binding affinities were investigated to provide a more accurate way for predicting the energetic change upon protein mutations. Overall, our result shows that a relatively long MD simulation (e.g. 100 ns) benefits the prediction accuracy for both MM/GBSA and MM/PBSA (with the best Pearson correlation coefficient between the predicted ∆∆G and the experimental data of ~ 0.44 for a challenging dataset). Further analyses shows that systems involving large perturbations (e.g. multiple mutations and large number of atoms change in the mutation site) are much easier to be accurately predicted since the algorithm works more sensitively to the large change of the systems. Besides, system-specific investigation reveals that conformational adjustment is needed to refine the micro-environment of the manually mutated systems and thus lead one to understand why longer MD simulation is necessary to improve the predicting result. The proposed strategy is expected to be applied in large-scale mutation effects investigation with interpretation.
Collapse
|
36
|
Zhai J, He X, Sun Y, Wan Z, Ji B, Liu S, Li S, Wang J. In silico binding affinity prediction for metabotropic glutamate receptors using both endpoint free energy methods and a machine learning-based scoring function. Phys Chem Chem Phys 2022; 24:18291-18305. [PMID: 35880533 PMCID: PMC9460939 DOI: 10.1039/d2cp01727j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Metabotropic glutamate receptors (mGluRs) play an important role in regulating glutamate signal pathways, which are involved in neuropathy and periphery homeostasis. mGluR4, which belongs to Group III mGluRs, is most widely distributed in the periphery among all the mGluRs. It has been proved that the regulation of this receptor is involved in diabetes, colorectal carcinoma and many other diseases. However, the application of structure-based drug design to identify small molecules to regulate the mGluR4 receptor is limited due to the absence of a resolved mGluR4 protein structure. In this work, we first built a homology model of mGluR4 based on a crystal structure of mGluR8, and then conducted hierarchical virtual screening (HVS) to identify possible active ligands for mGluR4. The HVS protocol consists of three hierarchical filters including Glide docking, molecular dynamic (MD) simulation and binding free energy calculation. We successfully prioritized active ligands of mGluR4 from a set of screening compounds using HVS. The predicted active ligands based on binding affinities can almost cover all the experiment-determined active ligands, with only one ligand missed. The correlation between the measured and predicted binding affinities is significantly improved for the MM-PB/GBSA-WSAS methods compared to the Glide docking method. More importantly, we have identified hotspots for ligand binding, and we found that SER157 and GLY158 tend to contribute to the selectivity of mGluR4 ligands, while ALA154 and ALA155 could account for the ligand selectivity to mGluR8. We also recognized other 5 key residues that are critical for ligand potency. The difference of the binding profiles between mGluR4 and mGluR8 can guide us to develop more potent and selective modulators. Moreover, we evaluated the performance of IPSF, a novel type of scoring function trained by a machine learning algorithm on residue-ligand interaction profiles, in guiding drug lead optimization. The cross-validation root-mean-square errors (RMSEs) are much smaller than those by the endpoint methods, and the correlation coefficients are comparable to the best endpoint methods for both mGluRs. Thus, machine learning-based IPSF can be applied to guide lead optimization, albeit the total number of actives/inactives are not big, a typical scenario in drug discovery projects.
Collapse
Affiliation(s)
- Jingchen Zhai
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Yuchen Sun
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Zhuoya Wan
- Department of Pharmaceutical Sciences and Center for Pharmacogenetics, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Beihong Ji
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Shuhan Liu
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| | - Song Li
- Department of Pharmaceutical Sciences and Center for Pharmacogenetics, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, PA 15261, USA.
| |
Collapse
|
37
|
Williams A, Zhan CG. Fast Prediction of Binding Affinities of SARS-CoV-2 Spike Protein and Its Mutants with Antibodies through Intermolecular Interaction Modeling-Based Machine Learning. J Phys Chem B 2022; 126:5194-5206. [PMID: 35817617 PMCID: PMC9301770 DOI: 10.1021/acs.jpcb.2c02123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 06/28/2022] [Indexed: 11/30/2022]
Abstract
Since the introduction of the novel SARS-CoV-2 virus (COVID-19) in late 2019, various new variants have appeared with mutations that confer resistance to the vaccines and monoclonal antibodies that were developed in response to the wild-type virus. As we continue through the pandemic, an accurate and efficient methodology is needed to help predict the effects certain mutations will have on both our currently produced therapeutics and those that are in development. Using published cryo-electron microscopy and X-ray crystallography structures of the spike receptor binding domain region with currently known antibodies, in the present study, we created and cross-validated an intermolecular interaction modeling-based multi-layer perceptron machine learning approach that can accurately predict the mutation-caused shifts in the binding affinity between the spike protein (wild-type or mutant) and various antibodies. This validated artificial intelligence (AI) model was used to predict the binding affinity (Kd) of reported SARS-CoV-2 antibodies with various variants of concern, including the most recently identified "Deltamicron" (or "Deltacron") variant. This AI model may be employed in the future to predict the Kd of developed novel antibody therapeutics to overcome the challenging antibody resistance issue and develop structural bases for the effects of both current and new mutants of the spike protein. In addition, the similar AI strategy and approach based on modeling of the intermolecular interactions may be useful in development of machine learning models predicting binding affinities for other protein-protein binding systems, including other antibodies binding with their antigens.
Collapse
Affiliation(s)
- Alexander
H. Williams
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| | - Chang-Guo Zhan
- Molecular
Modeling and Biopharmaceutical Center, University
of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
- Department
of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, Kentucky 40536, United States
| |
Collapse
|
38
|
Shin CH, Park SC, Park IG, Kim H, An B, Lee C, Kim SH, Lee J, Lee JM, Oh SJ. Cytosolic microRNA-inducible nuclear translocation of Cas9 protein for disease-specific genome modification. Nucleic Acids Res 2022; 50:5919-5933. [PMID: 35640600 PMCID: PMC9177975 DOI: 10.1093/nar/gkac431] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 05/03/2022] [Accepted: 05/10/2022] [Indexed: 11/13/2022] Open
Abstract
MicroRNA-dependent mRNA decay plays an important role in gene silencing by facilitating posttranscriptional and translational repression. Inspired by this intrinsic nature of microRNA-mediated mRNA cleavage, here, we describe a microRNA-targeting mRNA as a switch platform called mRNA bridge mimetics to regulate the translocation of proteins. We applied the mRNA bridge mimetics platform to Cas9 protein to confer it the ability to translocate into the nucleus via cleavage of the nuclear export signal. This system performed programmed gene editing in vitro and in vivo. Combinatorial treatment with cisplatin and miR-21-EZH2 axis-targeting CRISPR Self Check-In improved sensitivity to chemotherapeutic drugs in vivo. Using the endogenous microRNA-mediated mRNA decay mechanism, our platform is able to remodel a cell's natural biology to allow the entry of precise drugs into the nucleus, devoid of non-specific translocation. The mRNA bridge mimetics strategy is promising for applications in which the reaction must be controlled via intracellular stimuli and modulates Cas9 proteins to ensure safe genome modification in diseased conditions.
Collapse
Affiliation(s)
- Cheol-Hee Shin
- Center for Biomaterials, Biomedical Research Institute, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea
| | - Su Chan Park
- Graduate School of Medical Science & Engineering, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Il-Geun Park
- Graduate School of Medical Science & Engineering, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Hyerim Kim
- Program in Nanoscience and Technology, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, Republic of Korea
| | - Byoungha An
- Center for Biomaterials, Biomedical Research Institute, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea.,Division of Bio-Medical Science & Technology, Korea University of Science and Technology (UST), Seoul, Republic of Korea
| | - Choongil Lee
- Center for Biomaterials, Biomedical Research Institute, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea
| | - Sang-Heon Kim
- Center for Biomaterials, Biomedical Research Institute, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea.,Division of Bio-Medical Science & Technology, Korea University of Science and Technology (UST), Seoul, Republic of Korea
| | - Juyong Lee
- Department of Chemistry, College of Natural Science, Kangwon National University, Chuncheon 24341, Republic of Korea
| | - Ji Min Lee
- Graduate School of Medical Science & Engineering, Korea Advanced Institute of Science and Technology, 291 Daehak-ro, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Seung Ja Oh
- Center for Biomaterials, Biomedical Research Institute, Korea Institute of Science and Technology (KIST), Seoul 02792, Republic of Korea.,Division of Bio-Medical Science & Technology, Korea University of Science and Technology (UST), Seoul, Republic of Korea
| |
Collapse
|
39
|
da Fonseca AM, Soares NB, Meirú MIL, Colares RP, Neto MM, Sobrinho ACN, Dos Santos HS, Marinho ES. Combined study of docking and molecular dynamics against DNV-3 SN1 protein by bixinoids. J Biomol Struct Dyn 2022:1-11. [PMID: 35510585 DOI: 10.1080/07391102.2022.2070282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Dengue (DENV), Zica virus (ZIKV), and Chikungunya fever (CHIK) are tropical diseases that have caused a lot of problems in general worldwide. Transmitted by mosquitoes of the species Aedes aegypti and albopictus, they have not been completely eradicated in the country, and their proliferation has only increased in the Northeast region. Within the structure of the virus, it is possible to verify the presence of glycoprotein SN1, which is responsible for its replication. If this macromolecule is inhibited using a specific or complex linker, it can interrupt its replication activity. An alternative to this problem has been using structures derived from natural products that have pharmacological properties. A dynamic and molecular docking combined study used computational simulation in the four isomeric forms of bixin against the SN1 protein. The Z,E-bixin and E,E-bixin isomers, both with affinity energy -6.7 and -6.5 Kcal/mol, presented the best results. Thus, bixin and its isomers, found in annatto seeds, maybe an initial proposal in the search for prototype compounds to study to fight this lethal virus in the future.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Aluísio Marques da Fonseca
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Acarape, Brazil
| | - Neidelênio Baltazar Soares
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Acarape, Brazil
| | - Maria Imaculada Lourenço Meirú
- Academic Master in Sociobiodiversity and Sustainable Technologies, Institute of Engineering and Sustainable Development, University of International Integration of Afro-Brazilian Lusophony, Acarape, Brazil
| | - Regilany Paulo Colares
- Institute of Exact Sciences and Nature, University of International Integration of Afro-Brazilian Lusophony, Acarape, Brazil
| | | | | | | | - Emmanuel Silva Marinho
- Natural Sciences, Science and Technology Center, Ceará State University, Fortaleza, Brazil.,Theoretical and Electrochemical Chemistry Research Group/FAFIDAM, State University of Ceará, Limoeiro do Norte, Brazil
| |
Collapse
|
40
|
Zhou Y, Jiang Y, Chen SJ. RNA-ligand molecular docking: advances and challenges. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2022; 12:e1571. [PMID: 37293430 PMCID: PMC10250017 DOI: 10.1002/wcms.1571] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Accepted: 07/20/2021] [Indexed: 12/16/2022]
Abstract
With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.
Collapse
Affiliation(s)
- Yuanzhe Zhou
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Yangwei Jiang
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| | - Shi-Jie Chen
- Department of Physics and Astronomy, Department of Biochemistry, Institute of Data Sciences and Informatics, University of Missouri, Columbia, MO 65211-7010, USA
| |
Collapse
|
41
|
Williams AH, Zhan CG. Generalized Methodology for the Quick Prediction of Variant SARS-CoV-2 Spike Protein Binding Affinities with Human Angiotensin-Converting Enzyme II. J Phys Chem B 2022; 126:2353-2360. [PMID: 35315274 PMCID: PMC8982491 DOI: 10.1021/acs.jpcb.1c10718] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 02/24/2022] [Indexed: 12/25/2022]
Abstract
Variants of the SARS-CoV-2 virus continue to remain a threat 2 years from the beginning of the pandemic. As more variants arise, and the B.1.1.529 (Omicron) variant threatens to create another wave of infections, a method is needed to predict the binding affinity of the spike protein quickly and accurately with human angiotensin-converting enzyme II (ACE2). We present an accurate and convenient energy minimization/molecular mechanics Poisson-Boltzmann surface area methodology previously used with engineered ACE2 therapeutics to predict the binding affinity of the Omicron variant. Without any additional data from the variants discovered after the publication of our first model, the methodology can accurately predict the binding of the spike/ACE2 variant complexes. From this methodology, we predicted that the Omicron variant spike has a Kd of ∼22.69 nM (which is very close to the experimental Kd of 20.63 nM published during the review process of the current report) and that spike protein of the new "Stealth" Omicron variant (BA.2) will display a Kd of ∼12.9 nM with the wild-type ACE2 protein. This methodology can be used with as-yet discovered variants, allowing for quick determinations regarding the variant's infectivity versus either the wild-type virus or its variants.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
42
|
Xu Y, Huang SW, Ma YQ, Ding HM. Loading of DOX into a tetrahedral DNA nanostructure: the corner does matter. NANOSCALE ADVANCES 2022; 4:754-760. [PMID: 36131833 PMCID: PMC9416905 DOI: 10.1039/d1na00753j] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 12/05/2021] [Indexed: 06/03/2023]
Abstract
With the rapid development of nanotechnology, various DNA nanostructures have been synthesized and widely used in drug delivery. However, the underlying mechanisms of drug molecule loading into the DNA nanostructure are still elusive. In this work, we systematically investigate the interactions of a tetrahedral DNA nanostructure (TDN) with the anti-cancer drug doxorubicin (DOX) by combining molecular docking and all-atom molecular dynamics simulations. It is found that there are five possible binding modes in the single TDN-DOX interactions, namely the outside-corner mode, the inside-corner mode, the major-groove mode, the minor-groove mode, and the intercalation mode, where the van der Waals (VDW) interaction and the electrostatic (ELE) interaction dominate in the case of unionized DOX and ionized DOX, respectively. Moreover, with the increase of the DOX number, some of the interaction modes may disappear and the inside-corner mode is the most energy-favorable mode. The present study enhances the molecular understanding of the role of TDN as the drug carrier, which may provide a useful guideline for the future design of DNA nanostructures.
Collapse
Affiliation(s)
- Yao Xu
- National Laboratory of Solid State Microstructures and Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University Nanjing 210093 China
| | - Shu-Wei Huang
- National Laboratory of Solid State Microstructures and Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University Nanjing 210093 China
| | - Yu-Qiang Ma
- National Laboratory of Solid State Microstructures and Department of Physics, Collaborative Innovation Center of Advanced Microstructures, Nanjing University Nanjing 210093 China
| | - Hong-Ming Ding
- Center for Soft Condensed Matter Physics and Interdisciplinary Research, School of Physical Science and Technology, Soochow University Suzhou 215006 China
| |
Collapse
|
43
|
Mendis J, Kaya E, Kucukkal TG. Identification of Hotspot Residues in Binding of SARS-CoV-2 Spike and Human ACE2 Proteins. JOURNAL OF COMPUTATIONAL BIOPHYSICS AND CHEMISTRY 2021. [DOI: 10.1142/s2737416521500447] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Coronaviruses are a large family of viruses that can cause respiratory infections with varying severity from common cold to severe diseases such as novel coronavirus disease (COVID-19). COVID-19 has been declared as a global pandemic by the World Health Organization on March 11, 2020 and with the development of vaccines it slowed down as of this date. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uses its spike glycoprotein (Sgp) to bind human angiotensin-converting enzyme 2 (hACE2) receptor, and mediates membrane fusion and virus entry. The recognition of Sgp to human ACE2 and its high affinity for it has been of great importance since this provides the first step in viral entry to human cells. Therefore, it is crucial to identify key residues (hotspots) in this process. In this study, computational alanine scanning has been performed for Sgp and hACE2. The residues identified with significance in binding and other residues in close proximity were studied further through molecular mechanics-based protein binding free energy change prediction methods. Moreover, the interfacial residues in both proteins were investigated for their cooperative binding. Additionally, folding free energy changes upon mutation to Ala were calculated to assess their effect on stability of Sgp and hACE2. Our results taken together with findings from previous studies revealed the residues that are most significant and are relatively significant in binding of Sgp to human ACE2 protein.
Collapse
Affiliation(s)
- Jenny Mendis
- School of Science, Technology, Accessibility, Mathematics and Public Health, Gallaudet University, Washington, D.C. 20002, USA
| | - Ekrem Kaya
- Quest Student Research Institute, Chantilly, VA 20151, USA
- Freedom High School, Chantilly, VA 20152, USA
| | - Tugba G. Kucukkal
- School of Science, Technology, Accessibility, Mathematics and Public Health, Gallaudet University, Washington, D.C. 20002, USA
| |
Collapse
|
44
|
Feng Y, Yan Y, He J, Tao H, Wu Q, Huang SY. Docking and scoring for nucleic acid-ligand interactions: Principles and current status. Drug Discov Today 2021; 27:838-847. [PMID: 34718205 DOI: 10.1016/j.drudis.2021.10.013] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 09/06/2021] [Accepted: 10/20/2021] [Indexed: 12/24/2022]
Abstract
Nucleic acid (NA)-ligand interactions have crucial roles in many cellular processes and, thus, are increasingly attracting therapeutic interest in drug discovery. Molecular docking is a valuable tool for studying molecular interactions. However, because NAs differ significantly from proteins in both their physical and chemical properties, traditional docking algorithms and scoring functions for protein-ligand interactions might not be applicable to NA-ligand docking. Therefore, various sampling strategies and scoring functions for NA-ligand interactions have been developed. Here, we review the basic principles and current status of docking algorithms and scoring functions for DNA/RNA-ligand interactions. We also discuss challenges and limitations of current docking and scoring approaches.
Collapse
Affiliation(s)
- Yuyu Feng
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China
| | - Yumeng Yan
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China
| | - Jiahua He
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China
| | - Huanyu Tao
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China
| | - Qilong Wu
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China
| | - Sheng-You Huang
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, PR China.
| |
Collapse
|
45
|
King E, Aitchison E, Li H, Luo R. Recent Developments in Free Energy Calculations for Drug Discovery. Front Mol Biosci 2021; 8:712085. [PMID: 34458321 PMCID: PMC8387144 DOI: 10.3389/fmolb.2021.712085] [Citation(s) in RCA: 49] [Impact Index Per Article: 16.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 07/27/2021] [Indexed: 01/11/2023] Open
Abstract
The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.
Collapse
Affiliation(s)
- Edward King
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Erick Aitchison
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
| | - Han Li
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
| | - Ray Luo
- Department of Molecular Biology and Biochemistry, University of California, Irvine, CA, United States
- Department of Chemical and Biomolecular Engineering, University of California, Irvine, CA, United States
- Department of Materials Science and Engineering, University of California, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, CA, United States
| |
Collapse
|
46
|
Exploring the effect of temperature on inhibition of non-structural protease 3 of Chikungunya virus using molecular dynamics simulations and thermodynamics parameters. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116164] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|
47
|
Williams AH, Zhan CG. Fast Prediction of Binding Affinities of the SARS-CoV-2 Spike Protein Mutant N501Y (UK Variant) with ACE2 and Miniprotein Drug Candidates. J Phys Chem B 2021; 125:4330-4336. [PMID: 33881861 PMCID: PMC8084269 DOI: 10.1021/acs.jpcb.1c00869] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2021] [Revised: 03/30/2021] [Indexed: 12/28/2022]
Abstract
A recently identified variant of SARS-CoV-2 virus, known as the United Kingdom (UK) variant (lineage B.1.1.7), has an N501Y mutation on its spike protein. SARS-CoV-2 spike protein binds with angiotensin-converting enzyme 2 (ACE2), a key protein for the viral entry into the host cells. Here, we report an efficient computational approach, including the simple energy minimizations and binding free energy calculations, starting from an experimental structure of the binding complex along with experimental calibration of the calculated binding free energies, to rapidly and reliably predict the binding affinities of the N501Y mutant with human ACE2 (hACE2) and recently reported miniprotein and hACE2 decoy (CTC-445.2) drug candidates. It has been demonstrated that the N501Y mutation markedly increases the ACE2-spike protein binding affinity (Kd) from 22 to 0.44 nM, which could partially explain why the UK variant is more infectious. The miniproteins are predicted to have ∼10,000- to 100,000-fold diminished binding affinities with the N501Y mutant, creating a need for design of novel therapeutic candidates to overcome the N501Y mutation-induced drug resistance. The N501Y mutation is also predicted to decrease the binding affinity of a hACE2 decoy (CTC-445.2) binding with the spike protein by ∼200-fold. This convenient computational approach along with experimental calibration may be similarly used in the future to predict the binding affinities of potential new variants of the spike protein.
Collapse
Affiliation(s)
- Alexander H. Williams
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| | - Chang-Guo Zhan
- Molecular Modeling and Biopharmaceutical Center, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
- Department of Pharmaceutical Sciences, College of Pharmacy, University of Kentucky, 789 South Limestone Street, Lexington, KY 40536
| |
Collapse
|
48
|
Jiang Y, Liu HF, Liu R. Systematic comparison and prediction of the effects of missense mutations on protein-DNA and protein-RNA interactions. PLoS Comput Biol 2021; 17:e1008951. [PMID: 33872313 PMCID: PMC8084330 DOI: 10.1371/journal.pcbi.1008951] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Revised: 04/29/2021] [Accepted: 04/08/2021] [Indexed: 12/30/2022] Open
Abstract
The binding affinities of protein-nucleic acid interactions could be altered due to missense mutations occurring in DNA- or RNA-binding proteins, therefore resulting in various diseases. Unfortunately, a systematic comparison and prediction of the effects of mutations on protein-DNA and protein-RNA interactions (these two mutation classes are termed MPDs and MPRs, respectively) is still lacking. Here, we demonstrated that these two classes of mutations could generate similar or different tendencies for binding free energy changes in terms of the properties of mutated residues. We then developed regression algorithms separately for MPDs and MPRs by introducing novel geometric partition-based energy features and interface-based structural features. Through feature selection and ensemble learning, similar computational frameworks that integrated energy- and nonenergy-based models were established to estimate the binding affinity changes resulting from MPDs and MPRs, but the selected features for the final models were different and therefore reflected the specificity of these two mutation classes. Furthermore, the proposed methodology was extended to the identification of mutations that significantly decreased the binding affinities. Extensive validations indicated that our algorithm generally performed better than the state-of-the-art methods on both the regression and classification tasks. The webserver and software are freely available at http://liulab.hzau.edu.cn/PEMPNI and https://github.com/hzau-liulab/PEMPNI. Protein-nucleic acid interactions play important roles in various cellular processes. Missense mutations occurring in DNA- or RNA-binding proteins (termed MPDs and MPRs, respectively) could change the binding affinities of these interactions. Previous studies have compared protein-DNA and protein-RNA interactions from multifaceted viewpoints, but less attention has been given to the similarities and specific differences between the effects of MPDs and MPRs and between the methodologies for predicting the affinity changes induced by the two mutation classes. Therefore, we systematically compared their impacts and demonstrated that MPDs and MPRs could have specific preferences for binding affinity changes. These observations motivated us to construct regression models separately for MPDs and MPRs by introducing novel energy and nonenergy descriptors. Although similar frameworks were developed to estimate these two categories of mutation effects, different descriptors were selected in the regression models and further revealed the specificity of mutation classes. The interplay between the energy and nonenergy modules effectively improved prediction performance. Our algorithm can also be adopted to disentangle mutations significantly decreasing binding affinities from other mutations.
Collapse
Affiliation(s)
- Yao Jiang
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, P. R. China
| | - Hui-Fang Liu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, P. R. China
| | - Rong Liu
- Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan, P. R. China
| |
Collapse
|
49
|
Zhang X, Zhang Z, Guo J, Ma J, Xie S, Zhao Y, Wang C. Combination of multiple computational methods revealing specific sub-sectional recognition and hydrogen-bond dependent transportation of CKII peptide fragment in O-GlcNAc transferase. Comput Struct Biotechnol J 2021; 19:2045-2056. [PMID: 33995901 PMCID: PMC8085782 DOI: 10.1016/j.csbj.2021.04.009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 04/02/2021] [Accepted: 04/03/2021] [Indexed: 11/17/2022] Open
Abstract
Mechanism of CKII peptide recognition, transportation and binding in OGT is obtained. Peptide delivery is strong exothermic, highly dependent on hydrogen bond network. Typical ‘spread’ & ‘V’ conformation change noticed for peptide accompanies stable OGT. Specific subsection of peptide has diverse performance in its recognition and delivery. Multiple methods combination may be used in other bio-system with flexible substrate.
O-linked β-N-acetyl-D-glucosamine (O-GlcNAc) transferase (OGT) is an essential enzyme in many cellular physiological catalytic reactions that regulates protein O-GlcNAcylation. Aberrant O-GlcNAcylation is related to insulin resistance, diabetic complications, cancer and neurodegenerative diseases. Understanding the peptide delivery in OGT is significant in comprehending enzymatic catalytic process, target-protein recognition and pathogenic mechanism. Herein extensive molecular dynamics (MD) simulations combined with various techniques are utilized to study the recognizing and binding mechanism of peptide fragment extracted from casein kinase II by OGT from atomic level. The residues of His496, His558, Thr633, Lys634, and Pro897 are demonstrated to play a dominant role in the peptide stabilization via hydrogen bonds and σ-π interaction, whose van der Waals and non-polar solvent effects provide the main driving force. In addition, two channels are identified. The delivery mode, mechanism together with thermodynamic and dynamic characterizations for the most favorable channel are determined. The peptide is more inclined to be recognized by OGT through the cavity comprised of residues 799–812, 893–899, and 865–871, and Tyr13-terminal is prior recognized to Met26-terminal. The transportation process is accompanied with conformation changes between the “spread” and “V” shapes. The whole process is strong exothermic that is highly dependent on the variation of hydrogen bond interactions between peptide and OGT as well as the performance of different subsections of peptide. Besides that, multiple computational methods combinations may contribute meaningfully to calculation of similar bio-systems with long and flexible substrate.
Collapse
Affiliation(s)
- Xiao Zhang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng 475004, People's Republic of China
| | - Zhiyang Zhang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng 475004, People's Republic of China
| | - Jia Guo
- College of Chemistry and Chemical Engineering, Henan University, Kaifeng 475004, People's Republic of China
| | - Jing Ma
- School of Pharmacy, Henan University, Kaifeng 475004, People's Republic of China
| | - Songqiang Xie
- School of Pharmacy, Henan University, Kaifeng 475004, People's Republic of China
| | - Yuan Zhao
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng 475004, People's Republic of China
| | - Chaojie Wang
- The Key Laboratory of Natural Medicine and Immuno-Engineering, Henan University, Kaifeng 475004, People's Republic of China
| |
Collapse
|
50
|
Schizandrin A exhibits potent anticancer activity in colorectal cancer cells by inhibiting heat shock factor 1. Biosci Rep 2021; 40:222263. [PMID: 32110802 PMCID: PMC7069920 DOI: 10.1042/bsr20200203] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2020] [Revised: 02/18/2020] [Accepted: 02/19/2020] [Indexed: 02/06/2023] Open
Abstract
Heat shock factor 1 (HSF1) is a powerful multifaceted oncogenic modifier that plays a role in maintaining the protein balance of cancer cells under various stresses. In recent studies, there have been reports of increased expression of HSF1 in colorectal cancer (CRC) cells, and the depletion of the HSF1 gene knockdown has inhibited colon cancer growth both in vivo and in vitro. Therefore, HSF1 is a promising target for colon cancer treatment and chemoprevention. In the present study, we found that Schizandrin A (Sch A) significantly inhibited the growth of CRC cell lines by inducing cell cycle arrest, apoptosis and death. Through HSE luciferase reporter assay and quantitative PCR (qPCR), we identified Sch A as a novel HSF1 inhibitor. In addition, Sch A could effectively inhibit the induction of HSF1 target proteins such as heat-shock protein (HSP) 70 (HSP70) and HSP27, whether in heat shock or normal temperature culture. In the Surface Plasmon Resonance (SPR) experiment, Sch A showed moderate affinity with HSF1, further confirming that Sch A might be a direct HSF1 inhibitor. The molecular docking and molecular dynamic simulation results of HSF1/Sch A suggested that Sch A formed key hydrogen bond and hydrophobic interactions with HSF1, which may contribute to its potent HSF1 inhibition. These findings provide clues for the design of novel HSF1 inhibitors and drug candidates for colon cancer treatment.
Collapse
|