Efficient calculation of SAMPL4 hydration free energies using OMEGA, SZYBKI, QUACPAC, and Zap TK.

Several submissions for the SAMPL4 hydration free energy set were calculated using OpenEye tools, including many that were among the top performing submissions. All of our best submissions used AM1BCC charges and Poisson-Boltzmann solvation. Three submissions used a single conformer for calculating the hydration free energy and all performed very well with mean unsigned errors ranging from 0.94 to 1.08 kcal/mol. These calculations were very fast, only requiring 0.5-2.0 s per molecule. We observed that our two single-conformer methodologies have different types of failure cases and that these differences could be exploited for determining when the methods are likely to have substantial errors.

Mapping the binding of GluN2B-selective N-methyl-D-aspartate receptor negative allosteric modulators.

We have used recent structural advances in our understanding of the N-methyl-d-aspartate (NMDA) receptor amino terminal domain to explore the binding mode of multiple diaryl GluN2B-selective negative allosteric modulators at the interface between the GluN1 and GluN2B amino-terminal domains. We found that interaction of the A ring within the binding pocket seems largely invariant for a variety of structurally distinct ligands. In addition, a range of structurally diverse linkers between the two aryl rings can be accommodated by the binding site, providing a potential opportunity to tune interactions with the ligand binding pocket via changes in hydrogen bond donors, acceptors, as well as stereochemistry. The most diversity in atomic interactions between protein and ligand occur in the B ring, with functional groups that contain electron donors and acceptors providing additional atomic contacts within the pocket. A cluster of residues distant to the binding site also control ligand potency, the degree of inhibition, and show ligand-induced increases in motion during molecular dynamics simulations. Mutations at some of these residues seem to distinguish between structurally distinct ligands and raise the possibility that GluN2B-selective ligands can be divided into multiple classes. These results should help facilitate the development of well tolerated GluN2B subunit-selective antagonists.

The SAMPL3 blind prediction challenge: transfer energy overview.

Prediction of the free energy of solvation of a small molecule, or its transfer energy, is a necessary step along the path towards calculating the interactions between molecules that occur in an aqueous environment. A set of these transfer energies were gathered from the literature for series of chlorinated molecules with varying numbers of chlorines based on ethane, biphenyl, and dibenzo-p-dioxin. This focused set of molecules were then provided as a blinded challenge to assess the ability of current computational solvation methods to accurately model the interactions between water and increasingly chlorinated compounds. This was presented as part of the SAMPL3 challenge, which represented the fourth iterative blind prediction challenge involving transfer energies. The results of this exercise demonstrate that the field in general has difficulty predicting the transfer energies of more highly chlorinated compounds, and that methods seem to be erring in the same direction.

Blind prediction of host-guest binding affinities: a new SAMPL3 challenge.

The computational prediction of protein-ligand binding affinities is of central interest in early-stage drug-discovery, and there is a widely recognized need for improved methods. Low molecular weight receptors and their ligands--i.e., host-guest systems--represent valuable test-beds for such affinity prediction methods, because their small size makes for fast calculations and relatively facile numerical convergence. The SAMPL3 community exercise included the first ever blind prediction challenge for host-guest binding affinities, through the incorporation of 11 new host-guest complexes. Ten participating research groups addressed this challenge with a variety of approaches. Statistical assessment indicates that, although most methods performed well at predicting some general trends in binding affinity, overall accuracy was not high, as all the methods suffered from either poor correlation or high RMS errors or both. There was no clear advantage in using explicit versus implicit solvent models, any particular force field, or any particular approach to conformational sampling. In a few cases, predictions using very similar energy models but different sampling and/or free-energy methods resulted in significantly different results. The protonation states of one host and some guest molecules emerged as key uncertainties beyond the choice of computational approach. The present results have implications for methods development and future blind prediction exercises.

The SAMPL2 blind prediction challenge: introduction and overview.

The interactions between a molecule and the aqueous environment underpin any process that occurs in solution, from simple chemical reactions to protein-ligand binding to protein aggregation. Fundamental measures of the interaction between molecule and aqueous phase, such as the transfer energy between gas phase and water or the energetic difference between two tautomers of a molecule in solution, remain nontrivial to predict accurately using current computational methods. SAMPL2 represents the third annual blind prediction of transfer energies, and the first time tautomer ratios were included in the challenge. Over 60 sets of predictions were submitted, and each participant also attempted to estimate the error in their predictions, a task that proved difficult for most. The results of this blind assessment of the state of the field for transfer energy and tautomer ratio prediction both indicate where the field is performing well and point out flaws in current methods.

