Ligand selectivity in the recognition of protoberberine alkaloids by hybrid-2 human telomeric G-quadruplex: Binding free energy calculation, fluorescence binding, and NMR experiments

Nanjie Deng, Junchao Xia, Lauren Wickstrom, Clement Lin, Kaibo Wang, Peng He, Yunting Yin, Danzhou Yang

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

The human telomeric G-quadruplex (G4) is an attractive target for developing anticancer drugs. Natural products protoberberine alkaloids are known to bind human telomeric G4 and inhibit telomerase. Among several structurally similar protoberberine alkaloids, epiberberine (EPI) shows the greatest specificity in recognizing the human telomeric G4 over duplex DNA and other G4s. Recently, NMR study revealed that EPI recognizes specifically the hybrid-2 form human telomeric G4 by inducing large rearrangements in the 5 0 -flanking segment and loop regions to form a highly extensive four-layered binding pocket. Using the NMR structure of the EPI-human telomeric G4 complex, here we perform molecular dynamics free energy calculations to elucidate the ligand selectivity in the recognition of protoberberines by the human telomeric G4. The MM-PB(GB)SA (molecular mechanics-Poisson Boltzmann/Generalized Born) Surface Area) binding free energies calculated using the Amber force fields bsc0 and OL15 correlate well with the NMR titration and binding affinity measurements, with both calculations correctly identifying the EPI as the strongest binder to the hybrid-2 telomeric G4 wtTel26. The results demonstrated that accounting for the conformational flexibility of the DNA-ligand complexes is crucially important for explaining the ligand selectivity of the human telomeric G4. While the MD-simulated (molecular dynamics) structures of the G-quadruplex-alkaloid complexes help rationalize why the EPI-G4 interactions are optimal compared with the other protoberberines, structural deviations from the NMR structure near the binding site are observed in the MD simulations. We have also performed binding free energy calculation using the more rigorous double decoupling method (DDM); however, the results correlate less well with the experimental trend, likely due to the difficulty of adequately sampling the very large conformational reorganization in the G4 induced by the protoberberine binding.

Original languageEnglish (US)
Article number1574
JournalMolecules
DOIs
StatePublished - Apr 21 2019
Externally publishedYes

Fingerprint

G-Quadruplexes
alkaloids
Alkaloids
Free energy
Fluorescence
selectivity
free energy
Nuclear magnetic resonance
Ligands
fluorescence
nuclear magnetic resonance
ligands
Experiments
Molecular dynamics
deoxyribonucleic acid
Molecular Dynamics Simulation
molecular dynamics
Amber
Molecular mechanics
Telomerase

Keywords

  • Binding free energy
  • G-quadruplex
  • MM-PB(Gb)SA
  • Molecular dynamics simulation
  • NMR
  • NMR titration
  • Protoberberine

ASJC Scopus subject areas

  • Analytical Chemistry
  • Chemistry (miscellaneous)
  • Molecular Medicine
  • Pharmaceutical Science
  • Drug Discovery
  • Physical and Theoretical Chemistry
  • Organic Chemistry

Cite this

Ligand selectivity in the recognition of protoberberine alkaloids by hybrid-2 human telomeric G-quadruplex : Binding free energy calculation, fluorescence binding, and NMR experiments. / Deng, Nanjie; Xia, Junchao; Wickstrom, Lauren; Lin, Clement; Wang, Kaibo; He, Peng; Yin, Yunting; Yang, Danzhou.

In: Molecules, 21.04.2019.

Research output: Contribution to journalArticle

@article{c8ffe57342134193b3bc1f06b82c7e86,
title = "Ligand selectivity in the recognition of protoberberine alkaloids by hybrid-2 human telomeric G-quadruplex: Binding free energy calculation, fluorescence binding, and NMR experiments",
abstract = "The human telomeric G-quadruplex (G4) is an attractive target for developing anticancer drugs. Natural products protoberberine alkaloids are known to bind human telomeric G4 and inhibit telomerase. Among several structurally similar protoberberine alkaloids, epiberberine (EPI) shows the greatest specificity in recognizing the human telomeric G4 over duplex DNA and other G4s. Recently, NMR study revealed that EPI recognizes specifically the hybrid-2 form human telomeric G4 by inducing large rearrangements in the 5 0 -flanking segment and loop regions to form a highly extensive four-layered binding pocket. Using the NMR structure of the EPI-human telomeric G4 complex, here we perform molecular dynamics free energy calculations to elucidate the ligand selectivity in the recognition of protoberberines by the human telomeric G4. The MM-PB(GB)SA (molecular mechanics-Poisson Boltzmann/Generalized Born) Surface Area) binding free energies calculated using the Amber force fields bsc0 and OL15 correlate well with the NMR titration and binding affinity measurements, with both calculations correctly identifying the EPI as the strongest binder to the hybrid-2 telomeric G4 wtTel26. The results demonstrated that accounting for the conformational flexibility of the DNA-ligand complexes is crucially important for explaining the ligand selectivity of the human telomeric G4. While the MD-simulated (molecular dynamics) structures of the G-quadruplex-alkaloid complexes help rationalize why the EPI-G4 interactions are optimal compared with the other protoberberines, structural deviations from the NMR structure near the binding site are observed in the MD simulations. We have also performed binding free energy calculation using the more rigorous double decoupling method (DDM); however, the results correlate less well with the experimental trend, likely due to the difficulty of adequately sampling the very large conformational reorganization in the G4 induced by the protoberberine binding.",
keywords = "Binding free energy, G-quadruplex, MM-PB(Gb)SA, Molecular dynamics simulation, NMR, NMR titration, Protoberberine",
author = "Nanjie Deng and Junchao Xia and Lauren Wickstrom and Clement Lin and Kaibo Wang and Peng He and Yunting Yin and Danzhou Yang",
year = "2019",
month = "4",
day = "21",
doi = "10.3390/molecules24081574",
language = "English (US)",
journal = "Molecules",
issn = "1420-3049",
publisher = "Multidisciplinary Digital Publishing Institute (MDPI)",

}

