Use of a Molecular Switch Probe to Activate or Inhibit GIRK1 Heteromers In Silico Reveals a Novel Gating Mechanism.

G protein-gated inwardly rectifying K+ (GIRK) channels form highly active heterotetramers in the body, such as in neurons (GIRK1/GIRK2 or GIRK1/2) and heart (GIRK1/GIRK4 or GIRK1/4). Based on three-dimensional atomic resolution structures for GIRK2 homotetramers, we built heterotetrameric GIRK1/2 and GIRK1/4 models in a lipid bilayer environment. By employing a urea-based activator ML297 and its molecular switch, the inhibitor GAT1587, we captured channel gating transitions and K+ ion permeation in sub-microsecond molecular dynamics (MD) simulations. This allowed us to monitor the dynamics of the two channel gates (one transmembrane and one cytosolic) as well as their control by the required phosphatidylinositol bis 4-5-phosphate (PIP2). By comparing differences in the two trajectories, we identify three hydrophobic residues in the transmembrane domain 1 (TM1) of GIRK1, namely, F87, Y91, and W95, which form a hydrophobic wire induced by ML297 and de-induced by GAT1587 to orchestrate channel gating. This includes bending of the TM2 and alignment of a dipole of two acidic GIRK1 residues (E141 and D173) in the permeation pathway to facilitate K+ ion conduction. Moreover, the TM movements drive the movement of the Slide Helix relative to TM1 to adjust interactions of the CD-loop that controls the gating of the cytosolic gate. The simulations reveal that a key basic residue that coordinates PIP2 to stabilize the pre-open and open states of the transmembrane gate flips in the inhibited state to form a direct salt-bridge interaction with the cytosolic gate and destabilize its open state.

Design, synthesis and biological evaluation of renin inhibitors guided by simulated annealing of chemical potential simulations.

We have applied simulated annealing of chemical potential (SACP) to a diverse set of ∼150 very small molecules to provide insights into new interactions in the binding pocket of human renin, a historically difficult target for which to find low molecular weight (MW) inhibitors with good bioavailability. In one of its many uses in drug discovery, SACP provides an efficient, thermodynamically principled method of ranking chemotype replacements for scaffold hopping and manipulating physicochemical characteristics for drug development. We introduce the use of Constrained Fragment Analysis (CFA) to construct and analyze ligands composed of linking those fragments with predicted high affinity. This technique addresses the issue of effectively linking fragments together and provides a predictive mechanism to rank order prospective inhibitors for synthesis. The application of these techniques to the identification of novel inhibitors of human renin is described. Synthesis of a limited set of designed compounds provided potent, low MW analogs (IC50s<100nM) with good oral bioavailability (F>20-58%).

A human surfactant peptide-elastase inhibitor construct as a treatment for emphysema.

Two million Americans suffer from pulmonary emphysema, costing $2.5 billion/year and contributing to 100,000 deaths/year. Emphysema is thought to result from an imbalance between elastase and endogenous inhibitors of elastase, leading to tissue destruction and a loss of alveoli. Decades of research have still not resulted in an effective treatment other than stopping cigarette smoking, a highly addictive behavior. On the basis of our previous work, we hypothesize that small molecule inhibitors of human neutrophil elastase are ineffective because of rapid clearance from the lungs. To develop a long-acting elastase inhibitor with a lung pharmacodynamic profile that has minimal immunogenicity, we covalently linked an elastase inhibitor, similar to a trifluoro inhibitor that was used in clinical trials, to a 25-amino-acid fragment of human surfactant peptide B. We used this construct to prevent human neutrophil elastase-induced emphysema in a rodent model. The elastase inhibitor alone, although in a 70-fold molar excess to elastase in a mixture with <0.6% residual elastase activity, provided no protection from elastase-induced emphysema. Covalently combining an endogenous peptide from the target organ with a synthetic small molecule inhibitor is a unique way of endowing an active compound with the pharmacodynamic profile needed to create in vivo efficacy.

