Identification of Novel Influenza Polymerase PB2 Inhibitors Using a Cascade Docking Virtual Screening Approach

To discover novel inhibitors that target the influenza polymerase basic protein 2 (PB2) cap-binding domain (CBD), commercial ChemBridge compound libraries containing 384,796 compounds were screened using a cascade docking of LibDock–LigandFit–GOLD, and 60 compounds were selected for testing with cytopathic effect (CPE) inhibition assays and surface plasmon resonance (SPR) assay. Ten compounds were identified to rescue cells from H1N1 virus-mediated death at non-cytotoxic concentrations with EC₅₀ values ranging from 0.30 to 67.65 μM and could bind to the PB2 CBD of H1N1 with Kd values ranging from 0.21 to 6.77 μM. Among these, four compounds (11D4, 12C5, 21A5, and 21B1) showed inhibition of a broad spectrum of influenza virus strains, including oseltamivir-resistant ones, the PR/8-R292K mutant (H1N1, recombinant oseltamivir-resistant strain), the PR/8-I38T mutant (H1N1, recombinant baloxavir-resistant strain), and the influenza B/Lee/40 virus strain. These compounds have novel chemical scaffolds and relatively small molecular weights and are suitable for optimization as lead compounds. Based on sequence and structure comparisons of PB2 CBDs of various influenza virus subtypes, we propose that the Phe323/Gln325, Asn429/Ser431, and Arg355/Gly357 mutations, particularly the Arg355/Gly357 mutation, have a marked impact on the selectivities of PB2 CBD-targeted inhibitors of influenza A and influenza B.     

Temperature-gradient gel electrophoresis for the detection of polymorphic DNA and for quantitative polymerase chain reaction.

In order to detect mutations in a gene, either known mutations from human diseases or artificial ones in transgenic animals, or to screen for not yet identified mutations in patients, a method is required which guarantees detection of mutations which might occur in every single position of the whole open reading frame (ORF). It will be shown that a combination of polymerase chain reaction (PCR) and temperature gradient gel electrophoresis (TOGE) fulfills these requirements. By thermodynamic calculations the shift in the gel electrophoresis due to a mutation can be calculated in dependence on the position of the mutation. The theoretical results were tested with the mutations known so far. The quantitative determination of the copy number of a specific DNA or RNA sequence in a biological specimen (quantitative PCR) can be performed precisely and easily by combining PCR and TGGE. The system uses a quantification strategy of a new type of internal standardization. TGGE is applied to separate homo- and heteroduplexes which correspond respectively to standard and template sequences. The accuracy of this quantification strategy is very high, with a variability of < 15%. In addition to quantification, PCR/TGGE detects PCR artifacts and template mutants.     

Rapid mutation detection in complex genes by heteroduplex analysis with capillary array electrophoresis

Mutational analysis of large multiexon genes without prevalent mutations is a laborious undertaking that requires the use of a high-throughput scanning technique. The Human Genome Project has enabled the development of powerful techniques for mutation detection in large multiexon genes. We have transferred heteroduplex analysis (HA) by conformation-sensitive gel electrophoresis of the two major breast cancer (BC) predisposing genes, BRCA1 and BRCA2, to a multicapillary DNA sequencer in order to increase the throughput of this technique. This new method that we have called heteroduplex analysis by capillary array electrophoresis (HA-CAE) is based on the use of multiplex-polymerase chain reaction (PCR), different fluorescent labels and HA in a 16-capillary DNA sequencer. To date, a total of 114 different DNA sequence variants (19 insertions/deletions and 95 single-nucleotide substitutions - SNS) of BRCA1 and BRCA2 from 431 unrelated BC families have been successfully detected by HA-CAE. In addition, we have optimized the multiplex-PCR conditions for the colorectal cancer genes MLH1 and MSH2 in order to analyze them by HA-CAE. Both genes have been amplified in 13 multiplex groups, which contain the 35 exons, and their corresponding flanking intronic sequences. MLH1 and MSH2 have been analyzed in nine hereditary nonpolyposis colorectal cancer patients, and we have found six different DNA changes: one complex deletion/insertion mutation in MLH1 exon 19 and another five SNS. Only the complex mutation and one SNS may be classified as cancer-prone mutations. Our experience has revealed that HA-CAE is a simple, fast, reproducible and sensitive method to scan the sequences of complex genes.     

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.     

