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

S288T mutation altering MmpL3 periplasmic domain channel and H-bond network: a novel dual drug resistance mechanism

CONTEXT

Mycobacterial membrane proteins Large 3 (MmpL3) is responsible for the transport of mycobacterial acids out of cell membrane to form cell wall, which is essential for the survival of Mycobacterium tuberculosis (Mtb) and has become a potent anti-tuberculosis target. SQ109 is an ethambutol (EMB) analogue, as a novel anti-tuberculosis drug, can effectively inhibit MmpL3, and has completed phase 2b-3 clinical trials. Drug resistance has always been the bottleneck problem in clinical treatment of tuberculosis. The S288T mutant of MmpL3 shows significant resistance to the inhibitor SQ109, while the specific action mechanism remains unclear. The results show that MmpL3 S288T mutation causes local conformational change with little effect on the global structure. With MmpL3 bound by SQ109 inhibitor, the distance between D710 and R715 increases resulting in H-bond destruction, but their interactions and proton transfer function are still restored. In addition, the rotation of Y44 in the S288T mutant leads to an obvious bend in the periplasmic domain channel and an increased number of contact residues, reducing substrate transport efficiency. This work not only provides a possible dual drug resistance mechanism of MmpL3 S288T mutant but also aids the development of novel anti-tuberculosis inhibitors.

METHODS

In this work, molecular dynamics (MD) and quantum mechanics (QM) simulations both were performed to compare inhibitor (i.e., SQ109) recognition, motion characteristics, and H-bond energy change of MmpL3 after S288T mutation. In addition, the WT_SQ109 complex structure was obtained by molecular docking program (Autodock 4.2); Molecular Mechanics/ Poisson Boltzmann Surface Area (MM-PBSA) and Solvated Interaction Energy (SIE) methods were used to calculate the binding free energies (∆Gbind); Geometric criteria were used to analyze the changes of hydrogen bond networks.

Uncovering the interactions between PME and PMEI at the gene and protein levels: Implications for the design of specific PMEI

CONTEXT

Pectin methylesterase inhibitor (PMEI) can specifically bind and inhibit the activity of pectin methylesterase (PME), which has been widely used in fruit and vegetable juice processing. However, the limited three-dimensional structure, unclear action mechanism, low thermal stability and biological activity of PMEI severely limited its application. In this work, molecular recognition and conformational changes of PME and PMEI were analyzed by various molecular simulation methods. Then suggestions were proposed for improving thermal stability and affinity maturation of PMEI through semi-rational design.

METHODS

Phylogenetic trees of PME and PMEI were established using the Maximum likelihood (ML) method. The results show that PME and PMEI have good sequence and structure conservation in various plants, and the simulated data can be widely adopted. In this work, MD simulations were performed using AMBER20 package and ff14SB force field. Protein interaction analysis indicates that H-bonds, van der Waals forces, and the salt bridge formed of K224 with ID116 are the main driving forces for mutual molecular recognition of PME and PMEI. According to the analyses of free energy landscape (FEL), conformational cluster, and motion, the association with PMEI greatly disrupts PME's dispersed functional motion mode and biological function. By monitoring the changes of residue contact number and binding free energy, IG35M/ IG35R: IT93F and IT113W/ IT113W: ID116W mutations contribute to thermal stability and affinity maturation of the PME-PMEI complex system, respectively. This work reveals the interaction between PME and PMEI at the gene and protein levels and provides options for modifying specific PMEI.

Recognition and release of uridine and hCNT3: From multivariate interactions to molecular design

As a vital target for the development of novel anti-cancer drugs, human concentrative nucleoside transporter 3 (hCNT3) has been widely concerned. Nevertheless, the lack of a comprehensive understanding of molecular interactions and motion mechanism has greatly hindered the development of novel inhibitors against hCNT3. In this paper, molecular recognition of hCNT3 with uridine was investigated with molecular docking, conventional molecular dynamics (CMD) simulations and adaptive steered molecular dynamics (ASMD) simulations; and then, the uridine derivatives with possibly highly inhibitory activity were designed. The result of CMD showed that more water-mediated H-bonds and lower binding free energy both explained higher recognition ability and transported efficiency of hCNT3. While during the ASMD simulation, nucleoside transport process involved the significant side-chain flip of residues F321 and Q142, a typical substrate-induced conformational change. By considering electronegativity, atomic radius, functional group and key H-bonds factors, 25 novel uridine derivatives were constructed. Subsequently, the receptor-ligand binding free energy was predicted by solvated interaction energy (SIE) method to determine the inhibitor c8 with the best potential performance. This work not only revealed molecular recognition and release mechanism of uridine with hCNT3, but also designed a series of uridine derivatives to obtain lead compounds with potential high activity.

Omicron-included mutation-induced changes in epitopes of SARS-CoV-2 spike protein and effectiveness assessments of current antibodies