A fragment-based approach to the SAMPL3 Challenge.

The success of molecular fragment-based design depends critically on the ability to make predictions of binding poses and of affinity ranking for compounds assembled by linking fragments. The SAMPL3 Challenge provides a unique opportunity to evaluate the performance of a state-of-the-art fragment-based design methodology with respect to these requirements. In this article, we present results derived from linking fragments to predict affinity and pose in the SAMPL3 Challenge. The goal is to demonstrate how incorporating different aspects of modeling protein-ligand interactions impact the accuracy of the predictions, including protein dielectric models, charged versus neutral ligands, ΔΔGs solvation energies, and induced conformational stress. The core method is based on annealing of chemical potential in a Grand Canonical Monte Carlo (GC/MC) simulation. By imposing an initially very high chemical potential and then automatically running a sequence of simulations at successively decreasing chemical potentials, the GC/MC simulation efficiently discovers statistical distributions of bound fragment locations and orientations not found reliably without the annealing. This method accounts for configurational entropy, the role of bound water molecules, and results in a prediction of all the locations on the protein that have any affinity for the fragment. Disregarding any of these factors in affinity-rank prediction leads to significantly worse correlation with experimentally-determined free energies of binding. We relate three important conclusions from this challenge as applied to GC/MC: (1) modeling neutral ligands--regardless of the charged state in the active site--produced better affinity ranking than using charged ligands, although, in both cases, the poses were almost exactly overlaid; (2) simulating explicit water molecules in the GC/MC gave better affinity and pose predictions; and (3) applying a ΔΔGs solvation correction further improved the ranking of the neutral ligands. Using the GC/MC method under a variety of parameters in the blinded SAMPL3 Challenge provided important insights to the relevant parameters and boundaries in predicting binding affinities using simulated annealing of chemical potential calculations.

Diverse fragment clustering and water exclusion identify protein hot spots.

Simulated annealing of chemical potential located the highest affinity positions of eight organic probes and water on eight static structures of hen egg white lysozyme (HEWL) in various conformational states. In all HELW conformations, a diverse set of organic probes clustered in the known binding site (hot spot). Fragment clusters at other locations were excluded by tightly-bound waters so that only the hot-spot cluster remained in each case. The location of the hot spot was correctly predicted irrespective of the protein conformation and without accounting for protein flexibility during the simulations. Any one of the static structures could have been used to locate the hot spot. A site on a protein where a diversity of organic probes is calculated to cluster, but where water specifically does not bind, identifies a potential small-molecule binding site or protein-protein interaction hot spot.

Discovery of a novel class of non-ATP site DFG-out state p38 inhibitors utilizing computationally assisted virtual fragment-based drug design (vFBDD).

Discovery of a new class of DFG-out p38α kinase inhibitors with no hinge interaction is described. A computationally assisted, virtual fragment-based drug design (vFBDD) platform was utilized to identify novel non-aromatic fragments which make productive hydrogen bond interactions with Arg 70 on the αC-helix. Molecules incorporating these fragments were found to be potent inhibitors of p38 kinase. X-ray co-crystal structures confirmed the predicted binding modes. A lead compound was identified as a potent (p38α IC(50)=22 nM) and highly selective (≥ 150-fold against 150 kinase panel) DFG-out p38 kinase inhibitor.

Binding hot spots and amantadine orientation in the influenza a virus M2 proton channel.

Structures of truncated versions of the influenza A virus M2 proton channel have been determined recently by x-ray crystallography in the open conformation of the channel, and by NMR in the closed state. The structures differ in the position of the bound inhibitors. The x-ray structure shows a single amantadine molecule in the middle of the channel, whereas in the NMR structure four drug molecules bind at the channel's outer surface. To study this controversy we applied computational solvent mapping, a technique developed for the identification of the most druggable binding hot spots of proteins. The method moves molecular probes--small organic molecules containing various functional groups--around the protein surface, finds favorable positions using empirical free energy functions, clusters the conformations, and ranks the clusters on the basis of the average free energy. The results of the mapping show that in both structures the primary hot spot is an internal cavity overlapping the amantadine binding site seen in the x-ray structure. However, both structures also have weaker hot spots at the exterior locations that bind rimantadine in the NMR structure, although these sites are partially due to the favorable interactions with the interfacial region of the lipid bilayer. As confirmed by docking calculations, the open channel binds amantadine at the more favorable internal site, in good agreement with the x-ray structure. In contrast, the NMR structure is based on a peptide/micelle construct that is able to accommodate the small molecular probes used for the mapping, but has a too narrow pore for the rimantadine to access the internal hot spot, and hence the drug can bind only at the exterior sites.

