Unliganded Conformation Research Articles

On the function of TRAP substrate-binding proteins: Conformational variation of the sialic acid binding protein SiaP

Tripartite ATP-independent periplasmic (TRAP) transporters are analogous to ABC transporters in that they use a substrate-binding protein to scavenge metabolites (e.g., N-acetylneuraminate) and deliver them to the membrane components for import. TRAP substrate-binding proteins are thought to bind the substrate using a two-state (open and closed) induced-fit mechanism. We solved the structure of the TRAP N-acetylneuraminate substrate-binding protein from Aggregatibacter actinomycetemcomitans (AaSiaP) in both the open ligand-free and closed liganded conformations. Surprisingly, we also observed an intermediate conformation, where AaSiaP is mostly closed and is bound to a non-cognate ligand, acetate, which hints at how N-acetylneuraminate binding stabilises a fully closed state. AaSiaP preferentially binds N-acetylneuraminate (KD = 0.4 μM) compared to N-glycolylneuraminate (KD = 4.4 μM), which is explained by the closed-N-acetylneuraminate bound structure. Small-angle X-ray scattering data alongside molecular dynamics simulations suggest the AaSiaP adopts a more open state in solution than in crystal. However, the open unliganded conformation can also sample closed conformations. Molecular dynamics simulations also demonstrate the importance of water molecules for stabilising the closed conformation. Although our data is consistent with an induced fit model of binding, we suggest that the open unliganded conformation may sample multiple states capable of binding substrate. The mechanism by which the ligand is released for import remains to be determined.

Structural insights into the role and targeting of EGFRvIII

The epidermal growth factor receptor (EGFR) is a well-known oncogenic driver in lung and other cancers. In glioblastoma multiforme (GBM), the EGFR deletion variant III (EGFRvIII) is frequently found alongside EGFR amplification. Agents targeting the EGFR axis have shown limited clinical benefits in GBM and the role of EGFRvIII in GBM is poorly understood. To shed light on the role of EGFRvIII and its potential as a therapeutic target, we determined X-ray crystal structures of a monomeric EGFRvIII extracellular region (ECR). The EGFRvIII ECR resembles the unliganded conformation of EGFR, including the orientation of the C-terminal region of domain II. Domain II is mostly disordered, but the ECR structure is compact. We selected a nanobody with preferential binding to EGFRvIII relative to EGFR and structurally defined an epitope on domain IV that is occluded in the unliganded intact EGFR. These findings suggest new avenues for EGFRvIII targeting in GBM.

The role of C-terminal helix in the conformational transition of an arginine binding protein

The thermotoga maritima arginine binding protein (TmArgBP) is a periplasmic binding protein that has a short helix at the C-terminal end (CTH), which is swapped between the two chains. We apply a coarse-grained structure-based model (SBM) and all-atom MD simulation on this protein to understand the mechanism and the role of CTH in the conformational transition. When the results of SBM simulations of TmArgBP in the presence and absence of CTH are compared, we find that CTH is strategically located at the back of the binding pocket restraining the open-state conformation thereby disengaging access to the closed-state. We also ran all-atom MD simulations of open-state TmArgBP with and without CTH and discovered that in the absence of CTH the protein could reach the closed-state within 250 ns, while in its presence, the protein remained predominantly in its open-state conformation. In the simulation started from unliganded closed-state conformation without CTH, the protein exhibited multiple transitions between the two states, suggesting CTH as an essential structural element to stabilize the open-state conformation. In another simulation that began with an unliganded closed-state conformation with CTH, the protein was able to access the open-state. In this simulation the CTH was observed to reorient itself to interact with the protein emphasizing its role in assisting the conformational change. Based on our findings, we believe that CTH not only acts as a structural element that constraints the protein in its open-state but it may also guide the protein back to its open-state conformation upon ligand unbinding.

Mechanisms of activation and desensitization of full-length glycine receptor in lipid nanodiscs

Glycinergic synapses play a central role in motor control and pain processing in the central nervous system. Glycine receptors (GlyRs) are key players in mediating fast inhibitory neurotransmission at these synapses. While previous high-resolution structures have provided insights into the molecular architecture of GlyR, several mechanistic questions pertaining to channel function are still unanswered. Here, we present Cryo-EM structures of the full-length GlyR protein complex reconstituted into lipid nanodiscs that are captured in the unliganded (closed), glycine-bound (open and desensitized), and allosteric modulator-bound conformations. A comparison of these states reveals global conformational changes underlying GlyR channel gating and modulation. The functional state assignments were validated by molecular dynamics simulations, and the observed permeation events are in agreement with the anion selectivity and conductance of GlyR. These studies provide the structural basis for gating, ion selectivity, and single-channel conductance properties of GlyR in a lipid environment.

