Site Binding Mode Research Articles

Probing the pH-dependency of DC-SIGN/R multivalent lectin-glycan interactions using polyvalent glycan-gold nanoparticles.

The dendritic cell tetrameric lectin, DC-SIGN, and its closely related endothelial cell lectin, DC-SIGNR (collectively abbreviated as DC-SIGN/R) play a key role in the binding and transmission of deadly viruses, including Ebola, HIV, HCV, and SARS-CoV-2. Their virus binding/release processes involve a gradually acidifying environment following the natural intracellular trafficking pathways. Therefore, understanding DC-SIGN/R's pH-dependent binding properties with glycan ligands is of great importance. We have recently developed densely glycosylated gold nanoparticles (glycan-GNPs) as a powerful new tool for probing DC-SIGN/R multivalent lectin-glycan interaction (MLGI) mechanisms. They can provide not only quantitative MLGI affinities but also important structural information, such as binding site orientation and binding modes. Herein, we further employ the glycan-GNP probes to investigate the pH dependency of DC-SIGN/R MLGI properties. We find that DC-SIGN/R MLGIs exhibit distinct pH dependence over the normal physiological (7.4) to lysosomal (∼4.6) pH range. DC-SIGN binds glycan-GNPs strongly and stably from pH 7.4 to ∼5.8, but the binding is weakened significantly as pH decreases to ≤5.4 and may be fully dissociated at pH 4.6. This behaviour is fully consistent with DC-SIGN's role as an endocytic recycling receptor. In contrast, DC-SIGNR's affinity with glycan-GNPs is enhanced with the decreasing pH from 7.4 to 5.4, peaking at pH 5.4, and then reduced as pH is further lowered. Interestingly, both DC-SIGN/R binding with glycan-GNPs are found to be partially reversible in a pH-dependent manner.

Probing scaffold size effects on multivalent lectin-glycan binding affinity, thermodynamics and antiviral properties using polyvalent glycan-gold nanoparticles.

Multivalent lectin-glycan interactions (MLGIs) are pivotal for viral infections and immune regulation. Their structural and biophysical data are thus highly valuable, not only for understanding their basic mechanisms but also for designing potent glycoconjugate therapeutics against target MLGIs. However, such information for some important MGLIs remains poorly understood, greatly limiting research progress. We have recently developed densely glycosylated nanoparticles, e.g., ∼4 nm quantum dots (QDs) or ∼5 nm gold nanoparticles (GNPs), as mechanistic probes for MLGIs. Using two important model lectin viral receptors, DC-SIGN and DC-SIGNR, we have shown that these probes can not only offer sensitive fluorescence assays for quantifying MLGI affinities, but also reveal key structural information (e.g., binding site orientation and binding mode) useful for MLGI targeting. However, the small sizes of the previous scaffolds may not be optimal for maximising MLGI affinity and targeting specificity. Herein, using α-manno-α-1,2-biose (DiMan) functionalised GNP (GNP-DiMan) probes, we have systematically studied how GNP scaffold size (e.g., 5, 13, and 27 nm) and glycan density (e.g., 100, 75, 50 and 25%) determine their MLGI affinities, thermodynamics, and antiviral properties. We have developed a new GNP fluorescence quenching assay format to minimise the possible interference of GNP's strong inner filter effect in MLGI affinity quantification, revealing that increasing the GNP size is highly beneficial for enhancing MLGI affinity. We have further determined the MLGI thermodynamics by combining temperature-dependent affinity and Van't Hoff analyses, revealing that GNP-DiMan-DC-SIGN/R binding is enthalpy driven with favourable binding Gibbs free energy changes (ΔG°) being enhanced with increasing GNP size. Finally, we show that increasing the GNP size significantly enhances their antiviral potency. Notably, the DiMan coated 27 nm GNP potently and robustly blocks both DC-SIGN and DC-SIGNR mediated pseudo-Ebola virus cellular entry with an EC50 of ∼23 and ∼49 pM, respectively, making it the most potent glycoconjugate inhibitor against DC-SIGN/R-mediated Ebola cellular infections. Our results have established GNP-glycans as a new tool for quantifying MLGI biophysical parameters and revealed that increasing the GNP scaffold size significantly enhances their MLGI affinities and antiviral potencies.