Molecular Dynamics and MM-PBSA Analysis of the SARS-CoV-2 Gamma Variant in Complex with the hACE-2 Receptor

The SARS-CoV-2 virus, since its appearance in 2019, has caused millions of cases and deaths. To date, there is no effective treatment or a vaccine that is fully protective. Despite the efforts made by governments and health institutions around the globe to control its propagation, the evolution of the virus has accelerated, diverging into hundreds of variants. However, not all of them are variants of concern (VoC’s). VoC’s have appeared in different regions and throughout the two years of the pandemic they have spread around the world. Specifically, in South America, the gamma variant (previously known as P.1) appeared in early 2021, bringing with it a second wave of infections. This variant contains the N501Y, E484K and K417T mutations in the receptor binding domain (RBD) of the spike protein. Although these mutations have been described experimentally, there is still no clarity regarding their role in the stabilization of the complex with the human angiotensin converting enzyme 2 (hACE-2) receptor. In this article we dissect the influence of mutations on the interaction with the hACE-2 receptor using molecular dynamics and estimations of binding affinity through a screened version of the molecular mechanics Poisson Boltzmann surface area (MM-PBSA) and interaction entropy. Our results indicate that mutations E484K and K417T compensate each other in terms of binding affinity, while the mutation N501Y promotes a more convoluted effect. This effect consists in the adoption of a cis configuration in the backbone of residue Y495 within the RBD, which in turn promotes polar interactions with the hACE-2 receptor. These results not only correlate with experimental observations and complement previous knowledge, but also expose new features associated with the specific contribution of concerned mutations. Additionally, we propose a recipe to assess the residue-specific contribution to the interaction entropy.     

Drug Repurposing Identifies Inhibitors of Oseltamivir‐Resistant Influenza Viruses

The neuraminidase (NA) inhibitor, oseltamivir, is a widely used anti‐influenza drug. However, oseltamivir‐resistant H1N1 influenza viruses carrying the H275Y NA mutation spontaneously emerged as a result of natural genetic drift and drug treatment. Because H275Y and other potential mutations may generate a future pandemic influenza strain that is oseltamivir‐resistant, alternative therapy options are needed. Herein, we show that a structure‐based computational method can be used to identify existing drugs that inhibit resistant viruses, thereby providing a first line of pharmaceutical defense against this possible scenario. We identified two drugs, nalidixic acid and dorzolamide, that potently inhibit the NA activity of oseltamivir‐resistant H1N1 viruses with the H275Y NA mutation at very low concentrations, but have no effect on wild‐type H1N1 NA even at a much higher concentration, suggesting that the oseltamivir‐resistance mutation itself caused susceptibility to these drugs.     

CRISPR-Mediated Profiling of Viral RNA at Single-Nucleotide Resolution

Mass pathogen screening is critical to preventing the outbreaks and spread of infectious diseases. The large-scale epidemic of COVID-19 and the rapid mutation of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus have put forward new requirements for virus detection and identification techniques. Here, we report a CRISPR-based Amplification-free Viral RNA Electrical Detection platform (CAVRED) for the rapid detection and identification of SARS-CoV-2 variants. A series of CRISPR RNA assays were designed to amplify the CRISPR-Cas system‘s ability to discriminate between mutant and wild RNA genomes with a single-nucleotide difference. The identified viral RNA information was converted into readable electrical signals through field-effect transistor biosensors for the achievement of highly sensitive detection of single-base mutations. CAVRED can detect the SARS-CoV-2 virus genome as low as 1 cp μL⁻¹ within 20 mins without amplification, and this value is comparable to the detection limit of real-time quantitative polymerase chain reaction. Based on the excellent RNA mutation detection ability, an 8-in-1 CAVRED array was constructed and realized the rapid identification of 40 simulated throat swab samples of SARS-CoV-2 variants with a 95.0 % accuracy. The advantages of accuracy, sensitivity, and fast speed of CAVRED promise its application in rapid and large-scale epidemic screening.     

Capillary electrophoretic analysis of wild type and mutant Proteus penneri outer membrane proteins