Fragment-based design of small molecule PCSK9 inhibitors using simulated annealing of chemical potential simulations.

PCSK9 is a protein secreted by the liver that binds to the low-density lipoprotein receptor (LDLR), causing LDLR internalization, decreasing the clearance of circulating LDL particles. Mutations in PCSK9 that strengthen its interactions with LDLR result in familial hypercholesterolemia (FH) and early onset atherosclerosis, while nonsense mutations of PCSK9 result in cardio-protective hypocholesterolemia. These observations led to PCSK9 inhibition for cholesterol lowering becoming a high-interest therapeutic target, with antibody drugs reaching the market. An orally-available small molecule drug is highly desirable, but inhibiting the PCSK9/LDLR protein-protein interaction (PPI) has proven challenging. Alternate approaches to finding good lead candidates are needed. Motivated by the FH mutation data on PCSK9, we found that modeling the PCSK9/LDLR interface revealed extensive electron delocalization between and within the protein partners. Based on this, we hypothesized that compounds assembled from chemical fragments could achieve the affinity required to inhibit the PCSK9/LDLR PPI if they were selected to interact with PCSK9 in a way that, like LDLR, also involves significant fractional charge transfer to form partially covalent bonds. To identify such fragments, Simulated Annealing of Chemical Potential (SACP) fragment simulations were run on multiple PCSK9 structures, using optimized partial charges for the protein. We designed a small molecule, composed of several fragments, predicted to interact at two sites on the PCSK9. This compound inhibits the PPI with 1 µM affinity. Further, we designed two similar small molecules where one allows charge delocalization though a linker and the other doesn't. The first inhibitor with charge delocalization enhances LDLR surface expression by 60% at 10 nM, two orders of magnitude more potent than the EGF domain of LDLR. The other enhances LDLR expression by only 50% at 1 µM. This supports our conjecture that fragments can have surprisingly outsized efficacy in breaking PPI's by achieving fractional charge transfer leading to partially covalent bonding.

Conformational memories with variable bond angles.

Conformational Memories (CM) is a simulated annealing/Monte Carlo method that explores peptide and protein dihedral conformational space completely and efficiently, independent of the original conformation. Here we extend the CM method to include the variation of a randomly chosen bond angle, in addition to the standard variation of two or three randomly chosen dihedral angles, in each Monte Carlo trial of the CM exploratory and biased phases. We test the hypothesis that the inclusion of variable bond angles in CM leads to an improved sampling of conformational space. We compare the results with variable bond angles to CM with no bond angle variation for the following systems: (1) the pentapeptide Met-enkephalin, which is a standard test case for conformational search methods; (2) the proline ring pucker in a 17mer model peptide, (Ala)(8)Pro(Ala)(8); and (3) the conformations of the Ser 7.39 chi(1) in transmembrane helix 7 (TMH7) of the cannabinoid CB1 receptor, a 25-residue system. In each case, analysis of the CM results shows that the inclusion of variable bond angles results in sampling of regions of conformational space that are inaccessible to CM calculations with only variable dihedral angles, and/or a shift in conformational populations from those calculated when variable bond angles are not included. The incorporation of variable bond angles leads to an improved sampling of conformational space without loss of efficiency. Our examples show that this improved sampling leads to better exploration of biologically relevant conformations that have been experimentally validated.

Analysis of the Asymmetry of Activated EPO Receptor Enables Designing Small Molecule Agonists.