The effects of nature-inspired amino acid substitutions on structural and biochemical properties of the E. coli L-asparaginase EcAIII.

The Escherichia coli enzyme EcAIII catalyzes the hydrolysis of L-Asn to L-Asp and ammonia. Using a nature-inspired mutagenesis approach, we designed and produced five new EcAIII variants (M200I, M200L, M200K, M200T, M200W). The modified proteins were characterized by spectroscopic and crystallographic methods. All new variants were enzymatically active, confirming that the applied mutagenesis procedure has been successful. The determined crystal structures revealed new conformational states of the EcAIII molecule carrying the M200W mutation and allowed a high-resolution observation of an acyl-enzyme intermediate with the M200L mutant. In addition, we performed structure prediction, substrate docking, and molecular dynamics simulations for 25 selected bacterial orthologs of EcAIII, to gain insights into how mutations at the M200 residue affect the active site and substrate binding mode. This comprehensive strategy, including both experimental and computational methods, can be used to guide further enzyme engineering and can be applied to the study of other proteins of medicinal or biotechnological importance.

Identification of potential selective autophagy receptors from protein-content profiling of autophagosomes.

Selective autophagy receptors (SARs) are central to cellular homeostatic and organellar recycling pathways. Over the last two decades, more than 30 SARs have been discovered and validated using a variety of experimental approaches ranging from cell biology to biochemistry, including high-throughput imaging and screening methods. Yet, the extent of selective autophagy pathways operating under various cellular contexts, for example, under basal and starvation conditions, remains unresolved. Currently, our knowledge of all known SARs and their associated cargo components is fragmentary and limited by experimental data with varying degrees of resolution. Here, we use classical predictive and modeling approaches to integrate high-quality autophagosome content profiling data with disparate datasets. We identify a global set of potential SARs and their associated cargo components active under basal autophagy, starvation-induced, and proteasome-inhibition conditions. We provide a detailed account of cellular components, biochemical pathways, and molecular processes that are degraded via autophagy. Our analysis yields a catalog of new potential SARs that satisfy the characteristics of bonafide, well-characterized SARs. We categorize them by the subcellular compartments they emerge from and classify them based on their likely mode of action. Our structural modeling validates a large subset of predicted interactions with the human ATG8 family of proteins and shows characteristic, conserved LC3-interacting region (LIR)-LIR docking site (LDS) and ubiquitin-interacting motif (UIM)-UIM docking site (UDS) binding modes. Our analysis also revealed the most abundant cargo molecules targeted by these new SARs. Our findings expand the repertoire of SARs and provide unprecedented details into the global autophagic state of HeLa cells. Taken together, our findings provide motivation for the design of new experiments, testing the role of these novel factors in selective autophagy.

Ligand Dilution Analysis Facilitates Aptamer Binding Characterization at the Single-Molecule Level.

Cell-specific aptamers offer a powerful tool to study membrane receptors at the single-molecule level. Most target receptors of aptamers are highly expressed on the cell surface, but difficult to analyze in situ because of dense distribution and fast velocity. Therefore, we herein propose a random sampling-based analysis strategy termed ligand dilution analysis (LDA) for easily implemented aptamer-based receptor study. Receptor density on the cell surface can be calculated based on a regression model. By using a synergistic ligand dilution design, colocalization and differentiation of aptamer and monoclonal antibody (mAb) binding on a single receptor can be realized. Once this is accomplished, precise binding site and detailed aptamer-receptor binding mode can be further determined using molecular docking and molecular dynamics simulation. The ligand dilution strategy also sets the stage for an aptamer-based dynamics analysis of two- and three-dimensional motion and fluctuation of highly expressed receptors on the live cell membrane.

Computational analysis of non-competitive antagonist arylguanidine-α7 nAChR complexes

