On the detection of multiple-binding modes of ligands to proteins, from biological, structural, and modeling data

Revealing the graded activation mechanism of neurotensin receptor 1

Graded activation contributes to the precise regulation of GPCR activity, presenting new opportunities for drug design. In this work, a total of 10 μs enhanced-sampling simulations are performed to provide molecular insights into the binding dynamics differences of the neurotensin receptor 1 (NTSR1) to the full agonist SRI-9829, partial agonist RTI-3a and inverse agonist SR48692. The possible graded activation mechanism of NTSR1 is revealed by an integrated analysis utilizing the reweighted potential of mean force (PMF), deep learning (DL) and transfer entropy (TE). Specifically, the orthosteric pocket is observed to undergo expansion and contraction, with the G-protein-binding site experiencing interconversions among the inactive, intermediate and active-like states. Detailed structural comparisons capture subtle conformational differences arising from ligand binding in allosteric signaling, which can well explain the graded activation. Critical microswitches that contribute to graded activation are efficiently identified with the DL model. TE calculations enable the visualization of allosteric communication networks within the receptor, elucidating the driver-responder relationships associated with signal transduction. Fortunately, the dissociation of the full agonist from the orthosteric pocket is observed. The current findings systematically reveal the mechanism of NTSR1 graded activation, and also provide implications for structure-based drug design.

A Structure-Activity Investigation of the Fungal Metabolite (-)-TAN-2483B: Inhibition of Bruton's Tyrosine Kinase

The natural product (-)-TAN-2483B is a fungal secondary metabolite which displays promising anti-cancer and immunomodulatory activity. Our previous syntheses of (-)-TAN-2483B and sidechain analogues uncovered inhibitory activity against Bruton's tyrosine kinase (Btk), an established drug target for various leukaemia and immunological diseases. A structure-based computational study using ensemble docking and molecular dynamics was performed to determine plausible binding modes for (-)-TAN-2483B and analogues in the Btk binding site. These hypotheses guided the design of new analogues which were synthesised and their inhibitory activities determined, providing insights into the structural determinants of the furopyranone scaffold that confer both activity and selectivity for Btk. These findings offer new perspectives for generating optimised (-)-TAN-2483B-based kinase inhibitors for the treatment of leukaemia and immunological diseases.

Green approach to the synthesis of α-aminophosphonate-tetrahydroisoquinoline hybrids and their anti-cholinesterase activity

A series of 19 novel α-aminophosphonate-tetrahydroisoquinoline hybrids were synthesized through a cross dehydrogenative coupling reaction between N-aryl-tetrahydroisoquinolines and dialkylphosphites, using tert-butyl hydroperoxide as oxidazing agent. This simple procedure provided products with high atom economy and moderate to high yields. In vitro cholinesterase inhibitory activity of these compounds was evaluated. All the synthesized compounds showed good to excellent selective inhibition against butyrylcholinesterase. Compound 3bc was found to be the most active derivative with an IC50 of 9 nM. Molecular modelling studies suggested that the inhibitor is located in the peripheral anionic site (PAS) of the enzyme and interacts with some residue of the catalytic anionic site. Kinetic studies revealed that 3bc acts as a non-competitive inhibitor. Predicted ADME showed good pharmacokinetics and drug-likeness properties for most hybrids. Each newly synthesized compound was characterized by IR, 1H NMR, 13C NMR, 31P NMR spectral studies and also HRMS. The results of this study suggest that α-aminophosphonate-tetrahydroisoquinoline hybrids can be promising lead compounds in the discovery of new and improved drugs for the treatment of Alzheimer's disease and related neurodegenerative disorders.

Synthesis, DNA Binding, and Molecular Docking Studies of Dimethylaminobenzaldehyde-Based Bioactive Schiff Bases

A new series of p-dimethylaminobenzaldehyde derivatives were tested for therapeutic potential by exploring their properties through characterization. The derivatives were synthesized by 1 : 1 condensation reaction of p-dimethylaminobenzaldehyde and substituted amines. The synthesized compounds 1–8 were characterized by different characterization techniques including IR, mass, 1H NMR, and 13C NMR spectroscopy, elemental analysis, and mass spectrometry. Furthermore, binding of these Schiff bases to Ct-DNA was examined by absorption spectroscopy, fluorescence quenching, circular dichroic, viscosity measurement, molecular docking, and molecular dynamics simulation methods. Schiff bases were tested for antimicrobial activity against bacterial species Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa, and Staphylococcus aureus by the disc diffusion method. The pharmacological treatment of Schiff bases showed that 1–8 have promising potential against tested bacterial strains. The molecular docking study of the target compounds was also carried out against B-DNA dodecamer d(CGCGAATTCGCG)2, and it has been found that 1–8 can bind to Ct-DNA via an intercalative mode. DPPH free radical and hydrogen peroxide scavenging assays were employed to assess the antioxidant potential of synthesized Schiff bases.

Dual binding mode of "bitter sugars" to their human bitter taste receptor target