Amgen solved the high-resolution cocrystal structure of erythropoietin (EPO) bound to the extracellular part of the receptor (EPOR) in 1998, which reveals that the EPO-EPOR interaction surface is formed by 11 salt bridges, 17 H-bonds, and 2 hydrophobic clusters centered at a pair of crucial phenylalanines (F93). The EPOR has two domains, one that penetrates the membrane and a second extracellular domain that forms one arm of the binding site for the EPO ligand. The complete competent receptor-binding site is a homodimer of EPOR with the two arms forming a funnel-shaped cup where EPO binds. The two binding arms of the EPOR dimer meet at the membrane at a 120 degree angle, which Amgen characterizes as, "erythropoietin imposes a unique angular relationship and orientation that is responsible for optimal signaling." They come to this conclusion, because the EPOR cocrystallized with 2 equivalents of a 20 residue EPO mimetic peptide created at Robert Wood Johnson (RWJ) activates the receptor with a 3 order of magnitude reduction in potency, and the binding arms are forced to meet at the membrane with an angle of 180 degrees. The vast interaction surface between EPO and EPOR forms a singularly important three-dimensional structure responsible for hematopoietic stem cell proliferation and differentiation-this is Amgen's conclusion. This goal of this work is to present experimental and computational evidence that the Amgen structure is a postsignaling off-state and that the RWJ structure with the partially active peptide mimetics is an on-state. A detailed side-by-side comparison of the two structures will be presented along with literature evidence that calls into question the Amgen claim that their structure is a unique on-state. A computational fragment-based drug discovery method applied to the RWJ structure was used to locate and characterize a new predicted small molecule binding site and a fragment analysis was performed based on theories of asymmetry to create a proposed agonist with MW<300. When this molecule was experimentally tested, it displaced radiolabeled EPO with nanomolar potency and transformed human hematopoietic stem cells into red blood cells with subnanomolar potency. Obviously, this small molecule makes none of the EPO-EPOR interactions that Amgen stated were essential for fully turning on the receptor and provides strong evidence that stabilizing receptor asymmetry, not specific interactions, is the critical factor needed for activating signal transduction. Finally, when the agonist was altered to remove the asymmetric component, it still was able to displace radiolabeled EPO in competition binding experiments, but it no longer activated the receptor.

Hot-spot identification on a broad class of proteins and RNA suggest unifying principles of molecular recognition.

Chemically diverse fragments tend to collectively bind at localized sites on proteins, which is a cornerstone of fragment-based techniques. A central question is how general are these strategies for predicting a wide variety of molecular interactions such as small molecule-protein, protein-protein and protein-nucleic acid for both experimental and computational methods. To address this issue, we recently proposed three governing principles, (1) accurate prediction of fragment-macromolecule binding free energy, (2) accurate prediction of water-macromolecule binding free energy, and (3) locating sites on a macromolecule that have high affinity for a diversity of fragments and low affinity for water. To test the generality of these concepts we used the computational technique of Simulated Annealing of Chemical Potential to design one small fragment to break the RecA-RecA protein-protein interaction and three fragments that inhibit peptide-deformylase via water-mediated multi-body interactions. Experiments confirm the predictions that 6-hydroxydopamine potently inhibits RecA and that PDF inhibition quantitatively tracks the water-mediated binding predictions. Additionally, the principles correctly predict the essential bound waters in HIV Protease, the surprisingly extensive binding site of elastase, the pinpoint location of electron transfer in dihydrofolate reductase, the HIV TAT-TAR protein-RNA interactions, and the MDM2-MDM4 differential binding to p53. The experimental confirmations of highly non-obvious predictions combined with the precise characterization of a broad range of known phenomena lend strong support to the generality of fragment-based methods for characterizing molecular recognition.

Characterization of protein-ligand interaction sites using experimental and computational methods.