meta-Chlorophenylguanidine (1) is a non-competitive α7 nicotinic acetylcholine receptor (nAChR) antagonist. Here we examined the hydrogen bond donor role of the anilinic N1–H on the inhibitory effect of 1 by preparing its N1–CH3 counterpart 2. Analog 2 was found to be at least as potent as 1 as a non-competitive α7 nAChR antagonist in a patch-clamp assay.To establish a structural basis for the mode of interaction of guanidines 1 and 2, we generated 100 homology models of the hα7 nAChR. This was followed by Connolly surface (SYBYL-X2.1) and blind docking (AutoDock 4.1) studies to identify eight possible binding pockets, two of which were supported by empirical data and employed in our docking studies. The optimized model-ligand complexes were analyzed using a Hydropathic INTeractions (HINT) analysis in order to compare and contrast different binding pockets and modes. We identified a potential allosteric binding site and distinct rotameric binding modes for 1 and 2 at α7 nAChRs. These differences in the binding orientations minimized the importance of an anilinic NH function for the antagonist activity at nACh receptors.

Fragment Binding Pose Predictions Using Unbiased Simulations and Markov-State Models.

Predicting the costructure of small-molecule ligands and their respective target proteins has been a long-standing problem in drug discovery. For weak binding compounds typically identified in fragment-based screening (FBS) campaigns, determination of the correct binding site and correct binding mode is usually done experimentally via X-ray crystallography. For many targets of pharmaceutical interest, however, establishing an X-ray system which allows for sufficient throughput to support a drug discovery project is not possible. In this case, exploration of fragment hits becomes a very laborious and consequently slow process with the generation of protein/ligand cocrystal structures as the bottleneck of the entire process. In this work, we introduce a computational method which is able to reliably predict binding sites and binding modes of fragment-like small molecules using solely the structure of the apoprotein and the ligand's chemical structure as input information. The method is based on molecular dynamics simulations and Markov-state models and can be run as a fully automated protocol requiring minimal human intervention. We describe the application of the method to a representative subset of different target classes and fragments from historical FBS efforts at Boehringer Ingelheim and discuss its potential integration into the overall fragment-based drug discovery workflow.

Crystal structures of the closed form of Mycobacterium tuberculosis dihydrofolate reductase in complex with dihydrofolate and antifolates.

Tuberculosis is a disease caused by Mycobacterium tuberculosis and is the leading cause of death from a single infectious pathogen, with a high prevalence in developing countries in Africa and Asia. There still is a need for the development or repurposing of novel therapies to combat this disease owing to the long-term nature of current therapies and because of the number of reported resistant strains. Here, structures of dihydrofolate reductase from M.tuberculosis (MtDHFR), which is a key target of the folate pathway, are reported in complex with four antifolates, pyrimethamine, cycloguanil, diaverdine and pemetrexed, and its substrate dihydrofolate in order to understand their binding modes. The structures of all of these complexes were obtained in the closed-conformation state of the enzyme and a fine structural analysis indicated motion in key regions of the substrate-binding site and different binding modes of the ligands. In addition, the affinities, through Kd measurement, of diaverdine and methotrexate have been determined; MtDHFR has a lower affinity (highest Kd) for diaverdine than pyrimethamine and trimethoprim, and a very high affinity for methotrexate, as expected. The structural comparisons and analysis described in this work provide new information about the plasticity of MtDHFR and the binding effects of different antifolates.

Predicting Site-Binding Modes of Ions and Water to Nucleic Acids Using Molecular Solvation Theory.

Site binding of ions and water shapes nucleic acids folding, dynamics, and biological function, complementing the more diffuse, nonspecific "territorial" ion binding. Unlike territorial binding, prediction of site-specific binding to nucleic acids remains an unsolved challenge in computational biophysics. This work presents a new toolset based on the 3D-RISM molecular solvation theory and topological analysis that predicts cation and water site binding to nucleic acids. 3D-RISM is shown to accurately capture alkali cations and water binding to the central channel, transversal loops, and grooves of the Oxytricha nova's telomeres' G-quadruplex ( Oxy-GQ), in agreement with high-resolution crystallographic data. To improve the computed cation occupancy along the Oxy-GQ central channel, it was necessary to refine and validate new cation-oxygen parameters using structural and thermodynamic data available for crown ethers and ion channels. This single set of parameters that describes both localized and delocalized binding to various biological systems is used to gain insight into cation occupancy along the Oxy-GQ channel under various salt conditions. The paper concludes with prospects for extending the method to predict divalent cation binding to nucleic acids. This work advances the forefront of theoretical methods able to provide predictive insight into ion atmosphere effects on nucleic acids function.

