Relation between flexibility and positively selected HIV-1 protease mutants against inhibitors

Unveiling mutation effects on the structural dynamics of the main protease from SARS-CoV-2 with hybrid simulation methods

The main protease of SARS-CoV-2 (called M^pro or 3CL^pro) is essential for processing polyproteins encoded by viral RNA. Several M^pro mutations were found in SARS-CoV-2 variants, which are related to higher transmissibility, pathogenicity, and resistance to neutralization antibodies. Macromolecules adopt several favored conformations in solution depending on their structure and shape, determining their dynamics and function. In this study, we used a hybrid simulation method to generate intermediate structures along the six lowest frequency normal modes and sample the conformational space and characterize the structural dynamics and global motions of WT SARS-CoV-2 M^pro and 48 mutations, including mutations found in P.1, B.1.1.7, B.1.351, B.1.525 and B.1.429+B.1.427 variants. We tried to contribute to the elucidation of the effects of mutation in the structural dynamics of SARS-CoV-2 M^pro. A machine learning analysis was performed following the investigation regarding the influence of the K90R, P99L, P108S, and N151D mutations on the dimeric interface assembling of the SARS-CoV-2 M^pro. The parameters allowed the selection of potential structurally stable dimers, which demonstrated that some single surface aa substitutions not located at the dimeric interface (K90R, P99L, P108S, and N151D) are able to induce significant quaternary changes. Furthermore, our results demonstrated, by a Quantum Mechanics method, the influence of SARS-CoV-2 M^pro mutations on the catalytic mechanism, confirming that only one of the chains of the WT and mutant SARS-CoV-2 M^pros are prone to cleave substrates. Finally, it was also possible to identify the aa residue F140 as an important factor related to the increasing enzymatic reactivity of a significant number of SARS-CoV-2 M^pro conformations generated by the normal modes-based simulations.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach

Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?G_bind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.

Flexibility and extracellular opening determine the interaction between ligands and insect sulfakinin receptors

Despite their fundamental importance for growth, the mechanisms that regulate food intake are poorly understood. Our previous work demonstrated that insect sulfakinin (SK) signaling is involved in inhibiting feeding in an important model and pest insect, the red flour beetle Tribolium castaneum. Because the interaction of SK peptide and SK receptors (SKR) initiates the SK signaling, we have special interest on the structural factors that influence the SK-SKR interaction. First, the three-dimensional structures of the two T. castaneum SKRs (TcSKR1 and TcSKR2) were generated from molecular modeling and they displayed significance in terms of the outer opening of the cavity and protein flexibility. TcSKR1 contained a larger outer opening of the cavity than that in TcSKR2, which allows ligands a deep access into the cavity through cell membrane. Second, normal mode analysis revealed that TcSKR1 was more flexible than TcSKR2 during receptor-ligand interaction. Third, the sulfated SK (sSK) and sSK-related peptides were more potent than the nonsulfated SK, suggesting the importance of the sulfate moiety.

Recovery of the wild type atomic flexibility in the HIV-1 protease double mutants

The emergence of drug resistant mutations due to the selective pressure exerted by antiretrovirals, including protease inhibitors (PIs), remains a major problem in the treatment of AIDS. During PIs therapy, the occurrence of primary mutations in the wild type HIV-1 protease reduces both the affinity for the inhibitors and the viral replicative capacity compared to the wild type (WT) protein, but additional mutations compensate for this reduced viral fitness. To investigate this phenomenon from the structural point of view, we combined Molecular Dynamics and Normal Mode Analysis to analyze and compare the variations of the flexibility of C-alpha atoms and the differences in hydrogen bond (h-bond) network between the WT and double mutants. In most cases, the flexibility profile of the double mutants was more often similar to that of the WT than to that of the related single base mutants. All single mutants showed a significant alteration in h-bond formation compared to WT. Most of the significant changes occur in the border between the flap and cantilever regions. We found that all the considered double mutants have their h-bond pattern significantly altered in comparison to the respective single base mutants affecting their flexibility profile that becomes more similar to that of WT. This WT flexibility restoration in the double mutants appears as an important factor for the HIV-1 fitness recovery observed in patients.

Synthetic lethals in HIV: ways to avoid drug resistance : Running title: Preventing HIV resistance