The 25 human bitter taste receptors (hTAS2Rs) are responsible for detecting bitter molecules present in food, and they also play several physiological and pathological roles in extraoral compartments. Therefore, understanding their ligand specificity is important both for food research and for pharmacological applications. Here we provide a molecular insight into the exquisite molecular recognition of bitter β-glycopyranosides by one of the members of this receptor subclass, hTAS2R16. Most of its agonists have in common the presence of a β-glycopyranose unit along with an extremely structurally diverse aglycon moiety. This poses the question of how hTAS2R16 can recognize such a large number of "bitter sugars". By means of hybrid molecular mechanics/coarse grained molecular dynamics simulations, here we show that the three hTAS2R16 agonists salicin, arbutin and phenyl-β-D-glucopyranoside interact with the receptor through a previously unrecognized dual binding mode. Such mechanism may offer a seamless way to fit different aglycons inside the binding cavity, while maintaining the sugar bound, similar to the strategy used by several carbohydrate-binding lectins. Our prediction is validated a posteriori by comparison with mutagenesis data and also rationalizes a wealth of structure-activity relationship data. Therefore, our findings not only provide a deeper molecular characterization of the binding determinants for the three ligands studied here, but also give insights applicable to other hTAS2R16 agonists. Together with our results for other hTAS2Rs, this study paves the way to improve our overall understanding of the structural determinants of ligand specificity in bitter taste receptors.

GPCR drug discovery: integrating solution NMR data with crystal and cryo-EM structures

The 826 G protein-coupled receptors (GPCRs) in the human proteome regulate key physiological processes and thus have long been attractive drug targets. With the crystal structures of more than 50 different human GPCRs determined over the past decade, an initial platform for structure-based rational design has been established for drugs that target GPCRs, which is currently being augmented with cryo-electron microscopy (cryo-EM) structures of higher-order GPCR complexes. Nuclear magnetic resonance (NMR) spectroscopy in solution is one of the key approaches for expanding this platform with dynamic features, which can be accessed at physiological temperature and with minimal modification of the wild-type GPCR covalent structures. Here, we review strategies for the use of advanced biochemistry and NMR techniques with GPCRs, survey projects in which crystal or cryo-EM structures have been complemented with NMR investigations and discuss the impact of this integrative approach on GPCR biology and drug discovery.

Novel HIV-1 Non-nucleoside Reverse Transcriptase Inhibitors: A Combinatorial Approach

Highly active antiretroviral therapy (HAART) is one of the most effective means for fighting against HIV-infection. HAART primarily targets HIV-1 reverse transcriptase (RT), and 14 of 28 compounds approved by the FDA as anti-HIV drugs act on this enzyme. HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) hold a special place among HIV RT inhibitors owing to their high specificity and unique mode of action. Nonetheless, these drugs show a tendency to decrease their efficacy due to high HIV-1 variability and formation of resistant virus strains tolerant to clinically applied HIV NNRTIs. A combinatorial approach based on varying substituents within various fragments of the parent molecule that results in development of highly potent compounds is one of the approaches aimed at designing novel HIV NNRTIs. Generation of HIV NNRTIs based on pyrimidine derivatives explicitly exemplifies this approach, which is discussed in this review.

Exploring Multiple Binding Modes Using Confined Replica Exchange Molecular Dynamics

Molecular docking is extensively applied to determine the position of a ligand on its receptor despite the rather poor correspondence between docking scores and experimental binding affinities found in several studies, especially for systems structurally unrelated with those used in the scoring functions' training sets. Here, we present a method for the prediction of binding modes and binding free energies, which uses replica exchange molecular dynamics in combination with a receptor-shaped piecewise potential, confining the ligand in the proximity of the receptor surface and limiting the accessible conformational space of interest. We assess our methodology with a set of protein receptor-ligand test cases. In every case studied, the method is able to locate the ligand on the experimentally known receptor binding site, and it gives as output the binding free energy. The added value of our approach with respect to other available methods is that it quickly performs a conformational space search, providing a set of bound (or unbound) configurations, which can be used to determine phenomenological structural and energetic properties of an experimental binding state as a result of contributions provided by diversified multiple binding poses.

Focus on Chirality of HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors

Chiral HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are of great interest since one enantiomer is often more potent than the corresponding counterpart against the HIV-1 wild type (WT) and the HIV-1 drug resistant mutant strains. This review exemplifies the various studies made to investigate the effect of chirality on the antiretroviral activity of top HIV-1 NNRTI compounds, such as nevirapine (NVP), efavirenz (EFV), alkynyl- and alkenylquinazolinone DuPont compounds (DPC), diarylpyrimidine (DAPY), dihydroalkyloxybenzyloxopyrimidine (DABO), phenethylthiazolylthiourea (PETT), indolylarylsulfone (IAS), arylphosphoindole (API) and trifluoromethylated indole (TFMI) The chiral separation, the enantiosynthesis, along with the biological properties of these HIV-1 NNRTIs, are discussed.

Complex structures between the N-type calcium channel (CaV2.2) and ω-conotoxin GVIA predicted via molecular dynamics

The N-type voltage-gated Ca(2+) channel CaV2.2 is one of the important targets for pain management. ω-Conotoxins isolated from venoms of cone snails, which specifically inhibit CaV2.2, are promising scaffolds for novel analgesics. The inhibitory action of ω-conotoxins on CaV2.2 has been examined experimentally, but the modes of binding of the toxins to this and other related subfamilies of Ca(2+) channels are not understood in detail. Here molecular dynamics simulations are used to construct models of ω-conotoxin GVIA in complex with a homology model of the pore domain of CaV2.2. Three different binding modes in which the side chain of Lys2, Arg17, or Lys24 from the toxin protrudes into the selectivity filter of CaV2.2 are considered. In all the modes, the toxin forms a salt bridge with an aspartate residue of subunit II just above the EEEE ring of the selectivity filter. Using the umbrella sampling technique and potential of mean force calculations, the half-maximal inhibitory concentration (IC50) values are calculated to be 1.5 and 0.7 nM for the modes in which Lys2 and Arg17 occlude the ion conduction pathway, respectively. Both IC50 values compare favorably with the values of 0.04-1.0 nM determined experimentally. The similar IC50 values calculated for the different binding modes demonstrate that GVIA can inhibit CaV2.2 with alternative binding modes. Such a multiple-binding mode mechanism may be common for ω-conotoxins.