The synthesis and evaluation of phenoxyacylhydroxamic acids as potential agents for Helicobacter pylori infections

Two series of ω-phenoxy contained acylhydroxamic acids as novel urease inhibitors were designed and synthesized. Biological activity evaluations revealed that ω-phenoxypropinoylhydroxamic acids were more active than phenoxyacetohydroxamic acids. Out of these compounds, 3-(3,4-dichlorophenoxy)propionylhydroxamic acid c24 showed significant potency against urease in both cell free extract (IC50 = 0.061 ± 0.003 μM) and intact cell (IC50 = 0.89 ± 0.05 μM), being over 450- and 120-fold more potent than the clinically prescribed urease inhibitor AHA, repectively. Non-linear fitting of experimental data (V-[S]) suggested a mixed-type inhibition mechanism and a dual site binding mode of these compounds.

A High-Throughput Dose-Response Cellular Thermal Shift Assay for Rapid Screening of Drug Target Engagement in Living Cells, Exemplified Using SMYD3 and IDO1

A persistent problem in early small-molecule drug discovery is the frequent lack of rank-order correlation between biochemical potencies derived from initial screens using purified proteins and the diminished potency and efficacy observed in subsequent disease-relevant cellular phenotypic assays. The introduction of the cellular thermal shift assay (CETSA) has bridged this gap by enabling assessment of drug target engagement directly in live cells based on ligand-induced changes in protein thermal stability. Initial success in applying CETSA across multiple drug target classes motivated our investigation into replacing the low-throughput, manually intensive Western blot readout with a quantitative, automated higher-throughput assay that would provide sufficient capacity to use CETSA as a primary hit qualification strategy. We introduce a high-throughput dose-response cellular thermal shift assay (HTDR-CETSA), a single-pot homogenous assay adapted for high-density microtiter plate format. The assay features titratable BacMam expression of full-length target proteins fused to the DiscoverX 42 amino acid ePL tag in HeLa suspension cells, facilitating enzyme fragment complementation–based chemiluminescent quantification of ligand-stabilized soluble protein. This simplified format can accommodate determination of full-dose CETSA curves for hundreds of individual compounds/analyst/day in replicates. HTDR-CETSA data generated for substrate site and alternate binding mode inhibitors of the histone-lysine N-methyltransferase SMYD3 in HeLa suspension cells demonstrate excellent correlation with rank-order potencies observed in cellular mechanistic assays and direct translation to target engagement of endogenous Smyd3 in cancer-relevant cell lines. We envision this workflow to be generically applicable to HTDR-CETSA screening spanning a wide variety of soluble intracellular protein target classes.

Crystal structure of the Legionella pneumophila Lpg2936 in complex with the cofactor S-adenosyl-L-methionine reveals novel insights into the mechanism of RsmE family methyltransferases.

The methylation of U1498 located in the 16S ribosomal RNA of Escherichia coli is an important modification affecting ribosomal activity. RsmE methyltransferases methylate specifically this position in a mechanism that requires an S‐adenosyl‐L‐methionine (AdoMet) molecule as cofactor. Here we report the structure of Apo and AdoMet‐bound Lpg2936 from Legionella pneumophila at 1.5 and 2.3 Å, respectively. The protein comprises an N‐terminal PUA domain and a C‐terminal SPOUT domain. The latter is responsible for protein dimerization and cofactor binding. Comparison with similar structures suggests that Lpg2936 is an RsmE‐like enzyme that can target the equivalent of U1498 in the L. pneumophila ribosomal RNA, thereby potentially enhancing ribosomal activity during infection‐mediated effector production. The multiple copies of the enzyme found in both structures reveal a flexible conformation of the bound AdoMet ligand. Isothermal titration calorimetry measurements suggest an asymmetric two site binding mode. Our results therefore also provide unprecedented insights into AdoMet/RsmE interaction, furthering our understanding of the RsmE catalytic mechanism.