In this study the virulence factors, outer membrane proteins (OMP), lipopolysaccharides (LPS), hemolysin, and the in vivo and in vitro virulence of wild-type Proteus penneri 357 and its two isogenic mutant variants--a transposon and a spontaneous mutant--were examined. The OMPs of these variants were analyzed by a new and fast technique, "dynamic sieving" capillary electrophoresis (CE). The OMP profiles were dominated by two peaks (39 and 43 kDa). In the P. penneri clone examined, both the transposon and the spontaneous mutations induced significant changes in the OMP patterns (in the relative percentage of the dominant proteins). CE was suitable for the comparative analysis of bacterial protein patterns in the genetic variants of this strain, and provided valuable results in connection with the bacteriological virulence. The LPS composition of the genetic variants also showed alterations. The wild type of P. penneri 357 showed a typical ladder pattern, an "S" form, and the mutants possessed "R" LPS patterns (only few bands) in the gels. In the bacteriological virulence tests the wild type of P. penneri 357 was virulent in the in vivo, and toxic in the in vitro assays, while both mutants showed neither toxicity nor pathogenicity.     

Computational simulations of HIV-1 proteases--multi-drug resistance due to nonactive site mutation L90M

Human immunodeficiency virus type 1 protease (HIV-1 PR) is one of the proteins that currently available anti-HIV-1 drugs target. Inhibitors of HIV-1 PR have become available, and they have lowered the rate of mortality from acquired immune deficiency syndrome (AIDS) in advanced countries. However, the rate of emergence of drug-resistant HIV-1 variants is quite high because of their short retroviral life cycle and their high mutation rate. Serious drug-resistant mutations against HIV-1 PR inhibitors (PIs) frequently appear at the active site of PR. Exceptionally, some other mutations such as L90M cause drug resistance, although these appear at nonactive sites. The mechanism of resistance due to nonactive site mutations is difficult to explain. In this study, we carried out computational simulations of L90M PR in complex with each of three kinds of inhibitors and one typical substrate, and we clarified the mechanism of resistance. The L90M mutation causes changes in interaction between the side chain atoms of the 90th residue and the main chain atoms of the 25th residue, and a slight dislocation of the 25th residue causes rotation of the side chain at the 84th residue. The rotation of the 84th residue leads to displacement of the inhibitor from the appropriate binding location, resulting in a collision with the flap or loop region. The difference in levels of resistance to the three inhibitors has been explained from energetic and structural viewpoints, which provides the suggestion for promising drugs keeping its efficacy even for the L90M mutant.     

Efficient mutation detection in MEN1 gene using a combination of single-strand conformation polymorphism (MDGA) and heteroduplex analysis

For facilitated genotypic analysis of multiple endocrine neoplasia type 1 (MEN1), a familial syndrome associated with tumors of the parathyroid and neuroendocrine tissues, we developed two screening methods, heteroduplex mutation assay (HMA) and mutation detection gel analysis (MDGA), both based on electrophoretic discrimination of polymerase chain reaction (PCR) products, to detect the mutations. Forty-three genomic DNA samples were used for the evaluation of these techniques. The whole coding region of MEN1 was PCR-amplified with fluorescent primers and then denatured/renatured before electrophoresis on an automated sequencer. 100% of the mutations were detected, subsequently confirmed and identified by sequencing. "Negative" samples were used to evaluate the specificity and reproducibility of the two techniques. The combination of the two methods allows high throughput cost-effective mutation screening which is less laborious than systematic sequencing of the whole coding region of MEN1. Together, these methods provide an efficient screen for MEN1 mutations.     

Mechanisms of imidacloprid resistance in Nilaparvata lugens by molecular modeling

<正>Homology models of the ligand binding domain of the wild-type and Y151S mutant brown planthopper {Nilaparvata lugens)α1 and rat(Rattus norvegicus)β2 nicotinic acetylcholine receptor(nAChR) subunits were generated based on the crystal structure of acetylcholine binding protein of Lymnaea stagnalis.Neonicotinoid insecticide imidacloprid was docked into the putative binding site of wild-type and mutantα1β2 dimeric receptors by Surflex-docking,and the calculated docking energies were in agreement with experimental results.The resistance mechanisms and corresponding binding modes of imidacloprid on nAChRs containing the Y151S target-site mutation were discussed.     

BRCA1 mutation screening using restriction endonuclease fingerprinting-single-strand conformation polymorphism in an automated capillary electrophoresis system