The ability to identify the sites of a protein that can bind with high affinity to small, drug-like compounds has been an important goal in drug design. Accurate prediction of druggable sites and the identification of small compounds binding in those sites have provided the input for fragment-based combinatorial approaches that allow for a more thorough exploration of the chemical space, and that have the potential to yield molecules that are more lead-like than those found using traditional high-throughput screening. Current progress in experimental and computational methods for identifying and characterizing druggable ligand binding sites on protein targets is reviewed herein, including a discussion of successful nuclear magnetic resonance, X-ray crystallography and tethering technologies. Classical geometric and energy-based computational methods are also discussed, with particular focus on two powerful technologies, that is, computational solvent mapping and grand canonical Monte Carlo simulations (as used by Locus Pharmaceuticals Inc). Both methods can be used to reliably identify druggable sites on proteins and to facilitate the design of novel, low-nanomolar-affinity ligands.

Designing a small molecule erythropoietin mimetic.

Erythropoietin (EPO) is a protein made by the kidneys in response to low red blood cell count that is secreted into the bloodstream and binds to a receptor on hematopoietic stem cells in the bone marrow inducing them to become new red blood cells. EPO made with recombinant DNA technology was brought to market in the 1980s to treat anemia caused by kidney disease and cancer chemotherapy. Because EPO infusion was able to replace blood transfusions in many cases, it rapidly became a multibillion dollar per year drug and as the first biologic created with recombinant technology it launched the biotech industry. For many years intense research was focused on creating a small molecule orally available EPO mimetic. The Robert Wood Johnson (RWJ) group seemed to definitively establish that only large peptides with a minimum of 60 residues could replace EPO, as anything less was not a full agonist. An intense study of the published work led me to hypothesize that the size of the mimetic is not the real issue, but the symmetry making and breaking of the EPO receptor induced by the ligand is the key to activating the stem cells. This analysis meant that residues in the binding site of the receptor deemed absolutely essential for ligand binding and activation from mutagenesis experiments, were probably not really that important. My fundamental hypotheses were: (a) the symmetric state of the homodimeric receptor is the most stable state and thus must be the off-state, (b) a highly localized binding site exists at a pivot point where the two halves of the receptor meet, (c) small molecules can be created that have high potency for this site that will be competitive with EPO and thus can displace the protein-protein interaction, (d) small symmetric molecules will stabilize the symmetric off-state of the receptor, and (e) a key asymmetry in the small molecule will stabilize a mirror image asymmetry in the receptor resulting in the stabilization of the on-state and proliferation of the stem cells into red blood cells. Researchers at Amgen published a co-crystal structure of EPO bound to the EPO receptor, which has a beautiful twofold symmetry-it was argued that this is the active state of the receptor. Activating the EPO receptor with EPO induces an almost instantaneous shutdown mechanism to sharply curtail any proliferative signal transduction, and thus, my hypotheses lead to the conclusion that the Amgen co-crystal is actually the state after receptor downregulation and thus an off-state. To put these hypotheses to the test, my computational method of Simulated Annealing of Chemical Potential was run using the co-crystal created at RWJ, which is the receptor trapped in a partial agonist state. The simulations predicted a previously unknown high affinity binding site at the pivot point where the two halves of the dimeric receptor meet, and detailed analysis of the fragment patterns led to the prediction of a molecule less than 300 MW that is basically twofold symmetric with a chiral center on one side and not the other. Thus, to the degree that computer simulations can be taken seriously, these results support my hypotheses on small molecule receptor activation. When this small molecule was synthesized and tested it indeed induced human hematopoietic stems cells to become red blood cells. When the predicted chiral center of this molecule was removed eliminating its one asymmetric feature, the resulting molecule was an antagonist-it could potently displace hot EPO but could no longer induce stem cell proliferation and differentiation. These results provided strong support for my theories on how to create potent small molecule EPO agonists and were used to launch the new company Locus Pharmaceuticals. These molecules, however, required significant chemical changes in order to make them stable in other in vitro assays and to be in vivo active, but these alterations had to be done in a way that maintained the symmetry-asymmetry considerations that led to the creation of an in vitro active molecule. The combination of changing functional groups to enable good pharmacokinetics, while not changing the key intrinsic symmetry properties were never seriously pursued at Locus and the program died. Investigations into how red blood cells are created have occupied many prominent researchers for the entire twentieth century. In the second half of the century EPO was discovered and by the end of the century it became a blockbuster commercial product that launched the biotech revolution.

