Putative Allosteric Binding Site Research Articles

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation.

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Abstract 677: Discovery of a first-in-class allosteric inhibitor of the small GTPase Rac1

Abstract The Rac1 signaling pathway is known to be hyperactivated in numerous solid tumors. We developed a Rac1 AlphaLISA assay that is amenable to high-throughput screening (HTS) to screen large libraries of small molecules. Based on molecular modeling of a putative allosteric binding site for several of the hits, we optimized a series of compounds that lead to RP-102124. RP-102124 engages with Rac1 in both biochemical and cellular target engagement assays. RP-102124 has greater potency than previously reported Rac1 inhibitors in biochemical assays, with potent inhibition of phenotypic characteristics of Rac signaling such as cell migration and cytoskeletal rearrangements. RP-102124 is a Racmulti inhibitor with activity against oncogenic Rac1b. RP-102124 is orally bioavailable and shows dose dependent inhibition of tumor growth in xenograft models at doses that are well-tolerated. RP-102124 is a first-in-class allosteric Racmulti inhibitor that shows promising therapeutic potential. Citation Format: Robert Silverman, Alan Cooper, Ronodip Kar, Mayuri V. Raman, Ramu Ravirala, Mubarak Hasan Shah, Sandeep Mahankali, Erik T. Goka. Discovery of a first-in-class allosteric inhibitor of the small GTPase Rac1 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 677.

Determination of the Cannabinoid CB1 Receptor's Positive Allosteric Modulator Binding Site through Mutagenesis Studies.

Positive allosteric modulators (PAMs) of the cannabinoid CB1 receptor (CB1) offer potential therapeutic advantages in the treatment of neuropathic pain and addiction by avoiding the adverse effects associated with orthosteric CB1 activation. Here, molecular modeling and mutagenesis were used to identify residues central to PAM activity at CB1. Six putative allosteric binding sites were identified in silico, including novel sites previously associated with cholesterol binding, and key residues within each site were mutated to alanine. The recently determined ZCZ011 binding site was found to be essential for allosteric agonism, as GAT228, GAT229 and ZCZ011 all increased wild-type G protein dissociation in the absence of an orthosteric ligand; activity that was abolished in mutants F191A3.27 and I169A2.56. PAM activity was demonstrated for ZCZ011 in the presence of the orthosteric ligand CP55940, which was only abolished in I169A2.56. In contrast, the PAM activity of GAT229 was reduced for mutants R220A3.56, L404A8.50, F191A3.27 and I169A2.56. This indicates that allosteric modulation may represent the net effect of binding at multiple sites, and that allosteric agonism is likely to be mediated via the ZCZ011 site. This study underlines the need for detailed understanding of ligand receptor interactions in the search for pure CB1 allosteric modulators.

In silico identification of a β2-adrenoceptor allosteric site that selectively augments canonical β2AR-Gs signaling and function

Activation of β2-adrenoceptors (β2ARs) causes airway smooth muscle (ASM) relaxation and bronchodilation, and β2AR agonists (β-agonists) are front-line treatments for asthma and other obstructive lung diseases. However, the therapeutic efficacy of β-agonists is limited by agonist-induced β2AR desensitization and noncanonical β2AR signaling involving β-arrestin that is shown to promote asthma pathophysiology. Accordingly, we undertook the identification of an allosteric site on β2AR that could modulate the activity of β-agonists to overcome these limitations. We employed the site identification by ligand competitive saturation (SILCS) computational method to comprehensively map the entire 3D structure of in silico-generated β2AR intermediate conformations and identified a putative allosteric binding site. Subsequent database screening using SILCS identified drug-like molecules with the potential to bind to the site. Experimental assays in HEK293 cells (expressing recombinant wild-type human β2AR) and human ASM cells (expressing endogenous β2AR) identified positive and negative allosteric modulators (PAMs and NAMs) of β2AR as assessed by regulation of β-agonist-stimulation of cyclic AMP generation. PAMs/NAMs had no effect on β-agonist-induced recruitment of β-arrestin to β2AR- or β-agonist-induced loss of cell surface expression in HEK293 cells expressing β2AR. Mutagenesis analysis of β2AR confirmed the SILCS identified site based on mutants of amino acids R131, Y219, and F282. Finally, functional studies revealed augmentation of β-agonist-induced relaxation of contracted human ASM cells and bronchodilation of contracted airways. These findings identify a allosteric binding site on the β2AR, whose activation selectively augments β-agonist-induced Gs signaling, and increases relaxation of ASM cells, the principal therapeutic effect of β-agonists.