Structural determinants of D-cycloserine efficacy at the NR1/NR2C NMDA receptors.

We have studied relative efficacies of NR1 agonists glycine and d-cycloserine (DCS), and found efficacy to be dependent on the NR2 subunit. DCS shows partial agonism at NR1/NR2B but has higher relative efficacy than glycine at NR1/NR2C receptor. Molecular dynamics (MD) simulations of the NR1/NR2B and NR1/NR2C agonist binding domain dimer suggest only subtle differences in the interactions of DCS with NR1 binding site residues relative to glycine. The most pronounced differences were observed in the NR1/NR2C simulation between the orientation of helices F and G of the NR1 subunit. Interestingly, Helix F was previously proposed to influence receptor gating and to adopt an orientation depending on agonist efficacy. MD simulations and site-directed mutagenesis further suggest a role for residues at the agonist binding domain dimer interface in regulating DCS efficacy. To relate the structural rearrangements to receptor gating, we recorded single-channel currents from outside-out patches containing a single active NR1/NR2C receptor. DCS increased the mean open time and open probability of NR1/NR2C receptors compared with glycine. Maximum likelihood fitting of a gating model for NR1/NR2C receptor activation to the single-channel data suggests that DCS specifically accelerates the rate constant governing a fast gating step and reduces the closing rate. These changes appear to reflect a decreased activation energy for a pregating step and increased stability of the open states. We suggest that the higher efficacy of DCS at NR1/NR2C receptors involves structural rearrangements at the dimer interface and an effect on NR1/NR2C receptor pregating conformational changes.

Cyclostreptin and microtubules: is a low-affinity binding site required?

Cyclostreptin (CS) is a recently discovered natural product with cytotoxic activity caused by microtubule stabilization. It is the only known microtubule-stabilizing agent (MSA) that covalently binds to tubulin. It also exhibits the fast-binding kinetics seen for other MSAs. Through careful peptide digestion and mass spectrometry analysis, Buey et al. found that two amino acids are labeled by CS: Asn228, near the known taxane-binding site, and Thr220, in the type I microtubule pore. This led Buey et al. to propose Thr220 resides at the site previously predicted to be a way station or low-affinity site. By using molecular dynamics simulations and structural considerations of the microtubule pore and tubulin dimer, we conclude that postulation of a low-affinity site is unnecessary to explain the available experimental data. An alternative explanation views the microtubule pore as a structural entity that presents a substantial kinetic barrier to ligand passage to the known taxane-binding site-an entry point to the microtubule lumen that becomes completely blocked if cyclostreptin is bound at Thr220. Simulations of the free dimer also suggest a common mechanism of microtubule stabilization for taxane site MSAs through their conformational effect on the M-loop. Such an effect explains the low tubulin polymerization caused by cyclostreptin in vitro despite its covalent attachment.

The serine protease plasmin cleaves the amino-terminal domain of the NR2A subunit to relieve zinc inhibition of the N-methyl-D-aspartate receptors.

Zinc is hypothesized to be co-released with glutamate at synapses of the central nervous system. Zinc binds to NR1/NR2A N-methyl-d-aspartate (NMDA) receptors with high affinity and inhibits NMDAR function in a voltage-independent manner. The serine protease plasmin can cleave a number of substrates, including protease-activated receptors, and may play an important role in several disorders of the central nervous system, including ischemia and spinal cord injury. Here, we demonstrate that plasmin can cleave the native NR2A amino-terminal domain (NR2A(ATD)), removing the functional high affinity Zn(2+) binding site. Plasmin also cleaves recombinant NR2A(ATD) at lysine 317 (Lys(317)), thereby producing a approximately 40-kDa fragment, consistent with plasmin-induced NR2A cleavage fragments observed in rat brain membrane preparations. A homology model of the NR2A(ATD) predicts that Lys(317) is near the surface of the protein and is accessible to plasmin. Recombinant expression of NR2A with an amino-terminal deletion at Lys(317) is functional and Zn(2+) insensitive. Whole cell voltage-clamp recordings show that Zn(2+) inhibition of agonist-evoked NMDA receptor currents of NR1/NR2A-transfected HEK 293 cells and cultured cortical neurons is significantly reduced by plasmin treatment. Mutating the plasmin cleavage site Lys(317) on NR2A to alanine blocks the effect of plasmin on Zn(2+) inhibition. The relief of Zn(2+) inhibition by plasmin occurs in PAR1(-/-) cortical neurons and thus is independent of interaction with protease-activated receptors. These results suggest that plasmin can directly interact with NMDA receptors, and plasmin may increase NMDA receptor responses through disruption or removal of the amino-terminal domain and relief of Zn(2+) inhibition.

