Study on the molecular mechanism of inhibiting HIV-1 integrase by EBR28 peptide via molecular modeling approach

Inhibition Mechanism of Indoleamine 2, 3-Dioxygenase 1 (IDO1) by Amidoxime Derivatives and Its Revelation in Drug Design: Comparative Molecular Dynamics Simulations

For cancer treatment, in addition to the three standard therapies of surgery, chemotherapy, and radiotherapy, immunotherapy has become the fourth internationally-recognized alternative treatment. Indoleamine 2, 3-dioxygenase 1 (IDO1) catalyzes the conversion of tryptophan to kynurenine causing lysine depletion, which is an important target in the research and development of anticancer drugs. Epacadostat (INCB024360) is currently one of the most potent IDO1 inhibitors, nevertheless its inhibition mechanism still remains elusive. In this work, comparative molecular dynamics simulations were performed to reveal that the high inhibitory activity of INCB024360 mainly comes from two aspects: disturbing the ligand delivery tunnel and then preventing small molecules such as oxygen and water molecules from accessing the active site, as well as hindering the shuttle of substrate tryptophan with product kynurenine through the heme binding pocket. The scanning of key residues showed that L234 and R231 residues both were crucial to the catalytic activity of IDO1. With the association with INCB024360, L234 forms a stable hydrogen bond with G262, which significantly affects the spatial position of G262-A264 loop and then greatly disturbs the orderliness of ligand delivery tunnel. In addition, the cleavage of hydrogen bond between G380 and R231 increases the mobility of the GTGG conserved region, leading to the closure of the substrate tryptophan channel. This work provides new ideas for understanding action mechanism of amidoxime derivatives, improving its inhibitor activity and developing novel inhibitors of IDO1.

The global motion affecting electron transfer in Plasmodium falciparum type II NADH dehydrogenases: a novel non-competitive mechanism for quinoline ketone derivative inhibitors

With the emergence of drug-resistant Plasmodium falciparum, the treatment of malaria has become a significant challenge; therefore, the development of antimalarial drugs acting on new targets is extremely urgent. In Plasmodium falciparum, type II nicotinamide adenine dinucleotide (NADH) dehydrogenase (NDH-2) is responsible for catalyzing the transfer of two electrons from NADH to flavin adenine dinucleotide (FAD), which in turn transfers the electrons to coenzyme Q (CoQ). As an entry enzyme for oxidative phosphorylation, NDH-2 has become one of the popular targets for the development of new antimalarial drugs. In this study, reliable motion trajectories of the NDH-2 complex with its co-factors (NADH and FAD) and inhibitor, RYL-552, were obtained by comparative molecular dynamics simulations. The influence of cofactor binding on the global motion of NDH-2 was explored through conformational clustering, principal component analysis and free energy landscape. The molecular interactions of NDH-2 before and after its binding with the inhibitor RYL-552 were analyzed, and the key residues and important hydrogen bonds were also determined. The results show that the association of RYL-552 results in the weakening of intramolecular hydrogen bonds and large allosterism of NDH-2. There was a significant positive correlation between the angular change of the key pocket residues in the NADH-FAD-pockets that represents the global functional motion and the change in distance between NADH-C4 and FAD-N5 that represents the electron transfer efficiency. Finally, the possible non-competitive inhibitory mechanism of RYL-552 was proposed. Specifically, the association of inhibitors with NDH-2 significantly affects the global motion mode of NDH-2, leading to widening of the distance between NADH and FAD through cooperative motion induction; this reduces the electron transfer efficiency of the mitochondrial respiratory chain. The simulation results provide useful theoretical guidance for subsequent antimalarial drug design based on the NDH-2 structure and the respiratory chain electron transfer mechanism.

PD-L1 Nanobody Competitively Inhibits the Formation of the PD-1/PD-L1 Complex: Comparative Molecular Dynamics Simulations