Eltrombopag as an Allosteric Inhibitor of the METTL3-14 Complex Affecting the m6A Methylation of RNA in Acute Myeloid Leukemia Cells.

N6A-methyladenosine (m6A) post-transcriptional modification, the most abundant internal RNA modification, is catalyzed by the METTL3-14 methyltransferase complex. Recently, attention has been drawn to the METTL3-14 complex regarding its significant roles in the pathogenesis of acute myeloid leukemia (AML), attracting the potential of novel therapeutic targets for the disease. Herein, we report the identification and characterization of eltrombopag as a selective allosteric inhibitor of the METTL3-14 complex. Eltrombopag exhibited selective inhibitory activity in the most active catalytic form of the METTL3-14 complex by direct binding, and the mechanism of inhibition was confirmed as a noncompetitive inhibition by interacting at a putative allosteric binding site in METTL3, which was predicted by cavity search and molecular docking studies. At a cellular level, eltrombopag displayed anti-proliferative effects in the relevant AML cell line, MOLM-13, in correlation with a reduction in m6A levels. Molecular mechanism studies of eltrombopag using m6A-seq analysis provided further evidence of its cellular function by determining the hypomethylation of leukemogenic genes in eltrombopag-treated MOLM-13 cells and the overlapping of the pattern with those of METTL3-knockdown MOLM-13 cells. In conclusion, eltrombopag was first disclosed as a functional METTL3-14 allosteric inhibitor in AML cells, which could be utilized for the further development of novel anti-AML therapy.

Agonistic and potentiating effects of perfluoroalkyl substances (PFAS) on the Atlantic cod (Gadus morhua) peroxisome proliferator-activated receptors (Ppars)

Toxicity mediated by per- and polyfluoroalkyl substances (PFAS), and especially perfluoroalkyl acids (PFAAs), has been linked to activation of peroxisome proliferator-activated receptors (Ppar) in many vertebrates. Here, we present the primary structures, phylogeny, and tissue-specific distributions of the Atlantic cod (Gadus morhua) gmPpara1, gmPpara2, gmPparb, and gmPparg, and demonstrate that the carboxylic acids PFHxA, PFOA, PFNA, as well as the sulfonic acid PFHxS, activate gmPpara1 in vitro, which was also supported by in silico analyses. Intriguingly, a binary mixture of PFOA and the non-activating PFOS produced a higher activation of gmPpara1 compared to PFOA alone, suggesting that PFOS has a potentiating effect on receptor activation. Supporting the experimental data, docking and molecular dynamics simulations of single and double-ligand complexes led to the identification of a putative allosteric binding site, which upon binding of PFOS stabilizes an active conformation of gmPpara1. Notably, binary exposures of gmPpara1, gmPpara2, and gmPparb to model-agonists and PFAAs produced similar potentiating effects. This study provides novel mechanistic insights into how PFAAs may modulate the Ppar signaling pathway by either binding the canonical ligand-binding pocket or by interacting with an allosteric binding site. Thus, individual PFAAs, or mixtures, could potentially modulate the Ppar-signaling pathway in Atlantic cod by interfering with at least one gmPpar subtype.

Computer simulations of Piezo1 channel regulation by membrane and small molecule