Designing an orally available nontoxic p38 inhibitor with a fragment-based strategy.

The MAPK p38 became a focal point of inflammatory research when it was recognized that it played a key role in the production of the pro-inflammatory molecules TNF-alpha, IL-beta, and cyclooxygenase-2 (COX-2). The pharmaceutical industry devoted enormous efforts to creating p38 inhibitors, because blocking p38 had the potential of downregulating a group of pro-inflammatory mediators, and thus, one drug could have a cocktail effect. The market potential seemed to be clearly established (Bonafede et al., Clinicoecon Outcomes Res 6:381-388, 2014) with a multiplicity of TNF-alpha antibodies and a soluble receptor (Mewar and Wilson, Br J Pharmacol 162:785-791, 2011) already on the market, although the relationship between TNF-alpha production and p38 activation is a complicated two-way (Sabio and Davis, Semin Immunol 26:237-245, 2014) signal transduction process. With the discovery that activated p38 stabilizes (Mancini and Di Battista, Inflamm Res 60:1083-1092, 2011) COX-2 mRNA and upregulates expression of IL-beta (Bachstetter and Van Eldik, Aging Dis 1:199-211, 2010) probably in a similar manner, inhibiting p38 appeared to be a way of blocking TNF-alpha, COX-2, and IL-beta simultaneously. At Locus Pharmaceuticals we jumped on this opportunity, because we believed that our fragment-based drug discovery approach was ideally suited for making a potent small molecule p38 inhibitor that did not bind in the ATP site, but also had the solubility, lack of planarity, and low molecular weight required of a clinical candidate. Just to be clear, in our experience highly planar compounds often result in poor pharmacokinetics, because they tend to bind strongly to plasma proteins. At Locus we typically repeated assays by adding increasing amounts of plasma to check for potency degradation in the presence of blood. We found this to be a useful early indicator of pharmacokinetics and in vivo behavior. It became clear from our work and the work of others that binding to the ATP site resulted in nonspecific isoform toxicities, but binding in the adjacent allosteric DFG-site resulted in molecules that were too planar and too hydrophobic. Applying the computational method of Simulated Annealing of Chemical Potential (SACP) to this problem, we at Locus were able to come up with surprising fragment substitution patterns that led to potent non-ATP p38 inhibitors with the solubility and lack of planarity that resulted in potent in vivo efficacy in rodents with 33 % oral bioavailability. By using the simulations, we made only a small number of molecules and created a high quality clinical candidate. We also did extensive co-crystallography work, which demonstrated that the compounds bound in the mode predicted by the simulations. Unfortunately, all p38 programs ultimately shut down, because compelling evidence emerged that inhibiting p38 had no long-term clinical (Genovese, Arthritis Rheum 60:317-320, 2009) benefit. Devoting a large amount of limited resources to a target that ultimately turns out to be a mistake because it was not properly validated is a fatal error for a small company, and this is one of the reasons that Locus ultimately failed.

Conformational memories and the endocannabinoid binding site at the cannabinoid CB1 receptor.