The anti-PD-L1 monoclonal antibody (mAb) targeting PD-1/PD-L1 immune checkpoint has achieved outstanding results in clinical application and has become one of the most popular anti-cancer drugs. The mechanism of molecular recognition and inhibition of PD-L1 mAbs is not yet clear, which hinders the subsequent antibody design and modification. In this work, the trajectories of PD-1/PD-L1 and nanobody/PD-L1 complexes were obtained via comparative molecular dynamics simulations. Then, a series of physicochemical parameters including hydrogen bond, dihedral angle distribution, pKa value and binding free energy, and so forth, were all comparatively analyzed to investigate the recognition difference between PD-L1 and PD-1 and nanobody. Both LR113 (the amino acid residues in PD-L1 are represented by the lower left sign of L) and LR125 residues of PD-L1 undergo significant conformational change after association with mAbs, which dominates a strong electrostatic interaction. Solvation effect analysis revealed that solvent-water enhanced molecular recognition between PD-L1 and nanobody. By combining the analyses of the time-dependent root mean squared fluctuation (RMSF), free energy landscape, clustering and energy decomposition, the potential inhibition mechanism was proposed that the nanobody competitively and specifically bound to the β-sheet groups of PD-L1, reduced the PD-L1’s flexibility and finally blocked the formation of PD-1/PD-L1 complex. Based on the simulation results, site-directed mutagenesis of ND99 (the amino acid residues in Nano are displayed by the lower left sign of N) and NQ116 in the nanobody may be beneficial for improving antibody activity. This work offers some structural guidance for the design and modification of anticancer mAbs based on the structure of the PD-1/PD-L1 complex.

3D-QSAR, Molecular Docking and Molecular Dynamics Simulation of Pseudomonas aeruginosa LpxC Inhibitors

As an important target for the development of novel antibiotics, UDP-3-O-(R-3-hydroxymyristoyl)-N-acetylglucosamine deacetylase (LpxC) has been widely studied. Pyridone methylsulfone hydroxamate (PMH) compounds can effectively inhibit the catalytic activity of LpxC. In this work, the three-dimensional quantitative structure-activity relationships of PMH inhibitors were explored and models with good predictive ability were established using comparative molecular field analysis and comparative molecular similarity index analysis methods. The effect of PMH inhibitors' electrostatic potential on the inhibitory ability of Pseudomonas aeruginosa LpxC (PaLpxC) is revealed at the molecular level via molecular electrostatic potential analyses. Then, two molecular dynamics simulations for the PaLpxC and PaLpxC-inhibitor systems were also performed respectively to investigate the key residues of PaLpxC hydrolase binding to water molecules. The results indicate that orderly alternative water molecules can form stable hydrogen bonds with M62, E77, T190, and H264 in the catalytic center, and tetracoordinate to zinc ion along with H78, H237, and D241. It was found that the conformational transition space of PaLpxC changes after association with PMH inhibitors through free energy landscape and cluster analyses. Finally, a possible inhibitory mechanism of PMH inhibitors was proposed, based on our molecular simulation. This paper will provide a theoretical basis for the molecular design of LpxC-targeting antibiotics.

CPe-III-S Metabolism in Vitro and in Vivo and Molecular Simulation of Its Metabolites Using a p53-R273H Mutant

It was previously found that CPe-III-S, synthesized according to the chickpea peptide CPe-III (RQSHFANAQP), exhibited an antiproliferative effect. The aim of this study was to investigate the antiproliferative mechanism of CPe-III-S. CPe-III-S was treated by pepsin and trypsin in a simulated gastrointestinal digestion environment as well as in an animal experiment. With HPLC-ESI-MS analysis, three peptide fragments of Ser-His, His-Phe, and Ala-Asn-Ala-Gln were identified. Ser-His was the only common product from both in vitro and in vivo environments. The specific bindings between three peptides and p53-R273H were performed by molecular docking, and the molecular dynamic simulation between Ser-His and p53-R273H revealed the stability of the binding complex. The binding free energy of the complex was -12.56 ± 1.03 kcal/mol with a reliable hydrogen bond between the ligand and Thr284 of p53. Ser-His may restore mutant p53-R273H activity or inhibit its binding with a downstream signal. This metabolite is a potential anticancer factor for suppressing cell proliferation.

Substrate recognition and motion mode analyses of PFV integrase in complex with viral DNA via coarse-grained models

HIV-1 integrase (IN) is an important target in the development of drugs against the AIDS virus. Drug design based on the structure of IN was markedly hampered due to the lack of three-dimensional structure information of HIV-1 IN-viral DNA complex. The prototype foamy virus (PFV) IN has a highly functional and structural homology with HIV-1 IN. Recently, the X-ray crystal complex structure of PFV IN with its cognate viral DNA has been obtained. In this study, both Gaussian network model (GNM) and anisotropy network model (ANM) have been applied to comparatively investigate the motion modes of PFV DNA-free and DNA-bound IN. The results show that the motion mode of PFV IN has only a slight change after binding with DNA. The motion of this enzyme is in favor of association with DNA, and the binding ability is determined by its intrinsic structural topology. Molecular docking experiments were performed to gain the binding modes of a series of diketo acid (DKA) inhibitors with PFV IN obtained from ANM, from which the dependability of PFV IN-DNA used in the drug screen for strand transfer (ST) inhibitors was confirmed. It is also found that the functional groups of keto-enol, bis-diketo, tetrazole and azido play a key role in aiding the recognition of viral DNA, and thus finally increase the inhibition capability for the corresponding DKA inhibitor. Our study provides some theoretical information and helps to design anti-AIDS drug based on the structure of IN.