TY - JOUR

T1 - Ligand selectivity in the recognition of protoberberine alkaloids by hybrid-2 human telomeric G-quadruplex

T2 - Binding free energy calculation, fluorescence binding, and NMR experiments

AU - Deng, Nanjie

AU - Xia, Junchao

AU - Wickstrom, Lauren

AU - Lin, Clement

AU - Wang, Kaibo

AU - He, Peng

AU - Yin, Yunting

AU - Yang, Danzhou

PY - 2019/4/21

Y1 - 2019/4/21

N2 - The human telomeric G-quadruplex (G4) is an attractive target for developing anticancer drugs. Natural products protoberberine alkaloids are known to bind human telomeric G4 and inhibit telomerase. Among several structurally similar protoberberine alkaloids, epiberberine (EPI) shows the greatest specificity in recognizing the human telomeric G4 over duplex DNA and other G4s. Recently, NMR study revealed that EPI recognizes specifically the hybrid-2 form human telomeric G4 by inducing large rearrangements in the 5 0 -flanking segment and loop regions to form a highly extensive four-layered binding pocket. Using the NMR structure of the EPI-human telomeric G4 complex, here we perform molecular dynamics free energy calculations to elucidate the ligand selectivity in the recognition of protoberberines by the human telomeric G4. The MM-PB(GB)SA (molecular mechanics-Poisson Boltzmann/Generalized Born) Surface Area) binding free energies calculated using the Amber force fields bsc0 and OL15 correlate well with the NMR titration and binding affinity measurements, with both calculations correctly identifying the EPI as the strongest binder to the hybrid-2 telomeric G4 wtTel26. The results demonstrated that accounting for the conformational flexibility of the DNA-ligand complexes is crucially important for explaining the ligand selectivity of the human telomeric G4. While the MD-simulated (molecular dynamics) structures of the G-quadruplex-alkaloid complexes help rationalize why the EPI-G4 interactions are optimal compared with the other protoberberines, structural deviations from the NMR structure near the binding site are observed in the MD simulations. We have also performed binding free energy calculation using the more rigorous double decoupling method (DDM); however, the results correlate less well with the experimental trend, likely due to the difficulty of adequately sampling the very large conformational reorganization in the G4 induced by the protoberberine binding.

AB - The human telomeric G-quadruplex (G4) is an attractive target for developing anticancer drugs. Natural products protoberberine alkaloids are known to bind human telomeric G4 and inhibit telomerase. Among several structurally similar protoberberine alkaloids, epiberberine (EPI) shows the greatest specificity in recognizing the human telomeric G4 over duplex DNA and other G4s. Recently, NMR study revealed that EPI recognizes specifically the hybrid-2 form human telomeric G4 by inducing large rearrangements in the 5 0 -flanking segment and loop regions to form a highly extensive four-layered binding pocket. Using the NMR structure of the EPI-human telomeric G4 complex, here we perform molecular dynamics free energy calculations to elucidate the ligand selectivity in the recognition of protoberberines by the human telomeric G4. The MM-PB(GB)SA (molecular mechanics-Poisson Boltzmann/Generalized Born) Surface Area) binding free energies calculated using the Amber force fields bsc0 and OL15 correlate well with the NMR titration and binding affinity measurements, with both calculations correctly identifying the EPI as the strongest binder to the hybrid-2 telomeric G4 wtTel26. The results demonstrated that accounting for the conformational flexibility of the DNA-ligand complexes is crucially important for explaining the ligand selectivity of the human telomeric G4. While the MD-simulated (molecular dynamics) structures of the G-quadruplex-alkaloid complexes help rationalize why the EPI-G4 interactions are optimal compared with the other protoberberines, structural deviations from the NMR structure near the binding site are observed in the MD simulations. We have also performed binding free energy calculation using the more rigorous double decoupling method (DDM); however, the results correlate less well with the experimental trend, likely due to the difficulty of adequately sampling the very large conformational reorganization in the G4 induced by the protoberberine binding.

KW - Binding free energy

KW - G-quadruplex

KW - MM-PB(Gb)SA

KW - Molecular dynamics simulation

KW - NMR

KW - NMR titration

KW - Protoberberine

UR - http://www.scopus.com/inward/record.url?scp=85064815645&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85064815645&partnerID=8YFLogxK

U2 - 10.3390/molecules24081574

DO - 10.3390/molecules24081574

M3 - Article

C2 - 31010072

AN - SCOPUS:85064815645

JO - Molecules

JF - Molecules

SN - 1420-3049

M1 - 1574

ER -