Mechanisms of Activation and Desensitization of Full-length Glycine Receptor in Membranes

Glycinergic synapses play a central role in motor control and pain processing in the central nervoussystem. Glycine receptors (GlyR) are key players in mediating fast inhibitory neurotransmissionat these synapses. While previous high-resolution structural studies have provided insights into themolecular architecture of GlyR, several mechanistic questions pertaining to channel function arestill unknown. Here, we present Cryo-EM structures of the full-length GlyR protein reconstitutedinto lipid nanodiscs that are captured in the unliganded (closed) and glycine-bound (open anddesensitized) conformations. A comparison of the three states reveals global conformationalchanges underlying GlyR channel gating. The functional state assignments were validated bymolecular dynamics simulations of the structures incorporated in a lipid bilayer. Observedpermeation events are in agreement with the anion selectivity of the channel and the reportedsingle-channel conductance of GlyR. These studies establish the structural basis for gating,selectivity, and single-channel conductance of GlyR in a physiological environment.

Molecular dynamics simulations reveal distinct differences in conformational dynamics and thermodynamics between the unliganded and CD4-bound states of HIV-1 gp120.

The entry of human immunodeficiency virus type I (HIV-1) into host cells is initiated by binding to the cell-surface receptor CD4, which induces a conformational transition of the envelope (Env) glycoprotein gp120 from the closed, unliganded state to the open, CD4-bound state. Despite many available structures in these two states, detailed aspects on the dynamics and thermodynamics of gp120 remain elusive. Here, we performed microsecond-scale (μs-scale) multiple-replica molecular dynamics (MD) simulations to explore the differences in the conformational dynamics, protein motions, and thermodynamics between the unliganded and CD4-bound/complexed forms of gp120. Comparative analyses of MD trajectories reveal that CD4 binding promotes the structural deviations/changes and conformational flexibility, loosens the structural packing, and complicates the molecular motions of gp120. Comparison of the constructed free energy landscapes (FELs) reveals that the CD4-complexed gp120 has more conformational substates, larger conformational entropy, and lower thermostability than the unliganded form. Therefore, the unliganded conformation represents a structurally and energetically stable "ground state" for the full-length gp120. The observed great increase in the mobility of V1/V2 and V3 along with their more versatile movement directions in the CD4-bound gp120 compared to the unliganded form suggests that their orientations with respect to each other and to the structural core determine the differences in the conformational dynamics and thermodynamics between the two gp120 forms. The results presented here provide a basis by which to better understand the functional and immunological properties of gp120 and, furthermore, to deploy appropriate strategies for the development of anti-HIV-1 drugs or vaccines.

Zn2+ Modulates Hv1 Proton Channel Gating via Conformational Coupling to an Intracellular Coulombic Network

Gating of the intrinsic ‘aqueous’ H+ conductance (GAQ) in the Hv1 proton channel is sensitive to extracellular Zn2+ and other divalent cations in the micromolar range. Previous studies showed that the extracellularly-accessible residues H140 and H193 are required for Zn2+ to positively shift the GAQ-V relation. However, the molecular mechanism by which extracellular Zn2+ occupancy is coupled to voltage-dependent GAQ gating remains known. We used site-directed mutagenesis and electrophysiology to identify Zn2+-interacting residues and used the experimental data to build new models of the Hv1 voltage sensor (VS) domain in Zn2+-liganded (Hv1F·Zn2+) and unliganded (Hv1F + Zn2+) conformations. Although imidazole nitrogen atoms (H140-Nδ1 and H193-Nδ1) are close enough to simultaneously coordinate a single Zn2+ ion, we did not observe a stable Zn2+ coordination sphere during all-atom molecular dynamics (MD) simulations using the CHARMM36 force field. To simulate electronic interactions between Zn2+ and H140-Nδ1 and H193-Nδ1, we imposed harmonic constraints (5 kcal/mol/Å) and equilibrated the Hv1F·Zn2+ model during a 50 ns MD simulation. 20 ns after releasing constraints, ZN13 spontaneously dissociates from H140-Nδ1 and H193-Nδ1 but remains associated with the carboxyl group of E119. The VS also undergoes backbone reorganization involving the S1, S3 and S4 helices after ZN13 dissociation that suggest Zn2+-dependent conformational changes. Concomitant changes in the organization of an electrostatic network involving intracellular acidic residues (E153 and D174) and S4 ‘gating charge’ Arg side chains may reveal how voltage-dependent gating is altered by Zn2+. In support of this hypothesis, we find that E153 neutralization also alters Zn2+ potency. In summary, extracellular Zn2+ occupancy is coupled to reorganization of an intracellular Coulombic network that is required for voltage-dependent gating in Hv1.