Study on the interactions between diketo-acid inhibitors and prototype foamy virus integrase-DNA complex via molecular docking and comparative molecular dynamics simulation methods

Human immunodeficiency virus type 1 (HIV-1) integrase (IN) is an important drug target for anti-acquired immune deficiency disease (AIDS) treatment and diketo-acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN. Due to lack of three-dimensional structures including detail interactions between HIV-1 IN and its substrate viral DNA, the drug design and screening platform remains incompleteness and deficient. In addition, the action mechanism of DKA inhibitors with HIV-1 IN is not well understood. In view of the high homology between the structure of prototype foamy virus (PFV) IN and that of HIV-1 IN, we used PFV IN as a surrogate model for HIV-1 IN to investigate the inhibitory mechanism of raltegravir (RLV) and the binding modes with a series of DKA inhibitors. Firstly, molecular dynamics simulations of PFV IN, IN-RLV, IN-DNA, and IN-DNA-RLV systems were performed for 10 ns each. The interactions and inhibitory mechanism of RLV to PFV IN were explored through overall dynamics behaviors, catalytic loop conformation distribution, and hydrogen bond network analysis. The results show that the coordinated interactions of RLV with IN and viral DNA slightly reduce the flexibility of catalytic loop region of IN, and remarkably restrict the mobility of the CA end of viral DNA, which may lead to the partial loss of the inhibitory activity of IN. Then, we docked a series of DKA inhibitors into PFV IN-DNA receptor and obtained the IN-DNA-inhibitor complexes. The docking results between PFV IN-DNA and DKA inhibitors agree well with the corresponding complex of HIV-1 IN, which proves the dependability of PFV IN-DNA used for the anti-AIDS drug screening. Our study may help to make clear some theoretical questions and to design anti-AIDS drug based on the structure of IN.

Peptides that inhibit HIV-1 integrase by blocking its protein-protein interactions

HIV-1 integrase (IN) is one of the key enzymes in the viral replication cycle. It mediates the integration of viral cDNA into the host cell genome. IN activity requires interactions with several viral and cellular proteins, as well as IN oligomerization. Inhibition of IN is an important target for the development of anti-HIV therapies, but there is currently only one anti-HIV drug used in the clinic that targets IN. Several other small-molecule anti-IN drug leads are either undergoing clinical trials or in earlier stages of development. These molecules specifically inhibit one of the IN-mediated reactions necessary for successful integration. However, small-molecule inhibitors of protein-protein interactions are difficult to develop. In this review, we focus on peptides that inhibit IN. Peptides have advantages over small-molecule inhibitors of protein-protein interactions: they can mimic the structures of the binding domains within proteins, and are large enough to competitively inhibit protein-protein interactions. The development of peptides that bind IN and inhibit its protein-protein interactions will increase our understanding of the IN mode of action, and lead to the development of new drug leads, such as small molecules derived from these peptides, for better anti-HIV therapy.

STUDY ON THE INTERACTIONS OF Smac MIMETICS WITH XIAP-BIR3 DOMAIN BY DOCKING AND MOLECULAR DYNAMICS SIMULATIONS

Upon receiving an apoptotic stimulus, the mature mitochondrial protein second mitochondria-derived activator of caspases (Smac)/direct IAP-binding protein with low PI (DIABLO), which could be released from mitochondria into the cytosol together with cytochrome C , specifically binds to inhibitor of apoptosis proteins (IAPs) and relieves the inhibitory effect of caspase, thus promotes cell death. Some artificial small molecules (called Smac mimetics) can mimic the N-terminal four residues Ala1-Val2-Pro3-Ile4 (AVPI) sequence of mitochondrial protein Smac, and competitively bind to X-linked inhibitor of apoptosis protein baculoviral IAP repeats (XIAP-BIR3) domain with caspase-9, which leads to the removal of the inhibition of caspase-9 by XIAP and induce apoptosis. To gain an insight into the nature of XIAP-BIR3 domain recognizing Smac mimetics, we used docking and molecular dynamics simulations methods to study four representative Smac mimetics. The docking results show that the orientations of these backbones of ligands are identical with that of AVPI in the binding pocket. Each ligand corresponds to two competitive conformations, which are called extended and bended conformations. The results of molecular dynamics simulations show that the extended conformation is more stable, and the calculations of energy decomposition reveal that the residue Thr308 makes the strongest interaction with XIAP-BIR3. In addition, Asp309, Glu314, and Trp323 are indispensable for XIAP-BIR3 recognizing and binding Smac mimetics.

