Study of the Conformational Dynamics of the Catalytic Loop of WT and G140A/G149A HIV-1 Integrase Core Domain Using Reversible Digitally Filtered Molecular Dynamics

Scaffold-based selective SHP2 inhibitors design using core hopping, molecular docking, biological evaluation and molecular simulation

PTPN11 (coding the gene of SHP2), a classic non-receptor protein tyrosine phosphatase, is implicated in multiple cell signaling pathway. Abnormal activation of SHP2 has been shown to contribute to a variety of human diseases, including Juvenile myelomonocytic leukemia (JMML), Noonan syndrome and tumors. Thus, the SHP2 inhibitors have important therapeutic value. Here, based on the compound PubChem CID 8,478,960 (IC₅₀ = 45.01 μM), a series of thiophene [2,3-d] pyrimidine derivatives (IC₅₀ = 0.4-37.87 μM) were discovered as novel and efficient inhibitors of SHP2 through powerful "core hopping" and CDOCKER technology. Furthermore, the SHP2-PTP phosphatase activity assay indicated that Comp#5 (IC₅₀ = 0.4 μM) was the most active SHP2 inhibitor. Subsequently, the effects of Comp#5 on the structure and function of SHP2 were investigated through molecular dynamics (MD) simulation and post-kinetic analysis. The result indicated that Comp#5 enhanced the interaction of residues THR357, ARG362, LYS366, PRO424, CYS459, SER460, ALA461, ILE463, ARG465, THR507 and GLN510 with the surrounding residues, improving the stability of the catalytic active region and the entrance of catalytic active region. In particular, the Comp#5 conjugated with residue ARG362, elevating the efficient and selectivity of SHP2 protein. The study here may pave the way for discovering the novel SHP2 inhibitors for suffering cancer patients.

Evidence for Disruption of Mg2+ Pair as a Resistance Mechanism Against HIV-1 Integrase Strand Transfer Inhibitors

HIV-1 integrase is the enzyme responsible for integrating the viral DNA into the host genome and is one of the main targets for antiretroviral therapy; however, there are documented cases of resistance against all the currently used integrase strand transfer inhibitors (INSTIs). While some resistance-related mutations occur near the inhibitor’s binding site, the mutation N155H occurs on the opposite side of the drug-interacting Mg²⁺ ions, thus, not interacting directly with the drug molecules and currently lacking an explanation for its resistance mechanism. Moreover, mutation N155H and the resistance-related mutation Q148H are mutually exclusive for unknown reasons. In the present study, we use molecular dynamics simulations to understand the impact of the N155H mutation in the HIV-1 integrase structure and dynamics, when alone or in combination with Q148H. Our findings suggest that the Mg²⁺ ions of the active site adopt different orientations in each of the mutants, causing the catalytic triad residues involved in the ion coordination to adapt their side-chain configurations, completely changing the INSTIs binding site. The change in the ion coordination also seems to affect the flexibility of the terminal viral DNA nucleotide near the active site, potentially impairing the induced-fit mechanism of the drugs. The explanations obtained from our simulations corroborate previous hypotheses drawn from crystallographic studies. The proposed resistance mechanism can also explain the resistance caused by other mutations that take place in the same region of the integrase and help uncover the structural details of other HIV-1 resistance mechanisms.

Molecular dynamics simulation studies of the wild type and E92Q/N155H mutant of Elvitegravir-resistance HIV-1 integrase

Although Elvitegravir (EVG) is a newly developed antiretrovirals drug to treat the acquired immunodeficiency syndrome (AIDS), drug resistance has already been found in clinic, such as E92Q/N155H and Q148H/G140S. Several structural investigations have already been reported to reveal the molecular mechanism of the drug resistance. As full length crystal structure for HIV-1 integrase is still unsolved, we herein use the crystal structure of the full length prototype foamy virus (PFV) in complex with virus DNA and inhibitor Elvitegravir as a template to construct the wild type and E92Q/N155H mutant system of HIV-1 integrase. Molecular dynamic simulations was used to revel the binding mode and the drug resistance of the EVG ligand in E92Q/N155H. Several important interactions were discovered between the mutated residues and the residues in the active site of the E92Q/N155H double mutant pattern, and cross correlation and clustering methods were used for detailed analysis. The results from the MD simulation studies will be used to guide the experimental efforts of developing novel inhibitors against drug-resistant HIV integrase mutants.

Molecular dynamics simulation studies of the wild type and E92Q/N155H mutant of Elvitegravir-resistance HIV-1 integrase

Although Elvitegravir (EVG) is a newly developed antiretrovirals drug to treat the acquired immunodeficiency syndrome (AIDS), drug resistance has already been found in clinic, such as E92Q/N155H and Q148H/G140S. Several structural investigations have already been reported to reveal the molecular mechanism of the drug resistance. As full length crystal structure for HIV-1 integrase is still unsolved, we herein use the crystal structure of the full length prototype foamy virus (PFV) in complex with virus DNA and inhibitor Elvitegravir as a template to construct the wild type and E92Q/N155H mutant system of HIV-1 integrase. Molecular dynamic simulations was used to revel the binding mode and the drug resistance of the EVG ligand in E92Q/N155H. Several important interactions were discovered between the mutated residues and the residues in the active site of the E92Q/N155H double mutant pattern, and cross correlation and clustering methods were used for detailed analysis. The results from the MD simulation studies will be used to guide the experimental efforts of developing novel inhibitors against drug-resistant HIV integrase mutants.