ESR Resolves the C Terminus Structure of the Ligand-free Human Glutathione S-Transferase A1-1

Nitroxide- and Cu2+-based electron spin resonance (ESR) are combined to provide insight into the conformational states of the functionally important α-helix of the human glutathione S-transferase A1. Distance measurements on various spin-labeled dimeric human glutathione S-transferase A1-1 all result in bimodal distance distributions, indicating that the C-terminus exists in two distinct conformations in solution, one of which closely matches that found in the crystal structure of the ligand-bound enzyme. These measurements permit the generation of a model of the unliganded conformation. Room temperature ESR indicates that the second conformation has high mobility, potentially enabling the enzyme’s high degree of substrate promiscuity. This model is then validated using computational modeling and further Cu2+-based ESR distance measurements. Cu2+-based ESR also provides evidence that the secondary structure of the second conformation is of helical nature. Addition of S-hexyl glutathione results in a shift in relative populations, favoring the state that is similar to the previously known structure of the ligand-bound enzyme.

Improving the Immunogenicity of Native-like HIV-1 Envelope Trimers by Hyperstabilization

SummaryThe production of native-like recombinant versions of the HIV-1 envelope glycoprotein (Env) trimer requires overcoming the natural flexibility and instability of the complex. The engineered BG505 SOSIP.664 trimer mimics the structure and antigenicity of native Env. Here, we describe how the introduction of new disulfide bonds between the glycoprotein (gp)120 and gp41 subunits of SOSIP trimers of the BG505 and other genotypes improves their stability and antigenicity, reduces their conformational flexibility, and helps maintain them in the unliganded conformation. The resulting next-generation SOSIP.v5 trimers induce strong autologous tier-2 neutralizing antibody (NAb) responses in rabbits. In addition, the BG505 SOSIP.v6 trimers induced weak heterologous NAb responses against a subset of tier-2 viruses that were not elicited by the prototype BG505 SOSIP.664. These stabilization methods can be applied to trimers from multiple genotypes as components of multivalent vaccines aimed at inducing broadly NAbs (bNAbs).

Residues in the gp41 Ectodomain Regulate HIV-1 Envelope Glycoprotein Conformational Transitions Induced by gp120-Directed Inhibitors.

Interactions between the gp120 and gp41 subunits of the human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein (Env) trimer maintain the metastable unliganded form of the viral spike. Binding of gp120 to the receptor, CD4, changes the Env conformation to promote gp120 interaction with the second receptor, CCR5 or CXCR4. CD4 binding also induces the transformation of Env into the prehairpin intermediate, in which the gp41 heptad repeat 1 (HR1) coiled coil is assembled at the trimer axis. In nature, HIV-1 Envs must balance the requirements to maintain the noncovalent association of gp120 with gp41 and to evade the host antibody response with the need to respond to CD4 binding. Here we show that the gp41 HR1 region contributes to gp120 association with the unliganded Env trimer. Changes in particular amino acid residues in the gp41 HR1 region decreased the efficiency with which Env moved from the unliganded state. Thus, these gp41 changes decreased the sensitivity of HIV-1 to cold inactivation and ligands that require Env conformational changes to bind efficiently. Conversely, these gp41 changes increased HIV-1 sensitivity to small-molecule entry inhibitors that block Env conformational changes induced by CD4. Changes in particular gp41 HR1 amino acid residues can apparently affect the relative stability of the unliganded state and CD4-induced conformations. Thus, the gp41 HR1 region contributes to the association with gp120 and regulates Env transitions from the unliganded state to downstream conformations.IMPORTANCE The development of an efficient vaccine able to prevent HIV infection is a worldwide priority. Knowledge of the envelope glycoprotein structure and the conformational changes that occur after receptor engagement will help researchers to develop an immunogen able to elicit antibodies that block HIV-1 transmission. Here we identify residues in the HIV-1 transmembrane envelope glycoprotein that stabilize the unliganded state by modulating the transitions from the unliganded state to the CD4-bound state.