Structural insights into phenylethanolamines high-affinity binding site in NR2B from binding and molecular modeling studies.

BACKGROUND: Phenylethanolamines selectively bind to NR2B subunit-containing N-methyl-D-aspartate-subtype of ionotropic glutamate receptors and negatively modulate receptor activity. To investigate the structural and functional properties of the ifenprodil binding domain on the NR2B protein, we have purified a soluble recombinant rat NR2B protein fragment comprising the first ~400 amino acid amino-terminal domain (ATD2B) expressed in E. coli. Spectral measurements on refolded ATD2B protein demonstrated specific binding to ifenprodil. We have used site-directed mutagenesis, circular dichroism spectroscopy and molecular modeling to obtain structural information on the interactions between critical amino acid residues and ifenprodil of our soluble refolded ATD2B proteins. Ligand-induced changes in protein structure were inferred from changes in the circular dichroism spectrum, and the concentration dependence of these changes was used to determine binding constants for ifenprodil and its analogues. RESULTS: Ligand binding of ifenprodil, RO25,6981 and haloperidol on soluble recombinant ATD2B determined from circular dichroism spectroscopy yielded low-to-high micromolar equilibrium constants which concurred with functional IC50 measurement determined in heterologously expressed NR1/NR2B receptors in Xenopus oocytes. Amino acid residue substitutions of Asp101, Ile150 and Phe176 with alanine residue within the ATD2B protein altered the recombinant protein dissociation constants for ifenprodil, mirroring the pattern of their functional phenotypes. Molecular modeling of ATD2B as a clam-shell-like structure places these critical residues near a putative ligand binding site. CONCLUSION: We report for the first time biochemical measurements show that the functional measurements actually reflect binding to the ATD of NR2B subunit. Insights gained from this study help advance the theory that ifenprodil is a ligand for the ATD of NR2B subunit.

Enantiomeric propanolamines as selective N-methyl-D-aspartate 2B receptor antagonists.

Enantiomeric propanolamines have been identified as a new class of NR2B-selective NMDA receptor antagonists. The most effective agents are biaryl structures, synthesized in six steps with overall yields ranging from 11-64%. The compounds are potent and selective inhibitors of NR2B-containing recombinant NMDA receptors with IC 50 values between 30-100 nM. Potency is strongly controlled by substitution on both rings and the centrally located amine nitrogen. SAR analysis suggests that well-balanced polarity and chain-length factors provide the greatest inhibitory potency. Structural comparisons based on 3D shape analysis and electrostatic complementarity support this conclusion. The antagonists are neuroprotective in both in vitro and in vivo models of ischemic cell death. In addition, some compounds exhibit anticonvulsant properties. Unlike earlier generation NMDA receptor antagonists and some NR2B-selective antagonists, the present series of propanolamines does not cause increased locomotion in rodents. Thus, the NR2B-selective antagonists exhibit a range of therapeutically interesting properties.

Modulation of glycine potency in rat recombinant NMDA receptors containing chimeric NR2A/2D subunits expressed in Xenopus laevis oocytes.

Heteromeric NMDARs are composed of coagonist glycine-binding NR1 subunits and glutamate-binding NR2 subunits. The majority of functional NMDARs in the mammalian central nervous system (CNS) contain two NR1 subunits and two NR2 subunits of which there are four types (A-D). We show that the potency of a variety of endogenous and synthetic glycine-site coagonists varies between recombinant NMDARs such that the highest potency is seen at NR2D-containing and the lowest at NR2A-containing NMDARs. This heterogeneity is specified by the particular NR2 subunit within the NMDAR complex since the glycine-binding NR1 subunit is common to all NMDARs investigated. To identify the molecular determinants responsible for this heterogeneity, we generated chimeric NR2A/2D subunits where we exchanged the S1 and S2 regions that form the ligand-binding domains and coexpressed these with NR1 subunits in Xenopus laevis oocytes. Glycine concentration-response curves for NMDARs containing NR2A subunits including the NR2D S1 region gave mean glycine EC(50) values similar to NR2A(WT)-containing NMDARs. However, receptors containing NR2A subunits including the NR2D S2 region or both NR2D S1 and S2 regions gave glycine potencies similar to those seen in NR2D(WT)-containing NMDARs. In particular, two residues in the S2 region of the NR2A subunit (Lys719 and Tyr735) when mutated to the corresponding residues found in the NR2D subunit influence glycine potency. We conclude that the variation in glycine potency is caused by interactions between the NR1 and NR2 ligand-binding domains that occur following agonist binding and which may be involved in the initial conformation changes that determine channel gating.