Impact of resistance mutations on inhibitor binding to HIV-1 integrase

HIV-1 integrase (IN) is essential for HIV-1 replication, catalyzing two key reaction steps termed 3' processing and strand transfer. Therefore, IN has become an important target for antiviral drug discovery. However, mutants have emerged, such as E92Q/N155H and G140S/Q148H, which confer resistance to raltegravir (RAL), the first IN strand transfer inhibitor (INSTI) approved by the FDA, and to the recently approved elvitegravir (EVG). To gain insights into the molecular mechanisms of ligand binding and drug resistance, we performed molecular dynamics (MD) simulations of homology models of the HIV-1 IN and four relevant mutants complexed with viral DNA and RAL. The results show that the structure and dynamics of the 140s' loop, comprising residues 140 to 149, are strongly influenced by the IN mutations. In the simulation of the G140S/Q148H double mutant, we observe spontaneous dissociation of RAL from the active site, followed by an intrahelical swing-back of the 3'-OH group of nucleotide A17, consistent with the experimental observation that the G140S/Q148H mutant exhibits the highest resistance to RAL compared to other IN mutants. An important hydrogen bond between residues 145 and 148 is present in the wild-type IN but not in the G140S/Q148H mutant, accounting for the structural and dynamical differences of the 140s' loop and ultimately impairing RAL binding in the double mutant. End-point free energy calculations that broadly capture the experimentally known RAL binding profiles elucidate the contributions of the 140s' loop to RAL binding free energies and suggest possible approaches to overcoming drug resistance.

Molecular dynamics of HIV1-integrase in complex with 93del - a structural perspective on the mechanism of inhibition

HIV1 integrase is an important target for the antiviral therapy. Guanine-rich quadruplex, such as 93del, have been shown to be potent inhibitors of this enzyme and thus representing a new class of antiviral agents. Although X-ray and NMR structures of HIV1 integrase and 93del have been reported, there is no structural information of the complex and the mechanism of inhibition still remains unexplored. A number of computational methods including automated protein-DNA docking and molecular dynamics simulation in explicit solvent were used to model the binding of 93del to HIV1 integrase. Analysis of the dynamic behaviour of the complex using principal components analysis and elastic network modelling techniques allow us to understand how the binding of 93del aptamer and its interactions with key residues affect the intrinsic motions of the catalytic loops by stabilising them in catalytically inactive conformations. Such insights into the structural mechanism of inhibition can aid in improving the design of anti-HIV aptamers.

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.

Deconstructing honeybee vitellogenin: novel 40 kDa fragment assigned to its N terminus

Vitellogenin, an egg-yolk protein precursor common to oviparous animals, is found abundantly in honeybee workers - a caste of helpers that do not usually lay eggs. Instead, honeybee vitellogenin (180 kDa) participates in processes other than reproduction: it influences hormone signaling, food-related behavior, immunity, stress resistance and longevity. The molecular basis of these functions is largely unknown. Here, we establish and compare the molecular properties of vitellogenin from honeybee hemolymph (blood) and abdominal fat body, two compartments that are linked to vitellogenin functions. Our results reveal a novel 40 kDa vitellogenin fragment in abdominal fat body tissue, the main site for vitellogenin synthesis and storage. Using MALDI-TOF combined with MS/MS mass-spectroscopy, we assign the 40 kDa fragment to the N terminus of vitellogenin, whereas a previously observed 150 kDa fragment corresponded to the remainder of the protein. We show that both protein units are N glycosylated and phosphorylated. Focusing on the novel 40 kDa fragment, we present a homology model based on the structure of lamprey lipovitellin that includes a conserved β-barrel-like shape, with a lipophilic cavity in the interior and two insect-specific loops that have not been described before. Our data indicate that the honeybee fat body vitellogenin experiences cleavage unlike hemolymph vitellogenin, a pattern that can suggest a tissue-specific role. Our experiments advance the molecular understanding of vitellogenin, of which the multiple physiological and behavioral effects in honeybees are well established.

Conformational dynamics of the flexible catalytic loop in Mycobacterium tuberculosis 1-deoxy-D-xylulose 5-phosphate reductoisomerase

In mycobacteria, the biosynthesis of the precursors to the essential isoprenoids, isopentenyl diphosphate and dimethylallyl pyrophosphate is carried out by the methylerythritol phosphate pathway. This route of synthesis is absent in humans, who utilize the alternative mevalonate acid route, thus making the enzymes of the methylerythritol phosphate pathway of chemotherapeutic interest. One such identified target is the second enzyme of the pathway, 1-deoxy-D-xylulose 5-phosphate reductoisomerase. Only limited information is currently available concerning the catalytic mechanism and structural dynamics of this enzyme, and only recently has a crystal structure of Mycobacterium tuberculosis species of this enzyme been resolved including all factors required for binding. Here, the dynamics of the enzyme is studied in complex with NADPH, Mn2+, in the presence and absence of the fosmidomycin inhibitor using conventional molecular dynamics and an enhanced sampling technique, reversible digitally filtered molecular dynamics. The simulations reveal significant differences in the conformational dynamics of the vital catalytic loop between the inhibitor-free and inhibitor-bound enzyme complexes and highlight the contributions of conserved residues in this region. The substantial fluctuations observed suggest that 1-deoxy-D-xylulose 5-phosphate reductoisomerase may be a promising target for computer-aided drug discovery through the relaxed complex method.