Background

RNA viruses rapidly accumulate genetic variation, which can give rise to synthetic lethal (SL) and deleterious (SD) mutations. Synthetic lethal mutations (non-lethal when alone but lethal when combined in one genome) have been studied to develop cancer therapies. This principle can also be used against fast-evolving RNA-viruses. Indeed, targeting protein sites involved in SD + SL interactions with a drug would render any mutation of such sites, lethal.

Results

Here, we set up a strategy to detect intragenic pairs of SL and SD at the surface of the protein to predict less escapable drug target sites. For this, we detected SD + SL, studying HIV protease (PR) and reverse transcriptase (RT) sequence alignments from two groups of VIH⁺ individuals: treated with drugs (T) or not (NT). Using a series of statistical approaches, we were able to propose bona fide SD + SL couples. When focusing on spatially close co-variant SD + SL couples at the surface of the protein, we found 5 SD + SL groups (2 in the protease and 3 in the reverse transcriptase), which could be good candidates to form pockets to accommodate potential drugs.

Conclusions

Thus, designing drugs targeting these specific SD + SL groups would not allow the virus to mutate any residue involved in such groups without losing an essential function. Moreover, we also show that the selection pressure induced by the treatment leads to the appearance of new mutations, which change the mutational landscape of the protein. This drives the existence of differential SD + SL couples between the drug-treated and non-treated groups. Thus, new anti-viral drugs should be designed differently to target such groups.

Reviewers

This article was reviewed by Neil Greenspan Csaba Pal and István Simon.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0044-y) contains supplementary material, which is available to authorized users.

Chemical specificity and conformational flexibility in proteinase-inhibitor interaction: scaffolds for promiscuous binding

One of the most important roles of proteins in cellular milieu is recognition of other biomolecules including other proteins. Protein-protein complexes are involved in many essential cellular processes. Interfaces of protein-protein complexes are traditionally known to be conserved in evolution and less flexible than other solvent interacting tertiary structural surface. But many examples are emerging where these features do not hold good. An understanding of inter-play between flexibility and sequence conservation is emerging, providing a fresh dimension to the paradigm of sequence-structure-function relationship. The functional manifestation of the inter-relation between sequence conservation and flexibility of interface is exemplified in this review using proteinase-inhibitor protein complexes.

Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility

Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.

A cell-based reporter assay for screening for EcR agonist/antagonist activity of natural ecdysteroids in Lepidoptera (Bm5) and Diptera (S2) cell cultures, followed by modeling of ecdysteroid-EcR interactions and normal mode analysis

Ecdysteroid signal transduction is a key process in insect development and therefore an important target for insecticide development. We employed an in vitro cell-based reporter bioassay for the screening of potential ecdysone receptor (EcR) agonistic and antagonistic compounds. Natural ecdysteroids were assayed with ecdysteroid-responsive cell line cultures that were transiently transfected with the reporter plasmid ERE-b.act.luc. We used the dipteran Schneider S2 cells of Drosophila melanogaster and the lepidopteran Bm5 cells of Bombyx mori, representing important pest insects in medicine and agriculture. Measurements showed an EcR agonistic activity only for cyasterone both in S2 (EC50=3.3μM) and Bm5 cells (EC50=5.3μM), which was low compared to that of the commercial dibenzoylhydrazine-based insecticide tebufenozide (EC50=0.71μM and 0.00089μM, respectively). Interestingly, a strong antagonistic activity was found for castasterone in S2 cells with an IC50 of 0.039μM; in Bm5 cells this effect only became visible at much higher concentrations (IC50=18μM). To gain more insight in the EcR interaction, three-dimensional modeling of dipteran and lepidopteran EcR-LBD was performed. In conclusion, we showed that the EcR cell-based reporter bioassay tested here is a useful and practical tool for the screening of candidate EcR agonists and antagonists. The docking experiments as well as the normal mode analysis provided evidence that the antagonist activity of castasterone may be through direct binding with the receptor with specific changes in protein flexibility. The search for new ecdysteroid-like compounds may be particularly relevant for dipterans because the activity of dibenzoylhydrazines appears to be correlated with an extension of the EcR-LBD binding pocket that is prominent in lepidopteran receptors but less so in the modeled dipteran structure.