Design, synthesis and biological evaluations of N-Hydroxy thienopyrimidine-2,4-diones as inhibitors of HIV reverse transcriptase-associated RNase H

Human immunodeficiency virus (HIV) reverse transcriptase (RT) associated ribonuclease H (RNase H) is the only HIV enzymatic function not targeted by current antiviral drugs. Although various chemotypes have been reported to inhibit HIV RNase H, few have shown significant antiviral activities. We report herein the design, synthesis and biological evaluation of a novel N-hydroxy thienopyrimidine-2,3-dione chemotype (11) which potently and selectively inhibited RNase H with considerable potency against HIV-1 in cell culture. Current structure-activity-relationship (SAR) identified analogue 11d as a nanomolar inhibitor of RNase H (IC50 = 0.04 μM) with decent antiviral potency (EC50 = 7.4 μM) and no cytotoxicity (CC50 > 100 μM). In extended biochemical assays compound 11d did not inhibit RT polymerase (pol) while inhibiting integrase strand transfer (INST) with 53 fold lower potency (IC50 = 2.1 μM) than RNase H inhibition. Crystallographic and molecular modeling studies confirmed the RNase H active site binding mode.

Molecular and cellular mechanisms of HIF prolyl hydroxylase inhibitors in clinical trials.

Inhibition of the human 2-oxoglutarate (2OG) dependent hypoxia inducible factor (HIF) prolyl hydroxylases (human PHD1-3) causes upregulation of HIF, thus promoting erythropoiesis and is therefore of therapeutic interest. We describe cellular, biophysical, and biochemical studies comparing four PHD inhibitors currently in clinical trials for anaemia treatment, that describe their mechanisms of action, potency against isolated enzymes and in cells, and selectivities versus representatives of other human 2OG oxygenase subfamilies. The 'clinical' PHD inhibitors are potent inhibitors of PHD catalyzed hydroxylation of the HIF-α oxygen dependent degradation domains (ODDs), and selective against most, but not all, representatives of other human 2OG dependent dioxygenase subfamilies. Crystallographic and NMR studies provide insights into the different active site binding modes of the inhibitors. Cell-based results reveal the inhibitors have similar effects on the upregulation of HIF target genes, but differ in the kinetics of their effects and in extent of inhibition of hydroxylation of the N- and C-terminal ODDs; the latter differences correlate with the biophysical observations.

A dual substrate kinetic model for cytochrome P450BM3-F87G catalysis: simultaneous binding of long chain aldehydes and 4-fluorophenol.

To develop a model for binding and catalysis associated with the stimulation of 4-fluorophenol (4-FP) oxidation in the presence of long chain aldehydes by the enzymatic catalyst, cytochrome P450BM3-F87G. A variation of the Michaeli-Menten kinetic model was employed to describe interactions at the active site of the enzyme, along with computer aided modeling approaches. In addition to the hydroquinone product arising from de-fluorination of 4-FP, a second product (p-fluorocatechol) was also observed and, like the hydroquinone, its rate of formation increased in the presence of the aldehyde. When only aldehyde was present with the enzyme, BM3-F87G catalyzed its oxidation to the corresponding carboxylic acid; however, this activity was inhibited when 4-FP was added to the reaction. A 3D computer model of the active site containing both aldehyde and 4-FP was generated, guided by these kinetic observations. Finally, partitioning between the two phenolic products was examined with an emphasis on the conditions directing the initial epoxidation at either the 2,3- or 3,4-positions on the substrate. Temperature, reaction time, substrate concentration, and the structure of the aldehyde had no substantial effect on the overall product ratios, however the NADPH coupling efficiency decreased when unsaturated aldehydes were included, or when the temperature of the reaction was reduced. The unsaturated aldehyde, trans-2-decenal, stimulates BM3-F87G catalyzed oxidation of 4-fluorophenol through a cooperative active site binding mode that doesn't influence product distributions or coupling efficiencies, while 4-fluorophenol acts as a competitive inhibitor of aldehyde oxidation.