An Intrinsic Ligand Mediates N- and C-Terminal Interactions with the KCNH Gating Machinery

EAG, ERG and ELK are members of the KCNH potassium channel family important in excitability. Two conserved domains in the cytosolic N and C termini uniquely characterize members of this family. The C-terminal domain has homology to the cyclic nucleotide-binding domain of CNG/HCN channels but is unaffected by cyclic nucleotides. At the tail end of this cyclic nucleotide-binding homology domain (CNBhD) is an intrinsic ligand comprising two highly conserved residues that loop back to occupy the cyclic nucleotide-binding pocket. The CNBhD interacts with the N-terminal domain, with one interface also occupied by the intrinsic ligand. We probed the intrinsic ligand and its effect on C- and N-terminal modulation of gating in hEAG1 channels. Using two-electrode voltage clamp of channels expressed in Xenopus oocytes, we found that substitution of the intrinsic ligand residues, Y699 and L701, with AA, SS and GG progressively shifted the midpoint of the conductance-voltage relationship to more depolarized potentials. The substitution with G, a residue unable to support a ligand-receptor interaction and thus considered fully unliganded, resulted in the greatest steady-state shift of channels into the closed state. Interestingly, deleting the CNBhD created the same phenotype as the GG mutant, suggesting the unliganded conformation confers a loss of function of the CNBhD's effects on gating. Both AA and GG mutants abolished in vitro binding of the PAS domain to the CNBhD in anisotropy measurements. Functionally, when combining GG/CNBhD deletion with mutations at PAS-CNBhD interfaces, we observed a non-additive phenotype that is indistinguishable from the GG or CNBhD deletion mutant. These findings suggest that both the intrinsic ligand and CNBhD are required for the proper function of PAS domain, and the ligand and binding pocket can adopt multiple conformations, each conferring a different gating behavior.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study.

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal "glutaminase" (GAT) and C-terminal "synthetase" domain. The enzyme is identified as a potential target for anti-cancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

Effect of osmolytes on the binding of EGR1 transcription factor to DNA.

Osmolytes play a key role in maintaining protein stability and mediating macromolecular interactions within the intracellular environment of the cell. Herein, we show that osmolytes such as glycerol, sucrose, and polyethylene glycol 400 (PEG400) mitigate the binding of early growth response (protein) 1 (EGR1) transcription factor to DNA in a differential manner. Thus, while physiological concentrations of glycerol only moderately reduce the binding affinity, addition of sucrose and PEG400 is concomitant with a loss in the binding affinity by an order of magnitude. This salient observation suggests that EGR1 is most likely subject to conformational equilibrium and that the osmolytes exert their effect via favorable interactions with the unliganded conformation. Consistent with this notion, our analysis reveals that while EGR1 displays rather high structural stability in complex with DNA, the unliganded conformation becomes significantly destabilized in solution. In particular, while liganded EGR1 adopts a well-defined arc-like architecture, the unliganded protein samples a comparatively large conformational space between two distinct states that periodically interconvert between an elongated rod-like shape and an arc-like conformation on a submicrosecond time scale. Consequently, the ability of osmolytes to favorably interact with the unliganded conformation so as to stabilize it could account for the negative effect of osmotic stress on EGR1-DNA interaction observed here. Taken together, our study sheds new light on the role of osmolytes in modulating a key protein-DNA interaction.

CD4-Mimetic Small Molecules Sensitize Human Immunodeficiency Virus to Vaccine-Elicited Antibodies

Approaches to prevent human immunodeficiency virus (HIV-1) transmission are urgently needed. Difficulties in eliciting antibodies that bind conserved epitopes exposed on the unliganded conformation of the HIV-1 envelope glycoprotein (Env) trimer represent barriers to vaccine development. During HIV-1 entry, binding of the gp120 Env to the initial receptor, CD4, triggers conformational changes in Env that result in the formation and exposure of the highly conserved gp120 site for interaction with the coreceptors, CCR5 and CXCR4. The DMJ compounds (+)-DMJ-I-228 and (+)-DMJ-II-121 bind gp120 within the conserved Phe 43 cavity near the CD4-binding site, block CD4 binding, and inhibit HIV-1 infection. Here we show that the DMJ compounds sensitize primary HIV-1, including transmitted/founder viruses, to neutralization by monoclonal antibodies directed against CD4-induced (CD4i) epitopes and the V3 region, two gp120 elements involved in coreceptor binding. Importantly, the DMJ compounds rendered primary HIV-1 sensitive to neutralization by antisera elicited by immunization of rabbits with HIV-1 gp120 cores engineered to assume the CD4-bound state. Thus, small molecules like the DMJ compounds may be useful as microbicides to inhibit HIV-1 infection directly and to sensitize primary HIV-1 to neutralization by readily elicited antibodies. Preventing HIV-1 transmission is a priority for global health. Eliciting antibodies that can neutralize many different strains of HIV-1 is difficult, creating problems for the development of a vaccine. We found that certain small-molecule compounds can sensitize HIV-1 to particular antibodies. These antibodies can be elicited in rabbits. These results suggest an approach to prevent HIV-1 sexual transmission in which a virus-sensitizing microbicide is combined with a vaccine.