PIEZO proteins are homotrimeric cation-selective mechanosensitive ion channels that transduce various forms of mechanical forces into intracellular biochemical signals. Using multi-scale molecular dynamics simulations, we have generated an open state of the PIEOZ1 channel from the cryo-EM structure of a nonconducting conformation that reproduces unitary ionic conductance and selectivity in wild-type and mutants channels (Jiang, et al. Commun. Biol. 2021). This open state further allowed us to validate a putative allosteric agonist binding site of Yoda1 we previously identified using MD simulation and electrophysiology (Botello-Smith, et al.

Modeling and Structural Characterization of the Sweet Taste Receptor Heterodimer.

Sweet taste receptor, a heterodimer belonging to the class C G-protein coupled receptor (GPCR) family and composed of the T1R2 and T1R3 subunits, is responsible for the perception of natural sugars, sweet proteins, various d-amino acids, as well as artificial sweeteners. Despite the critical importance of the sweet receptor not only in mediating gustation but also in its role in the food industry, the architecture of the T1R2-T1R3 complex and the mechanism by which extracellular stimuli induce conformational changes that are propagated to the intracellular milieu, i.e., the signal transduction pathway, remain largely unknown. Here, we constructed and characterized a full-length structural model of the T1R2-T1R3 receptor, including both the transmembrane (TM) and extracellular (EC) domains of the heterodimer, using comparative modeling and extensive all-atom molecular dynamics simulations. Several heterodimer interfaces were first examined for the TM domain, and conformational changes occurring at the intracellular side and associated with the receptor's activation were characterized. From the analysis on the simulated data, putative allosteric binding sites for ligands, ions, and cholesterol were proposed. Also, insights into the protein interface of the TM domain upon activation are provided. The EC domain of the heterodimer, including both the Venus flytrap and cysteine-rich domains, was also investigated. Several important intersubunit interactions located at regions responsible for the receptor's proper function were observed, which resemble those recently identified in other class C GPCR members. Integration of the results from the TM and EC domains facilitates the generation of a full-length T1R2-T1R3 receptor. These findings along with the full-length structural model of the T1R2-T1R3 receptor provide a structural framework that may assist in understanding the mechanistic details associated with the receptor activation process for the sweet T1R2-T1R3 receptor as well as other members of the same family.

Cannabidiol binding and negative allosteric modulation at the cannabinoid type 1 receptor in the presence of delta-9-tetrahydrocannabinol: An In Silico study.

Recent evidence has raised in discussion the possibility that cannabidiol can act as a negative allosteric modulator of the cannabinoid type 1 receptor. Here we have used computational methods to study the modulation exerted by cannabidiol on the effects of delta-9-tetrahydrocannabinol in the cannabinoid receptor type 1 and the possibility of direct receptor blockade. We propose a putative allosteric binding site that is located in the N-terminal region of receptor, partially overlapping the orthosteric binding site. Molecular dynamics simulations reveled a coordinated movement involving the outward rotation of helixes 1 and 2 and subsequent expansion of the orthosteric binding site upon cannabidiol binding. Finally, changes in the cytoplasmic region and high helix 8 mobility were related to impaired receptor internalization. Together, these results offer a possible explanation to how cannabidiol can directly modulate effects of delta-9-tetrahydrocannabinol on the cannabinoid receptor type 1.

Could the presence of sodium ion influence the accuracy and precision of the ligand-posing in the human A2A adenosine receptor orthosteric binding site using a molecular docking approach? Insights from Dockbench.

The allosteric modulation of G protein-coupled receptors (GPCRs) by sodium ions has received considerable attention as crystal structures of several receptors, in their inactive conformation, show a Na+ ion bound to specific residues which, in the human A2A adenosine receptor (hA2A AR), are Ser913.39, Trp2466.48, Asn2807.45, and Asn2847.49. A cluster of water molecules completes the coordination of the sodium ion in the putative allosteric site. It is absolutely consolidated that the progress made in the field of GPCRs structural determination has increased the adoption of docking-driven approaches for the identification or the optimization of novel potent and selective ligands. Despite the extensive use of docking protocols in virtual screening approaches, to date, almost any of these studies have been carried out without taking into account the presence of the sodium cation and its first solvation shell in the putative allosteric binding site. In this study, we have focused our attention on determining how the presence of sodium ion binding and additionally its first hydration sphere, in hA2AAR could influence the ligand positioning accuracy during molecular docking simulations for most of the available resting and activated hA2A AR crystal structures, using DockBench as a comparative benchmarking tool and implementing a new correlation coefficient (EM). This work provides indications on the evidence that the posing performance (accuracy and/or precision) of the docking protocols in reproducing the crystallographic poses of different hA2A AR antagonists is generally increased in the presence of the sodium cation and its first solvation shell, in agreement with experimental observations. Consequently, the inclusion of sodium ion and its first solvation shell should be considered in order to facilitate the selection of new potential ligands in all molecular docking-based virtual screening protocols that aim to find novel GPCRs antagonists and inverse agonists.

Insights into the Allosteric Mechanism of Setmelanotide (RM-493) as a Potent and First-in-Class Melanocortin-4 Receptor (MC4R) Agonist To Treat Rare Genetic Disorders of Obesity through an in Silico Approach.

Human melanocortin-4 receptor (hMC4R) mutations have been implicated as the cause for about 6-8% of all severe obesity cases. Drug-like molecules that are able to rescue the functional activity of mutated receptors are highly desirable to combat genetic obesity among this population of patients. One such molecule is the selective MC4R agonist RM-493 (setmelanotide). While this molecule has been shown to activate mutated receptors with 20-fold higher potency over the endogenous agonist, little is known about its binding mode and how it effectively interacts with hMC4R despite the presence of mutations. In this study, a MC4R homology model was constructed based on the X-ray crystal structure of the adenosine A2A receptor in the active state. Four MC4R mutations commonly found in genetically obese patients and known to effect ligand binding in vitro were introduced into the constructed model. RM-493 was then docked into the wild-type and mutated models in order to better elucidate the possible binding modes for this promising drug candidate and assess how it may be interacting with MC4R to effectively activate receptor polymorphisms. The results reflected the orthosteric interactions of both the endogenous and synthetic ligands with the MC4R, which is supported by the site-directed mutagenesis studies. Meanwhile it helped explain the decremental affinity and potency of these ligands with the receptor polymorphisms. More significantly, our findings indicated that the structural characteristics of RM-493 may allow for enhanced receptor-ligand interactions, particularly through those with the putative allosteric binding sites, which facilitated the ligand to stabilize the active state of native and mutant MC4Rs to maintain reasonably high affinity and potency.

In silico mapping of allosteric ligand binding sites in type-1 cannabinoid receptor.

The recent resolution of the crystal structure of type-1 cannabinoid receptor (CB1 ) and the discovery of novel modulators for this target open the way to the possibility of elucidating the structural requirements for CB1 binding, and thereby facilitate a rational drug design. Compounds that target the orthosteric site of CB1 in some cases have shown side effects. Allosteric modulators could potentially avoid these side effects by influencing binding and/or efficacy of orthosteric ligands. Here, we summarize and compare previous data on different putative allosteric binding sites observed in CB1 homology models with an in silico docking study of the recently published crystal structure of the same receptor on endogenous and natural hydrophobic ligands that act as positive allosteric modulators and negative allosteric modulators of CB1 . In particular, a lipid-exposed pocket targeted by most of the tested molecules is reported and discussed.

Molecular tools for GABAA receptors: High affinity ligands for \u03b21-containing subtypes

γ-Aminobutyric acid type A (GABAA) receptors are pentameric GABA-gated chloride channels that are, in mammalians, drawn from a repertoire of 19 different genes, namely α1-6, β1-3, γ1-3, δ, ε, θ, π and ρ1-3. The existence of this wide variety of subunits as well as their diverse assembly into different subunit compositions result in miscellaneous receptor subtypes. In combination with the large number of known and putative allosteric binding sites, this leads to a highly complex pharmacology. Recently, a novel binding site at extracellular α+/β− interfaces was described as the site of modulatory action of several pyrazoloquinolinones. In this study we report a highly potent ligand from this class of compounds with pronounced β1-selectivity that mainly lacks α-subunit selectivity. It constitutes the most potent β1-selective positive allosteric modulatory ligand with known binding site. In addition, a proof of concept pyrazoloquinolinone ligand lacking the additional high affinity interaction with the benzodiazepine binding site is presented. Ultimately, such ligands can be used as invaluable molecular tools for the detection of β1-containing receptor subtypes and the investigation of their abundance and distribution.

New Anti‐Seizure (Arylalkyl)azole Derivatives: Synthesis, In Vivo and In Silico Studies

(Arylalkyl)azoles are a class of antiepileptic compounds including nafimidone, denzimol, and loreclezole (LRZ). Nafimidone and denzimol are thought to inhibit voltage-gated sodium channels (VGSCs) and enhance γ-aminobutyric acid (GABA)-mediated response. LRZ, a positive allosteric modulator of A-type GABA receptors (GABAA Rs), was reported to be sensitive to Asn265 of the β2/β3 subunit. Here, we report new N-[1-(4-chlorophenyl)-2-(1H-imidazol-1-yl)ethylidene]hydroxylamine esters showing anticonvulsant activity in animal models, including the 6-Hz psychomotor seizure test, a model for therapy-resistant partial seizure. We performed molecular docking studies for our active compounds using GABAA R and VGSC homology models. They predicted high affinity to the benzodiazepine binding site of GABAA R in line with the experimental results. Also, the binding mode and interactions of LRZ in its putative allosteric binding site of GABAA R is elucidated.

ATP Is an Allosteric Inhibitor of Coxsackievirus B3 Polymerase.

The RNA-dependent RNA polymerases from positive-strand RNA viruses, such as picornaviruses and flaviviruses, close their active sites for catalysis via a unique NTP-induced conformational change in the palm domain. Combined with a fully prepositioned templating nucleotide, this mechanism is error-prone and results in a distribution of random mutations in the viral progeny often described as a quasi-species. Here we examine the extent to which noncognate NTPs competitively inhibit single-cycle elongation by coxsackievirus B3 3D(pol), a polymerase that generates three to four mutations per 10 kb of RNA synthesized during viral infection. Using an RNA with a templating guanosine combined with 2-aminopurine fluorescence as a reporter for elongation, we find that the cognate CTP has a Km of 24 μM and the three noncognate nucleotides competitively inhibit the reaction with Kic values of 500 μM for GTP, 1300 μM for ATP, and 3000 μM for UTP. Unexpectedly, ATP also acted as an uncompetitive inhibitor with a Kiu of 1800 μM, resulting in allosteric modulation of 3D(pol) that slowed the polymerase elongation rate ≈4-fold. ATP uncompetitive inhibition required the β- and γ-phosphates, and its extent was significantly diminished in two previously characterized low-fidelity polymerases. This led to further mutational analysis and the identification of a putative allosteric binding site below the NTP entry channel at the interface of conserved motifs A and D, although cocrystallization failed to reveal any density for bound ATP in this pocket. The potential role of an ATP allosteric effect during the virus life cycle is discussed.

Benserazide, the first allosteric inhibitor of Coxsackievirus B3 3C protease

Coxsackievirus B3 is the main cause of human viral myocarditis and cardiomyopathy. Virally encoded Coxsackievirus 3C protease (3Cpro) plays an essential role in viral proliferation. Here, benserazide was discovered as a novel inhibitor from a drug library screen targeting Coxsackievirus 3Cpro using a FRET-based enzyme assay. Benserazide, whose chemical structure has no electrophilic functional groups, was characterized as a non-competitive inhibitor by enzyme kinetic studies. A molecular docking study with benserazide and its analogs indicated that a novel putative allosteric binding site was involved. Specifically, a 2,3,4-trihydroxybenzyl moiety was determined to be a key pharmacophore for the enzyme’s inhibitory activity. We suggest that the putative allosteric binding site may be a novel target for future therapeutic strategies.

The interaction of escitalopram and R-citalopram at the human serotonin transporter investigated in the mouse.

Escitalopram appears to be a superior antidepressant to racemic citalopram. It has been hypothesized that binding of R-citalopram to the serotonin transporter (SERT) antagonizes escitalopram binding to and inhibition of the SERT, there by curtailing the elevation of extracellular 5-hydroxytryptamine (5-HTExt), and hence anti-depressant efficacy. Further, it has been suggested that a putative allosteric binding site is important for binding of escitalopram to the primary, orthosteric, site, and for R-citalopram's inhibition here of. Primary: Investigate at the human (h)SERT, at clinical relevant doses, whether R-citalopram antagonizes escitalopram-induced 5-HTExt elevation. Secondary: Investigate whether abolishing the putative allosteric site affects escitalopram-induced 5-HTExt elevation and/or modulates the effect of R-citalopram. Recombinant generation of hSERT transgenic mice; in vivo microdialysis; SERT binding; pharmacokinetics; 5-HT sensitive behaviors (tail suspension, marble burying). We generated mice expressing either the wild-type human SERT (hSERT(WT)) or hSERT carrying amino acid substitutions (A505V, L506F, I507L, S574T and I575T) collectively abolishing the putative allosteric site (hSERT(ALI/VFL+SI/TT)). One mg/kg escitalopram yielded clinical relevant plasma levels and brain levels consistent with therapeutic SERT occupancy. The hSERT mice showed normal basal 5-HTExt levels. Escitalopram-induced 5-HTExt elevation was not decreased by R-citalopram co-treatment and was unaffected by loss of the allosteric site. The behavioral effects of the clinically relevant escitalopram dose were small and tended to be enhanced by R-citalopram co-administration. We find no evidence that R-citalopram directly antagonizes escitalopram or that the putative allosteric site is important for hSERT inhibition by escitalopram.

Toward a Global View of the Conformational Landscape of the Human Kinome

The human genome contains about 500 protein kinases, which play a central role in the regulation of the majority of cellular pathways. Mutations in kinase genes - often resulting in dysregulation of their phosphotransferase activity - are a frequent cause of disease, including many types of cancer. Kinases are especially flexible proteins, and undergo significant conformational changes during their catalytic and regulatory cycles. This conformational heterogeneity is also of fundamental importance in determining the binding affinity and selectivity of inhibitory ligands, their susceptibility to resistance mutations, and the discovery of putative allosteric binding sites. Markov state models have recently emerged as a practical computational approach to the enumeration of protein conformational states, and can be constructed by aggregating the data from multiple, independent, short molecular dynamics trajectories in a statistical fashion. We aim to apply this technique to the entire human kinome, simulating each protein kinase catalytic domain using a range of high performance compute resources, including the distributed simulation framework, Folding@Home. In combination with recent developments in GPU-accelerated simulation algorithms, this approach allows us to obtain aggregate trajectory lengths on the order of milliseconds. An automated software pipeline provides the ability to quickly generate multiple starting configurations for each kinase, while a central database of publicly available kinase data has been set up and used for tasks such as the selection of catalytic domain sequences and the assignment of relative priorities to each kinase. In parallel with our computational approach, we are working towards expressing a diverse range of kinases in bacterial systems, and scaling up a fluorescence-based assay to plate format for direct measurement of kinase inhibitor binding affinities. Our poster will present preliminary results from these efforts.

Assignments of backbone 1H, 13C and 15N resonances in H-Ras (1–166) complexed with GppNHp at physiological pH

The small GTPase Ras is an important signaling molecule acting as a molecular switch in eukaryotic cells. Recent findings of global conformational exchange and a putative allosteric binding site in the G domain of Ras opened an avenue to understanding novel aspects of Ras function. To facilitate detailed NMR studies of Ras in physiological solution conditions, we performed backbone resonance assignments of Ras bound to slowly hydrolysable GTP mimic, guanosine 5'-[ß, γ-imido]triphosphate at pH 7.2. Out of 163 non-proline residues of the G domain, signals from backbone amide proton, nitrogen and carbon spins of 127 residues were confidently assigned with the remaining unassigned residues mostly located at the exchange-broadened effectors interface.

The molecular binding interactions of inhibitors and activators of phosphoenolpyruvate carboxylase

We have performed molecular modelling studies of the binding to maize phosphoenolpyruvate carboxylase (PEPC) of phosphoenolpyruvate (PEP) and a number of representative competitive inhibitors. We predict that all these compounds share a common binding mode and that the differences in inhibitory activity of the various inhibitors arise mainly from either increased hydrophobic interactions of cis substituents or small but significant changes in their binding mode arising from steric clashes of trans substituents with the active site. We have also performed molecular modelling studies of glucose-6-phosphate (G6P) and a number of other allosteric activators in their putative allosteric binding site in this enzyme. We predict that these molecules share the same binding mode for their phosphate moiety while some of them engage in a variety of hydrogen bonding interactions with residues from different subunits of the enzyme, and others establish hydrophobic and van der Waals interactions with other regions of the allosteric binding site.