"Reverse" carbocyclic fleximers: synthesis of a new class of adenosine deaminase inhibitors

A series of flexible carbocyclic pyrimidine nucleosides has been designed and synthesized. In contrast to previously reported "fleximers" from our laboratory, these analogues have the connectivity of the heterocyclic base system "reversed", where the pyrimidine ring is attached to the sugar moiety, rather than the five membered imidazole ring. As was previously seen with the ribose fleximers, their inherent flexibility should allow them to adjust to enzyme binding site mutations, as well as increase the affinity for atypical enzymes. Preliminary biological screening has revealed surprising inhibition of adenosine deaminase, despite their lack of resemblance to adenosine.

Identifying and characterizing promiscuous targets: implications for virtual screening

INTRODUCTION

Ligand-based shape matching approaches have become established as important and popular virtual screening (VS) techniques. However, despite their relative success, the question of how to best choose the initial query compounds and their conformations remains largely unsolved. This issue gains importance when dealing with promiscuous targets, that is, proteins that bind multiple ligand scaffold families in one or more binding site. Conventional shape matching VS approaches assume that there is only one binding mode for a given protein target. This may be true for some targets, but it is certainly not true in all cases. Several recent studies have shown that some protein targets bind to different ligands in different ways.

AREAS COVERED

The authors discuss the concept of promiscuity in the context of virtual drug screening, and present and analyze several examples of promiscuous targets. The article also reports on the impact of the query conformation on the performance of shape-based VS and the potential to improve VS performance by using consensus shape clustering techniques.

EXPERT OPINION

The notion of polypharmacology is becoming highly relevant in drug discovery. Understanding and exploiting promiscuity present challenges and opportunities for drug discovery endeavors. The examples of promiscuity presented here suggest that promiscuous targets and ligands are much more common than previously assumed, and this should be taken into account in practical VS protocols. Although some progress has been made, there is a need to develop more sophisticated computational techniques and protocols that can identify and characterize promiscuous targets on a genomic scale.

Comparative thermodynamic studies on substrate and product binding of O-acetylserine sulfhydrylase reveals two different ligand recognition modes

BACKGROUND

The importance of understanding the detailed mechanism of cysteine biosynthesis in bacteria is underscored by the fact that cysteine is the only sulfur donor for all cellular components containing reduced sulfur. O-acetylserine sulfhydrylase (OASS) catalyzes this crucial last step in the cysteine biosynthesis and has been recognized as an important gene for the survival and virulence of pathogenic bacteria. Structural and kinetic studies have contributed to the understanding of mechanistic aspects of OASS, but details of ligand recognition features of OASS are not available. In the absence of any detailed study on the energetics of ligand binding, we have studied the thermodynamics of OASS from Salmonella typhimurium (StOASS), Haemophilus influenzae (HiOASS), and Mycobacterium tuberculosis (MtOASS) binding to their substrate O-acetylserine (OAS), substrate analogue (methionine), and product (cysteine).

RESULTS

Ligand binding properties of three OASS enzymes are studied under defined solution conditions. Both substrate and product binding is an exothermic reaction, but their thermodynamic signatures are very different. Cysteine binding to OASS shows that both enthalpy and entropy contribute significantly to the binding free energy at all temperatures (10-30°C) examined. The analyses of interaction between OASS with OAS (substrate) or methionine (substrate analogue) revealed a completely different mode of binding. Binding of both OAS and methionine to OASS is dominated by a favorable entropy change, with minor contribution from enthalpy change (ΔH(St-Met) = -1.5 ± 0.1 kJ/mol; TΔS(St-Met) = 8.2 kJ/mol) at 20°C. Our salt dependent ligand binding studies indicate that methionine binding affinity is more sensitive to [NaCl] as compared to cysteine affinity.

CONCLUSIONS

We show that OASS from three different pathogenic bacteria bind substrate and product through two different mechanisms. Results indicate that predominantly entropy driven methionine binding is not mediated through classical hydrophobic binding, instead, may involve desolvation of the polar active site. We speculate that OASS in general, may exhibit two different binding mechanisms for recognizing substrates and products.

Rigorous treatment of multispecies multimode ligand-receptor interactions in 3D-QSAR: CoMFA analysis of thyroxine analogs binding to transthyretin

