1
|
He Y, Zhu Y, Yin Z, Shi J, Shang K, Tian T, Shi H, Ding J, Zhang F. Design a novel of Brucellosis preventive vaccine based on IgV_CTLA-4 and multiple epitopes via immunoinformatics approach. Microb Pathog 2024:106909. [PMID: 39218373 DOI: 10.1016/j.micpath.2024.106909] [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: 04/13/2024] [Revised: 06/24/2024] [Accepted: 08/29/2024] [Indexed: 09/04/2024]
Abstract
Brucellosis is a zoonotic disease caused by Brucella, which is difficult to be eliminated by conventional drugs. Therefore, a novel multi-epitope vaccine (MEV) was designed to prevent human Brucella infection. Based on the method of "reverse vaccinology", cytotoxic T lymphocytes epitopes (CTLEs), helper T lymphocytes epitopes (HTLEs), linear B-cell epitopes (LBEs) and conformational B-cell epitopes (CBEs) of four Brucella proteins (VirB9, VirB10, Omp 19 and Omp 25) were obtained. In order to keep the correct protein folding, the multiple epitopes was constructed by connecting epitopes through linkers. In view of the significant connection between human leukocyte antigen CTLA-4 and B7 molecules found on antigen presenting cells (APCs), a new vaccine (V_C4MEV) for preventing brucellosis was created by combining CTLA-4 immunoglobulin variable region (IgV_CTLA-4) with MEV protein. Immunoinformatics analysis showed that V_C4MEV had a good secondary and tertiary structure. Additionally, molecular docking and molecular dynamics simulation (MD) revealed a robust binding affinity between IgV_ CTLA-4 and the B7 molecule. Notably, the vaccine V_C4MEV was demonstrated favorable immunogenicity and antigenicity in both in vitro and in vivo experiments. V_C4MEV had the potential to activate defensive cells and immune responses, offering a hopeful approach for developing vaccines against Brucella in the upcoming years.
Collapse
Affiliation(s)
- Yueyue He
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China; Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - YueJie Zhu
- Department of Reproductive Assistance, Center for Reproductive Medicine, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, Xinjiang, China; The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Zhengwei Yin
- The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Juan Shi
- The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Kaiyu Shang
- The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Tingting Tian
- The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Huidong Shi
- The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Jianbing Ding
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China; State Key Laboratory of Pathogenesis, Prevention, Treatment of Central Asian High Incidence Diseases, the First Affiliated Hospital of Xinjiang Medical University, Urumqi, China.
| | - Fengbo Zhang
- Department of Clinical Laboratory, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China; State Key Laboratory of Pathogenesis, Prevention, Treatment of Central Asian High Incidence Diseases, the First Affiliated Hospital of Xinjiang Medical University, Urumqi, China.
| |
Collapse
|
2
|
Roy A, Paul I, Paul T, Hazarika K, Dihidar A, Ray S. An in-silico receptor-pharmacophore based multistep molecular docking and simulation study to evaluate the inhibitory potentials against NS1 of DENV-2. J Biomol Struct Dyn 2024; 42:6136-6164. [PMID: 37517062 DOI: 10.1080/07391102.2023.2239925] [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/01/2023] [Accepted: 06/25/2023] [Indexed: 08/01/2023]
Abstract
DENV-2 strain is the most fatal and infectious of the five dengue virus serotypes. The non-structural protein NS1 encoded by its genome is the most significant protein required for viral pathogenesis and replication inside the host body. Thus, targeting the NS1 protein and designing an inhibitor to limit its stability and secretion is a propitious attempt in our fight against dengue. Four novel inhibitors are designed to target the conserved cysteine residues (C55, C313, C316, and C329) and glycosylation sites (N130 and N207) of the NS1 protein in an attempt to halt the spread of the dengue infection in the host body altogether. Numerous computer-aided drug designing techniques including molecular docking, molecular dynamics simulation, virtual screening, principal component analysis, and dynamic cross-correlation matrix were employed to determine the structural and functional activity of the NS1-inhibitor complexes. From our analysis, it was evident that the extent of structural and atomic level fluctuations of the ligand-bound protein exhibited a declining trend in contrast to unbound protein which was prominently noticeable through the RMSD, RMSF, Rg, and SASA graphs. The ADMET analysis of the four ligands revealed a promising pharmacokinetics and pharmacodynamic profile, along with good bioavailability and toxicity properties. The proposed drugs when bound to the targeted cavities resulted in stable conformations in comparison to their unbound state, implying they have good affinity promising effective drug action. Thus, they can be tested in vitro and used as potential anti-dengue drugs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Alankar Roy
- Amity Institute of Biotechnology, Amity University, Kolkata, India
| | - Ishani Paul
- Amity Institute of Biotechnology, Amity University, Kolkata, India
| | - Tanwi Paul
- Amity Institute of Biotechnology, Amity University, Kolkata, India
| | | | - Aritrika Dihidar
- Amity Institute of Biotechnology, Amity University, Kolkata, India
| | - Sujay Ray
- Amity Institute of Biotechnology, Amity University, Kolkata, India
| |
Collapse
|
3
|
Zhu Y, Yu M, Aisikaer M, Zhang C, He Y, Chen Z, Yang Y, Han R, Li Z, Zhang F, Ding J, Lu X. Contriving a novel of CHB therapeutic vaccine based on IgV_CTLA-4 and L protein via immunoinformatics approach. J Biomol Struct Dyn 2024; 42:6323-6341. [PMID: 37424209 DOI: 10.1080/07391102.2023.2234043] [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: 01/07/2023] [Accepted: 07/01/2023] [Indexed: 07/11/2023]
Abstract
Chronic infection induced by immune tolerance to hepatitis B virus (HBV) is one of the most common causes of hepatic cirrhosis and hepatoma. Fortunately, the application of therapeutic vaccine can not only reverse HBV-tolerance, but also serve a potentially effective therapeutic strategy for treating chronic hepatitis B (CHB). However, the clinical effect of the currently developed CHB therapeutic vaccine is not optimistic due to the weak immunogenicity. Given that the human leukocyte antigen CTLA-4 owns strong binding ability to the surface B7 molecules (CD80 and CD86) of antigen presenting cell (APCs), the immunoglobulin variable region of CTLA-4 (IgV_CTLA-4) was fused with the L protein of HBV to contrive a novel therapeutic vaccine (V_C4HBL) for CHB in this study. We found that the addition of IgV_CTLA-4 did not interfere with the formation of L protein T cell and B cell epitopes after analysis by means of immunoinformatics approaches. Meanwhile, we also found that the IgV_CTLA-4 had strong binding force to B7 molecules through molecular docking and molecular dynamics (MD) simulation. Notably, our vaccine V_C4HBL showed good immunogenicity and antigenicity by in vitro and in vivo experiments. Therefore, the V_C4HBL is promising to again effectively activate the cellular and humoral immunity of CHB patients, and provides a potentially effective therapeutic strategy for the treatment of CHB in the future.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Yuejie Zhu
- Reproductive Medicine Center, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
- Infectious Disease Center, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Mingkai Yu
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - Maierhaba Aisikaer
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - Chuntao Zhang
- Department of Microbiology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
| | - Yueyue He
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - Zhiqiang Chen
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - Yinyin Yang
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
- Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
| | - Rui Han
- Reproductive Medicine Center, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Zhiwei Li
- Clinical Laboratory Center, Xinjiang Uygur Autonomous Region People's Hospital, Urumqi, China
| | - Fengbo Zhang
- Department of Clinical Laboratory, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| | - Jianbing Ding
- Department of Immunology, School of Basic Medical Sciences, Xinjiang Medical University, Urumqi, China
- Xinjiang Key Molecular Biology Laboratory of Endemic Disease, Xinjiang Medical University, Urumqi, China
| | - Xiaobo Lu
- Infectious Disease Center, The First Affiliated Hospital of Xinjiang Medical University, Urumqi, China
| |
Collapse
|
4
|
Venanzi NE, Basciu A, Vargiu AV, Kiparissides A, Dalby PA, Dikicioglu D. Machine Learning Integrating Protein Structure, Sequence, and Dynamics to Predict the Enzyme Activity of Bovine Enterokinase Variants. J Chem Inf Model 2024; 64:2681-2694. [PMID: 38386417 PMCID: PMC11005043 DOI: 10.1021/acs.jcim.3c00999] [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] [Received: 07/03/2023] [Revised: 02/12/2024] [Accepted: 02/13/2024] [Indexed: 02/24/2024]
Abstract
Despite recent advances in computational protein science, the dynamic behavior of proteins, which directly governs their biological activity, cannot be gleaned from sequence information alone. To overcome this challenge, we propose a framework that integrates the peptide sequence, protein structure, and protein dynamics descriptors into machine learning algorithms to enhance their predictive capabilities and achieve improved prediction of the protein variant function. The resulting machine learning pipeline integrates traditional sequence and structure information with molecular dynamics simulation data to predict the effects of multiple point mutations on the fold improvement of the activity of bovine enterokinase variants. This study highlights how the combination of structural and dynamic data can provide predictive insights into protein functionality and address protein engineering challenges in industrial contexts.
Collapse
Affiliation(s)
| | - Andrea Basciu
- Department
of Physics, University of Cagliari, Cittadella
Universitaria, I-09042 Monserrato, Cagliari, Italy
| | - Attilio Vittorio Vargiu
- Department
of Physics, University of Cagliari, Cittadella
Universitaria, I-09042 Monserrato, Cagliari, Italy
| | - Alexandros Kiparissides
- Department
of Biochemical Engineering, University College
London, Gower Street, WC1E 6BT London, U.K.
- Department
of Chemical Engineering, Aristotle University
of Thessaloniki, 54 124 Thessaloniki, Greece
| | - Paul A. Dalby
- Department
of Biochemical Engineering, University College
London, Gower Street, WC1E 6BT London, U.K.
| | - Duygu Dikicioglu
- Department
of Biochemical Engineering, University College
London, Gower Street, WC1E 6BT London, U.K.
| |
Collapse
|
5
|
Alotayeq A, Ghannay S, Alhagri IA, Ahmed I, Hammami B, E. A. E. Albadri A, Patel H, Messaoudi S, Kadri A, M. Al-Hazmy S, Aouadi K. Synthesis, optical properties, DNA, β-cyclodextrin interaction, hydrogen isotope sensor and computational study of new enantiopure isoxazolidine derivative (ISoXD). Heliyon 2024; 10:e26341. [PMID: 38404822 PMCID: PMC10884473 DOI: 10.1016/j.heliyon.2024.e26341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 02/12/2024] [Accepted: 02/12/2024] [Indexed: 02/27/2024] Open
Abstract
A novel isoxazolidine derivative (ISoXD) dye was successfully synthesized and comprehensively characterized. In this study, we conducted a thorough examination of its various properties, including optical characteristics, interactions with DNA and β-cyclodextrin (β-CD), molecular docking, molecular dynamic simulation, and density functional theory (DFT) calculations. Our investigation encompassed a systematic analysis of the absorption and emission spectra of ISoXD in diverse solvents. The observed variations in the spectroscopic data were attributed to the specific solvent's capacity to engage in hydrogen bonding interactions. Remarkably, the most pronounced intensities were observed in glycol, which can establish many hydrogen bonds with ISoXD. Furthermore, our study revealed a significant distinction in the fluorescence behavior of ISoXD when subjected to different solvents, particularly between CHCl3 and CDCl3. Moreover, we explored the fluorescence intensity of the ISoXD complex in the presence of various metals, both in ethanol and water. The ISoXD complex exhibited a substantial increase of fluorescence upon interaction with different metal ions. The utilization of DFT calculations allowed us to propose an intramolecular charge transfer (ICT) mechanism as a plausible explanation for this quenching phenomenon. The interaction of ISoXD with DNA and β-CD was studied using absorption spectra. The binding constant (K) and the standard Gibbs free energy change (ΔGo) for the interaction between DNA and β-CD with ISoXD were determined. In docking study, ISoXD exhibited significant docking scores (-6.511) and MM-GBSA binding free energies (-66.27 kcal/mol) within the PARP-1 binding cavity. Its binding pattern closely resembles to the co-crystal ligand veliparib, and during a 100ns MD simulation, ISoXD displayed strong stability and formed robust hydrogen bonds with key amino acids. These findings suggest ISoXD's potential as a PARP-1 inhibitor for further investigation in therapeutic development.
Collapse
Affiliation(s)
- Afnan Alotayeq
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
| | - Siwar Ghannay
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
| | - Ibrahim A. Alhagri
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
- Department of Chemistry, Faculty of Sciences, Ibb University, Ibb, Yemen
| | - Iqrar Ahmed
- Department of Pharmaceutical Chemistry, Prof. Ravindra Nikam College of Pharmacy, Gondur, Dhule, 424002, India
- Division of Computer Aided Drug Design, Department of Pharmaceutical Chemistry, R. C. Patel Institute of Pharmaceutical Education and Research, Shirpur, 425405, Maharashtra, India
| | - Bechir Hammami
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
| | - Abuzar E. A. E. Albadri
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
| | - Harun Patel
- Division of Computer Aided Drug Design, Department of Pharmaceutical Chemistry, R. C. Patel Institute of Pharmaceutical Education and Research, Shirpur, 425405, Maharashtra, India
| | - Sabri Messaoudi
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
- Faculty of Sciences of Bizerte, Carthage University, Jarzouna, Bizerte 7021, Tunisia
| | - Adel Kadri
- Faculty of Science of Sfax, Department of Chemistry, University of Sfax, B.P. 1171, 3000 Sfax, Tunisia
- Department of Chemistry, Faculty of Science and Arts of Baljurashi, Al- Baha University, Saudi Arabia
| | - Sadeq M. Al-Hazmy
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
| | - Kaiss Aouadi
- Department of Chemistry, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
- Department of Chemistry, Laboratory of Heterocyclic Chemistry Natural Product and Reactivity/CHPNR, Faculty of Science of Monastir, University of Monastir, Avenue of the Environment, Monastir, 5019, Tunisia
| |
Collapse
|
6
|
Aulifa DL, Al Shofwan AA, Megantara S, Fakih TM, Budiman A. Elucidation of Molecular Interactions Between Drug-Polymer in Amorphous Solid Dispersion by a Computational Approach Using Molecular Dynamics Simulations. Adv Appl Bioinform Chem 2024; 17:1-19. [PMID: 38282640 PMCID: PMC10821732 DOI: 10.2147/aabc.s441628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 01/16/2024] [Indexed: 01/30/2024] Open
Abstract
Introduction Amorphous drug dispersion is frequently used to enhance the solubility and dissolution of poorly water-soluble drugs, thereby improving their oral bioavailability. The dispersion of these drugs into polymer matrix can inhibit their recrystallization. The inter-molecular interactions between drug and polymer plays a role in the improvement of the dissolution rate, solubility, and physical stability of drug. Aim This study aims to investigate the formation and interactions of ritonavir (RTV)/poloxamer (PLX) amorphous formulation using a computational approach via molecular dynamics (MD) simulations, which mimicked solvent evaporation and melt-quenching method. Methods TheRoot Mean Square Deviation (RMSD) value, Root Mean Square Fluctuation (RMSF), Radial Distribution Function (RDF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), and hydrogen bond interactions were analyzed to determine interaction mechanisms between RTV and PLX in amorphous solid dispersion. Results The pi-alkyl bonds between RTV and PLX were formed after simulations of solvent evaporation, while the hydrogen bond interactions of RTV-PLX was observed during melt method simulations. These results indicate the successful formulation of amorphous solid dispersion (ASD) from RTV and PLX. The RMSD values obtained from the solvent evaporation, melt-cooling-A, melt-cooling-B, and melt-cooling-C methods were 3.33 Å, 1.97 Å, 1.30 Å, and 1.29 Å, respectively, while the average RMSF values were 2.65 Å, 1.04 Å, 1.05 Å, and 1.07 Å, respectively. This indicates that the suppression of translational motion of RTV from the melt method can be stronger than solvent evaporation caused by the intermolecular interactions of RTV-PLX. Conclusion MD simulations helped in understanding the formation and interaction mechanisms of ASD formulations that were difficult to detect by experimental approaches.
Collapse
Affiliation(s)
- Diah Lia Aulifa
- Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
| | - Adnan Aly Al Shofwan
- Department of Pharmaceutics and Pharmaceutical Technology, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
| | - Sandra Megantara
- Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
| | - Taufik Muhammad Fakih
- Department of Pharmacy, Faculty of Mathematics and Natural Sciences, Universitas Islam Bandung, Bandung, Indonesia
| | - Arif Budiman
- Department of Pharmaceutics and Pharmaceutical Technology, Faculty of Pharmacy, Universitas Padjadjaran, Sumedang, Indonesia
| |
Collapse
|
7
|
Kumar A, Dutt M, Dehury B, Sganzerla Martinez G, Swan CL, Kelvin AA, Richardson CD, Kelvin DJ. Inhibition potential of natural flavonoids against selected omicron (B.1.19) mutations in the spike receptor binding domain of SARS-CoV-2: a molecular modeling approach. J Biomol Struct Dyn 2023:1-15. [PMID: 38115191 DOI: 10.1080/07391102.2023.2291165] [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: 05/29/2023] [Accepted: 09/09/2023] [Indexed: 12/21/2023]
Abstract
The omicron (B.1.19) variant of contagious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is considered a variant of concern (VOC) due to its increased transmissibility and highly infectious nature. The spike receptor-binding domain (RBD) is a hotspot of mutations and is regarded as a prominent target for screening drug candidates owing to its crucial role in viral entry and immune evasion. To date, no effective therapy or antivirals have been reported; therefore, there is an urgent need for rapid screening of antivirals. An extensive molecular modelling study has been performed with the primary goal to assess the inhibition potential of natural flavonoids as inhibitors against RBD from a manually curated library. Out of 40 natural flavonoids, five natural flavonoids, namely tomentin A (-8.7 kcal/mol), tomentin C (-8.6 kcal/mol), hyperoside (-8.4 kcal/mol), catechin gallate (-8.3 kcal/mol), and corylifol A (-8.2 kcal/mol), have been considered as the top-ranked compounds based on their binding affinity and molecular interaction profiling. The state-of-the-art molecular dynamics (MD) simulations of these top-ranked compounds in complex with RBD exhibited stable dynamics and structural compactness patterns on 200 nanoseconds. Additionally, complexes of these molecules demonstrated favorable free binding energies and affirmed the docking and simulation results. Moreover, the post-simulation validation of these interacted flavonoids using principal component analysis (PCA) revealed stable interaction patterns with RBD. The integrated results suggest that tomentin A, tomentin C, hyperoside, catechin gallate, and corylifol A might be effective against the emerging variants of SARS-CoV-2 and should be further evaluated using in-vitro and in-vivo experiments.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anuj Kumar
- Laboratory of Immunity, Shantou University Medical College, Shantou, China
- Department of Microbiology and Immunology, Faculty of Medicine, Dalhousie University, Halifax, Canada
- Department of Paediatrics, IWK Health Center, Canadian Centre for Vaccinology (CCfV), Halifax, Canada
| | - Mansi Dutt
- Laboratory of Immunity, Shantou University Medical College, Shantou, China
- Department of Microbiology and Immunology, Faculty of Medicine, Dalhousie University, Halifax, Canada
- Department of Paediatrics, IWK Health Center, Canadian Centre for Vaccinology (CCfV), Halifax, Canada
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, India
| | - Gustavo Sganzerla Martinez
- Laboratory of Immunity, Shantou University Medical College, Shantou, China
- Department of Microbiology and Immunology, Faculty of Medicine, Dalhousie University, Halifax, Canada
- Department of Paediatrics, IWK Health Center, Canadian Centre for Vaccinology (CCfV), Halifax, Canada
| | - Cynthia L Swan
- Vaccine and Infectious Disease Organization (VIDO), University of Saskatchewan, Saskatoon, Canada
| | - Alyson A Kelvin
- Vaccine and Infectious Disease Organization (VIDO), University of Saskatchewan, Saskatoon, Canada
- Department of Biochemistry, Microbiology, and Immunology, University of Saskatchewan, Saskatoon, Canada
| | - Christopher D Richardson
- Department of Microbiology and Immunology, Faculty of Medicine, Dalhousie University, Halifax, Canada
- Department of Paediatrics, IWK Health Center, Canadian Centre for Vaccinology (CCfV), Halifax, Canada
| | - David J Kelvin
- Laboratory of Immunity, Shantou University Medical College, Shantou, China
- Department of Microbiology and Immunology, Faculty of Medicine, Dalhousie University, Halifax, Canada
- Department of Paediatrics, IWK Health Center, Canadian Centre for Vaccinology (CCfV), Halifax, Canada
| |
Collapse
|
8
|
Adeosun IJ, Baloyi I, Aljoundi AK, Salifu EY, Ibrahim MA, Cosa S. Molecular modelling of SdiA protein by selected flavonoid and terpenes compounds to attenuate virulence in Klebsiella pneumoniae. J Biomol Struct Dyn 2023; 41:9938-9956. [PMID: 36416609 DOI: 10.1080/07391102.2022.2148753] [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: 02/15/2022] [Accepted: 11/12/2022] [Indexed: 11/25/2022]
Abstract
Klebsiella pneumoniae is one of the perturbing multidrug resistant (MDR) and ESKAPE pathogens contributing to the mounting morbidity, mortality and extended rate of hospitalization. Its virulence, often regulated by quorum sensing (QS) reinforces the need to explore alternative and prospective antivirulence agents, relatively from plants secondary metabolites. Computer aided drug discovery using molecular modelling techniques offers advantage to investigate prospective drugs to combat MDR pathogens. Thus, this study employed virtual screening of selected terpenes and flavonoids from medicinal plants to interrupt the QS associated SdiA protein in K. pneumoniae to attenuate its virulence. 4LFU was used as a template to model the structure of SdiA. ProCheck, Verify3D, Ramachandran plot scores, and ProSA-Web all attested to the model's good quality. Since SdiA protein in K. pneumoniae leads to the expression of virulence, 31 prospective bioactive compounds were docked for antagonistic potential. The stability of the protein-ligand complex, atomic motions and inter-atomic interactions were further investigated through molecular dynamics simulations (MDS) at 100 ns production runs. The binding free energy was estimated using the molecular mechanics/poisson-boltzmann surface area (MM/PB-SA). Furthermore, the drug-likeness properties of the studied compounds were validated. Docking studies showed phytol possesses the highest binding affinity (-9.205 kcal/mol) while glycitein had -9.752 kcal/mol highest docking score. The MDS of the protein in complex with the best-docked compounds revealed phytol with the highest binding energy of -44.2625 kcal/mol, a low root-mean-square deviation (RMSD) value of 1.54 Å and root-mean-square fluctuation (RMSF) score of 1.78 Å. Analysis of the drug-likeness properties prediction and bioavailability of these compounds revealed their conformed activity to lipinski's rules with bioavailability scores of 0.55 F. The studied terpenes and flavonoids compounds effectively thwart SdiA protein, therefore regulate inter- or intra cellular communication and associated in virulence Enterobacteriaceae, serving as prospective antivirulence drugs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Idowu Jesulayomi Adeosun
- Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Hatfield, South Africa
| | - Itumeleng Baloyi
- Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Hatfield, South Africa
| | - Aimen K Aljoundi
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Elliasu Y Salifu
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | | | - Sekelwa Cosa
- Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Hatfield, South Africa
| |
Collapse
|
9
|
Murali M, Ahmed F, Gowtham HG, Aribisala JO, Abdulsalam RA, Shati AA, Alfaifi MY, Sayyed RZ, Sabiu S, Amruthesh KN. Exploration of CviR-mediated quorum sensing inhibitors from Cladosporium spp. against Chromobacterium violaceum through computational studies. Sci Rep 2023; 13:15505. [PMID: 37726386 PMCID: PMC10509224 DOI: 10.1038/s41598-023-42833-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2023] [Accepted: 09/15/2023] [Indexed: 09/21/2023] Open
Abstract
An opportunistic human pathogenic bacterium, Chromobacterium violaceum resists the potency of most antibiotics by exploiting the quorum sensing system within their community to control virulence factor expression. Therefore, blocking the quorum sensing mechanism could help to treat several infectious caused by this organism. The quorum sensing receptor (CviR) of C. violaceum was used as a model target in the current investigation to identify potentially novel quorum sensing inhibitors from Cladosporium spp. through in silico computational approaches. The molecular docking results confirmed the anti-quorum sensing potential of bioactive compounds from Cladosporium spp. through binding to CviR with varying docking scores between - 5.2 and - 9.5 kcal/mol. Relative to the positive control [Azithromycin (- 7.4 kcal/mol)], the top six metabolites of Cladosporium spp. had higher docking scores and were generally greater than - 8.5 kcal/mol. The thermodynamic stability and binding affinity refinement of top-ranked CviR inhibitors were further studied through a 160 ns molecular dynamic (MD) simulation. The Post-MD simulation analysis confirmed the top-ranked compounds' affinity, stability, and biomolecular interactions with CviR at 50 ns, 100 ns, and 160 ns with Coniochaetone K of the Cladosporium spp. having the highest binding free energy (- 30.87 kcal/mol) and best interactions (two consistent hydrogen bond contact) following the 160 ns simulation. The predicted pharmacokinetics properties of top selected compounds point to their drug likeliness, potentiating their chance as a possible drug candidate. Overall, the top-ranked compounds from Cladosporium spp., especially Coniochaetone K, could be identified as potential C. violaceum CviR inhibitors. The development of these compounds as broad-spectrum antibacterial medicines is thus possible in the future following the completion of further preclinical and clinical research.
Collapse
Affiliation(s)
- Mahadevamurthy Murali
- Applied Plant Pathology Laboratory, Department of Studies in Botany, University of Mysore, Manasagangotri, Mysore, 570006, India
| | - Faiyaz Ahmed
- Department of Clinical Nutrition, College of Applied Health Sciences in Ar Rass, Qassim University, 51452, Buraydah, Saudi Arabia
| | | | - Jamiu Olaseni Aribisala
- Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, Durban, South Africa
| | - Rukayat Abiola Abdulsalam
- Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, Durban, South Africa
| | - Ali A Shati
- Faculty of Science, Biology Department, King Khalid University, 9004, Abha, Saudi Arabia
| | - Mohammad Y Alfaifi
- Faculty of Science, Biology Department, King Khalid University, 9004, Abha, Saudi Arabia
| | - R Z Sayyed
- Department of Microbiology, PSGVP Mandal's S I Patil Arts, G B Patel Science and STKV Sangh Commerce College, Shahada, 425409, India.
| | - Saheed Sabiu
- Department of Biotechnology and Food Science, Faculty of Applied Sciences, Durban University of Technology, Durban, South Africa.
| | - Kestur Nagaraj Amruthesh
- Applied Plant Pathology Laboratory, Department of Studies in Botany, University of Mysore, Manasagangotri, Mysore, 570006, India.
| |
Collapse
|
10
|
Prabhakaran P, Hebbani AV, Menon SV, Paital B, Murmu S, Kumar S, Singh MK, Sahoo DK, Desai PPD. Insilico generation of novel ligands for the inhibition of SARS-CoV-2 main protease (3CL pro) using deep learning. Front Microbiol 2023; 14:1194794. [PMID: 37448573 PMCID: PMC10338188 DOI: 10.3389/fmicb.2023.1194794] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Accepted: 06/05/2023] [Indexed: 07/15/2023] Open
Abstract
The recent emergence of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing the coronavirus disease (COVID-19) has become a global public health crisis, and a crucial need exists for rapid identification and development of novel therapeutic interventions. In this study, a recurrent neural network (RNN) is trained and optimized to produce novel ligands that could serve as potential inhibitors to the SARS-CoV-2 viral protease: 3 chymotrypsin-like protease (3CLpro). Structure-based virtual screening was performed through molecular docking, ADMET profiling, and predictions of various molecular properties were done to evaluate the toxicity and drug-likeness of the generated novel ligands. The properties of the generated ligands were also compared with current drugs under various phases of clinical trials to assess the efficacy of the novel ligands. Twenty novel ligands were selected that exhibited good drug-likeness properties, with most ligands conforming to Lipinski's rule of 5, high binding affinity (highest binding affinity: -9.4 kcal/mol), and promising ADMET profile. Additionally, the generated ligands complexed with 3CLpro were found to be stable based on the results of molecular dynamics simulation studies conducted over a 100 ns period. Overall, the findings offer a promising avenue for the rapid identification and development of effective therapeutic interventions to treat COVID-19.
Collapse
Affiliation(s)
- Prejwal Prabhakaran
- Department of Biotechnology, New Horizon College of Engineering, Bangalore, India
- Faculty of Biology, Albert-Ludwigs-Universität Freiburg, Freiburg im Breisgau, Germany
| | - Ananda Vardhan Hebbani
- Department of Biochemistry, Indian Academy Degree College (Autonomous), Bangalore, India
| | - Soumya V. Menon
- Department of Chemistry and Biochemistry, School of Sciences, Jain (Deemed-to-be) University, Bangalore, India
| | - Biswaranjan Paital
- Redox Regulation Laboratory, Department of Zoology, College of Basic Science and Humanities, Odisha University of Agriculture and Technology, Bhubaneswar, India
| | - Sneha Murmu
- ICAR-Indian Agricultural Statistics Research Institute, PUSA, New Delhi, India
| | - Sunil Kumar
- ICAR-Indian Agricultural Statistics Research Institute, PUSA, New Delhi, India
| | | | - Dipak Kumar Sahoo
- Department of Veterinary Clinical Sciences, College of Veterinary Medicine, Iowa State University, Ames, IA, United States
| | | |
Collapse
|
11
|
Minh Quang N, Tran Thai H, Le Thi H, Duc Cuong N, Hien NQ, Hoang D, Ngoc VTB, Ky Minh V, Van Tat P. Novel Thiosemicarbazone Quantum Dots in the Treatment of Alzheimer's Disease Combining In Silico Models Using Fingerprints and Physicochemical Descriptors. ACS OMEGA 2023; 8:11076-11099. [PMID: 37008140 PMCID: PMC10061515 DOI: 10.1021/acsomega.2c07934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/13/2022] [Accepted: 03/07/2023] [Indexed: 06/19/2023]
Abstract
Searching for thiosemicarbazone derivatives with the potential to inhibit acetylcholinesterase for the treatment of Alzheimer's disease (AD) is an important current goal. The QSARKPLS, QSARANN, and QSARSVR models were constructed using binary fingerprints and physicochemical (PC) descriptors of 129 thiosemicarbazone compounds screened from a database of 3791 derivatives. The R 2 and Q 2 values for the QSARKPLS, QSARANN, and QSARSVR models are greater than 0.925 and 0.713 using dendritic fingerprint (DF) and PC descriptors, respectively. The in vitro pIC50 activities of four new design-oriented compounds N1, N2, N3, and N4, from the QSARKPLS model using DFs, are consistent with the experimental results and those from the QSARANN and QSARSVR models. The designed compounds N1, N2, N3, and N4 do not violate Lipinski-5 and Veber rules using the ADME and BoiLED-Egg methods. The binding energy, kcal mol-1, of the novel compounds to the 1ACJ-PDB protein receptor of the AChE enzyme was also obtained by molecular docking and dynamics simulations consistent with those predicted from the QSARANN and QSARSVR models. New compounds N1, N2, N3, and N4 were synthesized, and the experimental in vitro pIC50 activity was determined in agreement with those obtained from in silico models. The newly synthesized thiosemicarbazones N1, N2, N3, and N4 can inhibit 1ACJ-PDB, which is predicted to be able to cross the barrier. The DFT B3LYP/def-SV(P)-ECP quantization calculation method was used to calculate E HOMO and E LUMO to account for the activities of compounds N1, N2, N3, and N4. The quantum calculation results explained are consistent with those obtained in in silico models. The successful results here may contribute to the search for new drugs for the treatment of AD.
Collapse
Affiliation(s)
- Nguyen Minh Quang
- Faculty
of Chemical Engineering, Industrial University
of Ho Chi Minh City, 12 Nguyen Van Bao, Dist. Go Vap, Ho Chi Minh 700000, Viet Nam
| | - Hoa Tran Thai
- Faculty
of Chemistry, Hue University of Sciences, Hue University, 77 Nguyen Hue, Hue City 530000, Viet Nam
| | - Hoa Le Thi
- Faculty
of Chemistry, Hue University of Sciences, Hue University, 77 Nguyen Hue, Hue City 530000, Viet Nam
| | - Nguyen Duc Cuong
- Faculty
of Chemistry, Hue University of Sciences, Hue University, 77 Nguyen Hue, Hue City 530000, Viet Nam
- School
of Hospitality and Tourism, Hue University, 22 Lam Hoang, Hue City 530000, Viet
Nam
| | - Nguyen Quoc Hien
- Vietnam
Atomic Energy Institute, 59 Ly Thuong Kiet, Dist. Hoan Kiem, Hanoi
City 100000, Viet Nam
| | - DongQuy Hoang
- Faculty
of
Materials Science and Technology, University of Science, Vietnam National University, Ho Chi Minh 700000, Viet Nam
- Vietnam
National University, Ho Chi Minh
City 700000, Viet Nam
| | - Vu Thi Bao Ngoc
- Faculty
of Chemistry and Environment, University
of Dalat, 01 Phu Dong Thien Vuong, Dalat City 660000, Viet Nam
| | - Vo Ky Minh
- Franklin
High School, 6400 Whitelock Pkwy, Elk Grove, California 95757, United States
| | - Pham Van Tat
- Department
of Sciences and Journal Management, Hoa
Sen University, 08 Nguyen Van Trang, Dist. 01, Ho Chi Minh 700000, Viet Nam
| |
Collapse
|
12
|
Saah SA, Sakyi PO, Adu-Poku D, Boadi NO, Djan G, Amponsah D, Devine RNOA, Ayittey K. Docking and Molecular Dynamics Identify Leads against 5 Alpha Reductase 2 for Benign Prostate Hyperplasia Treatment. J CHEM-NY 2023. [DOI: 10.1155/2023/8880213] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023] Open
Abstract
Steroid 5 alpha-reductase 2 (5αR-2) is a membrane-embedded protein that together with other isoforms plays a key role in the metabolism of steroids. This enzyme catalyzes the reduction of testosterone to the more potent ligand, dihydrotestosterone (DHT) in the prostate. Androgens, testosterone, and DHT play important roles in prostate growth, development, and function. At the same time, both testosterone and DHT have been implicated in the pathogenesis of benign prostate hyperplasia (BPH). Inhibition of the DHT formation, therefore, provides a therapeutic strategy that offers the possibility of preventing, delaying, or treating BPH. Currently, two steroidal drugs that inhibit 5αR-2, dutasteride and finasteride, have been approved for clinical use. These two come at a high cost and also portray undesirable sexual side effects which necessitate the need to find new chemotherapeutic alternatives for the disease. Based on the aforementioned, finasteride and dutasteride were subjected to scaffold hopping, fragment-based de novo design, molecular docking, and molecular dynamics simulations employing databases like ChEMBL, DrugBank, PubChem, ChemSpider, and Zinc15 in the identification of potential hits targeting 5αR-2. Altogether, ten novel compounds targeting 5αR-2 were identified with binding energies lower or comparable to finasteride and dutasteride, the main inhibitors for this target. Molecular docking and molecular dynamics simulations studies identify amino acid residues Glu57, Phe219, Phe223, and Leu224 to be critical for ligand binding and complex stability. The physicochemical and pharmacological profiling suggests the potential of the hit compounds to be drug-like and orally active. Similarly, the quality parameter assessments revealed the hits possess LELP greater than 3 implying their promise as lead-like molecules. The compounds A5, A9, and A10 were, respectively, predicted to treat prostate disorders with Pa (0.188, 0.361, and 0.270) and Pi (0.176, 0.050, and 0.093), while A8 and A9 were found to be associated with BPH treatment with Pa (0.09 and 0.127) and Pi (0.077 and 0.033), respectively. Structural similarity searches via DrugBank identified the drugs faropenem, acemetacin, estradiol valerate, and yohimbine to be useful for BPH treatment suggesting the de novo designed ligands as potential chemotherapeutic agents for treating this disease.
Collapse
|
13
|
Sakyi PO, Broni E, Amewu RK, Miller WA, Wilson MD, Kwofie SK. Homology Modeling, de Novo Design of Ligands, and Molecular Docking Identify Potential Inhibitors of Leishmania donovani 24-Sterol Methyltransferase. Front Cell Infect Microbiol 2022; 12:859981. [PMID: 35719359 PMCID: PMC9201040 DOI: 10.3389/fcimb.2022.859981] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2022] [Accepted: 04/28/2022] [Indexed: 11/13/2022] Open
Abstract
The therapeutic challenges pertaining to leishmaniasis due to reported chemoresistance and toxicity necessitate the need to explore novel pathways to identify plausible inhibitory molecules. Leishmania donovani 24-sterol methyltransferase (LdSMT) is vital for the synthesis of ergosterols, the main constituents of Leishmania cellular membranes. So far, mammals have not been shown to possess SMT or ergosterols, making the pathway a prime candidate for drug discovery. The structural model of LdSMT was elucidated using homology modeling to identify potential novel 24-SMT inhibitors via virtual screening, scaffold hopping, and de-novo fragment-based design. Altogether, six potential novel inhibitors were identified with binding energies ranging from −7.0 to −8.4 kcal/mol with e-LEA3D using 22,26-azasterol and S1–S4 obtained from scaffold hopping via the ChEMBL, DrugBank, PubChem, ChemSpider, and ZINC15 databases. These ligands showed comparable binding energy to 22,26-azasterol (−7.6 kcal/mol), the main inhibitor of LdSMT. Moreover, all the compounds had plausible ligand efficiency-dependent lipophilicity (LELP) scores above 3. The binding mechanism identified Tyr92 to be critical for binding, and this was corroborated via molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The ligand A1 was predicted to possess antileishmanial properties with a probability of activity (Pa) of 0.362 and a probability of inactivity (Pi) of 0.066, while A5 and A6 possessed dermatological properties with Pa values of 0.205 and 0.249 and Pi values of 0.162 and 0.120, respectively. Structural similarity search via DrugBank identified vabicaserin, daledalin, zanapezil, imipramine, and cefradine with antileishmanial properties suggesting that the de-novo compounds could be explored as potential antileishmanial agents.
Collapse
Affiliation(s)
- Patrick O. Sakyi
- Department of Chemistry, School of Physical and Mathematical Sciences, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
- Department of Chemical Sciences, School of Sciences, University of Energy and Natural Resources, Sunyani, Ghana
| | - Emmanuel Broni
- Department of Biomedical Engineering, School of Engineering Sciences, College of Basic & Applied Sciences, University of Ghana, Accra, Ghana
- Department of Parasitology, Noguchi Memorial Institute for Medical Research (NMIMR), College of Health Sciences (CHS), University of Ghana, Accra, Ghana
| | - Richard K. Amewu
- Department of Chemistry, School of Physical and Mathematical Sciences, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
| | - Whelton A. Miller
- Department of Medicine, Loyola University Medical Center, Maywood, IL, United States
- Department of Molecular Pharmacology and Neuroscience, Loyola University Medical Center, Maywood, IL, United States
- Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA, United States
| | - Michael D. Wilson
- Department of Parasitology, Noguchi Memorial Institute for Medical Research (NMIMR), College of Health Sciences (CHS), University of Ghana, Accra, Ghana
- Department of Medicine, Loyola University Medical Center, Maywood, IL, United States
| | - Samuel Kojo Kwofie
- Department of Biomedical Engineering, School of Engineering Sciences, College of Basic & Applied Sciences, University of Ghana, Accra, Ghana
- Department of Biochemistry, Cell and Molecular Biology, West African Centre for Cell Biology of Infectious Pathogens, College of Basic and Applied Sciences, University of Ghana, Accra, Ghana
- *Correspondence: Samuel Kojo Kwofie,
| |
Collapse
|
14
|
Avery C, Baker L, Jacobs DJ. Functional Dynamics of Substrate Recognition in TEM Beta-Lactamase. ENTROPY 2022; 24:e24050729. [PMID: 35626612 PMCID: PMC9140794 DOI: 10.3390/e24050729] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/05/2022] [Revised: 05/16/2022] [Accepted: 05/18/2022] [Indexed: 02/01/2023]
Abstract
The beta-lactamase enzyme provides effective resistance to beta-lactam antibiotics due to substrate recognition controlled by point mutations. Recently, extended-spectrum and inhibitor-resistant mutants have become a global health problem. Here, the functional dynamics that control substrate recognition in TEM beta-lactamase are investigated using all-atom molecular dynamics simulations. Comparisons are made between wild-type TEM-1 and TEM-2 and the extended-spectrum mutants TEM-10 and TEM-52, both in apo form and in complex with four different antibiotics (ampicillin, amoxicillin, cefotaxime and ceftazidime). Dynamic allostery is predicted based on a quasi-harmonic normal mode analysis using a perturbation scan. An allosteric mechanism known to inhibit enzymatic function in TEM beta-lactamase is identified, along with other allosteric binding targets. Mechanisms for substrate recognition are elucidated using multivariate comparative analysis of molecular dynamics trajectories to identify changes in dynamics resulting from point mutations and ligand binding, and the conserved dynamics, which are functionally important, are extracted as well. The results suggest that the H10-H11 loop (residues 214-221) is a secondary anchor for larger extended spectrum ligands, while the H9-H10 loop (residues 194-202) is distal from the active site and stabilizes the protein against structural changes. These secondary non-catalytically-active loops offer attractive targets for novel noncompetitive inhibitors of TEM beta-lactamase.
Collapse
Affiliation(s)
- Chris Avery
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28223, USA; (C.A.); (L.B.)
| | - Lonnie Baker
- Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28223, USA; (C.A.); (L.B.)
| | - Donald J. Jacobs
- Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28223, USA
- Correspondence:
| |
Collapse
|
15
|
JEDi: java essential dynamics inspector - a molecular trajectory analysis toolkit. BMC Bioinformatics 2021; 22:226. [PMID: 33932974 PMCID: PMC8088583 DOI: 10.1186/s12859-021-04140-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Accepted: 04/19/2021] [Indexed: 11/20/2022] Open
Abstract
Background Principal component analysis (PCA) is commonly applied to the atomic trajectories of biopolymers to extract essential dynamics that describe biologically relevant motions. Although application of PCA is straightforward, specialized software to facilitate workflows and analysis of molecular dynamics simulation data to fully harness the power of PCA is lacking. The Java Essential Dynamics inspector (JEDi) software is a major upgrade from the previous JED software. Results Employing multi-threading, JEDi features a user-friendly interface to control rapid workflows for interrogating conformational motions of biopolymers at various spatial resolutions and within subregions, including multiple chain proteins. JEDi has options for Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) to construct covariance (Q), correlation (R), and partial correlation (P) matrices. Shrinkage and outlier thresholding are implemented for the accurate estimation of covariance. The effect of rare events is quantified using outlier and inlier filters. Applying sparsity thresholds in statistical models identifies latent correlated motions. Within a hierarchical approach, small-scale atomic motion is first calculated with a separate local cPCA calculation per residue to obtain eigenresidues. Then PCA on the eigenresidues yields rapid and accurate description of large-scale motions. Local cPCA on all residue pairs creates a map of all residue-residue dynamical couplings. Additionally, kernel PCA is implemented. JEDi output gives high quality PNG images by default, with options for text files that include aligned coordinates, several metrics that quantify mobility, PCA modes with their eigenvalues, and displacement vector projections onto the top principal modes. JEDi provides PyMol scripts together with PDB files to visualize individual cPCA modes and the essential dynamics occurring within user-selected time scales. Subspace comparisons performed on the most relevant eigenvectors using several statistical metrics quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes are available for both cPCA and dpPCA. Conclusion JEDi is a convenient toolkit that applies best practices in multivariate statistics for comparative studies on the essential dynamics of similar biopolymers. JEDi helps identify functional mechanisms through many integrated tools and visual aids for inspecting and quantifying similarity/differences in mobility and dynamic correlations. Supplementary Information The online version contains supplementary material available at 10.1186/s12859-021-04140-5.
Collapse
|
16
|
Grear T, Avery C, Patterson J, Jacobs DJ. Molecular function recognition by supervised projection pursuit machine learning. Sci Rep 2021; 11:4247. [PMID: 33608593 PMCID: PMC7895977 DOI: 10.1038/s41598-021-83269-y] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2020] [Accepted: 01/28/2021] [Indexed: 01/31/2023] Open
Abstract
Identifying mechanisms that control molecular function is a significant challenge in pharmaceutical science and molecular engineering. Here, we present a novel projection pursuit recurrent neural network to identify functional mechanisms in the context of iterative supervised machine learning for discovery-based design optimization. Molecular function recognition is achieved by pairing experiments that categorize systems with digital twin molecular dynamics simulations to generate working hypotheses. Feature extraction decomposes emergent properties of a system into a complete set of basis vectors. Feature selection requires signal-to-noise, statistical significance, and clustering quality to concurrently surpass acceptance levels. Formulated as a multivariate description of differences and similarities between systems, the data-driven working hypothesis is refined by analyzing new systems prioritized by a discovery-likelihood. Utility and generality are demonstrated on several benchmarks, including the elucidation of antibiotic resistance in TEM-52 beta-lactamase. The software is freely available, enabling turnkey analysis of massive data streams found in computational biology and material science.
Collapse
Affiliation(s)
- Tyler Grear
- grid.266859.60000 0000 8598 2218Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262 USA
| | - Chris Avery
- grid.266859.60000 0000 8598 2218Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262 USA ,grid.266859.60000 0000 8598 2218Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28262 USA
| | - John Patterson
- grid.266859.60000 0000 8598 2218Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28262 USA
| | - Donald J. Jacobs
- grid.266859.60000 0000 8598 2218Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28262 USA ,grid.266859.60000 0000 8598 2218Center for Biomedical Engineering and Science, University of North Carolina at Charlotte, Charlotte, NC 28262 USA
| |
Collapse
|
17
|
Shruti SR, Rajasekaran R. Identification of therapeutic peptide scaffold from tritrpticin family for urinary tract infections using in silico techniques. J Biomol Struct Dyn 2019; 38:4407-4417. [DOI: 10.1080/07391102.2019.1680437] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- S. R. Shruti
- Department of Biotechnology, School of Biosciences and Technology, VIT (Deemed to Be University), Vellore, India
| | - R. Rajasekaran
- Department of Biotechnology, School of Biosciences and Technology, VIT (Deemed to Be University), Vellore, India
| |
Collapse
|
18
|
Wang F, Shen L, Zhou H, Wang S, Wang X, Tao P. Machine Learning Classification Model for Functional Binding Modes of TEM-1 β-Lactamase. Front Mol Biosci 2019; 6:47. [PMID: 31355207 PMCID: PMC6629954 DOI: 10.3389/fmolb.2019.00047] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Accepted: 06/11/2019] [Indexed: 11/13/2022] Open
Abstract
TEM family of enzymes is one of the most commonly encountered β-lactamases groups with different catalytic capabilities against various antibiotics. Despite the studies investigating the catalytic mechanism of TEM β-lactamases, the binding modes of these enzymes against ligands in different functional catalytic states have been largely overlooked. But the binding modes may play a critical role in the function and even the evolution of these proteins. In this work, a newly developed machine learning analysis approach to the recognition of protein dynamics states was applied to compare the binding modes of TEM-1 β-lactamase with regard to penicillin in different catalytic states. While conventional analysis methods, including principal components analysis (PCA), could not differentiate TEM-1 in different binding modes, the application of a machine learning method led to excellent classification models differentiating these states. It was also revealed that both reactant/product states and apo/product states are more differentiable than the apo/reactant states. The feature importance generated by the training procedure of the machine learning model was utilized to evaluate the contribution from residues at active sites and in different secondary structures. Key active site residues, Ser70 and Ser130, play a critical role in differentiating reactant/product states, while other active site residues are more important for differentiating apo/product states. Overall, this study provides new insights into the different dynamical function states of TEM-1 and may open a new venue for β-lactamases functional and evolutional studies in general.
Collapse
Affiliation(s)
- Feng Wang
- Department of Chemistry, Center for Scientific Computation, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX, United States
| | - Li Shen
- Department of Chemistry, Center for Scientific Computation, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX, United States
| | - Hongyu Zhou
- Department of Chemistry, Center for Scientific Computation, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX, United States
| | - Shouyi Wang
- Department of Industrial, Manufacturing, and Systems Engineering, University of Texas at Arlington, Arlington, TX, United States
| | - Xinlei Wang
- Department of Statistical Science, Southern Methodist University, Dallas, TX, United States
| | - Peng Tao
- Department of Chemistry, Center for Scientific Computation, Center for Drug Discovery, Design, and Delivery (CD4), Southern Methodist University, Dallas, TX, United States
| |
Collapse
|
19
|
Brown BP, Zhang YK, Westover D, Yan Y, Qiao H, Huang V, Du Z, Smith JA, Ross JS, Miller VA, Ali S, Bazhenova L, Schrock AB, Meiler J, Lovly CM. On-target Resistance to the Mutant-Selective EGFR Inhibitor Osimertinib Can Develop in an Allele-Specific Manner Dependent on the Original EGFR-Activating Mutation. Clin Cancer Res 2019; 25:3341-3351. [PMID: 30796031 PMCID: PMC6548651 DOI: 10.1158/1078-0432.ccr-18-3829] [Citation(s) in RCA: 68] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2018] [Revised: 01/18/2019] [Accepted: 02/14/2019] [Indexed: 01/10/2023]
Abstract
PURPOSE The third-generation EGFR inhibitor, osimertinib, is the first mutant-selective inhibitor that has received regulatory approval for the treatment of patients with EGFR-mutant lung cancer. Despite the development of highly selective third-generation inhibitors, acquired resistance remains a significant clinical challenge. Recently, we and others have identified a novel osimertinib resistance mutation, G724S, which was not predicted in in vitro screens. Here, we investigate how G724S confers resistance to osimertinib.Experimental Design: We combine structure-based predictive modeling of G724S in combination with the 2 most common EGFR-activating mutations, exon 19 deletion (Ex19Del) and L858R, with in vitro drug-response models and patient genomic profiling. RESULTS Our simulations suggest that the G724S mutation selectively reduces osimertinib-binding affinity in the context of Ex19Del. Consistent with our simulations, cell lines transduced with Ex19Del/G724S demonstrate resistance to osimertinib, whereas cells transduced with L858R/G724S are sensitive to osimertinib. Subsequent clinical genomic profiling data further suggest G724S occurs with Ex19Del but not L858R. Furthermore, we demonstrate that Ex19Del/G724S retains sensitivity to afatinib, but not to erlotinib, suggesting a possible therapy for patients at the time of disease relapse. CONCLUSIONS Altogether, these data suggest that G724S is an allele-specific resistance mutation emerging in the context of Ex19Del but not L858R. Our results fundamentally reframe the problem of targeted therapy resistance from one focused on the "drug-resistance mutation" pair to one focused on the "activating mutation-drug-resistance mutation" trio. This has broad implications across clinical oncology.
Collapse
Affiliation(s)
- Benjamin P. Brown
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA.,Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Yun-Kai Zhang
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - David Westover
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Yingjun Yan
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Huan Qiao
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Vincent Huang
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Zhenfang Du
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA
| | - Jarrod A. Smith
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA.,Department of Biochemistry, Vanderbilt University, Nashville, TN, USA
| | | | | | - Siraj Ali
- Foundation Medicine, Inc., Cambridge, MA, USA
| | | | | | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, Tennessee. .,Center for Structural Biology, Vanderbilt University, Nashville, Tennessee.,Vanderbilt Ingram Cancer Center, Vanderbilt University Medical Center, Nashville, Tennessee
| | - Christine M. Lovly
- Department of Medicine - Division of Hematology and Oncology, Vanderbilt University Medical Center, Nashville, TN 37232, USA.,Vanderbilt Ingram Cancer Center, Vanderbilt University Medical Center, Nashville, TN 37232, USA., Co-corresponding authors: Christine M. Lovly, MD, PhD, 2220 Pierce Avenue, 659 Preston Research Building, Nashville, TN 37232-6307, Phone: 615-936-3457, and Jens Meiler, PhD, Center for Structural Biology, 465 21st Ave South, BIOSCI/MRBIII, Room 5144B, Nashville, TN 37232-8725 USA, Phone: 615 936 5662,
| |
Collapse
|
20
|
High throughput nonparametric probability density estimation. PLoS One 2018; 13:e0196937. [PMID: 29750803 PMCID: PMC5947915 DOI: 10.1371/journal.pone.0196937] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2017] [Accepted: 04/23/2018] [Indexed: 12/04/2022] Open
Abstract
In high throughput applications, such as those found in bioinformatics and finance, it is important to determine accurate probability distribution functions despite only minimal information about data characteristics, and without using human subjectivity. Such an automated process for univariate data is implemented to achieve this goal by merging the maximum entropy method with single order statistics and maximum likelihood. The only required properties of the random variables are that they are continuous and that they are, or can be approximated as, independent and identically distributed. A quasi-log-likelihood function based on single order statistics for sampled uniform random data is used to empirically construct a sample size invariant universal scoring function. Then a probability density estimate is determined by iteratively improving trial cumulative distribution functions, where better estimates are quantified by the scoring function that identifies atypical fluctuations. This criterion resists under and over fitting data as an alternative to employing the Bayesian or Akaike information criterion. Multiple estimates for the probability density reflect uncertainties due to statistical fluctuations in random samples. Scaled quantile residual plots are also introduced as an effective diagnostic to visualize the quality of the estimated probability densities. Benchmark tests show that estimates for the probability density function (PDF) converge to the true PDF as sample size increases on particularly difficult test probability densities that include cases with discontinuities, multi-resolution scales, heavy tails, and singularities. These results indicate the method has general applicability for high throughput statistical inference.
Collapse
|
21
|
Molecular Dynamics vs. Stochastic Processes: Are We Heading Anywhere? ENTROPY 2018; 20:e20050348. [PMID: 33265438 PMCID: PMC7512867 DOI: 10.3390/e20050348] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 05/02/2018] [Revised: 05/02/2018] [Accepted: 05/02/2018] [Indexed: 01/18/2023]
|