The COVID-19 pandemic, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is spreading globally and continues to rage, posing a serious threat to human health and life quality. Antibody therapy and vaccines both have shown great efficacy in the prevention and treatment of COVID-19, whose development progress and adaptation range have attracted wide attention. However, with the emergence of variant strains of SARS-CoV-2, the neutralization activity of therapeutic or vaccine-induced antibodies may be reduced, requiring long-term virus monitoring and drug upgrade in response to its evolution. In this paper, conformational changes including continuous epitopes (CPs), discontinuous epitopes (DPs) and recognition interfaces of the three representative SARS-CoV-2 spike protein (SP) mutants (i.e., the Delta (B.1.617.2), Mu (B.1.621) and Omicron (B.1.1.529) strains), were analyzed to evaluate the effectiveness of current mainstream antibodies. The results showed that the conformation of SP wild type (WT) and mutants both remained stable, while the local antigenic epitopes underwent significant changes. Sufficient flexibility of SP CPs is critical for effective antibody recognition. The DPs of Delta, Mu and Omicron variants have showed stronger binding to human angiotensin converting enzyme-2 (hACE2) than WT; the possible drug resistance mechanisms of antibodies against three different epitopes (i.e., NTD_DP, RBD1_DP and RBD2_DP) were also proposed, respectively; the RBD2 of Delta, NTD of Mu, NTD and RBD2 of Omicron are deserve more attention in the subsequent design of next-generation vaccines. The simulation results not only revealed structural characteristics of SP antigenic epitopes, but also provided guidance for antibody modification, vaccine design and effectiveness evaluation.

Structural Insights on Retroviral DNA Integration: Learning from Foamy Viruses

Foamy viruses (FV) are retroviruses belonging to the Spumaretrovirinae subfamily. They are non-pathogenic viruses endemic in several mammalian hosts like non-human primates, felines, bovines, and equines. Retroviral DNA integration is a mandatory step and constitutes a prime target for antiretroviral therapy. This activity, conserved among retroviruses and long terminal repeat (LTR) retrotransposons, involves a viral nucleoprotein complex called intasome. In the last decade, a plethora of structural insights on retroviral DNA integration arose from the study of FV. Here, we review the biochemistry and the structural features of the FV integration apparatus and will also discuss the mechanism of action of strand transfer inhibitors.

An effective HIV-1 integrase inhibitor screening platform: Rationality validation of drug screening, conformational mobility and molecular recognition analysis for PFV integrase complex with viral DNA

As an important target for the development of novel anti-AIDS drugs, HIV-1 integrase (IN) has been widely concerned. However, the lack of a complete accurate crystal structure of HIV-1 IN greatly blocks the discovery of novel inhibitors. In this work, an effective HIV-1 IN inhibitor screening platform, namely PFV IN, was filtered from all species of INs. Next, the 40.8% similarity with HIV-1 IN, as well as the high efficiency of virtual screening and the good agreement between calculated binding free energies and experimental ones all proved PFV IN is a promising screening platform for HIV-1 IN inhibitors. Then, the molecular recognition mechanism of PFV IN by its substrate viral DNA and six naphthyridine derivatives (NRDs) inhibitors was investigated through molecular docking, molecular dynamics simulations and water-mediated interactions analyses. The functional partition of NRDs IN inhibitors could be divided into hydrophobic and hydrophilic ones, and the Mg²⁺ ions, water molecules and conserved DDE motif residues all interacted with the hydrophilic partition, while the bases in viral DNA and residues like Tyr212, Pro214 interacted with the hydrophobic one. Finally, the free energy landscape (FEL) and cluster analyses were performed to explore the molecular motion of PFV IN-DNA system. It is found that the association with NRDs inhibitors would obviously decrease the motion amplitude of PFV IN-DNA, which may be one of the most potential mechanisms of IN inhibitors. This work will provide a theoretical basis for the inhibitor design based on the structure of HIV-1 IN.

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.

Effect of HIV-1 Subtype C integrase mutations implied using molecular modeling and docking data

The degree of sequence variation in HIV-1 integrase genes among infected patients and their impact on clinical response to Anti retroviral therapy (ART) is of interest. Therefore, we collected plasma samples from 161 HIV-1 infected individuals for subsequent integrase gene amplification (1087 bp). Thus, 102 complete integrase gene sequences identified as HIV-1 subtype-C was assembled. This sequence data was further used for sequence analysis and multiple sequence alignment (MSA) to assess position specific frequency of mutations within pol gene among infected individuals. We also used biophysical geometric optimization technique based molecular modeling and docking (Schrodinger suite) methods to infer differential function caused by position specific sequence mutations towards improved inhibitor selection. We thus identified accessory mutations (usually reduce susceptibility) leading to the resistance of some known integrase inhibitors in 14% of sequences in this data set. The Stanford HIV-1 drug resistance database provided complementary information on integrase resistance mutations to deduce molecular basis for such observation. Modeling and docking analysis show reduced binding by mutants for known compounds. The predicted binding values further reduced for models with combination of mutations among subtype C clinical strains. Thus, the molecular basis implied for the consequence of mutations in different variants of integrase genes of HIV-1 subtype C clinical strains from South India is reported. This data finds utility in the design, modification and development of a representative yet an improved inhibitor for HIV-1 integrase.