Endocannabinoid stucture-activity relationships (SAR) indicate that the CB1 receptor recognizes ethanolamides whose fatty acid acyl chains have 20 or 22 carbons, with at least three homoallylic double bonds and saturation in at least the last five carbons of the acyl chain. To probe the molecular basis for these acyl chain requirements, the method of conformational memories (CM) was used to study the conformations available to an n-6 series of ethanolamide fatty acid acyl chain congeners: 22:4, n-6 (K(i) = 34.4 +/- 3.2 nM); 20:4, n-6 (K(i) = 39.2 +/- 5.7 nM); 20:3, n-6 (K(i) = 53.4 +/- 5.5 nM); and 20:2, n-6 (K(i) > 1500 nM). CM studies indicated that each analogue could form both extended and U/J-shaped families of conformers. However, for the low affinity 20:2, n-6 ethanolamide, the higher populated family was the extended conformer family, while for the other analogues in the series, the U/J-shaped family had the higher population. In addition, the 20:2, n-6 ethanolamide U-shaped family was not as tightly curved as were those of the other analogues studied. To quantitate this variation in curvature, the radius of curvature (in the C-3 to C-17 region) of each member of each U/J-shaped family was measured. The average radii of curvature (with their 95% confidence intervals) were found to be 5.8 A (5.3-6.2) for 20:2, n-6; 4.4 A (4.1-4.7) for 20:3, n-6; 4.0 A (3.7-4.2) for 20:4, n-6; and 4.0 A (3.6-4.5) for 22:4, n-6. These results suggest that higher CB1 affinity is associated with endocannabinoids that can form tightly curved structures. Endocannabinoid SAR also indicate that the CB1 receptor does not tolerate large endocannabinoid headgroups; however, it does recognize both polar and nonpolar moieties in the headgroup region. To identify a headgroup orientation that results in high CB1 affinity, a series of dimethyl anandamide analogues (R)-N-(1-methyl-2-hydroxyethyl)-2-(R)-methyl-arachidonamide (K(i) = 7.42 +/- 0.86 nM), (R)-N-(1-methyl-2-hydroxyethyl)-2-(S)-methyl-arachidonamide (K(i) = 185 +/- 12 nM), (S)-N-(1-methyl-2-hydroxyethyl)-2-(S)-methyl-arachidonamide (K(i) = 389 +/- 72 nM), and (S)-N-(1-methyl-2-hydroxyethyl)-2-(R)-methyl-arachidonamide (K(i) = 233 +/- 69 nM) were then studied using CM and computer receptor docking studies in an active state (R) model of CB1. These studies suggested that the high CB1 affinity of the R,R stereoisomer is due to the ability of the headgroup to form an intramolecular hydrogen bond between the carboxamide oxygen and the headgroup hydroxyl that orients the C2 and C1' methyl groups to have hydrophobic interactions with valine 3.32(196), while the carboxamide oxygen forms a hydrogen bond with lysine 3.28(192) at CB1. In this position in the CB1 binding pocket, the acyl chain has hydrophobic and C-H.pi interactions with residues in the transmembrane helix (TMH) 2-3-7 region. Taken together, the studies reported here suggest that anandamide and its congeners adopt tightly curved U/J-shaped conformations at CB1 and suggest that the TMH 2-3-7 region is the endocannabinoid binding region at CB1.

A critical role for a tyrosine residue in the cannabinoid receptors for ligand recognition.

Previous mutation and modeling studies have identified an aromatic cluster in the transmembrane helix (TMH) 3-4-5 region as important for ligand binding at the CB(1) and CB(2) cannabinoid receptors. Through novel mixed mode Monte Carlo/Stochastic Dynamics (MC/SD) calculations, we tested the importance of aromaticity at position 5.39(275) in CB(1). MC/SD calculations were performed on wild-type (WT) CB(1) and two mutants, Y5.39(275)F and Y5.39(275)I. Results indicated that while the CB(1) Y5.39(275)F mutant is very similar to WT, the Y5.39(275)I mutant shows pronounced topology changes in the TMH 3-4-5 region. Site-directed mutagenesis studies of tyrosine 5.39 to phenylalanine (Y-->F) or isoleucine (Y-->I) in both CB(1) and CB(2) were performed to determine the functional role of this amino acid in each receptor subtype. HEK 293 cells transfected with mutant receptor cDNAs were evaluated in radioligand binding and cyclic AMP assays. The CB(1) mutant and WT receptors were also co-expressed with G-protein-coupled inwardly rectifying channels (GIRK1 and GIRK4) in Xenopus oocytes to assess functional coupling. The Y-->F mutation resulted in cannnabinoid receptors with subtle differences in WT binding and signal transduction. In contrast, the Y-->I mutations produced receptors that could not produce signal transduction or bind to multiple cannabinoid compounds. However, immunofluorescence data indicate that the Y-->I mutation was compartmentalized and expressed at a level similar to that of the WT cannabinoid receptor. These results underscore the importance of aromaticity at position CB(1) 5.39(275) and CB(2) 5.39(191) for ligand recognition in the cannabinoid receptors.