Development, synthesis and structure–activity-relationships of inhibitors of the macrophage infectivity potentiator (Mip) proteins of Legionella pneumophila and Burkholderia pseudomallei

The bacteria Burkholderia pseudomallei and Legionella pneumophila cause severe diseases like melioidosis and Legionnaire’s disease with high mortality rates despite antibiotic treatment. Due to increasing antibiotic resistances against these and other Gram-negative bacteria, alternative therapeutical strategies are in urgent demand. As a virulence factor, the macrophage infectivity potentiator (Mip) protein constitutes an attractive target. The Mip proteins of B. pseudomallei and L. pneumophila exhibit peptidyl-prolyl cis/trans isomerase (PPIase) activity and belong to the PPIase superfamily. In previous studies, the pipecolic acid moiety proved to be a valuable scaffold for inhibiting this PPIase activity. Thus, a library of pipecolic acid derivatives was established guided by structural information and computational analyses of the binding site and possible binding modes. Stability and toxicity considerations were taken into account in iterative extensions of the library. Synthesis and evaluation of the compounds in PPIase assays resulted in highly active inhibitors. The activities can be interpreted in terms of a common binding mode obtained by docking calculations.

Active site binding modes of inhibitors of Staphylococcus aureus mevalonate diphosphate decarboxylase from docking and molecular dynamics simulations

Bacterial mevalonate diphosphate decarboxylase (MDD) is an attractive therapeutic target for antibacterial drug development. In this work, we discuss a combined docking and molecular dynamics strategy toward inhibitor binding to bacterial MDD. The docking parameters utilized in this study were first validated with observations for the inhibitors 6-fluoromevalonate diphosphate (FMVAPP) and diphosphoglycolylproline (DPGP) using existing structures for the Staphylococcus epidermidis enzyme. The validated docking protocol was then used to predict structures of the inhibitors bound to Staphylococcus aureus MDD using the unliganded crystal structure of Staphylococcus aureus MDD. We also investigated a possible interactions improvement by combining this docking method with molecular dynamics simulations. Thus, the predicted docking structures were analyzed in a molecular dynamics trajectory to generate dynamic models and reinforce the predicted binding modes. FMVAPP is predicted to make more extensive contacts with S. aureus MDD, forming stable hydrogen bonds with Arg144, Arg193, Lys21, Ser107, and Tyr18, as well as making stable hydrophobic interactions with Tyr18, Trp19, and Met196. The differences in predicted binding are supported by experimentally determined Ki values of 0.23 ± 0.02 and 34 ± 8 μM, for FMVAPP and DPGP, respectively. The structural information coupled with the kinetic characterization obtained from this study should be useful in defining the requirements for inhibition as well as in guiding the selection of active compounds for inhibitor optimization.

Hydrophobic and charged residues in the C-terminal arm of hepatitis C virus RNA-dependent RNA polymerase regulate initiation and elongation.

The RNA-dependent RNA polymerase (RdRp) of hepatitis C virus (HCV) is essential for viral genome replication. Crystal structures of the HCV RdRp reveal two C-terminal features, a β-loop and a C-terminal arm, suitably located for involvement in positioning components of the initiation complex. Here we show that these two elements intimately regulate template and nucleotide binding, initiation, and elongation. We constructed a series of β-loop and C-terminal arm mutants, which were used for in vitro analysis of RdRp de novo initiation and primer extension activities. All mutants showed a substantial decrease in initiation activities but a marked increase in primer extension activities, indicating an ability to form more stable elongation complexes with long primer-template RNAs. Structural studies of the mutants indicated that these enzyme properties might be attributed to an increased flexibility in the C-terminal features resulting in a more open polymerase cleft, which likely favors the elongation process but hampers the initiation steps. A UTP cocrystal structure of one mutant shows, in contrast to the wild-type protein, several alternate conformations of the substrate, confirming that even subtle changes in the C-terminal arm result in a more loosely organized active site and flexible binding modes of the nucleotide. We used a subgenomic replicon system to assess the effects of the same mutations on viral replication in cells. Even the subtlest mutations either severely impaired or completely abolished the ability of the replicon to replicate, further supporting the concept that the correct positioning of both the β-loop and C-terminal arm plays an essential role during initiation and in HCV replication in general. HCV RNA polymerase is a key target for the development of directly acting agents to cure HCV infections, which necessitates a thorough understanding of the functional roles of the various structural features of the RdRp. Here we show that even highly conservative changes, e.g., Tyr→Phe or Asp→Glu, in these seemingly peripheral structural features have profound effects on the initiation and elongation properties of the HCV polymerase.