Structures of adenosine kinase from Trypanosoma brucei brucei.

Trypanosoma brucei is a single-cellular parasite of the genus Kinetoplastida and is the causative agent of African sleeping sickness in humans. Adenosine kinase is a key enzyme in the purine-salvage pathway, phosphorylating adenosine to AMP, and also activates cytotoxic analogues such as cordycepin and Ara-A by their phosphorylation. The structures of T. brucei brucei adenosine kinase (TbAK) in its unliganded open conformation and complexed with adenosine and ADP in the closed conformation are both reported to 2.6 Å resolution. The structures give insight into the binding mode of the substrates and the conformational change induced upon substrate binding. This information can be used to guide the improvement of cytotoxic substrate analogues as potential antitrypanosomal drugs.

Molecular architecture of the uncleaved HIV-1 envelope glycoprotein trimer

The human immunodeficiency virus (HIV-1) envelope glycoprotein (Env) trimer, a membrane-fusing machine, mediates virus entry into host cells and is the sole virus-specific target for neutralizing antibodies. Binding the receptors, CD4 and CCR5/CXCR4, triggers Env conformational changes from the metastable unliganded state to the fusion-active state. We used cryo-electron microscopy to obtain a 6-A structure of the membrane-bound, heavily glycosylated HIV-1 Env trimer in its uncleaved and unliganded state. The spatial organization of secondary structure elements reveals that the unliganded conformations of both gp120 and gp41 subunits differ from those induced by receptor binding. The gp120 trimer association domains, which contribute to interprotomer contacts in the unliganded Env trimer, undergo rearrangement upon CD4 binding. In the unliganded Env, intersubunit interactions maintain the gp41 ectodomain helical bundles in a “spring-loaded” conformation distinct from the extended helical coils of the fusion-active state. Quaternary structure regulates the virus-neutralizing potency of antibodies targeting the conserved CD4-binding site on gp120. Recent studies that help validate the 3-D reconstruction of the unliganded HIV-1 Env precursor map will be presented. The Env trimer architecture provides mechanistic insights into the metastability of the unliganded state, receptor-induced conformational changes, and quaternary structure-based strategies for immune evasion.

Molecular architecture of the uncleaved HIV-1 envelope glycoprotein trimer

The human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein (Env) trimer, a membrane-fusing machine, mediates virus entry into host cells and is the sole virus-specific target for neutralizing antibodies. Binding the receptors, CD4 and CCR5/CXCR4, triggers Env conformational changes from the metastable unliganded state to the fusion-active state. We used cryo-electron microscopy to obtain a 6-Å structure of the membrane-bound, heavily glycosylated HIV-1 Env trimer in its uncleaved and unliganded state. The spatial organization of secondary structure elements reveals that the unliganded conformations of both glycoprotein (gp)120 and gp41 subunits differ from those induced by receptor binding. The gp120 trimer association domains, which contribute to interprotomer contacts in the unliganded Env trimer, undergo rearrangement upon CD4 binding. In the unliganded Env, intersubunit interactions maintain the gp41 ectodomain helical bundles in a "spring-loaded" conformation distinct from the extended helical coils of the fusion-active state. Quaternary structure regulates the virus-neutralizing potency of antibodies targeting the conserved CD4-binding site on gp120. The Env trimer architecture provides mechanistic insights into the metastability of the unliganded state, receptor-induced conformational changes, and quaternary structure-based strategies for immune evasion.

Structural, Biochemical, and Functional Characterization of the Cyclic Nucleotide Binding Homology Domain from the Mouse EAG1 Potassium Channel