Molecular modeling and functional confirmation of a predicted fatty acid binding site of mitochondrial aspartate aminotransferase.

Molecular interactions are necessary for proteins to perform their functions. The identification of a putative plasma membrane fatty acid transporter as mitochondrial aspartate aminotransferase (mAsp-AT) indicated that the protein must have a fatty acid binding site. Molecular modeling suggests that such a site exists in the form of a 500-Å(3) hydrophobic cleft on the surface of the molecule and identifies specific amino acid residues that are likely to be important for binding. The modeling and comparison with the cytosolic isoform indicated that two residues (Arg201 and Ala219) were likely to be important to the structure and function of the binding site. These residues were mutated to determine if they were essential to that function. Expression constructs with wild-type or mutated cDNAs were produced for bacteria and eukaryotic cells. Proteins expressed in Escherichia coli were tested for oleate binding affinity, which was decreased in the mutant proteins. 3T3 fibroblasts were transfected with expression constructs for both normal and mutated forms. Plasma membrane expression was documented by indirect immunofluorescence before [(3)H]oleic acid uptake kinetics were assayed. The V(max) for uptake was significantly increased by overexpression of the wild-type protein but changed little after transfection with mutated proteins, despite their presence on the plasma membrane. The hydrophobic cleft in mAsp-AT can serve as a fatty acid binding site. Specific residues are essential for normal fatty acid binding, without which fatty acid uptake is compromised. These results confirm the function of this protein as a fatty acid binding protein.

Grand canonical Monte Carlo simulation of ligand-protein binding.

A new application of the grand canonical thermodynamics ensemble to compute ligand-protein binding is described. The described method is sufficiently rapid that it is practical to compute ligand-protein binding free energies for a large number of poses over the entire protein surface, thus identifying multiple putative ligand binding sites. In addition, the method computes binding free energies for a large number of poses. The method is demonstrated by the simulation of two protein-ligand systems, thermolysin and T4 lysozyme, for which there is extensive thermodynamic and crystallographic data for the binding of small, rigid ligands. These low-molecular-weight ligands correspond to the molecular fragments used in computational fragment-based drug design. The simulations correctly identified the experimental binding poses and rank ordered the affinities of ligands in each of these systems.

Beta2 adrenergic receptor activation. Modulation of the proline kink in transmembrane 6 by a rotamer toggle switch.

In many rhodopsin-like G-protein-coupled receptors, agonist binding to a cluster of aromatic residues in TM6 may promote receptor activation by altering the configuration of the TM6 Pro-kink and by the subsequent movement of the cytoplasmic end of TM6 away from TM3. We hypothesized that the highly conserved Cys(6.47), in the vicinity of the conserved Pro(6.50), modulates the configuration of the aromatic cluster and the TM6 Pro-kink through specific interactions in its different rotamer configurations. In the beta(2) adrenergic receptor, mutation of Cys(6.47) to Thr, which in an alpha-helix has a different rotamer distribution from Cys and Ser, produced a constitutively active receptor, whereas the Ser mutant was similar to wild-type receptor. Use of the biased Monte Carlo technique of Conformational Memories showed that the rotamer changes among Cys/Ser/Thr(6.47), Trp(6.48), and Phe(6.52) are highly correlated, representing a rotamer "toggle switch" that may modulate the TM6 Pro-kink. Differential modulation of the accessibility of Cys(6.47) and an engineered Cys(6.52) in wild type and a constitutively active background provides experimental support for the association of this rotamer switch with receptor activation.