For a rigorous analysis of the receptor-ligand binding, speciation of the ligands caused by ionization, tautomerism, covalent hydration, and dynamic stereoisomerism needs to be considered. Each species may bind in several orientations or conformations (modes), especially for flexible ligands and receptors. A thermodynamic description of the multispecies (MS), multimode (MM) binding events shows that the overall association constant is equal to the weighted sum of the sums of microscopic association constants of individual modes for each species, with the weights given by the unbound fractions of individual species. This expression is a prerequisite for a precise quantitative characterization of the ligand-receptor interactions in both structure-based and ligand-based structure-activity analyses. We have implemented the MS-MM correlation expression into the comparative molecular field analysis (CoMFA), which deduces a map of the binding site from structures and binding affinities of a ligand set, in the absence of experimental structural information on the receptor. The MS-MM CoMFA approach was applied to published data for binding to transthyretin of 28 thyroxine analogs, each forming up to four ionization species under physiological conditions. The published X-ray structures of several analogs, exhibiting multiple binding modes, served as templates for the MS-MM superposition of thyroxine analogs. Additional modes were generated for compounds with flexible alkyl substituents, to identify bound conformations. The results demonstrate that the MS-MM modification improved predictive abilities of the CoMFA models, even for the standard procedure with MS-MM selected species and modes. The predicted prevalences of individual modes and the generated receptor site model are in reasonable agreement with the available X-ray data. The calibrated model can help in the design of inhibitors of transthyretin amyloid fibril formation.

Thiazolidinedione-based PI3Kα inhibitors: an analysis of biochemical and virtual screening methods

A series of synthesized and commercially available compounds were assessed against PI3Kα for in vitro inhibitory activity and the results compared to binding calculated in silico. Using published crystal structures of PI3Kγ and PI3Kδ co-crystallized with inhibitors as a template, docking was able to identify the majority of potent inhibitors from a decoy set of 1000 compounds. On the other hand, PI3Kα in the apo-form, modeled by induced fit docking, or built as a homology model gave only poor results. A PI3Kα homology model derived from a ligand-bound PI3Kδ crystal structure was developed that has a good ability to identify active compounds. The docking results identified binding poses for active compounds that differ from those identified to date and can contribute to our understanding of structure-activity relationships for PI3K inhibitors.

ConfGen: a conformational search method for efficient generation of bioactive conformers

We describe the methodology, parametrization, and application of a conformational search method, called ConfGen, designed to efficiently generate bioactive conformers. We define efficiency as the ability to generate a bioactive conformation within a small total number of conformations using a reasonable amount of computer time. The method combines physics-based force field calculations with empirically derived heuristics designed to achieve efficient searching and prioritization of the ligand's conformational space. While many parameter settings are supported, four modes spanning a range of speed and quality trades-offs are defined and characterized. The validation set used to test the method is composed of ligands from 667 crystal structures covering a broad array of target and ligand classes. With the fastest mode, ConfGen uses an average of 0.5 s per ligand and generates only 14.3 conformers per ligand, at least one of which lies within 2.0 A root-mean-squared deviation of the crystal structure for 96% of the ligands. The most computationally intensive mode raises this recovery rate to 99%, while taking 8 s per ligand. Combining multiple search modes to "fill-in" holes in the conformation space or energy minimizing using an all-atom force field each lead to improvements in the recovery rates at higher resolutions. Overall, ConfGen is at least as good as competing programs at high resolution and demonstrates higher efficiency at resolutions sufficient for many downstream applications, such as pharmacophore modeling.

Looking for an active conformation of the future HIV type-1 non-nucleoside reverse transcriptase inhibitors

HIV type-1 (HIV-1) non-nucleoside reverse transcriptase inhibitors (NNRTIs) are key drugs of highly active antiretroviral therapy (HAART) in the clinical management of AIDS/HIV infection. NNRTI-based HAART regimes effectively suppress viral reproduction, are not cytotoxic and show favourable pharmacokinetic properties. First-generation NNRTIs suffer the rapid selection of viral variants, hampering the binding of inhibitors into the reverse transcriptase (RT) non-nucleoside binding site (NNBS). Efforts to improve these first inhibitors led to the discovery of second-generation NNRTIs that proved to be effective against the drug-resistant mutant HIV-1 strains. The success of such agents launched a new season of NNRTI design and synthesis. This paper reviews the characteristics of second-generation NNRTIs, including etravirine, rilpivirine, RDEA-806, UK-453061, BIRL 355 BS, IDX 899, MK-4965 and HBY 097. In particular, the binding modes of these inhibitors into the NNBS of the HIV-1 RT and the most clinically relevant mutant RTs are analysed and discussed.

Molecular dynamics study of non-nucleoside reverse transcriptase inhibitor 4-[[4-[[4-[(E)-2-cyanoethenyl]-2,6-dimethylphenyl]amino]-2-pyrimidinyl]amino]benzonitrile (TMC278/rilpivirine) aggregates: correlation between amphiphilic properties of the drug and oral bioavailability

The non-nucleoside reverse transcriptase inhibitor (NNRTI) TMC278/rilpivirine is an anti-AIDS therapeutic agent with high oral bioavailability despite its high hydrophobicity. Previous studies established a correlation between ability of the drug molecule to form stable, homogeneous populations of spherical nanoparticles (approximately 100-120 nm in diameter) at low pH in surfactant-independent fashion and good oral bioavailability. Here, we hypothesize that the drug is able to assume surfactant-like properties under physiologically relevant conditions, thus facilitating formation of nanostructures in the absence of other surfactants. The results of all-atom molecular dynamics simulations indeed show that protonated drug molecules behave as surfactants at the water/aggregate interface while neutral drug molecules assist aggregate packing via conformational variability. Our simulation results suggest that amphiphilic behavior at low pH and intrinsic flexibility influence drug aggregation and are believed to play critical roles in the favorable oral bioavailability of hydrophobic drugs.