Efficient mutation scanning techniques are needed for the rapid detection of novel disease-associated mutations and rare-sequence variants of putative importance. The large size of the breast cancer 1 gene (BRCA1) and the many mutations found throughout its entire coding sequence make screening for mutations in this gene particularly challenging. We have developed a method for screening exon 11 of the BRCA1 gene based on restriction enzyme digestion of fluorescence-labeled polymerase chain reaction (PCR) products followed by single-strand conformation polymorphism (SSCP) using an automated capillary electrophoresis system, denoted capillary restriction endonuclease fingerprinting (REF)-SSCP electrophoresis. Using this strategy on a control set of samples, we were able to detect 17 of 18 known sequence alterations. The method was then applied to screen 73 Norwegian females with family histories of breast and/or ovarian cancer. A total of 172 sequence alterations were detected, including substitutions, insertions, and deletions. One novel substitution of unknown function was identified. Sequencing of all samples negative in the capillary REF-SSCP system gave no additional mutations confirming the high sensitivity of the described methodology. Capillary REF-SSCP electrophoresis appeared as a technically convenient technique, requiring amplification of fewer PCR fragments than traditional SSCP. The novel strategy allows high-throughput mutation scanning without radioactive labeling and polyacrylamide gel electrophoresis (PAGE).     

Impact of the R292K Mutation on Influenza A (H7N9) Virus Resistance towards Peramivir: A Molecular Dynamics Perspective

In March 2013, a novel avian influenza A (H7N9) virus emerged in China. By March 2021, it had infected more than 1500 people, raising concerns regarding its epidemic potential. Similar to the highly pathogenic H5N1 virus, the H7N9 virus causes severe pneumonia and acute respiratory distress syndrome in most patients. Moreover, genetic analysis showed that this avian H7N9 virus carries human adaptation markers in the hemagglutinin and polymerase basic 2 (PB2) genes associated with cross-species transmissibility. Clinical studies showed that a single mutation, neuraminidase (NA) R292K (N2 numbering), induces resistance to peramivir in the highly pathogenic H7N9 influenza A viruses. Therefore, to evaluate the risk for human public health and understand the possible source of drug resistance, we assessed the impact of the NA-R292K mutation on avian H7N9 virus resistance towards peramivir using various molecular dynamics approaches. We observed that the single point mutation led to a distorted peramivir orientation in the enzyme active site which, in turn, perturbed the inhibitor’s binding. The R292K mutation induced a decrease in the interaction among neighboring amino acid residues when compared to its wild-type counterpart, as shown by the high degree of fluctuations in the radius of gyration. MM/GBSA calculations revealed that the mutation caused a decrease in the drug binding affinity by 17.28 kcal/mol when compared to the that for the wild-type enzyme. The mutation caused a distortion of hydrogen bond-mediated interactions with peramivir and increased the accessibility of water molecules around the K292 mutated residue.     

In silico analysis of the thermodynamic stability changes of psychrophilic and mesophilic alpha-amylases upon exhaustive single-site mutations

Identifying sequence modifications that distinguish psychrophilic from mesophilic proteins is important for designing enzymes with different thermodynamic stabilities and to understand the underlying mechanisms. The PoPMuSiC algorithm is used to introduce, in silico, all the single-site mutations in four mesophilic and one psychrophilic chloride-dependent alpha-amylases and to evaluate the changes in thermodynamic stability. The analysis of the distribution of the sequence positions that could be stabilized upon mutation shows a clear difference between the three domains of psychrophilic and mesophilic alpha-amylases. Most of the mutations stabilizing the psychrophilic enzyme are found in domains B and C, contrary to the mesophilic proteins where they are preferentially situated in the catalytic domain A. Moreover, the calculations show that the environment of some residues responsible for the activity of the psychrophilic protein has evolved to reinforce favorable interactions with these residues. In the second part, these results are exploited to propose rationally designed mutations that are predicted to confer to the psychrophilic enzyme mesophilic-like thermodynamic properties. Interestingly, most of the mutations found in domain C strengthen the interactions with domain A, in agreement with suggestions made on the basis of structural analyses. Although this study focuses on single-site mutations, the thermodynamic effects of the recommended mutations should be additive if the mutated residues are not close in space.     

BRAF mutation detection and identification by cycling temperature capillary electrophoresis