Virtual screening with AutoDock Vina and the common pharmacophore engine of a low diversity library of fragments and hits against the three allosteric sites of HIV integrase: participation in the SAMPL4 protein–ligand binding challenge

To rigorously assess the tools and protocols that can be used to understand and predict macromolecular recognition, and to gain more structural insight into three newly discovered allosteric binding sites on a critical drug target involved in the treatment of HIV infections, the Olson and Levy labs collaborated on the SAMPL4 challenge. This computational blind challenge involved predicting protein-ligand binding against the three allosteric sites of HIV integrase (IN), a viral enzyme for which two drugs (that target the active site) have been approved by the FDA. Positive control cross-docking experiments were utilized to select 13 receptor models out of an initial ensemble of 41 different crystal structures of HIV IN. These 13 models of the targets were selected using our new "Rank Difference Ratio" metric. The first stage of SAMPL4 involved using virtual screens to identify 62 active, allosteric IN inhibitors out of a set of 321 compounds. The second stage involved predicting the binding site(s) and crystallographic binding mode(s) for 57 of these inhibitors. Our team submitted four entries for the first stage that utilized: (1) AutoDock Vina (AD Vina) plus visual inspection; (2) a new common pharmacophore engine; (3) BEDAM replica exchange free energy simulations, and a Consensus approach that combined the predictions of all three strategies. Even with the SAMPL4's very challenging compound library that displayed a significantly lower amount of structural diversity than most libraries that are conventionally employed in prospective virtual screens, these approaches produced hit rates of 24, 25, 34, and 27 %, respectively, on a set with 19 % declared binders. Our only entry for the second stage challenge was based on the results of AD Vina plus visual inspection, and it ranked third place overall according to several different metrics provided by the SAMPL4 organizers. The successful results displayed by these approaches highlight the utility of the computational structure-based drug discovery tools and strategies that are being developed to advance the goals of the newly created, multi-institution, NIH-funded center called the "HIV Interaction and Viral Evolution Center".

Prediction of the interaction between spermidine and the G-G mismatch containing acceptor stem in tRNAIle: Molecular modeling, density functional theory, and molecular dynamics study

Polyamines, putrescine, spermidine (SPD), and spermine are closely linked to cell growth, and highly regulate the levels of transcription, translation and protein turnover. We propose that SPD stimulates the formation of Ile-tRNAIle by inducing a selective structural change of the G-G mismatch containing acceptor stem in tRNAIle. Here, we provide insight into how SPD recognizes and stabilizes the G-G mismatch containing acceptor stem in tRNAIle with molecular modeling (MM), density functional theory (DFT) calculations, and molecular dynamics (MD) simulations. The results of the MM and DFT calculations indicate that the negatively charged region of the G-G mismatch containing acceptor stem in tRNAIle is preferentially recognized by positively charged SPD. In addition, MD simulations indicate that all of the positively charged amino groups of SPD under physiological conditions (N1(NH3+), N5(NH2+), and N10(NH3+) could form hydrogen bonds with tRNAIle and trigger the SPD-induced stabilization and structural change of the G-G mismatch containing acceptor stem in tRNAIle. Thus, this approach should be useful for determining the preferential binding site and appropriate binding mode of polyamines on tRNAIle.