Fragment-based discovery of BACE1 inhibitors using functional assays

Novel nonpeptidic inhibitors of beta-secretase (BACE1) have been discovered by employing a fragment-based biochemical screening approach. A diverse library of 20000 low-molecular weight compounds were screened and yielded 26 novel hits that were confirmed by biochemical and surface plasmon resonance secondary assays. We describe here fragment inhibitors cocrystallized with BACE1 in a flap open and flap closed conformation as determined by X-ray crystallography.

Conformational landscape of the human immunodeficiency virus type 1 reverse transcriptase non-nucleoside inhibitor binding pocket: lessons for inhibitor design from a cluster analysis of many crystal structures

Clustering of 99 available X-ray crystal structures of HIV-1 reverse transcriptase (RT) at the flexible non-nucleoside inhibitor binding pocket (NNIBP) provides information about features of the conformational landscape for binding non-nucleoside inhibitors (NNRTIs), including effects of mutation and crystal forms. The ensemble of NNIBP conformations is separated into eight discrete clusters based primarily on the position of the functionally important primer grip, the displacement of which is believed to be one of the mechanisms of inhibition of RT. Two of these clusters are populated by structures in which the primer grip exhibits novel conformations that differ from the predominant cluster by over 4 A and are induced by the unique inhibitors capravirine and rilpivirine/TMC278. This work identifies a new conformation of the NNIBP that may be used to design NNRTIs. It can also be used to guide more complete exploration of the NNIBP free energy landscape using advanced sampling techniques.

Efficacy and resistance of recently developed non-nucleoside reverse transcriptase inhibitors for HIV-1

Since the introduction of the HAART, non-nucleoside reverse transcriptase inhibitors (NNRTIs) have played an essential role in treating HIV: their strong antiviral potency, good metabolic profile and low pill burden make them an ideal option in the design of an optimized triple drug regimen. Nonetheless, the currently approved NNRTIs (efavirenz and nevirapine) are weighed by peculiar toxicities, while a low genetic barrier and the development of cross-resistance significantly limits their use in cases of suboptimal adherence. Many drugs are in development and they are all designed with the aim to overcome resistance problems. In this review we present data on virological efficacy and resistance profiles of some of the most promising new molecules: some (such as rilpivirine) are close to being marketed, others are in Phase II trials (IDX899 and RDEA806), others again have just completed preclinical studies and are having their first clinical evaluations (RO-5028, UK-453061 and BILR-355 BS); etravirine is already approved by the US FDA, but it is still not licensed in Europe. Other new molecules (Merck MK-4965, GlaxoSmithKline GW678284 and a pyridazinone derivative by Roche), which are currently in early-development phases, are also briefly described.

High quality binding modes in docking ligands to proteins

Multiple near-optimal conformations of protein-ligand complexes provide a better chance for accurate representation of biomolecular interactions, compared with a single structure. We present ISE-dock--a docking program which is based on the iterative stochastic elimination (ISE) algorithm. ISE eliminates values that consistently lead to the worst results, thus optimizing the search for docking poses. It constructs large sets of such poses with no additional computational cost compared with single poses. ISE-dock is validated using 81 protein-ligand complexes from the PDB and its performance was compared with those of Glide, GOLD, and AutoDock. ISE-dock has a better chance than the other three to find more than 60% top single poses under RMSD = 2.0 A and more than 80% under RMSD = 3.0 A from experimental. ISE alone produced at least one 3.0 A or better solutions among the top 20 poses in the entire test set. In 98% of the examined molecules, ISE produced solutions that are closer than 2.0 A from experimental. Paired t-tests (PTT) were used throughout to assess the significance of comparisons between the performances of the different programs. ISE-dock provides more than 100-fold docking solutions in a similar time frame as LGA in AutoDock. We demonstrate the usefulness of the large near optimal populations of ligand poses by showing a correlation between the docking results and experiments that support multiple binding modes in p38 MAP kinase (Pargellis et al., Nat Struct Biol 2002;9:268-272] and in Human Transthyretin (Hamilton, Benson, Cell Mol Life Sci 2001;58:1491-1521).

High-resolution structures of HIV-1 reverse transcriptase/TMC278 complexes: strategic flexibility explains potency against resistance mutations

TMC278 is a diarylpyrimidine (DAPY) nonnucleoside reverse transcriptase inhibitor (NNRTI) that is highly effective in treating wild-type and drug-resistant HIV-1 infections in clinical trials at relatively low doses ( approximately 25-75 mg/day). We have determined the structure of wild-type HIV-1 RT complexed with TMC278 at 1.8 A resolution, using an RT crystal form engineered by systematic RT mutagenesis. This high-resolution structure reveals that the cyanovinyl group of TMC278 is positioned in a hydrophobic tunnel connecting the NNRTI-binding pocket to the nucleic acid-binding cleft. The crystal structures of TMC278 in complexes with the double mutant K103N/Y181C (2.1 A) and L100I/K103N HIV-1 RTs (2.9 A) demonstrated that TMC278 adapts to bind mutant RTs. In the K103N/Y181C RT/TMC278 structure, loss of the aromatic ring interaction caused by the Y181C mutation is counterbalanced by interactions between the cyanovinyl group of TMC278 and the aromatic side chain of Y183, which is facilitated by an approximately 1.5 A shift of the conserved Y(183)MDD motif. In the L100I/K103N RT/TMC278 structure, the binding mode of TMC278 is significantly altered so that the drug conforms to changes in the binding pocket primarily caused by the L100I mutation. The flexible binding pocket acts as a molecular "shrink wrap" that makes a shape complementary to the optimized TMC278 in wild-type and drug-resistant forms of HIV-1 RT. The crystal structures provide a better understanding of how the flexibility of an inhibitor can compensate for drug-resistance mutations.