Enhanced microtubule binding and tubulin assembly properties of conformationally constrained paclitaxel derivatives.

Microtubule binding and tubulin assembly promotion by a series of conformationally restricted paclitaxel (PTX) derivatives was investigated. In these derivatives, the C-4 acetate of the taxane is tethered to the C-3' phenyl at ortho and meta positions with different length linkers. The apparent affinity of these derivatives for GMPCPP-stabilized microtubules was assessed by a competition assay, and their influence on microtubule polymerization was evaluated by measuring the critical concentration of GDP-tubulin in the presence of the respective molecule. In general, taxane derivatives with higher apparent affinity for microtubules induced tubulin assembly more efficiently. Among the derivatives, molecules with the shortest tether display the strongest affinity for microtubules. These derivatives exhibited enhanced microtubule stabilization properties and efficiently induced GDP-tubulin assembly into microtubules at low temperature of 12 degrees C and in the absence of Mg2+ ions in 0.1 M PIPES. Based on molecular dynamics simulations, we propose that the enhanced ability to assemble microtubules by these taxane derivatives is linked to their ability to effectively shape the conformation of the M-loop of tubulin for cross-protofilament interaction.

Improving solubility of NR2B amino-terminal domain of N-methyl-d-aspartate receptor expressed in Escherichia coli.

The amino-terminal domains (ATDs) of N-methyl-d-aspartate (NMDA) receptors contain binding sites for modulators and may serve as potential drug targets in neurological diseases. Here, three fusion tags (6xHis-, GST-, and MBP-) were fused to the ATD of NMDA receptor NR2B subunit (ATD2B) and expressed in Escherichia coli. Each tag's ability to confer enhanced solubility to ATD2B was assessed. Soluble ATD2B was successfully obtained as a MBP fusion protein. Dynamic light scattering revealed the protein (1mg/ml) exists as monodispersed species at 25 degrees C. Functional studies using circular dichroism showed that the soluble MBP-ATD2B bound ifenprodil in a dose-dependent manner. The dissociation constants obtained for ifenprodil were similar in the absence (64nM) and presence (116nM) of saturating concentration of maltose. Moreover, the yield of soluble MBP-ATD2B is 18 times higher than the refolded 6xHis-ATD2B. We have reported a systematic comparison of three different affinity tagging strategies and identified a rapid and efficient method to obtain large amount of ATD2B recombinant protein for biochemical and structural studies.

Subunit-specific agonist activity at NR2A-, NR2B-, NR2C-, and NR2D-containing N-methyl-D-aspartate glutamate receptors.

The four N-methyl-d-aspartate (NMDA) receptor NR2 subunits (NR2A-D) have different developmental, anatomical, and functional profiles that allow them to serve different roles in normal and neuropathological situations. Identification of subunit-selective NMDA receptor agonists, antagonists, or modulators could prove to be both valuable pharmacological tools as well as potential new therapeutic agents. We evaluated the potency and efficacy of a wide range of glutamate-like compounds at NR1/NR2A, NR1/NR2B, NR1/NR2C, and NR1/NR2D receptors. Twenty-five of 53 compounds examined exhibited agonist activity at the glutamate binding site of NMDA receptors. Concentration-response relationships were determined for these agonists at each NR2 subunit. We find consistently higher potency at the NR2D subunit for a wide range of dissimilar structures, with (2S,4R)-4-methylglutamate (SYM2081) showing the greatest differential potency between NR2A- and NR2D-containing receptors (46-fold). Analysis of chimeric NR2A/D receptors suggests that enhanced agonist potency for NR2D is controlled by residues in both of the domains (Domain1 and Domain2) that compose the bilobed agonist binding domain. Molecular dynamics (MD) simulations comparing a crystallography-based hydrated NR1/NR2A model with a homology-based NR1/NR2D hydrated model of the agonist binding domains suggest that glutamate exhibits a different binding mode in NR2D compared with NR2A that accommodates a 4-methyl substitution in SYM2081. Mutagenesis of functionally divergent residues supports the conclusions drawn based on the modeling studies. Despite high homology and conserved atomic contact residues within the agonist binding pocket of NR2A and NR2D, glutamate adopts a different binding orientation that could be exploited for the development of subunit selective agonists and competitive antagonists.

Glutamate receptors: variation in structure-function coupling.