BRAF mutations are found in many human tumors, namely melanomas ( approximately 70%) and colon carcinomas ( approximately 15%). This paper presents a method for identification of exon 15 BRAF mutations by denaturant capillary electrophoresis (CE), an analysis method that is sensitive, cost-effective (involving only polymerase chain reaction (PCR) and electrophoresis) and capable of high-throughput screening. In total, we found 21 (70%) out of 30 melanoma cell lines with BRAF mutations in exon 15: two of which were the p.Val600Asp (c.1799-800TG>AT) mutation, one cell line contained the p.Val600Arg (c.1798-99GT>AG) mutation, and 18 cell lines contained the p.Val600Glu (c.1799T>A) mutation. Of the nine cell lines that did not contain a BRAF mutation, five contained an NRAS mutation at exon 2, and no mutations were detected in NRAS exon 1. There was no overlap of NRAS and BRAF mutations in the same cell line. In addition, we looked at 221 colon biopsy samples and identified one further BRAF mutation, the p.Asp594Gly (c.1781A>G) mutation, in seven samples. The p.Val600Glu mutation was identified in 11 of the colon biopsy samples. Using the four mutations of BRAF exon 15, we then constructed a denaturing CE standard capable of distinguishing between each of the mutations; therefore, sequencing does not need to be performed to confirm the mutation. In conclusion, this sensitive, cost-effective mutation assay for BRAF (and RAS) will provide the opportunity to detect and determine mutations without the need to purify samples for sequencing. Future large-scale studies will provide the clinical usefulness of such mutations.     

In silico explanation for the causalities of deleterious RNF213 SNPs in Moyamoya disease and insulin resistance

Moyamoya disease (MMD), a cerebrovascular disorder caused by the RNF213 gene, is a cerebrovascular, neurological disorder leading to ischemic strokes. Our previous work suggested that RNF213 might be involved in the pro-inflammatory TNFα-mediated insulin-resistance pathway in adipocytes. Insulin resistance can lead to cerebrovascular diseases and ischemic strokes. Though p. R4810 K has been reported as the founder mutation for Asian population with this disease, there are several mutations continuously reported in clinical diagnosis. We are interested to know whether these mutations can modulate insulin resistance. Also, we are intended to understand the causalities of RNF213 and its associated mutations in MMD. For this, we have adopted a computational approach to characterize RNF213 and its naturally occurring SNPs. Clinically reported SNPs and the predicted SNPs were analyzed for their pathogenicity and effect on the biological function of the protein. To increase accuracy, this was performed through three different analysis software (PROVEAN, SIFT, and SNAP2). The mutations that were found to be deleterious in all the three platforms were further analyzed for their effect on the thermal stability of the protein through I-mutant and iStable. It was found that R4810 K and other mutations decreased the thermodynamic stability of the protein. Loss of function of RNF213 was suggested in some reports. Contrary to this, some studies reported a gain of function state due to the R4810K mutation. To understand this we have measured the ligand-binding ability of this mutated protein through COFACTOR and COACH. An increase in ligand binding is always related to the functional stability of a protein. We have observed that the R4810K mutation might increase the iron-binding efficiency of the amino acid residues. This increase in binding was further validated by analyzing the binding efficiencies by docking. Since RNF213 was previously reported as a target for Protein Tyrosine Phosphatase 1B (PTP1B), we have also analyzed whether PTP1B-binding positions are susceptible to mutations. We have re-analyzed our earlier report on the differential expression pattern of RNF213 in cancer and obese samples. We have provided a detailed analysis of the most deleterious SNPs related to RNF213. Also, we provide a prediction for the loss of function and gain of function attributes of RNF213 and its predicted causalities in MMD and insulin resistance.     

Use of denaturing gradient gel electrophoresis to detect mutation in VS2 of the 16S-23S rDNA spacer amplified from Staphylococcus aureus isolates

To develop a double gradient denaturing gradient gel electrophoresis (DG-DGGE) based typing method that rapidly and accurately types clinical isolates of Staphylococcus aureus, the VS2 region of the 16S-23S rRNA spacer region (ISR) was chosen because of its potential high variation. The VS2 region was amplified with a 40-mer GC-clamp attached to the 5'-end of the reverse primer. The 145 bp PCR product was then separated by DG-DGGE using denaturant concentrations of 25-40% and polyacrylamide concentrations of 6-12%. Of the five mutations identified in 336 S. aureus isolates, one mutation was found to be highly specific for 161/171 (94%) of methicillin-resistant S. aureus (MRSA) isolates from different geographic locations and isolation times. This same mutation was found in 15/160 (9%) of penicillin- or methicillin-sensitive S. aureus isolates. In some isolates two mutations occured together in the one genome suggesting some S. aureus isolates have two copies of VS2. In these 336 isolates nine genotypes with different combinations of the five mutations were identified. In 18 coagulase-negative staphylococci (CNS), the MRSA-specific mutation was found along with two other mutations in all isolates demonstrating consistent differences in the presence of these mutations between CNS and S. aureus. The marked differences in VS2 sequences found between MRSA, methicillin- or penicillin-sensitive S. aureus (SSA), and CNS by DGGE in the present study may be useful in evolutionary studies and in the development of a specific assay for MRSA from clinical specimens.     