Short-term antiviral activity of TMC278--a novel NNRTI--in treatment-naive HIV-1-infected subjects

OBJECTIVE

To evaluate antiviral activity, pharmacokinetics, tolerability and safety of TMC278, a non-nucleoside reverse transcriptase inhibitor (NNRTI), when given as a 25, 50, 100 or 150 mg once-daily dose for 7 days to antiretroviral-naive HIV-infected subjects.

DESIGN

Randomized, double-blind, placebo-controlled, phase IIa clinical trial.

METHODS

Participants were 47 antiretroviral naive HIV-infected subjects. Primary outcome was the change in plasma HIV-1 RNA viral load from baseline to day 8. Secondary outcomes were evaluation of pharmacokinetics of TMC278, immunologic changes, safety and tolerability, and evolution of viral genotypic and phenotypic patterns.

RESULTS

Patients treated with TMC278 achieved a median decrease in plasma viral load from baseline of 1.199 log10 copies/ml compared with a 0.002 log10 copies/ml gain in the placebo group (P < 0.01). A significantly higher proportion of subjects in the TMC278 groups obtained a viral load decrease of > 1.0 log10 compared with the placebo group (25/36 versus 0/11) (P < 0.01). No significant dose differences were noted in either antiviral effect or safety. No genotypic changes associated with antiretroviral resistance were detected between baseline and the end of the trial. Plasma concentrations of TMC278 were above the target concentration (13.5 ng/ml) at all time points for all TMC278-treated subjects. The most common reported adverse event was headache (TMC278 14%; placebo 18%).

CONCLUSIONS

TMC278 showed antiviral activity when given as monotherapy for 7 days at all doses studied and the drug was safe and well tolerated. Trials of longer treatment duration with TMC278, in combination with other antiretroviral drugs, are underway to assess the long-term durability of antiviral response, safety and tolerability.

"Molecular chameleons". Design and synthesis of C-4-substituted imidazole fleximers

The synthesis of two flexible nucleosides is presented. The "fleximers" feature the purine ring system split into its imidazole and pyrimidine components. This modification serves to introduce flexibility to the nucleoside while still retaining the elements essential for molecular recognition. As a result, these structurally innovative nucleosides can more readily adapt to capricious binding sites and, as such, should find use for investigating enzyme-coenzyme as well as nucleic acid-protein interactions.

Dual binding mode of a novel series of DHODH inhibitors

Human dihydroorotate dehydrogenase (DHODH) represents an important target for the treatment of hyperproliferative and inflammatory diseases. In the cell DHODH catalyzes the rate-limiting step of the de novo pyrimidine biosynthesis. DHODH inhibition results in beneficial immunosuppressant and antiproliferative effects in diseases such as rheumatoid arthritis. Here, we present high-resolution X-ray structures of human DHODH in complex with a novel class of low molecular weight compounds that inhibit the enzyme in the nanomolar range. Some compounds showed an interesting dual binding mode within the same cocrystal strongly depending on the nature of chemical substitution. Measured in vitro activity data correlated with the prevailing mode of binding and explained the observed structure-activity relationship. Additionally, the X-ray data confirmed the competitive nature of the inhibitors toward the putative ubiquinone binding site and will guide structure-based design and synthesis of molecules with higher activity.

Design, Synthesis, Biological Evaluation, and Molecular Modeling Studies of TIBO-Like Cyclic Sulfones as Non-Nucleoside HIV-1 Reverse Transcriptase Inhibitors

TIBO- and TBO-like sulfone derivatives 1 and 2 were designed, synthesized, and tested for their ability to block the replication cycle of HIV-1 in infected cells. The anti-HIV-1 activities of sulfones 3, which were intermediates in the syntheses of 1 and 2, were also evaluated. Surprisingly, the sulfone analogues of TIBO R82913 (compounds 1) were inactive, whereas interesting results were obtained for truncated derivatives 2. Compound 2 w was the most potent among this series in cell-based assays (EC50=0.07 microM, CC50>200 microM, SI>2857). It was twofold less potent than R82913, but more selective. An X-ray crystallographic analysis was carried out to establish the absolute configuration of 2 w and its enantiomer 2 x, which were obtained by semipreparative HPLC of 2 v, one of the most potent racemates. Compounds 1-3 were proven to target HIV-1 RT. In fact, representative derivatives inhibited recombinant HIV-1 RT in vitro at concentrations similar to those active in cell-based assays. 3D QSAR studies and docking simulations were developed on TIBO- and TBO-like sulfone derivatives to rationalize their anti-HIV-1 potencies and to predict the activity of novel untested sulfone derivatives. Predictive 3D QSAR models were obtained with a receptor-based alignment by docking of TIBO- and TBO-like derivatives into the NNBS of RT.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Crystallography and the design of anti-AIDS drugs: conformational flexibility and positional adaptability are important in the design of non-nucleoside HIV-1 reverse transcriptase inhibitors

