Dataset Information

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

ABSTRACT: Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp-ritonavir protein-ligand systems as cases, we observed that ligand association is a multistep process involving diffusion, localization, and conformational rearrangements of the protein, ligand, and water molecules. Moreover, these two ligands preferred different routes of binding, which reflect two well-known binding mechanisms: induced-fit and conformation selection models. Our study shows that xk263 has a stronger capacity for desolvating surrounding water molecules, thereby inducing a semiopen conformation of the HIVp flaps (induced-fit model). In contrast, the slow dehydration characteristic of ritonavir allows for gradual association with the binding pocket of HIVp when the protein's flap conformation is fully open (conformation selection model). By studying the mechanism of ligand association and understanding the role of solvent molecules during the binding event, we can obtain a different perspective on the mechanism of macromolecule recognition, providing insights into drug discovery.

SUBMITTER: Huang YM

PROVIDER: S-EPMC5499997 | biostudies-literature | 2017 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

Huang Yu-Ming M YM Raymundo Mark Anthony V MA Chen Wei W Chang Chia-En A CA

Biochemistry 20170221 9

Equilibrium constants, together with kinetic rate constants of binding, are key factors in the efficacy and safety of drug compounds, informing drug design. However, the association pathways of protein-ligand binding, which contribute to their kinetic behaviors, are little understood. In this work, we used unbiased all-atom molecular dynamics (MD) simulations with an explicit solvent model to study the association processes of protein-ligand binding. Using the HIV protease (HIVp)-xk263 and HIVp- ...[more]

PMID: 28060481

Similar Datasets

Project description:Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.

Project description:A new family of dicopper(I) complexes [CuI2RL](X)2 (R=H, 1X, R=tBu, 2X and R=NO2, 3X, X=CF3SO3, ClO4, SbF6, or BArF, BArF=[B{3,5-(CF3)2C6H3}4]-), where RL is a Schiff-base ligand containing two tridentate binding sites linked by a xylyl spacer, has been prepared and characterized, and its reaction with O2 has been studied. The complexes were designed with the aim of reproducing structural aspects of the active site of type 3 dicopper proteins; they contain two three-coordinate copper sites and a rather flexible podand ligand backbone. The solid-state structures of 1ClO4, 2CF3SO3, 2ClO4, and 3BArF.CH3CN have been established by single-crystal X-ray diffraction analysis. 1ClO4 adopts a polymeric structure in the solid state while 2CF3SO3, 2ClO4, and 3BArF.CH3CN are monomeric. The complexes have been studied in solution by means of 1H and 19F NMR spectroscopy, which put forward the presence of dynamic processes. 1-3BArF and 1-3CF3SO3 in acetone react rapidly with O2 to generate metaestable [CuIII2(mu-O)2(RL)]2+ 1-3(O2) and [CuIII2(mu-O)2(CF3SO3)(RL)]+ 1-3(O2)(CF3SO3) species, respectively, that have been characterized by UV-vis spectroscopy and resonance Raman analysis. Instead, reaction of 1-3BArF with O2 in CH2Cl2 results in intermolecular O2 binding. DFT methods have been used to study the chemical identities and structural parameters of the O2 adducts, and the relative stability of the CuIII2(mu-O)2 form with respect to the CuII2(mu-eta2:eta2-O2) isomer. The reaction of 1X, X = CF3SO3 and BArF, with O2 in acetone has been studied by stopped-flow UV-vis exhibiting an unexpected very fast reaction rate (k=3.82(4)x10(3) M-1 s-1, DeltaH=4.9+/-0.5 kJ.mol-1, DeltaS=-148+/-5 J.K-1.mol-1), nearly 3 orders of magnitude faster than in the parent [CuI2(m-XYLMeAN)]2+. Thermal decomposition of 1-3(O2) does not result in aromatic hydroxylation. The mechanism and kinetics of O2 binding to 1X (X=CF3SO3 and BArF) are discussed and compared with those associated with selected examples of reported models of O2-processing copper proteins. A synergistic role of the copper ions in O2 binding and activation is clearly established from this analysis.

Dataset Information

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

Publications

Mechanism of the Association Pathways for a Pair of Fast and Slow Binding Ligands of HIV-1 Protease.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets