1
|
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] [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
|
2
|
Singh A, Varadarajan A, Pant P, Singh TP, Vikram NK, Sharma S, Sharma P. Identification of potential anti-mucor agents by targeting endothelial cell receptor glucose-regulated protein-78 using in silico approach. J Biomol Struct Dyn 2024; 42:4344-4355. [PMID: 37288794 DOI: 10.1080/07391102.2023.2220809] [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/22/2022] [Accepted: 05/28/2023] [Indexed: 06/09/2023]
Abstract
Mucormycosis is a fungal infection of the sinuses, brain and lungs that is the cause of approximately 50% mortality rate despite the available first-line therapy. Glucose-Regulated Protein 78 (GRP78) is already reported to be a novel host receptor that mediates invasion and damage of human endothelial cells by Rhizopus oryzae and Rhizopus delemar, the most common etiologic species of Mucorales. The expression of GRP78 is also regulated by the levels of iron and glucose in the blood. There are several antifungal drugs in the market but they pose a serious side effect to the vital organs of the body. Therefore, there is an immediate need to discover effective drug molecules having increased efficacy with no side effects. With the help of various computational tools, the current study was attempted to determine potential antimucor agents against GRP78. The receptor molecule GRP78 was screened against 8820 known drugs deposited in DrugBank library using high-throughput virtual screening method. Total top 10 compounds were selected based on the binding energies greater than the reference co-crystal molecule. Furthermore, molecular dynamic (MD) simulations using AMBER were performed to calculate the stability of the top-ranked compounds in the active site of GRP78. After extensive computational studies, we propose that two compounds (CID439153 and CID5289104) have inhibitory potency against mucormycosis and can serve as potential drugs that can form the basis of treating mucormycosis disease.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Anamika Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Ashwin Varadarajan
- Department of Medicine, All India Institute of Medical Sciences, New Delhi, India
| | - Pradeep Pant
- Department of Chemistry, Indian Institute of Technology, New Delhi, India
| | - Tej P Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Naval K Vikram
- Department of Medicine, All India Institute of Medical Sciences, New Delhi, India
| | - Sujata Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Pradeep Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| |
Collapse
|
3
|
Wang Y, Wei S. Influence of hydrogen bonds on the reaction of guanine and hydroxyl radical: DFT calculations in C(H +)GC motif. Phys Chem Chem Phys 2024; 26:5683-5692. [PMID: 38288746 DOI: 10.1039/d3cp05885a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
A comprehensive theoretical investigation was performed to illuminate the influence of hydrogen bonds (H-bonds) on the obscure reaction of a hydroxyl radical (HO˙) and guanine (G) by selecting the building block of parallel triplex DNA, C(H+)GC, as the model. By mapping the energy profiles for addition and hydrogen abstraction reactions, the favorable pathway is predicted. The results reveal that in the C(H+)GC context, barrierless hydrogen abstraction from N2 of G leading to a neutral radical G(N2-H)˙ appears to become significant, but electrophilic attack by HO˙ on C8 of G resulting in 8-oxoG is the most thermodynamically favorable course. This shows a strong structural dependence due to the context constrained by the H-bond, which is dramatically different from the situation in unencumbered G. More interestingly, it proves that the stability order of resulting adduct radicals is not altered by H-bonding, but the activity for possible sites of the hydroxylation reaction changes. The significant influence of the H-bond on elementary reactions involved in the reaction is emphasized in the C(H+)GC context but is not restricted to the H-abstraction reaction. It is greatly anticipated that the present study could provide thoughtful insights into the vague hydroxyl radical-induced oxidation chemistry.
Collapse
Affiliation(s)
- Yinghui Wang
- College of Science, Chang'an University, Xi'an 710064, China.
| | - Simin Wei
- State Key Laboratory of Research & Development of Characteristic Qin Medicine Resources (Cultivation), Co-Construction Collaborative Innovation Center for Chinese Medicine Resources Industrialization by Shaanxi & Education Ministry, Shaanxi University of Chinese Medicine, Xianyang 712083, China.
| |
Collapse
|
4
|
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
|
5
|
Xie Z, Zhang S, Wu Y, Liang J, Yao W, Qu R, Tong X, Zhang G, Yang H. Interaction of isoquinoline alkaloids with pyrimidine motif triplex DNA by mass spectrometry and spectroscopies reveals diverse mechanisms. Heliyon 2023; 9:e14954. [PMID: 37082631 PMCID: PMC10112036 DOI: 10.1016/j.heliyon.2023.e14954] [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: 07/21/2022] [Revised: 03/21/2023] [Accepted: 03/22/2023] [Indexed: 04/22/2023] Open
Abstract
Isoquinoline alkaloids represent an important class of molecules due to their broad range of pharmacology and clinical utility. Prospective development and use of these alkaloids as effective anticancer agents have elicited great interest. In this study, in order to reveal structure-activity relationship, we present the characterization of bioactive isoquinoline alkaloid-DNA triplex interactions, with particular emphasis on the sequence selectivity and preference of binding to the two types of DNA triplexes, by electrospray ionization mass spectrometry (ESI-MS) and various spectroscopic techniques. The six alkaloids, including coptisine, columbamine, epiberberine, berberrubine, jateorhizine, and fangchinoline, were selected to explore their interactions with the TC and TTT triplex DNA structures. Berberrubine, fangchinoline, coptisine, columbamine, and epiberberine have preference for TC rich DNA sequences compared to TTT rich DNA triplex based on affinity values in MS. The experimental results from different fragmentation modes in tandem MS, subtractive and hyperchromic effects in UV absorption spectra, fluorescence quenching and enhancement in fluorescence spectra, and strong conformational changes in circular dichroism (CD) hinted that the interaction between isoquinoline alkaloid-TC/TTT DNA had diverse mechanisms including at least two different binding modes: the electrostatic binding and the intercalation binding. Interestingly, columbamine, berberrubine, and fangchinoline can stabilize TTT triplex as inferred from optical thermal melting profiles, while it was not the case in TC triplex. These results provide new insights into binding of isoquinoline alkaloids to pyrimidine motif triplex DNA.
Collapse
Affiliation(s)
- Zhaoyang Xie
- Northeast Asia Institute of Chinese Medicine, Changchun University of Chinese Medicine, Changchun, 130021, China
| | - Sunuo Zhang
- School of Materials Science and Engineering, Nankai University, Tianjin, 300350, China
| | - Yi Wu
- Northeast Asia Institute of Chinese Medicine, Changchun University of Chinese Medicine, Changchun, 130021, China
| | - Jinling Liang
- College of Pharmacy, Changchun University of Chinese Medicine, Changchun, 130117, China
| | - Wenbin Yao
- College of Pharmacy, Changchun University of Chinese Medicine, Changchun, 130117, China
| | - Ruoning Qu
- China-Japan Union Hospital of Jilin University, Changchun, 130033, China
| | - Xiaole Tong
- Jilin Jian Yisheng Pharmaceutical Co., Ltd., Jian, 134200, China
| | - Guang Zhang
- China-Japan Union Hospital of Jilin University, Changchun, 130033, China
- Corresponding author.
| | - Hongmei Yang
- Northeast Asia Institute of Chinese Medicine, Changchun University of Chinese Medicine, Changchun, 130021, China
- Corresponding author.
| |
Collapse
|
6
|
Rao GN, Jupudi S, Pant P, Palathoti N, Rajagopal K, Govindasamy R, Justin A. Ceftriaxone induces glial EAAT-2 promotor region via NF-kB conformational changes: An interaction analysis using HADDOCK. J Cell Biochem 2023; 124:359-372. [PMID: 36649127 DOI: 10.1002/jcb.30370] [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: 09/28/2022] [Revised: 12/30/2022] [Accepted: 01/03/2023] [Indexed: 01/18/2023]
Abstract
Excitotoxicity, depletion of energy metabolites, and ionic imbalance are the major factors involved in neurodegeneration mediated through excitatory amino acid transporter-2 (EAAT-2) dysfunction in ischemic insult. Recent studies have revealed that ceftriaxone expresses EAAT-2 via nuclear transcription factor kappa-B (NF-kB) signaling pathway, stimulation of EAAT-2 expression in the ischemic, and excitotoxic conditions that could provide potential benefits to control neurodegeneration. In this study, we have predicted the in silico model for interaction between NF-kB and EAAT-2 promoter region to rule out the conformational changes for the expression of EAAT-2 protein. Using homology-built model of NF-kB, we identified ceftriaxone-induced conformational changes in gene locus -272 of DNA where NF-kB binding with EAAT-2 promoter region through protein-DNA docking calculation. The interaction profile and conformational dynamics occurred between ceftriaxone predocked and postdocked conformations of NF-kB with DNA employing HADDOCK 2.2 web server followed by 250 ns long all atom explicit solvent molecular dynamics simulations. Both the protein and DNA exhibited modest conformational changes with respect to HADDOCK score, energy terms (desolvation energy [Edesolv ]), van der waal energy (Evdw ), electrostatic energy (Eelec ), restraints energy (Eair ), buried surface area, root mean square deviation, RMSF, radius of gyration, total hydrogen bonds when ceftriaxone pre- and postdocked NF-kB conformations were bound to DNA. Hence, the conformational changes in the C-terminal domain could be the reason for EAAT-2 expression through ceftriaxone specific binding pocket of -272 of DNA.
Collapse
Affiliation(s)
- Gaddam Narasimha Rao
- Department of Pharmacology, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Ooty, Tamilnadu, India
| | - Srikanth Jupudi
- Department of Pharmaceutical Chemistry, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Ooty, Tamilnadu, India
| | - Pradeep Pant
- Department of Chemistry, Indian Institute of Technology, HauzKhas, New Delhi, India
| | - Nagarjuna Palathoti
- Department of Pharmaceutical Chemistry, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Ooty, Tamilnadu, India
| | - Kalirajan Rajagopal
- Department of Pharmaceutical Chemistry, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Ooty, Tamilnadu, India
| | - Rathika Govindasamy
- Department of Chemistry, PSG College of Arts & Science, Coimbatore, Tamilnadu, India
| | - Antony Justin
- Department of Pharmacology, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Ooty, Tamilnadu, India
| |
Collapse
|
7
|
Gupta A, Vijayan V, Pant P, Kaur P, Singh TP, Sharma P, Sharma S. Structure prediction and discovery of inhibitors against phosphopantothenoyl cysteine synthetase of Acinetobacter baumannii. J Biomol Struct Dyn 2022; 40:11405-11417. [PMID: 34348086 DOI: 10.1080/07391102.2021.1958699] [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: 12/24/2022]
Abstract
Acinetobacter baumannii is an extremely dangerous multidrug-resistant (MDR) gram-negative pathogen which poses a serious life-threatening risk in immunocompromised patients. Phosphopantothenoyl cysteine synthetase (PPCS) catalyzes the formation of an amide bond between L-cysteine and phosphopantothenic acid (PPA) to form 4'- Phosphopantothenoylcysteine during Coenzyme A (CoA) biosynthesis. CoA is a crucial cofactor for cellular survival and inhibiting its synthesis will result in cell death. Bacterial PPCS differs from eukaryotic PPCS in a number of ways like it exists as a C-terminal domain of a PPCDC/PPCS fusion protein whereas eukaryotic PPCS exists as an independent protein. This difference makes it an attractive drug target. For which a conventional iterative approach of SBDD (structure-based drug design) was used, which began with three-dimensional structure prediction of AbPPCS using PHYRE 2.0. A database of FDA-approved compounds (Drug Bank) was then screened against the target of interest by means of docking score and glide energy, leading to the identification of 6 prominent drug candidates. The shortlisted 6 molecules were further subjected to all-atom MD simulation studies in explicit-solvent conditions (using AMBER force field). The MD simulation studies revealed that the ligands DB65103, DB449108 and DB443210, maintained several H-bonds with intense van der Waals contacts at the active site of the protein with high binding free energies: -11.42 kcal/mol, -10.49 kcal/mol and -10.98 kcal/mol, respectively, calculated via MM-PBSA method. Overall, binding of these compounds at the active site was found to be the most stable and robust highlighting the potential of these compounds to serve as antibacterials.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Akshita Gupta
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Viswanathan Vijayan
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi, India
| | - Punit Kaur
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Tej P Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Pradeep Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Sujata Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| |
Collapse
|
8
|
Pant P. Harmonizing Interstrand Electrostatic Repulsion by Conformational Rigidity in Counterion-Deprived Z-DNA: A Molecular Dynamics Study. J Phys Chem B 2022; 126:9956-9963. [PMID: 36412276 DOI: 10.1021/acs.jpcb.2c04527] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Abstract
Deoxyribonucleic acid (DNA) is a vital biomacromolecule. Although the right-handed B-DNA type helical structure is the most abundant and extensively studied form of DNA, several noncanonical forms, such as triplex, quadruplex, Z-DNA, A-DNA, and ss-DNA, have been probed from time to time to gain insights into the DNA's function. Z-DNA was recently found to be involved in cancer and several autoimmune diseases. In the present Article, we evaluated the conformational stability of locked-sugar-based Z-DNA via all-atom explicit-solvent molecular dynamics simulations and found that the modified DNA maintained the left-handed conformation even in the absence of counterions, wherein the structural rigidity dominates over the electrostatic repulsion between the complementary strands. The control Z-DNA without counterions, as expected, instantaneously resulted in unfolded states. The remarkable stability of the conformationally locked model system was thoroughly investigated via structural and energetic perspectives and was probably the result of the backbone widening in tandem with enhanced electrostatics between complementary strands. We believe that the design of the proposed modified Z-DNA construct could help understand the otherwise delicate Z-DNA conformation even in salt-deprived conditions. The design could also motivate the medicinal use of short segments of such modified nucleotides and could be utilized in more advanced modeling techniques, such as DNA origami which has gained popularity in recent years.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, New Delhi 110016, India
| |
Collapse
|
9
|
Pant P, Pathak A, Jayaram B. Bicyclo-DNA mimics with enhanced protein binding affinities: insights from molecular dynamics simulations. J Biomol Struct Dyn 2022; 41:4040-4047. [PMID: 35403569 DOI: 10.1080/07391102.2022.2061594] [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
DNA-protein interactions occur at all levels of DNA expression and replication and are crucial determinants for the survival of a cell. Several modified nucleotides have been utilized to manipulate these interactions and have implications in drug discovery. In the present article, we evaluated the binding of bicyclo-nucleotides (generated by forming a methylene bridge between C1' and C5' in sugar, leading to a bicyclo system with C2' axis of symmetry at the nucleotide level) to proteins. We utilized four ssDNA-protein complexes with experimentally known binding free energies and investigated the binding of modified nucleotides to proteins via all-atom explicit solvent molecular dynamics (MD) simulations (200 ns), and compared the binding with control ssDNA-protein systems. The modified ssDNA displayed enhanced binding to proteins as compared to the control ssDNA, as seen by means of MD simulations followed by MM-PBSA calculations. Further, the Delphi-based electrostatic estimation revealed that the high binding of modified ssDNA to protein might be related to the enhanced electrostatic complementarity displayed by the modified ssDNA molecules in all the four systems considered for the study. The improved binding achieved with modified nucleotides can be utilized to design and develop anticancer/antisense molecules capable of targeting proteins or ssRNAs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi, India
| | - Amita Pathak
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi, India
| | - B Jayaram
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi, India.,Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi, India
| |
Collapse
|
10
|
Ahmad N, Singh A, Gupta A, Pant P, Singh TP, Sharma S, Sharma P. Discovery of the Lead Molecules Targeting the First Step of the Histidine Biosynthesis Pathway of Acinetobacter baumannii. J Chem Inf Model 2022; 62:1744-1759. [PMID: 35333517 DOI: 10.1021/acs.jcim.1c01421] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Acinetobacter baumannii is a multidrug-resistant, opportunistic, nosocomial pathogen for which a new line of treatments is desperately needed. We have targeted the enzyme of the first step of the histidine biosynthesis pathway, viz., ATP-phosphoribosyltransferase (ATP-PRT). The three-dimensional structure of ATP-PRT was predicted on the template of the known three-dimensional structure of ATP-PRT from Psychrobacter arcticus (PaATPPRT) using a homology modeling approach. High-throughput virtual screening (HTVS) of the antibacterial library of Life Chemicals Inc., Ontario, Canada was carried out followed by molecular dynamics simulations of the top hit compounds. In silico results were then biochemically validated using surface plasmon resonance spectroscopy. We found that two compounds, namely, F0843-0019 and F0608-0626, were binding with micromolar affinities to the ATP-phosphoribosyltransferase from Acinetobacter baumannii (AbATPPRT). Both of these compounds were binding in the same way as AMP in PaATPPRT, and the important residues of the active site, viz., Val4, Ser72, Thr76, Tyr77, Glu95, Lys134, Val136, and Tyr156, were also interacting via hydrogen bonds. The calculated binding energies of these compounds were -10.5 kcal/mol and -11.1 kcal/mol, respectively. These two compounds can be used as the potential lead molecules for designing antibacterial compounds in the future, and this information will help in drug discovery programs against Acinetobacter worldwide.
Collapse
Affiliation(s)
- Nabeel Ahmad
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| | - Anamika Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| | - Akshita Gupta
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| | - Pradeep Pant
- Department of Chemistry, Indian Institute of Technology, Delhi 110016, India
| | - Tej P Singh
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| | - Sujata Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| | - Pradeep Sharma
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi 110029, India
| |
Collapse
|
11
|
Pant P, Aggarwal L. Assessing the DNA structural integrity via selective annihilation of Watson-Crick hydrogen bonds: Insights from molecular dynamics simulations. Biophys Chem 2022; 282:106758. [DOI: 10.1016/j.bpc.2021.106758] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 12/09/2021] [Accepted: 12/30/2021] [Indexed: 01/17/2023]
|
12
|
Wang J, Arantes PR, Bhattarai A, Hsu RV, Pawnikar S, Huang YMM, Palermo G, Miao Y. Gaussian accelerated molecular dynamics (GaMD): principles and applications. WILEY INTERDISCIPLINARY REVIEWS. COMPUTATIONAL MOLECULAR SCIENCE 2021; 11:e1521. [PMID: 34899998 PMCID: PMC8658739 DOI: 10.1002/wcms.1521] [Citation(s) in RCA: 92] [Impact Index Per Article: 30.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2020] [Accepted: 01/28/2021] [Indexed: 12/20/2022]
Abstract
Gaussian accelerated molecular dynamics (GaMD) is a robust computational method for simultaneous unconstrained enhanced sampling and free energy calculations of biomolecules. It works by adding a harmonic boost potential to smooth biomolecular potential energy surface and reduce energy barriers. GaMD greatly accelerates biomolecular simulations by orders of magnitude. Without the need to set predefined reaction coordinates or collective variables, GaMD provides unconstrained enhanced sampling and is advantageous for simulating complex biological processes. The GaMD boost potential exhibits a Gaussian distribution, thereby allowing for energetic reweighting via cumulant expansion to the second order (i.e., "Gaussian approximation"). This leads to accurate reconstruction of free energy landscapes of biomolecules. Hybrid schemes with other enhanced sampling methods, such as the replica exchange GaMD (rex-GaMD) and replica exchange umbrella sampling GaMD (GaREUS), have also been introduced, further improving sampling and free energy calculations. Recently, new "selective GaMD" algorithms including the ligand GaMD (LiGaMD) and peptide GaMD (Pep-GaMD) enabled microsecond simulations to capture repetitive dissociation and binding of small-molecule ligands and highly flexible peptides. The simulations then allowed highly efficient quantitative characterization of the ligand/peptide binding thermodynamics and kinetics. Taken together, GaMD and its innovative variants are applicable to simulate a wide variety of biomolecular dynamics, including protein folding, conformational changes and allostery, ligand binding, peptide binding, protein-protein/nucleic acid/carbohydrate interactions, and carbohydrate/nucleic acid interactions. In this review, we present principles of the GaMD algorithms and recent applications in biomolecular simulations and drug design.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Pablo R Arantes
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Apurba Bhattarai
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr, Lawrence, KS, 66047, United States
| | - Rohaine V Hsu
- Department of Bioengineering, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Shristi Pawnikar
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, KS, 66047, United States
| | - Yu-Ming M Huang
- Department of Physics & Astronomy, Wayne State University, 666 W Hancock St, Detroit, MI 48207, USA
| | - Giulia Palermo
- Department of Bioengineering and Department of Chemistry, University of California Riverside, 900 University Avenue, Riverside, CA 92512, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, 2030 Becker Dr., Lawrence, Kansas 66047, United States
| |
Collapse
|
13
|
Pant P, Pathak A, Jayaram B. Symmetric Nucleosides as Potent Purine Nucleoside Phosphorylase Inhibitors. J Phys Chem B 2021; 125:2856-2862. [PMID: 33715357 DOI: 10.1021/acs.jpcb.0c10553] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Nucleic acids are one of the most enigmatic biomolecules crucial to several biological processes. Nucleic acid-protein interactions are vital for the coordinated and controlled functioning of a cell, leading to the design of several nucleoside/nucleotide analogues capable of mimicking these interactions and hold paramount importance in the field of drug discovery. Purine nucleoside phosphorylase is a well-established drug target due to its association with numerous immunodeficiency diseases. Here, we study the binding of human purine nucleoside phosphorylase (PNP) to some bidirectional symmetric nucleosides, a class of nucleoside analogues that are more flexible due to the absence of sugar pucker restraints. We compared the binding energies of PNP-symmetric nucleosides to the binding energies of PNP-inosine/Imm-H (a transition-state analogue), by means of 200 ns long all-atom explicit-solvent Gaussian accelerated molecular dynamics simulations followed by energetics estimation using the MM-PBSA methodology. Quite interestingly, we observed that a few symmetric nucleosides, namely, ν3 and ν4, showed strong binding with PNP (-14.1 and -12.6 kcal/mol, respectively), higher than inosine (-6.3 kcal/mol) and Imm-H (-9.6 kcal/mol). This is rationalized by an enhanced hydrogen-bond network for symmetric nucleosides compared to inosine and Imm-H while maintaining similar van der Waals contacts. We note that the chemical structures of both ν3 and ν4, due to an additional unsaturation in them, resemble enzymatic transition states and fall in the category of transition-state analogues (TSAs), which are quite popular.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi 110016, India
| | - Amita Pathak
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi 110016, India
| | - B Jayaram
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.,Supercomputing Facility for Bioinformatics & Computational Biology, Hauz Khas, New Delhi 110016, India.,Kusuma School of Biological Sciences, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
14
|
Gupta A, Rani C, Pant P, Vijayan V, Vikram N, Kaur P, Singh TP, Sharma S, Sharma P. Structure-Based Virtual Screening and Biochemical Validation to Discover a Potential Inhibitor of the SARS-CoV-2 Main Protease. ACS OMEGA 2020; 5:33151-33161. [PMID: 33398250 PMCID: PMC7754785 DOI: 10.1021/acsomega.0c04808] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/01/2020] [Accepted: 12/09/2020] [Indexed: 05/27/2023]
Abstract
The recent pandemic caused by SARS-CoV-2 has led the world to a standstill, causing a medical and economic crisis worldwide. This crisis has triggered an urgent need to discover a possible treatment strategy against this novel virus using already-approved drugs. The main protease (Mpro) of this virus plays a critical role in cleaving the translated polypeptides that makes it a potential drug target against COVID-19. Taking advantage of the recently discovered three-dimensional structure of Mpro, we screened approved drugs from the Drug Bank to find a possible inhibitor against Mpro using computational methods and further validating them with biochemical studies. The docking and molecular dynamics study revealed that DB04983 (denufosol) showed the best glide docking score, -11.884 kcal/mol, and MM-PBSA binding free energy, -10.96 kcal/mol. Cobicistat, cangrelor (previous computational studies in our lab), and denufosol (current study) were tested for the in vitro inhibitory effects on Mpro. The IC50 values of these drugs were ∼6.7 μM, 0.9 mM, and 1.3 mM, respectively, while the values of dissociation constants calculated using surface plasmon resonance were ∼2.1 μM, 0.7 mM, and 1.4 mM, respectively. We found that cobicistat is the most efficient inhibitor of Mpro both in silico and in vitro. In conclusion, cobicistat, which is already an FDA-approved drug being used against HIV, may serve as a good inhibitor against the main protease of SARS-CoV-2 that, in turn, can help in combating COVID-19, and these results can also form the basis for the rational structure-based drug design against COVID-19.
Collapse
Affiliation(s)
- Akshita Gupta
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Chitra Rani
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Pradeep Pant
- Department
of Chemistry, Indian Institute of Technology
Delhi, New Delhi 110016, India
- Computational
Biochemistry, University of Duisburg Essen, Duisburg 47279, Germany
| | - Viswanathan Vijayan
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Naval Vikram
- Department
of Medicine, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Punit Kaur
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Tej Pal Singh
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Sujata Sharma
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| | - Pradeep Sharma
- Department
of Biophysics, All India Institute of Medical
Sciences, New Delhi 110029, India
| |
Collapse
|
15
|
Pant P, Fisher M. Marshall's nucleic acid: From double-helical structure to a potent intercalator. Biophys Chem 2020; 269:106525. [PMID: 33352335 DOI: 10.1016/j.bpc.2020.106525] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2020] [Revised: 12/09/2020] [Accepted: 12/09/2020] [Indexed: 02/06/2023]
Abstract
Deoxyribonucleic acid (DNA) not only stores genetic information but also emerged as a popular drug target. Modified nucleotides/nucleosides have been extensively studied in recent years wherein the sugar/nucleobase/phosphate-backbone has been altered. Several such molecules are FDA approved, capable of targeting nucleic acids and proteins. In this article, we modified negatively charged phosphate backbone to marshall's acid-based neutral backbone and analyzed the resultant structures by utilizing Gaussian accelerated molecular dynamics simulations (1 μs) in aqueous media at 150 mM salt concentration. We noted that the double-helical marshall's nucleic acid structure was partially denatured during the course of simulations, however, after using conformationally locked sugar, the marshall's nucleic acid (hereby called MNA) maintained the double-helical structure throughout the simulations. Despite the fact that MNA has a more extended backbone than the regular DNA, surprisingly, both showed similar helical rise (~3.4 Å) along with a comparable Watson-Crick hydrogen bond profile. The backbone difference was majorly compensated in terms of helical twist (~56° (MNA) and ~ 35° (control DNA)). Further, we examined a few MNA based ss-dinucleotides as intercalating ligands for a regular B-DNA. Quite strikingly, the ligands unwinded the DNA and showed intercalating properties with high DNA binding affinities. Hence, the use of small fragments of MNA based molecules in DNA targeted drug discovery is foreseen.
Collapse
Affiliation(s)
- Pradeep Pant
- Department of Chemistry, Indian Institute of Technology Delhi, India.
| | - Maria Fisher
- Department of Biosciences, University of Helsinki, Finland
| |
Collapse
|