Drug resistance is a key cause of failure for treatment of HIV infection. The efficacy of non-nucleoside reverse transcriptase inhibiting (NNRTI) drugs is impaired by rapid emergence of drug-resistance mutations. A multidisciplinary effort led to the discovery of the potent NNRTIs dapivirine and etravirine, both of which are diarylpyrimidine (DAPY) derivatives. Systematic structural and molecular modeling studies of HIV-1 reverse transcriptase (RT)/NNRTI complexes revealed different modes of inhibitor binding, and some of the DAPY inhibitors can bind to RT in different conformations. The torsional flexibility ("wiggling") of the inhibitors can generate numerous conformational variants and the compactness of the inhibitors permits significant repositioning and reorientation (translation and rotation) within the pocket ("jiggling"). Such adaptations appear to be critical for the ability of the diarylpyrimidine NNRTIs to retain their potency against a wide range of drug-resistant HIV-1 RTs. Exploitation of inhibitor conformational flexibility (such as torsional flexibility about strategically located chemical bonds) can be a powerful element of drug design, especially for the design of drugs that will be effective against rapidly mutating targets (which is a collection of related targets).

Comparison of Ligand-Based and Structure-Based 3D-QSAR Approaches: A Case Study on (Aryl-)Bridged 2-Aminobenzonitriles Inhibiting HIV-1 Reverse Transcriptase

Ligand- (GRIND) and structure-based (GLUE/GRIND) 3D-QSAR approaches were compared for 55 (aryl-)bridged 2-aminobenzonitriles inhibiting HIV-1 reverse transcriptase (HIV-1 RT). The ligand-based model was built from conformers selected by in vacuo minimization. The available X-ray structure of 3v in complex with HIV-1 RT allowed comparative structure-based calculations using the new docking software GLUE for conformer selection. Both models were validated via statistics and via virtual receptor sites (VRS) considering pharmacophoric regions and mutual distances, which were also compared with experimental evidence. The statistics show slight superiority of the structure-based approach in terms of fitting and prediction. By encoding relevant molecular interaction fields (MIF) into pharmacophoric regions, 10 such regions were derived from both models; they all fit the real receptor except HBD2. Also mutual distances highly agree between the real site and both VRS. Although distances from the structure-based approach are closer to the real receptor, present data prove the validity of the ligand-based GRIND approach.

Molecular modelling prediction of ligand binding site flexibility

We have investigated the efficacy of generating multiple sidechain conformations using a rotamer library in order to find the experimentally observed ligand binding site conformation of a protein in the presence of a bound ligand. We made use of a recently published algorithm that performs an exhaustive conformational search using a rotamer library to enumerate all possible sidechain conformations in a binding site. This approach was applied to a dataset of proteins whose structures were determined by X-ray and NMR methods. All chosen proteins had two or more structures, generally involving different bound ligands. By taking one of these structures as a reference, we were able in most cases to successfully reproduce the experimentally determined conformations of the other structures, as well as to suggest alternative low-energy conformations of the binding site. In those few cases where this procedure failed, we observed that the bound ligand had induced a high-energy conformation of the binding site. These results suggest that for most proteins that exhibit limited backbone motion, ligands tend to bind to low energy conformations of their binding sites. Our results also reveal that it is possible in most cases to use a rotamer search-based approach to predict alternative low-energy protein binding site conformations that can be used by different ligands. This opens the possibility of incorporating alternative binding site conformations to improve the efficacy of docking and structure-based drug design algorithms.

“Molecular Chameleons”. Design and Synthesis of a Second Series of Flexible Nucleosides

The second series of flexible shape-modified nucleosides is introduced. The "fleximers" feature the purine ring systems split into their individual imidazole and pyrimidine components. This structural modification serves to introduce flexibility to the nucleoside while still retaining the elements essential for recognition. As a consequence, these structurally innovative nucleosides can more readily adapt to their environment and should find use as bioprobes for investigating enzyme-coenzyme binding sites as well as nucleic acid and protein interactions. Their design and synthesis is described.

Pre-Incubation of Cell-Free HIV-1 Group M Isolates with Non-Nucleoside Reverse Transcriptase Inhibitors Blocks Subsequent Viral Replication in Co-Cultures of Dendritic Cells and T Cells

In order to study the inhibitory effect of various reverse transcriptase inhibitors (RTIs) on cell-free HIV, we adapted a recently described in vitro system, based on co-cultures of dendritic cells and resting CD4 T cells, modelling early target cells during sexual transmission. The compounds tested included the second-generation non-nucleoside RTI (NNRTI) TMC-120 (R147681, dapivirine) and TMC-125 (R165335, travertine), as well as the reference nucleoside RTI AZT (zidovudine), the nucleotide RTI PMPA (tenofovir) and the NNRTI UC-781. The virus strains included the reference strain HIV-1Ba-L and six primary isolates, representative of the HIV-1 group M pandemic. They all display the non-syncytium-inducing and CCR5 receptor-using (NSI/R5) phenotype, important in transmission. Cell-free virus was immobilized on a poly-L-lysine (PLL)-treated microwell plate and incubated with compound for 1h. Afterwards, the compound was thoroughly washed away; target cells were added and cultured for 2 weeks, followed by an extended culture with highly susceptible mitogen-activated T cells. Viral production in the cultures was measured on supernatant with HIV antigen ELISA. Negative results were confirmed by showing absence of proviral DNA in the cells.