Fast excitatory synaptic transmission in the CNS relies almost entirely on the neurotransmitter glutamate and its family of ion channel receptors. An appreciation of the coupling between agonist binding and channel opening has advanced rapidly during the past five years, largely as a result of new structural information about the agonist-binding site. Recent studies suggest that despite many structural similarities different family members use different mechanisms to translate agonist binding into channel opening.

A binding site tyrosine shapes desensitization kinetics and agonist potency at GluR2. A mutagenic, kinetic, and crystallographic study.

Binding of an agonist to the 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)-propionic acid (AMPA) receptor family of the glutamate receptors (GluRs) results in rapid activation of an ion channel. Continuous application results in a non-desensitizing response for agonists like kainate, whereas most other agonists, such as the endogenous agonist (S)-glutamate, induce desensitization. We demonstrate that a highly conserved tyrosine, forming a wedge between the agonist and the N-terminal part of the bi-lobed ligand-binding site, plays a key role in the receptor kinetics as well as agonist potency and selectivity. The AMPA receptor GluR2, with mutations in Tyr-450, were expressed in Xenopus laevis oocytes and characterized in a two-electrode voltage clamp setup. The mutation GluR2(Y450A) renders the receptor highly kainate selective, and rapid application of kainate to outside-out patches induced strongly desensitizing currents. When Tyr-450 was substituted with the larger tryptophan, the (S)-glutamate desensitization is attenuated with a 10-fold increase in steady-state/peak currents (19% compared with 1.9% at the wild type). Furthermore, the tryptophan mutant was introduced into the GluR2-S1S2J ligand binding core construct and co-crystallized with kainate, and the 2.1-A x-ray structure revealed a slightly more closed ligand binding core as compared with the wild-type complex. Through genetic manipulations combined with structural and electrophysiological analysis, we report that mutations in position 450 invert the potency of two central agonists while concurrently strongly shaping the agonist efficacy and the desensitization kinetics of the AMPA receptor GluR2.

Mechanism of partial agonism at NMDA receptors for a conformationally restricted glutamate analog.

The NMDA ionotropic glutamate receptor is ubiquitous in mammalian central neurons. Because partial agonists bind to the same site as glutamate but induce less channel activation, these compounds provide an opportunity to probe the mechanism of activation of NMDA-type glutamate receptors. Molecular dynamics simulations and site-directed mutagenesis demonstrate that the partial agonist homoquinolinate interacts differently with binding pocket residues than glutamate. Homoquinolinate and glutamate induce distinct changes in the binding pocket, and the binding pocket exhibits significantly more motion with homoquinolinate bound than with glutamate. Patch-clamp recording demonstrates that single-channel activity induced by glutamate or by homoquinolinate has identical single-channel current amplitude and mean open-channel duration but that homoquinolinate slows activation of channel opening relative to glutamate. We hypothesize that agonist-induced conformational changes in the binding pocket control the efficacy of a subunit-specific activation step that precedes the concerted global change in the receptor-channel complex associated with ion channel opening.

Structural features of the glutamate binding site in recombinant NR1/NR2A N-methyl-D-aspartate receptors determined by site-directed mutagenesis and molecular modeling.

We have used site-directed mutagenesis of amino acids located within the S1 and S2 ligand binding domains of the NR2A N-methyl-D-aspartate (NMDA) receptor subunit to explore the nature of ligand binding. Wild-type or mutated NR1/NR2A NMDA receptors were expressed in Xenopus laevis oocytes and studied using two electrode voltage clamp. We investigated the effects of mutations in the S1 and S2 regions on the potencies of the agonists L-glutamate, L-aspartate, (R,S)-tetrazol-5yl-glycine, and NMDA. Mutation of each of the corresponding residues found in the NR2A receptor subunit, suggested to be contact residues in the GluR2 alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor subunit, caused a rightward shift in the concentration-response curve for each agonist examined. None of the mutations examined altered the efficacy of glutamate as assessed by methanethiosulfonate ethylammonium potentiation of agonist-evoked currents. In addition, none of the mutations altered the potency of glycine. Homology modeling and molecular dynamics were used to evaluate molecular details of ligand binding of both wild-type and mutant receptors, as well as to explore potential explanations for agonist selectivity between glutamate receptor subtypes. The modeling studies support our interpretation of the mutagenesis data and indicate a similar binding strategy for L-glutamate and NMDA when they occupy the binding site in NMDA receptors, as has been proposed for glutamate binding to the GluR2 AMPA receptor subunit. Furthermore, we offer an explanation as to why "charge conserving" mutations of two residues in the binding pocket result in nonfunctional receptor channels and suggest a contributing molecular determinant for why NMDA is not an agonist at AMPA receptors.