KCNH channels are voltage-gated potassium channels with important physiological functions. In these channels, a C-terminal cytoplasmic region, known as the cyclic nucleotide binding homology (CNB-homology) domain displays strong sequence similarity to cyclic nucleotide binding (CNB) domains. However, the isolated domain does not bind cyclic nucleotides. Here, we report the X-ray structure of the CNB-homology domain from the mouse EAG1 channel. Through comparison with the recently determined structure of the CNB-homology domain from the zebrafish ELK (eag-like K(+)) channel and the CNB domains from the MlotiK1 and HCN (hyperpolarization-activated cyclic nucleotide-gated) potassium channels, we establish the structural features of CNB-homology domains that explain the low affinity for cyclic nucleotides. Our structure establishes that the "self-liganded" conformation, where two residues of the C-terminus of the domain are bound in an equivalent position to cyclic nucleotides in CNB domains, is a conserved feature of CNB-homology domains. Importantly, we provide biochemical evidence that suggests that there is also an unliganded conformation where the C-terminus of the domain peels away from its bound position. A functional characterization of this unliganded conformation reveals a role of the CNB-homology domain in channel gating.

Novel Inhibitor Discovery through Virtual Screening against Multiple Protein Conformations Generated via Ligand-Directed Modeling: A Maternal Embryonic Leucine Zipper Kinase Example

Kinase targets have been demonstrated to undergo major conformational reorganization upon ligand binding. Such protein conformational plasticity remains a significant challenge in structure-based virtual screening methodology and may be approximated by screening against an ensemble of diverse protein conformations. Maternal embryonic leucine zipper kinase (MELK), a member of serine-threonine kinase family, has been recently found to be involved in the tumerogenic state of glioblastoma, breast, ovarian, and colon cancers. We therefore modeled several conformers of MELK utilizing the available chemogenomic and crystallographic data of homologous kinases. We carried out docking pose prediction and virtual screening enrichment studies with these conformers. The performances of the ensembles were evaluated by their ability to reproduce known inhibitor bioactive conformations and to efficiently recover known active compounds early in the virtual screen when seeded with decoy sets. A few of the individual MELK conformers performed satisfactorily in reproducing the native protein-ligand pharmacophoric interactions up to 50% of the cases. By selecting an ensemble of a few representative conformational states, most of the known inhibitor binding poses could be rationalized. For example, a four conformer ensemble is able to recover 95% of the studied actives, especially with imperfect scoring function(s). The virtual screening enrichment varied considerably among different MELK conformers. Enrichment appears to improve by selection of a proper protein conformation. For example, several holo and unliganded active conformations are better to accommodate diverse chemotypes than ATP-bound conformer. These results prove that using an ensemble of diverse conformations could give a better performance. Applying this approach, we were able to screen a commercially available library of half a million compounds against three conformers to discover three novel inhibitors of MELK, one from each template. Among the three compounds validated via experimental enzyme inhibition assays, one is relatively potent (15; K(d) = 0.37 μM), one moderately active (12; K(d) = 3.2 μM), and one weak but very selective (9; K(d) = 18 μM). These novel hits may be utilized to assist in the development of small molecule therapeutic agents useful in diseases caused by deregulated MELK, and perhaps more importantly, the approach demonstrates the advantages of choosing an appropriate ensemble of a few conformers in pursuing compound potency, selectivity, and novel chemotypes over using single target conformation for structure-based drug design in general.

Crystallographic analysis of human hemoglobin elucidates the structural basis of the potent and dual antisickling activity of pyridyl derivatives of vanillin.

Vanillin has previously been studied clinically as an antisickling agent to treat sickle-cell disease. In vitro investigations with pyridyl derivatives of vanillin, including INN-312 and INN-298, showed as much as a 90-fold increase in antisickling activity compared with vanillin. The compounds preferentially bind to and modify sickle hemoglobin (Hb S) to increase the affinity of Hb for oxygen. INN-312 also led to a considerable increase in the solubility of deoxygenated Hb S under completely deoxygenated conditions. Crystallographic studies of normal human Hb with INN-312 and INN-298 showed that the compounds form Schiff-base adducts with the N-terminus of the α-subunits to constrain the liganded (or relaxed-state) Hb conformation relative to the unliganded (or tense-state) Hb conformation. Interestingly, while INN-298 binds and directs its meta-positioned pyridine-methoxy moiety (relative to the aldehyde moiety) further down the central water cavity of the protein, that of INN-312, which is ortho to the aldehyde, extends towards the surface of the protein. These studies suggest that these compounds may act to prevent sickling of SS cells by increasing the fraction of the soluble high-affinity Hb S and/or by stereospecific inhibition of deoxygenated Hb S polymerization.