TMC-120 and TMC-125 inhibited replication of HIV-1Ba-L with average EC50 values of 38 nM and 117 nM, respectively, whereas the EC50 of UC-781 was 517 nM. Complete suppression of virus and provirus was observed at compound concentrations of 100, 300 and 1000 nM, respectively. Inhibition of all primary isolates followed the same pattern as HIV-1Ba-L. In contrast, pre-treating the virus with the nucleotide RTI PMPA and AZT failed to inhibit infection even at a concentration of 100000 nM. These data clearly suggest that NNRTIs inactivate RT enzymatic activity of different viral clades (predominant in the epidemic) and might be proposed for further testing as a sterilizing microbicide worldwide.

Docking and 3-D QSAR Studies on Indolyl Aryl Sulfones. Binding Mode Exploration at the HIV-1 Reverse Transcriptase Non-Nucleoside Binding Site and Design of Highly Active N-(2-Hydroxyethyl)carboxamide and N-(2-Hydroxyethyl)carbohydrazide Derivatives

Three-dimensional quantitative structure-activity relationship (3-D QSAR) studies and docking simulations were developed on indolyl aryl sulfones (IASs), a class of novel HIV-1 non-nucleoside reverse transcriptase (RT) inhibitors (Silvestri, et al. J. Med. Chem. 2003, 46, 2482-2493) highly active against wild type and some clinically relevant resistant strains (Y181C, the double mutant K103N-Y181C, and the K103R-V179D-P225H strain, highly resistant to efavirenz). Predictive 3-D QSAR models using the combination of GRID and GOLPE programs were obtained using a receptor-based alignment by means of docking IASs into the non-nucleoside binding site (NNBS) of RT. The derived 3-D QSAR models showed conventional correlation (r(2)) and cross-validated (q(2)) coefficients values ranging from 0.79 to 0.93 and from 0.59 to 0.84, respectively. All described models were validated by an external test set compiled from previously reported pyrryl aryl sulfones (Artico, et al. J. Med. Chem. 1996, 39, 522-530). The most predictive 3-D QSAR model was then used to predict the activity of novel untested IASs. The synthesis of six designed derivatives (prediction set) allowed disclosure of new IASs endowed with high anti-HIV-1 activities.

Paul Janssen (1926-2003)

Paul Janssen (1926-2003) performed during fifty Years an intense pharmaceutical research activity. From 1953 to 2003, he discovered numerous new drugs in various fields of pharmacology. He widened the neuroleptic spectrum with haloperidol and risperidone, the opioid one with dextromoramide, fentanyl and its short-life derivates, constipating agents like loperamide, hypnotics, anaesthetics. In the field of anti-infectious agents, he discovered azole antifungals, parasiticides among which levamisole and mebendazole. Other therapeutic classes have been enriched by JANSSEN's works: vasodilatating agents, antihypertensive and anti-allergic drugs, etc. More recently, his research was oriented towards virology namely anti-HIV drugs. When Paul JANSSEN's life ended, his scientific production seemed to be one of the most eminent in the XXth century.

Correlations between Factors Determining the Pharmacokinetics and Antiviral Activity of HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors of the Diaryltriazine and Diarylpyrimidine Classes of Compounds

OBJECTIVE

To investigate the important factors that determine the bioavailability and the antiviral activity of the diaryltriazine (DATA) and diarylpyrimidine (DAPY) non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1 in animal species and humans using cell-based assays, physicochemical and computed parameters.

METHODS

This naturalistic study included 15 parameters ranging from molecular mechanics calculations to phase I clinical trials. The calculated parameters were solvent-accessible surface area (SASA), polar surface area and Gibbs free energy of solvation. Physicochemical parameters comprised lipophilicity (octanol/water partition coefficient [cLogP]), ionisation constant (pKa), solubility and aggregate radius. Cell-based assays included human colonic adenocarcinoma cell (Caco-2) permeability (transepithelial transport), drug metabolism and antiviral activity (negative logarithm of the molar effective concentration inhibiting viral replication by 50% [pEC50]). Exposure was tested in rats, dogs and human volunteers.

RESULTS

Of the 15 parameters, eight correlated consistently among one another. Exposure (area under the plasma concentration-time curve [AUC]) in humans correlated positively with that in rats (r = 1.00), with transepithelial transport (r = 0.83), lipophilicity (r = 0.60), ionisability (r = 0.89), hydrodynamic radius of aggregates (r = 0.66) and with antiviral activity (r = 0.61). Exposure in humans was also seen to correlate negatively with SASA (r = -0.89). No consistent correlation was found between exposure in dogs and the eight parameters. Of the 14 DATA/DAPY molecules, 11 form aggregates with radii between 34 and 100 nm.

CONCLUSIONS

We observed correlations between exposure in humans with exposure in rats, transepithelial transport (Caco-2 cells), ionisability, lipophilicity, aggregate radius and SASA in the class of DATA/DAPY NNRTI compounds. The lipophilic DATA/DAPY compounds form aggregates. It can be assumed that absorption in the intestinal tract and endocytosis in infected cells of these lipophilic compounds are governed by the common phenomenon of aggregate formation. As the lymphatic system offers a pathway for intestinal uptake of aggregates, this may offer a therapeutic advantage in the treatment of HIV-1 infection. Although it was not the objective of the study, we found that the rat was a better in vivo model than the dog for the prediction of systemic exposure in this particular set of compounds.