Global and Local Conformation of Human IgG Antibody Variants Rationalizes Loss of Thermodynamic Stability

Immunoglobulin G (IgG) monoclonal antibodies (mAbs) are a major class of medicines, with high specificity and affinity towards targets spanning many disease areas. The antibody Fc (fragment crystallizable) region is a vital component of existing antibody therapeutics, as well as many next generation biologic medicines. Thermodynamic stability is a critical property for the development of stable and effective therapeutic proteins. Herein, a combination of ion‐mobility mass spectrometry (IM‐MS) and hydrogen/deuterium exchange mass spectrometry (HDX‐MS) approaches have been used to inform on the global and local conformation and dynamics of engineered IgG Fc variants with reduced thermodynamic stability. The changes in conformation and dynamics have been correlated with their thermodynamic stability to better understand the destabilising effect of functional IgG Fc mutations and to inform engineering of future therapeutic proteins.     

Direct binding assay for the detection of type IV allosteric inhibitors of Abl

Abelson (Abl) tyrosine kinase is an important cellular enzyme that is rendered constitutively active in the breakpoint cluster region (BCR)-Abl fusion protein, contributing to several forms of leukemia. Although inhibiting BCR-Abl activity with imatinib shows great clinical success, many patients acquire secondary mutations that result in resistance to imatinib. Second-generation inhibitors such as dasatinib and nilotinib can overcome the majority of these mutations but fail to treat patients with an especially prevalent T315I mutation at the gatekeeper position of the kinase domain. However, a combination of nilotinib with an allosteric type IV inhibitor was recently shown to overcome this clinically relevant point mutation. In this study, we present the development of a direct binding assay that enables the straightforward detection of allosteric inhibitors which bind within the myristate pocket of Abl. The assay is amenable to high-throughput screening and exclusively detects the binding of ligands to this unique allosteric site.     

Aromatic interaction profile to understand the molecular basis of raltegravir resistance

Integrase (IN) is the enzyme of human immunodeficiency virus (HIV) which inserts the viral DNA (vDNA) into the host genome for successful viral replication leading to the infection. However, the chemical basis of HIV IN catalysis is speculative due to lack of complete co-crystal structure. Using the recently published prototype foamy virus IN crystal structure, we developed a model structure of HIV IN showing interaction of vDNA, the metal (Mg²⁺) cofactor, and raltegravir (RLT) in the active site. Molecular docking and dynamics simulations studies showed that RLT uses it core central ring with diketo motif for Mg²⁺ chelation and bridge interaction with DDE motif. The triple arene interactions mediated by RLT with neighboring molecular motifs (Y143, cytosine, and adenine) is maintained during long simulation in wild type (WT). The fluorobenzyl and oxadiazole moieties of RLT forms aromatic stacking with cytosine base (head stacking) aromatic side chain of Y143 (tail stacking), respectively, while central ring further establishes aromatic stacking with distorted adenine base of vDNA (central stacking). The novel triple stacking systems were further explored to understand the molecular basis of drug resistance by molecular simulation. The in silico mutation (N155H, Q148H, and Q148H + G140S) and simulation studies elucidated the structural mechanism of resistance to RLT. The simulation studies provided the molecular basis for interdependency observed for the primary and secondary (Q148H and G140S) mutations and also explained the mechanism of viral fitness regain. Our study reveals that triple stacking and its consequence in terms of VdW energetic profile acts as a critical point to understand the drug-resistance. Here, we demonstrate that the root mean square deviation of centroid system (aromatic stacking) can be used as a major determinant of RLT binding toward the fold resistance. This is first kind of report, which discloses a strategy to explore the molecular level of drug resistance profile using aromatic interactions.