Molecular dynamics simulations of the HIV-1 integrase dimerization interface: guidelines for the design of a novel class of integrase inhibitors

HIV-1 integrase (IN) is a validated target of anti-AIDS research. The classical approach of designing active-site directed ligands has largely been exploited. A promising alternative strategy to inactivate the enzyme is to prevent the formation of IN dimers. The rational design of dimerization inhibitors, however, is hampered by the lack of relevant structural data about the targeted monomeric form. Therefore, we performed molecular dynamics simulations and subsequent analyses to gain insight into the structural features of the IN catalytic-core-domain dimerization interface. As a result, the formation of a groove and a cavity along the dimerization interface of the IN monomer could be revealed. Both were shown to be suited for accommodating an inhibitory peptide. The results form a valuable basis for the design of ligands targeting the dimerization interface and, thus, of a whole new class of HIV-1 integrase inhibitors.

Communications: Electron polarization critically stabilizes the Mg2+ complex in the catalytic core domain of HIV-1 integrase

In this paper, we present a detailed dynamics study of the catalytic core domain (CCD) of HIV-1 integrase using both polarized and nonpolarized force fields. The numerical results reveal the critical role of protein polarization in stabilizing Mg(2+) coordination complex in CCD. Specifically, when nonpolarized force field is used, a remarkable drift of the Mg(2+) complex away from its equilibrium position is observed, which causes the binding site blocked by the Mg(2+) complex. In contrast, when polarized force field is employed in MD simulation, HIV-1 integrase CCD structure is stabilized and both the position of the Mg(2+) complex and the binding site are well preserved. The detailed analysis shows the transition of alpha-helix to 3(10)-helix adjacent to the catalytic loop (residues 139-147), which correlates with the dislocation of the Mg(2+) complex. The current study demonstrates the importance of electronic polarization of protein in stabilizing the metal complex in the catalytic core domain of HIV-1 integrase.

Raltegravir: molecular basis of its mechanism of action

Integration of the HIV-1 viral DNA generated by reverse transcription of the RNA genome into the host cell chromosomes is a key step of viral replication, catalyzed by the viral integrase. In October 2007, the first integrase inhibitor, raltegravir, was approved for clinical use under the name of Isentress™. The results of the various clinical trials that have evaluated raltegravir have been very encouraging with regard to the immunological and virological efficacy and tolerance. However, as observed for other anti-retrovirals, specific resistance mutations have been identified in patients failing to respond to treatment with raltegravir. Although knowledge of the integrase structural biology remains fragmentary, the structures and modeling data available might provide relevant clues on the origin of the emergence of these resistance mutations. In this review, we describe the mechanism of action of this drug and the main data relating to its use in vivo, together with recent structural data important to our understanding of the origin of viral resistance.

In search of second-generation HIV integrase inhibitors: targeting integration beyond strand transfer

Highly active antiretroviral therapy combines antiviral drugs targeting different steps in the HIV replication cycle in order to reduce viral loads in patients to undetectable levels. Since HIV readily develops resistance and can therefore escape the action of existing drugs, novel drugs with novel mechanisms of action must be developed. The integration of the viral genome into the human genome is an essential and critical replication step that is catalyzed by the viral integrase with the help of cellular cofactors. Although HIV-1 integrase has been studied for more than two decades, the first integrase inhibitor, raltegravir, was only recently approved for clinical use. A second compound, elvitegravir, is currently in advanced clinical trials. Both drugs interfere with the strand-transfer reaction of integrase. Due to the complexity and multistep nature of the integration reaction, several other functions of integrase can be exploited for drug discovery. In this review, we will describe these alternative strategies to inhibit integration. They have recently attracted considerable interest for the development of second-generation integrase inhibitors.

Dinâmica molecular: teoria e aplicações em planejamento de fármacos

Dinâmica Molecular (DM) é uma ferramenta computacional poderosa usada em Química Medicinal para o planejamento racional de fármacos. DM é uma extensão da Mecânica Molecular, onde o comportamento dinâmico de um sistema molecular é simulado através da integração numérica das equações de movimento. Esta técnica tem sido usada extensivamente para auxiliar e complementar o planejamento de novos ligantes de um alvo terapêutico, bem como estimar a sua potência. Este artigo enfoca a teoria básica da DM clássica e suas importantes aplicações no planejamento racional de potenciais compostos bioativos, particularmente compostos com atividade anti-HIV.