Combined 3D-QSAR, molecular docking, and molecular dynamics study of tacrine derivatives as potential acetylcholinesterase (AChE) inhibitors of Alzheimer's disease

Acetylcholinesterase (AChE) is one of the key targets of drugs for treating Alzheimer's disease (AD). Tacrine is an approved drug with AChE-inhibitory activity. In this paper, 3D-QSAR, molecular docking, and molecular dynamics were carried out in order to study 60 tacrine derivatives and their AChE-inhibitory activities. 3D-QSAR modeling resulted in an optimal CoMFA model with q(2) = 0.552 and r(2) = 0.983 and an optimal CoMSIA model with q(2) = 0.581 and r(2) = 0.989. These QSAR models also showed that the steric and H-bond fields of these compounds are important influences on their activities. The interactions between these inhibitors and AChE were further explored through molecular docking and molecular dynamics simulation. A few key residues (Tyr70, Trp84, Tyr121, Trp279, and Phe330) at the binding site of AChE were identified. The results of this study improve our understanding of the mechanisms of AChE inhibitors and afford valuable information that should aid the design of novel potential AChE inhibitors. Graphical Abstract Superposition of backbone atoms of the lowest-energy structure obtained from MD simulation (magenta) onto those of the structure of the initial molecular docking model (green).

Development of HIV-1 integrase inhibitors: recent molecular modeling perspectives

Of the three viral enzymes essential to HIV replication, HIV-1 integrase (IN) is gaining popularity as a target for the antiviral therapy of AIDS. Substantial work focusing on IN has been done over the past three decades, which has facilitated and led to the approval of three drugs. Here, we discuss in detail the development of IN inhibitors between January 2012 and May 2014, with a particular focus on molecular simulation. We highlight controversial aspects of computational drug design and refer to alternative practices where appropriate. The analysis of these computational approaches provides some useful clues to the possible future discovery of novel IN inhibitors.

Study on Molecular Recognition betweenEuphorbia FactorL713283 andβ-Tubulin via Molecular Simulation Methods

Euphorbia factorL713283 is a new lathyrane diterpene isolated fromEuphorbia lathyrisand shows strong anticancer activity. By using molecular similarity analysis,β-tubulin was identified as one of the possible targets of L713283. We further investigated the binding modes of L713283 withβ-tubulin using molecular docking and molecular dynamics (MD) simulation methods. The results indicated that the binding site betweenβ-tubulin and L713283 was composed of the four regions, that is, residues Phe20~Glu27, Leu225~Thr232, Phe270~Gly277, and Ile356~Met363. MM/GBSA method was used to calculate the binding free energy and determine the key residues for the association of L713283 withβ-tubulin. It was found that nonpolar interactions made the major contributions for the binding. In addition, we compared the binding pocket and motion modes of L713283-free and L713283-boundβ-tubulin systems. It is proposed that L713283 may bind toβ-tubulin and favor the formation ofαβ-tubulin dimmer. This work provides possible explanation for molecular mechanism of the anticancer agent L713283, and the strategy used here could benefit the investigation of possible target profile for those bioactive agents with unknown mechanisms.

Reduced HIV-1 integrase flexibility as a mechanism for raltegravir resistance

HIV-1 integrase is an essential enzyme necessary for the replication of the HIV virus as it catalyzes the insertion of the viral genome into the host chromosome. Raltegravir was the first integrase inhibitor approved by the FDA for antiretroviral treatment. HIV patients on raltegravir containing regimens often develop drug resistance mutations at residue 140 and 148 in the catalytic 140's loop resulting in a 5-10 fold decrease in susceptibility to raltegravir. Obtaining crystallographic structure information on the Q148H/R, G140S/A primary and secondary mutations has been elusive. Using 10 ns molecular dynamics simulations, we present a detailed analysis of the structural changes induced by these mutations. The formation frequency of a transient helix in the catalytic 140's loop is increased and the length of this helix is extended from 3-residues to 4 in the mutants relative to the wild type. This helix causes reduced flexibility in the protein active site and therefore serves as a gating mechanism restricting the access of raltegravir to the integrase binding pocket. These results suggest that resistance to raltegravir occurs through a common mechanism of altering the formation frequency of transient secondary structures such as α2 and β5 in addition to the conformational changes in the 140's loop therefore decreasing the flexibility of the HIV-1 integrase protein. The reduced integrase flexibility serves as a mechanism of resistance to raltegravir.

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.